A数学建模城市污染物.docx

资源描述

A数学建模城市污染物.docx

《A数学建模城市污染物.docx》由会员分享，可在线阅读，更多相关《A数学建模城市污染物.docx（45页珍藏版）》请在冰点文库上搜索。

A数学建模城市污染物.docx

A数学建模城市污染物

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

承诺书

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

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

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

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

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

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

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

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

参赛队员（打印并签名）：

指导教师或指导教师组负责人（打印并签名）：

日期：

年月日

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

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

编号专用页

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

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

评

阅

人

评

分

备

注

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

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

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

摘要

随着工农业的不断发展，我国土壤污染问题日益突出，进入土壤的持久性有机污染物、重金属污染等与日俱增，土壤环境面临着前所未有的压力。

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

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

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

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

工业区、交通区、生活区、公园绿地区、山区。

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

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

对于问题3，根据污染物扩散的特征，建立一维污染物扩散模型：

，以此推广建立污染物的二维空间扩散模型：

；再运用搜索算法在局部极值点附近搜索适量样本点；采用多元非线性拟合方法对二维空间扩散模型进行拟合，得到区域的最大值，即为污染源所在地,以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重金属的单因子指数

：

表示第i功能区的内梅罗综合指数

：

i从1到n分别表示n个一维观测点

：

表示污染源中心

：

表示污染源排污总量

：

表示第t年污染源的排污总量

：

表示观测点位置为x的污染物含量

：

表示二维污染源分布的中心位置

：

i从1到n分别表示n个二维观测点

四、模型的假设

假设1:

样本数据真实可靠，可以客观评价城市重金属的分布状况。

假设2:

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

假设3:

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

假设4:

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

假设5:

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

假设6:

污染源的的传播作用范围有限。

五、模型的建立与求解

数据的采集：

对采集点做散点图分布，观测采集点的地区分布，如下所示：

图1

采集点在地区空间上的分布相对均匀，基本覆盖城市的生活区，工业区，山区，交通区和公园绿地区，获得的信息能够较为全面的反映该城市各不同功能区的土壤污染状况。

5.1问题1模型的建立与求解

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

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

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

8种主要重金属元素在该城区的空间分布如下图所示（源代码见附录1）：

图1-1As污染分布图

图1-2Cd污染分布图

图1-3Cr污染分布图

图1-4Cu污染分布图

图1-5Hg污染分布图

图1-6Ni污染分布图

图1-7Pb污染分布图

图1-8Zn污染分布图

从各重金属空间分布图可大致看出：

1、山区和公园绿地区重金属污染程度较轻，生活区、交通区和工业区污染较为严重。

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

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

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

1）污染物评价模型的建立

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

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

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

1）单因子指数法：

（1.1）

为土壤中i种重金属的单项污染指数

为土壤中i种重金属的实测值

为i种重金属的土壤环境环境背景值

2）内梅罗综合污染指数：

（1.2）

j是城市功能区编号

i为重金属元素种类

n为重金属污染物种类

是j功能区单项污染因子

2）污染物评价模型的求解及分析

对数据预处理：

将附录1各观测点功能区分类B,C,D,E列复制黏贴到附录2中第J列，将数据复制黏贴到文件名为a.txt的txt文档中，保存路径如：

“C:

//a.txt”。

利用VisualC++6.0编程实现各功能区单因子指数和内梅罗综合指数的计算（源代码见附录2）。

单因子指数

功能区1

1.7418

2.2305

2.2264

3.7427

2.6583

1.4912

2.2292

3.4349

功能区2

2.0143

3.0239

1.7229

9.6618

18.353

1.6107

3.0013

4.0279

功能区3

1.1234

1.1717

1.2568

1.3119

1.1702

1.2564

1.1792

1.0622

功能区4

1.5856

2.7693

1.8727

4.7133

12.7664

1.4323

2.0495

3.5196

功能区5

1.7399

2.158

1.4076

2.2873

3.2855

1.2431

1.9583

2.2354

表1-1各功能区不同重金属的单因子指数

等级划分

单因子指数

≤1

1.0-2.0

2.0-3.0

〉3.0

污染等级

非污染

轻污染

中污染

重污染

表1-2单因子指数污染程度分级

注：

来源于中国绿色食品发展中心《绿色食品产地环境质量现状评价导则》2000年

‘

由数据发现，在工业区和交通区的重金属Hg含量严重超标，可推测该城市存在大型化工厂或电池制造厂等高排放量工厂。

内梅罗综合指数

功能1

功能2

功能3

功能4

功能5

Min

2.0904

9.2769

0.8242

6.4625

1.8689

Average

3.1706

13.533

1.2531

9.4264

2.7344

Max

6.6757

24.9653

2.8263

17.3667

5.0243

表1-3内梅罗综合指数

等级划分

综合污染指数

≤0.7

0.7-1.0

1.0-2.0

2.0-3.0

〉3.0

污染等级

安全

警戒级

轻污染

中污染

重污染

表1-4内梅罗综合指数污染程度分级

（注：

来源于中国绿色食品发展中心《绿色食品产地环境质量现状评价导则》2000年）

图1-9各不同功能区的内梅罗综合系数

5.2问题2分析

1）数据分析：

内梅罗综合污染指数受污染最严重的因子影响较大，一方面突出表现了高浓度的重金属元素Hg对工业区和交通区土壤环境质量的污染，另一方面也掩盖了其它重金属元素对该地区的影响。

因此，我们对剔除Hg后的数据再次进行计算，将两次结果进行对比，更加全面地反映出不同功能区的受污染情况。

内梅罗综合指数

功能区1

功能区2

功能区3

功能区4

功能区5

Min

2.0814

4.7344

0.8239

2.4765

1.4318

Average

3.1601

7.2859

1.2546

3.7937

2.0852

Max

6.662

15.7805

2.8331

8.0601

4.2902

表2-1排除Hg后的内梅罗综合指数

剔除Hg影响后，发现该城市工业区，交通区，生活区受其它重金属元素污染严重；公园绿地区污染情况也不容忽视，远离城市的山区受污染相对较轻。

2）通过数据特征，定性描述重金属污染原因

背景资料：

城市土壤重金属的来源存在两种途径：

一种是来自于成土母质；另一种是外源输入，重金属分布与城市化过程中人类的活动密切相关。

大量的城市土壤研究表明，城市重金属污染主要涉及As、Cd、Cr、Cu、Ni、Hg、Pb、Zn。

外源输入的重金属主要来源于：

1工业污染源。

主要是基础工业特别是采矿和冶炼业是向环境中释放重金属的主要污染源。

2生活污染源。

生活垃圾中含有的重金属能渗到土壤中;生活废水中的有机物可以吸收重金属,用这些水灌溉也可污染土壤。

3交通污染源。

机车的燃料中所含的重金属如对土壤的污染。

工业污染源

采矿冶金等行业三废排放：

Hg、Cr、Cd、As、Pb、Zn

交通污染源

汽车尾气排放：

Pb、Zn、Cd、Cr、Cu

生活污染源

生活废弃物：

Cr、Cu、Zn、Pb

表2-2各行业主要污染物类型

注:

中国环保信息网，城市土壤的重金属污染及生物治理，2010年[7]

对照表一，Cr,Ni,Cd,As含量较为接近相对背景值，由背景资料知道这几种元素受母岩控制较大，而人类活动对其影响较小。

而Hg,Cu,Zn,Pb等重金属含量在人类活动频繁的地势较为平坦的区域内含量明显增加，其中Hg，Cu含量明显幅度较大，工业区Hg平均值超标18.4倍,Cu超标9.7倍，交通区Hg平均值超标12.8倍，Cu超标4.7倍，说明该城市存在Hg，Cu排放严重超标的企业或工厂，并且其排放的污染物已经对该城市造成严重影响。

各重金属间的相关性分析：

各土壤重金属之间的依存关系，能帮助揭示重金属元素的污染来源，选用受人类活动影响较大的Hg,Cu,Pb和Zn为例。

分析元素之间的相关性，分析污染的来源。

由于Hg,Cu,Pb,Zn其数量值表现出等级特征，故采用spearman分析Hg,Cu,Pb,Zn之间的等级相关性。

使用spss中Correlate命令项spearman分别对功能区一，功能区二和功能区四的Hg，Cu，Pb,Zn做相关性分析，得到下表：