下载文档到电脑，查找使用更方便
20金贝
还剩页未读，继续阅读
地应力计算公式（一）、井中应力场的计算及其应用研究（秦绪英，陈有明，陆黄生 2003年6月）主应力计算根据泊松比、地层孔隙压力贡献系数、孔隙压力及密度测井值可以计算三个主应力值：相关系数计算：应用密度声波全波测井资料的纵波、横波时差（、）及测井的泥质含量可以计算泊松比、地层孔隙压力贡献系数、岩石弹性模量及岩石抗拉强度① 泊松比 ② 地层孔隙压力贡献系数 ③ 岩石弹性模量 ④ 岩石抗拉强度 注：分别为密度测井值，地层骨架密度，横波时差和纵波时差值为地区试验常数其它参数不同地区岩石抗压强度参数是参照岩石抗拉强度数值确定,一般是8～12倍,也可以通过岩心测试获得岩石内摩擦系数及岩石内聚力是岩石本身固有特性参数,可以通过测试分析获得地层孔隙压力由地层水密度针对深度积分求取,或者用重复地层测试器RFT测量也可以通过地层压裂测试获得,测试时,当井孔压力下降至不再变化时,为储层的孔隙压力二）、一种基于测井信息的山前挤压构造区地应力分析新方法（赵军 2005年4月）基于弹性力学的测井地应力分析以弹性力学理论为基础,经过一定的假设条件和边界条件可以推演出用于计算地下原地应力的数学模型,用地球物理测井信息(包括声波全波列和密度等)确定模型参数,对地应力进行连续计算与分析。
不同的研究者根据不同的条件提出了众多的地应力模型,如在油田得到较广泛应用的黄氏模型(黄荣樽等, 1995) 其模型如下：式中：分别为垂向、水平最大、水平最小地应力（单位：）；为密度（单位：）；为深度（单位：）；为有效应力系数；为孔隙压力（单位：）；为静泊松比；为地层厚度（单位：）；为重力加速度（单位：）；为构造校正量，必须分段考虑；可用依据阵列声波测井得到的纵、横波速度确定模型力学参数： 式中，为动态泊松比；为动态弹性模量(单位：)；为纵波速度（单位：）；为横波速度，（单位：）利用压实曲线计算估算附加构造应力大小基本原理：泥岩压实实验研究表明,在正常压实条件下,泥岩孔隙度随上覆压力(或埋深)呈指数递减规律(真柄钦茨,1981) ,即有:式中, 为任意深度处的泥岩孔隙度(单位: % ) ; 为泥岩初始孔隙度(单位: % ) ; 为埋藏深度(单位: m ) ; c为地区常数反映在单对数坐标图上为一条直线(图1) ,这就是通常的埋藏压实曲线注意: 这里的条件是正常压实条件, 受均匀的随深度呈线性增大的重力作用当岩石受到额外的侧向构造挤压应力的作用时, 泥岩的孔隙度会进一步减小,使得泥岩孔隙度偏离正常压实趋势线(图1) 。
这就是所谓的超压实作用这种额外的附加构造地应力可以通过偏移量的大小来估算图1中A点的孔隙度值明显偏离正常趋势线,A点在趋势线上的水平、垂直交点分别为B、C两点利用B、C两点与纵坐标的交点E、D之差(即△ED)可以估算附加构造应力大小(即深度差与岩石平均密度之积) 因此,可以通过正常趋势线方程与偏移量来求取附加构造应力电阻率、声波时差对地应力的敏感性研究发现,电阻率、声波时差对地应力的敏感性有不同特征当岩石承受的总地应力较弱时,岩石保持较高孔隙,电阻率对地应力响应不灵敏(图3) ,但声波时差能有效地反映地应力,因此只能采用声波时差建立地应力模型在岩石受到强地应力作用下,岩石致密,时波时差对地应力反映不灵敏,而电阻率能灵敏地反映地应力存在(图3) ,此时,应当采用电阻率建立地应力模型在中等地应力作用区,既可用声波时差建立模型,也可用电阻率建立模型三）、用测井资料计算地层应力（马建海 孙建孟 2002年）测井估算地层应力数值的方法测井计算的地层应力是原地层应力或扰动地层应力,从时间看主要是现今地层应力,文献[1 ,2 ,3 ]认为测井得出的是现今原地层应力计算的基本方法是首先应用密度测井积分估算出垂直应力,然后根据地层特点选择适当的模型计算水平地层应力。
1 应用密度测井估算垂直应力用密度测井资料计算垂直应力的公式为式中：：总垂直压力；：真垂直深度；：重力加速度；：偏移值；：体积密度，测量井段以上可用人工插值法获得连续的密度曲线，或借助垂直应力梯度反推 2 各种估算水平应力的模型方法各种模型基本是以垂直应力、孔隙应力和泊松比为基础，分别根据不同的理论假设来计算水平应力1） 多孔弹性水平应变模型法该模型为水平应力估算最常用的模型，它以三维弹性理论为基础最小水平主应力；：最大水平主应力；：总垂直应力；：垂直方向的有效应力系数（Biot系数）；：水平方向的有效应力系数（Biot系数）；：静态泊松比；：孔隙压力；：静态杨氏模量；：最小主应力方向的应变；：最大主应力方向的应变2） 双轴应变模型法双轴应变模型法是多孔弹性水平模型的一个特例，该特例以构造因子作输入参数，取代最大水平主应力方向的应变（）式中，为非平衡构造因子，反映的是构造力作用下最大水平应力和最小水平应力的地区经验关系3） 莫尔-库仑应力模型法此经验关系式以最大、最小主应力之间的关系给出其理论基础是莫尔-库仑破坏准则，即假设地层最大原地剪应力是由地层的抗剪强度决定的在假设地层处于剪切破坏临界状态的基础上，给出了地层应力经验关系式式中，；为三轴应力系数；为岩石内摩擦角；、为最大和最小主应力；为岩石单轴抗压强度。
当忽略地层强度时（认为破裂首先沿原有裂缝或断层发生），且垂向应力为最大主应力时，式为 进而有式中，；为岩石的内摩擦角此经验关系式有一定的物理基础,比较适合疏松砂岩地层,但其地层处于剪切破坏的临界状态的假定,没有普遍的意义该模型不考虑地层的形变机理和主应力方向,因此,它既可以用于拉张型盆地也可以用于挤压型盆地4） 一级压实模型一级压实模型通常用于表层地层,预测地层在一级压实过程中所产生的水平应力的关系（5） 单轴应变经验关系式这一类经验关系式发展最早,该经验关系式假设由于水平方向无限大,地层在沉积过程中只发生垂向变形,水平方向的变形受到限制,应变为0 ,水平方向的应力是由上覆岩层重量产生的主要有尼克经验关系式、Mattews & Kelly 经验关系式、Anderson 经验关系式、New-berry 经验关系式等近年来,有些研究者试图通过在单轴应变公式的基础上添加校正项来提高最小水平地层应力的预测精度,即式中，为Biot系数为构造应力作用的附加项，通过地层应力实测值与按上式计算得出的差来校正，且认为在一个断块内基本上为一常数，不随深度的不同而变化但实测数据来看，不同深度处是不同的6） 组合弹簧经验关系式该模型假设岩石为均质、各向同性的线弹性体,并假定在沉积和后期地质构造运动过程中,地层和地层之间不发生相对位移,所有地层两水平方向的应变均为常数。
由广义虎克定律得式中，分别为最小、最大水平主应力方向的应变,在同一断块内为常数此经验关系式把受力的地层比喻为2个平行板之间的一组弹簧,具有不同刚度的弹簧代表具有不同弹性参数的地层当两板受到力的作用时,只发生横向位移不发生偏转,从而使各弹簧的水平位移相等,刚度大的弹簧将受到较大的应力,即杨氏模量大的地层承受较高的应力该式有效地解释了砂岩地层比相邻的页岩层有更高的地层应力的现象该式假设各岩层水平方向应变相等,忽略了岩层的非线弹性特性,也没有考虑热应力7） 葛氏地层应力经验关系式葛洪魁提出了一组地层应力经验关系式,分别适用于水力压裂垂直缝和水平缝:水力压裂裂缝为垂直裂缝(最小地层应力在水平方向) 的经验关系式为式中，：热膨胀系数；：最小、最大水平地层应力方向的构造应力系数，在同一断块内可视为常数；：分别为考虑地层剥蚀的最小和最大水平地层应力附加量，在同一断块内可视为常数其中，水平应力的重力分量为，热应力分量为，构造应力分量为和,孔隙压力分量为，底层剥蚀的附加压力为和8） 应力-速度关系法X.M. Tang 等人基于理论和实验分析提出了应力与横波速度分裂的关系式,通过在实验室建立实验关系为 式中 分别为岩石受应力为0时平行和垂直于横波偏振方向的横波速度；分别为岩石受应力作用时平行和垂直于横波偏振方向的横波速度；分别为平行和垂直于横波偏振方向的应力-速度耦合系数，在实验室可有单轴应力实验测得的速度平方查核施加应力数据的线性拟合获得。
应用测井资料实际计算地层应力时,首先通过对岩石力学参数的动静态同步测试及岩石抗压和抗拉强度测试建立该地区动静态弹性参数经验关系,然后根据该地区地层应力特征选择上述模型之一应用密度和声波全波测井资料计算地层应力3 测井确定地层应力方向的方法（1）DSI（XMAC）各向异性处理结果估算地层应力方向在成岩期和成岩后,如果水平应力存在着较大的各向异性,岩石会表现出侧向差异压实现象此时,最大水平主应力方向上侧向压实程度较高,而在最小水平应力方向上侧向压实程度较低,从而造成了应力引起的岩石物理各向异性C. Esmersoy 等人研究表明在最大水平主应力方向上的横波传播速度大于最小水平主应力的横波传播速度具体方法见文章）T.J . Plona 等人还提出了应用频散特性识别岩石内在各向异性和应力引起各向异性,指出声波激发频率与测得的快慢横波时差曲线关系图上,当岩石为内在各向异性时,2 条频散曲线平行;当岩石为应力各向异性时2 条频散曲线交叉另外,岩石内在各向异性也可由地层微电阻率扫描成像、声波反射成像等直接识别2）双井径曲线估算地层应力方向理论和实验表明钻井过程中应力崩塌形成的椭圆井眼通常是由切向正应力作用于井壁形成的,椭圆井眼的长轴方向为最小应力方向,井眼表面上有拉应力,径向拉伸破坏岩石,造成在最大水平主应力方向上形成钻井诱导缝。
在实际应用时,需排除高士钧、储昭坦提出的非地层应力因素引起的椭圆井眼: (1) 溶蚀崩落变形井眼常发生于盐膏岩层,它是因岩盐、膏盐等岩层被钻井液溶蚀所形成,其基本形状一般为圆形,双井径读数均大于钻头直径2) 浸蚀崩落变形井眼,井壁周边岩石在经过长时间钻井液浸泡后,一些较软的岩石因吸收水份而使内部结构发生膨胀,强度降低,以致引起崩落由于岩石本身具有各向异性,这种崩落在井眼周边是不均匀的,也往往造成椭圆形井眼,在双井径曲线上表现为两井径读数不等,但都大于钻头直径3) 键槽形变形井眼,由于钻具偏心对井壁一侧反复碰撞磨损造成,多发生于井斜较大,岩石强度较低的井段,其特征为非对称椭圆井眼,在双井径曲线上表现为一条井径读数大于钻头直径,而另一条小于钻头直径4) 岩石弹塑变形井眼,有些柔性地层岩石具有弹塑流变特征,在水平压应力作用下发生缩径现象,形成对称椭圆井眼,两条井径读数均小于钻头直径,长轴指向最大水平主应力方向5) 高角度裂缝崩落变形井眼,一些与井壁相交的高角度裂缝,造成井壁邻近地层岩块强度降低,经过泥浆浸泡、冲刷及钻具反复碰撞振动,可能造成沿裂缝走向的井壁岩块崩落,形成椭圆井眼,在双井径曲线上表现为一条井径读数大于钻头直径,另一条等于钻头直径。
容易与地应力造成的椭圆井眼相混6) 井斜大造成的视椭圆形井眼:井眼并未变形,只是由于井斜大,地层倾角测量时仪器偏心,从而出现一条井径测量读数等于钻头直径,而另一条小于钻头直径四）用测井资料计算最大和最小水平应力剖面的新方法（谢刚2005年2月） （利用的是成像测井）具体见文章利用测井资料计算地层应力的大小和方向,目前国内多采用单轴应变模式确定地层应力,由此可以得到最小水平地层应力剖面,但不能得到最大水平地层应力剖面,因而无法计算地层破裂。