文档介绍：岩石物理参数计算及应力研究-llzlllo岩石物理参数计算及应力研究岩石物理参数计算地层岩石是地应力的载体,岩石物理性质对地应力的传递、衰减、集中、分散都会产生很大的影响,岩石物理参数与岩体赋存的地应力密切相关,岩石物理参数计算是地应力研究的必然步骤。通过纵、横波时差和密度等测井资料,能够计算地层条件下的岩石动态弹性模量,在此基础上,能够进行地应力分析、井眼稳定性分析、地层出砂分析、以及人工压裂设计等方面的研究。岩石物理参数包括岩石弹性参数和岩石机械强度参数。岩石弹性参数主要有泊松比、杨氏模量、剪切模量、体积模量、体压缩系数和、有效应力系数系数(比奥特系数);岩石机械强度主要有单轴抗压强度、岩石的抗剪强度和岩石抗张强度,以及内摩擦角等。1、岩石弹性参数对于各向同性均匀介质岩石来说,利用牛顿第二定律和三维虎克定律,经数学推导,可导出计算声波速度在岩石介质中的波动方程:(3-1-1)(3-1-2)根据上述的波动方程,能够得出各种弹性参数与声波时差的关系式。①泊松比定义为横向应变与纵向应变之比。(3-1-3)②切变模量定义为施加的应力与切应变之比。(3-1-4)③杨氏模量定义为施加的轴向应力与法向应变之比。(3-1-5)④体积模量定义为静水压力与体积应变之比。(3-1-6)⑤体积压缩系数定义为体积模量的倒数。即:(3-1-7)⑥有效应力系数(Boit)表示孔隙压力对岩石变形的影响,即:=Kb/Kma(3-1-8)式中:为岩石体积密度,;、为纵、横波时差,。公式中的a为单位转换系数。具体选择方法如下:如果密度单位为,时差单位为,弹性参数单位为psi,则;如果密度单位为,时差单位为,弹性参数单位为,则;如果密度单位为,时差单位为,弹性参数单位为,则。因此,利用阵列声波测井提供的纵、横波时差以及常规测井提供的密度资料就能够进行岩石弹性参数计算。但是由于费用等原因,并不是每口井都开展声波全波列或阵列声波测井,因而不能直接获取横波时差资料,在研究中则能够通过构造内某些井已有的横波时差曲线资料来建立横波时差曲线计算式。研究表明,横波时差与纵波时差、地层密度和纵波波阻抗之间有很好的相关性。通过对安棚地区4口井的纵横波时差曲线进行分析后,建立了纵横波时差经验关系式:(3-1-9)图3-1是由上式纵横波时差关系式得出的横波时差与实测横波时差的关系图,从图中能够看出,大部分点分布在斜率约为450的直线上,计算的横波时差与实测横波时差近似相等。图3-1-1合成横波时差与实测横波时差关系图当研究区内没有一口井具有横波时差资料时,则可用下面的公式来合成横波时差曲线:(3-1-10)式中:、为岩石骨架的横波时差和纵波时差,;、为流体的横波时差和纵波时差,。2、岩石动、静态弹性参数之间的转换方法岩石弹性参数的常用测定方法有动态法和静态法两种。静态法是通过对岩样进行静态加载测其变形得到,所得弹性参数称之为静态参数;动态法则是通过测定超声波在岩样中的传播速度转换得到,所得弹性参数称之为动态参数。因此,用测井资料计算得到的弹性参数是动态参数。根据地下岩层的应力形成、赋存和起作用的机理,特别是在应力幅值、加载速度和所引起的岩石变形等方面,更接近岩石静态测试的条件,另外,现有的力学本构关系一般是基于静态参数建立的,因此,在地应力计算和实际工程中应采用岩石的静态弹性参数。大量研究资料表明岩石的动态、静态弹性参数具有很好的相关性,且大部分情况下岩样的静态参数弹性模量小于其动态值。岩石动、静态弹性参数间存在较大差别,其原因主要是岩石中微裂缝和孔隙的存在。岩石这种孔隙的弹性材料有别于各向同性、均质的线弹性体。微裂缝的存在对岩石静态变形的影响较大,而超声波能够绕过一些微裂缝传播。在实际应用时,可通过岩石力学动、静态同步测试建立动、静态参数间的关系,从而把测井得到的动态参数转换为静态参数。由于研究区及其邻近区块没有条件做岩石力学试验,本次研究引用了辽河油田和大庆油田的实验结果:(3-1-11)(3-1-12)式中:、为静态杨氏模量和静态泊松比;、为动态杨氏模量和动态泊松比,即测井资料计算结果。式(3-11)和式(3-12)。3、岩石机械强度参数当前,岩石机械强度参数还没有理论计算式,一般通过岩石力学测试来确定。为了克服岩石力学试验存在的测试费用昂贵和数据量少等缺点,研究人员通过岩石力学试验建立了岩石强度参数的经验计算式:⑴单轴抗压强度Deer和Miller(1996)根据大量的室内试验结果建立了砂泥岩的单轴抗压强度与岩石动态杨氏模量和地层泥质含量之间的关系:①砂泥岩地层(3-1-13)②碳酸盐岩地层(3-1-14)式中:和的单位为;泥质含量为小数,由自然伽马测井资料确定。⑵岩石粘聚力(内聚强度或抗剪切强度)根据Brules和Coates的研究结果,粘聚力的计算公式为:(3-1-15)