地层破裂压力和地层坍塌压力预测新算法

1、地层破裂压力和地层坍塌压力预测新算法地层岩石作为一种多孔两相固体物质，其应力分析与普通单相固体物质是有区别的，但是，在我们目前所使用的地层岩石应力分析模型、理论中，都有意或无意地使用了单相固体应力分析的方法。为了分析两者的区别，在这里我们首先引入有效应力的概念。有效应力的概念是由李传亮老师首先提出来的，该理论认为岩石由两个有效应力：本体有效应力和结构有效应力。本体有效应力决定岩石的本体变形，结构有效应力决定岩石的结构变形。 (1) （2） (3) （4）式中：上覆地层压力；岩石骨架应力；岩石接触应力；岩石本体有效应力；岩石结构有效应力；岩石孔隙度；岩石触点孔隙度；（=）Pp岩石空隙流体压力。有

2、效应力通过孔隙度把普通材料和多孔介质统一起来了，有效应力计算公式中的孔隙度反映了孔隙压力对有效应力的贡献权值。在地应力分析中，我们所指的应力是结构有效应力。（1）借助结构有效应力公式，我们首先分析在非均匀地应力作用下井眼周围周向结构有效应力和径向结构有效应力分布规律。（5）式中：距井轴r距离并与按逆时针方向成角处的周向结构有效应力。 （6） （7）岩石泊松系数；A，B构造应力系数（构造应力系数对于不同的地质构造是不同的，但在统一构造断块内部，它是一个常数，且不随地层深度变化）；地层孔隙流体压力；上覆地层压力。需要说明的是，上述式（6）和式（7）式与现用的地应力计算公式（如下）是有区别的： （8

3、） （9）式中：有效应力系数；根据李传亮老师提出的岩石变形机制理论，岩石作为一种多孔胶结固体，其变形实际上是通过骨架颗粒的变形和排列方式的改变来实现的，骨架颗粒的变形即为岩石的本体变形，排列方式的改变即为岩石的结构变形。本体变形对应岩石的压缩变形过程，其孔隙度保持恒定不变，结构变形则对应岩石的压实变形过程，其孔隙度要发生变化。式（8）和式（9）中的，、分别为岩石骨架压缩率和岩石容积压缩率，通过实验模拟地层岩石应力状态测得，并认为就是地层岩石在地层中的孔隙度，其实这是没道理的，无论在地面还是在地下岩石的孔隙度是不变的，因为岩石仅仅发生了本体变形而没有发生结构变形。假设：=0.2，H=5000m，

4、A=1，B=0.5，地层流体压力梯度=12Mpa/km，上覆地层压力梯度=23MPA/km，=0.15.则在非均匀地应力作用下井眼周围周向结构有效应力和径向结构有效应力分布规律如下图。（2）我们再分析在井筒内液柱压力作用下井筒周围周向结构有效应力和径向结构有效应力分布规律。通过对厚壁圆筒在内压作用下周向应力和径向应力分布规律的分析： （10） 式中：，分别为厚壁圆筒所受内压和外压；，分别为厚壁圆筒内径和外径。从应力分布图（上图）我们可以看到，内压引起周围周向和径向应力变化主要集中在内壁附近，所以对于半径趋向于无穷大的井筒模型，可以忽略外边界应力变化，认为其周向应力和径向应力分布规律与厚壁圆筒受

5、内压情形相同。（3）对上述两种作用的周向结构有效应力分布进行叠加我们可以得到：（井筒液柱压力） （11）(式（10）在趋向于无穷大时可简化得到。)同样假设：=0.2，H=5000m，A=1，B=0.5，地层流体压力梯度=12Mpa/km，上覆地层压力梯度=23MPA/km，泥浆附加密度0.1g/cm3 ,=0.15.则在钻进过程中井眼周围周向应力和径向应力分布规律如下图。在=0°和180°方向，井壁上最大周向结构有效应力（拉正压负）： （12）在=90°和270°方向，井壁上最小周向结构有效应力（拉正压负）： （13）（4）再考虑井筒流体渗入地层在井壁上

6、产生的附加周向应力（压裂作业和钻井作业中泥饼形成之前）： （14）则在=0°和180°方向，井壁上最大周向结构有效应力（拉正压负）： （15）在=90°和270°方向，井壁上最小周向结构有效应力（拉正压负）： （16）（5）破裂压力和坍塌压力A破裂压力井壁上总的最大周向结构有效应力：当结构有效应力达到岩石的抗拉应力强度时，岩石即破裂，此时 （17）于是，结合式（17）和式（18）可得岩石破裂压力计算公式： （18）式中：。B坍塌压力由库尔摩尔强度准则可知，岩石破坏的条件是： （19）由上述井壁处周向结构有效应力分析可知，在=90°和270°方向，井壁上最小周向结构有效应力（拉正压负）： 即：在在=90°和270°方向，井壁上周向结构有效应力为最大正应力（取上式的负值），井筒液柱压力为最小正应力。 （20）经过瞬时滤失后形成优质泥饼，近似认为泥饼渗透率为零，则： （21）将（20）（21）式代入（19）式可以得到： （22）总结：结构有效应力的引入，使得对岩石这种多空介质的研究变得清晰明朗，通过以上的分析推倒，我们更正了最大水平地应力