地应力计算公式概要

地应力计算公式（一）、井中应力场的计算及其应用研究（秦绪英，陈有明，陆黄生2003年6月）主应力计算根据泊松比、地层孔隙压力贡献系数V、孔隙压力0P及密度测井值b可以计算三个主应力值：001HvAVPVP001hvBVPVP0Hvbdh相关系数计算：应用密度声波全波测井资料的纵波、横波时差（pt、st）及测井的泥质含量shV可以计算泊松比、地层孔隙压力贡献系数V、岩石弹性模量E及岩石抗拉强度TS。①泊松比22220.52()spsptttt②地层孔隙压力贡献系数22222(34)12()bsspmmsmptttVtt③岩石弹性模量2222234spbsspttEttt④岩石抗拉强度22(34)[(1)]TbspshshSattbEVcEV注：,,,mmsmptt分别为密度测井值，地层骨架密度，横波时差和纵波时差值。,,abc为地区试验常数。其它参数不同地区岩石抗压强度参数是参照岩石抗拉强度数值确定,一般是8～12倍,也可以通过岩心测试获得。岩石内摩擦系数及岩石内聚力是岩石本身固有特性参数,可以通过测试分析获得。地层孔隙压力由地层水密度针对深度积分求取,或者用重复地层测试器RFT测量。也可以通过地层压裂测试获得,测试时,当井孔压力下降至不再变化时,为储层的孔隙压力。（二）、一种基于测井信息的山前挤压构造区地应力分析新方法（赵军2005年4月）基于弹性力学的测井地应力分析以弹性力学理论为基础,经过一定的假设条件和边界条件可以推演出用于计算地下原地应力的数学模型,用地球物理测井信息(包括声波全波列和密度等)确定模型参数,对地应力进行连续计算与分析。不同的研究者根据不同的条件提出了众多的地应力模型,如在油田得到较广泛应用的黄氏模型(黄荣樽等,1995)。其模型如下：0()Hvhgdh1()1sHvppsustraPaPu2()1shvppsustraPaPu式中：,,vHh分别为垂向、水平最大、水平最小地应力（单位：MPa）；为密度（单位：3/gcm）；H为深度（单位：m）；a为有效应力系数；pP为孔隙压力（单位：MPa）；su为静泊松比；h为地层厚度（单位：m）；g为重力加速度（单位：2/Ns）；12,strstr为构造校正量，必须分段考虑；可用依据阵列声波测井得到的纵、横波速度确定模型力学参数：220.5(/)1(/)1psdpsVVUVV222322(34)10spsdpsVVVEVV式中，dU为动态泊松比；dE为动态弹性模量(单位：MPa)；pV为纵波速度（单位：/ms）；sV为横波速度，（单位：/ms）。利用压实曲线计算估算附加构造应力大小基本原理：泥岩压实实验研究表明,在正常压实条件下,泥岩孔隙度随上覆压力(或埋深)呈指数递减规律(真柄钦茨,1981),即有:0cZe式中,为任意深度处的泥岩孔隙度(单位:%);0为泥岩初始孔隙度(单位:%);Z为埋藏深度(单位:m);c为地区常数。反映在单对数坐标图上为一条直线(图1),这就是通常的埋藏压实曲线。注意:这里的条件是正常压实条件,受均匀的随深度呈线性增大的重力作用。当岩石受到额外的侧向构造挤压应力的作用时,泥岩的孔隙度会进一步减小,使得泥岩孔隙度偏离正常压实趋势线(图1)。这就是所谓的超压实作用。这种额外的附加构造地应力可以通过偏移量的大小来估算。图1中A点的孔隙度值明显偏离正常趋势线,A点在趋势线上的水平、垂直交点分别为B、C两点。利用B、C两点与纵坐标的交点E、D之差(即△ED)可以估算附加构造应力大小(即深度差与岩石平均密度之积)。因此,可以通过正常趋势线方程与偏移量来求取附加构造应力。电阻率、声波时差对地应力的敏感性研究发现,电阻率、声波时差对地应力的敏感性有不同特征。当岩石承受的总地应力较弱时,岩石保持较高孔隙,电阻率对地应力响应不灵敏(图3),但声波时差能有效地反映地应力,因此只能采用声波时差建立地应力模型。在岩石受到强地应力作用下,岩石致密,时波时差对地应力反映不灵敏,而电阻率能灵敏地反映地应力存在(图3),此时,应当采用电阻率建立地应力模型。在中等地应力作用区,既可用声波时差建立模型,也可用电阻率建立模型。（三）、用测井资料计算地层应力（马建海孙建孟2002年）测井估算地层应力数值的方法测井计算的地层应力是原地层应力或扰动地层应力,从时间看主要是现今地层应力,文献[1,2,3]认为测井得出的是现今原地层应力。计算的基本方法是首先应用密度测井积分估算出垂直应力,然后根据地层特点选择适当的模型计算水平地层应力。1应用密度测井估算垂直应力用密度测井资料计算垂直应力的公式为0()TVDvbghdhO式中：v：总垂直压力；TVD：真垂直深度；g：重力加速度；O：偏移值；b：体积密度，测量井段以上可用人工插值法获得连续的密度曲线，或借助垂直应力梯度反推。2各种估算水平应力的模型方法各种模型基本是以垂直应力、孔隙应力和泊松比为基础，分别根据不同的理论假设来计算水平应力。（1）多孔弹性水平应变模型法该模型为水平应力估算最常用的模型，它以三维弹性理论为基础。221111hvvertphorphHEEPP221111HvvertphorpHhEEPPh：最小水平主应力；H：最大水平主应力；v：总垂直应力；vert：垂直方向的有效应力系数（Biot系数）；hor：水平方向的有效应力系数（Biot系数）；：静态泊松比；pP：孔隙压力；E：静态杨氏模量；h：最小主应力方向的应变；H：最大主应力方向的应变。（2）双轴应变模型法双轴应变模型法是多孔弹性水平模型的一个特例，该特例以构造因子作输入参数，取代最大水平主应力方向的应变（H）()111hvvertphorphhhEPPKKHhhK式中，hK为非平衡构造因子，反映的是构造力作用下最大水平应力和最小水平应力的地区经验关系。（3）莫尔-库仑应力模型法此经验关系式以最大、最小主应力之间的关系给出。其理论基础是莫尔-库仑破坏准则，即假设地层最大原地剪应力是由地层的抗剪强度决定的。在假设地层处于剪切破坏临界状态的基础上，给出了地层应力经验关系式103()/ppPCPN式中，2(/4/2)Ntg；N为三轴应力系数；为岩石内摩擦角；1、3为最大和最小主应力；0C为岩石单轴抗压强度。当忽略地层强度0C时（认为破裂首先沿原有裂缝或断层发生），且垂向应力为最大主应力时，式103()/ppPCPN为13()/ppPPN进而有22111hvpPtgtgHhhK式中，/4/2；为岩石的内摩擦角。此经验关系式有一定的物理基础,比较适合疏松砂岩地层,但其地层处于剪切破坏的临界状态的假定,没有普遍的意义。该模型不考虑地层的形变机理和主应力方向,因此,它既可以用于拉张型盆地也可以用于挤压型盆地。（4）一级压实模型一级压实模型通常用于表层地层,预测地层在一级压实过程中所产生的水平应力的关系(1sin)hvHhhK（5）单轴应变经验关系式这一类经验关系式发展最早,该经验关系式假设由于水平方向无限大,地层在沉积过程中只发生垂向变形,水平方向的变形受到限制,应变为0,水平方向的应力是由上覆岩层重量产生的。主要有尼克经验关系式、Mattews&Kelly经验关系式、Anderson经验关系式、New-berry经验关系式等。近年来,有些研究者试图通过在单轴应变公式的基础上添加校正项来提高最小水平地层应力的预测精度,即()1hpvpTPP式中，为Biot系数。T为构造应力作用的附加项，通过地层应力实测值与按上式计算得出的差来校正，且认为在一个断块内T基本上为一常数，不随深度的不同而变化。但实测数据来看，不同深度处T是不同的。（6）组合弹簧经验关系式该模型假设岩石为均质、各向同性的线弹性体,并假定在沉积和后期地质构造运动过程中,地层和地层之间不发生相对位移,所有地层两水平方向的应变均为常数。由广义虎克定律得22()111hHhpvpEEPP22()111hHHpvpEEPP式中，,hH分别为最小、最大水平主应力方向的应变,在同一断块内为常数。此经验关系式把受力的地层比喻为2个平行板之间的一组弹簧,具有不同刚度的弹簧代表具有不同弹性参数的地层。当两板受到力的作用时,只发生横向位移不发生偏转,从而使各弹簧的水平位移相等,刚度大的弹簧将受到较大的应力,即杨氏模量大的地层承受较高的应力。该式有效地解释了砂岩地层比相邻的页岩层有更高的地层应力的现象。该式假设各岩层水平方向应变相等,忽略了岩层的非线弹性特性,也没有考虑热应力。（7）葛氏地层应力经验关系式葛洪魁提出了一组地层应力经验关系式,分别适用于水力压裂垂直缝和水平缝:水力压裂裂缝为垂直裂缝(最小地层应力在水平方向)的经验关系式为()()111vpThvphphEPETPKP()()111vpTHvpHpHEPETPKP式中，T：热膨胀系数；,hHKK：最小、最大水平地层应力方向的构造应力系数，在同一断块内可视为常数；,hH：分别为考虑地层剥蚀的最小和最大水平地层应力附加量，在同一断块内可视为常数。其中，水平应力的重力分量为1v，热应力分量为1TET，构造应力分量为()1vphEPK和()1vpHEPK,孔隙压力分量为121pP，底层剥蚀的附加压力为h和H。（8）应力-速度关系法X.M.Tang等人基于理论和实验分析提出了应力与横波速度分裂的关系式,通过在实验室建立实验关系为2222,,0,,0shshhhsHsHHHvvSvvS式中,0,0,shsHvv分别为岩石受应力为0时平行和垂直于横波偏振方向的横波速度；,,,shsHvv分别为岩石受应力作用时平行和垂直于横波偏振方向的横波速度；,hHSS分别为平行和垂直于横波偏振方向的应力-速度耦合系数，在实验室可有单轴应力实验测得的速度平方查核施加应力数据的线性拟合获得。应用测井资料实际计算地层应力时,首先通过对岩石力学参数的动静态同步测试及岩石抗压和抗拉强度测试建立该地区动静态弹性参数经验关系,然后根据该地区地层应力特征选择上述模型之一应用密度和声波全波测井资料计算地层应力。3测井确定地层应力方向的方法（1）DSI（XMAC）各向异性处理结果估算地层应力方向在成岩期和成岩后,如果水平应力存在着较大的各向异性,岩石会表现出侧向差异压实现象。此时,最大水平主应力方向上侧向压实程度较高,而在最小水平应力方向上侧向压实程度较低,从而造成了应力引起的岩石物理各向异性。C.Esmersoy等人研究表明在最大水平主应力方向上的横波传播速度大于最小水平主应力的横波传播速度。（具体方法见文章）T.J.Plona等人还提出了应用频散特性识别岩石内在各向异性和应力引起各向异性,指出声波激发频率与测得的快慢横波时差曲线关系图上,当岩石为内在各向异性时,2条频散曲线平行;当岩石为应力各向异性时2条频散曲线交叉。另外,岩石内在各向异性也可由地层微电阻率扫描成像、声波反射成像等直接识别。（2）双井径曲线估算地层应力方向理论和实验表明钻井过程中应力崩塌形成的椭圆井眼通常是由切向正应力作用于井壁形成的,椭圆井眼的长轴方向为最小应力方向