[論文レビュー] Analytical Solution for Wave Propagation in Stratified Poroelastic Medium. Part II: the 3D Case

Julien Diaz, Abdelaâziz Ezziani|ArXiv.org|Jul 25, 2008

Seismic Imaging and Inversion Techniques参考文献 8被引用数 8

ひとこと要約

本稿では、Cagniard-de Hoop法を用いて、二層構造の多孔質弾性体における三次元波動伝播の解析的解を提示する。反射波および透過波の明示的表現（高速P波、遅いP波、S波を含む）を導出し、数値シミュレーションコードの正確な検証および多孔質媒体内の波動挙動の詳細な分析を可能にする。

ABSTRACT

We are interested in the modeling of wave propagation in poroelastic media. We consider the biphasic Biot's model in an infinite bilayered medium with a plane interface. We adopt the Cagniard-De Hoop's technique. This report is devoted to the calculation of analytical solution in three dimension.

研究の動機と目的

• 三次元二層構造多孔質弾性体における一時的波動伝播の解析的解の開発を目的とする。
• Cagniard-de Hoop法を用いて3次元のBiotの多孔質弾性方程式を解き、波の種類を独立して計算可能にする。
• 多孔質弾性波問題の正確な基準解を提供することで、数値シミュレーションコード（例：Gar6more 3D）の検証を可能にする。
• 流体飽和固体界面における高速P波、遅いP波、S波の反射および透過特性を分析する。
• 地震探査に関連する不均質な多孔質媒体における波動分散および減衰の理解の基盤を提供する。

提案手法

• 固体および流体の変位場と細孔圧力を含む、多孔質弾性体の3次元Biot方程式を定式化する。
• Cagniard-de Hoop法を用いて問題をラプラス＝フォーリエ変換領域に変換し、経路積分を用いた解析的逆変換を可能にする。
• 周波数波数領域における点源のグリーン関数を導出し、それを時間領域解に逆変換する。
• 高速P波、遅いP波、S波のすべての波の寄与を、異なる伝播時間と振幅で計算する。
• 変位の直接微分を回避するため、源関数の導関数を用いて逆積分の数値評価を実装する。
• 実際の材料定数とさまざまな源タイプ（体積力および圧力源）を用いた数値例を通じて解の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 Cagniard-de Hoop法をどのように拡張することで、多孔質層状媒体における3次元波動伝播の解析的解を導出できるか？
• RQ2 高速P波、遅いP波、S波が、二層構造多孔質弾性系における全変位場に与える寄与はどのように異なるか？
• RQ3 3次元多孔質弾性構造において、源タイプ（体積力対圧力）および材料特性が波の振幅と伝播時間に与える影響は何か？
• RQ4 解析的解が、多孔質弾性波モデリングにおける3次元数値コードの検証基準としてどの程度有効に機能するか？
• RQ5 3次元において、多孔質弾性／多孔質弾性および音響／多孔質弾性層構造の間で波動挙動にどのような主要な差が生じるか？

主な発見

• 解析的解は、高速P波、遅いP波、S波をすべて別々にかつ正確に3次元で計算でき、それぞれ異なる伝播時間と振幅を持つ。
• 与えられた材料定数のもとで、上層では高速P波の位相速度が2692 m/s、下層では2535 m/sであり、遅いP波はそれぞれ1186 m/sおよび744 m/sである。
• S波の位相速度は上層で約1409 m/s、下層で1415 m/sであり、この構成では剪断波分散がほとんどないことが示された。
• 数値結果から、時間領域信号における波の到着が明確に分離されており、PfPf、PfS、PsPf、PfPs、PsS、PsPs波の種類に対応する明確なピークが観察された。
• 解は数値シミュレーションと優れた一致を示し、コード検証の基準としての有効性が確認された。
• ガウス源関数の4階微分を用いることで、主な15 Hz周波数が効果的に励起され、時間領域応答における波パッケージの明確な特定が可能になった。