4. 污染羽范围圈定方法与边界判定

污染羽范围圈定的本质,是把离散的监测井浓度数据转化为连续的空间浓度场,再按某一浓度阈值截取边界。因此边界位置并不唯一,它同时受三方面控制:插值方法如何"填补"井间空白、等值线如何绘制与闭合、以及用什么浓度值作为边界判据。理解这三者的耦合关系,才能解释为什么同一批数据在不同报告里会画出大小不同的污染羽。

一、插值方法:不同假设产生不同井间浓度

监测井是点状采样,井间浓度必须靠插值推断。常用方法各有隐含假设,直接影响污染羽的"胖瘦"。

反距离权重法(IDW)假设浓度随距离平滑衰减,幂指数 p 越大,高浓度点影响越局部。p=2 时井间过渡平缓,污染羽边缘较"鼓";p 增大到 4~5,等值线会更贴近高值井,羽状范围收缩。IDW 不产生超过实测最大值的极值,适合数据较密、梯度平缓的场地。

克里金法基于变异函数建模空间相关性,能给出估计方差。普通克里金在数据稀疏处会向局部均值回归,可能压低峰值、抬高边缘;指示克里金则直接对"是否超过阈值"建模,特别适合圈定超标范围。克里金的关键在于变异函数拟合是否合理,块金值过大意味着短距离随机性强,插值结果趋于平滑。

自然邻点法与最近邻法属于局部插值。最近邻法把每个网格点赋为最近井的值,结果呈阶梯状多边形(泰森多边形),边界生硬但对数据不做额外假设;自然邻点法过渡较自然,但外推能力弱。

实操中应至少用两种方法交叉验证:若 IDW 与克里金圈出的面积差异超过约 30%,说明井网对羽状形态控制不足,应补充井位或加密采样,而不是强行选一个"好看"的结果。

二、浓度等值线绘制:从网格到闭合边界

等值线绘制通常分三步:确定插值网格分辨率、生成浓度场、提取指定浓度等值线并闭合。网格过粗(如 50 m)会抹平羽状前缘的指状突起;网格过细(如 1 m)则在井间制造虚假精度。一般取平均井间距的 1/5~1/10 作为网格步长较为合理。

边界闭合需要处理外推问题。等值线在数据边缘会开口,常见做法包括:以场地边界或地下水下游方向截断、按最外侧井浓度趋势外推、或明确标注"边界不确定区"。切忌把开口等值线直接连成闭合圈,那等于凭空扩大了污染范围。

下面这张图对比了同一组监测数据在三种处理下的边界差异,并标出了检出限与监管标准两条阈值线对最终圈定范围的影响。

插值方法与阈值对污染羽边界的影响(示例) IDW(p=2) 井间平滑过渡,边缘较"鼓" 普通克里金 向局部均值回归,范围偏小 最近邻(泰森多边形) 阶梯状边界,不做平滑假设 阈值选择对边界的影响(同一浓度场,不同判据) 低 高 浓度 → 检出限 边界最外 监管标准 合规边界 修复目标值 治理边界 浓度分布曲线

三、检出限与监管标准:两条阈值线决定"圈多大"

边界判据通常有三档,含义完全不同:

  • 检出限(LOD):方法能可靠识别的最低浓度,常取信噪比 3 倍对应值。以 LOD 圈边界,得到的是"污染可识别范围",通常最大,用于判断污染是否到达某敏感受体。
  • 监管标准:如《地下水质量标准》Ⅲ类限值或场地风险管控目标,是合规判据。以它圈边界,得到的是"超标范围",是修复范围与监测布点的法定依据。
  • 修复目标值:通常严于监管标准,用于圈定治理工程量。

三档阈值画出的边界可能相差数倍。例如某示例场地苯的监管标准为 10 μg/L,检出限 0.5 μg/L,若浓度场在 0.5~10 μg/L 之间有大片缓坡区,则 LOD 边界会比标准边界外扩 1~2 个井距。因此报告必须明确标注所用阈值,并说明边界外是否仍有低于标准但高于检出的残留。

常见错误:用检出限圈出的范围直接作为修复范围,导致工程量虚高;或只画监管标准边界而不说明外推区,使下游受体评估缺失。正确做法是同时给出"检出边界"与"超标边界"两套图,并注明插值方法、网格步长与外推处理方式。

四、边界判定的实操要点

第一,先看数据密度。井间距大于污染羽尺度时,任何插值都是猜测,应优先补井而非调参。第二,做敏感性分析:改变插值方法、幂指数或变异函数块金值,观察边界摆动幅度,把摆动范围作为"边界不确定带"画在图上。第三,边界必须与水文地质条件对照——污染羽前缘通常沿地下水流向拉长,若等值线形态与流向明显矛盾,应检查数据或插值设置。第四,所有边界图应附浓度单位、阈值来源、插值软件与版本,保证可复现。

归根结底,污染羽边界不是一条客观存在的线,而是"数据 + 方法 + 判据"共同定义的推断结果。把这三项写清楚,边界才具备可比较、可审查、可用于决策的价值。