ARTICLE

Time-domain formulation of a perfectly matched layer for the second-order elastic wave equation with VTI media

JAEJOON LEE1 CHANGSOO SHIN2
Show Less
1 Computational Science and Technology Program, Seoul National University, Seoul, South Korea.,
2 Department of Energy System Engineering, Seoul National University, Seoul, South Korea.,
JSE 2015, 24(3), 231–257;
Submitted: 9 June 2025 | Revised: 9 June 2025 | Accepted: 9 June 2025 | Published: 9 June 2025
© 2025 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Lee, J. and Shin, C.S., 2015. Time-domain formulation of a perfectly matched layer for the second-order elastic wave equation with VTI media. Journal of Seismic Exploration, 24: 231-257. In this study, we introduce the unsplit Perfectly Matched Layer (PML) for the 2D and 3D second-order elastic wave equations with isotropic and transversely isotropic vertical axis of symmetry (VTI) media in the time domain. The introduced PML formulations are successfully applied to practical applications in terms of efficiency and stability. The PML formulations require less than or an equal number of auxiliary variables than other formulations, thereby decreasing the computational power necessary to calculate the solution in the PML zone. Derived directly from the second-order wave form, the PML formulation demonstrates an improved stability compared to first-order PMLs or second-order PMLs that are derived from first-order systems. Numerical examples demonstrate that the bulk waves and strong surface waves are perfectly damped out without introducing instability for an isotropic material in both 2D and 3D. The derived formulation also provides effective absorption with strong VTI materials, including zinc and apatite, that cause instability problems in other PML formulations.

Keywords
perfectly matched layer
time-domain seismic wave modeling
second-order elastic wave equation
VTI media
References
  1. Abarbanel, S. and Gottlieb, D., 1997. A mathematical analysis of the PML method. J. Computat.
  2. Phys., 134: 357-363.
  3. Abarbanel, S., Gottlieb, D. and Hesthaven, J.S., 2002. Long time behavior of the perfectly matched
  4. layer equations in computational electromagnetics. J. Scientif. Comput., 17: 405-422.
  5. Appelé, D. and Petersson, N.A., 2009. A stable finite difference method for the elastic wave
  6. equation on complex geometries with free surfaces. Communic. Computat. Phys. , 5: 84-107.
  7. Assi, H. and Cobbold, R.S., 2013. Perfectly matched layer for second-order time-domain elastic
  8. wave equation: formulation and stability. arXiv preprint arXiv: 1312.3722.
  9. Bécache, E., Fauqueux, S. and Joly, P., 2003. Stability of perfectly matched layers, group velocities
  10. and anisotropic waves. J. Computat. Phys., 188: 399-433.
  11. PML FOR SECOND-ORDER ELASTIC EQUATION 255
  12. Berenger, J.P., 1994. A perfectly matched layer for the absorption of electromagnetic waves. J.
  13. Computat. Phys., 114: 185-200.
  14. Carcione, J.M., Kosloff, D. and Kosloff, R., 1988. Wave-propagation simulation in an elastic
  15. anisotropic (transversely isotropic) solid. Quart. J. Mechan. Appl. Mathemat., 41: 319-346.
  16. Chang, W.F. and McMechan, G.A., 1987. Elastic reverse-time migration. Geophysics, 52:
  17. 1365-1375.
  18. Chew, W.C. and Liu, Q.H., 1996. Perfectly matched layers for elastodynamics: a new absorbing
  19. boundary condition. J. Computat. Acoust., 4: 341-359.
  20. Chew, W.C. and Weedon, W.H., 1994. A 3D perfectly matched medium from modified Maxwell’s
  21. equations with stretched coordinates. Microw. Opt. Technol. Lett., 7: 599-604.
  22. Collino, F., and Tsogka, C., 2001. Application of the perfectly matched absorbing layer model to
  23. the linear elastodynamic problem in anisotropic heterogeneous media. Geophysics, 66:
  24. 294-307.
  25. Crase, E., Pica, A., Noble, M., McDonald, J. and Tarantola, A., 1990. Robust elastic nonlinear
  26. waveform inversion: Application to real data. Geophysics, 55: 527-538.
  27. De Basabe, J.D., Sen, M.K. and Wheeler, M.F., 2008. The interior penalty discontinuous Galerkin
  28. method for elastic wave propagation: grid dispersion. Geophys. J. Internat., 175: 83-93.
  29. Duru, K., and Kreiss, G., 2012. A well-posed and discretely stable perfectly matched layer for
  30. elastic wave equations in second order formulation. Communic. Computat. Phys., 11: 1643.
  31. Grote, M.J. and Sim, I., 2010. Efficient PML for the wave equation, arXiv preprint arXiv:
  32. 0319.
  33. Hastings, F.D., Schneider, J.B. and Broschat, S.L., 1996. Application of the perfectly matched
  34. layer (PML) absorbing boundary condition to elastic wave propagation. J. Acoust. Soc. Am.,
  35. 100: 3061-3069.
  36. Kim, Y., Min, D.J. and Shin, C., 2011. Frequency-domain reverse-time migration with source
  37. estimation. Geophysics, 76: S41-S49.
  38. Komatitsch, D. and Martin, R., 2007. An unsplit convolutional perfectly matched layer improved
  39. at grazing incidence for the seismic wave equation. Geophysics, 72: SM155-SM167.
  40. Komatitsch, D. and Vilotte, J.P., 1998. The spectral element method: an efficient tool to simulate
  41. the seismic response of 2D and 3D geological structures. Bull. Seismol. Soc. Am., 88:
  42. 368-392.
  43. Komatitsch, D. and Tromp, J., 2003. A perfectly matched layer absorbing boundary condition for
  44. the second-order seismic wave equation. Geophys. J. Internat., 154: 146-153.
  45. Kosloff, D., Queiroz Filho, A., Tessmer, E. and Behle, A., 1989. Numerical solution of the
  46. acoustic and elastic wave equations by a new rapid expansion method. Geophys. Prosp., 37:
  47. 383-394.
  48. Ledbetter, H.M., 1977. Elastic properties of zinc: A compilation and a review. J. Phys. Chem. Ref.
  49. Data, 6: 1181-1203.
  50. Li, Y. and Matar, O.B., 2010. Convolutional perfectly matched layer for elastic second-order wave
  51. equation. J. Acoust. Soc. Am., 127: 1318-1327.
  52. Liu, Q.H. and Tao, J., 1997. The perfectly matched layer for acoustic waves in absorptive media.
  53. J. Acoust. Soc. Am., 102: 2072-2082.
  54. Marfurt, K.J., 1984. Accuracy of finite-difference and finite-element modeling of the scalar and
  55. elastic wave equations. Geophysics, 49: 533-549.
  56. Meza-Fajardo, K.C. and Papageorgiou, A.S., 2008. A nonconvolutional, split-field, perfectly
  57. matched layer for wave propagation in isotropic and anisotropic elastic media: Stability
  58. analysis. Bull. Seismol. Soc. Am., 98: 1811-1836.
  59. Mora, P., 1987. Nonlinear two-dimensional elastic inversion of multioffset seismic data. Geophysics,
  60. 52: 1211-1228.
  61. Pratt, R.G., Shin, C. and Hick, G.J., 1998. Gauss-Newton and full Newton methods in frequency-
  62. space seismic waveform inversion. Geophys. J. Internat., 133: 341-362.
  63. Qi, Q. and Geers, T.L., 1998. Evaluation of the perfectly matched layer for computational
  64. acoustics. J. Computat. Phys., 139: 166-183.
  65. 256 LEE & SHIN
  66. Shin, C. and Cha, Y.H., 2008. Waveform inversion in the Laplace domain. Geophys. J. Internat.,
  67. 173: 922-931.
  68. Shin, C. and Cha, Y.H., 2009. Waveform inversion in the Laplace-Fourier domain. Geophys. J.
  69. Internat., 177: 1067-1079.
  70. Sun, R. and McMechan, G.A., 2001. Scalar reverse-time depth migration of prestack elastic seismic
  71. data. Geophysics, 66: 1519-1527.
  72. Tarantola, A., 1986. A strategy for nonlinear elastic inversion of seismic reflection data.
  73. Geophysics, 51: 1893-1903.
  74. Thomsen, L., 1986. Weak elastic anisotropy. Geophysics, 51: 1954-1966.
  75. Yan, J. and Sava, P., 2008. Isotropic angle-domain elastic reverse-time migration. Geophysics, 73:
  76. $229-S239.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing