合理、高效地分配水资源成为水资源领域亟待解决的问题,水资源优化配置是解决该问题的重要手段,如何有效协调不同用水类别之间的竞争关系,并考虑用水需求和供水能力在时间尺度上的变化情况是难点。本文提出一种流域水资源配置的综合缺水量最小优化数学模型,即通过定义权重系数建立调度期综合缺水量最小目标函数,并在建立约束条件时建立水库自身蒸发渗漏损失以及退水流量的计算公式,最终联立目标函数与约束条件进行求解,其在南盘江流域的实际应用表明:1) 权重系数的设置有效协调了各种用水类别间的竞争关系,实现了水资源在用水配置上的最优分配,具有很强的适用性与应用价值。2) 供水条件时间变化特性会使得系统内出现结构性缺水,在水库蓄水期其影响较强。3) 受调度过程影响相同水库的年蒸发渗漏损失可能不同,且损失水量占比明显,不可忽略。 Rational and efficient allocation of water resources has become an urgent problem, and the optimal water resources allocation is an important means to solve this problem. How to effectively coordinate the competitive relationship between different water use categories and consider the changes of water demand and water supply capacity in time scale is the difficult point. This paper proposes an optimal water resources allocation model to minimize the comprehensive water shortage during the operation period. The weight coefficient, calculation formulas for the evaporation leakage loss, the return flow of the reservoir, and the constraint conditions are established. The application results in the Nanpanjiang River basin shows that: 1) The setting of weight coefficient effectively coordinates the competitive relationship among various water use categories, realizes the optimal allocation of water resources in water use allocation, and has strong applicability and application value. 2) The time-varying characteristics of water supply conditions will lead to structural water shortage in the system, which has a strong impact during the impounding period of the reservoir. 3) Affected by the operation process, the annual evaporation leakage loss of the same reservoir may be different, and the proportion of water loss is obvious, which cannot be ignored.
合理、高效地分配水资源成为水资源领域亟待解决的问题,水资源优化配置是解决该问题的重要手段,如何有效协调不同用水类别之间的竞争关系,并考虑用水需求和供水能力在时间尺度上的变化情况是难点。本文提出一种流域水资源配置的综合缺水量最小优化数学模型,即通过定义权重系数建立调度期综合缺水量最小目标函数,并在建立约束条件时建立水库自身蒸发渗漏损失以及退水流量的计算公式,最终联立目标函数与约束条件进行求解,其在南盘江流域的实际应用表明:1) 权重系数的设置有效协调了各种用水类别间的竞争关系,实现了水资源在用水配置上的最优分配,具有很强的适用性与应用价值。2) 供水条件时间变化特性会使得系统内出现结构性缺水,在水库蓄水期其影响较强。3) 受调度过程影响相同水库的年蒸发渗漏损失可能不同,且损失水量占比明显,不可忽略。
水资源配置,缺水量最小,权重系数,蒸发渗漏
Zhenyu Mu1, Anqi Tan1, Xueshan Ai1, Senlin Chen1, Tianqing Li2
1School of Water Resources and Hydropower Engineering, Wuhan University, Wuhan Hubei
2Kunming Engineering Corporation Limited, Kunming Yunnan
Received: May 1st, 2022; accepted: Jun. 14th, 2022; published: Jun. 29th, 2022
Rational and efficient allocation of water resources has become an urgent problem, and the optimal water resources allocation is an important means to solve this problem. How to effectively coordinate the competitive relationship between different water use categories and consider the changes of water demand and water supply capacity in time scale is the difficult point. This paper proposes an optimal water resources allocation model to minimize the comprehensive water shortage during the operation period. The weight coefficient, calculation formulas for the evaporation leakage loss, the return flow of the reservoir, and the constraint conditions are established. The application results in the Nanpanjiang River basin shows that: 1) The setting of weight coefficient effectively coordinates the competitive relationship among various water use categories, realizes the optimal allocation of water resources in water use allocation, and has strong applicability and application value. 2) The time-varying characteristics of water supply conditions will lead to structural water shortage in the system, which has a strong impact during the impounding period of the reservoir. 3) Affected by the operation process, the annual evaporation leakage loss of the same reservoir may be different, and the proportion of water loss is obvious, which cannot be ignored.
Keywords:Water Resources Allocation, Minimum Water Shortage, Weight Coefficient, Evaporation Leakage
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/
受人口增长以及气候变化影响,世界范围内发生的缺水以及水质安全问题与日俱增 [
水资源优化配置是指在流域或特定的区域范围内,以水资源的可持续利用和经济社会可持续发展为目标,对多种可利用水资源在区域间和各用水部门间进行的合理调配 [
基于上述背景,本文提出一种流域水资源配置的综合缺水量最小优化数学模型,即通过定义权重系数建立调度期综合缺水量最小目标函数,并在建立约束条件时建立水库自身蒸发渗漏损失以及退水流量的计算公式,最终联立目标函数与约束条件进行求解,其在南盘江流域的实际应用表明,本模型有效协调了流域内不同用水类别间的竞争关系,同时兼顾供水节点自身的供水能力变化情况,具有较强的应用价值。
设 y ( j ) 表示第j用水节点所属的用水类别, y ( j ) ∈ [ 1 , 6 ] 。对各用水类别缺水量设置等差权重系数(表1),优化准则为调度期综合缺水量最小,即目标函数为:
min ∑ t = 1 T ∑ j = 1 N Y S [ α y ( j ) ( D j , t − ∑ i = 1 N G S ( r i , j G Y S i , j , t ) Δ T t ) ] (1)
式中:T为调度期长度; N Y S 为水资源系统内的用水节点数; N G S 为水资源系统内的供水节点数; D j , t 为第j用水节点在第t时段的需水流量,m3/s; r i , j G Y 为第i供水节点与第j用水节点间的拓扑关系,当第i供水节点向第j用水节点供水时其值为1,反之为0; S i , j , t 为第i供水节点向第j用水节点的供水流量,m3/s; α y ( j ) 为第j用水节点所属类别对应第 y ( j ) 用水级别的重要性系数,其对应关系见表1。
用水节点所属类别 | 用水级别 | 权重系数 α n |
---|---|---|
重要生活 | 1_SH | 6 |
一般生活 | 0_SH | 5 |
重要工业 | 1_GY | 4 |
一般工业 | 0_GY | 3 |
重要农业 | 1_NY | 2 |
一般农业 | 0_NY | 1 |
表1. 用水类别的重要性系数
约束条件分为两类,一是工程水量平衡约束,二是决策变量边界约束。
由于不同节点的入流、出流等条件各不相同,本文针对不同节点类型分别建立水量平衡关系约束。
1) 水库、水闸的入库流量为:
Q i , t i n = Q i , t J + ∑ j ∈ Ω ¯ i G G Q j , t o u t + ∑ j ∈ Ω ¯ i Y G Q j , t U (2)
式中: Q i , t i n 为第i水库(闸) t时段的入库流量,m3/s; Q i , t J 为第i水库(闸) t时段上游的区间流量,m3/s; Q j , t o u t 为第i水库(闸) t时段上游第j ( j ∈ Ω ¯ i G G )供水节点的出流流量,m3/s; Q j , t U 为第j用水节点t时段的退水流量,m3/s。
如果给定第j用水节点的退水系数为 α j ,则退水流量 Q j , t U = α j Q j , t Y ;如果考虑滞时,比如2个时段滞时(分别为 β j , 1 、 β j , 2 ),则相邻2个时段的退水流量分别为 Q j , t U = α j β j , 1 Q j , t Y 、 Q j , t + 1 U = α j β j , 2 Q j , t Y 。
2) 供水节点水量平衡约束
水库水量平衡约束由五部分构成,一是水库入流量,二是水库出流量,三是水库供水量,四是水库蓄水变化量,五是渗漏及蒸发损失量,水量平衡约束如下:
{ V i , t = V i , t − 1 + ( Q i , t i n − Q i , t o u t − ∑ j ∈ Ω _ i G Y Q i , j , t D G − Q i , t S L − Q i , t Z F ) ⋅ Δ T t × 10 − 4 Q i , t S L = β i V i , t − 1 + V i , t 2 Q i , t Z F = γ i , t A i , t − 1 + A i , t 2 (3)
式中: V i , t − 1 、 V i , t 分别为第i水库t时段始、末兴利蓄水量,104 m3; Q i , t o u t 为第i水库t时段的出库流量,m3/s; Q i , j , t D G 为t时段第i水库供给第j用水节点的供水流量,m3/s; Q i , t S L 、 Q i , t Z F 为第i水库t时段的渗漏损失、蒸发损失,m3/s; β i 、 γ i , t 为第i水库的渗漏损失系数、蒸发损失系数,10−4 s−1; A i , t − 1 、 A i , t 分别为第i水库t时段始、末水库水面面积,104 m2。
水闸水量平衡约束:不需考虑蓄水变化量以及渗漏及蒸发损失量,流向下游的出闸流量加供水流量应该等于入闸总流量:
Q i o u t + ∑ j ∈ Ω _ i G Y Q i , j , t D G = Q i , t i n (4)
其中水回用工程水量平衡:流向下游河道的流量加供水流量应该等于用水流量乘以回收系数:
∑ j ∈ Ω _ i G G Q j , t o u t + ∑ j ∈ Ω _ i G Y Q i , j , t D G = ∑ j ∈ Ω ¯ i Y G λ j Q j , t Y (5)
式中: λ j 为第j用水节点的中水回用比例。
3) 用水节点供需水量平衡约束
用水节点水量平衡约束为:
Q i , t Y = ∑ j ∈ Ω ¯ i G Y Q j , i , t D G ≤ D i , t (6)
1) 水库库容上下边界约束为:
0 ≤ V i , t ≤ V i , t max (7)
式中: V i , t max 为第i水库t时段最大兴利库容,104 m3。
2) 水库及水闸下泄流量上下边界约束为:
Q i , t min ≤ Q i , t o u t ≤ Q i , t max (8)
式中: Q i , t min 为第i水库或水闸t时段的最小下泄流量,即生态流量,m3/s; Q i , t max 为第i水库t时段允许的最大下泄流量,m3/s。
3) 单一管道供给流量边界约束为:
0 ≤ Q i , j , t D G ≤ Q i , j D G M A X (9)
式中: Q i , j D G M A X 为第i供水节点供给第j用水节点的流量上限,m3/s。
4) 复合管道复合流量边界约束为:
0 ≤ ∑ i 、 j ∈ Ω k G D Q i , j , t D G ≤ Q k ¯ D G M A X (10)
式中: Ω k G D 为第k复合管道包含的所有供水–用水节点序号集合; Q k ¯ D G M A X 为第k复合管道的输水流量上限,m3/s。
式(1)~(10)所组成的优化数学模型属于典型的线性规划模型,可以利用相关商业软件进行求解,且能求得全局最优解。
南盘江发源于云南省曲靖市乌蒙山,经贵州省兴义市至广西红水河,是珠江的组成部分之一,其上游用水可分为农业、城镇、工业三类,为满足用水需求,需从临近的黑滩河流域进行调水。黑滩河为牛栏江右岸一级支流,属金沙江流域牛栏江水系,河道全长35.8 km,其中黑滩河水库为流域内最大水库,水库总库容13,500 m3,坝高52 m,输水工程全长约18.728 km,是主要供水节点之一,其蓄水对保障枯水年供水,减少供水破坏时间,降低破坏深度具有重要影响,对应流域网络节点图以及各类用水、供水节点如图1所示:
图1. 南盘江流域及水源调出区网络节点图
应用以上优化数学模型,采用线性规划方法进行求解,可得到各用水类别的空间分配结果,为便于比较,记无权重系数模型为A1,有权重系数(本方法)模型为A2,其计算结果见表2:
用水节点名称 | 用水级别 | 总需水量 (万m³) | 模型A1 | 模型A2 | ||
---|---|---|---|---|---|---|
总缺水量(万m³) | 缺水率(%) | 总缺水量(万m³) | 缺水率(%) | |||
黑滩河水库以上片–农业 | 0_NY | 784.913 | 49.02 | 6.25 | 364.18 | 46.40 |
块所坝片–农业 | 0_NY | 177.202 | 4.32 | 2.44 | 131.12 | 74.00 |
花山灌片–农业 | 0_NY | 2095.120 | 0 | 0 | 0 | 0 |
九龙灌左片–农业 | 0_NY | 346.031 | 0 | 0 | 0 | 0 |
九龙灌右片–农业 | 0_NY | 151.661 | 0 | 0 | 0 | 0 |
东风灌左片–农业 | 0_NY | 1320.689 | 0 | 0 | 0 | 0 |
东风灌右片–农业 | 0_NY | 482.587 | 0 | 0 | 0 | 0 |
菱角灌片–农业 | 1_NY | 202.955 | 14.29 | 7.04 | 35.82 | 17.65 |
白浪灌片–农业 | 1_NY | 1809.402 | 0 | 0 | 0 | 0 |
盘江集镇片–工业 | 0_GY | 197.376 | 26.71 | 13.53 | 13.36 | 6.77 |
花山街道片–工业 | 0_GY | 4433.029 | 369.72 | 8.34 | 369.73 | 8.34 |
沾益中心城区片–工业 | 1_GY | 2606.033 | 330.12 | 12.67 | 0 | 0 |
黑滩河水库以上片–农村 | 0_SH | 114.874 | 9.45 | 8.23 | 0 | 0 |
块所坝片–农村 | 0_SH | 114.874 | 9.17 | 7.99 | 0 | 0 |
盘江集镇片–城镇 | 1_SH | 403.568 | 61.53 | 15.25 | 0 | 0 |
花山街道片–城镇 | 1_SH | 301.291 | 22.71 | 7.54 | 0 | 0 |
沾益中心城区片–城镇 | 1_SH | 2774.470 | 0 | 0 | 0 | 0 |
总计 | - | 18316.075 | 897.04 | - | 914.21 | - |
表2. 有无权重系数时各模型计算结果
由上表可以得出,尽管在总缺水量上模型A2 (914.21万m³)比模型A1 (897.04万m³)多缺水17.17万m³(1.91%),但是,从破坏的用水类别数进行分析,模型A1的6个用水类别中有10个用水节点存在缺水,而模型A2只有3个用水类别中的5个用水节点存在缺水,且模型A2中权重系数较大的工业用水和生活用水较于模型A1缺水比例大幅度减少,代价为权重系数较小的农业用水缺水量增多。
结构性缺水指由水资源系统供需结构引起的,在用水级别较低的用水户还没全部破坏时,用水级别相对较高的用水户已经破坏的情况 [
时段 | 盘江集镇片–工业 | 花山街道片–工业 | 菱角灌片–农业 | |||||
---|---|---|---|---|---|---|---|---|
缺水流量 | 供水节点 | 缺水流量 | 供水节点 | 缺水流量 | 供水节点 | |||
白浪水库 | 黑滩河水库 | 花山水库 | 黑滩河水库 | 黑滩河水库 | ||||
1 | 0 | 0.0499 | 0 | 0 | 0.0557 | 1.3147 | 0 | 0.0804 |
2 | 0 | 0.0552 | 0 | 0.0111 | 0.0693 | 1.4368 | 0 | 0.11324 |
3 | 0 | 0.0499 | 0 | 0 | 0 | 1.3704 | 0 | 0.1347 |
4 | 0 | 0.0515 | 0 | 0 | 0 | 1.4161 | 0 | 0.1135 |
5 | 0 | 0 | 0.0499 | 0 | 0 | 1.3704 | 0 | 0.1425 |
6 | 0 | 0 | 0.0515 | 0 | 1.1268 | 0.2892 | 0 | 0.0377 |
7 | 0 | 0 | 0.0499 | 0 | 0.6124 | 0.7580 | 0 | 0.0084 |
8 | 0 | 0 | 0.0499 | 0 | 0.2794 | 1.0910 | 0 | 0.0080 |
9 | 0 | 0.0515 | 0 | 0 | 0.3739 | 1.0421 | 0 | 0.0014 |
10 | 0 | 0.0499 | 0 | 0.00002 | 0.2695 | 1.1008 | 0 | 0 |
11 | 0 | 0 | 0.0515 | 0 | 0.5662 | 0.8499 | 0.0831 | 0 |
12 | 0.0499 | 0 | 0 | 1.3704 | 0 | 0 | 0.05335 | 0 |
表3. 模型A2中存在缺水的一般工业和重要农业用水节点供水情况(m³/s)
由表3可得缺水时段主要集中在11月与12月,这是由于为了将水库续至目标库容,向上述用水节点供水的白浪水库、黑滩河水库与花山水库降低了供水流量而转向用于蓄水(水库详细调度过程图见图2~4),因此上述节点在11月、12月的缺水是合理的。
图2. 白浪水库模型A2的调度过程
图3. 花山水库模型A2调度过程
图4. 模型A2黑滩河水库调度过程
水资源配置系统内各水库的蒸发渗漏损失水量见表4:
模型A1 | 模型A2 | 差值 | |
---|---|---|---|
白浪 | 118.91 | 127.96 | 9.05 |
花山 | 463.56 | 463.26 | −0.30 |
水晶地 | 13.39 | 13.39 | 0.00 |
清水沟 | 45.92 | 45.92 | 0.00 |
九龙 | 0.00 | 0.00 | 0.00 |
姑娘桥 | 26.08 | 26.19 | 0.11 |
哈盆沟 | 5.28 | 6.28 | 1.00 |
鱼洞 | 40.12 | 42.11 | 1.99 |
三里桥 | 11.14 | 14.32 | 3.17 |
土地口 | 5.74 | 5.59 | −0.16 |
黑滩河 | 715.08 | 717.38 | 2.30 |
合计 | 1445.23 | 1462.40 | 17.17 |
表4. 水库年蒸发渗漏损失水量对比(104 m³)
可见模型A2 (1462.40万m³)比模型A1 (1445.23万m³)多损失17.17万m³,这是由水库蓄(或泄和供)水过程不同而产生的,其中白浪水库调度过程与三里桥水库调度过程见图5。
图5. 白浪水库与三里桥水库调度过程
对各用水类别缺水量设置不同权重系数(见表5),在三个模型下分别进行水资源优化配置的计算结果见表6。
用水级别 | 权重系数α | ||
---|---|---|---|
模型B1 | 模型B2 | 模型B3 | |
1_SH | 6 | 100 | 100 |
0_SH | 5 | 100 | 100 |
1_GY | 4 | 100 | 20 |
0_GY | 3 | 100 | 10 |
1_NY | 2 | 10 | 2 |
0_NY | 1 | 1 | 1 |
表5. 不同模型权重系数
用水节点编码单位 | 用水 级别 | 总需水量 (万m³) | 模型B1 | 模型B2 | 模型B3 | |||
---|---|---|---|---|---|---|---|---|
总缺水量 (万m³) | 缺水率(%) | 总缺水量 (万m³) | 缺水率(%) | 总缺水量 (万m³) | 缺水率(%) | |||
黑滩河水库以上片–农业 | 0_NY | 784.913 | 364.18 | 46.40 | 364.18 | 46.40 | 364.18 | 46.40 |
块所坝片–农业 | 0_NY | 177.202 | 131.12 | 74.00 | 131.12 | 74.00 | 131.12 | 74.00 |
花山灌片–农业 | 0_NY | 2095.120 | 0 | 0 | 0 | 0 | 0 | 0 |
九龙灌左片–农业 | 0_NY | 346.031 | 0 | 0 | 0 | 0 | 0 | 0 |
九龙灌右片–农业 | 0_NY | 151.661 | 0 | 0 | 0 | 0 | 0 | 0 |
东风灌左片–农业 | 0_NY | 1320.689 | 0 | 0 | 0 | 0 | 0 | 0 |
东风灌右片–农业 | 0_NY | 482.587 | 0 | 0 | 0 | 0 | 0 | 0 |
菱角灌片–农业 | 1_NY | 202.955 | 35.82 | 17.65 | 35.82 | 17.65 | 35.82 | 17.65 |
白浪灌片–农业 | 1_NY | 1809.402 | 0 | 0 | 0 | 0 | 0 | 0 |
盘江集镇片–工业 | 0_GY | 197.376 | 13.36 | 6.77 | 13.36 | 6.77 | 13.36 | 6.77 |
花山街道片–工业 | 0_GY | 4433.029 | 369.73 | 8.34 | 191.03 | 4.31 | 369.73 | 8.34 |
沾益中心城区片–工业 | 1_GY | 2606.033 | 0 | 0 | 162.25 | 6.23 | 0 | 0 |
黑滩河水库以上片–农村 | 0_SH | 114.874 | 0 | 0 | 7.24 | 6.30 | 0 | 0 |
块所坝片–农村 | 0_SH | 114.874 | 0 | 0 | 9.17 | 7.99 | 0 | 0 |
盘江集镇片–城镇 | 1_SH | 403.568 | 0 | 0 | 0 | 0 | 0 | 0 |
花山街道片–城镇 | 1_SH | 301.291 | 0 | 0 | 0 | 0 | 0 | 0 |
沾益中心城区片–城镇 | 1_SH | 2774.470 | 0 | 0 | 0 | 0 | 0 | 0 |
总计 | - | 18316.075 | 914.21 | - | 914.18 | 914.21 |
表6. 权重系数不同时各模型计算结果
由表6可以分析,模型B1 (914.21万m³)和模型B3缺水量相同,且两者均比模型B2 (914.18万m³)多缺水0.03万m³(0.003%),而模型B1和模型B3只有3个用水类别中的5个用水节点存在缺水。模型B2由于生活类和工业类用水节点的权重系数相同,从而将原本在模型B1和模型B3中只存在于一般工业节点(花山街道片–工业)的缺水量分摊到了重要工业和一般生活用水节点中。
从蒸发渗漏损失水量方面进行比较,模型B1和模型B3 (1462.40万m³)比模型B2 (1462.37万m³)多损失0.03万m³(见表7),这是由水库蓄水过程不同造成的水量损失。
模型B1 (B3) | 模型B2 | 差值 | |
---|---|---|---|
白浪 | 127.96 | 127.96 | 0 |
花山 | 463.26 | 463.58 | 0.32 |
水晶地 | 13.39 | 13.39 | 0 |
清水沟 | 45.92 | 45.92 | 0 |
九龙 | 0.00 | 0.00 | 0 |
姑娘桥 | 26.19 | 26.19 | 0 |
哈盆沟 | 6.28 | 6.28 | 0 |
鱼洞 | 42.11 | 42.11 | 0 |
三里桥 | 14.32 | 14.32 | 0 |
土地口 | 5.59 | 5.59 | 0 |
黑滩河 | 717.38 | 717.04 | −0.34 |
合计 | 1462.40 | 1462.37 | −0.03 |
表7. 不同权重系数条件下水库年蒸发渗流损失水量(104 m³)
本文首先定义权重系数建立调度期综合缺水量最小目标函数,并在建立约束条件时建立水库自身蒸发渗漏损失以及退水流量的计算公式,最终联立目标函数与约束条件进行求解,通过在南盘江流域的实际应用,可得出如下结论:
1) 权重系数的设置有效协调了各种用水类别间的竞争关系,实现了水资源在用水配置上的最优分配,具有很强的适用性与应用价值。
2) 供水条件时间变化特性会使得系统内出现结构性缺水,在水库蓄水期其影响较强。
3) 受调度过程影响相同水库的年蒸发渗漏损失可能不同,且损失水量占比明显,不可忽略。
云南省重点研发计划“水资源高效开发利用和水环境保护关键技术与管理决策支持系统研究”(2019BC002);国家自然科学基金项目(51779177)。
穆振宇,谭安琪,艾学山,陈森林,李天庆. 综合缺水量最小的水资源优化配置模型及应用Optimal Water Resources Allocation Model of Minimum Comprehensive Water Shortage and Its Application[J]. 水资源研究, 2022, 11(03): 274-285. https://doi.org/10.12677/JWRR.2022.113030