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

Zedong Wu; Tariq Alkhalifah

doi:10.1016/j.jcp.2018.03.046

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

Zedong Wu^*, Tariq Alkhalifah

^*Corresponding author for this work

Physical Sciences and Engineering

Research output: Contribution to journal › Article › peer-review

31 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 highly accurate and efficient.

Original language	English (US)
Pages (from-to)	350-361
Number of pages	12
Journal	Journal of Computational Physics
Volume	365
DOIs	https://doi.org/10.1016/j.jcp.2018.03.046
State	Published - Jul 15 2018

Bibliographical note

Funding Information:
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.

Publisher Copyright:
© 2018 Elsevier Inc.

Keywords

Acoustic wave equation
Anisotropic
Dispersion error
Finite difference

ASJC Scopus subject areas

Computational Mathematics
Physics and Astronomy(all)
Applied Mathematics
Numerical Analysis
Computer Science Applications
Modeling and Simulation
Physics and Astronomy (miscellaneous)

Access to Document

10.1016/j.jcp.2018.03.046

Cite this

@article{9c37953ed5a44b1b883c05648e660cfb,

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

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 highly accurate and efficient.",

keywords = "Acoustic wave equation, Anisotropic, Dispersion error, Finite difference",

author = "Zedong Wu and Tariq Alkhalifah",

note = "Funding Information: 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. Publisher Copyright: {\textcopyright} 2018 Elsevier Inc.",

year = "2018",

month = jul,

day = "15",

doi = "10.1016/j.jcp.2018.03.046",

language = "English (US)",

volume = "365",

pages = "350--361",

journal = "Journal of Computational Physics",

issn = "0021-9991",

publisher = "Elsevier",

}

TY - JOUR

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

AU - Wu, Zedong

AU - Alkhalifah, Tariq

N1 - Funding Information: 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. Publisher Copyright: © 2018 Elsevier Inc.

PY - 2018/7/15

Y1 - 2018/7/15

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 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 highly accurate and efficient.

KW - Acoustic wave equation

KW - Anisotropic

KW - Dispersion error

KW - Finite difference

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 -

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

Abstract

Bibliographical note

Keywords

ASJC Scopus subject areas

Access to Document

Other files and links

Fingerprint

Cite this