중간 규모 지진의 단층 파해 방향성 결정을 위한 새로운 주파수 영역 역산방법: 파쇄 전파 모델을 이용한 수치 시험

Anewwaveform inversion methodto determine the rupture directivity of moderate earthquakes: numerical tests for rupture models

  • 유승훈 (서울대학교 지구환경과학부) ;
  • 이준기 (서울대학교 지구환경과학부)
  • Yoo, Seung-Hoon (School of Earth and Environmental Science, Seoul National University) ;
  • Rhie, Jun-Kee (School of Earth and Environmental Science, Seoul National University)
  • 발행 : 2009.02.28

초록

단층 파쇄의 방향성은 지진 피해를 평가하는데 있어 매우 중요한 지진원 특성이다 하지만 기존의 모멘트 텐서 역산 방법으로는 단층의 파쇄 방향은 물론 실제 단층면의 방향도 정확하게 결정하기 어렵다. 본 연구에서는 중간 규모의 지진에 대하여 주파수 영역 파형 역산 방법을 이용하여 모멘트 텐서와 단층 파쇄의 방향성을 동시에 역산하는 방법을 제안하였다 여러 가지 다양한 파쇄 전파 모델을 가정한 수치 실험을 통해 역산 방법을 검증하였고, 실제 지진에 적응 가능성을 평가하기 위해 역산 해의 안정성에 가장 큰 영향을 주는 요소인 속도 구조 모델에 대한 민감도를 분석하였다. 민감도 분석 결과를 통해 속도 구조 모델이 실제 속도 구조와 크게 어긋나지 않을 경우 실제 지진에 대해서도 충분히 적응 가능하다는 것을 확인하였다. 향후 속도 구조가 비교적 잘 밝혀진 지역에 본 역산 방법을 적응 할 경우 중간 규모 지진의 단층 파쇄 효과를 효과적으로 추정할 수 있을 것으로 예상되며, 이를 통해 적용된 지역의 지진 발생 특성을 이해하는데 큰 도움을 줄 수 있을 것이다.

Rupture directivity is the important parameter in estimating damage due to earthquakes. However, the traditional moment tensor inversion technique cannot resolve the real fault plane or the rupture directivity. To overcome these limitations, we have developed a new inversion algorithm to determine the moment tensor solution and the rupture directivity for moderate earthquakes, using the waveform inversion technique in the frequency domain. Numerical experiments for unilateral and bilateral rupture models with various rupture velocities confirm that the method can resolve the ambiguity of the fault planes and the rupture directivity successfully. To verify the feasibility of the technique, we tested the sensitivity to velocity models, which must be the most critical factor in practice. The results of the sensitivity tests show that the method can be applied even though the velocity model is not perfect. If this method is applied in regions where the velocity model is well verified, we can estimate the rupture directivity of a moderate earthquake. This method makes a significant contribution to understanding the characteristics of earthquakes in those regions.

키워드

참고문헌

  1. Aki, K., and Richards, P. G., 2002, Quantitative Seismology, 2nd Edition, University Science Books, Sausalito, 700 pp
  2. Bertsekas, D. P., 1982, Enlarging the region of convergence of Newton’s method for constrained optimization: Journal of Optimization Theory and Applications, 36, 221–251. doi: 10.1007/BF00933831
  3. Dreger, D. S., and Helmberger, D. V., 1993, Determination of source parameters at regional distance with three-component sparse network data: Journal of Geophysical Research, 98, 8107–8125. doi: 10.1029/ 93JB00023
  4. Dziewonski, A. M., and Woodhouse, J. H., 1983,Anexperiment in systematicstudy of global seismicity: centroid moment solutions for 201 moderate and large earthquakes: Journal of Geophysical Research, 88, 3247–3271. doi: 10.1029/JB088iB04p03247
  5. Dziewonski, A. M., Chou, T. A., and Woodhouse, J. H., 1981, Determinationof earthquake source parameters from waveform data for studies of global and regional seismicity: Journal of Geophysical Research, 86, 2825–2852. doi: 10.1029/JB086iB04p02825
  6. Hartzell, S. H., and Helmberger, D. V., 1982, Strong-motion modeling of theImperial Valley earthquake of 1979: Bulletin of the Seismological Society of America, 72, 571–596
  7. Hartzell, S. H., and Heaton, T. H., 1983, Inversion of strong ground motionand teleseismic waveform data for the fault rupture history of the 1979 Imperial Valley, California, earthquake: Bulletin of the Seismological Society of America, 73, 1553–1583
  8. Haskell, N. A., 1964, Total energy and energy spectral density of elasticwave radiation from propagating faults: Bulletin of the Seismological Society of America, 54, 1811–1841
  9. Lay, T., and Wallace, T. C., 1995, Modern Global Seismology, Academic Press
  10. Patton, H. J., and Walter, W. R., 1993, Regional moment: magnitude relationsfor earthquakes and explosions: Geophysical Research Letters, 20, 277–280. doi: 10.1029/93GL00298
  11. Saikia, C. K., 1994, Modified frequency-wavenumber algorithm for regionalseismograms using Filon’s quadrature; modeling of Lg waves in Eastern North America: Geophysical Journal International, 118, 142–158. doi: 10.1111/j.1365-246X.1994.tb04680.x
  12. Shakal, A., Haddadi, H., Graizer, V., Lin, K., and Huang, M., 2006,Some KeyFeatures of the Strong-Motion Data from theM6.0 Parkfield, California, Earthquake of 28 September 2004: Bulletin of the Seismological Society of America, 96, s90–s118. doi: 10.1785/0120050817
  13. Shin, C. S., Pyun, S., and Bednar, J. B., 2007, Comparison of waveforminversion, part1: conventional wave field vs logarithmic wave field: Geophysical Prospecting, 55, 449–464. doi: 10.1111/j.1365-2478. 2007.00617.x
  14. Tarantola, A., 1987, Inverse Problem Theory: Methods for Data Fitting andParameter Estimation, Elsevier
  15. Zhao, L. S., and Helmberger, D. V., 1994, Source estimation from broadbandregional seismograms: Bulletin of the Seismological Society of America, 84, 91–104