用谱元法计算地震波场.pdfVIP

用谱元法计算地震波场 徐果明 沈玉松朱良保 翟永波 中国科学技术大学 地球和空间学院 安徽合肥230026 地球物理研究中的各种问题，可通过不同途径归结为不同形式的数学模型。在涉及地震波的正演计算领域， 则往往归结为波动方程的求解。对于简单的介质模型和地震源有解析解，而对于比较复杂的问题则要靠数值解法。 问题的求解或是表现为偏微分方程的边值问题，或是表现为区域上的变分问题，或是归结为边界上的积分方程。 这些不同的数学形式在理论上是等价的，但在实践中并不等效，它们分别导致有限差分法、有限元方法和边界元 方法、谱方法等不同的数值计算方法。在弹性动力学问题中，有限差分法是一种解决动力学问题的常用方法，它 的基本思想是基于将运动方程在时间和空间上均进行了均匀网格上的离散化。而有限元法是将复杂的弹性体看成 是由无数个质点组成的连续体i谱方法实际上是一种变分法，基本思想是将问题的解按某些整体光滑的函数系展 开，取前|7、r项部分和作为逼近解。逼近函数的选择是谱方法与有限元和有限差分相区别的特征之一。逼近函数的 versions中也是有区别的。配置方法是 选择在三种最常用的谱系(spectralschemes)：Galerkin，collocation和tau 谱方法中常用的。逼近函数的恰当选择和配置点的分配对解的准确性格外重要。谱方法最常用的函数系Chebyshev 多项式和Legendre多项式等。由函数逼近论的结果可知，函数的光滑性越好，则其谱逼近的精度越高．对于光滑 问题，谱方法解可获得类似的高精度，被称为具有“无穷阶”的收敛性。谱方法在例如流体力学问题的计算中， 已经取得很大成功，日益受到重视。 利用弹性波方程来解决不均匀介质中的地震波传播问题是地球物理学中的重要课题。地球物理学家们试图寻 找一种合适的方法来准确地解弹性波方程，同时还要考虑到实际应用中遇到的各种具体情况，如复杂的地形、2D 或3D的物质分界面、不均匀介质等等。前面提到的传统方法自然会被想到用来解决上述问题。但是，这些方法 或多或少都存在一些问题，使他们的使用有一定的局限性。 传播问题。其方法是基于Legendre多项式和Gauss·Lobatto 或Legendre)和配置点(涉及到数值积分)的选择不会显著的影响误差，但还是强烈关系到产生的代数多项式 形式的好坏和稀疏性，对平行迭代程序的效率影响也很大。在SEM方法中，模型被分隔为许多元，在每一个元中 函数在GLL积分点上取值。因此网格点取在元的内部、边和角点上。因此要区分阿格点。定义～个元的网格点称 为局部网格点。所有的网格点，其中有许多网格点被多于一个的元共享，称为全局网格点。因此需要确定局部网 格点与全局网格点之间的映射关系。在系统时间要前进一步时，要得到多个元共享的全局网格点上的函数值时， 应当将与此全局点相关元的贡献相加。这一步由于涉及到不同元信息之间的交换，对于平行计算可能会涉及不同 CPU之间信息的交换，因此是比较耗时的。谱元法对于网格点密度的要求要低一些，比如传统的二阶q，心差分中， 每个波长距离内至少要放置15个网格点。而谱元法只要5个网格点。它对复杂几何体系有很好的适应性。同时计 算精度高，非常有利于平行计算。 我们设计T--维空间中的谱元法的弹性波传播的计算程序。包括存在内界蕊和衰减的和速度随深度连续变化 的介质中激发的地震波，得到了地震波场波场和地震剖面图，包括纵波、横波和面波，直达波，反射波和转换波 等。对计算结果的仔细分析，得到一些感兴趣的结果。(国家自然科学基金资助项