정적 응력장 기하학으로 본 우츠 여진 생산법칙 해석

초록

우츠가 제시한 여진 생산법칙 K=K₀·e^{αM}을 고정된 정적 응력장을 기반으로 하는 ‘솔리드 세이시미티’ 이론으로 유도한다. 기하학적 연산을 통해 두 구간(소규모·대규모 지진)에서 서로 다른 α값(α=1.5 ln10, α=ln10)을 예측하고, 기존 카탈로그에서 α가 일정하게 보이는 것은 선택 방법의 인위적 연속 가정 때문일 수 있음을 제시한다.

상세 분석

본 논문은 우츠(1970)가 제시한 여진 생산법칙 K=K₀·exp(αM)이라는 경험적 관계를 물리적·수학적 기반 위에 놓고자 한다. 저자들은 ‘솔리드 세이시미티(Solid Seismicity)’라는 프레임워크를 도입한다. 이 이론은 지진 발생을 영구적인 정적 응력장(permanent static stress field) 위에 정의된 기하학적 영역으로 해석한다. 구체적으로, 주진(mainshock)의 파괴면을 중심으로 응력 변화가 일정 임계값을 초과하는 부피를 계산하고, 이 부피가 여진 발생 가능 영역이라고 가정한다.

정적 응력장은 지진 규모 M에 따라 파괴면의 반경 R∝10^{M/2} 로 확장된다. 여기서 응력 임계값을 초과하는 부피 V는 구형(또는 구면대칭) 가정하에 V∝R³∝10^{3M/2} 가 된다. 그러나 실제 여진은 이 부피 전체가 아니라, 응력 증가가 양의 방향으로만 작용하는 ‘양의 응력 영역(positive stress region)’에 국한된다. 이 영역은 구면의 일부(예: 상반구)로 제한될 수 있어 부피 스케일이 R³에서 R²·R₀ 형태로 변한다.

저자들은 두 가지 경우를 구분한다. 첫 번째는 작은 규모(M<M_c)에서 파괴면이 지표면에 가까워 응력 장이 반구형으로 제한되는 경우이다. 이때 여진 수 K∝V∝10^{1.5M} 이므로 α=1.5 ln10이 된다. 두 번째는 큰 규모(M>M_c)에서 파괴면이 충분히 깊어 전체 구형 응력 장이 작용하므로 K∝10^{M} 로 α=ln10이 된다. 이 두 스케일 사이에 이론적 ‘kink’(전이점)가 존재한다는 것이 핵심 주장이다.

논문은 또한 기존 여진 선택 방법—예를 들어, 일정한 시간·공간 윈도우를 적용하거나, Omori‑Utsu 법칙에 기반한 파라미터 추정—이 이러한 전이점을 매끄럽게 가리게 만든다고 지적한다. 시뮬레이션 실험에서 선택 기준을 바꾸면 α가 두 구간에서 각각 다른 값을 보이지만, 전통적인 방법을 사용하면 전체 구간에 대해 단일 α=ln10 값이 도출된다. 이는 관측 카탈로그에서 전이점이 사라진 것처럼 보이게 만든 ‘인위적 연속성’이다.

이러한 결과는 정적 응력 시뮬레이션에서 이미 보고된 α값(소규모에서 1.5 ln10, 대규모에서 ln10)과 일치한다. 그러나 실제 지진 데이터에서 전이점을 검증하려면 여진 선택 알고리즘을 재설계하고, 규모별로 별도 회귀 분석을 수행해야 한다. 저자들은 현재까지는 시뮬레이션 수준에서만 검증되었으며, 관측 데이터에서의 확증은 향후 연구 과제로 남겨두었다.

결론적으로, 본 연구는 여진 생산법칙을 경험적 관계에서 물리‑기하학적 원리로 전환시켰으며, 규모에 따른 α값 차이를 예측한다. 이는 지진 위험 평가와 여진 시나리오 구축에 새로운 이론적 토대를 제공한다.