Classical finite-difference method has some problems when dealing with irregular topography, including artifacts caused by staircase approximation to irregular free surface. Body-fitted grid is usually termed as “boundary conforming grid”, and it can solve the problem. On the other hand, Perfectly Matched Layer (PML) absorbing boundary is certified an efficient method to suppress spurious edge reflections. However, classical PML has the defects of its inherent instability and the effects of grazing incidence and some improved algorithms are developed. This paper models seismic wave in elastic media with irregular free surface by using the Lebedev Grid which is usually used in anisotropy media. In simulation process, the traction image method is used to implement the free surface conditions. And for other boundaries the improved perfectly matched layer technologies are chosen to absorb waves. Numerical tests on synthetic data demonstrate the validity of the proposed method. Finally we compare absorb effects and stability of some perfectly matched layer technologies.


Article metrics loading...

Loading full text...

Full text loading...


  1. ChongZeng, JianghaiXia, Richard D.Miller, et al.
    [2011] Application of the multiaxial perfectly matched layer(M-PML) to near-surface seismic modeling with Rayleigh waves. Geophysics, 76(3), T43–T52.
    [Google Scholar]
  2. HaiqiangLan, JingyiChen, YoushanLiu et al.
    [2013] Application of the perfectly matched layer in numerical modeling of wave propagation with an irregular free surface. 83rd Annual International Meeting, Expanded Abstracts, SEG, 3515–3520.
    [Google Scholar]
  3. Hestholm, S. and B.Ruud
    . [1994] 2D finite difference elastic wave modeling including surface topography. Geophysical Prospecting, 42(5), 371–390.
    [Google Scholar]
  4. Huang, J.P., Yang, Y., Li, Z.C. et al.
    [2014] Comparative study among implementations of several free-surface boundaries with perfectly matched layer conditions. Acta Seismologica Sinica, 36(5), 964–977.
    [Google Scholar]
  5. Komatitsch. D. and R.Martin
    , [2007] An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation. Geophysics, 72(5), SM155–SM167.
    [Google Scholar]
  6. Li, N., Huang, J.P., Li, Z.C. et al.
    [2014] The study on numerical simulation method of Lebedev Grid with dispersion improvement coefficients in TTI media. Chinese J. Geophys., 57(1), 261–269.
    [Google Scholar]
  7. Robertsson, J.O.A.
    [1996] A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography. Geophysics, 61(6), 1921–1934.
    [Google Scholar]
  8. Tian, K., Huang, J.P., Li, Z.C. et al.
    [2013] Simulation of the planar free surface in the finite-difference modeling of half-space viscoelastic medium. 75th EAGE Conference & Exhibition incorporating SPEEUROPEC, Tu-P13-03.
    [Google Scholar]
  9. Zhang, W. and Chen, X.
    [2006] Traction image method for irregular free surface boundaries in finite difference seismic wave simulation. Geophysical Journal International, 167(1), 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