由于不同水源汇入及分出等的交替作用与影响,多水源、多用户联通构成的多元汇入和分出的输调水系统运行是水量、水质混合及时空演变的过程,为量化不同节点及引出目标的水源构成和水质指标时变状态,研究建立了多元输调水系统动态水指标解析模型。模型以输水系统干网和各汇入、分出口及调蓄工程为研究范围,以沿程实时水量组成、水质因子为研究对象,提出系统多元动态水量与水质指标解析计算方法。以南水北调中线总干渠河南段为示例,应用模型得到干渠各节点、分出口的各水源水量构成、水质因子等指标成果,经分析可为制定科学合理的调度方案、实现预报预调或水量统计核算提供技术支撑,从而提高多水源联合调度的管理水平。 Due to the alternating effects and influences of multiple sources recirculation and flow cross-mixing, multi-source water, multi-user communication, multi-input and separate water transfer system opera-tion is the process of the amount of water, water quality mixing and water quality evolution. The water source composition of different nodes and the extraction targets and the dynamic water quality index need to quantify. An analytical dynamic water index model is developed for multi-source water delivery system in this study. The model takes the main stream, the import and export, storage projects as the research scope, the real time composition of water and the change of quality along the path as the research object. Taking the Henan section of the Middle Route of the South-to-North Water Transfer Project as an example, the application model is used to obtain the results of the water source flow composition and water quality factors of each node and branch of the main canal. It provides technical support for the realization of predictive pre-adjustment or statistical accounting, improving the management level of dispatch operations.
刘德波1,李程纯子1,程德虎2,程曦2,程伊文2
1河南省水利勘测设计研究有限公司,河南 郑州
2南水北调中线干线工程建设管理局,北京
收稿日期:2019年2月19日;录用日期:2019年4月1日;发布日期:2019年4月17日
由于不同水源汇入及分出等的交替作用与影响,多水源、多用户联通构成的多元汇入和分出的输调水系统运行是水量、水质混合及时空演变的过程,为量化不同节点及引出目标的水源构成和水质指标时变状态,研究建立了多元输调水系统动态水指标解析模型。模型以输水系统干网和各汇入、分出口及调蓄工程为研究范围,以沿程实时水量组成、水质因子为研究对象,提出系统多元动态水量与水质指标解析计算方法。以南水北调中线总干渠河南段为示例,应用模型得到干渠各节点、分出口的各水源水量构成、水质因子等指标成果,经分析可为制定科学合理的调度方案、实现预报预调或水量统计核算提供技术支撑,从而提高多水源联合调度的管理水平。
关键词 :多水源,联合调度,动态水指标,解析模型,预报预调
Copyright © 2019 by author(s) and Wuhan University.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
随着经济社会发展和人民生活水平的提高,对水资源质量、保证程度的要求越来越高,生活、生产、生态用水安全和环境改善的需求也日益增强。然而不同区域水资源时空分布不均,一些地区水资源 [
对于多水源、多用户联通构成的多元汇入和分出的输调水系统,各水源输入的水量,或可调入过程是有时空差异的,表现在输入水量、水质不同和过程不同步等方面。系统的需求可以是常态的、也可能是应急的,输出水量用途可能是生活、生产或生态等,以及不同的组合情况。在系统边界约束下,系统的运行就是多元汇入、演变分出的过程。如南水北调中线工程,总干渠输水线路长,极端情况下可能会出现输水中断 [
Hamalainen等,提出多水源调度需要有水量和水质综合分析作为前提 [
本文以跨区域水系干网及各汇入、分出口及调蓄工程为研究范围,以沿程实时水量组成、水质变化为研究对象,建立多元输调水系统动态水指标解析模型,关注各水源与各用户之间的关系,对系统各节点、各用户水质、水量演变进行量化。解决了多水源多用户之间的水量分配问题,明确各用水目标的需求满足程度、水源水量构成及调度的具体效果,为调度方案的分析论证、调整优化、实现预知预调或统计核算提供技术支撑,从而提高工程调度运行的管理水平。
多元输调水系统一般有输入、输出两类因子组成,随着不同节点、不同量级水体的汇入、引(输)出,水指标沿程随时间变化,干流不同节点及各引出口水源水量构成、水质不同。按照输入、输出的工程属性不同,可划分为纯输入、纯输出、既可输入又可输出的调蓄工程等类型,根据系统组成情况,可划分单元进行研究,单元输调水系统可概化为图1所示结构示意图。
图1. 输水系统概化结构示意图
图中1~9为输入、输出节点,Q(i)为第i个输入单元的流量,Oc(j)为沿程j节点的引出流量。 T表示调蓄水库连接。以典型断面或关键节点为计算目标,考虑水流能够充分自然混合,则沿程水量、水质指标主要计算公式如下:
Q u ( j , t + d t ( j ) ) = ∑ i = 1 j − 1 Q ( i , j − 1 , t ) (1)
Q d ( j , t + d t ( j ) ) = Q u ( j , t + d t ( j ) ) − Q c ( j , t ) (2)
Q c ( j , t ) = Q u ( j , t ) ,当 Q 引 ( j , t ) < Q u ( j , t ) (3)
Q c ( j , t ) = Q u ( j , t ) ,当 Q 引 ( j , t ) ≥ Q u ( j , t ) (4)
Q c ( i , j , t + d t ( j ) ) = Q ( i , j − 1 , t ) ⋅ Q c ( j , t ) / Q u ( j , t ) (5)
S u ( i , t + d t ( j ) , k ) = ∑ S ( i , t , k ) ⋅ Q ( i , j − 1 , t ) Q u ( j , t + d t ( j ) ) (6)
式中:i为输入单元;j为沿程节点或断面;t为时段(h);dt为时差(h),k为水质指标类别; S ( i , t , k ) 为t时刻第i个输入单元水质指标; Q u ( j , t ) 为t时刻节点j上游干流流量(m3/s); Q d ( j , t ) 为t时刻节点j下游干流流量(m3/s); Q c ( j , t ) 为t时刻在节点j处引出流量(m3/s); Q 引 ( j , t ) 为t时刻节点j处计划引出流量(m3/s); Q c ( i , j , t ) 为t时刻在节点j处引出的第i个输入单元的流量(m3/s); Q ( i , j , t ) 为t时刻在节点j处第i个输入单元的流量(m3/s); S u ( i , t , k ) 为t时刻节点j处第i个输入单元水质指标。
对于调蓄工程,其调度方式一般有两种类型:一是可以人工控制的,二是自由进出的。前者为可控制调节,即引出、输入水量皆可控制大小;自由进出的调蓄工程是干流水位高于调蓄水库水位时引出,反之则为干流补水。引出、输入流量计算如下:
1) 可控制调节型
按干流的流量Q(j)判别,可设定上边界Q1(j)、下边界Q2(j)条件,且 Q 1 ( j ) ≥ Q 2 ( j ) 。当 Q ( j ) > Q 1 ( j ) 时,引出;当 Q ( j ) < Q 2 ( j ) 时,输入;假定引出及输入能力分别为Q引限(j)、Q补限(j),则:
Q 引 出 ( j ) = ( Q ( j ) − Q 1 ( j ) ) ⋅ K c ( j ) ,按比例引出,且 Q 引 出 ( j ) ≤ Q 引 限 ( j ) (7)
Q 输 入 ( j ) = ( Q 2 ( j ) − Q ( j ) ) ⋅ K b ( j ) ,按比例输入,且 Q 输 入 ( j ) ≤ Q 补 限 ( j ) (6)
调度时可以调节Q1(j)、Q2(j)、Kc(j)、Kb(j)等参数,通过设置以上边界条件满足系统整体运行要求。
2) 自由调节型
自由调节与上、下游水位有关,由水位高的一边向水位低的一边供水。
Q 引 出 ( j ) = f 1 ( H 干 ( j ) , H 库 ( j ) ) ; Q 输 入 ( j ) = f 2 ( H 干 ( j ) , H 库 ( j ) ) (9)
式中:Q引出(j)、Q输入(j)分别为调蓄水库连接节点j处的引出、输入流量;H干(j)、H库(j)分别为节点j处干流和调蓄水库的水位。以上未考虑水体与边界之间水量、水质交换,蒸发、渗漏、降雨、水质挥发性等影响。
以南水北调中线总干渠河南段为示例,除引调丹江水外,考虑沿线周边水库作为供水水源,构建多水源多用户的多元耦合演变系统动态水指标解析模型,第一个输入水源为Q(1),假设流量为2 m3/s;第一个引出目标为Qc(2),计划引出流量为0.5 m3/s;第二个输入水源为Q(2),流量为1 m3/s;第二个引出目标为Qc(4),引出流量为1 m3/s,节点(断面) 1、2、3、4间水流传播时间分别为6、12、9小时。解析模型动态水指标计算主要参数摘录如表1,计算成果如表2、表3所示。
j | 间隔距离(km) | 属性 | 传播时间t(h) | Q(i) (m3/s) | Qc(j) (m3/s) | 控制条件 |
---|---|---|---|---|---|---|
1 | 20 | 控入 | 5 | 2 | 流量 | |
2 | 30 | 控出 | 6 | 0.5 | 流量 | |
3 | 80 | 控入 | 12 | 1 | 流量 | |
4 | 60 | 控出 | 9 | 1 | 流量 | |
5 | 50 | 控出 | 8 | 流量 | ||
6 | 70 | 调蓄 | 11 | 流量 | ||
7 | 40 | 控出 | 7 | 流量 | ||
8 | 90 | 控出 | 15 | 流量 | ||
9 | 15 | 控出 | 4 | 流量 | ||
10 | 110 | 分流 | 18 | 分水比例 |
表1. 典型输水模型计算参数
时间(t) | Q(1)下游水指标 | Qc(2)引水指标 | Qc(2)下游水指标 | ||||||
---|---|---|---|---|---|---|---|---|---|
Q(1)流量 | PH | 溶解氧(mg/L) | Qc(2)流量 | PH | 溶解氧(mg/L) | Qc(2)流量 | PH | 溶解氧(mg/L) | |
1 | 2 | 7 | 16 | 0 | 0 | 0 | 0 | 0 | 0 |
2 | 2 | 7 | 16 | 0 | 0 | 0 | 0 | 0 | 0 |
3 | 2 | 7 | 16 | 0 | 0 | 0 | 0 | 0 | 0 |
4 | 2 | 7 | 16 | 0 | 0 | 0 | 0 | 0 | 0 |
5 | 2 | 7 | 16 | 0 | 0 | 0 | 0 | 0 | 0 |
6 | 2 | 7 | 16 | 0 | 0 | 0 | 0 | 0 | 0 |
7 | 2 | 7 | 16 | 0.5 | 7 | 16 | 1.5 | 7 | 16 |
8 | 2 | 7 | 16 | 0.5 | 7 | 16 | 1.5 | 7 | 16 |
9 | 2 | 7 | 16 | 0.5 | 7 | 16 | 1.5 | 7 | 16 |
10 | 2 | 7 | 16 | 0.5 | 7 | 16 | 1.5 | 7 | 16 |
11 | 2 | 7 | 16 | 0.5 | 7 | 16 | 1.5 | 7 | 16 |
12 | 2 | 7 | 16 | 0.5 | 7 | 16 | 1.5 | 7 | 16 |
… | … | … | … | … | … | … | … | … | … |
表2. 水指标计算成果(一)
时间(t) | Q(2)下游水指标 | Qc(4)引水指标 | Qc(4)下游水指标 | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Q(1)流量 | Q(2)流量 | PH | 溶解氧(mg/L) | Q(1)流量 | Q(2)流量 | PH | 溶解氧(mg/L) | Q(1)流量 | Q(2)流量 | PH | 溶解氧(mg/L) | |
1 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
2 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
3 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
4 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
5 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
7 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
8 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
9 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
10 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
11 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
12 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
13 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
14 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
15 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
16 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
17 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
---|---|---|---|---|---|---|---|---|---|---|---|---|
18 | 0 | 1 | 7.5 | 15 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
19 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
20 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
21 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
22 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
23 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
24 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
25 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
26 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
27 | 1.5 | 1 | 7.2 | 15.6 | 0 | 1 | 7.5 | 15 | 0 | 0 | 0 | 0 |
28 | 1.5 | 1 | 7.2 | 15.6 | 0.6 | 0.4 | 7.2 | 15.6 | 0.9 | 0.6 | 7.2 | 15.6 |
29 | 1.5 | 1 | 7.2 | 15.6 | 0.6 | 0.4 | 7.2 | 15.6 | 0.9 | 0.6 | 7.2 | 15.6 |
30 | 1.5 | 1 | 7.2 | 15.6 | 0.6 | 0.4 | 7.2 | 15.6 | 0.9 | 0.6 | 7.2 | 15.6 |
… | … | … | … | … | … | … | … | … | … | … | … | … |
表3. 水指标计算成果(二)
考虑输调水系统结构、影响因素,根据水源数据和边界参数,采用本文研究模型进行计算,得到各节点(断面)干渠上、下游及分出口的各水源流量构成、水质因子变化等成果,对各节点(断面)指标进行分析是合理正确的。以下仅列出水源和节点流量、ph值、溶解氧等指标,部分成果摘入表1、表2,选取断面4实际引出Qc(4)的各水源构成及水质、引出Qc(4)后干渠下游水指标随时间变化成果,如图2、图3所示。
图2. Qc(4)实际引出水指标
图3. 断面4下游水指标
本文对多水源多用户输调水系统,构建了多水源动态水量、水质指标解析模型,通过解析计算,得到节点及用户不同水源水量构成、水质的时变状态。以南水北调中线联通沿线已建水库、纯调蓄工程构成的多水源输调水系统为例进行模拟验证。模型解决了多水源输水干渠、生态河道、综合利用河网水系等多水源多用户输调水系统水量分配、水价分担计算等问题,为进行水量和水质指标预测、调度效果评价、方案优化及调度管理提供技术支撑。
长距离输水及干流调节作用明显时,模型可加入相关变量及参数,考虑更多的影响因素。
刘德波,李程纯子,程德虎,程 曦,程伊文. 多水源联合调度动态水指标解析模型研究Analytical Dynamic Water Index Modeling for Dispatch Operations of Multi-Sources Water[J]. 水资源研究, 2019, 08(02): 177-184. https://doi.org/10.12677/JWRR.2019.82021