地震中の岩の挙動:シミュレーションからの知見
研究がコンピュータシミュレーションを使って、地震のストレス下で岩石がどう変わるかを探っている。
― 1 分で読む
地震が起こると岩石が変わるんだ。地震からのエネルギーである地震波が岩石を通ると、岩が柔らかくなることがよくあるんだ。この柔らかくなる現象は、実験室のテストでも自然界でもよく観察されてる。揺れが収まった後、こうした岩は一般的に時間が経つにつれてまた硬くなり始める。研究者たちは、これがどうして起こるのか、いくつかの理論を提案しているんだ。
この研究では、こうした現象の間に岩がどう振る舞うかについての2つの既存の理論をテストしてる。コンピュータシミュレーションを使って、地震波が岩をどのように移動するかをモデル化している。そして、これらのシミュレーションの結果を実際の実験室テストと比較すると、どちらのモデルも現実の観察に対して有望な精度を示しているんだ。
何千ものシミュレーションを通じて、面白いパターンが浮かび上がってくる。岩の初期状態や、岩の中にできては治る小さな亀裂が、この柔らかさと硬さのプロセスに重要な役割を果たしている。これらのコンピュータモデルは、地震中の地殻の振る舞いを予測する精度を向上させることを目的としているんだ。
非古典的非線形性の理解
岩石やセメントのような材料は、普通は見られない振る舞いをすることがある。この現象を非古典的非線形性って呼ぶんだ。つまり、これらの材料がストレスやひずみにどう反応するかが、弾性材料で見られる通常のパターンには当てはまらないってこと。一つの特定の非古典的非線形性のタイプは、ヒステリシスモデルを含んでいる。ヒステリシスモデルは、岩のストレスが現在の状態だけでなく、変形の歴史にも依存することに焦点を当てているんだ。
実験室の実験では、この非古典的非線形性の兆候が示されている。一種の実験は、非線形共鳴超音波分光法(NRUS)って呼ばれてる。NRUSでは、岩は同じ周波数で興奮したとき、周波数が上昇中か下降中かによって異なる振る舞いをする。
もう一つの実験、動的アクースティック弾性試験(DAET)では、動的ストレス下での岩の弾性特性の変化を調べている。観察によると、こうしたテスト中には、地震波の速度が予測可能な形で変わり、ストレスがかかるときと外されるときで異なるフェーズを経るんだ。
現場からの観察
フィールドスタディでも、地震中やその後に地震波が周囲の岩の速度を変えることが示されている。たとえば、大きな地震のときには地震波の速度が急激に減少することが記録されている。これは大きな地震だけでなく、小さな中間的なものでも起こるんだ。チリや日本の特定の地域の研究では、こうした地震速度の低下とその後の回復が記録されているんだ。
これらの観察にもかかわらず、実験室での発見が現実のシナリオにどのように関連するのかは不明なまま。研究者たちは実験室の実験とフィールド観察のスケールの違いによる課題に直面することが多い。しかし、数値モデルはこのギャップを埋める助けになり得るんだ。
岩の振る舞いの物理モデル
計算地震学では、岩が地震にどう反応するかを分析するための異なるモデルが存在している。一部のモデルは、非線形特性が地震の動態や地面の動きに与える影響に焦点を当てている。これらのモデルは、岩材料の複雑な振る舞いや相互作用を考慮する必要があるんだ。
この研究は、動的プロセス中に岩の損傷がどのように進化するかを記述する2つのモデルに焦点を当てている。内部変数モデル(IVM)は、材料の変化が岩内の損傷と相関する変数によって捉えられると提案している。連続損傷モデル(CDM)は、過剰なパラメータを必要とせずに損傷を全体の材料の振る舞いと結びつけているんだ。
モデルの実装
IVMとCDMモデルは、岩がストレスにどう反応するかを記述するための2つの異なる方法を提供している。両方のモデルは、実際の実験室データに対抗するための複雑なシミュレーションソフトウェア内に実装されたんだ。
シミュレーションでは、地震波が生成され、その結果のデータが両方のモデルで分析される。各モデルは、地震波が材料を通過する際の特定の損傷のシーケンスに焦点を当てて、岩がどう振る舞うかを解釈するようにセットされているんだ。
実験の設定
これらのモデルを評価するために、岩が地震活動を模した超音波信号にさらされる実験室のセットアップが作成される。2種類の超音波トランスデューサーが使用され、岩のサンプルに異なる周波数と振幅を適用する。これらの信号に対する岩の反応を測定するんだ。
実験中、研究者は岩の音響弾性率が時間とともにどう変化するかを追跡する。これは、岩の誘導ひずみ、結果として生じる損傷、そしてそれに続く回復フェーズの関係を研究することを含むんだ。
パラメータの評価
ベイジアン反転技術を利用して、研究者はモデルパラメータが実験室テストからの観察データとどれだけ一致しているかを定量化できる。この方法により、さまざまなパラメータに関連する不確実性を評価でき、結果の信頼性が向上するんだ。
モデルは、実世界の観察にどうフィットするかを発見するために、さまざまな入力パラメータに対してテストされる。これらの計算を通じて、モデルのさまざまな側面が、地震活動中およびその後の岩の挙動の理解にどう貢献するかが明らかになるんだ。
シミュレーションの結果
両方のモデルは、動的テスト中の岩の観察された挙動を再現できることを示した。損傷と回復の重要なフェーズを捉えたんだ。IVMとCDMモデルは実験測定と密接に相関していたが、CDMモデルは観察された損傷パターンのより包括的な説明を提供したんだ。
実験室からの測定結果は、適用されたストレスの周波数や材料の初期条件など、異なるパラメータが結果として生じる損傷と回復パターンに重要な役割を果たすことを示している。
結論と今後の方向性
結果は、特にCDMモデルが地震イベント中の岩の複雑な挙動についての重要な洞察を提供するかもしれないことを示唆している。自然界で見られるより複雑なシナリオを組み込めるように、大規模シミュレーションに向かう必要があることを強調しているんだ。
さらなる調査は、これらのモデルを洗練させ、変化する地形や地質の不均一性など、現実の要因を考慮できるようにすることを目指している。これらのシミュレーションの精度を向上させることで、研究者たちは地震中の地殻の挙動に関する予測を改善できることを望んでいるんだ。
これらのモデルを通じて岩の力学の理解が進むことで、広範な影響をもたらすことができる。地震に関連する損傷の軽減戦略や、活断層の近くにある都市部の安全対策が向上する可能性があるんだ。
今後は、実験室の実験とフィールドスタディの発見を数値モデルと統合することが重要になる。この包括的なアプローチを通じて、科学者たちは地震ストレス下での岩の複雑な挙動を完全に理解し、社会への地震の影響を予測する改善を目指すことができるんだ。
タイトル: Modeling and Quantifying Parameter Uncertainty of Co-seismic Non-classical Nonlinearity in Rocks
概要: Dynamic perturbations reveal unconventional nonlinear behavior in rocks, as evidenced by field and laboratory studies. During the passage of seismic waves, rocks exhibit a decrease in elastic moduli, slowly recovering after.Yet, comprehensive physical models describing these moduli alterations remain sparse and insufficiently validated against observations. Here, we demonstrate the applicability of two physical damage models - the internal variable model (IVM) and the continuum damage model (CDM) - to provide quantitative descriptions of nonlinear co-seismic elastic wave propagation observations. We recast the IVM and CDM models as nonlinear hyperbolic partial differential equations and implement 1D and 2D numerical simulations using an arbitrary high-order discontinuous Galerkin method. We verify the modeling results with co-propagating acousto-elastic experimental measurements. We find that the IVM time series of P-wave speed changes correlate slightly better with observations, while the CDM better explains the peak damage delay relative to peak strain. Subsequently, we infer the parameters for these nonlinear models from laboratory experiments using probabilistic Bayesian inversion and 2D simulations. By adopting the Adaptive Metropolis Markov Chain Monte Carlo method, we quantify the uncertainties of inferred parameters for both physical models, investigating their interplay in 70,000 simulations. We find that the damage variables can trade off with the stress-strain nonlinearity in discernible ways. We discuss physical interpretations of both damage models and that our CDM quantitatively captures an observed damage increase with perturbation frequency. Our results contribute to a more holistic understanding of non-classical non-linear damage with implications for co-seismic damage and post-seismic recovery after earthquakes.
著者: Zihua Niu, Alice-Agnes Gabriel, Linus Seelinger, Heiner Igel
最終更新: 2023-06-07 00:00:00
言語: English
ソースURL: https://arxiv.org/abs/2306.04197
ソースPDF: https://arxiv.org/pdf/2306.04197
ライセンス: https://creativecommons.org/licenses/by/4.0/
変更点: この要約はAIの助けを借りて作成されており、不正確な場合があります。正確な情報については、ここにリンクされている元のソース文書を参照してください。
オープンアクセスの相互運用性を利用させていただいた arxiv に感謝します。