Seismic modeling plays an irreplaceable fundamental role in the exploration of oil and gas resources. Due to its straightforward implementation, small memory and computational costs, finite-difference method (FDM) has been widely used in computational seismology. However, FDM is known to have difficulties reproducing properly the properties of the free-surface condition. Moreover, the topography is poorly represented by regular Cartesian grids. We present a novel FD scheme which can solve accurately seismic wave propagation in anisotropic medium with irregular free surface. The method is based upon the Fully Staggered Grid combined with body-fitted grid to honor the topography and a mimetic solution for the free-surface condition. Numerical tests verify that our method can obtain results with very few points per wavelength when compared to other similar approaches. Furthermore, no explicit meshing is required in order to set up the models, which adds to the robustness and simplicity of the scheme.


Article metrics loading...

Loading full text...

Full text loading...


  1. Castillo, J., Hyman, J., Shashkov, M. and Steinberg, S.
    [2001] Fourth- and sixth-order conservative finite difference approximations of the divergence and gradient. Appl. Numer. Math., 37(1–2), 171–187.
    [Google Scholar]
  2. Hestholm, S. and Ruud, B.
    [1994] 2-D finite-difference elastic wave modeling including surface topography. Geophys. Prosp., 42, 371–390.
    [Google Scholar]
  3. de la Puente, J. et al.
    [2013] Elastic Mimetic Finite-differences in the Presence of Topography. 75th EAGE Conference &Exhibition, Extended Abstracts, 2013, Tu-01-09.
    [Google Scholar]
  4. Lebedev, V.
    [1964] Difference analogies of orthogonal decompositions of basic differential operators and some boundary value problems: I. Soviet Computational Mathematics and Mathematical Physics, 4, 449–465.
    [Google Scholar]
  5. Lisitsa, V. and Vishnevskiy, D.
    [2010] Lebedev scheme for the numerical simulation of wave propagation in 3D anisotropic elasticity. Geophysical Prospecting, 58, 619–635.
    [Google Scholar]
  6. Virieux, J.
    [1986] P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method. Geophysics, 51, 889–901.
    [Google Scholar]
  7. Zhang, W. and Chen, X.
    [2006] Traction image method for irregular free surface boundaries in finite difference seismic wave simulation. Geophys. J. Int., 167, 337–353.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error