基于DEM等空间因子协克里格土壤环境评价模型研究.docVIP

基于DEM等空间因子协克里格土壤环境评价模型研究 　　[摘要]本文以地统计学方法为基础，综合考虑地形起伏、河流、污染源分布和气候等因素对土壤污染的影响，对克里格空间插值方法进行改进，提出了基于DEM等空间因子的协克里格模型，使之更符合土壤环境实际情况，能够更加准确反应土壤污染程度的空间变异性。并以南方某市的土壤监测数据为例，得出基于DEM等空间因子的协克里格模型评价结果。 　　[关键词]DEM TIN 克里格（Kriging） 协克里格 GIS 土壤环境评价模型 Pb 　　[中图分类号] Q938.1+3 [文献码] B [文章编号] 1000-405X（2014）-5-262-2 　　1引言 　　土壤染污指人类活动产生的污染物进入土壤并积累到一定程度， 引起土壤环境质量恶化， 对生物、水体、空气或（和）人体健康产生危害的现象[1]，它传统的评价模型通常采用内梅罗指数法[2]，该方法简单明了，易于建立模型。然而，该方法没有考虑土壤环境的非平稳性、非匀质性、空间变异性等特点，容易掩盖小尺度土壤污染严重的可能。近年来，克里格法在土壤科学和环境科学领域得到广泛的应用，然而，由于克里格方法仅用距离作为土壤采样点对未知点推定的影响权重，没有考虑到地形及其它因素土壤污染的影响，因而当监测区域内地形及其它因素变化较大时，使用克里格法进行土壤质量评估可能存在较大误差。 　　针对上述问题，本文提出了基于DEM等空间因子的协克里格土壤环境评价模型，该模型考虑了地形起伏变化、河流、污染源分布和气候等因素，得出了基于DEM等空间因子的协克里格土壤环境评价模型。 　　2基于DEM等空间因子的协克里格（Co-Kriging _DEM）评价模型 　　土壤环境影响因子非常复杂，评价模型如果只用普通克里格方法，则评价结果会在山脉河流两侧等天然屏障处存在较大误差。协克里格方法是建立在地质统计学理论基础上的克里格方法，是以少量稀疏不规则分布的监测数据作为约束条件，反演研究区内未知点数据。本文所采用的南方某市的土壤数据中，土壤监测点数据较少，且满足克里格要求的半变异函数和协方差要求，适合使用协克里格方法。本文以高程、河流、重点污染源企业和气候因素等作为协同克里格的条件和分区标准，对监测区域分为若干个分区分别应用克里格空间插值方法，从而消除地形起伏、河流等对克里格法的影响，使之能够更详细地反映土壤环境的实质。三种分区标准分别如下： 　　2.1以DEM为分区标准 　　土壤环境是一个相对比较稳定的自然环境，不像大气环境和水环境那样多变。因此，自然的大尺度的地物地貌很可能就是土壤污染分区的自然分界线，其中DEM是一个重要的分区标准，本文对高程数据进行处理再以一定标准建立不规则三角网（TIN，以下简称TIN），进行分区。 　　2.1.1高程数据预处理 　　由于本文所讨论的土壤监测区域的高程数据为二维dwg数据，即只有X，Y的数据，Z的值全部为0，所有高程都是以标注的形式给出，在没有经过处理之前，这样的高程数据是没有意义的。笔者所采取的处理方法是在CAD中用VBA对所有高程点编程，把标注的高程附给相应的Z，从而确定高程点的（X，Y，Z）。 　　2.1.2TIN模型及其分区 　　为了得到适用表达南方某市土壤监测数据地形数据，要求得到和生成三维TIN模型。笔者运用ARCGIS的ArcScene中的3D分析模块，对已经形成三维离散点的高程点的SHP文件生成TIN。 　　根据土壤污染扩散的实际情况，由于自然景观、成土过程和气候的空间连续性，使得土壤污染在空间分布上并非是彼此独立和均匀的，而是在某一特定空间区域内彼此相关，其变异存在一定的规律性。在宏观上，地形变化和土壤环境扩散是一对绝对负相关关系，地形差别越大越剧烈，土壤环境之间的影响就越小。文中采用133m和334m两个高程做为分区界线，一区：高程133m，的士平坦，污染物容易扩散。二区：133m高程334m，低山丘陵，污染物不易扩散。三区：334m高程1200，山地丘陵，污染物较难扩散。 　　2.2以主要河流为分区标准 　　由于土壤环境的相对稳定性，土壤污染物难以在河两岸相互扩散，大河流的河面成为土壤污染扩散的天然屏障。文中利用南方某市的8m分辨率的正射遥感影像，提取一级支流东江，二级支流西枝江以及龙白河等把研究区域分割开来，西北部主要以山地丘陵为主，东南主要以沿海平原和低山丘陵为主。因势而造，在仅考虑地形地貌因素下本文把研究区分为4大二级区。 　　2.3以重点污染源企业和气候因素为分区标准 　　从大尺度来说，南方某市位于北回归线以南，濒临南海，地处亚热带，属南亚热带海洋性气候。降水量丰富，多年平均降雨量为1897mm，对土壤污染的扩散具有一定的影响力。该市经济发展速度相当快，“十五”期间，石