基于gis的铜矿区土壤重金属空间分布研究.docxVIP

基于gis的铜矿区土壤重金属空间分布研究 采矿业对土地的严重破坏，已成为中国乃至世界上一个非常严重和突出的问题。在我国南方地区许多大型金属矿的开发利用过程中,部分矿区土、水、气、生环境质量明显下降,严重影响食品安全和人畜健康。土壤重金属污染因其长期性和难以修复的特点已受到广泛关注,而掌握重金属污染的空间变化规律则具有特别重要的意义。 我国在小尺度区域的土壤重金属空间变异研究中,将GIS和地统计模型集成技术应用于矿区土壤重金属污染方面的较少。Arrouays D等用地统计学方法研究了短距离土壤重金属,包括铅、铜和铬的空间变异,结果表明针对研究区域具体情况设置采样方案十分重要。本文采用地统计学方法,探讨了浙江省诸暨市浬浦铜矿区4号坑道周围农田、菜地、荒地土壤重金属含量的空间变异性,为改善该区环境质量以及今后更多相关区域空间变异研究提供参考。 1 研究领域和计划 1.1 矿区地貌和地貌 浙江浬浦铜矿位于浙江省中北部的诸暨市,29°43′23″N,119°59′09″E,采矿区总面积0.8 km2。矿区内地貌为单一丘陵地貌。矿区地属中亚热带季风气候,土壤类型主要为黄色砂页岩发育的红壤和山地黄红壤。采样区域位于矿区4号坑道以北,主要以4级阶梯形式东西向分布。 1.2 网格法调查布点 针对该铜矿区地形比较单一、研究区域属于较小范围的情况,为确保样品有代表性,在4号坑道北侧,主要依据土地利用类型和4级阶梯各层走势及长度状况,在各级阶梯上分别按照网格法调查布点,总共布点48个,点间距为10 m。矿区中发生过漫灌的4号坑道位于采样区域的南部偏东方向,距采样区南侧约300 m;且有灌溉沟渠沿大致由南向北方向流经采样区域中部。 1.3 土壤样品的采集 根据布设的采样点,按照《农田土壤环境质量监测技术规范(NY/T 395—2000)》要求采集、制备土壤样品。每个采样点在30 cm的直径范围内,土壤深度为0~20 cm的耕层混合土壤样品,风干后筛分备用。 1.4 自然干法、粉碎、整理、提取 样品分析:待测样品按《农田土壤环境质量监测技术规范(NY/T 395—2000)》标准程序进行自然风干、粉碎、过筛、消解(或提取)、测定。 测定项目包括铅、砷、铬、铜、锌、镍全量、p H值,测定参见《土壤理化分析》。 1.5 分析与试剂 铅、砷、铬、铜、锌、镍的全量测定:英国帕纳科公司Axios生产的波长色散型X射线荧光光谱仪(型号AFS-820)。 1.6 空间依赖关系的半方差函数 本文应用地统计学的空间分析方法,该方法一般以半方差图作为描述空间变异性的工具。它作为一种最优的空间插值方法,已被广泛地应用到具有区域化变量特征的土壤学、环境科学和生态学等领域。地统计学分析利用Arc GIS9.1进行,包括半方差函数(semi-variogram)的计算和拟合比较、克立格法(Kriging)空间插值和模拟误差分析。半方差函数也常称为变异函数,是地统计学中研究土壤变异性的关键函数。设在一维(二维或三维)空间中不同位置x1,x2,…xn上某土壤特性的观测值为Z(x1),Z(x2),…Z(xn),半方差[γ(h)]可反映区域化变量的空间依赖关系,计算可用下式估算: 式中:N(h)是以此为间距的所有观测点的成对数目。以h为横坐标,γ(h)为纵坐标作图,即为半方差图。空间局部插值(Kriging法)是建立在变异函数理论及其结构分析基础之上,在有限区域内对区域化变量的取值进行线性无偏最优估计的一种方法,即(λi是权重系数,Z(x)是样品值,Z是估计值): 2 结果与分析 2.1 土壤环境污染 由表1看出,研究区6种重金属含量变异最高的是As、Cu、Zn、Pb。说明外来污染影响较大,多年来矿区开采造成的矿山酸性废水污染(据实地调查,研究区即4号坑道在1967发生过矿山酸性废水外涌),同时当地村民大量施用化肥、农药也可能是重要原因,主要使用的化肥为过磷酸钙、氯化钾、尿素。Cr变异程度最低,Ni的变异程度也较低,说明外来Cr与Ni污染影响可能较少。 根据土壤环境质量标准(GB15618—1995),研究区土壤(p H6.5)采样点铜和锌100%超标,铅91.67%超标,砷79.17%超标,镍68.75%超标,铬未超标。研究区的土壤铜、锌、铅和砷含量达到严重污染水平,镍污染相对较轻,采样点镍的平均值为45.757 8 mg·kg-1(镍的二级标准为40 mg·kg-1),而铬尚未达到污染水平,因此研究区主要为Cu、Zn、Pb和As的复合污染。此外对偏度和峰度值的分析可以得知砷、铬、铜、锌、镍、铅为负偏,且除Cr外均为尖峰态。将Cr、As、Cu、Zn、Ni、Pb等土壤重金属含量的分析数据,做出不同土壤重金属含量的频数分布直方图(Histogram),得出研究区土壤6种重金属含量均为对数正