2-D elastic wave modeling with frequency-space 25-point finite-difference operators

作     者：Liao Jianping Wang Huazhong Ma Zaitian Liao Jianping;Wang Huazhong;Ma Zaitian

作者机构：School of Ocean and Earth Science Tongji University Shanghai China 

出 版 物：《Applied Geophysics》 (应用地球物理（英文版）)

年 卷 期：2009年第6卷第3期

页      面：259-266,300页

核心收录：

学科分类：0711[理学-系统科学] 081801[工学-矿产普查与勘探] 081802[工学-地球探测与信息技术] 07[理学] 08[工学] 0818[工学-地质资源与地质工程] 070105[理学-运筹学与控制论] 081101[工学-控制理论与控制工程] 071101[理学-系统理论] 0811[工学-控制科学与工程] 0701[理学-数学] 

基　　金：supported by the 863 Program (Grant no.2006AA09Z323) the 973 Program (Grant No.2006CB202402) 

主　　题：compressed storage frequency-space domain twenty-five point finite-difference optimal coefficients PML 

摘      要：Numerical simulation in the frequency-space domain has inherent advantages, such as： it is possible to simulate wave propagation from multiple sources simultaneously; there are no cumulative errors; only the interesting frequencies can be selected; and it is more suitable for wave propagation in viscoelastic media. The only obstacle to using the method is the requirement of huge computer storage. We extend the compressed format for storing the coefficient matrix. It can reduce the required computer storage dramatically. We get the optimal coefficients by least-squares method to suppress the numerical dispersion and adopt the perfectly matched layer （PML） boundary conditions to eliminate the artificial boundary reflections. Using larger grid intervals decreases computer storage requirements and provides high computational efficiency. Numerical experiments demonstrate that these means are economic and effective, providing a good basis for elastic wave imaging and inversion.