东北三省月降水量的时空克里金插值研究 Study on Spatial-temporal Kriging Interpolation of Monthly Precipitation in Three Provinces of Northeast China.pdfVIP

第31卷第3期 水文 Vol_31No．3 2011年6月 OFCHINA JOURNAL HYDROLOGY Jun．．2011 东北三省月降水量的时空克里金插值研究 李 莎12，舒红1，徐正全1 (1．武汉大学测绘遥感信息工程国家重点实验室，湖北武汉430079； 2．湖北第二师范学院机电系，湖北武汉430205) 摘要：为了实现对时空场任意点进行时空插值，以东北三省1999-2008年的月降水数据为研究对象，选 用了一类比较实用的积和式时空变异函数模型进行时空克里金插值。每个站点的降水量都是时间序 列，插值之前对各站点的降水进行了时序分解和去季节项处理。在分别得到纯空间和纯时间的变异函 数基础上建立时空变异函数。将二维的普通克里金插值扩展为三维，并对1999年10月各站点进行降水 量估计，同时与单纯空间克里金插值效果进行比较。结果表明时空插值效果理想，因为同时考虑了空 间、时间的相关性．插值精度较空间克里金更高。 关键词：克里金插值；时空相关；变异函数；降水量 中图分类号：P333 文献标识码：A 文章编号：1000—0852(2011)03—0031—05 降水量是气象数据中十分重要的元素之一。但是 流域进行降水插值分析．其结果证明对于月降水的插 由于人力、物力及地理条件的限制，气象站点的设置有 值克里金方法比较优越[31：孙然好等运用协克里金法 限且分布不均匀，同时站点所观测的降水量为点降水 对内蒙古降水量进行插值研究[41；魏义长等将经验正 数据，而要获得观测区域内面降水量和其它位置的降 交函数和地统计学方法结合起来对降水进行分析和空 水数据则需要通过插值来实现。关于空间数据插值的 间插值[51；孙鹏森等在大地坐标位置、海拔高度基础 方法较多，如反距离加权法(IDW)、泰森多边形法、样上，增加季风方向与坡向耦合效应的指标辅助进行区 条(Spline)插值法、径向基函数法(RBF)、克里金插值域降水的空间模拟和插值[61：刘智勇等基于DEM数据 选取反距离加权法、张力样条函数法、普通克里金法及 法(Kriging)等，这些方法各有特点也有各自的应用领 域．然而上述方法均没有考虑数据的时间维特性和时 协同克里金法对榆林市降水量进行空间插值同。 空间的相关性，只能进行空间插值却不能时间插值。本 l时空插值原理 文以东北三省1999～2008年的月降水数据(mm)为研 究对象，在地统计学的克里金方法基础上进行时空扩 许多自然现象的样本数据，如降水量的分布，在空 展，同时考虑空间域与时间域的相关特性，对研究区域 间或时问上表现出一定的相关性，即近距离的样本测 的降水进行时空插值。 量值之间的相似度往往比远距离的测量值之间更高。 对于降水数据的分析，国内外研究很多。Cowpert—克里金插值方法正是基于随机变量在空问或时间上的 wait运用Neyman—Scottrectangularpulses(NSRP)模相关结构而建立的。 型拟合降水数据并对缺失数据进行插值估计【1]： 1．1 平稳假设 Wheater等提出了基于泊松簇过程的单点模型用于降 平稳假设是进行统计分析的重要前提，地统计学 水模型拟合[21：刘登伟等选用了距离平方反比法 有两类平稳假设：二阶平稳假设和同有假设．两者均针 (IDS)、梯度距离平方比和普通克里金法分别对海河对于空间变量而言，将其扩展至时空场，如下： 收稿日期：2010—09—15 基金项目：半干旱气候变化教育部重点实验室(兰州大学)开放基金项目；国家高技术研究发展计划项目(2008AAl22201) 作者简介：李莎(