CN113326632A - 一种反推水库入库流量过程的优化修正方法 - Google Patents

一种反推水库入库流量过程的优化修正方法 Download PDF

Info

Publication number
CN113326632A
CN113326632A CN202110681901.1A CN202110681901A CN113326632A CN 113326632 A CN113326632 A CN 113326632A CN 202110681901 A CN202110681901 A CN 202110681901A CN 113326632 A CN113326632 A CN 113326632A
Authority
CN
China
Prior art keywords
water level
corrected
reservoir
flow process
flow
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
Application number
CN202110681901.1A
Other languages
English (en)
Other versions
CN113326632B (zh
Inventor
康传雄
刘轩萁
廖浩森
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang Institute of Technology
Original Assignee
Nanchang Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanchang Institute of Technology filed Critical Nanchang Institute of Technology
Priority to CN202110681901.1A priority Critical patent/CN113326632B/zh
Publication of CN113326632A publication Critical patent/CN113326632A/zh
Application granted granted Critical
Publication of CN113326632B publication Critical patent/CN113326632B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Feedback Control In General (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种反推水库入库流量过程的优化修正方法，本发明提供的方法允许在一定范围内修正观测水位值，以修正后的入库流量过程变幅最小为目标，建立优化修正模型，然后调用开源混合整数线性规划解算器Cbc进行求解，输出修正后的入库流量过程。本发明提供的优化修正模型物理意义明确，修正后的入库流量过程平滑，又能接近真实的天然入流情况，具有很好的实际应用价值。

Description

一种反推水库入库流量过程的优化修正方法
技术领域
本发明属于水库调度和水文测验领域，具体涉及一种反推水库入库流量过程的优化修正方法。
背景技术
水库入库流量是水库资料整编的一部分，对评估和提高水库运行效益有重要意义。对于没有入库流量观测的水库，根据水量平衡方程反推入库流量过程是常用的做法，即根据水库水位观测计算水库蓄水量变化，加上出库流量得到入库流量。该方法在实际运用中推算的入库流量过程常呈现“锯齿形”，甚至出现负值。即便计入蒸发和渗漏等水量损失，也不能改善上述问题。出现上述问题的根本原因是：受观测设备的精度、动库容、闸门启闭、风浪等因素的影响，观测的水位值不能准确反应水库的实际蓄水量。对于大中型水库，即便是0.01m的水位差，其对应的库容也是巨大的，因此导致直接反推的入库流量过程不合理。
发明内容
本发明针对直接反推入库流量过程不合理这一实际问题，基于观测水位值不能准确反应水库实际蓄水量这一客观原因，提供一种反推水库入库流量过程的优化修正方法，具体采用以下的技术方案：
一种反推水库入库流量过程的优化修正方法，其特征在于，包括以下步骤：将水库水位监测值和出库流量监测值输入优化修正模型，调用混合整数线性规划解算器Cbc求解后输出修正后的入库流量过程；
所述优化修正模型的目标函数为
Figure BDA0003122982370000011
其中，
Figure BDA0003122982370000012
和
Figure BDA0003122982370000013
分别表示修正后的相邻时段入库流量的增加量或减少量，在最小目标下，二者至少一个为0；
所述优化修正模型的约束条件包括：相邻时段的流量与变幅关系、允许的水位修正范围、水量平衡关系、水位-库容关系和非负约束；
所述相邻时段的流量与变幅关系具体为
Figure BDA0003122982370000014
其中，Qt′是修正后的时段平均入库流量；
所述允许的水位修正范围具体为-△z≤Zt′-Zt≤△z，其中，Zt和Zt′分别是修正前后的时刻水位，△z表示允许的修正范围，为0.01m～0.02m；
所述水量平衡关系具体为Vt′+1＝Vt′+(Qt′-qt)·△t-△wt，其中，Vt′表示修正后的水位Zt′对应的库容值，qt表示时段平均出库流量，△t表示时段时间长，△wt表示时段水量损失；
所述水位-库容关系具体为Vt′＝fZV(Zt′)，其中fZV表示水库的水位-库容曲线，以离散的水位-库容值表示，采用分段线性建模为混合整数线性规划模型；
所述非负约束具体为上述各变量均取非负值。
其中，△wt表示时段水量损失(实际中可忽略不计)。
本发明的有益效果为：本发明提供的优化修正模型物理意义明确，修正后的入库流量过程平滑，又能接近真实的天然入流情况，具有很好的实际应用价值。
附图说明
图1所示为入库流量优化修正方法的步骤图；
图2所示为优化修正前后的72个时段的三峡水库入库流量数据图。
具体实施方式
以下将结合实施例和附图对本发明的构思及产生的技术效果进行清楚、完整的描述，以充分地理解本发明的目的、方案和效果。
实施例1：
选取三峡水库2021/3/20 08:00～2021/3/23 08:00共72个时段(小时)的水位监测值和出库流量监测值，根据水量平衡关系直接反推的入库流量过程如图2实线所示，呈现明显的“锯齿形”，并且出现了5次负值，极为不合理。
运用本发明提供的方法修正入库流量过程的具体步骤如下(如图1所示)：
步骤一、导入2021/3/20 08:00～2021/3/23 08:00共72个时段的三峡水库数据，包括：水库水位监测值Zt和出库流量监测值qt。
步骤二、根据目标函数和约束条件建立优化修正模型，目标函数为
Figure BDA0003122982370000021
其中，
Figure BDA0003122982370000022
和
Figure BDA0003122982370000023
分别表示修正后的相邻时段入库流量的增加量或减少量，在最小目标下，二者至少一个为0；
约束条件包括：相邻时段的流量与变幅关系、允许的水位修正范围、水量平衡关系、水位-库容关系和非负约束；
相邻时段的流量与变幅关系具体为
Figure BDA0003122982370000024
其中，Qt′是修正后的时段平均入库流量；
允许的水位修正范围具体为-△z≤Zt′-Zt≤△z，其中，Zt和Zt′分别是修正前后的时刻水位，△z表示允许的修正范围，为0.02m；
水量平衡关系具体为Vt′+1＝Vt′+(Qt′-qt)·△t-△wt，其中，Vt′表示修正后的水位Zt′对应的库容值，qt表示时段平均出库流量，△t表示时段时间长，△wt表示时段水量损失；
水位-库容关系具体为Vt′＝fZV(Zt′)，其中fZV表示水库的水位-库容曲线，以离散的水位-库容值表示，采用分段线性建模为混合整数线性规划模型；
非负约束具体为上述各变量均取非负值；其中决策变量包括：修正后的水位Zt′和库容Vt′、修正后的入库流量Qt′、相邻时段的入库流量增加量
Figure BDA0003122982370000031
或减少量
Figure BDA0003122982370000032
步骤三、调用开源混合整数线性规划解算器Cbc求解优化修正模型，输出修正后的入库流量过程。
通过本实施例的方法修正得到的入库流量过程如图2虚线所示，该流量过程平滑，无不合理的极小值，变化趋势符合真实的入流情况。在无实测入库流量的情况下，运用本发明修正的入库流量过程可以作为水库的入库流量。
以上所述，只是本发明的一个实施例而已，本发明并不局限于上述实施方式，只要其以相同的手段达到本发明的技术效果，都应属于本发明的保护范围。在本发明的保护范围内其技术方案和/或实施方式可以有不同的修改和变化。

Claims (1)

1.一种反推水库入库流量过程的优化修正方法，其特征在于，包括以下步骤：将水库水位监测值和出库流量监测值输入优化修正模型，调用混合整数线性规划解算器Cbc求解后输出修正后的入库流量过程；
所述优化修正模型的目标函数为
Figure FDA0003122982360000011
其中，
Figure FDA0003122982360000012
和
Figure FDA0003122982360000013
分别表示修正后的相邻时段入库流量的增加量或减少量，在最小目标下，二者至少一个为0；
所述优化修正模型的约束条件包括：相邻时段的流量与变幅关系、允许的水位修正范围、水量平衡关系、水位-库容关系和非负约束；
所述相邻时段的流量与变幅关系具体为
Figure FDA0003122982360000014
其中，Q′t是修正后的时段平均入库流量；
所述允许的水位修正范围具体为-△z≤Z′t-Zt≤△z，其中，Zt和Z′t分别是修正前后的时刻水位，△z表示允许的修正范围，为0.01m～0.02m；
所述水量平衡关系具体为V′t+1＝Vt′+(Q′t-qt)·△t-△wt，其中，Vt′表示修正后的水位Z′t对应的库容值，qt表示时段平均出库流量，△t表示时段时间长，△wt表示时段水量损失；
所述水位-库容关系具体为Vt′＝fZV(Zt′)，其中fZV表示水库的水位-库容曲线，以离散的水位-库容值表示，采用分段线性建模为混合整数线性规划模型；
所述非负约束具体为上述各变量均取非负值。
CN202110681901.1A 2021-06-19 2021-06-19 一种反推水库入库流量过程的优化修正方法 Active CN113326632B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110681901.1A CN113326632B (zh) 2021-06-19 2021-06-19 一种反推水库入库流量过程的优化修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110681901.1A CN113326632B (zh) 2021-06-19 2021-06-19 一种反推水库入库流量过程的优化修正方法

Publications (2)

Publication Number Publication Date
CN113326632A true CN113326632A (zh) 2021-08-31
CN113326632B CN113326632B (zh) 2022-09-23

Family

ID=77423944

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110681901.1A Active CN113326632B (zh) 2021-06-19 2021-06-19 一种反推水库入库流量过程的优化修正方法

Country Status (1)

Country Link
CN (1) CN113326632B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024066902A1 (zh) * 2022-09-28 2024-04-04 中国长江三峡集团有限公司 入库径流修正优化方法及装置

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102155938A (zh) * 2011-04-07 2011-08-17 武汉大学 一种反推水库入库流量过程的测算方法
CN102968674A (zh) * 2012-12-05 2013-03-13 湖南省电力公司 一种基于水位控制的日调节水电站发电优化调度方法
CN103761583A (zh) * 2014-01-07 2014-04-30 中国长江三峡集团公司 一种水库排沙发电多目标优化调度方法
CN104021307A (zh) * 2014-06-20 2014-09-03 武汉大学 一种反推水库入库流量过程的测算方法
CN104123589A (zh) * 2014-06-24 2014-10-29 华中科技大学 一种梯级水电站短期优化调度方法
CN104268645A (zh) * 2014-09-25 2015-01-07 杨树滩 水资源配置中考虑供水优先序的行业供水量确定方法
CN105893672A (zh) * 2016-03-31 2016-08-24 中国水利水电科学研究院 一种狭长河道型水库全生命周期温度场研究方法
CN107180318A (zh) * 2017-07-28 2017-09-19 华中科技大学 基于变权重剩余防洪库容最大的水库群防洪库容分配方法
CN108681783A (zh) * 2018-04-04 2018-10-19 河海大学 一种水库实时多目标随机优化调度和风险评估方法
CN110188912A (zh) * 2019-02-28 2019-08-30 西安理工大学 基于改进花粉算法的地表水和地下水联合调度优化方法
CN111008790A (zh) * 2019-12-23 2020-04-14 华中科技大学 一种水电站群发电调度规则提取方法
CN111898253A (zh) * 2020-07-15 2020-11-06 武汉大学 一种水库调度及其下游河流生态环境保护的合作价值评估方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102155938A (zh) * 2011-04-07 2011-08-17 武汉大学 一种反推水库入库流量过程的测算方法
CN102968674A (zh) * 2012-12-05 2013-03-13 湖南省电力公司 一种基于水位控制的日调节水电站发电优化调度方法
CN103761583A (zh) * 2014-01-07 2014-04-30 中国长江三峡集团公司 一种水库排沙发电多目标优化调度方法
CN104021307A (zh) * 2014-06-20 2014-09-03 武汉大学 一种反推水库入库流量过程的测算方法
CN104123589A (zh) * 2014-06-24 2014-10-29 华中科技大学 一种梯级水电站短期优化调度方法
CN104268645A (zh) * 2014-09-25 2015-01-07 杨树滩 水资源配置中考虑供水优先序的行业供水量确定方法
CN105893672A (zh) * 2016-03-31 2016-08-24 中国水利水电科学研究院 一种狭长河道型水库全生命周期温度场研究方法
CN107180318A (zh) * 2017-07-28 2017-09-19 华中科技大学 基于变权重剩余防洪库容最大的水库群防洪库容分配方法
CN108681783A (zh) * 2018-04-04 2018-10-19 河海大学 一种水库实时多目标随机优化调度和风险评估方法
CN110188912A (zh) * 2019-02-28 2019-08-30 西安理工大学 基于改进花粉算法的地表水和地下水联合调度优化方法
CN111008790A (zh) * 2019-12-23 2020-04-14 华中科技大学 一种水电站群发电调度规则提取方法
CN111898253A (zh) * 2020-07-15 2020-11-06 武汉大学 一种水库调度及其下游河流生态环境保护的合作价值评估方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHUANXIONG KANG等: "Short-Term Hydrothermal Scheduling Using a Two-Stage Linear Programming with Special Ordered Sets Method", 《WATER RESOURCES MANAGEMENT》 *
TONGTIEGANG ZHAO等: "Optimizing Operation of Water Supply Reservoir: The Role of Constraints", 《MATHEMATICAL PROBLEMS IN ENGINEERING》 *
张睿: "流域大规模梯级电站群协同发电优化调度研究", 《中国博士学位论文电子期刊网 工程科技Ⅱ辑》 *
张高峰 等: "基于数学优化法绘制梯级水库水力剖面图的研究", 《水力发电》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024066902A1 (zh) * 2022-09-28 2024-04-04 中国长江三峡集团有限公司 入库径流修正优化方法及装置

Also Published As

Publication number Publication date
CN113326632B (zh) 2022-09-23

Similar Documents

Publication Publication Date Title
Liu et al. Deriving optimal refill rules for multi-purpose reservoir operation
CN107992989B (zh) 一种水库群洪水资源化预蓄预泄风险决策方法
Xu et al. Stochastic programming with a joint chance constraint model for reservoir refill operation considering flood risk
CN113326632B (zh) 一种反推水库入库流量过程的优化修正方法
Dittmann et al. Optimum multi-objective reservoir operation with emphasis on flood control and ecology
CN115099477B (zh) 一种水库旱限水位优化及抗旱调度方法
CN111898253B (zh) 水库调度及其下游河流生态环境保护的合作价值评估方法
Celeste et al. The role of spill and evaporation in reservoir optimization models
CN112149983A (zh) 一种耦合气象-水文不确定性的水库汛限水位动态控制风险分析方法
CN110223156A (zh) 基于逐步最优特征选择的自动化模型调优算法
CN115640881A (zh) 一种基于lstm的水库入库径流预报模型的水库水位修正方法
Wang et al. Scenario-based multi-objective optimization of reservoirs in silt-laden rivers: A case study in the Lower Yellow River
CN109706906A (zh) 一种三角门船闸开通闸运行的控制方法
CN110578317B (zh) 一种水文模型水库泄水量模拟方法
CN110647179B (zh) 一种水库实时防洪补偿调度的减压控制水位确定方法
Amirkhani et al. Multiobjective optimal operation of gated spillways
Sewagudde Lake Victoria’s water budget and the potential effects of climate change in the 21st century
CN110795688A (zh) 考虑遥相关因子的中长期水库调度方法以及自动控制***
CN115905776A (zh) 一种考虑上游水库调蓄影响的施工期设计洪水计算方法
CN112580230B (zh) 一种水库水温结构分析方法
Fisher et al. Modelling of hydrological persistence for hidden state Markov decision processes
Soler et al. Operation of an irrigation canal by means of the passive canal control
CN111754099A (zh) 一种基于三阶段风险对冲规则的实时防洪调度方法
JP5318081B2 (ja) 取水制御方法及び装置
Wu et al. Minimizing the impact of vacating instream storage of a multi-reservoir system: a trade-off study of water supply and empty flushing

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