1. 引言
Terzahi (1923)的有效应力定理揭示了孔隙压力在多孔弹性土材料中能引起土材料明显的体积增大,在数量级上类似于围压引起的体积变化,与围压的方向相反,并且土壤在应力条件下的应变和破坏取决于压力差(有效应力) [1]
(1)
式中:Pc是围压,且
,σ11、σ22和σ33分别为x、y、z三个方向的主应力,Pp是孔
隙流体压力,Pd是压力差(有效应力)。不过关于Terzahi (1923)的有效应力定理的精确性和有效性方面存在大量的分歧,部分异议偶尔来自多孔骨架应力和应变的含糊定义。
Biot和Willis (1957)根据实验室测量,在保持多孔材料孔隙压力不变的情况下,建立了样品应变与应力的关系,提出了基于各向同性材料的有效应力定理 [2]
(2)
式中:Pc是围压,且
,Pp是孔隙流体压力,Pe是有效应力。Kdry为多孔干燥岩石
的体积模量(排水状态下的体积模量),Kma为岩石基质的体积模量,n是有效应力系数(也称为Boit或Biot-Willis系数,有的文献中用α表示有效应力系数)。在两种极端情况下:1) 当孔隙度
时,Kdry = Kma,n = 0为固结良好的沉积岩、致密的沉积岩。2) 当孔隙度
≥
(
为临界孔隙度)时,Kdry = 0,n = 1为未固结的沉积物或流体,因此
(Hilterman, 2001) [3] 。Nur和Byerlee (1971)证实了Biot和Willis (1957)提出的有效应力定理的有效性 [4] ,事实上,Biot和Willis (1957)提出的有效应力定理式(2) [2] 包含了Terzaghi (1923)提出的有效应力定理式(1) [1] ,地面疏松干燥的多孔土壤的体积模量Kdry往往远小于土壤矿物颗粒的体积模量Kma。Sarker和Batzle (2008)的研究也证实在较浅的地层中
,n ≈ 1;因此在较浅的地层中上述式(2)可简化为式(1) [5] 。
Nur和Byerlee (1971) [4] 根据Terzahi (1923) [1] 和Biot (1941)的研究结果 [6] ,给出了一个基于统计学各向同性多孔集料的应力和孔隙压力的应变的张量通式
(3)
式中:δij是克罗内克尔符号,重复标记表示求和,即
,在式(3)中压缩为正方向。式(3)
右边的第一项是偏应力导致的应变,仅依赖于没有孔隙压力的有效剪切模量。第二项是由静水应力(平均应力σττ/3)单独导致的应变,仅依赖于无孔隙压力状态下岩石的有效体积模量K。最后一项是由孔隙压力导致的应变,并且依赖于由Biot (1941) [6] 介绍的一个新的有效模量H。有效模量K、μ和H是他们自己定义的,比较抽象,不是很明确,与实验室的测试是一致的。正如早期所定义的,它们是施加的平均应力或孔隙压力与被测量的平均应变之间的比值,方程(3)表明,在各向同性多孔集料中,孔隙压力影响法向应变,而不是剪切应变。Skempton (1954)的实验结果指出,对于各向同性多孔集料,孔隙压力仅仅依赖于平均主应力,但是在土壤中纯剪切应力对孔隙压力有显著的影响 [7] 。因此现有关于孔隙压力的预测模型都是基于各向同性多孔集料的,孔隙压力只与主应力有关,与剪切应力无关。
2. 多孔岩石可压缩性定义
为了解析Biot和Willis (1957)提出的有效应力定理 [2] ,在此我们引入Zimmerman (1991)给出的多孔岩石四个可压缩性的定义。对于多孔岩石而言,Pc是施加给岩石的围压,Pp是岩石孔隙流体压力;Vb为岩石的总体积,Vp为岩石的孔隙体积,Vm为岩石基质(岩石固体部分)的体积(图1) [8] 。
![](//html.hanspub.org/file/8-1770763x21_hanspub.png)
Figure 1. Confining pressure Pc and pore pressure Pp acting on a porous body [8]
图1. 围压Pc和孔隙压力Pp对多孔岩石的作用力示意图 [8]
(4)
(5)
(6)
(7)
并且 [8]
(8)
(9)
(10)
上述式(4)~(10)中:Cbc相当于多孔岩石在排水状态下的可压缩性,Cbc = 1/Kdry;Cbp被称为虚拟体积可压缩性(Fatt, 1958) [9] ,这个可压缩性被用着地层沉陷的计算(Geertsma, 1973) [10] ;Cpc和Cpp是孔隙的可压缩性,指压力的变化对岩石孔隙空间体积的影响;Hall (1953)指出Cpc作为“地层压实”系数 [11] ,孔隙可压缩性Cpp被用着油藏分析,因为它反映了超过孔隙流体的体积能被储存在孔隙空间中是由于孔隙压力的增加;Cm是岩石基质的可压缩性,Cm = 1/Kma;f是岩石的孔隙度。在这里围压Pc作用的方向为负方向(压缩),孔隙流体压力Pp作用的方向为正方向。
此外,我们还需要引入多孔岩石三个可压缩性定义(Green和Wang,1986) [12]
(11)
(12)
(13)
上述式(11)~(13)中:Cu是饱和岩石的可压缩性,Cu = 1/Ksat,Ksat是饱和岩石的体积模量;Cm是压力差Pd保持不变(
)时的岩石可压缩性,岩石形变是由内在的矿物模量来描述的,一般定义为岩石基质的可压缩性;
是围压变化量dPc与孔隙流体压力变化量dPp相等时的孔隙可压缩性,
是与之相对应的孔隙刚度(孔隙的体积模量),
。
在这里需要强调一下的是式(13)中岩石孔隙可压缩性
与Zimmerman (1991) [8] 定义的岩石孔隙可压缩性式(6)和式(7)是不同的,式(13)中岩石孔隙可压缩性
是围压变化量dPc与孔隙压力变化量dPp两者相等时岩石孔隙的可压缩性。
3. 有效应力定理的力学机制
对于均匀线性弹性固体内部的单个或多个任意形状的孔隙空间,有孔固体在干燥状态下的等效可压缩性(干燥体积模量的倒数)可以描述为(Walsh, 1965 [13] ; Zimmerman, 1991 [8] )
(14)
式中:Kdry为干燥岩石的等效体积模量;Kma为岩石基质的体积模量;
为岩石孔隙度;
为干燥孔隙空间的刚度。式(14)的物理意义是干燥有孔岩石的可压缩性由两部分构成,一是矿物内在的可压缩性;二是由孔隙空间的存在所引起的附加可压缩性。
根据上述式(2)式(14)得
(15)
式(15)描述了饱和岩石在围压的作用下,孔隙压力保持不变的情况下的三维力平衡模型,同时式(15)也说明了,围压作用于饱和岩石,保持孔隙压力不变的情况下,岩石的总应变等于岩石骨架(岩石排水状态)的总应变,岩石骨架(岩石排水状态)的总应变等于有效应力对岩石骨架引起的应变与岩石孔隙分担的孔隙压力引起的孔隙应变之和。对于固结的岩石,如果对多孔岩石施加围压Pc,当孔隙中没有流体时,孔隙压力Pp = 0,有效应力是作用在整个岩石骨架(干燥岩石)上的,而不是作用在单个岩石矿物颗粒上的(图2)。
![](//html.hanspub.org/file/8-1770763x42_hanspub.png)
Figure 2. The diagrammatic sketch of three-dimensional stress of the saturated rock
图2. 饱和岩石三维力平衡示意图
对于地下沉积岩而言,在弹性范围内,岩石的总应变受围压Pc和孔隙流体压力Pp控制
(16)
那么岩石总应变增量dεb为
(17)
岩石的微小形变可认为是线性的,根据式(8)和式(17)得
(18)
如果对饱和岩石施加围压Pc的过程中,孔隙流体压力Pp保持不变,假如岩石的形变是线性的,那么岩石的总应变εb等于岩石骨架(排水状态的岩石)的应变,由此可根据式(4)和式(18)得
(19)
由式(19)得
(20)
式(20)等价于Biot和Willis (1957) [2] 提出的有效应力定理式(2)。
如果对饱和岩石施加围压Pc的过程中,孔隙流体与外界不连通,孔隙流体的质量保持不变,那么在围压Pc增加的过程中,假如岩石的形变是线性的,岩石的总应变εb等于整个饱含流体岩石的应变,由此可根据式(11)和式(18)得
(21)
由式(21)得
(22)
式(22)描述了岩石在线性形变状态下,孔隙流体与外界不连通,并且孔隙压力增大时,孔隙压力Pp与施加围压Pc的关系。
4. Gassmann方程隐含的力学机制
沉积物由水动力搬运至湖盆中,在压实和成岩作用之前的最初沉积,矿物颗粒是疏松的,甚至是弱悬浮的,此时沉积物所具有的孔隙度称为临界孔隙度
。当孔隙度低于
时,矿物承载,而当孔隙度大于
时,岩石颗粒之间将失去接触和刚度,岩石“散开”并且变成承载液体中的悬浮物,岩石矿物颗粒所承受的有效应力与孔隙压力相等,
的值取决于沉积过程中颗粒的分选和磨圆(Mavko et al., 1998) [14] 。沉积岩在压实成岩过程中,埋藏深度逐渐增大,孔隙度逐渐减小,弹性刚度增加。当孔隙度小于
之后,岩石矿物颗粒要分担围压施加的应力,并且孔隙流体压力不断增大。事实上,沉积岩在压实成岩的过程中往往是一个非线性形变过程,这一非线性形变过程主要表现在孔隙度逐渐减小,饱含流体岩石的体积模量和与之对应的干岩石骨架的体积模量会逐渐增大(弹性刚度增加)。用式(17)描述沉积岩的孔隙流体力学机制要比式(18)更合理一些,对于饱和岩石而言,根据式(4)、式(11)和式(17)得
(23)
由式(23)得
(24)
式(24)是Skempton系数B (Skempton, 1954) 的 [7] 另外一种表达形式,描述了岩石在不排水的状态下围压增加所引起孔隙压力增加的关系,在数学上是孔隙压力Pp对围压Pc的导数,dPp/dPc = B。
来源于Zimmerman (1991, p. 62) [8] 和Jaeger 等人(2007, p. 177) [15] 推导的Skempton系数B表达式为
(25)
式(24)等价于式(25),并且式(24)和式(25)可还原于Gassmann方程 [16] 中
(26)
式(26)是Gassmann方程(Gassmann方程可以写成多种表达形式),由式(24)和(25)可知,Gassmann方程还可写成
(27)
式(27)揭示了Gassmann方程的力学和声学机制,是地震数据应用于油藏监测的理论基础,涉及到孔隙流体压力、地层应力和地震波反射信息。
除此之外,根据Gassmann方程还可以得到与式(14)相类似的结论。对于均匀线性弹性固体内部的单个或多个任意形状的孔隙空间,有孔固体在饱和状态下的等效可压缩性(饱和体积模量的倒数)可以描述为 [14]
(28)
式中:Ksat是饱和岩石的体积模量;Kma为岩石基质的体积模量;Kf是孔隙流体体积模量;
为岩石孔隙度;
为干燥孔隙空间的刚度。式(28)的物理意义是饱和岩石的可压缩性由两部分构成,一是矿物内在的可压缩性,二是由孔隙空间和孔隙流体耦合所引起的附加可压缩性。在干燥岩石中(岩石排水状态下),孔隙中无流体,流体体积模量Kf = 0,那么式(28)等价于式(14)。由此可见,同样矿物成分和孔隙度的干燥岩石,注入流体饱和后,岩石的体积模量会增大,通过比较式(14)和式(28),饱和岩石相对于干燥岩石增加的体积模量主要取决于孔隙流体与矿物本身的耦合效应。尽管Gassmann方程能较精确地描述了含流体岩石弹性形变状态下的力学机制,但是Brown和Korringa (1975) [17] 认为Gassmann方程建立在一个重要的理论假设基础之上。
Brown和Korringa (1975)放宽了Gassmann方程的各向同性的内在假设,并且推导出广义的Gassmann方程 [17]
(29)
式中:Cd是岩石排水状态下的可压缩性,Cd = 1/Kdry;Cu是岩石不排水状态下的可压缩性,Cu = 1/Ksat;Cm是岩石基质的可压缩性,Cm = 1/Kma;
是岩石孔隙可压缩性,与式(13)一致,
,
是孔隙体积模量;Cf是孔隙流体可压缩性,Cf = 1/Kf,Kf 是孔隙流体体积模量;
是岩石的孔隙度。
Brown和Korringa (1975)描述了岩石在不排水的状态下,Skempton系数B值(Skempton, 1954)的表达式为 [17]
(30)
在式(29)和式(30)中,如果
,那么式(29)就可简化为Gassmann方程
(31)
同时Gassmann方程还可以写成
(32)
式(31)等价于上述式(26),式(32)等价于上述式(24)、(27)。
从Brown和Korringa (1975) [17] 给出的广义Gassmann方程可知,Gassmann方程是基于
的假设条件,
很难直接地测量,
是多孔岩石的外部和内部承受同样的静水压力(dPc = dPp)时的岩石孔隙可压缩性。对于饱含流体岩石而言,围压dPc与孔隙压力dPp两者相等这种现象只存在于两种情况:第一种情况是假想岩石中的孔隙被基质充填,岩石是一个无孔的各向同性介质,对其施加围压dPc,岩石内部应力处处相等,并且等于dPc,此时
。在这里需要强调一点的是Zimmerman (1991) [8] 应用应力叠加原理推导式(9)时用了这个假设,式(9)实际为
;Nur和Byerlee [4] 在理论上推导Biot和 Willis (1957)的有效应力定理 [2] 时也用了这个假设(附录A)。第二种情况,对于孔隙饱含流体的岩石,对其施加围压dPc会产生孔隙压力dPp,因为岩石骨架要分担部分压力,也就是岩石骨架的有效应力,那么dPp必定是小于dPc,要保证dPp与dPc相等,必须要额外的增加孔隙流体压力,同时保持dPc不变,使dPp与dPc相等,此时岩石孔隙可压缩性
与第一种情况的物理模型是不一样的。对于海底疏松未固结的沉积物,孔隙完全与海水连通时,围压就是与之对应深度的海水压强Ph = ρwgh (ρw是海水密度,g是重力加速度,h是海水深度),此时海水深度增加产生的围压增量dPc与孔隙中的流体压力增量dPp是相等的,但是这种物理模型不适合固结的岩石。而对于地下的沉积岩而言,对其施加围压dPc所产生孔隙压力dPp两者是不相等的,因为岩石骨架要分担部分压力(有效应力),dPp必定是小于dPc。
因此根据以上的分析可知,
是一个具有理论意义的参数。Thomsen (2010 [21] ; 2017 [22] )讨论了
,认为Gassmann方程 [3] 应该修正为Brown 和 Korringa (1975)提出的广义Gassmann方程式(29) [16] ,但是在广义Gassmann方程式(29)中,
是一个不确定的参数,限制了其自身在实践中的应用价值。Hart和Wang (2010)根据Berea砂岩和Indiana灰岩的实验室测量结果发现,当岩石在较高的有效应力状态下(大于5 MPa),
是有效的,此时岩石呈现为线性弹性形变 [17] 。对于地下的沉积岩而言,埋藏深度超过600米,岩石的有效应力基本能达到5 MPa以上,这一结果可以根据地层的平均密度估算得到。虽然Thomsen (2010 [21] ; 2017 [22] )对Biot-Gassmann原理(Biot, 1956 [23] ; Gassmann, 1951 [3] )的有效性所持有的质疑略显犀利,但是Wang等人(2016) [19] 根据Pickett (1963) [20] 和他人的实验结果,干燥岩石的泊松比常常近似地等于矿物的泊松比,或Kdry/Kma = μdry/μma (Mavko et al., 1998, p. 244) [14] 的关系联立Biot-Gassmann原理推导出一种新的估算岩石基质体积模量Kma和干燥岩石体积模量Kdry的方法,其计算结果与理论非常吻合,这说明Biot-Gassmann原理在实践应用中具有重要的价值。笔者认为上述
是一个值得研究的参数,在这里之所以引入关于
的讨论,主要目的是力求本文所涉及到的理论来源更具有原滋原味,争取能起到抛砖引玉的作用。
5. 结论
Terzahi (1923)的有效应力定理 [1] 描述了地表饱和土壤的孔隙流体与土壤骨架(土壤骨架相当于排水状态下的多孔土壤)分担围压的力学机制,但是其应变与应力关系是模糊的,难以从土壤的三维力学模型解释其应变与应力的关系。Biot和Willis (1957) [2] 提出的有效应力定理能从理论上解释岩石围压增加的过程中,孔隙流体压力保持不变时岩石的应变与应力的关系,并且Biot和Willis (1957)提出的有效应力定理能解释Terzahi (1923)的有效应力定理成立的条件,对于地表的饱和土壤而言,土壤骨架的体积模量Kdry要远小于土壤矿物颗粒的体积模量Kma (Kdry/Kma ≈ 0),因此有效应力系数n ≈ 1,Biot和Willis (1957)的有效应力定理式(2)可简化为Terzahi(1923)的有效应力定理式(1)。
Biot和Willis (1957)提出的有效应力定理 [2] 描述了饱和多孔岩石在孔隙压力保持不变的情况下岩石应变与应力的关系,饱和岩石在围压Pc不断增加的过程中,要保持孔隙流体压力Pp不变,岩石孔隙流体必须与外界连通,排出部分流体才能保证孔隙流体压力不变(释放多余的压力)。从理论的严密性方面来讲,Biot和Willis (1957)的有效应力定理 [2] 是不足以描述孔隙不连通(非渗透)或者孔隙与外界不连通的饱和岩石的力学机制的,如果岩石孔隙不连通,意味着孔隙流体与外界不连通,当围压Pc不断增加时,孔隙内的部分流体无法排出,多余孔隙压力无法得到释放,随着围压Pc不断增加,孔隙流体压力Pp不断增大。
Gassmann方程 [3] 描述了饱和岩石在不排水的状态下(孔隙流体质量不变),围压增量dPc所引起孔隙压力增量dPp的这一力学机制,在数学上是孔隙压力Pp对围压Pc的导数。事实上,对于饱和岩石,孔隙流体与外界不连通时(孔隙流体质量不变),在围压Pc不断增加的过程中,孔隙压力Pp逐渐增大,岩石骨架(相当于排水状态下的岩石)和孔隙流体如何分担围压的问题并未得到完全地解决,Gassmann方程只是确定了孔隙压力增量dPp与围压增量dPc的关系dPp/dPc = B,而不是孔隙压力Pp与围压Pc的关系。
基金项目
中石化股份公司项目“古近系油藏地球物理识别技术研究”(P18051-2)和“储层流体的地震响应机理与油气判识”(P18051-3)联合资助。
附录A:
Nur和Byerlee (1957) [5] 对Biot和Willis (1957) [2] 提出的有效应力定理进行了理论推导,假定有一个具有任意形状、孔隙互相连通的固体材料各向同性集料,承受一个围压Pc和一个均匀的孔隙压力Pp,压缩方向为正方向。该应力状态可从概念上分两步完成:第一步,施加孔隙压力Pp和一个与之相等的围压
,即
;第二步,施加剩余的围压
,对孔隙压力不做任何进一步改变。如果现在假设集料在Pc的压力范围内呈线性弹性,那么可以将这两步导致的应变场进行叠加。
![](//html.hanspub.org/file/8-1770763x87_hanspub.png)
Figure A1. A homogeneous aggregate with pores (The lines S represent the surfaces of the pores that are subject to a pore pressure that is equal to the confining pressure)
图A1. 一个具有孔隙的各向同性集料(线S代表孔隙的表面,其承受的孔隙压力等于围压)
在分两步进行施加围压的时候,保持孔隙压力Pp不变,仅仅由围压
导致集料的体积应变
为
(A-1)
式中:β = 1/K为干燥集料的有效压缩系数。
同时假定有一个相等集料,所有的孔隙都被固体充填,并且假定内切表面为S (图A1),在多孔集料中相当于孔隙的边界。如果假定固体是各向同性的(如忽略粒间应力偏差),外部围压
在固体中导致一个相等的静水压力
。固体的体积应变θs为
(A-2)
式中:βs = 1/Ks为无孔隙固体的压缩系数(岩石基质的压缩系数)。正交于内切边界S的压力等于
。假如现在去掉内切区域内的材料,用流体替代,仍然保持孔隙压力
等于
。现在对应力边界值使用唯一性定理(Mushkelishvilli, 1963) [18] 。定理规定当法向应力被指定在所有边界上时唯一可确定的是弹性固体的形变。因为用流体代替固体夹杂物没有改变S上的应力,固体的形变仍然不变。因此,体积应变
等于
,即
(A-3)
结合以上公式,总应变为
(A-4)
或按照体积模量
(A-5)
式中:K为干燥集料的体积模量;Ks为固体本身内在的体积模量。式(A-5)支持所有呈线性弹性具有连通孔隙的各向同性集料。在上述推导过程中,饱和岩石在围压的作用下孔隙压力是保持不变,孔隙介质形变是由干燥集料的体积模量K来描述的,因此饱和岩石的总应变相当于干燥集料的总应变,干燥集料分担的是有效应力,因此有效应力Pe = Kθ,那么式(A-5)可以写成
(A-6)
式(A-6)等价于Biot和Willis (1957)提出的有效应力定理。