一种薄储层的预测方法及装置与流程

本申请属于石油勘探技术领域，尤其涉及一种薄储层的预测方法及装置。

背景技术：

薄储层是我国陆相盆地岩性地层油气藏的主要储集空间，储层厚度一般小于四分之一地震子波的波长，邻层干涉是常规地震资料不能识别薄储层的根本原因。

现有技术中，薄储层预测主要有两类技术，一类是：利用反演技术垂向分辨薄储层与邻层预测薄储层，由于反演技术不能直接提高分辨率，需要借助井震关系统计和测井约束进行外推，因此，该方法多解性较强，薄层预测精度有限。另外一种技术是：通过敏感地震属性或地层切片序列分析横向预测薄储层，其中敏感地震属性通过井震交汇分析寻找最佳地震属性，是一种统计关系，具有一定的局限性；通过分析地层切片序列，可以突出目标储层，但无法避开邻层干涉，在薄互层复杂情况下，邻层干涉甚至会掩盖目标储层，导致目标储层无法识别。上述方法都会导致薄储层的预测精度不高，因此，提供一种提高薄储层的预测精度的技术方案，是亟待解决的技术问题。

技术实现要素：

本申请目的在于提供一种薄储层的预测方法及装置，提高了薄储层预测的准确性。

一方面本申请提供了一种薄储层的预测方法，包括：

根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据；

根据所述目标井的时频谱、所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率；

根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据；

根据所述叠后地震数据构建地层时间模型；

沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据；

根据所述降频地层切片数据，预测薄储层的平面分布。

进一步地，所述方法的另一个实施例中，所述根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据，包括：

根据所述叠后地震数据、所述时深关系曲线、所述储层特征曲线，确定所述叠后地震数据的相位；

根据所述叠后地震数据的相位，利用希尔伯特变换生成所述-90度相位地震数据。

进一步地，所述方法的另一个实施例中，所述根据所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率，包括：

在所述目标井的时频谱中，进行频率扫描；

根据所述时深关系曲线、所述特征曲线，获取不同频率下，所述薄储层所述-90度相位地震数据中的位置；

将所述薄储层在所述-90度相位地震数据中的最大振幅位置处时的频率，作为所述目标干涉频率。

进一步地，所述方法的另一个实施例中，所述根据所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率，包括：

在所述目标井的时频谱中，进行频率扫描；

根据所述时深关系曲线、所述特征曲线，获取不同频率下，所述薄储层所述-90度相位地震数据中的位置；

将所述薄储层在所述-90度相位地震数据中的最大振幅位置处时的频率，作为所述目标干涉频率。

进一步地，所述方法的另一个实施例中，所述时频谱包括：小波时频谱、s变换时频谱、傅里叶变换时频谱、子波分解时频谱中的任意一个。

进一步地，所述方法的另一个实施例中，所述根据所述叠后地震数据构建地层时间模型，包括：

根据所述叠后地震数据，选择参考层，所述参考层包括地质时间相同的同相轴层位；

计算所述参考层之间的平均距离或最大距离；

根据所述平均距离或最大距离，计算所述参考层之间的内插地层切片个数；

根据所述内插地层切片个数和所述参考层，通过线性内插函数，构建所述地层时间模型。

进一步地，所述方法的另一个实施例中，所述根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，包括：

根据所述目标干涉频率，采用小波变换、广义s变换、傅里叶变换、子波分解中的任意一项技术对所述-90度相位地震数据进行时频信号分析。

另一方面，本申请提供了一种薄储层的预测装置，包括：

相位地震数据获取模块，用于根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据；

目标干涉频率确定模块，用于根据所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率；

时频信号分析模块，用于根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据；

地层时间模型构建模块，用于根据所述叠后地震数据构建地层时间模型；

降频处理模块，用于沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据；

薄储层预测模块，用于根据所述降频地层切片数据，预测薄储层的平面分布。

再一方面，本申请还提供了一种薄储层的预测装置，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据；

根据所述目标井的时频谱、所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率；

根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据；

根据所述叠后地震数据构建地层时间模型；

沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据；

根据所述降频地层切片数据，预测薄储层的平面分布。

又一方面，本申请还提供了一种地层品质因子的自动估算系统，包括：处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述薄储层的预测方法。

本申请提供的薄储层的预测方法及装置，通过-90度相位地震数据、目标井的时深关系曲线和储层特征曲线确定薄储层的位置；通过小波时频谱、目标井的时深关系曲线和薄储层的位置确定目标干涉频率；根据确定的目标干涉频率通过小波变换等方法生成目标分频数据；通过目标分频数据和参考层生成降频地层切片数据。降频地层切片数据清晰的刻画了薄储层的分布范围，解决了由于薄储层干涉引起的切片范围预测结果不准的问题，提高了薄储层预测的精度。提出了一种与传统的垂向分辨薄层完全不同的研究思路，提供了一种新的薄储层预测技术，是地震沉积学在薄层预测方面的重大突破。

附图说明

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

图1是本申请提供的一种薄储层的预测方法一个实施例的方法流程示意图；

图2(a)-2(b)是本申请一个实施例中不同储层特征曲线对应的波形相位变化的示意图；

图3是本申请一个实施例中目标井的合成地震记录示意图；

图4是本申请一个实施例中目标井的小波时频谱的示意图；

图5是现有技术中薄储层的地层切片示意图；

图6是本申请实施例中薄储层的降频地层切片示意图；

图7是本申请又一个实施例中薄储层预测流程示意图；

图8是本申请提供的薄储层的预测装置一个实施例的模块结构示意图；

图9是本申请提供的一种薄储层的预测系统的实施例的模块结构示意图。

具体实施方式

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

薄储层是我国陆相盆地岩性地层油气藏的主要储集空间，薄储层平面分布的预测是油气勘探开发的一项重要工作。薄储层指厚度小于地震分辨率的储层，薄储层叠加可以形成薄互层，且叠加的薄储层之间的厚度也小于地震分辨率。

本申请实施例以地层切片和小波分频为基础，通过降低频率来减弱邻层干涉，在此基础上再次生成地层切片。通过垂向降低邻层干涉、横向切片分析的薄层预测思路，是一种与传统的垂向分辨薄层完全不同的研究思路，可有效解决大部分薄互储层预测问题，且多解性弱，具有广阔的应用前景。

具体地，图1是本申请提供的一种薄储层的预测方法一个实施例的方法流程示意图，本申请提供的薄储层的预测方法包括：

s1、根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据。

具体地，本申请实施例可以根据叠后地震数据和测井曲线(如声波测井曲线)，获得目标井的地震数据的时深关系曲线。叠后地震数据可以表示对采集到的地震数据进行叠加处理后的地震数据，时深关系曲线可以包括目标井内记录的地震数据的时间和地下实际深度之间的关系。目标井可以表示钻穿目的层段，且钻遇目的薄储层的井，可以根据各个井的记录数据、测井曲线等确定该井是否属于目标井。地震数据包括叠后地震数据是以时间来记录，测井曲线以深度来记录，可以利用速度曲线(速度曲线是声波测井曲线的倒数)和地震子波褶积，将声波测井曲线转换成合成地震记录，将合成地震记录与实际地震数据对比，二者接近一致时，合成地震记录对应的测井曲线的深度与实际地震数据的时间具有一一对应关系，该关系就是时深关系，具体可以以曲线的形式表示。

储层特征曲线可以包括声波、密度、速度或波阻抗曲线等。利用储层特征曲线、时深关系曲线、叠后地震数据可以获得-90度相位地震数据。具体可以通过对比储层特征曲线、时深关系曲线，获得叠后地震数据的相位，通过对叠后地震数据进行相移，可以获得-90度相位地震数据。-90度相位地震数据也可以称为270度相位地震数据，表示地震数据的波形的相位是270度或者-90度。

本申请一个实施例中，所述根据叠后地震数据、储层特征曲线、地震记录的时深关系曲线，获得-90度相位地震数据，可以包括：

根据所述叠后地震数据、所述时深关系曲线、所述储层特征曲线，确定所述叠后地震数据的相位；

根据所述叠后地震数据的相位，利用希尔伯特变换生成所述-90度相位地震数据。

具体地，图2(a)-(b)是本申请一个实施例中不同储层特征曲线对应的波形相位变化的示意图，图2(a)可以表示储层特征曲线一对应的波形相位变化示意图，图2(b)可以表示储层特征曲线二对应的波形相位变化示意图。图2(a)-2(b)中左侧的特征曲线可以表示储层特征曲线，右侧为测井曲线生成的合成地震记录，其中横轴代表相位。可以看到，在特征曲线固定的情况下，不同相位得到的地震记录形态是不同的。

图2(a)中在特征曲线中界面以上的部分可以称为顶部，在界面以下的部分可以称为底部。当右侧的地震记录的曲线特征为顶部数值小，底部数值大且数值延续接近一个波长时(如图2(a))，若叠后地震数据对应为对称正极性波形，说明叠后地震数据为零相位(如图2(a)中左侧第一道波形曲线)；若叠后地震数据对应为对称反极性波形，说明叠后地震数据为180相位(如图2(a)中左侧第七道波形曲线)；若叠后地震数据对应为反对称正极性波形，说明叠后地震数据为90相位；若叠后地震数据对应为反对称负极性波形，说明叠后地震数据为-90相位。

当曲线特征为顶部数据小、底部数据大切片数值延续接近一个波长时，情况与上面相反。当无法找着这种曲线特征时，也可以寻找单一薄层情况下对应的曲线特征和叠后地震数据波形特征来确定数据体相位(如图2(b))，若叠后地震数据对应为对称正极性波形，说明叠后地震数据为270(也可以称为-90)相位；若叠后地震数据对应为对称反极性波形，说明叠后地震数据为90相位；若叠后地震数据对应为反对称正极性波形，说明叠后地震数据为零相位；若叠后地震数据对应为反对称负极性波形，说明叠后地震数据为180相位。

通过对比测井曲线生成的地震记录(即叠后地震数据的波形曲线)和储层特征曲线，可以确定叠后地震数据的相位。例如，图2(a)中左侧第1道波形是以中间水平线为基准上下对称的，黑色填充位置表示极性为正，可以称之为对称正极性波形，该道地震数据对应的相位为零相位。

获得叠后地震数据的相位后，可以通过相移例如：采用希尔伯特变换将叠后地震数据进行相移获得-90度相位地震数据。信号经希尔伯特变换后，在频域各频率分量的幅度保持不变，但相位将出现90°相移，采用希尔伯特变换可以确保叠后地震数据的其他能量保持不变，仅仅发生相位的移动。

图3是本申请一个实施例中目标井的合成地震记录示意图，如图3所示，图中左侧纵轴可以表示深度，右侧纵轴可以表示时间，well1可以表示目标井1，gr曲线可以表示自然伽玛曲线，图中oil表示油层，t3b、t2k2、t2k1、t1b1均可以表示地层。sp测井曲线(表示自然电位曲线)显示薄储层位于2210m左右，薄储层的厚度约6m，在其顶部和底部还各有一个约8m厚的储层，薄储层与顶部储层的夹层厚度约6m，与底部储层的夹层厚度约12m。薄储层远小于地震分辨率，且顶底储层都会对薄储层形成干涉，造成薄储层平面分布难以预测。通过对比速度曲线、叠后地震数据和合成地震记录确定叠后地震数据接近零相位，极性为反极性，相移角度为90度，可以通过希尔伯特变换生成-90度相位地震数据体。利用希尔伯特变换生成-90度相位地震数据体，方法简单，可以为后续薄储层的预测提供准确的数据基础。

s2、根据所述目标井的时频谱、所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率。

具体地，可以通过对比储层特征曲线和时深关系曲线，确定薄储层在-90度相位地震数据中的位置，这里可以确定是薄储层在-90度相位地震数据中的深度位置。可以根据在时频谱中薄储层在-90度相位地震数据中的位置，确定出薄储层是否存在邻层干涉，进一步确定出邻层干涉较小时时频谱中的频率，确定为目标干涉频率。

本申请一个实施例中，所述根据所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率，可以包括：

在所述目标井的时频谱中，进行频率扫描；

根据所述时深关系曲线、所述特征曲线，获取不同频率下，所述薄储层所述-90度相位地震数据中的位置；

将所述薄储层在所述-90度相位地震数据中的最大振幅位置处时的频率，作为所述目标干涉频率。

图4是本申请一个实施例中目标井的小波时频谱的示意图，图中well1可以表示目标井1，gr曲线可以表示自然伽玛曲线，sp测井曲线也可以表示自然电位曲线。如图4所示，叠后地震数据体主频约40hz，将目标井(即图中well1)按照时深关系曲线投影于40hz位置，通过-90度相位地震数据、目标井的时深关系曲线和储层特征曲线确定薄储层在-90度相位地震数据接近零点位置。薄储层位于零点位置，表明此时邻层干涉较大。计算目标井1的小波时频谱，即计算振幅随频率的变化，通过振幅和频率的变化分析发现当频率为20hz时，薄储层位于-90度相位地震数据中最大振幅位置处，表明在当前频率下，邻层干涉最小。因此，可以确定当前目标井的地震数据的目标干涉频率为20hz。

储层特征曲线是指对储层最为敏感的曲线，gr曲线和sp曲线是储层特征曲线其中两类曲线。gr曲线和sp曲线为深度域曲线，地震数据为时间域数据，二者都是地下地质情况的响应，本申请实施例中通过建立时间与深度的一一对应关系即时深关系曲线，将不同域的数据匹配起来。

本申请一个实施例中，所述时频谱可以包括：小波时频谱、s变换时频谱、傅里叶变换时频谱、子波分解时频谱中的任意一个。s变换是小波变换和短时傅里叶变换的继承和发展，是一种时频分析处理方法。

具体可以根据实际需要选择合适的时频谱进行目标干涉频率的确定，本申请实施例不作具体限定。

在本申请实施例实施过程中发现，在没有干涉的情况下，薄储层的中心应对应于-90度相位地震数据体的最大振幅处，当存在邻层干涉的情况下，薄层中心会偏离最大振幅，若能够使薄储层中心接近最大振幅，表明此时的邻层干涉最弱。本申请实施例通过降低频率可以减弱邻层干涉，而本申请实施例中的目标干涉频率为时频谱(如小波时频谱)上薄储层以最大振幅为中心对应的频率，此时邻层干涉达到最小。

目标干涉频率可以表明对地震数据进行多少赫兹的分频处理，是薄储层平面分布预测的理论基础。本申请实施例结合时频谱如小波时频谱、储层特征曲线、-90度相位地震数据等，确定出薄储层对应的目标干涉频率。

在本申请一个实施例中，所述根据所述储层特征曲线、所述地震记录的时深关系曲线、所述-90度相位地震数据，确定目标干涉频率，可以包括：

所述目标井包括测试井，获取预设数量口所述测试井的时频谱，在所述测试井的时频谱中进行频率扫描，获得所述测试井对应的目标干涉频率；

将所述预设数量口所述测试井对应的目标干涉频率的平均值，作为所述目标干涉频率。

具体地，实际应用中由于不同井的薄互层速度、厚度等参数都不相同，不同井点会得到不同的目标干涉频率，为保证所有井点的邻层干涉都达到最小，本申请实施例中可以选择一口钻穿薄储层的目标井来确定目标干涉频率。还可以选择多口井，采用上述实施例的方法，获得不同井对应的目标干涉频率。将不同井对应的目标干涉频率的平均值作为最终的整个地震数据的目标干涉频率。当然，此处在选择的井可以不是钻穿目的层段，且钻遇目的薄储层的目标井，可以是一般的钻井，本申请实施例中可以称为测试井。另外，选择测试井的数量也可以根据实际需要进行选择，本申请实施例不作具体限定。

本申请实施例利用目标井或多口测试井的时频谱，结合储层特征曲线、-90度相位地震数据等，通过频率扫描，可以确定出对薄储层干涉最小的目标干涉频率，为后续薄储层的预测提供了准确的数据基础。

s3、根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据。

具体地，可以以-90度相位地震数据为基础，根据确定的目标干涉频率通过小波变换方法进行时频信号分析，生成目标分频数据体。本申请实施例可以选择小波变换生成目标分频数据，还可以选择广义s变换、傅里叶变换、子波分解等技术对-90度相位地震数据进行时频信号分析，获得目标分频数据体，本申请实施例不作具体限定。

s4、根据所述叠后地震数据构建地层时间模型。

地层时间模型可以包括三维地震数据以及各个数据点拉平前对应的采集时间。

本申请一个实施例中，所述根据所述叠后地震数据构建地层时间模型，可以包括：

s40、根据所述叠后地震数据，选择参考层，所述参考层包括地质时间相同的同相轴层位；

s42、计算所述参考层之间的平均距离或最大距离；

s44、根据所述平均距离或最大距离，计算所述参考层之间的内插地层切片个数；

s46、根据所述内插地层切片个数和所述参考层，通过线性内插函数，构建所述地层时间模型。

具体地，可以以叠后地震数据为基础，选择地质时间相同的标志地震同相轴作为参考层。通常情况下，可选择最大洪(湖)泛面作为参考层，为了减少参考层的错误，可以尽可能多选择参考层。设定地层模型的起始相对地质时间和采样间隔，根据选择的参考层计算参考层间的平均距离或最大距离。相对地质时间就是地质体形成的相对年代，不是真实的年代，采样间隔表示采集地层模型的时间值。根据平均距离或最大距离及采样间隔计算参考层之间的内插地层切片个数，根据内插地层切片个数和参考层通过线性内插函数建立地层时间模型。根据未知函数在某区间内若干点的函数值，作出在该若干点的函数值与值相等的特定函数来近似原函数，进而可用此特定函数算出该区间内其他各点的原函数的近似值，这种方法可以称为内插法。按特定函数的性质分，有线性内插、非线性内插，本申请实施例中采用线性内插法构建地层时间模型，方法简单，操作方便。地层时间模型的建立为后续地层切片提供了数据基础。

s5、沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据。

具体地，沿地层时间模型中的每个地层层位从目标分频地震数据中提取振幅，形成一个降频地层切片数据。降频地层切片提取的仍是振幅数据，数据是分频数据中的振幅值组成的三维数据体。本申请实施例通过时频谱寻找目标层的最小干涉频率，而最小干涉频率一般都小于地震数据的主频，因此，这里属于降频。理论上降频会降低分辨率，更加不利于识别薄层，而本申请却可以利用降频来识别薄层。

s6、根据所述降频地层切片数据，预测薄储层的平面分布。

具体地，对降频地层切片数据进行分析，降频切片地层数据中，最大振幅值所处的范围即可以确定为薄储层的平面分布范围。

图5是现有技术中薄储层的地层切片示意图，如图5所示，图中对每口井的测井曲线进行了标记，图中1km可以表示比例尺，曲线中心的虚线表示薄储层所在深度。本实例中sp曲线对砂体变化较为敏感，从各井的sp曲线可以看出，仅有少数井与切片结果吻合，多数井如well5、well6、well7、well8等井点的钻井信息与切片结果不吻合。由于邻层干涉，传统地层切片更多的反映的是邻层的平面分布特征。

图6是本申请实施例中薄储层的降频地层切片示意图，如图6所示，图中1km可以表示比例尺，从各井的sp曲线可以看出，well1、well2、well3钻遇了薄储层，其它井未钻遇。降频地层切片预测结果与所有钻井结果相吻合，清晰了展示了目标储层的分布范围。

需要说明的是，本申请实施例中各步骤没有固定的执行顺序，具体步骤的执行顺序可以根据实际需要进行调整，例如：步骤s4可以与步骤s1或s2或s3同步进行，本申请实施例不作具体限定。

本申请实施例提供的薄储层的预测方法，通过-90度相位地震数据、目标井的时深关系曲线和储层特征曲线确定薄储层的位置；通过小波时频谱、目标井的时深关系曲线和薄储层的位置确定目标干涉频率；根据确定的目标干涉频率通过小波变换等方法生成目标分频数据；通过目标分频数据和参考层生成降频地层切片数据。降频地层切片数据清晰的刻画了薄储层的分布范围，解决了由于薄储层干涉引起的切片范围预测结果不准的问题，提高了薄储层预测的精度。提出了一种与传统的垂向分辨薄层完全不同的研究思路，提供了一种新的薄储层预测技术，是地震沉积学在薄层预测方面的重大突破。

图7是本申请又一个实施例中薄储层预测流程示意图，如图7所示，下面结合具体的示例，介绍本申请实施例的方案：

图7中主要包括两个部分，1目标分频数据的生成、2降频地层切片数据的生成。

其中，1目标分频数据的生成包括：

101、时深关系确定。以叠后地震数据和目标井的声波测井曲线为基础制作合成地震记录生成时深关系曲线。

102、-90度相位地震数据生成。通过叠后地震数据、目标井的时深关系曲线和储层特征曲线确定叠后地震数据的相位，通过希尔伯特变换生成-90度相位地震数据。相位的确定可根据地震数据的波形对称性来定性判断。

103、目标薄储层位置确定。通过-90度相位地震数据、目标井的时深关系曲线和储层特征曲线确定薄储层在-90度相位地震数据体的位置；

104、目标干涉频率确定。以-90度相位地震数据为基础，计算过目标井的小波时频谱，通过小波时频谱、目标井的时深关系曲线和薄储层的位置确定目标干涉频率。本申请实施例的目标干涉频率采用小波时频谱，但不作为本发明的限定，s变换时频谱、傅里叶变换时频谱、子波分解时频谱等方法同样适用于本申请。具体确定目标干涉频率的方法可以参考上述实施例，此处不再赘述。

105、目标分频数据体生成。以-90度相位地震数据为基础，根据确定的目标干涉频率通过小波变换方法生成目标分频数据。本申请实施例选择小波变换生成最佳分频数据体，但不作为本申请的限定。广义s变换、傅里叶变换、子波分解等技术同样适用于本申请。

2降频地层切片数据的生成具体包括：

201、参考层的选择。以叠后地震数据为基础，选择地质时间相同的标志地震同相轴作为参考层。

202、地层时间模型的建立。设定地层模型的起始相对地质时间和采样间隔，根据选择的参考层计算参考层间的平均距离或最大距离，根据平均距离或最大距离及采样间隔计算参考层之间的内插地层切片个数，根据内插地层切片个数和参考层通过线性内插函数建立地层时间模型。

203、降频地层切片数据的生成。沿地层时间模型中的每个地层层位从最佳分频数据体提取振幅形成一个降频地层切片数据。

利用降频地层切片数据，可以预测薄储层的分布范围。

本申请实施例适用于地震勘探或开发领域，特别是关于薄互储层中目标薄储层的预测方面，适合各个油田推广应用，特别适合岩性为碎屑岩、断层较少结构较为简单的区块，具有广阔的应用前景。

目前切片技术的基础数据体主要为-90度相位地震数据，该技术适合薄储层无干涉或干涉较小的地质情况，在海相盆地应用较广，然而在我国陆相含油气盆地中，薄互储层中的单一储层含油非常普遍，切片技术的应用主要集中在沉积相的预测，应用范围受到了较大限制。本申请解决了传统切片技术在薄互储层预测中的问题，主要适合薄互储层中目标薄储层的预测，可有效提高薄储层的预测精度，适合我国特殊的地质情况，因此，具有良好的应用前景。

基于上述所述的薄储层的预测方法，本说明书一个或多个实施例还提供一种薄储层的预测装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

具体地，图8是本申请提供的薄储层的预测装置一个实施例的模块结构示意图，如图8所示，本申请中提供的薄储层的预测装置包括：相位地震数据获取模块81、目标干涉频率确定模块82、时频信号分析模块83、地层时间模型构建模块84、降频处理模块85、薄储层预测模块86。

相位地震数据获取模块81，可以用于根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据；

目标干涉频率确定模块82，可以用于根据所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率；

时频信号分析模块83，可以用于根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据；

地层时间模型构建模块84，可以用于根据所述叠后地震数据构建地层时间模型；

降频处理模块85，可以用于沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据；

薄储层预测模块86，可以用于根据所述降频地层切片数据，预测薄储层的平面分布。

本申请提供的薄储层的预测装置，通过-90度相位地震数据、目标井的时深关系曲线和储层特征曲线确定薄储层的位置；通过小波时频谱、目标井的时深关系曲线和薄储层的位置确定目标干涉频率；根据确定的目标干涉频率通过小波变换等方法生成目标分频数据；通过目标分频数据和参考层生成降频地层切片数据。降频地层切片数据清晰的刻画了薄储层的分布范围，解决了由于薄储层干涉引起的切片范围预测结果不准的问题，提高了薄储层预测的精度。提出了一种与传统的垂向分辨薄层完全不同的研究思路，提供了一种新的薄储层预测技术，是地震沉积学在薄层预测方面的重大突破。

需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种薄储层的预测装置，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

存储器，所述指令被所述处理器执行时实现包括以下步骤：

根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据；

根据所述目标井的时频谱、所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率；

根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据；

根据所述叠后地震数据构建地层时间模型；

沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据；

根据所述降频地层切片数据，预测薄储层的平面分布。

所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如ram、rom等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、u盘；利用光学方式存储信息的装置如，cd或dvd。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。

需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。

上述实施例所述的一种薄储层的预测装置，提高了薄储层的预测精度。

本说明书实施例提供的上述薄储层的预测方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在pc端实现、linux系统实现，或其他例如使用android、ios系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种薄储层的预测系统的一个实施例中，图9是本申请提供的一种薄储层的预测系统的实施例的模块结构示意图，如图9所示，本申请另一实施例提供的薄储层的预测系统可以包括处理器91以及用于存储处理器可执行指令的存储器92，

处理器91和存储器92通过总线93完成相互间的通信；

所述处理器91用于调用所述存储器92中的程序指令，以执行上述各薄储层的预测方法实施例所提供的方法，例如包括：储器，所述指令被所述处理器执行时实现包括以下步骤：根据叠后地震数据、储层特征曲线、目标井的地震数据的时深关系曲线，获得-90度相位地震数据；根据所述目标井的时频谱、所述储层特征曲线、所述时深关系曲线、所述-90度相位地震数据，确定目标干涉频率；根据所述目标干涉频率对所述-90度相位地震数据进行时频信号分析，获得目标分频地震数据；根据所述叠后地震数据构建地层时间模型；沿所述地层时间模型中的地层层位从所述目标分频地震数据中提取振幅，获得降频地层切片数据；根据所述降频地层切片数据，预测薄储层的平面分布。

需要说明的是说明书上述所述的系统根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。

在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(programmablelogicdevice,pld)(例如现场可编程门阵列(fieldprogrammablegatearray，fpga))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片pld上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logiccompiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(hardwaredescriptionlanguage，hdl)，而hdl也并非仅有一种，而是有许多种，如abel(advancedbooleanexpressionlanguage)、ahdl(alterahardwaredescriptionlanguage)、confluence、cupl(cornelluniversityprogramminglanguage)、hdcal、jhdl(javahardwaredescriptionlanguage)、lava、lola、myhdl、palasm、rhdl(rubyhardwaredescriptionlanguage)等，目前最普遍使用的是vhdl(very-high-speedintegratedcircuithardwaredescriptionlanguage)与verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(applicationspecificintegratedcircuit，asic)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：arc625d、atmelat91sam、microchippic18f26k20以及siliconelabsc8051f320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

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

虽然本说明书一个或多个实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。第一，第二等词语用来表示名称，而并不表示任何特定的顺序。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

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

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

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

在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。

内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。

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

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

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

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

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