2. 地下水模型边界条件基础

边界条件是地下水数值模型与外界发生水量交换的“接口”。模型域内部的参数（渗透系数、给水度、储水率）决定水怎么在介质中流动，而边界条件决定有多少水能进来、能出去，以及模型域与外部环境如何耦合。边界条件设置错误，往往不会导致程序报错，却会让模型收敛良好、水位曲线平滑，而水量平衡完全失真。因此，理解三类基本边界的数学表达与物理意义，是后续水量平衡核查的前提。

三类边界条件的数学表达与物理意义

在地下水动力学中，边界条件按控制变量分为三类。第一类是定水头边界（Dirichlet 边界，第一类边界），即边界上水头已知：

H(x, y, t) |边界 = H_b(x, y, t)

其物理意义是边界与一个水位恒定或按已知过程变化的外部水体相连，例如与河流、湖泊、水库直接连通且水力联系极强的地段，或已知水位观测井控制的断面。模型通过调整穿过边界的流量来维持给定水头，因此这类边界的流量是计算结果而非输入。这一点非常关键：定水头边界会自动“补”或“排”任意多的水，若设置范围过大或水头值不合理，会掩盖模型内部参数错误。

第二类是定流量边界（Neumann 边界，第二类边界），即边界上流量（通量）已知：

K · ∂H/∂n |边界 = q_b(x, y, t)

其中 n 为边界外法线方向，q_b 为单位面积补给或排泄量。物理意义是边界处的水量交换已知但水头未知，典型情形包括：侧向径流流入/流出断面（用达西定律由区域水力坡度估算）、降雨入渗或蒸发蒸腾的面状通量、已知抽水量的定流量抽水井。零流量边界（隔水边界）是其特例，q_b = 0，表示不透水底板或对称面。定流量边界下，边界水头是未知量，由模型求解得到。

第三类是混合边界（Cauchy 边界，第三类边界），即水头与流量以线性关系耦合：

q_b = C · (H_b − H)  或  K · ∂H/∂n = C · (H_ext − H)

式中 C 为边界传导系数（量纲为 1/时间或长度/时间，取决于表达形式），H_ext 为外部水体水位，H 为模型边界处水头。物理意义是边界通过一层弱透水介质与外部水体相连，交换量既取决于两侧水位差，也取决于这层介质的渗透性。典型对象是河床淤积层、湖底沉积层、弱透水断层带，以及用“河流—含水层”交换系数刻画的河渠边界。当 C 很大时趋近定水头，C 趋近 0 时趋近隔水边界，因此混合边界是介于两者之间的连续过渡。

边界条件设置中的常见错误辨析

错误一：把定水头边界当“万能补水源”。初学者常把模型四周全部设为定水头，认为这样“稳定、好收敛”。但定水头边界会按需补水或排水，若模型内部渗透系数偏大、给水度偏小，定水头边界会以异常大的交换量补偿，使水位拟合看起来不错，而水量平衡表中边界交换项占比畸高。判断方法是查看边界流入/流出量占总体补排的比例，若侧向定水头交换量远超降雨入渗或开采量，就应怀疑边界范围过大或内部参数失真。

错误二：定流量边界的正负号与法线方向不一致。定流量边界的 q_b 是沿边界外法线方向的通量，流入模型域为负、流出为正（或按软件约定相反）。实际建模中常见把“流入”写成正值、把“流出”写成负值，导致整个边界补排方向颠倒。这类错误在稳定流模型中表现为水位整体偏高或偏低，在非稳定流中表现为水位持续单调漂移。规避方法是：在设置前明确写出“外法线指向模型外”，并逐个边界核对流入/流出符号，再用总量守恒做一次手算校核。

错误三：混合边界的传导系数 C 随意取值。混合边界的 C 综合反映河床淤积层厚度、垂向渗透系数和河道湿周，具有明确物理含义，不能为了“拟合好看”随意调大调小。C 过大时边界退化为定水头，失去弱透水层的调蓄作用；C 过小时边界接近隔水，河流对含水层的补给被严重低估。合理做法是先由河床沉积物厚度与垂向渗透系数估算 C 的量级，再在率定中有限调整，并记录调整依据。

错误四：边界类型与物理对象不匹配。例如把一条有明显淤积层的河流设为定水头边界，等于忽略了河床阻力，会高估河流对含水层的补给；反之把与含水层直接连通的岩溶大泉设为混合边界，会人为削弱其排泄能力。辨析时应回到对象本身：边界两侧是否有弱透水介质？水头是否被外部水体钳制？交换量是否已知？三个问题分别对应混合、定水头、定流量三类边界。