CN106338767A - 多参数地震属性分析方法及系统 - Google Patents

Abstract

公开了一种多参数地震属性分析方法及系统。该方法可以包括以下步骤：基于敏感地震属性的响应特征确定地质异常体；基于所述地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数；以及基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征。

Description

多参数地震属性分析方法及系统

技术领域

本发明涉及油气地球物理勘探领域，更具体地，涉及一种多参数地震属性分析方法及系统。

背景技术

在油气地球物理勘探领域，岩性隐蔽油气藏已成为现阶段油气勘探的重点目标之一，为此，利用地震资料开展薄储层的预测已成为越来越重要的研究课题，其中，如何分析各种地震属性的油气地质意义是至关重要的。

在现有技术中，早期的地震解释主要是解释反射层位并将其转换为深度剖面以确定地下的地质构造；1960年代下半叶产生的油气勘探亮点技术改变了地球物理工作者对地震解释的看法，是地震属性分析技术发展的基石。

1980年代中叶和1990年代初，出现了层位属性、层间属性以及从三维数据体中沿某一平面抽取出来的属性等技术方法。这些技术从计算的角度将地震属性分为两类：一类是单道计算的地震属性；一类是多道计算的地震属性。并由此归纳出瞬时振幅、瞬时相位、瞬时频率、相干数据体、波形聚类五种基本的地震属性。本世纪初，出现了由地震属性向储层参数转换的综合效果分析，就是将提取优化后的地震属性转换为各种储层参数如岩性、砂体厚度、泥质含量、孔隙度、渗透率、流体饱和度等等。但从岩石物理学的角度看，在储层参数和地震属性之间并不存在直接的解析关系，即不能用确定的函数表达式进行描述。因此针对多元逐步回归、相关滤波、神经网络、协克里金、非参数回归分析等方法，人们往往也仅凭经验来选择使用某一种预测方法。

发明人发现，目前的属性分析方法通常只表达了一维信息，比如振幅等，虽然波形聚类具有多维信息，但也只是以一个参数的方式来体现，具有地质模糊性。因此，尽管地震属性种类多样，但其本质只是对地震数据某个方面某种结构模式的一种刻画，并没有明确的地质含义。因此，有必要开发一种有效判识地震属性的油气地质意义的地震属性分析方法。

公开于本公开背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本公开提出了一种多参数地震属性分析方法及系统，该方法通过确定地质异常体，获得地质异常体的多个参数，并对所获得的多个参数进行地质解析，获得油气层空间展布特征，从而有效判识地震属性的油气地质意义。

根据本公开的一方面，提出了一种多参数地震属性分析方法，该方法可以包括以下步骤：基于敏感地震属性的响应特征确定地质异常体；基于地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数；以及基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征。

根据本公开的另一方面，提出了一种多参数地震属性分析系统，该系统可以包括以下单元：用于基于敏感地震属性的响应特征确定地质异常体的单元；用于基于地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数的单元；以及用于基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征的单元。

本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。

附图说明

通过结合附图对本公开示例性实施例进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施例中，相同的参考标号通常代表相同部件。

图1示出了根据本公开的示例性实施例的多参数地震属性分析方法的步骤的流程图。

图2示出了根据本公开的多参数地震属性分析方法的应用示例的示意图。

具体实施方式

下面将参照附图更详细地描述本公开的优选实施例。虽然附图中显示了本公开的优选实施例，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。

实施例1

图1示出了根据本公开的示例性实施例的多参数地震属性分析方法的步骤的流程图。该方法可以包括：步骤101，基于敏感地震属性的响应特征确定地质异常体；步骤102，基于地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数；以及步骤103，基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征。

本实施例通过确定地质异常体，获得地质异常体的多个参数，并对所获得的多个参数进行地质解析，获得油气层空间展布特征，从而有效判识地震属性的油气地质意义。

确定地质异常体

在一个示例中，可以基于敏感地震属性的响应特征确定地质异常体。

由地震反射特征分析可以得知，地震响应可以刻画油气层沉积特征。在地震层位精细解释的基础上，可以针对地质目标层段确定相对稳定的沉积界面(反射界面)。针对实钻测井分析并结合地震资料信息，可以确定油气藏的敏感地震属性。因此，通过井震精细标定，可以明确对某油气层敏感的敏感地震属性的响应特征，进而可以确定地质异常体。

获得地质异常体的参数

在一个示例中，可以基于地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数。

对于岩性隐蔽油气藏而言，由于地质体沉积的复杂性，往往油气层(特别是油气田开发阶段的含油气小层)地震响应并不严格与稳定地层界面平行，因此采用简单时窗分析技术分析该油气层的敏感地震属性十分困难。

本发明可以从三维空间出发刻画敏感地震属性的地质异常体。首先，可以确定目标地质异常体的分布时间范围，确定目标地质异常体的精确位置。在此基础上，通过大量实钻测井试气结果及试气层段，可以进一步确定地质异常体的三个参数：平均值参数、时间厚度参数及时差参数。

平均值参数是指能够表征地质异常体的强度平均值的参数，时间厚度参数是指能够以时间的形式表征地质异常体的厚度的参数，时差参数是指能够以时间差的形式表征地质异常体距离地质界面的距离的参数。

获得油气层空间展布特征

在一个示例中，可以基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征。

对于所获得的地质异常体的平均值参数、时间厚度参数及时差参数，可以从地震资料解释视角进行解析，针对地质异常体进行精细刻画。

具体地，平均值参数可以表征地质异常体的强度信息，也即地质异常体含油气性的强弱程度。时间厚度参数可以表征地质异常体的厚度信息，也即地质异常体的含油气层的相对厚度，可以反映地质异常体的大小。时差参数可以表征地质异常体距离地质界面的时差，也即可以确定该地质异常体代表的是哪套油气小层的响应特征。

如上所述，通过地质异常体的平均值参数、时间厚度参数及时差参数这三个参数的定量化解释，可以借助含油气敏感地震属性立体地有效刻画油气层的空间展布特征。相对于现有技术的定性化地震属性分析，本公开可以对地质异常体的多个参数进行定量化的地震属性分析。

应用示例

为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。

在鄂尔多斯盆地大牛地气田下石盒子组盒1段下亚段气层发育，通过研究发现盒1段气层具有负梯度、负截距的属性，针对盒1段气层可以应用AVO负截距、负梯度属性预测盒1段气层分布。

如何从异常剖面上区分出盒1段下亚段储层的有利区范围，是我们要解决的难题。为此我们引入了AVO三参数分析技术，通过大量井震信息精细分析，可以获得气层有利响应的平均值、时间厚度、时差等三参数，综合应用这三个参数可以更有效判识盒1段下亚段气层。

统计大量钻井结果分析，可以认为盒1段下亚段气层的响应为：AVO负截距负梯度属性要距离盒1段地质层位时差小于20ms，且负截距负梯度AVO异常平均值小于-3000，同时AVO负截距负梯度异常厚度要大于4ms。其中，AVO负截距负梯度属性距离盒1段地质层位时差可以作为时差参数的一个示例，负截距负梯度AVO异常平均值可以作为平均值参数的一个示例，AVO负截距负梯度异常厚度可以作为时间厚度参数的一个示例。利用此三参数，可以刻画大牛地气田盒1段下亚段气层发育区。

图2示出了根据本公开的多参数地震属性分析方法的应用示例的示意图。图2采用了盒1段AVO三参数分析技术，在图2中，过A井和B井的AVO负截距负梯度异常均大于4毫秒，且异常值小于-3000的剖面，剖面上AB井点处显示的为含气饱和度曲线，由B井可以看到在目的层段的异常分布范围较大(黑色线圈出范围)，距离盒1底部的T9d界面为10毫秒，含气饱和度曲线显示好；而A井异常距离T9d界面大于20毫秒，且测井含气饱和度曲线显示不好，异常刻画与实钻情况吻合。

本领域技术人员应理解，上面对本公开的实施例的描述的目的仅为了示例性地说明本公开的实施例的有益效果，并不意在将本公开的实施例限制于所给出的任何示例。

实施例2

在该实施例中，提供了一种多参数地震属性分析系统，该系统可以包括以下单元：用于基于敏感地震属性的响应特征确定地质异常体的单元；用于基于所述地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数的单元；以及用于基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征的单元。

本实施例通过确定地质异常体，获得地质异常体的多个参数，并对所获得的多个参数进行地质解析，获得油气层空间展布特征，从而有效判识地震属性的油气地质意义。

在一个示例中，平均值参数可以表征地质异常体的强度信息。

在一个示例中，时间厚度参数可以表征地质异常体的厚度信息。

在一个示例中，时差参数可以表征所述地质异常体距离地质界面的时差。

本领域技术人员应理解，上面对本公开的实施例的描述的目的仅为了示例性地说明本公开的实施例的有益效果，并不意在将本公开的实施例限制于所给出的任何示例。

本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。

计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。

这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。

用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。

这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。

这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。

也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。

附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (8)

1.一种多参数地震属性分析方法，所述方法包括以下步骤：

基于敏感地震属性的响应特征确定地质异常体；

基于所述地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数；以及

基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征。

2.根据权利要求1所述的多参数地震属性分析方法，其中，所述平均值参数表征所述地质异常体的强度信息。

3.根据权利要求1所述的多参数地震属性分析方法，其中，所述时间厚度参数表征所述地质异常体的厚度信息。

4.根据权利要求1所述的多参数地震属性分析方法，其中，所述时差参数表征所述地质异常体距离地质界面的时差。

5.一种多参数地震属性分析系统，所述系统包括以下单元：

用于基于敏感地震属性的响应特征确定地质异常体的单元；

用于基于所述地质异常体获得地质异常体的平均值参数、时间厚度参数及时差参数的单元；以及

用于基于所获得的平均值参数、时间厚度参数及时差参数进行地质解析，获得油气层空间展布特征的单元。

6.根据权利要求5所述的多参数地震属性分析系统，其中，所述平均值参数表征所述地质异常体的强度信息。

7.根据权利要求5所述的多参数地震属性分析系统，其中，所述时间厚度参数表征所述地质异常体的厚度信息。

8.根据权利要求5所述的多参数地震属性分析系统，其中，所述时差参数表征所述地质异常体距离地质界面的时差。