3. 污染羽空间插值的基本假设与工作流
空间插值不是把监测点连成漂亮色块的美化工具，它的本质是用有限的离散观测去推断未采样位置的浓度。任何插值结果都隐含了一组关于“空间相关性如何随距离变化”的假设。若假设与场地实际不符，等值线再平滑也可能是错误结论。本节围绕平稳性、各向异性、变异函数三个核心概念，梳理从数据准备到成果检验的完整工作流，并指出工程中最常见的误区。
一、平稳性：插值成立的前提
平稳性（stationarity）指污染羽的统计特征在空间上不随位置系统变化。最常用的是二阶平稳假设：区域内浓度的均值近似恒定，方差有限，且任意两点间的协方差只取决于它们的距离与方向，而与绝对位置无关。弱一些的是内蕴假设：只要求增量的方差平稳，即变异函数只依赖滞后距，不要求均值恒定。普通克里金采用内蕴假设，简单克里金则要求已知均值。
实际场地几乎不可能完全平稳。污染羽常因源区、地下水流向、地层界面而呈现趋势：靠近源区浓度高、下游逐渐降低。此时若直接对原始浓度做普通克里金，变异函数会被大尺度趋势“污染”，表现为远处点仍显示强相关，插值结果出现虚假的平滑拖尾。处理办法是先去趋势：用一阶或二阶多项式拟合区域趋势面，对残差做变异函数分析与克里金，最后把趋势面加回。判断是否需要去趋势，可先绘制浓度对坐标的散点图或分方向均值图，若沿地下水流向存在明显单调变化，就应警惕。
二、各向异性：方向不能忽略
污染物沿地下水流动方向迁移远、垂向扩散弱，因此空间相关性常呈几何各向异性：沿流向相关范围大，垂直流向相关范围小。变异函数模型需分别给出主变程（major range）与次变程（minor range），并用方位角描述主轴方向。若忽略各向异性、强行用各向同性模型，插值会把垂直流向的浓度“拉平”，低估羽状体长度或高估宽度。
确定主轴方向的实用做法是：按不同方位角（如每 15°或 30°）分组计算实验变异函数，比较各方向相关范围，相关范围最大的方向即主变程方向。该方向通常与地下水流向一致，但受地层透镜体影响也可能偏离，应以数据为准而非先验假定。示例中若主变程 800 m、次变程 200 m、方位角 60°，则各向异性比约 4:1，插值前必须写入模型参数。
三、变异函数：把“距离—差异”量化
实验变异函数 γ(h) 定义为滞后距 h 的两点浓度差平方的一半的平均值：
γ(h) = 1/(2N(h)) · Σ [Z(xi) − Z(xi+h)]²其中 N(h) 是距离落在 h 附近容差带内的点对数。随 h 增大，γ(h) 先上升，达到基台值（sill）后趋于稳定；基台对应的距离称变程（range），表示超出该距离后两点不再相关；h→0 时的截距为块金值（nugget），反映采样误差与小于最小采样间距的微尺度变异。常用模型有球状、指数、高斯模型，需用最小二乘或加权最小二乘拟合，并检查块金值与基台值之比：若块金比过大（如超过 50%），说明短尺度变异强、数据噪声大，插值可信度下降。
实验变异函数的点对数量随距离增大迅速减少，远距离段估计不稳定，拟合时应限制最大滞后距（常取最大采样距离的一半左右），并保证每个滞后距分组有足够点对（如不少于 30 对）。
四、从数据到成果的完整流程
第一步是数据清洗。统一坐标系统与高程基准，剔除采样记录错误、重复点与明显异常值；对低于检出限的数据，不能简单填 0，应记录为区间删失数据，或在变异函数分析中采用稳健估计。第二步做趋势与方向诊断，绘制浓度—坐标散点图、分方向均值图，判断是否需要去趋势、主轴大致方位。第三步去趋势，用多项式或广义最小二乘拟合趋势面，保留残差。第四步变异函数建模：按方位角分组计算实验变异函数，确定各向异性参数，拟合球状或指数模型，得到块金、基台、变程。第五步插值：选择普通克里金或泛克里金，在网格上估计浓度及克里金方差。第六步交叉验证：采用留一法，统计预测误差均值（应接近 0）、均方根误差与标准化误差（标准差应接近 1），若标准化误差方差明显偏离 1，说明变异函数参数或趋势设定不当，需回到第四步调整。
以某场地为例（示例参数，非实测）：监测井 24 口，主变程 800 m、次变程 200 m、方位角 60°，块金 0.15、基台 1.0（浓度取对数后），拟合球状模型。留一交叉验证得到均方根误差 0.42（对数单位），标准化误差方差 1.08，说明模型基本合理；若某口井残差持续偏大，应检查该点是否受局部源或采样误差影响。
五、常见误区
- 把插值当实测：克里金给出的是估计值与不确定性，远离开采区的位置克里金方差大，不能当作确定浓度用于风险判定。
- 忽略各向异性：直接用各向同性模型，导致羽状体形态失真，尤其影响下游边界判定。
- 变异函数“拟合好看”就通过：只看曲线贴合、不看块金比与交叉验证，参数可能物理意义不合理。
- 用全部数据外推：在数据覆盖范围外插值，方差迅速增大，结果不可信，应明确标注外推区。
- 浓度未变换：强偏态浓度直接插值易受高值点支配，通常先取对数再建模，成果再反变换。