水资源作为支撑区域经济社会发展的基础性自然资源与战略性经济资源,长期以来在保障粮食安全、能源开发、生态保护与城乡发展等方面发挥着关键作用
[1-2]。随着全球气候变化加剧和城镇化进程的快速推进,人口增长和社会发展进一步驱动水资源需求持续增长,致使水资源供需矛盾日益凸显
[3]。因此,如何科学规划与优化配置水资源,对于破解供需瓶颈、支撑区域可持续发展具有重要战略意义
[4]。而科学合理预测区域水资源供需平衡态势是制定科学决策的重要基础和依据
[5],当前主要的研究方法有时间序列法
[6]、定额法
[7]、多元线性回归法
[8]、WEAP模型法
[9]、TOPSIS法
[10]和系统动力学法等
[11]。影响水资源供需平衡的因素诸多,通常涉及产业规模调整、农业灌溉面积、人口规模、城镇化进程以及相关政策等多个方面,常规的预测方法难以系统解析上述因素间的非线性反馈机制、时滞效应与动态累积效应,无法捕获水资源供需之间系统行为。相较而言,系统动力学法既能考虑经济发展和水资源之间的复杂反馈关系,又能完整反映系统中各因素动态过程
[12-13],已被广泛用于流域水资源承载力评估
[14]、节水政策模拟
[15]与风险预警
[16],是模拟区域水资源供需状况的有效方法。
赣江流域是鄱阳湖流域乃至长江中下游地区的重要子流域之一
[17],其水资源的合理利用和可持续管理对区域经济社会发展和水安全保障具有重要意义。当前关于赣江流域内水资源供需关系的研究多以行政区划为分析单元进行分析,缺乏基于整个流域视角的系统性供需平衡分析。为此,本文采用系统动力学方法(System Dynamics, SD),构建赣江流域水资源供需预测模型,对不同情景下的水资源供需平衡进行模拟仿真,进而探索符合赣江流域发展的最优路径,为赣江流域水资源管理与可持续利用提供决策参考。
1 数据来源与方法
1.1 研究区概况
赣江流域位于长江中下游南岸,呈南北走向贯穿江西全境,经鄱阳湖与长江连通。流域地貌以山地丘陵与低丘岗地为主,平原及水域覆盖面积不足4%,气候类型为亚热带季风气候,多年平均气温为17.2~19.3 ℃,多年平均降水量为1400~1 800 mm。流域水系发育良好,长度超30 km的干支流超百条,水资源总量较为丰富。2022年末赣江流域水资源总量为689.42亿m
3,占江西省水资源量的44.3%。流域总面积约83 500 km
2,约占江西省国土面积的50.03%,涵盖赣州、吉安、抚州、宜春、南昌、新余和萍乡等7个地级市在内的44个县(市、区),是江西省人口最为密集的区域之一
[18]。流域内以农业为主要产业,广泛种植水稻、油菜和柑橘等作物,工业基础较好,已形成电子信息、有色金属、装备制造、食品加工等为主导的产业体系。
1.2 数据来源
本研究所用到的数据可分为社会经济数据、农业数据、水资源数据和污水回用数据,具体如
表1所示。
1.3 研究方法
1.3.1 系统动力学法
系统动力学(System Dynamics, SD)是一种用于理解复杂系统行为的建模方法和工具,通过构建动态模型来模拟系统随时间变化的行为,帮助决策者理解系统的结构、行为以及不同组成部分之间的相互作用
[21]。针对赣江流域用水系统的非线性关系,系统动力学在处理具有多变量、多回路、非线性和高阶系统的数据时具有独特的优势
[22],能更好地解析用水结构的动态演化过程。本研究使用Vensim-PLE仿真平台建立用水系统动力学模型,以确立社会经济发展、用水需求、水资源供给之间的关系。建模过程如下:
(1) 模型边界设定。以赣江流域江西境内的行政区为空间边界,涉及南昌市、吉安市、赣州市、宜春市、萍乡市、新余市和抚州市7个地级市。以2007—2035年为模拟的时间边界,选取2007年为基准年,时间步长为1 a。其中,2007—2022年为历史拟合年份,2023—2035年为仿真模拟年份。
(2) 模型变量及主要方程。本研究构建的赣江流域水资源供需平衡系统动力学模型共包括76个变量,其中变量参数58个,状态变量9个,速率变量9个。模型的主要方程详见
表2。
(3) 模型子系统划分及系统动力学流程图绘制。充分考虑赣江流域生活、工业、农业用水需求和水资源供给,将系统划分为社会经济子系统、用水需求子系统、水资源供给子系统和水资源供需子系统。其中,社会经济子系统提供用水需求外部驱动,用水需求子系统聚合农业、工业、生活、生态等所有用水需求,水资源供给子系统描述水资源的获取与供水能力,水资源供需子系统综合分析供需关系。根据各因素之间相互影响及其作用,绘制赣江流域水资源供需分析流图(见
图1)。
(4) 模型率定和验证。为提高模型的模拟精度与预测能力,按照7∶3的比例,划分2007—2017年为率定期,用于调整参数,使模拟结果尽可能贴近历史数据;2018—2022年为验证期,用于检验模型拟合效果与趋势再现能力。模型参数分为初始值、常数值和表函数等3类,其中初始值(总人口、有效灌溉面积等)直接来源于统计年鉴、政府公报等文件;常数值(工业废水排放系数、生活污水排放系数等)依据相关计算手册确定;表函数(地表水供水量、城市道路增长率等)通过历史数据,采用趋势外推、ARIMA模型
[23]、复利法等方法得到。
为评估模型的拟合精度与预测性能,本文采用相对误差(Relative Error, RE)对模型模拟结果与历史数据进行验证
[24]。公式如下:
式中:xi 为i年模型模拟值;yi 为i年实测值。
选取生活需水量、农业需水量、工业需水量和需水总量等关键变量进行误差分析,当模型模拟结果与实测数据的相对误差控制在±10%以内,则认为模型具备较好的拟合精度与可信度,适用于后续的未来情景模拟分析
[25-26]。经过参数优化后,如
表3所示,选取的4个关键变量模型模拟值在率定期内的误差绝对值均低于8%,而验证期误差绝对值均低于6%,满足精度要求,可进一步开展未来情景模拟。
1.3.2 情景设置
为评估赣江流域水资源供需关系在不同发展路径下的动态演变,本文在梳理历史时期国民经济与社会发展特征的基础上,结合相关规划,充分考虑未来影响水资源可持续利用的关键因素,并参考相关文献设置
[27-29],构建了4种发展情景,包括基准情景(S
1)、经济社会高增长情景(S
2)、节水优化情景(S
3)和协调发展情景(S
4),以模拟不同发展路径下的用水演变趋势及水资源供需状况。
(1) 现状延续情景(S1):以2007—2022年实测数据为基础,按照《江西省国民经济和社会发展第十四个五年规划和2035年远景目标纲要》及《江西省“十四五”节水型社会建设规划》要求,2020—2025年江西省城镇化率年均增长1.1%,万元工业增加值用水量年均下降4.3%,城镇居民用水定额年均下降1%。其他用水定额依据《生活及服务业用水定额第2部分:服务业、居民生活和建筑业》设置,同时假设未来无重大水资源管理政策出台,社会经济发展、用水习惯等因素延续历史发展趋势。
(2) 经济社会高增长情景(S2):在S1延续的基础上,以经济快速发展为主要导向,设定城镇化进程显著加快并带动人口集聚,同时提升工业发展水平,其中工业产业增加值增长率为4%,人口增长率、城镇化率分别每五年提高1.25‰、10%。
(3) 节水优化情景(S3):在S1延续的基础上,强化节水约束,社会经济发展适度放缓,万元工业产业增加值年增长率为1.8%,城镇化率每五年累计提升4%,万元工业产业增加值用水量每五年下降30%,其他定额每五年累计降低10%,非常规水资源利用水平显著提升,生活污水收集率及再生水利用率累计每五年分别提升6%和8%。
(4) 协调发展情景(S
4):在S
1,S
2和S
3基础上,融合“绿色发展”“高质量增长”理念,协同推进“多源水利用”和“生态优先”策略,设定万元工业产业增加值年增长率为3%,万元工业产业增加值用水量每五年下降25%,其他用水定额每五年累计降低6%,城镇化率、污水收集率及再生水利用率分别每五年累计提升6%,5%和5%。具体参数设定见
表4。
1.3.3 供需平衡指数
引入文献
[7]的供需平衡指数(RWRI)评价方法,如
表5所示,定性对水资源供需状况评价。
式中:RWRI为供需平衡指数,WD为区域水资源需水量,WS为区域水资源供水量。
根据
表5所示的评价标准,依据RWRI值所处区间,将当年的水资源供需状况划分为不同类型。
2 结果与分析
2.1 需水变化分析
基于4种情景设定,本文系统模拟赣江流域2023—2035年水资源供需的动态演变过程,重点分析生活、工业、农业以及总需水等不同用水类型的变化趋势结构演变特征。
2.1.1 需水量变化趋势
(1) 生活需水量(
图2A):4种情景的生活需水量呈现“两升两降”的变化趋势。其中,经济社会高增长情景(S
2)变化最为显著,2023年的11.74亿m
3增长至2035年的13.62亿m
3,增幅约16%,主要受到区域经济快速发展引发的人口规模扩大和城镇化进程加速的共同驱动;现状延续情景(S
1)增长相对平稳,由11.72亿m
3增至12.97亿m
3,增幅约11%;而节水优化情景(S
3)与协调发展情景(S
4)需水量呈现下降趋势,分别由11.59亿m
3降至10.87亿m
3,11.65亿m
3降至11.54亿m
3,下降幅度分别为6%和1%,表明节水与协调发展政策和人们节水意识提高有效抑制了生活用水需求的快速增长。
(2) 工业需水量(
图2B):4种情景下工业需水量均出现先小幅上升(2023—2025年),差距不大,主要源于初期技术改造及产能调整降低了单位产值耗水强度,随后随产业扩张速度和节水力度的不同快速分化(2026—2035)。其中经济社会高增长情景(S
2)总体增长最快,由2023年的25.49亿m
3增至2035年的42.44亿m
3,增幅约67%,主要为该情境下工业增值增长最快;现状延续情景(S
1)维持中等增速,2035年达到34.27亿m
3,较2023年增长约34%;协调发展情景(S
4)兼顾经济增长与资源节约,2035年增长至31.20亿m
3,较2023年增长约24%;相反,节水优化情景(S
3)总体呈现下降趋势,2035年达到14.66亿m
3,较2023年下降约6%,表明高效用水技术推广和产业结构优化对工业用水增长形成了显著抑制。
(3) 农业需水量(
图2C):与生活需水相似,农业需水在4种情景下亦呈现“两升两降”的分化趋势。其中现状延续情景(S
1)与经济高速发展情景(S
2)增长一致且需水量持续上升,从90.38亿m
3增长至96.72亿m
3,增幅约7%,主要是灌溉面积扩张和畜牧规模扩大,且缺乏强力节水手段所致;相反,节水优化情景(S
3)和协调发展情景(S
4)呈下降趋势,S
3借助高效灌溉与压减牲畜用水等措施将农业用水从86.94亿m
3压缩至59.55亿m
3,降幅约为32%;S
4在兼顾生态和粮食安全的前提下,农业用水由88.96亿m
3降至79.70亿m
3,降幅约10%。
(4) 总需水量(
图2D):4种情景下总用水量呈现不同的变化趋势。其中,现状延续情景(S
1)和经济高速发展情景(S
2)呈现增长趋势,情景S
1,情景S
2分别由128.31亿m
3升至144.76亿m
3(年均136.91亿m
3),128.32亿m
3升至153.57亿m
3(年均140.62亿m
3),增幅分别约为13%和20%,生活与工业需水的快速增长是两情景下总需水量上升的主要驱动因素;节水优化情景(S
3)和协调发展情景(S
4)呈现下降趋势,情景S
3有123.76亿m
3快速下降至85.88亿m
3(年均105.23亿m
3),降幅约为31%,农业大幅节水与工业用水增速受控共同促成了总需水量变化的主要原因;相较于情景S
3,情景S
4增长较为平缓,仅由126.41亿m
3升至123.23亿m
3(年均126.01亿m
3),降幅为3%,表明在适度经济增长与节水并举条件下,可有效控制需水总量。
2.1.2 需水结构变化分析
图3展现了2023—2035年4种情景下赣江流域生活、工业与农业需水在总需水量中占比结构演变趋势。现状延续情景(S
1)、经济高速发展情景(S
2)和协调发展情景(S
4)中均呈现工业占比上升、农业占比回落、生活和生态环境占比基本稳定的特征,其中,缺乏额外节水驱动的情景S
1中,农业需水比重由70%跌至67%,工业需水比重从20%稳步抬升至24%;工业扩张效应最强的情景S
2,工业需水占比由20%大幅攀升至28%,农业占比由70%跌至63%;情景S
4借助高效灌溉、产业升级等手段,工业需水占比由20%升至25%,农业需水占比70%降至65%,用水结构趋向均衡。相较之下,节水优化情景(S
3)通过强节水措施,将工业需水占比由20%压至17%,农业占比由70%减少至69%,同时生活需水比例则由9%上升至13%,生态环境用水占比保持基本稳定。
2.2 水资源供需平衡分析
如
图4所示,模拟期内4种情景年缺水量呈现不同的变化趋势。其中,现状延续情景(S
1)呈现出先下降后持续抬升,增长趋势逐渐放缓,从2023年7.66亿m
3至2035年23.11亿m
3(年均15.74亿m
3),增长了近3倍;经济高速发展情景(S
2)供需缺口加速扩大至2035年31.79亿m
3(年均19.41亿m
3),增长量是S
1的1.56倍,表现出高强度工业化和城镇化加大了水资源供需不平衡;相比之下,协调发展情景(S
4)曲线变化最为平缓,从2023年5.74亿m
3至2035年2.39亿m
3(年均1.50亿m
3),下降了0.26倍左右,是适度经济增长和节水力度共同作用的良性体现;最为显著的是节水优化情景(S
3),供需缺口持续下降,缺口于2024年后转负并在2035年盈余35.80亿m
3(年均-15.92亿m
3),累计削减了38.92亿m
3,是强节水力度的体现。
为进一步探讨4种情景下水资源供需状况,计算模拟期水资源供需平衡指数,结果如
图5所示。总体来看,4种情景呈现“二升二降”的分化格局:现状延续情景(S
1)、经济高速发展情景(S
2)的RWRI值总体上均保持上涨,协调发展情景(S
4)、节水优化情景(S
3)逐步下沉,且模拟期4种情景RWRI值均在1.4(轻度短缺范围)以下。经济高速发展情景(S
2)抬升最为迅猛,2035年RWRI值为1.26(年均1.16),缺水风险越过水资源平衡范围到达轻度短缺范围;现状延续情景(S
1)以较缓坡度上扬,模拟末期逼近临界短缺上边界,RWRI值为1.19左右(年均为1.13),表明若维持现有节水水平,缺水压力仍将逐年积累,存在一定的缺水风险;协调发展情景(S
4)曲线总体平稳下降,RWRI值由1.05降至1.01(年均1.04),临近平衡有余区域,表明适度增长与节水并举可将缺口控制在可调度范围;节水优化情景(S
3)RWRI值持续下滑,由2024年进入平衡有余区域,2028年进入轻度盈余区域,2032年跌破0.8进入中度盈余区域,2035年跌至0.71(年均0.87),成为4个情景中唯一落入“中度盈余”区的情景,缺水风险极小,但经济增长较慢。由此看来,未来赣江流域在维持一定经济的增长的同时,需强化节水措施,严控缺水风险。
2.3 情景结果比较
在未来发展过程中,若赣江流域按照经济高速发展情景(S2)路径推进,高耗水产业与人口集聚将使流域迅速逼近轻度短缺,需要“开源节流”,提前布局非常规水源并实施限水策略;从现状延续情景(S1)可以看出即便不追求高速增长,若缺乏新增节水措施,同样难以摆脱水资源临界短缺;节水优化情景(S3)提供了风险最低的参考路径,以节水为主导,提高节水水平的同时抑制工业增速,并大力开发非常规水源,最大限度缓解流域水资源紧张。协调发展情景(S4)在保障合理经济增长速率与民生用水需求的同时,综合考虑经济发展、污水处理技术、人民对生活水平的需要等多方面因素,以适当的速度保证未来经济发展,又能缓解水资源利用压力,符合可持续发展内涵。
综上所述,赣江流域未来水资源系统面临较大不确定性,特别是在经济高速发展情景下存在明显的供需矛盾。而协调发展策略在保证合理的经济增长的同时,可显著缓解缺水压力,提升系统韧性与适应能力,可作为流域水资源管理与政策制定中优先考虑。
3 结 论
以赣江流域为例,构建了基于SD的赣江流域水资源供需平衡仿真模型,并设定4种情景,对2023—2035年不同发展路径下的赣江流域用水演变趋势及水资源供需状况进行模拟分析,主要得出以下结论:
(1) 模拟期内,4种情景下总需水呈现不同的变化趋势。其中,情景S1,S2呈上升趋势,情景S3,S4呈下降趋势。情景S2上涨最快,情景S1次之,情景S3下降最快,情景S4下降较为缓慢,年均总需水量分别为140.62亿m3,136.91亿m3,105.23亿m3,126.01亿m3。就需水结构而言,情景S1,S2和S4均表现为工业需水占比上涨,同时农业需水占比下降,生活需水和生态需水占比基本稳定;而情景S3表现为工业和农业需水占比下降,生活需水占比提高,生态环境需水占比基本稳定,表明在维持适度经济增长的同时强化节水措施,不仅能有效控制需水总量增长,还能优化用水结构。
(2) 模拟期内,4种情景年均供需缺口依次为19.41亿m3(情景S2)>15.74亿m3(情景S1)>1.50亿m3(情景S4)>-15.92亿m3(情景S3),对应年均RWRI值分别为1.16,1.13,1.04,0.87,4种情景均处于临界短缺评价区间内及以下,但分化趋势明显:情景S2已超过临界短进入轻度缺范围,情景S1,情景S4均在临界短缺和平衡有余范围内,情景S3在由临界短缺进入中度盈余区域。若不加强节水和非常规水利用,赣江流域缺水风险将在高增长情景下显著加剧。
(3) 情景S4兼顾经济发展与节水目标,既保持适度增长,又将缺水量控制在可调度区间,明显缓解水资源压力,是最值得采纳的方案。未来应完善流域统筹规划,推广高效节水技术,提升污水处理及其他非常规水利用水平,以实现赣江流域水资源的可持续利用。