[论文解读] An introduction to linear poroelasticity

Andi Merxhani|arXiv (Cornell University)|Jul 14, 2016

Drilling and Well Engineering参考文献 17被引用 9

一句话总结

本文基于Verruijt的完全固结形式化方法，系统推导了线性孔弹性理论，并将其扩展至Rice与Cleary的完全不透水形式化方法，引入了一种广义流体质量守恒方程，以考虑未包裹孔隙的可压缩性。该方法直接从物理原理推导出流体含量增量、储水系数和完全不透水材料常数的本构关系，避免了人为定义和能量势能的使用，并提供了适用于有限元分析的变分形式，同时给出了混合离散化的稳定性条件。

ABSTRACT

This study is an introduction to the theory of three-dimensional consolidation. The point of departure in the description are the basic equations of elasticity (i.e. constitutive law, equations of equilibrium in terms of stresses, and the definition of strain), together with the principle of effective stress, and the law of Darcy for fluid flow in porous media. These equations, together with the principle of mass conservation, are the only premises used to derive a drained description of consolidation theory. Next, the increment of fluid content is defined and its constitutive law is derived - with its derivation following naturally from the equation of fluid mass balance. Pore pressure, storage, and undrained coefficients are also introduced and useful relations are proven. Where appropriate, the physical meaning of these coefficients is proven mathematically. Equations of equilibrium and fluid mass conservation are subsequently expressed in terms of the increment of fluid content and undrained coefficients, leading to an undrained description of consolidation theory. Thus the initial approach of using drained description is extended to the formalism of Rice and Cleary. This approach to three-dimensional consolidation is useful for its simplicity. It does not require the ad-hoc definition of constants and that of an elastic energy potential. Instead, it is a direct extension to linear elasticity that accounts for the coupling of skeletal deformations and fluid behaviour.

研究动机与目标

• 基于基本原理而非人为假设，提供一个自包含且物理解释清晰的线性孔弹性理论推导。
• 通过在流体质量守恒方程中引入未包裹孔隙可压缩性，将Verruijt的完全固结形式化方法加以扩展。
• 通过流体含量增量和孔隙压力系数，推导出孔弹性理论的完全不透水描述，与Rice和Cleary的形式化方法保持一致。
• 建立适用于有限元法（FEM）的变分形式，并施加适当的离散化约束。
• 证明关键系数（如储水率、Skempton系数B及完全不透水模量）的物理解释与数学意义。

提出的方法

• 从各向同性弹性出发，基于有效应力原理、达西定律和质量守恒，推导出平衡方程与流体质量守恒方程。
• 将流体含量增量引入为与孔隙压力共轭的变量，从流体质量守恒方程中推导其本构关系。
• 将储水系数推广至包含未包裹孔隙可压缩性，扩展了Verruijt原始形式化方法。
• 利用流体含量增量的本构关系，推导出孔隙压力系数B与B′，以及完全不透水泊松比νu。
• 将平衡方程与流体质量守恒方程用流体含量增量和完全不透水系数表示，从而实现完全不透水描述。
• 为FEM开发控制方程的弱形式，通过混合有限元离散化（如Taylor-Hood元）满足Babuška-Brezzi条件，确保稳定性。

实验结果

研究问题

• RQ1 如何将孔弹性中的流体质量守恒方程广义化，以超越Verruijt原始形式化方法，纳入未包裹孔隙的可压缩性？
• RQ2 在线性孔弹性背景下，Skempton孔隙压力系数B与完全不透水系数B′的物理解释与数学基础是什么？
• RQ3 如何直接从流体含量增量的本构关系推导出完全不透水孔弹性常数（如Eu、Ku、νu）？
• RQ4 在固结问题中，当储水率Sε趋近于零时会产生何种影响，这又如何影响FEM中的数值稳定性？
• RQ5 为确保孔弹性混合有限元形式化的稳定性，需要哪些离散化约束条件？

主要发现

• 流体质量守恒方程被扩展以包含未包裹孔隙可压缩性，从而导出更一般的储水率Sε表达式，该表达式同时考虑了固体相与内部空隙的可压缩性。
• 孔隙压力系数B被推导为B = α / M，其中α为Biot系数，M为Biot模量，其物理解释与体积应变及流体含量变化相关。
• 完全不透水系数B′被推导为B′ = α / (M + (1 - α)Ks)，表明其依赖于固体相的可压缩性，且在单轴应变条件下退化为B。
• 单轴储水系数S被推导为S = m_v + (1 - α)Cϕ，其中m_v为完全固结可压缩性，Cϕ为未包裹孔隙可压缩性，该式直接关联到可测量的材料参数。
• 完全不透水体积模量Ku被推导为Ku = K + (1 - α)²M / (1 - α + M/K)，表明完全不透水条件具有刚化效应。
• 变分形式表明，稳定的FEM离散化要求压力场的形函数阶次低于位移场，以满足Babuška-Brezzi条件。