岩石物理全参数计算及应力研究-llzlllo

第三章 岩石物理参数计算及应力研究

第一节 岩石物理参数计算

地层岩石是地应力的载体，岩石物理性质对地应力的传递、衰减、集中、分散都会产生很大的影响，岩石物理参数与岩体赋存的地应力密切相关，岩石物理参数计算是地应力研究的必然步骤。

通过纵、横波时差和密度等测井资料，可以计算地层条件下的岩石动态弹性模量，在此基础上，可以进行地应力分析、井眼稳定性分析、地层出砂分析、以及人工压裂设计等方面的研究。

岩石物理参数包括岩石弹性参数和岩石机械强度参数。岩石弹性参数主要有

泊松比μ、杨氏模量E 、剪切模量G 、体积模量K 、体压缩系数b C 和ma C 、有

效应力系数系数α（比奥特系数）；岩石机械强度主要有单轴抗压强度

c σ、岩石的抗剪强度0C 和岩石抗张强度t s ，以及内摩擦角ϕ等。

1、岩石弹性参数

对于各向同性均匀介质岩石来说，利用牛顿第二定律和三维虎克定律，经数学推导，可导出计算声波速度在岩石介质中的波动方程：

P ∆=-+-=+=t E G V 1)21)(1()1(2p μμρμρ

λ (3-1-1) s s t E G V ∆=+==1)1(2μρρ (3-1-2)

根据上述的波动方程，可以得出各种弹性参数与声波时差的关系式。

①泊松比

定义为横向应变与纵向应变之比。

22

2

25.0p s p

s t t t t ∆-∆∆-∆=μ (3-1-3)

②切变模量

定义为施加的应力与切应变之比。

a t G s b

⨯∆=2ρ (3-1-4)

③杨氏模量

定义为施加的轴向应力与法向应变之比。

)1(2μ+=G E

(3-1-5) ④体积模量

定义为静水压力与体积应变之比。

a t t K s p

b b ⨯⎪⎪⎭⎫ ⎝⎛∆-∆=22341ρ (3-1-6)

⑤体积压缩系数

定义为体积模量的倒数。即：

b b K C 1= (3-1-7)

⑥有效应力系数(Boit)

表示孔隙压力对岩石变形的影响，即：

b ma C C -=1α=K b /K ma (3-1-8)

式中：b ρ为岩石体积密度，3cm g ；s t ∆、p t ∆为纵、横波时差，ft s μ。公式中的a 为单位转换系数。具体选择方法如下：如果密度单位为3/cm g ，时差单位为ft s /μ，弹性参数单位为psi ，则101034.1⨯=a ；如果密度单位为3/cm g ，时差单位为ft s /μ，弹性参数单位为MPa ，则71029.9⨯=a ；如果密度单位为3/cm g ，时差单位为m s /μ，弹性参数单位为MPa ，则910=a 。

因此，利用阵列声波测井提供的纵、横波时差以及常规测井提供的密度资料就可以进行岩石弹性参数计算。

但是由于费用等原因,并不是每口井都开展声波全波列或阵列声波测井,因而不能直接获取横波时差资料，在研究中则可以通过构造内某些井已有的横波时差曲线资料来建立横波时差曲线计算式。研究表明，横波时差与纵波时差、地层密度和纵波波阻抗之间有很好的相关性。通过对安棚地区4口井的纵横波时差曲线进行分析后，建立了纵横波时差经验关系式：

34.1358.579.1+-∆=∆b p s t t ρ (3-1-9) 904.0=R

图3-1是由上式纵横波时差关系式得出的横波时差与实测横波时差的关系图，从图中可以看出，大部分点分布在斜率约为450的直线上，计算的横波时差

与实测横波时差近似相等。

图3-1-1 合成横波时差与实测横波时差关系图

当研究区内没有一口井具有横波时差资料时，则可用下面的公式来合成横波时差曲线：

()()()mas fp map p

mas fs mas s t t t t t t t t ∆-∆∆-∆∆-∆+∆=∆ (3-1-10)

式中：mas t ∆、map t ∆为岩石骨架的横波时差和纵波时差，ft s μ；fs t ∆、fp t ∆为流体的横波时差和纵波时差，ft s μ。

2、岩石动、静态弹性参数之间的转换方法

岩石弹性参数的常用测定方法有动态法和静态法两种。静态法是通过对岩样进行静态加载测其变形得到，所得弹性参数称之为静态参数；动态法则是通过测定超声波在岩样中的传播速度转换得到，所得弹性参数称之为动态参数。因此，用测井资料计算得到的弹性参数是动态参数。

根据地下岩层的应力形成、赋存和起作用的机理，特别是在应力幅值、加载速度和所引起的岩石变形等方面，更接近岩石静态测试的条件，另外，现有的力学本构关系一般是基于静态参数建立的，因此，在地应力计算和实际工程中应采用岩石的静态弹性参数。大量研究资料表明岩石的动态、静态弹性参数具有很好的相关性，且大部分情况下岩样的静态参数弹性模量小于其动态值。岩石动、静态弹性参数间存在较大差别，其原因主要是岩石中微裂缝和孔隙的存在。岩石这

种孔隙的弹性材料有别于各向同性、均质的线弹性体。微裂缝的存在对岩石静态变形的影响较大，而超声波可以绕过一些微裂缝传播。

在实际应用时，可通过岩石力学动、静态同步测试建立动、静态参数间的关系，从而把测井得到的动态参数转换为静态参数。

由于研究区及其邻近区块没有条件做岩石力学试验，本次研究引用了辽河油田和大庆油田的实验结果：

d s E E 7095.02526.0+=

(3-1-11) 37.036.0d s μμ⨯= (3-1-12)

式中：s E 、s μ为静态杨氏模量和静态泊松比；d E 、d μ为动态杨氏模量和动态泊松比，即测井资料计算结果。

式(3-11)和式(3-12)的相关系数分别为0.75和0.86。

3、岩石机械强度参数

目前，岩石机械强度参数还没有理论计算式，一般通过岩石力学测试来确定。为了克服岩石力学试验存在的测试费用昂贵和数据量少等缺点，研究人员通过岩石力学试验建立了岩石强度参数的经验计算式：

⑴单轴抗压强度C σ

Deer 和Miller （1996）根据大量的室内试验结果建立了砂泥岩的单轴抗压强度与岩石动态杨氏模量d E 和地层泥质含量sh V 之间的关系：

①砂泥岩地层

()sh d sh d c V E V E ⋅+-=008.010045.0σ

(3-1-13)

②碳酸盐岩地层 ()sh d sh d c V E V E ⋅+-=008.010026.0σ (3-1-14) 式中：c σ和d E 的单位为MPa ；泥质含量为小数，由自然伽马测井资料确定。

⑵岩石粘聚力0C （内聚强度或抗剪切强度）

根据Brules 和Coates 的研究结果，粘聚力0C 的计算公式为：

()()sh d d d b p V v C 78.0111211044.524150+⨯⎪⎪⎭⎫ ⎝⎛-+⨯-⨯=-μμμρ (3-1-15)