重金屬污染源的数学建模.docxVIP

2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话） 所属学校（请填写完整的全名）：中国人民解放军国防科学技术大学参赛队员 (打印并签名) ：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可将其归结为对所给数据