重金属积累在土壤表层原因【关于城市表层土壤重金属污染的数学模型分析】

1、重金属积累在土壤表层原因【关于城市表层土壤重金属污染的数学模型分析】 摘?要 本文以某城区为例，建立了分析土壤重金属的空间分布及各种重金属污染的主要原因的数学模型，并求出了地质环境演变模式。首先，运用surfer8.0软件中的克里格（kriging）插值模块对土壤重金属浓度进行空间数据插值并应用单项污染指数法对各重金属元素的污染情况进行估计，得到了各元素在不同功能区的分布特征及各功能区的污染程度。然后，利用多元统计分析中的因子分析法对8种元素进行相关性分析及主成分提取，得到该区域内重金属污染的主要原因是工业污染、交通污染和生活污染。最后，在前面分析结果的基础上，通过收集该地区历年来土壤中重金属

2、的分布数据，建立土壤重金属含量的动态变化模型：qt=q0k?t+qk-z，运用该模型可以描述该城区地质环境的演变模式。关键词 土壤重金属；空间分布；污染源位置；克里格插值；因子分析；地质演变模式中图分类号 tn914 文献标识码 a 文章编号 1673-9671-(x)071-0188-02随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等

3、。现对某城市城区土壤地质环境进行调查，为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（010厘米深度）进行取样、编号，并用gps记录采样点的位置。由于在地表各重金属浓度的分布是相互影响的，并且受多种因素的多重影响，因此，我们应用因子分析法来研究重金属污染的主要原因。地质环境是指由岩石圈、水圈和大气圈组成的环境系统。各种元素在土壤中都是处于一个动态的循环过程。一是土壤本身含有一定的量，即土壤背景值，这一值是自然形成的；二是元素的输入是多途径、多层次的，如工业、生活污染等；三是输入的元素会有一部分随着河流冲刷、地表侵蚀、植物吸收等因素流失。为了研究城市地质环

4、境的演变模式，应该首先研究土壤中重金属含量的输入和输出，这与该地区的地表河流分布，地势分布，风向及降雨等因素有关，因此还需要测定各种因素的叠加所导致的元素输入及输出后的累积系数，这些可以通过分析该地区历年的重金属浓度分布数据来求出。结合各方面因素，我们建立了土壤重金属含量的动态变化模型：qt=q0k?t+qk-z1 土壤重金属空间分布及各功能区污染程度由于测量得到的只是有限个采样点的重金属元素浓度值，不足以涵盖整个城市的重金属含量情况，因此，首先需要建立模型对已知数据进行空间插值，得到该城市内重金属元素含量的总体情况，在此基础上进一步求解出各种重金属元素的空间分布并绘制空间分布图，从而可以分析

5、不同功能区内重金属的污染情况。步骤1：各功能区的地形特征运用kriging插值对数据进行处理，并绘制出该城市的海拔高度图及各区域的地形特征图，从而得出各功能区所处的海拔范围，即居民区主要分布在海拔为0 m20 m的区域，工业区、主干道区以及公园绿地区主要分布在海拔20 m80 m范围内，而山区主要分布在海拔高于80 m的范围内。步骤2：土壤重金属元素的描述性统计根据已知各重金属的浓度数据求出土壤重金属元素浓度的数字特征，通过与该地区重金属元素的背景值比较，我们得出，该城市土壤中重金属的含量均超过了背景值。步骤3：土壤重金属元素的空间分布1）克里格（kriging）插值模型原理。克里格（krig

6、ing）插值法是地统计学中应用广泛的一种空间插值方法，也是精确度最高的一种方法。kriging插值方法是在考虑了信息样品的形状、大小及其与待估块段相互间的空间分布位置等几何特征及品位的空间结构之后，为了达到线性、无偏和最小方差的估计，而对每一样品值分别赋予一定的权系数，最后进行加权平均来估计块段品位的方法。2）运用surfer8.0软件建立kriging模型并绘制重金属空间分布图。我们在确定了克里格插值模型后，利用surfer8.0软件绘制出土壤重金属元素的空间区域分布等值线图（见图1）（按顺序依次为as、cd、cr、cu、hg、ni、pb、zn）：通过观察各等值线图可以得出，as、cd及pb

7、元素分布较为广泛，其中，as和cd元素主要分布在主干道区， pb元素主要分布在工业区和主干道区，较为集中；cr元素主要分布在工业区和主干道区，cu元素主要分布在工业区，hg元素主要分布在工业区，ni和pb元素在工业区的浓度值较高，zn元素主要分布在工业区和主干道区。步骤4：不同区域重金属元素的污染程度利用单项污染指数法公式pi=ci/si分析各区域的污染程度。其中，pi为区域重金属i的单项污染指数；ci为重金属i含量的实测值；si为重金属i含量的评价起始值，以重金属i的背景值加上2倍标准偏差的结果表示。若pi1，则表示该区域受到污染。pi1的样点数在样点总数中所占的比例称为超标率，以超标率作为衡量污染程度的评价指标，并规定，超标率小于30%为轻微污染，大于30%且小于80%为中等污染，大于80%为重度污染。2 重金属污染的主要原因分析首先，对8种重金属元素进行相关分析，找到各因子间的相互关系，相关程度比较大的元素可能来自同一污染源；其次，要说明重金属污染的主要原因，需要找到影响8种重金属浓度的主要因素，可以应用多元地统计中的主成分分析法对各重金属元素进行主成分分析，选取前三个最主要的因子进行因子分析，从而得到影响各重金