We propose a set of local operators to deal with internal discontinuities that do not coincide with collocated Cartesian grids in 2D transversely isotropic media. We use globally optimally accurate operators in order to obtain high accuracy. We derive modified operators by extrapolating wavefields from nearby grid points to the discontinuous point, introducing boundary conditions; we then distribute those conditions to the nearby grid points. Numerical examples suggest that the operators improve the coherency of the wavefront. We would like to optimise the local operators to control all the error propagation during the modelling.


Article metrics loading...

Loading full text...

Full text loading...


  1. Cuvilliez, C., Fuji, N.
    , 2015. OptimallyAccurateFD1Dhttps://github.com/IPGP/OptimallyAccurateFD1D
    [Google Scholar]
  2. Geller, R.J., Mizutani, H., Hirabayashi, N.
    , 2012. Existence of a second island of stability of predictor-corrector schemes for calculating synthetic seismograms. Geophys. J. Int., 188, 253–262.
    [Google Scholar]
  3. Geller, R.J., Takeuchi, N.
    , 1998. Optimally accurate second-order time-domain finite difference scheme for the elastic equation of motion: one-dimensional cases. Geophys. J. Int., 123, 449–470.
    [Google Scholar]
  4. Käser, M., Dumbser, M.,
    2006. An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes – I. The two-dimensional isotropic case with external source terms. Geophys. J. Int., 166, 855–877.
    [Google Scholar]
  5. Komatitsch, D., Vilotte, J.P.
    , 1998. The spectral-element method: an efficient tool to simulate the seismic response of 2D and 3D geological structures, Bull. Seism. Soc. Am., 88, 368–392.
    [Google Scholar]
  6. Lombard, B., Piraux, J.
    , Numerical treatment of two-dimensional interfaces for acoustic and elas-tic waves. J. Comp. Phys., 195, 90–116.
    [Google Scholar]
  7. Madariga, R.
    , 1976. Dynamics of an expanding circular fault, Bull. Seism. Soc. Am., 65, 163–182.
    [Google Scholar]
  8. Moczo, P., Robertsson, J.O.A., Eisner, L.
    , The finite-difference time-domain method for modeling of seismic wave propagation. Advances in Geophysics, 48, 421–516.
    [Google Scholar]
  9. Mizutani, H.
    2002. Accurate and efficient methods for calculating synthetic seismograms when elastic discontinuities do not coincide with the numerical grid. PhD thesis, Univ. Tokyo.
    [Google Scholar]
  10. Schönberg, M., Muir, F.
    , A calculus for finely layered anisotropic media. Geophysics, 54, 581–589.
    [Google Scholar]
  11. Virieux, J.
    , 1984. SH-wave propagation in heterogeneous media: velocity-stress finite-difference method., Geophysics, 49, 1933–1942.
    [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