2. 地下水模型参数体系梳理

在建立地下水数值模型之前，参数体系的梳理是决定模型能否反映真实水文地质条件的基础工作。参数敏感性分析与不确定性分析都建立在参数体系清晰、量纲统一、空间归属明确的前提之上。本节从参数分类、空间尺度与变异性、可识别性三个层面，系统梳理建模中常用的水文地质参数，并给出可操作的判断方法。

一、水文地质参数的分类

地下水模型中的参数可按其物理意义分为三大类：含水层传导与储水参数、溶质运移参数、源汇项参数。三类参数在模型中的角色不同，敏感性表现也差异明显。

1. 传导与储水参数。渗透系数 K（m/d）描述介质导水能力，是控制地下水流场形态最核心的参数；在各向异性介质中还需给出 Kx、Ky、Kz 三个方向分量。给水度 μ（无量纲）用于潜水含水层，表示单位面积水位下降时单位面积释放的重力水体积，典型砂砾层为 0.15～0.30，粉砂层为 0.05～0.15。储水系数 S（无量纲）用于承压含水层，表示单位水头下降时单位面积含水层释放的水量，典型值为 1×10⁻⁵～1×10⁻³。导水系数 T=K·M（M 为含水层厚度）常作为整体参数参与识别。

2. 溶质运移参数。弥散度 αL（纵向）、αT（横向）描述机械弥散强度，纵向弥散度常与尺度相关，室内试验为厘米级，区域尺度可达数十至数百米；有效孔隙度 ne 影响实际流速；阻滞因子 R 与分配系数 Kd 反映吸附作用。

3. 源汇项参数。包括降水入渗系数 α、河流渗漏系数、蒸发极限深度与蒸发系数、开采量 Q、越流系数 K′/M′ 等。这类参数往往具有强时空变化，且与边界条件耦合，是模型校准中最易出现补偿误差的部分。

二、参数的空间尺度与空间变异性

水文地质参数本质上是尺度依赖的。室内土柱测得的 K 代表厘米级均质体，抽水试验解析解得到的是影响半径范围内的等效 K，而区域模型网格（如 500 m×500 m）所需的 K 则是网格尺度上的等效值。三者不可直接混用。当介质存在明显非均质性时，等效渗透系数还取决于水流方向：平行层理方向为算术平均，垂直层理方向为调和平均，这解释了为何同一含水层水平与垂直渗透系数常相差 1～2 个数量级。

空间变异性的处理方式通常有三种：分区赋值（按水文地质单元划为若干参数区，区内取常数）、趋势面/插值赋值（由实测点插值到网格）、随机场方法（以均值、方差、相关长度描述，用于不确定性分析）。分区赋值简单但会人为制造参数突变界面；随机场方法能刻画渐变与非均质，但需要足够的实测样本支撑统计特征。实际建模中常采用“分区为主、区内插值为辅”的折中方案。

三、参数可识别性的初步判断

可识别性指仅凭现有观测数据能否唯一确定参数值。判断可识别性不必等到反演结束，建模初期即可通过以下三条线索做初步筛查。

第一，看观测对参数的敏感方向。若某参数变化 20% 而模拟水头变化小于观测误差（如 0.1 m），则该参数在当前观测条件下基本不可识别。例如给水度主要影响水位随时间的变化速率，若模型仅有稳态观测，给水度几乎不可识别；储水系数在承压含水层中影响瞬态响应，若无长序列水位动态，同样难以约束。

第二，看参数间是否存在结构性相关。导水系数 T 与给水度 μ 在抽水试验降深曲线中常同时影响曲线形态，二者存在相关；K 与含水层厚度 M 以乘积形式出现，单独识别困难；弥散度 αL 与孔隙度 ne 在穿透曲线中相互补偿。当两个参数在敏感矩阵中列向量近似平行时，应固定其中一个或引入额外观测（如示踪试验、分层水位）。

第三，看参数分区数与观测点数的匹配。若某参数划分了 8 个分区而区内仅有 2 个观测孔，则该参数的自由度过高，反演结果必然不唯一。经验上每个待识别参数分区至少应有 3～5 个有效观测约束，否则应合并分区或将该参数设为经验值固定。

四、参数体系梳理的落地做法

建议在建模准备阶段完成一张参数梳理表，字段包括：参数符号、物理意义、量纲、初值及来源、空间分区方案、是否参与反演、可识别性初判（高/中/低）、备注。对初判为“低”的参数，优先固定为经验值或文献值，仅对“高”“中”参数开展敏感性分析，可显著降低后续计算量并提高分析结论的可信度。参数体系梳理清楚之后，敏感性分析才有了明确的输入对象和合理的取值范围，这正是下一阶段工作的起点。