重金属污染

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

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

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

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

我们参赛选择的题号是（从A/B/C/D中选择一项填写）：

我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：

参赛队员(打印并签名) ：1.

2.

3.

指导教师或指导教师组负责人(打印并签名)：

日期： 2012 年 8 月15 日赛区评阅编号（由赛区组委会评阅前进行编号）：

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

编号专用页

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

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

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

摘要

本文运用Matlab、SPSS及Excel软件，综合分析该城区各区重金属污染程度及其重金属污染的主要原因，并根据重金属的传播特征建立数学模型，寻找出污染源，同时更进一步修改模型以更好研究城市地质环境的演变模式。

针对问题一，利用Matlab作出各重金属浓度空间分布图，可看出重金属浓度大多呈高浓度向外逐步扩散现象。同时，利用单因素污染指数公式及污染等级规则，通过模糊综合评价模型，得出各类区域重金属污染程度如下：

针对问题二，利用SPSS软件对各重金属变量进行R型聚类分析，将重金属污染分为四类：Cr、Ni、Cu为一类，Cd、Pb、Zn为一类，而As、Hg各为一类。并计算出各区域各单金属污染指数进行对照分析，说明该城区重金属污染主要来源于工业活动及其生产品和交通主干道污染。

针对问题三，我们根据重金属物质在土壤中扩散的特点，利用有衰减的扩散模型，构造出浓度与坐标的关系，并利用matlab编程，多次试验，合理筛选数据，通过回归分析的方法，求解参数，从而得到不同重金属的不同污染源。

针对问题四，由于问题三求解的局限性，我们在模型三的基础上，增加了时间、水文、人类活动等因素，得出一个更为具体并符合实际的模型。有利于计算土壤中重金属的浓度。

关键词单因素污染指数模糊综合评价 R型聚类分析有衰减的扩散模型扩散规律

一、问题重述

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

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

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

现根据采样点的位置、海拔高度及其所属功能区，以及其8种主要重金属元素在采样点处的浓度和背景值等信息，建立数学模型，解决以下问题：

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

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

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

(4) 分析所建立模型的优缺点，并收集更多信息，建立数学模型，以更好地研究城市地质环境的演变模式。

二、问题分析

对于问题一，可通过Matlab画图函数分别直观描述出该城区各类重金属的空间分布图，包括各采样点的平面位置、海拔以及其浓度分布，并大致观测各类区重金属浓度的变化。同时，对该城区内不同区域重金属的污染程度进行分析，首先需引入单因素污染指数及其计算公式，并给出重金属元素的污染等级评判标准。基于此，本文建立因素集和评价集，以简单模糊统计方法构造隶属函数，并在各重金属元素本身对土壤地质的作用是一样的假设下，即不考虑各重金属本身的毒性，通过模糊综合评价模型对该城区各区域的各类重金属总污染程度进行综合评定。

对于问题二，需对数据进行分析，来说明重金属污染的主要原因。考虑到重金属种类有8种，比较多，可通过R型聚类分析对相关性较高的重金属元素聚分成几类，基于此法，可大致判断出哪几类重金属元素很可能源于同一性质或同一类污染物或者相似人为性活动，同时计算出各区域各重金属的污染指数和污染等级，由此进行对照分析，综合推测出该城区重金属污染的主要原因。

对于问题三，考虑到土壤中重金属的传播特征与分子运动扩散性质很相似，即一般以污染源为起点向周围辐射状传播，并且浓度逐渐地衰减，传播方向是从高浓度到低浓度。针对这些传播特点，建立衰减扩散偏微分方程。首先根据问题一所做的图进行初步判断，观察高浓度的大约范围，由此选择数据。然后利用有衰减的扩散数学模型，和回归分析，由此进行参数估计，得到污染源的坐标。

对于问题四，我们考虑到问题三所建立的有衰减的扩散数学模型具有一定的局限性，

因此在问题三的基础上，研究考虑更多能够影响土壤重金属浓度的因素，并将其加入问题三的问题中，由此得到一个更全面的新模型。

三、模型假设

（1）问题一中，各重金属元素本身对土壤地质的作用是一样的，即综合模糊评价模型中的各重金属元素所占的权重一样。

（2）问题三中，土壤中重金属的传播服从扩散定律，并不考虑其它大气或水流等因素的影响。

（3）问题三中，假设重金属在水平方向上的扩散速度是一样的，土壤表面顺滑，不存在重金属物在非污染源处堆积现象。

四、符号说明