基于PER-Kriging插值方法的降水空间展布 Spatial interpolation of precipitation using the PER-Kriging method.pdfVIP

第22卷第6期 水科学进展 Voi．22，No．6 201 i年11月 ADVANCESINWATERSCIENCE Nov．．20II DOI：CNKI：32．1309．P．20ll1125．1638．019 基于PER—Kriging插值方法的降水空间展布 王 舒L2”，严登华2”，秦天玲2’3…，刘汝海1 (1．中国海洋大学海洋环境与生态教育部霞点实验室，山东青岛266100；2．流域水循环模拟与调控国家重点实验室， 摘要：将PRISM插值对高程因素的计算方法引入普通Kfiging插值中，提出了综合考虑观测点与插值点的位置、距 离和高程关系的PER—Kriging插值方法。通过引入计算参考平面并进行数学变换，使该方法具有较好的可操作性。 在对高程～降水量采用线性拟合的基础上，加入非线性拟合，进一步提高PER-Kriging插值的表达能力。在澜沧江 mm以上；不 地形差异太大造成的插值异常，并且能较好地消除插值运算的平滑效应，平均误差比后者减少了20 J『司曲线拟合插值的相对误差均小于12％。 关键词：PER-Kriging插值；降水；高程；澜沧江 中图分类号：P426．6；TVl25文献标志码：A 文章编号：1001-6791(2011)06-0756-08 由于全球气候变化对自然环境与社会经济的影响不断加剧，水文水资源、生态与环境、资源科学等现代 科学研究中经常需要使用高精度的降水空间分布数字信息¨…。然而由于研究区地形地貌、技术手段等多方 面因素的综合限制，观测点的设置较为有限，其布局也往往不能满足研究的需要。利用分散的观测点数据进 行空间插值展布已经成为获取高精度降水空间分布数字信息的主要方法之一。现阶段，常用的插值方法包括 响，优于其它线性插值方法…，但是仍然存在平滑效应，在地形起伏比较大的区域效果并不理想”引；而 PRISM插值则根据高程与插值要素之间的关系进行插值，当高程与插值要素相关系数较高时能取得较好的 插值效果，然而其输入参数较多，各要素的权蓖计算方法不确定…。 Regres- sionsofKriging)插值方法，并在澜沧江流域对该方法进行了验证。 1 PER-Kriging插值方法的建立 1．1普通Kriging插值 普通Kriging方法将地统计学运用到插值运算中，利用半变异函数来进行无偏估计内插，应用范围非常 广泛。普通Kriging方法对每一个插值点_『的插值计算公式如下h3川： 彳=∑AiZ。 (1) i=l 收稿日期：201l-03·17；网络出版时间：201l-11．25 网络出版地址：http：／／www．enki．net／kems／detail／32．1309．P．20111i25．1638．019．html ■ 目 作者简介：王舒(1983一)，男，云南昆明人。硕士研究生，主要从事生态水文学研究。 C11．COrn E-mail：wangshul944@21 万方数据 第6期 王舒，等：基于PER-Kriging插值方法的降水空间展布 757 式中z?为插值点．『降水量的计算值；n为参与插值计算的观测点个数；Z。为第i个观测点的观测值；A。 为第i个观测点的权重。 普通Kriging插值采用半方差函数作为确定权苇参数Aj的依据，同时考虑了距离和方向因素，因此可以