王美琪 李 建

(西南石油大学计算机科学学院 四川 成都 610500)

0 引 言

克里金插值法是一种在地统计学中常用的空间插值方法[1]。该插值方法最优无偏，效果优于距离加权法和趋势面法，目前已应用于许多领域，几乎是20世纪中后期以来应用最为广泛的地理空间插值方法[2]，如地质学、土壤学、生态学、气象学、遥感甚至其他研究“时空变量”的领域[3-5]。克里金插值法主要分为两大部分：(1) 量化已知数据的空间相关性即变差函数，一般来说，距离较近的点属性值越接近，距离较远的点属性值相差较多；(2) 确定邻域范围并搜索邻域点，确定克里金方程组，求解克里金方程组得到权重系数，利用权重系数加权求和得到待插值点的属性值。

在石油行业中，通过提取三维地震目的层的属性数据[6]，计算得到属性数据的分布平面图，如振幅、频率、属性聚类等含油气预测平面图，为找油找气提供目标方向。三维地震属性平面图的形成需要插值，由于三维地震属性数据巨大，在实际应用中难以全区域计算，一般采用局部邻域为每一个待插值点搜索出合适的邻域点作为输入。目前邻域点搜索多采用指定半径的圆或椭圆，将在圆或椭圆内的点作为邻域点输入，该方法须遍历，并计算所有已知点到待插值点的距离。文献[7]漂移模式和搜索范围参数的影响较小,可忽略不计；文献[8]改进了邻域克里金插值法；文献[9]采用Delaunay三角划分搜索和固定距离的方式搜索点；文献[10]引入KD-Tree的二维索引结构，提出改进J邻近点搜索方法，效率有所提高，但该类搜索仍需计算已知点到待插值点的距离，这无形中增加了计算量。……