重金属污染源地数学建模.docx

重金属污染源地数学建模.docx

《重金属污染源地数学建模.docx》由会员分享，可在线阅读，更多相关《重金属污染源地数学建模.docx（12页珍藏版）》请在冰豆网上搜索。

重金属污染源地数学建模

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

承诺书

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

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

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

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

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

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

A

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

20001032

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

中国人民解放军国防科学技术大学

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

1.腾骥

2.袁雪强

3.泉朴

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

文强

日期：

2011年9月9日

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

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

编号专用页

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

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

评

阅

人

评

分

备

注

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

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

基于MCMC的城市表层土壤重金属污染分析

摘要：

本文针对某城市城区表层土壤重金属污染分析的问题，从各分区重金属元素污染的空间分布、污染程度、污染原因及传播途径等多方面进行了全面的统计与研究，同时采用贝叶斯推理和马尔科夫链蒙特卡洛法（MCMC），利用MATLAB进行编程，对所设计的模型进行求解，得到了各重金属污染源的地理位置坐标，最后对城市地质环境演变模式的问题进行了探究，得出了初步的结论。

首先，根据所给数据信息，利用MATLAB及EXCEL作图，得到了该地区的地势图及各重金属污染程度的分布图。

对于问题1，通过统计各重金属的重污染点在各区域的分布和计算不同区域各中金属污染程度之间的相关系数，得到了各重金属的大致空间分布，并采用改进的梅罗指数法，得出不同区域重金属的污染程度。

对于问题2，在问题1所得数据及结论的基础上，结合实际情况，对各重金属污染的主要原因进行了较为全面的分析。

对于问题3，通过对污染的传播特征的研究，基于城市的地质特性，得到空气传播为污染主要传播方式的结论，并由此建立关于气体传播的概率模型，采用贝叶斯推理和马尔科夫链蒙特卡洛法，利用MATLAB进行50000次抽样，得到收敛的坐标值，从而计算出各重金属污染源的地理位置坐标分别为：

As（2742,7294）、Cd（2382,3693）、Cr（3470,2309）、Cu（2707,2295）、Hg（2708,2294）、Ni（1647,2729）、Pb（2882,3617）、Zn（4152,2299）。

对于问题4，以题中城市土壤重金属污染的演变规律为例进行了一定深度的探究，并对所需信息及模型设计方案进行了初步的讨论。

关键词：

梅罗指数法相关系数贝叶斯推理MCMC法MH算法

一、问题的提出与分析

（一）背景介绍

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

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

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

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

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

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

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

（二）要解决的问题

现给出采样点的位置、海拔高度及其所属功能区等信息，8种主要重金属元素在采样点处的浓度及8种主要重金属元素的背景值，解答以下问题：

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

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

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

4.分析所建立模型的优缺点，并说明为更好地研究城市地质环境的演变模式，还应收集什么信息，有了这些信息，如何建立模型解决问题。

（三）问题分析

整个问题可以分为三类分别进行求解：

问题1和问题2

可将其归结为对所给数据的统计与分析，可以通过作图或列表等直观的形式针对不同金属和区域进行相关统计，再对统计的结果加以比对与分析，即可得出结论。

问题3

确定污染源的位置实质就是要找到使周围污染达到最大的点，这样可以转化为概率模型，在统计数据的基础上，运用贝叶斯推理马尔科夫链蒙特卡洛法可加以计算和确定。

问题4

开放类问题，可根据现有资料进行回答。

二、模型假设

1.题目中所给出的样本均具有普遍性及代表性，而非所测地区中的特殊情况。

2.该地区的空气流动没有显著的方向性特征。

3.污染物浓度在取样前后较长一段时间不会有显著的变化。

三、符号说明

:

土壤污染元素综合污染指数

:

土壤中污染元素i的实测值

:

土壤中污染元素i的评价标准值

:

任一闭合曲面

:

指定所围的区域

:

通过S流入Ω的质量

:

的质量

:

Ω由于浓度的变化引起的质量增加

:

衰减系数

:

贝叶斯公式模型参数

m:

模型参数的个数

y:

贝叶斯公式观测数据

：

马尔可夫蒙特链卡洛法中的分布

：

接受概率

:

目标概率密度函数

s:

距离污染源的距离

p：

沉降系数

h:

海拔高度

：

临界海拔

：

重金属元素吸附系数

四、模型建立与求解

（一）空间分布的特征

根据所给数据，对采样点按区域进行重新的排序，由此作出各采样点的地势分布图（图一），又由各种金属在各区域的浓度值，作出金属污染程度在各区域的分布图（以Pb（铅）为例，见图二，其它金属元素分布图见附件），由此分析可得到各金属元素在该城区的空间分布。

以Pb为例进行分析，经统计，Pb的重度污染区域主要集中于图的左下方，且Pb的重度污染点在工业区占62.8%，另有20.9%重度污染点则分布在左下方的交通区中，由此可得Pb的污染区域主要分布于工业区和左下方的部分交通区，对于其它金属元素，可以用相同的方法，还可以通过计算不同区域各金属元素污染程度的相关系数，找出各金属元素污染之间的联系，这样可以有效提高分析效率和准确性。

例如对表一中工业区的各金属元素污染程度相关系数（其它区域见附件）的比较可以得出，Cr与Cu、Cr与Hg、Cu与Hg、Cd与Pb、Cd与Zn、Pb与Zn相关系数较大，可以认为它们具有相同或相近的空间分布。

表一工业区各金属元素污染程度相关系数

Cd

Cr

Cu

Hg

Ni

Pb

Zn

As

0.328601

0.379628

0.152932

0.181295

0.689725

0.394663

0.517731

Cd

-

0.540985

0.566458

0.533171

0.488666

0.829216

0.753606

Cr

-

-

0.919726

0.902139

0.69825

0.675399

0.69507

Cu

-

-

-

0.98348

0.50278

0.66966

0.621675

Hg

-

-

-

-

0.479139

0.612465

0.590405

Ni

-

-

-

-

-

0.57762

0.634098

Pb

-

-

-

-

-

-

0.738756

根据以上数据分析，可以得到各重金属元素大致的空间分布，分别为：

1.Hg、Cd、Pb、Zn、Cr均主要分布于工业区和左下方的部分交通区；

2.As、Cu主要分布于工业区；

3.Ni的分布比较分散，没有明显的富集现象。

（二）重金属污染程度的划分

为了全面、综合地反映污染物的整体污染水平，需要一种同时考虑污染物综合污染水平的多因子评价方法，这里采用梅罗指数法（《梅罗指数评价法的修正及其应用》），其计算公式如下：

式中，P为土壤污染元素综合污染指数；为土壤中污染元素i的实测值；为土壤中污染元素i的评价标准值。

梅罗指数土壤污染评价标准如表二。

表二梅罗指数土壤污染评价标准

等级

梅罗污染指数

污染等级

I

P0.8

清洁（安全）

II

0.7P≤2.50

尚清洁（警戒线）

III

2.50≤4.25

轻度污染

IV

4.25≤7.20

中度污染

V

P7.20

重污染

但由于梅罗指数法过分突出极大值对水质污染的影响，会出现评价项目中只有一项指标值偏高，而其他指标值均较低也会使综合评分值偏高的情况，故对梅罗法进行改进如下：

其中

式中：

为权重值前n项的平均值，n值根据评价数据确定；为前n项的值；m为前n项中的项数。

对于以上8种重金属，取权重值如表三：

表三各重金属权重值

As

Cd

Cr

Cu

Hg

Ni

Pb

Zn

权重值

0.070

0.140

0.014

0.020

0.698

0.030

0.070

0.020

计算得到域不同重金属的污染程度评价结果如表四：

表四各区域重金属污染程度

As

Cd

Cr

Cu

Hg

Ni

Pb

Zn

生活区

II

III

III

III

III

II

III

IV

工业区

III

III

III

V

V

II

III

IV

山区

II

II

II

II

II

III

II

II

交通区

III

III

IV

V

V

III

III

IV

公园绿地区

II

III

II

III

IV

II

III

IV

（三）重金属污染原因分析

根据

（一）、

（二）中所得数据及结论，将各重金属污染的原因做如下分析：

1.对于除Ni以外的其他重金属元素，均有60%以上的重污染点分布在工业区，由此可认为重金属污染主要源自工业区，这主要是由于工业区中的工厂在生产过程中产生了含该种金属的废弃物，在排放时也没有进行必要的净化处理而直接排入水体或大气中，从而造成了污染的形成与扩散，

2.Hg、Cd、Pb、Zn四种金属元素的均有约20%重污染点分布在交通区，这些污染主要来自于公路上行驶车辆燃烧汽油或柴油所产生的尾气沉降及轮胎的磨损。

3.金属元素Zn在生活区和公园绿地区也有一定的分布，这可能是由于该区域所使用的肥料中含有Zn元素。

4.金属元素Ni的污染点较为分散，可认为该污染与人为因素关系不大，主要来自于土壤母质本身。

（四）重金属污染物的传播特征

重金属一般不会直接通过土壤进行传播，它一般通过如下方式进行传播：

1.空气传播

工业区中的工厂会排放