整理重金属污染源的数学建模.docx
- 文档编号:5401206
- 上传时间:2022-12-15
- 格式:DOCX
- 页数:17
- 大小:549.40KB
整理重金属污染源的数学建模.docx
《整理重金属污染源的数学建模.docx》由会员分享,可在线阅读,更多相关《整理重金属污染源的数学建模.docx(17页珍藏版)》请在冰豆网上搜索。
整理重金属污染源的数学建模
2)规划实施可能对环境和人群健康产生的长远影响。
一、安全评价
(二)环境影响经济损益分析的步骤
(2)区域、流域、海域的建设、开发利用规划。
环境影响篇章或说明
2.环境影响评价技术导则
①主体是人类;
2.环境影响评价技术导则
(5)污染防止措施能否达到要求。
发现规划环境影响报告书质量存在重大问题的,审查时应当提出对环境影响报告书进行修改并重新审查的意见。
2.环境影响评价工作等级的划分依据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
P
0.8
清洁(安全)
II
0.7
P≤2.50
尚清洁(警戒线)
III
2.50
≤4.25
轻度污染
IV
4.25
≤7.20
中度污染
V
P
7.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.空气传播
工业区中的工厂会排放出大量的工业废气,这些废气中可能含有大量的重金属颗粒,这些颗粒会随空气的流动而分散到城市各处,从而造成污染,另外,废气也可能会随降雨渗入土壤。
2.水源传播
工厂排出的工业废水中含有大量的重金属离子及沉淀物,这些物质均可随水渗入土壤,造成污染。
3.工业废渣和生活垃圾
废渣及生活垃圾的堆积会造成重金属元素的大量积聚。
另外,重金属的传播受地形因素的影响较为明显,从数据来看,该城市整体海拔较低,地势起伏不大。
相对1000米的取样间隔,300米的高度差不会因高空吸附或雨水冲刷造成污染物沉积异常。
但山势对风向的影响势必导致污染物浓度区域性变化。
如上图的标注,山脚①、②、③较易形成污染物富集;由于Ⅰ、Ⅱ山谷的存在,经过山谷的大风在A、B区域减速后极易产生污染物富集;此外,迎向污染源的山坡受到的污染会明显大于背向污染源的山坡。
从Hg、Cd的污染分布情况可以验证这几点(如图四):
(五)建立模型确定污染源
在进行污染源的确定之前先进行筛选,如果污染是由点污染源引起的,那么污染源就一定处于污染程度最大,即重金属浓度最高的采样点的附近。
由于城市中建筑物较多,人类活动频繁,导致局部气流复杂,因此通过空气传播占主要地位,又由于城市缺少地表水系,雨水及其它污水会通过下水管排出城外,故影响相对较小,因此只讨论通过空气传播的模型:
设
是t时刻点
处某种重金属污染物的浓度。
任取一个闭合曲面S,它所围的区域是Ω,由于气体扩散,由t到t+
时刻的这段时间内,通过S流入Ω的质量为:
由高斯公式得:
其中
分别是沿
方向的扩散系数。
由于衰减(例如吸收、代谢等),
的质量减少为:
其中
为衰减系数。
由于重金属污染物在土壤中很难降解,因此衰减系数可以取为0。
由质量守恒定律,在Ω内由于扩散和衰减的合作用,积存于Ω内物质的总量为
从另外的一个角度来看,Ω内由于浓度的变化引起的质量增加为:
显然
,即
由
设扩散源在点(
由假设可知,如果认为经过相当长的时间后,扩散已经停止,则物质处于平衡状态,有:
即:
由于污染物在空中传播时,地面上污染物的沉积量与z方向上的传播无关,故上述方程可以改写成:
因此问题转化成二维含源反演问题,可以用贝叶斯推理来估计参数
和
。
贝叶斯推理的基础是贝叶斯定理,它可以表述如下:
其中,
为模型参数,y为观测数据,
为参数的先验概率密度数,
为似然函数,
为参数的后验概率密度函数。
在贝叶斯推理中,p(
)表示在未获得测量数据之前,对模型参数分布的认识,主要来源于以往数据、经验和主观判断等。
代表模型参数拟合测量数据的程度,越大表示拟和效果越好,反之则差。
表示了获得测量数据以后模型参数的分布规律,即为在统计反演意义下的反问题的解。
其步骤通常分为以下三步:
1.确定先验概率密度数
一般情况下,对于参数
,如果仅知道它分布的区间
,那么可以假设先验概率为均匀分布:
=
假设参数之间相互独立则:
其中m为模型参数的个数。
2.确定似然函数
似然函数的构造对于贝叶斯推理的结果具有很大影响,一般可以认为测量误差服从正态分布
,此时似然函数可以表示为:
其中n为测量个数,
可以用解析解表示,一般情况下,
为数值解。
3.对后验概率密度函数抽样,进而获得参数的估计值
选取马尔科夫链蒙特卡洛法(MarkovchainMontCarlo,简记为MCMC)。
算法为MH算法,表示如下:
(1)设定参数初始值
并设定
;
(2)proposal分布
中生成
,计算接受概率:
(3)生成[0,1]之间均匀分布的随机数
;
(4)如果
,设置
;否则
。
当
(迭代次数)时,停止;否则
,继续第
(2)步。
其中,
为目标概率密度函数,在贝叶斯推理中为似然函数。
(六)模型的求解
通过MATLAB编程,使用MH算法进行50000次抽样,使目标坐标值趋于平稳,最终得到各重金属元素污染源的地理位置坐标如表五。
(图五为Hg元素进行抽样后的迭代曲线图,可以看到图线最终趋于一条平稳的直线。
)
表五各重金属元素污染源地理位置坐标
As
Cd
Cr
Cu
Hg
Ni
Pb
Zn
X值
2742
2382
3470
2707
2708
1647
2882
4152
Y值
7294
3693
2309
2295
2294
2729
3617
2299
(七)对城市地质环境的演变模式的探究
地质环境的演变模式是一个地域在一段较长时期内所表现出来的地质变化规律与趋势,因此,要对城市地质环境的演变模式进行研究,就必须掌握所研究城市在长期的地质变化过程中形成的规律,以研究题中城市土壤重金属污染的演变规律为例,就需要收集该城市每年生活、工业等重要污染源的废弃物排放量,以及每年的生物降解量,另外,该城市的水文特征及气候特征也会对重金属元素扩散产生重要的影响,因此还应该考虑降雨量及气候条件的信息。
拥有这些数据后,可以建立因子分析法,回归分析,曲线拟合等模型进行分析,进一步还可使用马尔科夫链等方法对该城市未来的污染状况进行预测与评估。
五、模型评价
优点:
1.通过改进后的内梅罗指数法引入了权值,避免了极大值对结果的较大影响,提高了数据的准确性和可信度。
2.通过马尔科夫链蒙特卡洛法(MCMC)简化了微分方程的求解,降低了计算时的复杂度。
3.通过对地形的综合分析,剔除了地形影响较大的区域,简化了模型。
缺点:
1.没有兼顾降雨、降解等其他对污染产生作用的因素的影响。
2.对污染的传播规律还有进一步探索的余地。
六、参考文献
[1]朱嵩、刘国华、毛根海、程伟平、黄跃飞,《利用贝叶斯推理估计二维含源对流扩散方程参数》,《四川大学学报》,第40卷第2期,2008年5月。
[2]李尚志、程继新,《数学建模教程》,2011年
[3]杨忠平、卢文喜、辛欣,《长春市城区表层土壤重金属污染来源解析》,
[4]赵彦峰、郭恒亮、孙志英,《基于土壤学知识的主成分分析判断土壤重金属来源》,《地理科学》,第28卷第1期,2008年2月
[5]李亚松,张兆吉,费宇虹,王昭,《内梅罗指数评价法的修正及其应用》,《水资源保护》,第25卷第六期,2009年11月
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 整理 重金属 污染源 数学 建模