CN112698382B - 小尺度断控岩溶型储层地震预测方法及装置 - Google Patents

Abstract

本发明提供一种小尺度断控岩溶型储层地震预测方法及装置，该方法包括：获取研究区的地震数据体以及测井数据；将该地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据；根据该构造域弱地震反射特征数据以及该测井数据得到表征小尺度断控岩溶型储层的反演数据，通过将该地震数据体转化到wheeler域后进行特征提取，能消除地震数据体或其他属性体中由构造运动产生的变化，恢复古沉积环境，凸显小断距断裂的地震响应特征，以便有效识别弱振幅地震反射特征储层，识别小尺度断控岩溶型储层内部结构特征。

Description

小尺度断控岩溶型储层地震预测方法及装置

技术领域

本发明涉及石油勘探开发技术领域，尤其涉及一种小尺度断控岩溶型储层地震预测方法及装置。

背景技术

石油勘探开发成果表明，碳酸盐岩主要油气富集区大多与断裂有密切关系，油气富集于连通油源的断裂破碎带附近得岩溶发育区，成为目前探寻岩溶型油藏的重要方向。随着近年勘探开发的深入，在一些开发成熟区，尺度较大的岩溶孔洞储层已经钻探殆尽，剩余的多是储层空间较小的储层，其储层空间多以裂缝、溶蚀孔洞为主，储层与非储层的波阻抗差异小，呈弱能量地震反射特征，受走滑小断距断裂控制，并且由于上覆地层界面的地震强能量反射特征，掩盖了这类储层的弱振幅响应特征，因此具有较强的隐蔽性，这类弱振幅响应特征岩溶体也钻遇高产油气流。

现有小尺度断控岩溶型储层地震预测技术主要包括：1、利用三维地震数据体相干、曲率、蚂蚁追踪算法等常规属性断裂检测方法，表征断裂分布特征；2、利用三维地震数据体均方根属性或波阻抗表征岩溶储层分布特征；3、将断裂分布特征与岩溶储层分布特征叠合，从而得到断裂和溶蚀洞穴储层的分布。

现有技术的缺点包括：识别弱振幅地震反射特征储层的能力差、难以识别小尺度断控岩溶型储层内部结构特征，制约着储量潜力的充分挖掘、油藏的高效开发以及新井投产成功率的提高。

发明内容

针对现有技术中的问题，本发明提供一种小尺度断控岩溶型储层地震预测方法及装置、电子设备以及计算机可读存储介质，能够至少部分地解决现有技术中存在的问题。

为了实现上述目的，本发明采用如下技术方案：

第一方面，提供一种小尺度断控岩溶型储层地震预测方法，包括：

获取研究区的地震数据体以及测井数据；

将该地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；

对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据；

根据该构造域弱地震反射特征数据以及该测井数据得到表征小尺度断控岩溶型储层的反演数据。

进一步地，该将该地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据，包括：

将该地震数据体转化到wheeler域；

提取wheeler域的地震数据体中具有地震波相位一致性的波组；

利用该wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据。

进一步地，该提取wheeler域的地震数据体中具有地震波相位一致性的波组，包括：

对wheeler域的地震数据体进行核主成分分析得到具有地震波相位一致性的波组。

进一步地，该利用该wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据，包括：

利用该wheeler域的地震数据体减去提取到的波组得到wheeler域弱地震反射特征数据。

进一步地，该根据该构造域弱地震反射特征数据以及该测井数据得到表征小尺度断控岩溶型储层的反演数据，包括：

根据该构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体；

根据该断裂检测属性体、该测井数据以及该构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据。

进一步地，该根据该构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体，包括：

提取该构造域弱地震反射特征数据中目的层相干属性得到断裂发育的方位；

根据该方位旋转该构造域弱地震反射特征数据的测网信息；

对旋转后的数据进行频谱分解得到不同频率的振幅体和相位体；

对不同频率的振幅体和相位体进行边沿增强；

根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体。

进一步地，该根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体，包括：

对边缘增强后的不同频率的振幅体和相位体进行自适应的核主成分分析得到反映小尺度的断裂检测属性体。

进一步地，该测井数据包括：已知井的阻抗；

该根据该断裂检测属性体、该测井数据以及该构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据，包括：

对该断裂检测属性体与该已知井的阻抗进行交汇分析得到断裂检测属性体与波阻抗体的转换关系；

利用该转换关系将该断裂检测属性体转换成对应的波阻抗体；

将该波阻抗体作为初始模型对该构造域弱地震反射特征数据进行反演得到表征小尺度断控岩溶型储层的反演数据。

第二方面，提供一种小尺度断控岩溶型储层地震预测装置，包括：

数据获取模块，获取研究区的地震数据体以及测井数据；

数据转换提取模块，将该地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；

数据反变换模块，对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据；

数据处理模块，根据该构造域弱地震反射特征数据以及该测井数据得到表征小尺度断控岩溶型储层的反演数据。

进一步地，该数据转换提取模块，包括：

数据转换子模块，将该地震数据体转化到wheeler域；

波组提取子模块，提取wheeler域的地震数据体中具有地震波相位一致性的波组；

特征增强子模块，利用该wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据。

进一步地，该波组提取子模块包括：

核主成分分析单元，对wheeler域的地震数据体进行核主成分分析得到具有地震波相位一致性的波组。

进一步地，该波组提取子模块包括：

计算单元，利用该wheeler域的地震数据体减去提取到的波组得到wheeler域弱地震反射特征数据。

进一步地，该数据处理模块包括：

属性体获取子模块，根据该构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体；

反演子模块，根据该断裂检测属性体、该测井数据以及该构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据。

进一步地，该属性体获取子模块包括：

方位提取单元，提取该构造域弱地震反射特征数据中目的层相干属性得到断裂发育的方位；

测网信息旋转单元，根据该方位旋转该构造域弱地震反射特征数据的测网信息；

频谱分解单元，对旋转后的数据进行频谱分解得到不同频率的振幅体和相位体；

边沿增强单元，对不同频率的振幅体和相位体进行边沿增强；

属性体获取单元，根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体。

进一步地，该属性体获取单元包括：

自适应核主成分分析子单元，对边缘增强后的不同频率的振幅体和相位体进行自适应的核主成分分析得到反映小尺度的断裂检测属性体。

进一步地，该测井数据包括：已知井的阻抗；

该反演子模块包括：

交汇分析单元，对该断裂检测属性体与该已知井的阻抗进行交汇分析得到断裂检测属性体与波阻抗体的转换关系；

转换单元，利用该转换关系将该断裂检测属性体转换成对应的波阻抗体；

地震数据反演单元，将该波阻抗体作为初始模型对该构造域弱地震反射特征数据进行反演得到表征小尺度断控岩溶型储层的反演数据。

第三方面，提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，该处理器执行该程序时实现上述的小尺度断控岩溶型储层地震预测方法的步骤。

第四方面，提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述的小尺度断控岩溶型储层地震预测方法的步骤。

本发明提供的小尺度断控岩溶型储层地震预测方法及装置，该方法包括：获取研究区的地震数据体以及测井数据；将该地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据；根据该构造域弱地震反射特征数据以及该测井数据得到表征小尺度断控岩溶型储层的反演数据，通过将该地震数据体转化到wheeler域后进行特征提取，能够消除地震数据体或者其他属性体中由构造运动产生的变化，恢复古沉积环境，直观了解同时期沉积环境下地层的情况，断溶体上覆连续强反射同相轴被恢复到水平状态，地震波的相位具有较好的一致性，wheeler域弱地震反射特征数据凸显了小断距断裂的地震响应特征，以便有效识别弱振幅地震反射特征储层，识别小尺度断控岩溶型储层内部结构特征，利于储量潜力的充分挖掘、油藏的高效开发以及新井投产成功率的提高。

为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。

附图说明

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

图1为本发明实施例中的服务器S1与客户端设备B1之间的架构示意图；

图2是本发明实施例中的小尺度断控岩溶型储层地震预测方法的流程示意图；

图3示出了本发明实施例中步骤S200的具体步骤；

图4示出了某研究区的原始地震剖面；

图5示出了采用本发明实施例提供的技术处理后的某研究区的去除地层地震反射特征的地震剖面；

图6示出了本发明实施例中的步骤S400的具体步骤；

图7示出了本发明实施例中的步骤S410的具体步骤；

图8示出了本发明实施例中的步骤S420的具体步骤；

图9示出了本发明实施例中的HA6油田的地震原始剖面、去地层强轴后反演剖面以及断控岩溶储层内部结构剖面。

图10示出了本发明实施例中的HA6油田的断溶体油藏边界划分及评价图；

图11是本发明实施例中的小尺度断控岩溶型储层地震预测装置的结构框图；

图12示出了本发明实施例中的数据转换提取模块的结构框图；

图13示出了本发明实施例中的数据处理模块的结构框图；

图14示出了本发明实施例中的属性体获取子模块的结构框图；

图15示出了本发明实施例中的反演子模块的结构框图；

图16为本发明实施例电子设备的结构图。

具体实施方式

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

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

需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。

小尺度断控岩溶型储层地震预测主要存在两个难点：①振幅能量弱、断层附近储层识别难度大；②同时由于一些走滑断层断距小，表征断层和溶洞(断溶体系统)难度大；③内部岩溶结构识别难。

由于小尺度断控岩溶型储层上覆地层界面往往是一连续地震强振幅反射特征，这种区域性连续强反射往往会对下伏储层地震相应特征有屏蔽作用，特别是孔洞-裂缝型储层地震相应特征相对较弱，易受上覆地层强反射屏蔽，而难以识别，本发明实施例将地震反射同相轴看成有地层界面地震响应特征以及地层界面之下的储层地震响应特征的叠加，提取隐蔽地震反射特征，去除影响小断距断裂识别的地层地震相应特征，识别小尺度断裂轮廓，利用断裂的轮廓作为约束条件参与断控岩溶储层内部结构反演，既考虑的断裂的控制因素又能表征岩溶作用对储层的影响，精细刻画小尺度断控岩溶型储层轮廓及其内部结构。

有鉴于此，本申请提供了一种小尺度断控岩溶型储层地震预测装置，该装置可以为一种服务器S1，参见图1，该服务器S1可以与至少一个客户端设备B1通信连接，所述客户端设备B1可以将研究区的地震数据体以及测井数据发送至所述服务器S1，所述服务器S1可以在线接收所述研究区的地震数据体以及测井数据。所述服务器S1可以在线或者离线对获取的研究区的地震数据体以及测井数据进行预处理，将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；根据所述wheeler域弱地震反射特征数据以及所述测井数据得到表征小尺度断控岩溶型储层的反演数据。而后，所述服务器S1可以将表征小尺度断控岩溶型储层的反演数据在线发送至所述客户端设备B1。所述客户端设备B1可以在线接收所述表征小尺度断控岩溶型储层的反演数据。

可以理解的是，所述客户端设备B1可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。

在实际应用中，进行小尺度断控岩溶型储层地震预测的部分可以在如上述内容所述的服务器S1侧执行，即，如图1所示的架构，也可以所有的操作都在所述客户端设备B1中完成，且该所述客户端设备B1可以直接与数据库服务器S2进行通信连接。具体可以根据所述客户端设备B1的处理能力，以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备B1中完成，所述客户端设备B1还可以包括处理器，用于进行小尺度断控岩溶型储层地震预测的具体处理。

所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信，包括在本申请提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然，所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol，远程过程调用协议)、REST协议(Representational State Transfer，表述性状态转移协议)等。

图2是本发明实施例中的小尺度断控岩溶型储层地震预测方法的流程示意图；如图2所示，该小尺度断控岩溶型储层地震预测方法可以包括以下内容：

步骤S100：获取研究区的地震数据体以及测井数据；

具体地，研究区的地震数据体由地震数据处理后得到。

步骤S200：将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；

其中，将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据，以便去除上覆地层强反射的屏蔽作用，实现小尺度断裂及孔洞储层地震响应特征提取。

具体地，将地震数据体转化到wheeler域实现地层体转换，目的是消除地震数据体或者其他属性体中由构造运动产生的变化，恢复古沉积环境，直观了解同时期沉积环境下地层的情况。

步骤S300：对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据。

步骤S400：根据所述构造域弱地震反射特征数据以及所述测井数据得到表征小尺度断控岩溶型储层的反演数据。

通过采用上述技术方案，实现对弱地震响应的断裂识别，能更好、更精确的提供断裂信息；弱振幅隐蔽岩溶储层地震相应特征，对断控岩溶储层发育带边界更清晰；能够刻画断控岩溶缝洞储层内部结构。

在一个可选的实施例中，参见图3，该步骤S200可以包括以下内容：

步骤S210：将所述地震数据体转化到wheeler域；

其中，wheeler域的地震数据体的地震反射同相轴的产状、相位等信息与地层、沉积等信息的相关性较强，断溶体上覆连续强反射同相轴被恢复到水平状态，地震波的相位此时具有较好的一致性。

步骤S220：提取wheeler域的地震数据体中具有地震波相位一致性的波组；

具体地，对wheeler域的地震数据体进行核主成分分析得到具有地震波相位一致性的波组。

步骤S230：利用所述wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据。

具体地，利用所述wheeler域的地震数据体减去提取到的波组，得到wheeler域弱地震反射特征数据。

通过采用上述技术方案，能够实现数据滤波过程，滤除低频的特征，凸显小断距断裂的地震响应特征，实现了波形分离，参见图4和图5，去除断溶体上覆地层强反射地震波形，剩余地震波形将原始弱反射地震波凸显出来增强，反映出断溶体内部弱反射储层发育特征，同时可以使地震剖面上断裂系统更加清晰。

在一个可选的实施例中，参见图6，该步骤S400可以包括以下内容：

步骤S410：根据所述构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体；

通过对构造域弱地震反射特征数据进行处理，实现小尺度断控岩溶储层边界识别。

步骤S420：根据所述断裂检测属性体、所述测井数据以及所述构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据。

其中，通过采用断裂检测属性体、所述测井数据以及所述构造域弱地震反射特征数据，能够实现小尺度断控岩溶储层内部结构反演。

在一个可选的实施例中，参见图7，该步骤S410可以包括以下内容：

步骤S411：提取所述构造域弱地震反射特征数据中目的层相干属性得到断裂发育的方位；

步骤S412：根据所述方位旋转所述构造域弱地震反射特征数据的测网信息；

其中，通过旋转所述构造域弱地震反射特征数据的测网信息，使地震测网主测线方向与断裂发育方位垂直或近似垂直。

步骤S413：对旋转后的数据进行频谱分解得到不同频率的振幅体和相位体；

其中，对旋转后的数据进行高分辨率的频谱分解。

步骤S414：对不同频率的振幅体和相位体进行边沿增强；

具体地，对不同频率的振幅体和相位体进行边沿增强，识别不同频率上波形、振幅和相位等多种优选的不连续性属性，如相干、曲率、最大似然概率体、梯度结构张量等。

步骤S415：根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体。

具体地，对边缘增强后的不同频率的振幅体和相位体进行自适应的核主成分分析得到反映小尺度的断裂检测属性体。

通过采用上述技术方案，可以实现小尺度边缘检测功能，可以得到反映沉积结构的属性体，相比常规断裂识别方法，具有更强的抗噪性；相比断裂识别属性，突出横向地层内的不连续性，有利于复杂沉积体边界刻画，异常地质体的边界更加清晰。

在一个可选的实施例中，测井数据包括：已知井的阻抗，参见图8，该步骤S420可以包括以下内容：

步骤S421：对所述断裂检测属性体与所述已知井的阻抗进行交汇分析得到断裂检测属性体与波阻抗体的转换关系；

步骤S422：利用所述转换关系将所述断裂检测属性体转换成对应的波阻抗体；

步骤S423：将所述波阻抗体作为初始模型对所述构造域弱地震反射特征数据进行反演得到表征小尺度断控岩溶型储层的反演数据。

具体地，在去除断控岩溶储层上覆地层界面强反射的基础上，将断裂检测属性体与已知井的阻抗进行交汇分析，得到断裂检测属性体转换成波阻抗体的关系式，利用该关系式把断裂检测属性体转换成波组抗体，作为波阻抗反演的初始模型，再将该初始模型作为约束条件，应用到反演迭代流程中，得到既能表征断溶体轮廓特征又能刻画内部结构特征的反演数据。

其中，在储层定量或半定量描述中，常用到波阻抗反演技术，但常规波阻抗反演只能反映大型的溶洞型储层，对纵向上发育的断控岩溶体刻画效果较差，本发明实施例基于去除地层信息的断控岩溶储层地震波形指示反演，能够得到既能表征断溶体轮廓特征又能刻画内部结构特征的反演数据。

为了验证本发明实施例提供的小尺度断控岩溶型储层地震预测技术的效果，在HA6油田实施应用，HA6油田位于塔北隆起斜坡区，发育中下奥陶统主要发育断裂裂缝、洞穴和孔洞三种类型储集体，储层发育规模受控于多期活动的走滑断裂体系及顺断裂溶蚀作用，为典型的断溶体油藏，埋藏较深(>5 500m)，发育玉尔吐斯组烃源岩，油源断裂发育，油气成藏条件优越。目的层各组地层较平缓，表现为向东南方向倾斜的单斜构造,现今构造面貌是多期构造运动叠加改造的结果，断裂展布及发育受控于多期次的构造运动。该研究区小尺度断控岩溶型储层上覆地层界面是一连续地震强振幅反射特征。这种区域性连续强反射会对下伏储层地震相应特征有屏蔽作用，特别是孔洞-裂缝型储层地震相应特征相对较弱，易受上覆地层强反射屏蔽，而难以识别。图9的(a)是该研究区的地震原始剖面，图9的(b)为基于去除地层反射特征地震剖面，图9的(c)为按照单井岩相岩石物理分析将反演结果划分不同的岩相带,钻遇红色区带GR值升高说明该区被泥质充填，从该区带的形态看，应为断裂发育区，也体现断溶体纵向长度远大于横向宽度，横向上分出主断裂区、破碎带、围岩带的地质结构特征。

采用断溶体内部结构识别技术，进行断溶体内部结构划分，为断裂分段评价提供基础资料，解释了不同断裂带之间、同一断溶体不同部位含油气规模的差异。在断溶体识别的基础上进行断溶体发育带的划分和评价，分小层统计各储层时间厚度，分别乘以各储层的平均速度以求得各层段储层的厚度，最后将各层段储层的厚度相加即可得到相应的储层厚度，将预测储层厚度D1与钻井储层厚度D2对比，预测误差e小于4m，参见表1。依据储层厚度及是否有油源断裂匹配，评价出断溶体油气富集带，如图10，等值线为储层厚度趋势线，黑色实线为断溶体油藏单元边界。断溶体划分和评价为断溶体油藏勘探和高效开发井布署、剩余油挖潜发挥了重要的作用，近年来新钻井26口钻井成功率由原来的78％提高到86％。

表1反演预测储层厚度与实钻厚度对比

基于同一发明构思，本申请实施例还提供了一种小尺度断控岩溶型储层地震预测装置，可以用于实现上述实施例所描述的方法，如下面的实施例所述。由于小尺度断控岩溶型储层地震预测装置解决问题的原理与上述方法相似，因此小尺度断控岩溶型储层地震预测装置的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图11是本发明实施例中的小尺度断控岩溶型储层地震预测装置的结构框图。如图11所示，该小尺度断控岩溶型储层地震预测装置具体包括：数据获取模块10、数据转换提取模块20、数据反变换模块30以及数据处理模块40。

数据获取模块10获取研究区的地震数据体以及测井数据；

数据转换提取模块20将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；

数据反变换模块30对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据。

数据处理模块40根据所述构造域弱地震反射特征数据以及所述测井数据得到表征小尺度断控岩溶型储层的反演数据。

通过采用上述技术方案，实现对弱地震响应的断裂识别，能更好、更精确的提供断裂信息；弱振幅隐蔽岩溶储层地震相应特征，对断控岩溶储层发育带边界更清晰；能够刻画断控岩溶缝洞储层内部结构。

在一个可选的实施例中，参见图12，该数据转换提取模块，包括：数据转换子模块21、波组提取子模块22以及特征增强子模块23。

数据转换子模块21将所述地震数据体转化到wheeler域；

波组提取子模块22提取wheeler域的地震数据体中具有地震波相位一致性的波组；

特征增强子模块23利用所述wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据。

在一个可选的实施例中，波组提取子模块包括：核主成分分析单元，对wheeler域的地震数据体进行核主成分分析得到具有地震波相位一致性的波组。

在一个可选的实施例中，波组提取子模块包括：计算单元，利用所述wheeler域的地震数据体减去提取到的波组得到wheeler域弱地震反射特征数据。

在一个可选的实施例中，参见图13，数据处理模块包括：属性体获取子模块31以及反演子模块32。

属性体获取子模块31根据所述构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体；

反演子模块32根据所述断裂检测属性体、所述测井数据以及所述构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据。

在一个可选的实施例中，参见图14，属性体获取子模块包括：方位提取单元311、测网信息旋转单元312、频谱分解单元313、边沿增强单元314以及属性体获取单元315。

方位提取单元311提取所述构造域弱地震反射特征数据中目的层相干属性得到断裂发育的方位；

测网信息旋转单元312根据所述方位旋转所述构造域弱地震反射特征数据的测网信息；

频谱分解单元313对旋转后的数据进行频谱分解得到不同频率的振幅体和相位体；

边沿增强单元314对不同频率的振幅体和相位体进行边沿增强；

属性体获取单元315根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体。

在一个可选的实施例中，属性提获取单元包括：自适应核主成分分析子单元，对边缘增强后的不同频率的振幅体和相位体进行自适应的核主成分分析得到反映小尺度的断裂检测属性体。

在一个可选的实施例中，测井数据包括：已知井的阻抗；参见图15，该反演子模块包括：交汇分析单元321、转换单元322以及地震数据反演单元323。

交汇分析单元321对所述断裂检测属性体与所述已知井的阻抗进行交汇分析得到断裂检测属性体与波阻抗体的转换关系；

转换单元322利用所述转换关系将所述断裂检测属性体转换成对应的波阻抗体；

地震数据反演单元323将所述波阻抗体作为初始模型对所述构造域弱地震反射特征数据进行反演得到表征小尺度断控岩溶型储层的反演数据。

上述实施例阐明的装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为电子设备，具体的，电子设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。

在一个典型的实例中电子设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述的小尺度断控岩溶型储层地震预测方法的步骤。

下面参考图16，其示出了适于用来实现本申请实施例的电子设备600的结构示意图。

如图16所示，电子设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。

以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。

特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述的小尺度断控岩溶型储层地震预测方法的步骤。

在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。

为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

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

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

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

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。

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

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

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (4)

1.一种小尺度断控岩溶型储层地震预测方法，其特征在于，包括：

获取研究区的地震数据体以及测井数据；

将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；

对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据；

根据所述构造域弱地震反射特征数据以及所述测井数据得到表征小尺度断控岩溶型储层的反演数据；

其中，所述将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据，包括：

将所述地震数据体转化到wheeler域；

提取wheeler域的地震数据体中具有地震波相位一致性的波组；

利用所述wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据；

其中，所述提取wheeler域的地震数据体中具有地震波相位一致性的波组，包括：

对wheeler域的地震数据体进行核主成分分析得到具有地震波相位一致性的波组；

其中，所述利用所述wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据，包括：

利用所述wheeler域的地震数据体减去提取到的波组得到wheeler域弱地震反射特征数据；

其中，所述根据所述构造域弱地震反射特征数据以及所述测井数据得到表征小尺度断控岩溶型储层的反演数据，包括：

根据所述构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体；

根据所述断裂检测属性体、所述测井数据以及所述构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据；

其中，所述根据所述构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体，包括：

提取所述构造域弱地震反射特征数据中目的层相干属性得到断裂发育的方位；

根据所述方位旋转所述构造域弱地震反射特征数据的测网信息；

对旋转后的数据进行频谱分解得到不同频率的振幅体和相位体；

对不同频率的振幅体和相位体进行边缘增强；

根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体；

其中，所述根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体，包括：

对边缘增强后的不同频率的振幅体和相位体进行自适应的核主成分分析得到反映小尺度的断裂检测属性体；

其中，所述测井数据包括：已知井的阻抗；

其中，所述根据所述断裂检测属性体、所述测井数据以及所述构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据，包括：

对所述断裂检测属性体与所述已知井的阻抗进行交汇分析得到断裂检测属性体与波阻抗体的转换关系；

利用所述转换关系将所述断裂检测属性体转换成对应的波阻抗体；

将所述波阻抗体作为初始模型对所述构造域弱地震反射特征数据进行反演得到表征小尺度断控岩溶型储层的反演数据。

2.一种小尺度断控岩溶型储层地震预测装置，其特征在于，包括：

数据获取模块，获取研究区的地震数据体以及测井数据；

数据转换提取模块，将所述地震数据体转化到wheeler域后进行特征提取得到wheeler域弱地震反射特征数据；

数据反变换模块，对所述wheeler域弱地震反射特征数据进行反wheeler域变换得到构造域弱地震反射特征数据；

数据处理模块，根据所述构造域弱地震反射特征数据以及所述测井数据得到表征小尺度断控岩溶型储层的反演数据；

其中，所述数据转换提取模块，包括：

数据转换子模块，将所述地震数据体转化到wheeler域；

波组提取子模块，提取wheeler域的地震数据体中具有地震波相位一致性的波组；

特征增强子模块，利用所述wheeler域的地震数据体与提取到的波组得到wheeler域弱地震反射特征数据；

其中，所述波组提取子模块包括：

核主成分分析单元，对wheeler域的地震数据体进行核主成分分析得到具有地震波相位一致性的波组；

其中，所述波组提取子模块包括：

计算单元，利用所述wheeler域的地震数据体减去提取到的波组得到wheeler域弱地震反射特征数据；

其中，所述数据处理模块包括：

属性体获取子模块，根据所述构造域弱地震反射特征数据得到反映小尺度的断裂检测属性体；

反演子模块，根据所述断裂检测属性体、所述测井数据以及所述构造域弱地震反射特征数据得到表征小尺度断控岩溶型储层的反演数据；

其中，所述属性体获取子模块包括：

方位提取单元，提取所述构造域弱地震反射特征数据中目的层相干属性得到断裂发育的方位；

测网信息旋转单元，根据所述方位旋转所述构造域弱地震反射特征数据的测网信息；

频谱分解单元，对旋转后的数据进行频谱分解得到不同频率的振幅体和相位体；

边缘增强单元，对不同频率的振幅体和相位体进行边缘增强；

属性体获取单元，根据边缘增强后的不同频率的振幅体和相位体得到反映小尺度的断裂检测属性体；

其中，所述属性体获取单元包括：

自适应核主成分分析子单元，对边缘增强后的不同频率的振幅体和相位体进行自适应的核主成分分析得到反映小尺度的断裂检测属性体；

其中，所述测井数据包括：已知井的阻抗；

其中，所述反演子模块包括：

交汇分析单元，对所述断裂检测属性体与所述已知井的阻抗进行交汇分析得到断裂检测属性体与波阻抗体的转换关系；

转换单元，利用所述转换关系将所述断裂检测属性体转换成对应的波阻抗体；

地震数据反演单元，将所述波阻抗体作为初始模型对所述构造域弱地震反射特征数据进行反演得到表征小尺度断控岩溶型储层的反演数据。

3.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1所述的小尺度断控岩溶型储层地震预测方法的步骤。

4.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1所述的小尺度断控岩溶型储层地震预测方法的步骤。