지진 현재상태 추정에서 자연시간 예측까지 DIY 방법

초록

본 논문은 기존의 지진 현재상태(nowcasting) 기법인 “카운팅 방법”에 ROC 분석을 결합해, 큰 지진이 발생하기 전까지 관측된 작은 지진의 수를 자연시간으로 활용하여 확률을 직접 계산한다. 확률은 ROC 곡선의 양성예측값(PPV)으로 정의되며, 별도의 확률 모델을 가정하지 않는다. 로스앤젤레스 지역의 1994년 노스리지를 지진 사례에 적용해, 작은 지진이 누적될수록 PPV가 처음에는 감소하고 이후 급격히 상승하는 패턴을 확인하였다.

상세 요약

이 연구는 지진 위험을 실시간으로 평가하는 ‘nowcasting’ 접근법을 한 단계 발전시킨다. 핵심은 Gutenberg‑Richter(GR) 관계, 즉 대형 지진(M > M_T) 하나당 평균적으로 N_GR개의 소형 지진(M > M_S)이 발생한다는 통계적 사실을 이용하는 것이다. 기존 카운팅 방법은 마지막 대형 지진 이후 발생한 소형 지진의 누적 수 n(t)를 단순히 비례 관계에 대입해 위험도를 추정했지만, 확률 모델을 가정해야 하는 한계가 있었다. 본 논문은 이를 극복하기 위해 ROC(수신자 동작 특성) 분석을 도입한다. 관측된 n값을 임계치로 삼아 ‘대형 지진 발생 여부’를 이진 분류 문제로 전환하고, 다양한 임계값에 대한 민감도와 특이도를 계산해 ROC 곡선을 만든다. 그 곡선 아래 면적(AUC) 자체는 모델 성능을 나타내지만, 여기서는 곡선상의 각 점에 대응하는 양성예측값(PPV)을 직접 확률로 해석한다. PPV는 “현재까지 관측된 소형 지진 수가 특정 값일 때, 다음 대형 지진이 곧 발생할 확률”을 의미한다.

자연시간 개념은 물리적 시간 대신 ‘소형 지진의 발생 횟수’를 시간 척도로 삼아, 지역별 지진 활동률 차이를 자연스럽게 보정한다. 이는 특히 관측망의 감도 차이나 지역별 누적 스트레스 차이를 반영하는 데 유리하다. 실험에서는 로스앤젤레스 주변을 연구구역으로 설정하고, M ≥ 6.7인 1994년 노스리지를 기준 사건 이후의 M ≥ 2.5 소형 지진 데이터를 사용했다. 결과는 흥미로운 두 단계 패턴을 보인다. 초기에는 소형 지진이 증가함에 따라 PPV가 감소하는데, 이는 대형 지진이 발생한 직후 잔류 응력이 감소해 일시적으로 큰 지진 위험이 낮아지는 현상을 반영한다. 이후 일정 수치에 도달하면 PPV가 다시 상승하기 시작하고, 어느 시점에서 급격히 가파르게 증가한다. 이 전환점은 누적된 소형 지진이 충분히 스트레스를 재축적했음을 시사한다.

이 방법의 장점은 (1) 확률 모델을 사전에 가정하지 않아 데이터에 대한 과잉적합 위험이 적다, (2) ROC 기반 PPV가 직관적인 위험도 지표를 제공한다, (3) 자연시간을 사용해 지역별 비교가 용이하다는 점이다. 반면 한계점도 존재한다. 첫째, GR 관계의 기울기(b값)와 완전성 임계값(M_S)이 정확히 설정되지 않으면 n(t) 자체가 편향될 수 있다. 둘째, ROC 곡선을 신뢰성 있게 만들기 위해서는 충분히 많은 대형 지진 사례가 필요하지만, 실제로는 데이터가 희소하다. 셋째, PPV는 사전 확률(대형 지진 발생 빈도)에 크게 의존하므로, 장기적인 평균 발생률이 변할 경우 재조정이 필요하다. 마지막으로, 자연시간이 물리적 시간과 완전히 일치하지 않기 때문에, 인프라 대비 위험 관리(예: 긴급 대응 계획)에는 보완적인 물리시간 분석이 병행되어야 한다.

전반적으로 이 논문은 지진 위험을 실시간으로 정량화하는 새로운 DIY 프레임워크를 제시하며, 특히 데이터가 풍부한 지역에서 빠르게 적용 가능하다는 실용성을 강조한다. 향후 연구에서는 다중 규모(다중 M_S) 카운팅, 베이지안 업데이트와 결합한 동적 ROC, 그리고 다른 지진대(예: 일본, 터키)에서의 검증을 통해 방법론을 확장할 여지가 크다.