북한 지하 핵실험과 그 여진의 P/S 스펙트럼 비율을 이용한 정밀 구분 연구

초록

본 논문은 2006년 이후 북한이 실시한 6건의 지하 핵실험과 그에 따른 6건의 여진을 P‑S 스펙트럼 진폭 비율(Pg/Lg)으로 구분하는 방법을 평가한다. IMS 배열국(USRK, KSRS)과 비 IMS 3‑C 관측소(SEHB, MDJ)를 활용했으며, 주파수가 높아질수록 실험과 여진 사이의 비율 차이가 크게 나타난다. 단일 관측소와 다중 관측소(마할라노비스 거리) 모두에서 99.996% 이상의 확신도로 두 집단을 구분할 수 있음을 보였다. 다만, 이 결과는 북한 실험 여진에 한정되며, 동일 지역의 자연 지진에 바로 적용하기는 어렵다.

상세 분석

이 연구는 핵실험 감시 체계에서 가장 핵심적인 문제 중 하나인 인공 지진(핵실험)과 자연 여진을 구분하는 방법론을 정량적으로 검증한다. 저자들은 P‑S 스펙트럼 진폭 비율, 즉 Pg(지표면 파)와 Lg(지구 내부 전파)의 진폭을 주파수별로 비교함으로써 두 현상의 물리적 차이를 드러냈다. 핵실험은 일반적으로 고주파 성분이 강하고, Pg/Lg 비율이 높은 반면, 여진은 저주파 성분이 상대적으로 강해 비율이 낮다. 논문에서는 0.5 Hz부터 12 Hz까지 0.5 Hz 간격으로 스펙트럼을 분석했으며, 특히 4 Hz 이상에서 두 집단 간 차이가 눈에 띄게 확대됨을 확인했다.

관측망은 IMS 배열국인 USRK와 KSRS, 그리고 비 IMS 3‑C 관측소인 SEHB(한국)와 MDJ(중국) 네 곳을 사용했다. 각 관측소에서 얻은 Pg/Lg 비율의 평균과 표준편차를 기반으로 단일 관측소 수준에서 통계적 검정을 수행했는데, KSRS에서는 여진 6건 중 어느 하나도 폭발 집단의 표본이 될 확률이 99.996% 이상이며, USRK에서는 이 확률이 더욱 높았다. 이는 단일 관측소만으로도 충분히 높은 구분 신뢰도를 제공한다는 의미다.

다중 관측소 결합을 위해 저자들은 마할라노비스 거리(D²)를 적용했다. 이는 각 관측소의 Pg/Lg 비율을 다변량 정규분포로 가정하고, 폭발군 평균과 여진군 평균 사이의 거리(제곱)를 계산하는 방법이다. 결과는 4 Hz 이상에서 D²가 100을 초과하고, 특히 USRK의 6–12 Hz 대역에서는 D²가 21 000을 넘는 등, 두 집단 간 차이가 통계적으로 거의 불가능할 정도로 크게 나타났다.

또한, 여진 탐지는 기존 IDC(International Data Centre) 자동 처리와 교차 상관 프로시저를 병행했다. IDC 자동 처리에서는 2건, 교차 상관에서는 3건, 그리고 여진 자체를 마스터 이벤트로 사용한 경우 1건을 추가로 검출했다. 이는 교차 상관 기법이 미약한 신호에도 높은 감도와 정확성을 제공함을 시사한다.

하지만 논문은 몇 가지 제한점을 명시한다. 첫째, 연구 대상이 북한 실험 여진에 국한되므로, 동일 지역의 일반 지진(예: 전단파 지진)과의 구분 능력은 검증되지 않았다. 둘째, Pg/Lg 비율은 지역 구조와 경로 효과에 민감하므로, 다른 지리적 환경에서는 동일한 임계값을 적용하기 어려울 수 있다. 셋째, 데이터 수가 비교적 적은(폭발 6건, 여진 6건) 점은 통계적 일반화에 제약을 준다.

종합적으로, 이 논문은 고주파 영역에서 P‑S 스펙트럼 비율을 활용한 정량적 구분 방법이 핵실험 감시 체계에 매우 유용함을 입증한다. 특히 마할라노비스 거리와 같은 다변량 통계 기법을 결합하면, 단일 관측소보다 더 강력한 구분 신뢰도를 확보할 수 있다. 향후 연구에서는 다양한 지진 유형과 지역을 포함한 대규모 데이터셋을 활용해 일반화 가능한 임계값을 도출하고, 실시간 자동 구분 알고리즘에 적용하는 것이 필요하다.