海底岩性薄砂体储层预测方法和装置、存储介质及终端与流程

本发明涉及海洋石油勘探领域，更具体地说，涉及一种海底岩性薄砂体储层预测方法和装置、存储介质及终端。

背景技术：

近年来，我国石油对外依存度已经超越警戒线的情况频频出现，出于国家石油安全考虑，海洋油气的勘探和开发急需增储上产以满足国内经济发展的需要。珠江口盆地(东部)的惠州凹陷已经被证实为最好的富生烃凹陷之一。自珠江口盆地1973年开展勘探工作以来，惠州凹陷一直是该区勘探工作的主战场。但随着勘探的不断深入，经过近40年的勘探开发，新近系及古近系珠海组这些埋深浅的地层的大型构造圈闭所剩无几，传统方式的勘探面临着严峻的形式。为了扩宽勘探视野，近几年已将勘探方向转移到深层的古近系勘探等新领域的复式勘探上。但仍面临一系列问题：如地震资料的品质有待提高、采集处理解释有进一步优化和提升的空间、深层开发作业成本高、施工难度大、技术要求高，等等。因此，继续拓展新的勘探领域，不仅能助力惠州凹陷油田群的持续开发，更能利用老油区生产设施齐全且便捷的特点更快地开发动用。其中的新勘探领域之一就是中浅层地层岩性圈闭的勘探。

目前，国内外岩性勘探技术可以主要分为两类和7项技术方法，分别为：1、两类。建立在调谐厚度理论上的应用及其延伸；利用波形分类等分析手段识别砂体。2、7项技术方法。(1)基于经典调谐理论的技术；(2)特色属性提取与分析方法；(3)基于地震沉积学的技术；(4)拟声波反演技术；(5)去薄层强屏蔽的目标处理技术；(6)基于压缩感知的薄层识别方法；(7)基于模式识别的薄砂体储层检测技术。这些技术和方法各擅其场，百家争鸣，但是又有各自的局限和不足，主要表现在经验参数的选取，特定地质条件的约束等，很难比较大范围的推广应用。

技术实现要素：

本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种海底岩性薄砂体储层预测方法和装置、存储介质及终端。

本发明解决其技术问题所采用的技术方案是：构造一种海底岩性薄砂体储层预测方法，包括：

s1、对薄砂体储层进行探测，获取所述薄砂体储层中一探测点的均方根振幅a和主频值f；

s2、根据所述均方根振幅a和主频值f获取所述探测点对应的fa因子；

s3、由所述fa因子得到所述探测点对应的储层厚度。

进一步，本发明所述的海底岩性薄砂体储层预测方法，所述步骤s2包括：

通过下述公式得到所述fa因子：

其中所述f(x)为所述fa因子。

进一步，本发明所述的海底岩性薄砂体储层预测方法，在所述步骤s3后所述方法还包括：

s4、重复执行所述步骤s1至步骤s3，获取所述薄砂体储层中多个探测点的储层厚度；

s5、由所有所述测量点的储层厚度得到所述薄砂体储层的厚度分布。

进一步，本发明所述的海底岩性薄砂体储层预测方法，在所述步骤s5后所述方法还包括：

s6、根据所述薄砂体储层的厚度分布预测所述薄砂体储层岩性尖灭线的位置，其中所述岩性尖灭线的尖灭点位于所述薄砂体储层的厚度最小处。

另，本发明还提供一种海底岩性薄砂体储层预测装置，包括：

第一探测单元，用于对薄砂体储层进行探测，获取所述薄砂体储层中一探测点的均方根振幅a和主频值f；

fa因子获取单元，用于根据所述均方根振幅a和主频值f获取所述探测点对应的fa因子；

第一厚度获取单元，用于由所述fa因子得到所述探测点对应的储层厚度。

进一步，本发明所述的海底岩性薄砂体储层预测装置，所述fa因子获取单元中所述fa因子为：

其中所述f(x)为所述fa因子。

进一步，本发明所述的海底岩性薄砂体储层预测装置，所述装置还包括：

第二探测单元，用于重复执行所述步骤s1至步骤s3，获取所述薄砂体储层中多个探测点的储层厚度；

第二厚度获取单元，用于由所有所述测量点的储层厚度得到所述薄砂体储层的厚度分布。

进一步，本发明所述的海底岩性薄砂体储层预测装置，所述装置还包括：

尖灭点获取单元，用于根据所述薄砂体储层的厚度分布预测所述薄砂体储层岩性尖灭线的位置，其中所述岩性尖灭线的尖灭点位于所述薄砂体储层的厚度最小处。

另，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如上述的海底岩性薄砂体储层预测方法。

另，本发明还提供一种终端，所述终端包括处理器，所述处理器用于执行存储器中存储的计算机程序时实现如上述海底岩性薄砂体储层预测方法的步骤。

实施本发明的一种海底岩性薄砂体储层预测方法和装置、存储介质及终端，具有以下有益效果：该方法包括：s1、对薄砂体储层进行探测，获取薄砂体储层中一探测点的均方根振幅a和主频值f；s2、根据均方根振幅a和主频值f获取探测点对应的fa因子；s3、由fa因子得到探测点对应的储层厚度。s4、重复执行步骤s1至步骤s3，获取薄砂体储层中多个探测点的储层厚度；s5、由所有测量点的储层厚度得到薄砂体储层的厚度分布。s6、根据薄砂体储层的厚度分布预测薄砂体储层岩性尖灭线的位置。本发明中fa因子综合了振幅、频率属性与砂体厚度的关系，通过算法使fa因子与砂体厚度呈现一一对应的关系，不存在调谐效应，从而在预测薄砂体储层时避免误差。

附图说明

下面将结合附图及实施例对本发明作进一步说明，附图中：

图1是本发明中振幅随砂体厚度变化示意图；

图2是本发明中频率随砂体厚度变化示意图；

图3是本发明一实施例中海底岩性薄砂体储层预测方法的流程图；

图4是本发明中fa因子值随砂体厚度变化示意图；

图5是本发明一实施例中海底岩性薄砂体储层预测方法的流程图；

图6是本发明一实施例中海底岩性薄砂体储层预测方法的流程图；

图7是本发明一种终端的结构示意图。

具体实施方式

为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。

实施例

本方法的目的首先在于提升海底地层薄砂体储层的预测精度和岩性尖灭点的识别，提高地层岩性圈闭的钻探成功率，达到降本增效的效果；其次，本方法在现有技术和方法的基础上提出对算法的优化，创新性提出fa因子的概念和算法，使得砂体储层厚度与fa因子呈现一一对应的关系，有利于储层厚度的准确预测。

算法的原理如图1所示，简要介绍如下：

根据经典的调谐厚度理论，在楔状体正演模型中，砂体厚度与振幅呈现先正相关，再负相关，再波动的特征。其中，在1/4λ的砂体厚度处振幅达到最大值。但相同的振幅值就会对应两个砂体厚度。这就导致仅使用单一振幅属性预测砂体厚度存在多解性。

同理，如图2所示，频率随砂体厚度的变化呈现先负相关，再正相关，再波动的特征。其中，在3/8λ的砂体厚度处频率达到最小值。也就造成了同一频率对应两个或多个砂体厚度。这就导致仅使用单一频率属性预测砂体厚度也存在多解性。

通过观察振幅随砂体厚度变化的关系和频率随砂体厚度变化的关系，总结两者属性的互补，创新性提出fa因子的概念和算法，使得fa因子值与砂体厚度值呈现一一对应的关系。如图3所示，本实施例的种海底岩性薄砂体储层预测方法包括下述步骤：

s1、对薄砂体储层进行探测，获取薄砂体储层中一探测点的均方根振幅a和主频值f；

s2、根据均方根振幅a和主频值f获取探测点对应的fa因子；进一步，本实施例的海底岩性薄砂体储层预测方法，步骤s2包括：

通过下述公式得到fa因子：

其中f(x)为fa因子。

s3、由fa因子得到探测点对应的储层厚度。

如图4所示。从而通过fa因子值预测砂体厚度，消除多解性。图中的红色散点为fa因子值，蓝色曲线为通过散点拟合的曲线，通过线性回归系数r2＝0.9956可知，曲线拟合精度较高，可以认为预测较为准确。

本实施例中fa因子综合了振幅、频率属性与砂体厚度的关系，通过算法使fa因子与砂体厚度呈现一一对应的关系，不存在调谐效应，通过fa因子获取探测点对应的储层厚度。

实施例

如图5所示，在上述实施例的基础上，本实施例的海底岩性薄砂体储层预测方法在步骤s3后方法还包括：

s4、重复执行步骤s1至步骤s3，获取薄砂体储层中多个探测点的储层厚度；

s5、由所有测量点的储层厚度得到薄砂体储层的厚度分布。

本实施例中fa因子综合了振幅、频率属性与砂体厚度的关系，通过算法使fa因子与砂体厚度呈现一一对应的关系，不存在调谐效应，通过fa因子获取探测点对应的储层厚度，通过采集多个探测点的储层厚度数据，得到整个探测区的薄砂体储层的厚度分布。进而可根据薄砂体储层的厚度分布进行井位部署。

实施例

如图6所示，在上述实施例的基础上，本实施例的海底岩性薄砂体储层预测方法在步骤s5后方法还包括：

s6、根据薄砂体储层的厚度分布预测薄砂体储层岩性尖灭线的位置，其中岩性尖灭线的尖灭点位于薄砂体储层的厚度最小处。优选地，在实际操作过程中，岩性尖灭线中的尖灭点位于薄砂体储层的预测厚度接近零处。

本实施例对上述实施例获取的薄砂体储层的厚度分布进行分析，得到探测区中薄砂体储层的储层厚度变化规律，通过变化规律即可得到薄砂体储层厚度逐渐变薄的地方，从而预测薄砂体储层岩性尖灭线的位置。

实施例

本实施例还提供一种海底岩性薄砂体储层预测装置，包括：

第一探测单元，用于对薄砂体储层进行探测，获取薄砂体储层中一探测点的均方根振幅a和主频值f；

fa因子获取单元，用于根据均方根振幅a和主频值f获取探测点对应的fa因子；进一步，本实施例的海底岩性薄砂体储层预测装置，fa因子获取单元中fa因子为：

其中f(x)为fa因子。

第一厚度获取单元，用于由fa因子得到探测点对应的储层厚度。

实施例

在上述实施例的基础上，本实施例的海底岩性薄砂体储层预测装置还包括：

第二探测单元，用于重复执行步骤s1至步骤s3，获取薄砂体储层中多个探测点的储层厚度；

第二厚度获取单元，用于由所有测量点的储层厚度得到薄砂体储层的厚度分布。

实施例

在上述实施例的基础上，本实施例的海底岩性薄砂体储层预测装置还包括：

尖灭点获取单元，用于根据薄砂体储层的厚度分布预测薄砂体储层岩性尖灭线的位置，其中岩性尖灭线的尖灭点位于薄砂体储层的厚度最小处。优选地，在实际操作过程中，岩性尖灭线中的尖灭点位于薄砂体储层的预测厚度接近零处。

实施例

本实施例的计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现如上述的海底岩性薄砂体储层预测方法。本实施例中海底岩性薄砂体储层预测方法使用算法简单，所需处理数据量小，所以在计算机程序安装时耗时少；在计算机程序执行时处理数据量小，处理速度快。

实施例

如图7所示，本实施例的终端包括处理器，处理器用于执行存储器中存储的计算机程序时实现如上述海底岩性薄砂体储层预测方法的步骤。作为选择，终端包括但不限于智能手机、平板电脑、笔记本电脑、台式电脑、服务器等。本实施例中海底岩性薄砂体储层预测方法使用算法简单，所需输入的参数少，操作简单，对操作人员的技能要求低，可快速上手工作。

本实施例应用珠江口盆地惠州凹陷实钻薄砂体储层的5口井进行验证后，模型井处砂体厚度预测误差0.15m，二次三维采集工区内井位砂体厚度预测误差0.13m，其余一次三维工区三口井位砂体厚度预测误差分别为-1.31m、1.35m和-0.90m，总体达到米级误差。应用于目标评价的勘探实践后，将构造油田a升级为岩性油田a-plus，岩性圈闭面积较构造圈闭面积增大约40％，新增单层储量两千多万方，提高约198.8％，油田整体储量增加约203％。

本发明中fa因子综合了振幅、频率属性与砂体厚度的关系，通过算法使fa因子与砂体厚度呈现一一对应的关系，不存在调谐效应，从而在预测薄砂体储层时避免误差。

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

专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(ram)、内存、只读存储器(rom)、电可编程rom、电可擦除可编程rom、寄存器、硬盘、可移动磁盘、cd-rom、或技术领域内所公知的任意其它形式的存储介质中。

以上实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并据此实施，并不能限制本发明的保护范围。凡跟本发明权利要求范围所做的均等变化与修饰，均应属于本发明权利要求的涵盖范围。