엇격자 유한차분법을 이용한 극지해역 지진파 모델링 = Earthquake Wave Propagation Using Staggered-grid Finite-difference Method in the Model of the Antarctic Region

다국어 초록 (Multilingual Abstract)

We simulate the propagation of earthquake waves in the continental margin of Antarctica using the elastic wave modeling algorithm, which is modified to be suitable for acoustic-elastic coupled media and earthquake source. To simulate the various types of earthquake source, the staggered-grid finite-difference method, which is composed of velocity-stress formulae, can be more appropriate to use than the conventional, displacement-based, finite-difference method. We simulate the elastic wave propagation generated by earthquakes combining 3D staggered-grid finite-difference algorithm composed of displacement-velocity-stress formulae with double couple mechanisms for earthquake source.
Through numerical tests for left-lateral strike-slip fault, normal fault and reverse fault, we could confirm that the first arrival of P waves at the surface is in a good agreement with the theoretically-predicted results based on the focal mechanism of an earthquake. Numerical results for a model made after the subduction zone in the continental margin of Antarctica showed that earthquake waves, generated by the reverse fault and propagating through the continental crust, the oceanic crust and the ocean, are accurately described.

국문 초록 (Abstract)

이 연구에서는 기존의 탄성파 모델링 알고리즘에 지진 송신원을 적용하고, 음향-탄성파 결합 매질을 구현하여남극대륙 주변과 같은 극지해역에서 발생할 수 있는 지진파의 거동을 모사한다. 기존의 변위근사 유한차분법과 달리 속도-응력 식으로 구성되는 엇격자 유한차분법의 경우 다양한 송신원을 구현하는데 적합하므로 변위-속도-응력 식에 기초하여 개발된 3차원 엇격자 유한차분법 알고리즘과 이중 우력(Double Couple Forces)을 이용하여 구현한 지진 송신원을접목시켜 지진파의 거동을 모사한다. 좌수향 주향이동단층, 정단층, 역단층 형태의 지진 송신원에 대해서 개발된 알고리즘을 검증한 결과 이론적으로 예측되는 P파의 초동을 정확히 모사할 수 있었고, 섭입대 모델에 대한 수치모형실험 결과 섭입대에서 역단층에 의해 발생된 후 대륙지각, 해양지각 및 해양에서 전파되는 지진파의 거동양상이 정확하게 모사되는 것을 확인할 수 있었다.

참고문헌 (Reference)

1 조창수, "회전된 엇갈린 격자를 이용한 탄성파 모델링에의 CPML 경계조건 적용" 한국지구물리.물리탐사학회 12 (12): 183-191, 2009

2 최윤석, "탄성파 및 음향-탄성파 결합 파동방정식 유한요소 모델링 기법을 이용한 주파수 영역 파형역산" 서울대학교 2007

3 백진주, "최근 5년간 한반도 중서부 지역에서 발생한 지진의 진원 특성 분석" 한국지구과학회 32 (32): 161-169, 2011

4 양준모, "청양 지자기관측소에서 획득된 지자기전달함수와 분극값의 시간변동성에 대한 연구" 한국지구과학회 30 (30): 824-833, 2009

5 오석훈, "안동지진 (2009년 5월 2일) 발생 기간 지자기장 자료의 변동성 분석" 한국지구과학회 30 (30): 638-691, 2009

6 Fornberg, B, "The pseudospectral method: Comparisons with finite differences for the elastic wave equation" 52 : 483-501, 1987

7 Fulbert, F, "The Darmstadt VLF radiophysical station VADAR for the detection of possible earthquake precursors" 13 : 2713-, 2011

8 Kelly, K.R, "Synthetic seismograms: A finite-difference approach" 53 : 1045-1055, 1976

9 신성렬, "Staggered Grid를 이용한 유 한차분법 탄성파 모델링" 34 : 168-174, 1997

10 Moczo, P, "Stability and grid dispersion of the P-SV 4th-order staggered-grid finite-difference schemes" 44 : 381-402, 2000

11 Graves, W.R, "Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences" 86 : 1091-1106, 1996

12 Silva, H.G, "Seismo-electromagnetic phenomena in the westerm part of the Eurasia-Nubia plate boundary" 11 : 241-248, 2011

13 Virieux, J, "SH-wave propagation in heterogeneous media: Velocity-stress finite-difference method" 49 : 1933-1942, 1984

14 Aki, K, "Quantitative seismology" University Science Books Sausalito 700-, 1980

15 Lin, T, "Probabilistic seismic hazard deaggregation of ground motion prediction models" 1-12, 2011

16 Virieux, J, "P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method" 51 : 889-899, 1986

17 Adamova, P, "Non-double-couple earthquake mechanism as an artofact of the point-source approach applied to a finite-extent focus" 100 : 447-557, 2010

18 Yunga, S, "Non double couple seismic sources, faults interaction and hypothesis of self-organized criticality" 5 : 11-15, 2005

19 Mussett, E.A, "Looking into the Earth: An introduction to geological geophysics" Cambridge University Press 470-, 2000

20 Min, D.J, "Free-surface boundary condition in finite-difference elastic wave modeling" 94 : 237-250, 2004

21 Levander, A.R, "Fourth-order finite-difference P-SV seismograms" 53 : 1425-1436, 1988

22 Kosloff, D.D, "Forward modeling by a Fourier method" 47 : 1402-1412, 1982

23 Gaspar-Escibano, J.M, "Communicating earthquake risk: Mapped parameters and cartographic representation" 11 : 359-366, 2011

24 Reynolds, C.A, "Boundary conditions for the numerical solution of wave propagation problems" 43 : 1099-1100, 1978

25 Stein, S, "An introduction to seismology, earthquakes, and earth structure" Blackwell Publishing 498-, 2003

26 Marfurt, K.J, "Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equation" 49 : 533-549, 1984

27 Clayton, R, "Absorbing boundary conditions for acoustic and elastic wave equations" 67 : 1529-1540, 1977

28 Narayan, J.P, "A fourth order accurate SH-wave staggered grid finite-difference algorithm with variable grid size and VGR-stress imaging technique" 165 : 271-294, 2008

29 Symons, N.P, "3d acoustic and elastic modeling with marmousi2" 2171-2175, 2006

30 Ho-Yong Lee, "2D time-domain acoustic-elastic coupled modeling: a cell-based finite-difference method" 한국지질과학협의회 13 (13): 407-414, 2009