깊은 비화산 진동의 메커니즘: 에피소드 트레머와 슬립의 새로운 해석

초록

본 논문은 사인‑고든 방정식으로 기술되는 단층 내부 마찰 현상에서 슬립 펄스(솔리톤)와 진동(포논)이 서로 다른 해로 존재함을 제시한다. 이동하는 슬립 펄스가 구조적 이질성 및 미소 어스퍼리와 상호작용하면서 공명형 진동을 유발하고, 이는 ETS 사건에서 관측되는 비화산 진동(NVT)의 원인으로 설명된다. 모델은 진동의 중심 주파수와 감쇠 특성을 단층의 탄성계수, 유효 정상응력, 관통 경도, 마찰계수와 연결시키며, 진동 전파와 역방향 이동 현상도 재현한다.

상세 분석

본 연구는 비화산 진동(NVT)을 사인‑고든(Sine‑Gordon) 방정식의 두 종류 해, 즉 솔리톤과 포논으로 해석한다. 솔리톤은 단층면을 따라 이동하는 슬립 펄스로, 비선형 파동의 특성을 지니며 에너지 손실이 적어 장거리 전파가 가능하다. 반면 포논은 작은 진동 모드, 즉 공명형 진동으로서 단층 내부의 특정 깊이 구간에서 고유 주파수를 갖는다. 슬립 펄스가 단층 내 구조적 이질성(예: 강도 변화, 물성 불균일)이나 미소 어스퍼리와 충돌하면, 에너지의 일부가 포논 모드로 전환되어 진동이 발생한다. 이 과정은 ‘진동‑슬립 상호작용’이라고 부를 수 있으며, 사인‑고든 방정식의 비선형 항이 이러한 에너지 교환을 수학적으로 기술한다.

모델은 진동의 중심 주파수 f₀가 단층의 전단 탄성계수 μ, 유효 정상응력 σₙ′, 관통 경도 H, 마찰계수 μ_f와 직접 연관된다고 제시한다. 구체적으로 f₀ ∝ √(μ/ρ)·(σₙ′/H)·(1‑μ_f) 형태로 표현되며, 여기서 ρ는 주변 암석의 밀도이다. 따라서 정상응력이 낮고 관통 경도가 높은 깊은 영역일수록 저주파 진동이 우세하고, 반대로 높은 정상응력·낮은 경도 구역에서는 고주파 성분이 강화된다. 이러한 파라미터 의존성은 관측된 ETS 진동이 주변(배경) 진동이나 지진에 의해 유발된 진동과 미세하게 다른 스펙트럼을 보이는 이유를 설명한다.

진동 전파 양상에 대해서는 두 가지 메커니즘을 제시한다. 첫째, 슬립 펄스 자체가 이동함에 따라 펄스 앞뒤로 포논이 지속적으로 방출되므로, 진동은 슬립 진행 방향(스트라이크 방향)과 거의 일치하는 속도로 전파한다. 둘째, 2차원 슬립 펄스 내부에 존재하는 ‘킥(kink)’이라 불리는 국소적인 변형이 이동하면서 포논을 빠르게 전파시킬 수 있다. 이 경우 진동 전파 속도는 슬립 펄스 자체보다 크게 증가하여 관측된 ‘급속 진동 전파’를 설명한다. 또한 슬립 펄스가 특정 지점에서 고정(pinning)되거나 파편화될 경우, 파편화된 소규모 펄스가 역방향으로 이동하면서 역방향 진동 전파가 발생한다. 이러한 현상은 ETS 사건에서 보고된 진동의 역방향 이동 패턴을 자연스럽게 재현한다.

수치 시뮬레이션 결과는 모델이 제시한 주파수‑감쇠 곡선이 실제 관측된 ETS 진동 스펙트럼과 높은 일치도를 보임을 확인한다. 특히, 진동의 중심 주파수가 2–8 Hz 범위에 머물면서 깊이별로 약간씩 변하는 경향은 모델 파라미터(σₙ′, H, μ_f)의 변동에 의해 정량적으로 설명된다. 또한, 진동 발생 빈도와 슬립 펄스 속도 사이의 상관관계도 시뮬레이션에서 재현되어, 슬립 속도가 증가할수록 진동 발생률이 상승하고, 이는 관측된 ETS 사건의 ‘진동 폭발’ 현상과 일치한다.

결론적으로, 사인‑고든 방정식을 기반으로 한 이 모델은 비화산 진동의 발생 메커니즘, 주파수 특성, 전파 및 역전파 현상을 통합적으로 설명한다. 이는 기존의 마찰‑탄성‑유체 상호작용 모델이 설명하기 어려웠던 복잡한 진동 패턴을 물리적으로 해석할 수 있는 새로운 틀을 제공한다.