Some Aspects of Numerical Simulation of Seismic Wave Propagation
Journal of Geology & Geophysics

Journal of Geology & Geophysics
Open Access

ISSN: 2381-8719

+44 7868 792050

Editorial - (2013) Volume 2, Issue 2

Some Aspects of Numerical Simulation of Seismic Wave Propagation

Jingyi Chen*
Department of Geosciences, The University of Tulsa, Tulsa, OK 74104, USA
*Corresponding Author: Jingyi Chen, Department of Geosciences, The University of Tulsa, Tulsa, OK 74104, USA Email:

Over the last decades, numerical simulation of seismic wave propagation has become an important tool to better understand the characterization of wave propagation and rock properties. It also plays a significant role in seismic data processing and interpretation [1]. Here, several important aspects of seismic numerical modeling of wave propagation will be mentioned, such as numerical algorithms for solving partial differential wave equations, absorbing boundary condition and high performance computing techniques.

Numerical modeling of wave propagation has been used in the most types of wave equations in time and frequency domains [2,3], for example, acoustic, elastic, poroelastic and viscoelastic equations. Several numerical algorithms such as finite difference, finite element, spectral element and pseudo-spectral methods have been used to solve those wave equations [2,4-6]. These numerical algorithms have different numerical precisions due to grid dispersion and numerical approximation.

To model seismic wave propagation in an unbounded medium, many absorbing boundary conditions have been developed for suppressing the artificial reflections at the edges of the truncated medium, for instance, paraxial conditions [7], continued fraction conditions [8] and sponge absorbing conditions [9]. Currently, widely applied boundary conditions are the perfectly matched layer (PML) [10,11] and modified PMLs such as Convolutional perfectly matched layer (CPML) and nearly perfectly matched layer (NPML) [2,12]. These PMLs, originally proposed to the first-order electromagnetic (EM) wave equations, were applied to first-order seismic wave equations with velocity-stress forms. Soon after, some PMLs are developed to second-order displacement wave equations in acoustic and elastic media [13,14].

Due to needs of large model space and computational memory, high performance techniques have received more and more attention on the development of seismic numerical modeling. After the early 90s, parallel computing techniques were fast developed based on conventional message passing interface (MPI) and open multiprocessing (OpenMP) [15]. Currently, modern graphics processing unit (GPU) and GPU cluster have the tendency toward replace the conventional parallel computing techniques [16].


  1. Virieux J and Operto S (2009) An overview of full-waveform inversion in exploration geophysics. Geophysics 74: WCC127-WCC152.
  2. Chen J and Zhao J (2011) Application of the nearly perfectly matched layer to seismic wave propagation modeling in elastic anisotropic media. Bulletin of the Seismological Society of America 101: 2866-2871.
  3. Dupuy B, Barros L D, Garambois S and Virieux J (2011) Wave propagation in heterogeneous porous media formulated in the frequency-space domain using a discountinuous Galerkin method. Geophysics 76: N13-N28.
  4. Takenaka H, Wang Y and Furumura T (1999) An efficient approach of the pseudospectral method for modeling of geometrically symmetric seismic wavefield. Earth Planets Space 51: 73-79.
  5. Komatitsch D, Barnes C and Tromp J (2000) Simulation of anisotropic wave propagation based upon a spectral element method. Geophysics 65: 1251-1260.
  6. De Basabe J D and Sen M K (2007) Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations. Geophysics 72: T81-T95.
  7. Clayton R and Engquist B (1977) Absorbing boundary conditions for acoustic and elastic wave equations. Bulletin of the Seismological Society of America 67: 1529-1540.
  8. Guddati M N and Lim K W (2006) Continued fraction absorbing boundary conditions for convex polygonal domains. International Journal for Numerical Methods in Engineering 66: 949-977.
  9. Cerjan C, Kosloff D, Kosloff R and Reshef M (1985) A nonreflecting boundary condition for discrete acoustic and elastic wave equations. Geophysics 50: 705-708.
  10. Berenger J (1994) A perfectly matched layer for the absorption of electromagnetic waves. Journal of Computational Physics 114: 185-200.
  11. Collino F and Tsogka C (2001) Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics 66: 294-307.
  12. Martin R, Komatitsch D and Ezziani A (2008) An unsplit convolutional perfectly matched layer improved at grazing incidence for seismic wave propagation in poroelastic media. Geophysics 73: T51-T61.
  13. Li Y and Matar O B (2010) Convolutional perfectly matched layer for elastic second-order wave equation. Journal of the Acoustical Society of America 127: 1318-1327.
  14. Metin G, Chen J and Ozsoy C (2013) Application of nearly perfectly matched layer with second-order acoustic equations in seismic numerical modeling. Journal of Geology & Geosciences.
  15. Bao H, Bielak J, Ghattas O, Kallivokas L F, O’Hallaron D R, et al. (1998) Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers. Computer Methods in Applied Mechanics and Engineering 152: 85-102.
  16. Komatitsch D, Erlebacher G, Goddeke D and Michea D (2010) High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster. Journal of Computational Physics 229: 7692-7714.
Citation: Chen J (2013) Some Aspects of Numerical Simulation of Seismic Wave Propagation. J Geol Geosci 2:e109.

Copyright: © 2013 Chen J. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.