技术领域
[0001] 本
发明涉及地球物理勘探的储层预测领域，更具体地讲，涉及一种扩展弹性阻抗反演含流体性概率的预测方法。
背景技术
[0002] 叠前弹性阻抗(EI)反演是利用不同炮检距(或入射
角)道集数据及横波、纵波、
密度等
测井资料，联合反演出与岩性、含油气性相关的多种弹性参数，综合判断储层物性及含油气性。常规的弹性阻抗方法存在以下问题：求取的值随角度变化在尺度上发生很大变化，求取的反射系数也不太稳定。由于
地震振幅随偏移距的变化(AVO，Amplitude versus Offset)效应的影响，反演的EI值随入射角变化而发生剧烈变化。近、中、远不同角度的弹性波阻抗值的量级不在一个级别上，虽然弹性阻抗值有变化，但相对于量级的变化是极其微小的，这给不同入射角下的弹性阻抗的交会对比带来困难。
[0003] 扩展弹性阻抗(EEI)的提出，消除了入射角对EI的影响，提高了入射角的适用范围，提高了对气
水差异反映的灵敏度。但目前对反演结果的分析方法主要通过交会图的方式，存在人为性大、分析结果
稳定性差等问题。
[0004] 图1是示出
现有技术的扩展弹性阻抗反演含流体性概率的预测方法的
流程图。图1中示出的方法是现有技术中的对反演结果的分析方法，主要通过交会图的方式，下面将参照图1描述现有技术的扩展弹性阻抗反演含流体性概率的预测方法。
[0005] 如图1所示，在步骤S101，通过测井资料计算气层、水层不同角度EEI值。
[0006] 在步骤S102，对计算的气层、水层的不同角度的EEI值制作交会图。由于气层、水层随入射角的变化，其EEI值的变化规律有差异，反映在交会图上的分布
位置不同，分析气层(或水层)主要分布区域，在交会图上圈出有利区，建立解释模板。
[0007] 在步骤S103，对地震资料进行反演。通过反演得到不同角度的扩展弹性阻抗反演数据体，并去除非储层段，仅保留储层段的反演值。
[0008] 在步骤S104，利用测井资料制作交会图解释模板，对地震资料的反演结果进行解释，获得含气性预测结果。
[0009] 然而，通过交会图建立解释模板的方式，气层(或水层)分布区域的划分受人为因素干扰较大，难以建立统一的标准，解释结果的稳定性差，对解释结果的准确性影响较大。
[0010] 因此，需要一种不受人为因素干扰的扩展弹性阻抗反演含流体性概率的预测方法。
发明内容
[0011] 针对现有技术中存在的气层或水层分布区域的划分受人为因素干扰较大的问题，本发明提出了一种扩展弹性阻抗反演含流体性概率的预测方法。所述扩展弹性阻抗反演含流体性概率的预测方法可包括以下步骤：通过测井资料计算各个样本的不同角度的扩展弹性阻抗EEI值，来获得EEI(θ)j，其中，θ表示样本的角度，j表示样本号；计算出每个样本的不同角度的EEI差值ΔEEIj，其中，j表示样本号；应用直方图分析ΔEEIj的分布范围，获得气层和水层的概率分布范围，以建立气层概率分布的量化分析模板，或者建立气层和水层概率分布的量化分析模板；对地震资料进行反演，得到不同角度的扩展弹性阻抗反演剖面，并去除非储层段，仅保留储层段的反演数据；利用储层段的反演数据计算扩展弹性阻抗差值数据体；将建立的量化分析模板应用于扩展弹性阻抗差值数据体，获得含气概率数据体，通过对目的层进行时窗分析来获得含气概率平面分析结果。
[0012] 获得EEI(θ)j的步骤可包括：在目标区域内选取典型井，在选取的典型井内针对目的层段的气层、水层确定分析样本，针对确定的分析样本计算不同角度的扩展弹性阻抗EEI值，来获得EEI(θ)j。
[0013] 根据本发明示例性
实施例的扩展弹性阻抗反演含流体性概率的预测方法主要通过EEI差值计算和直方图建立量化分析模板，将量化分析模板应用于实际数据体(扩展弹性阻抗差值数据体)，最终获得流体(含气或含水)的概率分布图。应用本定量分析扩展弹性阻抗反演结果的方法后，可以降低采用传统交会图圈定有利区域分析方法的人为性，提高分析结果的
精度和效率。
附图说明
[0014] 通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会变得更加清楚，其中：
[0015] 图1是示出现有技术的扩展弹性阻抗反演含流体性概率的预测方法的流程图；
[0016] 图2是示出根据本发明示例性实施例的扩展弹性阻抗反演含流体性概率的预测方法的流程图；
[0017] 图3是由测井资料计算的气层(用方形表示)、水层(用五角星表示)的入射角为5°和25°的EEI值交会图；
[0018] 图4是统计气层和水层的远、近入射角扩展弹性阻抗差分布的直方图；
[0019] 图5是特定井去除了非储层的须二段扩展弹性阻抗反演差值的剖面图；
[0020] 图6是利用量化分析模板获得的含流体性概率的剖面图；
[0021] 图7为通过概率统计获得的目的层段的含气概率的平面图。
具体实施方式
[0022] 现在，将参照附图更充分地描述根据本发明的示例实施例。
[0023] 图2是示出根据本发明示例性实施例的扩展弹性阻抗反演含流体性概率的预测方法的流程图。
[0024] 参照图2，在步骤S201，通过测井资料计算各个样本的不同角度的扩展弹性阻抗(EEI)值，来获得EEI(θ)j，其中，θ表示样本的角度，j表示样本号。具体地讲，可在目标区域内选取典型井，在选取的典型井内针对目的层段的气层、水层确定分析样本，针对确定的分析样本计算不同角度的扩展弹性阻抗EEI值，来获得EEI(θ)j。
[0025] 由于气层、水层随入射角的变化，其EEI值的变化规律有差异。因此，在步骤S202，可计算出每个样本的不同角度的EEI差值ΔEEIj，j表示样本号；
[0026] 在步骤S203，建立气层(或水层)概率分布的量化分析模板。具体地讲，应用直方图分析ΔEEIj的分布范围，获得气层和水层的概率分布范围，以建立气层概率分布的量化分析模板，或者建立气层和水层概率分布的量化分析模板。
[0027] 在步骤S204，对地震资料进行反演，得到不同角度的扩展弹性阻抗反演剖面，并去除非储层段，仅保留储层段的反演数据。
[0028] 在步骤S205，利用储层段的反演数据计算扩展弹性阻抗差值数据体。
[0029] 在步骤S206，将在步骤S203建立的量化分析模板应用于扩展弹性阻抗差值数据体，获得含气概率数据体，通过对目的层进行时窗分析来获得含气概率平面分析结果。
[0030] 下面将通过图3至图7来描述通过图2中示出的根据本发明示例性实施例的扩展弹性阻抗反演含流体性概率的预测方法的各种示图。
[0031] 图3是由测井资料计算的气层(用方形表示)、水层(用五角星表示)的入射角为5°和25°的EEI值交会图。参照图3，在小入射角度弹性阻抗相同的情况下，水层在25°的EEI值明显大于气层，由此可以看出，通过分析远、近入射角扩展弹性阻抗差异，可以达到区分水层和气层的目的。
[0032] 图4是统计气层和水层的远、近入射角扩展弹性阻抗差分布的直方图，其中，ΔEEIj＝EEI(25°)j-EEI(5°)j，横轴为ΔEEI值，图4的上部的纵轴为样点出现累积频次，即分布概率，图4的下部的纵轴为统计的样点出现的频次。图4中的气层的ΔEEIj主要分布在650以下，80％以上分布在400以下，而80％以上的水层分布在450～1100，450～650是气层和水层均出现的区间，但出现的概率有所不同。根据图4中示出的统计结果可以获得气层、水层的分布概率关系。
[0033] 图5是特定井(例如，蓬莱2井)去除了非储层的须二段扩展弹性阻抗反演差值的剖面图，图5中的井附近须二上段弹性阻抗差值相对小，大多在700以下，而须二下段弹性阻抗差值大，以800以上的高值为主，该井实际测井解释储层主要发育在须二上段和下段，其中上段主要为气层，下段为水层。
[0034] 图6是利用量化分析模板获得的含流体性概率的剖面图，图6的上部示图表示含气概率的剖面图，图6的下部示图表示含水概率的剖面图。从图6可以看出，须二上段含气概率较高，大部分在60％以上，须二下段含气概率较低，大部分小于20％，图6的上部示图中示出的概率与下部示图中示出的概率呈互补的趋势。
[0035] 图7是通过概率统计获得的目的层段的含气概率的平面图，
颜色最浅的区域表示含气概率高，颜色最深的区域表示含气概率较低，颜色居中的区域表示含气概率低，从颜色居中的区域到颜色最深的区域再到颜色最浅的区域，含气概率逐渐增高到1，由此可定量地识别出含气有利区。
[0036] 从以上描述可以看出，根据本发明示例性实施例的扩展弹性阻抗反演含流体性概率的预测方法主要通过EEI差值计算和直方图建立量化分析模板，将量化分析模板应用于实际数据体(扩展弹性阻抗差值数据体)，最终获得含流体(含气或含水)的概率分布图。应用本定量分析扩展弹性阻抗反演结果的方法后，可以降低采用传统交会图圈定有利区域分析方法的人为性，提高分析结果的精度和效率。另外，根据本发明示例性实施例的扩展弹性阻抗反演含流体性概率的预测方法，预测了陆相碎屑岩气藏的含气概率分布，预测结果表明高产气井和气水同产的工业气井大多位于高、中含气概率区，其含气性预测符合率达到81％，具有较高的预测精度，对陆相碎屑岩气藏具有良好的应用前景。
[0037] 尽管已参照本发明的特定示例性实施例显示和描述了本发明，但本领域的技术人员应该理解，在不脱离
权利要求及其等同物限定的本发明的精神和范围的情况下，可在形式和细节上进行各种改变。
