多孔介质有效应力的分形模型及计算.pdf

多孔介质有效应力的分形模型及计算 李玮 闫铁 苏鹏 孙维国 邢云龙 （大庆石油学院 黑龙江大庆 163318） 摘要：本文就多孔介质应力关系展开研究。在分析了有效应力方程，尤其是泰尔扎吉方程存在的问题之后,根据分形方 法从二维到三维推导了更为合理的多孔介质应力关系的分形方程，该方程能反应任意孔隙结构下应力作用关系。为了解 决现场应用方便，并进一步给出多孔介质应力关系统计模型。新方程的特点是引入了孔隙度和分形维数两个能反映多孔 介质结构特性的重要参数，使有效应力方程更加完善。实验数据计算结果表明：新的基于分形的有效应力模型与原方程 的计算结果存在着巨大差异。 关键词： 多孔介质；有效应力；Terzaghi 方程；分形几何；分形维数 岩石是一种矿物组成上非均匀的，孔隙 维的地层内部应力处理方法后来被Biot所 结构上复杂多变的多孔介质。其结构上是由 总结，并用相容性理论阐明了扩散与变形耦 相互连接在一起的固体颖粒(骨架)所组成 合作用的一体过程[2,3] 。目前式（1）已经广 的，孔隙中通常被一种或几种流体填充。通 泛应用于地基处理、岩土工程、地质、石油 常情况下，多孔介质同时受到外部应力和内 工程等许多领域。 部应力的共同作用。Terzaghi对深入地分析 1959 年Habbert 等人[4]借鉴了泰尔扎 多孔介质上的总应力、有效应力、孔隙压力， 吉方程原理。认为上覆岩层总体重力所形成 并提出著名的“有效应力原理”。该方程阐 的压力σ 即为总应力σ ，它由沉积物颗 OV O 述了地下储层岩石骨架及孔隙组成、骨架起 粒之间接触点上的应力σ （基岩应力）与孔 到支撑作用。目前已经广泛应用于地基处 隙流体压力P 共同支撑，即 p 理、岩土工程、地质、石油工程等许多领域 σ σ+P （2） OV p [1-3] 。 式中：σ 上覆岩层压力，MPa。在简单修 OV 本文应用分形理论，从多孔介质结构上 改变量的概念后，式（2）通过形式上基本 的随机性、非均匀性入手，重新建立了有效 照搬式（1），但内涵有所不同。由于方程式 应力的分形计算模型，并进行了相应的参数 （2）中只有一个未知量且和上覆岩层压力 分析和简化计算，为进一步应用分形方法研 联系起来，因此在石油工程中得到了广泛的 究多孔介质的应力关系奠定基础。 应用。 σ 和P 两个量 1、有效应力的基本模型 由以上两式的结构可知， P 均衡地分担总应力的作用，而忽略了岩石本 1923年泰尔扎吉在研究浅层岩土多孔 体结构特征在σ 和P 之间协调所起的作 P 介质时，首先引入了一维压实有效应力的概 用。 [1] 念，并给出如下关系 [5] 1963年Hand