《岩石物理参数计算及应力研究llzlllo》由会员分享，可在线阅读，更多相关《岩石物理参数计算及应力研究llzlllo（37页珍藏版）》请在金锄头文库上搜索。
1、第三章 岩石物理参数计算及应力研究第一节 岩石物理参数计算地层岩石是地应力的载体，岩石物理性质对地应力的传递、衰减、集中、分散都会产生很大的影响，岩石物理参数与岩体赋存的地应力密切有关，岩石物理参数计算是地应力研究的必然环节。通过纵、横波时差和密度等测井资料，可以计算地层条件下的岩石动态弹性模量，在此基本上，可以进行地应力分析、井眼稳定性分析、地层出砂分析、以及人工压裂设计等方面的研究。岩石物理参数涉及岩石弹性参数和岩石机械强度参数。岩石弹性参数重要有泊松比、杨氏模量、剪切模量、体积模量、体压缩系数和、有效应力系数系数（比奥特系数）；岩石机械强度重要有单轴抗压强度、岩石的抗剪强度和岩石抗张强度
2、，以及内摩擦角等。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)式中：为
3、岩石体积密度，；、为纵、横波时差，。公式中的a为单位转换系数。具体选择措施如下：如果密度单位为，时差单位为，弹性参数单位为psi，则；如果密度单位为，时差单位为，弹性参数单位为，则；如果密度单位为，时差单位为，弹性参数单位为，则。因此，运用阵列声波测井提供的纵、横波时差以及常规测井提供的密度资料就可以进行岩石弹性参数计算。但是由于费用等因素,并不是每口井都开展声波全波列或阵列声波测井,因而不能直接获取横波时差资料，在研究中则可以通过构造内某些井已有的横波时差曲线资料来建立横波时差曲线计算式。研究表白，横波时差与纵波时差、地层密度和纵波波阻抗之间有较好的有关性。通过对安棚地区4口井的纵横波时差曲
4、线进行分析后，建立了纵横波时差经验关系式：(3-1-9)图3-1是由上式纵横波时差关系式得出的横波时差与实测横波时差的关系图，从图中可以看出，大部分点分布在斜率约为450的直线上，计算的横波时差与实测横波时差近似相等。图3-1-1 合成横波时差与实测横波时差关系图当研究区内没有一口井具有横波时差资料时，则可用下面的公式来合成横波时差曲线：(3-1-10)式中：、为岩石骨架的横波时差和纵波时差，；、为流体的横波时差和纵波时差，。2、岩石动、静态弹性参数之间的转换措施岩石弹性参数的常用测定措施有动态法和静态法两种。静态法是通过对岩样进行静态加载测其变形得到，所得弹性参数称之为静态参数；动态法则是通
5、过测定超声波在岩样中的传播速度转换得到，所得弹性参数称之为动态参数。因此，用测井资料计算得到的弹性参数是动态参数。根据地下岩层的应力形成、赋存和起作用的机理，特别是在应力幅值、加载速度和所引起的岩石变形等方面，更接近岩石静态测试的条件，此外，既有的力学本构关系一般是基于静态参数建立的，因此，在地应力计算和实际工程中应采用岩石的静态弹性参数。大量研究资料表白岩石的动态、静态弹性参数具有较好的有关性，且大部分状况下岩样的静态参数弹性模量不不小于其动态值。岩石动、静态弹性参数间存在较大差别，其因素重要是岩石中微裂缝和孔隙的存在。岩石这种孔隙的弹性材料有别于各向同性、均质的线弹性体。微裂缝的存在对岩石
6、静态变形的影响较大，而超声波可以绕过某些微裂缝传播。在实际应用时，可通过岩石力学动、静态同步测试建立动、静态参数间的关系，从而把测井得到的动态参数转换为静态参数。由于研究区及其邻近区块没有条件做岩石力学实验，本次研究引用了辽河油田和大庆油田的实验成果：(3-1-11)(3-1-12)式中：、为静态杨氏模量和静态泊松比；、为动态杨氏模量和动态泊松比，即测井资料计算成果。式(3-11)和式(3-12)的有关系数分别为0.75和0.86。3、岩石机械强度参数 目前，岩石机械强度参数还没有理论计算式，一般通过岩石力学测试来拟定。为了克服岩石力学实验存在的测试费用昂贵和数据量少等缺陷，研究人员通过岩石力
7、学实验建立了岩石强度参数的经验计算式：单轴抗压强度Deer和Miller（1996）根据大量的室内实验成果建立了砂泥岩的单轴抗压强度与岩石动态杨氏模量和地层泥质含量之间的关系：砂泥岩地层(3-1-13)碳酸盐岩地层(3-1-14)式中：和的单位为；泥质含量为小数，由自然伽马测井资料拟定。岩石粘聚力（内聚强度或抗剪切强度）根据Brules和Coates的研究成果，粘聚力的计算公式为：(3-1-15)式中：为岩石动态泊松比，无量纲；为岩石的体积密度，单位为；为岩石的纵波速度，；为地层的泥质含量，小数。岩石的抗张强度(3-1-16)岩石的内摩擦角对于岩石内摩擦角的拟定，斯伦贝谢和西方阿特拉斯公司在计
8、算时假定岩石的所有内摩擦角都为，这与实际状况有一定的误差，岩石类型和构成岩石的颗粒的相对大小对内摩擦角有影响。一般岩石的摩擦角在150450之间。根据Brie强度公式，摩擦角与泊松比之间关系有：(3-1-17)此外，内摩擦角也可以按石油大学提出的经验关系式进行计算：(3-1-18)其中：。具体采用哪个公式，视实际状况而定。4、岩石物理参数计算实例分析根据上述参数计算公式，运用的相应测井资料对安塞油田的几口井进行岩石机械特性进行分析。从安2051和泌356井的岩石物理参数成果图（图3-1-2、图3-1-3）可以看出，这两口井地层岩石的多种弹性模量参数(动态)比较大，杨氏模量在40000MPa-6
9、5000MPa之间，即岩石的抗破坏能力比较强。从图中还可以看出，泥岩层的抗剪切强度和单轴抗拉强度比砂岩层的大，这是在储层改造时为什么泥岩层能作为遮挡层的因素之一。固然，井壁岩石的破坏不仅与岩石自身的强度有关，还与地应力有关。一般状况下，若不考虑地应力的影响，泥岩抗破坏的能力比砂岩强。图3-1-2 安2051井岩石机械特性成果图图3-1-3 泌356井岩石机械特性成果图图3-1-4 安2051井实测横波、合成横波及力学参数对比图计算多种弹性模量参数需要同步具有有纵、横波时差，并且弹性模量是计算地应力及多种破裂压力的基本参数，因此，纵、横波时差是地应力研究中的重要原始数据。在缺少横波时差资料的状况
10、下，可由纵波时差以及密度来合成横波时差。下面给出了安2051井采用实测横波时差和合成横波时差计算弹性力学参数的对比图。从图3-1-4看出，合成的横波时差与实测横波时差几乎相等，并且分别由这两种时差计算得出的多种弹性参数曲线基本重叠，只有少数深度点相差较大，从总体上来讲，通过合成的横波来计算弹性力学参数的精度还是比较高的。因此，在缺少横波时差的状况下，可以由纵波时差与密度资料来合成横波时差，并且能获得比较精确的弹性力学参数，进而进行地应力以及井壁稳定性等方面的分析。第二节 地应力大小及方位1、地应力地应力是存在于地层中的未受工程扰动的天然应力，也成岩体初始应力、绝相应力或原岩应力。它是由于地壳内
11、部的垂直运动和水平运动及其她因素而引起介质内部单位面积上的作用力。地壳中的不同地区，不同深度地层中的地应力的大小和方向随空间和时间的变化而变化构成应力场。地层中每一种质点的地应力均有其大小和方向，其中涉及最大水平主应力、最小水平主应力、垂向应力的大小和方向。地应力是石油勘探与开发中的一种重要基本参数。含油气盆地构造的形成和演化是在一定的地应力场作用下的产物，只有弄清含油气盆地、含油气区块的地应力场分布，才干对的结识古构造行迹的发生演化历史，才干有效地分析和解决油气勘探开发的有关问题。目前已经结识到地应力对油气勘探开发的作用和影响越来越多地从各个方面体现出来，如：地质构造形成玉演化是构造应力作用
12、及变化的成果；储层中油气运移和汇集与地应力有关，油气总是由强应力区向弱应力区运移；天然裂缝和裂隙面与最大主应力方向平行；油田应力场状态决定着断层的形态和分布；在渗入率各向异性、低渗入率油田中，主渗入率方向与最大目前最大水平主应力方向一致；在钻井过程中，井壁的稳定性与地层岩石的力学性质、地层剖面的地应力状态有密切关系；油井采油过程中的出砂与地层的岩石力学性质、油层的应力环境、出砂指数有关；油层改造过程中，地应力场状态、地层岩石的力学性质决定着水力压裂的裂缝的形态、方位、高度和宽度，影响着压裂的增产效果等等。地应力的大小、方向、分布规律及其演化史是油气勘探开发中地应力研究的重要内容，而岩石的力学性
13、质、储层的孔隙压力、地层温度、构造应力、重力及地层剥蚀等是影响油田应力场状态的重要因素。2、地应力测量的途径及措施地应力测量从原理上可分为直接测量与间接测量两类。前者通过测量岩石的破裂，直接拟定应力，例如20世纪70年代发展起来的水压裂法就是一种直接测量的措施；后者通过测井岩石的变形和物性变化来拟定介质的受力状态，如根据岩石受力时的变形特性、弹性波速度变化、电阻率变化、声发射特性和矿物颗粒的显微构造变化拟定介质的受力状态，20世纪50年代发展起来的应力解除法就是一种间接测量措施。地应力测量从内容上可分为绝对值测量与相对值测量。前者是测量岩石所受的应力数值与方向，后者是测量固定点上随时间变化的应
14、力状态。地应力测量从仪器安装的形式上可分为钻孔法和非钻孔法。前者将测量仪牢固地贴附于钻孔底部，测量钻孔附近的应力状态；非钻孔法能理解较大空间岩体的应力状态，能避免由于钻孔开挖而来的应力变化的影响，但由于仪器读数换算为岩石应力时，影响因素较多。因此，较常用的地应力测量措施大多是钻孔法。3、运用测井资料计算地应力地层间或层内的不同岩性岩石的物理特性、力学特性和地层孔隙压力异常等方面的差别导致了层间或者层内地应力分布的非均质性。地应力大小是随处层性质变化的：山前构造带地应力重要来源于上覆地层压力机地质构造运动产生的构造力，在不同性质的地层由于其抵御外力的变形性质不同，因而其承受构造力也不相似。若依托
15、实测找寻层内或层间地应力的分布规律，这是不切实际的。因此，可以结合测井资料和分层地应力解释模型，可分析层内或层间地应力的大小。水力压裂法或声发射凯塞尔效应法只可以测试岩心点的地应力值，而一般我们需要理解层内或层间地应力的分布规律，但用这两种措施是不切实际的。测井资料具有持续、来源广、成本低的特点，因而结合分层地应力理论，建立分层地应力剖面测井解释技术，具有非常重要的意义。地应力计算模型中对于垂向应力的拟定，普遍采用了垂向应力为一主应力并且等于上覆岩层压力的假设，即：(3-2-1)如果用式(3-2-1)来计算垂向应力时，需要有岩石体积密度随深度变化的函数，但是，实际岩层的密度随深度的变化关系无法用一种简朴的函数来体现，必须按深度用分段求和的措施来计算上覆岩层压力的值。因此，可以用下式计算垂向应力。(3-2-2)式中：为目的层以上第段地层的平均密度，；为目的层以上第地层厚度，；为重力加速度，；为目的层以上的地层层数，整数；垂向应力，。人们在地应力成因和分布规律研究的基本上，提出了如下几种水平应力经验计算模型：双井径模型该模型在单轴应变的基本上，引入了一种应力不平衡系数来表征水平应力的非均质性：(3-2-3)
