傅里叶有限差分法三维波动方程正演模拟.pdf

第 50 卷 第 6 期 地 　球 　物 　理 　学 　报 Vol. 50 , No. 6 2007 年 11 月 CHINESE 　JOURNAL 　OF 　GEOPHYSICS Nov. , 2007 张金海 ,王卫民 ,赵连锋等. 傅里叶有限差分法三维波动方程正演模拟. 地球物理学报 ,2007 ,50 (6) :1854～1862 Zhang J H ,Wang W M ,Zhao L F ,et al. Modeling 3D scalar waves using the Fourier finitedifference method. Chinese J . Geophys . (in Chinese) , 2007 , 50 (6) :1854～1862 傅里叶有限差分法三维波动方程正演模拟 张金海 ,王卫民 ,赵连锋 ,姚振兴 中国科学院地质与地球物理研究所,北京　100029 ( ) 摘 　要 　傅里叶有限差分 FFD 法兼有相位屏法和隐式有限差分法二者的优势 ,能够处理复杂地质构造中的波传 播问题 ,但在三维情形下 ,算子的双向分裂会引起明显的方位各向异性误差. 本文用 Fourier 变换计算双向分裂过程 中的高阶交叉项 ,消除了方位各向异性误差. 该方法充分利用了 FFD 法在双域实现的算法结构 , 明显减少了由于引 入误差校正所带来的计算量. 将该方法应用于修改后的三维 French 模型的地震正演问题 ,并将得到的叠后记录、单 炮记录同全波有限差分法的模拟结果进行对比 ,结果证实了该方法对一次反射波具有较高的模拟精度 ,在内存需 求和计算效率方面则具有更大的优势. 关键词 　三维正演模拟 , 傅里叶有限差分 , 双向分裂 , 波动方程 ( ) 文章编号　0001 - 5733 2007 06 - 1854 - 09 中图分类号　P631 收稿日期　2007 - 04 - 04 , 2007 - 05 - 15 收修定稿 Modeling 3D scalar waves using the Fourier finitedifference method ZHANGJinHai ,WANG WeiMin ,ZHAO LianFeng ,YAO ZhenXing Institute of Geology and Geophysics , Chinese Academy of Sciences ,Beijing 100029 , China Abstract 　Fourier finitedifference ( FFD) operator has both of the advantages of Fourier and finitedifference methods, it can handle the wave propagations in complex media with large velocity contrasts and wide propagation angles. In 3D case ,however ,twoway splitting may cause artificial azimuthal anisotropy. In this paper ,we handle the remnants of the conventional twoway splitting technique by Fourier tran