TY - JOUR
T1 - A highly accurate finite-difference method with minimum dispersion error for solving the Helmholtz equation
AU - Wu, Zedong
AU - Alkhalifah, Tariq Ali
N1 - KAUST Repository Item: Exported on 2021-02-19
Acknowledgements: We thank KAUST for its support and the SWAG group for the collaborative environment. We also thank BP for providing the benchmark dataset. The research reported in this publication is supported by funding from King Abdullah University of Science and Technology (KAUST). For computer time, this research used the resources of the Supercomputing Laboratory at King Abdullah University of Science and Technology (KAUST) in Thuwal, Saudi Arabia. We also thank the associate editor Eli Turkel and another anonymous reviewer for their fruitful suggestions and comments.
PY - 2018/4/5
Y1 - 2018/4/5
N2 - Numerical simulation of the acoustic wave equation in either isotropic or anisotropic media is crucial to seismic modeling, imaging and inversion. Actually, it represents the core computation cost of these highly advanced seismic processing methods. However, the conventional finite-difference method suffers from severe numerical dispersion errors and S-wave artifacts when solving the acoustic wave equation for anisotropic media. We propose a method to obtain the finite-difference coefficients by comparing its numerical dispersion with the exact form. We find the optimal finite difference coefficients that share the dispersion characteristics of the exact equation with minimal dispersion error. The method is extended to solve the acoustic wave equation in transversely isotropic (TI) media without S-wave artifacts. Numerical examples show that the method is is highly accurate and efficient.
AB - Numerical simulation of the acoustic wave equation in either isotropic or anisotropic media is crucial to seismic modeling, imaging and inversion. Actually, it represents the core computation cost of these highly advanced seismic processing methods. However, the conventional finite-difference method suffers from severe numerical dispersion errors and S-wave artifacts when solving the acoustic wave equation for anisotropic media. We propose a method to obtain the finite-difference coefficients by comparing its numerical dispersion with the exact form. We find the optimal finite difference coefficients that share the dispersion characteristics of the exact equation with minimal dispersion error. The method is extended to solve the acoustic wave equation in transversely isotropic (TI) media without S-wave artifacts. Numerical examples show that the method is is highly accurate and efficient.
UR - http://hdl.handle.net/10754/627484
UR - http://www.sciencedirect.com/science/article/pii/S0021999118302134
UR - http://www.scopus.com/inward/record.url?scp=85045045568&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2018.03.046
DO - 10.1016/j.jcp.2018.03.046
M3 - Article
AN - SCOPUS:85045045568
SN - 0021-9991
VL - 365
SP - 350
EP - 361
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -