2011数学建模土壤污染问题.doc

城市表层土壤重金属污染分析摘要随着我国城市化进程发展，城市土壤重金属污染日益成为人们关注的焦点。本文通过对城市表层土壤重金属污染的详细分析，探究城市表层土壤重金属污染的程度，研究其污染的原因，寻找污染源，运用数学的方法来解决这一问题。对于问题一，本文将附件中给出的数据，通过matlab采用插值法建立三维模型，利用Excel建立数据图表，通过模型图我们可以清楚的看到各种元素不同的空间分布。然后通过单因子污染指数法和内梅罗综合指数法，分析了不同元素对不同功能区的污染程度，算出不同区域内各种重金属元素的污染程度。对于问题二，本文对文献中的重金属污染的主要原因进行了分析，结合问题一中得出的数据，对每种重金属的空间位置及污染情况分别进行了Excel画图处理，得出了污染的主要原因，然后利用层次分析法对结论进行了正确性的分析与讨论，最终证实结论的正确性，即该城区的重金属污染主要原因在于工业释放及交通排放。对于问题三，本文根据前两问的结论分析重金属的传播特征，建立了偏微分方程模型，将各个元素的浓度与坐标带人方程求出各参量的值，再结合已利用插值法画出的三维立体图从而确定污染源的位置，并对已经得出的数据进行了误差性的分析误差分析。对于问题四，我们仔细分析了模型的优缺点。为更好地研究城市地质环境的演变模式，还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，以及每年的生物降解量，降雨量对重金属元素扩散的影响，空气污染也应该考虑进去。有了这些数据以后建立因子分析法，回归分析，曲线拟合等模型解决问题。关键词位置不对单因子指数法、内梅洛综合污染指数法、层次分析法、模糊综合评判、灰色理论一 问题重述按照重述可以简洁，但研究八种金属的污染问题要提出来功能划分，城区一般可分为生活区、工业区、山区、主要干道区及公园绿地区等，分别记为1类区、2类区、5类区，不同的区域环境受人类活动影响的程度不同。现要求你们通过数学建模来完成一下任务： 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 通过数据分析，说明重金属污染的主要原因。 分析重金属污染的传播特征，由此建立模型，确定污染源的位置。 分析你所建模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二 模型假设 假设污染源采集过程中不丢失，能正确反映具体情况； 假设调查过程中，不发生自然灾害等恶劣天气，不考虑风力等其他因素对重金属传播过程的影响，污染源浓度不变； 假设所选择的采样点能正确反映该城区污染的浓度； 假设各个采集点进行取样时间差对问题研究没有影响。三 符号说明Pi-污染物i的污染指数；Ci-污染物实测值（质量分数）；Si-污染物评价标准；P综-综合污染指数；(Ci/Si)max-土壤重金属元素中污染指数最大值；(Ci/Si)ave-为土壤各污染指数的平均值；Ni-污染物累计指数；Mi-各元素背景浓度平均值四 问题分析 随着社会生产力的不断提高和科学技术的不断进步，人类在开发环境、利用环境创造和丰富物质文明的同时，也在污染和破坏着自身的生活环境。现根据该城区的土壤地质环境的调查结果，分析问题。问题一的分析：要描述8种金属在该城区的空间分布，根据附件提供的数据，结合数据特点，运用matlab软件编程，分别建立8种金属在该城区的空间分布三维图，就可以形象地表达出主要重金属在该城区的空间分布。 根据该城区不同地区重金属的污染程度的分析，在现在的社会生活中，本文结合土壤重金属污染的特点，选用单因子分析法和内梅罗综合污染评价法来评价该区不同地区的污染程度。再根据我国土壤环境质量标准，以背景值作为参比值，得到5个区的污染指数，再运用excel软件绘制条形图来形象展现各个城区的不同金属的污染指数，加以分析其污染程度。问题二的分析：要分析该城区的污染原因，首先要知道土壤中的重金属来源多种多样，但人类工农业生产活动中也会产生一定的重金属。我们了解8种金属污染的主要来源，然后分析各个区的各个元素间的相关性，同时考虑各元素主要存在于哪个功能区，判断相关系数高的元素是否由该功能区的某种污染物造成，从而确定重金属的污染原因。问题三的分析：根据重金属污染物的传播特征，建立模型确定污染源。首先我们了解了这8种重金属的主要传播特征，但是由于这些金属传播特征各异，无法建立一个统一的模型确定出污染源的位置。但考虑到大部分金属颗粒的状态可以通过大气传播，于是依据烟羽模型，通过测量样本点的浓度和位置，推出污染源所在的轨迹，将多个样本反推出的轨迹相交，得以确定污染源的位置。问题四的分析：我们仔细分析了模型的优缺点。为更好地研究城市地质环境的演变模式，还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，以及每年的生物降解量，降雨量对重金属元素扩散的影响，空气污染也应该考虑进去。有了这些数据以后建立因子分析法，回归分析，曲线拟合等模型解决问题。五 模型的建立与求解5.1问题一金属的空间分布 对于要展现8种重金属元素的空间分布，本文将附表中的数据，运用matlab软件编程（见附录），建立三维立体图，首先描绘了各个取样点的空间位置，其中不同的颜色表示不同的区域，如图5.1.1图5.1.1采样点的位置空间分布 为了能清晰的在空间分布图中展现8中金属元素的在城区的分布情况，我们运用matlab软件编程（见附录）分别建立每种元素在空间分布的三维图，并用不同的颜色表示重金属的不同浓度，既体现每种元素在城区的分布，又体现出各个地区的重金属的污染程度。如图.9 图5.1.2As元素的空间分布 图5.1.3Cd元素的空间分布 图5.1.4Cr元素的空间分布 图5.1.5Cu元素的空间分布 图5.1.6Hg元素的空间分布 图5.1.7Ni元素的空间分布 图5.1.8Pb元素的空间分布 图5.1.9Zn元素的空间分布为了能够通过浓度数据直观的展现在不同的功能区内不同的重金属的浓度分布情况将附表中的数据通过Excel作图，如图5.1.10图5.1.10在整个城区的均值浓度分布为了能够论述各元素的空间分布情况，我们对各元素在整个城区浓度均值柱状图进行了客观的分析，如表1表1 8种主要金属元素空间分布分析As (砷)该元素在城区各个功能区内都分布较少Cd该元素主要分布在工业区及城市内的主干道附近，是由于工业废物及汽车尾气的排放造成的Cr该元素在城区的各个功能区都分布较为均匀，或许由于该地的气候等多方面的影响。Cu该元素主要分布在工业区，山区较少，人类的生产生活的影响工业废物排放都能影响其的浓度Hg该元素主要分布在工业区及城市内的主干道附近，是由于工业废物及汽车尾气的排放造成的Ni该元素在城区的各个功能区都分布较为均匀，或许由于该地的气候等多方面的影响。Pd该元素在城区的各个功能区都分布较为均匀，或许由于该地的气候等多方面的影响，但是山区较少人类活动对其有显著影响。Zn该元素在城区的各个功能区都分布都较多，或许由于该地的气候，水源等多方面的影响。5.2该城区各个地区重金属污染程度的评价 综合各种污染评价方法的适用范围、评价目的、优缺点，我们选用单因子指数法和内梅罗综合污染指数法两种方法来评价重金属的污染程度。（2）单因子指数法该方法是对单一的变化因子进行数据的处理得到相应的指数，再与指数区间比较得到相应的等级评价 式中：Pi为污染物i的污染指数；Ci为污染物实测值（质量分数）；Si为污染物评价标准；i代表某种污染物。当Pi 1时，表示土壤未受污染；1Pi2时，表示土壤已受轻微污染；2Pi3时，表示土壤已受中度污染；Pi3时，表示土壤已受严重污染。表2各元素在功能区污染指数值单因子污染指数一类区二类区三类区四类区五类区As0.418030.483430.269610.380.42Cd1.449811.965560.761.811.40271Cr0.766870.590.430.646160.48484Cu1.413.640.494781.786260.86262Hg0.624.282370.273042.987620.76661Ni0.458560.495290.386350.440320.38224Pb1.974472.658311.044461.821141.73453Zn2.370092.779280.732942.440991.54最大单因子污染指数2.370094.282371.044462.987621.73453平均单因子污染指数1.18372.112690.549461.53870.9492根据表2可得：一类区生活区未受As、Cr、Hg的污染，受Cd、Cu、Pb轻微污染，受Zn的中度污染;二类区工业区未受As、Cr、Ni得污染，受Cd的轻度污染，受Pb、Zn的的中度污染，受Cu、Hg的严重污染;三类区山区未受As、Cd、Cr、Cu、Hg、Ni、Zn的污染，受Pb的轻微污染。四类区交通区未受As、Cr、Ni的污染，受Cd、Cu、Pb的轻微污染，受Hg、Zn的中度污染；五类区公园绿地区未受As、Cr、Cu、Hg、Ni的污染，受Cd、Pb、Zn的轻微污染。我们有根据表1的数据做了各个区域不同重金属的单因子指数条形图，形象地体现出各个地区受每种重金属的污染程度。如图5.2.1图5.2.1各个区域不同重金属的单因子指数由表2图5.2.1可知：功能区一各元素污染程度 ：Zn PbCd Cu Hg Cr Ni As功能区二各元素污染程度: Cd Cu Zn PbHgCr Ni As 功能区三各元素污染程度：Pb CdZn CuCrNi HgAs 功能区四各元素污染程度：HgZn Pb Cd Cu C r Ni As 功能区五各元素污染程度：PbZn CdCu Hg Cr As Ni （2）内梅罗综合污染指数法内梅罗指数是一种兼顾极值或称突出最大值的计权型多因子环境质量指数。内梅罗指数特别考虑了污染最严重的因子，内梅罗环境质量指数在加权过程中避免了权系数中主观因素的影响，是目前仍然应用较多的一种环境质量指数。内梅罗综合污染指数法计算公式： 式中：P综为综合污染指数；(Ci/Si)max为土壤重金属元素中污染指数最大值；(Ci/Si)ave为土壤各污染指数的平均值。 将附表中的数据通过Excel处理我们得到了各个区域不同重金属浓度均值，如表3表3 各个区域不同重金属浓度均值As(g/g)Cd(ng/g)Cr(g/g)Cu(g/g)Hg(ng/g)Ni(g/g)Pb(g/g)Zn(g/g)1区均值6.2704545289.9613669.01840949.40318293.04068218.34227369.106364237.008642区均值7.2513889393.1111153.409167127.53583642.3552819.81166793.040833277.92753区均值4.0440909152.319738.95969717.31727340.95606115.45378836.55590973.2942424区均值5.7080435360.0144958.05391362.214928446.8225417.61710163.534203242.854865区均值6.2637143280.5428643.63630.191714114.9917115.28971460.708571154.24229 8种主要金属元素的背景值（见附表）元素平均值标准偏差范围As (g/g)5.4Cd (ng/g)1303070190Cr (g/g)3191349Cu (g/g)20.4Hg (ng/g)3581951Ni (g/g)19.9Pb (g/g)3161943Zn (g/g)69144197将已经求出的5大功能区8种重金属的环境质量指数用内梅罗综合指数法的表达式进行处理，即可求出五大功能区土壤综合污染指数。结果如表5：表5 五大功能区土壤综合污染指数功能区12345P综合3.04245.89271.272413.25782.9818表6 污染等级评价等级一P综=0.7清洁（安全）等级二0.7P综=1.0尚清洁（警戒线）等级三1.0P综=2.0轻度污染等级四2.0P综=3.0中度污染等级五P综3.0重度污染由五大功能区土壤综合污染指数与污染等级对比可知：1:工业区，交通区和生活区属于重度污染且有交通区工业区生活区2:公园绿地区属于中度污染3：山区属于轻度污染5.3问题二重金属污染的主要原因（1）土壤重金属污染原因主要有以下几个方面：工业三废引起的重金属污染由于部分矿产开发中选矿、冶炼工艺水平落后，个别矿区没有环保治理设备，废水、废气排放出的大量废弃物未经处理直接投放环境，而其中的重金属随着自然的沉降、雨水的淋溶等途径进入土壤，进入正常循环的生态系统，造成重金属污染严重危害人们的生产生活。化肥农药的过度使用 重金属元素是过量肥料中最多的污染物质，化肥中的过磷酸钙和磷矿粉中含有微量的As、Cd重金属元素（WILLIAMS C H，1973）。含铅及有机汞的农药发挥作用的同时也为土壤重金属污染埋下了祸根，造成土壤的胶质结构改变，营养流失，对农作物的产量及品质都造成极大的不良影响。目前的饲料添加剂中也常含有高含量的Cu和Zn（夏家淇，1996），这使得有机肥料中的Cu、Zn含量也明显增加并随着肥料施入农田。汽车尾气的排放以公路、铁路为中心成条带状分布的重金属污染土壤主要是由于汽车尾气的排放、汽车轮胎磨损产生的大量含重金属的有害气体和粉尘的沉降所引起的，污染元素中主要为Pb、Cu、Zn等元素（李波，2005）。这些物质随风飘落，进入土壤中引起重金属污染。实验证明，道路两旁土壤中重金属的污染比较严重，并随着离公路距离的由近到远，土壤的污染程度渐轻。 （2）土壤累积污染指数法为了能够展现污染物累积度进而分析污染原因我们利用土壤累积污染指数计算公式：（Ni为污染物累积指数，Ci为实测浓度平均值，Mi为各元素背景浓度平均值）：根据附件中的数据，运用土壤累积污染指数计算公式，计算出土壤污染累积指数表，见表7表7 土壤污染累积指数表全区平均值背景平均值标准偏差范围污染物累积指数As (g/g)5.905.41.64Cd (ng/g)295.4913030701902.27Cr (g/g)52.6431913491.70Cu (g/g)57.3920.44.35Hg (ng/g)267.9035819517.65Ni (g/g)17.3019.91.41Pb (g/g)64.6331619432.08Zn (g/g)197.31691441972.86图5.3.1元素在不同功能区的浓度均值及背景值根据表7给出的污染累积指数，及图5.2.2可以推断，城区的重金属污染，与人类的生产、生活有密切关系，通过上图将全区金属浓度平均值与背景值做比较，可以初步得出此城区的重金属污染物主要为Cd、Cu、Hg、Pb、Zn。，该城区的污染可能为工业排放和交通尾气排放。基于对问题一中所得数据的分析，我们绘制出了各元素在各区的单因子污染指数柱状图，见图5.2.3 图5.3.2各元素在不同功能区的污染指数 图5.3.3区单因子污染指数分布图据图5.2.3分析，该城区不同功能区土壤中重金属污染贡献不同，从图中可以清晰的得出，Pb的污染范围最广，基本对各个区都有影响，Cd、Zn影响4个区， Cu影响3个区，Hg影响2个区，As、Cr、Ni对各个区基本没有影响。由问题一得出的结论绘制图如表8：表8 污染贡献值 污染贡献元素主要污染贡献元素一类区Cd、Cu、Pb、ZnZn二类区Cd、Pb、Cu、Hg、ZnHg三类区PbPb四类区Cd、Pb、Cu、Hg、ZnHg五类区Cd、Pb、ZnPb据图5.2.4和表8分析，该城区不同功能区土壤中重金属污染贡献具有一定的差别，在一类区中，Zn、Pb的污染贡献最大，Cd、Cu居其次，根据查阅到的文献资料，Cd、Cu、Pb、Zn来源主要为工业废水的排放和农药的使用，由于一类区属于生活区，故可判断其污染的主要原因可能为生产活动中化肥施用过量，或是城市垃圾的焚烧，燃煤等。在二类区中，Hg、Cu、Zn、Pb的污染贡献最大，Cd居次，Hg、Cu、Zn、Pb的主要来源于工业废水，污泥，农药等，由于二类区为工业区，可基本判断为工业污染，主要原因可能为颜料厂，冶金厂，电镀厂等工业废水、废渣处理不合格，乱排乱放。在三类区中，仅Pb具有一定的污染贡献，且其含量很低，与背景值最为接近，鉴于三类区为山区，其Pb污染可能为气流运动，沉降累计。在四类区中，Cd、Pb、Cu、Hg、Zn的污染贡献普遍很大，由于四类区为主干道路区，其污染原因最为可能是汽车尾气的排放，以及汽车轮胎的磨损。在五类区中，Cd、Pb、Zn的污染贡献大，五类区为公园绿地区，其污染原因最为可能是通过气流运动，降水等沉降积累。综合以上分析，该城区的重金属污染主要原因在于工业释放及交通排放。5.4问题三 确定污染源的位置（1）通过一些资料的查询我们了解到重金属的主要的传播特征有：重金属的沉降大气中重金属主要来源于工业生产、汽车尾气排放、汽车轮胎磨损产生的大量含有重金属的有害气体和粉尘。它们主要分布在工矿的周围及公路、铁路两侧。大气中大多数重金属通过自然沉降和雨淋沉降进入土壤圈。经过自然沉降和雨淋沉降进入到土壤中的重金属污染，主要以工矿烟囱、废物堆和公路为中心向四周及两侧扩散：由城市郊区农区，污染随着距城市距离的加大而降低，城市的郊区污染较为严重；此外，还与城市的人口密度、城市土地利用率、机动车密度呈正相关；重工业越发达，污染相对就越严重。农药、化肥和塑料薄膜的使用使用含有Pb、Cd、Hg、As等的农药和不合理地施用化肥都可以导致土壤中重金属的污染。一般过磷酸盐中含有高量的重金属Hg、Cd、As、Zn、Pb，磷肥次之，氮肥和钾肥含量较低，但氮肥中Pb、As和Cd含量较高17。农用薄膜生产中应用到的热稳定剂中含有Cd、Pb，在大量使用塑料大棚和地膜过程中都可以造成土壤重金属污染。污泥施肥污泥中含有大量的有机质和氮、磷、钾等营养元素，但同时污泥中也含有大量的重金属元素，随着市政污水处理产生的大量污泥被施加于农田，农田中的重金属含量也会不断增高。污泥施肥可导致土壤中Cd、Hg、Cr、Cu、Zn、Ni、Pb含量的增加，且污泥施用越多，污染就越严重。污水灌溉污水灌溉一般是指用经过一定处理的城市生活污水来灌溉农用土地、森林及草地。由于大量工业废水同生活污水一起进入市政污水网，使得城市污水中含有的大量重金属离子，随着污水灌溉而进入土壤。近年来，污水灌溉已经成分我国农业灌溉的重要组成部分。北方旱作地区污灌最为普遍，占全国污灌面积的90%以上，所以土壤重金属污染也比较严重。含有重金属的废物的堆积废弃物堆中重金属含量一般比较高，污染的范围一般以废弃物堆为中心向四周扩散。重金属在土壤中的含量和形态分布特征受废弃物种类和释放率的影响，如铬渣堆放区的Cd、Hg、Pb为重污染，Zn为中度污染，Cr、Cu为轻污染。金属矿山的酸性废水污染金属矿山的开采、冶炼、重金属尾矿、冶炼废渣和矿渣堆放等，以及可以酸溶出含重金属的矿山酸性废水，随着矿山排水和降雨使之带入水环境或直接进入土壤，都可以直接或间接地造成土壤重金属污染。矿山酸性废水重金属污染范围一般在矿山的周围或河流的下游，在河流不同河段的重金属污染程度往往受污染源(矿山)的控制。河流同一污染源的下段自上游到下游，由于金属元素迁移能力的减弱和水体自净化能力的适度恢复，重金属化学污染的程度逐渐降低。（2）模型的建立此模型的主要目的是根据测得的污染数据,得到一个关于重金属污染物传播的偏微分方程,进而确定污染源位置的计算方法。我们根据给出的各个观测点的位置坐标和污染数据利用对重金属污染的线性插值，画出了重金属污染图，将附件三给出的的八种重金属元素的值作为重金属污染物的背景浓度（），观测点的处理如下：第一，光滑插值后如果有不光滑数据出现，考虑到计算机处理数据的误差，予以忽略。第二，分别对八种元素处理，得到各元素的污染坐标，对于相对误差不超过5%的坐标予以合并，取其平均值。第三，单个元素分析时，只考虑权重达到0.8以上的污染数据，其他数据因其影响太小，予以忽略。考虑到污染物在土壤中的运动状态，我们忽略了重金属污染的垂直扩散，假设重金属均匀扩散且是一维的，以污染源为原点，观测坐标X轴方向作为X轴，设污染物在原点的污染为连续过程，选用一维场中污染物质的二维扩散微分方程为：初始条件为：边界条件为：当处瞬时点源的解析解：其中是hantush函数，且常量参数：以污染源位置、 、为变量，结合格点搜索算法，反复调整变量的值，直到收敛函数的值最小，其中为个测量点的污染值；为背景值。（3）模型的求解此模型预测结果表明：As的污染源在第84个和第178个测量点附近；Cd的污染源在第22个和第95个测量点附近；Cr的污染源在第22个测量点附近；Cu的污染源在第8个和第22个测量点附近；Hg的污染源在第9个和第182个和第257个测量点附近；Ni的污染源在第22个测量点附近；Pb的污染源在第6个和第16个测量点附近；Zn的污染源在第61个测量点附近；分析结果：可看出：主要污染来自于工业区和交通区，且在交通区测得的第22个测量点污染最严重，说明22号附近的区域是主要污染源，与前文分析结果一致，说明模型正确。（4）模型验证我们使用matlab软件的griddata命令下的cubic（以三角形为基础的三次方程内插）插值的方式对海量数据进行了处理和绘图，得到了各元素的分布图，经过matlab的计算，可得到污染最严重处的坐标，从图中的位置可判断污染源的位置如下： 图5.4.1 As分布图可得和两个点，对比前面的模型，误差为 图5.4.2 Cd分布图可得和两个点，对比前面的模型，误差为 图5.4.3 Cr分布图可得一个点，对比前面的模型，误差为 图5.4.4 Cu分布图可得和两个点，对比前面的模型，误差为 图5.4.5 Hg分布图可得 、和三个点，对比前面的模型，误差为 图5.4.6 Ni分布图可得一个点，对比前面的模型，误差为 图5.4.7 Pb分布图可得和两个点，对比前面的模型，误差为 图5.4.8 Zn分布图可得一个点，对比前面的模型，误差为由此可看出，模型的准确性稍差。5.5问题四模型的优点、不足及改进随着全球城市化的加快，城市经济的快速发展，城市人口的不断增加，人类活动对城市环境质量的负面影响也日益凸显。城市土壤的污染就是人类健康的第一大杀手。土壤重金属污染的覆盖范围广、持续时间长、污染隐蔽性强还无法被微生物降解并且还可能通过食物链不断地在生物体内富集，甚至可转化为毒性更大的甲基化合物，对生物链中某些生物产生危害，最终页会在人体内蓄积而危害人体健康因此对城市土壤地质环境进行及时查证，并通过查证获得的海量数据资料开展城市环境质量评价，建立重金属污染的传播模型，研究人类活动下城市地质环境的演变模式，显得尤为重要。（1） 模型的优点本文在解决问题过程中，从各种影响因素和对不同评价方法可能产生的结果进行了有效的验证，结果表明，与单因子评价相比较，内梅罗增和指数法更加全面，但考虑到内梅罗指数法也有其缺陷，我们综合考虑后，觉得这两种方法都适合于研究重金属污染评价，我们决定了同时用这两种方法。让这两种方法的相结合，进行综和评价，使的结果更加完美，这是本文的一大特点。本文第二大特点就是，在问题一中我们运用matlab软件编程，分别建立了8种重金属元素在该城区的空间分布三维图，更加形象地表示出了元素空间分布，有效研究了金属的污染程度问题。也同时为解决后面问题做好基础。（2） 模型的缺点与改进本文研究的是整个城区的重金属污染状况，由于采样点数据个数有限，所得重金属含量空间分布一定是不连续的，只能采用以点带面的方式来表达土壤重金属的空间分布，所以对这些重金属含量进行空间展示结果会有一定的误差，该问题还有待改进。可以结合变程布设采样间距，这样以点带面更加合理。本文虽然用了单因子评价法和内梅罗综合指数法相结合，但是内梅罗污染指数法存在一定的缺陷，如有时过分突出极大值的影响，未考虑次大值对评价结果的影响。为了使得综合评价结果更加准确，可以采用一些数学方法如模糊综合评判或灰色理论进行评价研究。（3）模型的结果分析总体来说，模型基本上可以表示出污染源的传播特征，并且可以找出污染源，从而可以各污染物的总体分布图，数据吻合度也不错，不过由于时间仓促，无法更深入探讨模型的特征，是本文一点小小的遗憾。综合我们所建立末模型的优缺点，为更好研究地质环境的演变模式，我们还应该收集一些关于该城区地貌地质特征的知识，了解一下当地的土壤环境，各地区所含微生物的量，研究是否影响重金属含量，当地岩石侵蚀状况，了解一下当地环境温度，各个小区的人口变化以及人类活动，各地矿物质的分布状况等等。六 参考文献1宋新山，matlab在环境科学中的应用，/view/5c0e48cdalc7aa00b52acb76.html,2011-09-102蔡建安、张文艺，环境质量评价与系统分析，合肥工业大学出版社3刘来福、杨淳、黄海洋，数学建模方法与分析（原书第三版），机械工业出版社（2009-05出版）4黄建国，福建主要滨海湿地生态系统重金属污染特征及评价c2007年4月5谷超豪等编，数学物理方程（第三版）m，北京高等教育出版社，2007，P56-606张从，环境评价教程，北京：中国环境科学出版社，2004。7朱建平、范霄文Excel在统计工作中的应用，北京：清华大学出版社，2007.8李凌云、阎子健，兰州市西固区土壤重金属污染及空间分布特征，2011.49曾传祚，环境科学丛书第十二卷第四期七 附录附录一：As的空间分布Matlab程序a=xlsread(d:data,sheet1,b2:m320);T=a(:,10);D=a(:,11);K=a(:,12);c=a(:,1);X,Y,Z=griddata(T,D,K,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值X,Y,C=griddata(T,D,c,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值figure,surf(X,Y,Z,C)附录二：Cd的空间分布Matlab程序a=xlsread(d:data,sheet1,b2:m320);T=a(:,10);D=a(:,11);K=a(:,12);c=a(:,2);X,Y,Z=griddata(T,D,K,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值X,Y,C=griddata(T,D,c,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值figure,surf(X,Y,Z,C)附录三：Cr的空间分布Matlab程序a=xlsread(d:data,sheet1,b2:m320);T=a(:,10);D=a(:,11);K=a(:,12);c=a(:,3;X,Y,Z=griddata(T,D,K,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值X,Y,C=griddata(T,D,c,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值figure,surf(X,Y,Z,C)附录四：Cu的空间分布Matlab程序a=xlsread(d:data,sheet1,b2:m320);T=a(:,10);D=a(:,11);K=a(:,12);c=a(:,4);X,Y,Z=griddata(T,D,K,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值X,Y,C=griddata(T,D,c,linspace(min(T),max(T),linspace(min(D),max(D),v4);%插值figure,surf(