3. 水量平衡原理与核查指标
地下水模型无论采用有限差分还是有限元方法，其数值解在物理上必须服从质量守恒。水量平衡核查就是把这个物理约束变成可量化、可比较的验收指标：模型算出来的每一口抽水、每一条河流的渗漏、每一个边界的流入流出，最终都要能和储量变化对上账。如果对不上，说明模型存在数值误差、边界条件设置错误或参数不合理，必须先排查再用于预测。
一、水量平衡方程的基本形式
对任意一个均衡期（一个应力期或整个模拟期），把模拟域视为一个控制体，水量平衡方程可写为：
ΣQin − ΣQout = ΔS
即：总流入量 − 总流出量 = 储量变化量
这里需要先明确符号约定：本文统一采用流入为正、流出为负的记账方式，储量增加为正、减少为负。这样方程可以写成更通用的形式：
ΣQ_i = ΔS
其中 Q_i 为各均衡项（流入为正、流出为负），ΔS 为均衡期内储量变化在 MODFLOW 类模型中，储量变化项由 STORAGE 子程序包计算，各边界项分别由 WEL、RIV、CHD、GHB、DRN、RCH、EVT 等子程序包累计。模型输出文件（如 MODFLOW 的 .lst 列表文件）会按应力期和整个模拟期分别给出各项的体积量（单位通常为 m³ 或 ft³），核查时直接读取这些累计值即可。
二、各均衡项的物理定义
均衡项不是抽象符号，每一项都对应模型中一类具体的边界或源汇机制，理解其定义是正确核查的前提：
- 定水头边界（CHD）：维持指定水头的边界，流入或流出量由模型求解得到，是模型与外部环境交换水量的主要通道，符号取决于水头梯度方向。
- 河流（RIV）：通过河床弱透水层与含水层发生渗漏交换，量值由河水位与地下水位的差、河床导水系数共同决定，河水位高于地下水位时为流入。
- 排水沟（DRN）：仅在地下水头高于沟底高程时排水，是单向流出项，不产生流入。
- 通用水头边界（GHB）：以外部水头和导水系数描述的半透水边界，可正可负，常用于侧向补给或排泄。
- 抽水井（WEL）：人为开采，通常为流出项；回灌井则为流入项。
- 补给（RCH）与蒸发（EVT）：面状源汇，补给为流入，蒸发为流出且受埋深控制。
- 储量变化（ΔS）：均衡期内含水层释放或储存的水量，由给水度或储水率与水位变化共同决定。
需要强调的是，同一个物理过程在不同模型中可能被归入不同项，例如河流既可用 RIV 也可用 GHB 近似，核查时必须按模型实际启用的子程序包逐项对应，不能凭物理直觉想当然。
三、相对误差的计算方法
绝对误差（ΣQin − ΣQout − ΔS）的量级与模型规模有关，不便于跨模型比较，因此工程上普遍采用相对误差（也称水量平衡误差百分比）：
Discrepancy (%) = |ΣQ_in − ΣQ_out − ΔS| / ( (|ΣQ_in| + |ΣQ_out|) / 2 ) × 100%分母取流入绝对值与流出绝对值的平均值，是为了避免流入流出相互抵消导致分母过小、误差被虚假放大。MODFLOW 列表文件中给出的 "Percent Discrepancy" 即按此思路计算。实际核查时，应分别检查每个应力期和整个模拟期的相对误差，因为某些应力期（如抽水启动初期）可能出现较大的瞬时不平衡，而全期累计值可能掩盖局部问题。
四、允许阈值的设定
阈值设定没有唯一标准，需结合模型规模、非线性程度和用途综合确定。一般经验是：
| 核查对象 | 建议阈值 | 说明 |
|---|---|---|
| 单个应力期相对误差 | ≤ 1% | 常规稳定流或缓变流问题 |
| 整个模拟期累计相对误差 | ≤ 0.5% | 作为模型验收的基本门槛 |
| 强非线性工况（干湿交替、承压转无压） | 可放宽至 1%～2% | 需在报告中说明原因 |
| 预测情景 | 应比率定期更严 | 避免误差在长序列中累积 |
阈值一旦设定，应写入模型验收标准并在报告中固定下来，不能因为结果不达标就临时放宽。若相对误差超限，排查顺序通常是：先看收敛判据是否过松（如 solver 的 head change 和 residual 设置），再看时间步长是否过大导致储量项计算失真，然后检查边界条件是否存在物理上矛盾（如 CHD 与 WEL 在同一网格相互抵消），最后核对各子程序包的量纲与单位换算是否一致。
五、核查的实际操作要点
在具体执行时，建议按以下顺序读取和核对数据。以 MODFLOW 为例，打开 .lst 文件后先定位 "VOLUMETRIC BUDGET FOR ENTIRE MODEL" 段落，逐项抄录 IN、OUT 和 STORAGE 的累计值，再按上面的公式手算一遍相对误差，与程序输出的 Percent Discrepancy 对照，确认自己理解的口径与程序一致。这一步看似繁琐，但能有效发现单位换算错误——例如补给量按 mm/d 输入而模型长度单位是 m，量级会差三个数量级，相对误差会异常放大。
其次要区分"数值误差"和"物理不平衡"。数值误差通常随收敛判据收紧而减小；如果收紧 solver 设置后相对误差几乎不变，说明问题出在边界条件或参数本身，例如同一位置既有定水头又有抽水井、河流参数导致水量在河与含水层之间反复振荡。此时应回到边界条件章节的方法逐项复核，而不是继续调求解器。
实用建议：把每个应力期的流入、流出、储量变化和相对误差整理成一张核查表，随模型文件一同归档。这样在后续参数调整或情景预测时，可以快速对比哪一项发生了变化、变化是否合理，也便于向审查方说明模型的水量闭合情况。
最后需要提醒，水量平衡通过只是模型可用的必要条件，不是充分条件。一个水量完全闭合的模型仍可能因为边界位置错误、参数分区不合理而给出错误预测。因此平衡核查应与水位拟合、参数敏感性分析结合使用，共同构成模型验收的证据链。