改进的高精度拟谱微分方法及其在三维非均匀介质地震传播研究中的应用 Time-staggered difference for pseudospectral differentiation openrator and its application in 3-D heterogeneous medium.pdfVIP

第51卷第6期 地 球 物 理 学 报 v01．51，No．6 2008年11月 CHINESEJOURNAl。OFGE()PHYSICS Nov．，2008 武哗，李小凡，顾观文等．改进的高精度拟谱微分方法及其在■维非均匀介质地震传播研究中的应用．地球物理学报．2008， 5l(6)：1868～1875 WuY。I。iXF，GuGW，et differencefor differentiationandits in3一 a1．Time-staggeredpseudospectral openratorapplication D Chinese)，2008，51(6)：1868--1875 heterogeneousmedium．ChineseJ．Geophys．(in 改进的高精度拟谱微分方法及其在三维非均匀 介质地震传播研究中的应用 武 哗1’2，李小凡2，顾观文3，张美根2，龙桂华2，赵晓燕 1防灾科技学院，北京燕郊06520l 2中圃科学院地质与地球物理研究所．北京 100029 3中国地质科学院物化探研究所．廊坊065000 摘 要 改进了．用于模拟地震波场的傅里叶拟谱微分方法，改进后的方法精度足常规拟谱方法的4倍，称为改进 的傅咀叶拟谱方法．在较高数值精度的一阶应力一速度弹性波动方程的基础上，采用该方法和常规拟谱方法对 Marmousi模型进行数值求解，结果表明，该方法的数值频散效应明显比常规拟谱方法弱．将该方法与有限元方法 在各向异性介质中进行模拟比较，发现该方法的精度接近有限元方法，数值频散效应比有限元方法明显减小，而且 可在较大空l日j网格间距下进行汁算，从而提高计算效率．在3一D非均匀介质中的地震波传播数值模拟结果表明，该 方法足一种研究复杂非均匀介质中地震波传播问题的高效方法． 关键词 弹性波．傅咀叶变换，数值模拟，拟谱，三维 08 文章编号0001—5733(2008)06—1868中图分类号P631 收稿日期2007—11-06，2008-07—21收修定稿 differencefor differentiation Time-staggered pseudospectral openrator andits in3一D medium applicationheterogeneous WUYel．。，LIXiao—Fan2，GUGuan—Wen3，ZHANGMei—Gen2， LONGGui—Hua2．ZHA0Xiao—Yanl 1 InstituteDisasterPreventionScienceand 065201．China of Technology．BeijingYanjiao 2 Institute and 100029．China ofGeologyGeophysics，ChineseAcademyofSciences，Beijing 3 Institute andGeochemical 065000。China ofGeophysical Exploration．Langfang AbstractAnewfastdifferentiation difference