随着堵河水库群完建并进入正常运行期运用，堵河流域形成了以潘口水库为核心的控制性水库群。为减少堵河流域水电站在汛末蓄水期无效弃水，提高水资源综合利用率和兴利效益，本文开展了水库汛末蓄水调度方式优化研究。首先采用K判别系数法，分析了流域梯级水库蓄水次序；进而针对来水较枯、潘口水库汛后无法蓄满的年份，提出了基于改进K判别系数法的梯级水库蓄水次序优化模型；并以蓄满率和发电指标对蓄水方案进行了评价。研究表明，改进后的梯级水库联合蓄水方案，可有效增加总发电效益约0.266亿kW·h。 With the completion and being in practice of Duhe reservoir group, the operation of the key reservoirs represented by Pankou has begun to take shape. In order to reduce the ineffective water abandonment of cascade hydropower stations at the end of flood season and improve the comprehensive utilization rate of water resources, the optimal operation model of cascade reservoirs at the end of flood season is established. Firstly, the storage inlet order of cascade reservoirs is analyzed by K discrimination coefficient method. Then, for years when the incoming water flow is relatively less and Pankou Reservoir cannot be fully filled after flood season, the water storage sequence of cascade reservoirs is improved based on the K discrimination coefficient method, and the impound scheme is evaluated according to the full storage rate and power generation index. The joint operation of reservoir group can increase 26.6 million kW·h power generation annually.
随着堵河水库群完建并进入正常运行期运用，堵河流域形成了以潘口水库为核心的控制性水库群。为减少堵河流域水电站在汛末蓄水期无效弃水，提高水资源综合利用率和兴利效益，本文开展了水库汛末蓄水调度方式优化研究。首先采用K判别系数法，分析了流域梯级水库蓄水次序；进而针对来水较枯、潘口水库汛后无法蓄满的年份，提出了基于改进K判别系数法的梯级水库蓄水次序优化模型；并以蓄满率和发电指标对蓄水方案进行了评价。研究表明，改进后的梯级水库联合蓄水方案，可有效增加总发电效益约0.266亿kW·h。
梯级水库，汛末蓄水，K判别系数，堵河流域
Wei Wang1, Xin Zhang2, Weiqing Yu3, Shaofeng Gu4
1Hanjiang Water Resources & Hydropower Group Co., Ltd., Wuhan Hubei
2National Energy Changyuan Shiyan Hydropower Development Co., Ltd., Shiyan Hubei
3Changjiang Institute of Survey, Planning, Design and Research, Wuhan Hubei
4Hubei Electric Power Co., Ltd. Huanglongtan Hydropower Plant, Shiyan Hubei
Received: Nov. 7th, 2022; accepted: Nov. 20th, 2022; published: Dec. 31st, 2022
With the completion and being in practice of Duhe reservoir group, the operation of the key reservoirs represented by Pankou has begun to take shape. In order to reduce the ineffective water abandonment of cascade hydropower stations at the end of flood season and improve the comprehensive utilization rate of water resources, the optimal operation model of cascade reservoirs at the end of flood season is established. Firstly, the storage inlet order of cascade reservoirs is analyzed by K discrimination coefficient method. Then, for years when the incoming water flow is relatively less and Pankou Reservoir cannot be fully filled after flood season, the water storage sequence of cascade reservoirs is improved based on the K discrimination coefficient method, and the impound scheme is evaluated according to the full storage rate and power generation index. The joint operation of reservoir group can increase 26.6 million kW·h power generation annually.
Keywords:Cascade Reservoirs, Impounding at the End of Flooding Season, K Discrimination Coefficient, Duhe River Basin
Copyright © 2022 by author(s) and Wuhan University.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
堵河流域位于鄂西北汉江右岸，为汉江一级支流，干流全长约354 km，流域面积约12,430 km2。随着堵河梯级水库完建并进入正常运行期，水库上下游之间、干支流之间径流特性处于动态变化；由于潘口水库汛后需蓄水库容较大，汛末蓄水压力大，特别是遇流域枯水年份，需采用梯级联合优化调度方式提高梯级水库汛末有效蓄水。
许多学者对汛末梯级水库蓄水方式开展了大量研究。易灵等 [
K判别系数法 [
堵河干支流已建成主要水库16座，总库容约50亿m3、总调节库容约23亿m3；包括干流泗河河段的鄂坪、汇湾、白果坪，干流两河口以下的潘口、小漩、黄龙滩；支流官渡河龙背湾、松树岭，支流杨家河石庙子水库，支流潭口河的冯家湾水库，支流竹溪河水库，支流谭家河水库，支流泉河红岩一级、大峡、白沙河，支流霍河上的霍河水库。综合考虑梯级水库防洪任务及现有水文资料情况，本文重点讨论汛期预留有防洪库容或库容较大的鄂坪、龙背湾、潘口、霍河、黄龙滩水库的汛后蓄水方式。堵河流域水库群示意图见图1。
图1. 堵河流域水库群示意图
本次研究采用1956~2018年长系列旬径流资料，梯级水库主要调度参数及调度图采用设计阶段成果。堵河流域梯级水库防洪库容释放时间为：潘口防洪库容释放时间为8月21日。黄龙滩虽然没有预留防洪库容，但是在实际运行中，汛期结合预报预泄提前将水位预泄至243 m迎洪，结合上游潘口防洪库容释放时机，本文假设黄龙滩8月21日可以开始蓄水。
1) 目标函数
在满足综合利用要求的基础上，以梯级电站发电效益最大为目标。
2) 约束条件
以水库坝前水位、入库流量等为输入条件，以最下一级水库的下泄流量或者蓄水量作为约束条件。
上下游水库之间水量平衡计算公式：
Q i n , t i + 1 = Q o u t , t i + Q t i , i + 1 (1)
Q o u t , t i = f ( h t i , Q i n , t i ) (2)
W t i + 1 = f 2 i + 1 ( h t i + 1 ) (3)
W t + 1 i + 1 = W t i + 1 + ( Q i n , t i + 1 − Q s , t i + 1 ) ⋅ Δ t (4)
h t + 1 i + 1 = f 3 i + 1 ( W t + 1 i + 1 ) (5)
式中： Q i n , t i + 1 为t时段下游水库的平均入库流量； Q o u t , t i 为t时段上游水库的平均出库流量； Q t i , i + 1 为t时段上游水库和下游水库之间的平均区间流量； f ( h t i , Q i n , t i ) 为水库i应用调度图计算出库流量函数； W t i + 1 为下游水库时段初的库容； f 2 i + 1 为下游水库的水位库容计算函数，由水位计算库容； h t i + 1 为下游水库t时段初的坝前水位； W t + 1 i + 1 为下游水库t + 1时段初(即t时段末)的库容； Δ t 为计算时段； Q s , t i + 1 为下游水库t时段的下泄流量，根据蓄水调度需求计算下泄流量； f 3 i + 1 为下游水库的水位库容计算函数，由库容计算水位； h t + 1 i + 1 为下游水库t + 1时段初的坝前水位。
3) 决策变量
堵河梯级水库中，鄂坪位于最上一级，黄龙滩位于干流最后一级，龙背湾、霍河等水库均位于支流。根据梯级水库位置关系，比较有相互水力联系水库的K值，选择K值大的水库 i max 为当前计算时段的蓄水水库。
K值计算公式：
K i = W i + ∑ V i F i ∑ H i (6)
式中， W i 表示梯级水库群中第i个水库的入库总水量； ∑ V i 为第i个水库及其上游梯级各水库可供发电的总蓄水量； F i 表示第i个水库的水面面积； ∑ H i 表示第i个水库及其所有下游水库的总水头。
堵河梯级水库是典型的山区型水库，发电水头是影响发电效益敏感性因素，因此尽量提高水库蓄满率，是提高梯级发电效益的手段之一。按设计方案调度梯级水库长系列蓄水过程显示，一般情况下，黄龙滩水库汛末基本上都可以蓄满。对于枯水年份，由于蓄水期9~11月天然来水偏少，导致在这些年份潘口、鄂坪、龙背湾三个水库均有汛后无法蓄满的情况发生。
为此，针对来水较枯，梯级水库汛后无法蓄满的年份，开展提高梯级水库汛后蓄满率的调度方式研究。采用K判别系数法判别堵河梯级汛末蓄水方式，并结合对潘口、黄龙滩水库蓄水条件加以约束的尝试，寻求更有利的汛末梯级水库蓄水方式。不同蓄水条件约束方案见表1。
| 黄龙滩 | 潘口 | 鄂坪、龙背湾 | 霍河 | |
|---|---|---|---|---|
| 方案1 | K值最大 | K值最大，或H黄≥ 246.0 m | K值最大 | K值最大 |
| 方案2 | K值最大 | K值最大，或H黄≥ 246.5 m | K值最大，或H潘≥ 354.5 m | K值最大 |
| 方案3 | K值最大 | K值最大，或H黄≥ 246.0 m | K值最大，或H潘≥ 354.0 m | K值最大 |
| 方案4 | K值最大 | K值最大，或H黄≥ 245.0 m | K值最大，或H潘≥ 354.5 m | K值最大 |
表1. 堵河梯级水库汛末蓄水优化方案
注：表中H黄、H潘分别为当前时段黄龙滩、潘口水库坝前水位。
研究结果显示，K判别系数法判别堵河梯级汛末蓄水方式，基本上没有发生由于K值较大的水库因过多承担蓄水储能任务，而导致水库可能提前达到正常蓄水位，产生弃水的情况，提高了梯级水库汛末蓄满率。尤其是对于来水较枯、潘口水库汛后无法蓄满的年份，潘口库水位抬高的较为明显。这主要是因为梯级水库单独调度时，仅从当前时段自身水库来水和库水位来判断电站运行方式；K判别系数法则结合梯级水库蓄水情势，当
图2. 潘口水库运行过程示意图
下游水库蓄至较高水位时，上游水库可以提前减小出力蓄水，这样既可以提高上游水库运行水位，增加发电水头；也可以尽量避免下游水库提前达到正常蓄水位而产生无效弃水的情况。
以枯水年份较多的1990~2018年为例，绘出潘口水库运行过程示意图见图2。
K判别系数法汛后蓄水较快，对汛后发电略有影响，但水库一旦蓄满，后期运行水位抬高，发电效益相应抬高；总体来看，汛末按K判别式法蓄水增加了梯级电站的发电效益，并且主要增加的是枯水期发电效益，对于来水偏枯、年发电量小于平均值的年份，K判别系数法对发电量的增加尤为明显。梯级水库汛末蓄水情况统计见表2，梯级水库旬发电效益对比示意图见图3，潘口水库年发电效益对比示意图见图4。
与设计阶段拟定的各电站调度方式发电效益相比，梯级水库蓄满率增加1.6%~36.5%，尤其是潘口和霍河水库蓄满率增加明显；梯级总发电效益增加0.240~0.266亿kW·h。随着潘口水库汛末蓄水条件放宽，潘口水库的蓄满率和发电效益随之提高；虽然方案3潘口水库蓄满率和发电效益较方案1、方案4略有下降，从梯级水库联合运行总发电效益来看，方案3优于其他几个方案。因此，本次研究推荐方案3为堵河梯级水库汛末蓄水方式。
| 电站 | 鄂坪 | 龙背湾 | 潘口 | 霍河 | 黄龙滩 | 合计 | |
|---|---|---|---|---|---|---|---|
| 设计阶段方案 | 蓄满率 | 69.80% | 60.30% | 63.50% | 47.60% | 96.80% | |
| 年均发电量(亿kW·h) | 2.01 | 3.39 | 9.65 | 0.46 | 9.74 | 25.248 | |
| K判别系数法 | 蓄满率 | 73% | 66.70% | 88.90% | 84.10% | 98.40% | |
| 年均发电量(亿kW·h) | 2.07 | 3.4 | 9.69 | 0.46 | 9.86 | 25.488 | |
| 优化方案1 | 蓄满率 | 74.60% | 66.70% | 93.70% | 84.10% | 98.40% | |
| 年均发电量(亿kW·h) | 2.08 | 3.4 | 9.7 | 0.46 | 9.87 | 25.511 | |
| 优化方案2 | 蓄满率 | 74.60% | 66.70% | 85.70% | 84.10% | 98.40% | |
| 年均发电量(亿kW·h) | 2.08 | 3.4 | 9.69 | 0.46 | 9.86 | 25.497 | |
| 优化方案3 | 蓄满率 | 74.60% | 68.30% | 92.10% | 84.10% | 98.40% | |
| 年均发电量(亿kW·h) | 2.08 | 3.4 | 9.7 | 0.46 | 9.87 | 25.514 | |
| 优化方案4 | 蓄满率 | 74.60% | 68.30% | 93.70% | 84.10% | 98.40% | |
| 年均发电量(亿kW·h) | 2.08 | 3.4 | 9.7 | 0.46 | 9.87 | 25.512 | |
表2. 堵河梯级水库汛末蓄满率和年均发电量统计
图3. 梯级电站旬发电效益对比示意图
图4. 潘口水库年发电效益对比示意图
本文以梯级水库发电量最大化为目标，对堵河流域联合调度汛末蓄水方式进行初步探索，提出了基于K判别系数法结合潘口、黄龙滩水库蓄水位的堵河梯级汛末蓄水次序判别方式。主要结论如下：
1) 基本上没有发生由于K值较大的水库因过多承担蓄水储能任务，而导致水库可能提前达到正常蓄水位，产生弃水的情况，可减少无效弃水；
2) 提高了梯级水库汛末蓄满率，尤其是对于来水较枯、潘口水库汛后无法蓄满的年份，潘口水库汛末水位抬高较为明显；
3) 增加梯级电站枯水期发电量，年均增加梯级电站发电量约0.266亿kW·h；
4) 判别条件简单，适用于堵河梯级水库汛后蓄水次序判别。
王 伟,张 欣,余蔚卿,顾绍峰. 堵河流域水库群汛末蓄水方式研究Study on Impounding Method of Reservoir Group in Duhe River Basin at the End of Flooding Season[J]. 水资源研究, 2022, 11(06): 657-662. https://doi.org/10.12677/JWRR.2022.116071