下载文档原格式
城市表层土壤重金属污染分析摘要对于问题一，首先采用克里格插值法根据附件中给的采样点的浓度数据对城区内每种重金属浓度值进行插值，进而绘制每种重金属的浓度的空间分布图，对此为了方便计算我们借助于Surfer软件进行绘图，然后我们采用单因子指数法和内梅罗综合指数法对各功能区的污染程度进行评价，但是这两种方法无法从自然异常中分离出人为异常，为了弥补其不足，采用地累积指数法[1]对土壤重金属污染做进一步评价，继而得到不同区域重金属的污染程度。
评价结果如表9所示。
对于问题二采用主成分分析法，旨在利用降维的思想，把多指标转化为少数几个综合指标，根据主成分分析法的一般步骤，首先对附件中给出的重金属浓度进行标准化处理，然后得到各金属之间的相关系数矩阵，求出相关系数矩阵的特征值和向量值，再得到因子的成分矩阵，确定出主成分的个数，计算出各因子的成分得分矩阵，最后通过算综合主成分中各个因子的权重系数得到污染性较大的因子，最后分析该污染主要原因。
结果如表14和评价结果。
问题三我们用对流-扩散偏微分方程来进行描述，对流扩散方程是描述粘性流体运动的非线性偏微分方程模型，我们将对流扩散方程进行简化，即变为二阶椭圆形偏微分方程，利用有限插值数值法估计出污染物的浓度分布，并与实际检测值相比较，偏差较小的即为污染源的位置及源强，即将问题转化为非线性最优化问题，结果如表所示。
对于问题四在分析问题三模型的优缺点基础上，为更好的研究城市地质环境的演变模式，将问题三的模型中二维对流扩散偏微分方程扩展到三维，这样就可以全面考虑影响模型准确性的各相关参数，通过收集不同的地理，天气条件下地质元素的空间分布信息，进而根据这些信息，建立优化模型，即三维方程的参数，在此基础上通过模拟仿真进行分析。
关键字：单因子指数法内梅罗综合指数法地累积指数法主成分分析法偏微分模型一问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
城市表层土壤重金属污染分析摘要城市表层土壤重金属污染严重影响人的生活和生产，在已知319个样本的坐标以及各重金属污染物浓度情况下，按照各问题的要求，我们建立相应的模型来对具体问题进行求解。
首先为了了解该城市的地形，我们使用了MATLAB软件绘制各分区分布地形图,根据区域化变量理论，运用克里金插值法对未采样点的区域化变量的值进行线性无偏最优估计，绘制出各污染元素的拟合分布图。
在分析不同区域受污染程度时，以题中所给的元素背景值以及国家土壤环境质量准则为基础制定污染等级评判标准，据此进行单因子指数评价以及模糊数学评价分析，得出每个区域内主要污染元素以及污染程度。
接下来，为了分析污染的原因，我们将单因子指数的结果与聚类分析算法的结果进行比较，得出生活区内铬、铜、镉、锌含量较高；工业区内的主要污染元素为镉、铬、汞、镍、铅、锌；山区各重金属含量均较低；交通区内镉、铜、汞污染较重；在公园绿地区内，镉、汞、铅、锌的污染较为严重。
更进一步，为了确定每一种重金属元素的来源区域，我们结合重金属污染物在土壤内易沉积，不易传播的特性以及当地海拔高度，通过主成分分析法和相关分析法计算出污染物的所在区域，并与surfer的绘图点进行比较，发现两者有很好的吻合度，找到了Cd的污染主要来自工业区和交通区；Cr污染主要来自交通区和生活区；Cu污染主要来自交通区和工业区；Hg污染主要来自交通区；Ni污染主要来自交通区；Pb污染主要来自生活区和工业区；Zn污染主要来自交通区，生活区，工业区。
对于污染物传播特性的分析，我们首先用SPSS对各区域内各污染因子进行正态检验，然后对于正态分布和非正态分布的因子我们各给出一种函数描述其传播特性。
对于问题四，我们提出建立城市污染物数据库，记录城市污染物随时间和空间的变化情况，有利于对未知点进行时间和空间的预测。
关键字：重金属污染、克里金插值法、多因素分析法、聚类分析、主成分分析法一、问题重述城市是人类社会不断发展的产物,是人类重要的生活环境,也是人类生产生活的主要场所. 随着经济和社会发展速度的不断加快,城市环境问题,尤其是重金属作为一个威胁人类生存和城市发展的重要的城市环境污染物,引起了学者、专家们的关注.随着工农业发展和城市化进程的加剧, 通过农业灌溉、工业排放、交通运输、市政建设和大气沉降等造成城市土壤重金属的污染日趋严重。
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A我们的参赛报名号为（如果赛区设置报名号的话）：A0504所属学校（请填写完整的全名）：集美大学诚毅学院参赛队员(打印并签名) ：1. 张建栋2. 蔡艺鑫3. 杨榕榕指导教师或指导教师组负责人(打印并签名)：日期： 2011 年 9月 11 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要：本文通过对测量数据的描述性统计分析，运用评价体系、单因子分析法、内梅罗综合污染指数法对各功能区的污染程度进行综合评价；其后，运用主成分分析法, 计算主成分贡献率及累计贡献率,主成分载荷确定各区污染物的主要构成成分，并且结合相关系数分析，根据传播过程中各元素的结伴组合，得出各重金属污染物的传播特征，找到污染的原因和污染源方位；在确定了污染源的区域后，为进一步找出污染源距离勘测点的距离，根据传质机理，运用等浓度线密集程度确定半径距离，得出勘测点与污染源的距离，以实际勘测点为圆心画半径，依据浓度差梯度方向，从而反推算出近相关的各元素污染源位置。
最后，为进一步完善模型，考虑到了污染源的流入速率、随着重力影响而流入土壤深处的金属元素比率、以离子状态扩散到周围的比率、被植物吸收等被大自然降解掉的比率，运用化学知识、微分方程，元素守恒综合考虑城市地质环境的演变模式。
2011高教社杯全国大学生数学建模竞赛城市表层土壤重金属污染分析摘要本文主要研究重金属对城市表层土壤污染的问题,我们根据题目所给定的一些数据和信息分析并建立了扩散传播模型、权重分配模型、对比模型和转换模型解决问题。
首先，我们利用Matlab 软件拟出该城区地势图（图1），根据所给数据绘出该地区的三维地势及采样点在其上的综合空间分布图。
之后将8种重金属的浓度等高线投影到该地区三维地形图曲面上，接着分别计算8种重金属在五个区域的平均值，立体图和平面图（图1附件）相结合便可得出8种重金属元素在该城区的空间分布。
其次，在确定该城区内不同区域重金属的污染程度时，我们运用两种方法进行解答。
先假设各重金属毒性及其它性质相同，运用公式ijij P C P ='求出各区域各金属相对于背景平均值的比值作为金属污染程度，再运用1ji ij j C C ==∑求出各区域重金属污染程度，并将各区进行比较。
之后，我们加上各重金属的毒性，对各重金属求出权数，再结合国标重金属污染等级和已知的各组数据来确定金属的污染程度。
由上述两种方法的对比，更准确地得出重金属对各区的影响程度。
即： 工业区>交通区>生活区>公园绿地区>山区 并根据第一个模型的数据来说明重金属污染的主要原因。
再次，对重金属污染物的传播特征进行了分析，判断出重金属污染物主要是通过大气、土壤和水流进行传播。
在分析之中，我们得出这三种状态的传播并不是孤立存在的，而是可以相互影响和叠加的，因此，我们分别建立三个传播模型，再对这三个传播模型进行了时间和空间上的拟合，得出重金属浓度最高的区域图，并结合各重金属的分布图（图6）来确定各污染源的位置。
最后，本题中只给出了重金属对土壤的污染，对于研究城市地质环境的演变模式，还需要搜集一些信息（图7）。
根据每种因素对地质环境的影响程度进行由定性到定量的转化。
建立同一地质时期地质环境中各因素的正影响和负影响的权重分配模型，再对这些权重进行验算和修正。
A题城市表层土壤重金属污染分析摘要一．问题重述随着人类活动对城市环境影响程度的逐步加深，研究城市地质环境的演变模式，日益成为人们关注的焦点。
由于不同的区域环境受人类活动影响的程度不同，所以通过对城市划分的五个功能区所采取样本点的分析，建立数学模型，并对以下问题进行求解：(1) 求解出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。
(2) 通过数据分析，说明重金属污染的主要原因。
(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。
(4) 分析所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二．问题分析2.1．问题重要性分析土壤重金属污染是指由于人类活动，土壤中的微量有害元素在土壤中的含量超过背景值，过量沉积而引起的含量过高，统称为土壤重金属污染。
土壤受重金属或类金属毒物污染后，常常通过农作物和水进入人体，造成人体中毒。
如铅、汞、镉、砷、铬、铊等重金属污染土壤后，都会对居民健康造成各种危害。
随着全球经济化的迅速发展，含重金属的污染物通过各种途径进入土壤，造成土壤严重污染。
土壤重金属污染可影响农作物产量和质量的下降，并可通过食物链危害人类的健康，也可以导致大气和水环境质量的进一步恶化。
而且随着今后中国经济社会的发展如生态退耕、农业结构调整及自然灾害损毁等，土壤资源将进一步减少。
因而如何有效地得出重金属元素在城区的分布情况，并及时控制,治理土壤重金属的污染，改良土壤质量,将成为生态环境保护工作中十分重要的一项内容。
2.2. 问题思路分析三．模型假设（1）忽略海拔对该城区重金属元素的空间分布的影响；（2）假设该城区表层土壤中元素的背景值可以作为土壤质量评价标准；四．符号说明五．模型的建立及求解5.11. 模型准备：由于每个样本所含的多种化学元素的浓度数据是按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS 记录采样点的位置。
2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则 .我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等与队外的任何人(包括指导教师研究、讨论与赛题有关的问题。
我们知道, 抄袭别人的成果是违反竞赛规则的 , 如果引用别人的成果或其他公开的资料(包括网上查到的资料 ,必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们将受到严肃处理。
我们参赛选择的题号是(从 A/B/C/D中选择一项填写 : A 我们的参赛报名号为(如果赛区设置报名号的话 :所属学校(请填写完整的全名 :重庆交通大学参赛队员 (打印并签名 :1. 陈训教2. 范雷3. 陈芮指导教师或指导教师组负责人 (打印并签名 :胡小虎日期:2011 年 9 月 12日赛区评阅编号(由赛区组委会评阅前进行编号:2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号(由赛区组委会评阅前进行编号 :全国统一编号(由赛区组委会送交全国前编号 :全国评阅编号(由全国组委会评阅前进行编号 :城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染做出了详细的分析, 对于本题中所提出的问题一,我们利用 MATLAB 软件对所给的数值进行空间作图,然后分别作出了八种重金属元素的空间分布特征, 然后, 我们利用综合指数(内梅罗指数评价的方法, 对五个区域进行了综合评价, 得出结果令人满意。
对于问题二, 我们根据第一问和题目所给的数据进行综合分析, 得出了重金属污染的主要原因来自于交通区含铅为主的大量排放, 和工业区污水的大量排放等等。
对于问题三, 我们通过对问题一中的八张重金属元素空间分布的图可以看出, 发现大多数金属都呈中心发散性传播, 同时经过分析, 我们发现, 如果考虑大气传播和固态传播, 很难得出结论, 在交通区, 由于是汽车尾气造成的传播, 发现重金属的传播无规律可循等,所以,我们考虑液态形式的传播,以针对地表水污染物的物理运动过程, 以偏微分方程为建模基础, 通过和假设和模型参数的估计, 得出了可能污染源位置, 最后, 我们对模型进行了稳定性检验即灵敏性分析和拟合检验, 发现在参数变化在 10%左右,模型的稳定性良好。
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从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（低）。
题目： A参赛队员：队员1蔡新星队员2古振炎队员3黄祥振指导教师：指导教练组单位：赣南师学院承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括、电子、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：赣南师学院参赛队员 (打印并签名) ：1. 蔡新星2. 古振炎3. 黄祥振指导教师或指导教师组负责人 (打印并签名)：指导教练组日期： 2011 年 9月 12 日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析【摘要】在经济快速发展的今天，人类活动对城市环境质量的影响日显突出。
研究人类活动影响下城市地质环境的演变模式显得尤为重要，本论文所建模型即为了研究某城区八种重要重金属形成污染作出评价并探究其成因。
对于问题一，本论文先用matlab插值法画出8种元素的空间分布图，对元素污染程度做出定性分析，然后采用单因子指数法，梅罗综合指数法和地积累指数法（Muller指数法）对土壤进行重金属污染评价，再对不同区域重金属的污染程度做出定量分析。
通过对三个模型的综合分析得出以下结论（没参加污染程度比较的元素即认为在此功能区无该金属污染）：1.生活区中各元素污染程度为：As < Cr < Pb< Cd < Hg < Zn< Cu。
2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：重庆交通大学参赛队员 (打印并签名) ：1. 陈训教2. 范雷3. 陈芮指导教师或指导教师组负责人 (打印并签名)：胡小虎日期：2011 年9 月 12日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要本文针对城市表层土壤重金属污染做出了详细的分析，对于本题中所提出的问题一，我们利用MATLAB软件对所给的数值进行空间作图，然后分别作出了八种重金属元素的空间分布特征，然后，我们利用综合指数（内梅罗指数）评价的方法，对五个区域进行了综合评价，得出结果令人满意。
对于问题二，我们根据第一问和题目所给的数据进行综合分析，得出了重金属污染的主要原因来自于交通区含铅为主的大量排放，和工业区污水的大量排放等等。
对于问题三，我们通过对问题一中的八张重金属元素空间分布的图可以看出，发现大多数金属都呈中心发散性传播，同时经过分析，我们发现，如果考虑大气传播和固态传播，很难得出结论，在交通区，由于是汽车尾气造成的传播，发现重金属的传播无规律可循等，所以，我们考虑液态形式的传播，以针对地表水污染物的物理运动过程，以偏微分方程为建模基础，通过和假设和模型参数的估计，得出了可能污染源位置，最后，我们对模型进行了稳定性检验即灵敏性分析和拟合检验，发现在参数变化在10%左右，模型的稳定性良好。
最后我们全面分析了模型的优缺点，，最后可以用MATLAB软件得出相应的结果。
为更好地研究城市地质环境的演变模式，测定污染源范围还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，地下水流动方向以及每年的生物降解量，降雨量对重金属元素扩散的影响。
一但有污染证据，我们可以在该污染源附近沿地下水流动方向设定更多采样点，由此，我们可以构造一个三维公式来计算污染物质浓度的浮动就可以模拟三维空间内的重金属分布影响。
关键字：表层土壤重金属污染 MATLAB 内梅罗指数偏微分方程稳定性检验灵敏性分析地质演变生物降解量一、问题重述随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点，对城区重金属污染分析以及传播越来越有其必要性。
对于本题中所提出的问题一，我们利用MATLAB软件对所给的数值进行空间作图，分别作出了八种重金属元素的空间分布特征图，然后，我们运用综合指数（内梅罗指数）评价的方法，对五个区域进行了综合评价，得出结果令人满意。
对于问题二，我们根据第一问和题目所给的数据进行综合分析，得出了重金属污染的主要原因来自于交通区大量排放的含铅为主污染物，和工业区污水的大量排放等等。
对问题三，通过分析第一问中所给定各种元素空间分布规律，和查阅大量资料，我们发现，如果考虑大气传播和固态传播，很难得出结论，在交通区，由于是汽车尾气造成的传播，发现重金属的传播无规律可循等，所以，我们考虑液态形式的传播，以针对地表水污染物的物理运动过程，以偏微分方程为建模基础，通过和假设和模型参数的估计，得出了可能污染源位置，最后，我们对模型进行了稳定性检验即灵敏性分析，发现在参数变化在10%左右，模型的稳定性良好。
为更好地研究城市地质环境的演变模式，测定污染源范围还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，地下水流动方向以及每年的生物降解量，降雨量对重金属元素扩散的影响。
一旦有污染证据，我们可以在该污染源附近沿地下水流动方向设定更多采样点，由此，我们可以构造一个三维公式来计算污染物质浓度的浮动就可以模拟三维空间内的重金属分布影响。
二、问题分析针对本题所提出的问题一，我们利用题给的数据运用MATLAB软件作出相应的八中重金属元素空间分布图，同时考虑到分析和评价城区累不同金属的污染度，我们根据数据应用内梅罗指数综合评价城区内不同区域的重金属污染度。
得到的结果令人满意。
针对问题二，通过对第一问和原始数据的分析，我们得出了重金属污染的主要原因还是交通区汽车含铅气体的大量排放，工业区不达标的污水大量排放等等。
对问题三，通过分析第一问中所给定各种元素空间分布规律，和查阅大量资料，我们发现，如果考虑大气传播和固态传播，很难得出结论，在交通区，由于是汽车尾气造成的传播，发现重金属的传播无规律可循等，所以，我们考虑液态形式的传播，以针对地表水污染物的物理运动过程，以偏微分方程为建模基础，通过和假设和模型参数的估计，得出了可能污染源位置，最后，我们对模型进行了稳定性检验即灵敏性分析，发现在参数变化在10%左右，模型的稳定性良好。
证明了我们模型的正确性。
对于问题四，三、模型假设假设一：问题中附件给出的原始数据真实，有效。
假设二：（0～10cm）的地表是由土和沙砾所构成。
假设三：所有的土和蓄水层的性质在浸透到水层和未浸透水层都是均匀的。
假设四：稳定、均匀内的水流只发生在整个未浸透水层的垂直方向，以及在浸透水层水平（纵向）平面中。
沿地表水流动方向。
假设五：物理过程起着重要作用，在此模型中只考虑物理过程（扩散），不考虑在其过程中的化学反应。
假设六：所有的污染源均为点源。
四、符号说明i P ：单项内梅罗指标数iave P ：i P 的平均值max i P ：iP 的最大值 P ：内梅罗综合指数I C ：污染物实测值aX ：污染金属背景值的最小值 cX ：污染金属背景值的平均值 p X ：污染金属背景值最大值l a ：水平扩散系数（m ）t a ：垂直扩散系数（m ）C ：污染物质的浓度（liter mg） bC ：背景浓度（liter mg ） 0C :污染物源处的浓度（liter mg） D ：渗透系数（m 2）K ：水力传导性（2ft day gal ）I ：水力梯度m ：污染物质的排放率（day mg ）n ：有效多孔性q :污染物质的排放率（day mg ）d R ：阻滞因子S :复合参数0t ：污染开始的时间（yr ）d V ：地下水的流速（day ft ）W ：hantush 函数),(00y x ：污染源的坐标五、模型的建立与求解5.1.1问题一的模型建立从附件-1中得到的数据，运用MATLAB软件，绘制了下面的三维地表图形，图形中的黑点就是取样的地点。
（z代表的是海拔）。
下面建立该城区内不同区域重金属的污染程度，为了更好地评价该城区不同区域的重金属污染程度，我们引入目前比较广泛应用的评价指标——内梅罗指数，进行综合评定，单项污染指数法能够比较直观地反映环境中各项污染指标的情况；内梅罗综合指数法不仅考虑到了所有评价因子单项污染程度的平均水平，而且还考虑到了最大污染指数，因此能够更为科学、综合的反映评价区域内总体环境质量状况。
评价公式如下：P：单项指标数i⎪⎪⎩⎪⎪⎨⎧>--+=<<--+=<<--+=≤=pi c p p i i p i c c p a i i c i a a c a i i a i ii i X C X X X C P X C X X X X C P X C X X X X C P X C X C P )/()(3)/()(2)/()(1/（1）其中：i C 为污染物实测值，p c a X X X ,,是依据附件给出的平均值，偏差，以及范围所定出的界限值。
综合污染指数（P ）采用内梅罗综合指数法：222max iavei P P P +=（2）其中：max i P 为单项指标数i P 的最大值，iave P 为单项指标数i P 的平均值。
给定污染综合指数等级的划分评定表，见表一：5.1.2问题一的模型求解依据原始数据和公式（1）（2），利用MATLAB 编程求得八种重金属元素的不同区域的单项内梅罗指数，见表二：种重金属元素的P 值的平均值，最终评价出污染等级，其果见表二，最后得出五大区域的污染等级见（表三）5.2分析说明重金属污染的主要原因5.2．1 该市表层土壤重金属含量基本状况分析将实际测得的该市五大区域土壤样品重金属含量与自然区土壤重金属含量背景值进行比较发现，五大区8大重金属含量平均值均高于背景值（表4）并且8大金属含量大部分样点超过自然区表层土壤均值，占总样点数的百分比最小为67.71%，最大为88.09%。
这说明表土层8大重金属均有外源物质的进入，并有了一定的积累。
表四5.2．2 表层土壤重金属污染总体评价根据五大区域8种重金属内梅罗指数（表二）求出各元素的单项内梅罗平均指数（见表五），再对照土壤的等级评价标准（表一）对表层土壤污染进行评价，评价得出 Cd、Cu、Hg、Zn这四种重金属指数达到了污染等级指数，污染级污染样点分别占79.94%、88.08%、67.71%、78.68%，均属于大面积轻度污染，其中Pb临近警戒线。
这就可以看出，Cd、Cu、Hg、Zn是构成污染的的主要原因，需要加强控制并采取相应措施进行降排，而Pb需要预防性控制，其余的As、Cr、Ni均较为安全。
5.2．3地区对四种轻度污染元素的影响经过对样点金属含量的详细分析，可以看出各种重金属在不同的区域的分不存在着很大的差异性。
于是针对达到污染等级的四种重金属在各个区域的单项内梅罗指数作出其分布图（图9）。
通过对各金属污染指数的分析得到这四中重金属元素的主要污染来源区为工业区和交通区。
工业区内和交通区内的采样点的单项污染指数明显高于其他地区。
工业排放和增加了重金属的含量，通过图9知道在工业区内的四种污染等级的重金属元素污染指数均超过了2且为两种以上的重金属复合污染，其四种重金属元素含量均达到了较高的水平。
交通区域污染指数也较高，也为多种重金属复合污染，这说明该市的交通表较发达，车辆排污水较高，大量排放的尾气，和车胎摩擦所带来的重金属是造成交通区污染的主要原因。
