CN105319582A - 一种地震属性参数的选择方法和装置 - Google Patents

Abstract

本公开是关于一种地震属性参数的选择方法和装置，属于勘探开发领域。所述方法包括：根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合；丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二集合中与所述储层参数相关度最高的多个地震属性参数。采用本公开，可以提高对储层参数特性判断的准确性。

Description

一种地震属性参数的选择方法和装置

技术领域

本公开涉及勘探开发领域，特别涉及一种地震属性参数的选择方法和装置。

背景技术

当前，随着经济的发展，对油气资源即石油和天然气的需求迅速增长。在开采效率不能明显提高的前提下，只有从大量的地震属性参数中选出最能体现储层参数特性的少量地震属性参数，才能提高确定钻井位置的准确性，最终提高油气资源的开采量。

在相关技术中，对于每个目标储层，都会采集数百个地震属性参数(如振幅类、层序统计类)，该数百个地震属性参数与目标储层的储层参数(如饱和度、孔隙度)具有不同的相关度，相关度越高，对储层参数的影响就越大。通过选出多个与目标储层的储层参数相关度最高的地震属性参数，得到每个目标储层的储层参数，进而确定钻井的位置。现有的地震属性参数的选择方法主要是专家选择法，即邀请具有丰富经验的专家，根据经验对大量的地震属性参数进行筛选，得到少量的、能体现储层结构的地震属性参数。

在实现本公开的过程中，发明人发现相关技术至少存在以下问题：

通过专家选择地震属性参数会带有专家的主观性，不能完全客观的选择地震属性参数，这样带有主观性选出的地震属性参数，会降低对储层参数特性判断的准确性。

发明内容

为了解决相关技术的问题，本公开实施例提供了一种地震属性参数的选择方法和装置，所述技术方案如下：

一方面，提供了一种地震属性参数的选择方法，所述方法包括：

根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，所述第一参数集合包含多个地震属性参数；

丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；

根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二集合中与所述储层参数相关度最高的多个地震属性参数；

根据所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置。

可选的，所述丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，包括：

根据专家选择法，从所述第一参数集合中选出与所述预设的目标储层的储层参数相关度最低的多个地震属性参数；

丢弃所述与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，将所述第一集合中剩余的与地震参数作为所述第二参数集合。

可选的，所述根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二参数集合中与所述储层参数相关度最高的多个地震属性参数，包括：

根据所述预设的目标储层的储层参数，结合所述第二参数集合中与所述储层参数相关度最高的多个地震属性参数，建立回归方程；

根据所述回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表所述预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

可选的，所述方法还包括：

根据与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置。

可选的，所述根据与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置，包括：

确定所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值；

根据所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从所述不同位置中确定所述钻井的位置。

可选的，所述钻井的位置，包括：

所述钻井位置的经纬度。

另一方面，提供了一种地震属性参数的选择装置，所述装置包括：

获取单元，用于根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，所述第一参数集合包含多个地震属性参数；

选择单元，用于丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；

处理单元，用于根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二集合中与所述储层参数相关度最高的多个地震属性参数；

所述处理单元，还用于根据所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置。

可选的，所述选择单元，用于：

根据专家选择法，从所述第一参数集合中选出与所述预设的目标储层的储层参数相关度最低的多个地震属性参数；

丢弃所述与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，将所述第一集合中剩余的与地震参数作为所述第二参数集合。

可选的，所述处理单元，用于：

根据所述预设的目标储层的储层参数，结合所述第二参数集合中与所述储层参数相关度最高的多个地震属性参数，建立回归方程；

根据所述回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表所述预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

可选的，所述处理单元还用于：

根据与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置。

可选的，所述处理单元，用于：

确定所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值；

根据所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从所述不同位置中确定所述钻井的位置。

可选的，所述钻井的位置，包括：

所述钻井位置的经纬度。

本公开的一些有益效果可以包括：

本公开实施例中，通过根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数，从而避免与储层参数相关度较低的地震属性参数对储层参数的影响，提高对储层参数特性判断的准确性。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性的，并不能限制本公开。

附图说明

此处所说明的附图用来提供对本公开的进一步理解，构成本申请的一部分，并不构成对本公开的限定。在附图中：

图1是根据一示例性实施例示出的一种地震属性参数的选择方法的流程图。

图2是根据一示例性实施例示出的确定地震属性参数的方法的示意图。

图3是根据一示例性实施例示出的一种地震属性参数的选择装置的结构示意图。

图4是根据一示例性实施例示出的一种装置的框图。

通过上述附图，已示出本公开明确的实施例，后文中将有更详细的描述。这些附图和文字描述并不是为了通过任何方式限制本公开构思的范围，而是通过参考特定实施例为本领域技术人员说明本公开的概念。

具体实施方式

为使本公开的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本公开做进一步详细说明。在此，本公开的示意性实施方式及其说明用于解释本公开，但并不作为对本公开的限定。

本公开实施例提供一种地震属性参数的选择方法和装置，以下结合附图对本公开进行详细说明。

实施例一

本公开实施例提供了一种地震属性参数的选择方法，如图1所示，该方法的处理流程可以包括如下步骤：

在步骤101中，根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，第一参数集合包含多个地震属性参数。

在步骤102中，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合。

在步骤103中，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数。

本公开实施例中，通过根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数，从而避免与储层参数相关度较低的地震属性参数对储层参数的影响，提高对储层参数特性判断的准确性。

实施例二

本公开实施例提供了一种地震属性参数的选择方法，下面将结合具体实施方式，对图1所示的处理流程进行详细的说明，内容可以如下：

在步骤101中，根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，第一参数集合包含多个地震属性参数。

其中，在预先存储的储层与参数集合的对应关系中，对应一个储层，有多个参数，每个参数包含多个地震参数。

在实施中，第一参数集合中包括孔隙度、饱和度等参数，本实施例中以饱和度为例，进行说明。

表1饱和度与地震属性参数的对应关系

表1中详细列出了多个储层中，每个储层中不同位置饱和度的数据，每一行代表一个储层，其中hd代表每个储层的单气层厚度，∑hd代表每个储层的累计气层厚度，sw代表每个储层的单气层饱和度，加权sw代表每个储层的加权平均饱和度。

加权平均饱和度的计算公式为： $\mathrm{sw}\_\mathrm{hd}=\frac{\mathrm{\Σ}{\mathrm{hd}}_i\×{\mathrm{sw}}_i}{\mathrm{\Σ}{\mathrm{hd}}_i}---\left(1\right)$

在公式(1)中，sw_hd代表的是厚度加权平均饱和度，hd_i代表的是单储层厚度，sw_i代表的是单储层饱和度，根据公式(1)进行运算，得到每个储层中的加权平均饱和度，并将该加权平均饱和度作为该储层的饱和度。

在步骤102中，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合。

在实施中，丢弃第一集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数可采用多种方法，由于步骤102的实现方式多种多样，本处实施例仅列举其中一种实现方式，具体方法流程如下：

步骤一，根据专家选择法，从第一参数集合中选出与预设的目标储层的储层参数相关度最低的多个地震属性参数。

其中，专家选择法是邀请具有丰富勘探领域知识的专家，通过自身掌握的理论知识，并结合自身的经验，从众多地震属性参数中选取与预设的目标储层的储层参数相关度最高的地震属性参数。

在实施中，应用专家选择法从第一参数集合中选出与预设的目标储层的储层参数相关度最低的多个地震属性参数，即选出第一参数集合中与目标储层的储层参数的最不相关的多个地震参数。

步骤二，丢弃与预设的目标储层的储层参数相关度最低的多个地震属性参数，将第一集合中剩余的与地震参数作为第二参数集合。

在实施中，将步骤一中选出的第一集合中与预设的目标储层的储层参数相关度最低的多个地震属性项丢弃。其中，第一集合中与饱和度相关的地震属性参数分为五大类，具体为：

(1)振幅统计类属性，包含均方根振幅，平均绝对振幅，最大波峰振幅，平均波峰振幅，最大波谷振幅，平均波谷振幅，最大绝对振幅，绝对振幅总量，振幅总量，平均能量，能量总体，平均振幅，振幅变化，振幅变化的不对称性，振幅的峰态等。

(2)复数道统计类属性，包含平均反射强度，平均瞬时频率，平均瞬时相位，反射强度斜率，瞬时频率斜率。它们是通过振幅包络，相位、频率及它们在纵向上的变化来反映波组特征的变化。

(3)谱统计类属性，包含有效带宽，弧线长度，平均零交叉点频率，三个频谱的主频率分量，谱峰值频率，由主频峰值到最大频率处的谱斜率。

(4)层序统计类属性，包含能量半衰时，能量半衰处斜率，正负采样点比，高于振幅门限的百分比，低于振幅门限的百分比，波峰数，波谷数。它反映了波级能量分布变化情况，能量分布是能量居前，能量居中，还是能量居后。

(5)相关统计类属性，包括：与下一个CDP的协方差，到下一个CDP相关时窗的滑动时间，平均信噪比，相关长度，相关分量，K-L信号复杂度等。

经过步骤一的专家选择法，从上述五大类的地震属性参数中选出与储层参数饱和度的相关度最低的多个地震属性参数，剩余的地震属性参数构成第二参数集合，如表2所示。

表2剩余的地震属性参数构成的第二参数集合

在步骤103中，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数。

在实施中，由于步骤103的实现方式多种多样，本实施例中仅列举其中一种实现方式，具体的方法流程如图2所示：

在步骤201中，根据预设的目标储层的储层参数，结合第二参数集合中与储层参数相关度最高的多个地震属性参数，建立回归方程。

其中，回归方程是根据第二参数集合中与储层参数相关度最高的多个地震属性参数，得到反映储层参数的数学关系表达式。

在实施中，根据表2中序号1的均方根振幅、平均瞬时频率、谱斜率、有效带宽、远近道振幅比、能量半衰时6个地震属性参数值，得到与序号1对应的饱和度的回归方程。

在步骤202中，根据回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

其中，逐步递归是在建立多元回归方程的过程中，按偏相关系数的大小次序将自变量逐个引入方程，对引入方程中的每个自变量偏相关系数进行统计检验，效应显著的自变量留在回归方程内，循此继续遴选下一个自变量。

在实施中，根据步骤201建立的回归方程，依次引入变量均方根振幅、平均瞬时频率、谱斜率、有效带宽、远近道振幅比、能量半衰时6个地震属性参数值，每当引入一个地震属性参数后，需要对引入的该地震属性参数的偏相关系数进行检验，如果检验效果不显著，停止引入新地震属性参数。由于新地震属性参数的引入，原已引入方程中的地震属性参数由于变量之间的相互作用其效应有可能变得不显著者，经统计检验确证后要随时从方程中剔除，只保留效应显著的地震属性参数。直至不再引入和剔除地震属性参数为止，从而得到最优的回归方程，格式例如S_w＝C₀+∑C_iA_i，其中，S_w代表预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

通过逐个引入地震属相参数的方法，根据表2中的实际地震属性参数值，最终得到的与饱和度对应的逐步回归方程为：

sw＝73.18-0.2879Inf+0.2127Slopf+0.0235EB-0.0205FN+0.0005Halft(2)

在公式(2)，sw代表预设的目标储层的储层参数，Inf代表平均瞬时频率，Slopf代表谱斜率，EB代表有效带宽，FN代表远近道振幅比，Half代表能量半衰时。各地震属性参数前的系数为各个地震属性参数，相对于饱和度这一储层参数的权重系数。

在确定第二集合中与储层参数相关度最高的多个地震属性参数后，还可以进行多种后续处理，下面以确定钻井位置为例进行说明：

根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值，从不同位置中确定钻井的位置。

其中，钻井的位置包括：钻井位置的经纬度。

在实施中，确定钻井位置的方法流程如下所示：

步骤一，确定与储层参数相关度最高的多个地震属性参数在不同位置的参数值。

步骤二，根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从不同位置中确定钻井的位置。

根据与每个井号6对应的平均瞬时频率、谱斜率、有效带宽、远近道振幅比、能量半衰时，确定该井号的饱和度的数值，按照该方法，依次确定表2中剩余井号的饱和度的数值，

在之前的步骤中根据公式2确定了与饱和度最相关的平均瞬时频率、谱斜率、有效带宽、远近道振幅比、能量半衰时5个地震属性参数，每个地震参数在不同的地理位置有不同的地震属性参数值，根据这5个地震属性参数的在不同的地理位置的地震属性参数值，确定能准确开采油气资源的钻井的地理位置，其中，该地理位置可以为钻井位置的经纬度。

本公开实施例中，通过根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数，从而避免与储层参数相关度较低的地震属性参数对储层参数的影响，提高对储层参数特性判断的准确性。

实施例三

基于相同的技术构思，本公开实施例还提供了一种地震属性参数的选择装置，如图3所示，装置包括：

获取单元310，用于根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，第一参数集合包含多个地震属性参数；

选择单元320，用于丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；

处理单元330，用于根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数。

可选的，选择单元320，用于：

根据专家选择法，从第一参数集合中选出与预设的目标储层的储层参数相关度最低的多个地震属性参数；

丢弃与预设的目标储层的储层参数相关度最低的多个地震属性参数，将第一集合中剩余的与地震参数作为第二参数集合。

可选的，处理单元330，用于：

根据预设的目标储层的储层参数，结合第二参数集合中与储层参数相关度最高的多个地震属性参数，建立回归方程；

根据回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

可选的，处理单元330还用于：

根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值，从不同位置中确定钻井的位置。

可选的，处理单元330，用于：

确定与储层参数相关度最高的多个地震属性参数在不同位置的参数值；

根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从不同位置中确定钻井的位置。

可选的，处理单元330确定的钻井的位置，包括：

所述钻井位置的经纬度。

本公开实施例中，通过根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数，从而避免与储层参数相关度较低的地震属性参数对储层参数的影响，提高对储层参数特性判断的准确性。

实施例四

图4是根据一示例性实施例示出的一种用于地震属性参数的选择的终端800的框图。例如，终端800可以是移动电话，计算机，数字广播终端，消息收发设备，游戏控制台，平板设备，医疗设备，健身设备，个人数字助理等。

参照图4，终端800可以包括以下一个或多个组件：处理组件802，存储器804，电源组件806，多媒体组件808，音频组件810，输入/输出(I/O)的接口812，传感器组件814，以及通信组件816。

处理组件802通常控制终端800的整体操作，诸如与显示，电话呼叫，数据通信，相机操作和记录操作相关联的操作。处理元件802可以包括一个或多个处理器820来执行指令，以完成上述的方法的全部或部分步骤。此外，处理组件802可以包括一个或多个模块，便于处理组件802和其他组件之间的交互。例如，处理部件802可以包括多媒体模块，以方便多媒体组件808和处理组件802之间的交互。

存储器804被配置为存储各种类型的数据以支持在终端800的操作。这些数据的示例包括用于在终端800上操作的任何应用程序或方法的指令，联系人数据，电话簿数据，消息，图片，视频等。存储器804可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。

电力组件806为终端800的各种组件提供电力。电力组件806可以包括电源管理系统，一个或多个电源，及其他与为终端800生成、管理和分配电力相关联的组件。

多媒体组件808包括在所述终端800和用户之间的提供一个输出接口的屏幕。在一些实施例中，屏幕可以包括液晶显示器(LCD)和触摸面板(TP)。如果屏幕包括触摸面板，屏幕可以被实现为触摸屏，以接收来自用户的输入信号。触摸面板包括一个或多个触摸传感器以感测触摸、滑动和触摸面板上的手势。所述触摸传感器可以不仅感测触摸或滑动动作的边界，而且还检测与所述触摸或滑动操作相关的持续时间和压力。在一些实施例中，多媒体组件808包括一个前置摄像头和/或后置摄像头。当终端800处于操作模式，如拍摄模式或视频模式时，前置摄像头和/或后置摄像头可以接收外部的多媒体数据。每个前置摄像头和后置摄像头可以是一个固定的光学透镜系统或具有焦距和光学变焦能力。

音频组件810被配置为输出和/或输入音频信号。例如，音频组件810包括一个麦克风(MIC)，当终端800处于操作模式，如呼叫模式、记录模式和语音识别模式时，麦克风被配置为接收外部音频信号。所接收的音频信号可以被进一步存储在存储器804或经由通信组件816发送。在一些实施例中，音频组件810还包括一个扬声器，用于输出音频信号。

I/O接口812为处理组件802和外围接口模块之间提供接口，上述外围接口模块可以是键盘，点击轮，按钮等。这些按钮可包括但不限于：主页按钮、音量按钮、启动按钮和锁定按钮。

传感器组件814包括一个或多个传感器，用于为终端800提供各个方面的状态评估。例如，传感器组件814可以检测到终端800的打开/关闭状态，组件的相对定位，例如所述组件为终端800的显示器和小键盘，传感器组件814还可以检测终端800或终端800一个组件的位置改变，用户与终端800接触的存在或不存在，终端800方位或加速/减速和终端800的温度变化。传感器组件814可以包括接近传感器，被配置用来在没有任何的物理接触时检测附近物体的存在。传感器组件814还可以包括光传感器，如CMOS或CCD图像传感器，用于在成像应用中使用。在一些实施例中，该传感器组件814还可以包括加速度传感器，陀螺仪传感器，磁传感器，压力传感器或温度传感器。

通信组件816被配置为便于终端800和其他设备之间有线或无线方式的通信。终端800可以接入基于通信标准的无线网络，如WiFi，2G或3G，或它们的组合。在一个示例性实施例中，通信部件816经由广播信道接收来自外部广播管理系统的广播信号或广播相关信息。在一个示例性实施例中，所述通信部件816还包括近场通信(NFC)模块，以促进短程通信。例如，在NFC模块可基于射频识别(RFID)技术，红外数据协会(IrDA)技术，超宽带(UWB)技术，蓝牙(BT)技术和其他技术来实现。

在示例性实施例中，终端800可以被一个或多个应用专用集成电路(ASIC)、数字信号处理器(DSP)、数字信号处理设备(DSPD)、可编程逻辑器件(PLD)、现场可编程门阵列(FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述方法。

在示例性实施例中，还提供了一种包括指令的非临时性计算机可读存储介质，例如包括指令的存储器804，上述指令可由终端800的处理器820执行以完成上述方法。例如，所述非临时性计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。

一种非临时性计算机可读存储介质，当所述存储介质中的指令由移动终端的处理器执行时，使得移动终端能够执行一种地震属性参数的选择方法，所述方法包括：

根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，第一参数集合包含多个地震属性参数；

丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；

根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数。

可选的，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，包括：

根据专家选择法，从第一参数集合中选出与预设的目标储层的储层参数相关度最低的多个地震属性参数；

丢弃与预设的目标储层的储层参数相关度最低的多个地震属性参数，将第一集合中剩余的与地震参数作为第二参数集合。

可选的，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二参数集合中与储层参数相关度最高的多个地震属性参数，包括：

根据预设的目标储层的储层参数，结合第二参数集合中与储层参数相关度最高的多个地震属性参数，建立回归方程；

根据回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

可选的，该方法还包括：

根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值，从不同位置中确定钻井的位置。

可选的，根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值，从不同位置中确定钻井的位置包括：

确定与储层参数相关度最高的多个地震属性参数在不同位置的参数值；

根据与储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从不同位置中确定钻井的位置。

可选的，钻井的位置，包括：

钻井位置的经纬度。

本公开实施例中，通过根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，丢弃第一参数集合中与预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，根据预设的目标储层的储层参数，通过逐步回归算法，确定第二集合中与储层参数相关度最高的多个地震属性参数，从而避免与储层参数相关度较低的地震属性参数对储层参数的影响，提高对储层参数特性判断的准确性。

本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由下面的权利要求指出。

应当理解的是，本发明并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。

Claims (12)

1.一种地震属性参数的选择方法，其特征在于，所述方法包括：

根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，所述第一参数集合包含多个地震属性参数；

丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；

根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二集合中与所述储层参数相关度最高的多个地震属性参数。

2.根据权利要求1所述的方法，其特征在于，所述丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合，包括：

根据专家选择法，从所述第一参数集合中选出与所述预设的目标储层的储层参数相关度最低的多个地震属性参数；

丢弃所述与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，将所述第一集合中剩余的与地震参数作为所述第二参数集合。

3.根据权利要求1所述的方法，其特征在于，所述根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二参数集合中与所述储层参数相关度最高的多个地震属性参数，包括：

根据所述预设的目标储层的储层参数，结合所述第二参数集合中与所述储层参数相关度最高的多个地震属性参数，建立回归方程；

根据所述回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表所述预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

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

根据与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置。

5.根据权利要求4所述的方法，其特征在于，所述根据与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置包括：

确定所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值；

根据所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从所述不同位置中确定所述钻井的位置。

6.根据权利要求5所述的方法，其特征在于，所述钻井的位置，包括：

所述钻井位置的经纬度。

7.一种地震属性参数的选择装置，其特征在于，所述装置包括：

获取单元，用于根据预先存储的储层与参数集合的对应关系，获取目标储层对应的第一参数集合，所述第一参数集合包含多个地震属性参数；

选择单元，用于丢弃所述第一参数集合中与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，得到第二参数集合；

处理单元，用于根据所述预设的目标储层的储层参数，通过逐步回归算法，确定所述第二集合中与所述储层参数相关度最高的多个地震属性参数。

8.根据权利要求7所述的装置，其特征在于，所述选择单元，用于：

根据专家选择法，从所述第一参数集合中选出与所述预设的目标储层的储层参数相关度最低的多个地震属性参数；

丢弃所述与所述预设的目标储层的储层参数相关度最低的多个地震属性参数，将所述第一集合中剩余的与地震参数作为所述第二参数集合。

9.根据权利要求7所述的装置，其特征在于，所述处理单元，用于：

根据所述预设的目标储层的储层参数，结合所述第二参数集合中与所述储层参数相关度最高的多个地震属性参数，建立回归方程；

根据所述回归方程，通过逐步回归的算法，得到如下回归方程S_w＝C₀+∑C_iA_i，其中，S_w代表所述预设的目标储层的储层参数，C₀代表回归方程中的初始值，C_i代表回归方程中的回归系数，A_i代表回归方程中优选出的地震属性参数，i为不为零的自然数。

10.根据权利要求7所述的装置，其特征在于，所述处理单元还用于：

根据与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值，从所述不同位置中确定钻井的位置。

11.根据权利要求10所述的装置，其特征在于，所述处理单元，用于：

确定所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值；

根据所述与所述储层参数相关度最高的多个地震属性参数在不同位置的参数值对应的地理位置参数，从所述不同位置中确定所述钻井的位置。

12.根据权利要求11所述的装置，其特征在于，所述钻井的位置，包括：

所述钻井位置的经纬度。