• 지구물리와물리탐사
• /
• 제10권3호
• /
• pp.191-202
• /
• 2007

한반도 남부지역의 세부적인 지진파 감쇠특성 규명을 위해, 기존에는 적용이 불가능하였던 Q 토모그래피 역산을 위한 사전 수치검증 연구를 수행하였다. 특히 강지진동모사를 위해 일반적으로 사용되고 있는 추계학적 점지진원 지진동 모델(stochastic point-source ground-motion model; Boore, 2003)에서 사용되는 Q 값에 대한 2차원(2D; 2 Dimensional) 토모그래피 역산을 시도함으로써 역산 결과가 강지진동모사에 직접 활용될 수 있도록 하였다. 수치검증 방법으로는 Q 토모그래피 checkerboard 시험방법이 사용되었는데, 이를 위해 광역 단일 Q 모델의 추계학적 지진동모델 파라미터 역산결과의 지진원과 부지효과 모델 파라미터 값을 이용해서 관측자료와 지진규모-거리-주파수-오차 분포가 동일한 스펙트럼 합성자료를 생성하였다. 수치검증을 위한 Q 블록 격자의 총 개수는 75개(내륙지역=69개(약 $35{\times}44km^2$의 격자크기); 해양지역=6개)로 설정하였으며, $Q_0f^{\eta}$ 함수형태의 Q 블록 값은 $Q_0$=100, 500, ${\eta}=0.0{\sim}1.0$의 분포를 갖도록 하고, 파선의 깊이는 별도로 고려하지 않았다. 스펙트럼 합성자료 생성에 이용된 모델파라미터의 정해와 모델파라미터의 역산결과를 비교하기 위한 checkerboard 수치검증은 3단계에 걸쳐 수행되었는데, 1단계는 블록별 Q의 초기값 추정 단계이며, 2단계는 관측소별 부지증폭함수를 추정하는 단계, 마지막 3단계에서는 최종적인 Q를 도출하는 단계이다. 관측소별 부지증폭함수의 초기 추정값으로는 기 분류된 관측소 등급에 대한 평균 부지증폭함수(연관희, 서정희, 2007)가 사용되었으며, 3단계의 checkerboard 수치검증 결과 최종적으로 추정된 부지효과 모델에는 오차가 발생하였으나 블록별 Q의 정해는 만족할 정도로 추정할 수 있었다.

• 지구물리와물리탐사
• /
• 제9권3호
• /
• pp.241-249
• /
• 2006

한국의 지질 환경은 암석 분포가 매우 다양하고 복잡한 구조 활동의 영향을 받아 지하매질의 이방성 특성이 국부적으로 심하게 변화한다. 기존의 이방성 주시 모델링의 경우 지질 모델을 2차원으로 단순화시킴으로써 이러한 복잡한 지질 환경을 제대로 고려할 수 없었다. 또한 약 이방성 가정을 사용하여 실제로 나타날 수 있는 지하 매질의 심각한 이방성 영향을 주시 모델링에서 고려할 수 없었다. 이에 이 연구에서는 보다 실제적이고 복잡한 3차원 횡등방성 매질(transversely isotropic media)에서 q-P파의 초동 주시 양상을 모사할 수 있는 주시 모델링 알고리듬을 개발하였다. 이 알고리듬에서는 2차원 비선형 주시 내삽(2D nonlinear traveltime interpolation) 기법과 주시의 3차원 격자 채움법(mapping)을 이용한 직접 전파법(direct calculation)을 통해 급격한 물성의 변화에도 주시 계산이 가능하도록 하였다. 또한, 최소 주시 계산과정에서 수치 미분을 통한 최대 경사법(steepest descent method)을 사용하여 약 이방성 가정을 극복하였다. 개발된 알고리듬은 해석해와 비교하여 그 타당성을 검증하였고 3차원 2층구조에 대한 주시 계산을 수행하여 물성이 급격히 변화하는 모델에 대해서도 안정적으로 주시 계산이 이루어짐을 확인하였다. 이 연구에서 개발한 3차원 주시 모델링 알고리듬은 향후 구조보정이나 토모그래피 알고리듬 개발에 사용될 수 있을 것으로 기대한다.