利用DEM作为辅助信息推定大区域地下水初始流场.PDF

水 　　利 　　学 　　报 　2004 年 11 月 SHUILI 　　XUEBAO 第 11 期 ( ) 文章编号 2004 利用 DEM 作为辅助信息推定大区域地下水初始流场 1 ,3 ,4 2 1 ,3 3 朱奎 ,张祥伟 ,夏军 ,牛存稳 ( 1 中国科学院地理科学与资源研究所 陆地水循环与地表过程重点实验室,北京　100101 ; 2 水利部规划计划司 ,北京　100053 ; 　3 武汉大学 水资源与水电工程科学国家重点实验室 ,湖北 武汉　430072 ; 4 中国科学院 研究生院,北京　100101) 摘要 : 由于大区域地下水位实测资料比较少,信息不足 ,影响水位空间分布的推定精度 , 因此 ,充分利用辅助信息弥 ( ) 补大区域实测水位信息不足 ,提高地下水位的推定精度十分必要。本文提出了利用 DEM 数字地形模型 中的地形 ( ) 标高作为辅助信息修正地下水位趋势面的 RK Residual Kriging 法 ,并将此方法运用到华北平原滏阳河平原区。推 定地下水位空间分布的过程中 ,利用 DEM 信息的趋势方程 ,复相关系数较高 ,残差值范围小、更好的满足正态分布 要求。对推定的水位结果进行分析 :单个验证点的地形、位置对趋势方程的影响很大 ,但验证点的三种平均误差则 表明 ,利用 DEM 信息的修正 RK 法得到的推定结果误差较小 ,且三类方程采用一次趋势面效果更佳。 关键词 :大区域地下水 ;初始水位 ; Kriging 法 ;DEM 中图分类号:P6412 文献标识码 :A 研究大区域地下水资源 ,首要问题是对区域地下水流动进行符合实际的模拟 ,掌握区域地下水的运 动变化规律。采用数值计算方法对地下水流动模拟时 ,需要根据实测地下水位值来推定初始地下水流 场 ,初始流场的推定直接影响着整个地下水资源的评价精度。初始流场的推定通常采用 Universal ( ) [1～2 ] [3 ] [4 ] Kriging UK 方法 ,但由于 UK 法的计算行列式较大 , 比较复杂;Neuman 等 和 Sun 提出了 Residual ( ) Kriging RK 法推定区域地下水位 ,其基本思想是将地下水位的趋势成分过滤后得到实测点的水位残差 , ( ) 然后运用 Ordinary Kriging OK 法推定水位残差的空间分布 ,趋势成分加推定的残差分布得到地下水位 的空间分布。RK 法计算简单 ,便于理解 ,但是该方法所需实测资料较多 ,当实测水位点密度较小时 ,计 算的趋势面方程就会与实际地下水位面相差很大 ,从而影响水位的推定精度。为弥补实测资料不足 ,陈 [5 ] ( ) 家军 将地形标高作为辅助信息。运用 CoUniversal Kriging CK 法推定大区域地下水位空间分布 ,提高 了推定精度。然而 ,该法需要识别地下水位和地形标高的交叉实验变差函数 ,经过多