《重金属污染源数学建模》由会员分享，可在线阅读，更多相关《重金属污染源数学建模（14页珍藏版）》请在金锄头文库上搜索。
1、2021高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规那么.我们完全明白，在竞赛开始后参赛队员不能以任何方式包括 、电子邮件、网上咨询等与队外的任何人包括指导教师研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规那么的, 如果引用别人的成果或其他公开的资料包括网上查到的资料，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规那么，以保证竞赛的公正、公平性。如有违反竞赛规那么的行为，我们将受到严肃处理。我们参赛选择的题号是从A/B/C/D中选择一项填写： A 我们的参赛报名号为如果赛区设置报名号的话：
2、 20001032 所属学校请填写完整的全名： 中国人民解放军国防科学技术大学 参赛队员 (打印并签名) ：1. 李腾骥 2. 袁雪强 3. 赵泉朴 指导教师或指导教师组负责人 (打印并签名)： 杨文强 日期： 2021 年 9 月 9 日赛区评阅编号由赛区组委会评阅前进行编号：2021高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号由赛区组委会评阅前进行编号：赛区评阅记录可供赛区评阅时使用：评阅人评分备注全国统一编号由赛区组委会送交全国前编号：全国评阅编号由全国组委会评阅前进行编号：基于MCMC的城市表层土壤重金属污染分析摘 要：本文针对某城市城区表层土壤重金属污染分析的问题，从
3、各分区重金属元素污染的空间分布、污染程度、污染原因及传播途径等多方面进行了全面的统计与研究，同时采用贝叶斯推理和马尔科夫链蒙特卡洛法MCMC，利用MATLAB进行编程，对所设计的模型进行求解，得到了各重金属污染源的地理位置坐标，最后对城市地质环境演变模式的问题进行了探究，得出了初步的结论。 首先，根据所给数据信息，利用MATLAB及EXCEL作图，得到了该地区的地势图及各重金属污染程度的分布图。 对于问题1，通过统计各重金属的重污染点在各区域的分布和计算不同区域各中金属污染程度之间的相关系数，得到了各重金属的大致空间分布，并采用改良的内梅罗指数法，得出不同区域重金属的污染程度。 对于问题2，在
4、问题1所得数据及结论的根底上，结合实际情况，对各重金属污染的主要原因进行了较为全面的分析。 对于问题3，通过对污染的传播特征的研究，基于城市的地质特性，得到空气传播为污染主要传播方式的结论，并由此建立关于气体传播的概率模型，采用贝叶斯推理和马尔科夫链蒙特卡洛法，利用MATLAB进行50000次抽样，得到收敛的坐标值，从而计算出各重金属污染源的地理位置坐标分别为：As2742,7294、Cd2382,3693、Cr3470,2309、Cu2707,2295、Hg(2708,2294)、Ni1647,2729、Pb2882,3617、Zn4152,2299。 对于问题4，以题中城市土壤重金属污染的
5、演变规律为例进行了一定深度的探究，并对所需信息及模型设计方案进行了初步的讨论。关键词：内梅罗指数法 相关系数 贝叶斯推理 MCMC法 MH算法 一、问题的提出与分析一背景介绍随着城市经济的快速开展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将
6、所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土010 厘米深度进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。二要解决的问题现给出采样点的位置、海拔高度及其所属功能区等信息，8种主要重金属元素在采样点处的浓度及8种主要重金属元素的背景值，解答以下问题：1. 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。2. 通过数据分析，说明重金属污染的主要原因。3. 分析重金属污
7、染物的传播特征，由此建立模型，确定污染源的位置。4. 分析所建立模型的优缺点，并说明为更好地研究城市地质环境的演变模式，还应收集什么信息，有了这些信息，如何建立模型解决问题。三问题分析整个问题可以分为三类分别进行求解：问题1和问题2可将其归结为对所给数据的统计与分析，可以通过作图或列表等直观的形式针对不同金属和区域进行相关统计，再对统计的结果加以比对与分析，即可得出结论。问题3确定污染源的位置实质就是要找到使周围污染到达最大的点，这样可以转化为概率模型，在统计数据的根底上，运用贝叶斯推理马尔科夫链蒙特卡洛法可加以计算和确定。问题4开放类问题，可根据现有资料进行答复。二、模型假设1. 题目中所给
8、出的样本均具有普遍性及代表性，而非所测地区中的特殊情况。2. 该地区的空气流动没有显著的方向性特征。3. 污染物浓度在取样前后较长一段时间内不会有显著的变化。三、符号说明: 土壤污染元素综合污染指数Ci: 土壤中污染元素i的实测值Si: 土壤中污染元素i的评价标准值: 任一闭合曲面: 内指定所围的区域: 通过S流入的质量: 内减少的质量: 内由于浓度的变化引起的质量增加: 衰减系数: 贝叶斯公式模型参数m: 模型参数的个数y: 贝叶斯公式观测数据：马尔可夫蒙特链卡洛法中的分布：接受概率: 目标概率密度函数s : 距离污染源的距离p： 沉降系数h: 海拔高度：临界海拔： 重金属元素吸附系数四、模
9、型建立与求解一空间分布的特征根据所给数据，对采样点按区域进行重新的排序，由此作出各采样点的地势分布图图一，又由各种金属在各区域的浓度值，作出金属污染程度在各区域的分布图以Pb铅为例，见图二，其它金属元素分布图见附件，由此分析可得到各金属元素在该城区的空间分布。3图一 采样点地势分布图海拔m图二 Pb污染程度分布图以Pb为例进行分析，经统计，Pb的重度污染区域主要集中于图的左下方，且Pb的重度污染点在工业区占62.8%，另有20.9%重度污染点那么分布在左下方的交通区中，由此可得Pb的污染区域主要分布于工业区和左下方的局部交通区，对于其它金属元素，可以用相同的方法，还可以通过计算不同区域各金属元
10、素污染程度的相关系数，找出各金属元素污染之间的联系，这样可以有效提高分析效率和准确性。例如对表一中工业区的各金属元素污染程度相关系数其它区域见附件的比拟可以得出，Cr与Cu、Cr与Hg、Cu与Hg、Cd与Pb、Cd与Zn、Pb与Zn相关系数较大，可以认为它们具有相同或相近的空间分布。表一 工业区各金属元素污染程度相关系数CdCrCuHgNiPbZnAsCd-Cr-Cu-Hg-Ni-Pb-根据以上数据分析，可以得到各重金属元素大致的空间分布，分别为：1. Hg、Cd、Pb、Zn、Cr均主要分布于工业区和左下方的局部交通区；2. As、Cu主要分布于工业区；3. Ni的分布比拟分散，没有明显的富集
11、现象。二重金属污染程度的划分为了全面、综合地反映污染物的整体污染水平，需要一种同时考虑污染物综合污染水平的多因子评价方法，这里采用内梅罗指数法?内梅罗指数评价法的修正及其应用?，其计算公式如下：式中，P为土壤污染元素综合污染指数；Ci为土壤中污染元素i的实测值；Si为土壤中污染元素i的评价标准值。内梅罗指数土壤污染评价标准如表二。表二 内梅罗指数土壤污染评价标准等级内梅罗污染指数污染等级IP0.8清洁平安IIP2.50尚清洁警戒线III0P轻度污染IV重污染但由于内梅罗指数法过分突出极大值对水质污染的影响，会出现评价工程中只有一项指标值偏高，而其他指标值均较低也会使综合评分值偏高的情况，故对内
12、梅罗法进行改良如下：其中 式中：为权重值前n项的平均值，n值根据评价数据确定；为前n项的值；m为前n项中的项数。对于以上8种重金属，取权重值如表三：表三 各重金属权重值AsCdCrCuHgNiPbZn权重值计算得到域内不同重金属的污染程度评价结果如表四：表四 各区域内重金属污染程度AsCdCrCuHgNiPbZn生活区IIIIIIIIIIIIIIIIIIIIV工业区IIIIIIIIIVVIIIIIIV山区IIIIIIIIIIIIIIIII交通区IIIIIIIVVVIIIIIIIV公园绿地区IIIIIIIIIIIVIIIIIIV三重金属污染原因分析根据一、二中所得数据及结论，将各重金属污染的原因做如下分析：1. 对于除Ni以外的其他重金属元素，均有60%以上的重污染点分布在工业区，由此可认为重金属污染主要源自工业区，这主要是由于工业区中的工厂在生产过程中产生了含该种金属的废弃物，在排放时也没有进行必要的净化处理而直接排入水体或大气中，从而造成了污染的形成与扩散，2. Hg、Cd、Pb、Zn四种金属元素的均有约20%重污染点分布在交通区，这些污染主要来自于公路上行驶车辆燃烧汽油或柴油所产生的尾气沉降及轮胎的磨
