等级秩序克里格法在位山灌区土壤墒情分析中的应用.pdfVIP

第 18卷第 12期 水利科技与经济 V01．18 No．12 2012年 12月 WaterConservancyScienceandTechnologyandEconomy Dec．，2012 等 级 秩 序 克 里 格 法 在 位 山 灌 区 土 壤 墒 ，睛 分 析 中 的 应 用 薄宏波 ，胡 健 ，曹洪博 ，许月勇 (1．位山灌区管理处，山东 聊城 252000；2．中国水利水电科学研究院，北京 100048) [摘 要] 通过大量的位山灌区田间实测土壤含水量资料，采用等级秩序克里格插值方法研究 土壤水分的空间分布格局。由于该方法存在平滑效应，使得等级秩序克里格的估计结果无法 与样本数据所表现出来的统计特征和空间分布结构相一致。采用一种后处理算法来修正等级 秩序克里格的平滑效应，并将其应用于位山灌区农田土壤水分的空间估计。结果显示，该方法 能够还原含水量样本数据的真实空间变异特征 。 [关键词] 土壤水分；空间变异；等级秩序克里格法；平滑效应；位山灌区 [中图分类号] $27 [文献标识码] B [文章编号] 1006—7175(2012)12—0092—02 1 概 述 2 等级秩序克里格法及其平滑修正 土壤水分的分布具有高度的空间异质性 ，了解和认 等级秩序克里格和普通克里格类似，只不过所有的 识土壤水的空间变异规律和分布格局对许多水文问题、 分析和计算都是在经过均匀变换后 的尺度下进行的。 农业管理问题的解决有着重要的实际意义。地统计学为 Yamamoto在 2005年初次提出了普通克里格后处理算 研究诸如土壤水分、土壤盐分、污染物含量，这类区域化 法 j，前者在确定修正因子f时采用的是数值优化方法 ， 变量的空间变异和空间分布提供了一套非常有效的工具 后者在统计推导的基础上给出了最优修正因子f的理论 集，在诸多不 同种类的克里格法中又以普通克里格法的 计算公式 ，因而更加简洁直观。由于等级秩序克里格实 实践应用最为广泛。为了提高基于不对称样本克里格估 际上是在均匀分布尺度下的普通克里格插值过程，因而 计结果的精度 ，除了借助一些变异函数的稳健估计量之 可以将Yamam~o的后处理算法应用到等级秩序克里格 外，各种非线性克里格法也应运而生，如指示克里格、多 平滑效应修正中。将二者结合 ，等级秩序克里格法可以 元高斯克里格、对数正态克里格和等级秩序克里格等。 简单描述为：①将初始样本观测值按照升序进行排列；② 等级秩序克里格是一种非线性克里格插值技术，它首先 将初始样本观测值转换成相应的均匀分布值 ，实际上是 把不对称的样本数据变换成[0，1]上的均匀分布，然后对 经验累计频率的计算过程；③根据公式计算均匀尺度下 变换后的值进行变异函数的构建、普通克里格估计 ，最后 的样本变异函数，并构建相应的理论模型；④选择合适的 再经逆变换将估计值变换成原始尺度下的值…。该法可 邻域搜索策略，保证参与估计的样本数 目既不能过多也 以有效地克服样本的不对称分布对估计结果的不利影 不能过少，过多或过少的样点都会引起平滑效应的增大； 响，已经被应用到土壤重金属污染物以及土壤机械组成 ⑤通过交叉证实检验计算每个样本点的普通克里格估计 空间变异的研究中。等级秩序克里格法作为克里格方法 值，从而可以计算对应于每个样本点的插值标准差、插值 的一种同样也受到平滑效应的困扰 ，如何使其估计结果 标准差个数；⑥采用跟交叉证实检验过程相同的邻域搜 再现样本特征也是值得探讨的问题。本文的目的在于将 索策略进行普通克里格插值 ，计算每个待估点处的插值 Yamamoto的后处理算法应用到等级秩序克里格平滑效应 标准差个数估计值、均匀值估计值以及相应的插值标准 的修正中，结合位山灌区实际的农