4. 污染羽范围圈定方法与边界判定
污染羽范围圈定的本质,是把离散的监测井浓度数据转化为连续的空间浓度场,再按某一浓度阈值截取边界。因此边界位置并不唯一,它同时受三方面控制:插值方法如何"填补"井间空白、等值线如何绘制与闭合、以及用什么浓度值作为边界判据。理解这三者的耦合关系,才能解释为什么同一批数据在不同报告里会画出大小不同的污染羽。
一、插值方法:不同假设产生不同井间浓度
监测井是点状采样,井间浓度必须靠插值推断。常用方法各有隐含假设,直接影响污染羽的"胖瘦"。
反距离权重法(IDW)假设浓度随距离平滑衰减,幂指数 p 越大,高浓度点影响越局部。p=2 时井间过渡平缓,污染羽边缘较"鼓";p 增大到 4~5,等值线会更贴近高值井,羽状范围收缩。IDW 不产生超过实测最大值的极值,适合数据较密、梯度平缓的场地。
克里金法基于变异函数建模空间相关性,能给出估计方差。普通克里金在数据稀疏处会向局部均值回归,可能压低峰值、抬高边缘;指示克里金则直接对"是否超过阈值"建模,特别适合圈定超标范围。克里金的关键在于变异函数拟合是否合理,块金值过大意味着短距离随机性强,插值结果趋于平滑。
自然邻点法与最近邻法属于局部插值。最近邻法把每个网格点赋为最近井的值,结果呈阶梯状多边形(泰森多边形),边界生硬但对数据不做额外假设;自然邻点法过渡较自然,但外推能力弱。
实操中应至少用两种方法交叉验证:若 IDW 与克里金圈出的面积差异超过约 30%,说明井网对羽状形态控制不足,应补充井位或加密采样,而不是强行选一个"好看"的结果。
二、浓度等值线绘制:从网格到闭合边界
等值线绘制通常分三步:确定插值网格分辨率、生成浓度场、提取指定浓度等值线并闭合。网格过粗(如 50 m)会抹平羽状前缘的指状突起;网格过细(如 1 m)则在井间制造虚假精度。一般取平均井间距的 1/5~1/10 作为网格步长较为合理。
边界闭合需要处理外推问题。等值线在数据边缘会开口,常见做法包括:以场地边界或地下水下游方向截断、按最外侧井浓度趋势外推、或明确标注"边界不确定区"。切忌把开口等值线直接连成闭合圈,那等于凭空扩大了污染范围。
下面这张图对比了同一组监测数据在三种处理下的边界差异,并标出了检出限与监管标准两条阈值线对最终圈定范围的影响。
三、检出限与监管标准:两条阈值线决定"圈多大"
边界判据通常有三档,含义完全不同:
- 检出限(LOD):方法能可靠识别的最低浓度,常取信噪比 3 倍对应值。以 LOD 圈边界,得到的是"污染可识别范围",通常最大,用于判断污染是否到达某敏感受体。
- 监管标准:如《地下水质量标准》Ⅲ类限值或场地风险管控目标,是合规判据。以它圈边界,得到的是"超标范围",是修复范围与监测布点的法定依据。
- 修复目标值:通常严于监管标准,用于圈定治理工程量。
三档阈值画出的边界可能相差数倍。例如某示例场地苯的监管标准为 10 μg/L,检出限 0.5 μg/L,若浓度场在 0.5~10 μg/L 之间有大片缓坡区,则 LOD 边界会比标准边界外扩 1~2 个井距。因此报告必须明确标注所用阈值,并说明边界外是否仍有低于标准但高于检出的残留。
四、边界判定的实操要点
第一,先看数据密度。井间距大于污染羽尺度时,任何插值都是猜测,应优先补井而非调参。第二,做敏感性分析:改变插值方法、幂指数或变异函数块金值,观察边界摆动幅度,把摆动范围作为"边界不确定带"画在图上。第三,边界必须与水文地质条件对照——污染羽前缘通常沿地下水流向拉长,若等值线形态与流向明显矛盾,应检查数据或插值设置。第四,所有边界图应附浓度单位、阈值来源、插值软件与版本,保证可复现。
归根结底,污染羽边界不是一条客观存在的线,而是"数据 + 方法 + 判据"共同定义的推断结果。把这三项写清楚,边界才具备可比较、可审查、可用于决策的价值。