Tricubic插值在三维波前构建法中的应用.pdfVIP

·1174 · 中国地球科学联合学术年会 2014 Tricubic 插值在三维波前构建法中的应用 * 高正辉 孙建国 吉林大学地球探测科学与技术学院 长春 130026 三维波前构建法在 Kirchhoff 型偏移，层析成像，观测装置设计等方面扮演着重要的角色。相比于 其他只计算初至走时的射线追踪方法，该方法有一个明显的优势，即能计算多值走时。但在三维波前构 建法中，模型的参数给定在规则的立方体网格点上。在计算的过程中，由射线端点组成的三角波前网格 是不规则分布的，所以射线端点不太可能落在规则的立方体网格点上。因此，需要插值来得到射线端点 处的速度及其导数。 目前主要有两种插值方法来处理上述问题：三线性插值和基数样条插值。前者具有非常快的计算速 度，但是精度相对较差。后者在一定程度上提高了计算精度。但是，这两种方法都是由一维的插值公式 拓展而来，并不是本质上的三维和各向同性（插值结果与插值方向的先后次序无关）。 本文在三维波前构建法中首次引进 Tricubic 插值方法来计算插值点的速度及其导数，它是一个本质 上 三 维 的 各 向 同 性 局 部 插 值 方 法 。 采 用 插 值 点 周 围 八 个 网 格 点 的 函 数 值 及 其 导 数 ¶f ¶f ¶f ¶ 2f ¶ 2f ¶ 2f ¶ 3f {f , , , , , , , } 来作为限制条件。看起来需要太多的参数，但是近些年来， ¶ x ¶y ¶ z 抖x y 抖x z 抖y z 抖x y z 振幅信息越来越重要。为了获得振幅信息，必须进行动力学射线追踪，这就需要计算速度的二阶导数。 ¶ 3f 由此，只需要在每个立方网格点处额外地计算 。 抖x y z 3 i j k Tricubic f x , y ,z = a x y z f x , y ,z 插值由 Lekien 和 Marsden （2005 ）给出，其形式为： ( ) å ijk 。式中， ( ) , , 0 i j k = 表示插值结果，aijk 表示系数矩阵的元素，x ,y ,z 表示在一个单位立方网格内的相对位置。因为需要确定 2 2 2