地震波数值模拟方法研究综述

地震波数值模拟方法研究综述

在地学领域，对于许多地球物理问题，人们已经得到了它应遵循的基本方程(常微分方程或偏微分方程)和相应的定解条件，但能用解析方法求得精确解的只是少数方程性质比较简单，且几何形状相当规则的问题。对于大多数问题，由于方程的非线性性质，或由于求解区域的几何形状比较复杂，贝y 不能得到解析解。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但这种方法只是在有限的

情况下是可行的，过多的简化可能导致很大的误差甚至错误的解答。因此人们多年来寻找和发展了另一种求解方法一- 数值模拟方法。

地震数值模拟(SeismicNumericalModeling) 是地震勘探和地震学的基础，同时也是地震反演的基础。所谓地震数值模拟，就是在假定地下介质结构模型和相应的物理参数已知的情况下，模拟研究地震波在地下各种介质中的传播规律，并计算在地面或地下各观测点所观测到的数值地震记录的一种地震模拟方法。地震波场数值模拟是研究复杂地区地震资料采集、处理和解释的有效辅助手段，这种地震数值模拟方法已经在地震勘探和天然地震领域中得到广泛应用。

地震数值模拟的发展非常迅速，现在已经有各种各样的

地震数值模拟方法在地震勘探和地震学中得到广泛而有效

的应用。这些地震波场数值模拟方法可以归纳为三大类，即几何射线法、积分方程法和波动方程法。波动方程数值模拟方法实质上是求解地震波动方程，因此模拟的地震波场包含

了地震波传播的所有信息，但其计算速度相对于几何射线法要慢。几何射线法也就是射线追踪法，属于几何地震学方法，由于它将地震波波动理论简化为射线理论，主要考虑的是地震波传播的运动学特征，缺少地震波的动力学信息，因此该方法计算速度快。因为波动方程模拟包含了丰富的波动信息，为研究地震波的传播机理和复杂地层的解释提供了更多的佐证，所以波动方程数值模拟方法一直在地震模拟中占有重要地位。

1地震波数值模拟的理论基础

地震波数值模拟是在已知地下介质结构的情况下，研究地震波在地下各种介质中传播规律的一种地震模拟方法，其理论基础就是表征地震波在地下各种介质中传播的地震波传播理论。上述三类地震波数值模拟方法相应的地震波传播理论的数学物理表达方式不尽相同。射线追踪法是建立在以射线

理论为基础的波动方程高频近似理论基础上的，其数学表形式为程函方程和传输方程。积分方程法是建立在以惠更斯原理为基础的波叠加原理基础上的，其数学表达形式为波动方程的格林函数域积分方程表达式和边界积分方程表达式动方程数值解法是建立在以弹性或粘弹性理论和牛顿力学

为基础的双曲型偏微分方程一波动方程的理论基础上的。由于地下介质性质不同，其相应的地震波传播方程也不同。

由于地震波波动方程在复杂介质中地震波传播研究的

广泛适应性及地震波方程数值解法在地震波数值模拟中应

用的广泛性和有效性，本文将重点研究地震波波动方程数值解法，同时对几何射线法和积分方程法也作适当的讨论。

2地震波数值模拟的内容及特点

2 . 1地震波数值模拟的内容

地震波数值模拟是以地震波传播理论为基础的。描述地震波在各种介质中传播的波动方程是一个变系数的偏微分方程。地震波波动方程的定解问题(即正演问题)包括微分算子、算子系数、震源、初始条件和边界条件等。地震波正演过程中，求解微分方程可以计算出系统中表示状态的参数随时间的变化。地震波的正演过程数学上可描述为

d=A(m)

式中d为合成地震数据向量；

以为正演算子；

m为模型向量。

d的精度受m的离散化精度和正演算子以计

算精度的影响。

地震波波动方程正演问题的内容主要包括：

(1)地震波数值模拟的基本原理；

(2)地震波数值模拟的算法；

(3)程序设计及其质量，它主要受计算精度、计

算效率和计算稳定性的影响。

2. 2模型的离散化

研究的目的不同，构成地球物理模型的物理量

也不同。对于均匀各向同性介质中的声波方程而言，地球物理模型小可以表示为m=(p,v P,)；而对于均匀各向同性介质

中的弹性波方程而言，其地球物理模型则可以表示为m —m(p , vp,vs)。

地球物理模型的离散化是通过对模型的空间剖分实现的。地球物理模型的空间剖分方法目前主要有两种，即正交网格剖分和非正交网格剖分。所谓的正交网格就是在平面上

是矩形网格，而非正交网格在平面上是三角形网格和不规则四边形网格。对于地下介质进行非正交网格剖分可以差分考虑地下介质分布的几何形状，并且不受边界几何形态的限制。基于这一点，非正交网格数值模拟方法要优于正交网格的数值模拟方法。为了准确刻画地下介质物理性质的空间变化，网格剖分必须要足够精细，但是模型剖分得越细，空间网格点的数目就越多，这必然会占用大量的计算内存，加大计算量，降低效率，从而增加计

算的成本，同时引起误差的

积累。因此，模型离散化时必须考虑数值模拟的分辨率（或网

格大小）和计算成本。

2 . 3波动方程的离散化

模型空间的网格化必然带来波场的网格化。由于这种网格化把一个连续的地震波动问题转化成一个离散的地震波动问题，因此必然涉及到波场逼近，并且在空间网格化以后，尽可能以较小的逼近误差表示离散波场的空间微分。有限差分法通过有限差分算子将波动方程离散化，以差分代替微分，将微分方程问题转化为代数方程问题，然后求解相关的线性代数方程组以获得微分方程问题的数值解。差分算子是一个空间局部的算子，在空间域具有较高的分辨率，可以较好地适应剧烈变化的地下介质。但是在频率域中，有限差分算子的分辨率就非常低了。算法的稳定性和收敛性受空间采样率和时间采样率的影响，但算法的速度较快。

基于变分原理和网格插值的有限元法比较适合几何条

件和物理条件都较复杂的问题。但是算法复杂，计算速度慢，一般要求插值基函数是分段线性函数，不具有正交性，算子也是空间局部算子，空间分辨率高，但是频率域中分辨率却很低。

另外一种逼近空间微分的方法是伪谱法，它是利用傅立叶变换将波场函数表示为傅立叶级数的展开形式，将时间域的波动方程在频率域中求解。伪谱法对微分算子的逼近程度

可以达到尼奎斯特频率，并且收敛速度快。但由于傅立叶变换是