A highly accurate finite-difference method with minimum dispersion error for solving the Helmholtz equation

Research output: Contribution to journalArticlepeer-review

18 Scopus citations

Abstract

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.
Original languageEnglish (US)
Pages (from-to)350-361
Number of pages12
JournalJournal of Computational Physics
Volume365
DOIs
StatePublished - Apr 5 2018

Fingerprint

Dive into the research topics of 'A highly accurate finite-difference method with minimum dispersion error for solving the Helmholtz equation'. Together they form a unique fingerprint.

Cite this