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

技術要約：有効法線応力の振動がいかに断層ガウジの破壊を促進するか

問題提起
地殻内の断層の安定性は、有効法線応力（$\sigma'_n = \sigma_n - P$、ここで $\sigma_n$ は全法線応力、$P$ は間隙流体圧）によって支配されている。間隙圧の単調な変化が誘発地震に及ぼす影響については十分に研究されているが、粒状の断層ガウジに対する「周期的な」有効法線応力の変化への応答については、粒子スケールでの制約が不十分である。自然現象（潮汐、通過する地震波）や工学的システム（パルス注入）は、断層に振動荷重を課す。先行研究では、周期的な法線応力が静的なモール・クーロンの閾値を下回るレベルで断層を弱体化させることが確立されているが、広い周波数スペクトルにわたる挙動、安定した周波数ウィンドウの存在、およびサブクリティカル（閾値未満）な破壊を駆動する具体的な粒子スケールのメカニズムについては、体系的に解明されていなかった。さらに、飽和したガウジにおける間隙圧の振動が、乾燥したガウジにおける直接的な法線応力の振動と同じ応答を生み出すのかについても不明であった。

手法
著者らは、離散要素法（DEM）と連続体流体力学ソルバーを結合させ、二次元の流体飽和粒状断層ガウジ層をシミュレートした。このモデルは、粒子レベルの微細力学と間隙流体圧の進化を追跡する。

• セットアップ: 硬質な壁の間に閉じ込められた粒状層（空隙率 $\phi_0 \approx 0.185$ の高密度・過圧密状態）を用い、静的な破壊閾値を下回る一定の剪断応力（$\tau = 0.2\sigma_0$）を負荷した。
• 強制プロトコル: 0.5から$10^4$ Hzの周波数範囲において、2つの補完的なプロトコルをテストした。これらは、ピーク荷重を静的なモール・クーロン破壊閾値未満に厳格に維持している。
  1. 乾燥コントロール: 乾燥した層に対する法線応力（$\sigma_n$）の直接的な周期振動。
  2. ウェット・プロトコル: 一定の $\sigma_n$ を保持したまま、流体飽和層の孔隙圧（$P$）を周期的に振動させるもの。これにより、テルツァーギの有効応力の原理の等価性を検証した。
• パラメータ: シミュレーションには現実的な粒子特性（ヤング率、摩擦係数）を用い、粒子の運動量拡散に対して間隙圧の拡散が迅速に進むよう、高い浸透率（$k = 10^{-11}$ m$^2$）を設定した。

主要な貢献と結果
本研究は、断層の安定性における非単調な周波数依存性を特定し、3つの異なる破壊モードと、一つの堅牢な「非破壊ウィンドウ」を明らかにした。

1. 3つの破壊モードと非破壊ウィンドウ:

   • モード I (低周波数, $f \sim 0.5\text{--}1$ Hz): 周期的な破壊と停止。システムは負荷サイクルと同期して滑り、そして停止する。
   • モード II (中低周波数, $f \sim 5\text{--}10$ Hz): 停止を伴わない周期的な滑り。層はサイクルによって変調されながら連続的に滑り続け、完全に再ジャミング（固結）することはない。
   • 非破壊ウィンドウ ($f \sim 30\text{--}200$ Hz): システムが小さな振幅の周期的なクリープのみを示す、明確な領域。サブクリティカルな荷重下であるにもかかわらず、マクロな破壊は発生しない。
   • モード III (高周波数, $f \ge 500$ Hz): 連続的な滑りに至る破壊。システムは、個々の荷重ピークとは無関係に、一見安定しているように見える期間を経て、持続的な滑りへと移行する。

2. 乾燥およびウェット・プロトコルの等価性:
   両方のプロトコルは、同じ周波数において同一のレジーム構造（3つのモードと非破壊ウィンドウ）を生じさせた。これは、間隙圧の振動と直接的な法線応力の振動が、テルツァーギの原理を介して等価な強制力として作用することを裏付けている。主な違いは、膨張硬化（dilatant hardening）（膨張に伴う一時的な間隙圧の低下）による時間遅れであったが、根本的な安定性の境界は変わらなかった。

3. サブ・クーロン破壊のメカニズム:
   論文では、剪断誘起の膨張によって駆動される摩擦係数 $\mu(\phi)$ の漸進的な劣化に破壊の原因を求めている。メカニズムは周波数によって異なる。

   • 低周波数 (剪断駆動型ラチェット現象): 振動周期 $T$ が粒子の運動量拡散時間スケール（$\zeta_v$）よりもはるかに長い場合、各サイクルにおいて、剪断駆動のマイクロスリップが不可逆的な膨張を引き起こすのに十分な時間が与えられる。層は次のサイクルが来る前に完全に再圧密されないため、空隙率がサイクルごとに蓄積していく「ラチェット」効果が生じ、臨界的な破壊空隙率に達するまで進行する。
   • 高周波数 (振動駆動型／音響流動化): $T \ll \zeta_v$ の場合、剪断駆動のラチェット現象は抑制される。代わりに、急速な振動が（ウェットの場合に）増幅された浸透力、応力勾配、および慣性力を生成する。これらの動的な力が粒子集合体を撹乱し、剪断応力とは無関係に、膨張と再配列（音響流動化）を誘発して、最終的に破壊に至る。
   • 非破壊ギャップ: この中間的な窓は、剪断駆動のラチェットには早すぎ（膨張の蓄積に不十分）、かつ振動駆動の力が粒子を撹乱するほど大きくならないほど遅い、というレジームを表している。

意義と主張
著者らは、これらの結果が以下の現象に対する粒子的メカニズムの基礎を提供するとしている。

• サブ・クーロン破壊: 荷重が静的な強度閾値を決して超えない振動荷重の下で、いかにして断層が破壊し得るかを説明する。
• 動的な地震誘発: 地震波による動的なトリガーの効率は、波の周波数と断層ガウジ固有の時間スケール（$\zeta_v$）との相互作用に依存することを示唆している。非破壊ウィンドウの存在は、すべての周波数が滑りを誘発するのに等しく効果的であるわけではないことを意味する。
• 周期的な注入プロトコル: 地熱発電やCO2回収・貯留プロジェクトにおける流体注入を最適化するための理論的枠組みを提供する。本結果は、「非破壊ウィンドウ」内の振動周波数を選択することが、誘発地震を抑制する手段となり得る一方で、他の周波数ではパラドキシカルに（逆説的に）サブクリティカルな圧力下での破壊を促進する可能性があることを示唆している。

本研究は、これらのレジームの境界が普遍的なものではなく、システムサイズ（具体的には、$\zeta_v \propto h^2$ であることからガウジの厚さに依存する）に応じてスケールすることを強調している。したがって、実験室の結果を、断層帯がはるかに厚いフィールド条件に適用するには、周波数ウィンドウをより低い値へとシフトさせる必要がある。なお、これらの知見は、今回研究された高浸透率の条件下でのものであり、低浸透率のガウジでは、間隙圧の拡散がより遅いため、異なる挙動を示す可能性がある。