CN105974473B - 一种确定薄储层的地震属性的方法和装置 - Google Patents

Abstract

本申请实施例提供一种确定薄储层的地震属性的方法和装置。该方法包括：根据获取的目的层的主频和层速度计算得到所述目的层的波长；获取所述目的层的层厚度；判断所述目的层的层厚度是否小于所述目的层的波长的四分之一；当所述判断的结果为是时，确定所述目的层为薄储层；计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。利用本申请实施例提供的技术方案有效的确定储层的发育区。与现有技术相比，利用本申请提供的技术方案可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测出储层的分布特征。

Description

一种确定薄储层的地震属性的方法和装置

技术领域

本发明涉及油气勘探技术领域，尤其涉及一种确定薄储层的地震属性的方法和装置。

背景技术

地震属性指的是那些由叠前或叠后地震数据，经过数学变换而得出的有关地震波的几何学、运动学、动力学或统计学特征。地震属性主要有时间、振幅、频率、相位、相干和衰减等几大类。在油气勘探技术领域，地震属性主要应用于油藏储层厚度、孔隙度、渗透率等储层的分布特征预测，其中地震资料的振幅、频率等单项地震属性的应用比较广泛。

现有技术中为了准确反映出储层的分布特征，需要放大储层的振幅、频率的特征；目前应用较广泛的甜心体属性是振幅与频率平方根的商，但甜心体属性只能适应于具有强振幅、低频的特征的储层。对于薄储层(厚度小于四分之一波长)情况下，由于干涉作用，储层表现为强振幅、高频的特征，频率越高，振幅与频率平方根的商越小，导致甜心体属性(振幅与频率平方根的商)无法有效反映薄储层的分布特征。

因此，现有技术亟需一种确定薄储层的地震属性的方法，可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测储层的分布特征。

发明内容

本申请的目的是提供一种确定薄储层的地震属性的方法和装置，可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测储层的分布特征。

本申请提供的确定薄储层的地震属性的方法和装置是这样实现的：

一种确定薄储层的地震属性的方法，所述方法包括：

根据获取的目的层的主频和层速度计算得到所述目的层的波长；

获取所述目的层的层厚度；

判断所述目的层的层厚度是否小于所述目的层的波长的四分之一；

当所述判断的结果为是时，确定所述目的层为薄储层；

计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。

在一个优选的实施例中，所述方法还包括：

在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；

利用所述均方根振幅对所述薄储层进行储层分布特征分析处理。

在一个优选的实施例中，所述获取的目的层的主频和层速度包括：

利用提取的所述目的层的振幅频谱确定出的所述目的层的主频，以及利用声波时差测井曲线确定出的所述目的层的层速度。

在一个优选的实施例中，至少采用下述方式中的一种获取所述目的层的层厚度：

利用测井曲线确定出所述目的层的层厚度；

利用录井图确定出所述目的层的层厚度。

在一个优选的实施例中，所述获取的所述目的层瞬时振幅数据体和瞬时主频数据体包括：

利用所述目的层的地震数据计算得到的瞬时振幅数据体和瞬时主频数据体。

一种确定薄储层的地震属性的装置，所述装置包括：

波长计算模块，用于根据获取的目的层的主频和层速度计算得到所述目的层的波长；

厚度获取模块，用于获取所述目的层的层厚度；

判断模块，用于判断所述目的层的层厚度是否小于所述目的层的波长的四分之一；

薄储层确定模块，当所述判断模块判断的结果为是时，确定所述目的层为薄储层；

地震属性计算模块，用于计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。

在一个优选的实施例中，所述装置还包括：

均方根振幅提取模块，用于在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；

储层分布特征分析处理模块，用于利用所述均方根振幅对所述薄储层进行储层分布特征分析处理。

在一个优选的实施例中，所述获取的目的层的主频和层速度包括：

利用提取的所述目的层的振幅频谱确定出的所述目的层的主频，以及利用声波时差测井曲线确定出的所述目的层的层速度。

在一个优选的实施例中，所述厚度获取模块至少包括下述单元之一：

第一厚度获取单元，用于利用测井曲线确定出所述目的层的层厚度；

第二厚度获取单元，用于利用录井图确定出所述目的层的层厚度。

在一个优选的实施例中，所述获取的所述目的层瞬时振幅数据体和瞬时主频数据体包括：

利用所述目的层的地震数据计算得到的瞬时振幅数据体和瞬时主频数据体。

本申请通过判断所述目的层的层厚度与所述目的层的波长的四分之一之间的大小，可以将目的层的层厚度小于目的层的波长的四分之一所对应的所述目的层确定为薄储层；由于薄储层表现为强振幅、高频的特征，可以计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果。与现有技术相比，利用本申请提供的技术方案可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测储层的分布特征。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本申请提供的确定薄储层的地震属性的方法的一种实施例的流程图；

图2是本申请提供的一种瞬时振幅数据体的剖面示意图；

图3是本申请提供的一种瞬时主频数据体的剖面示意图；

图4是本申请提供的一种时振幅数据体和瞬时主频数据体的乘积的结果数据体的剖面示意图；

图5是本申请提供的确定薄储层的地震属性的方法的另一种实施例的流程图；

图6是本申请提供的一种在所述结果数据体中提取的均方根振幅的平面示意图；

图7是本申请提供的确定薄储层的地震属性的装置的一种实施例中的结构示意图；

图8是本申请提供的确定薄储层的地震属性的装置的另一种实施例中的结构示意图。

具体实施方式

为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。

下面以几个具体的例子详细说明本申请实施例的具体实现。

以下首先介绍本申请一种确定薄储层的地震属性的方法的一种实施例。图1是本申请提供的确定薄储层的地震属性的方法的一种实施例的流程图，本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。具体的如图1所示，所述方法可以包括：

S110：根据获取的目的层的主频和层速度计算得到所述目的层的波长。

本申请实施例中，可以根据获取的目的层的主频和层速度计算得到所述目的层的波长。具体的，可以将所述目的层的层速度除以所述目的层的主频，将商作为所述目的层的波长。

具体的，所述获取的目的层的主频和层速度可以包括：

利用提取的所述目的层的振幅频谱确定出的所述目的层的主频，以及利用声波时差测井曲线确定出的所述目的层的层速度。

具体的，在实际应用中，所述目的层可以采用下述方式确定：利用测井资料中的声波时差测井曲线、密度曲线、地震子波制作合成地震记录；然后，利用合成地震记录、钻井分层数据、地震数据进行层位标定，将钻井的目的层的分层在地震数据上标定的位置确定为目的层。

具体的，所述目的层的振幅频谱可以包括在地震数据中，沿着预先拾取的层位以至少一个子波长度为时窗提取的振幅频谱。相应的，所述预先拾取的层位可以包括在地震数据中拾取所述目的层的主要层位。在实际应用中，拾取密度可以为8x8至4x4，本申请实施例并不以此为限。

在一个具体的实施例中，假设所述目的层的主频为30HZ，所述目的层的层速度为5000m/s，相应的，可以计算出所述目的层的波长为167m。

S120：获取所述目的层的层厚度。

本申请实施例中，可以获取所述目的层的层厚度。具体的，可以至少采用下述方式中的一种获取所述目的层的层厚度：

利用测井曲线确定出所述目的层的层厚度；

利用录井图确定出所述目的层的层厚度。

具体的，所述测井曲线和录井图可以从测井过程中的测井资料中获取。

S130：判断所述目的层的层厚度是否小于所述目的层的波长的四分之一。

本申请实施例中，在获取所述目的层的波长和厚度之后，可以判断所述目的层的层厚度是否小于所述目的层的波长的四分之一。

S140：当所述判断的结果为是时，确定所述目的层为薄储层。

本申请实施例中，当步骤S130判断的结果为是时，可以确定所述目的层为薄储层。在一个具体的实施例中，假设所述目的层的层厚度为25m，结合上述计算出的所述目的层的波长167m，可以计算出所述目的层的波长167m的四分之一为41.75m，判断出25m<41.75m。相应的，可以判断出所述目的层的层厚度小于所述目的层的波长的四分之一。因此，可以确定所述目的层为薄储层。

S150：计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。

本申请实施例中，当确定所述目的层为薄储层之后，可以计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。由于薄储层表现为强振幅、高频的特征，这里将目的层瞬时振幅数据体和瞬时主频数据体的乘积的结果数据体作为所述薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果。

具体的，所述获取的所述目的层瞬时振幅数据体和瞬时主频数据体可以包括：

利用所述目的层的地震数据计算得到的瞬时振幅数据体和瞬时主频数据体。

具体的，在实际应用中，所述利用所述目的层的地震数据计算得到的瞬时振幅数据体和瞬时主频数据体可以结合相应的计算软件，本申请实施例中所述计算软件可以包括但不限于GeoEast。具体的，可以选择目的层的地震数据为输入的数据体，解释目的层对应的层位的顶、底，将层位的顶、底作为计算时窗，根据预设的瞬时振幅算法和预设的瞬时主频算法分别计算得到瞬时振幅数据体和瞬时主频数据体。

具体的，如图2所示，图2是本申请提供的一种瞬时振幅数据体的剖面示意图。其中，横坐标表示地震线道号，其中424、413表示道号，539、639表示线号；纵坐标表示时间，单位为s；Lanea SE-1和Lanea SE-2表示测井的井号；色标210表示地震记录的振幅。

具体的，如图3所示，图3是本申请提供的一种瞬时主频数据体的剖面示意图。其中，横坐标表示地震线道号，其中424、413表示道号，539、639表示线号；纵坐标表示时间，单位为s；Lanea SE-1和Lanea SE-2表示测井的井号；色标310表示地震记录的频率，单位为Hz。

具体的，如图4所示，结合上述图2对应的瞬时振幅数据体和图3对应的瞬时主频数据体，图4是本申请提供的一种时振幅数据体和瞬时主频数据体的乘积的结果数据体的剖面示意图。其中，横坐标表示地震线道号，其中424、413表示道号，539、639表示线号；纵坐标表示时间，单位为s；Lanea SE-1和Lanea SE-2表示测井的井号；色标410表示地震记录的频率与振幅的乘积。

由以上本申请一种确定薄储层的地震属性的方法的实施例可见，本申请通过判断所述目的层的层厚度与所述目的层的波长的四分之一之间的大小，可以将目的层的层厚度小于目的层的波长的四分之一所对应的所述目的层确定为薄储层；由于薄储层表现为强振幅、高频的特征，可以计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果。与现有技术相比，利用本申请提供的技术方案可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测储层的分布特征。

以下介绍本申请一种确定薄储层的地震属性的方法的另一种实施例。图5是本申请提供的确定薄储层的地震属性的方法的另一种实施例的流程图，本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的系统或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。具体的如图5所示，所述方法可以包括：

S510：根据获取的目的层的主频和层速度计算得到所述目的层的波长。

S520：获取所述目的层的层厚度。

S530：判断所述目的层的层厚度是否小于所述目的层的波长的四分之一。

S540：当所述判断的结果为是时，确定所述目的层为薄储层。

S550：计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。

S560：在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅。

S570：利用所述均方根振幅对所述薄储层进行储层分布特征分析处理。

本申请实施例中，在步骤S560提取均方根振幅之后，可以利用所述均方根振幅对所述薄储层进行储层分布特征分析处理，确定出储层的发育区。具体的，结合上述图4中对应的结果数据体，如图6所示，图6是本申请提供的一种在所述结果数据体中提取的均方根振幅的平面示意图。图6中包括井号为Lanea E-1、Lanea E-2、Lanea E-3、Lanea E-4、Laneajy-1、Lanea SE-1和Lanea SE-2的测井对应的区域，从图中可见，610和620所对应的区域为储层的发育区，其中610对应井号为Lanea E-2的测井所对应的区域，620对应井号为LaneaSE-1的测井所对应的区域；色标630表示均方根振幅。

由以上本申请一种确定薄储层的地震属性的方法的实施例可见，本申请通过判断所述目的层的层厚度与所述目的层的波长的四分之一之间的大小，可以将目的层的层厚度小于目的层的波长的四分之一所对应的所述目的层确定为薄储层；由于薄储层表现为强振幅、高频的特征，可以计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果；接着，可以在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；利用所述均方根振幅对所述薄储层进行储层分布特征分析处理，可以有效的确定储层的发育区。与现有技术相比，利用本申请提供的技术方案可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测出储层的分布特征。

基于本申请所述的一种确定薄储层的地震属性的方法，本申请另一方面还提供一种确定薄储层的地震属性的装置，图7是本申请提供的确定薄储层的地震属性的装置的一种实施例中的结构示意图；如图7所示，所述装置700可以包括：

波长计算模块710，可以用于根据获取的目的层的主频和层速度计算得到所述目的层的波长；

厚度获取模块720，可以用于获取所述目的层的层厚度；

判断模块730，可以用于判断所述目的层的层厚度是否小于所述目的层的波长的四分之一；

薄储层确定模块740，可以当所述判断模块730判断的结果为是时，确定所述目的层为薄储层；

地震属性计算模块750，可以用于计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性。

本申请实施例中，判断模块730可以根据波长计算模块710计算得到的目的层的波长和厚度获取模块720获取的目的层的层厚度，可以判断出所述目的层的层厚度是否小于所述目的层的波长的四分之一；当所述判断模块730判断的结果为是时，薄储层确定模块740可以确定所述目的层为薄储层；最后，地震属性计算模块750可以计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性，可以实现对薄储层的强振幅、高频的特征放大效果。

图8是本申请提供的确定薄储层的地震属性的装置的另一种实施例中的结构示意图；另一个实施例中，如图8所示，所述装置700还可以包括：

均方根振幅提取模块760，可以用于在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；

储层分布特征分析处理模块770，可以用于利用所述均方根振幅对所述薄储层进行储层分布特征分析处理。

在一个优选的实施例中，所述获取的目的层的主频和层速度可以包括：

利用提取的所述目的层的振幅频谱确定出的所述目的层的主频，以及利用声波时差测井曲线确定出的所述目的层的层速度。

在一个优选的实施例中，所述厚度获取模块可以至少包括下述单元之一：

第一厚度获取单元，可以用于利用测井曲线确定出所述目的层的层厚度；

第二厚度获取单元，可以用于利用录井图确定出所述目的层的层厚度。

在一个优选的实施例中，所述获取的所述目的层瞬时振幅数据体和瞬时主频数据体可以包括：

利用所述目的层的地震数据计算得到的瞬时振幅数据体和瞬时主频数据体。

由以上本申请一种确定薄储层的地震属性的方法或装置的实施例可见，本申请通过判断所述目的层的层厚度与所述目的层的波长的四分之一之间的大小，可以将目的层的层厚度小于目的层的波长的四分之一所对应的所述目的层确定为薄储层；由于薄储层表现为强振幅、高频的特征，可以计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果；接着，可以在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；利用所述均方根振幅对所述薄储层进行储层分布特征分析处理，可以有效的确定储层的发育区。与现有技术相比，利用本申请提供的技术方案可以有效的确定出薄储层的地震属性，实现对薄储层的强振幅、高频的特征放大效果，从而更好地预测出储层的分布特征。

本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (8)

1.一种确定薄储层的地震属性的方法，其特征在于，所述方法包括：

根据获取的目的层的主频和层速度计算得到所述目的层的波长；

获取所述目的层的层厚度；

判断所述目的层的层厚度是否小于所述目的层的波长的四分之一；

当所述判断的结果为是时，确定所述目的层为薄储层；

计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性；

其中，所述获取的所述目的层瞬时振幅数据体和瞬时主频数据体包括采用下述方式获取：

结合相应的计算软件，选择目的层的地震数据为输入的数据体，解释目的层对应的层位的顶、底，将层位的顶、底作为计算时窗，根据预设的瞬时振幅算法和预设的瞬时主频算法分别计算得到瞬时振幅数据体和瞬时主频数据体，所述计算软件包括GeoEast。

2.根据权利要求1所述的方法，其特征在于，所述方法还包括：

在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；

利用所述均方根振幅对所述薄储层进行储层分布特征分析处理。

3.根据权利要求1或2任意一项所述的方法，其特征在于，所述获取的目的层的主频和层速度包括：

利用提取的所述目的层的振幅频谱确定出的所述目的层的主频，以及利用声波时差测井曲线确定出的所述目的层的层速度。

4.根据权利要求1或2任意一项所述的方法，其特征在于，至少采用下述方式中的一种获取所述目的层的层厚度：

利用测井曲线确定出所述目的层的层厚度；

利用录井图确定出所述目的层的层厚度。

5.一种确定薄储层的地震属性的装置，其特征在于，所述装置包括：

波长计算模块，用于根据获取的目的层的主频和层速度计算得到所述目的层的波长；

厚度获取模块，用于获取所述目的层的层厚度；

判断模块，用于判断所述目的层的层厚度是否小于所述目的层的波长的四分之一；

薄储层确定模块，当所述判断模块判断的结果为是时，确定所述目的层为薄储层；

地震属性计算模块，用于计算获取的所述目的层瞬时振幅数据体和瞬时主频数据体的乘积，将所述乘积的结果数据体作为所述薄储层的地震属性；

其中，所述获取的所述目的层瞬时振幅数据体和瞬时主频数据体包括采用下述方式获取：

结合相应的计算软件，选择目的层的地震数据为输入的数据体，解释目的层对应的层位的顶、底，将层位的顶、底作为计算时窗，根据预设的瞬时振幅算法和预设的瞬时主频算法分别计算得到瞬时振幅数据体和瞬时主频数据体，所述计算软件包括GeoEast。

6.根据权利要求5所述的装置，其特征在于，所述装置还包括：

均方根振幅提取模块，用于在所述结果数据体中沿预先拾取的层位以至少一个子波长度为时窗提取均方根振幅；

储层分布特征分析处理模块，用于利用所述均方根振幅对所述薄储层进行储层分布特征分析处理。

7.根据权利要求5或6任意一项所述的装置，其特征在于，所述获取的目的层的主频和层速度包括：

利用提取的所述目的层的振幅频谱确定出的所述目的层的主频，以及利用声波时差测井曲线确定出的所述目的层的层速度。

8.根据权利要求5或6任意一项所述的装置，其特征在于，所述厚度获取模块至少包括下述单元之一：

第一厚度获取单元，用于利用测井曲线确定出所述目的层的层厚度；

第二厚度获取单元，用于利用录井图确定出所述目的层的层厚度。