Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

地球の深部、石油、ガス、そして地熱を蓄える多孔質な岩石の中では、固体の岩石と、その微細な孔の中に閉じ込められた流体との間で、絶え間ない静かな戦いが繰り広げられています。エンジニアがこれらの資源を抽出しようとする際、しばしば岩石を押し広げるために流体を注入しますが、これは水圧破砕（ハイドリック・フラクチャリング）として知られるプロセスです。数十年にわたり、科学者たちはこのプロセスを、岩石を単なる弾性のあるスポンジのように、伸びては元に戻るものとしてモデル化してきました。しかし、実際の岩石はより複雑です。特に高い圧力で押しつぶされるとき、岩石は永久に変形したり、崩れたり、内部の亀裂に沿って滑ったりすることがあります。固体骨格と内部の流体の両方が押し引きを行う中で、岩石がどのように破壊されるかを正確に理解することは、井戸の安全な掘削、廃棄物貯蔵の管理、さらには氷山が割れたりマグマが地殻を押し上げたりといった自然界の地質学的イベントを予測する上でも極めて重要です。課題は、計算が不可能になるような複雑さに陥ることなく、岩石の破壊という乱雑な現実を捉える数学的記述を作り出すことにあります。

研究チームは、この破壊プロセスを記述するために用いられる基本的な方程式を再検討し、従来のモデルで使用されていた一般的なショートカット（近道）が重大な誤差を招いていることを発見しました。岩石力学の世界において、科学者は「フェーズフィールド」アプローチを用いています。これは、亀裂を鋭い目に見えない線としてではなく、材料が徐々に弱まっていく「ぼやけた領域」として扱う手法です。これにより、コンピュータは亀裂がどのように成長し、形状を変えるかをシミュレートできます。これを正確に行うためには、モデルは2つの異なる挙動を考慮しなければなりません。それは、岩石が引き離されるとき（引張）と、押しつぶされるとき（圧縮）です。引張では、微細な亀裂が開きます。圧縮では、それらは閉じ、内部の粗い表面が互いに摩擦を起こしながら滑ります。研究者たちは、既存の多くのモデル（産業界やアカデミアで広く使用されているもの）が、系のエネルギーを記述するための変数を混同していることを発見しました。具体的には、流体圧をエネルギー方程式における直接的な入力として扱うのではなく、孔内の流体量として扱うことで、数学的な「跳躍（ジャンプ）」または不連続性が生じることを突き止めました。この跳躍により、モデルは、岩石が引張から圧縮へと切り替わる際に、岩石の強度が不自然に急激に変化すると予測してしまい、亀裂がどのように、どこで形成されるかについての予測を不正確にするのです。

これを解決するために、研究者たちは、岩石の変形と流体含有量の関係を注意深く追跡する、新しい一貫した定式化を開発しました。彼らは、岩石を流体で満たされた微細な小銭型の亀裂の集合体として想定し、微視的な視点に基づいてモデルを構築しました。これらのマイクロクラックがストレスの下でどのように開閉するかを分析することで、圧力の変化に関わらず、岩石の強度表面（材料がいつ破壊されるかを定義する数学的な境界）が滑らかで連続的であることを保証する一連の規則を導き出しました。彼らは、特定のタイプの流動則を用いることで、岩石の永久変形を作用している力に直接結びつけ、さらにエネルギー方程式における欠落していた結合項を正しく含めることで、モデルが引張から圧縮への遷移を正確に記述できることを証明しました。この修正は単なる微調整ではありません。これは、モデルが亀裂を前進させる駆動力の計算方法を根本的に変えるものです。

チームは、新開発のモデルが通用するかどうかを確認するため、2つの全く異なるシナリオに対してテストを行いました。第一に、岩石に流体を注入して亀裂を作るという、典型的な水圧破法ベンチマークのシミュレーションを行いました。彼らは、理想的な条件下で導き出された精密な数学的解である既知の解析解と比較しました。新モデルは、以前の混合された定式化よりもこの解にずっと近く、亀裂を開始させるために必要な圧力と、成長に伴う亀裂の幅を正確に予測しました。古いモデルは、必要な流体圧を過大評価する傾向があり、これは非効率または不安全な掘削作業につながる欠陥でした。第二のテストでは、岩石が両側から押しつぶされる状況、つまり岩石が甚大な圧力にさらされている深部の地下環境を模したシミュレーションを行いました。ここでは、岩石は単にパチンと割れるのではなく、塑性変形を起こし、材料が自らの前を滑るような剪断帯を形成しました。新モデルは、機械的な圧縮によって引き起こされる剪断破壊と、流体注入によって駆動される引張破壊の両方を捉え、この複雑な挙動を再現することに成功しました。

これらのシミュレーションの結果は、岩石の塑性変形が亀裂の成長に対する重大な障壁として機能することを明らかにしています。岩石が永久に変形することを許容すると、岩石はより多くのエネルギーを吸収するため、亀裂を開いたまま維持するためには、より高い流体圧が必要になります。この追加の抵抗は亀裂のプロファイルを完全に変え、純粋な弾性モデルが示唆するよりも、亀裂をより狭くし、伝播させるためにより大きな力を必要とさせます。さらに、圧縮テストにおいて、モデルは岩石内の孔隙圧が亀裂形成に伴ってどのように変化するかを示しました。岩石が変形して流体が移動できるスペースが生まれるにつれ、局所的に圧力が低下し、損傷領域へとより多くの流体を引き込みます。この相互作用は、流体の流れと岩石の機械的破壊が密接に結合したフィードバックループを生み出しており、新モデルはこのダイナミクスを高精度に捉えています。

数学的基礎を修正することにより、研究者たちは、多孔質媒体を扱うエンジニアや科学者にとって、より信頼性の高いツールを提供しました。この新しい定式化は、異なる種類の岩石破壊の間の遷移を滑らかかつ物理的に現実的なものにし、これまでの試みが抱えていた人工的な跳躍を取り除きました。この正確さは、石油・ガス生産の最適化から、二酸化炭素回収のための安全な地下貯蔵のデザインに至るまで、幅広いアプリケーションにおいて不可欠です。本研究は、固体と流体が相互作用する複雑なシステムにおいても、基礎となる物理学への注意深く一貫したアプローチが、材料の破壊の真の性質を明らかにし得ることを示しています。このモデルは、理論的な力学と、足元にある地面の乱雑で塑性的な現実との間の溝を埋め、より高い確信を持って地殻の挙動を予測する準備が整っています。