CN103454686B - 基于地层切片的小尺度沉积相进行储层预测的方法及系统 - Google Patents

Abstract

本发明提供一种基于地层切片的小尺度沉积相进行储层预测的方法及系统，该方法包括：采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；根据测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相；根据测井曲线以及地震资料确定出区分所述的沉积微相的地震属性；根据沉积微相以及地震属性提取井点地震属性；根据沉积微相以及井点地震属性确定数学映射关系；根据地震资料生成地层切片；根据沉积微相与特征地震属性的数学映射关系将地层切片转换为沉积相平面分布图；根据沉积相平面分布图进行储层预测。实现了井点控制下的地层切片数据自动转为沉积相平面图，进而可根据沉积相平面图进行储层预测。

Description

基于地层切片的小尺度沉积相进行储层预测的方法及系统

技术领域

本发明关于石油勘探技术领域，特别是关于基于地震沉积学的储层预测技术，具体的讲是一种基于地层切片的小尺度沉积相进行储层预测的方法及系统。

背景技术

沉积相的研究一直是油气勘探开发的基础，随着油气勘探开发程度的不断深入，勘探对象日益复杂，逐渐由构造油气藏转向岩性油气藏，由斜坡区转向深洼区，这无疑对沉积相的研究提出了更精细的要求。以往研究中以组、段(大于1/4波长)等大尺度级别为单位的沉积相图远不能满足岩性油气藏的勘探需求。高精度三维地震资料、测井资料、钻井资料为小尺度(单砂层，小于1/4波长)的沉积相研究提供重要的资料基础，但是，现有技术中对于沉积相研究的作图尺度依然比较大，无法满足精细油气藏勘探的需求。另一方面，研究人员目前编绘沉积相图的方法主要是沉积相模式指导、测井相控制、地震属性约束、手工编绘，精度低、效率低，无法满足研究需求。

发明人通过大量调研发现，随着地震沉积学研究技术的广泛应用，特别是等时地层切片技术应用于沉积体系和沉积相平面展布研究中，使得沉积相的研究尺度精细了很多，但在以往的研究中，基于地层切片的沉积相分析主要通过沉积模式指导下、定性程度、手动编辑的工作方式完成小尺度沉积相平面图的编绘，作图精度低、效率低，没有充分挖掘沉积特征与地震属性特征之间的内在联系，没有将二者很好的结合起来，不能实现地震数据到地质成果的快速自动转化，因此未能形成一种将地层切片数据定量、快速、有效、自动转为沉积相平面图进而进行储层预测的研究技术方案。

发明内容

为了解决上述技术问题，本发明提供了一种地层切片的小尺度沉积相进行储层预测的方法及系统，借助井点沉积相的可靠解释成果和等时地层切片清晰反映的沉积特征，在小尺度研究范围内，充分挖掘沉积特征与地震属性之间的内在联系，实现了井点控制下的地层切片数据定量、快速、有效、自动转为沉积相平面图，进而根据沉积相平面图进行储层预测。

本发明的目的之一是，提供一种基于地层切片的小尺度沉积相进行储层预测的方法，包括：采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；根据所述的测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相；根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性；根据所述的沉积微相以及所述的地震资料提取井点地震属性；根据所述的沉积微相以及井点地震属性确定数学映射关系；根据所述的地震资料生成地层切片；根据所述的沉积微相与特征地震属性的数学映射关系将所述的地层切片转换为沉积相平面分布图；根据所述的沉积相平面分布图进行储层预测。

本发明的目的之一是，提供了一种基于地层切片的小尺度沉积相进行储层预测的系统，包括：采集装置，用于采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；沉积微相解释划分装置，用于根据所述的测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相；地震属性确定装置，用于根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性；井点地震属性提取装置，用于根据所述的沉积微相以及所述的地震资料提取井点地震属性；数学映射关系确定装置，用于根据所述的沉积微相以及井点地震属性确定数学映射关系；地层切片生成装置，用于根据所述的地震资料生成地层切片；沉积分布图转换装置，用于根据所述的沉积微相与特征地震属性的数学映射关系将所述的地层切片转换为沉积相平面分布图；储层预测装置，用于根据所述的沉积相平面分布图进行储层预测。

本发明的有益效果在于，提供了一种地层切片的小尺度沉积相进行储层预测的方法及系统，实现了依据地层切片进行沉积相自动解释，极大提高了工作效率；井点地震属性快速提取技术极大提高了测井-地震-地质相结合的定量化研究程度；该技术方案原理简单、操作方便、技术上易于实现。

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

附图说明

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

图1为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测的方法的流程图；

图2为图1中的步骤S103的具体流程图；

图3为图1中的步骤S104的具体流程图；

图4为图1中的步骤S105的具体流程图；

图5为图1中的步骤S106的具体流程图；

图6为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测系统的结构框图；

图7为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测系统中地震属性确定装置的结构框图；

图8为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测系统中井点地震属性提取装置的结构框图；

图9为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测的系统中数学映射关系确定装置的结构框图；

图10为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测的系统中地层切片生成装置的结构框图；

图11为qn6井沙一中井沙一中重力流水道测井相特征示意图；

图12为ch52井沙一中深湖测井相特征示意图；

图13为qn6井精细合成记录示意图；

图14为ch52井精细合成记录示意图；

图15为QN地区井点平均振幅属性与伽马、沉积微相交汇图示意图；

图16为wheeler域沙一中某振幅地层切片示意图；

图17为数据处理后的地层切片示意图；

图18为自动沉积相解释结果示意图。

具体实施方式

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

本发明的主要目的在于借助井点沉积相的可靠解释成果和等时地层切片清晰反映的沉积特征，在小尺度研究范围内，充分挖掘沉积特征与地震属性之间的内在联系，提供一种井点控制下的地层切片数据定量、快速、有效、自动转为沉积相平面图的研究技术方法。

本发明提及的小尺度沉积相的研究范围是指小于通常沉积相图的组、段编制尺度，定量来说小于1/4波长厚度，达到单砂层级别，甚至更小，主要依据地层切片采样率而定。

图1为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测的方法的流程图，由图1可知，该方法具体包括：

S101：采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；

S102：根据所述的测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相。该步骤在进行沉积微相的划分时，主要应用到测井曲线中的自然伽马曲线、自然电位曲线、微电极电阻率曲线的幅度特征、形态特征、光滑度三个参数，结合岩心资料中的岩心描述、钻井取心等地质对沉积微相进行划分。得到各个类型的沉积微相可按照微相类别依次定名为1、2、3、4、……。

S103：根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性。步骤S103的具体流程图如图2所示。

S104：根据所述的沉积微相以及所述的地震资料提取井点地震属性。步骤S104的具体流程图如图3所示。

S105：根据所述的沉积微相以及井点地震属性确定数学映射关系。步骤S105的具体流程图如图4所示。

S106：根据所述的地震资料生成地层切片。步骤S106的具体流程图如图5所示。

S107：根据所述的数学映射关系将所述的地层切片转换为沉积相平面分布图。将沉积微相与特征地震属性的数学映射关系应用到地层切片上，并以该地层切片的地震属性值为约束条件，从井点出发，外推到无井区，实现地层切片的自动沉积相解释。

S108：根据所述的沉积相平面分布图进行储层预测。

图2为图1中的步骤S103的具体流程图，由图2可知，该步骤具体包括：

S201：从所述的测井曲线中提取声波曲线、密度曲线；

S202：根据所述的声波曲线、密度曲线以及地震资料制作当前钻井的精细合成记录；

S203：根据所述的精细合成记录以及地震资料分析所述目的层的地震反射特征，所述地震反射特征包括振幅、频率、相位。

S204：根据所述的地震反射特征确定出区分所述的沉积微相的地震属性。

该步骤主要分析目的层的各种沉积微相的振幅、频率等基本特征，初步确定出区分沉积微相的地震属性。所述的地震属性属于体属性范畴，不包括时窗属性。

图3为图1中的步骤S104的具体流程图，由图3可知，该步骤具体包括：

S301：获取所述的沉积微相对应的顶深度、底深度；

S302：获取当前钻井的精细合成记录，此处的精细合成记录即为步骤S202中确定出的合成记录。

S303：根据所述的精细合成记录将所述沉积微相的顶深度转换为地震时间域的顶时间范围；

S304：根据所述的精细合成记录将所述沉积微相的底深度转换为地震时间域的底时间范围；

S305：根据地震资料在所述的顶时间范围、所述的底时间范围内提取井点地震属性。

在具体的实施方式中，诸如在井旁道范围依据默认提取井眼为中心的7个地震道，也可在井旁道范围依据实际地层倾斜程度而定，地层越倾斜，井旁道范围越小，地层越平缓，井旁道范围可以适当变大，一般为1～10；若微相纵向时间范围较小，不能满足提取某种地震属性的需要，则默认为提取该地震属性需要的最小时间范围。井点地震属性的提取是在井旁道区域内，根据距井眼的距离加权平均的方法，距离井眼越远权重越小，涉及的计算公式有：

距离井眼距离的计算公式：

$R_i=\sqrt{{(x_0-x_i)}^2+{(y_0-y_i)}^2}$

其中，R_i为地震道距离井眼的距离；x₀为0井的井眼位置横坐标；y₀为0井的井眼位置纵坐标；x_i为0井旁第i个地震道位置的横坐标；y_i为0井旁第i个地震道位置的纵坐标；

距离加权平均的井点地震属性计算公式：

$A_0=\frac{\underset{i=1}{\overset{j}{\Σ}}{W}_i{R}_i}{\underset{i=1}{\overset{j}{\Σ}}{R}_i}$

其中，A₀为0井井眼处距离加权平均后的地震属性值；W_i为0井井旁第i个地震道的属性值，j为井旁道范围(1～10)。

图4为图1中的步骤S105的具体流程图，由图4可知，该步骤具体包括：

S401：根据所述的沉积微相以及所述的井点地震属性确定二维散点交汇关系。此处的二维散点交汇关系在具体的实际使用中，诸如为砂岩厚度与均方根振幅属性、弧长属性的数学定量关系为：y＝2.3731×10^-4x₁-3.35×10^-5x₂，y为计算的砂岩厚度，x₁为均方根振幅属性，x₂为弧长属性，相关系数为0.847。

S402：根据所述的二维散点交汇关系从所述的井点地震属性中选出将所述的沉积微相区分开的井点地震属性；

S403：将所述的沉积微相区分开的井点地震属性，称为特征地震属性；

S404：根据所述的二维散点交汇关系建立所述的沉积微相与特征地震属性的数学映射关系。

优选的特征地震属性主要通过二维散点交汇图的线性拟合、区间分布法进行，每个沉积微相对应一个地震属性值，获得微相类别与地震属性的数学映射关系。

优选的特征地震属性在具体的实例中主要包括优选的单一地震属性或优选的融合地震属性两个方面。其中，优先选择单一地震属性与沉积微相的数学映射关系，若不能满足研究需求，则进行融合地震属性与沉积微相的数学映射关系研究。

优选的单一地震属性的具体做法是：

1)利用所述的沉积微相与所述井眼处的地震属性计算结果做二维散点交汇；

2)根据微相的不同交汇图中的数据点进行类型区分，分析是否能将沉积微相清晰区分，若不能，则选择其它地震属性如此分析；

优选的融合地震属性的具体做法是：

1)根据单一地震属性与沉积微相的交汇趋势关系，从地震属性计算方法角度考虑各个地震属性之间的关系，不选择相似度较大的地震属性进行融合，以免带来冗余信息的影响，融合方法采用多元非线性回归的方法。

2)根据微相的不同将交汇中的数据点进行类型区分，分析是否能将沉积微相清晰区分，若不能，则更换地震属性类别继续进行属性融合。

图5为图1中的步骤S106的具体流程图，由图6可知，该步骤具体包括：

S501：根据所述的地震资料确定三维地震数据体；

S502：将所述的三维地震数据体转换到Wheeler域，得到Wheeler域数据体，Wheeler域数据体异于常规叠后三维地震数据体，是常规叠后三维地震数据体经过wheeler变换后的wheeler域地震数据体。

S503：根据所述的特征地震属性在所述的Wheeler域数据体上提取三维特征属性体。在Wheeler域数据体上提取地震属性的计算方法与常规叠后三维地震数据体的属性计算方法一样，根据特征地震属性在Wheeler域数据体提取的地震属性成为三维特征属性体。S504：根据提取的三维特征属性体生成地层切片。该步骤主要包括在提取的三维地震属性体上通过等比例内插法、非线性插值法提取地层切片。

以上本发明提供的一项基于地层切片的小尺度沉积相进行储层预测的方法包括的技术主要有两种:井点处各沉积微相层段的地震属性提取、地层切片-沉积相快速计算机自动转化。已有的资料中没有提及井点处各沉积微相层段的地震属性提取技术，未见利用地层切片自动进行沉积相解释的方法及相关技术报道。

图6为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测系统的结构框图，由图6可知，该系统具体包括：

采集装置100，用于采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；

沉积微相解释划分装置200，用于根据所述的测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相。该装置在进行沉积微相的划分时，主要应用到测井曲线中的自然伽马曲线、自然电位曲线、微电极电阻率曲线的幅度特征、形态特征、光滑度三个参数，结合岩心资料中的岩心描述、钻井取心等地质对沉积微相进行划分。得到各个类型的沉积微相可按照微相类别依次定名为1、2、3、4、……。

地震属性确定装置300，用于根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性。地震属性确定装置300的具体结构如图7所示。

井点地震属性提取装置400，用于根据所述的沉积微相以及所述的地震资料提取井点地震属性。井点地震属性提取装置400的具体结构如图8所示。

数学映射关系确定装置500，用于根据所述的沉积微相以及井点地震属性确定数学映射关系。数学映射关系确定装置500的具体结构如图9所示。

地层切片生成装置600，用于根据所述的地震资料生成地层切片。地层切片生成装置600的具体结构如图10所示。

沉积分布图转换装置700，用于根据所述的数学映射关系将所述的地层切片转换为沉积相平面分布图。将沉积微相与特征地震属性的数学映射关系应用到地层切片上，并以该地层切片的地震属性值为约束条件，从井点出发，外推到无井区，实现地层切片的自动沉积相解释。

储层预测装置800，用于根据所述的沉积相平面分布图进行储层预测。

图7为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测系统中地震属性确定装置的结构框图，由图7可知，地震属性确定装置具体包括：

曲线提取模块201，用于从所述的测井曲线中提取声波曲线、密度曲线；

合成记录制作模块202，用于根据所述的声波曲线、密度曲线以及地震资料制作当前钻井的精细合成记录；

地震反射特征分析模块203，用于根据所述的精细合成记录以及地震资料分析所述目的层的地震反射特征，所述地震反射特征包括振幅、频率、相位；

地震属性确定模块204，用于根据所述的地震反射特征确定出区分所述的沉积微相的地震属性。

该装置主要分析目的层的各种沉积微相的振幅、频率等基本特征，初步确定出区分沉积微相的地震属性。所述的地震属性属于体属性范畴，不包括时窗属性。

图8为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测系统中井点地震属性提取装置的结构框图，由图8可知，井点地震属性提取装置具体包括：

深度获取模块301，用于获取所述的沉积微相对应的顶深度、底深度；

合成记录获取模块302，用于获取当前钻井的精细合成记录，此处的精细合成记录即为合成记录制作模块202中确定出的合成记录。

顶时间范围转换模块303，用于根据所述的精细合成记录将所述沉积微相的顶深度转换为地震时间域的顶时间范围；

底时间范围转换模块304，用于根据所述的精细合成记录将所述沉积微相的底深度转换为地震时间域的底时间范围；

井点地震属性提取模块305，用于根据地震资料在所述的顶时间范围、所述的底时间范围内提取井点地震属性。

在具体的实施方式中，诸如在井旁道范围依据默认提取井眼为中心的7个地震道，也可在井旁道范围依据实际地层倾斜程度而定，地层越倾斜，井旁道范围越小，地层越平缓，井旁道范围可以适当变大，一般为1～10；若微相纵向时间范围较小，不能满足提取某种地震属性的需要，则默认为提取该地震属性需要的最小时间范围。井点地震属性的提取是在井旁道区域内，根据距井眼的距离加权平均的方法，距离井眼越远权重越小，涉及的计算公式有：

距离井眼距离的计算公式：

$R_i=\sqrt{{(x_0-x_i)}^2+{(y_0-y_i)}^2}$

其中，R_i为地震道距离井眼的距离；x₀为0井的井眼位置横坐标；y₀为0井的井眼位置纵坐标；x_i为0井旁第i个地震道位置的横坐标；y_i为0井旁第i个地震道位置的纵坐标；

距离加权平均的井点地震属性计算公式：

$A_0=\frac{\underset{i=1}{\overset{j}{\Σ}}{W}_i{R}_i}{\underset{i=1}{\overset{j}{\Σ}}{R}_i}$

其中，A₀为0井井眼处距离加权平均后的地震属性值；W_i为0井井旁第i个地震道的属性值，j为井旁道范围(1～10)。

图9为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测的系统中数学映射关系确定装置的结构框图，由图9可知，数学映射关系确定装置具体包括：

交汇关系确定模块401，用于根据所述的沉积微相以及所述的井点地震属性确定二维散点交汇关系。此处的二维散点交汇关系在具体的实际使用中，诸如为砂岩厚度与均方根振幅属性、弧长属性的数学定量关系为：y＝2.3731×10^-4x₁-3.35×10^-5x₂，y为计算的砂岩厚度，x₁为均方根振幅属性，x₂为弧长属性，相关系数为0.847。

特征地震属性确定模块402，用于根据所述的二维散点交汇关系从所述的井点地震属性中选出将所述的沉积微相区分开的井点地震属性，称为特征地震属性；

数学映射关系确定模块403，用于根据所述的二维散点交汇关系建立所述的沉积微相与特征地震属性的数学映射关系。

优选的特征地震属性主要通过二维散点交汇图的线性拟合、区间分布法进行，每个沉积微相对应一个地震属性值，获得微相类别与地震属性的数学映射关系。

优选的特征地震属性在具体的实例中主要包括优选的单一地震属性或优选的融合地震属性两个方面。其中，优先选择单一地震属性与沉积微相的数学映射关系，若不能满足研究需求，则进行融合地震属性与沉积微相的数学映射关系研究。

优选的单一地震属性的具体做法是：

1)利用所述的沉积微相与所述井眼处的地震属性计算结果做二维散点交汇；

2)根据微相的不同交汇图中的数据点进行类型区分，分析是否能将沉积微相清晰区分，若不能，则选择其它地震属性如此分析；

优选的融合地震属性的具体做法是：

1)根据单一地震属性与沉积微相的交汇趋势关系，从地震属性计算方法角度考虑各个地震属性之间的关系，不选择相似度较大的地震属性进行融合，以免带来冗余信息的影响，融合方法采用多元非线性回归的方法。

2)根据微相的不同将交汇中的数据点进行类型区分，分析是否能将沉积微相清晰区分，若不能，则更换地震属性类别继续进行属性融合。

图10为本发明实施例提供的一种基于地层切片的小尺度沉积相进行储层预测的系统中地层切片生成装置的结构框图，由图10可知，地层切片生成装置具体包括：

三维地震数据体确定模块501，用于根据所述的地震资料确定三维地震数据体；

地层切片数据体确定模块502，用于将所述的三维地震数据体转换到Wheeler域，得到Wheeler域数据体，Wheeler域数据体异于常规叠后三维地震数据体，是常规叠后三维地震数据体经过wheeler变换后的wheeler域地震数据体。

三维属性体提取模块503，用于根据所述的特征地震属性在所述的Wheeler域数据体上提取三维特征属性体。在Wheeler域数据体上提取地震属性的计算方法与常规叠后三维地震数据体的属性计算方法一样，根据特征地震属性在Wheeler域数据体上提取的地震属性成为三维特征属性体。地层切片生成模块504，用于根据提取的三维特征属性体生成地层切片。该模块主要包括在提取的三维地震属性体上通过等比例内插法、非线性插值法提取地层切片。

以上本发明提供的一项基于地层切片的小尺度沉积相进行储层预测的系统包括的技术主要有两种:井点处各沉积微相层段的地震属性提取、地层切片-沉积相快速计算机自动转化。已有的资料中没有提及井点处各沉积微相层段的地震属性提取技术，未见利用地层切片自动进行沉积相解释的方法及相关技术报道。

下面结合具体的实施例，详细介绍本发明的技术方案。本发明提供的一种基于地层切片的小尺度沉积相进行储层预测方案主要包括：

S1、目标区内重点井的沉积微相划分。

利用常规测井曲线(自然伽马、自然电位、微电极电阻率曲线)、录井数据及岩心资料，根据曲线形态特征、光滑度、幅度特征进行目的层的沉积微相单井解释。如图11所示的qn6井沙一中重力流水道测井相特征，以及图12所示的ch52井沙一中深湖测井相特征。

S2、制作高精度地震合成记录，进行精细井震标定，将单井沉积微相解释结果与地震资料做到时深相互匹配，初步分析地震反射振幅、频率等特征，确定要提取的属性类型。如图13所示的qn6井精细合成记录示意图，以及图14所示的ch52井精细合成记录示意图。

S3、在地震资料中提取各微相层段对应的井点地震属性，多为体属性类型，采用交汇图的方式优选出能够将各沉积微相区分开的特征地震属性类型，建立沉积微相与特征地震属性的数学映射关系。本实施例通过多个交汇图的对比，优选出瞬时振幅属性，如图15所示的QN地区井点平均振幅属性与伽马、沉积微相交汇图示意图。

S4、将原始地震数据转为wheeler域，然后依据wheeler数据体提取选择好的特征地震属性，本实施例为瞬时振幅属性，最后在此地震属性数据体制作地层切片并选择出合适的地层切片，如图16所示的wheeler域沙一中某振幅地层切片示意图。

S5、对选择好的地层切片进行初步预处理，主要包括断层处数据校正、数据滤波，使得地层切片图像特征更加能够体现沉积现象。如图17所示的数据处理后的地层切片示意图；

S6、在处理后的地层切片上，利用沉积微相与特征地震属性的数学映射关系，快速将此地层切片转化为沉积相平面展布图。如图18所示的自动沉积相解释结果示意图。由图18可知，重力流水道填充区域为储层发育区，是油气勘探重点区，而深湖填充区域主要为泥岩发育区，不发育储层。

综上所述，本发明提供了一种基于地层切片的小尺度沉积相进行储层预测的方法及系统，借助井点沉积相的可靠解释成果和等时地层切片清晰反映的沉积特征，在小尺度研究范围内，充分挖掘沉积特征与地震属性之间的内在联系，实现了井点控制下的地层切片数据定量、快速、有效、自动转为沉积相平面图，进而根据沉积相平面图进行储层预测。

与前人研究的沉积相解释的技术方法对比，本发明具备以下几个显著特点：1)沉积相研究精度显著提高，主要依据提取的地层切片采样率而定；2)计算机自动实现了依据地层切片进行沉积相自动解释，极大提高了工作效率；3)井点地震属性快速提取技术极大提高了测井-地震-地质相结合的定量化研究程度；4)该技术方法原理简单、操作方便、技术上易于实现。利用该发明对QN地区沙一中下部的沉积相展布特征进行自动解释，解释成果符合地质规律，通过与人工解释成果的对比，统计发现，二者的相关程度高达95％，能够满足科研生产需要。

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

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

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

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

Claims (6)

1.一种基于地层切片的小尺度沉积相进行储层预测的方法，其特征是，所述的方法包括：

采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；

根据所述的测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相；

根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性，所述的地震属性属于体属性范畴，不包括时窗属性；

根据所述的沉积微相以及所述的地震资料提取井点地震属性；

根据所述的沉积微相以及井点地震属性确定数学映射关系；

根据所述的地震资料生成地层切片；

根据所述的数学映射关系将所述的地层切片转换为沉积相平面分布图；

根据所述的沉积相平面分布图进行储层预测；

其中，根据所述的沉积微相以及井点地震属性确定数学映射关系包括：根据所述的沉积微相以及所述的井点地震属性确定二维散点交汇关系；根据所述的二维散点交汇关系从所述的井点地震属性中选出将所述的沉积微相区分开的井点地震属性，称为特征地震属性；根据所述的二维散点交汇关系建立所述的沉积微相与特征地震属性的数学映射关系；

根据所述的地震资料生成地层切片包括：根据所述的地震资料确定三维地震数据体；将所述的三维地震数据体转换到Wheeler域，得到Wheeler域数据体；根据所述的特征地震属性在所述的Wheeler域数据体上提取三维特征属性体；根据提取的三维特征属性体生成地层切片。

2.根据权利要求1所述的方法，其特征是，根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性具体包括：

从所述的测井曲线中提取声波曲线、密度曲线；

根据所述的声波曲线、密度曲线以及地震资料制作当前钻井的精细合成记录；

根据所述的精细合成记录以及地震资料分析所述目的层的地震反射特征，所述地震反射特征包括振幅、频率、相位；

根据所述的地震反射特征确定出区分所述的沉积微相的地震属性。

3.根据权利要求2所述的方法，其特征是，根据所述的沉积微相以及所述的地震资料提取井点地震属性具体包括：

获取所述的沉积微相对应的顶深度、底深度；

获取当前钻井的精细合成记录；

根据所述的精细合成记录将所述沉积微相的顶深度转换为地震时间域的顶时间范围；

根据所述的精细合成记录将所述沉积微相的底深度转换为地震时间域的底时间范围；

根据地震资料在所述的顶时间范围、所述的底时间范围内提取井点地震属性。

4.一种基于地层切片的小尺度沉积相进行储层预测的系统，其特征是，所述的系统包括：

采集装置，用于采集当前钻井的测井曲线、录井剖面、岩心资料以及地震资料；

沉积微相解释划分装置，用于根据所述的测井曲线、录井剖面以及岩心资料对目的层的沉积微相进行解释、划分，得到各个类型的沉积微相；

地震属性确定装置，用于根据所述的测井曲线以及地震资料确定出区分所述的沉积微相的地震属性，所述的地震属性属于体属性范畴，不包括时窗属性；

井点地震属性提取装置，用于根据所述的沉积微相以及所述的地震资料提取井点地震属性；

数学映射关系确定装置，用于根据所述的沉积微相以及井点地震属性确定数学映射关系；

地层切片生成装置，用于根据所述的地震资料生成地层切片；

沉积分布图转换装置，用于根据所述的数学映射关系将所述的地层切片转换为沉积相平面分布图；

储层预测装置，用于根据所述的沉积相平面分布图进行储层预测；

其中，所述的数学映射关系确定装置具体包括：交汇关系确定模块，用于根据所述的沉积微相以及所述的井点地震属性确定二维散点交汇关系；特征地震属性确定模块，用于根据所述的二维散点交汇关系从所述的井点地震属性中选出将所述的沉积微相区分开的井点地震属性，称为特征地震属性；数学映射关系确定模块，用于根据所述的二维散点交汇关系建立所述的沉积微相与特征地震属性的数学映射关系；

所述的地层切片生成装置具体包括：三维地震数据体确定模块，用于根据所述的地震资料确定三维地震数据体；地层切片数据体确定模块，用于将所述的三维地震数据体转换到Wheeler域，得到Wheeler域数据体；三维属性体提取模块，用于根据所述的特征地震属性在所述的Wheeler域数据体上提取三维特征属性体；地层切片生成模块，用于根据提取的三维特征属性体生成地层切片。

5.根据权利要求4所述的系统，其特征是，所述的地震属性确定装置具体包括：

曲线提取模块，用于从所述的测井曲线中提取声波曲线、密度曲线；

合成记录制作模块，用于根据所述的声波曲线、密度曲线以及地震资料制作当前钻井的精细合成记录；

地震反射特征分析模块，用于根据所述的精细合成记录以及地震资料分析所述目的层的地震反射特征，所述地震反射特征包括振幅、频率、相位；

地震属性确定模块，用于根据所述的地震反射特征确定出区分所述的沉积微相的地震属性。

6.根据权利要求5所述的系统，其特征是，所述的井点地震属性提取装置具体包括：

深度获取模块，用于获取所述的沉积微相对应的顶深度、底深度；

合成记录获取模块，用于获取当前钻井的精细合成记录；

顶时间范围转换模块，用于根据所述的精细合成记录将所述沉积微相的顶深度转换为地震时间域的顶时间范围；

底时间范围转换模块，用于根据所述的精细合成记录将所述沉积微相的底深度转换为地震时间域的底时间范围；

井点地震属性提取模块，用于根据地震资料在所述的顶时间范围、所述的底时间范围内提取井点地震属性。