How effective normal stress oscillations advance failure in fault gouge: frequency dependence, non-failure window, and the role of dilation

기술 요약: 유효 수직 응력 진동이 단층 구유(Fault Gouge)의 파쇄를 촉진하는 방식

문제 제기
지각 내 단층의 안정성은 유효 수직 응력($\sigma'_n = \sigma_n - P$, 여기서 $\sigma_n$은 총 수직 응력, $P$는 간극 유체 압력)에 의해 결정된다. 단조적인 간극 압력 변화가 유발 지진에 미치는 영향은 잘 연구되어 있으나, 입자 규모에서 주기적인 유효 수직 응력 변화에 대한 과립형 단층 구유의 반응은 여전히 규명되지 않은 부분이 많다. 자연 현상(조석, 통과하는 지진파)과 공학적 시스템(펄스형 주입)은 단층에 진동 하중을 가한다. 선행 연구들은 주기적인 수직 응력이 정적 모어-쿨롱(Mohr–Coulomb) 임계값 미만에서도 단층을 약화시킬 수 있음을 입증했으나, 넓은 주파수 스펙트럼에 걸친 거동, 안정적인 주파수 창(window)의 존재 여부, 그리고 임계 미달 파쇄를 유도하는 구체적인 입자 규모 메커니즘은 체계적으로 다뤄지지 않았다. 또한, 포화된 구유 내의 간극 압력 진동이 건조한 구유에서의 직접적인 수직 응력 진동과 동일한 반응을 생성하는지도 불분명했다.

연구 방법
저자들은 결합된 이산 요소법(DEM)과 연속체 유체 역학 솔버를 사용하여 2차원 유체 포화 과립형 단층 구유 층을 시뮬레이션하였다. 이 모델은 입자 수준의 미시 역학과 간극 유체 압력의 진화를 추적한다.

• 설정: 강성 벽 사이에 갇힌 과립 층을 모델링하였으며, 초기 상태는 밀도가 높고 과압밀된 상태(공극률 $\phi_0 \approx 0.185$)이다. 이 층에는 정적 파쇄 임계값 미만의 일정한 전단 응력($\tau = 0.2\sigma_0$)이 가해진다.
• 강제 프로토콜: 0.5에서 $10^4$ Hz 사이의 주파수 범위에서 두 가지 상호 보완적인 프로토콜을 테스트하였다. 이때 피크 하중은 준정적 모어-쿨롱 파쇄 임계값보다 엄격하게 낮게 유지되었다:
  1. 건조 대조군 (Dry Control): 건조한 층에 수직 응력($\sigma_n$)을 직접 주기적으로 진동시킴.
  2. 습윤 프로토콜 (Wet Protocol): 일정한 $\sigma_n$ 하에서 간극 압력($P$)을 주기적으로 진동시켜 테르자기(Terzaghi)의 유효 응력 등가 원리를 테스트함.
• 매개변수: 시뮬레이션에는 실제적인 입자 특성(영률, 마찰 계수)과 높은 투과율($k = 10^{-11}$ m$^2$)을 사용하여, 입자의 운동량 확산 속도에 비해 간극 압력 확산이 빠르게 일어나도록 하였다.

주요 기여 및 결과
본 연구는 단층 안정성에 대한 비단조적 주파수 의존성을 식별하여, 세 가지 뚜렷한 파쇄 모드와 하나의 견고한 "비파쇄 창(non-failure window)"을 밝혀냈다.

1. 세 가지 파쇄 모드와 비파쇄 창:

   • 모드 I (저주파, $f \sim 0.5\text{--}1$ Hz): 주기적 파쇄 및 정지. 시스템이 하중 주기와 동기화되어 미끄러지고 멈춘다.
   • 모드 II (중저주파, $f \sim 5\text{--}10$ Hz): 정지 없는 주기적 슬라이딩. 층이 지속적으로 미끄러지되 주기에 의해 조절되며, 완전히 재재킹(re-jamming)되지 않는다.
   • 비파쇄 창 (Non-Failure Window, $f \sim 30\text{--}200$ Hz): 시스템이 작은 진폭의 주기적 크리프(creep)만을 보이는 뚜렷한 영역. 하위 임계 하중에도 불구하고 거시적인 파쇄는 발생하지 않는다.
   • 모드 III (고주파, $f \ge 500$ Hz): 지속적인 슬라이딩으로 이어지는 파쇄. 시스템이 개별 하중 피크와 무관하게 일정 기간의 안정기 이후 지속적인 슬라이딩으로 전이된다.

2. 건조 및 습윤 프로토콜의 등가성:
   두 프로토콜 모두 동일한 주파수에서 동일한 영역 구조(세 가지 모드 및 비파쇄 창)를 생성하였다. 이는 간극 압력 진동과 직접적인 수직 응력 진동이 테르자기의 원리에 따라 등가적인 강제력을 제공함을 확인해 준다. 주요 차이점은 팽창 경화(dilatant hardening)(팽창 중 발생하는 일시적인 간극 압력 저하)로 인한 습윤 사례에서의 시간 지연이었으나, 근본적인 안정성 경계는 변하지 않았다.

3. 임계 미달 쿨롱 파쇄(Sub-Coulomb Failure)의 메커니즘:
   논문은 파쇄의 원인을 전단 유도 팽창에 의한 마찰 계수 $\mu(\phi)$의 점진적인 저하로 돌린다. 메커니즘은 주파수에 따라 달라진다:

   • 저주파 (전단 주도 래칫 효과, Shear-Driven Ratcheting): 진동 주기 $T$가 입자 운동량 확산 시간($\zeta_v$)보다 훨씬 길 때, 각 주기는 전단 유도 미세 미끄러짐이 비가역적 팽창을 일으키기에 충분한 시간을 제공한다. 층이 다음 주기 전에 완전히 재압착되지 않으므로, 공극률이 주기마다 축적되어 임계 파쇄 공극률에 도달할 때까지 "래칫(ratchet)" 효과가 발생한다.
   • 고주파 (진동 주도/음향 유동화, Vibration-Driven/Acoustic Fluidization): $T \ll \zeta_v$일 때, 전단 주도 래칫 효과는 억제된다. 대신, 빠른 진동이 증폭된 침투력(wet case의 경우), 응력 구배 및 관성력을 생성한다. 이러한 동적 힘들이 입자 집합체를 교란하여, 전단 응력과 무관하게 팽창과 재배열(음향 유동화)을 유도하고 결국 파쇄로 이어진다.
   • 비파쇄 갭 (The Non-Failure Gap): 중간 영역은 진동이 전단 주도 래칫을 일으키기에는 너무 빠르고(팽창 축적 시간이 부족), 진동 주도 힘이 입자를 교란할 만큼 커지기에는 너무 느린(진동력이 충분히 발달하지 못함) 영역을 나타낸다.

의의 및 주장
저자들은 이러한 결과가 다음 사항에 대한 입자 규모의 기계적 근거를 제공한다고 주장한다:

• 임계 미달 파쇄 (Sub-Coulomb Failure): 정적 강도를 초과하지 않는 진동 하중 하에서도 단층이 어떻게 파쇄될 수 있는지를 설명한다.
• 동적 지진 트리거링 (Dynamic Earthquake Triggering): 지진파에 의한 동적 트리거링의 효율성이 파동 주파수와 단층 구유의 고유 시간 척도($\zeta_v$) 사이의 상호작용에 달려 있음을 시사한다. 비파쇄 창의 존재는 모든 주파수가 슬립을 유발하는 데 동일하게 효과적인 것은 아님을 의미한다.
• 주기적 주입 프로토콜: 지열 및 이산화탄소 격리 프로젝트에서 유체 주입을 최적화하기 위한 이론적 틀을 제공한다. 연구 결과는 진동 주파수를 "비파쇄 창" 내에서 선택함으로써 유발 지진을 완화할 수 있는 반면, 다른 주파수는 역설적으로 하위 임계 압력에서 파쇄를 촉진할 수 있음을 시사한다.

본 연구는 이러한 영역의 경계가 보편적인 것이 아니라 시스템 크기(특히 $\zeta_v \propto h^2$인 구유 두께)에 따라 스케일링된다는 점을 강조한다. 따라서 실험실 결과는 단층대가 훨씬 더 두꺼운 현장 조건으로 스케일링되어야 하며, 이에 따라 주파수 창이 낮은 값으로 이동하게 된다. 본 연구의 결과는 연구된 고투과율 영역에 국한된 것이며, 저투과율 구유는 간극 압력 확산이 더 느리기 때문에 다른 거동을 보일 수 있다.