城市表层土壤重金属污染的空间分布特征分析
城市表层土壤重金属污染分析
城市表层土壤重金属污染分析摘要：文章利用单项污染指数衡量各区域内每种重金属元素对各监测点的污染程度，由尼梅罗算法得到8种重金属元素对各区域的污染程度；用因子分析法得到各种重金属元素污染的主要原因；由重金属元素的传播特征利用优化方法确定了污染源位置。
关键词:重金属污染尼梅罗算法因子分析法1 引言在以经济建设为一切工作重心的今天，工业化进程突飞猛进的同时重金属污染问题日趋严重。
重金属一旦进入土壤很难在生物循环过程中分解，当重金属在土壤中累积量超过土壤本身的承受能力时，不仅会影响土壤动植物的生长发育，而且还会通过植物的吸收、富集，并最终通过食物链进入人体，给人体健康带来巨大的危害。
目前，关于土壤重金属污染的研究已成为一个热点问题。
本文以2011年全国大学生数学建模竞赛题为背景，就某城区As、Cd、Cr、Cu、Hg、Ni、Pb、Zn八种主要重金属对土壤的污染状况展开研究。
考虑到不同的区域环境受人类活动影响的程度不同，所以按照功能，将城区划分为生活区、工业区、山区、主干道路区及公园绿地区。
研究过程中主要采用标点检测取样的办法获得各重金属的浓度数据，在此基础上给出了土壤重金属污染的研究办法。
具体做法是先由尼梅罗算法确定各区域的污染程度，同时利用因子分析法寻求污染原因，而后依据重金属的传播特征进行回溯，这样即可确定污染源的位置。
2 各区域重金属的污染程度对于重金属对土壤环境的污染程度，由于涉及多种元素，可用单项污染指数来衡量某一监测点某种元素对该点的污染程度，并用综合污染指数来衡量这八种重金属元素对该点的综合污染程度。
研究过程中，监测取样的方法获得的只是各金属在某一监测点的浓度，而通过这些数据很难直接评价污染程度，所以可选取一个统一的标准，将这些元素的浓度进行转化。
将各金属元素浓度背景值的上限作为标准，以浓度值在背景上限值中所占的比重作为污染程度。
可定义单项污染指数为：参照国家GB15618-1995《土壤环境质量标准》中对土壤质量等级给出的标准，就能得到重金属元素对各功能区的污染程度。
城市表层土壤重金属污染分析
一 、问题重述土壤是人类赖以生存的主要自然资源之一，也是人类生态环境的重要组成部分。
然而随着工业、城市污染的加剧和农用化学物质种类、数量的增加以及人类随着经济和社会及科学的发展逐渐向原始生态环境的扩进，土壤重金属污染日益严重。
目前，全世界各类重金属的排放量居高不下，其中Ni 的排放量大约100万吨、Mn 的排放量约在1500万吨、Pb 大约500万吨、Cu 约340万吨、Hg 大约在1.5万吨。
另据我国农业部进行的全国污灌区调查显示，土壤重金属污染具有污染物在土壤中移动性差、滞留时间长、不能被微生物降解的特点，并可经水、植物等介质最终影响人类的健康，总体上治理和恢复的难度较大。
随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
本文针对题目提出的几个问题，就以下四个方面展开讨论：（1） 应用点模式空间分析概念给出8种主要重金属元素在该城区的空间分布，这里不仅考虑每种重金属元素在该城区的空间分布，还考虑了不同区域中8中不同重金属元素的空间分布，从而结合不同的视角分析该城区内不同区域重金属的污染程度；（2） 重金属污染源主要来自随着大气沉降进入土壤的重金属、随污水进入土壤的重金属、随固体废弃物进入土壤的重金属和随农用物资进入土壤的重金属4个主要方面，本文结合主成分分析，给出该城区主要的污染源以及不同类型区域的污染源，进而结合实际讨论重金属污染的主要原因；（3） 针对现有数据的分布特征，包括该城区8种重金属空间分布和不同类型区域的重金属空间分布，建立数学规划模型，讨论了重金属扩散的中心位置和扩散方向，确定了污染源的位置；（4） 讨论了模型的优缺点，并分析了各类重金属污染对地质变化的前瞻性后果预测，具体给出了不同重金属对于环境污染的危害程度，提出了可能的解决方案，主要是针对预测结果的土壤重金属污染修复的可能性规划方案。
城市表层土壤重金属污染分析模型
城市表层土壤重金属污染分析模型摘要：针对经济的快速发展，城市人口的不断增加和人类活动对城市环境质量的影响也日益加剧的现状，该文对某城市城区表层土壤重金属进行了分析评价。
针在单因子指数评价基础上采用内梅罗综合污染指数评价土壤的综合污染，比较该城区的各个功能区重金属的污染程度。
基于重金属在大气、水体中传播特性的不同，利用高斯扩散推广模型确定重金属污染程度较大的污染源位置。
为更好地研究城市地质环境的演变模式，还应收集的信息有该城市常年的风速、冲洗系数、亨利系数；通过对以上数据的分析，建立重金属污染物在气体和土壤中扩散模型。
关键词：指数法因子分析重金属污染高斯扩散改进模型中图分类号：tu2 文献标识码：a 文章编号：1674-098x（2013）03（a）-0-021 问题分析针对海量数据，应从整体上对污染程度进行评价。
而内梅罗综合污染指数法评价土壤的综合污染，以突出最高一项污染指数的作用。
在土壤中有很多重金属元素有相似的存在形式和传播途径，并且有相同的污染源，因此在进行通过数据分析，说明重金属污染的主要原因时，基于统计原理建立起来正态模型，不同的重金属有不同的传播方式，其大体分为大气传播、水体传播、固体传播，因金属元素在土壤中大部分以稳定形态存在，故忽略重金属元素在固体土壤中的传播。
根据收集的信息和题目中的有关资料对重金属污染物的传播特征的分析，可将8种重金属污染物分为两类。
一类是在大气中传播，而大气传播的污染物最终经空气沉降进入土壤；一类是在土壤中传播。
对于在大气中传播的重金属污染物，文章建立重金属污染物在气体中扩散模型，根据所在的空间任意位置土壤表面的重金属污染物浓度的多少来确立污染源的位置，函数的最大值即为污染源的位置；同理建立了重金属污染物在土壤中的传播模型。
2 模型建立及求解2.1 土壤的环境质量评价与分级2.1.1 单因子指数法2.1.3 评价分级标准该文采用gb15618-1995《土壤环境质量标准》。
城市表层重金属污染分析
关 于 Ａ ：从空 间三维 图图 ２中可 以看到 ，１处的波峰很高说明该处污染情况很严重 ；有 ２４处明显波 ｓ －
峰 ，说 明相应位 置污 染情 况 比较严重 ：还有 １ 处面积较广且所处 高度稍低 ：表 明该 处所受污 染情 况相对严 重且污 染的范围较广 ； 同样分析二维等高线 图 １ ，图中有 ２处等 高线之 间的间距越来越 密集 且颜色 很深， 表 明该处受污 染情 况很严重，有 ２４处等高线 比较密集颜色相对较深表 明这两 处的污 染情况 相对严 重，还 ． 有 １ 等高线 间的距 离较密集但是所涉范 围较广说 明该处 的污染 也较严 重且污 染的面积很广 ．结合 有关背 处 景数据 可知 中心污染源 集中在 （ ２ １ ７ ９ ，（ ２ ９ ，３２ ） １８ ４ ０ ０ ）处 ．污染源 主要集 中在工 ５ ９ ，５ ３ ） １６ ６ ０ ４ ，（ ７ １ ，１７ ７ 业 区，还有 １ 处污染源位 于山林密集区 ，其污染级别不是特别严重 ． 重 金属 ｃ 、ｃ 、Ｃ 、Ｈｇ ｉ ｂ ｎ的污染程度 可作类似 分析 ． ｄ ｒ ｕ 、Ｎ 、Ｐ 、Ｚ
重金属 Ａ 分布平 面图 ｓ
据 ．表 中数值 ０表示在该污染级别下不存在观测样
嚣
本点 ．这是个大样本事件 ，可 以认为该级别污染很 轻微 ， 甚至不存在 这种级别 的污染， 而百分 比越大 ， 就说 明在该污染级别下涉及 的样本 点比较多，污染
波及范围较广 ．
傅
， ａ
５
；
中图分类号：ｘ ． １１ ３ ３ 文献标识码 ：Ａ 文章编 号：１ ０ — １ ５（ ０ ２ ３ ０ — ５ ０ ９ ８ ２ １ ）０ — １ １ ０ ３ ０
随着工业 发展和 城市化进程的加剧 ，通过交通运输 、工业排放和市政建设等造成城 市重 金属污 染越 来
城市表层土壤重金属污染分析
城市表层土壤重金属污染分析摘要随着近代工农业的发展，重金属通过各种途径进入土壤中，由于其不会被生物降解、迁移性小，很难从土壤中去除，容易在土壤中积累。
重金属污染物往往具有很强的生物毒性，当重金属积累超出土壤的承受能力或环境条件发生变化时，有毒物质可能会突然活化，导致严重的环境危害，因而土壤重金属污染有“化学定时炸弹”的说法。
另一方面，土壤重金属污染物可通过摄取、吸入、皮肤接触等多种途径危害人体康。
1 .问题重述现对某城市城区土壤地质环境进行调查。
为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。
应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。
另一方面，按照2公里的间距在那些远离人随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。
对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。
按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。
群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。
附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。
现要求你们通过数学建模来完成以下任务：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。
(2) 通过数据分析，说明重金属污染的主要原因。
(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。
(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？2.问题分析问题一首先要求得出8种主要重金属元素在该城区的空间分布，此问题主要是利用附录中所给的数据，研究给定采集点处的重金属含量。
城市表层土壤重金属污染分析
城市表层土壤重金属污染分析一、引言随着城市化进程的加快，城市土壤受到重金属等污染物的威胁问题日益凸显。
城市表层土壤是城市生态环境中的重要组成部分，受到重金属污染的影响会对人类健康和生态系统造成重大影响。
因此，对城市表层土壤中重金属污染的分析具有重要意义。
二、重金属在城市表层土壤中的来源城市表层土壤中重金属主要来源于工业排放、交通尾气、生活垃圾填埋和农药施用等活动。
这些活动导致了土壤中重金属含量的逐渐积累，从而引发了土壤污染问题。
三、常见的城市表层土壤重金属污染物种城市表层土壤中常见的重金属污染物种包括铅（Pb）、镉（Cd）、铬（Cr）、汞（Hg）等。
这些重金属对人体健康和环境造成严重危害，需要引起重视。
四、城市表层土壤重金属污染的影响1.对人体健康的影响–长期暴露于重金属污染土壤中会导致慢性中毒，严重影响身体健康。
–儿童和孕妇更容易受到重金属污染的影响，引起神经系统和生殖系统的损伤。
2.对生态系统的影响–土壤中的重金属会影响土壤微生物的活性，破坏土壤生态系统平衡。
–重金属还会进一步污染地下水，威胁周围生态环境的稳定性。
五、城市表层土壤重金属污染分析方法1.采样方法–选择合适的采样点位，并采用土壤钻孔或其它方法获取土壤样品。
2.实验分析–利用化学分析方法，对土壤样品中的重金属进行检测和分析，包括原子吸收光谱等技术手段。
3.数据处理–对实验数据进行统计分析和处理，得出城市表层土壤中重金属的含量及分布情况。
六、城市表层土壤重金属污染治理建议1.减少污染源–减少工业废气排放、加强交通管理，从源头减少重金属排放。
2.土壤修复–利用植物吸收、土壤修复技术等手段，对污染土壤进行修复和改良。
3.加强监测–定期对城市表层土壤进行监测，及时发现并处理重金属污染问题。
结论城市表层土壤中的重金属污染是一个严重的环境问题，对人类健康和生态系统造成威胁。
因此，开展城市表层土壤重金属污染的分析研究具有重要意义，可以为环境保护和城市可持续发展提供科学依据。
城市表层土壤重金属污染分析-2011高教社杯全国大学生数学建模竞赛全国一等奖A题
2011高教社杯全国大学生数学建模竞赛城市表层土壤重金属污染分析摘要本文主要研究重金属对城市表层土壤污染的问题,我们根据题目所给定的一些数据和信息分析并建立了扩散传播模型、权重分配模型、对比模型和转换模型解决问题。
首先，我们利用Matlab 软件拟出该城区地势图（图1），根据所给数据绘出该地区的三维地势及采样点在其上的综合空间分布图。
之后将8种重金属的浓度等高线投影到该地区三维地形图曲面上，接着分别计算8种重金属在五个区域的平均值，立体图和平面图（图1附件）相结合便可得出8种重金属元素在该城区的空间分布。
其次，在确定该城区内不同区域重金属的污染程度时，我们运用两种方法进行解答。
先假设各重金属毒性及其它性质相同，运用公式ijij P C P ='求出各区域各金属相对于背景平均值的比值作为金属污染程度，再运用1ji ij j C C ==∑求出各区域重金属污染程度，并将各区进行比较。
之后，我们加上各重金属的毒性，对各重金属求出权数，再结合国标重金属污染等级和已知的各组数据来确定金属的污染程度。
由上述两种方法的对比，更准确地得出重金属对各区的影响程度。
即： 工业区>交通区>生活区>公园绿地区>山区 并根据第一个模型的数据来说明重金属污染的主要原因。
再次，对重金属污染物的传播特征进行了分析，判断出重金属污染物主要是通过大气、土壤和水流进行传播。
在分析之中，我们得出这三种状态的传播并不是孤立存在的，而是可以相互影响和叠加的，因此，我们分别建立三个传播模型，再对这三个传播模型进行了时间和空间上的拟合，得出重金属浓度最高的区域图，并结合各重金属的分布图（图6）来确定各污染源的位置。
最后，本题中只给出了重金属对土壤的污染，对于研究城市地质环境的演变模式，还需要搜集一些信息（图7）。
根据每种因素对地质环境的影响程度进行由定性到定量的转化。
建立同一地质时期地质环境中各因素的正影响和负影响的权重分配模型，再对这些权重进行验算和修正。
城市表层土壤重金属污染分析模型
描述 8 种重金属元素的空间分布，应该将其用图像表示出来。为了更好地体现空间 性，首先画出该城区的地形图，参考 MATLAB 中等高线绘制函数 contour，将元素的浓度 表示为不同的颜色，绘制等浓度图。分别画出不同元素在整个城区的分布图和分别在各 个区的分布图。但这不能体现海拔与浓度的关系，因此下面还应画出海拔和浓度的关系 图。对于污染程度，首先考虑不同区域各个元素的污染程度，再综合考虑每个区域的综 合污染程度，引入内梅罗指数和不同元素的权重，计算不同区域的综合污染程度。
表 1 城区土壤重金属元素含量分析
元素 极小值
极大值
平均值
背景值
变异系数
As
1.61
Cd
0.04
Cr
15.32
Cu
2.29
Hg
0.0857
Ni
4.27
Pb
19.68
Zn
32.86
30.13 1.6198 920.84 2528.48
16 142.5 472.48 3760.82
5.6765 0.3023 53.5097 55.0167 0.2997 17.2618 61.7409 201.2026
3
图 2 8 种元素在该城区的分布情况图
图 3 8 种元素在一类区的分布情况图
4
图 4 各个元素浓度随海拔变化图 从该图可以看出浓度随海拔的降低而升高，这可能是人类在海拔低处的活动较多， 对重金属元素的浓度产生了较大的影响。 1.2 污染程度分析 下面对该城区中不同区域重金属的污染程度进行具体分析： i）整个城区元素含量分析 根据所给的数据，对该城区土壤重金属元素的含量进行分析，可得到表 1 的结果。
同样，可以得到其他四个区中 8 种元素的权重ω（k) , k = 2,3,⋯,5 .
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
城市表层土壤重金属污染的空间分布特征分析摘要：分析了某城市城区表层土壤中的 as、cd 等8种重金属在生活区、工业区、山区、主干道路区及公园绿地区的含量水平，得出了不同区域重金属的污染程度；运用污染负荷指数法对影响土壤各重金属主要因子进行分析，确定不同区域重金属污染的主要特征；建立重金属污染物的传播模型，运用kriging插值法对重金属含量进行最优无偏估计插值，对重金属污染的空间分布进行分析，揭示了城市表层土壤中重金属含量的空间分布特征。
关键词：城市城区；表层土壤；重金属污染；空间分布特征中图分类号：x530.2 文献标识码：a 文章编号：0439-8114（2013）06-1287-05随着城市经济的快速发展和城市人口的不断增加，大量工业“三废”、城市生活垃圾和污泥等污染物的排放与不恰当处置使重金属在土壤中不断积累，加重了土壤重金属的污染负荷，导致我国城市表层土壤的重金属污染日趋严重。
而城市土壤重金属污染是能有效反映城市环境污染状况的重要指标之一。
因此，对城市土壤环境异常的查证并应用查证数据开展城市环境质量评价、研究人类活动影响城市土壤环境的演变模式日益成为人们关注的焦点。
1 数据来源与研究方法以某城市城区为研究区，将其划分为间距1 km左右的网格子区域，按照每平方公里1个采样点对表层土壤（0～10 cm土层）进行取样、编号，并用gps记录采样点的位置。
应用专门仪器测试分析，获得了每个样本所含的多种重金属元素的浓度数据。
另外，按照2 km 的间距在那些远离人群及工业活动的自然区取样，将分析数据作为该城区表层土壤中重金属元素的背景值。
研究以2011年高教社杯全国大学生数学建模竞赛a题[1]所列的数据为数据来源，文献[1]列出了采样点的位置、海拔高度及其所属功能区的信息、8种主要重金属元素在采样点处的浓度和8种主要重金属元素的背景值。
按照功能划分，现代城市整个城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、3类区、4类区、5类区，不同的区域环境受人类活动影响的程度不同。
作者对城市表层土壤中的 as、cd、cr、cu、hg、ni、pb、zn 8种重金属元素在不同区域的含量水平进行了研究。
采用统计学原理和污染负荷指数法[2，3]，根据国家污染负荷指数与污染程度关系分级标准进行评价，得出该城市不同区域的土壤重金属污染程度。
对采样点进行最优估计，借助空间变异分析理论与方法[4]，确定不同区域重金属污染的主要污染源：工业区重金属污染主要来源于工业生产产生的工业垃圾；生活区的重金属污染主要是由于人们造成的生活垃圾形成的；在交通要道主要是汽车尾气的排放；山区人口分布比较分散，自然绿地的覆盖面积广，环境自净能力比较强，所以土壤重金属含量相对较低。
建立重金属污染物的传播模型，运用kriging插值法[5]对重金属含量进行最优无偏估计插值，对重金属污染的空间分布进行检验，从而揭示城市表层土壤中重金属含量的空间分布特征。
2 结果与分析2.1 重金属污染的空间分布根据统计学的基本原理，采用excel 2003进行数据分析，得出所有数据的平均值、最小值、最大值、标准差、变异系数、偏度、峰度等统计信息，同时采用区域识别法[6]对原始数据进行异常值检验。
对于已处理的数据，采用spss中的k-s方法[7]进行数据的正态分布检验，然后采用surfer 8.0软件绘出了8种重金属污染在该区的空间分布图，图中纵横坐标表示采样点的位置，图例表示采样点所在的功能区（图1a至图1h）。
2.2 不同区域重金属污染程度分析土壤重金属污染程度的评价分别采用研究区土壤重金属含量的自然背景值和国家一级标准值（保护区域自然生态、维持自然背景的土壤环境质量的限值）为参考标准，当评价标准为背景值时，相应的评价结果称之为元素富集[8，9]，当以国家一级标准为评价标准时，称之为元素污染。
2.2.1 研究区重金属污染程度的污染负荷指数法分析污染负荷指数法是在从事重金属污染水平的分级研究中提出来的一种评价方法，被广泛应用于土壤重金属污染的评价。
污染负荷指数法优点是该指数由评价区域所包含的多种重金属成分共同构成，并使用了求积的统计方法，通过这种方法对整个区域各个点位各种重金属进行了定量评价，并对各点的污染程度进行分级，能直观反映对环境污染最严重的元素和各种元素对环境污染的贡献程度以及重金属在时间、空间上的变化趋势[8]。
重金属污染负荷指数与污染程度关系见表1。
根据重金属污染负荷指数与污染程度关系，不同的污染负荷指数可将污染程度分为无污染、中等污染、强污染和极强污染4个等级，分别用0、ⅰ、ⅱ、ⅲ来表示[10，11]。
根据污染负荷指数法计算得到研究区不同区域的污染负荷指数，参照污染负荷指数与污染程度关系分级标准，得到研究区不同区域的重金属污染程度见表2。
2.2.2 研究区重金属单元素污染程度的单项污染指数法分析单项污染指数法是以土壤单项污染物的实测值与评价标准的比值，表示某一污染物影响下的环境污染指数，可以反映各污染物的污染程度。
这里污染分级标准采用国家一级标准为参考标准。
选用国家土壤环境质量标准中的一、二、三级标准分别作为土壤重金属污染评价的污染积累起始值、中度污染起始值和重度污染起始值，把土壤重金属污染分为4级（表3）。
利用单项污染指数计算公式[12]求出五类区中采样点的各重金属元素的单项污染指数，并求其平均值，将平均值作为此元素在该类区的单项污染指数，按单项污染分级标准对重金属单元素的污染程度进行分级。
根据上述方法，得到研究区表层土壤中各重金属元素的单项污染指数见表4。
从全区各元素的平均污染指数可以看出研究区表层土壤受as、cd、cr和 ni污染相对较轻，污染程度低；其中cu、hg、pb和zn的污染指数较大，表明研究区受这4种重金属元素的轻度污染。
2.3 重金属污染的传播模型地统计学是建立在区域化变量、平稳假设和内蕴假设等概念基础上，以区域化变量理论为基础，以半方差函数为主要工具，研究那些在空间分布上既有随机性又有结构性，或空间相关性和依赖性的自然现象的科学。
由于重金属污染的传播途径受很多因素的影响，比如土壤养分、物理特性、土壤重金属等性质的变异以及时间上的变异，基于地统计学的空间变异理论[13]，将地统计学方法与gis 相结合来研究变量的空间变异特征。
半方差函数又称半变异函数，能同时描述区域化变量的随机性和结构性[14，15]。
利用半方差分析时将这8种重金属元素的浓度进行对数转换，连同采样点的地理坐标输入地统计软件gs+，拟合半方差函数，并选择最佳拟合模型及其参数。
球形模型：上式中，c0为块金值；c0+c为基台值。
以as为例，在gs+软件中拟合半方差函数，得到4种拟合模型如图2，该4种拟合模型的参数见表5。
在半方差函数中，决定系数（r2）较大、残差（rss）较小回归效果越好，由表5可知，as元素的最佳拟合模型为线性模型。
对8种元素用相同的方法确定其最佳拟合模型，运用gs+7.0地统计软件系统完成半方差函数理论模型的自动拟合、最优选择、参数计算及半方差图绘制，结果如表6及图2（a～d）所示。
将基于地统计学的空间变异理论与gis的空间分析功能相结合分析研究区表层土壤重金属含量的空间变异特征和空间分布特征。
按照区域化变量空间相关性程度分级标准，当块金系数c0/（c0+c）0.75时，表明变量以随机性变异为主，变量的空间相关性很弱，变量受人为因素影响较大。
由表6可知，研究区表层土壤中cu和pb 的c0/（c0+c）较小，表现出强烈的空间相关性，其空间变异以结构性变异为主；as、cd、cr、zn和ni的c0/（c0+c）为0.25～0.75，具有中等程度空间相关，同时注意到，as、cd、cr和ni的c0/（c0+c）大于0.50，表明随机因素影响超过了结构性因素的影响，以随机因素占主导地位；hg的c0/（c0+c）值大于0.75，具有较弱的空间相关性，表明随机因素是其主要影响因素。
研究区表层土壤中cu和pb的变程较大，表明其空间相关性的范围较大；其余元素的变程相对较小，这可能说明灰尘中cu和pb的含量更容易受到随机因素干扰。
结果表明，研究区表层土壤重金属含量具有较为明显的空间变异性，均表现出很好的迁移性、空间结构性。
3 小结从块金值与基台值的比值可以看出，目前城市表层土壤中重金属含量受随机因素（主要是人为源输入）影响显著，其空间分布与工业生产等人类活动密切相关。
这主要是因为城市土壤接纳了来自工业生产与居民生活所产生的污（废）水、固体废弃物等各方面的90%的污染物，这些污染物中所携带的重金属也随之大量进入到城市土壤中，造成城市土壤重金属污染，与此同时，交通频繁，汽车尾气的排放，使得大气干湿沉降和城市灰尘作为城市环境中除土壤外重金属的主要赋存介质，吸附着大量重金属污染物，这些吸附在大气干湿沉降和城市灰尘中的重金属污染物通过大气沉降等途径持续地大量输入到地表环境。
参考文献：[1] 教育部高等教育司，中国工业与应用数学学会.2011高教社杯全国大学生数学建模竞赛赛题[eb ol].http：///html_cn/node/a1ffc4c5587c8a6f96eacefb8db cc34e.html，2011-09-09.[2] 范拴喜，甘卓婷，李美娟，等.土壤重金属污染评价方法进展[j].中国农学通报，2010，26（17）：310-315.[3] 张尧庭，方开泰.多元统计分析引论[m]. 北京：科学出版社，1982.[4] 刘静，蔡国学，刘洪斌.西南丘陵地区土壤有机质含量的空间插值法研究[j].西南大学学报（自然科学版），2008，30（3）：107-111.[5] 吴学文，晏路明.普通kriging法的参数设置及变异函数模型选择方法——以福建省一月均温空间内插为例[j].地球信息科学，2007，9（3）：104-108.[6] 张朝生，章申，等.长江水系河流沉积物重金属元素含量的计算方法研究[j].环境科学学报，1995，15（3）：258-264.[7] 刘付程，史学正，于东升，等.基于地统计学和gis的太湖典型地区土壤属性制图研究——以土壤全氮制图为例[j].土壤学报，2004，41（1）：63-70.[8] 国家环境保护局.中国土壤元素背景值[m].北京：中国环境科学出版社，1990.[9] 刘凤枝.农业环境监测实用手册[m].北京：中国标准出版社，2001.[10] gb15618-1995，土壤环境质量标准[s].[11] 尹君.基于gis 绿色食品基地土壤环境质量评价方法研究[j].农业环境保护，2001，20（6）：10-11.[12] 陈俊坚，张会华.广东省区域地质背景下土壤表层重金属元素空间分布特征及其影响因子分析[j].生态环境学报，2011， 20（4）：646-651.[13] 张长波，李志博，姚春霞，等.污染场地土壤重金属含量的空间变异特征及其污染源识别指示意义[j].土壤学报，2006， 38（5）：525-533.[14] 蒋炳言，方凤满，汪琳琳.芜湖市区土壤重金属污染评价及来源分析[j].城市环境与城市生态学报，2010，23（4）：36-40.[15] 郭旭东，傅伯杰，马克明，等.基于gis 和地统计学的土壤养分空间变异特征研究——以河北省遵化市为例[j].应用生态学报，2000，11（4）：557-563.。