2011高教社杯全国大学生数学建模竞赛

承 诺 书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)：

日期： 年 月 日

赛区评阅编号（由赛区组委会评阅前进行编号）：

2011高教社杯全国大学生数学建模竞赛

编 号 专 用 页

评 阅 人 评 分 备 注 赛区评阅编号（由赛区组委会评阅前进行编号）：

赛区评阅记录（可供赛区评阅时使用）：

全国统一编号（由赛区组委会送交全国前编号）：

全国评阅编号（由全国组委会评阅前进行编号）：

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

摘要

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

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

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

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

?c?2c?2c?c?c的扩散和迁移，并由此建立污染物二维传播模型?Dx2?Dy2?Vx。随?Vy?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 P 意义 第i个监测点重金属浓度的测量值 重金属浓度的评价标准值 单因子污染指数 尼梅罗污染指数 x，y方向上的弥散系数 x，y方向上的对流速度 符号 c(x,y) F M 意义 某个时刻坐标（x，y）下重金属污染浓度 污染级别 重金属污染物的质量 土壤密度 分配系数 土壤含水率 ?b Kd Dx,Dy Vx,Vy?

四、问题分析

4.1问题1分析

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

4.2问题2分析

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

4.2问题3分析

3