徐 源， 师华定,*， 王 超， 费 杨， 于靖靖， 舒 密

1.中国环境科学研究院土壤与固体废物环境研究所， 北京 100012

2.生态环境部土壤与农业农村生态环境监管技术中心， 北京 100012

3.北京信息科技大学， 北京 100101

近年来，随着城市化进程的加快，重金属通过污水灌溉、大气干湿沉降、工业废渣等途径进入土壤[1-2]，重金属污染已成为当今土壤污染中污染面积最广、危害最大的环境问题之一[3-4]. Cd、Hg、As、Pb、Cr作为重金属“五毒”元素，在土壤中具有移动性大、毒性高、无法降解等特点，在生产活动中容易被作物吸收富集，不仅严重影响作物的产量和品质，还可以通过食物链在人体积累，危害人体健康[5-7]. 研究土壤“五毒”元素污染的来源以及污染源的主要影响区域，对维持区域环境健康具有重要意义[8-9].

目前，土壤重金属污染来源解析主要利用基于重金属总量的受体模型，常用的受体模型有化学质量平衡(CMB)模型、主成分分析/因子分析-多元线性回归(PCA/FA-MLR)模型、正定矩阵因子分解(PMF)模型、UNMIX模型等[10-11]. 其中，PMF模型具有适用性广、不需要测量源成分谱、能对因子贡献作非负约束等优点，得到广泛的应用. 例如，黄华斌等[12]利用PMF模型解析出研究区土壤重金属的来源；董騄睿等[13]通过PMF模型解析结果与Pb稳定同位素比值结果进行对比分析，指出PMF模型可以较好地应用于土壤重金属源解析研究.

然而，单一的PMF模型的解析结果往往比较笼统，对污染源的判别中主观性较强，且缺乏直观视觉效果，因此较多学者将其与地统计……