1887

Abstract

Summary

We extend the Robertson’s image method ( ) to a 3D finite-difference time-domain (FDTD) wave propagation in the presence of surface topography. Robertson method is accurate for a horizontal free surface and easily implemented on a standard staggered grid (SSG), but has some artefact for a topographic surface if it is staircase represented. We use a surface-conforming mapping to transfer the physical topographic domain (x,y,z) to a regular Cartesian domain (ξ,κ,η), which is equivalent to a curvilinear coordinate system. We take the velocity-stress first-order equation form for the (anisotropic) elastic wave equation, and discretize the transferred Cartesian domain on a fully staggered grid (FSG). Following the Robertsson procedure, we put all the particle velocity zero above the free surface, the velocity derivative with respect to η at the free surface is replaced with velocity derivatives with respect to ξ and κ, and anti-symmetrize the traction component normal to the curvilinear surface. We present the numerical results to show the efficiency of this extension.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201412694
2015-06-01
2024-04-28
Loading full text...

Full text loading...

References

  1. Graves, Robert W.
    , 1996, Simulating Seismic Wave Propagation in 3D Elastic Media Using Staggered-Grid Finite Differences, Bulletin of the Seismological Society of America, Vol. 86, No. 4, 1091–1106.
    [Google Scholar]
  2. Hestholm, S.
    , 1999, Three-dimensional finite difference viscoelastic wave modelling including surface topography: Geophysical Journal International, 139, 852–878.
    [Google Scholar]
  3. Hestholm, S., and B.Ruud
    , 2002, 3D free-boundary conditions for coordinate-transform finite-difference seismic modelling: Geophysical Prospecting, 50, 463–474.
    [Google Scholar]
  4. Komatitsch, D., and R.Martin
    , 2007, An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation: Geophysics, 72, no. 4, SM155–SM167.
    [Google Scholar]
  5. Kristek, J., P.Moczo, and R. J.Archuleta
    , 2002, Efficient methods to simulate planar free surface in the 3D 4th-order staggered-grid finite difference schemes: Studia Geophysica et Geodaetica, 46, 355–381.
    [Google Scholar]
  6. Levander, A. R.
    , 1988, Fourth-order finite-difference P-SV seismograms: Geophysics, 53, 1425–1436.
    [Google Scholar]
  7. Moczo, P., E.Bystricky, J.Kristek, J.Carcione, and M.Bouchon
    , 1997, Hybrid modeling of P-SV seismic motion at inhomogeneous viscoelastic topographic structures: Bulletin of the Seismological Society of America, 87, 1305–1323.
    [Google Scholar]
  8. Josepde la Puente, MiguelFerrer, MauricioHanzich, José E.Castillo, and José M.Cela
    , 2014, Mimetic seismic wave modeling including topography on deformed staggered grids: GEOPHYSICS, VOL. 79, T125–T141.
    [Google Scholar]
  9. Robertsson, O. A. Johan
    , 1996, A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography: Geophysics61, 1921–1934.
    [Google Scholar]
  10. Zhang, W., and X.Chen
    , 2006, Traction image method for irregular free surface boundaries in finite difference seismic wave simulation: Geophysical Journal International, 167, 337–353
    [Google Scholar]
  11. Zhang, W., Z.Zhang, and X.Chen
    , 2012, Three-dimensional elastic wave numerical modelling in the presence of surface topography by a collocated-grid finite-difference method on curvilinear grids: Geophysical Journal International, 190, 358–378.
    [Google Scholar]
  12. Zeng, C., J.Xia, R. D.Miller, and G. P.Tsoflias
    , 2012, An improved vacuum formulation for 2D finite-difference modeling of Rayleigh waves including surface topography and internal discontinuities: Geophysics, 77, no. 1, T1–T9.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201412694
Loading
/content/papers/10.3997/2214-4609.201412694
Loading

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