一种储层特征参数的预测方法及装置与流程

本申请涉及地球物理勘探技术领域，特别涉及一种储层特征参数的预测方法及装置。

背景技术：

储层特征参数的预测，是对储层进行预测的基础。准确而全面的储层特征参数数据，是能够精确预测储层的保证。因此，在储层预测过程中，储层特征参数的预测是一个重要环节。

现有技术中，对储层特征参数的预测，主要是采用交会图技术。利用所述交会图技术，可以较为准确地预测出储层的孔隙度和饱和度，进而可以较为准确地预测出储层的结构等特性。但是，对于一些致密储层来说，由于储层的孔隙结构比较复杂，可能具有两种甚至两种以上的孔隙结构。这个时候，仅仅依靠储层的孔隙度和饱和度两种参数，就无法更准确地预测储层的真实结构了。而现有技术中，没有同时预测储层的孔隙度、饱和度和储层孔隙结构特征参数的方法。

现有技术中至少存在如下问题：

现有技术只预测得到了储层的孔隙度和饱和度两种储层特征参数，没有得到反映储层的孔隙结构的特征参数，而且现有技术中也没有同时预测储层的孔隙度、饱和度和储层孔隙结构特征参数的方法。这样就导致了在利用得到的储层特征参数进行储层预测时，不能准确真实地预测出储层的真实结构。

技术实现要素：

本申请实施例的目的是提供一种储层特征参数的预测方法及装置，以构建一种三维岩石物理模板，提供更全面更准确的用于储层预测的特征参数数据。

本申请实施例提供一种储层特征参数的预测方法及装置是这样实现的：

一种储层特征参数的预测方法，所述方法包括：

获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，所述三种地震响应特征参数数据通过地震叠前反演得到；

将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测，所述三维岩石物理模版包括采用下述构建方式生成：确定储层特征参数的分布范围，对所述分布范围进行网格剖分，生成三维网格，计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度，利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的节点确定出所述三维岩石物理模版；

输出预测得到的待预测位置的储层特征参数数据。

优选实施例中，所述对所述分布范围进行网格剖分，生成三维网格的方式包括：

分别以三种预定储层特征参数作为三个坐标轴，建立第一个三维坐标系；

分别在所述三个坐标轴上对三种储层特征参数的分布范围进行划分，得到所述三维网格。

优选实施例中，所述利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的节点确定出所述三维岩石物理模版的方式，包括：

建立分别以纵波阻抗、横波阻抗和密度为三个坐标轴的第二个三维坐标系；

以所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度为坐标，将所述三维网格的各节点定位于所述第二个三维坐标系中，得到定位后的节点；

以所述定位后的节点确定出所述三维岩石物理模版。

优选实施例中，所述储层特征参数包括：

储层的孔隙扁度，储层的孔隙度，储层的饱和度。

优选实施例中，所述计算所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗的方式包括：

获取储层的岩石基质的体积模量、剪切模量和密度；

获取储层的孔隙中流体的体积模量和密度；

利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量；

利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量，进而计算得到储层孔隙中流体的密度；

利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量；

利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度；

利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。

优选实施例中，采用包括下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：

式中，Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

Gm表示储层的岩石基质的剪切模量；

φ表示各节点对应的孔隙度；

α表示各节点的孔隙扁度；

a和b满足：

和表示干岩石骨架的极化因子；

K^*(y)和G^*(y)表示岩石等效体积模量和剪切模量。

优选实施例中，采用包括下述计算公式计算得到储层孔隙中的体积模量：

Kf＝SWKW+(1-SW)Kg

式中，Kf表示所述储层孔隙中流体的体积模量；

Sw表示所述各节点对应的饱和度；

Kw表示储层孔隙中水的体积模量；

Kg表示储层孔隙中气体的体积模量。

优选实施例中，采用包括下述计算公式计算得到储层孔隙中流体的密度：

ρf＝Swρw+(1-Sw)ρg

式中，ρf表示所述储层孔隙中流体的密度；

Sw表示所述各节点对应的饱和度；

ρw表示储层孔隙中水的密度；

ρg表示储层孔隙中气体的密度。

优选实施例中，采用包括下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：

Gsat＝Gdry

式中，Ksat表示所述储层的饱和流体岩石的体积模量；

Gsat表示所述储层的饱和流体岩石的体积模量；

Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

φ表示各节点对应的孔隙度；

Kf表示所述储层孔隙中流体的体积模量。

优选实施例中，采用包括下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：

ρsat＝(1-φ)ρm+φρf；

式中，ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度；

Gsat表示所述储层的饱和流体岩石的体积模量；

Ksat表示所述储层的饱和流体岩石的体积模量；

ρm表示储层的岩石基质的密度；

ρf表示所述储层孔隙中流体的密度；

φ表示各节点对应的孔隙度。

优选实施例中，采用包括下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：

Ip＝ρVp；

Is＝ρVs；

ρ＝ρsat

式中，Ip表示所述各节点对应的纵波阻抗；

Is表示所述各节点对应的横波阻抗；

ρ表示所述各节点对应的密度；

ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度。

优选实施例中，所述将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置进行储层特征参数预测的方式，包括：

利用所述待预测位置的三种地震响应特征参数数据，在所述三维岩石物理模版中，确定出所述待预测位置的储层特征参数数据。

优选实施例中，所述装置包括：

数据获取模块，用于获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据；

模版构建模块，用于构建预先构建好的三维岩石物理模版；

参数预测模块，用于将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测；

数据输出模块，用于输出预测得到的待预测位置的储层特征参数数据。

优选实施例中，所述模版构建模块包括：

分布范围确定模块，用于确定储层特征参数的分布范围；

网格生成模块，用于对所述储层特征参数的分布范围进行网格剖分，生成三维网格；

数据计算模块，用于计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度；

模版生成模块，用于利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，构建所述三维岩石物理模版。

优选实施例中，所述数据计算模块包括：

模量计算模块，用于利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量，还用于利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量，还用于利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量；

地震响应特征参数计算模块，用于利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度，还用于计算得到储层孔隙中流体的密度；

波阻抗计算模块，用于利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。

优选实施例中，采用包括下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：

式中，Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

Gm表示储层的岩石基质的剪切模量；

φ表示各节点对应的孔隙度；

α表示各节点的孔隙扁度；

a和b满足：

和表示干岩石骨架的极化因子；

K^*(y)和G^*(y)表示岩石等效体积模量和剪切模量。

优选实施例中，采用包括下述计算公式计算得到储层孔隙中的体积模量：

Kf＝SwKw+(1-Sw)Kg

式中，Kf表示所述储层孔隙中流体的体积模量；

Sw表示所述各节点对应的饱和度；

Kw表示储层孔隙中水的体积模量；

Kg表示储层孔隙中气体的体积模量。

优选实施例中，采用包括下述计算公式计算得到储层孔隙中流体的密度：

ρf＝Swρw+(1-Sw)ρg

式中，ρf表示所述储层孔隙中流体的密度；

Sw表示所述各节点对应的饱和度；

ρw表示储层孔隙中水的密度；

ρg表示储层孔隙中气体的密度。

优选实施例中，采用包括下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：

Gsat＝Gdry

式中，Ksat表示所述储层的饱和流体岩石的体积模量；

Gsat表示所述储层的饱和流体岩石的体积模量；

Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

φ表示各节点对应的孔隙度；

Kf表示所述储层孔隙中流体的体积模量。

优选实施例中，采用包括下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：

ρsat＝(1-φ)ρm+φρf；

式中，ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度；

Gsat表示所述储层的饱和流体岩石的体积模量；

Ksat表示所述储层的饱和流体岩石的体积模量；

ρm表示储层的岩石基质的密度；

ρf表示所述储层孔隙中流体的密度；

φ表示各节点对应的孔隙度。

优选实施例中，采用包括下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：

Ip＝ρVp；

Is＝ρVs；

ρ＝ρsat

式中，Ip表示所述各节点对应的纵波阻抗；

Is表示所述各节点对应的横波阻抗；

ρ表示所述各节点对应的密度；

ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度。

利用本申请实施例提供的一种储层特征参数的预测方法，可以通过构建的三维岩石物理模版，准确预测出储层的孔隙扁度、饱和度和孔隙度。其中，所述储层的孔隙扁度可以反映储层的孔隙结构。利用所述三维岩石物理模版预测得到的三种储层特征参数，不仅可以准确地反映储层的孔隙大小以及储层对流体的储存特性，还可以准确预测出储层的孔隙结构，这样就可以更全面地反映出储层的结构特性。进而有效提高储层的结构分析结果的真实性和准确度。利用本申请实施例提供的一种储层特征参数的预测装置，可以自动实施所述储层特征参数的预测方法，并自动输出预测结果，操作简单快捷，不需要实施人员的过多参与，有效提高了用户体验。

附图说明

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

图1是本申请一个实施例中提供的一种储层特征参数的预测方法的方法流程图；

图2是本申请一个实施例中提供的一种三维岩石物理模版的构建方式的方法流程图；

图3是本申请一个实施例中提供的一种确定出所述三维岩石物理模版的方法流程图；

图4是本申请一个实施例中提供的一种储层特征参数的预测装置的模块结构示意图；

图5是本申请一个实施例中提供的模版构建模块的模块结构示意图；

图6(a)是本申请一个实施例中得到的预先构建好的三维岩石物理模版；

图6(b)是本申请一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比；

图7(a)是本申请另一个实施例中得到的预先构建好的三维岩石物理模版；

图7(b)是本申请另一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比；

图8(a)是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙度连井剖面；

图8(b)是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的饱和度连井剖面；

图8(c)是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙扁度连井剖面；

图9是本申请一个实施例中提供的计算处理模块的模块结构示意图。

具体实施方式

本申请实施例提供一种储层特征参数的预测方法及装置。

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

图1是本申请所述一种储层特征参数的预测方法一种实施例的方法流程图。虽然本申请提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本申请实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理的实施环境)。

具体的，如图1所示，本申请提供的一种储层特征参数的预测方法的一种实施例可以包括：

S1：获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，所述三种地震响应特征参数数据包括通过地震叠前反演得到。

所述储层中待预测位置可以包括储层中的任意位置。

所述待预测位置可以只是一个位置，也可以是N个位置，N≥2。

在本申请一个实施例中，所述待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据，可以通过地震叠前反演得到。

当然，具体的获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据的方式不必限定，在本申请其他实施例中，也可以通过其他方式获取，比如在本申请另一个实施例中，可以通过直接测量的方式得到所述三种地震响应特征数据。

S2：将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测。

所述三维岩石物理模版主要由三维网格中的各节点确定出来。

所述储层特征参数至少包括储层的孔隙扁度、储层的孔隙度、储层的饱和度。

图2是本申请一个实施例中提供的一种三维岩石物理模版的构建方式的方法流程图，具体的，如图2所示，所述预先构建好的三维岩石物理模版的构建方式可以包括：

S201：确定储层特征参数的分布范围，对所述分布范围进行网格剖分，生成三维网格。

所述储层特征参数的分布范围，可以认为是所述储层特征参数的变化区间，每个储层特征参数的变化区间都是可以确定的，比如，储层的孔隙扁度的分布范围可以是大于等于0且小于等于1。所述确定储层特征参数的分布范围属于本领域技术人员可以掌握的常规技术，具体的确定方式不必限定。

所述对所述分布范围进行网格剖分，可以包括对每种储层特征参数的分布范围进行划分，所述划分的方式，可以是等间隔划分，也可以是不等间隔划分，划分的间隔的宽度也不必限定。比如，对于上述的储层的孔隙扁度，划分的间隔的宽度可以均取为0.1，比如所述储层的孔隙扁度的分布范围为[0,1]，划分之后，每一个小间隔可以依次为[0,0.1)、(0.1,0.2)……(0.9,1.0]。当然，划分的间隔的宽度越小，最终得到的三维岩石物理模版的预测精度就越高。具体的，在本申请其他实施例中，可以根据所需要的预测精度以及实际实施情形，确定所述对每种储层特征参数的分布范围进行划分的方式以及划分的间隔的宽度。

最终由每种储层特征参数的分布范围划分出的间隔组成三维网格，所述三维网格的每个节点对应一组储层特征参数数据，所述一组储层特征参数可以包括储层的孔隙扁度、孔隙度和饱和度。

S202：计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度。

所述三维网格中各节点的孔隙扁度、孔隙度和饱和度是已经确定的。

利用所述各节点的孔隙扁度、孔隙度和饱和度三种储层特征参数，结合储层的岩石基质的体积模量、剪切模量和密度，储层的孔隙中流体的体积模量和密度，就可以计算出所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度。

S203：利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的各节点确定出所述三维岩石物理模版。

所述对所述三维网格中各节点进行定位，可以包括：将计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度分别作为所述各节点的三个坐标，利用所述坐标，将所述各节点在同一个分别以纵波阻抗、横波阻抗和密度为坐标轴的三维坐标系中定位。

利用定位后的各节点，就可以确定出所述三维岩石物理模版。

所述确定出所述三维岩石物理模版的方式属于常规的数据处理手段，可以采用拟合等方式确定出所述三维岩石物理模版，具体的确定方式不必限定。

利用上述实施例中得到的所述三维岩石物理模版，可以预测出储层的三种储层特征参数，可以包括：储层的孔隙扁度，储层的孔隙度，储层的饱和度。

S3：输出预测得到的待预测位置的储层特征参数数据。

所述待预测位置的储层特征参数可以包括：储层的孔隙扁度，储层的孔隙度和储层的饱和度。

具体的，输出所述预测得到的待预测位置的储层特征参数数据，可以根据实施人员的需要选择性输出，比如可以同时输出待预测位置的孔隙扁度、孔隙度和饱和度的数据，也可以只输出孔隙度和饱和度的数据，当然，也可以只输出孔隙扁度一种参数的数据。

同样的，在数据输出阶段，所述待预测位置可以只是一个位置，也可以是N个位置，N≥2。也就是说，在本申请其他实施例中，可以输出储层中一个位置的储层特征参数数据，也可以同时输出储层中N个位置的储层特征参数数据，N≥2。

利用上述各实施例提供的一种储层特征参数的预测方法可以利用所述三维岩石物理模版预测得到的储层的孔隙扁度、储层的孔隙度和储层的饱和度三种储层特征参数，不仅可以准确地反映储层的孔隙大小以及储层对流体的储存特性，还可以准确预测出储层的孔隙结构，这样就可以更全面地反映出储层的结构特性。进而有效提高储层的结构分析结果的真实性和准确度。

在本申请另一个实施例中，对所述分布范围进行网格剖分，生成三维网格的方式可以包括：

分别以三种预定储层特征参数作为三个坐标轴，建立第一个三维坐标系。

本申请一个实施例中，所述三种预定储层特征参数可以包括储层的孔隙扁度、储层的孔隙度、储层的饱和度，选择上述三种储层特征参数，可以更全面地反映储层的结构特征，尤其是储层的孔隙特征。

当然，在本申请其他实施例中，实施人员也可以根据实际作业情况和作业需要，选择其他的储层特征参数作为三种预定储层特征参数，比如三种预定储层特征参数可以是储层的渗透率、储层的孔隙度和储层的饱和度。

分别在所述三个坐标轴上对三种储层特征参数的分布范围进行划分，得到所述三维网格。

所述分别在所述三个坐标轴上对三种储层特征参数的分布范围进行划分的方式，可以包括对在所述每种储层特征参数对应的坐标轴上，对所述每种储层特征参数的分布范围进行划分，所述划分的方式，可以是等间隔划分，也可以是不等间隔划分，划分的间隔的宽度也不必限定。

利用上述实施例，可以得到一个三维网格，并且三维网格中的每个节点对应着一组三种储层特征参数数据。

在本申请又一个实施例中，利用计算得到的三维网格中各节点对应的纵波阻抗、横波阻抗和密度，对所述三维网格中各节点进行定位，以定位后的各节点确定出所述三维岩石物理模版的方式如图3所示，具体的，可以包括：

S2031：建立分别以纵波阻抗、横波阻抗和密度为三个坐标轴的第二个三维坐标系。

S2032：以所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度为坐标，将所述三维网格的各节点定位于所述第二个三维坐标系中，得到定位后的各节点。

S2033：以所述定位后的各节点确定出所述三维岩石物理模版。

利用所述三维岩石物理模版，准确预测出储层的孔隙度、饱和度和孔隙扁度。

在本申请再一个实施例中，将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置进行储层特征参数预测的方式，可以包括：

利用所述待预测位置的三种地震响应特征参数数据，在所述三维岩石物理模版中，确定出所述待预测位置的储层特征参数数据。

具体的，确定出所述待预测位置的储层特征参数数据的方式，属于常规数学方法，可以以三种地震响应特征参数数据为坐标，在所述三维岩石物理模版中，确定出所述坐标对应的点，所述坐标对应的点对应的储层特征参数数据就是所述待预测位置的储层特征参数数据。

在本申请另一个实施例中，计算所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗的方式，可以包括：

获取储层的岩石基质的体积模量、剪切模量和密度。

获取储层的孔隙中流体的体积模量和密度。

利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量。

在本申请一个实施例中，可以包括采用下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：

式中，Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

Gm表示储层的岩石基质的剪切模量；

φ表示各节点对应的孔隙度；

α表示各节点的孔隙扁度；

a和b满足：

和表示干岩石骨架的极化因子；

K^*(y)和G^*(y)表示岩石等效体积模量和剪切模量。

利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量，进而计算得到储层孔隙中流体的密度。

在本申请一个实施例中，可以包括采用下述计算公式计算得到储层孔隙中流体的体积模量：

Kf＝SwKw+(1-Sw)Kg

式中，Kf表示所述储层孔隙中流体的体积模量；

Sw表示所述各节点对应的饱和度；

Kw表示储层孔隙中水的体积模量；

Kg表示储层孔隙中气体的体积模量。

在本申请一个实施例中，所述计算得到储层孔隙中流体的密度的计算公式可以包括：

ρf＝Swρw+(1-Sw)ρg

式中，ρf表示所述储层孔隙中流体的密度；

Sw表示所述各节点对应的饱和度；

ρw表示储层孔隙中水的密度；

ρg表示储层孔隙中气体的密度。

利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量。

在本申请一个实施例中，可以包括采用下述计算公式计算得到储层的流体饱和岩石的体积模量和剪切模量：

Gsat＝Gdry

式中，Ksat表示所述储层的饱和流体岩石的体积模量；

Gsat表示所述储层的饱和流体岩石的体积模量；

Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

φ表示各节点对应的孔隙度；

Kf表示所述储层孔隙中流体的体积模量。

利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度。

在本申请一个实施例中，可以包括采用下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：

ρsat＝(1-φ)ρm+φρf；

式中，ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度；

Gsat表示所述储层的饱和流体岩石的体积模量；

Ksat表示所述储层的饱和流体岩石的体积模量；

ρm表示储层的岩石基质的密度；

ρf表示所述储层孔隙中流体的密度。

利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。

在本申请一个实施例中，可以包括采用下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：

Ip＝ρVp；

Is＝ρVs；

ρ＝ρsat

式中，Ip表示所述各节点对应的纵波阻抗；

Is表示所述各节点对应的横波阻抗；

ρ表示所述各节点对应的密度；

ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度。

利用上述实施例提供的计算方式，可以计算出所述三维网格中各节点对应的密度、纵波阻抗和横波阻抗。

利用上述各实施例提供的一种储层特征参数的预测方法，可以得到的所述三维岩石物理模版，并利用所述三维岩石物理模版，准确预测出待预测位置的储层特征参数数据。

图6(a)是本申请一个实施例中得到的预先构建好的三维岩石物理模版，图6(b)是本申请一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。

图7(a)是本申请另一个实施例中得到的预先构建好的三维岩石物理模版，图7(b)是本申请另一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。从图6和图7中，可以看出，利用构建的三维岩石物理模版，可以准确预测出储层的孔隙扁度、饱和度和孔隙度。并且预测得到的结果的真实性和准确度得到了有效提高。

图6(b)是本申请一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。图中，纵坐标表示所述待预测位置的深度，单位是米(m)。图中黑色曲线表示的是测量得到的纵横波速度和密度曲线以及孔隙度和饱和度曲线，灰色曲线表示的是根据所述三维岩石物理模版预测得到的孔隙扁度及孔隙度和饱和度。可以看出，无论是储层的孔隙度还是储层的饱和度，测井解释和利用所述三维岩石物理模版预测得到的结果吻合度较高，这也证明了利用所述三维岩石物理模版进行储层特征参数预测具有较高的可靠性和准确度。从反演的结果看，储层段的孔隙扁度变化范围大，从0.1到1.0都存在分布，从图6(a)的三维岩石物理模版上也可以发现，这说明碳酸岩盐储层的纵横波速度不仅受孔隙度和饱和度的影响，也受孔隙扁度的影响，利用纵横波密度数据反演孔隙度和饱和度需要考虑孔隙扁度这一参数。

图7(b)是本申请另一个实施例中预测得到的孔隙扁度、孔隙度和饱和度与测井解释结果对比。图中，纵坐标表示所述待预测位置的深度，单位是米(m)。可以看到，利用所述三维岩石物理模版预测得到的孔隙度和饱和度与测井解释结果具有较高的吻合度。证明了利用所述三维岩石物理模版进行储层特征参数预测具有较高的可靠性和准确度。

图8是本申请一个实例中根据预测得到的储层特征参数数据反演得出的连井剖面。

其中，图8(a)对应的是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙度连井剖面。

图8(b)对应的是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的饱和度连井剖面。

图8(c)对应的是本申请一个实例中，反演得到的四川高石梯MX203井和MX8井的孔隙扁度连井剖面。

利用储层的孔隙扁度、饱和度和孔隙度，可以得到相应的连井剖面。

基于本申请所述的储层特征参数的预测方法，本申请提供一种储层特征参数的预测装置，所述装置可以集成在储层预测的功能组件中，进行储层预测。图4是本申请一个实施例中提供的储层特征参数的预测装置的模块结构示意图。如图4所示，所述装置可以包括：

数据获取模块101，可以用于获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据。

模版构建模块102，可以用于构建预先构建好的三维岩石物理模版。

参数预测模块103，可以用于将待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中，对所述待预测位置的储层特征参数进行预测。

数据输出模块104，可以用于输出预测得到的待预测位置的储层特征参数数据。利用上述实施例提供的一种储层特征参数的预测装置，可以自动实施所述储层特征参数的预测方法，并自动输出预测结果，操作简单快捷，不需要实施人员的过多参与，有效提高了用户体验。

在本申请另一个实施例中，所述参数预测模块的模块结构示意图如图5所示，具体的，可以包括：

分布范围确定模块1031，可以用于确定储层特征参数的分布范围。

网格生成模块1032，可以用于对所述储层特征参数的分布范围进行网格剖分，生成三维网格。

数据计算模块1033，可以用于计算所述三维网格中各节点对应的纵波阻抗、横波阻抗和密度。

模版生成模块1034，可以用于利用计算得到的三维网格节点对应的纵波阻抗、横波阻抗和密度，构建所述三维岩石物理模版。

利用上述实施例提供的参数预测模块，可以自动构建所述三维岩石物理模版，操作简单快捷，提高了用户体验。

在本申请又一个实施例中，所述数据计算模块的模块结构示意图如图9所示，具体的，可以包括：

模量计算模块10331，可以用于利用所述岩石基质的体积模量、剪切模量、所述各节点对应的孔隙扁度和孔隙度，计算得到储层的干岩石骨架的体积模量和剪切模量。

还可以用于利用所述干岩石骨架的体积模量和剪切模量，计算得到储层的流体饱和岩石的体积模量和剪切模量。

还可以用于利用所述各节点对应的饱和度，计算得到储层孔隙中流体的体积模量。

地震响应特征参数计算模块10332，可以用于利用所述流体饱和岩石的体积模量和剪切模量、所述岩石基质的密度、所述流体的密度，计算得到所述饱和流体岩石的密度、纵波速度和横波速度。还可以用于计算得到储层孔隙中流体的密度。

波阻抗计算模块10333，可以用于利用所述饱和流体岩石的纵波速度和横波速度，计算得到所述各节点对应的纵波阻抗和横波阻抗。

利用上述实施例提供的一种数据计算模块，可以直接计算出所述各节点对应的密度、纵波阻抗和横波阻抗。

本申请另一个实施例中，所述计算得到储层的干岩石骨架的体积模量和剪切模量，可以包括采用下述公式计算得到储层的干岩石骨架的体积模量和剪切模量：

式中，Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

Gm表示储层的岩石基质的剪切模量；

φ表示各节点对应的孔隙度；

α表示各节点的孔隙扁度；

a和b满足：

和表示干岩石骨架的极化因子；

K^*(y)和G^*(y)表示岩石等效体积模量和剪切模量。

本申请另一个实施例中，可以包括采用下述计算公式计算得到储层孔隙中的体积模量：

Kf＝SwKw+(1-Sw)Kg

式中，Kf表示所述储层孔隙中流体的体积模量；

Sw表示所述各节点对应的饱和度；

Kw表示储层孔隙中水的体积模量；

Kg表示储层孔隙中气体的体积模量。

本申请另一个实施例中，可以包括采用下述计算公式计算得到储层孔隙中流体的密度：

ρf＝Swρw+(1-Sw)ρg

式中，ρf表示所述储层孔隙中流体的密度；

Sw表示所述各节点对应的饱和度；

ρw表示储层孔隙中水的密度；

ρg表示储层孔隙中气体的密度。

本申请另一个实施例中，可以包括采用下述计算公式得到储层的流体饱和岩石的体积模量和剪切模量：

Gsat＝Gdry

式中，Ksat表示所述储层的饱和流体岩石的体积模量；

Gsat表示所述储层的饱和流体岩石的体积模量；

Kdry表示所述储层的干岩石骨架的体积模量；

Gdry表示所述储层的干岩石骨架的剪切模量；

Km表示储层的岩石基质的体积模量；

φ表示各节点对应的孔隙度；

Kf表示所述储层孔隙中流体的体积模量。

本申请另一个实施例中，可以包括采用下述计算公式计算得到所述饱和流体岩石的密度、纵波速度和横波速度：

ρsat＝(1-φ)ρm+φρf；

式中，ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度；

Gsat表示所述储层的饱和流体岩石的体积模量；

Ksat表示所述储层的饱和流体岩石的体积模量；

ρm表示储层的岩石基质的密度；

ρf表示所述储层孔隙中流体的密度。

本申请另一个实施例中，可以包括采用下述计算公式计算得到所述各节点对应的纵波阻抗和横波阻抗：

Ip＝ρVp；

Is＝ρVs；

ρ＝ρsat

式中，Ip表示所述各节点对应的纵波阻抗；

Is表示所述各节点对应的横波阻抗；

ρ表示所述各节点对应的密度；

ρsat表示所述饱和流体岩石的密度；

Vp表示所述饱和流体岩石的纵波速度；

Vs表示所述饱和流体岩石的横波速度。

所述储层特征参数的预测装置中，所述构建三维岩石物理模版、将所述待预测位置的三种地震响应特征参数数据输入预先构建好的三维岩石物理模版中对所述待预测位置的储层特征参数进行预测、输出预测得到的待预测位置的储层特征参数数据的实施方式的扩展可以参照前述方法的相关描述。

尽管本申请内容中提到不同的储层特征参数的预测方式，从获取储层中待预测位置的纵波阻抗、横波阻抗和密度三种地震响应特征参数数据、构建三维岩石物理模版、对所述待预测位置的储层特征参数进行预测到输出预测得到的待预测位置的储层特征参数数据的各种时序方式、数据获取/处理/输出方式等的描述，但是，本申请并不局限于必须是行业标准或实施例所描述的情况等，某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例，仍然可以属于本申请的可选实施方案范围之内。

虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。

上述实施例阐明的装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。

本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

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