城市表层土壤重金属污染的空间分布特征分析

城市表层土壤重金属污染的空间分布特征分析

摘要：分析了某城市城区表层土壤中的 as、cd 等8种重金属在生活区、工业区、山区、主干道路区及公园绿地区的含量水平，得出了不同区域重金属的污染程度；运用污染负荷指数法对影响土壤各重金属主要因子进行分析，确定不同区域重金属污染的主要特征；建立重金属污染物的传播模型，运用kriging插值法对重金属含量进行最优无偏估计插值，对重金属污染的空间分布进行分析，揭示了城市表层土壤中重金属含量的空间分布特征。

关键词：城市城区；表层土壤；重金属污染；空间分布特征

中图分类号：x530.2 文献标识码：a 文章编号：0439-8114（2013）06-1287-05

随着城市经济的快速发展和城市人口的不断增加，大量工业“三废”、城市生活垃圾和污泥等污染物的排放与不恰当处置使重金属在土壤中不断积累，加重了土壤重金属的污染负荷，导致我国城市表层土壤的重金属污染日趋严重。而城市土壤重金属污染是能有效反映城市环境污染状况的重要指标之一。因此，对城市土壤环境异常的查证并应用查证数据开展城市环境质量评价、研究人类活动影响城市土壤环境的演变模式日益成为人们关注的焦点。

1 数据来源与研究方法

以某城市城区为研究区，将其划分为间距1 km左右的网格子区域，按照每平方公里1个采样点对表层土壤（0～10 cm土层）进行取样、编号，并用gps记录采样点的位置。应用专门仪器测试分析，

获得了每个样本所含的多种重金属元素的浓度数据。另外，按照2 km 的间距在那些远离人群及工业活动的自然区取样，将分析数据作为该城区表层土壤中重金属元素的背景值。研究以2011年高教社杯全国大学生数学建模竞赛a题[1]所列的数据为数据来源，文献[1]列出了采样点的位置、海拔高度及其所属功能区的信息、8种主要重金属元素在采样点处的浓度和8种主要重金属元素的背景值。按照功能划分，现代城市整个城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、3类区、4类区、5类区，不同的区域环境受人类活动影响的程度不同。作者对城市表层土壤中的 as、cd、cr、cu、hg、ni、pb、zn 8种重金属元素在不同区域的含量水平进行了研究。采用统计学原理和污染负荷指数法[2，3]，根据国家污染负荷指数与污染程度关系分级标准进行评价，得出该城市不同区域的土壤重金属污染程度。对采样点进行最优估计，借助空间变异分析理论与方法[4]，确定不同区域重金属污染的主要污染源：工业区重金属污染主要来源于工业生产产生的工业垃圾；生活区的重金属污染主要是由于人们造成的生活垃圾形成的；在交通要道主要是汽车尾气的排放；山区人口分布比较分散，自然绿地的覆盖面积广，环境自净能力比较强，所以土壤重金属含量相对较低。建立重金属污染物的传播模型，运用kriging插值法[5]对重金属含量进行最优无偏估计插值，对重金属污染的空间分布进行检验，从而揭示城市表层土壤中重金属含量的空间分布特征。

2 结果与分析

2.1 重金属污染的空间分布

根据统计学的基本原理，采用excel 2003进行数据分析，得出所有数据的平均值、最小值、最大值、标准差、变异系数、偏度、峰度等统计信息，同时采用区域识别法[6]对原始数据进行异常值检验。对于已处理的数据，采用spss中的k-s方法[7]进行数据的正态分布检验，然后采用surfer 8.0软件绘出了8种重金属污染在该区的空间分布图，图中纵横坐标表示采样点的位置，图例表示采样点所在的功能区（图1a至图1h）。

2.2 不同区域重金属污染程度分析

土壤重金属污染程度的评价分别采用研究区土壤重金属含量的自然背景值和国家一级标准值（保护区域自然生态、维持自然背景的土壤环境质量的限值）为参考标准，当评价标准为背景值时，相应的评价结果称之为元素富集[8，9]，当以国家一级标准为评价标准时，称之为元素污染。

2.2.1 研究区重金属污染程度的污染负荷指数法分析污染负荷指数法是在从事重金属污染水平的分级研究中提出来的一种评价方法，被广泛应用于土壤重金属污染的评价。污染负荷指数法优点是该指数由评价区域所包含的多种重金属成分共同构成，并使用了求积的统计方法，通过这种方法对整个区域各个点位各种重金属进行了定量评价，并对各点的污染程度进行分级，能直观反映对环境污染最严重的元素和各种元素对环境污染的贡献程度以及重金属

在时间、空间上的变化趋势[8]。重金属污染负荷指数与污染程度关系见表1。

根据重金属污染负荷指数与污染程度关系，不同的污染负荷指数可将污染程度分为无污染、中等污染、强污染和极强污染4个等级，分别用0、ⅰ、ⅱ、ⅲ来表示[10，11]。根据污染负荷指数法计算得到研究区不同区域的污染负荷指数，参照污染负荷指数与污染程度关系分级标准，得到研究区不同区域的重金属污染程度见表2。

2.2.2 研究区重金属单元素污染程度的单项污染指数法分析单项污染指数法是以土壤单项污染物的实测值与评价标准的比值，表示某一污染物影响下的环境污染指数，可以反映各污染物的污染程度。这里污染分级标准采用国家一级标准为参考标准。选用国家土壤环境质量标准中的一、二、三级标准分别作为土壤重金属污染评价的污染积累起始值、中度污染起始值和重度污染起始值，把土壤重金属污染分为4级（表3）。

利用单项污染指数计算公式[12]求出五类区中采样点的各重金

属元素的单项污染指数，并求其平均值，将平均值作为此元素在该类区的单项污染指数，按单项污染分级标准对重金属单元素的污染程度进行分级。根据上述方法，得到研究区表层土壤中各重金属元素的单项污染指数见表4。从全区各元素的平均污染指数可以看出研究区表层土壤受as、cd、cr和 ni污染相对较轻，污染程度低；其中cu、hg、pb和zn的污染指数较大，表明研究区受这4种重金属元素的轻度污染。

2.3 重金属污染的传播模型

地统计学是建立在区域化变量、平稳假设和内蕴假设等概念基础上，以区域化变量理论为基础，以半方差函数为主要工具，研究那些在空间分布上既有随机性又有结构性，或空间相关性和依赖性的自然现象的科学。由于重金属污染的传播途径受很多因素的影响，比如土壤养分、物理特性、土壤重金属等性质的变异以及时间上的变异，基于地统计学的空间变异理论[13]，将地统计学方法与gis 相结合来研究变量的空间变异特征。

半方差函数又称半变异函数，能同时描述区域化变量的随机性和结构性[14，15]。利用半方差分析时将这8种重金属元素的浓度进行对数转换，连同采样点的地理坐标输入地统计软件gs+，拟合半方差函数，并选择最佳拟合模型及其参数。

球形模型：

上式中，c0为块金值；c0+c为基台值。以as为例，在gs+软件中拟合半方差函数，得到4种拟合模型如图2，该4种拟合模型的参数见表5。在半方差函数中，决定系数（r2）较大、残差（rss）较小回归效果越好，由表5可知，as元素的最佳拟合模型为线性模型。对8种元素用相同的方法确定其最佳拟合模型，运用gs+7.0

地统计软件系统完成半方差函数理论模型的自动拟合、最优选择、参数计算及半方差图绘制，结果如表6及图2（a～d）所示。

将基于地统计学的空间变异理论与gis的空间分析功能相结合分析研究区表层土壤重金属含量的空间变异特征和空间分布特征。按