典型城市单元的土壤重金属溯源方法与实证研究.docxVIP

典型城市单元的土壤重金属溯源方法与实证研究

刘胜然,王铁宇,汤洁,孟晶,何博,赵慧,肖荣波

1吉林大学新能源与环境学院,长春1300122中国科学院生态环境研究中心,

北京1000853中国科学院大学,北京1000494广东工业大学环境科学与工程

学院,广州510006

土壤是生物圈中最活跃的部分,极易受到各种人类活动(工业、农业、交通)所产生的污染的影响。而在土壤化学组成中重金属含量对土壤的理化性质影响最大,当重金属在土壤中大量累积时会导致土壤性质发生变化,使土壤肥力下降,更为严重的会导致重金属在食物链中富集,最终影响人类的健康[1]。中国正处于快速城市化阶段,快速城市化不可避免的给土壤保护带来了极大的压力。根据全国土壤污染状况调查,土壤环境状况不容乐观。全国土壤总的点位超标率为16.1%,其中无机污染物超标点位数占全部超标点位的82.8%。土壤重金属的污染程度在不断加剧,但目前对土壤中重金属污染来源的认知尚不清晰。因此,探明土壤中重金属污染来源,对土壤重金属污染的防治、控制及保障农业生产安全

具有重要意义。因此,应尽早开展各尺度上的土壤重金属污染的源解析研究,以明确土壤重金属中不同污染来源及其相对贡献,为重金属污染的风险监控与防治管理提供决策依据。

污染源解析方法概述

污染源解析技术是一种针对污染区域的对污染来源进行定性及定量研究的方法。污染源解析技术最初主要应用于大气颗粒物,目前,已经在沉积物、水、土壤等多种介质中得到了广泛应用[2-5]。

污染源解析技术按研究对象的不同可以分为扩散模型与受体模型。其中扩散模型是以污染源为研究对象,其根据各污染源的污染物排放量、污染物的浓度、与环境受体空间距离,以及风向风速、温度等环境因素综合叠加各源对环境受体的影响程度,但由于排放清单的不确定性以及环境的多变性,因此模型预测的误差往往较大[6]。而受体模型的研究对象是基于受污染的环境介质,通过分析环境介质的污染特征及污染源的特征,从而定性识别污染源并定量确定每一类污染源的贡献率。由于这种方法不受环境因素多变的限制,并且不依赖于排放源的排放条件,也不用考虑污染物的迁移过程,因此得到了更为广泛的应用[7]。

受体模型按分析对象的不同可分为三类。第一类是以单一化学元素浓度为基础的单因子统计法,具有代表性的有元素浓度比值法、富集因子法。第二类是将多种化学元素浓度相结合的多参数统计法,具有代表性的有因子分析/多元线性回归分析法、主成分分析、正定矩阵因子法等。第三类是基于同位素测定的同位素比值分析法。

源解析方法体系及基本算法

现存的土壤重金属污染源解析的模型与方法有很多种,常用的方法包括同位素比值分析法、化学质量平衡法、主成分分析法、正定矩阵因子分析、地统计及污染空间分析法、富集因子法等。但由于各种方法都在存在各自的优势及局限性,而土壤环境污染问题往往是多种污染源共同导致的。多源并存普遍存在的情况给源解析带来诸多困难,单一的源解析方法往往难以保证污染源解析的准确性,因此需要多种方法联合使用,以提高源解析的准确性,增加源解析结果的可信

度。

本文应用正定矩阵因子分析法、同位素比值分析法及物质流分析法对特定区域重金属污染进行源解析及通量核算,建立了一套以三种方法为基础的源解析方法体系,为后续的案例实践及方法应用提供技术参考。具体步骤见图1。

图1结合3种方法的源解析方法体系Fig.1SourceappointmentsystemofcombiningthreemethodsPMF:正定矩阵因子分析法PositiveMatrixFactorization；IRA:同位素比值分析法IsotopeRatioAnalysis；MFA:物质流分析法MaterialFlowAnalysis

正定矩阵因子分析模型

正定矩阵因子分析法是美国环境保护署(USEPA)推荐使用的源解析方法之一,是一种基于因子分析原理的数据分析方法。其首先利用权重确定受体化学组分中

的误差,再通过最小二乘法迭代确定污染源及其贡献比例,是目前广泛应用于大气、土壤和沉积物污染物源解析的受体模型。其将原始浓度矩阵Eik分解成两个因子矩阵即Aij和Bjk以及一个残差矩阵εik,方程如下：

(1)

式中,Eik代表第i个样品第k个污染物的浓度,即浓度矩阵；Aij代表第i个样品在第j个源中贡献,即源分担率矩阵；Bjk为第k个污染物在第j个源中贡献浓度,即源成分谱矩阵；εik为随机误差。

PMF主要通过多次迭代计算分解原始矩阵,得到最优矩阵A和B,使目标函数达到最小值,目标函数Q为：

(2)

式中,σik表示Eik的不确定度。

当重金属浓度小于或等于相应的M