基于三维弹电岩石物理量版确定储层参数的方法及装置与流程
本发明涉及石油天然气地震勘探技术领域，特别地，涉及一种基于三维弹电岩石物理量版确定储层参数的方法及装置。
背景技术
岩石物理模版是通过建立适当的岩石物理模型，模拟储层物性参数与地震响应特征之间的联系，建立储层参数和地震属性之间的定量解释关系。这种关系的确立对储层物性及其状态研究是非常重要，可以增强储层预测是结果的可靠性，提高储层流体和岩性识别的能力，降低定量解释的风险。
在岩石物理理论研究的基础上，国内外部分学者对使用岩石物理量板进行储层参数和岩性预测开展了探索和研究。目前，基于二维岩石物理量板来预测储层参数或者识别岩性比较成熟，但二维岩石物理量板在预测孔隙度、含油气饱和度和泥质含量中的两个参数时需要固定其中一个参数，这样才能得到另外两个参数的预测结果。对于泥质含量比较稳定的地区比较适用，而对于泥质含量变化比较大的地方预测的含油气饱和度和泥质含量与实际预测结果误差较大。因此，业内需要一种可以更加准确的预测储层参数的方法。
技术实现要素：
本申请实施例的目的在于提供一种基于三维弹电岩石物理量版确定储层参数的方法及装置，可以更加准确的定量预测储层参数，进一步提高用于油气检测的物性参数的精度。
本申请提供的一种基于三维弹电岩石物理量版确定储层参数的方法及装置是通过包括以下方式实现的：
一种基于三维弹电岩石物理量版确定储层参数的方法，包括：
获取目标地层的三维弹性以及电性参数实测数据；
将所述三维弹性以及电性参数实测数据向预先构建的三维弹电岩石物理量版进行投影，获得投影点，其中，三维弹电岩石物理量版的构建包括：将三维储层参数空间进行网格化剖分处理，计算网格节点对应的弹性以及电性参数数据，获得三维弹电岩石物理量版；
计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据；
根据所述投影点的储层参数数据确定所述目标地层的储层参数数据。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述将三维储层参数空间进行网格化剖分处理，包括：
根据测井数据确定目标工区的储层参数数据取值区间；
将三维储层参数空间中的所述储层参数数据取值区间基于预设精度值进行划分，获得网格化剖分处理后的三维储层参数空间。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述计算网格节点对应的弹性参数数据，包括：
基于椭球颗粒接触模型构建弹性岩石计算模型，根据所述弹性岩石计算模型计算确定网格节点对应的弹性参数数据。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述计算网格节点对应的电性参数数据，包括：
根据电性岩石计算模型计算网格节点对应的电性参数数据，其中，所述电性岩石计算模型包括改进的西门度公式：
其中，rw、rsh、rt分别表示地层水电阻率、黏土电阻率、地层电阻率，φ、sw、vsh分别表示孔隙度、含水饱和度、泥质含量，a、m和n为常数。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述计算确定网格节点对应的弹性参数数据，包括：
所述弹性参数包括纵波速度；
基于椭球颗粒接触模型计算干岩石的体积模量和剪切模量；
根据所述干岩石的体积模量和剪切模量、流体饱和条件下的岩石体积模量和剪切模量、岩石密度计算获得纵波速度数据。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述基于椭球颗粒接触模型计算干岩石的体积模量和剪切模量，包括：
其中，kdry、gdry分别表示干岩石的体积模量、和剪切模量，f1＝1-[α-0.0903-1]1.456，f2＝1-[α-0.1026-1]1.531，
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述计算网格节点对应的弹性以及电性参数数据，包括：
根据测井数据或者岩心数据对所述弹性岩石计算模型以及电性岩石计算模型进行校正；
利用校正后的述弹性岩石计算模型以及电性岩石计算模型，计算网格节点对应的弹性以及电性参数数据。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据，包括：
筛选出所述投影点与各网格节点的距离小于预设阈值的网格节点；
将距离小于预设阈值的网格节点按照距离的大小进行排序，获得距离最小值对应的网格节点。本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述获取目标地层的三维弹性以及电性参数实测数据，包括：
获取目标地层的叠前地震数据以及电磁数据；
根据所述叠前地震数据反演获得弹性参数实测数据，以及根据所述电磁数据反演获得电性参数实测数据。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的方法，所述获取目标地层的三维弹性以及电性参数实测数据，包括：
获取预设类型的弹性参数以及电性参数，根据对预设类型的弹性参数以及电性参数相对储层参数的敏感性分析结果，优选出敏感弹性参数以及敏感电性参数；
将所述敏感弹性参数以及敏感电性参数作为目标地层的三维弹性以及电性参数，获取目标地层的三维弹性以及电性参数实测数据。
另一方面，本申请实施例还提供一种基于三维弹电岩石物理量版确定储层参数的装置，包括：
获取模块，用于获取目标地层的三维弹性以及电性参数实测数据；
投影模块，用于将所述三维弹性以及电性参数实测数据向预先构建的三维弹电岩石物理量版进行投影，获得投影点，其中，三维弹电岩石物理量版的构建包括：将三维储层参数空间进行网格化剖分处理，计算网格节点对应的弹性以及电性参数数据，获得三维弹电岩石物理量版；
储层参数确定模块，用于计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据，根据所述投影点的储层参数数据确定所述目标地层的储层参数数据。
本申请实施例的基于三维弹电岩石物理量版确定储层参数的装置，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取目标地层的三维弹性以及电性参数实测数据；
将所述三维弹性以及电性参数实测数据向预先构建的三维弹电岩石物理量版进行投影，获得投影点，其中，三维弹电岩石物理量版的构建包括：将三维储层参数空间进行网格化剖分处理，计算网格节点对应的弹性以及电性参数数据，获得三维弹电岩石物理量版；
计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据，根据所述投影点的储层参数数据确定所述目标地层的储层参数数据。
本说明书一个或多个实施例提供的一种基于三维弹电岩石物理量版确定储层参数的方法及装置，可以通过预先构建三维弹电岩石物理量版，在三维空间中确定储层参数与弹-电参数之间的对应关系；然后将目标地层的三维弹性以及电性参数实测数据投影到预先构建的三维弹电岩石物理量版中，以定量确定目标地层对应的储层参数数据。利用本申请各个实施例，可以根据三维弹电岩石物理量版来准确的刻画弹-电参数与储层参数间的三维定量解释关系，用之于三维弹-电参数数据中反演获得储层参数数据，从而进一步提高用于油气检测的物性参数的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本说明书提供的一种基于三维弹电岩石物理量版确定储层参数的方法实施例的流程示意图；
图2为本说明书提供的一个实例中基于三维弹电物理量板定量预测储层参数的流程图；
图3为本说明书提供的一个实例中基于纵波速度、密度和对数电阻率构建的三维弹电物理量板预测含气饱和度的示意图；
图4为本说明书提供的一个实例中基于纵波速度、密度和对数电阻率构建的三维弹电物理量板预测孔隙度的示意图；
图5为本说明书提供的一个实例中基于纵波速度、密度和对数电阻率构建的三维弹电物理量板预测泥质含量的示意图；
图6为本说明书提供的另一个实例中测井测量的纵波阻抗、密度和对数电阻率数据示意图；
图7为本说明书提供的另一个实例中基于测井测量的纵波阻抗、密度和对数电阻率数据预测的含气饱和度、孔隙度和泥质含量数据示意图；
图8为本说明书提供的一种基于三维弹电岩石物理量版确定储层参数的装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是说明书一部分实施例，而不是全部的实施例。基于说明书一个或多个实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书实施例方案保护的范围。
针对目前二维岩石物理量版来预测储层参数或者识别岩性时，只能利用二维测量数据来定量预测二维储层参数，无法准确的完成多维储层参数的预测工作的问题。本申请提供了一种三维弹电岩石物理量版来定量刻画含油气储层。例如，可以通过构建三维弹电岩石物理量版建立一种能准确刻画弹性参数、电性参数与储层参数(如孔隙度-饱和度-泥质含量)间的三维定量解释关系，用之于三维弹-电参数数据中反演获得储层参数数据。从而，利用本申请提供的方案可以进一步提高用于油气检测的物性参数的精度。
本申请各个实施中，弹性参数、电性参数以及待预测的储层参数的类型可以根据实际情况确定。例如，所述储层参数可以包括储层孔隙度、含气饱和度、泥质含量等，所述弹性参数可以包括纵波(或者横波)速度、密度、纵波(横波)阻抗、泊松比、纵横波速度比等，所述电性参数可以包括电阻率、对数电阻率等。
图1是本说明书提供的所述一种基于三维弹电岩石物理量版确定储层参数的方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
具体的一个实施例如图1所示，本说明书提供的基于三维弹电岩石量版确定储层参数的方法的一个实施例中，所述方法可以包括：
s2、获取目标地层的三维弹性以及电性参数实测数据。
本实施例中，可以针对目标工区的不同目标地层，获取相应的三维弹性以及电性参数实测数据。其中，弹性参数可以包括纵波(或者横波)速度、岩石密度、纵波(横波)阻抗、泊松比、纵横波速度比等，所述电性参数可以包括地层的电阻率或者对数电阻率等。具体实施时，可以根据预先构建的三维弹电岩石物理量版中的弹性参数以及电性参数类型确定需要获取的三维弹性以及电性参数实测数据。
例如，预先构建的三维弹电岩石物理量版中的弹性参数以及电性参数类型分别为：纵波速度、岩石密度以及对数电阻率，相应的，可以获取目标地层的纵波速度、岩石密度以及对数电阻率实测数据，组成目标地层的三维弹性以及电性参数实测数据。本说明书的一个实施例中，可以获取目标工区的叠前地震数据以及电磁数据，然后，利用叠前地震数据反演获得弹性参数数据，以及利用电磁数据反演获得电性参数数据，从而提高获得的弹性以及电性参数数据的准确性。
s4、将所述三维弹性以及电性参数实测数据向预先构建的三维弹电岩石物理量版进行投影，获得投影点。
本实施例中，可以将上述步骤中获取的弹性以及电性参数实测数据构成的三维数据，向预先构建的三维弹电岩石物理量版进行投影，获得该三维数据对应的投影点。例如，假设上述三维数据为(vp,ρ,rt)，其中，vp表示纵波速度，ρ表示岩石密度，rt表示电阻率。则可以将(vp,ρ,rt)向三维弹电岩石物理量版进行投影，获得(vp,ρ,rt)在的三维弹电岩石物理量版空间中的投影点位置。
本说明书的一个实施例中，所述三维弹电岩石物理量版的构建可以包括：将三维储层参数空间进行网格化剖分处理，然后计算各网格节点对应的弹性以及电性参数数据，从而将三维储层参数空间中各网格节点对应的储层参数数据与三维弹电参数空间中的弹-电参数数据一一关联起来，获得三维弹电岩石物理量版。
进一步的，本说明书的一个实施例中，可以首先确定目标工区的储层参数数据变化范围，根据所述储层参数数据变化范围确定各储层参数数据取值区间；然后，将三维储层参数空间中的所述储层参数数据取值区间基于预设精度值进行划分，从而将三维储层参数空间进行网格化剖分，获得网格化剖分处理后的三维储层参数空间。具体实施时，所述预设精度值可以在兼顾精度与效率的情况下根据实际情况进行预先设定。所述目标工区的储层参数数据变化范围可以根据如测井数据等实测数据进行确定。然后，计算网格化剖分处理后的三维储层参数空间中各网格节点对应的弹性以及电性参数数据，从而构建三维弹电岩石物理量版。
本实施例中，通过首先根据目标工区的实测数据确定各储层参数的数据变化范围，然后，在该范围内对三维储层参数空间进行网格化剖分处理，可以提高网格化处理的效率；同时，还可以减少网格节点数，降低网格节点对应的弹-电参数的数据计算量，大幅度提高数据处理的效率。进一步的，还可以利用实测数据有效限定数据筛选范围，降低筛选误差，提高后续基于弹电岩石物理量版确定目标工区的储层参数数据的准确性。
例如，本说明书的一个或者多个实施例中，假设所述储层参数类型包括：孔隙度、含水饱和度以及泥质含量。则可以首先根据岩心和测井解释获得的数据确定目标工区的孔隙度、含水饱和度以及泥质含量三个储层参数的变化范围，并根据所述储层参数数据变化范围确定各储层参数数据取值区间。然后，将三维储层参数空间中的所述储层参数数据取值区间基于相应的预设精度值进行划分。相应的，孔隙度对应数值相等的点连成线即为孔隙度对应的等势线，同理，含水饱和度以及泥质含量对应数值相等的点连成线即为含水饱和度以及泥质含量对应的等势线。各参数等势线的交点为网格节点，从而形成由孔隙度、含水饱和度以及泥质含量组成的三维网格空间，从而获得网格化剖分处理后的三维储层参数空间。
然后，可以计算网格化后的三维储层参数空间中各网格节点对应的弹性以及电性参数数据。本说明书的一些实施方式中，可以根据预先构建的弹性岩石计算模型以及电性岩石计算模型，计算各网格节点对应的弹性以及电性参数数据。本说明书的一个实施例中，可以基于椭球颗粒接触模型构建弹性岩石计算模型，根据所述弹性岩石计算模型计算确定弹性参数数据。其中，所述椭球颗粒接触模型包括：可以将地层岩石等效为椭球状颗粒，即地层孔隙纵横比不等于1。将地层岩石等效为椭球状颗粒更符合地层实际情况，从而基于椭球颗粒接触模型可以进一步提高弹性参数数据计算的准确性。
本说明书的一个或者多个实施例中，假设三维弹电岩石物理量版对应的弹性以及电性参数分别为vp,ρ,rt，待预测的储层参数为孔隙度φ、含水饱和度sw和泥质含量vsh。则在上述由孔隙度、含水饱和度以及泥质含量组成的三维网格空间内，对于各网格节点的φ、sw、vsh的三维数据，可以由岩石矿物成分及其含量利用沃伊特-罗伊斯-伊尔(voigt-reuss-hill)公式平均计算岩石基质的体积模量和剪切模量。然后，可以基于椭球颗粒接触模型计算干岩石的体积模量kdry和剪切模量gdry，其中，对应的计算公式可以表达为：
式中，f1＝1-[α-0.0903-1]1.456，f2＝1-[α-0.1026-1]1.531，
然后，可以采用伍德(wood)公式计算混合流体的等效体积模量：
式中，kfl为流体的体积模量，kg为油气的体积模量，kw为水的体积模量。
之后，可以采用加斯曼(gassmann)公式计算流体饱和情况下的岩石体积模量和剪切模量，公式如下：
gsat＝gdry(5)
式中，ksat和gsat分别为流体饱和的体积模量和剪切模量。
之后，可以由岩石物理体积模型计算出岩石的密度，然后再计算出纵波速度：
式中，ρ为饱和岩石的密度。
本说明书的另一个实施例中，可以采用改进的西门度公式计算电阻率:：
式中，rw、rsh和rt分别为地层水电阻率、黏土电阻率和地层电阻率，a、m和n均为阿尔奇公式的常数。
重复上述过程，可以计算获得所有网格节点的弹性参数数据和电性参数数据，将储层参数对应的三维空间网格与弹-电参数对应的三维空间一一关联起来，形成三维弹电岩石物理量版。
本说明书的一个实施例中，还可以根据测井数据或者岩心数据对所述弹性岩石计算模型以及电性岩石计算模型进行校正，具体实施时，可以对弹性岩石计算模型以及电性岩石计算模型中的常数数据进行校正。然后，利用校正后的弹性岩石计算模型以及电性岩石计算模型，计算各网格节点对应的弹性以及电性参数数据，从而进一步提高三维弹电岩石物理量版构建的精确度。
s6、计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据，根据所述投影点的储层参数数据确定所述目标地层的储层参数数据。
本实施例中，对于s4中获得的投影点，可以在预先确定的储层参数数据变化范围内，计算投影点与三维弹电岩石物理量版中各网格节点之间的距离(例如欧式距离)。然后，筛选出投影点与各网格节点的距离最小值对应的网格节点，获取所述距离最小值对应的网格节点的储层参数数据(如：孔隙度、含水饱和度和泥质含量)作为该投影点的储层参数数据。而投影点与各目标地层的实测弹-电参数数据一一对应，则确定投影点的储层参数数据也即确定了各目标地层的储层参数数据。从而利用上述步骤定量的预测出了目标工区中各地层的储层参数数据。
本说明书的一个实施例中，可以首先筛选出投影点与各网格节点的距离小于预设阈值的网格节点，然后将距离小于预设阈值的网格节点按照距离的大小进行排序，从而确定出距离最小值对应的网格节点。之后，再获取所述距离最小值对应的网格节点的储层参数数据作为该投影点的储层参数数据。其中，所述预设阈值的大小可以根据实际情况进行调整，例如，如果筛选出的网格节点较多，则可以缩小预设阈值的大小；反之，如果没有筛选出网格节点，则可以适当的扩大预设阈值的大小。利用本实施例上述方案，可以进一步减少运算量，提高数据处理效率。
重复上述步骤，可以定量的确定目标工区各地层对应的储层参数数据。同时，本说明书各个实施例还通过引入电性参数(如电阻率)信息弥补了弹性参数对饱和度不敏感的问题，从而进一步提高含油气饱和度的预测精度。
本说明书的一些实施方式中，还可以在构建三维弹电岩石物理量版之前，根据不同的目标工区，优选出相应的敏感弹性参数以及敏感电性参数，或者相适应于该目标工区的待预测储层参数。然后，利用所述敏感弹(电)性参数以及带预测储层参数构建更适用于该目标工区的三维弹电岩石物理量版，以使得构建的三维弹电岩石物理量版更加适应不同地区的复杂地形状况。
本说明书的一个实施例中，可以获取预设类型的弹性参数以及电性参数，然后对预设类型的弹性参数以及电性参数进行相对待预测储层参数的敏感性分析，根据分析结果优选出敏感弹性或者电性参数。将优选出的敏感弹性或者电性参数用于构建三维弹电岩石物理量版，以及获取目标地层的所述敏感弹性或者电性参数的实测数据，用于确定目标地层的储层参数数据。例如，可以对比分析多个类型的弹(电)性参数数据变化时，相应的待预测储层参数的变化程度。将待预测储层参数变化相对明显时对应的弹(电)性参数类型作为敏感弹(电)性参数。利用本实施例的方案，可以提高建立的弹电岩石物理量版对目标地层的适应性，从而进一步提高储层参数预测的精确性。
当然，本说明书的一些实施方式中，也可以根据目标地区的实际情况，选择可以较易获得的弹性参数或者电性参数构建三维弹电岩石物理量版。然后，针对三维弹电岩石物理量版中的弹性参数以及电性参数，获取目标地层的三维弹性以及电性参数实测数据，进一步确定所述目标地层的储层参数。
从而使得上述实施例中的方案，不仅可以应用于常规的泥沙岩，也可以通过选用合适的弹性、电性参数以及储层参数构建不同的三维弹电岩石物理量版，以满足复杂储层参数预测的应用环境。并且，本说明书各个实施例提供的方案还可以应用于多个尺度(岩心、测井和地震-电磁)的数据预测储层参数，从而为定量预测储层提供了一种非常简单实用的方法。
为了进一步验证本方法的实用性和可行性，本说明书还提供了应用上述方案的一个具体实例，图2表示利用三维弹电岩石物理量版定量预测储层参数的流程图，根据目标工区的实际情况，本实例中构建三维弹电岩石物理量版采用的弹性参数以及电性参数为：纵波速度、岩石密度以及对数电阻率，待预测的储层参数为：含气饱和度、孔隙度以及泥质含量。如图2所示，具体实施方案可以如下：
(1)获取弹性参数以及电性参数实测数据。
获取目标工区的叠前地震数据、电磁数据，并反演获得弹性参数(纵波速度vp、岩石密度ρ)以及电性参数(对数电阻率log10(rt))实测数据。
(2)将实测数据向预先构建的三维弹电岩石物理量版进行投影。
图3、图4、图5分别表示采用预先构建的三维弹电岩石物理量版预测孔隙度φ、含水饱和度sw以及泥质含量vsh的示意图，将纵波速度(p-wavevelocity)、岩石密度(density)以及对数电阻率(log10(rt))组成的三维实测数据向三维弹电岩石物理量版中投影，如图3-图5所示，图中的点表示相应的投影点。其中，纵波速度的单位为：千米/秒(km/s)，岩石密度的单位为：克/立方厘米(g/cm3)，对数电阻率的单位为：欧姆米(ohm-m)。
相应的，可以根据如下步骤预先构建三维弹电岩石物理量版：
首先，可以根据岩心和测井解释成果确定孔隙度、含水饱和度以及泥质含量的变化范围，将相应数据变化范围内的孔隙度、含水饱和度以及泥质含量按照预设的要求进行等分。从而将三维储层参数空间进行网格化剖分，形成由孔隙度、含水饱和度以及泥质含量组成的三维空间网格。
如图3-图5所示，孔隙度(porosity)的变化范围为0.1-0.3，以0.1为步长进行等分；含气饱和度(gassaturation)的变化范围为0.1-0.6，以0.1为步长进行等分；泥质含量(claycontent)的变化范围为0-0.8，以0.4为步长进行等分。其中，孔隙度为0.1、0.2、0.3的线即为孔隙度对应的等势线，含气饱和度为0.1、0.2、0.3、0.4、0.5、0.6的线即为含气饱和度对应的等势线，泥质含量为0、0.4、0.8的线即为泥质含量对应的等势线。孔隙度、含气饱和度以及泥质含量的各条等势线形成的交点即为网格节点。从而形成由孔隙度、含水饱和度以及泥质含量组成的三维空间网格。
其次，在三维空间网格内，对于各网格节点对应的孔隙度、含水饱和度和泥质含量，基于椭球颗粒接触模型，可以利用公式(1)-(6)计算确定纵波速度，以及根据岩石物理体积模型计算确定岩石密度；并利用改进的西门度公式(7)计算确定电阻率，从而获得对数电阻率数据。从而，对于每一网格节点对应的储层参数数据(φ,sw,vsh)，可以利用预先构建的弹性以及电性岩石计算模型计算确定出相应的三维弹-电参数数据(vp,ρ,log10(rt))。从而将储层参数对应的三维空间网格与弹-电参数对应的三维空间一一关联起来，获得三维弹电岩石物理量版。
(3)确定目标地层的储层参数。
对获得的投影点，可以在预先确定的储层参数数据变化范围内，计算投影点与三维弹电岩石物理量版中各网格节点之间的距离。并计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据。投影点与各目标地层的实测弹-电参数数据一一对应，则确定投影点的储层参数数据也即确定了各目标地层的储层参数数据。从而利用上述步骤定量的预测出了目标工区中各地层的储层参数数据。
图6和图7分别表示基于三维弹电岩石物理量版定量预测某未固结储层测井数据的实例示意图。其中，弹性参数为：纵波阻抗(p-waveimpedance)、岩石密度(density)，电性参数为：对数电阻率(log10(rt))，其中，纵波阻抗的单位为：千米/秒*克/立方厘米(km/s*g/cm3)。待预测的储层参数为：含气饱和度(gassaturation)、孔隙度(porosity)以及泥质含量(claycontent)。
图6(a)、图6(b)、图6(c)分别表示实测以及反演获得的纵波阻抗、岩石密度以及对数电阻率随地层深度(depth)，其中实线表示实测的纵波阻抗、岩石密度以及对数电阻率数据，虚线表示根据上述弹性岩石计算模型以及电性岩石计算模型反演获得的纵波阻抗、岩石密度以及对数电阻率数据。图7(a)、图7(b)、图7(c)分别表示预测的含气饱和度、孔隙度以及泥质含量的数据。图7(a)、图7(b)、图7(c)中的圆点表示测井实际测得的含气饱和度、孔隙度以及泥质含量数据，实线表示利用三维弹电岩石物理量版预测的含气饱和度、孔隙度以及泥质含量数据。其中，地层深度的单位为：海底向下深度、单位为米(mbsf)。
分析图6(a)、图6(b)、图6(c)可知，仅图6(b)中的虚线与实线没有完全重合。即利用本说明书实施例校正后的计算模型反演获得的纵波阻抗和对数电阻率数据与实测数据基本吻合，仅岩石密度实测与反演数据存在一些偏差。从而表明基于校正后的弹性或者电性岩石计算模型反演获得的电-弹性数据精确度较高，进一步提高了后续储层参数数据预测的准确度。
分析图7(a)、图7(b)、图7(c)可知，基于图6(a)、图6(b)、图6(c)反演的弹-电参数数据预测出的孔隙度、含气饱和度和泥质含量与测井解释结果吻合较好。从而表明，利用本说明书上述一个或者多个实施例提供的方案，可以大大提高储层参数(特别是含水饱和度)预测精度，而且也可以较好地完成多源数据和多维储层参数的预测工作。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书一个或多个实施例提供的一种基于三维弹电岩石物理量版确定储层参数的方法，可以通过预先构建三维弹电岩石物理量版，在三维空间中确定储层参数与弹-电参数之间的对应关系；然后将目标地层的三维弹性以及电性参数实测数据投影到预先构建的三维弹电岩石物理量版中，以定量确定目标地层对应的储层参数数据。利用本申请各个实施例，可以根据三维弹电岩石物理量版来准确的刻画弹-电参数与储层参数间的三维定量解释关系，用之于三维弹-电参数数据中反演获得储层参数数据，从而进一步提高用于油气检测的物性参数的精度。
基于上述所述的基于三维弹电岩石物理量版确定储层参数的方法，本说明书一个或多个实施例还提供一种基于三维弹电岩石物理量版确定储层参数的装置。所述的装置可以包括使用了本说明书实施例所述方法的系统、软件(应用)、模块、组件、服务器等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。具体的，图8说明书提供的一种基于三维弹电岩石物理量版确定储层参数的装置实施例的模块结构示意图，如图8，所述装置可以包括：
获取模块102，可以用于获取目标地层的三维弹性以及电性参数实测数据；
投影模块104，可以用于将所述三维弹性以及电性参数实测数据向预先构建的三维弹电岩石物理量版进行投影，获得投影点，其中，三维弹电岩石物理量版的构建包括：将三维储层参数空间进行网格化剖分处理，计算网格节点对应的弹性以及电性参数数据，获得三维弹电岩石物理量版；
储层参数确定模块106，可以用于计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据；根据所述投影点的储层参数数据确定所述目标地层的储层参数数据。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书一个或多个实施例提供的一种基于三维弹电岩石物理量版确定储层参数的装置，可以通过预先构建三维弹电岩石物理量版，在三维空间中确定储层参数与弹-电参数之间的对应关系；然后将目标地层的三维弹性以及电性参数实测数据投影到预先构建的三维弹电岩石物理量版中，以定量确定目标地层对应的储层参数数据。利用本申请各个实施例，可以根据三维弹电岩石物理量版来准确的刻画弹-电参数与储层参数间的三维定量解释关系，用之于三维弹-电参数数据中反演获得储层参数数据，从而进一步提高用于油气检测的物性参数的精度。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种基于三维弹电岩石物理量版确定储层参数的装置，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
获取目标地层的三维弹性以及电性参数实测数据；
将所述三维弹性以及电性参数实测数据向预先构建的三维弹电岩石物理量版进行投影，获得投影点，其中，三维弹电岩石物理量版的构建包括：将三维储层参数空间进行网格化剖分处理，计算网格节点对应的弹性以及电性参数数据，获得三维弹电岩石物理量版；
计算所述投影点与各网格节点的距离，根据距离最小值对应的网格节点的储层参数数据确定所述投影点的储层参数数据；根据所述投影点的储层参数数据确定所述目标地层的储层参数数据。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如ram、rom等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、u盘；利用光学方式存储信息的装置如，cd或dvd。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种基于三维弹电岩石物理量版确定储层参数的装置，可以通过预先构建三维弹电岩石物理量版，在三维空间中确定储层参数与弹-电参数之间的对应关系；然后将目标地层的三维弹性以及电性参数实测数据投影到预先构建的三维弹电岩石物理量版中，以定量确定目标地层对应的储层参数数据。利用本申请各个实施例，可以根据三维弹电岩石物理量版来准确的刻画弹-电参数与储层参数间的三维定量解释关系，用之于三维弹-电参数数据中反演获得储层参数数据，从而进一步提高用于油气检测的物性参数的精度。
需要说明的是，本说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类、存储介质+程序实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、平板计算机或者这些设备中的任何设备的组合。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述并不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
- 还没有人留言评论。精彩留言会获得点赞！