CN115705451A - 岩性敏感弹性参数的反演方法、装置及计算机存储介质 - Google Patents

Abstract

本申请实施例公开了一种岩性敏感弹性参数的反演方法、装置及计算机存储介质，属于地震勘探技术领域。该方法包括：获取目标区储层的地震相编码、地层格架编码和测井资料；根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，所述弹性阻抗模型用于描述所述目标区储层的地震反射特征空间变化规律；根据所述弹性阻抗模型进行反演处理，得到所述目标区的岩性敏感弹性参数，所述岩性敏感弹性参数用于模拟所述目标区储层的岩性三维空间分布特征。本申请实施例通过将地震反射特征的空间变化规律与测井资料的纵向特征结合开展建模和约束反演，提高了确定岩性敏感弹性参数的准确性，提高了岩性储层预测的精度。

Description

岩性敏感弹性参数的反演方法、装置及计算机存储介质

技术领域

本申请实施例涉及地震勘探技术领域，特别涉及一种岩性敏感弹性参数的反演方法、装置及计算机存储介质。

背景技术

随着地震勘探技术的发展，通常能够通过地震资料对储层的岩性进行预测。由于叠前地震数据中通常包含了地震波振幅随偏移距离变化的丰富信息，因此，能够通过叠前地震数据进行反演来获得地层纵波速度、横波速度、密度、纵横波速度比、泊松比等对岩性特征比较敏感的弹性参数。

目前，在通过叠前地震数据对岩性敏感弹性参数进行反演时，能够建立地震资料和测井资料之间的映射关系，并利用叠前地震资料反演出波阻抗、纵横波速度比、泊松比、泊松阻抗等岩性敏感弹性参数，实现对岩性储层空间分布特征的表征。

但是，常规叠前地震反演中采用测井资料内插外推建立的反演约束模型，在空间上的横向变化关系难以反映地震反射模式在空间上的差异性，使得弹性参数反演结果难以精确描述非均质性较强的复杂岩性储层空间分布特征。

发明内容

本申请实施例提供了一种岩性敏感弹性参数的反演方法、装置及计算机存储介质，可以用于解决相关技术中弹性参数反演结果不准确，导致不能准确描述复杂岩性储层空间分布特征的问题。所述技术方案如下：

一方面，提供了一种岩性敏感弹性参数的反演方法，所述方法包括：

获取目标区储层的地震相编码、地层格架编码和测井资料；

根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，所述弹性阻抗模型用于描述所述目标区储层的地震反射特征空间变化规律；

根据所述弹性阻抗模型进行反演处理，得到所述目标区的岩性敏感弹性参数，所述岩性敏感弹性参数用于模拟所述目标区储层的岩性三维空间分布特征。

在一些实施例中，获取目标区储层的地震相编码和地层格架编码，包括：

获取所述目标区储层的三维地震数据和地质信息；

根据所述三维地震数据，构建所述地震相编码；

根据所述三维地震数据和所述地质信息，构建所述地层格架编码。

在一些实施例中，所述根据所述三维地震数据，构建所述地震相编码，包括：

从所述目标区储层的三维地震数据中提取地震属性数据体，所述地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；

通过聚类算法对所述地震属性数据体进行聚类处理，得到所述目标区储层的地震相类型以及每种地震相的分布范围；

根据所述地震相类型和所述每种地震相的分布范围，对每种地震相类型进行编码，得到所述地震相编码。

在一些实施例中，所述根据所述三维地震数据和所述地质信息，构建所述地层格架编码，包括：

根据所述三维地震数据和所述地质信息，获取所述目标区储层的地震层位、断层分布特征和空间组合关系；

根据所述目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；

确定所述地震层位与所述断层趋势面之间的接触关系；

根据所述地震层位与所述断层趋势面之间的接触关系，对所述目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；

对相邻两个层位之间的离散化网格用数字进行编码，得到所述目标区储层的地层格架编码。

在一些实施例中，所述根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，包括：

从所述测井资料中获取所述目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，所述第一角度弹性阻抗曲线的第一角度范围中最小角度大于所述第二角度弹性阻抗曲线的第二角度范围中最大角度，所述第二角度弹性阻抗曲线的第三角度范围中最小角度大于所述第三角度弹性阻抗曲线的角度范围中最大角度；

从所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；

根据所述每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；

通过所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线进行空间插值外推，得到所述目标位置处的弹性阻抗数值，所述目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；

根据所述同一地震相类型对应的参考弹性阻抗数值和所述目标位置的弹性阻抗数值，确定所述目标区储层的弹性阻抗模型。

在一些实施例中，所述根据所述弹性阻抗模型进行反演处理，得到所述目标区储层的岩性敏感弹性参数，包括：

通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据；

根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量；

根据所述目标区储层的体积模量和剪切模量，确定所述目标区储层的岩性敏感弹性参数。

在一些实施例中，所述通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据，包括：

根据所述测井资料，将所述目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；

根据所述角度部分叠加地震数据和所述测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；

根据所述角度部分叠加地震数据、所述三个角度的地震子波和所述弹性阻抗模型，建立弹性阻抗反演目标函数；

对所述弹性阻抗反演目标函数通过迭代算法进行求解处理，得到所述目标区储层的弹性阻抗反演数据。

在一些实施例中，所述根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量，包括：

根据所述测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出所述测井资料中的弹性阻抗曲线与所述目标测井的体积模量和剪切模量之间的非线性模量转换系数，所述目标测井为所述目标区储层范围内的任一测井；

根据所述非线性模量转换系数和所述目标区储层的弹性阻抗反演数据，确定所述目标区储层的体积模量和剪切模量。

在一些实施例中，所述根据所述目标区储层的体积模量和剪切模量数据，确定所述目标区储层的岩性敏感弹性参数，包括：

根据所述测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定所述目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，所述数据点为所述目标区存储范围内具有测井资料的位置，所述参考体积模量为对所述所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，所述参考剪切模量为对所述所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；

根据所述目标区储层的体积模量和剪切模量，以及所述参考体积模量和所述参考剪切模量，确定所述目标区储层的岩性敏感弹性参数。

另一方面，提供了一种岩性敏感弹性参数的反演装置，所述装置包括：

获取模块，用于获取目标区储层的地震相编码、地层格架编码和测井资料；

确定模块，用于根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，所述弹性阻抗模型用于描述所述目标区储层的地震反射特征空间变化规律；

反演模块，用于根据所述弹性阻抗模型进行反演处理，得到所述目标区的岩性敏感弹性参数，所述岩性敏感弹性参数用于模拟所述目标区储层的岩性三维空间分布特征。

在一些实施例中，所述获取模块包括：

第一获取子模块，用于获取所述目标区储层的三维地震数据和地质信息；

第一构建子模块，用于根据所述三维地震数据，构建所述地震相编码；

第二构建子模块，用于根据所述三维地震数据和所述地质信息，构建所述地层格架编码。

在一些实施例中，所述第一构建子模块用于：

从所述目标区储层的三维地震数据中提取地震属性数据体，所述地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；

通过聚类算法对所述地震属性数据体进行聚类处理，得到所述目标区储层的地震相类型以及每种地震相的分布范围；

根据所述地震相类型和所述每种地震相的分布范围，对每种地震相类型进行编码，得到所述地震相编码。

在一些实施例中，所述第二构建子模块用于：

根据所述三维地震数据和所述地质信息，获取所述目标区储层的地震层位、断层分布特征和空间组合关系；

根据所述目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；

确定所述地震层位与所述断层趋势面之间的接触关系；

根据所述地震层位与所述断层趋势面之间的接触关系，对所述目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；

对相邻两个层位之间的离散化网格用数字进行编码，得到所述目标区储层的地层格架编码。

在一些实施例中，所述确定模块包括：

第二获取子模块，用于从所述测井资料中获取所述目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，所述第一角度弹性阻抗曲线的第一角度范围中最小角度大于所述第二角度弹性阻抗曲线的第二角度范围中最大角度，所述第二角度弹性阻抗曲线的第三角度范围中最小角度大于所述第三角度弹性阻抗曲线的角度范围中最大角度；

第三获取子模块，用于从所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；

第一确定子模块，用于根据所述每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；

插值子模块，用于通过所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线进行空间插值外推，得到所述目标位置处的弹性阻抗数值，所述目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；

第二确定子模块，用于根据所述同一地震相类型对应的参考弹性阻抗数值和所述目标位置的弹性阻抗数值，确定所述目标区储层的弹性阻抗模型。

在一些实施例中，所述反演模块包括：

反演子模块，用于通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据；

第三确定子模块，用于根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量；

第四确定子模块，用于根据所述目标区储层的体积模量和剪切模量，确定所述目标区储层的岩性敏感弹性参数。

在一些实施例中，所述反演子模块用于：

根据所述测井资料，将所述目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；

根据所述角度部分叠加地震数据和所述测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；

根据所述角度部分叠加地震数据、所述三个角度的地震子波和所述弹性阻抗模型，建立弹性阻抗反演目标函数；

对所述弹性阻抗反演目标函数通过迭代算法进行求解处理，得到所述目标区储层的弹性阻抗反演数据。

在一些实施例中，所述第三确定子模块用于：

根据所述测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出所述测井资料中的弹性阻抗曲线与所述目标测井的体积模量和剪切模量之间的非线性模量转换系数，所述目标测井为所述目标区储层范围内的任一测井；

根据所述非线性模量转换系数和所述目标区储层的弹性阻抗反演数据，确定所述目标区储层的体积模量和剪切模量。

在一些实施例中，所述第四确定子模块用于：

根据所述测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定所述目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，所述数据点为所述目标区存储范围内具有测井资料的位置，所述参考体积模量为对所述所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，所述参考剪切模量为对所述所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；

根据所述目标区储层的体积模量和剪切模量，以及所述参考体积模量和所述参考剪切模量，确定所述目标区储层的岩性敏感弹性参数。

另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，所述指令被处理器执行时实现上述岩性敏感弹性参数的反演方法中的任一步骤。

本申请实施例提供的技术方案带来的有益效果至少包括：

在本申请实施例中，能够通过地震相编码、地层格架编码和测井资料进行弹性阻抗模型的建模，并根据弹性阻抗模型进行反演，得到岩性敏感弹性参数，即通过将地震反射特征的空间变化规律与测井资料的纵向特征结合开展建模和约束反演，提高了确定岩性敏感弹性参数的准确性，提高了岩性储层预测的精度。

附图说明

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

图1是本申请实施例提供的一种岩性敏感弹性参数的反演方法流程图；

图2是本申请实施例提供的一种岩性敏感弹性参数的反演方法流程图；

图3是本申请实施例提供的一种地震属性切片图的示意图；

图4是本申请实施例提供的一种地震相类型的示意图；

图5是本申请实施例提供的一种岩性敏感弹性参数的反演装置的结构示意图；

图6是本申请实施例提供的一种获取模块的结构示意图；

图7是本申请实施例提供的一种确定模块的结构示意图；

图8是本申请实施例提供的一种反演模块的结构示意图；

图9是本申请实施例提供的一种终端的结构示意图。

具体实施方式

为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。

在对本申请实施例提供的一种岩性敏感弹性参数的反演方法进行详细地解释说明之前，先对本申请实施例提供的一种应用场景进行解释说明。

地震波在实际地球介质中传播时携带了大量与储层性质密切相关的特征信息，采用地震资料开展储层岩性预测已经在油气勘探中发挥着非常重要的作用。由于叠前地震数据中包含了地震波振幅随偏移距变化的丰富信息，通过叠前地震反演能够获得地层纵波速度、横波速度、密度、纵横波速度比、泊松比等对岩性特征比较敏感的弹性参数，因此，为了进行岩性储层预测，通常需要利用叠前地震资料反演岩性敏感弹性参数。

基于这样的应用场景，本申请实施例提供了一种提高确定岩性敏感弹性参数准确性的岩性敏感弹性参数的反演方法。

图1是本申请实施例提供的一种岩性敏感弹性参数的反演方法流程图，该岩性敏感弹性参数的反演方法可以包括如下几个步骤：

步骤101：获取目标区储层的地震相编码、地层格架编码和测井资料。

步骤102：根据该目标区储层的地震相编码、地层格架编码和测井资料，确定该目标区储层的弹性阻抗模型，该弹性阻抗模型用于描述该目标区储层的地震反射特征空间变化规律。

步骤103：根据该弹性阻抗模型进行反演处理，得到该目标区的岩性敏感弹性参数，该岩性敏感弹性参数用于模拟该目标区储层的岩性三维空间分布特征。

在本申请实施例中，能够通过地震相编码、地层格架编码和测井资料进行弹性阻抗模型的建模，并根据弹性阻抗模型进行反演，得到岩性敏感弹性参数，即通过将地震反射特征的空间变化规律与测井资料的纵向特征结合开展建模和约束反演，提高了确定岩性敏感弹性参数的准确性，提高了岩性储层预测的精度。

在一些实施例中，获取目标区储层的地震相编码和地层格架编码，包括：

获取该目标区储层的三维地震数据和地质信息；

根据该三维地震数据，构建该地震相编码；

根据该三维地震数据和该地质信息，构建该地层格架编码。

在一些实施例中，根据该三维地震数据，构建该地震相编码，包括：

从该目标区储层的三维地震数据中提取地震属性数据体，该地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；

通过聚类算法对该地震属性数据体进行聚类处理，得到该目标区储层的地震相类型以及每种地震相的分布范围；

根据该地震相类型和该每种地震相的分布范围，对每种地震相类型进行编码，得到该地震相编码。

在一些实施例中，根据该三维地震数据和该地质信息，构建该地层格架编码，包括：

根据该三维地震数据和该地质信息，获取该目标区储层的地震层位、断层分布特征和空间组合关系；

根据该目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；

确定该地震层位与该断层趋势面之间的接触关系；

根据该地震层位与该断层趋势面之间的接触关系，对该目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；

对相邻两个层位之间的离散化网格用数字进行编码，得到该目标区储层的地层格架编码。

在一些实施例中，根据该目标区储层的地震相编码、地层格架编码和测井资料，确定该目标区储层的弹性阻抗模型，包括：

从该测井资料中获取该目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，该第一角度弹性阻抗曲线的第一角度范围中最小角度大于该第二角度弹性阻抗曲线的第二角度范围中最大角度，该第二角度弹性阻抗曲线的第三角度范围中最小角度大于该第三角度弹性阻抗曲线的角度范围中最大角度；

从该第一角度弹性阻抗曲线、该第二角度弹性阻抗曲线和该第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；

根据该每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；

通过该第一角度弹性阻抗曲线、该第二角度弹性阻抗曲线和该第三角度弹性阻抗曲线进行空间插值外推，得到该目标位置处的弹性阻抗数值，该目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；

根据该同一地震相类型对应的参考弹性阻抗数值和该目标位置的弹性阻抗数值，确定该目标区储层的弹性阻抗模型。

在一些实施例中，根据该弹性阻抗模型进行反演处理，得到该目标区储层的岩性敏感弹性参数，包括：

通过该弹性阻抗模型进行弹性阻抗反演，得到该目标区储层的弹性阻抗反演数据；

根据该弹性阻抗反演数据和该目标区储层的测井资料，确定该目标区储层的体积模量和剪切模量；

根据该目标区储层的体积模量和剪切模量，确定该目标区储层的岩性敏感弹性参数。

在一些实施例中，通过该弹性阻抗模型进行弹性阻抗反演，得到该目标区储层的弹性阻抗反演数据，包括：

根据该测井资料，将该目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；

根据该角度部分叠加地震数据和该测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；

根据该角度部分叠加地震数据、该三个角度的地震子波和该弹性阻抗模型，建立弹性阻抗反演目标函数；

对该弹性阻抗反演目标函数通过迭代算法进行求解处理，得到该目标区储层的弹性阻抗反演数据。

在一些实施例中，根据该弹性阻抗反演数据和该目标区储层的测井资料，确定该目标区储层的体积模量和剪切模量，包括：

根据该测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出该测井资料中的弹性阻抗曲线与该目标测井的体积模量和剪切模量之间的非线性模量转换系数，该目标测井为该目标区储层范围内的任一测井；

根据该非线性模量转换系数和该目标区储层的弹性阻抗反演数据，确定该目标区储层的体积模量和剪切模量。

在一些实施例中，根据该目标区储层的体积模量和剪切模量数据，确定该目标区储层的岩性敏感弹性参数，包括：

根据该测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定该目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，该数据点为该目标区存储范围内具有测井资料的位置，该参考体积模量为对该所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，该参考剪切模量为对该所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；

根据该目标区储层的体积模量和剪切模量，以及该参考体积模量和该参考剪切模量，确定该目标区储层的岩性敏感弹性参数。

上述所有可选技术方案，均可按照任意结合形成本申请的可选实施例，本申请实施例对此不再一一赘述。

图2是本申请实施例提供的一种岩性敏感弹性参数的反演方法流程图，本实施例以该岩性敏感弹性参数的反演方法应用于终端中进行举例说明，该岩性敏感弹性参数的反演方法可以包括如下几个步骤：

步骤201：终端获取目标区储层的地震相编码、地层格架编码和测井资料。

由于地震相能够反应地层的地震特征，进而反映出地层特征，且不同地震相能够反应出不同的地震特征及地层特征，地层格架是指广泛应用于地层序列中各类地层或岩石单位的区域性时空有序排列型式，即地层格架能够反应地层的时空分布特征，因此，为了准确确定目标区储层的岩性敏感弹性参数，终端能够获取目标区储层的地震相编码、地层格架编码和测井资料。该测井资料为针对目标区储层的所有测井的测井资料或部分测井的测井资料。

作为一种示例，终端能够在接收到获取指令时，从存储的文件中获取目标区储层的地震相编码、地层格架编码和测井资料。也即是，终端能够事先获取并存储目标区储层的地震相编码、地层格架编码和测井资料，并在接收到获取指令时，从存储的文件中获取目标区储层的地震相编码、地层格架编码和测井资料。或者，终端还能够在接收到获取指令时，向其他设备发送获取请求，该获取请求中携带目标区储层的标识；其他设备在接收到终端获取请求时，根据目标区储层的标识向终端返回目标区储层的地震相编码、地层格架编码和测井资料，从而使终端获取到目标区储层的地震相编码、地层格架编码和测井资料。

需要说明的是，获取指令为用户通过指定操作作用在终端上时触发，该指定操作包括语音操作、点击操作、滑动操作等等。

作为一种示例，终端不仅能够在接收到获取指令时，按照上述方式获取目标区储层的地震相编码、地层格架编码和测井资料，还能够通过其他方式获取目标区储层的地震相编码和地层格架编码。

作为一种示例，终端获取目标区储层的地震相编码和地层格架编码的操作至少包括：获取目标区储层的三维地震数据和地质信息；根据三维地震数据，构建地震相编码；根据三维地震数据和地质信息，构建地层格架编码。

由于目标区储层的地震相编码和地层格架编码与目标区储层的三维地震数据和所处地质的地质信息相关，因此，终端能够获取目标区储层的三维地震数据和地质信息。

在一些实施例中，终端能够从存储的文件中获取三维地震数据和地质信息。

在一些实施例中，终端根据三维地震数据，构建地震相编码的操作包括：从目标区储层的三维地震数据中提取地震属性数据体，该地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；通过聚类算法对地震属性数据体进行聚类处理，得到目标区储层的地震相类型以及每种地震相的分布范围；根据地震相类型和每种地震相的分布范围，对每种地震相类型进行编码，得到地震相编码。

由于地震属性数据体能够反应底下岩性特征，比如，振幅类属性能够用于检测地层变薄带，频率信息能够用于检测断层等等，因此，终端能够从目标区储层的三维地震数据中提取地震属性数据体，该地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据。

由于目标区储层中可能会因地质分布范围不同，导致存在多种类型的地震相，因此，终端能够通过聚类算法对地震属性数据体进行聚类处理，得到目标区储层的地震相类型以及每种地震相的分布范围。

需要说明的是，终端从目标区储层的三维地震数据中提取地震属性数据体，以及通过聚类算法对地震属性数据体进行聚类处理，得到目标区储层的地震相类型以及每种地震相的分布范围的操作能够参考相关技术，本申请实施例对此不再进行一一赘述。

作为一种示例，终端对每种地震相类型进行编码，得到地震相编码的操作包括：从任一整数开始，按照从小到大顺序，为每种地震相类型设置对应的数字标识，得到每种地震相类型对应的地震相编码，任意两个类型不相同的地震相对应不同的数字标识。

在一种实施环境中，终端从目标区储层的三维地震数据中提取地震属性数据体，该三种地震属性数据体对应的地震属性切片图如图3所示。终端通过聚类算法对目标区储层的地震属性数据体进行聚类处理后，能够将目标区储层的地震相类型划分为N类，N的取值与目标区储层的地质特征有关。示例性地，终端能够得到如图4所示的地震相。在确定目标区储层的地震相类型后，终端能够从1开始为每种地震相类型设置对应的数字标识，得到地震相编码：1，2，……，N。

作为一种示例，终端根据三维地震数据和地质信息，构建地层格架编码的操作包括：根据三维地震数据和地质信息，获取目标区储层的地震层位、断层分布特征和空间组合关系；根据目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；确定地震层位与断层趋势面之间的接触关系；根据地震层位与断层趋势面之间的接触关系，对目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；对相邻两个层位之间的离散化网格用数字进行编码，得到目标区储层的地层格架编码。

由于三维地震数据能够反应地层的相关信息，因此，终端能够根据三维地震数据对目标区储层进行地震层位和断层解释。对地震层位与断层的解释包括根据三维地震数据和地质信息，获取目标区储层的地震层位、断层分布特征和空间组合关系；根据目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；确定地震层位与断层趋势面之间的接触关系。

在一些实施例中，断层分布特征中包括多个断层特征点，终端在根据目标区储层的断层分布特征和空间组合关系，模拟断层趋势面之前，能够将多个断层特征点中突兀的断层特征点从三维空间中剔除，该突兀的断层特征点是指与相邻的参考特征点之间的距离大于距离阈值的断层特征点，该参考特征点为多个断层特征点中除突兀的断层特征点之外的其他断层特征点。

在一些实施例中，终端能够将多个断层特征点中相邻断层特征点连接后建立平滑的断层趋势面。之后将地震层位与断层趋势面相结合，以建立地层格架，并由深层到浅层逐一确定每个地震层位与断层趋势面之间的接触点，每个地震层位与断层趋势面之间的接触关系一致，以确保目标区储层的地层格架中每个磁层之间的关系合理。

在一些实施例中，终端对相邻两个层位之间的离散化网格用数字进行编码，得到目标区储层的地层格架编码时，按照从深层到浅层的顺序，对两个相邻层位之间的离散化网格通过数字进行编码，且在通过数字编码时，按照从小到大的顺序依次进行编码，以得到目标区储层的地层格架编码。比如，终端能够从数字1开始按照从小到大的顺序依次进行编码，得到目标区储层的地层格架编码。

步骤202：终端根据目标区储层的地震相编码、地层格架编码和测井资料，确定目标区储层的弹性阻抗模型，该弹性阻抗模型用于描述目标区储层的地震反射特征空间变化规律。

作为一种示例，终端根据目标区储层的地震相编码、地层格架编码和测井资料，确定目标区储层的弹性阻抗模型的操作包括：从测井资料中获取目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，该第一角度弹性阻抗曲线的角度范围中最小角度大于第二角度弹性阻抗曲线的角度范围中最大角度，第二角度弹性阻抗曲线的角度范围中最小角度大于第三角度弹性阻抗曲线的角度范围中最大角度；从第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；根据每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；通过第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线进行空间插值外推，得到目标位置处的弹性阻抗数值，该目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；根据同一地震相类型对应的参考弹性阻抗数值和目标位置的弹性阻抗数值，确定目标区储层的弹性阻抗模型。

由于测井资料中包括所有测井相关的信息，因此，终端能够从测井资料中获取目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线。终端能够将角度范围位于第一角度范围的弹性阻抗曲线确定为第一角度弹性阻抗曲线，将角度范围位于第二角度范围的弹性阻抗曲线确定为第二角度弹性阻抗曲线，将角度范围位于第三角度范围的弹性阻抗曲线确定为第三角度弹性阻抗曲线，且第一角度范围的最小角度大于第二角度范围的最大角度，第二角度范围的最小角度大于第三角度范围的最大角度。

在一些实施例中，终端根据每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值的操作包括：当目标地震相类型对应的相同地震相编码的数量只有一个时，将目标地震相类型的地震相编码对应的弹性阻抗数值确定为目标地震相类型的参考弹性阻抗数值，目标地震相类型为目标区储层中包括的多个地震相类型中的任一地震相类型；当目标地震相类型对应的相同地震相编码的数量有多个时，将目标地震相类型对应的多个弹性阻抗数值相加后除以目标地震相类型对应的相同地震相编码的数量，并将计算得到的均值确定为目标地震相类型的参考弹性阻抗数值。

由于目标区储层内对于同一个地震相类型，该地震相可能并不是分布在同一个位置，因此，同一个地震相类型可能对应有多个相同的地震相编码，比如，填充相类型的地震相在目标区储层范围内的东边和西边同时存在，且没有连接在一起，因此，类型为填充相的地震相对应有2个相同的地震相编码。

需要说明的是，终端通过第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线进行空间插值外推时，能够按照地震相编码从小到大的顺序，依次在具有同一个地震相编码的离散化网格内进行空间插值外推，得到目标位置处的弹性阻抗数值。

需要说明的是，终端通过第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线进行空间插值外推的操作能够参考相关技术，本申请实施例对此不再进行一一赘述。

在一些实施例中，终端根据同一地震相类型对应的参考弹性阻抗数值和目标位置的弹性阻抗数值，通过下述第一公式确定目标区储层的弹性阻抗模型。

EI_model＝λ₁EI_well+λ₂EI_seismic (1)

需要说明的是，在上述第一公式(1)中，EI_model为地层网格内任一点位置处估计的弹性阻抗模型的值，EI_well为目标位置的弹性阻抗数值，EI_seismic为参考弹性阻抗数值，λ₁为目标位置的弹性阻抗数值在目标位置处对应的权重系数，λ₂为参考弹性阻抗数值在对应位置处的权重系数。

步骤203：终端根据弹性阻抗模型进行反演处理，得到目标区的岩性敏感弹性参数，该岩性敏感弹性参数用于模拟目标区储层的岩性三维空间分布特征。

作为一种示例，终端根据弹性阻抗模型进行反演处理，得到目标区储层的岩性敏感弹性参数的操作包括：通过弹性阻抗模型进行弹性阻抗反演，得到目标区储层的弹性阻抗反演数据；根据弹性阻抗反演数据和目标区储层的测井资料，确定目标区储层的体积模量和剪切模量；根据目标区储层的体积模量和剪切模量，确定目标区储层的岩性敏感弹性参数。

在一些实施例中，终端通过弹性阻抗模型进行弹性阻抗反演，得到目标区储层的弹性阻抗反演数据的操作至少包括：根据测井资料，将目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；根据角度部分叠加地震数据和测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；根据角度部分叠加地震数据、三个角度的地震子波和弹性阻抗模型，建立弹性阻抗反演目标函数；对弹性阻抗反演目标函数通过迭代算法进行求解处理，得到目标区储层的弹性阻抗反演数据。

在一些实施例中，终端根据测井资料，将目标区储层的三维地震数据进行角度叠加是指将第一角度范围的地震数据、第二角度范围的地震数据和第三角度范围的地震数据进行角度部分叠加。

需要说明的是，角度部分叠加中，该第一角度范围与第一角度弹性阻抗曲线的第一角度范围相同，第二角度范围与第二角度弹性阻抗曲线的第二角度范围相同，第三角度范围与第三角度弹性阻抗曲线的第三角度范围相同。

在一些实施例中，终端根据角度部分叠加地震数据和测井资料中的弹性阻抗曲线进行井震标定，并提取三个角度的地震子波，该三个角度的地震子波包括第一角度范围的地震子波、第二角度范围的地震子波和第三角度范围的地震子波。

需要说明的是，终端根据角度部分叠加地震数据和测井资料中的弹性阻抗曲线进行井震标定，并提取地震子波的操作能够参考相关技术，本申请实施例对此不再进行一一赘述。

在一些实施例中，终端根据角度部分叠加地震数据、三个角度的地震子波和弹性阻抗模型，建立的弹性阻抗反演目标函数F₁能够如下述第二公式所示。

需要说明的是，d为角度部分叠加地震数据，E为待确定(或为待求取)的弹性阻抗反演数据，G为三个角度的地震子波建立的矩阵，EI_model为目标区储层的弹性阻抗模型，α为弹性阻抗模型在反演过程中的约束权重系数，取值在0-1之间，取值越大代表弹性阻抗模型在弹性阻抗反演目标函数中的约束作用越大。β为稀疏约束权重系数，取值在0-1之间，取值越大代表反演结果的稀疏性越强。

在一些实施例中，在对弹性阻抗反演目标函数通过迭代算法进行求解处理时，当迭代结果不发生变化，或发生的变化小于预设变化阈值时，确定此时得到的迭代结果为目标区储层的弹性阻抗反演数据。或者，当迭代此时达到次数阈值时，确定得到的迭代结果为目标区储层的弹性阻抗反演数据。

需要说明的是，该预设变化阈值与次数阈值均能够根据需求事先进行设置，比如，该预设变化阈值能够为1、2等等，次数阈值能够为10次、15次等等。

还需要说明的是，终端在对弹性阻抗反演目标函数通过迭代算法进行求解处理的过程能够参考相关技术，本申请实施例对此不再进行一一赘述。

由于目标区储层的体积模量和剪切模量能够反应目标区储层的岩性敏感弹性特征，因此，终端需要根据弹性阻抗反演数据和目标区储层的测井资料，确定目标区储层的体积模量和剪切模量。

作为一种示例，终端根据弹性阻抗反演数据和目标区储层的测井资料，确定目标区储层的体积模量和剪切模量的操作至少包括：根据测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出测井资料中的弹性阻抗曲线与目标测井的体积模量和剪切模量之间的非线性模量转换系数，目标测井为目标区储层范围内的任一测井；根据非线性模量转换系数和目标区储层的弹性阻抗反演数据，确定目标区储层的体积模量和剪切模量。

由于测井资料中记载目标测井的体积模量和剪切模量可能不能准确描述目标区储层的岩性敏感岩性特征，因此，终端需要通过测井资料中目标测井的体积模量和剪切模量进行地震数据的转换，以确定目标区储层的体积模量和剪切模量。

在一种实施例中，该模量转换系数反演函数F₂能够为如下述第三公式所示的函数。

需要说明的是，在上述第三公式(3)中，K_well、μ_well和ρ_well分别为测井资料中目标测井的体积模量、剪切模量和密度信息(该密度信息能够为密度曲线)，EI_1well、EI_2well和EI_3well分别为测井资料中第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，

和

均为模量转换系数。

在一些实施例中，终端能够在模量转换系数反演函数F₂的值最小时，确定对应的模量转换系数为测井资料中的弹性阻抗曲线与体积模量和剪切模量之间的非线性模量转换系数。

由于弹性阻抗曲线能够描述弹性阻抗数据，且弹性阻抗曲线与体积模量和剪切模量之间的关系能够通过非线性模量转换系数进行转换，因此，终端在确定非线性模量转换系数后，能够根据非线性模量转换系数和目标区储层的弹性阻抗反演数据，确定出目标区储层的体积模量和剪切模量。

在一些实施例中，终端能够根据非线性模量转换系数和目标区储层的弹性阻抗反演数据，通过下述第四公式确定目标区储层的体积模量和剪切模量。

第四公式(4)：

需要说明的是，在上述第四公式(4)中，EI₁、EI₂和EI₃分别为弹性阻抗反演数据中的第一角度弹性角度反演数据、第二角度弹性阻抗反演数据和第三角度弹性阻抗反演数据，K为目标区储层的体积模量，μ为目标区储层的剪切模量。

在一些实施例中，终端根据目标区储层的体积模量和剪切模量，确定目标区储层的岩性敏感弹性参数的操作包括：根据测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，数据点为目标区存储范围内具有测井资料的位置，参考体积模量为对所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，参考剪切模量为对所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；根据目标区储层的体积模量和剪切模量，以及参考体积模量和参考剪切模量，确定目标区储层的岩性敏感弹性参数。

由于针对目标区储层可能会通过多个测井进行多次测井操作，从而目标区储层的范围内会包括多个测井，不同测井所在位置体积模量和参考剪切模量的能相同也可能不同，因此，为了准确的确定目标区储层的岩性敏感弹性参数，终端能够根据测井资料中的岩石孔隙度、各种矿物组分的弹性模量及体积比例等信息，对目标区储层范围内的所有数据点的体积模量和剪切模量分别进行统计规律分析与加权平均，从而得到参考体积模量和参考剪切模量。

在一些实施例中，终端根据目标区储层的体积模量和剪切模量，以及参考体积模量和参考剪切模量，能够通过下述第五公式确定目标区储层的岩性敏感弹性参数。

需要说明的是，在上述第五公式(5)中，K为目标区储层的体积模量，μ为目标区储层的剪切模量，K_s为参考体积模量，μ_s为参考剪切模量，C为岩性敏感弹性参数，即为目标区储层的泊松比。

步骤204：终端通过岩性敏感弹性参数模拟目标区储层的岩性三维空间分布特征。

由于岩性敏感弹性参数能够反应目标区储层的岩性三维空间分布特征，因此，终端能够通过岩性敏感弹性参数模拟目标区储层的岩性三维空间分布特征。

需要说明的是，终端通过岩性敏感弹性参数模拟目标区储层的岩性三维空间分布特征的操作能够参考相关技术，本申请实施例对此不再进行一一赘述。

在本申请实施例中，终端能够通过地震相编码、地层格架编码和测井资料进行弹性阻抗模型的建模，并根据弹性阻抗模型进行反演，得到岩性敏感弹性参数，即通过将地震反射特征的空间变化规律与测井资料的纵向特征结合开展建模和约束反演，提高了确定岩性敏感弹性参数的准确性，能够改善复杂储层中岩性特征横向变化快、非均质性严重的问题，提高了岩性储层预测的精度。

图5是本申请实施例提供的一种岩性敏感弹性参数的反演装置的结构示意图，该岩性敏感弹性参数的反演装置可以由软件、硬件或者两者的结合实现。该岩性敏感弹性参数的反演装置可以包括：获取模块501、确定模块502和反演模块503。

获取模块501，用于获取目标区储层的地震相编码、地层格架编码和测井资料；

确定模块502，用于根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，所述弹性阻抗模型用于描述所述目标区储层的地震反射特征空间变化规律；

反演模块503，用于根据所述弹性阻抗模型进行反演处理，得到所述目标区的岩性敏感弹性参数，所述岩性敏感弹性参数用于模拟所述目标区储层的岩性三维空间分布特征。

在一些实施例中，参见图6，所述获取模块501包括：

第一获取子模块5011，用于获取所述目标区储层的三维地震数据和地质信息；

第一构建子模块5012，用于根据所述三维地震数据，构建所述地震相编码；

第二构建子模块5013，用于根据所述三维地震数据和所述地质信息，构建所述地层格架编码。

在一些实施例中，所述第一构建子模块5012用于：

从所述目标区储层的三维地震数据中提取地震属性数据体，所述地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；

通过聚类算法对所述地震属性数据体进行聚类处理，得到所述目标区储层的地震相类型以及每种地震相的分布范围；

根据所述地震相类型和所述每种地震相的分布范围，对每种地震相类型进行编码，得到所述地震相编码。

在一些实施例中，所述第二构建子模块5013用于：

根据所述三维地震数据和所述地质信息，获取所述目标区储层的地震层位、断层分布特征和空间组合关系；

根据所述目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；

确定所述地震层位与所述断层趋势面之间的接触关系；

根据所述地震层位与所述断层趋势面之间的接触关系，对所述目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；

对相邻两个层位之间的离散化网格用数字进行编码，得到所述目标区储层的地层格架编码。

在一些实施例中，参见图7，所述确定模块502包括：

第二获取子模块5021，用于从所述测井资料中获取所述目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，所述第一角度弹性阻抗曲线的第一角度范围中最小角度大于所述第二角度弹性阻抗曲线的第二角度范围中最大角度，所述第二角度弹性阻抗曲线的第三角度范围中最小角度大于所述第三角度弹性阻抗曲线的角度范围中最大角度；

第三获取子模块5022，用于从所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；

第一确定子模块5023，用于根据所述每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；

插值子模块5024，用于通过所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线进行空间插值外推，得到所述目标位置处的弹性阻抗数值，所述目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；

第二确定子模块5025，用于根据所述同一地震相类型对应的参考弹性阻抗数值和所述目标位置的弹性阻抗数值，确定所述目标区储层的弹性阻抗模型。

在一些实施例中，参见图8，所述反演模块503包括：

反演子模块5031，用于通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据；

第三确定子模块5032，用于根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量；

第四确定子模块5033，用于根据所述目标区储层的体积模量和剪切模量，确定所述目标区储层的岩性敏感弹性参数。

在一些实施例中，所述反演子模块5031用于：

根据所述测井资料，将所述目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；

根据所述角度部分叠加地震数据和所述测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；

根据所述角度部分叠加地震数据、所述三个角度的地震子波和所述弹性阻抗模型，建立弹性阻抗反演目标函数；

对所述弹性阻抗反演目标函数通过迭代算法进行求解处理，得到所述目标区储层的弹性阻抗反演数据。

在一些实施例中，所述第三确定子模块5032用于：

根据所述测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出所述测井资料中的弹性阻抗曲线与所述目标测井的体积模量和剪切模量之间的非线性模量转换系数，所述目标测井为所述目标区储层范围内的任一测井；

根据所述非线性模量转换系数和所述目标区储层的弹性阻抗反演数据，确定所述目标区储层的体积模量和剪切模量。

在一些实施例中，所述第四确定子模块5033用于：

根据所述测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定所述目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，所述数据点为所述目标区存储范围内具有测井资料的位置，所述参考体积模量为对所述所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，所述参考剪切模量为对所述所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；

根据所述目标区储层的体积模量和剪切模量，以及所述参考体积模量和所述参考剪切模量，确定所述目标区储层的岩性敏感弹性参数。

在本申请实施例中，终端能够通过地震相编码、地层格架编码和测井资料进行弹性阻抗模型的建模，并根据弹性阻抗模型进行反演，得到岩性敏感弹性参数，即通过将地震反射特征的空间变化规律与测井资料的纵向特征结合开展建模和约束反演，提高了确定岩性敏感弹性参数的准确性，能够改善复杂储层中岩性特征横向变化快、非均质性严重的问题，提高了岩性储层预测的精度。

需要说明的是：上述实施例提供的岩性敏感弹性参数的反演装置在进行岩性敏感弹性参数的反演时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的岩性敏感弹性参数的反演装置与岩性敏感弹性参数的反演方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。

图9示出了本申请一个示例性实施例提供的终端900的结构框图。该终端900可以是：智能手机、平板电脑、MP3播放器(Moving Picture Experts Group Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts Group Audio LayerIV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端900还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。

通常，终端900包括有：处理器901和存储器902。

处理器901可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器901可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器901也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器901可以集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器901还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。

存储器902可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器902还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器902中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器901所执行以实现本申请中方法实施例提供的岩性敏感弹性参数的反演方法。

在一些实施例中，终端900还可选包括有：外围设备接口903和至少一个外围设备。处理器901、存储器902和外围设备接口903之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口903相连。具体地，外围设备包括：射频电路904、显示屏905、摄像头组件906、音频电路907、定位组件908和电源909中的至少一种。

外围设备接口903可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器901和存储器902。在一些实施例中，处理器901、存储器902和外围设备接口903被集成在同一芯片或电路板上；在一些其他实施例中，处理器901、存储器902和外围设备接口903中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。

射频电路904用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路904通过电磁信号与通信网络以及其他通信设备进行通信。射频电路904将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路904包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路904可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路904还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。

显示屏905用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏905是触摸显示屏时，显示屏905还具有采集在显示屏905的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器901进行处理。此时，显示屏905还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏905可以为一个，设置终端900的前面板；在另一些实施例中，显示屏905可以为至少两个，分别设置在终端900的不同表面或呈折叠设计；在另一些实施例中，显示屏905可以是柔性显示屏，设置在终端900的弯曲表面上或折叠面上。甚至，显示屏905还可以设置成非矩形的不规则图形，也即异形屏。显示屏905可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。

摄像头组件906用于采集图像或视频。可选地，摄像头组件906包括前置摄像头和后置摄像头。通常，前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件906还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。

音频电路907可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器901进行处理，或者输入至射频电路904以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端900的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器901或射频电路904的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路907还可以包括耳机插孔。

定位组件908用于定位终端900的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件908可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。

电源909用于为终端900中的各个组件进行供电。电源909可以是交流电、直流电、一次性电池或可充电电池。当电源909包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。

在一些实施例中，终端900还包括有一个或多个传感器910。该一个或多个传感器910包括但不限于：加速度传感器911、陀螺仪传感器912、压力传感器913、指纹传感器914、光学传感器915以及接近传感器916。

加速度传感器911可以检测以终端900建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器911可以用于检测重力加速度在三个坐标轴上的分量。处理器901可以根据加速度传感器911采集的重力加速度信号，控制显示屏905以横向视图或纵向视图进行用户界面的显示。加速度传感器911还可以用于游戏或者用户的运动数据的采集。

陀螺仪传感器912可以检测终端900的机体方向及转动角度，陀螺仪传感器912可以与加速度传感器911协同采集用户对终端900的3D动作。处理器901根据陀螺仪传感器912采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。

压力传感器913可以设置在终端900的侧边框和/或显示屏905的下层。当压力传感器913设置在终端900的侧边框时，可以检测用户对终端900的握持信号，由处理器901根据压力传感器913采集的握持信号进行左右手识别或快捷操作。当压力传感器913设置在显示屏905的下层时，由处理器901根据用户对显示屏905的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。

指纹传感器914用于采集用户的指纹，由处理器901根据指纹传感器914采集到的指纹识别用户的身份，或者，由指纹传感器914根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器901授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器914可以被设置终端900的正面、背面或侧面。当终端900上设置有物理按键或厂商Logo时，指纹传感器914可以与物理按键或厂商Logo集成在一起。

光学传感器915用于采集环境光强度。在一个实施例中，处理器901可以根据光学传感器915采集的环境光强度，控制显示屏905的显示亮度。具体地，当环境光强度较高时，调高显示屏905的显示亮度；当环境光强度较低时，调低显示屏905的显示亮度。在另一个实施例中，处理器901还可以根据光学传感器915采集的环境光强度，动态调整摄像头组件906的拍摄参数。

接近传感器916，也称距离传感器，通常设置在终端900的前面板。接近传感器916用于采集用户与终端900的正面之间的距离。在一个实施例中，当接近传感器916检测到用户与终端900的正面之间的距离逐渐变小时，由处理器901控制显示屏905从亮屏状态切换为息屏状态；当接近传感器916检测到用户与终端900的正面之间的距离逐渐变大时，由处理器901控制显示屏905从息屏状态切换为亮屏状态。

本领域技术人员可以理解，图9中示出的结构并不构成对终端900的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。

本申请实施例还提供了一种非临时性计算机可读存储介质，当所述存储介质中的指令由终端的处理器执行时，使得终端能够执行上实施例提供的岩性敏感弹性参数的反演方法。

本申请实施例还提供了一种包含指令的计算机程序产品，当其在终端上运行时，使得终端执行上述实施例提供的岩性敏感弹性参数的反演方法。

本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

以上所述仅为本申请实施例的较佳实施例，并不用以限制本申请实施例，凡在本申请实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (19)

1.一种岩性敏感弹性参数的反演方法，其特征在于，所述方法包括：

获取目标区储层的地震相编码、地层格架编码和测井资料；

根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，所述弹性阻抗模型用于描述所述目标区储层的地震反射特征空间变化规律；

根据所述弹性阻抗模型进行反演处理，得到所述目标区的岩性敏感弹性参数，所述岩性敏感弹性参数用于模拟所述目标区储层的岩性三维空间分布特征。

2.如权利要求1所述的方法，其特征在于，获取目标区储层的地震相编码和地层格架编码，包括：

获取所述目标区储层的三维地震数据和地质信息；

根据所述三维地震数据，构建所述地震相编码；

根据所述三维地震数据和所述地质信息，构建所述地层格架编码。

3.如权利要求2所述的方法，其特征在于，所述根据所述三维地震数据，构建所述地震相编码，包括：

从所述目标区储层的三维地震数据中提取地震属性数据体，所述地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；

通过聚类算法对所述地震属性数据体进行聚类处理，得到所述目标区储层的地震相类型以及每种地震相的分布范围；

根据所述地震相类型和所述每种地震相的分布范围，对每种地震相类型进行编码，得到所述地震相编码。

4.如权利要求2所述的方法，其特征在于，所述根据所述三维地震数据和所述地质信息，构建所述地层格架编码，包括：

根据所述三维地震数据和所述地质信息，获取所述目标区储层的地震层位、断层分布特征和空间组合关系；

根据所述目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；

确定所述地震层位与所述断层趋势面之间的接触关系；

根据所述地震层位与所述断层趋势面之间的接触关系，对所述目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；

对相邻两个层位之间的离散化网格用数字进行编码，得到所述目标区储层的地层格架编码。

5.如权利要求1-4任一权利要求所述的方法，其特征在于，所述根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，包括：

从所述测井资料中获取所述目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，所述第一角度弹性阻抗曲线的第一角度范围中最小角度大于所述第二角度弹性阻抗曲线的第二角度范围中最大角度，所述第二角度弹性阻抗曲线的第三角度范围中最小角度大于所述第三角度弹性阻抗曲线的角度范围中最大角度；

从所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；

根据所述每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；

通过所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线进行空间插值外推，得到所述目标位置处的弹性阻抗数值，所述目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；

根据所述同一地震相类型对应的参考弹性阻抗数值和所述目标位置的弹性阻抗数值，确定所述目标区储层的弹性阻抗模型。

6.如权利要求1所述的方法，其特征在于，所述根据所述弹性阻抗模型进行反演处理，得到所述目标区储层的岩性敏感弹性参数，包括：

通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据；

根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量；

根据所述目标区储层的体积模量和剪切模量，确定所述目标区储层的岩性敏感弹性参数。

7.如权利要求6所述的方法，其特征在于，所述通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据，包括：

根据所述测井资料，将所述目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；

根据所述角度部分叠加地震数据和所述测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；

根据所述角度部分叠加地震数据、所述三个角度的地震子波和所述弹性阻抗模型，建立弹性阻抗反演目标函数；

对所述弹性阻抗反演目标函数通过迭代算法进行求解处理，得到所述目标区储层的弹性阻抗反演数据。

8.如权利要求6所述的方法，其特征在于，所述根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量，包括：

根据所述测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出所述测井资料中的弹性阻抗曲线与所述目标测井的体积模量和剪切模量之间的非线性模量转换系数，所述目标测井为所述目标区储层范围内的任一测井；

根据所述非线性模量转换系数和所述目标区储层的弹性阻抗反演数据，确定所述目标区储层的体积模量和剪切模量。

9.如权利要求6-8任一权利要求所述的方法，其特征在于，所述根据所述目标区储层的体积模量和剪切模量数据，确定所述目标区储层的岩性敏感弹性参数，包括：

根据所述测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定所述目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，所述数据点为所述目标区存储范围内具有测井资料的位置，所述参考体积模量为对所述所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，所述参考剪切模量为对所述所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；

根据所述目标区储层的体积模量和剪切模量，以及所述参考体积模量和所述参考剪切模量，确定所述目标区储层的岩性敏感弹性参数。

10.一种岩性敏感弹性参数的反演装置，其特征在于，所述装置包括：

获取模块，用于获取目标区储层的地震相编码、地层格架编码和测井资料；

确定模块，用于根据所述目标区储层的地震相编码、地层格架编码和测井资料，确定所述目标区储层的弹性阻抗模型，所述弹性阻抗模型用于描述所述目标区储层的地震反射特征空间变化规律；

反演模块，用于根据所述弹性阻抗模型进行反演处理，得到所述目标区的岩性敏感弹性参数，所述岩性敏感弹性参数用于模拟所述目标区储层的岩性三维空间分布特征。

11.如权利要求10所述的装置，其特征在于，所述获取模块包括：

第一获取子模块，用于获取所述目标区储层的三维地震数据和地质信息；

第一构建子模块，用于根据所述三维地震数据，构建所述地震相编码；

第二构建子模块，用于根据所述三维地震数据和所述地质信息，构建所述地层格架编码。

12.如权利要求11所述的装置，其特征在于，所述第一构建子模块用于：

从所述目标区储层的三维地震数据中提取地震属性数据体，所述地震属性数据体包括均方根振幅数据、相对阻抗数据和绝对振幅数据；

通过聚类算法对所述地震属性数据体进行聚类处理，得到所述目标区储层的地震相类型以及每种地震相的分布范围；

根据所述地震相类型和所述每种地震相的分布范围，对每种地震相类型进行编码，得到所述地震相编码。

13.如权利要求11所述的装置，其特征在于，所述第二构建子模块用于：

根据所述三维地震数据和所述地质信息，获取所述目标区储层的地震层位、断层分布特征和空间组合关系；

根据所述目标区储层的断层分布特征和空间组合关系，模拟断层趋势面；

确定所述地震层位与所述断层趋势面之间的接触关系；

根据所述地震层位与所述断层趋势面之间的接触关系，对所述目标区储层对应的三维空间采用均匀网格进行离散化划分，得到多个离散化网格；

对相邻两个层位之间的离散化网格用数字进行编码，得到所述目标区储层的地层格架编码。

14.如权利要求10-13任一权利要求所述的装置，其特征在于，所述确定模块包括：

第二获取子模块，用于从所述测井资料中获取所述目标区储层的第一角度弹性阻抗曲线、第二角度弹性阻抗曲线和第三角度弹性阻抗曲线，所述第一角度弹性阻抗曲线的第一角度范围中最小角度大于所述第二角度弹性阻抗曲线的第二角度范围中最大角度，所述第二角度弹性阻抗曲线的第三角度范围中最小角度大于所述第三角度弹性阻抗曲线的角度范围中最大角度；

第三获取子模块，用于从所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线获取每个地震相编码对应的弹性阻抗数值；

第一确定子模块，用于根据所述每个地震相编码对应的弹性阻抗数值，确定同一地震相类型对应的参考弹性阻抗数值；

插值子模块，用于通过所述第一角度弹性阻抗曲线、所述第二角度弹性阻抗曲线和所述第三角度弹性阻抗曲线进行空间插值外推，得到所述目标位置处的弹性阻抗数值，所述目标位置为同时具有同一个地震相编码和同一个地层格架编码的离散化地层网格对应的位置；

第二确定子模块，用于根据所述同一地震相类型对应的参考弹性阻抗数值和所述目标位置的弹性阻抗数值，确定所述目标区储层的弹性阻抗模型。

15.如权利要求10所述的装置，其特征在于，所述反演模块包括：

反演子模块，用于通过所述弹性阻抗模型进行弹性阻抗反演，得到所述目标区储层的弹性阻抗反演数据；

第三确定子模块，用于根据所述弹性阻抗反演数据和所述目标区储层的测井资料，确定所述目标区储层的体积模量和剪切模量；

第四确定子模块，用于根据所述目标区储层的体积模量和剪切模量，确定所述目标区储层的岩性敏感弹性参数。

16.如权利要求15所述的装置，其特征在于，所述反演子模块用于：

根据所述测井资料，将所述目标区储层的三维地震数据进行角度叠加，得到角度部分叠加的地震数据；

根据所述角度部分叠加地震数据和所述测井资料中的弹性阻抗曲线进行井震标定，得到三个角度的地震子波；

根据所述角度部分叠加地震数据、所述三个角度的地震子波和所述弹性阻抗模型，建立弹性阻抗反演目标函数；

对所述弹性阻抗反演目标函数通过迭代算法进行求解处理，得到所述目标区储层的弹性阻抗反演数据。

17.如权利要求15所述的装置，其特征在于，所述第三确定子模块用于：

根据所述测井资料中目标测井的体积模量、剪切模量和密度信息，通过模量转换系数反演函数，反演出所述测井资料中的弹性阻抗曲线与所述目标测井的体积模量和剪切模量之间的非线性模量转换系数，所述目标测井为所述目标区储层范围内的任一测井；

根据所述非线性模量转换系数和所述目标区储层的弹性阻抗反演数据，确定所述目标区储层的体积模量和剪切模量。

18.如权利要求15-17任一权利要求所述的装置，其特征在于，所述第四确定子模块用于：

根据所述测井资料中的岩石孔隙度、矿物组分的弹性模量及体积比例，确定所述目标区储层范围内的所有数据点的参考体积模量和参考剪切模量，所述数据点为所述目标区存储范围内具有测井资料的位置，所述参考体积模量为对所述所有数据点的体积模量进行加权平均后得到的体积模量统计平均值，所述参考剪切模量为对所述所有数据点的剪切模量进行加权平均后得到的剪切模量统计平均值；

根据所述目标区储层的体积模量和剪切模量，以及所述参考体积模量和所述参考剪切模量，确定所述目标区储层的岩性敏感弹性参数。

19.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有指令，所述指令被处理器执行时实现上述权利要求1至权利要求9中的任一项权利要求所述的方法的步骤。

Images