区域地下水位插值的整体-局部组合方法.PDFVIP

第 27 卷 第 8 期 农 业 工 程 学 报 Vol.27 No.8 2011 年 8 月 Transactions of the CSAE Aug. 2011 63 区域地下水位插值的整体-局部组合方法 1 2 1※ 3 王春颖 ，尚松浩 ，毛晓敏 ，Jorge Kazuo Yamamoto （1. 中国农业大学水利与土木工程学院 中国农业水问题研究中心，北京 100083 ； 2. 清华大学水利水电工程系 水沙科学与水利水电工程国家重点实验室，北京 100084； 3. Department of Environmental and Sedimentary Geology, University of Sao Paulo, Sao Paulo 05508-900, Brazil ） 摘 要：利用趋势面分析法进行空间插值时估值趋势面光滑，对局部数据起伏较大的区域插值精度不理想。针对这一问 题，该文基于非平稳地下水位数据，在用趋势面分析法进行空间插值的基础上，运用 Yamamoto 提出的修正的普通克里 金法对区域残差进行了空间插值。将趋势面分析法的估计值与残差估计值相加，得到区域地下水位的空间分布规律。结 果表明，采用这种整体与局部结合的组合插值方法，地下水位插值的结果比趋势面分析结果更接近实际情况，在保证全 局最优的同时较好地解决了估值趋势面光滑、局部插值精度低的问题，更能真实地反映区域地下水位空间分布特征。 关键词：地下水,插值，模型，趋势面分析法，普通克里金法，平滑效应 doi ：10.3969/j.issn.1002-6819.2011.08.011 中图分类号：P641.72 ；P628+.2 文献标志码：A 文章编号：1002-6819(2011)-08-0063-06 王春颖，尚松浩，毛晓敏，等. 区域地下水位插值的整体-局部组合方法 [J]. 农业工程学报 ，2011，27(8)：63－ 68. Wang Chunying, Shang Songhao, Mao Xiaomin, et al. Combined global-local interpolation method for regional groundwater level[J]. Transactions of the CSAE, 2011, 27(8): 63 －68. (in Chinese with English abstract) 波动较大的数据[4-5] 。 0 引 言 趋势面分析法（Trend Surface Analysis，TSA ）是分 空间属性（如地下水位）的插值指利用已知点的属 析地理系统、资源、疾病、人口、环境、农业等要素的 性值估计未知点的属性值，常用的空间插值方法包括反 空间分布规律与变化趋势的常用方法，通常可以取得较 距离加权法、克里金法、径向基函数法、趋势面分析法 好的插值效果[6-10] 。趋势面分析法利用不同阶次的趋势面 等。克里金法是一种线性无偏估计法，与其他的估计方 方程研究多个影响因素对因变量的作用，趋势图可以直 法相比，方差最小是克里金法的显著优点，其中简单克 观地分析因变量随自变量变动的趋势，可应用