Application of Convolutional Perfectly Matched Layer Method to Numerical Elastic Modeling Using Rotated Staggered Grid

회전된 엇갈린 격자를 이용한 탄성파 모델링에의 CPML 경계조건 적용

  • Cho, Chang-Soo (Earthquake Research Center, Korea Institute of Geoscience and Mineral Resources) ;
  • Lee, Hee-Il (Earthquake Research Center, Korea Institute of Geoscience and Mineral Resources)
  • 조창수 (한국지질자원연구원 지진연구센터) ;
  • 이희일 (한국지질자원연구원 지진연구센터)
  • Published : 2009.05.28

Abstract

Finite difference method using not general SSG (standard staggered grid) but RSG (rotated staggered grid) was applied to simulation of elastic wave propagation. Special free surface boundary condition such as imaging method is needed in finite difference method using SSG in elastic wave propagation. But free surface boundary condition in finite difference method using RSG is easily solved with adding air layer or vacuum layer. Recently PML (Perfectly Matched layer) is widely used to eliminate artificial reflection waves from finite boundary because of its' greate efficiency. Absorbing ability of CPML (convolutional Perfectly Matched Layer) that is more efficient than that of PML and CPML that don't use splitting of wave equation that should be adapted to PML was applied to FDM using RSG in this study. Frequency absorbing characteristic and energy absorbing ability in CPML layer were investigated and CPML eliminated artificial boundary waves very effectively in FDM using RSG in being compared with that of Cerjan's absorbing method. CPML method also diminished amplitude of waves in boundary layer of solid-liquid model very well.

탄성파 수치 모형 계산에 있어서 널리 사용되는 엇갈린 격자 방법이 아니라 회전된 엇갈린 격자 방법을 사용하여 탄성파 수치 모사를 수행하였다. 표준 엇갈린 격자 방법에서는 특별한 자유 경계조건을 적용하여야 하는 단점이 있지만 회전된 격자 방법에서는 물성으로 진공 또는 공기층을 부여함으로써 자유 경계조건을 실현가능하다는 것을 확인할 수 있었다. 파동전파에 있어서 유한 경계 조건에서 발생하는 인공 반사파를 제거하기 위해 PML (Perfectly Matched Layer)의 파동식 분해라는 단점을 극복할 수 있고 좋은 성능을 보이는 CPML (Convolutional Perfectly Matched Layer)법을 회전된 엇갈린 격자법(RSG: Rotatged Staggered Grid)에 적용하였다. 회전된 격자 유한 차분법에서 CPML의 고주파수 흡수 특성과 에너지 흡수율 조사, Cerjan법의 감쇠를 비교한 결과 흡수경계조건으로 좋은 성능을 확인하였다. 유체와 고체의 모형에 대한 경계에 대하여서도 매우 효과적으로 경계면에서 발생하는 반사파를 제거할 수 있음을 알 수 있었다.

Keywords

References

  1. Bohlen, T. and Saenger, E. H., 2006, Accuracy of Heterogeneous Staggered-grid finite-difference modeling of Rayleigh Waves, Geophysics, 71, T109-T115
  2. Bohlen, T. and Saenger, E. H., 2006, Accuracy of Heterogene-ous Staggered-grid finite-difference modeling of Rayleigh Waves, Geophysics, 71, T109-T115
  3. Cerjan, C., Kosloff, D., Kosloff, R. and Reshef, M., 1985, A nonreflecting boundary condition for discrete acoustic and elastic wave equations, Geophysics, 50, 705-708
  4. Chen, H., Wang, X. and Zhao, H., 2006, A rotated staggered grid finite difference with the aborbing boundary condition of a perfectly matched layer, Chinese Science Bulletin, 51, 2304-2314
  5. Correia, D. and Jin, J.-M., 2005, On the Development of a Higher-Order PML, IEEE Trans. Antennas Propag., 53, 4157-4163
  6. Graves, R. W. 1996, Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences, BSSA, 86, 1091-1106
  7. Ewing, W. M., Tardetzky, W. S., and Press, F., 1957, Elastic Waves in Layered Media: McGraw-Hill
  8. Levander, A., 1988, Fourth-order finite-difference P-SV seismo-grams, Geophysics, 53, 1425-1436
  9. 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
  10. Roden, J. A. and Gedney, S. D., 2000, Convolution PML(CPML): An efficcient FDTD implementation of the CFS-PML for arbitrary media, Microwave and Optical Technology Letters, 27, 334-339
  11. Saenger, E. H., Gold, N. and Shapiro, S. A., 2000, Modeling the propagation of elastic waves using a modified finite-difference grid, Wave Motion, 31, 77-92