错格实数傅里叶变换微分算子及其在非均匀介质波动传播研究中的应用.pdf

第 46 卷 第 2 期 地 球 物 理 学 报 Vol . 46 , No . 2 　　　　　　　　　　 　　　　　　　　　　 　2003 年 3 月 CHINESE JOURNAL OF GEOPHYSICS 　 Mar . , 2 0 0 3 错格实数傅里叶变换微分算子及其在 非均匀介质波动传播研究中的应用 1 1 2 赵志新 　徐纪人 　堀内 茂木 1 中国地质科学院地质研究所 ,北京 　100037 2 日本科技厅防灾科学技术研究所 ,筑波 　305 - 0006 提出一种新的数值微分运算方法 ,即错格实数傅里叶变换微分法. 该方法的运算速度比错格复数傅里 摘 　要 　 叶变换数值微分解法快 0 . 33 倍 ; 因为该微分算法在整个微分运算过程中保留了奈奎斯特分量 ,使得它比普通分格 的实数傅里叶变换数值微分算法的精度高 ,稳定性好. 将该方法和 CagniardDe Hoop 解析法在求解半无限空间地震 波动的问题中进行比较 ,结果表明 ,新微分法的精度和解析方法的精度相同. 在非均匀介质中的地震波传播数值模 拟的结果表明 ,该方法是一种研究非均匀介质中地震波传播问题的有效的数值微分方法. 关键词 　错格分割法 　错格实数傅里叶变换微分 　数值模拟 　非均匀介质 　地震波传播 ( ) 文章编号 　0001 - 5733 2003 02 - 0234 - 07 　　中图分类号 　P631 　　收稿日期 　2002 - 01 - 15 ,2002 - 08 - 16收修定稿 STA GGERED GRID REAL VAL UE FFT DIFFERENTIATION OPERATOR AND ITS APPL ICATION ON WAVE PROPA GATION SIMULATION IN THE HETERO GENEOUS MEDIUM ZHAO ZHIXIN 1 　XU J IREN1 　SHIGEKI HORIUCHI2 1 Institute of Geology , Chinese Academy of Geological Science ,Beij ing 100037 , China 2 N ational Research Institute f or Earth Science and Disaster Prevention , Science and Technology Agency , Tsukuba 305 - 0006 ,J ap an Abstract 　A new fast differentiation operator , “staggered grid real value FFT differentiation operator ”, is proposed in this paper . The met