Elastic wave propagation in fractured media with spring-type and frictional contact deformation laws

私たちの足元の地面は、めったに単一の、途切れることのない固形ブロックではありません。深い石油貯留層から地球の地殻そのものに至るまで、多くの場所で岩石は、数キロメートルにわたって伸びることもある薄い平面的な亀裂、すなわち割れ目によって交差しています。地震による地震波であれ、材料を検査するために使用される音響パルスであれ、エネルギーがこの割れた景観の中を伝わる際、これらの亀裂は物語（展開）を変えてしまいます。それらはエネルギーを散乱させ、信号を弱め、時には岩の面同士を擦り合わせたり、滑らせたりすることさえあります。波がどのようにこれらの亀裂と相互作用するかを正確に理解することは、地下に何があるかをマッピングしようとする地質学者にとっても、構造物の安全性を確保しようとする技術者にとっても極めて重要です。しかし、この相互作用の物理学を突き止めることは、非常に困難であることで知られています。単純な亀裂は単なる隙間ではありません。それは、波がほんの一瞬のうちに通過する間に、岩の面が押し付け合ったり、摩擦を伴って滑ったり、あるいはパカッと開閉したりする、複雑な境界なのです。

数十年にわたり、科学者たちはこれらの挙動を予測するために簡略化されたモデルに頼ってきました。それらはしばれて、亀裂をあたかも引き伸ばされたり圧縮されたりする目に見えないバネでできているかのように扱うことがよくありました。これらは有用ではありますが、波が十分に強く、岩の面を接触させたり、滑らせたり、あるいは完全に分離させたりする場合には、機能しなくなります。これらのモデルは、摩擦が岩をその場に留めておく力が強すぎるまで保持する「スティック・スリップ（固着と滑り）」運動を捉えることができず、また、亀裂が押しつぶされる際に起こる非線形な硬化を正確に記述することもできません。ベルゲン大学、ノルウェー研究センター、およびミュンヘン工科大学の研究者による新しい研究は、破砕された岩石の中を移動する弾性波を前例のないリアリズムでシミュレートする包括的な計算フレームワークを構築することで、これらのギャップに対処しています。簡略化された仮定に頼る代わりに、チームは、基本的なバネのような挙動から、摩擦や亀裂の物理的な開閉がリアルタイムで計算される完全な接触力学に至るまで、4つの異なるレベルの複雑さを扱える手法を開発しました。

研究者たちは、岩石と亀裂が単一の、相互に連結したシステムとして扱われる統一されたデジタル環境を構築しました。彼らのシミュレーションにおいて、岩石は三次元の体積として表現され、一方で亀裂は、その中に埋め込まれた二次元の表面としてモデル化されています。これにより、コンピュータは、岩石を通って伝わる波が亀裂にどのように当たり、亀裂がどのように変形し、その変形がいかにして岩石へと波紋を送り戻すかを追跡することができます。チームは、亀裂がどのように振る舞うかについて、4つの特定の方法をテストしました。最初の2つのモデルは、亀裂を横切る力が岩の動きに比例するというバネのような法則を用いたもので、一方は単純な線形関係を許容し、もう一方は、亀裂が閉じるときに硬くなるという事実を考慮した、より複雑な非線形ルールを用いています。他の2つのモデルは、摩擦と接触という現実を導入しました。これらのシナリオでは、亀裂の面は互いに接触し、車のタイヤがスリップするのを防ぐ法則に似た摩擦の法則に従って滑ることが許容されます。同時に、岩の面が互いに通り抜けることはできないというルールも課されています。

この新しいフレームワークの正確性を確保するために、チームは一連の厳格なテストを実施しました。彼らは単一の亀裂に波が当たる様子をシミュレートし、そのコンピュータの結果を既知の数学的解と比較しました。その結果、彼らの手法は、波の速度と、透過または反射するエネルギーの量の両方について、理論的な期待値と高い精度で一致することを発見しました。次に、彼らはシステムをさらに押し進め、亀裂を開いたり、固着させたり、滑らせたりするのに十分な強さを持つ波をシミュレートしました。亀裂の面が摩擦係数1.0で互いに擦り合うという、より複雑なシナリオにおいても、シミュレーションは安定しており、一貫した結果を示しました。研究者たちは、亀裂が「固着」から「滑り」へと切り替わる際に起こる物理学の急激な変化を、彼らのアプローチが処理できることを検証しました。この遷移は、単純なモデルをしばしば躓かせます。

チームが4つの異なるモデルを並べて比較したとき、結果として生じる波のパターンの違いは驚くべきものでした。シミュレーションは、モデルの選択が極めて重要であることを明らかにしました。亀裂が非線形ルールに従って閉じ、硬化したとき、単純な線形バネでモデル化された亀裂よりもはるかに多くの波エネルギーを伝達しました。これは、非線形モデルが、亀裂が閉じるにつれて岩の面がより密接に押し付けられ、亀裂が実質的に硬くなり、波に対してより透明（透過的）になるという事実を考慮しているためです。さらに、摩擦の導入は波面の形状を変化させました。摩擦モデルにおいては、亀裂の反対側に現れる波は、摩擦のないバージョンよりも幅広く、歪んでいました。研究者たちは、摩擦によって亀裂の面がより長い時間、互いに固着し、エネルギーの放出を遅らせ、波の散乱の仕方を変えることを観察しました。

研究は単一の亀裂に留まりませんでした。研究者たちは、彼らのフレームワークを、異なる特性を持つ6つの交差する亀裂を含む三次元領域に適用し、小さな亀裂のネットワークを作成しました。この複雑な設定において、彼らは異なる亀裂のセットに対して異なる最大閉鎖限界を割り当てました。つまり、ある亀裂は他のものよりもきつく閉じることができるように設定しました。シミュレーションによれば、波の場は高度に不均質になりました。より硬く、よりきつく閉じることができる亀裂は、より多くのエネルギーを反射し、より少なく伝達し、波のパターンの中に明確な影を作り出しました。研究者たちは、各亀裂がどの程度閉じ、どの程度滑ったかを正確に追跡することができ、亀裂の物理的特性のわずかな違いが、目に見えて異なる変形反応をもたらすことを明らかにしました。例えば、よりきつく閉じることが許容された亀裂は、より多くの「固着（スティック）」挙動を示し、滑る前にその面をより長く保持し、それがそれを通過する波に影響を与えました。

これらの発見の含意は、直接的な結果を超えて広がっています。単純なバネのような挙動から複雑な摩擦接触まで、あらゆる事象を扱うことができる単一の計算フレームワークを示すことで、研究者たちは、現実世界の地震データをより高い忠実度で解釈するために使用できるツールを提供しました。摩擦、開口、および閉鎖を伴う亀裂の中で波がどのように散乱し減衰するかをシミュレートできるということは、地球科学者が異なる種類の岩石形成や流体で満たされた亀裂をより正確に区別できることを意味します。本研究は、亀裂力学の非線形および摩擦的な側面を無視することが、地下におけるエネルギーの移動を予測する上で重大な誤差につながる可能性があることを裏付けています。現在の研究は純粋に機械的な相互作用に焦点を当てていますが、このフレームワークは流体の流れのモデルと互換性のある基礎の上に構築されており、水や石油が亀裂の中を移動しながら通過する波と相互作用する、より現実的なシミュレーションへの道を示唆しています。この研究は、岩石の亀裂の振る舞いにおける完全な複雑さを受け入れることで、足元に隠された世界をより鮮明に描き出すことができるということを示す実証となっています。