弹性波理论.docVIP

地震波交错网格高阶差分数值模拟研究 摘要: 地震波数值模拟技术是勘探地球物理学中的重要组成部分，研究通过弹性波一阶速度——应力方程，采用交错网格高阶有限差分法实现了地震波在各向同性介质中的高精度的数值模拟，并采用完全匹配层( PML) 吸收边界来消除边界反射，可取得较好的效果。通过模型的正演计算和复杂模型的处理结果表明，交错网格高阶有限差分法数值模拟是一种快速有效的地震波数值模拟方法。 关键词: 地震勘探; 交错网格; 有限差分; 数值模拟 引 言 地震数值模拟是模拟地震波在介质中传播的一种数值模拟技术，随着地震波理论在天然地震和地震勘探中的应用，地震模拟技术便应运而生，并随着地震波理论和计算机技术的发展，地震数值模拟技术自20世纪60年代以来也得到了飞速发展，形成了目前具有有限差分法、有限元法、虚谱法和积分方程法等各种数值模拟方法的现代地震数值模拟技术。 有限差分法是偏微分方程的主要数值解法之一。在各种地震数值模拟方法中，最早出现的数值模拟方法是有限差分法。Ａlterman和Karal(1968)首先将有限差分法应用于层状介质弹性波传播的数值模拟中。此后，Boore(1972)又将有限差分法用于非均匀介质地震波传播的模拟。Alford等（1974）研究了声波方程有限差分法模拟的精确性。Kelly等（1976）研究了用有限差分法制作人工合成地震记录的方法。Virieux(1986)提出了应用速度——应力一阶方程交错网格有限差分法模拟P——SV波在非均匀介质中的传播。交错网格方法提高了地震模拟的精度和稳定性，并消除了部分假想。 有限元法也是偏微分方程的数值解法之一。Lysmer和Drake（1972）最早将有限元法应用于地震数值模拟。Marfurt(1984)研究对比了模拟弹性波传播的有限差分法和有限元法的精度。Seron等（1990,1996）给出了弹性波传播有限元模拟方法。Padovani等（1994）研究了地震波模拟的低阶和高阶有限元法。Sarma等（1998）给出了三维声波模拟的虚谱法。 积分方程法是建立在波动方程的积分表达式的基础上的，其理论基础是惠更斯原理。积分方程法也是有限元法之后发展起来的一种地震数值模拟方法。Pao和Varatharajulu（1976）提出了弹性波散射的积分表达式。Bennett和Mieras(1981)给出了流体目标声波散射的时间域积分方程解。Bouchon(1987)给出了裂隙或孔洞弹性波绕射的离散波数法模拟方法。Bouchon等（1989）研究了具有不规则界面的多层介质中波传播的边界积分方程——离散波数法。Bakamjian（1992）给出了三维地震波传播模拟的边界积分方程法。符力耘和牟永光（1994）提出了弹性波正演模拟的边界元法。符力耘等（1997）提出了非线性Fredholm积分方程的正演问题。符力耘（2003）给出了含起伏地表的广义Lipmann—Schwinger积分方程的数值模拟方法。 射线追踪方法是建立在波动方程的高频近似基础上的一种地震数值模拟方法（cerveny等，1977）。这种方法实际只计算了最奇异部分的解，即旅行时和振幅函数的特征曲线，它们分别是程函方程和传播方程的解。这种方法计算效率高。但是，一些复杂的本构方程由于积分方程法和射线追踪法不满足假设条件而限制了这些方法的应用。 上述这些地震数值模拟方法各有优缺点。对于复杂构造、复杂地质体和复杂岩性地震模拟而言，交错网格高阶有限差分法其综合性能（占内存大小、模拟精度、计算效率和并行算法实现）最好，是实用性最好的方法。 在声波方程正演的数值模拟中,由于有限的计算空间区域无形之中引入了人为边界,不可避免的需要对在数值网格边界上产生的反射或回绕能量进行合适的处理,否则这些人为产生的反射或回绕能量会在很大程度上扭曲真正的波动传播信号,使模拟剖面变得模糊不清,不利于对地层构造信息进行解释。 为了消除这些人为产生的边界反射或回绕能量,人们发展了多种方法,其中最常见的主要有以下几种:一是最简单的扩展边界法,即在需要计算的数值网格外增加一些额外的网格数目,这样可以使人为边界反射效应远离所需要的计算网格,但是这种方法带来的负面效应是所需要的网格数目大大增加,因而也大大增加了计算量,对计算机的计算速度和存储能力提出了更高的要求,所以这种方法并没有得到很好的推广;二是海绵吸收法,即在计算网格的边界区域设置一定宽度的阻尼带,利用某些衰减函数对数值模拟波场进行逐步衰减;三是反周期扩展法,即利用正反周期函数极性相反的特点消除回绕波场;四是傍轴近似法,即利用波动方程的傍轴近似条件来消除计算边界上的反射。如何选择合适的边界吸收是一个值得研究的问题。 数值模拟基本原理 各向同性介质是最基本的一种介质模型，目前地震勘探中大多都是基于这种介质模型，根据