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 22222)(1)(5.0)(1)(5.0)(2AI SI AI SI Vp Vs Vp Vs --=--=+=μλλσ上式表明泊松比σ相对来说对流体检测因子拉梅系数λ的较敏感，是一个较好的流体检测弹性参数。

垂直入射（自激自收）时，反射系数为11221122Vp Vp Vp Vp R pp ρρρρ+-=其中：R pp 为纵波反射系数；ρ1，ρ2对应为上、下介质密度；Vp 1，Vp 2分别为上、下层介质的纵波速度。

而非垂直入射（炮检距不为零）时，纵、横波的反射和透射系数是以佐布里兹（Zoeppritz ）方程的矩阵形式⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡11112112221122111122111222221112221111221122112cos 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=θ时，纵波反射系数为1212)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)()(222211⎥⎦⎤⎢⎣⎡+=P S P S V V V V K 根据上述公式定义的弹性波阻抗，可得入射角为θ时的反射系数可近似为：1212)(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 +-+≈图1 弹性波阻抗反演技术路线式中有Vp、Vs、 三个未知数，利用不同入射角数据进行反演，就得到多个入射角弹性阻抗，由此建立方程组可求取其它弹性参数，用于岩性及油气预测。

弹性波阻抗反演的基本思路如图1所示。

基于流体置换模型技术，应用纵波声波时差、密度、泥质含量、孔隙度、含水饱和度和骨架、流体的各种弹性参量，反演井中横波速度（图2）。

EPT的流体置换模型考虑了流体的饱和度、孔隙形状和流体压力的影响。

流体压力改变会引起孔隙尤其是微细孔隙的闭合与开启。

因此，考虑形状不同的孔隙在受载应力下的变化和响应非常重要。

孔隙的形状及取向分布是影响介质弹性刚度系数的主要参量，在不同的流体压力及应力场作用下，孔隙的形状会发生变形，从而改变了孔隙的形状与大小。

实际上，扁率小于0.01或更小的细微孔隙，其密度在不足1%的情况下就会发生速度20%的变化。

图3 某油田井实测横波(兰色) 和反演横波(红色)对比图图3所示为车660井实测横波(兰色)和反演横波(红色)对比图，有图可见井中实测横波和反演横波之间吻合很好，证明了基于流体置换模型技术反演横波速度的方法技术是切实的。

根据井中纵波速度、横波速度和密度计算井中弹性波阻抗，在复杂构造框架和多种储层沉积模式的约束下，采用地震分形插值技术建立可保留复杂构造和地层沉积学特征的弹性波阻抗模型，使反演结果符合研究区的构造、沉积和异常体特征。

采用广义线性反演技术反演各个角度的地震子波，得到与入射角有关的地震子波。

在每一个角道集上，采用宽带约束反演方法反演弹性波阻抗，得到与入射角有关的弹性波阻抗。

最后对不同角度的弹性波阻抗反演纵横波阻抗，进而获得泊松比等弹性参数。

图4为两口井的连井线泊松比剖面。

图2 流体置换模型3参数和使用说明3.1 软件流程弹性波阻抗反演的软件使用基本流程如下：A. 数据准备B. 层位标定C. 井中横波反演D. 井中弹性波阻抗反演E. 子波提取E. 地震弹性波阻抗反演F. 地震弹性参数反演G. 地质解释3.2 层位标定和子波提取使用提示：层位标定只对入射角道集中的一个入射角数据进行，一般对最小入射角的那一道。

若信号质量不佳，可选其它入射角的。

这样保证了所有入射角数据有同样的时深关系。

3.3 井中横波反演( Well VS Inversion )参数描述和提示：1)在进行横波速度反演时需要使用下面的五种曲线，若曲线不全，反演无法进行。

b)，c)，d)三种曲线可由测井解释（GMbase: Reservoir Interpretation）获得，或者拷贝油田已解释的数据。

a)P-wave sonic log纵波声波测井曲线b)Shale content (%)解释的泥质含量曲线c)Porosity解释的孔隙度曲线d)Water Saturation解释的含水饱和度曲线e)Density岩石密度曲线2)Number of Segments 由于井中深层和浅层的处理参数不同，不同的沉积环境也需要不同的处理参数，用户可以根据岩性、地质分层等数据，将测井曲线分成几段来处理。

指定处理的段数。

最多可以分16段。

3)Depth From … To …当前段的深度范围4)P-wave Velocity of Shale当前段泥岩的纵波速度5)S-wave Velocity of Shale当前段泥岩的横波速度6)Shale Density 当前段泥岩的密度7)P-wave Velocity of Sand当前段砂岩的纵波速度8)S-wave Velocity of Sand当前段砂岩的横波速度9)Sand Density当前段砂岩的密度10)Matrix Bulk Modulus当前段基质的体积模量11)Matrix Density当前段基质的密度12)Fluid Density当前段流体的密度13)Fluid Bulk Modulus当前段流体的体积模量14)Statistic Parameters 设置深度间隔后，按“Statistic Parameter”按钮，系统根据该井的测井曲线自动统计各个深度段内的速度和密度。