含浊沸石砂砾岩预测方法、系统、装置及存储介质

摘要

本发明提供了一种含浊沸石砂砾岩预测方法、系统、装置及存储介质。该方法包括：利用工区薄片岩性分析数据及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；利用三维地震数据体反演得到电阻率反演数据体；利用三维地震数据体反演得到密度反演数据体；利用电阻率门槛值约束电阻率反演数据体得到含浊沸石砂砾岩的有效电阻率反演数据体；利用密度门槛值约束密度反演数据体得到含浊沸石砂砾岩的有效密度反演数据体；将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体；将融合数据体利用地震解释的目标层位顶底界进行约束，得到含浊沸石砂砾岩的平面预测分布图。

说明书

技术领域

本发明涉及石油勘探技术领域，特别涉及一种含浊沸石砂砾岩储层预测方法、系统、装置及存储介质。

背景技术

砂砾岩储层是指以砾岩、砂砾岩为主的油气藏储集层，在国内以准噶尔盆地西北缘的玛湖凹陷二三叠系为典型代表。砂砾岩储层具有物性横向变化快，低孔低渗等特点。而含浊沸石砂砾岩指粒间填隙物以普遍含浊沸石为特征的一类砂砾岩储层。它具有砂砾岩储层的基本特征，又有其本身的特殊性。其本身的特殊性是识别和预测含浊沸石砂砾岩的关键。

目前，学者的研究主要集中在利用测井曲线对砂砾岩储层单井定量识别(例如，CN108303752A公开的砂砾岩有效储层常规测井定量识别方法)、利用波阻抗反演及地震属性等方法来对砂砾岩储层进行平面预测(例如《中拐凸起二叠系致密砂砾岩储层的成因及有利区块预测》，万敏，2011年，新疆石油天然气)、通过单井储层样品的X衍射和薄片镜下观察结果来预测浊沸石纵向分布规律(《吴起-志丹地区延长组下组合浊沸石的纵向分布特征与成因》，王朝等，2019年，吉林大学学报(地球科学版))。这些研究只是关注了砂砾岩储层测井识别方法、单因素地震平面预测技术及储层粒间浊沸石矿物的分布规律等方面的研究，针对含浊沸石砂砾岩储层的分布预测没有涉及和形成系统的预测方法。

发明内容

针对现有技术的不足，本发明的目的在于提供一种快速有效的地质-测井-地震一体化的含浊沸石砂砾岩储层预测方法。

为了实现上述目的，本发明提供了一种含浊沸石砂砾岩预测方法，其中，该方法包括：

利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；

利用三维地震数据体反演得到电阻率反演数据体；

利用三维地震数据体反演得到密度反演数据体；

利用确定的电阻率门槛值对电阻率反演数据体进行约束得到含浊沸石砂砾岩的有效电阻率反演数据体；

利用确定的密度门槛值对密度反演数据体进行约束得到含浊沸石砂砾岩的有效密度反演数据体；

将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体；

将得到的融合数据体，利用地震解释的目标层位顶底界进行约束，得到含浊沸石砂砾岩的平面预测分布图；

完成含浊沸石砂砾岩预测。

在上述含浊沸石砂砾岩预测方法中，优选地，所述利用工区薄片分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值包括：

获取工区薄片的岩性分析数据确定薄片的岩性；所述岩性包括含浊沸石砂砾岩以及非含浊沸石砂砾岩；

获取工区薄片对应的电阻率测井曲线、密度测井曲线；

基于确定的岩性以及获取的电阻率测井曲线、密度测井曲线，建立工区薄片含浊沸石砂砾岩以及非含浊沸石砂砾岩的电阻率-密度交会图版；

基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

在上述含浊沸石砂砾岩预测方法中，优选地，所述利用工区薄片分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值包括：

获取工区薄片的岩性分析数据确定岩性为含浊沸石砂砾岩的薄片；

获取确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线；

根据确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线，建立含浊沸石砂砾岩的电阻率-密度交会图版；

基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

在上述含浊沸石砂砾岩预测方法中，优选地，该方法进一步包括：利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的测井曲线响应特征。

在上述含浊沸石砂砾岩预测方法中，优选地，利用确定的电阻率门槛值对电阻率反演数据体进行约束得到含浊沸石砂砾岩的有效电阻率反演数据体可以通过下述方式实现：利用确定的电阻率门槛值对电阻率反演数据体进行约束：在电阻率反演数据体中，电阻率值分布在电阻率门槛值之外的电阻率反演数据体赋值为零，电阻率值分布在电阻率门槛值之内的电阻率反演数据体不为零，从而得到经电阻率门槛值约束后的有效电阻率反演数据体，以此作为含浊沸石砂砾岩的有效电阻率反演数据体。

在上述含浊沸石砂砾岩预测方法中，优选地，利用确定的密度门槛值对密度反演数据体进行约束得到含浊沸石砂砾岩的有效密度反演数据体可以通过下述方式实现：利用确定的密度门槛值对密度反演数据体进行约束：在密度反演数据体中，密度值分布在密度门槛值之外的密度反演数据体赋值为零，密度值分布在密度门槛值之内的密度反演数据体不为零，从而得到经密度门槛值约束后的有效密度反演数据体，以此作为含浊沸石砂砾岩的有效密度反演数据体。

在上述含浊沸石砂砾岩预测方法中，优选地，将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体可以通过下述方式实现：任意选定密度反演数据体和电阻率反演数据体中的一种作为基准有效反演数据体，并将另一种有效反演数据体中的全部非零数据赋值为1得到第一有效反演数据体；然后将基准有效反演数据体与第一有效反演数据体相乘，得到能预测含浊沸石砂砾岩的融合数据体。

本发明还提供了一种含浊沸石砂砾岩预测系统，其中，该系统包括：

门槛值确定模块：用于利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；

电阻率反演数据体建立模块：用于利用三维地震数据体反演得到电阻率反演数据体；

密度反演数据体建立模块：用于利用三维地震数据体反演得到密度反演数据体；

有效电阻率反演数据体建立模块：用于利用确定的电阻率门槛值对电阻率反演数据体进行约束得到含浊沸石砂砾岩的有效电阻率反演数据体；

有效密度反演数据体建立模块：用于利用确定的密度门槛值对密度反演数据体进行约束得到含浊沸石砂砾岩的有效密度反演数据体；

融合数据体建立模块：用于将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体；

浊沸石砂砾岩预测模块：用于将得到的融合数据体，利用地震解释的目标层位顶底界进行约束，得到含浊沸石砂砾岩的平面预测分布图。

在上述含浊沸石砂砾岩预测系统中，优选地，所述门槛值确定模块包括：

第一岩性获取子模块：用于获取工区薄片的岩性分析数据确定薄片的岩性；所述岩性包括含浊沸石砂砾岩以及非含浊沸石砂砾岩；

第一测井曲线获取子模块：用于获取工区薄片对应的电阻率测井曲线、密度测井曲线；

第一交会图版建立模块：用于基于确定的岩性以及获取的电阻率测井曲线、密度测井曲线，建立工区薄片含浊沸石砂砾岩以及非含浊沸石砂砾岩的电阻率-密度交会图版；

第一电阻率门槛值和密度门槛值确定子模块：用于基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

在上述含浊沸石砂砾岩预测系统中，优选地，所述门槛值确定模块包括：

第二岩性获取子模块：用于获取工区薄片的岩性分析数据确定岩性为含浊沸石砂砾岩的薄片；

第二测井曲线获取子模块：用于获取确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线；

第二交会图版建立模块：用于根据确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线，建立含浊沸石砂砾岩的电阻率-密度交会图版；

第二电阻率门槛值和密度门槛值确定子模块：用于基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

在上述含浊沸石砂砾岩预测系统中，优选地，该系统进一步包括：测井曲线响应特征确定模块：用于利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的测井曲线响应特征。

在上述含浊沸石砂砾岩预测系统中，优选地，有效电阻率反演数据体建立模块用于在电阻率反演数据体中，将电阻率值分布在电阻率门槛值之外的电阻率反演数据体赋值为零，电阻率值分布在电阻率门槛值之内的电阻率反演数据体不为零，从而得到经电阻率门槛值约束后的有效电阻率反演数据体，以此作为含浊沸石砂砾岩的有效电阻率反演数据体。

在上述含浊沸石砂砾岩预测系统中，优选地，有效电阻率反演数据体建立模块用于在密度反演数据体中，将密度值分布在密度门槛值之外的密度反演数据体赋值为零，密度值分布在密度门槛值之内的密度反演数据体不为零，从而得到经密度门槛值约束后的有效密度反演数据体，以此作为含浊沸石砂砾岩的有效密度反演数据体。

在上述含浊沸石砂砾岩预测系统中，优选地，融合数据体建立模块包括：

第一赋值子模块：用于任意选定密度反演数据体和电阻率反演数据体中的一种作为基准有效反演数据体，并将另一种有效反演数据体中的全部非零数据赋值为1得到第一有效反演数据体；

第一融合子模块：用于将基准有效反演数据体与第一有效反演数据体相乘，得到能预测含浊沸石砂砾岩的融合数据体。

本发明还提供了一种含浊沸石砂砾岩预测装置，包括处理器及存储器；其中，

存储器，用于存放计算机程序；

处理器，用于执行存储器上所存放的程序时，实现上述含浊沸石砂砾岩预测方法的步骤。

本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述含浊沸石砂砾岩预测方法的步骤。

本发明提供的上述技术方案，基于薄片与电阻率、密度测井曲线相结合，明确含浊沸石砂砾岩的电阻率门槛值、密度门槛值；利用该门槛值分别对电阻率反演数据体和密度反演数据体进行约束，将约束后的有效电阻率反演数据体与有效密度反演数据体进行融合，最终实现含浊沸石砂砾岩的平面预测。该技术方案将地质、测井与地震有效结合，能够广泛应用于含浊沸石砂砾岩预测，实现对含浊沸石砂砾岩进行快速有效的预测。

附图说明

图1为本发明一实施例提供的含浊沸石砂砾岩预测方法的流程示意图。

图2为本发明一实施例提供的含浊沸石砂砾岩预测方法中确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值步骤的优化示意图。

图3为本发明一实施例提供的含浊沸石砂砾岩预测方法中确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值步骤的优化示意图。

图4为本发明一实施例提供的含浊沸石砂砾岩预测方法的流程示意图。

图5为本发明一实施例中建立的电阻率-密度交会图版。

图6A为本发明一实施例中电阻率反演数据体图。

图6B为本发明一实施例中密度反演数据体图。

图6C为本发明一实施例中有效电阻率反演数据体图。

图6D为本发明一实施例中有效密度反演数据体图。

图6E为本发明一实施例中融合数据体图。

图7为本发明一实施例中含浊沸石砂砾岩的平面预测分布图。

图8为本发明一实施例提供的含浊沸石砂砾岩预测系统的结构示意图。

图9为本发明一实施例提供的含浊沸石砂砾岩预测系统中门槛值确定模块的结构示意图。

图10为本发明一实施例提供的含浊沸石砂砾岩预测系统中门槛值确定模块的结构示意图。

图11为本发明一实施例提供的含浊沸石砂砾岩预测系统中融合数据体建立模块的结构示意图。

图12为本发明一实施例提供的含浊沸石砂砾岩预测装置的结构示意图。

具体实施方式

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

下面参考本发明的若干代表性实施方式，详细阐述本发明的原理和精神。

参见图1，为了实现上述目的，本发明提供了一种含浊沸石砂砾岩预测方法，其中，该方法包括：

步骤S1：利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；

步骤S2：利用三维地震数据体反演得到电阻率反演数据体；

步骤S3：利用三维地震数据体反演得到密度反演数据体；

步骤S4：利用确定的电阻率门槛值对电阻率反演数据体进行约束得到含浊沸石砂砾岩的有效电阻率反演数据体；

步骤S5：利用确定的密度门槛值对密度反演数据体进行约束得到含浊沸石砂砾岩的有效密度反演数据体；

步骤S6：将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体；

步骤S7：得到的融合数据体，利用地震解释的目标层位顶底界进行约束，得到含浊沸石砂砾岩的平面预测分布图。

参见图2，在一个实施方式中，步骤S1可以包括：

步骤S11：获取工区薄片的岩性分析数据确定各薄片的岩性；所述岩性包括含浊沸石砂砾岩以及非含浊沸石砂砾岩；

步骤S12：获取工区薄片对应的电阻率测井曲线、密度测井曲线；

步骤S13：基于确定的岩性以及获取的电阻率测井曲线、密度测井曲线，建立工区薄片含浊沸石砂砾岩以及非含浊沸石砂砾岩的电阻率-密度交会图版；

步骤S14：基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

参见图3，在一个实施方式中，步骤S1可以包括：

步骤S111：获取工区薄片的岩性分析数据确定岩性为含浊沸石砂砾岩的薄片；

步骤S121：获取确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线；

步骤S131：根据确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线，建立含浊沸石砂砾岩的电阻率-密度交会图版；

步骤S141：基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

其中，工区薄片的岩性分析数据采用常规方法获取即可；例如，采用岩石薄片镜进行常规的薄片岩性观察或者薄片岩性鉴定等等。

其中，获取薄片对应的电阻率测井曲线、密度测井曲线采用常规方法即可；例如，可以确定薄片对应的深度，获取对应深度的电阻率测井曲线、密度测井曲线。

其中，建立电阻率-密度交会图版采用常规方法即可；例如根据各个薄片对应的电阻率响应值、密度响应值，建立以电阻率、密度为坐标轴的电阻率-密度交会图版即可，其中电阻率响应值、密度响应值通过电阻率测井曲线、密度测井曲线即可获知。

在一个实施方式中，该方法进一步可以包括：利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的测井曲线响应特征。

在一个实施方式中，步骤S2利用三维地震数据体反演得到电阻率反演数据体可以采用地质统计学反演方法，如地震波形指示法。在一示例中，利用单井电阻率测井曲线作为敏感参数曲线，建立样本点，然后利用多井样本拟合建立概率密度函数，其次进行空间变差函数分析，确定空间变异程度，在三维地震数据的约束下，对井间储层进行统计评估，得出电阻率反演数据体。

在一个实施例中，步骤S3利用三维地震数据体反演得到密度反演数据体可以采用地质统计学反演方法，如地震波形指示法。在一示例中，利用单井密度测井曲线作为敏感参数曲线，建立样本点，然后利用多井样本拟合建立概率密度函数，其次进行空间变差函数分析，确定空间变异程度，在三维地震数据的约束下，对井间储层进行统计评估，得出密度反演数据体。

在一个实施方式中，步骤S4利用确定的电阻率门槛值对电阻率反演数据体进行约束得到含浊沸石砂砾岩的有效电阻率反演数据体可以通过下述方式实现：利用确定的电阻率门槛值对电阻率反演数据体进行约束：在电阻率反演数据体中，电阻率值分布在电阻率门槛值之外的电阻率反演数据体赋值为零，电阻率值分布在电阻率门槛值之内的电阻率反演数据体不为零，从而得到经电阻率门槛值约束后的有效电阻率反演数据体，以此作为含浊沸石砂砾岩的有效电阻率反演数据体。其中，电阻率分布在电阻率门槛值之内的电阻率反演数据体保留(即保持不变)或者修改为其他任意不为零的数均可。

在一个实施方式中，步骤S5利用确定的密度门槛值对密度反演数据体进行约束得到含浊沸石砂砾岩的有效密度反演数据体可以通过下述方式实现：利用确定的密度门槛值对密度反演数据体进行约束：在密度反演数据体中，密度值分布在密度门槛值之外的密度反演数据体赋值为零，密度分布在密度门槛值之内的密度反演数据体不为零，从而得到经密度门槛值约束后的有效密度反演数据体，以此作为含浊沸石砂砾岩的有效密度反演数据体。其中，密度分布在电阻率门槛值之内的密度反演数据体保留(即保持不变)或者修改为其他任意不为零的数均可。

在一个实施方式中，步骤S6将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体可以通过下述方式实现：任意选定密度反演数据体和电阻率反演数据体中的一种作为基准有效反演数据体，并将另一种有效反演数据体中的全部非零数据赋值为1得到第一有效反演数据体；然后将基准有效反演数据体与第一有效反演数据体相乘，得到能预测含浊沸石砂砾岩的融合数据体。利用该融合数据体能够有效识别哪些为含浊沸石砂砾岩，具体而言：在同一个线道号，同一时间深度点上，如果基准有效反演数据体与第一有效反演数据体都有非零值的情况下，最终得到的值也非零，也就是融合后的表征含浊沸石砂砾岩的有效电阻率-密度数据；而在同一个线道号，同一时间深度点上，只要有一个值为零，则最终得到的值也为零，也就是融合后的表征非含浊沸石砂砾岩的非有效电阻率-密度数据，基于此，该融合数据体能预测含浊沸石砂砾岩。

在一个实施方式中，步骤S7得到的融合数据体，利用地震解释的目标层位顶底界进行约束，得到含浊沸石砂砾岩的平面预测分布图可以采用常规方式进行；例如，将得到的含浊沸石砂砾岩储层融合数据体，利用地震解释的目标层位顶底界进行约束，可以运用均方根、平均值、最大值、求和等数学运算方法，计算得出含浊沸石砂砾岩的平面预测分布图。

参见图4，本发明的又一实施例提供了一种对某工区进行含浊沸石砂砾岩预测的方法，其中，该方法包括：

步骤101、确定含浊沸石砂砾岩的测井相应特征：利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的测井曲线响应特征；具体而言：

利用工区薄片的岩性分析数据确定各薄片的岩性；所述岩性包括含浊沸石砂砾岩以及非含浊沸石砂砾岩，其中，非含浊沸石砂砾岩包括致密砂砾岩、钙质胶结砂砾岩、含泥含砾砂砾岩和泥岩；

获取工区薄片对应的电阻率测井曲线、密度测井曲线(具体而言：确定薄片对应的深度，获取对应深度的电阻率测井曲线、密度测井曲线)；

分析含浊沸石砂砾岩的测井曲线响应特征；

得出含浊沸石砂砾岩具有高电阻率、低密度的测井曲线响应特征；

步骤102、确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值：利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；具体而言：

利用工区薄片的岩性分析数据确定各薄片的岩性；所述岩性包括含浊沸石砂砾岩以及非含浊沸石砂砾岩，其中，非含浊沸石砂砾岩包括致密砂砾岩、钙质胶结砂砾岩、含泥含砾砂砾岩和泥岩；

获取工区薄片对应的电阻率测井曲线、密度测井曲线(具体而言：确定薄片对应的深度，获取对应深度的电阻率测井曲线、密度测井曲线)；

基于确定的岩性以及获取的电阻率测井曲线、密度测井曲线，建立工区薄片含浊沸石砂砾岩以及非含浊沸石砂砾岩的电阻率-密度交会图版(结果如图5所示)；

基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；

由图5可以明确含浊沸石砂砾岩的值域范围为电阻率大于12Ω·m，密度值为2.36g/cm

步骤103、三维地震数据体经过地质统计学反演方法得出电阻率反演数据体：利用单井电阻率测井曲线作为敏感参数曲线，建立样本点，然后利用多井样本拟合建立概率密度函数，其次进行空间变差函数分析，确定空间变异程度，在三维地震数据的约束下，对井间储层进行统计评估，得出电阻率反演数据体，结果如图6A所示；

步骤104、三维地震数据体经过地质统计学反演方法得出密度反演数据体：利用单井密度测井曲线作为敏感参数曲线，建立样本点，然后利用多井样本拟合建立概率密度函数，其次进行空间变差函数分析，确定空间变异程度，在三维地震数据的约束下，对井间储层进行统计评估，得出密度反演数据体，结果如图6B所示；

步骤105、门槛值约束得到有效电阻率反演数据体：利用确定的电阻率门槛值对电阻率反演数据体进行约束：在电阻率反演数据体中，电阻率值小于等于12Ω·m的电阻率反演数据体赋值为零，电阻率值大于12Ω·m的电阻率反演数据体保留(即保持不变)，从而得到经电阻率门槛值约束后的有效电阻率反演数据体，以此作为含浊沸石砂砾岩的有效电阻率反演数据体，结果如图6C所示。

步骤106、门槛值约束得到有效密度反演数据体：利用确定的密度门槛值对密度反演数据体进行约束：在密度反演数据体中，密度小于2.36g/cm

步骤107、有效电阻率反演数据体与有效密度反演数据体融合得到能预测含浊沸石砂砾岩的融合数据体：以电阻率反演数据体作为基准有效反演数据体，并将密度反演数据体中的全部非零数据赋值为1得到第一有效反演数据体；然后将基准有效反演数据体与第一有效反演数据体相乘，得到能预测含浊沸石砂砾岩的融合数据体，结果如图6E所示；

利用该融合数据体能够有效识别哪些为含浊沸石砂砾岩，具体而言：在同一个线道号，同一时间深度点上，如果基准有效反演数据体与第一有效反演数据体都有非零值的情况下，最终得到的值也非零，也就是融合后的表征含浊沸石砂砾岩的有效电阻率-密度数据；而在同一个线道号，同一时间深度点上，只要有一个值为零，则最终得到的值也为零，也就是融合后的表征非含浊沸石砂砾岩的非有效电阻率-密度数据，基于此，该融合数据体能预测含浊沸石砂砾岩；

步骤108、利用层间数学方法得出含浊沸石砂砾岩的平面预测分布图：得到的融合数据体，利用地震解释的目标层位顶底界进行约束，运用均方根运算方法计算得到含浊沸石砂砾岩的平面预测分布图，结果如图7所示。

本发明实施例还提供了一种含浊沸石砂砾岩预测系统，优选地，该系统用于实现上述的方法实施例。

图8是根据本发明实施例的含浊沸石砂砾岩预测系统的结构框图，如图8所示，该装置包括：门槛值确定模块81、电阻率反演数据体建立模块82、密度反演数据体建立模块83、有效电阻率反演数据体建立模块84、有效密度反演数据体建立模块85、融合数据体建立模块86、浊沸石砂砾岩预测模块87；其中，

门槛值确定模块81：用于利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值；

电阻率反演数据体建立模块82：用于利用三维地震数据体反演得到电阻率反演数据体；

密度反演数据体建立模块83：用于利用三维地震数据体反演得到密度反演数据体；

有效电阻率反演数据体建立模块84：用于利用确定的电阻率门槛值对电阻率反演数据体进行约束得到含浊沸石砂砾岩的有效电阻率反演数据体；

有效密度反演数据体建立模块85：用于利用确定的密度门槛值对密度反演数据体进行约束得到含浊沸石砂砾岩的有效密度反演数据体；

融合数据体建立模块86：用于将有效电阻率反演数据体与有效密度反演数据体进行融合得到能预测含浊沸石砂砾岩的融合数据体；

浊沸石砂砾岩预测模块87：用于将得到的融合数据体，利用地震解释的目标层位顶底界进行约束，得到含浊沸石砂砾岩的平面预测分布图。

继续参见图9，上述门槛值确定模块81可以包括：

第一岩性获取子模块811、第一测井曲线获取子模块812、第一交会图版建立模块813、第一电阻率门槛值和密度门槛值确定子模块814；其中，

第一岩性获取子模块811：用于获取工区薄片的岩性分析数据确定薄片的岩性；所述岩性包括含浊沸石砂砾岩以及非含浊沸石砂砾岩；

第一测井曲线获取子模块812：用于获取工区薄片对应的电阻率测井曲线、密度测井曲线；

第一交会图版建立模块813：用于基于确定的岩性以及获取的电阻率测井曲线、密度测井曲线，建立工区薄片含浊沸石砂砾岩以及非含浊沸石砂砾岩的电阻率-密度交会图版；

第一电阻率门槛值和密度门槛值确定子模块814：用于基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

继续参见图10，上述门槛值确定模块81可以包括：第二岩性获取子模块8111、第二测井曲线获取子模块8121、第二交会图版建立模块8131、第二电阻率门槛值和密度门槛值确定子模块8141；

第二岩性获取子模块8111：用于获取工区薄片的岩性分析数据确定岩性为含浊沸石砂砾岩的薄片；

第二测井曲线获取子模块8121：用于获取确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线；

第二交会图版建立模块8131：用于根据确定的岩性为含浊沸石砂砾岩的薄片对应的电阻率测井曲线、密度测井曲线，建立含浊沸石砂砾岩的电阻率-密度交会图版；

第二电阻率门槛值和密度门槛值确定子模块8141：用于基于建立的电阻率-密度交会图版确定含浊沸石砂砾岩的电阻率门槛值和密度门槛值。

进一步，该系统可以进一步包括：测井曲线响应特征确定模块：用于利用工区薄片岩性分析数据以及对应电阻率测井曲线、密度测井曲线确定含浊沸石砂砾岩的测井曲线响应特征。

进一步，有效电阻率反演数据体建立模块84可以用于在电阻率反演数据体中，将电阻率值分布在电阻率门槛值之外的电阻率反演数据体赋值为零，电阻率值分布在电阻率门槛值之内的电阻率反演数据体不为零，从而得到经电阻率门槛值约束后的有效电阻率反演数据体，以此作为含浊沸石砂砾岩的有效电阻率反演数据体。其中，电阻率分布在电阻率门槛值之内的电阻率反演数据体保留(即保持不变)或者修改为其他任意不为零的数均可。

进一步，有效电阻率反演数据体建立模块85可以用于在密度反演数据体中，将密度值分布在密度门槛值之外的密度反演数据体赋值为零，密度值分布在密度门槛值之内的密度反演数据体不为零，从而得到经密度门槛值约束后的有效密度反演数据体，以此作为含浊沸石砂砾岩的有效密度反演数据体。其中，密度分布在电阻率门槛值之内的密度反演数据体保留(即保持不变)或者修改为其他任意不为零的数均可。

继续参见图11，上述融合数据体建立模块86可以包括：第一赋值子模块861、第一融合子模块862；

第一赋值子模块861：用于任意选定密度反演数据体和电阻率反演数据体中的一种作为基准有效反演数据体，并将另一种有效反演数据体中的全部非零数据赋值为1得到第一有效反演数据体；

第一融合子模块862：用于将基准有效反演数据体与第一有效反演数据体相乘，得到能预测含浊沸石砂砾岩的融合数据体。

上述各模块、子模块的具体实施过程，可以参见上述方法实施例的描述，在此不再赘述。

图12是根据本发明实施例的含浊沸石砂砾岩预测装置的示意图。图12所示的含浊沸石砂砾岩预测装置为通用数据处理装置，其包含通用的计算机硬件结构，其至少包含处理器1000、存储器1111；所述处理器1000用于执行所述存储器中存储的含浊沸石砂砾岩预测程序，以实现各方法实施例所述的含浊沸石砂砾岩预测方法(具体方法参见上述方法实施例的描述，在此不再赘述)。

本发明实施例还提供了一种计算机可读存储介质，所述存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现各方法实施例所述的含浊沸石砂砾岩预测方法(具体方法参见上述方法实施例的描述，在此不再赘述)。

以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的，因此权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外，由于本领域的技术人员容易想到很多修改和改变，因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作，而是可以涵盖落入其范围内的所有合适修改和等同物。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。