정규 지진과 느린 미끄럼을 아우르는 단순 1차원 모델

초록

본 논문은 초기 전단‑정규 응력 분포와 마찰 법칙을 보다 현실적으로 개선한 1차원 모델을 제시한다. 모델은 단층 강도, 전단‑정규 응력 비율 및 그 기울기에 따라 정규 지진과 느린 미끄럼 사건(SSE)을 구분하고, 파열 속도·슬립 속도·응력 강하 범위를 정량화한다. 또한 파열 모드가 균열형인지 펄스형인지도 응력 이질성 형태에 따라 결정된다는 점을 밝힌다. 중요한 결과는 불안정 발생 후 파열 역학이 전통적인 rate‑state 파라미터 a·b 비에 의존하지 않으며, 기존 rate‑state 기반 모델의 한계를 제시한다는 것이다.

상세 분석

이 연구는 기존의 지진‑슬립 모델이 지나치게 이상화된 초기 응력 조건과 단순한 속도‑의존 마찰식에 의존한다는 점을 비판하고, 보다 물리적으로 타당한 초기 전단‑정규 응력 분포와 복합 마찰 법칙을 도입한다는 점에서 의미가 크다. 모델은 1차원 연속체를 가정하고, 전단 응력 τ와 정규 응력 σ의 비율 r = τ/σ를 핵심 변수로 삼는다. r과 그 공간적 기울기 ∂r/∂x가 일정 임계값을 초과하면 전단‑정규 응력 비가 급격히 변하면서 파열이 시작된다. 여기서 파열 유형은 두 가지로 구분된다. 첫째, r과 ∇r가 충분히 크면 파열은 균열형(crack‑like)으로 전파되며, 파열 전선은 뒤따르는 회복 구역 없이 지속적으로 전진한다. 둘째, r·∇r가 낮은 경우 파열은 자체 치유 펄스(self‑healing pulse) 형태로 전파되며, 파열 전선 뒤에 슬립이 멈추는 구역이 형성된다. 특히, 응력 이질성이 계단형(step‑like)이라면 파열은 무조건 균열형이 되며, 국소적인 이질성에서는 위의 두 모드가 전환될 수 있다.

정규 지진과 느린 미끄럼 사건을 구분하는 기준은 단층 강도와 r, ∇r의 조합이다. 모델은 정규 지진의 경우 전단 파열 속도가 전단파 속도 Vs의 0.1 ~ 1배, 경우에 따라 초과(초음속)까지 도달할 수 있음을 보여준다. 슬립 속도는 보통 1 m/s 수준이며, 응력 강하는 초기 전단 응력의 1 ~ 10 %에 해당한다. 반면 SSE는 슬립 속도가 연간 수 cm에서 0.1 m/s까지 광범위하게 변하고, 응력 강하는 0.1 ~ 10 % 수준으로 정규 지진보다 낮지만, 여전히 비선형적인 응답을 보인다. 이러한 정량적 범위는 관측된 ETS(episodic tremor and slip) 현상과도 일치한다.

가장 파격적인 결과는 파열 후 역학이 전통적인 rate‑state 마찰 파라미터 a와 b의 비율에 무관하다는 점이다. 즉, 마찰이 속도 약화(velocity‑weakening)든 속도 강화(velocity‑strengthening)든 파열 전개 속도·형태는 동일하게 결정된다. 이는 기존 rate‑state 모델이 파열 전개를 설명하는 데 근본적인 한계가 있음을 시사한다. 따라서 향후 모델링에서는 a·b 비보다 전단‑정규 응력 비와 그 기울기에 초점을 맞추어야 할 필요가 있다.

이 모델은 단순함에도 불구하고 복잡한 지진‑슬립 현상을 포괄적으로 설명한다는 점에서 학계와 실무 현장 모두에 큰 영향을 미칠 것으로 기대된다. 특히, 초기 응력 이질성의 측정과 단층 강도 추정이 가능한 경우, 이 모델을 통해 지진 위험도와 슬립 이벤트의 발생 가능성을 정량적으로 예측할 수 있다.