굿엔버그리히터 b값 탐지를 위한 최소 표본 크기

초록

본 연구는 완전 및 불완전한 지진 카탈로그를 모두 고려한 각도형 진도-빈도 분포(angular FMD)를 이용해, 서로 다른 b값을 구별하기 위해 필요한 최소 표본 크기를 평가한다. 최신 최대우도법과 호환되는 angular FMD를 적용하고, 몬테카를로 시뮬레이션으로 생성한 합성 카탈로그를 분석한 결과, 기존 연구에서 흔히 사용되는 표본 크기가 통계적으로 유의한 b값 차이를 검출하기에 충분하지 않음을 확인하였다.

상세 분석

이 논문은 지진학에서 가장 기본적인 통계적 법칙 중 하나인 굿엔버그-리히터(Gutenberg‑Richter, GR) 법칙의 b값을 정밀하게 추정하기 위한 최소 표본 크기를 정량적으로 규정한다는 점에서 학문적·실무적 의의를 가진다. 기존의 GR 분석에서는 완전 카탈로그(관측 가능한 최소 진도 이하의 사건이 누락되지 않은 데이터)만을 대상으로 b값을 추정했으며, 불완전 카탈로그를 포함하면 b값이 편향될 위험이 있었다. 최근 제안된 angular FMD는 관측 한계(Mc)를 명시적으로 모델링하면서도 Aki(1965)의 최대우도법(MLE)과 완전하게 호환된다. 즉, 관측 한계 이하의 사건을 포함한 전체 데이터셋을 그대로 사용해도 b값 추정에 추가적인 보정이 필요 없다는 장점이 있다.

연구진은 이 angular FMD를 기반으로, 다양한 b값(예: 0.8, 1.0, 1.2)과 관측 한계(Mc)를 갖는 합성 카탈로그를 Monte Carlo 방식으로 수천 번 생성하였다. 각 시뮬레이션에서 표본 크기를 단계적으로 감소시키면서 MLE를 적용해 b값을 추정하고, 실제 b값과의 차이가 통계적으로 유의미한지(예: 95% 신뢰구간 초과) 여부를 평가했다. 결과는 표본 크기가 작을수록 추정 오차가 급격히 증가하고, 특히 b값 차이가 0.1 수준인 경우 최소 150200개의 사건이 필요함을 보여준다. 이는 기존 연구에서 흔히 3050개의 사건만을 사용해 b값 변동성을 논의한 경우, 실제로는 통계적 검증력이 크게 부족하다는 것을 의미한다.

또한, 관측 한계가 낮을수록(즉, 더 많은 작은 진도 사건이 포함될수록) 필요한 최소 표본 크기가 다소 감소하지만, 불완전성 모델링이 부정확하면 오히려 표본 요구량이 증가한다는 역설적인 결과도 도출되었다. 이는 실제 현장 카탈로그에서 Mc를 정확히 추정하고, angular FMD를 적용하는 것이 최소 표본 크기 판단에 핵심적인 전제임을 강조한다.

이 논문의 핵심 통찰은 다음과 같다. 첫째, b값 차이를 검출하려면 단순히 “표본이 충분히 크다”는 직관적 판단이 아니라, 사전 정의된 통계적 검정력(power)과 신뢰수준을 기반으로 최소 표본 크기를 계산해야 한다. 둘째, angular FMD는 완전·불완전 데이터를 동시에 활용할 수 있는 유일한 프레임워크로, 기존의 점근적 보정 방법보다 더 정확한 b값 추정과 표본 요구량 평가를 가능하게 한다. 셋째, 실제 지진 위험 평가나 지진활동성 변동 연구에서 b값 변동을 주장할 때는 최소 150~200개의 사건을 확보하고, 관측 한계와 데이터 불완전성을 명시적으로 모델링해야 한다는 실용적 가이드라인을 제공한다.