面向大范围降水空间插值普通克里金模型开发与实例分析.pdfVIP

第22卷 1期 应用基础与工程科学学报 Vo1．22，No．1 2014年2月 JOURNALOFBASICSCIENCEANDENGINEERING February2014 文章编号：1005-0930(2014)01-0106-012 中图分类号：P333．6 文献标识码：A doi：10．3969／j．issn．1005-0930．2014．O1．Ol1 面向大范围降水空间插值的普通克里金 模型开发与实例分析 胡庆芳 ，胡 艳3，杨大文 ，王银堂2，葛 慧2 (I．清华大学水利水电工程系，北京 100084；2．南京水利科学研究院水文水资源与水利工程科学国家重点实验 室，南京210029；3．水利部太湖流域管理局，上海200434) 摘要：基于Matlab平台，开发了面向大范围降水空间插值的普通克里金模型一 MatlabbasedordinaryKriging(MatOK)．与已有模型相比，MatOK的主要特点是： (1)在降水空间变异函数计算环节，采用SCE—UA算法拟合理论变异函数；(2)在 OK方程组估值环节，引入降水空间发生概率的计算，完善了日等短时间尺度降 水空间估值方法；(3)通过对OK方程组的标准化处理，有效提高了模型数值计 算的稳定性．同时，将MatOK初步应用于集水面积为83374km 的赣江流域年、 月、日降水空间插值中，并重点讨论了MatOK的计算稳定性和计算效率，结果初 步说明了采用MatOK进行大范围区域降水空间插值是可行的． 关键词：降水空间插值；普通克里金；变异函数；数值稳定性；计算效率；赣江流域 降水是最基本的气象水文变量之一，其空间异质性深刻影响着地表产汇流及其伴生 过程．因此精确获取降水空间分布一直是水文领域一项重要的科学 目标．而基于雨量站实 测数据进行空间插值是降水由离散的点信息扩展为连续空间信息的一种基本方法．普通 克立金法(OrdinaryKriging，OK)是一种常用的降水空间插值方法，同时也广泛应用于土 壤水分、矿产储量等众多变量的空间估计中．该方法的思想精髓是 “全局的空间变异函数 计算，局域的线性无偏最优估计”．尽管OK提出已有较长时间，但 目前针对该方法的研究 和应用仍比较活跃．在OK主要思路的基础上，近年来一些学者又进行了改进和发展．如 Ma和Yao提出了利用快速傅立叶变换拟合空间变异函数的方法 J，Harris等提出了基于 地理加权变异函数的移动窗口克里金方法 ；Hengl和Guan各 自将回归分析和OK结合 起来发展了回归．克里金方法 ，Yamamoto提出了基于 “插值方差”修正OK平滑效应的 途径 等． OK的计算过程比较复杂，为便于应用，国外学者开发了一系列模型或程序，其中比 较著名的有GSLIB ]、DACE 等．然而降水在空间上具有连续性与间断性、确定性与随机 性的双重属性，不仅具有量的不同，而且可能具有 “干湿”状态的不同，但已有模型在插值 收稿 日期 ：2012-08-28；修订Et期：2013-0729 基金项目：国家自然科学基金项 目；水利部公益性行业科研专项经费项 目(201001002) 作者简介：胡庆芳(198l一)，男，博士，工程师．E-mail：hqfwork@163．corn 应用基础与工程科学学报 是否大于某一临界降水发生概率(criticalprecipitationocculTenceprobability，POPcr)，再估 算降水量．这样的优点是在兼顾了降水空间插值的定量估计精度和分类估计精度．与 Thornton的方法不同的是，MatOK无论是在降水发生概率的估计还是在降水量的估计两 个环节均采用OK，而不是采用距离反比方法．具体而言，MatOK在空间变异函数计算模 块，同时计算降水量 自身的理论空间变异函数 (称之为 “一般理论变异函数”)和降水量对 应的指示变量的变异函数(称为 “指示理论变异函数”，指示变量的定义如式 (1)，其中Pi 指某位置的降水量)．在OK方程组估值模块，先根据指示理论变异函数估计POP，再根据 一 般理论变异函数估计降水量．另外，在求解两个 OK方程组前，对MatOK原始方程组进 行标准化处理，以保