流体识别方法和装置的制造方法
【技术领域】
[0001] 本发明涉及地质勘探技术领域，特别涉及一种流体识别方法和装置。
【背景技术】
[0002] 随着勘探开发程度的不断提高，勘探目标逐渐由最初的构造油气藏转变为岩性油 气藏和隐蔽性油气藏，勘探策略也由最开始的寻找构造圈闭转变为寻找岩性异常体或者直 接进行流体检测。近年来，基于地震资料反演储层弹性参数的方法已被广泛应用于储层预 测当中，并且储层预测的目的也不再局限于构造、岩性圈闭识别，而是要对储层含流体性 (包括流体类型以及流体饱和度)进行检测，这就对流体识别技术提出了更高的要求。
[0003] 上世纪70年代，基于地震资料振幅信息的流体识别技术作为首个直接进行地下含 油气性检测技术被广泛应用在实际的勘探开发中。随后，依据地震波在地下介质中传播时 会在储层油气水分界面处产生沿水平方向振幅强度稳定的反射波同相轴的平点勘探方式 也开始被使用。80年代提出的AV0技术，主要是利用地震波反射振幅随偏移距或者入射角的 变化规律判断含气砂岩储层，使得利用AV0分析技术检测流体成为可能。流体因子是利用纵 横波速度的加权组合来预测储层流体，随后将流体因子引入到反射系数域，计算纵横波反 射系数的加权组合。1999年之后，对流体因子的研究逐渐从反射系数域转变到波阻抗域，其 中λρ、μρ属性是阻抗域最主要且应用效果最显著的流体因子，λρ_μρ属性交会图要比常规的 纵横波阻抗交会图对储层流体有更好的识别能力。
[0004] 随后，Russell等对流体因子进行了进一步的改进得到了适用性更强的流体因子。 近年来，对于流体因子的研究不断深入，应用领域逐渐从叠后发展到叠前、从各向同性介质 发展到各向异性介质、从时间域发展到频率域。
[0005] 然而，如何准确有效地对流体进行识别，目前尚未提出有效的解决方案。
【发明内容】
[0006] 本发明实施例提供了一种流体识别方法，以达到准确识别流体的目的，该方法包 括：
[0007] 获取研究区的叠前角道集地震数据体；
[0008] 对所述叠前角道集地震数据进行分角度叠加，得到分角度叠加数据体；
[0009] 采用所述分角度叠加数据体进行叠前AV0反演，得到纵波阻抗体和密度体；
[0010] 对所述纵波阻抗体进行时频分解处理，得到分频波阻抗体；
[0011] 结合所述密度体，对所述分频波阻抗体进行反演运算，得到所述研究区的储层地 震频散属性；
[0012] 根据所述储层地震频散属性进行流体识别。
[0013] 在一个实施方式中，按照以下公式结合所述密度体，对所述分频波阻抗体进行反 演运算，得到所述研究区的储层地震频散属性：
[0014]
[0015] 其中，z (t，f)表示特定频率时的波阻抗体，z (t，fo)表示参考频率时的波阻抗体，P 表示所述密度体，D表示反演得到的速度频散值。
[0016] 在一个实施方式中，所述参考频率为所述叠前角道集地震数据体的主频率。
[0017] 在一个实施方式中，对所述叠前角道集地震数据进行分角度叠加，得到分角度叠 加数据体，包括：
[0018] 确定分角度叠加的多个角度区间；
[0019] 根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠加，得到所述多 个角度区间中各个角度区间对应的部分叠加数据体。
[0020] 在一个实施方式中，对所述纵波阻抗体进行时频分解处理，得到分频波阻抗体，包 括：
[0021 ]通过S变换对所述纵波阻抗体进行时频分解处理，得到分频波阻抗体。
[0022] 本发明实施例还提供了一种流体识别装置，以达到准确识别流体的目的，该装置 包括：
[0023] 获取模块，用于获取研究区的叠前角道集地震数据体；
[0024] 叠加模块，用于对所述叠前角道集地震数据进行分角度叠加，得到分角度叠加数 据体；
[0025] 反演模块，用于采用所述分角度叠加数据体进行叠前AV0反演，得到纵波阻抗体和 密度体；
[0026] 时频分解模块，用于对所述纵波阻抗体进行时频分解处理，得到分频波阻抗体；
[0027] 运算模块，用于结合所述密度体，对所述分频波阻抗体进行反演运算，得到所述研 究区的储层地震频散属性；
[0028] 识别模块，用于根据所述储层地震频散属性进行流体识别。
[0029] 在一个实施方式中，所述运算模块具体用于按照以下公式结合所述密度体，对所 述分频波阻抗体进行反演运算，得到所述研究区的储层地震频散属性：
[0030]
[0031 ]其中，z (t，f)表示特定频率时的波阻抗体，z (t，f 〇)表示参考频率时的波阻抗体，P 表示所述密度体，D表示反演得到的速度频散值。
[0032] 在一个实施方式中，所述参考频率为所述叠前角道集地震数据体的主频率。
[0033] 在一个实施方式中，所述叠加模块包括：
[0034]角度区间确定单元，用于确定分角度叠加的多个角度区间；
[0035] 叠加单元，用于根据所述多个角度区间对所述叠前角道集地震数据进行分角度叠 加，得到所述多个角度区间中各个角度区间对应的部分叠加数据体。
[0036] 在一个实施方式中，所述时频分解模块具体用于通过S变换对所述纵波阻抗体进 行时频分解处理，得到分频波阻抗体。
[0037] 在本发明实施例中，基于研究区叠前AV0反演得到的纵波阻抗数据体与密度体，利 用时频分析技术对波阻抗数据体进行时频分解得到不同瞬时频率的阻抗体，然后结合地震 波阻抗依赖分频反演公式提取地震频散数据体，最后根据频散属性实现流体识别，从而实 现对储层流体分布规律的有效预测，提高了流体识别的准确性。
【附图说明】
[0038] 此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不 构成对本发明的限定。在附图中：
[0039] 图1是根据本发明实施例的流体识别方法流程图；
[0040] 图2是根据本发明实施例的储层饱和不同流体时纵波速度值随频率变化关系图；
[0041] 图3是根据本发明实施例的储层饱和不同流体时纵波频散值随频率变化关系图；
[0042] 图4是根据本发明实施例的利用地震波阻抗依赖分频反演得到的地震频散属性剖 面示意图；
[0043] 图5是根据本发明实施例的利用地震波阻抗依赖分频反演得到的地震频散属性沿 目的层的切片示意图；
[0044]图6是根据本发明实施例的流体识别装置的结构框图。
【具体实施方式】
[0045] 为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对 本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并 不作为对本发明的限定。
[0046] 鉴于地震波在地下含流体介质中传播时会在不同频带范围和不同尺度上发生反 射波能量衰减与速度频散，而这种现象主要是由储层孔隙流体的存在导致的。频散属性作 为常规储层弹性信息的一个延伸和补充，从另一个角度展示了地震数据中被主频分量所掩 盖的高频和低频信息，对于岩性及隐蔽油气藏的识别尤其是对流体检测具有独特的优势。 与频率有关的速度变化主要受到储层孔隙流体的影响，并且流体的流动性（即岩石渗透率 与流体粘滞性的比值)是影响地震波速度频变特征的主要因素。
[0047] 发明人考虑到可以基于波阻抗依赖频散属性反演进行流体识别，具体的，可以综 合利用波