1. 引言
随着社会经济的高速发展,人口数量的不断增长,城镇化进程的加快,用水需求量日益增加且供水保证率越高,导致水资源的供需矛盾愈发尖锐,使得生态环境保护和经济发展面临挑战 [1] [2] 。合理分配有限的水资源,提高水资源利用效率,实现生态与社会经济协调、可持续发展已变得十分重要。水资源配置可对水资源进行重新分配与布局,是目前解决水资源供需矛盾问题的主要调控手段之一 [3],已成为国内外大量学者重点研究的课题。为与水资源管理部门的日常运行管理相适应,现有水资源配置研究常以日、月、年等单一时间尺度进行。甘治国等在年时间尺度条件下,基于北京市水资源利用现状,提出并开发了基于规则控制的北京市水资源配置模拟模型 [4];刘德地等在月时间尺度条件下,在分析咸潮影响区水资源特点的基础之上,以珠海市为例,建立了成潮影响区的水资源优化配置模型 [5];王维平等在月时间尺度条件下,结合山东省莱州滨海平原地区海水入侵问题,以生态需水为约束条件、最大水资源开发能力为目标,构建起生态型水资源优化配置模型 [6];田晶等在月时间尺度条件下,以社会、经济和生态效益为目标,构建了汉江流域中下游地区水资源多目标优化配置模型 [7];樊尔兰等在旬时间尺度条件下,分析了西安市黑河水库水温、水质与水库取水建筑物之间的关系,建立了分层型水库水量水质综合优化调度数学模型 [8];张守平等在月时间尺度条件下,以改进的“三次平衡”思想为指导,在湟水干流地区进行了水质型缺水识别,构建了水量水质联合模拟与配置模型 [9] 。但现有水资源配置研究中时间尺度相对单一,较少考虑到时间尺度的变化给配置方案带来的风险,缺少不同时间尺度条件下配置方案的对比分析。Vörösmarty等与Mekonnen等分别基于年尺度与月尺度条件下的来水与需水,预测了全球范围的缺水状况,发现月时间尺度条件下的缺水量超出年时间尺度条件下缺水量的40%,全球范围内面临缺水威胁的人口数目增长了将近一倍 [10] [11],不同时间尺度条件下水资源短缺效应有所不同。而由美国康奈尔大学开发研制的Interactive River-Aquifer Simulation (IRAS) [12] 水资源配置模型具有灵活调整计算时间步长的功能,可适应多种时间尺度的输入(如日、旬、月、季、年等),能满足不同时间尺度条件下水资源配置效应的研究。
本文以汉江流域中下游地区为例,选择年、季、月三种时间尺度,在分析水资源系统中供、用、耗、排等主要环节所涉及的要素和相互连接关系的基础上,采用IRAS模型研究了水资源配置结果的时间尺度效应,以期为保障供水安全,提高配置效率提供支持。
2. IRAS模型结构与原理
IRAS是一款能够灵活调整计算时间步长,适用于多种时间尺度条件的水资源配置模拟模型,目前该模型已经在英国泰晤士河流域、墨西哥瓜纳华托河流域与葡萄牙瓜迪亚纳河流域等地区得到了成功应用。模型以物理水网为基础,通过分析水资源系统中供、用、耗、排等主要环节中所涉及的要素和相互连接关系,结合水资源系统的来水序列、需水序列、各类参数以及约束条件,在不同的调度运行规则下,模拟系统内水资源的配置过程。其原理如图1所示:
![](//html.hanspub.org/file/2-2410780x10_hanspub.png)
Figure 1. Framework of the IRAS model
图1. IRAS模型原理图
2.1. 系统概化
水资源系统规模庞大、涉及因素众多,结构纷繁复杂,而节点网络图是系统概化的具体表现,因此提炼出具有代表意义的节点网络图是建立系统框架与通用性模型的前提。
IRAS模型水资源系统中所涉及的实体元素可以分为节点与有向线段两种类型。节点用于描述研究区域中的工程设施、控制结构与用水单元,根据实体元素的属性可将其分为水源节点、来水节点、用水节点、分流节点、汇流节点。有向线段用于描述节点之间存在的水量传输与影响关系,根据水流是否能够双向流动可将其分为双向线段与单向线段。其中单向线段又可进一步分为需水线段、分流线段、天然线段,参与到地表水的输送与分配过程之中;双向线段通常结合达西定律用于地下水的模拟。
2.2. 参数设置
水资源配置的最终目的通常是为宏观决策提供理论参考,而宏观结果的可靠性离不开微观过程的精准性与合理性。模型中各类参数均能以“日”为单位实时变化,因此可结合工程实际状况与人工经验对各类参数进行灵活地调整与修正,以提高模拟过程的精准程度,为配置方案的可靠性提供保障。
2.3. 运行规则
1) 时间步长划分
IRAS模型以年为周期进行模拟计算,其时间步长可灵活调整,以适应不同时间尺度数据的输入,为在日、旬、月、季、年等时间尺度条件下水资源配置模拟计算提供了支持。每个时间步长可进一步划分为不同子时间步长(最小为1天),子时间步长的数目越多,模拟过程越精细,配置结果精度越高,所需计算时间越长。
2) 缺水计算
当用水户的用水需求高于来水条件时,采用公式(1)、(2)计算用水户的缺水量:
(1)
(2)
式中,t表示当前时间步长;tst表示时间步长包含子时间步长的总数;st表示当前子时间步长;
表示节点外推天然来水量;
表示节点总入流量;
表示节点来源于水库的入流量;
表示节点总需水量;f表示需求抑制系数;
表示节点缺水量;水量单位均为万m3。
按照水资源配置公平性、有效性与可持续性原则,根据水库与计算单元之间的水力联系,确定水库对用水户的供水比例,分解用水户缺水量,并加以累积,得到水库的拟放水量:
(3)
式中,
表示第k水库的拟放水量;
表示第i单元的缺水量;
表示第k水库与第i计算分区之间的水力联系;水量单位均为万m3。
3) 供水计算
水库群的供水流量根据水库运行调度规则确定。其中龙头水库的供水流量根据水库的蓄水线与供水流量线确定,其供水流量计算图如图2所示。根据模拟子时间步长在运行调度规则表中所处时段的初始最大、最小蓄水量,结束最大、最小蓄水量与初始最大、最小供水流量,结束最大、最小供水流量,并结合当前时刻水库群的总蓄水量,插值计算得该子时间步长内龙头水库的供水流量:
![](//html.hanspub.org/file/2-2410780x22_hanspub.png)
Figure 2. Calculation diagram of water supply discharge of the main reservoir
图2. 龙头水库供水流量计算图
(4)
(5)
(6)
(7)
(8)
(9)
(10)
式中,t表示模拟当前时刻;T表示模拟时段的总长度;Pt表示已模拟时段占模拟时段总长度比例;
、
表示水库群在t时刻的最大与最小蓄水量;
、
、
、
表示水库群在时段起始与结束时的最大、最小蓄水量;
、
表示龙头水库在t时刻的最大与最小供水流量;
、
、
、
表示龙头水库在时段起始与结束时的最大、最小供水流量;
表示水库群在t时刻的蓄水量;
表示水库群蓄水量的插值比例;
表示龙头水库在t时刻的供水流量;水量单位均为万m3;流量单位均为m3/s。
水库群中的其他水库的供水流量根据水量平衡关系表确定。当水库的蓄水量大于水量平衡关系表中的临界蓄水量时,超过临界值部分的蓄水量将排向下游。
4) 配水与耗水计算
水资源在各个用水户之间的分配按照优先权从高到低的顺序进行。各类用水户用水过程乘以耗水系数,即可得到耗水过程。
2.4. 约束条件
1) 计算分区水量平衡约束
(11)
式中:
表示第i计算分区的水量;
表示第m计算分区的水量;
表示第i与第m计算分区之间的水力联系;
表示第k水库的下泄量;
表示第k水库与第i计算分区之间的水力联系;
表示第i计算分区的天然来水量;
表示第i计算分区第j用水户的配水量;
表示第i计算分区第j用水户的耗水系数;
表示第i计算分区的损失水量;
表示第i计算分区的外调水量;水量单位均为万m3。
2) 水库约束
① 水库水量平衡约束
(12)
式中:
与
分别表示第k水库在第st + 1与第st子时间步长内的库容;
与
分别表示第k水库的来水量与下泄量;
表示第k水库的损失水量;水量单位均为万m3。
② 水库库容约束
(13)
式中:
与
与分别表示第n水库在第st子时间步长内最小与最大库容;水量单位均为万m3。
3) 非负约束
(14)
4) 需水约束
(15)
式中:
表示第i计算分区第j用水户在第st子时间步长内的需水量,万m3。
5) 供水能力约束
(16)
式中:
表示第i计算分区在第st子时间步内的可供水量,万m3。
3. 汉江流域中下游地区水资源配置模型构建
3.1. 研究区——汉江流域中下游地区概况
根据研究需要及掌握资料,本文以汉江中下游湖北省境内范围为研究对象,区域内国土面积约为4.11万km2,涉及神农架、十堰、襄阳、荆门、潜江、天门、仙桃、孝感和武汉等重要城市及其所辖的市、县、区以及“五三”、“沙洋”、“沉湖”等农场的全部或部分范围。汉江流域中下游地区包含以丹江口水库至汉口龙王庙的汉江干流及其分支东荆河为主要水源以及补充水源的供水范围,是实施最严格水资源管理制度的试点流域之一,且位于其上游侧的丹江口水库是南水北调中线工程的战略水源地,在全国水资源配置中具有重要地位。研究区域位于亚热带季风气候区,水资源总量丰富,多年平均降水量达800~1100 mm,但时空分布不均,呈现出南多北少,水量集中在汛期的格局。并且随着人口数目增长与社会经济发展,需水量不断增大,供水保证率越来越高,区域内水资源的供需矛盾日益凸显,水资源条件已经成为限制流域内经济发展与生态稳定的瓶颈。
3.2. 计算单元划分
为便于与行政区的来水、需水信息结合,以水资源四级区为基础,结合县(市)级行政区划分计算单元并将研究区域划分为28个计算单元,计算分区空间分布与系统节点网络概化分别如图3与图4所示。
3.3. 模型输入
1) 水平年及计算步长
以2010年为基准年,2030年为规划水平年,计算步长取为“年”、“季”、“月”。
![](//html.hanspub.org/file/2-2410780x70_hanspub.png)
Figure 3. Spatial distribution of computing units and reservoirs in the middle and lower reaches of Hanjiang River basin
图3. 汉江中下游各个计算单元与水库的空间分布
![](//html.hanspub.org/file/2-2410780x71_hanspub.png)
Figure 4. Schematic diagram of the middle and lower reaches of Hanjiang River basin
图4. 汉江流域中下游水资源系统节点网络概化图
2) 来水条件与需水条件
根据长江水利委员会编写的《长江流域水资源综合规划》成果,输入各个计算单元的水资源量及水库入库流量采用历史时期(1956~2011年)长序列水文资料逐年、季、月长系列径流过程和水库监测入库流量过程进行调算;需水基础数据主要来自长江水利委员会编写的《汉江流域综合规划》(下称《规划》)及《2010年湖北省水资源公报》(下称《公报》)相关成果。各计算单元的用水部门包含:生活用水、工业生产用水、农业生产用水3种河道外需水和河道内生态用水共4类,以《公报》相关数据对基准年的用水指标进行校核和修正,进而对规划水平年的用水指标及定额进行预测,表1展示了汉江流域中下游地区在2010基准年与2030水平年河道外各用水户的平均需水状况;河道内生态需水量的估算采用Tennant法 [13],其逐月需水数据通过年需水量乘以各个月份的分配系数得到,分配系数以《规划》报告提供数据为准。
![](Images/Table_Tmp.jpg)
Table 1. Annual off-stream water demand in the present and planning years (100 million m3)
表1. 现状与规划水平年河道外需水状况(亿m3)
3) 参数设置
根据优先保证生活和生态基本用水的原则,供水优先权从高到低依次为:河道内生态用水、生活用水、城镇生产用水与农村生产用水,以此分析汉江流域中下游地区的水资源供需情况;汉江流域中下游地区的河道内生态用水、生活用水、城镇生产用水与农村生产用水的回归系取值见表2;研究区域内共有17座大中型水库,水库的特征参数参见文献 [7] 。
![](Images/Table_Tmp.jpg)
Table 2. Regression coefficient of every type of water use in the Hanjiang River basin
表2. 汉江流域各类用水的回归系数
4. 结果分析
根据构建的水资源配置模型,对汉江流域中下游地区的用水、供水进行了模拟计算,分析了各计算单元2010基准年与2030水平年的三种时间尺度下的供需情况。汉江流域中下游地区在2010基准年与2030水平年中,年、季和月三种时间尺度下,各类用水户的总体缺水状况如表3所示:由表3可知:在相同时间尺度下,2030规划水平年的总体缺水率与各类用水户的缺水率均大于或等于(不缺水时) 2010基准年的总体缺水率与各类用水户的缺水率;在“年”、“季”、“月”三种时间尺度条件下,2010基准年与2030规划水平年中汉江流域中下游地区总体缺水率分别为0.7%、2.59%、5.39%与1.19%、3.27%、6.34%,均随着时间尺度的缩短而上升。
![](Images/Table_Tmp.jpg)
Table 3. Total water scarcity in the middle and lower reaches of Hanjiang River basin (10,000 m3)
表3. 汉江流域中下游地区总体缺水状况(万m3)
由于水资源管理常以地级市为单位,将各个计算单元的缺水状况汇总统计到地级市之中,以进一步分析研究区域缺水分布状况与不同时间尺度条件下各地级市的缺水状况。表4展示了2030规划水平年中汉江流域中下游地区各地级市在三种时间尺度下的缺水状况。天门、孝感、潜江、仙桃与武汉由于紧邻汉江下游干流,来水条件良好,水利工程调控能力强,供水充沛,在三种时间尺度条件下的缺水率均很小,水资源配置的时间尺度效应并不明显,故不列入表4之中。
![](Images/Table_Tmp.jpg)
Table 4. Water scarcity in every city in the middle and lower reaches of Hanjiang River basin (10,000 m3)
表4. 汉江流域中下游地区各地级市缺水状况(万m3)
十堰、神农架、襄阳与荆门在三种时间尺度条件下的缺水率均随着时间尺度的缩短而增加,水资源配置的时间尺度效应比较明显。以荆门市为例:荆门市需水量大,超过了天然来水条件,且水库调控能力有限,缺乏有效调节,缺水相对严重。图5展示了荆门市在三种时间尺度下的来水—需水关系,可以发现,长时间尺度配置条件平滑了来水与需水过程,需水与来水被均匀到年内其他时段,用水缺口有一定的减小:荆门市在月、季和年尺度下的用水缺口分别为136,674、132,527和132,527万m3,但是由于其年内各个时段的需水均处于较高水平,三种时间尺度下的缺水量的减少并不明显。而在长时间尺度条件平滑来水与需水过程同时,用水缺口在年内的分布也变得更加均匀,其变差系数分别为0.55、0.45和0,从而减小了水库在集中供水时期的供水压力,提升了区域的调控能力,因此长时间尺度条件下的配置结果缺水率低于短时间尺度条件下的配置结果缺水率。
5. 结论
本文以2010年为基准年,2030年为规划水平年,在分析水资源系统中供、用、耗、排等主要环节中所涉及的要素和相互连接关系的基础上,运用IRAS模型在汉江流域中下游地区建立了水资源配置模拟模型,并利用模型可灵活调整计算时间步长,时间尺度适应性好的特点,研究了汉江流域中下游流域在年、季、月时间尺度条件下的水资源配置状况,得到的主要研究结论如下:水资源配置结果的时间尺度效应在不同地区有所差异,
![](//html.hanspub.org/file/2-2410780x72_hanspub.png)
Figure 5. The process of water inflow and water demand in Jingmen city under seasonal and monthly scales
图5. 荆门市在季、月尺度下来水–需水过程
这是由于不同地区由于水资源条件和系统调节能力不同;具有良好来水条件区域,在不同时间尺度条件下均不发生缺水,配置结果明显无变化,水资源配置时间尺度效应不明显;在水资源短缺且系统的调节能力较弱地区,长时间尺度配置条件将平滑需水变量与来水变量在时间上的不均匀性,进而提高需水变量与来水变量的匹配程度,减小总体用水缺口,均匀缺水在年内的分布,使得长时间尺度条件下的配置结果缺水率低于短时间尺度条件下的缺水率,水资源配置时间尺度效应比较敏感。因此水资源配置中各个计算单元的时间尺度效应与区域水资源条件与系统调节能力相关,选择合适的时间尺度不仅是反映出水资源配置空间区域性特点,也是新时期水资源的精细化配置和风险管理的必然要求。
基金项目
国家自然科学基金(51879194, 91647106, 51579183)资助项目。