多孔介质弹性力学.docxVIP

多孔弹性力学 迄今为止，我们一直把岩石当作均质性固体物质处理。然而，岩石通常是复杂的物质，因此，在微观规模上， 它是非均质性的。在很大程度上，岩石表现的弹性反应、破裂应力等都依赖于其非固体部分。在本章中，我们将 考虑岩石的孔隙体积，其不仅在储层的石油开采中是必须的，而且在岩石的机械特性方面也起到重要的作用。首 先，我们对多孔和渗透性介质进行宏观描述，其允许我们去研究静态和动态的机械特性。这个处理方法基于Maurice A.Biot 的理论。 液体中固体颗粒的悬浮 首先，让我们看一个非常简单的多孔介质；即在这个多孔介质中，固体和流体部分相互之间分别产生形变。 实际上，我们可以把这个介质看成是流体中固体颗粒的悬浮，或例如一个饱水的、完全未固结的砂岩。如果我们 把这个混合物放进一个容器中，由于外挤压力的作用而产生的体积应变是： ? ? ? p v K eff （1.79） ? 是混合物的体积应变。总变形必须等于每一部分变形的总和，其值来自于每一分量的体积部分。 v ? ? Vs ?  v?s Vf ? v? f （1.80） v V tot V tot 下标 s 和 f 分别指固体和流体，V 指总体积。现在，我们定义孔隙度为流体所占体积与总体积的比值，即： tot V ? ? f V （1.81） 因为V ? V ? V V ＝V ＋V ，故而固体所占体积为： tot s f tot f s tot Vs ? 1 ? ? =1－Ф V  （1.82） tot 应变? 和? 分别由固体的体积模量K 和流体的体积模量K 得到，根据公式（1.49），公式（1.81）可写成： v,s v,f s f （1.83） 结合公式（1.79）和（1.83），此时我们发现悬浮的有效模量是： 1 1 ?? ? = + K K K eff s f （1.84） 这是一个特别简单的多孔物质的实例。现在，我们要通过考虑岩石由固体框架和孔隙流体组成的事实进行总结， 固体框架和孔隙流体不能被分别处理。 毕奥静力特性多孔弹性力学理论 现在，我们将考虑一个各向同性，多孔和渗透性的介质，该介质由两部分组成：固体和流体。固体的位移用 u s us 表示，流体的位移用u f uf 表示。对于一个附着在固体上的体积单元，应变由us 的组分的导数得到。利用公 式（1.41），我们得到体积应变：  ? ? ? ? ? ? u v s  （1.85） 对于流体部分，我们定义一个应变参数? ，其描述了流体相对于固体的体积变形 ? ?? ? ? ? ? ? (u s ? u ) （1.86） f 应力张量? 表示作用在附着于固体骨架上一个体积单元的总外应力。体积单元由固体骨架上的应力和流体中的 静水压力即孔隙压力 p f 来平衡这个应力。按照符号惯例，所有的应力（包括孔隙压力）在压方向上为正。 孔隙流体的存在使材料的应变能增加了额外的限制条件。因此，应力应变关系（公式（1.48））也将被修正。 毕奥（1962）证明了二物相系统的线性应力应变关系，这个二物相系统可以用应变参数? 和? 来表示，即应力 v 张量单元和孔隙压力 p ： f  ? ? ?? x v ? ? ?? y v  ? 2G? x ? 2G? y  ? C? （1.87） ? C? （1.88） ? ? ?? z v ? 2G? z ? C? （1.89） ? ? 2G? xy xy ? ? 2G? yz yz ? ? 2G? zx zx （1.90） （1.91） （1.92） P ? C? f v ? M? （1.93） ? 和 G 是多孔介质的Lame 参数，C 和 M 是为描述二物相系统而增加的弹性模量。如我们在公式（1.87）－（1.89） 和公式（1.93）看到的那样，C 在固体和流体形变中起作用，M 粗略形象地描绘了多孔流体的弹性特性。关于固体和流体弹性模量的确定表达式在以后的叙述中会给出。由公式（1.87）－（1.89）可以累加得出： ? ＝K? －C? （1.94） v ? 是由公式（1.9）定义的。K＝ ? ＋2 G/3 是不排水条件下多孔岩石的体积模量，不排水条件指的是流体不允许漏泄的地方。 现在，我们对弹性模量K，C 和 M 怎样与岩石组分模量相联系做一个试验。首先，设想我们做一个夹套实验 （见图 1.20a）：一多孔介质置于一不透水的夹套内，其承受一外部静水压力? 。加载过程中允许孔隙流体泄漏， p 因此孔隙压力保持不变，从而应力完全由固体骨架来传递。从公式（1.93）和（1.94），我们得到： p? ? K ? C2 ? K p ? M fr v  （1.95） 因为这个实验描述了岩石固体部分的硬度，K 称为骨架体积模量。因为没有与流