伊朗伊斯法罕水-食物-能源纽带方法用于地面沉降管理

《Journal of Hydrology: Regional Studies》：A water–food–energy nexus approach to land subsidence management in Isfahan, Iran

　　

研究区域：伊朗中部伊斯法罕-博尔哈尔平原（Isfahan–Borkhar Plain）是一个干旱、依赖地下水的流域。密集的农业开采、反复干旱以及地表水流入量减少已导致含水层枯竭并引发严重地面沉降。 研究焦点：本研究开发了一个集成的水-食物-能源（WFE）纽带模拟-优化框架，以量化地面沉降减缓、农业收入与地下水抽能需求之间的权衡。一个系统动力学（SD）模型区分了子系统内相互关系（子系统内过程）与跨部门相互作用（跨子系统资源交换），并与基于非支配排序遗传算法II（NSGA-II）的多目标优化方法耦合，以推导帕累托有效管理组合。能源维度通过地下水抽能需求和电价假设来表示，而非内生的电力供应动态。流域尺度地面沉降被建模为对地下水位下降的不可逆响应，并用作筛选指标。 区域新水文见解：在一切照旧（BAU）情景下，预计到2060年地下水枯竭量将超过400亿立方米，而累积地面沉降预计将达到约3.5米。结果表明，仅提高灌溉效率是不够的，因为减少的回归流会降低含水层补给。模拟表明，大幅减缓需要将需求管理与人工补给含水层（MAR）相结合。一个平衡组合，包括农业地下水使用量减少40%和每年1.5亿立方米的人工补给含水层，可将地面沉降限制在小于0.5米，并将抽能需求降低50%以上。该框架揭示了纽带权衡并支持流域尺度筛选。

论文解读：伊朗伊斯法罕平原水-食物-能源纽带框架下的地面沉降管理研究

研究背景与意义

伊朗中部伊斯法罕-博尔哈尔平原（Isfahan–Borkhar Plain）是一个干旱且高度依赖地下水的流域。大都市伊斯法罕是世界上沉降速率最快的城市区域之一，合成孔径雷达干涉测量（InSAR）观测到的城市区域变形速率为5至100毫米/年，邻近平原官方报告甚至接近40厘米/年，远超其他地区采用的规划阈值。导致这一危机的根本原因在于数十年来基本不受监管的地下水开采，以支持灌溉农业和城市化，同时叠加了反复干旱以及扎扬德鲁德河（Zayandeh-Rud River）等地表水资源的长期管理不善。长期地下水抽取不仅降低了孔隙压力、压实易变形沉积物并产生不可逆的地表形变，还会破坏基础设施并减少含水层储水能力。由于驱动沉降的抽水行为受灌溉农业和抽水所需能源的维持，仅从地下水单一视角无法评估减缓措施的有效性。水-食物-能源（WFE）纽带框架为这一问题提供了合适视角，因为水资源、能源和粮食安全是相互依存的政策领域。然而，现有文献中存在明显空白：地面沉降诊断研究多为描述性，未提供集成政策模拟或优化框架；WFE纽带研究通常未将地面沉降作为主要输出或目标函数；地下水管理优化研究则往往未充分整合食物和能源部门的反馈回路。为填补这些空白，研究人员在《Journal of Hydrology: Regional Studies》上发表论文，开发了一个集成模拟-优化框架，旨在筛选技术上可实现的地下水管理组合，以减缓伊斯法罕平原的地面沉降。

主要技术方法

研究人员采用系统动力学（System Dynamics, SD）建模方法，在Vensim决策支持系统（DSS）软件中构建了包含水、食物（农业）和能源三个核心子系统及独立地面沉降模块的因果反馈结构，明确区分了子系统内相互关系与跨部门相互作用。模型利用2005年至2019年的年尺度外生强迫数据进行参数化，以2005年至2014年数据校准，2015年至2019年数据验证。验证后的SD模拟器通过Python接口与基于非支配排序遗传算法II（NSGA-II）的多目标优化（MOO）引擎耦合，以最小化累积地面沉降、最大化农业收入和最小化抽能需求为目标函数，生成收敛的帕累托有效管理组合。决策变量涵盖农业地下水配额、灌溉效率目标、作物种植结构转变、人工补给含水层（MAR）量和用水价格。

研究结果

1. 1.
   引言
   引言部分通过文献综述指出，地面沉降不仅是地表形变现象，更是地下水系统压力和潜在储水能力丧失的重要指标。研究人员强调，有效的沉降减缓需要解决地下水枯竭的根本驱动因素，而非仅针对形变症状。伊斯法罕平原的案例凸显了平衡水文地质恢复与经济连续性的必要性。论文明确了四个具体目标：构建WFE纽带SD模型、验证模型并将地面沉降作为主要输出、集成MOO算法、生成并分析帕累托有效前沿以量化环境恢复与社会经济稳定之间的权衡。
2. 2.
   材料与方法
   2.1 方法框架概述
   方法开发分为两个阶段。阶段1开发并验证包含四个耦合模块（水、食物、能源、地面沉降）的SD模型，并区分系统内相互关系与跨系统相互作用。阶段2将验证后的模拟器与NSGA-II耦合，识别帕累托有效管理组合。

2.2 研究区域：伊斯法罕-博尔哈尔平原

该平原位于扎扬德鲁德河集水区内，是一个无承压冲积含水层，历史上接受降雨入渗、灌溉回归流和河流渗漏补给。近年来，持续干旱、上游水库调节和跨流域调水减少了下游河流流量，使区域水平衡转向地下水开采，加速了地下水位下降和地面沉降。

2.3 SD模型开发

2.3.1 概念模型与因果回路图

概念模型识别了关键反馈回路：一个主导的增强回路“抽水-压力回路”（R1）描述了农业需水增加导致地下水开采量上升、水位下降、单位提水能耗增加，而补贴能源价格削弱了价格信号，使开采持续；一个平衡回路“沉降-成本回路”（B1）反映了水位下降引发沉降、基础设施受损并产生财务成本，进而形成减少抽水的压力。

2.3.2 存量与流量模型公式

水子系统以地下水储量（S_GW）为核心，根据质量平衡原理，总补给量（R_total）包括降水补给（R_precip）、河流渗漏补给（R_river）、灌溉回归流补给（R_ret）和市政补给（R_muni），总开采量（E_total）包括农业（E_agri）、市政（E_muni）和工业（E_ind）开采。灌溉回归流补给（R_ret）与灌溉效率（η）直接相关，效率提升会减少回归流，进而降低补给，揭示了效率-补给权衡。食物子系统计算总农业需水量（AWD），由作物组面积（A_i）、作物需水量（CWR_i）和灌溉效率（η_i）决定，并引入单位净作物需水量毛利（GM_water）指导作物种植结构优化。能源子系统计算抽能消耗（EC_pump），与总开采量（E_total）和提水扬程（H_lift）成正比，扬程随地下水位（H_GW）下降而增加。地面沉降模块将流域尺度沉降建模为对地下水位下降的不可逆响应，用作筛选指标。

2.4 数据与参数化

研究使用了2005年至2019年的年尺度外生强迫数据，包括降水、河流流量、作物面积和人口统计等，参数化参考了FAO-56 Penman-Monteith方法估算的潜在蒸散发。

2.5 模型校准与验证

模型以2005年至2014年为校准期，2015年至2019年为验证期，采用标准统计指标评估模拟性能，确保模型在历史观测下的可靠性。

2.6 多目标优化

2.6.1 目标函数

三个冲突目标分别为最小化累积地面沉降、最大化农业收入和最小化抽能需求。

2.6.2 决策变量

包括农业地下水配额、灌溉效率目标、作物模式转变、人工补给含水层（MAR）量和用水价格。

2.6.3 NSGA-II优化引擎

通过Python接口将SD模型与NSGA-II耦合，生成帕累托有效前沿。

1. 3.
   结果与讨论
   3.1 模型校准与验证结果
   验证结果表明SD模型能够可靠地模拟历史地下水水位和地面沉降趋势，为后续优化提供可信的模拟平台。

3.2 一切照旧（BAU）情景

BAU情景下，预计到2060年地下水枯竭量将超过400亿立方米，累积地面沉降将达到约3.5米，突显了现行管理方式的不可持续性。

3.3 帕累托有效前沿与权衡分析

优化结果显示，帕累托有效前沿揭示了三个目标之间的非线性权衡。仅提高灌溉效率的方案因减少灌溉回归流补给而效果有限。一个平衡管理组合，即农业地下水使用量减少40%配合每年1.5亿立方米的人工补给含水层（MAR），可将累积沉降限制在小于0.5米，同时使抽能需求降低50%以上，但会对农业收入产生一定影响。

3.4 管理启示

研究强调，有效的沉降减缓必须结合需求管理与人工补给含水层（MAR），以应对灌溉效率与含水层补给之间的耦合关系。能源共效益显著，抽水量减少直接降低了能源消耗。框架为决策者提供了技术可行的政策组合及其权衡的清晰图谱。

结论

研究得出结论，在伊斯法罕-博尔哈尔平原，一切照旧（BAU）管理将导致严重的地下水枯竭和约3.5米的累积地面沉降。仅提高灌溉效率不足以缓解沉降，因为减少的回归流会降低含水层补给。模拟表明，大幅减缓需要将需求管理与人工补给含水层（MAR）相结合。一个平衡组合，包括农业地下水使用量减少40%和每年1.5亿立方米的人工补给含水层（MAR），可将沉降限制在小于0.5米，并将抽能需求降低50%以上。该水-食物-能源（WFE）纽带模拟-优化框架成功量化了地面沉降减缓、农业收入与抽能需求之间的权衡，为干旱区地下水管理提供了集成评估工具。