CN110929929A - 一种考虑下游河道水力安全的水电站泄洪优化调控方法 - Google Patents

一种考虑下游河道水力安全的水电站泄洪优化调控方法 Download PDF

Info

Publication number
CN110929929A
CN110929929A CN201911136971.8A CN201911136971A CN110929929A CN 110929929 A CN110929929 A CN 110929929A CN 201911136971 A CN201911136971 A CN 201911136971A CN 110929929 A CN110929929 A CN 110929929A
Authority
CN
China
Prior art keywords
hydropower station
regulation
daily
flow
representing
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.)
Pending
Application number
CN201911136971.8A
Other languages
English (en)
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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201911136971.8A priority Critical patent/CN110929929A/zh
Publication of CN110929929A publication Critical patent/CN110929929A/zh
Pending legal-status Critical Current

Links

Images

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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067—Enterprise or organisation modelling
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use

Abstract

本发明属于水利水电工程的技术领域，具体涉及一种考虑下游河道水力安全的水电站泄洪优化调控方法，包括：步骤一，选取枢纽水电站下游的具有防冲安全问题的日调节水电站作为调控对象；步骤二，建立优化调控模型，并设置优化调控模型的目标函数和约束条件；步骤三，选取决策变量进行优化，获得优化结果。本发明显著地提高了日调节水电站的防冲安全度，通过上游枢纽水电站的调蓄库容削弱洪水流量峰值，并合理配置日调节水电站的下游水位与各个闸门的泄流量，有效地提高了下游易冲刷位置的水力安全度。

Description

一种考虑下游河道水力安全的水电站泄洪优化调控方法
技术领域
本发明属于水利水电工程的技术领域，具体涉及一种考虑下游河道水力安全的水电站泄洪优化调控方法。
背景技术
中低水头、大引用流量水电站下游是容易发生冲蚀淘刷破坏的重灾区。我国因为泄洪引发冲蚀淘刷问题造成泄水建筑物破坏的水电站比比皆是：如马迹塘水电站闸坝下游河床出现较大的冲刷坑，南津渡水电站在下游回流区发生了明显冲蚀破坏，长湖水电站下游两岸发生了严重的淘刷破坏。
申请人发现：导致这些问题的原因在于泄洪流量峰值高、波动大，使得下游水流无法正常消能；并且，闸门间开度组合不当，容易造成下游流速分布与流态不好。因此，亟需一种新型的水电站泄洪优化调控方法来解决上述问题。
发明内容
本发明的目的在于：针对现有技术的不足，提供一种考虑下游河道水力安全的水电站泄洪优化调控方法，提高了日调节水电站的防冲安全度，通过上游枢纽水电站的调蓄库容削弱洪水流量峰值，并合理配置日调节水电站的下游水位与各个闸门的泄流量，有效地提高了下游易冲刷位置的水力安全度。
为了实现上述目的，本发明采用如下技术方案：
一种考虑下游河道水力安全的水电站泄洪优化调控方法，包括：
步骤一：选取枢纽水电站下游的具有防冲安全问题的日调节水电站作为调控对象，获取所述枢纽水电站与所述日调节水电站的区间来流量
Figure BDA0002279833290000021
获取所述枢纽水电站调蓄前的所述日调节水电站的上游来流量
Figure BDA0002279833290000022
获取所述枢纽水电站调蓄前第t时段的下泄流量
Figure BDA0002279833290000023
步骤二：分别以所述日调节水电站的洪水流量峰值最小及以所述日调节水电站的闸门调控工作量最小为目标函数建立优化调控模型，并设置所述优化调控模型的约束条件；
步骤三：生成下泄流量序列Qzxxt,所述Qzxxt由若干个所述
Figure BDA0002279833290000024
组成，以所述Qzxxt为决策变量进行优化，获得以所述日调节水电站的洪水流量峰值最小为优化导向的所述日调节水电站的洪水流量峰值及对应的闸门调控总工作量，或获得以所述日调节水电站的闸门调控工作量最小为优化导向的闸门调控总工作量及对应的所述日调节水电站的洪水流量峰值。
进一步地，所述步骤二中，还包括：以所述日调节水电站的洪水流量峰值最小建立所述优化调控模型的目标函数K1，所述目标函数K1满足关系式：
Figure BDA0002279833290000025
其中，K1表示所述日调节水电站的洪水流量峰值最小目标，T表示所述优化调控模型的总时段数，t表示所述优化调控模型的第t时段，t＝1,2,…,T，
Figure BDA0002279833290000031
表示经过所述枢纽水电站调蓄后的所述日调节水电站的洪水流量，所述优化调控模型的第t时段和所述
Figure BDA0002279833290000032
的第t时段相同，所述
Figure BDA0002279833290000033
的单位为m3/s。
进一步地，所述步骤二中，还包括：以所述日调节水电站的闸门调控工作量最小建立所述优化调控模型的目标函数K2，所述目标函数K2满足关系式：
Figure BDA0002279833290000034
其中，K2表示所述日调节水电站的闸门调控工作量最小目标，T表示所述优化调控模型的总时段数，t表示所述优化调控模型的第t时段，t＝1,2,…,T，所述优化调控模型的第t时段和所述
Figure BDA0002279833290000035
的第t时段相同，k表示所述优化调控模型的第k号闸门，ngate表示所述优化调控模型的总闸门数；
Figure BDA0002279833290000036
表示第k号闸门在第t时段初的开度，
Figure BDA0002279833290000037
表示第k号闸门在第t时段末的开度；其中，开度的单位为m。
进一步地，所述步骤二中，还包括：将水量平衡设置为所述优化调控模型的约束条件，所述水量平衡满足关系式：
Figure BDA0002279833290000038
Figure BDA0002279833290000039
Figure BDA00022798332900000310
其中，
Figure BDA00022798332900000311
表示所述枢纽水电站在第t时段末的库容，
Figure BDA00022798332900000312
表示所述枢纽水电站在第t时段初的库容；
Figure BDA00022798332900000313
表示所述枢纽水电站第t时段的上游来流量，其单位为m3/s；
Figure BDA0002279833290000041
表示所述日调节水电站在第t时段末的库容，
Figure BDA0002279833290000042
表示所述日调节水电站在第t时段初的库容，
Figure BDA0002279833290000043
表示所述日调节水电站的发电流量，
Figure BDA0002279833290000044
表示所述日调节水电站的闸门下泄流量，Δt表示时间步长；其中，库容的单位为亿立方米。
进一步地，所述
Figure BDA0002279833290000045
的上限为所述枢纽水电站的发电流量与溢弃流量之和，所述
Figure BDA0002279833290000046
的下限为所述枢纽水电站的发电流量。
进一步地，所述步骤二中，还包括：将所述枢纽水电站的调蓄限度设置为所述优化调控模型的约束条件，所述枢纽水电站的调蓄限度满足关系式：
Figure BDA0002279833290000047
其中，
Figure BDA0002279833290000048
表示所述枢纽水电站对所述
Figure BDA0002279833290000049
削峰所动用的调蓄库容，其单位为亿立方米；
Figure BDA00022798332900000410
分别表示所述枢纽水电站能够动用的调蓄库容的上限和下限。
进一步地，所述步骤二中，还包括：对所述
Figure BDA00022798332900000411
设置N个(N>1)流量等级，将所述日调节水电站的闸门开度限度设置为所述优化调控模型的约束条件，所述日调节水电站的闸门开度限度满足关系式：
Figure BDA00022798332900000412
其中，
Figure BDA00022798332900000413
表示所述
Figure BDA00022798332900000414
属于第n个(n∈N)流量等级时的第k号闸门所对应的开度；
Figure BDA00022798332900000415
表示所述
Figure BDA00022798332900000416
属于第n个(n∈N)流量等级时的下限；
Figure BDA00022798332900000417
表示所述
Figure BDA00022798332900000418
属于第n个(n∈N)流量等级时的上限，下限和上限的单位均为m3/s。
进一步地，所述步骤二中，还包括：将所述枢纽水电站的水位和库容、所述日调节水电站的水位和库容设置为所述优化调控模型的约束条件，该约束条件满足关系式：
Figure BDA0002279833290000051
Figure BDA0002279833290000052
Figure BDA0002279833290000053
Figure BDA0002279833290000054
Figure BDA0002279833290000055
Figure BDA0002279833290000056
其中，
Figure BDA0002279833290000057
分别表示所述枢纽水电站的库容的上限和下限；
Figure BDA0002279833290000058
分别表示所述日调节水电站的库容的上限和下限；
Figure BDA0002279833290000059
表示所述枢纽水电站在第t时段末的水位，
Figure BDA00022798332900000510
表示所述日调节水电站在第t时段末的水位，水位的单位为m；
Figure BDA00022798332900000511
分别表示所述枢纽水电站的水位的上限和下限；
Figure BDA00022798332900000512
分别表示所述日调节水电站的水位的上限和下限；fA v～z( )表示所述枢纽水电站的库容与水位之间所满足的正相关函数；
Figure BDA00022798332900000513
表示所述日调节水电站的库容与水位之间所满足的正相关函数，fv～z为水位库容的关系曲线。
进一步地，所述步骤二中，还包括：将所述枢纽水电站的尾水位和所述日调节水电站的尾水位设置为所述优化调控模型的约束条件，该约束条件满足关系式：
Figure BDA00022798332900000514
Figure BDA00022798332900000515
其中，
Figure BDA00022798332900000516
表示所述枢纽水电站在第t时段的尾水位，
Figure BDA00022798332900000517
表示所述日调节水电站在第t时段的尾水位，尾水位的单位为m；
Figure BDA0002279833290000061
表示所述枢纽水电站的尾水位与总下泄流量之间所满足的正相关函数；
Figure BDA0002279833290000062
表示所述日调节水电站的发电流量及闸门下泄流量与尾水位之间所满足的正相关函数。
进一步地，所述步骤三中，还包括：生成若干组满足所述约束条件的所述Qzxxt，运用遗传算法对若干组所述Qzxxt进行迭代寻优计算。
进一步地，所述遗传算法中的交叉和变异步骤的交叉概率和变异概率分别设置为0.1。
本发明的有益效果在于：本发明分别以所述日调节水电站的洪水流量峰值最小及所述日调节水电站的闸门调控工作量最小为目标函数建立优化调控模型，生成下泄流量序列Qzxxt，以所述Qzxxt为决策变量进行优化，最大化地对洪水削峰，保障了下游水力防冲安全，并显著地降低了闸门调控的总工作量，有效地提高了日调节水电站的防冲安全度和下游易冲刷位置的水力安全度。
附图说明
图1为本发明实施例1的流程示意图。
图2为本发明实施例1日调节水电站上游来流量和上下游两电站区间来流量的示意图。
图3为本发明实施例1闸门开度限度的示意图。
图4为本发明实施例1遗传算法的流程示意图。
图5为本发明实施例2遗传算法的流程示意图。
图6为本发明实施例1和实施例2的解集图。
具体实施方式
如在说明书及权利要求当中使用了某些词汇来指称特定组件，本领域技术人员应可理解，制造商可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名称的差异来作为区分组件的方式，而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”为一开放式用语，故应解释成“包含但不限定于”。“大致”是指在可接受的误差范围内，本领域技术人员能够在一定误差范围内解决技术问题，基本达到技术效果。
在本发明的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、水平”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
在发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
以下结合附图和具体实施例对本发明作进一步详细说明，但不作为对本发明的限定。
实施例1
如图1～2所示，一种考虑下游河道水力安全的水电站泄洪优化调控方法，包括：
步骤一：选取枢纽水电站下游的具有防冲安全问题的日调节水电站作为调控对象，其中，以二滩水电站为枢纽水电站，以桐子林水电站为日调节水电站，桐子林水电站在初始时段的来流量约为11000m3/s，以15天为调控周期，根据水文预报获取枢纽水电站与日调节水电站的区间来流量
Figure BDA0002279833290000081
以及枢纽水电站调蓄前的日调节水电站的上游来流量
Figure BDA0002279833290000082
获取枢纽水电站调蓄前第t时段的下泄流量
Figure BDA0002279833290000083
t的取值范围为1～120，每个时段包含3个小时。
步骤二：分别以日调节水电站的洪水流量峰值最小及以日调节水电站的闸门调控工作量最小为目标函数建立优化调控模型，并设置优化调控模型的约束条件。
以日调节水电站的洪水流量峰值最小建立优化调控模型的目标函数K1，目标函数K1满足关系式：
Figure BDA0002279833290000084
其中，K1表示日调节水电站的洪水流量峰值最小目标，T表示优化调控模型的总时段数，总时段数为120；t表示优化调控模型的第t时段，t＝1,2,…,T，每个时段包含3个小时；
Figure BDA0002279833290000085
表示经过枢纽水电站调蓄后的日调节水电站的洪水流量。
以日调节水电站的闸门调控工作量最小建立优化调控模型的目标函数K2，目标函数K2满足关系式：
Figure BDA0002279833290000091
其中，K2表示日调节水电站的闸门调控工作量最小目标，T表示优化调控模型的总时段数，总时段数为120；t表示优化调控模型的第t时段，t＝1,2,…,120，每个时段包含3个小时；k表示优化调控模型的第k号闸门，ngate表示优化调控模型的总闸门数；
Figure BDA0002279833290000092
表示第k号闸门在第t时段初的开度，
Figure BDA0002279833290000093
表示第k号闸门在第t时段末的开度。
在步骤二中，首先将水量平衡设置为优化调控模型的约束条件，水量平衡满足关系式：
Figure BDA0002279833290000094
Figure BDA0002279833290000095
Figure BDA0002279833290000096
其中，
Figure BDA0002279833290000097
表示枢纽水电站在第t时段末的库容，
Figure BDA0002279833290000098
表示枢纽水电站在第t时段初的库容；
Figure BDA0002279833290000099
表示枢纽水电站第t时段的上游来流量；
Figure BDA00022798332900000910
表示日调节水电站在第t时段末的库容，
Figure BDA00022798332900000911
表示日调节水电站在第t时段初的库容，
Figure BDA00022798332900000912
表示日调节水电站的发电流量，
Figure BDA00022798332900000913
表示日调节水电站的闸门下泄流量，Δt表示时间步长。
然后，将枢纽水电站的调蓄限度设置为优化调控模型的约束条件，枢纽水电站的调蓄限度满足关系式：
Figure BDA00022798332900000914
其中，
Figure BDA00022798332900000915
表示枢纽水电站对
Figure BDA00022798332900000916
削峰所动用的调蓄库容；
Figure BDA0002279833290000101
分别表示枢纽水电站能够动用的调蓄库容的上限和下限。
接着，将枢纽水电站的水位和库容、日调节水电站的水位和库容设置为优化调控模型的约束条件，该约束条件满足关系式：
Figure BDA0002279833290000102
Figure BDA0002279833290000103
Figure BDA0002279833290000104
Figure BDA0002279833290000105
Figure BDA0002279833290000106
Figure BDA0002279833290000107
其中，
Figure BDA0002279833290000108
分别表示枢纽水电站的库容的上限和下限；
Figure BDA0002279833290000109
分别表示日调节水电站的库容的上限和下限；
Figure BDA00022798332900001010
表示枢纽水电站在第t时段末的水位，
Figure BDA00022798332900001011
表示日调节水电站在第t时段末的水位；
Figure BDA00022798332900001012
分别表示枢纽水电站的水位的上限和下限；
Figure BDA00022798332900001013
分别表示日调节水电站的水位的上限和下限；fA v～z( )表示枢纽水电站的库容与水位之间所满足的正相关函数；
Figure BDA00022798332900001014
表示日调节水电站的库容与水位之间所满足的正相关函数。
最后，将枢纽水电站的尾水位和日调节水电站的尾水位设置为优化调控模型的约束条件，该约束条件满足关系式：
Figure BDA00022798332900001015
Figure BDA00022798332900001016
其中，
Figure BDA00022798332900001017
表示枢纽水电站在第t时段的尾水位，
Figure BDA00022798332900001018
表示日调节水电站在第t时段的尾水位；
Figure BDA00022798332900001019
表示枢纽水电站的尾水位与总下泄流量之间所满足的正相关函数；
Figure BDA0002279833290000111
表示日调节水电站的发电流量及闸门下泄流量与尾水位之间所满足的正相关函数。
在步骤二中，还对
Figure BDA0002279833290000112
设置N个(N>1)流量等级，以桐子林水电站的第5号、第6号和第7号闸门为例，记录桐子林水电站的明渠段在各泄洪流量等级下的能够保障下游水力安全闸门开度，桐子林水电站明渠段泄洪闸开度情况如表1所示。
表1桐子林水电站明渠段泄洪闸开度情况
Figure BDA0002279833290000113
将日调节水电站的闸门开度限度设置为优化调控模型的约束条件，日调节水电站的闸门开度限度满足关系式：
Figure BDA0002279833290000114
其中，
Figure BDA0002279833290000115
表示
Figure BDA0002279833290000116
属于第n个(n∈N)流量等级时的第k号闸门所对应的开度；
Figure BDA0002279833290000117
表示
Figure BDA0002279833290000118
属于第n个(n∈N)流量等级时的下限；
Figure BDA0002279833290000121
表示
Figure BDA0002279833290000122
属于第n个(n∈N)流量等级时的上限，因此，桐子林水电站的闸门开度限度设置为如图3所示，当设置好闸门开度限度之后，可以根据目标函数K2计算满足日调节水电站的闸门开度限度约束的闸门调控总工作量。
步骤三：如图4所示，生成100组满足约束条件的下泄流量序列Qzxxt，100组Qzxxt构成初始群体，每组Qzxxt为初始群体的初始个体，每组Qzxxt包含100个
Figure BDA0002279833290000123
运用遗传算法对100组Qzxxt进行迭代寻优计算，并设置100个进化代，每代群体均包含100个个体，其中，
Figure BDA0002279833290000124
的上限为枢纽水电站的发电流量与溢弃流量之和，
Figure BDA0002279833290000125
的下限为枢纽水电站的发电流量，遗传算法中的交叉和变异步骤的交叉概率和变异概率分别设置为0.1；再以Qzxxt为决策变量进行优化，本实施例以目标函数K2为优化导向，获得以日调节水电站的闸门调控工作量最小为优化导向的闸门调控总工作量，及对应的日调节水电站的洪水流量峰值，对应的日调节水电站的洪水流量峰值由目标函数K1求得，优化结果如表2所示。
表2以目标函数K2为优化导向的优化结果
Figure BDA0002279833290000126
从上述步骤可以发现，本实施例通过模型试验手段保证每个泄洪流量等级下的泄洪闸下游形成淹没式水跃，并寻求满足流速流态安全要求的明渠段的最大泄流量，由此确定明渠段在各泄洪流量等级下的闸门开度；在明渠段按照上述要求最大化分担河床段泄洪压力后，再寻求河床段河道导墙内侧临底流速的闸门开度最优组合，从而得到各泄洪流量等级下满足水力安全要求的闸门开度最优组合。
本实施例中，还可以运用优化调控模型获得以日调节水电站的闸门调控工作量最小为优化导向的枢纽水电站的下泄流量序列。
实施例2
如图5所示，本实施例与实施例1不同的是，步骤三中，以目标函数K1为优化导向，获得优化后的枢纽水电站的下泄流量序列，并获得以日调节水电站的洪水流量峰值最小为优化导向的日调节水电站的洪水流量峰值，及对应的闸门调控总工作量，对应的闸门调控总工作量由目标函数K2求得，优化结果如表3所示。
表3以目标函数K1为优化导向的优化结果
Figure BDA0002279833290000131
如图6所示，以15天为调控周期，获得以目标函数K1为优化导向和以目标函数K2为优化导向的解集，其中，以目标函数K1为优化导向的枢纽水电站的下泄流量的波动程度要小于以目标函数K2为优化导向的枢纽水电站的下泄流量。
由上述实施例可知，本发明可以通过目标函数K1为导向或通过目标函数K2为导向获得优化结果，本发明的泄洪优化方法最大程度上对洪水进行了削峰，保障了下游水力防冲安全，并显著地降低了闸门调控的总工作量，提高了日调节水电站的防冲安全性能，提高了下游易冲刷位置的水力安全度。
根据上述说明书的揭示和教导，本发明所属领域的技术人员还能够对上述实施方式进行变更和修改。因此，本发明并不局限于上述的具体实施方式，凡是本领域技术人员在本发明的基础上所作出的任何显而易见的改进、替换或变型均属于本发明的保护范围。此外，尽管本说明书中使用了一些特定的术语，但这些术语只是为了方便说明，并不对本发明构成任何限制。

Claims (10)

1.一种考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，包括：
步骤一：选取枢纽水电站下游的具有防冲安全问题的日调节水电站作为调控对象，获取所述枢纽水电站与所述日调节水电站的区间来流量
Figure FDA0002279833280000011
获取所述枢纽水电站调蓄前的所述日调节水电站的上游来流量
Figure FDA0002279833280000012
获取所述枢纽水电站调蓄前第t时段的下泄流量
Figure FDA0002279833280000013
步骤二：分别以所述日调节水电站的洪水流量峰值最小及以所述日调节水电站的闸门调控工作量最小为目标函数建立优化调控模型，并设置所述优化调控模型的约束条件；
步骤三：生成下泄流量序列Qzxxt,所述Qzxxt由若干个所述
Figure FDA0002279833280000014
组成，以所述Qzxxt为决策变量进行优化，获得以所述日调节水电站的洪水流量峰值最小为优化导向的所述日调节水电站的洪水流量峰值及对应的闸门调控总工作量，或获得以所述日调节水电站的闸门调控工作量最小为优化导向的闸门调控总工作量及对应的所述日调节水电站的洪水流量峰值。
2.如权利要求1所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
以所述日调节水电站的洪水流量峰值最小建立所述优化调控模型的目标函数K1，所述目标函数K1满足关系式：
Figure FDA0002279833280000015
其中，K1表示所述日调节水电站的洪水流量峰值最小目标，T表示所述优化调控模型的总时段数，t表示所述优化调控模型的第t时段，t＝1,2,…,T，
Figure FDA0002279833280000021
表示经过所述枢纽水电站调蓄后的所述日调节水电站的洪水流量。
3.如权利要求1所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
以所述日调节水电站的闸门调控工作量最小建立所述优化调控模型的目标函数K2，所述目标函数K2满足关系式：
Figure FDA0002279833280000022
其中，K2表示所述日调节水电站的闸门调控工作量最小目标，T表示所述优化调控模型的总时段数，t表示所述优化调控模型的第t时段，t＝1,2,…,T，k表示所述优化调控模型的第k号闸门，ngate表示所述优化调控模型的总闸门数；
Figure FDA0002279833280000023
表示第k号闸门在第t时段初的开度，
Figure FDA0002279833280000024
表示第k号闸门在第t时段末的开度。
4.如权利要求2或3所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
将水量平衡设置为所述优化调控模型的约束条件，所述水量平衡满足关系式：
Figure FDA0002279833280000025
Figure FDA0002279833280000026
Figure FDA0002279833280000027
其中，
Figure FDA0002279833280000028
表示所述枢纽水电站在第t时段末的库容，
Figure FDA0002279833280000029
表示所述枢纽水电站在第t时段初的库容；
Figure FDA0002279833280000031
表示所述枢纽水电站第t时段的上游来流量；
Figure FDA0002279833280000032
表示所述日调节水电站在第t时段末的库容，
Figure FDA0002279833280000033
表示所述日调节水电站在第t时段初的库容，
Figure FDA0002279833280000034
表示所述日调节水电站的发电流量，
Figure FDA0002279833280000035
表示所述日调节水电站的闸门下泄流量，Δt表示时间步长。
5.如权利要求4所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
将所述枢纽水电站的调蓄限度设置为所述优化调控模型的约束条件，所述枢纽水电站的调蓄限度满足关系式：
Figure FDA0002279833280000036
其中，
Figure FDA0002279833280000037
表示所述枢纽水电站对所述
Figure FDA0002279833280000038
削峰所动用的调蓄库容；
Figure FDA0002279833280000039
分别表示所述枢纽水电站能够动用的调蓄库容的上限和下限。
6.如权利要求4所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
对所述
Figure FDA00022798332800000310
设置N个(N>1)流量等级，将所述日调节水电站的闸门开度限度设置为所述优化调控模型的约束条件，所述日调节水电站的闸门开度限度满足关系式：
Figure FDA00022798332800000311
其中，
Figure FDA00022798332800000312
表示所述
Figure FDA00022798332800000313
属于第n个(n∈N)流量等级时的第k号闸门所对应的开度；
Figure FDA00022798332800000314
表示所述
Figure FDA00022798332800000315
属于第n个(n∈N)流量等级时的下限；
Figure FDA00022798332800000316
表示所述
Figure FDA00022798332800000317
属于第n个(n∈N)流量等级时的上限。
7.如权利要求4所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
将所述枢纽水电站的水位和库容、所述日调节水电站的水位和库容设置为所述优化调控模型的约束条件，该约束条件满足关系式：
Figure FDA0002279833280000041
Figure FDA0002279833280000042
Figure FDA0002279833280000043
Figure FDA0002279833280000044
Figure FDA0002279833280000045
Figure FDA0002279833280000046
其中，
Figure FDA0002279833280000047
分别表示所述枢纽水电站的库容的上限和下限；
Figure FDA0002279833280000048
分别表示所述日调节水电站的库容的上限和下限；
Figure FDA0002279833280000049
表示所述枢纽水电站在第t时段末的水位，
Figure FDA00022798332800000410
表示所述日调节水电站在第t时段末的水位；
Figure FDA00022798332800000411
分别表示所述枢纽水电站的水位的上限和下限；
Figure FDA00022798332800000412
分别表示所述日调节水电站的水位的上限和下限；fA v～z()表示所述枢纽水电站的库容与水位之间所满足的正相关函数；
Figure FDA00022798332800000413
表示所述日调节水电站的库容与水位之间所满足的正相关函数。
8.如权利要求7所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤二中，还包括：
将所述枢纽水电站的尾水位和所述日调节水电站的尾水位设置为所述优化调控模型的约束条件，该约束条件满足关系式：
Figure FDA00022798332800000414
Figure FDA0002279833280000051
其中，
Figure FDA0002279833280000052
表示所述枢纽水电站在第t时段的尾水位，
Figure FDA0002279833280000053
表示所述日调节水电站在第t时段的尾水位；
Figure FDA0002279833280000054
表示所述枢纽水电站的尾水位与总下泄流量之间所满足的正相关函数；
Figure FDA0002279833280000055
表示所述日调节水电站的发电流量及闸门下泄流量与尾水位之间所满足的正相关函数。
9.如权利要求1所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于，所述步骤三中，还包括：
生成若干组满足所述约束条件的所述Qzxxt，运用遗传算法对若干组所述Qzxxt进行迭代寻优计算。
10.如权利要求9所述的考虑下游河道水力安全的水电站泄洪优化调控方法，其特征在于：所述遗传算法中的交叉和变异步骤的交叉概率和变异概率分别设置为0.1。
CN201911136971.8A 2019-11-19 2019-11-19 一种考虑下游河道水力安全的水电站泄洪优化调控方法 Pending CN110929929A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911136971.8A CN110929929A (zh) 2019-11-19 2019-11-19 一种考虑下游河道水力安全的水电站泄洪优化调控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911136971.8A CN110929929A (zh) 2019-11-19 2019-11-19 一种考虑下游河道水力安全的水电站泄洪优化调控方法

Publications (1)

Publication Number Publication Date
CN110929929A true CN110929929A (zh) 2020-03-27

Family

ID=69851180

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911136971.8A Pending CN110929929A (zh) 2019-11-19 2019-11-19 一种考虑下游河道水力安全的水电站泄洪优化调控方法

Country Status (1)

Country Link
CN (1) CN110929929A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117094144A (zh) * 2023-08-16 2023-11-21 郑州大学 一种考虑泄洪设施闸门控制的水库防洪调度方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646294A (zh) * 2013-11-12 2014-03-19 国网电力科学研究院 一种考虑闸门操作规则的水电站洪水优化调度方法
CN104182634A (zh) * 2014-08-21 2014-12-03 华中科技大学 一种梯级水电站联合运行水位控制断面优选方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646294A (zh) * 2013-11-12 2014-03-19 国网电力科学研究院 一种考虑闸门操作规则的水电站洪水优化调度方法
CN104182634A (zh) * 2014-08-21 2014-12-03 华中科技大学 一种梯级水电站联合运行水位控制断面优选方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘广宇等: "基于优化组合方案的梯级水库泄洪闸门数字化及其应用" *
王方方等: "考虑水电调蓄的西江水库群应急防洪调度研究" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117094144A (zh) * 2023-08-16 2023-11-21 郑州大学 一种考虑泄洪设施闸门控制的水库防洪调度方法
CN117094144B (zh) * 2023-08-16 2024-04-26 郑州大学 一种考虑泄洪设施闸门控制的水库防洪调度方法

Similar Documents

Publication Publication Date Title
WO2021073192A1 (zh) 一种考虑预报误差降低水库洪水起调水位的预报调度方法
CN107992989B (zh) 一种水库群洪水资源化预蓄预泄风险决策方法
Harby et al. Pumped storage hydropower
CN110210109B (zh) 一种河网中堰闸工程反向水流的数值模拟方法及系统
CN111445080A (zh) 一种基于随机森林的西北内陆河生态基流预测方法
CN101714193B (zh) 面向河流生态系统保护的发电型水库调度函数优化方法
CN110929929A (zh) 一种考虑下游河道水力安全的水电站泄洪优化调控方法
CN115271304A (zh) 基于气候变化的跨流域调水工程优化调度方法
CN109919397A (zh) 平原河网地区排涝容量校验优化方法
CN114239992A (zh) 一种水库动态汛限水位制定方法
CN108133331B (zh) 一种黄河下游宽滩区滞洪沉沙功能和减灾效应的评价方法
Kasperek et al. Hydropower generation on the Nysa Klodzka river
CN109933892B (zh) 一种水温模型电站下泄流量边界条件的修正方法
CN102776872B (zh) 水电站高效防洪发电的优化方法
CN112766565B (zh) 梯级水库群防洪优化调度分数阶逐步优化方法
CN112580230B (zh) 一种水库水温结构分析方法
CN111754362A (zh) 一种多级放空坝放空能力评估方法
CN111126693B (zh) 基于上游水库运行对下游电站发电能力影响的调度方法
CN111144621B (zh) 一种基于趋势判别的水库群防洪发电竞争关系判别方法
Nguyen et al. Short-term reservoir system operation for flood mitigation with 1D hydraulic model
CN114186877A (zh) 一种考虑水库调节库容的日内水光互补计算方法
US20200048124A1 (en) System and method for closed-loop dissolved oxygen monitoring and control
Baczyński et al. Possibility of power generation control in small hydroinstallations for a period of several days
CN107844871A (zh) 一种电站发电效益优化方法
CN105256774B (zh) 一种利用下游水库回水减小上游水电站截流难度的方法

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