1. 引言
随着经济社会发展和人民生活水平的提高,对水资源质量、保证程度的要求越来越高,生活、生产、生态用水安全和环境改善的需求也日益增强。然而不同区域水资源时空分布不均,一些地区水资源 [1] [2] 供需矛盾突出,总量不够、水质不满足要求,或过程不匹配等,成为制约社会发展的主要因素。为实现空间均衡、人水和谐、健康持续发展,改善供求关系,提高供水保障,通过水系联通 [3] 、多水源联合调度进行水源调配是未来解决用水安全问题、提高抗风险能力的有效途径。
对于多水源、多用户联通构成的多元汇入和分出的输调水系统,各水源输入的水量,或可调入过程是有时空差异的,表现在输入水量、水质不同和过程不同步等方面。系统的需求可以是常态的、也可能是应急的,输出水量用途可能是生活、生产或生态等,以及不同的组合情况。在系统边界约束下,系统的运行就是多元汇入、演变分出的过程。如南水北调中线工程,总干渠输水线路长,极端情况下可能会出现输水中断 [4] [5] ,为保障供水安全、平稳,沿线已建水库、调蓄工程联通总干渠供水,即构成多水源、多用户输调水系统。系统调度方案选择、效果评价、调度的有效性、用户的满足程度、水源分配及水价机制等,是需要解决的重要问题。
Hamalainen等,提出多水源调度需要有水量和水质综合分析作为前提 [6] 。张静等以供水调度系统成本最小为目标构建了地表水源、地下水源、外来水源联合调度模型 [7] 。章燕喃等以北京市典型年供水安全和弃水量最小为目标,建立了多水源联合调度模型 [8] 。潘莉等也对北京市多水源水量调度模型进行了研究 [9] 。蒋任飞等建立了南水北调中线调水多目标供水水量分配优化模型 [10] 。宋健峰和殷建军提出了南水北调供水断面水价成本折算分摊法 [11] 。彭辉等建立了岗南、黄壁庄、王快和西大洋水库向南水北调中线应急供水联合调度预报模型 [12] 。目前,大多数文献对于多水源输调水系统的研究,主要是针对一定区域水资源条件,从水源、需求两方面进行的配置模型研究,假定水源条件分析得到典型供水过程,或是针对单水源多用水户供水系统的研究。总体上是设置需求,研究供水方案,是基于统计分析、典型计算提出可能性配置成果。而针对多水源多用户跨区域输调水系统的不同水源与各用户之间动态水量分配关系、水质指标演变的研究较少。
本文以跨区域水系干网及各汇入、分出口及调蓄工程为研究范围,以沿程实时水量组成、水质变化为研究对象,建立多元输调水系统动态水指标解析模型,关注各水源与各用户之间的关系,对系统各节点、各用户水质、水量演变进行量化。解决了多水源多用户之间的水量分配问题,明确各用水目标的需求满足程度、水源水量构成及调度的具体效果,为调度方案的分析论证、调整优化、实现预知预调或统计核算提供技术支撑,从而提高工程调度运行的管理水平。
2. 多水源联合调度动态水指标解析模型
多元输调水系统一般有输入、输出两类因子组成,随着不同节点、不同量级水体的汇入、引(输)出,水指标沿程随时间变化,干流不同节点及各引出口水源水量构成、水质不同。按照输入、输出的工程属性不同,可划分为纯输入、纯输出、既可输入又可输出的调蓄工程等类型,根据系统组成情况,可划分单元进行研究,单元输调水系统可概化为图1所示结构示意图。
![](//html.hanspub.org/file/9-2410740x10_hanspub.png)
Figure 1. Schematic diagram of generalized structure of water conveyance system
图1. 输水系统概化结构示意图
图中1~9为输入、输出节点,Q(i)为第i个输入单元的流量,Oc(j)为沿程j节点的引出流量。 T表示调蓄水库连接。以典型断面或关键节点为计算目标,考虑水流能够充分自然混合,则沿程水量、水质指标主要计算公式如下:
(1)
(2)
,当
(3)
,当
(4)
(5)
(6)
式中:i为输入单元;j为沿程节点或断面;t为时段(h);dt为时差(h),k为水质指标类别;
为t时刻第i个输入单元水质指标;
为t时刻节点j上游干流流量(m3/s);
为t时刻节点j下游干流流量(m3/s);
为t时刻在节点j处引出流量(m3/s);
为t时刻节点j处计划引出流量(m3/s);
为t时刻在节点j处引出的第i个输入单元的流量(m3/s);
为t时刻在节点j处第i个输入单元的流量(m3/s);
为t时刻节点j处第i个输入单元水质指标。
对于调蓄工程,其调度方式一般有两种类型:一是可以人工控制的,二是自由进出的。前者为可控制调节,即引出、输入水量皆可控制大小;自由进出的调蓄工程是干流水位高于调蓄水库水位时引出,反之则为干流补水。引出、输入流量计算如下:
1) 可控制调节型
按干流的流量Q(j)判别,可设定上边界Q1(j)、下边界Q2(j)条件,且
。当
时,引出;当
时,输入;假定引出及输入能力分别为Q引限(j)、Q补限(j),则:
,按比例引出,且
(7)
,按比例输入,且
(6)
调度时可以调节Q1(j)、Q2(j)、Kc(j)、Kb(j)等参数,通过设置以上边界条件满足系统整体运行要求。
2) 自由调节型
自由调节与上、下游水位有关,由水位高的一边向水位低的一边供水。
;
(9)
式中:Q引出(j)、Q输入(j)分别为调蓄水库连接节点j处的引出、输入流量;H干(j)、H库(j)分别为节点j处干流和调蓄水库的水位。以上未考虑水体与边界之间水量、水质交换,蒸发、渗漏、降雨、水质挥发性等影响。
3. 案例模拟
3.1. 解析模拟示例
以南水北调中线总干渠河南段为示例,除引调丹江水外,考虑沿线周边水库作为供水水源,构建多水源多用户的多元耦合演变系统动态水指标解析模型,第一个输入水源为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所示。
![](Images/Table_Tmp.jpg)
Table 1. Calculation index of typical water conveyance models
表1. 典型输水模型计算参数
![](Images/Table_Tmp.jpg)
Table 2. The results of water index (1)
表2. 水指标计算成果(一)
![](Images/Table_Tmp.jpg)
Table 3. The results of water index (2)
表3. 水指标计算成果(二)
3.2. 计算结果及分析
考虑输调水系统结构、影响因素,根据水源数据和边界参数,采用本文研究模型进行计算,得到各节点(断面)干渠上、下游及分出口的各水源流量构成、水质因子变化等成果,对各节点(断面)指标进行分析是合理正确的。以下仅列出水源和节点流量、ph值、溶解氧等指标,部分成果摘入表1、表2,选取断面4实际引出Qc(4)的各水源构成及水质、引出Qc(4)后干渠下游水指标随时间变化成果,如图2、图3所示。
![](//html.hanspub.org/file/9-2410740x36_hanspub.png)
Figure 2. The real time water index of Qc(4)
图2. Qc(4)实际引出水指标
![](//html.hanspub.org/file/9-2410740x37_hanspub.png)
Figure 3. The downstream water index of Section 4
图3. 断面4下游水指标
4. 结语
本文对多水源多用户输调水系统,构建了多水源动态水量、水质指标解析模型,通过解析计算,得到节点及用户不同水源水量构成、水质的时变状态。以南水北调中线联通沿线已建水库、纯调蓄工程构成的多水源输调水系统为例进行模拟验证。模型解决了多水源输水干渠、生态河道、综合利用河网水系等多水源多用户输调水系统水量分配、水价分担计算等问题,为进行水量和水质指标预测、调度效果评价、方案优化及调度管理提供技术支撑。
长距离输水及干流调节作用明显时,模型可加入相关变量及参数,考虑更多的影响因素。
参考文献