[论文解读] 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 方法求解三维 Biot 多孔弹性方程，实现对各类波的独立计算。
• 通过提供多孔弹性波问题的精确参考解，验证数值模拟代码（如 Gar6more 3D）的可靠性。
• 分析快 P 波、慢 P 波和 S 波在流体饱和固-固界面处的反射与透射特性。
• 为理解非均质多孔介质中波的频散与衰减特性提供理论基础，相关研究对地震勘探具有重要意义。

提出的方法

• 建立三维 Biot 多孔弹性介质方程，包含固体位移场、流体位移场及孔隙压力。
• 应用 Cagniard-de Hoop 方法将问题转换至拉普拉斯-傅里叶域，通过围线积分实现解析反演。
• 推导点源在频率-波数域的格林函数，并进一步反演为时域解。
• 利用该方法分别计算快 P 波、慢 P 波和 S 波的贡献，其具有不同的传播时间和振幅。
• 通过源函数的导数实现反演积分的数值计算，避免对位移场进行直接微分。
• 通过包含真实材料参数和多种源类型（体源与压力源）的数值算例验证解的正确性。

实验结果

研究问题

• RQ1 Cagniard-de Hoop 方法如何扩展以推导多孔弹性分层介质中三维波传播的解析解？
• RQ2 快 P 波、慢 P 波和 S 波对双层多孔弹性系统总位移场的贡献有何差异？
• RQ3 在三维多孔弹性构型中，波的振幅与传播时间如何随源类型（体源 vs. 压力源）和材料参数变化？
• RQ4 该解析解在多孔弹性波模拟中作为三维数值代码验证基准的适用程度如何？
• RQ5 在三维条件下，多孔弹性/多孔弹性与声学/多孔弹性层状结构中的波行为存在哪些关键差异？

主要发现

• 该解析解能够分别且精确地计算出所有波型——快 P 波、慢 P 波和 S 波，其具有不同的传播时间与振幅。
• 在给定材料参数下，上层介质中快 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 等波型的独立峰值。
• 该解与数值模拟结果高度一致，证实其作为代码验证基准的实用性。
• 采用高斯源函数的四阶导数可有效激发主导频率 15 Hz，从而在时域响应中清晰识别出波包。