改进的高精度拟谱微分方法及其在三维非均匀介质地震传播研究中的应用  

Time-staggered difference for pseudospectral differentiation openrator and its application in 3-D heterogeneous medium

作　　者：武晔[1,2] 李小凡[2] 顾观文[3] 张美根[2] 龙桂华[2] 赵晓燕[1] 

机构地区：[1]防灾科技学院,北京燕郊065201 [2]中国科学院地质与地球物理研究所,北京100029 [3]中国地质科学院物化探研究所,廊坊065000

出　　处：《地球物理学报》2008年第6期1868-1875,共8页Chinese Journal of Geophysics

基　　金：国家自然科学基金重点项目(40437018);国家自然科学基金青年基金项目(40304004);国家杰出青年科学基金项目(40025412);防灾科技学院防灾减灾青年科技基金项目联合资助

摘　　要：改进了用于模拟地震波场的傅里叶拟谱微分方法,改进后的方法精度是常规拟谱方法的4倍,称为改进的傅里叶拟谱方法.在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用该方法和常规拟谱方法对Marmousi模型进行数值求解,结果表明,该方法的数值频散效应明显比常规拟谱方法弱.将该方法与有限元方法在各向异性介质中进行模拟比较,发现该方法的精度接近有限元方法,数值频散效应比有限元方法明显减小,而且可在较大空间网格间距下进行计算,从而提高计算效率.在3-D非均匀介质中的地震波传播数值模拟结果表明,该方法是一种研究复杂非均匀介质中地震波传播问题的高效方法.A new fast differentiation operator, ＂ pseudospectral time-staggered difference operator＂, is proposed in this paper. The accuracy of this method is four time of the normal pseudospectral method. It is confirmed that the numerical dispersion of this method is less than that of normal pseudospectral method by comparing the calculation result of Marmousi model. The new operation is applied to simulate seismic wave propagation in 3-D heterogeneous medium, and the results show that the time-staggered difference operator for the pseudospectral method is valid for the seismic wave simulation in the complex heterogeneous medium.

关 键 词：弹性波 傅里叶变换 数值模拟 拟谱 三维 

分 类 号：P631[天文地球—地质矿产勘探]