Application and Improvement of Complex Frequency Shifted Perfectly Matched Layers for Elastic Wave Modeling in the Frequency-domain

주파수영역 탄성파모델링에 대한 CFS-PML경계조건의 적용 및 개선

  • Son, Min-Kyung (Earthquake Research Center, Korea Institute of Geoscience and Mineral Resources) ;
  • Cho, Chang-Soo (Earthquake Research Center, Korea Institute of Geoscience and Mineral Resources)
  • 손민경 (한국지질자원연구원 지진연구센터) ;
  • 조창수 (한국지질자원연구원 지진연구센터)
  • Received : 2012.05.17
  • Accepted : 2012.07.23
  • Published : 2012.08.31


Absorbing boundary conditions are used to mitigate undesired reflections that can arise at the model's truncation boundaries. We apply a complex frequency shifted perfectly matched layer (CFS-PML) to elastic wave modeling in the frequency domain. Modeling results show that the performance of our implementation is superior to other absorbing boundaries. We consider the coefficients of CFS-PML to be optimal when the kinetic energy becomes to the minimum, and propose the modified CFS-PML that has the CFS-PML coefficient ${\alpha}_{max}$ defined as a function of frequency. Results with CFS-PML and modified CFS-PML are significantly improved compared with those of the classical PML technique suffering from large spurious reflections at grazing incidence.


Grant : 한반도 동부지역 잠재적 위험 지진원 평가

Supported by : 한국지질자원연구원


  1. Bao, H., Bielak, J., Ghattas, O., Kallivokas, L. F., O'Hallaron, D. R., Shewchuk, J. R., and Xu, J., 1998, Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers, Computer Methods in Applied Mechanics and Engineering, 152, 85-102.
  2. Berenger, J., 1994, A perfectly matched layer for the absorption of electromagnetic waves, Journal of Computational Physics, 114, 185-200.
  3. Berenger, J., 2002, Application of the CFS PML to the absorption of evanescent waves in waveguides, IEEE Microwave and Wireless Components Letters, 12, 218-220.
  4. Berg, P., If, F., Nielsen, P., and Skovgaard, O., 1994, Analytical reference solutions, in Helbig, K., Ed., Modeling the Earth for Oil Exploration, Pergamon Press, 421-427.
  5. Cerjan, C., Kosloff, D., Ksloff, R., and Reshef, M., 1985, A nonreflecting boundary condition ofr discrete acoustic and elastic wave equation, Geophysics, 50, 705-708.
  6. Chew, W., and Liu, Q., 1996, Perfeclty matched layers for elastodynamics: A new absorbing boundary condition, Journal of Computational Acoustics, 4, 341-359.
  7. Cho, C. S., and Lee, H. I., 2009, Application of convolutional perfectly matched layer method to numerical elastic modeling using rotated staggered grid, Jigu-Mulli-wa-Mulli-Tamsa, 12, 183-191.
  8. Clayton, R., and Engquist, B., 1977, Absorbing boundary conditions for acoustic and elastic wave equations, Bulletin of the Seismological Society of America, 67, 1529-1540.
  9. Collino, F., and Monk, P., 1998, Optimizing the perfectly matched layer, Computer Methods in Applied Mechanics and Engineering, 164, 157-171.
  10. Collino, F., and Tsogka, C., 2001, Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media, Geophysics, 66, 294-307.
  11. Correia, D. and Jin, J.-M., 2005, On the Development of a Higher-Order PML, IEEE Trans. Antennas Propag., 53, 4157- 4163.
  12. Drossaert, F. H., and Giannopoulos, A., 2007a, A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves, Geophysics, 72, T9-T17.
  13. Drossaert, F. H., and Giannopoulos, A., 2007b, Complex frequency shifted convolution PML for FDTD modelling of elastic waves, Wave Motion, 44, 593-604.
  14. Higdon, R., 1991, Absorbing boundary conditions for elastic waves, Geophysics, 56, 231-241.
  15. Komatitsch, D., and Martin, R., 2007, An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation, Geophysics, 72, SM155-SM167.
  16. Komatitsch, D., and Tromp, J., 1999, Introduction to the spectral element method for three-dimensional seismic wave propagation, Geophysical Journal International, 139, 806-822.
  17. Pratt, R. G., and Worthington, M. H., 1990, Inverse theory applied to multi-source cross-hole tomography. Part 1: Acoustic wave-equation method, Geophysical Prospecting, 38, 287- 310.
  18. Roden, J. A., and Gedney, S. D., 2000, Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML for arbitrary media, Microwave and Optical Technology Letters, 27, 334-339.<334::AID-MOP14>3.0.CO;2-A
  19. Shin, C., 1995, Sponge boundary condition for frequency-domain modeling, Geophysics, 60, 1870-1874.
  20. Tarantola, A., 1984, Inversion of seismic reflection data in the acoustic approximation, Geophysics, 49, 1259-1266.
  21. Virieux, J., 1986, P-SV wave propagation in heterogeneous media: velocity-stress finite-difference method, Geophysics, 51, 889-901.
  22. Virieux, J., and Operto, S., 2009, An overview of full-waveform inversion in exploration geophysics, Geophysics, 74, WCC127-WCC152.
  23. Zhang, W., and Shen, Y., 2010, Unsplit complex frequencyshifted PML implementation using auxiliary differential equations for seismic wave modeling, Geophysics, 75, T141-T154.