气候变化对南水北调西线工程的影响及建议

关键词：气候变化；ＲＣＣＣ－ＷＢＭ 模型；水资源；南水北调西线工程

中图分类号：Ｐ４６７；ＴＶ６８；Ｐ３３３ 文献标志码：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２４．０９．００９

引用格式：王国庆，虞畅，金君良，等．气候变化对南水北调西线工程的影响及建议［Ｊ］．人民黄河，２０２４，４６（９）：６０－６６．

０引言

我国水资源空间分布不均衡，呈“南多北少” 格局［１］ ，长江流域年降水量超过１ １００ ｍｍ，而黄河流域大部分地区年降水量只有２００～６００ ｍｍ。黄淮海流域作为重要的工业和粮食生产基地，以全国７％的水资源量支撑着全国２４．２％的人口、四分之一的耕地和国内生产总值［２－３］ 。北方水资源短缺已成为严重影响工业、农业和经济社会发展的一个制约性因素［４］ 。２００２年，水利部批准了南水北调国家战略性工程，旨在缓解北方水资源短缺问题［５－６］ 。南水北调工程通过东、中、西３ 条调水线路将长江、黄河、淮河、海河四大流域相联系，构成了以“四横三纵”为主体的国家水资源配置总体布局，是实现水资源南北调配、东西互济，提高水资源与人口、经济时空匹配性的“国之重器”［７－９］ ，计划每年通过３ 条线路调水１４０ 亿ｍ^３［１０］ 。截至２０２１ 年年底，南水北调东线、中线一期工程已为华北地区输送了４９４ 亿ｍ^３的水量，相当于长江流域年径流量的５％，黄河流域年径流量的８０％。西线工程作为我国水资源配置“四横三纵”总体格局中重要“一纵”，规划将从长江上游金沙江、雅砻江、大渡河等河流水系引水入黄河上游，以缓解我国北方地区干旱缺水问题，目前尚处于前期论证阶段［１１］ 。

以全球气候变暖为特征的气候变化将改变区域水循环，进而导致降水、蒸发、径流等要素在时空分布上发生复杂变化［１２－１３］ 。联合国政府间气候变化专门委员会（ Ｉｎｔｅｒｇｏｖｅｒｎｍｅｎｔａｌ Ｐａｎｅｌ ｏｎ Ｃｌｉｍａｔｅ Ｃｈａｎｇｅ，ＩＰＣＣ）第六次评估报告指出，２０１１—２０２０ 年全球平均气温较工业革命前升高１．１ ℃，未来全球气温将继续升高，预计大多数陆地区域的降水和地表水资源季节和年际变化将更加剧烈［１４］ 。利用合适的水文模型耦合全球气候模式（Ｇｌｏｂａｌ Ｃｌｉｍａｔｅ Ｍｏｄｅｌｓ，ＧＣＭｓ） 并模拟不同气候情景下流域的水文响应，是评估气候变化对水资源影响的主要手段之一［１５－１７］ 。Ｂａｏ 等［１８］ 将ＶＩＣ 模型与ＣＭＩＰ５ 的气候模式输出相结合，预测了海河流域未来流量的变化趋势，得出流量有增大的趋势，北部流域的增大趋势较南部流域的明显。Ｚｈａｎｇ等［１９］ 利用ＳＷＡＴ 模型分析了不同排放情景下黄河上游径流量的变化，发现中等排放情景和低排放情景下唐乃亥水文站流量具有增大趋势，而高排放情景下则呈减小趋势。Ｗａｎｇ 等［２０］ 利用８ 种全球气候模式情景和ＶＩＣ 模型评估了未来气候变化对长江流域的潜在影响，结果表明，气温升高将导致年平均径流量小幅减小。Ｓｕ 等［２１］ 评估了气候变化对长江上游流域径流量的影响，发现到２１ 世纪末，年径流量可能有所增加。尽管多数研究结果表明长江上游未来降水可能增多，但考虑到气温升高和寒区冰雪的变化，一些研究认为未来汛期流量可能会减小，径流趋势预估的不确定性仍然不可忽视［２２－２４］ 。

气候变化将使南水北调工程所涉长江和黄淮海流域来水条件发生变化，直接改变调水区可调水量和受水区水资源需求，进而影响调水规模和工程运行方式。本研究基于１９６１—２０２０ 年西线工程调水区和受水区实测资料以及未来气候模式情景，通过数理统计与水文模拟，分析未来气温、降水变化对西线工程调水区和受水区水资源的影响，以期为南水北调西线工程的规划、建设和运行管理提供借鉴。

１资料与方法

１．１研究区域概况与数据来源

根据《南水北调工程总体规划》，南水北调西线工程在长江上游通天河、支流雅砻江和大渡河上游筑坝建库，开凿穿过长江与黄河分水岭巴颜喀拉山的输水隧洞，调长江水入黄河上游。规划调水规模为１７０ 亿ｍ３，主要目标是解决涉及青海、甘肃、宁夏、内蒙古、陕西、山西等六省（自治区）黄河上中游地区和渭河关中平原缺水问题。根据黄河勘测规划设计研究院有限公司完成的《南水北调西线工程规划方案比选论证》（２０２０年１０ 月），规划的西线工程调水区为金沙江巴塘站以上区域、雅砻江雅江站以上区域和大渡河大金站以上区域，将黄河中上游流域作为工程的受水区。图１ 为南水北调西线工程调水区和受水区范围。

本研究采用的数据资料包括南水北调西线工程调水区和受水区的实测气温、降水量、径流量，以及全球气候模式预估的未来气温、降水量。其中：０．２５°格点历史再分析气象资料来源于中国气象局国家气象信息中心，实测水文资料来源于长江水利委员会水文局和黄河水利委员会水文局，资料序列为１９６１—２０２０ 年；花园口站天然径流量资料来源于第三次全国水资源评价和黄河流域水资源公报。未来气候情景采用１７ 个ＣＭＩＰ６ 全球气候模式在ＳＳＰ２－４．５ 排放情景下降尺度后的气温和降水资料，来源于国家气候中心，资料序列为２０２５—２０９９年，空间分辨率为０．２５°。

根据《南水北调工程总体规划》，将１９６１—２０００年作为基准期，２００１—２０２０ 年作为重点比较期；未来中期关注２０３５年水平，以２０３１—２０４０年均值代表；未来远期关注２０５０年水平，以２０４６—２０５５年均值代表。

１．２研究方法

采用Ｍａｎｎ－Ｋｅｎｄａｌｌ 检验方法诊断气温、降水和径流的演变趋势。该方法是一种分析数据序列随时间变化趋势的非参数统计检验方法［２５］ 。其优点是不需数据服从特定的分布， 同时检验范围较宽， 被广泛应用于水文气象要素的趋势分析和显著性检验［２６］ 。

采用ＲＣＣＣ－ＷＢＭ 模型，模拟南水北调西线工程调水区和受水区的天然径流量过程。该模型是水利部应对气候变化研究中心研发的大尺度流域水文模型［２７－ ３０］ ，考虑了地面径流、地下径流和融雪径流３ 种径流成分。针对大尺度流域，为解决下垫面和气候空间分布不均匀问题，该模型将流域划分为若干个正交网格，对每个格点分别进行产流计算，汇流计算则借鉴ＶＩＣ 模型汇流计算方案；模型输入包括每个格点逐时段的气温、降水量和蒸发能力。

２结果与讨论

２．１西线工程调水区和受水区气候与径流演变特征

西线工程调水区和受水区多年平均气温分别为－２．４２ ℃和５．７９ ℃，多年平均年降水量分别为５６２、４７２ ｍｍ。图２ 为西线工程调水区和受水区１９６１—２０２０年平均气温和年降水量的演变过程。

由图２ 可以看出：

１）西线工程调水区和受水区年均气温均呈现显著性升高趋势，线性上升速率分别为０．３３４、０．３１７ ℃ ／１０ ａ，均高于全球和全国平均升温速率；相比而言，调水区年均气温较低，但其升温速率高于受水区。

２）西线工程调水区年降水量为非显著性增加趋势，线性增长速率为１２．２ ｍｍ／１０ ａ；受水区年降水量呈现先减少后增加趋势，１９６１—２０２０年总体为非显著性增加趋势，线性增长速率为４．８ ｍｍ／１０ａ。

调水区由金沙江、雅砻江和大渡河３ 个流域的子区域组成，多年平均年降水量分别为４８５．３、８４４．４、９５０．２ ｍｍ；经９５％置信水平的统计检验，金沙江流域子区和雅砻江流域子区年降水量均为非显著性增加趋势，大渡河流域子区年降水量为非显著性减少趋势。

２００１—２０２０年调水区和受水区平均气温和降水量均较基准期（１９６１—２０００ 年）有所升高和增加（见表１），两个区域气温分别升高１．１９、１．０６ ℃，年降水量分别增加３４．６、２７．２ ｍｍ（占基准期的６．２８％和５．８７％）。

长江上游位于青藏高原东部，海拔高，生态环境脆弱，河川径流量受人类活动影响相对较小，调水区的巴塘站、雅江站和大金站实测径流量接近天然状态。受水区黄河中上游流域气候干旱，水资源短缺。调水区和受水区年径流量演变过程见图３。

由图３ 可以看出：调水区巴塘站和雅江站实测径流量均呈现先减少后增加总体略增的态势，线性增长速率分别为０．７６０ 亿、０．５３３ 亿ｍ^３／ａ，而大金站实测径流量呈现非显著性减少趋势，线性增长速率为－０．０４５亿ｍ^３／ａ；受水区（花园口站）天然径流量和实测径流量的线性增长速率分别为－２．９８２ 亿、－４．４５８ 亿ｍ^３／ａ，经９５％置信水平统计检验，均呈现显著性减少趋势。

统计结果表明：

１）受水区花园口站多年平均天然径流量和实测径流量分别为５３８ 亿、３５２ 亿ｍ^３；而调水区巴塘站、雅江站和大金站１９６１—２０２０ 年平均年径流量分别为２９１ 亿、２１７ 亿、１６５ 亿ｍ^３；调水区３ 个站的实测年径流量接近受水区花园口站实测年径流量的２ 倍。

２） ２００１—２０２０ 年调水区年径流量较基准期（１９６１—２０００ 年）增加９．０％，其中：巴塘站和雅江站年径流量分别增加１６．８％和９．３％，大金站年径流量减少４．３％；受水区花园口站实测年径流量和天然年径流量较基准期分别减少２６．９％和１２．９％。

２．２ 西线工程调水区和受水区逐月径流量模拟

将调水区和受水区划分为分辨率为０．２５°的网格，调水区和受水区分别包括３５８和８１４ 个格点。采用基准期１９６１—２０００ 年的格点气象资料驱动基于０．２５°格点的ＲＣＣＣ－ＷＢＭ 模型，以Ｎａｓｈ－Ｓｕｔｃｌｉｆｆｅ模型效率系数和相对误差为目标函数，以１９６１—１９９０年为率定期率定模型参数，以１９９１—２０００年为验证期检验模型的稳定性和模拟效果，分别模拟巴塘站、雅江站、大金站和花园口站月径流量过程。４ 个站月径流量的模拟效果统计见表２。为直观起见，将调水区雅江站和受水区花园口站１９６１—１９６９ 年实测与模拟月径流深变化过程绘于图４。

由表２ 可以看出，ＲＣＣＣ－ＷＢＭ 模型对调水区和受水区天然径流量具有较好的模拟效果，率定期和检验期Ｎａｓｈ－Ｓｕｔｃｌｉｆｆｅ 模型效率系数均超过０．７５，平均相对误差在±３％以内。由图４ 可以看出，两个代表性水文站实测与模拟月径流量过程总体拟合较好，只有个别峰值径流量模拟误差较大。总体而言，ＲＣＣＣ－ＷＢＭ模型能够较好地模拟南水北调西线工程调水区和受水区天然径流量过程，可以用来模拟未来气候变化下两个区域的水资源情势。

２．３西线工程调水区和受水区未来气候变化趋势

鉴于不同模式预估的气温和降水变化存在差异，参考ＩＰＣＣ处理不确定性的方式，采用四分位分析方法描述未来气温升降、降水增减的可能性及变化幅度。图５ 和图６ 分别给出了南水北调西线调水区和受水区未来中期、长期气温、降水较基准期的变化。

由图５ 和图６ 可以看出：

１）１７ 个气候模式均预估调水区和受水区未来中期和长期气温较基准期有所升高，多数模式预估未来降水量有所增加。

２）调水区在未来中期和长期气温较基准期分别升高１．９２ ℃［１．４３ ℃， ２．９４ ℃］和２．３６ ℃［１．６４ ℃， ３．６１ ℃］；降水变化３．３７％ ［－ ０． １５％， １３． ４％］ 和３． ３９％ ［－０．８４％，１７．５０％］；调水区未来气候以暖湿化为主。

３）受水区在未来中期和长期气温较基准期分别升高１．７７ ℃［１．０８ ℃， ２．５９ ℃］和２．３４ ℃ ［１．１７ ℃，３．３４ ℃］；降水变化３．２４％［－２．２４％， ７．３２％］和４．３４％［－４．４８％， ９．１３％］；受水区未来气候变化趋势与调水区总体一致，以暖湿化为主。

２．４未来气候变化下西线工程调水区和受水区径流演变趋势

基于不同模式情景下模拟的未来年径流量变化过程，分析南水北调西线工程调水区和受水区未来不同阶段年径流量较基准期的变化，见图７。

从图７可以看出：调水区未来年径流量以略增为主，未来中期和长期年径流量较基准期分别增加１．４２％ ［－６．６９％， ２２．９５％］ 和２．０８％ ［－５．８４％，３１．７７％］；受水区未来年径流量以正常略减少为主，未来中期和长期年径流量较基准期变幅分别为－１．０２％［－１４．７９％，４．３４％］和－０．２８％［－２３．３１％， ６．２１％］。

２．５讨论与建议

２．５．１西线工程水资源条件的影响分析

１）根据降水、径流量的历史变化规律和未来中长期变化态势，西线调水区水源条件可能会有所改善。２００１—２０２０年西线调水区年降水量和径流量总体以增加为主，但存在明显空间差异性，其中金沙江巴塘水文qWrISiRsNow/8C5n5NoZP8PNJnU07QaAnYEQngFgIAQ=站、雅砻江雅江水文站天然年径流量分别增加１６．３％和９．８％，大渡河大金水文站天然年径流量减少４．３％；预估未来中期和长期天然年径流量较基准期也有所增加，但增幅仅２％左右。

２）西线工程受水区天然年径流量不足的局面在２０００年之后有所加剧，且在中长期可能无法扭转，将增加西线受水区跨流域调水需求。西线受水区为黄河流域，２００１—２０２０ 年实测年径流量和天然年径流量均较基准期明显衰减，分别减少２６．９％和１２．９％；预估未来中期和长期天然年径流量较基准期偏少１％左右。因此，受水区天然径流量不足局面可能难以改变。

２．５．２主要认识与建议

１）根据水资源变化空间差异性，优化南水北调西线调水方案。西线调水区２００１—２０２０ 年天然径流量较基准期虽总体有所增加，但存在明显的空间差异性，因此应充分考虑和主动适应气候变化背景下西线调水区有关江河水资源条件演变差异性，优化调水方案，确保西线工程科学规划建设和长效运行。

２）深化气候变化对流域水资源条件不确定性影响的研究。尽管全球气候模式是目前预估气候变化最有力的工具，但未来气候变化情景依然存在不确定性，特别是对降水预估结果不确定性更大；另外，随着气温继续升高，极端气候和水事件增多趋强，降雨时空分布发生变化，更不利于水资源的管理。

３结论

南水北调西线工程是国家水安全保障体系的重要组成部分，对于提升国家水资源调配能力、保障国家水安全具有重要意义。气候变化对南水北调西线工程调水区和受水区水资源具有一定影响，研究得出主要结论如下：

１）１９６１—２０２０年，西线工程调水区和受水区气温呈显著升高趋势，年降水量呈非显著性增加趋势，调水区和受水区气候以暖湿化趋势为主；调水区年径流量呈现增加趋势，但空间上存在一定差异，受水区年径流量显著减少。

２）未来中期和长期，西线工程调水区和受水区气温继续显著升高，年降水量多呈非显著性增加趋势；受气候变化影响，西线工程调水区年径流量较基准期略增，受水区年径流量与基准期基本持平略有减少。

３）总体而言，未来气候变化下西线调水区年径流量相对稳定，能够满足调水的需要，受水区径流量不足局面可能难以改变，跨流域调水需求有所增加。气候变化及其影响预估结果具有一定的不确定性，未来仍需加强气候变化影响的不确定性研究。