《GIS技术在城市表层土壤重金属污染分析中的应用-毕业论文》由会员分享，可在线阅读，更多相关《GIS技术在城市表层土壤重金属污染分析中的应用-毕业论文（11页珍藏版）》请在金锄头文库上搜索。
1、1 南阳师范学院20XX届毕业生毕业论文（设计） 题 目：GIS技术在城市表层土壤重金属污染 分析中的应用 完 成 人： 班 级： 学 制： 专 业： 地理信息系统 指导教师： 完成日期： 目 录摘要(1)1 污染特征分析(1) 1.1 创建分析区(2)1.2 数据转换(4)1.3 空间插值(5)1.4 重金属污染的模型分析(4)2 结论与讨论(4)参考文献(14)Abstract(15)第1页(共9页)GIS技术在城市表层土壤重金属污染分析中的应用 摘要：通过使用地理信息系统主流软件ArcGIS，在处理和分析已有数据的基础上，对重金属污染特征及与功能区之间的关联等展开研究。本文重点探讨了在城
2、市表层土壤重金属污染分析中所运用的关键GIS技术，及GIS技术在解决该类问题中的优势和前景。关键词：重金属污染；空间分析；空间插值；GIS应用 随着工业化、城市化的发展，土壤环境污染问题越来越突出。土壤中的重金属因不被微生物降解，不易移动，故会不断积累，造成严重污染，并可通过植物吸收进入食物链，造成农产品质量安全隐患，危害人类健康，本文通过对城市表层土壤重金属污染进行分析，总结了此研究中使用的关键GIS技术，为相关课题研究作参考。本文主要使用ArcGIS处理相关数据、分析重金属污染与功能区中间的关系，即，各功能区内各种金属污染的具体特征及与环境的相互作用关系，并重点阐述所应用的关键GIS技术，
3、GIS技术在解决该类问题中的优势和前景。 由于地理数据的纷繁复杂，利用GIS处理空间数据，在此基础上利用GIS强大的分析功能对空间数据研究挖掘出内在的信息、规律已经成为提取地理信息的有效手段。GIS技术在资源与环境应用领域中，它发挥着技术先导的作用。GIS技术不仅可以有效地管理具有空间属性的各种资源环境信息，对资源环境管理和实践模式进行快速和重复的分析测试，便于制定决策、进行科学和政策的标准评价，而且可以有效地对多时期的资源环境状况及生产活动变化进行动态监测和分析比较，也可将数据收集、空间分析和决策过程综合为一个共同的信息流，明显地提高工作效率和经济效益1。1 污染特征分析1.1 创建分析区对
4、某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（010 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据2。只给了采样点，没限定分析区边界，要确定分析区只能根据现有条件，最好的方式是做现有采样点的凸包，根据1000m随机采样的前提，然后做凸包500m的缓冲区作为研究区域范围。图1 为样本点建立分析区本文首先的前提假设是城市重金属污染和城市的功能区有莫大的联系。基于此，要根据已得到的采样点所属功能区划分城市功能区。由此可以采样点图层创建泰森多边形。将每个采
5、样点的重金属浓度分配给与本采样点的距离最小的区域。根据地学第一定律：地理位置上越是相近的区域，其所表现出的特性就越相近3。通过泰森多边形的创建，每个采样点所影响的大致范围就初具模型。图2 建立泰森多边形 泰森多边形具有广泛的用途。如，可以用离散点的性质来描述泰森多边形区域的性质；可用离散点的数据来计算泰森多边形区域的数据；判断一个离散点与其它哪些离散点相邻时，可根据泰森多边形直接得出，且若泰森多边形是n边形，则就与n个离散点相邻；当某一数据点落入某一泰森多边形中时，它与相应的离散点最邻近，无需计算距离4。1.2 数据转换 在分析区的范围中建立泰森多边形后，就需要根据每个多边形的功能区属性分成功
6、能区块，方便观察各功能区的分布，更为了以后的分析做基础。经过字段融合技术处理后得到5种功能区。接下来就要将功能区的矢量数据转化成栅格数据。 栅格和矢量数据结构都有各自的优缺点。在本课题中，就需要用到栅格和矢量数据结构相互转换算法。把矢量数据转化成栅格数据，利用栅格数据的单元格的值进行空间数据统计分析。观察各种重金属浓度在各个地域的变化，研究其分布特征。 矢量数据转换成栅格数据算法也是GIS技术的重要成分。由于矢量数据的基本要素是点、线、面，因而矢量数据向栅格数据的转换实质上是解决点、线和面数据的转换5。1.3 空间插值 利用空间插值可以得到每一个单元格的重金属浓度值。经过ArcGIS强大的渲染
7、功能处理后，可以很容易的观察出重金属污染的分布特征。图3重金属As污染图空间位置上越靠近的点，越可能具有相似的特征值，而距离越远的点，其特征值相似的可能性越小，这是空间插值最基本的理论假设。一般地，以下情况会用到空间插值：(1)现有的离散曲面的分辨率，像元大小或方向与所要求的不符，需要重新插值。(2)现有的连续曲面的数据模型与所需的数据模型不符，需要重新插值。(3)现有的数据不能完全覆盖所要求的区域范围，需要插值空间插值方法的主要目标：(1)对不足或缺失数据的估计，由于观测台站分布的密度及分布位置的原因，不可能任何空间地点的数据都能实测得到，需要用到插值，以了解区域内观测变量的完整空间分布；(
8、2)数据的网格化：规则格网能够更好地反映连续分布的空间现象，并对他们的变化作出模拟，对已知观测台站的观测数据进行空间内插，便可得到格网化数据；(3)内插等值线：以等值线的形式直观地显示数据的空间分布；(4)对不同分区未知数据的推求。 进行空间插值时，一般包括以下过程：(1)空间插值数据源的获取；(2)对数据进行分析，找出源数据的分布特性、统计特性，便于选择最恰当的插值方法；(3)插值方法的选择并进行插值计算；(4)对插值结果的评价；(5)运用多种插值方法进行计算，对各种方法的插值结果进行比较、分析并选择最佳的插值方法6。 本研究主要使用克立格插值。插值是一项复杂的运算，具体要忽略掉相关性很小的
9、采样点的影响。此时首先就要确定采样点之间具有相关性的距离值。本内容是重金属污染，样本是从土壤中取得，由于可能或受到河流、地势等因素的影响，还要排除掉在方向上的走势相关性。这涉及到探索性空间数据分析7。在本课题中主要用到了直方图分析，Q-Q概率图，趋势图，半变异/协方差函数云图。直方图是一种适用于对大量样点数据进行整理加工，找出其统计规律，即分析数据分布状态，以便对其总体分布特征进行推断的方法。本课题检验数据是否符合正态分布以及发现离群值，对于离群值特殊关注，排除错误样本。Q-Q概率图是根据变量分布的分位数对所指定的理论分布分位数绘制图形，是一种用来检验样点分布的统计图，如果被检验的样点数据符合
10、所指定的分布，则代表样点的点簇在一条直线上。本课题用它与标准正态分布（直线）对比，检验数据是否符合正态分布以及发现离群值。趋势分析是将样点的位置在X,Y平面上表示，对于感兴趣的属性值，通过垂直方向上的Z轴来表示，构成三维图。如果经过投影后曲线是平直，表明没有趋势，如果多项式有确定形式，就表明数据中存在全局趋势。半变异/协方差函数云图是数据集中所有样本点对的理论半变异值和协方差，把他们用两点间的距离函数来表示绘制图形。半变异/协方差函数云图技术是用来探索和量化空间自相关性。确定相关性所影响的范围，找出整体样本点的影响的临界范围。 经过优化后生成插值。具体用的插值方法要视所给样本点的重金属污染溶度
11、而定。此课题中主要应用到了普通克吕格插值。当区域化变量Z(x)的数学期望EZ(x)=m为未知数的时候，常用到普通克吕格法进行局部统计。普通克吕格模型为：Z(s)=+(s).其中，是未知的常量；(s)为随机误差。普通克吕格作为一种估计方法，有独特的灵活性。在运用普通克吕格法进行局部估计时，设待估块段为V，中心为x,其平均值为ZV。则 （1）在待估块段V的领域内，存在一组N个已知样点，其观测值为Z(x)，其数学期望也为M，令为ZV的线性估计量，由几个已知点观测值构成线性组合，即公式为， （2）在满足无偏性条件，最优性条件两个条件时，得到线性无偏，最优估计量。（1）无偏性条件 （2）最优性条件 普通
12、克吕格可用于估计那些看起来有某种趋势的数据。 1.4 重金属污染的模型分析在研究地区重金属分析时，可能遇到统计数学的方法来对数据进行定性与定量的分析。具体的评价方法有：土壤单项污染指数，该评价指标属于GB156181995中的二级标准，评价土壤单项污染情况。变异系数，将标准差作为算数平均数的百分率来表示，以说明样本的分散程度。本文利用内梅罗法计算重金属的综合污染指数。土壤单项污染指数计算公式为：土壤单项污染指数=土壤污染物实测值/污染物质量标准。变异系数计算公式为： （3） 内梅罗法的计算公式为： （4）这些算法需要通过栅格图层统计的运算，分析结果见表1。表1 各功能区8种重金属污染样品所占比
13、例元素各功能区污染样品所占百分比（%）1生活区2工业区3山区4交通区5公园绿地区As (g/g)0.618180.611110.166670.50.74286Cd (ng/g)0.727270.861110.257580.768160.54286Cr (g/g)0.454550.416670.227270.362320.25714Cu (g/g)0.772730.944440.257580.826090.68751Hg (ng/g)0.568180.694440.227270.521740.51286Ni (g/g)0.272730.416670.169700.195620.14826Pb (
14、g/g)0.590910.805560.196970.666670.4Zn (g/g)0.613640.777780.136360.717390.48571 注：污染样品即污染指数1的样品图4 综合污染指数从上表中可以看出各个功能区的综合污染情况，从大到小依次为：工业区生活区交通区公园绿地区山区 工业区几乎为5，属于重度污染。生活区大约比工业区小2，比山区大1.5，属于中等水平。交通区综合污染指数比生活区小了0.16，比公园绿地区大了0.06可以看出交通区的污染情况是次于生活区的。公园绿地区比交通区降了0.06，可见比交通区的污染情形要好。山区的综合污染指数值为1.5，重金属污染影响最小。从该图可以最直接的看出各种功能区所受的重金属污染是有很大差别的,所以功能区与重金属污染有莫大关系8。2 结论与讨论GIS技术在空间
