ARTICLE

A non-split perfectly matched layer absorbing boundary condition for the second-order wave equation modeling

C. ZHANG1,3 B. SUN1 H. YANG1 J. MA2
Show Less
1 School of Aerospace, Tsinghua University, Beijing, P.R. China. bbsunok@126.com,
2 Department of Mathematics, Harbin Institute of Technology, Harbin, P.R. China.,
3 Department of Geophysics, Research Institute of Petroleum Exploration & Development, Beijing, P.R. China.,
JSE 2016, 25(6), 513–525;
Submitted: 2 March 2015 | Accepted: 29 September 2016 | Published: 1 December 2016
© 2016 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

Zhang, C., Sun, B., Yang, H. and Ma, J., 2016. A non-split perfectly matched layer absorbing boundary condition for the second-order wave equation modeling. Journal of Seismic Exploration, 25: 513-525. The perfectly matched layer (PML) absorbing boundary conditions (ABC) have been well studied for seismic wavefield modeling. However, existing approaches are either based on a wavefield variable-split or are only applicable to first-order wave equations. In this paper, we present a non-split PML ABC for the second order wave equations in displacement. The principle of the proposed method lies in introducing a series of auxiliary variables to represent the partial derivatives associated with the stretching axis. We derive the non-split PML formula by using basic tensor algebra. The derived equations are in a compact form which makes the ABC condition easier for implementation in practice. Furthermore, as no extra splitting wavefield variables are introduced, the computer memory usage of wave equation modeling can be reduced accordingly. Numerical results for both acoustic and elastic examples show the quality of the performance of the proposed new method.

Keywords
perfectly matched layer
absorbing boundary condition
second-order
wave equation
seismic modeling
finite difference
References
  1. Abarbanel, S. and Gottlieb, D., 1997. A mathematical analysis of the PML method. J. Comput.Phys., 134: 357-363.
  2. Aki, K. and Richards, P., 1980. Quantitative Seismology: Theory and Methods. W.H. Freeman andCo., San Francisco.
  3. Basu, U. and Chopra, A., 2003. Perfectly matched layers for time-harmonic elastodynamics ofunbounded domains: theory and finite-element implementation. Comput. Meth. Appl. Mech.Engin., 192: 1337-1375.
  4. Baysal, E., Kosloff, D. and Sherwood, J., 1983. Reverse time migration. Geophysics, 48:1514-1524.
  5. Berenger, J.P., 1994. A perfectly matched layer for the absorption of electro-magnetic waves. J.Comput. Phys., 114: 185-200.
  6. Cerjan, C., Kosloff, D., Kosloff, R. and Resef, M., 1985. A non reflecting boundary condition fordiscrete acoustic and elastic wave equations. Geophysics, 50: 705-708.
  7. Clayton, R. and Engquist, B., 1977. Absorbing boundary conditions for acoustic and elastic waveequations. Bull. Seismol. Soc. Am., 67: 1529-1540.
  8. Chew, W. and Liu, Q., 1996. Perfectly matched layers for elastodynamics: a new absorbingboundary condition. J. Comput. Acoust., 4: 341-359.
  9. Collino, F. and Monk, P., 1998. Optimizing the perfectly matched layer. Comput. Meth. Appl.Mech. Eng., 164: 157-171.
  10. Collino, F. and Tsogka, C., 2001. Application of the PML absorbing layer model to the linearelastodynamic problem in anisotropic heterogeneous media. Geophysics, 66: 294-307.
  11. Dablain, M., 1986. The application of high-order differencing to the scalar wave equationGeophysics, 51: 54-66.
  12. Drossaert, F. and Giannopoulos, A., 2007. A nonsplit complex frequency shifted PML based onrecursive integration for FDTD modeling of elastic waves. Geophysics, 72(2): T9-T17.LAYER ABSORBING BOUNDARY CONDITION 523
  13. Festa, G. and Vilotte, J., 2005. The Newmark scheme as velocity-stress time staggering: Anefficient PML implementation for spectral-element simulations of elastodynamics. Geophys.J. Internat., 161: 789-812.
  14. Gedney, S., 1996. An anisotropic PML absorbing media for the FDTD simulation of fields in flossyand dispersive media. Electromagnet., 16: 399-415.
  15. Hastings, F., Schneider, J. and Broschat, S., 1996. Application of the perfectly matched layer(PML) absorbing boundary condition to elastic wave propagation. J. Acoust. Soc. Am., 100:3061-3069.
  16. Hagstrom, T. and Hariharan, S., 1998. A formulation of asymptotic and exact boundary conditionsusing local operators. Appl. Numer. Math., 27: 403-416.
  17. Kelly, K., Ward, R., Treitel, S. and Alford, R., 1976. Synthetic seismograms: A finite-differenceapproach. Geophysics, 41: 2-27.
  18. Komatitsch, D. and Martin, R., 2007. An unsplit convolutional perfectly matched layer improvedat grazing incidence for the seismic wave equation. Geophysics, 72(5): SM155-SM167.
  19. Komatitsch, D. and Tromp, J., 2003. A perfectly matched layer absorbing boundary condition forthe second-order seismic wave equation. Geophys. J. Internat., 154: 146-153.
  20. Liu, J., Ma, J. and Yang, H., 2009. The study of perfectly matched layer absorbing boundaries forSH wave fields. Appl. Geophys., 6: 267-274.
  21. Liu, Q. and Tao, J., 1997. The perfectly matched layer for acoustic waves in absorptive media. J.Acoust. Soc. Am., 102: 2072-2082.
  22. Liu, Y. and Sen, M.K., 2010. A hybrid scheme for absorbing edge reflections in numericalmodeling of wave propagation. Geophysics, 75(2): Al1-A6.
  23. Loewenthal, D. and Mufti, I., 1983. Reverse time migration in spatial frequency domain.Geophysics, 48: 627-635.
  24. McMechan, G.A., 1983. Migration by extrapolation of time-dependent boundary values. Geophys.Prosp., 31: 413-420.
  25. Mur, G., 1981. Absorbing boundary conditions for the finite-difference approximation of thetime-domain electromagnetic field equations. IEEE Transact. Electromagn. Compat., 23:377-382.
  26. Qi. Q. and Geers, T., 1998. Evaluation of the perfectly matched layer for computational acoustics.J. Comput. Phys., 139: 166-183.
  27. Ramadan, O., 2003. Auxiliary differential equation formulation: An efficient implementation of theperfectly matched layer. IEEE Microwave Wireless Compon. Lett., 13: 69-71.
  28. Sacks, Z., Kingsland, D., Lee, R. and Lee, J., 1995. A perfectly matched anisotropic absorber foruse as an absorbing boundary condition. IEEE Transact. Antenn. Propag., 43: 1460-1463.
  29. Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation. Geophysics,49: 1259-1266.
  30. Teixeira, F. and Chew, W., 2000. Complex space approach to perfectly matched layers: A reviewand some new developments. Internat. J. Numer. Model. - Electron. Netw., Devices Fak13: 441-455.
  31. Virieux, J., 1984. SH-wave propagation in heterogeneous media: Velocity stress finite-difference: method. Geophysics, 49: 1933-1957.
  32. Wang, L. and Liang, C., 2006. A new implementation of CFS-PML for ADI-FDTD method.Microw. Optic. Technol. Lett., 48: 1924-1928.
  33. Wang, T. and Tang, X., 2003. Finite-difference modeling of elastic wave propagation: Anon-splitting perfectly matched layer approach. Geophysics, 68: 1749-1755.
  34. Whitmore, D. 1983. Iterative depth imaging by back time propagation. Expanded Abstr., 53rd Ann.Internat. SEG Mtg., Las Vegas: 382-385.
  35. Zeng, Y. and Liu, Q., 2004. A multidomain PSTD method for 3D elastic wave equations. Bull.Seismol. Soc. Am., 94: 1002-1015.
  36. Zhang, W. and Shen, Y., 2010. Unsplit complex frequency-shifted PML implementation usingauxiliary differential equations for seismic wave modeling. Geophysics, 75(4): T141-T154.524 ZHANG, SUN, YANG & MA
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing