4. 边界条件设定前的资料准备

边界条件不是凭空“画”在地图上的一条线,它必须由水文地质概念模型支撑。概念模型回答三个问题:含水层系统是什么样子、水从哪里来又到哪里去、模型范围四周与外界如何交换水量。只有这三件事在资料层面站得住,边界位置与性质(定水头、定流量、通用水头、隔水、河流、潜水面等)才有依据。本节按“先建骨架、再理收支、后定边界”的顺序,梳理设定边界前应准备的数据,并给出可直接对照使用的数据清单。

一、概念模型构建:先搭出含水层系统的骨架

概念模型的核心是把地层、构造与水力特征组织成可计算的结构。需要准备的第一类资料是空间结构数据:含水层与弱透水层的顶底板高程、厚度、垂向分层关系,以及基岩面起伏。这些数据通常来自钻孔柱状图、剖面图、区域地质图与物探解译成果。第二类是水文地质参数:渗透系数(或导水系数)、给水度与储水系数、垂向越流系数,以及参数的分区边界。参数不必一开始就精确,但必须给出初值与合理区间,并注明来源(抽水试验、经验值、类比区)。

第三类是水位与流场信息:初始水头分布、等水位线图、地下水流向与水力坡度。初始水头既是模型的初始条件,也是判断边界性质的重要参照——例如区域流场中水头等值线近似垂直模型边界时,该边界往往适合处理为隔水边界或零流量边界;等值线近似平行且穿过边界时,则更可能是流量交换边界。第四类是系统概化范围:模型边界应尽量落在天然水文地质边界(分水岭、隔水断层、河流、海岸线)上,若不得不采用人为边界,则需说明其合理性并准备敏感性分析。

概念模型不是一次性成果。资料补充后应回头修正分层、参数分区与边界位置,形成“资料—概化—检验”的迭代记录,这也是后续水量平衡核查能追溯依据的前提。

二、源汇项识别:把每一项收支落到数据上

源汇项识别要求逐项列出进入和离开含水层系统的水量,并为每一项准备“位置、强度、时间变化”三类数据。常见项包括:降水入渗补给、灌溉回归补给、河流与湖泊的渗漏或排泄、侧向流入流出、人工开采、蒸发蒸腾排泄、越流交换。以降水入渗为例,需要降雨量系列、入渗系数分区及其随土地利用类型的变化;以开采为例,需要井位坐标、开采层位、逐月或逐年开采量,以及开采井的过滤段深度。

识别过程中最容易出问题的是“量级判断”:某些项虽然存在但数量级很小,可以暂缓细化;某些项看似次要,却可能是边界附近的主要交换量。建议为每一项标注数据可靠度(实测/估算/经验)与时间分辨率(年/月/日),并做一次粗略的收支排序。下面的图示把概念模型、源汇项与边界判定所需的数据组织成一条从资料到边界决策的关系链,便于对照检查缺项。

边界设定前的资料组织与判定关系 空间结构资料 钻孔柱状图 / 剖面图 含水层顶底板高程、厚度 基岩面起伏、断层位置 参数与水位资料 渗透系数 / 导水系数分区 给水度、储水系数、越流系数 初始水头、等水位线、流向 源汇项资料 降水入渗、灌溉回归 河流渗漏 / 排泄、蒸发 开采井位、层位、开采量 侧向流入流出估算 水文地质概念模型 分层结构 + 参数分区 初始流场 + 系统范围 收支项清单与量级排序 数据可靠度标注 (迭代修正) 边界位置判定 优先落在天然边界 分水岭 / 隔水断层 / 河流 人为边界需说明合理性 边界性质判定 定水头 / 定流量 / 通用水头 隔水 / 河流 / 潜水面 依据:流场方向与交换量 数据清单核对 位置、强度、时间变化 可靠度:实测/估算/经验 缺项标注与补充计划 为水量平衡核查留接口

三、边界位置与性质判定:用数据说话

边界位置判定优先遵循“天然边界优先”原则。若模型东侧存在一条常年有水的河流且河水位与地下水位存在稳定水力联系,可将其处理为定水头边界或河流边界;若边界处为隔水断层,则处理为零流量边界。对于没有天然边界可依的情况,人为边界应尽量远离研究重点区,并用区域等水位线图验证边界处流场方向是否与假定一致。

边界性质判定需要三类证据:一是流场证据,即边界附近地下水流向与水力坡度;二是交换量证据,即通过达西定律估算的侧向流量量级;三是动态证据,即边界附近观测孔水位随时间的变化是否与假定边界一致。例如假定某侧为定水头边界,则边界附近观测孔水位应相对稳定;若观测到明显趋势性下降,则说明该边界可能并非定水头,或存在未识别的开采排泄。下面给出一份可直接填写的数据清单示例,字段按本节任务设置,参数为教学示例值。

资料类别 具体数据项 示例值 / 说明 可靠度
空间结构 含水层顶板高程 示例:28.5 m(钻孔 ZK03) 实测
空间结构 含水层厚度 示例:12~18 m,东薄西厚 实测/推断
参数 渗透系数分区 示例:K=8.5 m/d(冲积层)、1.2 m/d(坡积层) 抽水试验
水位 初始水头 示例:西部 35.2 m,东部 29.8 m 统测
源汇项 降水入渗系数 示例:0.18(耕地)、0.09(建设用地) 经验/估算
源汇项 开采量 示例:单井 800 m³/d,共 14 眼 计量
边界 东侧河流水位 示例:29.6 m,年变幅 0.4 m 水文站
边界 西侧分水岭位置 示例:地形分水岭与地下水分水岭基本重合 推断

四、资料准备的常见缺口与处理

实际项目中最常见的缺口有三类。第一类是参数分区缺少依据,此时应先用区域经验值建立初值,并在模型中设置参数敏感性分析,而不是反复“调参”掩盖资料不足。第二类是源汇项时间分辨率不足,例如只有年开采量而无月分配,可按灌溉季节或用水规律做合理拆分,并明确标注为估算。第三类是边界附近缺少观测孔,导致边界性质无法验证,此时应补充临时水位统测,或在模型中把该边界设为通用水头边界并做情景对比。

不要把“资料不足”直接转化为“边界随便设”。边界性质一旦错误,后续水量平衡核查会出现系统性偏差,且难以通过参数调整纠正。资料缺口应显式记录在案,作为模型不确定性的组成部分。

完成上述准备后,应形成一份可追溯的资料台账:每项数据注明来源、时间、精度与可靠度,每个边界注明位置依据与性质依据。这份台账既是边界条件设定的输入,也是后续水量平衡核查中解释“为什么收支不平衡”的第一手依据。资料准备做得扎实,边界设定与平衡核查才不是两张皮,而是同一套概念模型下的前后衔接。