农田大尺度下土壤电导率的空间分布模拟方法.docxVIP

农田大尺度下土壤电导率的空间分布模拟方法 0 地质统计模拟 在区域空间特征的研究中，克立格法的空间变化分布相对平坦，因此很难重建区域变量的空间结构，并且很难反映极端值。而在很多情况下,再现区域化变量的空间结构,发现空间变异中极端值的分布类型和规律比局部值的精确估计或大范围的平均值的估计更重要。地质统计模拟就是根据已知的均值、方差和空间结构特征,用地质统计学的方法产生随机变量的空间分布,它克服了克立格法的缺点,较好的再现区域化变量的空间结构性。 地质统计模拟研究的开端是以转向带法的提出为标志的,转向带法是将多维(二、三维)的模拟转换为一系列对等的一维模拟。运用一维的协方差或变异函数沿着不同的直线来进行模拟,然后在二维或三维空间的每个点模拟值是由一维模拟过程中相对应值的加权平均而获得。然而,用转向带法进行条件模拟的速度慢,步骤烦琐, 且通常得到的模拟实现显示出呈带状的人为影响(这是由于最多只能用15 条线,对三维空间只是规则但很粗糙的划分)。另外,它也不能简单地处理与模拟网格的坐标轴方向不一致的各向异性的方向。 而在转向带法基础上发展起来的序贯高斯模拟方法只需要一个有效的克立格算法,计算类型单一,且能自动处理各向异性的问题。此方法主要根据现有数据计算待模拟点值的条件概率分布,从该分布中随机取一值作为模拟现实。每得出一个模拟值,就把它连同原始数据、此前得到的模拟数据一起作为条件数据,进入下一点的模拟。 土壤盐碱化一直是制约河套灌区农业生产发展的主要障碍因素,弄清土壤中盐分的空间变异及其分布规律对河套灌区盐土资源的合理开发、利用和改良有重要意义。由于上述的原因,本文选用了序贯高斯模拟方法对河套灌区部分农田的土壤电导率进行了条件模拟。 1 土壤样品采集与分析 试验区设在巴彦淖尔市五原县,位于河套灌区的中部,面积约8 000 hm2,土壤以壤土为主,该地以种植玉米、向日葵为主,2007年6月26日在试验区用剖面法采集土壤样本。方案如下:将试验区划分若干1.5 km×1.5 km的网格,南北、东西基线上加密布点,如图1所示,共设置了60个取样点。用剖面法采集表层0～20 cm深度内的土样, 用便携式快速测定仪器(WalklAB),共测得60个土壤电导率值。 在数据处理过程中采用DPS Version8.01和GS+Version7软件分别进行常规统计和地统计的条件随机模拟分析。 2 土壤电导率的空间变异性 2.1 表层土壤电导率与变异系数 按照经典统计学方法,对田间实测数据进行了基本参数的计算和χ2拟合检验,得到的结果列于表1。由表1可以看出,表层土壤电导率服从对数正态分布,其变化范围很大,最大值是最小值的20多倍。变异系数(Cv)反映的是相对变异,即随机变量的离散程度。根据相关研究,Cv≤0.1为弱变异性;0.1 Cv1为中等变异性;Cv≥1为强变异性。试验区土壤表层电导率变异系数为0.977属中等变异性。 2.2 空间分布函数模型 半方差图是研究土壤特性空间变异的一种有效方法。它表示了土壤区域化变量对于在不同位差分隔和方向上的所有成对点之间的观测值的空间相关性。表层土壤电导率半方差模拟模型及其拟合参数见表2,表层土壤电导率半方差图及其拟合结果见图2。从中可以看出表层土壤电导率的半变异函数符合指数模型: γ(h)={0C0+C(1?e?ha)h=0h＞0(1)γ(h)={0h=0C0+C(1-e-ha)h＞0(1) 其决定系数为0.809。在半方差函数模型中,块金值C0反映了区域化变量内部的随机程度,一般来说,C0为微观结构与采样及化验误差之和;C为结构方差,由土壤母质、地形、地貌、气候等非人为的结构性因素引起的变异;(C0+C)为基台值,表示系统内总的变异。C0/(C0+C)可以表示空间变异性程度(由随机性因素引起的空间变异性占系统总变异的比例)。如果C0/(C0+C)小于25%,表现为强空间相关性; 在25% ～75%之间的空间相关性中等;大于75%,空间相关性很弱;若比值接近于1,说明在整个尺度上具有恒定的变异。由表2可知,表层土壤电导率表现为强空间相关性,说明表层土壤电导率的空间分布主要是由结构性因素(如气候、母质、地形、土壤类型等)作用的结果。3a表示变程,即某变量观测值之间的距离大于该值时,则说明它们之间是相互独立的。若小于该值时,则说明它们之间存在着一定的相关关系, 由表2可知,表层土壤电导率的空间自相关距为5.91 km。 3 条件模拟和克立格插值 3.1 条件数据集合的模拟 序贯高斯条件模拟值法是贝叶斯理论的一个应用,是条件模拟值法中应用较广,理论成熟的方法之一. 此方法根据现有数据计算待模拟点值的条件概率分布,从该分布中随机取一值作为模拟现实.每得出一个模拟值,就把它连同原始数据、此前得到的模拟数据一起作为条件数据,