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

错格实数傅里叶变换微分算子及其在 非均匀介质波动传播研究中的应用 赵志新1 徐纪人1 堀内 茂木2 1 中国地质科学院地质研究所 ,北京 100037 2 日本科技厅防灾科学技术研究所 ,筑波 305 - 0006 摘 要 提出一种新的数值微分运算方法 ,即错格实数傅里叶变换微分法 . 该方法的运算速度比错格复数傅里 叶变换数值微分解法快 0. 33 倍 ;因为该微分算法在整个微分运算过程中保留了奈奎斯特分量 ,使得它比普通分格 的实数傅里叶变换数值微分算法的精度高 ,稳定性好 . 将该方法和 Cagniard2De Hoop 解析法在求解半无限空间地震 波动的问题中进行比较 ,结果表明 ,新微分法的精度和解析方法的精度相同 . 在非均匀介质中的地震波传播数值模 拟的结果表明 ,该方法是一种研究非均匀介质中地震波传播问题的有效的数值微分方法 . 关键词 错格分割法 错格实数傅里叶变换微分 数值模拟 非均匀介质 地震波传播 文章编号 0001 - 5733 (2003) 02 - 0234 - 07 中图分类号 P631 收稿日期 2002 - 01 - 15 ,2002 - 08 - 16收修定稿 STA GGERED GRID REAL VAL UE FFT DIFFERENTIATIO N OPERATO R AND ITS APPL ICATIO N O N WAVE PROPA GATIO N SIMUL ATIO N IN THE HETERO GENEO US MEDIUM ZHAO ZHIXIN1 XU J IREN1 SHIGEKI HORIUCHI2 1 Institute of Geology , Chinese Academy of Geological Science , Beijing 100037 , China 2 National Research Institute f or Earth Science and Disaster Prevention , Science and Technology Agency , Tsukuba 305 - 0006 , J apan Abstract A new fast differentiation operator “, staggered grid real value FFT differentiation operator”, is proposed in this paper . The method is one third faster on average than the traditional differentiation method by using the complex FFT differentiation , and is very efficient in economy. Furthermore , the new method has higher accuracy and higher stability than the traditional grid real FFT method because the Nyquist component is remained throughout calculation procedure . The validity of staggered grid real value FFT differentiation operator is confirmed by comparing with the analytic Cagniard2De Hoop method in the half space SH problem. The new operation is applied to simulate seismic wave propagation in inhomogeneous medium ,and the results show that the staggered grid real value FFT differentiation operator for the p seudospectral method is valid for the seismic wave simulation in the heterogeneous medium. Staggered grid , Staggered grid real value FFT differentiation , Numerical