1、1城市表层土壤重金属污染分析摘 要随着全球经济化的迅速发展，人类活动对城市地质环境影响日益显著，重金属污染物通过各种途径进入土壤，土壤重金属污染造成了生态系统严重破坏。模型一：将采样点的空间坐标与重金属污染物浓度综合起来，绘制出元素的空间分布图（三维图和等高图） ，综合使用空间分布图，可直观的看出各功能区的密集程度以及各种重金属污染物浓度。通过地质累积指数、内梅罗指数可求出每种元素在每个区域内的综合污染指数，利用各种重金属元素的综合污染指数，建立重金属污染物对环境污染程度的模糊综合评判矩阵，并运用贴近度算法，求出各功能区的相对污染程度，得到各区域污染程度大小序列：交通区工业区生活区公园绿地区山
2、区。模型二：利用因子分析法和主成分分析法，求出各金属元素的相关系数矩阵，求得两个主成分：对第一个主成分贡献较大的元素有 Hg、Cu，对第二主成分贡献较大的元素有 As 、Cr 、Ni。最终得出结论：该城区的污染主要是由人为活动影响的。交通污染和工业污染是主要原因。 根据各功能区各种重金属元素的综合污染指数，绘制出重金属元素污染指数直方图，可直观地观察出在不同区域内其主要污染元素，并寻找其源头及造成因素。 模型三：综合考虑重金属污染物的传播特征，根据重金属污染的不同传播特征，建立两个不同的污染源确定模型。一、根据污染物传播特征（高浓度向低浓度扩散） ，判断污染源可能在浓度较高的位置。利用各重金属
3、元素的污染指数生成归一化矩阵，按照各污染物的综合污染指数，将污染区域的位置分成三个模块。利用 fminsearch 函数最优解算法和初始模糊猜测值求目标函数最优解对应的空间坐标。从而得到模糊最优污染源位置。二、利用城区地貌特征和污染物传播特征（高海拔向低海拔扩散） ，判断原始污染源可能在海拔较高的位置，运用曲线拟合和微分方程的性质，可得到多个近似污染源；对重金属污染源位置分类讨论，得出原始污染源：（15600，15286，118） 。模型四：通过对模型三的综合评价，可看出该模型缺乏重金属污染源的分布随时间而产生的变化，及其对污染源位置确定的影响。采取按一定时间梯度在整个城区进行取样考察的方法，
4、对不同时间城市地貌进行插值逼近，求出城市地质环境中重金属污染物以及污染源在不同时间内的分布规律，进一步的探索人类活动影响下城市地质环境的演变模式。关键词：内梅罗指数、模糊综合评判、主成分分析2一、问题的重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为 1 类区、2 类区、5 类区，不同的区域环境受人类活动影响的程度不同。现对
5、某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距 1 公里左右的网格子区域，按照每平方公里 1 个采样点对表层土（010 厘米深度）进行取样、编号，并用 GPS 记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照 2 公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息，附件 2 列出了 8种主要重金属元素在采样点处的浓度，附件 3 列出了 8 种主要重金属元素的背景值。下面我们通过建立数学模型解决以下问题：(1) 给出 8 种主要重金属元素
6、在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的内在因素和人为因素。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析问题（ 3）所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集哪方面信息，以及如何利用这些信息，建立更优化的模型来确定污染源的位置。二、问题的分析由样本点的坐标，各种重金属元素的浓度，可以绘制出三维图形。通过内梅罗指数方程，确定各浓度污染物的权重，结合数据可确定各重金属污染物对环境影响的综合污染指数，并对内梅罗综合污染指数贴近度模糊评判，从而确定不同区域的重金属污染物对环境的污染程度。
7、对各种重金属元素进行因子分析和主成分分析，确定各城区重金属污染主要由某几种元素造成，具体分析该元素浓度大的原因。观察地貌图结合重金属污染物的传播特征，将其分为两方面进行考虑（即高浓度向低浓度扩散：求浓度较高点为污染源；高海拔向低海拔扩散：求海拔较高点） 。可以对数据进行拟合，或者运用 fminsearch 函数求最优解求出污染源的大概最优位置。研究城市地质环境的演变模式，需要结合时间问题来综合判断污染源的位置。一定程度上反映出人类活动对城市地质环境演变的影响。三、模型的假设与符号说明3.1、 模型的假设3（1）假设一：测量数据准确，不考虑测量仪器的误差及人为因素；（2）假设二: 采集数据在同一
8、时间点采集；（3）假设三：地质结构稳定，在短时间内不发生改变；（4）假设四：各种重金属元素的传播特征相同；（5）假设五：忽略生物因素对重金属元素传播造成的影响；3.2、 符号说明第i种元素在j功能区的地质累积指；ij第i种元素在j功能区的浓度；A第i种元素的背景浓度；jBE第 i 种元素在第 j 功能区的综合污染指数；ip地质累积指数算术平均值；()jave地质累积指数最大值;mxia归一化后各元素污染指数矩阵；y归一化后同样本点各元素污染指数和；x分别表示 X,Y,海拔高度；四、模型的建立与求解4.1.问题（1）:为了描绘出 8 种主要重金素元素在该城区的空间分布，分析了每个采样点各种重金属
9、元素浓度，并利用多元统计分析方法研究不同功能区与污染物浓度之间的相互联系,从而得出空间累积规律。将附件 1 与附件 2 整合到一起，为 shuju 表。对于该问题我们采用了画不同元素在整个区域的空间分布图，将采取样点的 x,y 作为 x 轴 y 轴坐标，z 轴为重金属元素在各个样点的浓度大小，绘制三维立体图；由于每种元素对于五个功能区的影响因素不同，我们可以根据 shuju 表，将五个区用不同的符号表示在图中，由每种元素在各功能区的浓度及地理位置，可以画出一个二维图。利用 MATLAB 软件（程序代码见附录表 1） ，绘制出 8 种主要重金属元素在该城区的空间分布图，如下：456根据 MATL
10、AB 软件画出的重金属元素在该城区内的空间分布图，可以分析出在某种重金属元素高对应的地理位置功能区分类的密集程度，判断出单元素对某个功能区的污染最大。我们将数据定量化分析，建立重金属对环境的模糊综合评判模型。 模型一：重金属对环境污染的模糊综合评判模型每种元素在每个功能区的综合污染指数：首先，我们为了研究沉积物及其他物质中重金属污染程度的定量指标，将 shuju表按功能区序号进行排序，令 Aj 为第 j 区样点各元素的浓度数据表，利用重金属污染物的地质累积指数(GeoaccumulationIndex 又称为 Muller 指数)，其表达公式如下：ij2=log1.5ijABE其中， 第i种元
11、素在j功能区的地质累积指数；ij第i种元素在j功能区的浓度；A元素i的背景浓度；jBE计算出 As、Cd、Cr、Cu、Hg、Ni、Pb、Zn 八种重金属元素分别对每个功能区的地质积累指数。其次，突出环境要素中浓度最大的污染物对环境质量的影响，采用地质累积指数与内梅罗指数相结合的综合指数法进行污染评价，得到新的综合指数，公式如下： 22max()ijaveijijp其中，第 i 种元素在第 j 功能区的综合污染指数；ijp地质累积指数算术平均值；()ave地质累积指数最大值;mxijMATLAB 运行（程序代码见附录 2） ，求得每种元素在每个功能区的综合污染指数 Cj,元素指数 As Cd C
12、r Cu Hg Ni Pb Zn C1 0.772 1.7287 2.831 2.6555 2.4005 0.58955 2.3712 3.4138C2 1.4348 1.8466 1.8517 5.0423 5.7374 0.83216 2.3298 2.8744C3 0.81694 0.83919 1.3788 1.3047 1.4612 1.4539 0.96855 0.91522C4 1.7542 2.2005 3.046 4.3847 5.8439 2.0895 1.4034 3.69417C5 0.79299 1.6977 0.75168 2.0384 3.3072 0.52323
13、 1.6217 2.6496采用贴近度的方法对五个功能区进行污染程度排序：以绝对污染程度（由同一元素五个功能区的污染程度的最大值作为该元素对环境的污染程度）和五个不同地区的不同元素的污染程度作为指标，建立模糊评价矩阵: 1.7542.03.465.02.8439.052.713.694818.7.7.696151.7542R.03.46.85.4392.08.43.69112072 采用贴近度 来计算五个功能区的相对污染程度：1()mikjijkrn生活区 工业区 山区 交通区 公园绿地区 0.6115 0.8515 0.3601 0.9721 0.4789 所得数据恰好与元素空间分布图的分布
14、状况相吻合，对八种元素在五个功能区的综合污染从强到弱排序：交通区、工业区、生活区、公园绿地区、山区。4.2.问题（2）重金属的主要原因：重金属以各种化学状态或化学形态存在于自然界中，随着大气、水、土壤、生物体的流动进入生态系统循环，然而重金属循环中不能被生物降解，可在土壤中富集，所以生活中的重金属含量会越来越多。重金属过量存在严重的损坏了人们的身体健康，因此分析某种重金属污染的主要原因十分重要。本文对重金属污染的主要原因从两个标准来进行评判。首先，对于整个城区而言，各种重金属对整个城区的危害程度不同。综合罗综合指数图表以及重金元素在各功能区的空间分布，对数据进行主成分分析，每种重金属元素的综合
15、污染指数与每个功能区具有线性关系，可求出每种元素在整个区的主成分贡献率：As Cd Cr Cu Hg Ni Pb Zn 0.0210 0.0258 0.0959 0.2551 0.3988 0.0453 0.0376 0.1206对各重金属元素在整个城区的贡献率进行排序：Hg、Cu、Zn、Cr、Ni、Pb、Cd、As。经贡献率排序可得，在整个城区内，重金属Hg、Cu 两种元素对环境的危害程度较大，而 Cd、As 元素对环境的危害程度较小。其次，对五个功能区而言，各种重金属对于每个功能区的危害程度不同。由模型一，可绘制出内梅罗综合指数图表：8根据图表，可以直观的观察出环境污染严重由哪种重金属元素
16、造成。在 C1（生活区）内，C r、Zn 含量较多，这些元素主要由生活中的化妆品、染发剂、燃煤、皮革、塑料门窗、蓄电池等生活用品产生。在 C2（工业区）内，As、Hg 含量较多，这些元素主要由工厂排放的废渣、废气、废水进入环境，也告诉我们在工业生产中应注意“三废”的先处理后排放产生。在 C3（山区）内，Ni 含量较多，这是因为含镍废水灌溉、植物生长吸收土壤中的镍、动植物的残体腐烂、岩石风化等原因造成山区重金属的积累。在 C4（交通区）内重金属元素含量比重都较大，这些元素的主要由机械中所含的重金属、交通工具尾气排放及轮胎磨损产生的大量含重金属的有害气体和粉尘等产生。在 C5（公园绿地区）内，Hg
17、、Zn 含量较多，这些元素主要由人类在公园留下的废弃物、大型娱乐机械产生。4.3.问题（3）模型三：归一化数值分析与污染源浓度目标函数求解模型污染源位置的确定与污染物的浓度有关，根据污染物的扩散方式可以知道，污染物可能是由高浓度向低浓度扩散，也可能是高地势向低地势扩散。因此污染源可能在（正常范围内）污染物浓度较高的地方，也可能在地势较高的地方。 按照同种浓度不同重金属对环境影响程度标准，将所给的数据进行归一化处理，可以得到矩阵 a（见附录） 。从得到的数据画直方图，可以看出图中有几处污染物的含量明显偏高，污染源的位置可能就在其附近，污染物含量偏高的点可能为特殊点（可能是人为原因导致，如有人向此
18、处堆积污染物） ，将此类点看为特殊点，对它们进行剔除。将污染源大概化为三个且分布在三个地理位置，然后按坐标位置将其划分为三个地理区域905101520251 19 37 55 73 91109127145163181199217235253271289307As ( g/g)Cd (ng/g)Cr ( g/g)Cu ( g/g)Hg (ng/g)Ni ( g/g)Pb ( g/g)Zn ( g/g)数据直方图对上图中的样本点 9,182,257 对应的地理坐标，由图 2 分析，对这三个点地理位置分别分析可知：样本点 9 附近的生活区高度密集，受人类活动影响因素较大，重金属元素 Cr、Zn 浓度
19、呈现较高的现象；样本点 182 附近交通区与生活区交叉分布，重金属元素 Hg、As 浓度呈现较高的现象；样本点 257 附近的工业区分布集中，重金属元素Hg、Cu 浓度呈现较高的现象；这些地区可能是由附近的交通要道、生活垃圾、工厂排放“三废”导致的重金属污染物富集，致使重金属污染物浓度较高。同时综合考虑到上述三处地方地势较低，也可能是山上的重金属元素通过雨水冲刷、岩石风化、大气沉降等原因导致的。因此，可以分两情况考虑。假如是上述第一种情况导致的，我们可根据函数最优规划对此地区的污染源进行较精确判断，跟据此地区的地貌特征和污染物综合指数的关系进行数据多元拟合建立方程。设 分别用数据进行拟合得到1
20、23(,)(,),xyzqxyzC(t)=f(x),C(t)=f(y),C(t)=f(z)。则可得到污染物浓度合坐标的对应方程 C（t）=P（q1,q2,q3）;然后，我们可根据无约束条件目标函数求解方法，用 matlab 中的 fminsearch 函数让目标函数在出事猜测点 8 点所对应的地理位置附近求污染物浓度最大点，即求 min-C(t)（也即目标函数在点 8 对应坐标的位置附近取最大值） ，当 C（t）最大时求出的点坐标就是理想的污染源所在点。10图 2上面我们大致确定了现在污染比较严重的区域的位置，通过对地貌的分析发现这些位置大部分在山脚下，而在高层的重金属会随着雨水的冲刷流向低处
21、，在相对低洼的地方富集，使得这些地方的重金属污染浓度较高，因此估计这些地区的污染可能是有山上的某个地方堆积了大量的重金属元素引起的，即可确定污染源的位置。利用归一化和数据，分别与 X,Y,海拔高度建立拟合模型；经分析得到的方程如下：y1=1.027e-012x3-2.71e-008x2 -0.0001086x+1.375y2=-7.061e-013x3+3.706e-008x2-0.000641x+4.305y3=-9.117e-007x3+0.0004927x2-0.07837x+4.185再利用微分方程的基本知识就可确定污染源大致在（15600，15286，118） 。4.4.问题（4）：
22、对模型三进行评价：优点：该模型综合考虑了城区的地貌、重金属污染物的传播特征与其分布状况，并对污染源的位置进行分类思考，利用本方案来确定污染源的位置更具科学性。 缺点：对城区土壤地质环境进行调查时，采样点的数据是在某个时刻测量出的重金属元素的浓度，该模型缺乏对时间因素的考虑。不能够很好的反映出城市地质环境的演变过程。为了对模型三进行优化，还应按一定梯度时间，在相同采样点进行测量各种重金属元素的浓度大小和对地表特征进行逼近。设 为扩散域D的地表曲面方程，其定义域为 .剖分 为由凸图形z=G (x,y) 2R构成的有限多子域 之和, 对所有的 、 ，有:ij，i如果，令 那么，必然有 1pk而且 其中 为 的边界, 为域 的近似 中子域 与 有公共点( l) 的子域边界， 为指标集。.
