基于聚合–分解框架,构建了并联水库群防洪调度函数:在聚合步骤中,将总出库流量表示为聚合水库库容的分段线性函数;在分解步骤中,按各水库来水量比例分解总出库流量,得到各水库出库流量。从而构建了并联水库群防洪优化调度模拟模型,采用遗传算法进行模拟优化,识别分段线性函数的关键参数。以广西西江流域百色、龙滩、青狮潭三个并联水库为实例开展研究,计算结果表明,较常规调度,优化调度在率定期(共10场次洪水)和检验期(共5场次洪水)分别能使下游防洪控制点梧州站平均多削峰970 m3/s和1320 m3/s。因此,提出的方法对指导水库群防洪优化调度有一定价值。<br/>The optimal operating rules for parallel reservoirs were built based on the aggregation-decom- position method. The relationship between the total release of the reservoirs and the aggregation reservoir storage could be expressed by piecewise-linear operating rules in the aggregation step, and the release was decomposed to each parallel reservoir according to its ratio of inflow in the decomposing procedure. Then the optimization model of parallel reservoirs was built, and its key parameters were calibrated by the Genetic Algorithm. Three parallel reservoirs (the Baise, Longtan and Qingshitan reservoirs) of the Xijiang River in China were selected for a case study. The optimal flood operating rules were compared with the conventional operating rules, and the flood peak of reference station for flood control could be decreased by 970 m3/s and 1320 m3/s in calibrated and validation periods, respectively. Thus, this method is beneficial for the reservoir’s real- time operation to decrease the flood disasters.
陈西臻1,2,刘攀1,2,何素明3,张旺1,2,李泽君1,2,张靖文1,2
1武汉大学水资源与水电工程科学国家重点实验室,湖北 武汉
2水资源安全保障湖北省协同创新中心,湖北 武汉
3广西壮族自治区水利电力勘测设计研究院,广西 南宁
Email: chenxizhen15@163.com
收稿日期:2015年2月2日;录用日期:2015年2月9日;发布日期:2015年2月16日
基于聚合–分解框架,构建了并联水库群防洪调度函数:在聚合步骤中,将总出库流量表示为聚合水库库容的分段线性函数;在分解步骤中,按各水库来水量比例分解总出库流量,得到各水库出库流量。从而构建了并联水库群防洪优化调度模拟模型,采用遗传算法进行模拟优化,识别分段线性函数的关键参数。以广西西江流域百色、龙滩、青狮潭三个并联水库为实例开展研究,计算结果表明,较常规调度,优化调度在率定期(共10场次洪水)和检验期(共5场次洪水)分别能使下游防洪控制点梧州站平均多削峰970 m3/s和1320 m3/s。因此,提出的方法对指导水库群防洪优化调度有一定价值。
关键词 :分段线性,并联水库,聚合–分解,遗传算法,调度函数
洪水是我国最严重的自然灾害之一,通过科学的水库防洪优化调度可以有效地减轻灾害损失。对于大流域并联水库群的防洪优化调度问题,由于洪水的时滞性增加了时间维数,水库个数增加了空间维数,优化调度的求解存在维数灾难题。梅亚东[
在水库群隐随机调度中,可采用“优化–拟合–验证”框架编制水库群联合调度函数,即建立水库群防洪优化调度模型,通过确定性优化调度得到最优样本,再采用多元线性回归或其他方法分析最优样本得到调度函数,并基于模拟调度验证调度函数的合理性[
聚合–分解模型的示意图如图1所示,用分段线性表示聚合水库的调度函数,并因此得到总出库流量,再通过分解模块得到各水库出库流量。
基于大系统聚合分解思想,建立水库群联合优化调度模型[
式中:
由于线性调度函数[
图1. 水库群聚合–分解示意图
图2. 分段线性调度函数
式中:
如图2所示,当调度函数的四个参数A、B、C、D横纵坐标确定时,调度函数随之确定。水库的预报库容
出库
式中:
在分解阶段,可根据各水库来水量比例分解出库流量。即
式中:
聚合水库
采用最大削峰准则[
式中:
1) 水量平衡约束:
式中:
2) 库容约束:
式中:
3) 库容曲线约束:
式中:
4) 水库出库流量约束:
式中:
5) 河道汇流约束:
式中:
遗传算法(Genetic Algorithm)是一种通过模拟自然进化过程搜索最优解的方法。该方法按照择优汰劣的法则,将适者生存与自然界基因变异、繁衍等规律相结合,采用随机搜索,以种群为单位,根据个体的适应度进行选择、交叉及变异等操作,最终可收敛于全局最优解[
1) 确定模型的目标函数、个体的适应值、适应度表达方式;
2) 随机产生初始群体,每个个体包含点B、C、D的横纵坐标,共6个优化变量;
3) 计算每个个体的适应值、适应度,记录最佳个体,按精英保留策略复制个体基因;
4) 产生随机数,当达到交叉、变异概率后个体之间实行交叉或变异,从而产生新一代种群;
5) 计算新一代种群个体的适应值、适应度;
6) 重复步骤3)、4)、5)直到满足终止条件为止。
综上所述,可以得出模型[
图3. 模型优化求解流程图
西江属珠江流域干流,跨越滇、黔、桂、粤四省及越南,广西属珠江流域西江水系的面积为202,082 km2,占全区总面积的85.39%。目前广西境内的西江流域面临着比较严峻的防洪形势。西江流域属于亚热带华南季风气候区,受季风气候影响,降雨时空分布不均匀,降雨量年际变化大,汛期(4~9月)降雨量占全年总降雨量的70%~85%,个别年份最多达到90%。修建水库及提防是有效的防洪措施。如图4所示,西江流域目前投入运行的防洪水库主要有郁江上游的百色水库、红水河上游的龙滩水库和柳江上游的青狮潭水库。三个水库的主要参数如表1所示。
百色水库防洪库容为16.4亿m3,按现有调度规则运行的防洪任务为:使下游南宁市的防洪能力从20年一遇提高至50年一遇,使郁江中下游区防洪能力从10年一遇提高至20年一遇,同时减轻右江沿岸洪水灾害。龙滩水库防洪库容为50亿m3,现有防洪调度规则为:1) 在梧州涨水期,龙滩水库下泄流量不大于6000 m3/s,当其入库流量超过25,000 m3/s时,龙滩泄量不超过4000 m3/s;2) 在梧州退水期,当其流量在42,000 m3/s以上时,龙滩水库仍按不大于4000 m3/s下泄;如梧州站流量小于42,000 m3/s,则龙滩按入库流量下泄;3) 当龙滩水库蓄满时,龙滩水库按入库流量下泄。青狮潭水库防洪库容0.541亿m3/s,青狮潭水库的防洪任务为:将桂林市30年一遇洪水降至20年一遇(防洪堤安全泄量4870 m3/s)以下。
对于1962、1966、1968、1970、1974、1976、1978、1983、1988、1994、1996、1998、2001、2002、2005共15场经放大后得到的百年一遇洪水,经百色、龙滩、青狮潭三个水库常规调度后梧州站调洪结果如表2和表3所示,梧州站百年一遇洪水洪峰流量平均削减753 m3/s,其中最高削峰流量3800 m3/s,能将梧州站百年一遇洪水降低至44年一遇,削峰效果最差的年份是1996年,梧州站洪峰流量比设计值增大600 m3/s。
输入1970、1974、1976、1978、1983、1988、1994、1996、1998、2001共10场放大得到的1%设计洪水作为率定期识别调度函数参数,选取1962、1966、1968、2002、2005共5场放大得到的1%设计洪
图4. 西江流域节点概化图
水库 | 汛限水位(m) | 防洪高水位(m) | 防洪库容(亿m3) | 洪水传播至梧州时间(h) |
---|---|---|---|---|
百色 | 214 | 228 | 16.4 | 132 |
龙滩 | 359.3 | 376 | 50 | 120 |
青狮潭 | 224.2 | 226 | 0.514 | 55 |
表1. 三库主要参数
年份 | 百色流量(m3/s) | 龙滩流量(m3/s) | 青狮潭流量(m3/s) | 梧州流量(m3/s) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
设计 | 常规 | 优化 | 设计 | 常规 | 优化 | 设计 | 常规 | 优化 | 设计 | 常规 | 优化 | |
1970 | 4310 | 3000 | 3630 | 18,800 | 13,300 | 12,900 | 1070 | 1090 | 1040 | 52,700 | 52,000 | 50,900 |
1974 | 4560 | 3000 | 2410 | 17,700 | 12,700 | 13,000 | 1980 | 1860 | 1650 | 52,700 | 50,200 | 52,000 |
1976 | 4010 | 2830 | 1060 | 14,900 | 12,000 | 13,000 | 1990 | 1990 | 1650 | 52,700 | 52,600 | 53,500 |
1978 | 6310 | 3000 | 4190 | 12,600 | 10,600 | 11,400 | 1940 | 1940 | 1650 | 52,700 | 52,700 | 52,500 |
1983 | 2980 | 2660 | 1010 | 19,100 | 11,500 | 18,800 | 1400 | 1400 | 1410 | 52,700 | 52,700 | 52,300 |
1988 | 1960 | 1960 | 830 | 20,900 | 12,000 | 11,700 | 1020 | 936 | 836 | 52,700 | 51,200 | 48,800 |
1994 | 2670 | 2610 | 1120 | 12,300 | 8300 | 4500 | 1660 | 1660 | 1550 | 52,700 | 52,700 | 49,800 |
1996 | 2930 | 2690 | 1440 | 15,600 | 11,600 | 9900 | 1440 | 1440 | 1350 | 52,700 | 53,300 | 52,800 |
1998 | 1140 | 1140 | 544 | 8400 | 5200 | 4230 | 2010 | 2010 | 1650 | 52,700 | 52,700 | 50,300 |
2001 | 11100 | 3000 | 4140 | 20,800 | 15,500 | 20,800 | 349 | 349 | 222 | 52,700 | 48,900 | 46,400 |
平均 | 4200 | 2590 | 2040 | 16,100 | 11,300 | 12,000 | 1490 | 1470 | 1300 | 52,700 | 51,900 | 50,900 |
表2. 率定期百年一遇洪水三库及梧州站调度结果
年份 | 百色流量(m3/s) | 龙滩流量(m3/s) | 青狮潭流量(m3/s) | 梧州流量(m3/s) | ||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
设计 | 常规 | 优化 | 设计 | 常规 | 优化 | 设计 | 常规 | 优化 | 设计 | 常规 | 优化 | |
1962 | 3630 | 2880 | 1310 | 11,500 | 6000 | 8310 | 963 | 1060 | 792 | 52,700 | 51,300 | 49,700 |
1966 | 6310 | 3000 | 2400 | 18,200 | 18,200 | 18,700 | 727 | 727 | 636 | 52,700 | 52,600 | 51,800 |
1968 | 7030 | 3000 | 3370 | 21,400 | 21,400 | 20,200 | 1500 | 1100 | 1130 | 52,700 | 51,600 | 48,800 |
2002 | 4370 | 3130 | 3450 | 20,500 | 12,700 | 12,900 | 1470 | 1220 | 1310 | 52,700 | 52,700 | 52,600 |
2005 | 1830 | 1830 | 1360 | 7860 | 5500 | 3810 | 1027 | 935 | 694 | 52,700 | 52,000 | 50,700 |
平均 | 4630 | 2770 | 2380 | 15,900 | 12,800 | 12,800 | 1140 | 1010 | 912 | 52,700 | 52,000 | 50,700 |
表3. 检验期百年一遇洪水三库及梧州站调度结果
水作为检验期检验调度函数的合理性。用遗传算法独立运算10次,每次计算100代,取交叉概率0.6,变异概率0.08,得到最优解,即四点坐标为:A (151.342, 0);B (187.887, 4591);C (188.237, 4620);D (189.178, 4625)。率定期和检验期梧州站百年一遇设计洪水调度结果如表2和3所示。
如表2和表3所示,经三库优化调度后:
1) 在率定期和检验期共十五场洪水中,梧州站洪峰削峰流量最高可达6300 m3/s,能将百年一遇洪水降到24年一遇左右。
2) 在率定期平均削减梧州站洪峰流量1770 m3/s,较常规调度削减970 m3/s;能将百年一遇洪水平均降到63年一遇左右;检验期平均削减洪峰流量1980 m3/s,较常规调度削减1320 m3/s,使洪水平均降到60年一遇左右。
青狮潭入库流量和防洪库容较小,对梧州站洪峰流量影响不大,龙滩入库流量较大,对梧州站影响较大,2005年百年一遇洪水龙滩和梧州调度结果如图5和图6所示。由于龙滩水库入流洪峰与梧州洪峰出现时间不同步,为了起到削峰作用,在入库流量较小时,龙滩水库出库流量比入库流量小。为了找出使梧州削峰效果最明显的调度函数参数,聚合水库的出流量越小越好,这很容易导致出现水库水位上升,防洪库容提前用完,洪峰过后水库水位没有回落的情况。可在洪水结束之后腾空水库。
本文以广西西江流域百色、龙滩、青狮潭三个水库为例,建立了基于聚合-分解的并联水库群防洪优化调度模型,采用分段线性调度函数模拟聚合水库调度函数型式,并按各水库来水量比例分解聚合水库出库流量。输入10场百年一遇设计洪水并用遗传算法识别函数参数,用5场百年一遇洪水验证调度函数的合理性。结果表明,该方法较能有效削减下游防洪控制点的洪峰流量。本次研究采用的是3段线性调度函数模拟聚合调度函数型式,下一步研究将考虑其他形式的调度函数型式[
图5. 龙滩水库1%设计洪水调度过程(2005年)
图6. 梧州站1%设计洪水调度结果(2005年)
优秀青年科学基金(51422907);教育部新世纪优秀人才支持计划(NCET-11-0401);水利部公益性行业科研专项经费项目(201201051);中央分成水资源费项目(项目编号:1261430210028)。
陈西臻,刘 攀,何素明,张 旺,李泽君,张靖文, (2015) 基于聚合–分解的并联水库群防洪优化调度研究Optimization of Parallel Reservoirs for Flood Control Based on the Aggregation-Decomposition Method. 水资源研究,01,23-31. doi: 10.12677/JWRR.2015.41003