ept软件模块培训教材_弹性波阻抗反演

EPT软件/功能模块

培训系列教材

GMAX v1.0 – Elastic Impedance Inversion

弹性波阻抗反演

1.模块功能

2.原理和方法

3.参数和使用说明

4.应用注意事项

EPT公司

,

1.模块功能

弹性波阻抗反演（EI）是叠前地震反演重要方法之一。基于流体置换模型技术，应用纵波声波时差、密度、泥质含量、孔隙度、含水饱和度和骨架、流体的各种弹性参量，反演井中横波速度。根据井中纵波速度、横波速度和密度计算井中弹性波阻抗，在复杂构造框架和多种储层沉积模式的约束下，采用地震分形插值技术建立可保留复杂构造和地层沉积学特征的弹性波阻抗模型，使反演结果符合研究区的构造、沉积和异常体特征。采用广义线性反演技术反演各个角度的地震子波，得到与入射角有关的地震子波。在每一个角道集上，采用宽带约束反演方法反演弹性波阻抗，得到与入射角有关的弹性波阻抗。最后对不同角度的弹性波阻抗反演纵横波阻抗，进而获得泊松比等弹性参数, 对储层的几何、物性和含流体特性进行精细描述。叠前地震弹性参数反演的关键技术包括:

◆基于流体替换模型的井中横波速度反演技术

◆与偏移距有关的子波反演

◆复杂地质构造情况下弹性阻抗建模

◆纵横波阻抗、泊松比、拉梅系数和剪切模量反演

2.原理和方法

地震反射振幅不仅与分界面两侧介质的地震弹性参数有关，而且随入射角变化而变化。叠前弹性波阻抗反演技术利用不同炮检距的地震数据及横波、纵波、密度等测井资料，联合反演出与岩性、含油气性相关的多种弹性参数，综合判别储层物性及含油气性。正是由于叠前弹性波阻抗反演利用了大量地震及测井信息，所以进行多参数分析的结果较叠后声阻抗反演在可信度方面有很大提高，可对含油气性进行半定量—定量描述。

传统的A VO 和岩石物理分析是提取和分析纵横波速度的异常变化来确定孔隙流体和岩性的变化。纵横波速度和密度对反射系数的重要性，可以从平面波的Zoeppritz方程中看出。但是，在波动方程中，Md2U/dX2= d2U/dX2，（U是位移），其表达式并不与地震波速度直接相关，而与岩石密度和弹性模量相关。因此，直接考虑泊松比、拉梅系数和岩石剪切模量比采用地震波速度能更好地反映岩石物理特征。地震的纵波速度与含孔隙流体岩石特征的关系是靠体变模量K

联系在一起的，体变模量K 和纵波速度都包含了最敏感的流体检测因子拉梅系数(λ)，但都因纵波速度和体变模量中包含μ而减弱了λ的敏感性，对这可以由关系式，Vp 2=(λ+2μ)/ρ和Vs 2=μ/ρ看出。最近，A VO 的反演试图包含密度参数，以获取准确的弹性模量参数。对于反演的准确性而言，其精度随未知量的增多而降低，使方程的解变得不稳健，提取的参数也就更不准确。因此，在反演中将考虑用弹性模量/密度关系或阻抗参数，具体为：

AI 2=(Vp*ρ)2=(λ+2μ)*ρ

SI 2=(Vs*ρ)2=μ*ρ

其中，AI 为纵波波阻抗，SI 为横波波阻抗。通过叠前地震资料反演得到的纵横波波阻抗通过如下变换，我们可以得到拉梅系数和岩石密度的乘积剖面和剪切模量和岩石密度的乘积剖面，即λρ和μρ。

λ= Vp 2ρ-2 Vs 2ρ

μ= Vs 2ρ

λρ= PI 2-2 SI 2 μρ= SI 2

2222)(1)(5.0)(1)(5.0)(2AI SI AI SI Vp Vs Vp Vs --=--=+=μλλσ

上式表明泊松比σ相对来说对流体检测因子拉梅系数λ的较敏感，是一个较好的流体检测弹性参数。垂直入射（自激自收）时，反射系数为

1

1221122Vp Vp Vp Vp R pp ρρρρ+-=

其中：R pp 为纵波反射系数；ρ1，ρ2对应为上、下介质密度；Vp 1，Vp 2分别为

上、下层介质的纵波速度。而非垂直入射（炮检距不为零）时，纵、横波的反射和透射系数是以佐布里兹（Zoeppritz ）方程的矩阵形式

⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡11112112221

122111122111222221112221111221122112cos 2sin cos sin 2sin 2cos 2sin 2cos 2cos 2cos 2cos 2sin sin cos sin cos cos sin cos sin ϕθθθϕρρϕρρϕϕϕρρθρρϕθϕθϕθϕθϕθp s p p p s s p s p s p s s p ps pp ps pp V V V V V V V V V V V V V V V T T R R

表示的。其中：R pp 和R ps 分别为纵、横波反射系数，T pp 和T ps 分别为纵、横波透射系

数。但该式并未直观表述纵、横波速度及密度对反射系数的贡献。Connolly 等学者对上述反射系数表达式作出近似。Connolly 定义P 波入射角θ的弹性波阻抗EI(θ)为

)sin 41()sin 8()tan 1(222θθθρK K S P V V EI --+=

弹性波阻抗的基本作用是代替与入射角相关的P 波反射率，就象AI 代表零偏移距反射率一样。当 0=θ时，纵波反射系数为

1

212)0(AI AI AI AI R pp +-= 此时，弹性阻抗与声阻抗相等，即Vp AI EI ρ==。

如果我们定义反射界面上下介质的弹性波阻抗EI 1和EI 2的数学表达式为：

)sin 41(1)sin 8(1

)tan 1(11222θθθρK K S P V V EI --+= )

sin 41(2)sin 8(2)tan 1(22222θθθρK K S P V V EI --+=

式中下标1、2分别表示界面上、下介质，K 的表达方式为：

2

)()(22221

1⎥⎦⎤⎢⎣⎡+=P S P S V V V V K 根据上述公式定义的弹性波阻抗，可得入射角为θ时的反射系数可近似为：

1

212)(EI EI EI EI R pp +-≈θ 由上式可见，非垂直入射时反射系数表达式与垂直入射时反射系数表达式一样，这样我们就可以借用传统相对成熟的叠后波阻抗反演方法反演弹性波阻抗，这也是Connolly 定义弹性波阻抗的原因。

由弹性波阻抗的表达式可得：

()()[](){}()()[]θρθρρρθ422sin ln 2ln 4sin ln ln 2ln 4)ln()ln(cos S S P V K V K V EI ++++-= 当入射角小于30度时，tan 2θ≈ sin 2θ，sin 4θ≈0, 上式可简化为

()()[]{}θρρ2sin ln 2ln 4)ln()ln()ln(S P P V K V V EI +-+≈