Geophysics一种不连续变网格有限差分模拟技术.PDFVIP

Geophysics: 一种不连续变网格有限差分模拟技术 有限差分模拟计算是理解地震波传播规律的主要手段。区域网格剖分 需要考虑实际地球介质的速度变化。因为地下介质非常复杂，存在较大 的速度差异，通常采用不同尺度的网格进行剖分。对背景速度场使用较 大的网格剖分，对相对较小的速度异常体采用较小的网格。这样，既可 以控制模拟精度，又能够提高计算效率。但是，粗细网格过渡区需要使 用插值技术获得网格节点上的波场值，因此产生的插值误差较大，严重 干扰了真实地震波场的传播。 地质与地球物理研究所地球深部结构与过程研究室博士研究生范娜 和导师赵连锋副研究员与其合作者提出了一种高精度变网格有限差分方 法，在粗细网格过渡区采用旋转有限差分算子代替插值算子，完全消除 了插值误差。图1 是使用旋转有限差分算子的不连续变网格方法示意图 。这种方法兼具连续变网格方法的空间高精度和不连续变网格方法较高 的计算效率。图2 比较了基于统一的连续网格、使用旋转有限差分算子 和插值算子的不连续变网格模拟的波场中波形上的虚假反射波。在使用 旋转有限差分算子的不连续变网格的有限差分模拟的波形中看不到虚假 反射波。