地应力计算公式.docx

2003年6月）地应力计算公式 2003年6月） （一）、井中应力场的计算及其应用研究（秦绪英，陈有明，陆黄生 主应力计算 根据泊松比J、地层孔隙压力贡献系数 V、孔隙压力P0及密度测井值 心可以计算三个 主应力值: -VR VR -VR VR -? 1 6 B - v -VF0 VPo H ?二订 dh 0 相关系数计算： 应用密度声波全波测井资料的纵波、 横波时差（.Vtp、?迸s）及测井的泥质含量 Vsh可以 计算泊松比 」、地层孔隙压力贡献系数 V、岩石弹性模量 E及岩石抗拉强度 ST。 ①泊松比0?砂2-比 ①泊松比 2（疋-就） ②地层孔隙压力贡献系数 、，PbAgt—Atp） V =1 2 2 2咋幽伫—^p） ③岩石弹性模量 厂Pb 3尿-4比 ④岩石抗拉强度 St 二 a 订（3 ￡ 一4弋）[b 厂（1-冷）c 厂 V^] 注：;?, ;?m^tms^tmp分别为密度测井值，地层骨架密度，横波时差和纵波时差值。 a,b, C为 地区试验常数。 其它参数 不同地区岩石抗压强度参数是参照岩石抗拉强度数值确定 ，一般是8?12倍，也可以通过 岩心测试获得。岩石内摩擦系数及岩石内聚力是岩石本身固有特性参数 ，可以通过测试分析 获得。地层孔隙压力由地层水密度针对深度积分求取 ，或者用重复地层测试器 RFT测量。也 可以通过地层压裂测试获得 ，测试时，当井孔压力下降至不再变化时 ，为储层的孔隙压力。 （二）、一种基于测井信息的山前挤压构造区地应力分析新方法（赵军 2005年 2005年4月) 以弹性力学理论为基础，经过一定的假设条件和边界条件可以推演出用于计算地下原地 应力的数学模型，用地球物理测井信息（包括声波全波列和密度等）确定模型参数，对地应力进 行连续计算与分析。不同的研究者根据不同的条件提出了众多的地应力模型 ，如在油田得到 较广泛应用的黄氏模型（黄荣樽等，1995）。其模型如下： H J 二」(h)gdh -HUs+ strjJ -H Us + strj J _Us v - aPp) aPp 1 -Us str2 (j-aPp)+aPp 式中：v^H^h分别为垂向、水平最大、水平最小地应力（单位： MPa）； P为密度（单 位：g / cm3）； H为深度（单位：m）； a为有效应力系数；Pp为孔隙压力（单位：MPa）； 2 Us为静泊松比；h为地层厚度（单位： m）； g为重力加速度（单位： N/s ）； stnstD为 构造校正量，必须分段考虑； 可用依据阵列声波测井得到的纵、横波速度确定模型力学参数: 0.5(Vp/Vs)2 -1 (Vp/Vs)2 -1 Ed「 Ed 「Vs2(3Vp2-4Vs2) 103 Vp2 -Vs2 式中，Ud为动态泊松比；Ed为动态弹性模量（单位：MPa）； Vp为纵波速度（单位：m/s）； Vs为横波速度，（单位：m/s ）。 利用压实曲线计算估算附加构造应力大小 基本原理：泥岩压实实验研究表明，在正常压实条件下，泥岩孔隙度随上覆压力（或埋深） 呈指数递减规律（真柄钦茨,1981）,即有： 式中，??为任意深度处的泥岩孔隙度 伸位：% ） ; 0为泥岩初始孔隙度（单位：% ） ; Z 为埋藏深度（单位:m ） ; c为地区常数。反映在单对数坐标图上为一条直线 （图1），这就是通常 的埋藏压实曲线。注意：这里的条件是正常压实条件，受均匀的随深度呈线性增大的重力作 用。当岩石受到额外的侧向构造挤压应力的作用时 ，泥岩的孔隙度会进一步减小 ，使得泥岩孔 隙度偏离正常压实趋势线（图1）。这就是所谓的超压实作用。这种额外的附加构造地应力可 以通过偏移量的大小来估算。图 1中A点的孔隙度值明显偏离正常趋势线 ,A点在趋势线上的 水平、垂直交点分别为 B、C两点。利用B、C两点与纵坐标的交点 E、D之差（即厶ED）可以估 算附加构造应力大小（即深度差与岩石平均密度之积）。因此，可以通过正常趋势线方程与偏 移量来求取附加构造应力。 电阻率、声波时差对地应力的敏感性 研究发现，电阻率、声波时差对地应力的敏感性有不同特征。当岩石承受的总地应力较 弱时，岩石保持较高孔隙，电阻率对地应力响应不灵敏 （图3），但声波时差能有效地反映地应力 因此只能采用声波时差建立地应力模型。 在岩石受到强地应力作用下，岩石致密，时波时差对地应力反映不灵敏，而电阻率能灵敏 地反映地应力存在（图3），此时，应当采用电阻率建立地应力模型。在中等地应力作用区 ，既可 用声波时差建立模型，也可用电阻率建立模型。  （三）、用测井资料计算地层应力（马建海 孙建孟 2002年） 测井估算地层应力数值的方法 测井计算的地层应力是原地层应力或扰动地层应力 ，从时间看主要是现今地层应力，文献 [1 ,2 ,3 ]认为测井得出