CN113625340B - 一种基于属性敏感度选择储层敏感地震属性组合的方法 - Google Patents

Abstract

本发明涉及一种基于属性敏感度选择敏感属性组合的方法，包括步骤：提取储层中潜在的多种类型储层地震属性；计算所述第一交会椭圆得到地震属性敏感度S1；根据所述地震属性敏感度S1从大到小选出若干种储层地震属性作为储层敏感地震属性集；计算对应的两种储层敏感地震属性之间的地震属性敏感度S2，并根据所述的地震属性敏感度S2得到地震属性敏感度矩阵；计算所述地震属性敏感度矩阵每行的敏感度平均值，作为该行所对应的储层敏感地震属性与其他储层敏感属性之间的平均相关程度，选择平均相关程度较小的多个储层敏感地震属性作为最终的储层敏感地震属性组合。述选择敏感属性组合的方法能够提高对储层结构预测的准确性。

Description

一种基于属性敏感度选择储层敏感地震属性组合的方法

技术领域

本发明涉及油气田开发技术领域，尤其涉及一种基于属性敏感度选择储层敏感地震属性组合的方法。

背景技术

地震属性分析是油气田开发阶段储层预测常用的一种技术手段，通过沿着解释的砂体顶、底层面提取振幅、频率、波形等各种类别的地震属性，就能根据属性值在平面上呈现出的分布形态，实现储层展布范围刻画、砂体优势发育区域厘定以及沉积微相划分等储层研究工作。目前，常采用多种地震属性作为组合，通过多元线性回归或非线性神经网络算法形成一种最终的复合属性来预测储层，即多属性储层预测。而选用哪些储层敏感地震属性形成储层敏感属性组合，则直接决定了最终预测结果是否与井上实钻结果相吻合，以及属性趋势是否符合地质沉积模式认识。

而现有技术并没有可行的方法能够从若干种属性组合中选择出与储层参数的相关程度最高，同时包含冗余信息最少的敏感属性组合，因此导致储层预测结果存在可靠性不高的问题。

发明内容

针对上述问题，本发明的目的是提供一种基于属性敏感度选择储层敏感地震属性组合的方法、装置、介质及设备，能够选出的与储层参数的相关程度高，且包含冗余信息少的多种敏感地震属性组合，进一步提高储层预测结果的准确性。

为实现上述目的，本发明采取以下技术方案：

一种基于属性敏感度选择储层敏感地震属性组合的方法，包括步骤：

提取储层中潜在的多种类型储层地震属性；

计算所述多种类型储层地震属性的地震属性敏感度S1；

根据所述地震属性敏感度S1从大到小选出若干种储层地震属性作为储层敏感地震属性集；

计算所述储层敏感地震属性集中的任意两种储层敏感地震属性之间的地震属性敏感度S2，并根据所述的地震属性敏感度S2得到地震属性敏感度矩阵；

计算所述地震属性敏感度矩阵每行的敏感度平均值，作为该行所对应的储层敏感地震属性与其他储层敏感属性之间的平均相关程度，选择平均相关程度较小的多个储层敏感地震属性作为最终的储层敏感地震属性组合。

进一步地，计算所述多种类型储层地震属性的地震属性敏感度S1包括：制作所述多种类型储层地震属性与待预测的储层参数的交会图，并根据所述交会图制作第一交会椭圆，测量所述第一交会椭圆的长轴长度R1与短轴长度R2，通过所述长轴长度R1除以短轴长度R2得到地震属性敏感度S1。

进一步地，所述第一交会椭圆的制作方法具体为：

先绘制标准圆形以能够尽可能包含储层地震属性与待预测的储层参数之间的交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第一交会椭圆。

进一步地，计算对应的任意两种储层敏感地震属性之间的地震属性敏感度S2包括：对所述储层敏感地震属性集中的任意两种储层敏感地震属性进行两两交会得到第二交会椭圆，并将第二交会椭圆的长轴长度r1除以短轴长度r2得到的值作为两种储层敏感地震属性之间的属性敏感度S2。

进一步地，所述第二交会椭圆的制作方法包括：

先绘制标准圆形以能够尽可能包含两种储层敏感地震属性之间的交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第二交会椭圆。

进一步地，采用如下方式计算属性敏感度矩阵：

对所述的若干种储层敏感地震属性进行编号，具体为1,2,…,n，其中n为储层敏感地震属性总数，制作n行n列的矩阵，其中第i行第j列所对应的值为第i个储层敏感地震属性与第j个储层敏感地震属性组合的属性敏感度，获得属性敏感度矩阵，其中，1<＝i、j<＝n；

该属性敏感度矩阵中每行元素的敏感度平均值代表该行对应储层敏感地震属性与其他行对应的储层敏感地震属性的平均相关程度，所述平均相关程度越小说明该行所对应的储层敏感属性包含的冗余信息越少。

进一步地，提取储层中潜在的多种类型储层地震属性，具体为：

提取总负振幅、最小振幅、弧长、平均瞬时相位、平均瞬时频率共5种类型储层地震属性。

一种基于属性敏感度选择储层敏感地震属性组合的装置，包括：

第一处理装置，用于提取储层中潜在的多种类型储层地震属性；

第二处理装置，用于计算所述多种类型储层地震属性的地震属性敏感度S1；

第三处理装置，用于根据所述地震属性敏感度S1从大到小选出若干种储层地震属性作为储层敏感地震属性集；

第四处理装置，用于计算任意两种储层敏感地震属性之间的地震属性敏感度S2，并根据所述的地震属性敏感度S2得到地震属性敏感度矩阵；

第五处理装置，用于计算所述地震属性敏感度矩阵每行的敏感度平均值，作为该行所对应的储层敏感地震属性与其他储层敏感属性之间的平均相关程度，选择平均相关程度较小的多个储层敏感地震属性作为最终的储层敏感地震属性组合。

一种计算机可读存储介质，存储有计算机程序，所述计算机程序用于被处理器执行时实现所述的选择储层敏感地震属性组合的方法的步骤。

一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述的选择储层敏感地震属性组合的方法的步骤。

本发明由于采取以上技术方案，其具有以下优点：采用本发明提供的选择敏感属性组合的方法能够优选出与储层参数(如砂体厚度、砂地比等)相关程度最高的多种地震属性，同时各属性之间保持相互独立、即尽可能少包含冗余信息，进而提高对储层结构预测的准确性。

附图说明

图1是本发明实施的提供的选择储层敏感地震属性组合的方法的流程图；

图2(a)-(e)分别表示实施例中提取的某海上油田某砂体顶、底层面之间总负振幅、最小振幅、弧长、平均瞬时相位、平均瞬时频率共5种地震属性；

图3(a)-(e)分别表示总负振幅、最小振幅、弧长、平均瞬时相位、平均瞬时频率分别5中地震属性与井上钻遇砂体厚度的交会图与交会椭圆；

图4(a)-(f)分别表示总负振幅、最小振幅、弧长、平均瞬时相位4种储层敏感属性两两交会的得到的交会椭圆，以及任意两种属性之间的属性敏感度；

图5是基于属性敏感度计算的属性敏感度矩阵，其中第i行第j列元素值代表第i种储层敏感属性与第j种储层敏感属性组合的属性敏感度；

具体实施方式

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

在本发明的描述中，需要说明的是，术语“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的系统或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，使用术语“第一”、“第二”等词语来限定零部件，仅仅是为了便于对上述零部件进行区别，如没有另行声明，上述词语并没有特殊含义，不能理解为指示或暗示相对重要性。

在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

如图1所示，以渤海某油田某砂体储层预测为例，采用本发明方法优选出3种储层敏感地震属性，形成储层敏感地震属性组合，具体包括以下步骤：

1)提取总负振幅、最小振幅、弧长、平均瞬时相位、平均瞬时频率共5种储层地震属性，如图2(a)-(e)所示；

2)制作上述5种储层地震属性与井上钻遇砂体厚度的交会图，并根据所述交会图绘制第一交会椭圆，如图3(a)-(e)所示，根据第一交会椭圆上长轴长度R1与短轴长度R2之比求得地震属性敏感度S1＝R1/R2，即长轴长度除以短轴长度；

其中，所述第一交会椭圆的制作方法具体为：

先绘制标准圆形以能够尽可能包含主要交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第一交会椭圆。

3)根据地震属性敏感度从大到小的原则选出了4种储层敏感地震属性：总负振幅、最小振幅、弧长、平均瞬时频率作为储层敏感地震属性；

4)对所述储层敏感地震属性两两进行交会，得到共6种交会图，分别为：最小振幅-总负振幅交会、弧长-总负振幅交会、平均瞬时频率-总负振幅交会、弧长-最小振幅交会、平均瞬时频率-最小振幅交会、平均瞬时频率-弧长交会，计算任意两种敏感地震属性组合的属性敏感度S2，如图4(a)-(f)所示；

该步骤中计算任意两种敏感地震属性的属性敏感度S2与步骤2)基本相同，具体为：

对储层敏感地震属性两两进行交会，得到第二交会椭圆，并将第二交会椭圆的长轴长度r1除以短轴长度r2得到的值为两种地震敏感属性组合的属性敏感度S2。

5)根据上述任意两种敏感地震属性组合的属性敏感度S2得到属性敏感度矩阵，如图5所示，矩阵中第i行第j列元素代表的是总负振幅、最小振幅、弧长、平均瞬时频率这4种属性中第i种储层敏感属性与第j种储层敏感地震属性组合的属性敏感度。

6)计算属性敏感度矩阵每行的平均值，其中第i行的平均值作为总负振幅、最小振幅、弧长、平均瞬时频率这4种属性中第i个储层敏感地震属性与其它敏感地震属性的平均相关程度。由此计算出总负振幅属性与其他敏感地震属性的平均相关程度为2.2，最小振幅属性与其他敏感地震属性的平均相关程度为2.5，弧长属性与其他敏感地震属性的平均相关程度为2.4，平均瞬时频率属性与其它敏感地震属性的平均相关程度为1.9。该矩阵每行元素的平均值代表了该行的储层敏感地震属性与其他敏感地震属性的平均相关程度，相关程度越小说明该储层敏感地震属性包含的冗余信息越少。因此，选取平均相关程度最小的三种属性作为最终优选出的储层敏感属性组合，即总负振幅、弧长、平均瞬时频率。

采用本发明提供的选择敏感属性组合的方法能够优选出与储层参数(如砂体厚度、砂地比等)相关程度最高的多种地震属性，同时各属性之间保持相互独立、即尽可能少包含冗余信息，进而提高对储层结构预测的准确性。

最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (5)

1.一种基于属性敏感度选择储层敏感地震属性组合的方法，其特征在于，包括步骤：

提取储层中潜在的多种类型储层地震属性；

计算所述多种类型储层地震属性的地震属性敏感度S1；

根据所述地震属性敏感度S1从大到小选出若干种储层地震属性作为储层敏感地震属性集；

计算所述储层敏感地震属性集中的任意两种储层敏感地震属性之间的地震属性敏感度S2，并根据所述的地震属性敏感度S2得到地震属性敏感度矩阵；

计算所述地震属性敏感度矩阵每行的敏感度平均值，作为该行所对应的储层敏感地震属性与其他储层敏感属性之间的平均相关程度，选择平均相关程度较小的多个储层敏感地震属性作为最终的储层敏感地震属性组合；

计算所述多种类型储层地震属性的地震属性敏感度S1包括：制作所述多种类型储层地震属性与待预测的储层参数的交会图，并根据所述交会图制作第一交会椭圆，测量所述第一交会椭圆的长轴长度R1与短轴长度R2，通过所述长轴长度R1除以短轴长度R2得到地震属性敏感度S1；

所述第一交会椭圆的制作方法具体为：

先绘制标准圆形以能够尽可能包含储层地震属性与待预测的储层参数之间的交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第一交会椭圆；

计算对应的任意两种储层敏感地震属性之间的地震属性敏感度S2包括：对所述储层敏感地震属性集中的任意两种储层敏感地震属性进行两两交会得到第二交会椭圆，并将第二交会椭圆的长轴长度r1除以短轴长度r2得到的值作为两种储层敏感地震属性之间的属性敏感度S2；

所述第二交会椭圆的制作方法包括：

先绘制标准圆形以能够尽可能包含两种储层敏感地震属性之间的交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第二交会椭圆；

提取储层中潜在的多种类型储层地震属性，具体为：

提取总负振幅、最小振幅、弧长、平均瞬时相位、平均瞬时频率共5种类型储层地震属性。

2.根据权利要求1所述的方法，其特征在于，采用如下方式计算属性敏感度矩阵：

对所述的若干种储层敏感地震属性进行编号，具体为1,2,…,n，其中n为储层敏感地震属性总数，制作n行n列的矩阵，其中第i行第j列所对应的值为第i个储层敏感地震属性与第j个储层敏感地震属性组合的属性敏感度，获得属性敏感度矩阵，其中，1<＝i、j<＝n；

该属性敏感度矩阵中每行元素的敏感度平均值代表该行对应储层敏感地震属性与其他行对应的储层敏感地震属性的平均相关程度，所述平均相关程度越小说明该行所对应的储层敏感属性包含的冗余信息越少。

3.一种基于属性敏感度选择储层敏感地震属性组合的装置，其特征在于，包括：

第一处理装置，用于提取储层中潜在的多种类型储层地震属性；

第二处理装置，用于计算所述多种类型储层地震属性的地震属性敏感度S1；

第三处理装置，用于根据所述地震属性敏感度S1从大到小选出若干种储层地震属性作为储层敏感地震属性集；

第四处理装置，用于计算任意两种储层敏感地震属性之间的地震属性敏感度S2，并根据所述的地震属性敏感度S2得到地震属性敏感度矩阵；

第五处理装置，用于计算所述地震属性敏感度矩阵每行的敏感度平均值，作为该行所对应的储层敏感地震属性与其他储层敏感属性之间的平均相关程度，选择平均相关程度较小的多个储层敏感地震属性作为最终的储层敏感地震属性组合；

计算所述多种类型储层地震属性的地震属性敏感度S1包括：制作所述多种类型储层地震属性与待预测的储层参数的交会图，并根据所述交会图制作第一交会椭圆，测量所述第一交会椭圆的长轴长度R1与短轴长度R2，通过所述长轴长度R1除以短轴长度R2得到地震属性敏感度S1；

所述第一交会椭圆的制作方法具体为：

先绘制标准圆形以能够尽可能包含储层地震属性与待预测的储层参数之间的交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第一交会椭圆；

计算对应的任意两种储层敏感地震属性之间的地震属性敏感度S2包括：对所述储层敏感地震属性集中的任意两种储层敏感地震属性进行两两交会得到第二交会椭圆，并将第二交会椭圆的长轴长度r1除以短轴长度r2得到的值作为两种储层敏感地震属性之间的属性敏感度S2；

所述第二交会椭圆的制作方法包括：

先绘制标准圆形以能够尽可能包含两种储层敏感地震属性之间的交会点为原则，之后对标准圆形分别沿着长轴和短轴两个轴向方向进行压缩与拉伸将标准圆形改造为椭圆，以椭圆尽可能少包含非交会点区域为原则，得到第二交会椭圆；

提取储层中潜在的多种类型储层地震属性，具体为：

提取总负振幅、最小振幅、弧长、平均瞬时相位、平均瞬时频率共5种类型储层地震属性。

4.一种计算机可读存储介质，其特征在于，存储有计算机程序，所述计算机程序用于被处理器执行时实现如权利要求1或2所述的选择储层敏感地震属性组合的方法的步骤。

5.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1或2所述的选择储层敏感地震属性组合的方法的步骤。