2011国赛数学建模优秀论文A题
2011高教社杯全国大学生数学建模竞赛A题
2011高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A
我们的参赛报名号为（如果赛区设置报名号的话）：
所属学校（请填写完整的全名）：_______________
参赛队员(打印并签名) ：1.
2.
3.
指导教师或指导教师组负责人(打印并签名)：
日期：2011年9月12日赛区评阅编号（由赛区组委会评阅前进行编号）：
2011高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号（由赛区组委会评阅前进行编号）：
全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：。
2011年全国大学生数学建模竞赛A题论文
城市表层土壤重金属污染分析摘要土壤作为城市环境的重要组成部分，不仅提供人类生存所需的各种营养物质，而且接受来自工业和生活废水、固体废物、农药化肥、及大气降尘等物质的污染．很容易导致金属元素的蓄积，从而造成土壤重金属的污染．本文讨论了某城市表层土壤重金属污染的空间分布分布状况，地区污染程度，以及污染传播特征，污染源等，建立了相应的几何与数学模型或算法，得到了较好的结果，为防治城市表层重金属污染，保护和提高土壤资源和生态环境，提供参考．对于问题一：通过给定数据的相关分析，不考虑地形高低对污染浓度变化的影响，用Matlab 软件编程绘制个重金属元素污染浓度空间分布三维网格图和二维等高线图，综合研究该城市各功能区的空间分布以及污染程度分布．建立了Muller 地积累指数模分析模型：)]/([log 2Bn C Fn ⨯=ℜ，确定污染程度水平分级标准，通过统计计算，分析了各重金属在不同功能区的污染状况及程度．结论是：主干道路区和工业区的重金属元素的污染最严重，其他次之．对于问题二，为说明重金属元素污染的主要原因，采用单因子指数模型和内梅罗综合指数模型进行综合指标评价分析，结合问题一中统计数据进行综合分析，得到个重金属元素在各功能区及城区的综合污染程度指标．污染最严重的功能区是主干道路区，其次按照污染程度从大小的顺序依次为：工业区、生活区、公园绿地区、山区．主干道路区土壤表层重金属元素含量很高，且种类多．根据地区的差异性和元素的特殊性，分析出重金属污染Hg 和Cu 污染是最严重的污染源，且污染最严重的地区在主干道路区和工业区．这些污染主要由于含铅汽油的燃烧、汽车轮胎磨损产生的含锌粉尘、工业污水的排放、生活废水的排放、化肥农药的多度使用、金属矿山的开采．详细情况见正文．对于问题三，为了找出该城区的污染源，在分析出重金属元素的主要传播特征之后，考虑大气空间传播情况，建立了微分方程模型，通过模型求解分析，用其等效的向内（向污染浓度较高的方向搜索）搜索算法，计算确定了重金属元素主要污染源的位置，其中As 较严重的中心污染源坐标分别为：（5291，7349，10）、（12696,3024，27）、（18134、10046、41）、（17814,10707,64）、（27700,11609,165）．这五个污染源主要分布在主干道路区．（5291，5739，10），（12696,3024，27），（17814,10707，64）分布在工业区，其它两种污染源分布在生活区．其余元素的中心污染源见正文．对于问题四，需对前面所建立的模型进行分析与评价并进行模型的优化，在详细分析了前三个问题的求解模型及过程之后，评价出所建立模型的优缺点．在问题三中，重金属元素除了在大气中传播以外，还通过水土流动传播．另外，前几个模型都是静态的，但污染物传播的过程与时间有关，是一个动态的过程．最后建立了一个扩算方程模型进行优化，能够为更好的研究城市地理环境的演变模式做贡献．关键词：重金属污染 地积累指数模型 单因子指数模型 内梅罗综合污染指数 微分方程模型一、问题重述1．1 问题背景随着工业发展和城市化进程的加剧，通过交通运输、工业排放、市政建设和大气沉降等造成城市重金属污染越来越严重．对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究城市不同功能区表层土壤重金属污染特征和污染空间分布性，以便更好的研究人类活动影响下城市地质环境的演变模式．本文就如何应用查证获得的海量数据资料展开城市环境质量评价，研究地质环境的演变模式建立数学模型．附录1列出了采样点的位置、海拔高度及其所属功能区等信息，附录2列出了8种列出了8种主要重金属元素在采样点处的浓度，附录3列出了8种主要重金属元素的背景值．1．2 需要解决的问题有(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度．(2) 通过数据分析，说明重金属污染的主要原因．(3) 分析重金属污染物的传播特征，由此建立数学模型，确定污染源的位置．(4) 分析所建立模型的优缺点，为更好的研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析该题目一方面通过GPS记录了该城市大量样本点的位置以及所属功能区，再应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据，通过这两个表的数据就大致可以提取出一些对于解决问题的重要信息，另一方面，题目给出了自然区各样本点的重金属元素的背景值，作为重金属污染情况的指标．对于分析研究各个样本点的污染程度至关重要．利用Matlab软件进行三维网格图和等高线图的制作并结合相关的数据统计分析，可以分析该城区不同区域重金属的污染程度．后面利用地积累指数法和内梅尔综合评价指数对数据进行处理，分析污染严重的功能区和重金属．结合图形的分析以及模型的建立综合分析重金属污染物的传播特征．接着对模型进行一定的优化处理，使得处理的结果更准确．三、模型假设1、假设题目所给的数据合理正确．2、该区域的划分是稳定的，不会出现大的变动．3、不考虑观测误差、随机误差和其他外在因素所产生的误差．4、重金属在大气中无穷空间扩散，不受风的影响，其扩散服从热传导定律．5、重金属污染程度连续变化，大气中重金属元素浓度连续变化．6、界限不明显区域有扩大、缩小、消失的过程，穿过大气进入仪器的重金属含量只有浓度大小之分，浓度大小由仪器灵敏度确定．四、变量与符号说明eo lg地积累指数n ()8,7,6,5,4,3,2,1=n 分别表示As,Cd,Cr,Cu,Hg,Ni,Pb,Zn 元素Fn 污染物重金属元素n 的浓度 Bn第n 种重金属元素的背景值上限P 综综合污染指数 n C重金属n 的实测值（ug/g ） max (/)n n C S 重金属污染物中污染指数最大值 (/)n n wg C S重金属污染物中污染指数平均值 n χ 重金属污染物n 的环境质量指数；n α 重金属污染物n 的实测值 n β 重金属污染物的评价标准． Ω 重金属元素通过的平面t 时间 h 海拔高度 V体积五、模型建立与求解针对问题一，首先想到的是用Matlab 软件编程，进行三维网格图、三维曲面图、等高线图和散点图的制作．5．1 问题(1)的分析、模型建立与求解： 5．1．1 问题（1）的分析对于问题一，首先来分析一下， 要给出8种主要重金属在该城区的空间分布， 就必须确定每个重金属元素与他们所对应的地区之间的联系．刚好题目给出了每个样本点的各元素浓度，那么 是不是可以将每种重金属元素含量浓度含量与该目标点所在的功能区建立联系？由此 想到利用Matlab 软件画出每种元素在该城区的三维曲面和空间曲面图．同时 在分析不同区域重金属的污染程度时，考虑到这个污染程度是否可以量化．并且是否能够建立一种模型将这种指标量化．这道问题还要求考虑每个功能区的污染程度， 知道每个功能区的每种重金属污染程度又是不一样的．那 通过什么指标来判断每个功能区的污染程度大小，这也是 为什么用权重作为评价每个功能区的污染程度的指标．5．1．2 问题（1）的模型建立该城区受这八种重金属元素As 、Cd 、Cr 、Cu 、Hg 、Ni 、Pb 、Zn 污染程度不一样．题目提供了每种重金属元素的背景值，那么 怎么样利用这些背景值和每种元素相关的浓度确定不同区域重金属的污染程度？所以 需要找出一种方法来准确的分析该城区内不同区域重金属的污染程度，并且最好能够量化．建立8种主要重金属元素在该城区的空间分布模型如下：引入了一种用于研究沉积物及其他物质中重金属污染程度的该区内不同地域重金属的污染程度的定量指标——地积累指数又称Muller 指数法，Muller 指数法表达式为:)]/([log 2Bn C Fn ⨯=ℜ式中Fn 表示污染物重金属元素n 的浓度；Bn 表示第n 种重金属元素的背景值上限，C 为考虑各地岩石差异可能会引起背景值的变动而取得一系列系数(一般取值为1．5)，用来表征沉积特征、岩石地质及其他影响．Muller 地积累指数评价和分级标准分级标准具体详见表1表1：地积累指数分级标准地积累指数ℜ 分级污染程度105≤ℜ<6及严重污染 54≤ℜ< 5强-及严重污染 43≤ℜ< 4强污染 32≤ℜ< 3中等-强污染 21≤ℜ< 2中等污染 10≤ℜ< 1轻度-中等污染 0≤ℜ 0无污染 该方法指标主要是通过每种重金属元素测得的实际浓度以及他们的相关背景值，计算出每种元素的地积累指数．然后根据上面这张表 就可以判断出每种元素的污染级别，这样就可以对每种元素的污染情况进行分析．然后 再利用Matlab 软件对题目所给数据进行处理，画出相应的网格曲面图和等高线曲线图．这里需要对Matlab 进行编程，首先利用每个样本点的横坐标、纵坐标、海拔高度建立等高线图，程序语句见附录一．通过该图，可以直观的看到该城区各功能区的空间分布．但是这张图不能反映出8种主要元素在城区的污染情况， 需要借助于各种主要元素的浓度．所以 需要再建立一张等高曲线图以及相应的网格曲面图，将主要元素的浓度作为第三坐标，命令语句见附录一．5．1．3 问题（1）的求解过程首先通过Matlab 软件，调用每个样本点的位置相关数据．就是以海拔为第三坐标，并且对每个功能区进行颜色区分，画出该城区每个功能区的二维等高线图．最后把每个样本点显示在图上．得到如下这张图：图一：重金属As空间二维等高线分布图这张图只反映出了该城区各功能区的空间分布，还不能看出每种重金属污染的情况．将每种重金属元素的浓度在图上反应出来，做出该城区重金属污染的二维等高线图．具体程序语句见附录二，得到如下这张图：图二：重金属As分布平面图同时为了对应这张As含量分布平面图，也画出了三维网格曲面图（图三）．图三：重金属As含量分布的空间三维图从空间三维图三中可以看到，有一处的波峰很高说明该处污染情况很严重，有二处处于波峰说明污染情况比较严重的主要有二处，还有一处面积比较广且所处高度稍微低一点这表明该处所受污染情况相对严重且污染的范围较广；同样分析二维等高线图二，图中有一处等高线之间的间距越来越密集且颜色很深表明该处受污染情况很严重，有二处等高线比较密集颜色相对较深表明这二处的污染情况相对严重，还有一处等高线间的距离较密集但是所包围的面积较广说明该处的污染也较严重且污染的面积很广．再结合前面的数据他们中心污染源的坐标分别为：（5291，5739），（12696,3024），（17814,10707）．都是分布在工业区，还有一处污染级别不是特别严重,但是在该处存在着污染源，此处刚好是山林密集区．通过观察图三，会发现刚好有三个点处于波峰，还有个点波峰稍微偏低，但还是能很直观的看出来．再来看一下，Cd这种重金属的城区各功能区的二维等高线图，分布平面图，空间分布图（图四、图五）：图四：重金属Cd空间二维等高线分布图图五：重金属Cd含量空间分布平面图以及相应的三维网格曲面图（图六）：图六：重金属Cd含量空间分布图从空间分布图六中可以看到，污染情况比较严重且面积比较广的主要有一处，还有五处污染也相对严重．以及几处小的污染；同样从二维分布图五可以看出等高曲线所谓面积有一处颜色很深，说明该区域污染情况很严重，同时也观察到又五处等高曲线所围的面积颜色比较深，这说明了这五处区域污染情况相对严重，很明显的是有一处等高曲线所围成的面积比较广且颜色较深，表明了该区域有一处污染情况较严重且污染面积比较广，由此可见不管是从二维还是三维图形进行分析的结果是相吻合的．再结合前面的数据它金属Cd中心污染源的坐标为：（22304,10527）．分布在主干道路区，还有一处污染级别不是特别严重．再观察图三，会发现刚好有三个点处于波峰．如此，通过同样的方法，都能够得到对其它六种种重金属在该城区的空间分布以及污染情况的了解（参见附录三）通过观察每种元素的三维曲面图以及等高曲线图．很容易观察到，每种重金属对该城区都存在或大或小的污染．其中有些地区是存在多种重金属污染，并且污染很严重，通过观察这8张图会发现这六种元素Cd,Cr,Cu,Hg,Ni,Pb 在横坐标在[3000,4000],纵坐标在[3000,6000]这个区域内含量都非常高，大致可以判定这段区域属于重度污染区．下面将题目中所给的数据用excel进行分类处理，得到样本点的地积累指数．然后运用数学统计法得到各种元素污染程度数据分布表，通过这些表就可以确定该城区内不同区域重金属的污染程度．统计该表时，是通过统计每个功能区的总样本点个数，然后通过地积累指数法分别计算出每种样本点的地积累指数，并判断他们的所在的污染级别．然后统计每种污染级别下，各功能区的污染点数占总点数的百分比也就是说的权重，通过该权重就能够分析出每种重金属元素的污染程度大小，以及污染所波及的范围．从而得到每种重金属元素污染最严重的地区．通过Excel对数据运算，得到重金属元素As 污染情况分布表：表二：As污染程度分布数据表下面通过同样的数据处理，得到Cd污染程度数据分布表：表三：Cd污染程度数据分布表其它六种元素的污染程度数据分布表见附录三．表中数值0的意义是在该污染级别下不存在观测的样本点．这是个大样本事件，可以认为该级别污染很轻微，甚至不存在这种级别的污染．而百分比越大，就说明在该污染级别下涉及的样本点比较多，污染波及范围较广．5．1．4问题（1）的结果分析5．1．4．1 As这种重金属污染情况分析由该表可以看出各个区域受As的污染程度，其中一类区即是生活区31．82%无污染，63．64%轻度—中度污染，4．55%为中等污染，无强污染和及严重污染的情况；二类区即是工业区38．89%不受重金属污染，52．78%受轻度—中度污染，5．56%受中等污染，2．78%受中等—强污染；三类区即是山区大多数不受污染，只有15．15%受轻度—中度污染，1．51%受中等污染；四类区即是主干路区47．83%不受污染，50．00%受轻度—中度污染，0．72%受中等污染，1．45%受中等—强污染；五类区即是公园绿地区大多数受轻度—中度污染，25．71%不受污染，2．86%受中等污染．再结合相应的几何图形，会发现在四区存在三个很明显的污染源，在污染源附近会看到，有很多二区的样本点．有个别一区的点，说明这种元素对一区的影响相对来说轻点．所以由分析可知工业区受污染最严重，污染面积达到了61．11%，其次是生活区、主干道路区，生活区污染面积都达到了50%以上，也就是说这三个区有至少一半的土壤受到该元素的不同程度的污染．其余功能区受污染程度就次之．5．1．4．2 Cd这种重金属污染情况分析由该表可以看出各个区域受Cd的污染程度，其中一类区即是生活区29．55%无污染，54．55%轻度—中度污染，13．64%为中等污染，无强污染和及严重污染的情况；二类区即是工业区16．77%不受重金属污染，44．44%受轻度—中度污染，30．56%受中等污染，8．33%受中等—强污染；三类区即是山区大多数不受污染，只有75．76%受轻度—中度污染，21．21%受中等污染；四类区即是主干路区23．91%不受污染，44．2%受轻度—中度污染，26．09%受中等污染，5．07%受中等—强污染；五类区即是公园绿地区大多数受轻度—中度污染，48．57%不受污染，31．43%受轻度-重度污染，11．43%受中等污染，8．57%受中等-强污染．再结合相应的几何图形，会发现在四区存在三个很明显的污染源，在污染源附近会看到，有很多二区的样本点．有个别一区的点，说明这种元素对一区的影响相对来说轻点．所以由分析可知工业区受污染最严重，污染面积达到了61．11%，其次是生活区、主干道路区，生活区污染面积都达到了50%以上，也就是说这三个区有至少一半的土壤受到该元素的不同程度的污染．其余功能区受污染程度就次之．5．1 这六种重金属Cr、Cu、Hg、Ni、Pb、Zn污染情况分析由于重金属含量越多，说明该地区的重金属污染程度越严重．Cr污染最严重的有一处，该中心污染源的坐标为：（3299，6018），所在地区为主干道路区，一定程度上波及到了生活区．一区和四区存在强-及严重污染，一区波及面积达到了52．27%，四区波及面积达到了41．3%，该元素污染最严重的就是生活区．Cu污染最严重的有一处，该中心污染源的坐标为：（2427，3971），所在地区为生活区，一定程度上波及到了工业区和主干道路区．一区和四区存在及严重污染，一区污染波及范围达到了84．09%，四区污染波及范围达到了84．06%，该元素污染最严重的就是生活区和主干道路区．Hg污染最严重的有一处，中心污染源的坐标为：（3299，6018），所在地区为主干道路区，一定程度上波及到了生活区．一区和四区存在及严重污染，一区污染波及范围达到了54．55%，四区污染波及范围达到了50．74%，该元素污染最严重的就是主干道路区．Ni污染最严重的有一处，中心污染源的坐标为：（3299，6018），所在地区为主干道路区，一定程度上波及到了生活区．一区、二区和四区存在及严重污染，一区污染波及范围达到了90．91%，二区污染波及范围达到了94．44%，四区污染波及范围达到了93．48%，该元素污染最严重的就是主干道路区和生活区．Pb污染最严重的有二处，中心污染源的坐标为：（2383，3692）、（5062，4339），所在地区为生活区和主干道路区，一定程度上波及到了工业区．一区和四区存在及严重污染，一区污染波及范围达到了52．73%，四区污染波及范围达到了80．87%，该元素污染最严重的就是主干道路区．Zn污染最严重的有一处，中心污染源的坐标为：（14065，10987），所在地区为主干道路区，一定程度上波及到了工业区．四区存在及严重污染，四区污染波及范围达到了67．39%，该元素污染最严重的就是主干道路区．所以，该城区不同区域重金属污染最严重的区域是主干道路区和工业区，其次是生活区、公园绿地区、山区．5．2 问题（2）的求解：5．2．1问题（2）的分析通过问题一的分析，可粗劣的判断哪几种元素污染比较大，哪个功能区污染比较严重，但是怎么样才能具体到哪个功能区污染最严重，被污染的功能区的土壤哪种重金属污染最严重？所以，针对问题二给出的数据分析，不能简单的进行数据处理．为了使得所寻找出来的原因更有说服力，用两种方法分别进行说明和验证，还要进行综合指标评价．最后确定了最严重的污染地区以及污染最严重的相关元素，根据地区的差异性和元素的特殊性，才能说明重金属污染的主要原因．5．2．2数据的统计分析首先通过数据的处理，建立每个功能区各重金属元素的污染程度样本所占的百分比表．一功能区的相关百分比数据如下：表四：一功能区各重金属污染程度所占百分比在此功能区从总体来看，重金属污染程度处于中等-强污染，其中主要污染来自重金属元素Ni，另外在该区域有少数地方Cu污染及严重．表五： 二功能区各重金属污染程度所占百分比在该功能区重金属Hg 和重金属Ni 的污染极为严重，尤其是在该区域的某些地方．由此可见，在此功能区照成重金属污染的罪魁祸首为重金属元素Hg 和重金属元素Ni ． 通过这两张表， 会发现有些地区之所以污染严重，主要是因为个别元素污染所导致的．所以 要分析重金属污染的原因，就得分析该重金属在该功能区为什么会产生污染．其它三个功能区各重金属污染程度百分比见附录三．通过该附录表 可以看到在该功能区里，重金属污染程度较轻，污染等级集中在轻度污染及以下． 再观察功能区四，重金属污染十分严重，大多数重金属污染元素都集中在在各个功能区，但是在这个功能区，Pb 污染级别比较轻，没有中度甚至以上级别的污染． 再看功能区五，从总体上分析，该地区重金属污染中等、强污染几乎没有，正因如此造成重金属污染的少数种类重金属元素就凸显出来了——Ni 元素和Hg 元素．纵观整体，分析所有的功能区， 很容易发现造成重金属污染的主要重金属元素，他们就是Ni 元素和Hg 元素．知道前面的数据分析理由不充分，只是一个粗劣的判断．为了综合前面处理的数据，准确找出各个功能区污染的主要元素． 需要利用单因子指数法和内梅罗综合污染指数法进行综合评价．5．2．3 单因子指数法和内梅罗综合污染指数法的建立与求解单因子指数法是目前国内土壤重金属的单项污染指数评价方法之一，其计算公式为：n n n βαχ=，式中n χ为重金属污染物n 的环境质量指数；n α为重金属污染物i 的实测值；n β为重金属污染物的评价标准．n χ﹥1表示污染；n χ=1或n χ﹤1表示无污染；且n χ值越大，则污染物越严重．为了更全面的反应各重金属对土壤的不同作用．突出高浓度重金属对环境质量的影响， 采用内梅罗综合污染指数法．其计算公式为:2)/(/22max n wgn n n S C S C P +=）（综,式中max )(n n βα表示重金属污染物种污染指数nn βα的最大值；(/)n n wg C S 表示重金属污染物中污染指数的平均值．土壤污染水平分级标准采用国家土壤环境二级标准．土壤污染综合污染指数分级标准为综合污染指数>3为重污染，2~3为中污染，1~2为轻污染，0．7~1为警戒级，≤ 0．7为安全级．下面为了找到每种元素在该城区的综合污染指数，借助于Matlab 循环计算．编写如下系列命令见附录七．运行程序结果为As 综合污染指数：p=4．0093，分别运行另外几种程序，得到每种重金属元素的综合评价指标，简单结果如下表：。
2011年数学建模A题优秀论文
2011年数学建模A题优秀论文城市表层土壤重金属污染分析摘要随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
对于问题一我们首先用EXCEL 对数据进行处理，然后用MATLAB 等软件对所给的数值进行空间作图，然后分别做出了八种重金属元素的空间分布特征，我们利用综合指数（内梅罗指数）评价的方法，建立模型：ij j j P C S =22,,max ()2N j ave j P P P =+区域生活区 工业区 山区 交通区 公园绿地区 污染程度重污染重污染轻度污染重污染中度污染析，并作出了不同重金属浓度与海拔的分布图；然后结合第一问给出的空间分布图和区域散点图，参照主要重金属含量土壤单项污染的指数，分析得出各重金属污染的主要原因主要来自工业区、交通区和生活区。
对于问题三我们建立模型，建立目标函数；=jmk H P C e-⨯⨯综应用MATLAB 软件对数据处理，作出可能为污染源的三个位置；然后用MATLAB 进行三次拟合后，得到污染源的位置。
对于问题四，我们在已有信息的基础上，还应收集不同时间内的样点对应的浓度以及各污染源重金属的产生率。
根据高斯浓度模型建立高斯修正模型，得到浓度关于时间和空间的表达式ut e C C -⋅=0。
关键词：重金属污染 内梅罗污染指数 相关性分析 污染源 高斯浓度一.问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。
2011国赛数学建模优秀论文A题
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：04505所属学校（请填写完整的全名）：德州学院参赛队员(打印并签名) ：1. 翟明元2. 张小凤3. 邹菲指导教师或指导教师组负责人(打印并签名)：赵学杰高秀莲日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于系统综合评价的城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染问题，首先对各重金属元素进行分析，然后对各种重金属元素的基本数据进行统计分析及无量纲化处理，再对各金属元素进行相关性分析，最后针对各个问题建立模型并求解。
针对问题一，我们首先利用EXCEL 和 SPSS 统计软件对各金属元素的数据进行处理，再利用Matlab 软件绘制出该城区内8种重金属元素的空间分布图最后通过内梅罗污染模型：2/12max22⎪⎪⎭⎫ ⎝⎛+=P P P 平均综，其中平均P 为所有单项污染指数的平均值，m ax P 为土壤环境中针对问题二,我们首先利用EXCELL 软件画出8种元素在各个区内相对含量的柱状图，由图可以明显地看出各个区内各种元素的污染情况，然后再根据重金属元素污染来源及传播特征进行分析，可以得出工业区及生活区重金属的堆积和迁移是造成污染的主要原因，Cu 、Hg 、Zn 主要在工业区和交通区如公路、铁路等交通设施的两侧富集，随时间的推移，工业区、交通区的土壤重金属具有很强的叠加性，受人类活动的影响较大。
数模2011国赛A题优秀论文重金属污染
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：13257006所属学校（请填写完整的全名）：武汉纺织大学参赛队员(打印并签名) ：1. 刘浩2. 郭子雷3. 房旭指导教师或指导教师组负责人(打印并签名)：日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析与预测摘要本文通过对城市表层土壤重金属各指标的分析，建立各种适宜的数学模型，并通过matlab软件计算，得出所需要的结果，达到题目的要求。
对于问题一运用matlab软件将绘制了城区土壤中的8种主要重金属的浓度的三维坐标图，动态直观显示城区的重金属空间分布图。
运用潜在生态风险指数来反映各个类区的各种重金属的污染程度。
首先导出不同类区的样本数据（包括样点数目N，As等八种重金属的浓度C(i)以及背景浓度Cn（i），取各类区样点的各种重金属的平均浓度作为该类区的浓度值，最后运用hakanson指数的计算公式算出各个类区的综合生态风险指数RI和单因子风险参数Er(i)。
反映了不同类区的不同重金属的污染程度，得出RI：交通区=575.47（重）〉工业区=523（重）〉公园绿地区=248.65（中）〉生活区=233.7（中）〉山区=113.18（低）。
2011数学建模国赛A题
城市表层土壤重金属污染分析与预测摘要近几十年，由于工业化和城市化的进程加快，土壤重金属污染日益严重。
本文针对重金属污染问题，主要包括四个相关问题，分别为确定重金属元素的空间分布，并分析该城区内不同区域重金属的污染程度，说明重金属污染的主要原因，确定污染源的位置，对模型进行评价与扩展并分析该城市地质环境演变模式。
我们在合理的假设下，先绘制出重金属元素空间分布图，再根据单因子指数与内梅罗综合污染指数分析方法对各功能区污染程度分析，说明重金属污染的主要原因。
接着，在此基础上建立二维污染物空间扩散模型，采用搜索算法在污染较为严重的区域搜索出适量的样本点，再采用多元非线性回归方法对二维扩散模型进行拟合，求解相应的参数。
最后通过对第三问模型优缺点的分析，建立时空扩散模型，更加精确反映城市地质环境的演变情况。
针对问题一，首先用软件Surfer绘制出该城区的空间分布图，得出该城区大致的地形地貌，再将纵坐标改为八种重金属元素的含量，采用克里格插值，得到各种重金属在城区的分布图，得出各重金属在不同功能区的污染物空间分布特点。
再选用单因子指数和内梅罗综合指数评价污染状况，由数据可以反映出工业区和交通区Hg和Cu严重超标，各功能区污染程度分别为：生活区、山区、交通区及公园绿地区受到了轻度污染，工业区受到了重度污染。
针对问题二，对各功能区的污染物含量进行分析，分析推测该城市可能存在规模较大或者数量较多的Hg,Cu,Zn，Pb严重超标的工厂。
为进一步揭示表层土壤重金属污染来源，我们对三个污染相对严重的功能区的超标较严重的重金属元素Hg,Cu,Zn,Pb进行相关性分析，运用SPSS软件计算出各相关系数，揭示出不同功能区内土壤表层重金属污染的不同来源，分析可得，交通区污染主要来源于汽车尾气排排放，工业区污染主要来源于工厂的超标排放。
针对问题三，根据污染物扩散的特征，建立一维污染物扩散模型：Q(x)=Q0eαx−12x2x m，以此推广成污染物的高斯二维空间扩散模型：Q(x,y)=Q0exp[−(x−x o)2−(y−y0)22σ2] 。
2011全国大学生数学建模竞赛A题获奖论文——一篇
城市表层土壤重金属污染分析的数学模型摘要为研究城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式。
本文通过处理和分析已给数据，给出金属的空间分布说明污染程度和主要原因；建立数学模型确定污染源位置；最后收集其他信息讨论城市地质环境的演变模式。
问题一，利用matlab软件作出位置坐标x、y与八种总金属元素浓度的空间分布图；分析采集的重金属元素浓度所在区域的大致情形。
对采集的重金属元素浓度的数据进行分析，并计算单因子和多因子污染指数，根据土壤污染分级标准判断出不同重金属元素在各功能区的污染程度和各功能区的综合污染程度，其中工业区中铜是所有元素在不同功能区中污染程度最严重的，而工业区和交通区的综合污染程度是最严重的。
问题二，首先利用SAS软件对八种重金属元素在五个城区的含量进行主成分分析，得到八种重金属对各功能区的贡献率，可初步推断出工业生产、交通设施和生活垃圾造成重金属污染。
再利用SAS软件对各城区的重金属进行因子分析，进一步判断八种不同重金属污染的原因，如汞污染的原因为工业生产中三废的排放、交通运输业中汽油的燃烧和汽车轮胎磨损产生的粉尘等。
问题三，根据所给数据，分析重金属污染传播特征，即分别是介质的迁移运动、污染物的分散运动、污染物的累积与转化、污染物被环境介质吸收或吸附、污染物的沉淀，然后利用Matlab软件，采用多元纯二次二项式回归分析方法，分别得到每种重金属元素浓度与坐标的回归方程，并根据该方程利用多元函数求极值的方法确定出污染源的可能位置分别为：As（1878.2634，6003.7263，4.5846），Cd(970.5835，3946.7518，6.5891)，Cr（1235.1956，2658.3427，8.5402），Cu（138.4682，6223.4521，3.2461），Hg （1231.5782，2561.5483，5.2478），Ni（12234.2587，5865.1656，23.2461），Pb （2310.68914145.2674，3.2651），Zn（3015.43418642.2365 5.0543）；问题四，基于前三问，分析所建模型的优缺点。
2011全国数模竞赛A题
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）： 20111663 所属学校（请填写完整的全名）：南京邮电大学参赛队员 (打印并签名) ：1. 刘超舟2. 刘立亿3. 董睿指导教师或指导教师组负责人 (打印并签名)：日期： 2011 年 9 月12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：A题城市表层土壤重金属污染分析一．摘要二．问题重述三．符号约定四．基本假设五．建模与求解5.1问题（1）的建模与求解5.1.1问题分析根据采样的坐标，我们可以得到一个大概的区域分布散点图，如下图：图5.1.1 功能区采样点分布散点图图中各功能区得分布无规律，采样间隔比较大，很难建立各重金属测量值与样本采集点坐标的平滑函数，故采用等高线图来描述各种重金属在该城区的空间分布。
附件中给了采样点的三维坐标和各重金属元素的含量值，由于我们只需要考虑该区域的污染程度，因此不需要考虑坐标，只用知道采样地点位于某一区即可。
我们从采集样本的样本值中可以分析出每一区域的污染程度。
由于重金属分子的相对分子质量较大，取样时只对表层土壤0~10厘米深度进行取样，从收集的资料来看，可对取样值污染程度采用地积累指数分析法进行评级分析。
2011年高教杯A题数学建模论文
关键词：土壤重金属污染 数据标准化 插值分析 梯度场分析 污染评价
目
录
第一部分 问题重述…………………………………………………………() 第二部分 问题分析…………………………………………………………() 第三部分 模型的假设………………………………………………………() 第四部分 定义与符号说明…………………………………………………() 第五部分 模型的建立与求解………………………………………………() 1．数据处理及图像的给出 ……………………………………………() 2．题 1 的解答…………………………………………………………() 3．题 2 的解答……………………………………………………………() 4．模型的建立………………………………………………………() 5．模型的求解………………………………………………………() 第六部分 结果分析…………………………………………………………() 第七部分 模型的评价与改进……………………………………………（） 1．优点…………………………………………………………………(). 2．缺点…………………………………………………………（） 3．改进措施……………………………………………………(). 第八部分 参考资料…………………………………………… ……………() 第九部分 附录…………………………………………………… …………()
ቤተ መጻሕፍቲ ባይዱ 公园绿地区
山区
生活区
工业区
Hu 主干道路区
问题 1 的分析 问题 1 属于概率分布问题，运用插值拟合的方法分析。我们通过 matlab 拟合曲线, 得到概率模型为 p=λ e−αρ ，ρ = x 2 + y 2 （λ 为扩散系数，α 为实系数，p 为某元素的 密度） 。通过对对附件中所给数据分析可得到该城区内不同区域重金属的污染程度。
全国大学生数学建模大赛2011 A题
城市表层土壤重金属污染分析摘要通过分析城市城区土壤地质环境数据，选取采集样本为研究对象，建立模型综合评价城区污染状况，分析确立重金属污染源及其污染的传播方式。
模型一：采用模糊综合评价模型通过建立各重金属浓度的隶属度函数进而求出模糊关系矩阵和模糊权重矩阵，最后求出模糊评价向量，按照最大隶属度原则计算得出各区污染程度的等级，其中生活区、山区、公园绿地属于一级污染，工业区和主干道路区属于二级污染。
模型二：采用改进单因子污染指数模型计算各种重金属污染指数，再根据标准确定各区主要污染因子及其污染程度。
本文得到山区和公园绿地区各种金属污染程度都较低。
生活区主要污染因子Cd和Zn，工业区主要污染因子Cd、Cu、Hg、、Zn，主干道路区主要污染因子Cd和Zn，最后由这些不同区域的不同属性判断重金属污染的原因主要包括生活垃圾，工业排放，汽车尾气等。
模型三：利用地统计学中的半方差函数公式描述八种重金属元素分布在整个城市区域上的浓度变化，分析出这八个变量的空间变异方式，即确定变量的传播方式。
利用克里克法对重金属含量进行插值，模拟出各种重金属的分布扩散区域图，其中As,Cd,Hg 在东北部的污染源全都为点状以小范围的方式向周围辐射，Cd,Cr,Cu,Hg,Ni,Pb,Zn等元素都城市东南部以大范围的方式向四周面状和带状扩散。
污染源中心为22号、8号、20号等许多采样点所处区域。
最后再收集数据，对模型三进行改进，重金属元素通过大气、河流和沉降进入土壤，但都需要经过植物根系主导的土壤水循环，以此建立微分方程模型研究重金属元素的迁移规律，分析地质环境演变。
关键字模糊综合评价单因子污染指数半方差函数微分方程1 问题重述土壤是人类赖以生存的最基本的自然资源之一，也是人类环境的重要组成成分。
随着工业化、城市化进程的不断加快，废弃物排放、农业化肥使用量增加，土壤重金属污染越来越严重。
土壤重金属污染是由于人类活动使重金属在土壤中的累计量明显高于土壤环境质量标准或土壤环境背景值，致使土壤环境质量下降和生态环境恶化的现象。
2011数学建模A题参考答案
2011高教社杯全国大学生数学建模竞赛城市表层土壤重金属污染模型初探摘要土壤是人类赖以生存和发展的物质基础。
对我国这个人口大国而言，土壤问题显得尤为重要，但随着我国社会经济的飞速发展，土壤污染总体上却呈加剧趋势。
本文通过查阅众多文献资料，分析大量数据，探究城市表层土壤重金属污染状况。
通过使用Eviews6.0软件，计算数据的统计特征值，测算重金属在土壤表层的污染指数和变异系数，从而得出各重金属元素的空间分布特征。
然后，分别采用单因子污染指数法和多因子污染指数法，找出同一区域内的主要污染元素，并给出出主要污染原因。
由重金属的空间分布特征、统计特征值及污染原因，分析得出在8种重金属元素中，Cu、Hg、Zn 三种元素污染最为明显，故本文在第三部分重点分析三者的传播特征。
通过分析Matlab 画出的重金属元素污染的分布图，得出重金属元素污染以同心圆的方式向四周传播扩散，而且各点处的重金属元素污染物浓度与该点距污染源的距离成负趋势函数关系的结论。
为了避开由于多污染源叠加影响造成的误差，本文按照元素的污染高浓度聚集点，对整个城区进行分块处理，使得每一个污染块中仅含有一个污染源，从而简化模型，构造出污染物浓度与坐标值间的函数关系式(,,)，并通过Cu、Hg、Zn三种元素的Q f x y z数据对函数进行拟合得出参数值，由参数值得出的模型通过了统计中的T检验与F检验，模型拟合度合理，验证了上述模型的正确性及有效性。
为了研究地质环境的演变模式，本文提出引入定期收集同一采样点各金属元素浓度的信息，构建浓度与时间、坐标值之间的函数关系，从而缩小之前模型的误差，使其更符合实际规律。
关键词：统计特征值，污染指数，分块，拟合。
一问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
2011数学建模A题
Columns 1 through 7 详细数据见附录一列示部分如下： 0 96 137 190 245 298 248 96 0 163 94 149 202 166 137 163 0 143 198 251 201 190 94 143 0 57 110 94 245 149 198 57 0 53 151 298 202 251 110 53 0 204 248 166 201 94 151 204 0 328 232 281 138 97 150 80 372 276 325 182 127 86 124 308 256 261 184 241 294 90 397 301 350 207 152 179 149 453 357 406 263 208 179 205 505 409 458 315 260 207 275 490 394 443 300 245 192 242 596 500 549 406 351 298 348 698 602 651 508 453 400 450 723 627 676 533 478 425 475 599 503 552 409 354 301 351 420 324 373 230 273 326 172 706 610 659 516 461 408 458 750 654 703 560 505 452 502 690 594 643 500 445 392 442 698 602 651 508 453 400 450 574 478 527 384 329 332 326 669 573 622 479 424 371 421 639 543 592 449 394 431 391
由垃圾转运站垃圾转运量等情况统计表（南山）统计出所有小区垃圾总量为422.63104吨， 结合深圳市南山区垃圾转运站分布图，导入MATLAB可以确定各个小区及转运站的位置坐标， 及分布。建立模型如下：用spss将所有的转运站分为五类，会得到五个中心，五个中心中 离得最远的两个点的连线的中点作为圆心作为大小橱余设备的位置。位置模型如下图1，其 中一个圆表示聚合后的一类；如图紫色区域为大小设橱余备的位置
cumcm2011A
2011高教社杯全国大学生数学建模竞赛题目（请先阅读“全国大学生数学建模竞赛论文格式规范”）
A题城市表层土壤重金属污染分析
随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。
现对某城市城区土壤地质环境进行调查。
为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。
应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。
另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。
附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。
现要求你们通过数学建模来完成以下任务：
(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。
(2) 通过数据分析，说明重金属污染的主要原因。
(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。
(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？
附件1. 取样点位置及其所属功能区。
2011年数学建模A题国家一等奖
2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：2011 年09 月12 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要本文针对城市表层土壤受重金属污染问题进行综合分析。
首先运用Matlab软件求解8种主要重金属元素在城区空间分布，土壤受污染主要原因，传播特征以及为今后如何更好研究地质演变问题分别建立了相应的数学模型，并对其求解结果作出了分析。
针对问题1，根据各种污染物浓度在不同区域内分布的随机性，利用空间内插值法，以城市位置为平面“横纵向”，污染物浓度为“竖向”建立三维空间模型，得到这8种重金属元素的空间分布图，直观反映出污染物浓度的等值线，得到污染物浓度的分布规律和各种重金属元素的污染浓度范围区块。
然后通过建立污染负荷指数法模型算出各区的)(PLI为1.0602，(PLI为2.1573，山区) (PLI值，生活区)(PLI为1.8336，工业区)交通区)(PLI为1.5780；结果表明工业区的污染程度最高(PLI为1.9209，公园绿地区)为2级强污染，其他区为1级中等污染。
2011年全国数学建模A获奖论文
2011 高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的 资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规 则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从 A/B/C/D 中选择一项填写） ： 我们的参赛报名号为（如果赛区设置报名号的话） ： 所属学校（请填写完整的全名） ： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 文涛 北京大学A日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：2011 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：摘要随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响 日显突出。
本文主要对 Pb、Cd、Cr、Cu、Ni、Zn、As 和 Hg 这 8 大重金属在不同城区 的污染程度进行分析。
对于问题一，本文采用地统计学软件 A rcG IS ，建立8种主要重金属元素在该城区的 空间分布图，并从图中直观地看出8种主要重金属元素在不同区的污染情况。
为了更客 观、准确地对不同区域重金属的污染程度分析，采用目前国内外普遍采用的单因子指数 法和内梅罗综合指数法建立模型，算出不同区域内各种重金属的污染程度。
2011全国数模竞赛A题
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）： 20111663 所属学校（请填写完整的全名）：南京邮电大学参赛队员 (打印并签名) ：1. 刘超舟2. 刘立亿3. 董睿指导教师或指导教师组负责人 (打印并签名)：日期： 2011 年 9 月12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：A题城市表层土壤重金属污染分析一．摘要二．问题重述三．符号约定四．基本假设五．建模与求解5.1问题（1）的建模与求解5.1.1问题分析根据采样的坐标，我们可以得到一个大概的区域分布散点图，如下图：图5.1.1 功能区采样点分布散点图图中各功能区得分布无规律，采样间隔比较大，很难建立各重金属测量值与样本采集点坐标的平滑函数，故采用等高线图来描述各种重金属在该城区的空间分布。
附件中给了采样点的三维坐标和各重金属元素的含量值，由于我们只需要考虑该区域的污染程度，因此不需要考虑坐标，只用知道采样地点位于某一区即可。
我们从采集样本的样本值中可以分析出每一区域的污染程度。
由于重金属分子的相对分子质量较大，取样时只对表层土壤0~10厘米深度进行取样，从收集的资料来看，可对取样值污染程度采用地积累指数分析法进行评级分析。
全国大学生数学建模竞赛2011A题评阅要点（精选多篇）
全国大学生数学建模竞赛2011A题评阅要点（精选多篇）第一篇：全国大学生数学建模竞赛2011A题评阅要点2011高教社杯全国大学生数学建模竞赛A题评阅要点[说明]本要点仅供参考，各赛区评阅组应根据对题目的理解及学生的解答，自主地进行评阅。
本问题的数据来源于某城市对土壤环境的实地监测。
评阅时，应着重注意数学模型的建立、计算方法（或所选软件的程序语句）及选择该方法的理由。
(1)可用插值拟合的方法获得各重金属污染物浓度的空间分布。
再参考由背景值确定的阈值，定量分析城区各区域的污染程度。
由于空间数据是不规则的，较好的方法是用散乱数据插值，例如Kriging插值、Shepard插值等。
也可以用其他方法插值拟合，但应明确所使用的方法，并作出分析，不能只简单套用软件。
各个污染元素浓度的最大值与插值后浓度的最大值距离不会太远。
(2)分析污染产生的原因，必须有充分的数据分析以及明确的结论。
例如，可以根据各区域的污染浓度信息进行聚类，考察污染物出现的相关性，发现某些污染物结伴出现（如Cr与Ni，Cd与Pb的相关性较高），这与污染物产生的原因是密切相关的，由此可大致确定出产生这些污染的原因。
(3)本小题可以在不同的假设下建立相应的模型，但必须有合理的假设、建立明确的数学模型，并根据模型和所给的数据进行数值计算。
例如，由于雨水的作用是重金属在土壤表层中传播的主要原因之一，可以假设传播以对流形式为主，由此建立对流方程，并以给出的重金属污染物浓度数据作为初始值（实际上是终值），从而得到偏微分方程的定解问题。
类似于(1)，采用插值拟合的方法，可以得到地形高度函数。
利用特征线法，可以得到各区域在各个时间点上的重金属污染物浓度数据，从而可以得到各时间的污染范围，由此确定出污染源的位置。
(4)本问题只给出一个时间点上的数据，信息量明显不足，需要补充更多的信息。
如果学生考虑到多个时间点上的采样信息，给出更好的演化模式，应予以鼓励。
第二篇：全国大学生数学建模竞赛2011D题评阅要点2011高教社杯全国大学生数学建模竞赛D题评阅要点[说明]本要点仅供参考，各赛区评阅组应根据对题目的理解及学生的解答，自主地进行评阅。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：04505所属学校（请填写完整的全名）：德州学院参赛队员(打印并签名) ：1. 翟明元2. 张小凤3. 邹菲指导教师或指导教师组负责人(打印并签名)：赵学杰高秀莲日期： 2011 年 9 月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于系统综合评价的城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染问题，首先对各重金属元素进行分析，然后对各种重金属元素的基本数据进行统计分析及无量纲化处理，再对各金属元素进行相关性分析，最后针对各个问题建立模型并求解。
针对问题一，我们首先利用EXCEL 和 SPSS 统计软件对各金属元素的数据进行处理，再利用Matlab 软件绘制出该城区内8种重金属元素的空间分布图最后通过内梅罗污染模型：2/12max22⎪⎪⎭⎫ ⎝⎛+=P P P 平均综，其中平均P 为所有单项污染指数的平均值，m ax P 为土壤环境中针对问题二,我们首先利用EXCELL 软件画出8种元素在各个区内相对含量的柱状图，由图可以明显地看出各个区内各种元素的污染情况，然后再根据重金属元素污染来源及传播特征进行分析，可以得出工业区及生活区重金属的堆积和迁移是造成污染的主要原因，Cu 、Hg 、Zn 主要在工业区和交通区如公路、铁路等交通设施的两侧富集，随时间的推移，工业区、交通区的土壤重金属具有很强的叠加性，受人类活动的影响较大。
同时城市人口密度，土地利用率，机动车密度也是造成重金属污染的原因。
针对问题三，我们从两个方面考虑建模即以点为传染源和以线为传染源。
针对以点为传染源我们建立了两个模型：无约束优化模型()[]()[]()22y i y x i x m D -+-=，得到污染源的位置坐标()6782,5567；有衰减的扩散过程模型得位置坐标（8500,5500），模型为：u k zu c y u b x u a h u 2222222222-∂∂+∂∂+∂∂=∂∂， 针对以线为传染源我们建立了l c be u Y ∆-+=0模型，并通过线性拟合分析线性污染源的位置。
针对问题四，我们在已有信息的基础上，还应收集不同时间内的样点对应的浓度以及各污染源重金属的产生率。
根据高斯浓度模型建立高斯修正模型，得到浓度关于时间和空间的表达式ut e C C -⋅=0。
在本题求解过程中，我们所建立的模型与实际紧密联系，有很好的通用性和推广性。
但在求点污染源时，我们假设只有一个污染源，而实际上可能有多个点污染源，从而使得误差增大，或者使污染源的位置够不准确。
关键词 内梅罗污染模型 无量纲化 相关性 回归模型 高斯浓度模型一、问题重述俗话说：“一方水土养一方人”。
城市是人类活动最密集的地区，但在废物处理设施仍不发达的绝大多数地区 ,城市及其周边土壤依然发挥着重要的容纳和净化污染物的功能,在强烈的环境负荷冲击下,土壤的服务功能面临极大的威胁,换言之,土壤的缓冲净化功能将接近极限并有被超过的危险,因而将导致严重土壤污染的产生,而其结果将是长远和危险的。
随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
我们将城区分为生活区、工业区、山区、主干道路区及公园绿地区五个部分，分别进行土壤地质环境的调查，对城市环境质量做出评价，希望能有效控制重金属污染物的排放及扩散，制定相关措施保护好我们赖以生存的周边环境，根据题意，本文需要解决的问题有：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。
(2) 通过数据分析，说明重金属污染的主要原因。
(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。
(4)为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、模型假设⑴不考虑元素间的相互作用的影响⑵短期内重金属元素的物理、化学变化及迁移对周围环境影响不大 ⑶假设附录中所给8种重金属元素的背景值真实 ⑷不考虑历史沉积的重金属的影响三、符号说明i x 第i 种元素在第j 个采样点的浓度（8,...2,1=i ）；x 第i 种元素浓度的平均值（8,...2,1=i ）； i x ' 第i 种元素在第j 个采样点无量纲化后的数值（5,...2,1=j ）； i P 第j 个功能区重金属i 的单项污染指数（5,...2,1=j ）； i C 第j 个功能区重金属i 含量的实测值（5,...2,1=j ）；D 污染距离积；h 污染源位置与已知采样点的距离；()()()i yi x,给定采样点的坐标；四、数据处理4.1 对重元素的分析城市工业“三废”排放，金属采矿和冶炼，家庭燃煤，生活垃圾，汽车尾气排放都增加了城市土壤重金属的负荷。
重金属污染环境的主要有汞、铅、铬、锌镉、铜等。
其中汞的毒性最大，铬、铅、锌等也有相当大毒性。
此外还有砷，砷虽不属于金属.但它的毒性与重金属相似，因此归于重金属一类阐述，称为类金属。
目前对我国土壤污染比较普遍的重金属有汞、铬、砷。
根据该城区重金属污染的情况，下面对重金属在土壤污染中的来源及传播途径作简要介绍。
4.1.1砷元素该元素毒性很低，水体中含砷污染物主要来自砷和含砷金属矿的开采、冶炼，以及和砷化物为原料的玻璃、颜料、药物、纸张的生产都可产生含砷的废水，造成水体的砷污染。
砷及砷化物在水中会在水生物体内累积，但累积程度比其他重金属要低。
砷和砷化物，一般可通过水、大气和食物进入人体。
4.1.2镉元素当环境受到镉污染后，镉可在生物体内富集，通过食物链进入人体引起慢性中毒。
镉的主要污染源是电镀、采矿、冶炼、染料、电池和化学工业等排放的废水。
相当数量的镉通过废气、废水、废渣排入环境，造成污染。
镉对土壤的污染主要有气型和水型两种。
气型污染主要来自工业废气。
镉随废气扩散到工厂周围并自然沉降，蓄积于工厂周围的土壤中，可使土壤中的镉浓度达到40ppm。
水型污染主要是铅锌矿的选矿废水和有关工业（电镀、碱性电池等）废水排入地面水或渗入地下水引起。
4.1.3铬元素对水体污染的铬主要来源于电镀、制革、铝盐生产以及铬矿石开采所排放的废水。
是我国水体中一种普遍的污染物。
水体中铬污染主要是三价铬和六价铬，它们在水体中的迁移转化有一定的规律性。
4.1.4铜元素铜(Cu)及其化合物在环境中所造成的污染称为铜污染。
主要污染来源是铜锌矿的开采和冶炼、金属加工、机械制造、钢铁生产等。
冶炼排放的烟尘是大气铜污染的主要来源。
世界铜的年迁移量为:岩石风化20万吨,河流输送11万吨4.1.5汞元素汞是在常温下唯一呈液态的金属元素。
人类活动造成水体汞污染，主要来自氯碱、塑料、电池、电子等工业排放的废水。
由于天然本底情况下汞在大气、土壤和水体中均有分布，所以汞的迁移转化也在陆、水、空之间发生。
4.1.6镍元素镍污染是由镍及其化合物所引起的环境污染。
大部分煤含有微量镍,通过燃烧过程被释放出来,这是大气中镍的主要来源。
镍可以在土壤中富集。
土壤中的镍主要来源于岩石风化，大气降尘，灌溉用水（包括含镍废水），农田施肥，植物和动物残体的腐烂等。
全世界每年镍的迁移状况是：岩石风化量为320 000吨，河流输送量为19 000吨,开采量为560 000吨，矿物燃料燃烧排放5 600吨。
4.1.7铅元素铅对环境的污染，一是由冶炼、制造和使用铅制品的工矿企业，尤其是来自有色金属冶炼过程中所排出的含铅废水、废气和废渣造成的。
二是由汽车排出的含铅废气造成的，汽油中用四乙基铅作为抗爆剂（每公斤汽油用1～3克），在汽油燃烧过程中，铅便随汽车排出的废气进入大气，成为大气的主要铅污染源4.1.8锌元素锌在土壤中富集，会使植物体中也富集而导致食用这种植物的人和动物受害。
金属锌本身无毒，但在焙烧硫化锌矿石、熔锌、冶炼其他含有锌杂质的金属的过程中,以及在铸铜过程中产生的大量氧化锌等金属烟尘，对人有直接的危害。
其他如橡胶轮胎的磨损以及煤的燃烧也是大气锌污染的原因。
各种工业废水的排放是引起水体锌污染的主要原因。
4.2 对基本数据的分析用 EXCELL软件和 SPSS统计软件处理数据如表1所示：表14.3 元素浓度的无量纲化处理在利用SPSS 统计软件数据进行聚类分析的时候，因为单位不统一需要进行无量纲化处理，我们采用均值化方法，即每一个变量除以该变量的平均值，即xx x ii ='， (1) 标准化以后各变量的平均值都为1，标准差为原始变量的变异系数。
该方法在消除量纲和数量级影响的同时，保留了各变量取值差异程度上的信息，差异程度越大的变量对综合分析的影响也越大。
4.4 重金属元素间的相关性分析研究土壤中重金属的相关性可以推测重金属的来源是否相同,若重金属含量有显著的相关性,说明有相同来源的可能性较大,否则来源可能不止一个.我们用积差法来计算各重金属之间的相关系数，所谓积差法就是用两个变量的协方差与两个变量的标准差的乘积之比：()()()()()()2222222yyxxyx xy y y x x y y x x r yx xy---=----==∑∑∑σσσ表2 重金属元素间的相关系数A 组：Cr ，Ni ，CuB 组：Pb ，Cd ，Zn而对于As 、Hg 由相关系数表可见，其相关系数较小，我们认为相关参数小的元素间没有关系，所以将其各自单独一组。
下面我们建立回归模型图像验证它们之间的函数关系：Ni图1 Cr与Ni0.51.5Pb图2 Cd与Pb00.51 1.52Zn图3 Pb与Zn0.51.5Cu图4 Cr与Cu显然Cr，Ni和Cu及Pb，Cd和Zn显示属于适度空间相关性，反映区域因素（土壤母质）对其含量的影响较大，而As、Hg元素则属于低空间相关性，说明其受到人为因素（工业布局施肥灌溉和土地利用方式等）作用较强。
五、模型的分析、建立与求解5.1 问题一5.1.1重金属的空间分布由附件中所给的数据，我们考虑将各采样点的坐标和重金属的浓度建立对应关系，利用Matble软件画出等高线来体现该城区8种重金属的空间分布。