지진 시계열의 난류 유사 행동

초록

본 연구는 지진 전후의 수직 지진속도 시계열을 난류 흐름에서 사용되는 계층적 확률 모델로 분석한다. 평균을 제거한 증분을 detrend한 뒤 확률밀도함수(PDF)의 형태 변화를 추적했으며, 중·대형 지진 발생 5~10시간 전부터 Gaussian에서 비Gaussian(꼬리가 두꺼운) 형태로 급격히 전이함을 확인했다. 이 전이는 장거리 상관성 감소와 동시에 나타나며, 향후 지진 전조 탐지에 활용될 수 있다.

상세 분석

본 논문은 지진계측 데이터, 특히 지표면 수직 속도(Vz) 시계열을 난류 현상에서 도입된 계층적 확률 모델에 적용함으로써 기존의 선형 혹은 단순 통계 분석을 넘어서는 새로운 통찰을 제공한다. 먼저, 원시 시계열에 대해 일정 구간(Δt)마다 2차 다항식으로 추세를 제거(detrending)하고, 그 차분 ΔV(t,Δt)=V(t+Δt)−V(t)−trend를 구한다. 이러한 증분은 시간 스케일에 따라 서로 다른 통계적 특성을 보이며, 특히 큰 스케일에서는 장거리 상관성을, 작은 스케일에서는 급격한 변동성을 반영한다.

PDF 분석에서는 증분의 분포를 로그정규(LN) 형태의 곱셈적 카스케이드 모델로 설명한다. 즉, ΔV(t,Δt)=W(t)·ε(t)이며, W(t)=exp(σ·G(t)) (G는 평균 0, 분산 1인 가우시안)이고 ε(t)는 단위 분산의 가우시안 백색 잡음이다. σ는 카스케이드 강도를 나타내는 파라미터로, σ=0이면 순수 가우시안, σ>0이면 꼬리가 두꺼운 비가우시안이 된다. 논문에서는 σ를 이동 윈도우(예: 6시간) 내에서 추정했으며, 지진 전후에 σ가 급격히 증가하는 현상을 발견했다.

특히, 중·대형 지진(예: Mw≥6.0) 발생 510시간 전부터 σ가 0.20.3 수준으로 상승하면서 PDF가 가우시안에서 뚜렷한 레비-스틸러 형태로 변한다. 동시에, Detrended Fluctuation Analysis(DFA)를 적용한 결과, 지표면 속도 증분의 장거리 상관 지수 α가 0.9 수준에서 0.6 이하로 감소한다. 이는 시스템이 기존의 장거리 연관성을 잃고, 보다 무작위적인 변동성으로 전이함을 의미한다.

이러한 두 현상—PDF 꼬리 두꺼워짐과 상관 지수 감소—가 동시에 나타나는 시점을 ‘전이 시점’으로 정의하고, 이를 지진 전조로 활용한다는 점이 논문의 핵심 기여이다. 또한, 난류 흐름에서 관측되는 에너지 카스케이드와 유사하게, 지진 전조 단계에서는 지구 내부의 응력 전달이 다중 스케일에서 비선형적으로 재분배되는 것으로 해석한다.

실험적으로는 일본과 캘리포니아 지역의 20여 건의 지진 데이터를 사용했으며, 각 사건에 대해 σ와 α의 시간적 변화를 비교하였다. 결과는 일관되게 전이 시점이 지진 발생 전 5~10시간 내에 존재함을 보여준다. 또한, 소규모 지진(Mw<5)에서는 σ 변화가 미미하거나 전이 시점이 불분명해, 이 방법이 중·대형 지진에 특히 유효함을 시사한다.

이와 같은 접근은 기존의 지진 전조 탐지(예: 전자기 방출, 가스 배출 등)와는 독립적인 물리적 메커니즘을 제공하므로, 다중 센서 융합 시스템에 통합될 경우 조기 경보 정확도를 크게 향상시킬 가능성이 있다. 다만, σ 추정에 사용되는 윈도우 길이와 데이터 품질에 민감하므로, 실시간 적용을 위해서는 최적화된 신호 처리 파이프라인과 노이즈 억제 기법이 필요하다.