地震数值模拟.doc

地震数值模拟 1.1 地震数值模拟 地震数值模拟是地震勘探和地震学的重要基础。所谓地震数值模拟就是在假定地下介质结构模型和相应物理参数已知的情况下，模拟研究地震波在地下各种介质中的传播规律，并计算在地面或地下各观测点所应观测到的数值地震记录的一种地震模拟方法。这种地震数值模拟方法已在地震勘探和天然地震领域中得到广泛的应用。它不但在石油、天然气、煤、金属和非金属等矿产资源及工程和环境地球物理中得到普遍的应用，而且在地震灾害预测、地震区带划分以及地壳构造和地球内部结构研究中，也得到相当广泛的应用。地震数值模拟在地震勘探和地震学各工作阶段中都有重要的作用。在地震数据采集设计中，地震数值模拟可用于野外地震观测系统的设计和评估，并进行地震观测系统的优化。在地震数据处理中，地震数值模拟可以检验各种反演方法的正确性。在地震数据处理结果的解释中，地震数值模拟又可以对地震解释结果的正确性进行检验。 地震数值模拟与地震物理模拟同属于地震正演模拟，即在已知地下介质结构和物理参数模型情况下，预测在地面或井中各观测点所应观测到的地震记录。而地震勘探的目的则是根据地面或井中各观测点所观测的地震记录来刻画地下介质结构模型，并描述其状态或岩性。这是一个反演过程，但是，这个反演过程是建立在地震正演模型的基础上的。因此地震数值模拟不仅可以进行地震正演模拟研究，同时也是地震反演的基础。 地震数值模拟的发展非常迅速，现在已经研制出了各种各样的地震数值模拟方法，并均在地震勘探和地震学中得到广泛而有效的应用，这些地震数值模拟方法可以归纳为地震波方程数值解法、积分方程法和射线追踪法三大类（Carcione,2002）。 1.2 地震数值模拟的理论基础 地震数值模拟是在已知地下介质结构情况下，研究地震波在地下各种介质中传播规律的一种地震模拟方法，其理论基础就是表征地震波在地下各种介质中传播的地震波传播理论。 上述三类地震数值模拟方法相应的地震波传播理论的数学物理表达方式也不尽相同。地震波方程数值解法是建立在以弹性或黏弹性理论和Newton力学为基础的双曲型偏微分方程—地震波传播方程的理论基础上的。由于地下介质性质不同，其相应的地震波传播方程也不同。如声学介质中的声波波动方程；弹性介质中的弹性波波动方程；黏弹性介质中的黏弹性波波动方程；孔隙弹性介质（双相或多相介质）中的双相（或多相）介质弹性波方程；各向异性介质中的各向异性弹性波波动方程等。积分方程法是建立在以Huygens原理为基础的波叠加原理基础上的，其数学表达形式为波动方程的Green函数域积分方程式和边界积分方程式。射线追踪法是建立在以射线理论为基础的波动方程高频近似理论基础上的，其数学表达形式为程函方程和传输方程。 由于地震波传播方程在三维复杂介质中地震波传播的广泛适应性及地震波方程数值解法在地震波数值模拟中应用的广泛性和有效性，下面，我们将着重研究地震波方程数值解法，同时，对积分方程和射线追踪法也进行适当的讨论。 2 地震数值模拟方法及其特点 2.1 地震数值模拟方法 地震数值模拟是以地震波传播理论为基础的。描述地震波在各种介质中传播的波动方程属于双曲型变系数偏微分方程。地震波波动方程的定解问题（或正演方程）包括：微分算子；微分算子的系数（可变或常系数）；震源项；边界条件；初始条件等。 地震波正演过程是由“因”推“果”的确定过程，求解微分方程可以算出系统中表示状态的参量随时间的演化。地震波的正演过程数学上可描写为： 其中， —合成地震数据向量； —正演算子； —模型向量。 的精度受 的离散化精度和正演算子 计算精度的影响。 地震波传播方程正演问题的研究内容: （1） 地震波数值模拟原理 （2） 地震波数值模拟算法 （3） 计算程序的质量。它受三个主要因素影响，即计算精度，计算效率和计算的稳定性。 另外，现代计算机的飞速发展极大地促进了地震数值模拟研究的进展。 现代地震波数值模拟是利用先进的计算机技术将野外一定规模的复杂地质构造、地质体和复杂岩性通过可视化建模技术建立二维、三维地质-地球物理模型，并用先进的数值模拟方法对野外地震勘探方法进行模拟（如图1所示）。 2.2 地质—地球物理模型离散化 不同的研究目的，构成地质—地球物理模型的地球物理量也不同。对于均匀各向同性介质声波方程，地质—地球物理模型 可表示为 而对均匀各向同性介质弹性波方程而言，其地质—地球物理模型则表示为： 即模型 离散化为包含有限个模型物理参数的向量 。 地质—地球物理模型的离散化是通过对模型的空间剖分来实现的。目前地质—地球物理模型的空间剖分方法可分为两种，即正交网格剖分和非正交网格剖分。正交网格在平面上是矩形网格，而非正交网格剖分在平面上是三角形网格和不规则四边形网格。对于地下介质进行非正交网格剖