重金属污染溯源技术指南(试行)

重金属污染溯源技术指南(试行)1适用范围本指南规定了重金属污染溯源的工作流程、技术方法、质量控制与结果应用要求，适用于我国境内农用地、建设用地及流域水体沉积物、地下水等介质中镉（Cd）、铅（Pb）、汞（Hg）、砷（As）、铬（Cr）、铜（Cu）、锌（Zn）、镍（Ni）等重金属污染来源识别与解析工作。其他重金属污染溯源可参照执行。2术语与定义2.1重金属污染溯源通过资料收集、现场调查、样品检测、模型分析等技术手段，识别区域重金属污染来源，定量计算不同污染源贡献比例的过程。2.2源解析对区域或点位重金属污染的各类来源进行定性识别与定量贡献计算的技术环节。2.3源清单研究区域内各类潜在重金属污染源的位置、排放强度、排放方式、历史排放累积量等信息的汇编列表。3工作流程重金属污染溯源工作分为五个阶段，依次为：前期准备与源清单建立、点位布设与样品采集、样品检测与质量控制、源识别与源解析、结果验证与报告编制，具体流程见图1。4技术方法4.1前期准备与源清单建立4.1.1资料收集收集研究区域基础地理信息、地质背景、土地利用历史、污染源分布、环境质量监测历史数据五类信息，具体要求如下：基础地理与地质信息：比例尺不低于1:50000的地形图、地质图、土壤类型图，区域成土母质重金属背景值，区域地壳元素丰度数据。依据《中国土壤元素背景值》，我国土壤Cd背景值范围为0.0170.312mg/kg，中位值0.079mg/kg；Pb背景值范围10.656.3mg/kg，中位值23.4mg/kg，实际工作需采用本地实测背景值。土地利用与开发信息：近50年土地利用变更记录、历史工矿企业厂区分布、污染事故记录、农业投入品使用记录等。环境质量数据：近10年土壤、水、沉积物重金属监测数据，重点污染源监督性监测数据。4.1.2潜在污染源分类根据排放途径将重金属污染源分为四类：自然源：成土母质风化、区域地质高背景等原生来源。我国西南喀斯特地区、华南多金属成矿带普遍存在Cd、Pb等重金属原生高背景，部分区域成土母质Cd含量可达0.8mg/kg以上，超出农用地土壤污染风险筛选值。工业源：矿山开采与选矿、有色金属冶炼、化工、电镀、皮革鞣制、火力发电等工矿排放，主要通过废气沉降、废水排放、固体废物渗滤进入环境。单位产品产污系数参考：铅锌矿开采每吨矿石Cd排放强度为0.52.0g，有色金属冶炼每吨粗铅Cd排放强度为310g。农业源：磷肥与复混肥施用、畜禽粪便还田、农药使用、污灌等。我国磷肥Cd平均含量约为0.52.5mg/kg，长期施用磷肥区域农田土壤Cd年累积量可达0.0030.01mg/kg；规模化养殖场畜禽粪便Cu、Zn平均含量分别为50300mg/kg、1501000mg/kg。生活源与交通源：生活垃圾填埋渗滤、汽车尾气沉降、轮胎刹车片磨损等。无铅汽油推广后，交通源对Pb的贡献下降约70%，但轮胎磨损仍每年向环境排放约400吨Zn、50吨Cd进入周边土壤。

4.1.3源清单编制

对每类潜在污染源统计以下参数：污染源名称、地理位置、投产与停运时间、生产规模、历年重金属排放量、排放方式，计算1950年至今的累积排放量，建立网格化源清单，网格分辨率不低于1km×1km。4.2点位布设与样品采集4.2.1点位布设原则溯源针对性原则：覆盖潜在污染源影响范围，同时设置背景对照点，背景对照点需布置在未受人为活动影响、与研究区域成土母质一致的区域，数量不低于总点位的10%。网格布点与加密布点结合：流域尺度溯源采用网格布点，网格大小为1km×1km~5km×5km；工矿周边污染地块采用加密布点，以污染源为中心向外梯度布设，污染核心区点位间距不超过50m，影响区不超过200m。总点位数量满足解析精度要求，模型解析有效样本量不低于30个。4.2.2不同介质样品采集要求土壤样品：农用地采集020cm表层土壤，建设用地采集060cm表层土壤，污染地块需采集不同深度剖面样品，垂向每10cm分层采集。每个样品由5~10个子样混合而成，采样重量不低于2kg，采用聚乙烯塑料袋密封，避免金属工具交叉污染。水体与沉积物样品：水体采集表层下0.5m处水样，每个样品采集1L，加硝酸酸化至pH<2保存；表层沉积物采集0~10cm层，每个混合样品重量不低于1.5kg，密封冷藏保存。污染源排放样品：废气采集总悬浮颗粒物样品，采用大流量采样器过滤收集；废水采集瞬时混合样；固废、肥料、畜禽粪便等固体源采集具有代表性的样品，每个样品重量不低于1kg。4.3样品检测与质量控制4.3.1前处理要求土壤、沉积物固体样品经室温自然风干，去除碎石与植物残体，用玛瑙研磨机研磨，过100目尼龙筛，避免使用金属器具引入污染。水样经0.45μm滤膜过滤后酸化保存。4.3.2检测方法不同重金属检测方法选择符合下表要求，检测方法检出限需满足定量要求：重金属推荐检测方法方法检出限（mg/kg，土壤）Cd石墨炉原子吸收分光光度法（GFAAS）、电感耦合等离子体质谱法（ICP-MS）0.01Pb石墨炉原子吸收分光光度法、ICP-MS0.1As、Hg原子荧光光谱法（AFS）、ICP-MSAs：0.01；Hg：0.001Cr、Cu、Zn、Ni电感耦合等离子体发射光谱法（ICP-OES）、ICP-MSCr：0.4；Cu：0.3；Zn：1.0；Ni：0.24.3.3质量控制要求每批样品插入10%平行样，平行样相对偏差不超过±10%；每20个样品插入1个国家标准物质（GBW系列），测定值偏差在标准值允许范围内，回收率控制在90%~110%之间。实验器具全部采用聚乙烯或聚四氟乙烯材料，提前用10%硝酸浸泡24小时以上，去离子水冲洗干净，避免交叉污染。空白实验：每批样品做1个全程空白，空白浓度低于方法检出限视为合格。4.4源识别与源解析源解析分为定性识别与定量解析两个环节，根据研究尺度和精度要求选择适配技术方法。4.4.1定性源识别技术空间插值法：采用克里金插值、反距离加权插值绘制重金属含量空间分布栅格图，对比污染源位置，识别污染扩散梯度，判断主导源类型。若高值区沿工矿企业向外递减，或沿灌溉渠系分布，可初步判定工业源或农业污灌为主要来源。多元统计法：包括相关性分析、主成分分析（PCA）、聚类分析（CA）：将不同重金属含量数据标准化后，通过因子提取识别同源性，高载荷的重金属组合对应典型来源：如Cd-Pb-As组合通常对应有色金属矿冶源，Cu-Zn组合通常对应畜禽养殖与肥料施用源，Cr-Ni组合通常对应成土母质或电镀工业源。铅同位素示踪法：利用不同来源Pb同位素组成（²⁰⁸Pb/²⁰⁶Pb、²⁰⁷Pb/²⁰⁶Pb）的特征差异识别来源：自然源Pb的²⁰⁶Pb/²⁰⁷Pb比值通常>1.20，汽油燃烧源Pb比值约1.081.15，矿冶源Pb比值约1.151.18，可精准区分自然与人为污染源。磁性参数识别法：工矿与交通排放的颗粒物通常具有高磁性，通过测量土壤低频磁化率、频率磁化率，可快速筛选污染范围，识别工业与交通源影响。4.4.2定量源解析技术源清单法：基于污染源排放强度、污染物迁移扩散系数，计算各污染源对研究区域的累积贡献量，计算公式为：

式中，Ci为污染源i对目标区域的重金属贡献浓度，Ei为污染源i累积排放量，Fi为迁移到目标区域的比例，S为影响区域面积。

该方法适用于大尺度区域污染溯源，精度依赖源清单准确性，不确定度通常在20%~40%之间。

2.化学质量平衡模型（

式中，Cj为环境样品中元素j的总浓度，aij为源i中元素j的占比，Si为源i对受体的贡献浓度。

该方法要求明确列出所有排放源，源成分谱已知，适合污染源类型清晰的中小尺度区域解析，要求源成分谱之间共线性相关系数不超过0.8，否则解析结果偏差超过30%。

3.正定矩阵因子分解模型（PMF）：属于受体模型，不需要提前获取源成分谱，通过对受体样品数据的因子分解，约束因子得分非负，输出不同源的成分谱与贡献比例，是目前土壤重金属溯源应用最广泛的模型。模型运行要求：信号噪声比（S/N）<0.5的变量予以剔除，0.5<S/N<1的变量做弱信号处理，运行次数不低于20次，选择Q值最小的运算结果。PMF对样本量要求不低于30个，当污染源数量在3~6个之间时，解析结果相对偏差通常低于20%。

4.同位素源解析模型：采用Pb、Cd

式中，f为源1的贡献比例，Rs为样品同位素比值，R1、R2分别为两个端源的同位素比值。

该方法准确性最高，不确定度通常低于10%，适合小范围高精度溯源，受同位素测试成本限制，大尺度应用较少。

5.随机森林模型（RF4.5结果验证采用以下至少一种方法对解析结果进行验证：比值验证：定量解析的各源贡献浓度总和与样品实测总浓度的相对偏差控制在±15%以内，超过该范围需重新解析。同位素验证：采用同位素示踪结果对比模型解析结果，偏差超过20%需核查源分类与模型参数。独立样品验证：选取预留的10%~20%独立样品，对比预测浓度与实测浓度，决定系数R2≥0.7视为合格。

5不同类型区域溯源技术要点

5.1重点区分地质高背景与人为活动贡献：当土壤重金属含量高于区域背景值但低于筛选值时，首先开展成土母质贡献解析，若成土母质贡献比例超过60%，可判定为自然来源为主。优先识别工业源与农业源：工矿企业周边10km范围内农田，重点排查矿冶、工业排放的大气沉降与废水灌溉贡献；远郊平原农田，重点排查磷肥、畜禽粪便施用的累积贡献。统计数据显示我国农田Cd污染中，工业源贡献占比约42%，农业源贡献约38%，自然源贡献约20%。点位布设沿污染源-农田梯度设置，采集代表性肥料、灌溉水样品作为源样品，优先采用PMF模型结合铅同位素验证开展解析。5.2建设用地污染地块溯源针对地块内历史生产活动，逐一排查现有和历史工矿企业的排放，重点关注涉重金属生产车间、废渣堆放区等核心污染区。采用高密度加密布点，垂向分层采集土壤与地下水样品，明确污染扩散范围，结合污染源分布的空间叠置分析识别排放源。对责任主体不明确的污染地块，采用铅同位素指纹、重金属形态组成匹配污染源特征，追溯污染责任主体，解析结果可作为环境损害鉴定的参考依据。5.3流域重金属污染溯源按子流域单元划分溯源单元，沿干流、支流梯度布设沉积物与水体点位，同步采集上游背景点、污染源入汇口下游点位样品。优先采用源清单法结合PMF模型，计算不同行政区、不同污染源对流域断面的污染贡献，输出流域污染源排放贡献清单。我国流域重金属污染中，矿冶源贡献占比可达55%以上，是流域重金属污染的首要来源。水文条件变化大的流域，需分别采集丰水期、枯水期样品，分析不同水文条件下污染源贡献的变化。6质量控制与质量保证全流程质量控制：从资料收集、点位布设、采样、前处理、检测、解析全环节建立质量控制记录，所有原始记录留存时间不低于10年。不确定度评估：报告中明确溯源结果的不确定度范围，基于模型的定量解析不确定度不得超过30%，同位素解析不确定度不得超过15%。数据审核：检测数据需经三级审核，源解析结果需由两名以上相关领域技术人员校核，确保源分类合理、模型参数选择正确、结