지진원천 파라미터 추정을 위한 확률적 스펙트럼 반전 방법 및 이탈리아 중앙지진 사례 적용

초록

본 논문은 변위 진폭 스펙트럼을 이용해 지진의 모멘트, 코너 주파수, 고주파 감쇠 지수 및 품질인자(Q)를 동시에 추정하는 확률론적 프레임워크를 제시한다. 베이시안 접근으로 사후 확률밀도함수(PDF)를 정의하고, 전역 탐색‑지역 최적화를 결합한 Basin‑Hopping 기법으로 최대 사후 확률점을 찾는다. 합성시험을 통해 신호‑대‑잡음비와 대역폭이 파라미터 상관관계에 미치는 영향을 검증했으며, 2016‑2017년 이탈리아 중앙지진 연속에 적용해 평균 전단강하 2.1 MPa, 주요 사건은 평균보다 큰 전단강하를 보였음을 보고한다.

상세 분석

이 연구는 지진원천 파라미터 추정 문제를 “데이터와 모델 사이의 정보 결합”이라는 베이시안 관점에서 재정의한다. 먼저 변위 스펙트럼을 일반화된 Brune 모델(모멘트 M, 코너 주파수 f_c, 고주파 감쇠 지수 γ)과 거리·전파 손실을 포함한 Green 함수(품질인자 Q)로 표현하고, 로그 변환을 통해 선형화된 관측 연산자를 도출한다. 사전 확률분포는 균일하게 가정하고, 관측오차와 모델오차를 각각 정규분포로 모델링함으로써 전체 사후 PDF는 다변량 정규분포 형태의 우도함수와 결합된다. 핵심은 이 사후 PDF의 최대값을 찾는 최적화 단계이다. 저자는 전역 탐색을 위한 마코프 체인 기반의 Basin‑Hopping 알고리즘을 선택했는데, 이는 전역적인 탐색 능력과 지역적인 quasi‑Newton 최적화의 빠른 수렴을 동시에 제공한다. 탐색 과정에서 세 개의 파라미터(온도, 스텝 크기, 재시작 횟수)만을 조정해 문제 특성에 맞게 튜닝할 수 있다.

최적점 탐색 후에는 최대 가능도 주변에서 PDF를 샘플링해 주변 평균, 분산, 상관계수를 계산한다. 여기서 저자는 마진 PDF가 가우시안 형태와 얼마나 일치하는지를 평가하는 품질 지표와, 신호‑대‑잡음비(SNR) 기반의 주파수 대역 선택 기준을 도입해 불안정하거나 제약이 약한 해를 자동으로 배제한다. 합성시험에서는 SNR을 10 이하로 낮추면 코너 주파수와 γ 사이의 상관관계가 급격히 증가하고, 최소 코너 주파수 주변 ±1 옥타브 이상의 대역폭이 확보돼야 파라미터가 독립적으로 회복된다는 것을 확인했다.

실제 데이터 적용에서는 2016‑2017년 이탈리아 중앙지진 연속의 30여 건 사건을 분석했다. 모든 사건에 대해 모멘트‑코너 주파수 관계가 Brune 스케일링을 따르는 것으로 나타났으며, 평균 전단강하(Δσ)는 2.1 ± 0.3 MPa였다. 그러나 Mw ≥ 5.5 수준의 주요 사건들은 3 ~ 4 MPa 수준의 높은 전단강하를 보였으며, 이는 전형적인 규모‑전단강하 스케일링과의 편차를 시사한다. 또한, 추정된 Q값은 150 ~ 300 범위에 머물렀고, 지진 효율성(에너지 대비 방출된 파동 에너지 비율)은 낮은 편으로, 이는 동적 효과에 의한 스트레스 오버슈트 혹은 마찰 효율성 증가와 연관될 가능성을 제시한다.

이와 같이 확률적 프레임워크는 파라미터 간 상관관계를 정량적으로 파악하고, 불확실성을 체계적으로 전파함으로써 기존 단일 최적화 방식보다 더 신뢰성 있는 지진원천 해석을 제공한다. 특히 대규모 자동화 분석에 적합한 사전‑사후 확률 구조와 품질 기준은 향후 지진 카탈로그 구축 및 지진 위험 평가에 직접적인 활용 가치를 가진다.