《2.5维各向异性介质中地震波场数值模拟》.pdfVIP

第28卷第4期 物探化探计算技术 2006年1t月 文章编号：1001—1749(2006)04—0289—05 2．5维各向异性介质中地震波场数值模拟 武 晔1，李小凡1，赵玉莲2…，李信富1 (1．中科院地质与地球物理研究所，北京100029； 2．中国石油勘探开发研究院西北分院，兰州730020； 3．成都理工大学，成都610059) 摘要：这里实现了在2．5维各向异性介质中地震波传播的数值模拟。首先从2．5维一阶速度 一应力弹性各向异性波动方程出发，得出了该方程的拟谱法数值解法；然后通过对Marmousi速 度模型进行数值模拟，表明了拟谱法的有效性；最后计算并分析了均匀和混合各向异性介质中波 的传播快照和理论模拟的地震记录，进一步认识了波在各向异性介质中的传播规律。 关键词：各向异性；弹性波；傅立叶变换；数值模拟 中图分类号：P631．4 文献标识码：A 甲的传播规律，计算f坷匀和琨合各同异性哆}．庾中 0前言 的波场快照和合成记录，并观测到了波的反射和透 射等现象。 在进行波动方程数值模拟中，通常采用二阶位 1 方法原理 移方程，也可采用一阶速度一应力方程，一阶方程 不用对弹性参数进行空间微分…，更能适应较复 1．1 各向异性介质中的弹性波动方程建立 杂的地质构造模型。作者在本文中采用2．5维一 由弹性动力学的基本方程和弹性本构方程，可 阶速度一应力波动方程进行数值模拟。 推导出TI介质中2．5维一阶双曲型波动方程的表 地震波动方程的数值模拟方法主要为有限差 达形式‘43为 分法(FD)、有限元法(FE)、边界元法(BE)、傅立 p以2i+i+^纵=鲁+鲁+正 叶拟谱法(PSM)、谱元法(SEM)等。有限差分法 的算法简单，但存在频散效应，要解决频散问题，需 pry=鲁+警+‘2i+i+^ 要加密网格，这样就导致计算量增加，效率下降。 有限元法是一种比较精确的数值模拟方法，但它需 纵=鲁+警+z纵2i+i+^ 要很高的计算机内存。拟谱法最先由Kosloff和 Baysal【21应用于地震波动问题。由于拟谱法可以 以=cll警+C13O以v__!, (1) 等效为高阶网格差分法，当其空间差分精度的阶数 达到无穷时的极限情况，故拟谱法的精度很高。作 以=c13鲁+C33O以v__i 者在本文中为了进一步验证拟谱法具有较高的计 算效率和精度，对Marmousi速度模型口1进行了数 以=c44警+c44鲁 值模拟。 为了进一步认识弹性波在复杂各向异性介质 吒=C66鲁x癖=c44誓 基金项目：国家自然科学基金项目；国家自然科学基金资助项目 收稿日期：2005—12—30 万方数据 290 物探化探计算技术 28卷 1．2 汲动方程的拟谱法数值模拟公式 弹性各向同性Marmousi模型已经成为估计速 J 傅立叶变换求空间导数的原理‘5 度精度的标准模型。Marmousi模型在地质构造上 氅掣