지진 군집 탐지 검정의 검정력 평가

초록

전 세계 대형 지진 카탈로그가 포아송 과정으로 설명될 수 있는지를 검정하는 기존 연구는 귀무가설을 기각하지 못할 경우를 주로 다루었다. 본 논문은 그 반대, 즉 비포아송적 군집성을 실제로 탐지할 능력, 즉 제2종 오류(검정력)를 평가한다. 저자는 인위적으로 군집성을 갖는 확률 과정을 설계하고, 현재 사용되는 여러 통계 검정에 대해 시뮬레이션을 수행한다. 결과는 대형 사건의 발생 빈도가 낮고 관측 기간이 짧아 검정력이 크게 감소함을 보여준다. 따라서 “포아송 가정이 옳다”는 결론은 데이터 부족으로 인한 검정력 부족을 간과한 것일 수 있다.

상세 분석

본 연구는 지진 발생을 포아송 과정으로 모델링하는 전통적 접근법에 내재된 검정력(power) 문제를 정량적으로 분석한다. 검정력은 귀무가설이 거짓일 때 이를 올바르게 기각할 확률이며, 제2종 오류율(β)과 직접적으로 연결된다. 저자는 먼저 전 세계 M≥8.5 규모 지진 카탈로그가 약 100년 정도의 관측 기간에 불과하다는 사실을 강조한다. 이러한 데이터 희소성은 사건 간 평균 간격이 수십 년에 달함을 의미하며, 통계적 검정이 충분한 자유도를 확보하기 어렵게 만든다.

문제의 핵심은 “군집성”을 어떻게 정의하고 시뮬레이션에 반영하느냐이다. 저자는 포아송 과정에 시간적 변동성을 추가한 혼합 모델을 제안한다. 구체적으로, 기본 포아송 레이트 λ₀에 더해 일정 기간(예: 10년) 동안 레이트가 λ₁>λ₀로 급증하는 클러스터 구간을 무작위로 삽입한다. 클러스터 발생 자체도 포아송 과정으로 가정해, 전체 과정은 “포아송 + 클러스터” 형태가 된다. 이렇게 생성된 인공 카탈로그는 명백히 비포아송적이지만, 실제 관측 데이터와 동일한 길이와 사건 수를 갖는다.

다음 단계에서는 기존 문헌에서 널리 사용되는 세 가지 검정을 적용한다. 첫째는 사건 간 간격이 지수분포를 따르는지를 검정하는 Kolmogorov‑Smirnov(KS) 검정, 둘째는 연속적인 시간 구간 내 사건 수가 포아송 분포를 따르는지를 검정하는 χ² 적합도 검정, 셋째는 변동 레이트 모델과 고정 레이트 모델 간의 우도비(Likelihood Ratio Test, LRT)이다. 각 검정에 대해 10⁴개의 시뮬레이션을 수행하고, p‑값 분포를 수집한다.

시뮬레이션 결과는 놀라울 정도로 낮은 검정력을 보여준다. 예를 들어, 클러스터 레이트 λ₁을 기본 레이트의 5배로 설정했음에도 KS 검정의 평균 p‑값은 0.42에 머물렀으며, 5% 유의수준에서 귀무가설을 기각한 비율은 12%에 불과했다. χ² 검정은 약간 더 민감했지만 여전히 20% 이하의 검정력만을 보였다. LRT는 가장 높은 검정력을 기록했지만, 전체 시뮬레이션 중 30% 정도만이 유의수준을 초과했다. 이는 실제 데이터가 동일한 조건이라면, 비포아송적 군집성을 가지고 있더라도 기존 검정으로는 거의 탐지되지 않을 가능성이 높다는 것을 의미한다.

또한 저자는 검정력에 영향을 미치는 주요 요인을 정량화한다. 사건 수(N)가 증가할수록 검정력은 기하급수적으로 상승하지만, N이 30~40 이하인 현재의 대형 지진 카탈로그에서는 검정력이 0.5 이하에 머문다. 관측 기간을 200년으로 연장하거나, M≥7.5 규모 이하의 사건을 포함하면 검정력이 0.8 이상으로 개선되지만, 이는 데이터 품질과 기록의 일관성 문제를 야기한다.

결론적으로, 본 논문은 “포아송 가정이 맞다”는 결론이 데이터 부족에 기인한 검정력 결핍을 간과한 결과일 수 있음을 경고한다. 향후 연구에서는 더 긴 관측 기록, 다중 지역 통합, 혹은 베이지안 사전 정보를 활용한 검정 설계가 필요함을 제안한다.