【文献分享】多空间尺度下水资源供需平衡对驱动因素的敏感性研究

摘要

黄河流域日益严峻的水资源不安全问题，本质上源于气候驱动的产水量（WY）与农业、工业、城市用水需求（WC）之间不断加深的结构性失衡。本研究结合GIS制图、随机森林模型与Geodetector—OLS分析，构建了综合水足迹核算框架，从栅格、城市和子流域三个尺度系统评估了2000—2024年黄河流域水资源供需平衡及其驱动机制。研究发现，黄河流域水资源压力呈现显著的“上游相对宽松—中下游长期赤字”梯度格局。进一步分析表明，降水季节性和土地利用转型是WY的主导因素，灌溉面积扩张与工业用水强度则是用水需求增长的关键驱动因素；而严重赤字区往往源于耕地扩张与地下水依赖等因素的耦合作用。结果表明，流域水资源压力并非单一自然约束的结果，需要基于赤字类型实施差异化治理，以实现水资源安全与流域可持续发展的协同推进。

淡水资源是支撑生态系统稳定与人类社会运行的关键基础，但在人口增长、工业扩张、农业集约化与气候变化等多重因素作用下，越来越多流域正面临水资源供需失衡加剧的问题。黄河流域作为典型的水资源约束型流域，内部降水梯度显著、年际波动强烈，叠加灌溉农业、能源开发与快速城镇化等高强度人类活动，水资源供需矛盾尤为突出。现有研究仍难以充分揭示气候驱动的水资源供给与部门用水需求在嵌套空间单元中的耦合关系及其空间分异特征，尤其流域尺度的总体评估往往掩盖子流域和城市尺度上具有重要治理意义的局部赤字。同时，传统体积法难以区分蓝水、绿水和灰水的差异化形成机制与治理路径。相较之下，多组分水足迹分解框架能够更精准地识别水资源压力来源，为多尺度供需失衡诊断、部门需求归因与差异化治理提供了更具针对性的分析工具。

本研究为气候异质性显著的大型流域水资源治理与多尺度水文过程认知做出了重要贡献：（1）构建了一个全面的多尺度水足迹核算框架，综合刻画栅格、城市和子流域尺度上的水资源供需动态，能够同时识别黄河流域农业、工业与生态用水压力的空间分异特征。（2）基于2000—2024年长时间序列数据，揭示了在研究单一行政或水文单元时可能被忽视的关键治理盲点。（3）通过整合25年的环境和社会经济数据，本研究提出了具有较强可操作性和数据支撑的分层分类治理路径，为黄河流域不同尺度下的水资源优化配置提供了决策依据，也为其他面临水资源、生态与社会经济复合压力的流域治理提供了可借鉴的分析框架。

1.蓝水和绿水足迹的校准与验证

蓝水与绿水足迹模型基于2000—2012年黄河水资源公报的流量与分水数据进行校准，并在2013—2024年开展独立样本验证。全流域尺度校准的Kling-Gupta效率（KGE）为0.88，相对体积误差（RVE）为+2.3%，优于KGE≥0.75、RVE≤±8%的接受阈值。子流域尺度上，KGE值在0.76（下游灌区）至0.84（上游）之间，所有子流域RVE均保持在±7%以内。黄土高原河段均方根误差（RMSE）最大，为5.93毫米，主要受对流性降水导致的高年际径流变率影响。2013—2024年验证期的流域尺度KGE为0.85，相对体积误差为+3.2%，说明模型在未参与校准的时段中依然具有较好的预测能力。各空间尺度的校准与验证统计量汇总于表1。

模拟与观测的年均蓝水足迹在子流域尺度上整体贴近1:1线（见图1），在50–70毫米的动态区间内异方差性最小，表明模型在上游高供水与下游低供水条件下均表现一致。残差分布近似正态、均值为零，未出现系统性高估或低估，验证了水量平衡校准的有效性（见图2）。

2.分行业蓝水和灰水足迹动态变化特征

2022—2024年间，全流域年均蓝水足迹为39.0±4.8km³/年，其中灌溉农业占比最大（58.3%，24.7±3.1km³/年），其次是工业和能源开采（18.9%）、生活和城市用水（9.8%）以及生态用水（7.3%）。研究期间总蓝水需求量净增长19%，主要由工业用水强度增长（31%）和城市生活用水增长（44%）驱动，而生态用水分配则下降了8%。绿水足迹在黄土高原雨养农业区呈小幅下降（5%），与退耕还林工程实施后的植被恢复一致。年均灰水足迹为16.2±2.1km³/年，其中农业非点源污染占52%，工业废水占31%，生活污水占17%。灰水足迹在2000—2012年间的增速超过蓝水足迹，反映出化肥施用强度增加与工业废水排放的影响。各部门足迹的量级、占比、趋势与空间分布汇总于表2。

3.综合水压力指数：空间分布与时间轨迹

综合水压力指数（WSI）显示，黄河流域呈现出清晰的上游至下游加剧梯度，并叠加了西向东的社会经济需求放大模式。流域平均WSI从2000—2005年基准期的0.46增至2019—2024年的0.63，于2016年首次跨过高压力阈值（>0.6）。子流域层面，下游灌区WSI最高且上升最快，从0.72增至0.87，处于严重压力类别；关中平原走廊次之，从0.61升至0.79，反映城市工业需求在稳定蓝水供应基础上的快速增长；黄土高原子流域WSI绝对增幅最大（+0.21），主要由灰水压力上升驱动，凸显了灰水分指数的诊断价值；上游河段WSI最低，保持在0.18—0.27之间。表3展示了子流域WSI的分阶段变化，图3则展示了基准期与近期的网格尺度空间分布。图4展示了2000—2024年WSI的时空动态分布。图5展示了2000—2024年黄河流域WSI随机森林变量重要性评分。

4.随机森林驱动因素重要性与供需影响因素

随机森林模型在整个网格级年度面板数据上实现了R²=0.87和0.042WSI单位的RMSE，显示出较强的预测能力。重要性排序显示，灌溉面积分数是WSI空间变异性的最重要预测因子，其次为年降水量、地下水依赖性指数、土地利用类型和人口密度。供给侧驱动因素共同解释了41%的WSI方差，而需求侧的社会经济与土地管理因素共同解释了53%的方差。在下游灌区和关中平原，灌溉面积与地下水依赖性指数呈现高度非线性关系，前者在灌溉面积超过约0.45的阈值后对WSI的影响急剧上升；而在上游和黄土高原河段，年降水量主导WSI方差，反映出不同区域供给侧与需求侧控制作用的显著空间分异。

5.地理探测器因子交互作用分析与OLS残差的影响因素

地理探测器分析显示，灌溉面积分数和年降水量是单因子解释力最强的变量，q值分别为0.41和0.44。但交互作用探测结果更为关键：所有主要因子对均呈现双线性或非线性增强交互作用，其中灌溉面积分数与地下水依赖性指数的联合q值达到0.83，超过了二者单因子q值之和，说明地表水取用与地下水超采的叠加，是形成严重水资源压力区的核心耦合机制。人口密度与工业用水强度的交互作用（q=0.74）则集中在山西—陕西能源走廊与郑州—济南城市轴线，反映出家庭需求与工业取水在特定区域内形成了显著的叠加放大效应。表4总结了六个最具诊断信息价值的因子对的Geodetector q-值和交互作用分类。

Moran特征向量空间滤波的OLS回归进一步验证了上述发现。滤波后残差的Moran’s I从0.41降至0.06，空间自相关被有效吸收。标准化系数显示，灌溉面积分数是需求侧最强的预测因子（β=0.38），年降水量是供给侧最强的压力抑制因子（β=-0.31），地下水依赖性指数、人口密度与工业用水强度亦呈显著正相关。所有VIF值低于3.8，表明多重共线性控制良好。OLS残差图（图6）显示，观测值超过模型预测的正残差集中在河南与山东的下游灌区，提示现有变量集未能完全捕捉这些地区的地下水枯竭过程与跨省调水依赖等治理因素。

6.水足迹聚合中的尺度依赖偏差

为进一步表征聚合敏感性，进行了一项随机空间聚合实验。结果如图7所示，流域平均WSI的估计值在各层级均保持稳定，变异系数均低于0.06。然而，WSI的空间方差及其与蓝、绿、灰水组分之间的相关性对聚合层级更为敏感：单位级WSI标准差从最高分辨率层级降至最低层级时下降了38%，蓝水与灰水足迹的相关性也从0.61降至0.44。这表明，虽然平均水压力估计对空间聚合不敏感，但聚合过程会系统性地压缩压力分布的结构异质性，并削弱反映体积稀缺性与水质胁迫之间相对贡献的关键组分关联。

7.多组分水足迹动态与时空压力轨迹

图8展示了2000—2024年间全流域WSI及其组分分指数的时间演变，并按河段进行了分解。蓝—绿—灰水足迹分解揭示了传统供需比无法捕捉的流域动态。全流域WSI从2000—2005年的0.46升至2019—2024年的0.63，于2016年跨过高压力阈值。这一轨迹由三个并发过程驱动：蓝水农业开采量持续增长、城市与能源走廊工业及生活用水加速扩张，以及灰水同化能力逐步恶化—后者代表的基于水质的水资源需求在官方供水账户中未被记录，因而长期从政策诊断中缺失。上下游WSI差距从0.54扩大至0.60，反映出上游供应改善（源头区降水增加与生态恢复）与下游消费区需求增长之间的空间脱钩。黄土高原子流域呈现出独特轨迹，水压力上升主要来自灰水足迹扩张而非蓝水开采，源于农业化学投入集约化带来的面源污染压力。这种以灰水为主导的压力路径，与下游体积开采主导型压力在机制上本质不同，需要差异化的治理响应——从水量配额转向污染源控制与面源管理。

8.尺度差异化的驱动因素归因与治理启示

图9展示了所有12个预测因子的完整地理探测器单因子和交互作用q值矩阵。结果显示，最严重的赤字区域并非由单一因素决定。灌溉面积分数与地下水依赖性的非线性增强交互作用联合解释力达0.83，揭示地表水取用与地下水超采的共存是极端水压力的核心协同机制，尤其集中于下游宁夏、内蒙古灌区及河南—山东平原。这种耦合动态在缺乏明确地下水模块的传统水量平衡模型中容易被系统性低估。另一类交互作用体现在人口密度与工业用水强度之间，其双线性增强效应集中在山西—陕西能源走廊，表明人口驱动的家庭需求与能源部门工业取水在同一水文流域内相互强化，亟需整合性的城市—工业水治理。采用空间滤波的OLS残差分析进一步发现，河南与山东下游灌区存在系统性的正残差，即观测水压力超过模型预测，最合理的解释是南水北调等跨省调水在分配账户中降低了表观稀缺性，却未相应减少实际消耗压力，反映出单流域分析框架难以处理的核算错配问题。

9.水足迹组分交互作用与土地利用调节

灰水足迹对土地利用类型尤为敏感，图10展示了整个研究期间，六种主要的CLUD土地覆盖类别单位面积的蓝、绿、灰水足迹量级比较；图11展示了按主导土地覆盖类型、黄河流域的单位面积水足迹组成量大小：灌溉耕地单位面积灰水足迹是雨养耕地的3.8倍，是恢复草地的6.4倍；建设用地虽人均蓝水足迹最高，但单位面积绿水消耗最低，反映出土地利用将稀缺路径从绿水充足性转向蓝水取水强度。这些交互作用表明，黄河流域水压力管理不能仅局限于水文或社会经济维度，而需构建综合的土地—水治理框架，系统考虑土地利用转变对蓝、绿、灰三个足迹组分的协同影响。

2000年至2024年间，黄河流域水安全持续恶化，综合水压力指数从0.46升至0.63，于2016年首次跨过高压力阈值，下游灌区至2019年已达严重压力水平。这一趋势并非由单一因素驱动，而是源于社会—水文耦合交互作用，其中灌溉面积分数与地下水依赖性的非线性增强（联合q=0.83）最为关键，产生了任何单一因素都无法形成的压力协同效应，且持续存在于降水变率之外，反映出结构性而非气候性的水不安全根源。由此得出三点主要结论：其一，在异质性流域中，基于统一网格的蓝—绿—灰水足迹分解是实现可靠水压力诊断的前提——传统体积供需比无法识别灰水同化维度，会低估黄土高原等以灰水为主导的子流域压力，也难以区分性质不同的稀缺机制。其二，供给侧驱动因素在上中游主导水压力的方差，需求侧驱动因素在下游主导，这种尺度依赖的驱动结构要求分析框架能够同时解析两个尺度，避免在空间分辨率与因素归因之间取舍。其三，最严重的赤字区域由耦合驱动因子对共同控制，意味着针对孤立部门的治理框架效果将系统性低于解决多重压力源相互作用的综合干预措施。