1. 课程导入与目标设定

任何一个地下水数值模型，本质上都在求解同一件事：在给定的含水层结构与参数下，水从哪里进来、从哪里出去、水位如何响应。而“从哪里进来、从哪里出去”正是边界条件要回答的问题；“进来多少、出去多少、储变量变化多少”则是水量平衡核查要回答的问题。前者决定模型的物理合理性，后者决定模型数值结果的可靠性。二者缺一，模型就可能出现“看起来收敛、实际上错误”的危险状态。

在实际项目中，常见两类典型失败。第一类是边界条件设置不当：把定水头边界设在离研究区过近的位置，人为固定了本应自由波动的水位，导致降落漏斗被“压平”；或者把河流边界简化为定水头，忽略了河床淤积层带来的水力阻力，使河流补给量被严重高估。第二类是水量平衡未核查：模型虽然收敛，但某应力期总补给与总排泄之差与储变量变化不匹配，误差被隐藏在某个未察觉的边界项中。这两类问题往往不会让软件报错，却会让后续的预测、方案比选和资源评价全部失真。

为建立整体认识，先看下面这张课程知识结构图。它把本课程要讲的内容组织成一条从“认识边界”到“核查平衡”再到“定位误差”的主线，而不是零散的知识点堆叠。

上图把课程内容压缩为三层递进关系：先识别和设定边界条件，再用逐期水量平衡表检验模型是否“账目清楚”，最后在账目不平或单项异常时反查边界与参数。三层之间不是并列关系，而是前者为后者提供输入、后者为前者提供验证的闭环关系。

为什么先讲边界，再讲平衡

边界条件是模型的“外部世界接口”。以常见的 MODFLOW 类模型为例，定水头边界（CHD 或常水头单元）会直接向模型注入或抽出水量，其流量不是事先给定的，而是求解后由水头差和导水系数计算得到；定流量边界（WEL）则相反，流量给定、水头自由。这两类边界的性质完全不同：前者是“水位受控”，后者是“流量受控”。如果把一个实际受区域水位控制的断面误设为定流量，模型就失去了对外部水位变化的响应能力；反之，把一条流量已知的排水沟设为定水头，则会人为制造一个无限大的“水库”。

正因为边界流量是求解结果而非输入，水量平衡核查才显得关键。模型算完后，软件通常会输出各应力期的平衡项：补给、排泄、各边界进出流量、储变量变化。核查的逻辑是：总流入 − 总流出 = 储变量变化。若等式两边不闭合，说明存在数值误差或某项被遗漏；若闭合但某一项数值明显违背物理常识（例如干旱期河流补给量反而增大），则说明边界或参数设置有问题。前者是“账不平”，后者是“账平但账目可疑”，两类问题都需要通过核查发现。

本课程的学习路径与实战任务

课程将沿着“概念—类型—设置—核查—排障—复盘”的顺序推进，但不会把每一节都做成同样的步骤清单。具体而言，你会先掌握各类边界条件的物理含义与适用场景，再学习如何在模型中正确表达它们；随后进入水量平衡核查，学会编制逐应力期平衡表、设定误差判据、识别异常项；最后通过一个完整案例，把边界审查与平衡核查串成一次可复现的实战。

本课程的实战任务设定如下（示例参数，用于教学演示，非实测结果）：某三层含水层模型，研究区面积约 120 km²，模拟期 10 年，划分为 20 个应力期。模型边界包括：北侧定水头边界（区域水位控制）、南侧河流边界（考虑河床阻力）、东侧隔水边界、西侧定流量开采边界。任务要求你完成三件事：第一，逐条审查四类边界的类型、位置与参数依据是否合理；第二，编制 20 个应力期的水量平衡表，计算相对误差；第三，对相对误差超过 1% 或单项流量突变的应力期进行定位，给出修正建议。

为帮助你建立“什么样的结果算通过”的直观标准，下面给出一个简化的平衡核查示例表（数值为教学假设，仅示意格式与判读方式）：

这张表的价值不在于数字本身，而在于它体现的核查习惯：每个应力期都要有闭合差和相对误差两列，且要对“数值通过但物理可疑”的情形单独标注。例如第 3 期储变量变化接近零，意味着该期补给与排泄几乎相等，若此时河流边界流量却出现明显跳变，就应回头检查河流边界的水位或河床参数是否被误设。