CN110555602A - 一种径流资料短缺情况下的旬径流分配方法和*** - Google Patents
一种径流资料短缺情况下的旬径流分配方法和*** Download PDFInfo
- Publication number
- CN110555602A CN110555602A CN201910759536.4A CN201910759536A CN110555602A CN 110555602 A CN110555602 A CN 110555602A CN 201910759536 A CN201910759536 A CN 201910759536A CN 110555602 A CN110555602 A CN 110555602A
- Authority
- CN
- China
- Prior art keywords
- runoff
- day
- data
- division
- distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000011156 evaluation Methods 0.000 claims abstract description 23
- 238000011160 research Methods 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 abstract description 5
- 239000013589 supplement Substances 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种径流资料短缺情况下的旬径流分配方法和***,属于径流资料短缺或无资料地区水文研究领域,所述方法包括:根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分;根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则继续进行划分。本发明方法能够充分利用已有的月径流(年径流资料)资料,对更小时间单位的径流数据进行补充,有效减少现有分配方法的固有误差,无需参数限制,使分配结果更加贴合实际情况,可以应用于无资料水文研究领域中解决相关径流模拟分配问题,简单高效、便于实施,具有广泛的应用价值。
Description
技术领域
本发明属于径流资料短缺或无资料地区水文研究领域,更具体地,涉及一种径流资料短缺情况下的旬径流分配方法和***。
背景技术
无资料或资料短缺地区的水文研究是近代国际水文水资源研究的热点和难点问题之一,就我国而言,很多水库、水电站建设之初都没有***的径流量数据记载,大多数径流资料都是以月作为最小单位。但在实际工程应用或者科学研究中,往往需要更为详细的数据支撑,如旬、日数据,此时便需要将月径流资料进行科学合理的划分。
但目前的划分方法比较粗略,常见的方法是将月径流均匀地分到每一旬或者每一日(等分),这种方法虽然简单,但存在着很明显的不足,地表径流的形成受到气候、土地利用类型等多种因素的影响,因此径流的时间序列,以及降雨-径流的关系通常具有明显的非平稳性,均等划分使旬与旬(日与日)间的径流值跳跃常常过大或者跳跃为零,与实际相悖。
当然也有很多用随机模拟的方法生成随机数代替实际径流数据,但涉及到对统计特征的模拟、随机模型的建立(含定阶、参数估计等),需要考虑到趋势项、跳跃项、周期项、随机项等一系列因素,步骤上相对复杂,增加了计算负担,没有充分利用到已有的月径流资料。
总体而言,现有径流分配方法存在分配旬径流量时有明显误差,使得模拟数据与实际情况不符的问题。
发明内容
针对现有技术的缺陷,本发明的目的在于提供一种径流资料短缺情况下的旬径流分配方法和***,旨在解决采用现有方法分配旬径流量时有明显误差,使得模拟数据与实际情况不符的问题。
为实现上述目的,本发明一方面提供了一种径流资料短缺情况下的旬径流分配方法,包括:
(1)根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分;
(2)根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则返回执行步骤(1)。
进一步地,步骤(1)具体为,采用大规模线性规划求解器根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分。
进一步地,所述二次规划模型为:
其中,i表示旬径流编号,N表示径流序列中所划分旬的个数,Qi为时段i的旬径流量。
进一步地,所述二次规划模型满足以下约束条件:
总量平衡约束:Q3j+Q3j+1+Q3j+2=3qj,其中,j表示月径流编号,qj为时段j的月径流量;非负约束:Qi≥0。
进一步地,所述评价参数为:
其中,Cvm为月径流序列离差系数,Cvx为优化分配后旬径流序列离差系数,K为评价参数,反映旬径流序列与月径流序列离差参数的变化量。
本发明另一方面提供了一种径流资料短缺情况下的旬径流分配***,包括:
旬径流分配模块,用于根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分;
分配结果评价模块,用于根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则继续采用旬径流分配模块进行旬径流分配。
进一步地,所述旬径流分配模块采用大规模线性规划求解器根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分。
进一步地,所述二次规划模型为:
其中,i表示旬径流编号,N表示径流序列中所划分旬的个数,Qi为时段i的旬径流量。
进一步地,所述二次规划模型满足以下约束条件:
总量平衡约束:Q3j+Q3j+1+Q3j+2=3qj,其中,j表示月径流编号,qj为时段j的月径流量;非负约束:Qi≥0。
进一步地,所述评价参数为:
其中,Cvm为月径流序列离差系数,Cvx为优化分配后旬径流序列离差系数,K为评价参数,反映旬径流序列与月径流序列离差参数的变化量。
通过本发明所构思的以上技术方案,与现有技术相比,能够取得以下有益效果:
(1)相比现有无资料情形的年内逐旬径流分配方法,本发明方法能够充分利用已有的月径流(年径流资料)资料,对更小时间单位的径流数据进行补充,有效减少了现有分配方法的固有误差,无需参数限制,使分配结果更加贴合实际情况。
(2)本发明通过调用大规模线性规划求解器,能够快速准确地求解得到二次规划模型的全局最优解,从而实现年内逐旬流量的最优分配,简单高效、便于实施,具有广泛的应用价值。
附图说明
图1为本发明提出的一种径流资料短缺情况下的旬径流分配方法流程图;
图2(a)为应用本发明方法得到的逐旬流量序列分配示意图;
图2(b)为应用现有分配方法得到的逐旬流量序列分配示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
参考图1,本发明实施例一方面提供了一种径流资料短缺情况下的旬径流分配方法,包括:
(1)根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分。
具体地,根据已知的月径流或者部分已知的旬径流资料,以旬流量序列中相邻旬流量的波动值的平方和最小为目标建立数学规划模型,由于线性规划求解器gurobi采用先进的优化技术,充分利用多核处理器的优势,实现高效稳定的求解,并且提供方便轻巧的接口,便于直接调用,因此,本发明采用大规模线性规划求解器gurobi对已知的径流资料进行以旬为单位的径流划分,为使数学规划模型适应大规模线性规划求解器以求得最优解,将数学规划模型重构为严格的二次规划模型:其中,i表示旬径流编号,N表示径流序列中所划分旬的个数,Qi为时段i的旬径流量。本发明中二次规划模型满足以下约束条件:总量平衡约束:Q3j+Q3j+1+Q3j+2=3qj,其中,j表示月径流编号,qj为时段j的月径流量;非负约束:Qi≥0。
(2)根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则返回执行步骤(1)。
具体地,本发明对分配得到的旬径流量评价参数为:
其中,Cvm为月径流序列离差系数,Cvx为优化分配后旬径流序列离差系数,K为评价参数,反映旬径流序列与月径流序列离差参数的变化量,K值越小,说明分配后的旬径流序列与之前的月径流序列越相似,根据先前的分配经验,一般认为当K<0.01时,旬径流序列分配与月径流相似度满足要求。
为了验证本发明方法的有效性,选取某资料短缺地区水文站的某年径流序列进行分配,图2(a)、图2(b)分别为采用本发明分配方法及现有分配方法所产生的分配序列,可以看出,图2(a)的旬径流序列相比图2(b)而言,其径流量波动更加平稳,波动趋势更加符合实际,其中,图2(a)K值为0.06%,图2(b)K值为2.90%,旬径流量序列分配的相似度相比之下提升近两个数量级,且满足评价标准的要求,进一步佐证本发明分配方法所生成的旬径流序列更具有代表性及实际意义。
本发明实施例另一方面提供了一种径流资料短缺情况下的旬径流分配***,包括:旬径流分配模块,用于根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分;分配结果评价模块,用于根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则继续采用旬径流分配模块进行旬径流分配。本发明旬径流分配模块采用大规模线性规划求解器根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分。采用的二次规划模型为:其中,i表示旬径流编号,N表示径流序列中所划分旬的个数,Qi为时段i的旬径流量。二次规划模型满足以下约束条件:总量平衡约束:Q3j+Q3j+1+Q3j+2=3qj,其中,j表示月径流编号,qj为时段j的月径流量;非负约束:Qi≥0。采用的评价参数为:其中,Cvm为月径流序列离差系数,Cvx为优化分配后旬径流序列离差系数,K为评价参数,反映旬径流序列与月径流序列离差参数的变化量。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种径流资料短缺情况下的旬径流分配方法,其特征在于,包括:
(1)根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分;
(2)根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则返回执行步骤(1)。
2.根据权利要求1所述的一种径流资料短缺情况下的旬径流分配方法,其特征在于,步骤(1)具体为,采用大规模线性规划求解器根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分。
3.根据权利要求1或2所述的一种径流资料短缺情况下的旬径流分配方法,其特征在于,所述二次规划模型为:
其中,i表示旬径流编号,N表示径流序列中所划分旬的个数,Qi为时段i的旬径流量。
4.根据权利要求3所述的一种径流资料短缺情况下的旬径流分配方法,其特征在于,所述二次规划模型满足以下约束条件:
总量平衡约束:Q3j+Q3j+1+Q3j+2=3qj,其中,j表示月径流编号,qj为时段j的月径流量;非负约束:Qi≥0。
5.根据权利要求1-4任一项所述的一种径流资料短缺情况下的旬径流分配方法,其特征在于,所述评价参数为:
其中,Cvm为月径流序列离差系数,Cvx为优化分配后旬径流序列离差系数,K为评价参数,反映旬径流序列与月径流序列离差参数的变化量。
6.一种径流资料短缺情况下的旬径流分配***,其特征在于,包括:
旬径流分配模块,用于根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分;
分配结果评价模块,用于根据评价参数判断划分得到的旬径流量是否达到设定标准,如是,则划分结束;如否,则继续采用旬径流分配模块进行旬径流分配。
7.根据权利要求6所述的一种径流资料短缺情况下的旬径流分配***,其特征在于,所述旬径流分配模块采用大规模线性规划求解器根据构建的二次规划模型对已知的径流资料进行以旬为单位的径流划分。
8.根据权利要求6或7所述的一种径流资料短缺情况下的旬径流分配***,其特征在于,所述二次规划模型为:
其中,i表示旬径流编号,N表示径流序列中所划分旬的个数,Qi为时段i的旬径流量。
9.根据权利要求8所述的一种径流资料短缺情况下的旬径流分配***,其特征在于,所述二次规划模型满足以下约束条件:
总量平衡约束:Q3j+Q3j+1+Q3j+2=3qj,其中,j表示月径流编号,qj为时段j的月径流量;非负约束:Qi≥0。
10.根据权利要求1-9任一项所述的一种径流资料短缺情况下的旬径流分配***,其特征在于,所述评价参数为:
其中,Cvm为月径流序列离差系数,Cvx为优化分配后旬径流序列离差系数,K为评价参数,反映旬径流序列与月径流序列离差参数的变化量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910759536.4A CN110555602B (zh) | 2019-08-16 | 2019-08-16 | 一种径流资料短缺情况下的旬径流分配方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910759536.4A CN110555602B (zh) | 2019-08-16 | 2019-08-16 | 一种径流资料短缺情况下的旬径流分配方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110555602A true CN110555602A (zh) | 2019-12-10 |
CN110555602B CN110555602B (zh) | 2023-04-07 |
Family
ID=68737741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910759536.4A Active CN110555602B (zh) | 2019-08-16 | 2019-08-16 | 一种径流资料短缺情况下的旬径流分配方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110555602B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103577895A (zh) * | 2013-11-08 | 2014-02-12 | 河海大学 | 一种资料短缺情形下的二次耦合月径流预报方法 |
US20180128939A1 (en) * | 2016-11-09 | 2018-05-10 | Guizhou Normal University | Adjustable karst underground water and soil loss simulation apparatus |
CN108053083A (zh) * | 2018-01-16 | 2018-05-18 | 河南创辉水利水电工程有限公司 | 一种水库水电站非汛期联合优化发电调度方法 |
-
2019
- 2019-08-16 CN CN201910759536.4A patent/CN110555602B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103577895A (zh) * | 2013-11-08 | 2014-02-12 | 河海大学 | 一种资料短缺情形下的二次耦合月径流预报方法 |
US20180128939A1 (en) * | 2016-11-09 | 2018-05-10 | Guizhou Normal University | Adjustable karst underground water and soil loss simulation apparatus |
CN108053083A (zh) * | 2018-01-16 | 2018-05-18 | 河南创辉水利水电工程有限公司 | 一种水库水电站非汛期联合优化发电调度方法 |
Non-Patent Citations (6)
Title |
---|
周婷等: "基于小波分解的优化支持向量机模型在水库年径流预测中的应用", 《华北水利水电大学学报(自然科学版)》 * |
徐长江等: "资料短缺地区径流模拟研究" * |
洪家宁等: "反褶积模型在日径流模拟中的应用" * |
罗吉忠等: "基于SWAT模型的缺资料流域径流模拟研究", 《西南民族大学学报(自然科学版)》 * |
胡彩虹等: "汾河水库上游流域径流年内分配时程变化规律研究", 《干旱区资源与环境》 * |
靳晟等: "设计年径流及其年内分配计算软件的开发研究", 《人民黄河》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110555602B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104881582B (zh) | 空气质量的预测方法和装置 | |
CN103186566A (zh) | 一种数据分级存储方法、装置及*** | |
CN110501733B (zh) | 一种自适应格网vrs生成与服务方法 | |
CN104036327A (zh) | 一种用于智能配电网的快速母线负荷预测方法 | |
CN109766349A (zh) | 任务防重方法、装置、计算机设备及存储介质 | |
Veselov et al. | Methodological approach for harmonization of the investment and pricing policy options in the electric power industry | |
CN104951903A (zh) | 一种基于金融资产管理***架构的自定义资产评估方法 | |
CN110555602B (zh) | 一种径流资料短缺情况下的旬径流分配方法和*** | |
CN113902160A (zh) | 一种应用于水利枢纽工程的水调方法 | |
CN103246539B (zh) | 基于用户使用习惯的gis***模块加载方法 | |
Zagona et al. | RiverWare | |
CN113849884A (zh) | 一种基于改进长系列法的水库多级调度图编制方法 | |
CN111090710B (zh) | 一种配网多时间尺度图模异动自动维护方法 | |
CN112965665B (zh) | 一种基于sas和ssd的gp数据库数据存储方法 | |
CN109872250B (zh) | 一种电力市场初期基数偏差电量结算处理方法及*** | |
CN112036604A (zh) | 一种考虑多时序过程因子的中期径流预报方法 | |
CN108599976B (zh) | 网络质量评估模型构建方法、装置及计算机设备 | |
CN111080115A (zh) | 一种基于供水水力模型的用水量分配方法及装置 | |
CN111596908B (zh) | 一种新安江模型云计算方法 | |
CN115953251B (zh) | 一种保险佣金核算方法和*** | |
CN110162839A (zh) | 一种流域水沙变化影响因素贡献率的辨识方法及*** | |
CN115965056A (zh) | 一种中长期电量变权组合预测方法及*** | |
CN115269564B (zh) | 一种用于大规模***的集中信创化迁移方法 | |
CN103559266A (zh) | 多模式匹配方法及装置 | |
Chunye | Research on Measuring Method of Function Size of Information System |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230506 Address after: No.73 Tuodong Road, Guandu District, Kunming, Yunnan 650011 Patentee after: KUNMING ELECTRIC POWER TRANSACTION CENTER CO.,LTD. Address before: No. 73 Tuodong Road, Panlong District, Kunming City, Yunnan Province, 650011 Patentee before: KUNMING ELECTRIC POWER TRANSACTION CENTER CO.,LTD. Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY |