城市表层土壤重金属污染分析.doc

2011高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： A0509 所属学校（请填写完整的全名）： 集美大学诚毅学院 参赛队员 (打印并签名) ：1. 唐翠清 2. 林超 3. 李睿 指导教师或指导教师组负责人 (打印并签名)： 日期： 2011 年 9 月 12 日33目录摘要1一、问题的问题重述2二、问题的假设和符号说明21、问题的假设22、符号说明2三、问题分析与模型建立、求解：31、问题（1）32、问题（2）83、问题（3）84、问题（4）14四、模型的改进15五、模型的推广15六、参考文献15七、附录15城市表层土壤重金属污染分析摘要：该问题研究的是城市表层土壤重金属污染分析问题。运用Matlab与Excel工具，结合“附件1列出了采样点的位置、海拔高度及其所属功能区等信息，建立对应的数学模型，解决问题。 问题1：通过Matlab和Excel二者结合应用，做出8种主要重金属污染物在城区的空间分布图，采用单因子污染指数法和土壤综合污染指数（内梅罗污染指数）建立基本数学模型，求解出城市表层土壤中8种主要重金属元素的环境质量指数，进行该城市不同区域重金属污染程度评价，得出该城市内不同区域重金属的土壤综合污染指数为1类区：1.85; 2类区:2.01; 3类区:1.03; 4类区:1.88; 5类区:1.50，因此污染程度从严重到轻度为2类区（工业区）4类区（交通区）1类区（生活区）5类区（公园绿地区）3类区（山区）。问题2：通过已知数据以及图表，观测出在该城市重金属污染程度高的城区，从而分析重金属污染的主要原因为2类区（工业区）的工业污染。问题3: 通过分析重金属污染物的传播特征建立高斯扩散模型，借助查找稳定度级别划分表和帕斯奎尔曲线的和的值的表格进行参数估计，求解出重金属的最大浓度，从而根据附件2所给的数据得出污染源的位置。问题4：通过分析问题3所建立的模型的优缺点，从而加以改进，使之更好地应用于研究城市地质环境的演变模式，加以推广。关键词：单因子污染指数法 土壤综合污染指数 高斯扩散模型 一、 问题重述： 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。根据题目所给的数据，通过数学建模来完成以下任务：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、 问题的假设和符号说明：1、问题的假设：1）土壤是由土壤中矿物质、有机质、土壤水分及土壤空气四个组成部分的合称。 这四个组成部分的分配和比例是不同的，按容积来说，土壤固体部分约有38的矿物质和12的有机质。而土壤含50的孔隙，在这部分孔隙中，土壤水分及土壤空气各占一半。在自然条件下，空气和水分的比例经常变动。现假设空气和水的比例是。2） 在重金属污染物传播的过程中，有生物化学反应，尤其是水中，会影响计算，对此在建模的过程中，我们忽略生物化学反应。3） 假设污水流入水中还未来得及沉淀就随水流运动，此时，即使有沉淀物，但重金属的含量可以忽略。4） 水中含有水体悬浮沉积物等它们是天然水体中存在的主要胶体物质。而水体中胶体物质对重金属有吸附作用，这不是重金属污染物主要来源，而且这会对我们求解模型造成影响，因此我们对此忽略不计。5） 假设表层土壤与空气是相连接，则表层土壤中的空气属于大气层。因此我们可以近似认为在求地面气体浓度，即是在求表层土壤含有重金属的气体浓度。6） 表层土壤的水分主要来自河流、灌溉，忽略天气变化，如降雨导致土壤的重金属含量的增加或晴天水分蒸发，使得所得的数据不精确。7） 附件1、2、3所给的数据均是精确的，且所测得的城市表层土壤重金属污染的浓度应保持不变。2、符号说明： 城市表层土壤中8种主要重金属元素，即时，依次表示为 城市表层土壤中8种主要重金属元素的环境质量指数 城市表层土壤中污染中8种主要重金属元素的样本实测值的平均值 城市表层土壤中8种主要重金属元素的评价标准的平均值 城市表层土壤重金属污染物某个位置（对应的坐标为：的浓度 该位置的空气重金属污染物的浓度 该位置的水中的重金属浓度 该位置的植物吸收的浓度 该位置的还有可能产生重金属的浓度 该位置的空气吸收占的比例系数 该位置的水吸收占的比例系数 该位置的植物吸收占的比例系数 该位置的还有可能产生重金属的比例系数 三、 问题分析与模型建立、求解：1、问题（1）：首先，我们画出8种主要重金属元素在采样点的浓度在城市采样点位置的分布图：图1-1 8种主要重金属元素在采样点的浓度在城市采样点位置的分布图 由图我们可以观察出，8种主要重金属元素在图中的左下角区域分布密集。由此，我们做出8种主要重金属元素在该城区的空间分布图，如下所示。图1-2 8种主要重金属元素在该城区的空间分布图 三维曲线表示地理三维，有颜色的三维曲面表示各金属的浓度，颜色越靠近红色说明浓度越高。2 分析该城区内不同区域重金属的污染程度：我们先做出城市表层土壤中8种主要重金属元素在该城市五个区域的饼图，如下图所示。图1-3 城市表层土壤中8种主要重金属元素在该城市五个区域的饼图 从上图以及通过前面所做的8种主要重金属元素在该城区的空间分布图，我们可以得知8种主要重金属元素主要分布在2类区（工业区），也就是说2类区（工业区）的污染程度最大。于是，我们采用污染指数法对该城市表层土壤的重金属污染程度进行现状分析。由于地域差异等各种因素的影响,污染指数法没有统一的评价标准。通常进行土壤环境质量评价时，参照中国土壤环境质量标准（GB15618-1995）来确定研究的评价标准。采用土壤综合污染指数（内梅罗污染指数）进行土壤环境质量综合评价，划定污染等级。土壤综合污染指数计算公式如下： 单因子污染指数法：土壤单项污染指数= 即 内梅罗污染指数：土壤综合污染指数=因此，我们根据附件1、附件2、附件3所给的数据以及上述公式，得出如下表格：表1-1 各个城区8种主要重金属元素的平均值城市功能区As 平均值(g/g)Cd平均值(g/g)Cr平均值 (g/g)Cu平均值 (g/g)Hg平均值(g/g)Ni平均值 (g/g)Pb平均值(g/g)Zn平均值 (g/g)1类区6.27 0.28996 69.02 49.40 0.93041 18.34 69.11 237.01 2类区7.25 0.39311 53.41 127.54 0.64236 19.81 93.04 277.93 3类区4.04 0.15232 38.96 17.32 0.40956 15.45 36.56 73.29 4类区5.71 0.36001 58.05 62.21 0.44682 17.62 63.53 242.85 5类区6.26 0.28054 43.64 30.19 0.11499 15.29 60.71 154.24 表1-2 各个城区8种主要重金属元素的土壤单项污染指数及综合指数城市功能区土壤单项污染指数土壤综合污染指数As(g/g)Cd(g/g)Cr(g/g)Cu(g/g)Hg(g/g)Ni(g/g)Pb(g/g)Zn(g/g)1类区1.74 2.23 2.23 3.74 2.66 1.49 2.23 3.43 1.85 2类区2.01 3.02 1.72 9.66 18.35 1.61 3.00 4.03 2.01 3类区1.12 1.17 1.26 1.31 1.17 1.26 1.18 1.06 1.03 4类区1.59 2.77 1.87 4.71 12.77 1.43 2.05 3.52 1.88 5类区1.74 2.16 1.41 2.29 3.29 1.24 1.96 2.24 1.50 或表示未污染，为轻污染状态；为中度污染状态；为重污染状态，且值越大，则污染越严重。由此我们可以知道均大于1，即各个城区8种主要重金属元素都造成了污染。而由土壤综合污染指数可得1类区、3类区、4类区、5类区均为轻污染状态，2类区为中度污染状态且该城区内不同区域重金属的污染程度为2类区（工业区）4类区（交通区）1类区（生活区）5类区（公园绿地区）3类区（山区）。2、问题（2）：重金属污染是由重金属或其化合物造成的环境污染。重金属污染与其他有机化合物的污染不同。不少有机化合物可以通过自然界本身物理的、化学的或生物的净化，使有害性降低或解除。而重金属具有富集性，很难在环境中降解。由第（1）问的数据分析结果，我们知道该城区表层土壤重金属的污染主要在工业区。因此，我们得出重金属污染的主要原因是工业污染，其次是交通污染和生活垃圾污染。工业污染大多通过废渣、废水、废气排入环境，在人和动物、植物中富集，从而对环境和人的健康造成很大的危害，工业污染的治理可以通过一些技术方法、管理措施来降低它的污染，最终达到国家的污染物排放标准；交通污染主要是汽车尾气的排放，国家制定了一系列的管理办法，例如：使用乙醇汽油、安装汽车尾气净化器等；生活污染主要是一些生活垃圾的污染，废旧电池、破碎的照明灯、没有用完的化妆品、上彩釉的碗碟等，对于重金属的污染只要我们从其来源加以控制，就多多少少可以减少重金属污染。目前我国塑料生产企业的工艺、设备、技术研发较落后，是造成污染严重的主要原因，而管理不善、地方保护及人们环保意识淡薄，加剧了污染，强化治理迫在眉睫。生产企业应放眼未来，倡导环保，使用环保型助剂才能使PVC行业健康长远发展。3、问题（3）：分析重金属污染物的传播特征：a) 重金属污染物通过空气进入土壤：空气中的重金属是存在于大气颗粒物(TSP)当中的，而粒径小于10m的大气颗粒物(可吸入颗粒物PM10)更含有高含量的重金属。据报道，大约75%90%的重金属分布在PM10中，且粒径越小，重金属含量越高。同时，风、湍流是决定污染物在大气中稀释扩散的最直接因素。b) 重金属污染物通过水进入土壤：水中的重金属主要可以分为两大部分传播：首先是含于水中，随着水的流动而流动。而主要影响因素有河流量、污水中的污染物浓度；其次是重金属进入水中一段时间会沉淀，成为沉淀物存在水中。c) 随含重金属废弃物堆积进入土壤：含重金属废弃物以其堆积为中心，一般以废弃堆为中心向四周扩散。d) 通过农用物质进入土壤：农用物质进入土壤主要是以农药、化肥和塑料薄膜使用为主。不合理地施用含有铅、汞、镉、砷等的农药和化肥，都可以导致土壤中重金属的污染。还有农用塑料薄膜生产应用的热稳定剂中含有Cd、Pb，在大量使用塑料大棚和地膜过程中都可以造成土壤重金属的污染。模型的建立与求解： 我们通过分析重金属污染物的传播特征，建立模型： 但是结合实际来研究，会发现：1)植物吸收重金属的含量远小于产生重金属的含量，故可以近似忽略；2）通过其他方法产生重金属对土壤产生影响，这些方法是在特殊的地方进行的，故我们可以的忽略其影响；3）计算重金属在水中的浓度，虽然是一种重要的影响，但是我们并不知道需要求的城区水源分布状况等具体情况，而空气是到处都有，没有以上的问题，我们可以主要考虑空气重金属污染物的浓度这一问题；4）对于土壤的影响，空气和水所占的比例重，根据相关资料，我们可以将、约为0，而空气和水的比例是。现在我们主要讨论空气重金属污染物的浓度：.对于该位置的空气重金属污染物的浓度的计算：1)符号补充说明： 平均风速 源强（即源释放速率） 风所吹向的目的地的距离 横向距离 纵向距离 泄漏源的有效高度 烟气抬升的高度 释放烟气的物体自身的高度 水平扩散参数，污染物在y方向分布的标准偏差 铅直扩散参数，污染物在z方向分布的标准偏差2）模型的求解：a、 假设污染物排到空气中有一定的高度，则以某个位置为原点的坐标系中的垂直坐标为，根据高斯扩散模式，可得该位置形成的污染物浓度为： b、 取z=0时，即可以得到地面气体浓度： (1)c、 首先，我们可以分析：污染源高除了烟囱的高度以外，在地面上也可以为烟气抬升的高度，即。根据大气稳定度可知我们的测量数据是在距离地面10m处，而实际上，在制作大气稳定度表时，这距离地面10m处是个大约值，根据有关资料，大气稳定度只要在10m高度以上，就能用的，还有就是烟气抬升高度是可以达到10m，所以我们可以取h=10m，即可求得表层土壤金属污染物的浓度： d、 令h=10m，即可得到对应于地面轴线气体浓度的表层土壤金属污染物浓度的计算公式： e、 连续点源地面最大浓度： 由于，则 由于、，则求得 ，3）参数的估计：扩散参数，：扩散参数是表征湍流扩散剧烈程度的物理量，是影响污染物浓度的重要参数。扩散参数，的大小与观测到的风速、云的密集程度、离地高度、气体泄漏持续时间和日照等等因素有关。大气的湍流结构和风速受大气稳定度的影响。大气稳定度是在10m高度以上，才能用的。为了计算大气稳定度，我们将大气扩散稀释能力分为6个等级：A 极不稳定，B 不稳定，C 弱不稳定，D 中性，E 弱稳定，F 稳定。若稳定级别为AB，则表示按A 、B级的数据内插。根据表3-1确定稳定度级别：表3-1 稳定度级别划分表地面风速 (距地面10m处)/白天太阳辐射阴天的白天或夜间有云的夜间强中弱薄云这天或低云5/10云量4/106CDDDDD我们先规定：A、白天太阳辐射符号：强为、中为、弱为；阴天的白天或夜间为；有云的夜间薄云遮天或低云5/10为和有云的夜间云量4/10为；B、地面平均风速。C、为选取的天气状况，根据天气情况查表即可得到大气稳定度，例如，我们取，则得大气稳定度为A。从表中查得对应的扩散参数，：要确定扩散参数，主要有两种方法：利用扩散曲线确定，和根据英国伦敦气象局给出的表3-2来确定。扩散曲线是求扩散参数，的基础，但是查找扩散曲线来确定，难懂且误差大，若我们需要查找20km距离内的扩散参数，那么用根据英国伦敦气象局给出的表3-2，更好看懂和理解。我们就选用英国伦敦气象局给出的表3-2来确定扩散参数，。表3-2 帕斯奎尔曲线的和值帕斯奎尔曲线的和值稳定度ABCDEF标准差距离r/km0.127.014.019.110.712.67.448.374.656.053.724.192.330.249.829.335.820.523.314.015.38.3711.66.057.914.190.371.647.451.630.233.520.521.9410.75.580.440.543.326.528.815.321.410.714.46.980.511210581.451.253.532.635.318.126.513.017.78.370.613215395.862.862.838.640.920.931.214.920.59.770.817027912384.680.950.753.527.040.018.626.512.11.020745615110999.161.465.632.148.821.432.614.01.224367417813311673.076.737.257.724.727893020315713383.787.941.965.627.043.317.21.6313123022818114995.398.647.073.529.348.819.11.825320716610710952.182.331.654.520.52.027823318211612156.785.633.560.521.93.039536326916717379.112941.986.527.04.050849333521922110016648.610231.26.072377747431631514023760.915637.78.060340940517730670.720742.81073549848821236679.124246.51256924442787.428550.216 72930754410036555.820 88437265911143760.5确定风所吹向的目的地的距离和扩散参数，估计：根据前两题我们可以得到该城区内不同区域重金属的污染程度，既而可以得出不同重金属元素主要分布哪个区域。As主要分布在2类区（工业区），Cd 主要分布在2类区（工业区），Cr主要分布在1类区（生活区），Cu主要分布在2类区（工业区），Hg主要分布在2类区（工业区），Ni主要分布在2类区（工业区），Pb主要分布在2类区（工业区），Zn主要分布在2类区（工业区）。从而我们可以整理知道重金属主要分布在2类区（工业区）有As、Cd、Cu、Hg、Ni、Pb、Zn，而重金属主要分布在1类区（生活区）只有Cr。为此我们结合所举的例子和的图表，由于我们是为了确定污染源的具体位置，则我们就选取下风向距离最近的对应的和，从而求出扩散参数的值。1） 平均风速：我们根据中的例子得知。2） 源强（即源释放速率）：假设重金属As、Cd、Cu、Hg、Ni、Pb、Zn、Cr对应的源强分别。源强主要是通过计算或实测得到的，但实际操作上我们遇到参考资料有限和无法进行实测，这导致于对于重金属的源强，我们难以求出准确的数据。对于这个问题的处理，我们采用与其有所近似的数据，即我们选取“现有污染源大气污染物排放限值”中的数据作为8种重金属的源强，理由如下：首先，我们通过前两题的分析，我们得出8种重金属中的7种重金属主要分布在工业区；其次，在工业区对于污染物排放有一部分是来自于排气筒，而我们选取“现有污染源大气污染物排放限值”中的数据就是针对排气筒，因而可以用“现有污染源大气污染物排放限值”中的数据计算8种重金属的源强；此外，根据大气稳定度可知我们的测量数据是在距离地面处，但是实际上，在制作大气稳定度表时，这距离地面处是个大约值。根据有关资料，大气稳定度只要在高度以上，就能用的，根据“现有污染源大气污染物排放限值”中的数据。因此我们可知排气筒高度最低为，但为了计算和查表的方便，我们在计算的过程中大胆的统一假设泄漏源的有效高度。对于现有污染源大气污染物排放限值”中的数据，题中所给的重金属主要是以化合物的形式存在的，这就导致排放的气体并不纯，但是我们进行的是参数估计，进而我们可以将其认为代表某种重金属的源强。现有污染源大气污染物排放限值(排烟筒高15m)（节选）序号污染元素最高允许排放浓度/mgm-1最高允许排放速率/kgh-1一级二级三级1Cr0.08禁排0.0090.0142Pb0.9禁排0.0050.0073Cd1禁排0.060.094Ni1禁排0.180.285Hg1禁排0.00180.0028其中，上表是我们从“现有污染源大气污染物排放限值”节选出有用数据，表中主要体现以上5种数据。根据上表，我们选取重金属的二级最高允许排放速率为例，进行接下来的计算。查表得。1） 浓度的估算：将确定的参数、代入（2）求得Cr的地面最大浓度为由于参数源强（即源释放速率）的数据有限，我们无法得知该区实际的源强，既而得出的Cr的地面最大浓度与题目所给的数据相比差很多。若计算出准确的Cr的地面最大浓度，那么接下来就是求解具体位置：我们先计算Cr在工业区的浓度差值，然后，在比较36组Cr在工业区的浓度差值，取浓度差值最小的那一组，即可作为污染源位置。若得出不只一组浓度差值最小，那么我们要进行分析：倘若这几组数据位置相近，那么，我们可以推知污染源是一片区域；若几组数据位置相差大，那么，我们可以推知污染源不只一个地方。根据以上方法，我们在有准确的除Cr以外的重金属的地面最大浓度的基础上，可以算出其它重金属的污染源位置。4、问题（4）:分析你所建立模型的优缺点:优点：我们是通过建立表层土壤重金属污染物某个位置的浓度与其空气重金属污染物的浓度、水中的重金属浓度、植物吸收的浓度、还有可能产生重金属的浓度的比例关系，这种模型将所有的情况考虑进去，计算过程少，主要是查找图表，而且还浅显易懂。缺点：通过此模型找出污染源的位置，但是只是一个大致的位置，并不是最精确的，与实际有误差。为了更好地研究城市地质环境的演变模式还应收集该地区的每年生活、工业的排污量；该城市处理污染物的能力；水土，地貌状况，以及环境保护力度；认为改造环境的能力；还有空气污染，自然灾害的发生频率等等。有了这些数据以后建立因子分析法，回归分析，曲线拟合，模糊数学等模型解决问题。 四、 模型的改进：问题三的模型在参数的估计中，对于源强（即源释放速率）的估计不够准确，现在，我们用另一种方法求解源强（即源释放速率）。我们对于某种重金属元素在各个地方所对应的浓度进行拟合求解出拟合的浓度；通过式子，其中为选取的体积、为单位时间。我们可以令。剩下求解浓度的过程与上面相同。对于问题三的模型，我还可以改成先用matlab程序求解每个重金属拟合曲线图，研究每个重金属拟合曲线图，结合excel对找出在最接近拟合曲线的浓度所对应的位置，即为污染源地。对于问题三的模型，我们可以添加其他参数计算除了空气重金属污染浓度对于表层土壤重金属浓度的关系比例，即根据实际情况计算，进而使表层土壤重金属浓度更加精确。五、 模型的推广：本文所述的高斯扩散模型是来解决城市表层土壤重金属污染分析过程中的实际问题，应用过程中充分考虑了现实环境的多样性，使得模型贴近实际，通用性和推广性强。这种方法还可以应用于其他类的多目标评价决策与分类问题，如雨量评价，多个目标的决策问题等。六、 参考文献：1马莉,MATLAB数学实验与建模, 清华大学出版社, 2010年1月2郝吉明,马广大,大气污染控制工程,高教出版社,1989年3杨若明、金军，环境监测，化学工业出版社，2009年1月4重金属污染的主要原因，/question/317414370.html，2011年9月10日5 空气中重金属处理，/question/199687769.html ，2011年9月11日七、 附录：1、重金属元素在采样点的浓度在城市采样点位置的分布图的程序 A=747817.84 13737315.93 132117914.90 017876.56 104921276.35 1647272814.08 288336178.94 238336929.62 270822957.41 293317678.72 42338955.93 404318959.17 242739715.72 352643574.49 506243395.51 4777489711.45 586849046.14 653456417.84 548160047.41 459246038.50 248659995.51 329960189.84 357362139.39 474164343.30 537586434.09 563579656.14 539486315.31 529173493.69 4742729321.87 4948729318.38 5567678210.53 700462263.50 730452306.35 704846005.51 818044964.49 932843113.50 909053655.51 804954394.29 807764014.29 801772106.56 6869728616.58 705683487.41 774782605.93 845789914.69 946083114.90 906276395.31 931967994.29 1063164725.51 1068555284.69 1064344727.20 1170244805.31 1173055324.90 1148263544.90 1070081843.89 1063087743.69 1167886183.11 1190277093.89 1324470563.89 1274684502.91 1285589453.30 1379796214.90 1432586664.09 1546786582.91 1244243292.72 1309343393.11 1392053543.30 1484455193.30 1656960556.14 1638766093.69 1606173524.49 1565875943.11 1429874188.06 1417766843.69 1509269363.69 1277857993.507250775369 19007114886.14 187381092110.9935 181341004630.13 1719898103.89 1714490812.91 1839391831.96 1976788106.98 2100688192.91 2109194825.93 2284691495.93 2366497907.41 22304105274.29 21418107214.69 21439113835.72 20554112286.77 20101107744.49 21072104043.11 2021599512.9129 19968129617.63 21766123485.93 22674121732.34 22535112932.91 2522157955.72 2645355772.34 2641665086.56 2781655814.69 2536164236.35 2406573535.10 2599870324.69 2717777713.50 2642486394.69 2607388074.49 2463194223.30 2470295222.91 2546198344.09 24813107995.72 26086110944.90 26015120784.90 27700116094.09 27696116215.93 27346133312.91 26591137152.72 27823147372.34 27232144822.53 24580133193.89 24153124502.34 22965135352.34 23198135231.96 24685142782.72 2865487553.11 24003152863.50 21684131012.72 22193121851.77 1707958943.50 1525551102.53 1500755359.17 351825717.41 346923085.72 376221708.06 392721105.72 415322998.72 32677936.77 468413646.56 549512056.77 566416536.98 554120936.14 545127576.69 402029908.00 402639136.69 510140808.23 543839949.35 538230128.90 531420603.77 550311275.41 56361337.78 66053745.62 709313815.41 710024494.58 683734906.91 790639785.00 804530525.62 839420356.91 840310757.78 807906.26 966312887.56 946922864.79 917832992.77 909539756.26 1022538217.34 1021027894.17 1034017645.00 1155715817.56 1141525855.62 1164935156.05 1273440154.79 12696302423.72 1240020606.47 1259110634.17 1376513535.00 1369423574.17 1385533456.26 1486225245.00 1489616031.61 153877294.58 1581023075.41 1603230613.57 1580139662.00 1508735123.37 1687227986.47 1773436293.37 1682342076.91 1700847752.38 1720362187.56 1700572128.67 1694774874.79 1630182995.41 1790482875.83 1830373855.62 1843865394.38 1855655882.77 1895448745.41 1801244145.62 1907285196.69 2028285907.56 2147585402.77 2145075551.61 2026175865.62 1956973482.97 1941169344.58 1950160915.20 2058265483.17 1990953005.00 2101857647.56 2217654925.20 2335953256.05 2323865022.77 2262448182.38 2170365915.20 500688468.67 573496595.41 6395104437.12 7405109814.58 8446112008.23 7612119384.58 7912128406.05 8866131438.23 9296131024.17 9475120003.97 92121130510.74 8629120863.77 77761061310.27 8622106385.20 923798726.47 830797266.47 710694678.90 642388313.37 745889206.69 890488684.17 1054795918.23 10398103604.38 10395112035.20 11529112435.20 11563102985.41 1164693815.00 1264195609.58 1400089703.17 1420799803.7713 12734103447.34 1272776915.6279203734 1524891063.57 1642890694.38 16289100725.41 16267110584.38 16440120685.62 16440132326.0526209158 12153123368.67 11958133136.47 10800132827.12 10022122043.77 9333146315.41 92771614810.97 11121164329.81 10856147278.23 12644149435.41 12625162592.77 9036175387.78 10599179806.47 12632179495.6277573891 16629144817.1238 20591135493.97 20983158628.00 20177176423.57 19041157693.371797126283 16607173653.9717 22605143012.77 23146153823.17 22046176342.97 23785176433.17 25981180511.80 27380182023.17 25021162903.77 23325167014.79 26852161143.57473423 103521713310.74 90951641411.68 10510153147.34 1395456156.05 1014216625.41 1776535616.26 692456966.47 467837656.47 618220054.79 598525677.56 765319529.35 ;x=A(:,1);y=A(:,2);z=A(:,3);scatter(x,y,5,z)%散点图figureX,Y,Z=griddata(x,y,z,linspace(0,28654),linspace(0,18449),v4);%插值pcolor(X,Y,Z);shading interp%伪彩色图figure,contourf(X,Y,Z) %等高线图figure,surf(X,Y,Z)%三维曲面2、重金属元素在该城区的空间分布图的程序 A=7478157.84 13737311