石文插值算法参数使用说明1.docVIP

石文插值算法参数使用说明 本软件提供有反距离、可里金、多面函数、移动拟合内插法、Shepard方法和最小曲率等插值方法。为了方便广大用户使用，现将各方法适用情况及参数设置对插值影响做以说明。 一、 每一种插值方法都有的参数说明： 参估点数、最大距离： 插值就是在一个未知结点周围找已知的数据点，并用找到的已知数据点估计出这个未知节点的值。 图1 用所有数据点插值结果 图2 用参估点数为30时插值结果图 这就涉及到要在未知节点周围找多少个已知数据点，“参估点数”是给定的最大参估点的个数，在位置点周围最多选这么多点数（可能会少选）。 一般，插值的权值是两点距离的函数。已知点和待估点之间距离越大，已知点对待估点的贡献越少。“最大距离”表示当已知点和待估点之间距离大于或等于这个距离时，已知点对待估点的贡献为0。因此超过“最大距离”的数据点不会被选入。 图1和图2为同一组数据插值的结果，图1使用所有数据，图2使用了部分数据，两图趋势差别不大，但用所有数据时插值结果比较光滑。 在软件中，若已知数据点个数小于800，参估点数等于数据点数；若大于800，参估点数选20~50之间的数（为了提高插值速度）。当数据量特别大时，如地震数据，参估点数可选8~20之间的数，若数据点分布均匀，不会影响插值结果。“最大距离”软件默认为1.0e+6，几乎包含所有数据点，若用户知道实际数据范围，可自行设置，但要尽量保证选够参估点数。 各向异性参数 各向异性参数包含各向异性比和方位角两个参数。该参数在软件的可里金插值方法与其他插值方法中规定不尽相同，但不影响使用。最新版本已统一参数设置标准，使其符合地质研究要求。 图3 方位角为0，各项异性比为3的插值结果 各向异性比是以纵轴为标准，即以南北方向为标准。若各向异性比为R，那么在南北方向一个单位距离引起的属性的变化量，相当于在东西方向（横轴）经过R个单位距离，属性所产生的变化量。 方位角以北方（纵轴正向）为基准，顺时针旋转为正方向。 图4 方位角为20０，各项异性比为3的插值结果 　　　　　图5 方位角为60０，各项异性比为3的插值结果 若各向异性比取值为1时，调节方位角插值结果一般不会变化。可用简单数据给不同参数试验几次比较一下结果，了解具体的参数设置。 二、插值方法说明 1、反距离，Shepard方法 反距离是一种常用的距离插值方法，直接以距离远近控制已知数据对待估节点的贡献率（确定权值），但插值结果对距离较敏感，易产生“牛眼睛”现象。Shepard方法是该方法的一种改进，在某种程度上克服了“牛眼睛”现象。 反距离方法有参数“指数”，取值为1，以欧氏距离为标准度量；取值为2，以欧氏距离平方为标准度量；取值为3，以欧氏距离立方为标准度量； 克可里金方法（简单克可里金和普通克可里金） 克可里金法是地质上最常用的插值方法，因而在本系统中设置该方法为系统默认方法。 两种方法参数相同，有块金值、拱高、变程和稳定参数等。一般块金值与拱高之和等于数据方差。大小变程之比等于各向异性比。调整这些参数相当于改变变差函数，本软件在克里金方法中给了7个不同的变差函数，通常这些参数不需要调整。 由于插值时需要求解克里金方程，有可能方程病态或是矛盾方程，通过设置稳定参数可解决这个问题，也可通过改变变差函数或修改参数（大变程，拱高等）解决这个问题。 软件默认的变差函数为线性变差函数，有断层时默认的变差函数为指数变差函数。有断层时用指数变差函数效果比较稳定。同样的数据使用不同变差函数，插值结果不尽相同，使用有些变差函数可能会出现错误结果。遇到这种情况，可改变变差函数，调节参数等以改善插值结果。系统所给的默认参数是经过大量试验，在大多数情况下能得出正常结果的。 多面函数 该方法有一个参数“核值”，核值取值在0~1之间。一般随着核值的增大，等值线更加平滑，偏向于描绘趋势。 移动拟合内插法 该方法是一种拟合方法，注重描绘属性的变化趋势。该方法有一参数“指数”，它表示采用的拟合多项式最高次幂，它的取值一般为1~5，系统默认为2。用户可根据需要调节该参数。 最小曲率 最小曲率方法所得等值线光滑，变化自然。该方法有几个特殊的参数：松弛因子{纵向松弛因子，横向松弛因子}，误差和迭代次数。“松弛因子”控制等值线的平滑程度，系统默认为3。松弛因子值越大，等值线越平滑。 “误差”是前后两次迭代结果的差值，当前后两次迭代结果的