오모리 법칙 150주년 지진 여진의 미분 방정식 재조명

초록

본 논문은 오모리의 여진 법칙을 미분 방정식 형태로 재구성하고, 주진 후 암석이 ‘냉각’되는 비정상성을 반영한다. 이를 통해 전통적인 옴리‑우츠 법칙의 한계를 보완하고, 여진 시계열의 물리적 해석 가능성을 제시한다.

상세 분석

오모리 법칙은 여진 발생 빈도가 시간에 따라 t⁻¹ 형태로 감소한다는 경험적 관계로, 19세기 말부터 지진학의 기본 모형으로 활용돼 왔다. 그러나 기존 식은 암석 물성의 시간 의존성을 무시하고, 파라미터 p 와 c 를 단순 상수로 가정한다는 근본적인 한계가 있다. 저자들은 이 문제를 해결하기 위해 여진 발생률 n(t) 에 대한 미분 방정식
dn/dt = ‑α n / (t + τ) + β · f(t)
형태를 제시한다. 여기서 α 는 초기 여진 강도를, τ 는 초기 지연 시간을, β 는 암석 ‘냉각’ 효과를 기술하는 함수 f(t) 를 통해 비정상성을 도입한다. f(t) 는 온도 혹은 응력 완화에 비례하는 지수 감쇠 함수로 가정되며, 이는 실험실에서 관찰된 암석의 점탄성 회복 현상을 반영한다.

수식 유도 과정에서는 에너지 방출률과 응력 재분배 메커니즘을 결합해, 여진 발생률이 응력 감소와 동시에 물성 변화에 의해 추가 감쇠된다는 물리적 해석을 제공한다. 파라미터 추정은 비선형 최소제곱법과 베이지안 사전분포를 이용해 실제 지진 기록(예: 1995년 쿠마모토, 2011년 동일본)과 비교 검증하였다. 결과는 전통적인 옴리‑우츠 모델보다 R² 값이 평균 15 % 향상되고, 특히 초기 t < 1 일에서 과대예측을 크게 감소시켰다.

또한, 저자는 ‘암석 냉각’ 개념을 지진 발생 후 미세균열 네트워크의 재구성 속도와 연계시켜, f(t) 의 시간 상수 τ_c 가 지역 지질조건에 따라 변동함을 실증하였다. 이는 여진 위험도 평가 시 지역 맞춤형 파라미터 설정이 필요함을 시사한다. 한계점으로는 f(t) 의 형태가 단순 지수 감쇠에 국한되어 있어, 복합적인 열‑기계‑화학 상호작용을 완전히 포착하지 못한다는 점이다. 향후 연구에서는 다중 스케일 모델링과 실시간 관측 데이터를 결합해 f(t) 를 동적으로 업데이트하는 방안을 제시한다.