2011A数学建模城市污染物

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

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

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

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。

如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员（打印并签名）：1.2.3.指导教师或指导教师组负责人（打印并签名）：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘 要随着工农业的不断发展，我国土壤污染问题日益突出，进入土壤的持久性有机污染物、重金属污染等与日俱增，土壤环境面临着前所未有的压力。

因此，了解土壤污染物来源是切实有效地控制土壤污染，保障环境安全和农业可持续发展的重要前提。

本文针对城市表层土壤重金属污染问题，在合理的假设下，先绘制出重金属元素空间分布图；再根据单因子指数与内梅罗综合污染指数分析方法对各功能区污染程度分析，说明重金属污染的主要原因；接着，在此基础上建立二维污染物空间扩散模型，采用搜索算法在污染较为严重的区域搜索出适量的样本点，再采用多元非线性回归方法对二维扩散模型进行拟合，求解相应的参数；最后通过对第三问模型优缺点的分析，建立时空扩散模型，更加精确反映城市地质环境的演变情况。

对于问题1，根据附件数据，调用Matlab 软件中的surf 函数绘制出重金属元素空间分布图，采用三维图像表示城区地貌特征，用灰度描绘该地受污染程度，清晰反映出在不同地貌，不同功能区的重金属污染物空间分布特点。

选用单因子指数和内梅罗指数模型评价污染状况，运用Visual C++ 6.0编程计算出不同功能区单因子指数与内梅罗综合指数，由数据反映出工业区和交通区Hg 和Cu 严重超标，各功能区染程度由重到轻依次为：工业区、交通区、生活区、公园绿地区、山区。

对于问题2，对各功能区的污染物含量进行分析，分析推测该城市存在规模较大或者数量较多的Hg,Cu,Zn 严重超标的工厂。

为进一步揭示表层土壤重金属污染来源，对各个功能区的不同重金属污染元素进行相关性分析，运用spss 软件计算出各相关系数，揭示了不同功能区内土壤表层重金属污染的不同来源，分析可得，交通区污染主要来源于汽车尾气排排放，工业区污染主要来源于工厂的超标排放。

对于问题3，根据污染物扩散的特征，建立一维污染物扩散模型：22()0()x u Q x Q e σ--=，以此推广建立污染物的二维空间扩散模型：222000()()(,)exp[]x x y y Q x y Q σ-+-=-；再运用搜索算法在局部极值点附近搜索适量样本点；采用多元非线性拟合方法对二维空间扩散模型进行拟合，得到区域的最大值，即为污染源所在地,以Hg 的分布为例说明。

对于问题4，通过对上述模型优缺点的分析，进一步搜集相关数据，依据重金属传播所具有的富集性和不易分解性等特征，在空间分布模型的基础上引进时间变量，突出时间对于污染物扩散的影响，建立时空扩散模型，以更好地研究地质环境的演化模式。

本文描述了城市土壤重金属在城市不同功能区的含量分布、不同功能区各元素间分布格局的差异性、给出定量描述城区不同区域重金属的污染程度评价模型，通过对数据的分析以及数据间相关性的分析推测不同污染物的来源，并根据重金属的污染物的传播特征，建立确定污染源位置的模型，为合理规划和利用城市土壤、改造和提高城市环境质量、保障人类健康等提供重要依据。

关键词：城市土壤 重金属污染 单因子指数 内梅罗综合指数 扩散模型 非线性拟合一、问题重述1．1背景资料背景值：未受人类活动明显影响的土壤本身的化学元素组成和含量称为土壤背景值，或土壤环境背景值。

影响土壤背景值的因素很复杂。

包括数万年以来人类活动的综合影响，风化、淋溶、淀积等地球化学作用的影响，生物小循环的影响，母质成因、质地与有机物含量的影响等等。

因此，土壤背景值是—个范围值，而不是—个确定值。

1.2问题描述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。

对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。

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

现对某城市城区土壤地质环境进行调查。

为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。

应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。

另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。

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

1.3需要解决的问题：问题（1）：给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。

问题（2）：通过数据分析，说明重金属污染的主要原因。

问题（3）：分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。

问题（4）：分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析2.1 问题1：先对采集点做分布进行分析，若分布较不均匀，则需做空间插值处理，随后根据附录数据并运用matlab的surf函数进行绘图，得到用三维坐标反映四维信息图像，其中以白色到黑色的渐变过程描绘污染程度的由轻到重的变化。

再采用最为常用的土地污染评价方法单因子指数法和内梅罗综合污染指数法对五个不同功能区的污染状况进行评价。

2.2 问题2：由问题1得到各功能区不同种类的重金属含量的单因子指数，并由此分析其重金属污染物的主要类型，推测污染物来源的性质；因为对不同污染物的之间的相关性分析可以帮助分析污染源的性质，通过选择污染较为严重的几种元素进行相关性分析，放映了不同功能区污染来源的差异，并推测各功能区重金属污染的主要原因。

2.3 问题3：通过污染物的传播特点，建立二维污染物空间分布模型；考虑到污染源的位于污染严重的附近区域，仅在污染严重的区域搜索样本点，极大地简化算法，且能准确反映污染源的位置；采用多元非线性拟合方法求解二维高斯分布函数，确定局部极值所对应的坐标轴，即污染源所在地；对相关参数进行解释，说明重金属污染物的传播特点。

2.4问题4：对模型一和模型三存在的优缺点进行评价，通过对模型优缺点的分析，了解模型存在的问题，收集更多的信息，借用这些信息对模型进行改善，更好地研究城市地质环境的演变模式。

思路的流程图如下：三、 符号说明i P ：表示第i 重金属的单因子指数j p ：表示第i 功能区的内梅罗综合指数12,,,n q q q ：i 从1到n 分别表示n 个一维观测点Q ：表示污染源中心0Q ：表示污染源排污总量0()t Q ：表示第t 年污染源的排污总量()Q x ：表示观测点位置为x 的污染物含量00(,)x y ：表示二维污染源分布的中心位置1122(,),(,),,(,)n n q x y q x y q x y ：i 从1到n 分别表示n 个二维观测点四、 模型的假设假设1:样本数据真实可靠，可以客观评价城市重金属的分布状况。

假设2:各重金属污染物之间没有相互反应，传输中没有引起某种物质的消除和形成。

假设3:各重金属污染物在传播的过程中总含量没有衰减。

假设4:该城市重金属均由本城市市民生产和生活活动产生，无外来重金属的影响。

假设5:污染严重的局部地区仅存在唯一的污染源。

假设6:污染源的的传播作用范围有限。

五、 模型的建立与求解数据的采集：对采集点做散点图分布，观测采集点的地区分布，如下所示：图1采集点在地区空间上的分布相对均匀，基本覆盖城市的生活区，工业区，山区，交通区和公园绿地区，获得的信息能够较为全面的反映该城市各不同功能区的土壤污染状况。

5.1问题1模型的建立与求解为了分析8种重金属污染的分布情况，首先做出这8种重金属污染的分布图，以对重金属污染分布有总体认识；然后根据所给出的数据，用内梅罗综合污染指数法，定量的分析8种重金属的污染分布。

考虑到海拔因素对人类活动存在重大影响，为直观表现在不同地貌下，人类活动对城市不同功能区土壤重金属的影响。

以三维图像表示城市的地貌特征，用灰度描绘各种重金属污染物的污染严重程度（其中纯白代表无污染）。

8种主要重金属元素在该城区的空间分布如下图所示（源代码见附录1）：图1-1 As污染分布图图1-2 Cd污染分布图图1-3 Cr污染分布图图1-4 Cu污染分布图图1-5 Hg污染分布图图1-6 Ni污染分布图图1-7 Pb污染分布图图1-8 Zn污染分布图从各重金属空间分布图可大致看出：1、 山区和公园绿地区重金属污染程度较轻，生活区、交通区和工业区污染较为严重。

2、 Hg 分布较广，各功能区污染程度不一；各区As 、Cd 、Ni 含量普遍较低。

3、 重金属污染集中分布在海拔较低处，污染程度随着海拔的升高而减轻。

4、 工业区呈现混合污染态势。

1）污染物评价模型的建立土壤重金属评价方法包括生态危害指数法、地积累指数法、单因子指数法和内梅罗综合污染指数法、模糊贴近度方法、潜在生态危害指数法等。

因为内梅罗污染指数反映了各污染物对土壤的作用，并且特别考虑了污染最严重的因子，突出了高浓度污染物对土壤环境质量的影响，其环境质量指数在加权过程中避免了权系数中主观因素的影响。

故本次研究所采用目前在土壤重金属研究领域被广泛运用的单因子指数法和内梅罗综合指数法来进行重金属污染等级划分。

1）单因子指数法：i i iC P S =(1.1)i P 为土壤中i 种重金属的单项污染指数 i C 为土壤中i 种重金属的实测值 i S 为i 种重金属的土壤环境环境背景值2）内梅罗综合污染指数：2211(())()2nij iji j M ax P P np =+=∑ (1.2)j 是城市功能区编号 i 为重金属元素种类 n 为重金属污染物种类ij P 是j 功能区单项污染因子2）污染物评价模型的求解及分析对数据预处理：将附录1各观测点功能区分类B,C,D,E 列复制黏贴到附录2中第J 列，将数据复制黏贴到文件名为a.txt 的txt 文档中，保存路径如：“C://a.txt ”。