承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

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

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

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

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：中国人民解放军国防科学技术大学参赛队员(打印并签名) ：1. 李腾骥2. 袁雪强3. 赵泉朴指导教师或指导教师组负责人(打印并签名)：杨文强日期：2011 年9 月9 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于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. 污染物浓度在取样前后较长一段时间内不会有显著的变化。

三、符号说明P: 土壤污染元素综合污染指数C i: 土壤中污染元素i的实测值S i: 土壤中污染元素i的评价标准值S: 任一闭合曲面M: S内指定所围的区域3M: 通过S流入Ω的质量1M: Ω内减少的质量2M: Ω内由于浓度的变化引起的质量增加32k: 衰减系数θ: 贝叶斯公式模型参数m : 模型参数的个数 y : 贝叶斯公式观测数据(),k T x x ：马尔可夫蒙特链卡洛法中的分布 (,)k x x α：接受概率()x π: 目标概率密度函数s : 距离污染源的距离 p ： 沉降系数 h : 海拔高度0h ：临界海拔r ： 重金属元素吸附系数四、模型建立与求解（一）空间分布的特征根据所给数据，对采样点按区域进行重新的排序，由此作出各采样点的地势分布图（图一），又由各种金属在各区域的浓度值，作出金属污染程度在各区域的分布图（以Pb （铅）为例，见图二，其它金属元素分布图见附件），由此分析可得到各金属元素在该城区的空间分布。

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

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

表一工业区各金属元素污染程度相关系数根据以上数据分析，可以得到各重金属元素大致的空间分布，分别为： 1. Hg 、Cd 、Pb 、Zn 、Cr 均主要分布于工业区和左下方的部分交通区；Cd Cr Cu Hg Ni Pb Zn As 0. 0. 0. 0. 0. 0. 0. Cd - 0. 0. 0. 0. 0. 0. Cr - - 0. 0. 0.69825 0. 0.69507 Cu - - - 0.983480.50278 0.66966 0. Hg - - - - 0. 0. 0. Ni - - - - - 0.577620. Pb ------0.图二 Pb 污染程度分布图2. As 、Cu 主要分布于工业区；3. Ni 的分布比较分散，没有明显的富集现象。

（二）重金属污染程度的划分为了全面、综合地反映污染物的整体污染水平，需要一种同时考虑污染物综合污染水平的多因子评价方法，这里采用内梅罗指数法（《内梅罗指数评价法的修正及其应用》），其计算公式如下：式中，P 为土壤污染元素综合污染指数；C i 为土壤中污染元素i 的实测值；S i 为土壤中污染元素i 的评价标准值。

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

但由于内梅罗指数法过分突出极大值对水质污染的影响，会出现评价项目中只有一项指标值偏高，而其他指标值均较低也会使综合评分值偏高的情况，故对内梅罗法进行改进如下：'max max2w C C C+=其中 1max nii CC m==∑式中：w C 为权重值前n 项的平均值，n 值根据评价数据确定；i C 为前n 项的值；m 为前n项中0h h =的项数。

对于以上8种重金属，取权重值如表三：P =计算得到域内不同重金属的污染程度评价结果如表四：（三）重金属污染原因分析根据（一）、（二）中所得数据及结论，将各重金属污染的原因做如下分析：1. 对于除Ni以外的其他重金属元素，均有60%以上的重污染点分布在工业区，由此可认为重金属污染主要源自工业区，这主要是由于工业区中的工厂在生产过程中产生了含该种金属的废弃物，在排放时也没有进行必要的净化处理而直接排入水体或大气中，从而造成了污染的形成与扩散，2. Hg、Cd、Pb、Zn四种金属元素的均有约20%重污染点分布在交通区，这些污染主要来自于公路上行驶车辆燃烧汽油或柴油所产生的尾气沉降及轮胎的磨损。

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

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

（四）重金属污染物的传播特征重金属一般不会直接通过土壤进行传播，它一般通过如下方式进行传播：1. 空气传播工业区中的工厂会排放出大量的工业废气，这些废气中可能含有大量的重金属颗粒，这些颗粒会随空气的流动而分散到城市各处，从而造成污染，另外，废气也可能会随降雨渗入土壤。

2. 水源传播工厂排出的工业废水中含有大量的重金属离子及沉淀物，这些物质均可随水渗入土壤，造成污染。

3. 工业废渣和生活垃圾废渣及生活垃圾的堆积会造成重金属元素的大量积聚。

另外，重金属的传播受地形因素的影响较为明显，从数据来看，该城市整体海拔较低，地势起伏不大。

相对1000米的取样间隔，300米的高度差不会因高空吸附或雨水冲刷造成污染物沉积异常。

但山势对风向的影响势必导致污染物浓度区域性变化。

图三城区地势图如上图的标注，山脚①、②、③较易形成污染物富集；由于Ⅰ、Ⅱ山谷的存在，经过山谷的大风在A、B区域减速后极易产生污染物富集；此外，迎向污染源的山坡受到的污染会明显大于背向污染源的山坡。

从Hg、Cd的污染分布情况可以验证这几点（如图四）：图四Cd污染程度分布图（五）建立模型确定污染源在进行污染源的确定之前先进行筛选，如果污染是由点污染源引起的，那么污染源就一定处于污染程度最大，即重金属浓度最高的采样点的附近。