内容发布更新时间 : 2025/3/15 5:05:13星期一 下面是文章的全部内容请认真阅读。

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

摘要

随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出，城市土壤表层的重金属污染问题日益严重。本文采用克里格插值法绘制污染物空间分布图，内梅罗评价法对污染状况进行评价，相关分析和主成分分析去探究污染原因，建立二维传播模型确定污染源位置。

针对问题1，为了解决污染物空间分布问题，我们首先采用统计分析的方法大致分析了污染物在各个研究区的分布离散度等情况，然后根据克里格插值法用MAPGIS软件绘制了重金属污染物的空间等值线图并对其进行分析。为了解决污染程度问题，我们采用尼罗德评价法，分析出生活区Zn污染最为严重，整体一般清洁，工业区和交通区污染较为严重，山区和公园绿地大体环境清洁。

针对问题2，为分析重金属污染主要原因，我们首先对重金属污染进行Pearson相关分析，相关系数表明大部分重金属污染物之间都有较强的正相关关系，又因信息的重叠度较高，我们采用主成分分析法解析污染来源，发现Cd、Cr、Cu、Pb、Zn这几种重金属污染元素很可能来自农药施用，工业“三废”排放，As、Ni主要来自工业区的污染物排放。

针对问题3，为确定污染源位置，首先我们查阅了污染物的传播特征，根据污染物

?c?2c?2c?c?c?Vy的扩散和迁移，并由此建立污染物二维传播模型?Dx2?Dy2?Vx。随

?t?x?y?x?y后用1stopt进行多项式拟合确定传播系数，最后用MATLAB编程求解得: 元素 污染源个数

针对问题4，我们先根据模型分析其优缺点，然后考虑运移方程中的土壤密度，分配系数，土壤含水率，重金属污染的积累等条件，对模型进行优化并建立方程。

As 5 Cd 7 Cr 4 Cu 2 Hg 5 Ni 2 Pb 5 Zn 7

关键词：空间插值 尼梅罗指数评价 主成分分析 相关系数分析 二维传播微分方程

1

一、问题重述

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

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

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

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

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

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

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

(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。 (4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集的信息，以及有了这些信息，建立模型解决问题的方法。

二、模型假设

1、重金属在土壤中的传播方向仅仅是从高浓度向低浓度传播，不考虑其他的传播方向； 2、重金属在土壤中传播忽略高程在其中的作用，重金属污染只看做是迁移、扩散；

2

3、不考虑重金属在土壤中的自我降解能力；

4、该城市某个区域土壤是均匀的，其中的水的对流速度和扩散系数都是恒定的。 5、题目提供的相关数据皆是真实、可信，在没有特殊情况下测得的；

三、符号说明

符号 Ci Si Pi 意义 第i个监测点重金属浓度的测量值 重金属浓度的评价标准值 单因子污染指数 尼梅罗污染指数 x，y方向上的弥散系数 x，y方向上的对流速度 符号 c(x,y) F M 意义 某个时刻坐标（x，y）下重金属污染浓度 污染级别 重金属污染物的质量 土壤密度 分配系数 土壤含水率 P Dx,Dy?b Kd Vx,Vy?

四、问题分析

4.1问题1分析

问题1 题目要求我们给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。我们对题目给出的附件进行了数据的预处理和正态分布检验，消除了异常值带来的影响和比例效应。根据附件1、2给出的样本位置和样本点污染物浓度我们先对数据进行统计分析，然后用MAPGIS软件绘制空间分布等值线图更加直观的观察重金属元素的空间分布。我们采用内梅罗环境污染指数评价的方法结合土壤环境质量标准，做出了污染评价。

4.2问题2分析

问题2题目要求我们通过数据分析，说明重金属污染的主要原因。我们查阅资料，查出了我国境内8种重金属污染物的主要来源。对8种重金属污染物先采用Pearson相关分析法求出了它们的相关系数矩阵，求出的结果表明部分重金属之间表现了较强的正相关关系，信息存在一定的重叠，故我们后又采用主成分分析的方法，求出影响的综合指标，判断不同的重金属污染原因。

4.2问题3分析

3

问题3题目要求我们分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。我们查阅资料，分析了重金属污染物在土壤中的传播特征，据此建立重金属污染物的二维传播模型。在建立模型的过程中我们先根据重金属污染的空间分布图确定预备污染源位置，然后采用1stpot软件确定了参数扩散系数和流动系数，最后用MATLAB软件对方程进行求解得出污染源的位置，并对结果进行分析。

4.4问题4分析

问题四，要求我们分析我们建立模型的优缺点并且对模型进行优化，我们发现在上述问题进行求解时，我们并没有考虑土壤密度，分配系数，土壤含水率等条件，所以我们应该对系数进行优化，然后根据二维传播方程进行求解。

五、模型建立与求解

5.1问题一

5.1.1数据预处理

（1）对异常数据的处理

由于测量的值具有不稳定性，且所产生的异常值会对变异函数有很重要的影响，会影响变异函数模型的精度。所以我们将A?3s所表示的范围认为是正常测量值，而超过A?3s的测量值则用A?3s进行替代，将小于A?3s的测量值用A?3s替代。其中A表示某一重金属元素所有测得浓度的算术平均值，而s则表示标准差。所处理得到的数据如附录一所示。[2]

（2）正态分布检验

本文对八种重金属含量的分布进行了Klmogorov一smirnov检验,由检验结果发现As、Pb、Hg、Cd、Zn、Cr、Cu、Ni这八种重金属元素全部拒绝正态分布零假设，需要进行对数转换，以满足统计学对数据的要求，减少因异常值带来的影响。 5.1.2模型分析

（1）主要重金属元素的空间分布

问题1要求我们给出8种主要重金属元素在该城区的空间分布。根据附件1、2给出的污染物浓度和采样点坐标，对其先进行数据处理，然后采用统计分析的方法求出各个重金属污染物的离散程度，再采用克里格空间插值法结合MAPGIS软件绘制金属污染物空间等值线图并对求出的结果进行分析。

（2）分析不同区域重金属的污染程度——尼梅罗污染指数评价

在分析不同区域金属污染物浓度的时候，我们应该先筛选出每个区域所对应的站点，然后先用单因子评价对每个站点的每种污染物做评价，然后用尼梅罗指数法对该区域的所有站点进行分析。[2]

5.1.3模型的建立与求解

4