城市土壤表层重金属污染分析.doc

城市表层土壤重金属污染分析 摘 要：本文主要分析了城市表层土壤重金属污染情况。基于所给数据，画出重金属元素在该城区空间分布。通过熵值法、内梅罗污染指数法，得到城区不同功能区的污染程度。并利用微分方程的数学分析方法，确定污染源的位置。最后利用灰色模型理论，分析该地区的地质演变模式。 问题一 首先对所给数据进行整理，筛选出不同功能区的各种重金属元素的通过Matlab编程，并通过等高线法表示海拔高度。并通过熵值模型、内梅罗模型计算得到，工业区污染最为严重，交通区次之，接下来依次是生活区、公园绿地区，山区的污染程度最小。 问题二 由分布图可以知道，工业区的主要污染金属为Hg、Cd交通区内的主要污染元素为Pb，生活区内的主要污染元素是As，公园绿地区和山区的污染较小元素种类 As (μg/g) Cd (ng/g) Cr (μg/g) Cu (μg/g) Hg (ng/g) Ni (μg/g) Pb (μg/g) Zn (μg/g) max 20 1000 300 250 1900 45 180 900 min 0 0 0 0 0 0 0 0 其中：为某对象在指标下的平均值，为第个指标的平均值，为标准差 计算第项指标下第功能区指标值的比重: 5.2.3 计算信息熵值及权重 式中: 计算信息效用值: 5.2.4 计算评价指标权重 指被评价的对象个数，这里就是5 计算样本评价值如下： 最终结果如表2所示： 表2 不同功能区熵值法评价值 功能区 1 2 3 4 5 评价值 7.4728 33.6662 0.7793 21.5877 6.0467 式中，p综为某样点( 或某区域) 所有元素的内梅罗综合污染指数; pmax为监测点土壤重金属最大单项污染指数；pave为监测点所有重金属单项污染指数平均值。 得到的不同功能区内梅罗综合污染指数如表3所示（具体过程见附件三）。 表3 各个功能区综合污染指数 功能区 1 2 3 4 5 内梅罗综合污染指数 3.1706 13.5330 1.2531 9.4264 2.7344 重金属污染主要原因 通过数据以及问题中各个金属浓度的分布图，可以看出在工业区的主要污染金属为Hg、Cd，可能原因是工业区内有大量的含汞废水、冶炼排放物； 在交通区内的主要污染元素为Pb，因为汽车排放的废气中含有大量的铅元素； 生活区内的主要污染元素是As，可能是砷被大量用作杀虫剂、杀菌剂、杀鼠剂和除草剂； 公园绿地区和山区的污染较小。处，则此扩散问题满足Cauchy问题[3] （1） （2） 对（1）（2）进行傅立叶变换，且令 ， 由于 故得常微分方程Cauchy问题 得唯一解 .　 　（3） 对（3）求逆变换，由于 ， ， 故得 （4） 对（4）式两边同时取对数得： （5） 令： 则式（5）变形为： （6） 下面主要用Matlab软件对式（6）进行非线性拟合（见附件四），拟合的过程是：首先根据问题一中画出的金属浓度分布图，大致地判断污染源的区间并找出这个区间所有取样点的坐标。这样就有很多组W、x、y、z的值，用Matlab进行非线性规划，即可求出各个参数的值。 现在以As为例来求其污染源。 由As的浓度分布图，污染源所在的区间为，。找出这个区间内所有取样点得坐标，见表4。 表4 所取取样点的坐标 x/m y/m z/m 15467 8658 17 17044 10691 93 17087 11933 43 18413 11721 88 19007 11488 84 18738 10921 53 17814 10707 64 18134 10046 41 17198 9810 37 17144 9081 20 18393 9183 26 19767 8810 46 20554 11228 43 20101 10774 40 20215 9951 31 16301 8299 24 17904 8287 25 19072 8519 36 20282 8590 57 15140 11101 30 15198 10100 28 15248 9106 16 16428 9069 20 16289 10072 43 16267 11058 60 用Matlab程序[见附件 ]对（6）进行非线性拟合，结果如下： 置信度为95%的置信区间，残差范围是[-0.263,0.335]，得到上述各参数的拟合值见表5