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

《地震波数值模拟方法研究综述new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、万方数据第27卷第2期2007年5月防灾减灾工程学报JournalofDisasterPreventionandMitigationEngineeringV01．27No．2May2007地震波数值模拟方法研究综述+李信富，李小凡，张美根(中国科学院地质与地球物理研究所岩石圈演化国家重点实验室，北京100029)摘要：地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段，而且是研究地球深部精细构造的有效工具。地震波场数值模拟的主要方法包括3大类，即几何射线法、积分方程法和波动方程法。首先介绍了地震波数值模拟方法的理论基础和内容及特点，然后概括介绍

2、了20世纪90年代以来几何射线法和各种波动方程数值模拟方法的新进展。因为波动方程数值模拟包含了丰富的波动信息，可为研究地震波的传播机理和复杂地层的解释提供更多的佐证，所以重点介绍波动方程数值模拟方法。最后，对各种数值模拟方法进行了评述，并给出了地震波传播理论和数值模拟方法发展的展望。关键词：地震波；数值模拟；射线追踪；积分方程；波动方程中图分类号：TU315．3“文献标识码：A文章编号：1672．2132(2007)02一0241一080引言在地学领域，对于许多地球物理问题，人们已经得到了它们应遵循的基本方程(常微分方程或偏微分方程)和相应的定解条件，但能用解析方法

3、求得精确解的只是少数方程性质比较简单，且几何形状相当规则的问题。对于大多数问题，由于方程的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析解。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但这种方法只是在有限的情况下是可行的，过多的简化可能导致很大的误差甚至错误的解答。因此人们多年来寻找和发展了另一种求解方法——数值模拟方法。地震数值模拟(SeismicNumericalModeling)是地震勘探和地震学的基础，同时也是地震反演的基础。所谓地震数值模拟，就是在假定地下介质结构模型和

4、相应的物理参数已知的情况下，模拟研究地震波在地下各种介质中的传播规律，并计算在地面或地下各观测点所观测到的数值地震记录的一种地震模拟方法。地震波场数值模拟是研究复杂地区地震资料采集、处理和解释的有效辅助手段，这种地震数值模拟方法已经在地震勘探和天然地震领域中得到广泛应用。地震数值模拟的发展非常迅速，现在已经有各种各样的地震数值模拟方法在地震勘探和地震学中得到广泛而有效的应用。这些地震波场数值模拟方法可以归纳为三大类，即几何射线法、积分方程法和波动方程法。波动方程数值模拟方法实质上是求解地震波动方程，因此模拟的地震波场包含了地震波传播的所有信息，但其计算速度相对于几何

5、射线法要慢。几何射线法也就是射线追踪法，属于几何地震学方法，由于它将地震波波动理论简化为射线理论，主要考虑的是地震波传播的运动学特征，缺少地震波的动力学信息，因此该方法计算速度快。因为波动方程模拟包含了丰富的波动信息，为研究地震波的传播机理和复杂地层的解释提供了更多的佐证，所以波动方程数值模拟方法一直在地震模拟中占有重要地位。1地震波数值模拟的理论基础地震波数值模拟是在已知地下介质结构的情况下，研究地震波在地下各种介质中传播规律的一种地震模拟方法，其理论基础就是表征地震波在地下各种介质中传播的地震波传播理论。上述三类地震·收稿日期：2006一08一08；修回日期：2

6、006一09—20基金项目：国家自然科学基金重点项目(40437018)；国家自然科学基金项目(40304004)作者简介：李信富(1976一)，男，博士研究生。主要从事地震波数值模拟和非线形插值方法的研究。Email：xinfulee@mail．iggcas．ac．cn万方数据242防灾减灾工程学报第27卷波数值模拟方法相应的地震波传播理论的数学物理表达方式不尽相同。射线追踪法是建立在以射线理论为基础的波动方程高频近似理论基础上的，其数学表达形式为程函方程和传输方程。积分方程法是建立在以惠更斯原理为基础的波叠加原理基础上的，其数学表达形式为波动方程的格林函数域积分

7、方程表达式和边界积分方程表达式。波动方程数值解法是建立在以弹性或粘弹性理论和牛顿力学为基础的双曲型偏微分方程一波动方程的理论基础上的。由于地下介质性质不同，其相应的地震波传播方程也不同。由于地震波波动方程在复杂介质中地震波传播研究的广泛适应性及地震波方程数值解法在地震波数值模拟中应用的广泛性和有效性，本文将重点研究地震波波动方程数值解法，同时对几何射线法和积分方程法也作适当的讨论。2地震波数值模拟的内容及特点2．1地震波数值模拟的内容地震波数值模拟是以地震波传播理论为基础的。描述地震波在各种介质中传播的波动方程是一个变系数的偏微分方程。地震波波动方程的定解问题(