1. 引言
近年来,随着经济的快速发展和人们生活水平的不断提高,水资源需求量日益增加 [1],跨地区河流上下游、左右岸用水矛盾也日渐突出,研究流域水库群优化调度是十分必要的。
目前关于流域水库群优化调度的研究较多。陈立华等 [2] 根据梯级水电站优化调度特点,建立并行遗传算法求解多阶段最优化问题数学模型。马立亚等 [3] 以水库运行模拟模型为基础,建立了基于DPSA的梯级水库群中长期优化调度模型。冯仲恺等 [4] 将蛛群优化方法(SSO)这一新型群体智能寻优方法引入梯级水库群优化调度领域,并提出精英集聚蛛群优化方法(ESSO)。张培等 [5] 从风险管理的角度出发,构建了考虑径流预报误差,滞时的不确定性两类风险因子耦合作用的水库群短期优化调度风险分析模型。陈悦云等 [6] 建立面向发电、供水、生态要求的赣江流域水库群优化调度模型。万东辉等 [7] 提出了一种考虑水库调度影响的梯级水库群汛期分期方法。
维持河流健康,实现人水和谐是新时期的治水理念,不仅需要考虑人类社会发展对水资源的需求,同时也要考虑维持河流健康所需要的生态水量。本文以泗河流域作为算例,建立水库闸坝群交互式优化调度模型,以人水和谐理念为指导思想,以实现全流域水资源的统一管理和统一调度为基本准则,以全流域水库群闸坝群的科学合理调度为基本手段,以实现泗河流域水资源系统在灌溉、工业供水、生态补水等方面的综合利用效益最佳为优化目标,探索了实现人与河流和谐发展的有效途径。
2. 研究区概况
泗河发源于泰安新泰市太平顶,流经新泰、泗水、曲阜、兖州、邹城、微山、任城,于任城区辛闸村入南阳湖,是南四湖东区最大的入湖河流。干流全长159 km,流域面积2357 km2,其中济宁市内2030 km2,泰安市327 km2。泗河流域属鲁南泰沂低丘陵与山前冲洪积平原交接地带,东部山峦绵亘,丘陵起伏;西部为泰沂山前冲洪积平原,地形平坦开阔,山区921 km2,丘陵区566 km2,平原区850 km2。泗河流域地理位置及水系如图1所示。
![](//html.hanspub.org/file/10-2411066x10_hanspub.png?20211118082424387)
Figure 1. Geographical location and drainage map of the Sihe River basin
图1. 泗河流域地理位置及水系图
![](//html.hanspub.org/file/10-2411066x11_hanspub.png?20211118082424387)
Figure 2. The present situation of Sihe River basin and the layout of planned water conservancy projects
图2. 泗河流域现状及规划水利工程布局示意图
泗河流域属于暖温带半湿润气候区,四季分明,气候变化显著,夏季炎热,冬季寒冷。流域内年平均气温为13.7℃,多年平均降水量为715 mm。降水分布的年际变化和季节变化都很大,最大降水量为1192 mm (1964年),年最小降水量为375 mm (2002年);年内雨季多集中在7、8月份,占年平均降水量的51.7%。根据实测资料分析计算,泗河流域多年平均径流量44,352万m3,其中济宁泗水、曲阜、邹城、兖州、济宁市区及以下分别为21,935、8764、5736、125、107万m3,泰安7684万m3。
3. 水库闸坝群交互式优化调度模型构建
针对泗河流域上中下游生产、生态用水矛盾尖锐的特点,以全流域河道外生产供水量最大、河道内水面面积最大为优化目标,以县际控制断面下泄水量为约束条件,构建泗河流域水库闸坝群交互式优化调度模型;采用将遗传算法嵌套于长系列变动时历法中的算法进行联合兴利调算;建立评价模型论证方案的合理性。
3.1. 泗河流域水库闸坝群概述
泗河流域现状共建有大型水库1座,中型水库4座,拦河闸坝11座。规划拟在泗河干流新建橡胶坝14座。另外,在泗河下游左岸新建煤矿塌陷区补源水库1座。流域内将形成31座闸坝水库群。泗河流域现状及规划水利工程布局如图2所示。
3.2. 水库闸坝群交互式优化调度模型构建
1) 优化目标
目标函数1:全流域河道外生产供水量最大:
(1)
目标函数2:河道内水面面积最大:
(2)
式中:
为第i个水库(闸坝)第t时段的生产供水量;
为第i个水库(闸坝)第t时段的水面面积;n为水库群及闸坝群数量;T为联合调算中总时段数。
2) 约束条件
本模型的约束条件主要有县际控制断面下泄水量约束、水量平衡约束、水量连续约束等。
3) 决策变量
将各闸坝的灌溉面积作为决策变量。其中,决策变量最小值为各闸坝近5年来平均实灌面积,决策变量最大值为有效灌溉面积。
4) 求解算法
首先利用线性加权和方法将基于流域协同用水的水库闸坝群交互式优化调度模型转化为单目标优化问题,如下式(3)所示:
(3)
再采用将遗传算法嵌套于长系列变动时历法中的算法进行联合供水调算。优化调度模型结构以及泗河流域31座水库群闸坝群水量联合优化调度计算流程,如图3和图4所示。
5) 合理性评价
为评价水库群闸坝群水量联合调度方案的合理性,基于耗散结构理论建立了流域内各地区间用水相对协同度的评价模型,通过流域内各地区间用水协同度的变化关系,进行水库群闸坝群水量联合调度方案的合理性评价,流程见图5。
![](//html.hanspub.org/file/10-2411066x17_hanspub.png?20211118082424387)
Figure 3. Optimization scheduling model structure diagram
图3. 优化调度模型结构图
![](//html.hanspub.org/file/10-2411066x18_hanspub.png?20211118082424387)
Figure 4. Flow chart of joint optimal operation of water volume of dam and sluice gate of 31 reservoirs in Sihe River basin
图4. 泗河流域31座水库群闸坝群水量联合优化调度计算流程图
![](//html.hanspub.org/file/10-2411066x19_hanspub.png?20211118082424387)
Figure 5. Flow chart of rationality evaluation model
图5. 合理性评价模型流程图
水资源系统的用水相对协调度ui,可由公式(4)计算得出:
(4)
水资源系统用水相对协同度E(p),可由公式(5)计算得出:
(5)
4. 拦河闸坝优化效果评价及建设方案
4.1. 泗河流域水库闸坝群优化效果评价
利用水库闸坝群长系列联合兴利调算成果,以流域内各地区间用水相对协同度、拦河闸坝的蓄水状态及建设位置为评价指标,优化干流闸坝群建设方案。总体上看,流域上游产流量大,拦河闸坝的蓄水状况最佳,中游、下游略逊,特别是兖州以下,区间几乎无产水面积,闸坝来水仅为上游弃水,但部分闸坝因有中水汇入及塌陷区水库相机补源,蓄水状况有所改善,如图6所示。
![](//html.hanspub.org/file/10-2411066x22_hanspub.png?20211118082424387)
Figure 6. Annual average water storage before the optimization of river dam and sluice gate
图6. 拦河闸坝优化前多年平均蓄水情况
4.2. 拦河闸坝优化建设方案
在规划工况下,计算得出泗河流域各地区间用水相对协同度为0.62。同时结合闸坝的控制流域面积、上级蓄水工程状况等因素,初步提出规划拦河闸坝优化方案(如图7所示)。
![](//html.hanspub.org/file/10-2411066x23_hanspub.png?20211118082424387)
Figure 7. Schematic diagram of optimized construction scheme for river sluice and dam
图7. 拦河闸坝优化建设方案示意图
1) 泗源橡胶坝来水量仅为贺庄水库弃水量,来水条件不佳,距离上游中型水库贺庄水库较近,建成后供水功能不强,且距离泗水城区较远生态景观功能弱,综上各因素本次建议将泗源坝去掉。
2) 官园大桥橡胶坝因临泗坝及红旗闸,所以均属泗河干流上调蓄能力较强的水利枢纽,来水条件较为优越,但其建成后供水功能不强,建议去掉。
3) 颜家村橡胶坝来水量仅为陈寨坝弃水量,来水条件不佳,且距离上级已建成的陈寨坝及规划建设张家村坝距离近,建议去掉。
4) 崇文大道橡胶坝下移1 km后,距离塌陷区相机补源水库更近,具有节省工程投资、蓄水量更大的优势,更有利于为塌陷区相机补源水库补水;并且,下游1 km处河底高程相差不大,回水长度仍能至崇文大道处,景观效果更好。
4.3. 拦河闸坝优化后模型成果
4.3.1. 拦河闸坝优化后蓄水情况
结合拦河闸坝优化后多年平均蓄水情况(图8所示)以及梯级闸坝蓄水状况(表1所示),对拦河闸坝优化后多年平均总体蓄水状况分析如下:年内完全蓄满的月份有5.3个月,占44.5%;部分蓄满的有5.3个月,占44.5%;完全蓄不满的有1.4个月,占11%。各闸坝的蓄水状况分析如下:黄阴集闸蓄水状态最佳,完全蓄满的月份有8.3个月;书院坝由于非汛期向外调出水,造成完全蓄满的月份有1.8个月;金口坝因有城市供水任务,虽然年内完全蓄满的月份仅有0.4个月,但是部分蓄满的月份有9.9个,蓄水状态尚可。
![](//html.hanspub.org/file/10-2411066x24_hanspub.png?20211118082424387)
Figure 8. Annual average water storage after the optimization of river dam and sluice gate
图8. 拦河闸坝优化后多年平均蓄水情况
![](Images/Table_Tmp.jpg)
Table 1. Water storage status of cascade dam
表1. 梯级闸坝蓄水状况
4.3.2. 河道内生态景观效果评价
泗河干流闸坝群优化结果如表2所示,河道内多年月均水面面积达2364万m2。泗水县境内、曲阜市境内、兖州市及济宁市区境内干流河道内多年月均水面面积分别为753、463、1148万m2,较现状工况新增水面面积分别为365、233、790万m2,全河道新增水面面积1388万m2。
![](Images/Table_Tmp.jpg)
Table 2. Evaluation of ecological landscape effect in river channel (104 m2)
表2. 河道内生态景观效果评价(万m2)
4.3.3. 地表水资源开发利用率
地表水资源开发利用率结果,如表3所示。对结果分析可知,规划工程布局下泗河干流水库群及闸坝群多年平均向河道外生产用户供水12,454万m3,向外流域调出水量313万m3,计算干流地表水资源开发利用率为35%。其中,泗水县境内、曲阜市境内、兖州市及济宁市区境内地表水资源开发利用率分别为23%、44%、28%。与现状相比,泗水县境内水资源开发利用率略有下降,为中下游地区提供了更多的水资源量,有利于地区间协调发展;兖州市及济宁市区水资源开发利用率增加,增加利用的水资源量主要为煤矿塌陷区水库存蓄的雨洪水资源,以及上中游地区的生态泄水量。
![](Images/Table_Tmp.jpg)
Table 3. Utilization rate of surface water resources development
表3. 地表水资源开发利用率
4.3.4. 县际控制断面下泄水量
规划工程布局下,经优化调度后,得出不同县际控制断面下泄流量如表4所示。泗水曲阜边界、曲阜兖州边界、兖州济宁市区边界多年平均下泄水量分别为9398、6515、5570万m3,较现状分别新增下泄水量938、1212、1121万m3。
![](Images/Table_Tmp.jpg)
Table 4. Intercounty control section drainage
表4. 县际控制断面下泄水量
5. 基于县际控制断面下泄水量的干流水量分配方案研究
水量分配问题涉及到地区与部门多个决策主体,近期与远期多个决策时段,社会、经济、生态多个决策目标,是一个多维、复杂问题。水量分配应遵循以流域为基础,满足生态、公平性、尊重现状、效率、丰增枯减、市场、产水量优先的基本原则。课题组创新性提出基于水库闸坝群交互式优化调度模型,研究制定了流域内县际控制断面下泄水量分配的控制指标,如表5所示。在不同行政地区水量分配方案合理可行的基础上,秉承人水和谐的治水理念,实现区域间雨洪资源的充分协调可持续利用。
![](Images/Table_Tmp.jpg)
Table 5. Control index of water discharge at intercounty section of main stream
表5. 河干流县际断面下泄水量控制指标
6. 结语
本文以流域协同用水为目标,研究构建水库闸坝群交互式优化调度模型,以流域内各地区间用水相对协同度、拦河闸坝蓄水状态、河道生态用水满足程度等为评价指标,研究河道拦河闸坝建设优化方案;研究制定流域内不同行政地区的水量分配方案,对水权水市场建设和最严格水资源管理制度的顺利实施具有重要意义。
参考文献