应用Kriging算法实现气象资料空问内插 - Read.PDF

维普资讯 第 32卷第 2期 气 象 科 技 Vo1．32，No．2 2004年 4月 METEOROLOGICALSCIENCE AND TECHNOLOGY Apr．，2004 应用Kriging算法实现气象资料空问内插 刘 峰 (北京大学物理学院大气科学系 ，北京 100871；民航中南空管局气象中心，广州 510406) 摘要 使用 Kriging算法对海平面气压和气温等气象观测资料进行空间内插，实现数据网格化 。通过 GRADS绘制 等值线图与人工分析的天气 图以及 Surfer软件的Kriging方案的结果进行对 比，并对 比分析多个统计量，验证 自制 软件结果的好坏 。结果表明：对于气象物理量场 ，变异 函数应采用球面模型；对量级较大的海平面气压值应进行缩 小化处理 ；对于常规观测资料 ，球面模型半径应取 10个经纬距以上。 关键词 Kriging算法 空间内插 数据网格化 通过对各个观测值赋予不同的权重 w 后累加得 引言 到，即Y =∑WU。如果该点上的真实值为y， 通过观测场可以获得一定数量的气象观测值， i= l 这些观测值反映了气象要素空间分布的全部或部分 则误差大小是 =(y 一y)，即 =(∑w ，【 特征。根据观测值来估计未知点上的要素值，需要 i： l 通过空间内插来实现。空间内插对于观测台站稀少 — y)，该式展开为 = ∑∑Ww ，【 — 或台站分布不合理的地区具有重要 的实际意义 ，并 i= lJ= l 使得利用网格数据进行客观分析成为可能。郑 良 2 ：w ，【y+y。如果该点上没有布置观测场，则 i= l 杰…等在华东中尺度天气试验 中采用逐步订正客 没法知道真正的误差大小，所 以上式不可解。 观分析方法和逐次最优插值客观分析方法对观测资 1．2 通过距离相关量求权重 料进行网格化。这两种方法需要先给出一个初始估 空间统计方法认为估计值与观测值的相似程度 计场 ，然后利用观测场与初始估计场进行订正。由 是通过点对的距离相关量来度量的。点对距离相关 于可能需要订正多次才能得到较好 的精度 ，因而计 量只与观测场间的相互距离有关 ，而与它们 的绝对 算量较大。Kriging(克里格法) 内插方法属于空 位置无关。气象观测场以及估计点的位置坐标是已 间统计方法(Geostatistics)，是由南非矿山地质工程 知的，因而可得到观测场之间、观测场与估计点之间 师 D．G．Krige于 1951年提出，并由法国地质学家 的距离相关量。距离相关量 c(【，，【，，)表示观测场 Dr．Matheron于 1962年加 以完善。Kriging是一种 i和J的距离相关关系，c(【，，y)表示观测场 i和 求最优 、线形、无偏的空间内插方法。在充分考虑观 估计点y的距离相关关系。 测资料之间的相互关系后 ，对每一个观测资料赋予 Kriging算法中使用距离相关量 c( ，【，)和 一 定的权重系数，加权平均得到估计值。 c( ，y)来代