基于拟谱区域分解法地震波数值模拟.pdfVIP

基于拟谱区域分解法地震波数值模拟 武晔 李小凡 (中科院地质与地球物理研究所 北京 100029) 地震波场数值模拟是对地下介质结构和相应物理参数已知的情况下。模拟研究地震波在地下各种介质中的 传播规律。并计算在地面和地下各观测点所对应的地震记录的一种地震模拟方法。波场数值模拟技术现已发展成 为勘探地球物理学的一个重要领域，形成一套完整的理论体系和一系列的方法技术，主要体现在以下几个方面： 模拟方法种类多样主要包括积分法、射线追踪法、波动方程法：模拟维数由二维走向三维：模拟介质由简单走向 复杂比如由水平层状介质发展到各向异性、双相、多相介质。地震波场数值模拟方法已在地震勘探和天然地震 领域中得到广泛的应用，而且在地震灾害预测、地震区带划分以及地壳构造和地球内部结构研究中，也得到广泛 的应用。 地震波动方程的数值模拟方法主要有：有限差分法(FD)、有限元法(FE)、边界元法(BE)、傅立时拟谱法、 谱元法(sEM)等。有限差分法的算法简单快速，但存在频散效应，要解决频散问题，需要加密数值计算的网格，这 样导致计算量增加，效率下降。有限元法是一种比较精确的数值模拟方法，但是它对计算机内存要求很高，计算量 用予地震波动问题。该算法通过对空间坐标的偏微分实施快速傅立叶变换，只在时间域作差分计算。由于拟谱法 可以等效为高阶网格差分法当其空间差分精度的阶数达到无穷时的极限情况，故拟谱法的精度是无限阶的。 区域分解算法(DomainDecompositionMethod)是一类求解偏微分方程数值解的新技术，由于近年来并行计 算机和并行算法的发展，该方法越来越受到人们的重视，并且已经成为计算数学的热门领域。区域分解算法是把 计算区域分解为若干个子域，子域的形状尽可能规则，这样原阀题就转化为子域上的求解。在各予域内采用最有 效的方法求解，能充分发挥各种数值模拟方法的优点，大大提高求解的有效性。区域分解算法的优越性在于： 1． 子区域如果形状规则，其上可以使用各种快速算法。 2． 可以在不同子区域上选用不同的数学模型，以便整体模型更适合具体实际情况。 3． 可以局部采用一致网格，不需要整体采用一致网格。各子区域上可以采用不同的离散方法计算，故区域 分解算法具有很大的灵活性。 4． 把大问题分解为若干小问题，减小计算规模。 5． 算法高度并行，而且对于对称区域问题有更简单的区域分解算法。 本文所涉及的工作将区域分解法应用于地震波场数值模拟计算，进一步研究地震波在各向异性介质中的传 播特征和传播规律．由于傅立叶变换微分算子是一个全局算予，如果把该算子直接应用到整个区域，将平滑掉大 部分局部信息。为了克服这一缺点，本研究将拟傅氏谱算子局域化，以突出局部信息。此外，本工作将把该方法 应用到声波方程和各向异性介质波动方程的数值模拟计算，以证明该方法的有效性。本工作还将比较拟谱区域分 解法和拟谱法在计算速度、计算精度方面的差异，以论证本方法是速度快，精度高，所占内存少的一种可行有效的 数值模拟方法。由于模拟中在处理非均匀介质目题时需要大量计算资源，因此与本模拟有关的数值计算，在浪 潮天梭10000大型并行计算集群服务器系统上进行。 参考文献 a1．ELASTICWAVECALCULATlONSBYTHEFOURIER [1]DanKosloff,．MosheReeshefet METHOD， BSSA，1984，74：875～891 [2]吕涛，石济民，林振宝，区域分解算法一偏微分方程数值解新技术，科学出版社，1991 [3]杨绍国等，区域分解法地震正演模拟方法研究，‘地球物理学报》，1996(39)：690～697． ·528·