3. 空间分布判别数据基础

土壤重金属空间分布判别,本质上是用有限点位的实测数据去推断连续空间上的污染格局。判别结果是否可信,首先取决于数据本身是否具有代表性、可比性和可追溯性。本节围绕“采样设计—样品采集与制备—含量测定—质量控制”这条数据链展开,说明每个环节的关键参数与操作要点,以及它们如何共同决定后续空间插值与分布制图的质量。

一、布点方法:让点位“代表”空间

布点方法决定样本在空间上的覆盖方式,直接影响能否捕捉污染的空间结构。常见方法有三类:

  • 系统网格布点:按固定间距划分网格,每格取一点(或格内随机取一点)。优点是覆盖均匀、便于计算采样密度、适合大范围区域调查与地统计插值;缺点是可能遗漏局部热点,网格间距应小于预期污染斑块尺度。
  • 随机与分层随机布点:简单随机布点无偏但可能局部聚集;分层随机先按土地利用类型、母质、距污染源距离等划分 strata,再在各层内随机布点,兼顾代表性与重点区域加密。
  • 判断性布点:依据污染源、风向、地形、水流方向等先验信息在疑似高值区布点,常用于污染源周边调查或加密验证,但不宜单独用于区域均值估计。

实际调查常采用“系统网格 + 重点区加密 + 判断性验证”的组合策略。例如某工业园区周边农田调查,以 500 m×500 m 网格覆盖全域,在主导风向下风向 1 km 范围内加密至 250 m×250 m,并在历史污灌区增设判断性点位。

二、采样密度:精度与成本的平衡

采样密度通常以“点/km²”或“网格间距”表示。密度过低会丢失空间变异信息,过高则成本陡增。确定密度时可参考:区域面积与变异程度、判别目标尺度(区域尺度还是地块尺度)、相关技术规范要求。一般区域环境调查常用 1~4 点/km²;农田地块尺度调查可达 8~16 点/km² 甚至更高;污染源周边加密区可提高到数十点/km²。一个实用原则是:网格间距应不大于预期污染斑块特征尺度的 1/2,否则插值结果易出现“伪平滑”。

三、采样深度:区分表层与剖面

重金属在土壤剖面中常呈表层富集。判别空间分布时,通常采集 0~20 cm 耕层或表层土壤;林地、草地可按 0~10 cm 或凋落物层下采集。若需评估垂向迁移,则按 0~20、20~40、40~60 cm 分层采样。采样深度必须在全区域统一,并在记录中明确,否则不同深度样品混用会引入系统性偏差。采样时去除表面砾石与植物残体,使用木铲、竹片或不锈钢铲(测定铬、镍时注意避免不锈钢工具污染)。

四、样品采集与制备:防止交叉污染

每个点位按“梅花形”或“蛇形”采集 5~8 个子样,混合后四分法缩分至约 1 kg,装入聚乙烯自封袋或棕色广口瓶,贴标签记录点位编号、经纬度、深度、土地利用、采样日期与采样人。样品在阴凉处自然风干,剔除石块与根系,用木棒碾碎后过 2 mm 尼龙筛;用于重金属全量测定的样品再经玛瑙研钵研磨至过 100 目(0.149 mm)筛。制备过程中应避免使用金属筛和金属研磨器,防止引入铁、铬、镍等污染。

五、测定方法:按元素与目的选择

常用方法各有适用场景:

方法适用元素/场景特点
原子吸收光谱(AAS)Cd、Pb、Cu、Zn 等单元素设备普及、成本低;火焰法适合 mg/kg 级以上,石墨炉法适合痕量 Cd
电感耦合等离子体质谱(ICP-MS)多元素同时测定,痕量至超痕量灵敏度高、线性范围宽,适合大批量多元素调查
X 射线荧光光谱(XRF)现场快速筛查、多元素半定量无损、快速,适合布点前预调查与热点初筛,精确定量需实验室方法验证

无论采用何种方法,样品前处理须与仪器匹配:AAS 与 ICP-MS 通常采用王水或 HNO₃-HF-HClO₄ 消解体系,XRF 则多用压片法或熔片法。测定时应同步做空白、平行样与标准物质对照。

六、数据质量控制:让数据“可判别”

质量控制贯穿全流程,核心指标包括:

  • 空白值:每批样品至少 2 个全程序空白,空白测定值应低于方法检出限或规定限值。
  • 平行样:每 20 个样品至少 1 组平行,相对偏差一般控制在 10%~20% 以内(痕量元素可放宽)。
  • 标准物质:使用有证标准物质(CRM)核查准确度,回收率一般要求在 90%~110%。
  • 加标回收:按批次做加标回收,回收率宜在 80%~120%。
  • 检出限与数据修约:低于检出限的数据应统一标注处理方式(如按 1/2 检出限参与统计),不可随意剔除。

此外,应建立点位—样品—检测数据的唯一编码对应关系,保留原始记录与仪器谱图,确保数据可追溯。只有通过质量控制的数据,才适合进入空间插值与分布判别环节。

土壤重金属空间分布判别数据链 ① 采样设计 系统网格 / 分层随机 密度:1~16 点/km² 深度:0~20 cm 统一 重点区加密 + 判断性验证 ② 采集与制备 5~8 子样混合缩分 风干 → 过 2 mm 筛 研磨过 100 目筛 避免金属工具污染 ③ 含量测定 AAS:单元素痕量 ICP-MS:多元素同时 XRF:现场快速筛查 消解 / 压片前处理匹配 ④ 质量控制 空白 / 平行样 标准物质回收率 加标回收 80~120% 数据可追溯编码 质控不合格 → 复采 / 复测 / 修正布点 合格数据集 点位坐标 + 元素含量 + 深度 + 质控记录 空间插值与分布判别 半变异函数 → 克里金 / IDW → 污染等级图 关键判断依据 网格间距 ≤ 斑块尺度 1/2 平行样相对偏差 ≤ 20% CRM 回收率 90~110% 常见数据缺陷 深度不统一 → 系统偏差 金属工具 → 交叉污染 无质控 → 插值结果失真

把上述四个环节串起来看,采样设计回答“在哪里采、采多密、采多深”,采集制备回答“样品是否干净、是否均质”,测定方法回答“用什么手段测、测到什么精度”,质量控制回答“数据能不能用”。四者缺一不可:布点再合理,若制备引入污染或质控缺失,空间插值仍会给出错误的热点位置。因此,在进入半变异函数拟合与克里金插值之前,应先完成数据集的完整性与合规性检查,确认点位坐标、深度、含量与质控记录一一对应,再开展后续空间分布判别。