文档介绍：岩石物理参数计算及应力研究-llzlllo
第三章岩石物理参数计算及应力研究
第一节岩石物理参数计算
地层岩石是地应力的载体,岩石物理性质对地应力的传递、衰减、集中、分散都会产生很大的影响,岩石物理参数与岩体赋存的地应力密切相关,岩石物理参数计算是地应力研究的必然步骤。
通过纵、横波时差和密度等测井资料,可以计算地层条件下的岩石动态弹性模量,在此基础上,可以进行地应力分析、井眼稳定性分析、地层出砂分析、以及人工压裂设计等方面的研究。
岩石物理参数包括岩石弹性参数和岩石机械强度参数。岩石弹性参数主要有泊松比μ、杨氏模量E、剪切模量G、体积模量K、体压缩系数Cb和Cma、有σ效应力系数系数α(比奥特系数);岩石机械强度主要有单轴抗压强度c、岩石
s的抗剪强度C0和岩石抗张强度t,以及内摩擦角?等。
1、岩石弹性参数
对于各向同性均匀介质岩石来说,利用牛顿第二定律和三维虎克定律,经数学推导,可导出计算声波速度在岩石介质中的波动方程:
Vp=+2G=ρ
G=E(1-)1 =ρ(1+μ)(1-2μ)?tP(3-1-1) Vs=ρE1 =ρ2(1+μ)?ts(3-1-2)
根据上述的波动方程,可以得出各种弹性参数与声波时差的关系式。①泊松比
定义为横向应变与纵向应变之比。
μ=?ts2-?t2
p
?ts2-?t2
p (3-1-3)
②切变模量
定义为施加的应力与切应变之比。
G=ρb
?t2
s?a (3-1-4)
③杨氏模量
定义为施加的轴向应力与法向应变之比。
E=2G(1+μ)
④体积模量
定义为静水压力与体积应变之比。
?14? ??a Kb=ρb- ?t23?t2?s??p(3-1-5) (3-1-6)
⑤体积压缩系数
定义为体积模量的倒数。即:
Cb=1 Kb(3-1-7)
⑥有效应力系数(Boit)
表示孔隙压力对岩石变形的影响,即:
α=1-Cma=Kb/Kma Cb(3-1-8)
式中:ρb为岩石体积密度,gcm3;?ts、?tp为纵、横波时差,μsft。公式中的a为单位转换系数。具体选择方法如下:如果密度单位为g/cm3,时差单位为μs/ft,弹性参数单位为psi,则a=?1010;如果密度单位为g/cm3,时差单位为μs/ft,弹性参数单位为MPa,则a=?107;如果密度单位为g/cm3,时差单位为μs/m,弹性参数单位为MPa,则a=109。
因此,利用阵列声波测井提供的纵、横波时差以及常规测井提供的密度资料就可以进行岩石弹性参数计算。
但是由于费用等原因,并不是每口井都开展声波全波列或阵列声波测井,因而不能直接获取横波时差资料,在研究中则可以通过构造内某些井已有的横波时差曲线资料来建立横波时差曲线计算式。研究表明,横波时差与纵波时差、地层密度和纵波波阻抗之间有很好的相关性。通过对安棚地区4口井的纵横波时差曲线进行分析后,建立了纵横波时差经验关系式:
?ts=?tp-+
R= (3-1-9)
图3-1是由上式纵横波时差关系式得出的横波时差与实测横波时差的关系图,从图中可以看出,大部分点分布在斜率约为450的直线上,计算的横波时差
与实测横波时差近似相等。
图3-1-1 合成横波时差与实测横波时差关系图
当研究区内没有一口井具有横波时差资料时,则可用下面的公式来合成横波时差曲线:
?ts=?tmas+(?tfs-?tmas)(?t?tp
fp-?tmap)-?tmas (3-1-10)
式中:?tmas、?tmap为岩石骨架的横波时差和纵波时差,μsft;?tfs、?tfp为流体的横波时差和纵波时差,μsft。
2、岩石动、静态弹性参数之间的转换方法
岩石弹性参数的常用测定方法有动态法和静态法两种。静态法是通过对岩样进行静态加载测其变形得到,所得弹性参数称之为静态参数;动态法则是通过测定超声波在岩样中的传播速度转换得到,所得弹性参数称之为动态参数。因此,用测井资料计算得到的弹性参数是动态参数。
根据地下岩层的应力形成、赋存和起作用的机理,特别是在应力幅值、加载速度和所引起的岩石变形等方面,更接近岩石静态测试的条件,另外,现有的力学本构关系一般是基于静态参数建立的,因此,在地应力计算和实际工程中应采用岩石的静态弹性参数。大量研究资料表明岩石的动态、静态弹性参数具有很好的相关性,且大部分情况下岩样的静态参数弹性模量小于其动态值。岩石动、静态弹性参数间存在较大差别,其原因主要是岩石中微裂缝和孔隙的存在。岩石这
种孔隙的弹性材料有别于各向同性、均质的线弹性体。微裂缝