区域土壤质地层次三维空间分布的地统计模拟方法比较.pdfVIP

土壤 (Soils)， 2010, 42 (3): 429-437 区域土壤质地层次三维空间分布的地统计模拟方法比较① 贺 勇 l ， 胡克林 1· ， 李保国 1 ， 黄元仿， , De1i Chen2 (1中国农业大学土壤和水科学系，教育部槌物土壤相宜作用重点实验室，北京 100193; 2 澳大利亚锺尔本大学土地和食品学院，澳大利亚维多利亚 3052) 摘要2 为了~I:化描述冲权平原阪土攘倒回原地层次的复杂变异性，本研究根据华.:f~)fft积平原…块面积约 IS km2 阻域 内的 139 个土攘剖面 (0 - 2 m) 观测敛掘，分别应用顺序指示模拟和多点地统计学模拟方法分析了区域内土攘质地层次的三维 费问分布特征.前者的结泉显示土填质地居在水平方向上的分布具有很强的随机特性，而后者的结果表明土壤质地层在水平方 向上的分布具有一定的规律性.将转移概率函数得到的各向阅伎分布比例和各向异做平均氏度(厚度〉用米线成立维训练、倒像， 然后将生成的训练图像用于多点地统计学模拟，并与传统的顺序指示模拟结果进行了对比.结果表明，多点地统计学模拟在估 计精度以及法直方向上土壤质地层次连贯性的再现性上优于传统顺序指示模拟，所得纺果能更好地反映土壤质地层次的空间分 布特征. 关键词: rK域尺度:工t编质地层次;二三维:顺序指示模拟t 多点地统计学模拟 中图分类号 : SIS2.7 土壤质地层次的空间分布十分复杂。典型的土壤 制面只能代表有限的眩域，很难将其扩展到3巨大的研 究区域[11. 在传统的土壤分类制图过程中，土壤制阁 者常根据其经验将有限的土壤观测剖面扩展到整个研 究眩，然而，当外界成土环境条件复杂多变导致研究 区的土壤质地层次空间变异较大时，土壤制图者需要 消耗大量的人力物力未获取更多的土壤剖面资料，这 常常是不切实陈的.由于土壤质地层次在土壤科学和 农业生产实践中具有很霆耍的意义[21，因此，有必要 探求精确定量化的方法来描述土壤质地层次的空间分 布. 目前，用于模拟指示变量〈如岩石屈就质地摆〉 的主要方法是基于变异函敬的两点地统计模拟(如顺 序指示模拟) [3η. .这种方法不仅考虑到了样点数据 的空间结构特征，而且能体现未来样点的空间变异性 [81. 因此，广泛地用于 2-D 土壤制罔(9)和13-。土墙性 质模拟(10- 1 1].然而，对于具有较高连贯性的土壤变最 (如土壤质地层)，基于变异函数的条件模拟结果常表 现出高捕或低连贯性[12-13).而且，当数据稀疏时，变 异函数常表现为纯块金效应成具有很高的随机特性 (8， 13- 141. 这种情况下进行条件模拟，无法体现变异函数 的优点，此外，传统的变异踊数撤然可以反映空间上 网点之间的相关性，但难于表征复杂的空间结构和再 现复杂目标的几何形态(15) ? 鉴于上述传统建模方法的不足，近年米一种用于 模拟地质结构的新方法一一多点地统计学模拟，得以 研究和发展(15 叫。这种方法在一定程度上能够克服传 统条件模拟中出现的高精或低连贯性问题I叫。但由于 其算法上的缺陷，未能在拱际中推广应用 [15， 191. 加ebelle 等(161在前人研究的基础上将算法加以改进， 提出了多点模拟的 SNES (single normal equation simulation) 算法.国内外很多学者应用该方法进行了 实际建模，井证明了多点地质统计模拟优于传统地质 统计学建模方法[15吨 16， 19- 201. 然而至IJ 目前为止，将此方法 用于土壤剖面质地层次模拟的研究还少见报道。 因此，本文选取了华北冲积平腺~块面积的 15 km2 的区域，分别运用指示变异函数和转移楠率函数 分析土壤质地层次在不同方向上的空间变异特性;在 此基础上，分别运用顺序指示模拟和多点地统计学模 拟方法对土壤质地层次进行王维模拟:通过对模拟结 果的检验，评价这两种方法在区域土壤质地层次模拟 的适用性， ①J事金项目 z 阂署~.点荔础研究发展规划项目。∞9CB118ω川、国家自然科学基金项目 (40401归川和教育部新世纪优秀人才交待计划项目 (NCET.()7-0809 )资助. * lSH孔作者 (hukel@ωu.edu.cn) 作者简介s 贺勇 (1980…).男，内盖世宙人.博士研究生，主要从事土地利用与倍息技术的研究. E-mail: hegaol叭@163.ωm 430 土壤 第 42 卷 1 材料与为法 1. 1 数据采集 研究区位于河北省曲周县…块面积约 15 旷的 区域，采取平行四边形闷格化布点 (GPS ~伎)，四jiJ. 形边长 350m，共布设 139 个点(圈。。采取开挖剖 面和钻孔制面相结合的方法，观测土壤削面质地层次 的变化。观测深度为 2 m，质