大学生数学建模 城市表层土壤重金属污染分析.doc

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

摘要

本文我们就城市土壤重金属污染程度的区域性和污染原因进行了定量的制图、说明和分析，并根据重金属污染传播特征中的地域性特点采用地统计学软件GS+对其进行空间地理位置拟合，结合重金属污染含量的地理分布和密集程度确定污染源。

为消除由于偶然因素造成的奇异选样数据对整体数据规律的影响，我们首先对所给的319个样本8种重金属含量浓度值的数据进行预处理、正态检验，对采样数据超出范围的异常数据进行区间极值赋值。采用地统计学软件GS+7.0对8种重金属含量浓度值在城市各个区域的分布规律做了空间3D图，我们可以更直观地看出每种重金属含量浓度在城市各个区域的分布规律。

我们以背景值为参考值采用单因子指数评价法和尼梅罗综合污染指数法评价土壤重金属各要素的污染程度，确定土壤尼梅罗污染指数的5个区间，再将经过单因子指数处理后的数据进行归类，进而确定城市每个区域的每种重金属污染的级别。之后，基于预处理后的采样数据，我们采用SPSS软件对其进行均值和方差等常规量统计，并对其变异系数和自相关系数进行了统计，以此来反映重金属污染的离散程度和各种重金属之间的相互作用规律，进而更科学地分析其污染原因。而后我们又根据重金属元素的传播特征，用地统计学软件GS+7.0对每种重金属污染元素的含量与其相对应的地理位置的坐标进行空间拟合，得出一个置信度较高的2D平面，再利用Photoshop软件对8种重金属元素含量分布2D图进行线性加深合成，得到一个8种重金属元素含量分布的2D合成图，再用GS+7.0确定污染源的位置在工业区,具体位置为X坐标范围，Y坐标范围，用假设检验对污染源所在地分析，在置信水平的情况下，否定原假设，选择备择假设，从而验证污染源所在地。

最后通过模型的优缺点，我们提供了相应的处理建议和方法，从而更精确获得数据，并能较准确模拟相应图形。

关键词：变异函数地统计软件GS+7.0尼梅罗综合污染指数

Persons相关系数

[正文页]

问题的提出与分析

1.1问题重述

随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。

按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。

现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其

附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。

现要求你们通过数学建模来完成以下任务：

(1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。

(2)通过数据分析，说明重金属污染的主要原因。

(3)分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。

(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？

1.2问题分析

由题目知我们可以分四步对城市表层土壤重金属污染进行分析。

针对问题（1）由附件1中的编号我们可以确定8种重金属元素浓度值及与其相对应的维度，经度和海拔值。再由功能区对每个编号进行分组处理，结合地统计学软件GS+7.0可以分别绘出8种重金属元素的浓度值在该城市不同区域内的不同地理位置的3D图，便于我们定性地对该城市不同区域重金属的污染程度进行定性分析。至于对该城区内不同区域重金属的污染程度的评价，我们可以考虑使用较为广泛的单因子指数评价法和尼梅罗综合污染指数法评价土壤重金属各要素的污染程度。首先为了消除由于数据选取的偶然性造成的对污染状况规律分析的影响，我们对数据进行预处理，考虑以采样数据超出范围的数据为异常数据，其中和为所处理数据的均值和标准差，当数据大于时，以采样点的最大值代替，当数据小于时，以采样点的最小值代替。对8种重金属元素浓度值依次进行过滤异常值处理，并且对于不同区的数据也要处理。以附件1中所给的背景值为参考值处理分区后的数据可以得到每种重金属元素浓度值在不同编号下的尼梅罗污染指数，由尼梅罗污染指数的评价标准区间，我们可以得到该城市每区每种重金属元素浓度的尼梅罗污染