무작위 응력과 오모리 법칙 파워리히와 인버스가우시안 혼합 모델

초록

본 논문은 지진 후 응력 변동을 브라운 운동으로 모델링한 레비와 역가우시안 분포가 오모리 법칙의 즉시 높은 여진 발생률을 설명하지 못한다는 점을 지적한다. 저자는 초점부 근처의 파워‑법칙 응력 분포와 기존 확률분포를 결합해 새로운 레비와 혼합 역가우시안 분포를 도출하고, 이들 분포가 시간 간격에 따라 오모리 법칙과 동일한 지수‑거듭제곱 형태를 나타냄을 보인다. 또한 장기적인 배경 지진율을 판구조 변형과 연결시켜 계산하고, 실제 지진 카탈로그를 통해 이론을 검증한다.

상세 분석

이 연구는 지진 후 응력의 시간적 진화를 확률 과정으로 보는 전통적 접근에 근본적인 수정점을 제시한다. 기존에 레비 분포와 역가우시안 분포가 브라운 운동에 의해 발생하는 응력 변동을 기술한다는 가정은, 시간 간격이 0에 가까워질수록 확률밀도가 0으로 수렴한다는 수학적 특성을 갖는다. 이는 오모리 법칙이 제시하는 즉각적인 높은 여진 발생률과 정면으로 충돌한다. 저자는 이 모순을 해소하기 위해, 지진 파열 영역 내에서 응력이 파워‑법칙 형태로 분포한다는 물리적 근거를 도입한다. 구체적으로 응력 크기 σ가 σ‑ψ 형태의 확률밀도를 가진다고 가정하고, 이를 레비 혹은 역가우시안 확률밀도와 곱셈적 혼합(Mixture) 형태로 결합한다. 수학적 전개 결과, 새로운 레비 혼합 분포는 순수한 파워‑법칙 형태를 띠며 지수는 -(1+ψ/2) 로 나타난다. 이는 작은 시간 간격에서 여진 발생률이 무한히 커지는 오모리 법칙의 초기 기울기와 일치한다. 반면 혼합 역가우시안 분포는 두 개의 지수 구간을 보인다. 짧은 시간에서는 레비와 동일한 -(1+ψ/2) 지수를, 보다 긴 시간에서는 -(1+ψ) 지수를 갖는다. 이는 관측된 여진률이 초기 급격히 감소한 뒤, 보다 완만한 감쇠를 보이는 현상을 자연스럽게 재현한다. 더욱 긴 시간 스케일에서는 이러한 파워‑법칙 거동이 일정한 배경 지진율로 전이된다고 가정한다. 배경율은 저자가 이전 연구에서 제시한 지진 규모 분포와 판구조 변형 속도와의 연계식을 이용해 정량화한다. 실증 분석에서는 전 세계 주요 지진 카탈로그를 대상으로 파라미터 ψ와 확산 계수 D 등을 최적화했으며, 제안된 혼합 모델이 기존 레비·역가우시안 모델보다 잔차가 현저히 작고, 오모리 법칙의 로그‑선형 적합도와도 높은 일치를 보였다. 특히 대형 지진 직후의 초단위부터 수일까지의 여진 데이터에서 새로운 모델이 초기 급증을 정확히 포착함을 확인하였다. 마지막으로 저자는 실험실 규모의 암석 변형 실험에서 응력 변동의 파워‑법칙 지수를 직접 측정함으로써 모델 파라미터를 독립적으로 검증할 수 있는 방안을 제시한다. 전체적으로 이 논문은 무작위 응력 동역학과 파워‑법칙 응력 분포의 결합을 통해 오모리 법칙의 미세구조적 메커니즘을 통계적으로 설명하는 새로운 이론적 틀을 제공한다.