火成岩储层边界的刻画方法、装置及存储介质与流程

本发明涉及油气藏勘探开发技术领域，特别涉及一种火成岩储层边界的刻画的方法、装置及存储介质。

背景技术：

在火成岩储层的勘探开发过程中，由于不同类型的火成岩储层的孔隙和裂缝的发育特征差异较大，造成不同类型的火成岩储层的含油性或含气性也存在差异，因此，为了便于对火成岩储层的勘探开发，可以预先对不同类型的火成岩储层的边界进行刻画，以确定不同类型的火成岩储层的边界，从而基于不同类型的火成岩储层的边界，对孔隙和裂缝的发育相对较明显的火成岩储层进行勘探开发。比如，因火山活动频繁导致不同时期形成的喷出岩和侵入岩的孔隙和裂缝的发育差异较大，其中，喷出岩是熔浆由火山通道喷溢出地表凝固形成的岩石，该类岩石的类型主要有玄武岩和安山岩等，该类岩石在构造、风化和淋滤等因素的影响下孔隙和裂缝的发育相对较明显，侵入岩是熔浆在地壳的一定深度凝结而形成的岩石，孔隙和裂缝的发育相对较少。因此，亟需一种不同类型的火成岩储层边界的刻画方法以对火成岩储层的边界进行刻画，从而明确孔隙和裂缝的发育相对较明显的火成岩储层的分布范围，进而有效指导火成岩储层的井位部署。

技术实现要素：

为了确定火成岩储层中不同类型火成岩储层的边界，从而有效指导火成岩储层的井位部署问题，本发明提供了一种火成岩储层边界的刻画方法、装置及存储介质。所述技术方案如下：

第一方面，提供了一种火成岩储层边界的刻画方法，所述方法包括：

确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，所述目标区块是指待研究开发的区块，所述测井响应特征是指火成岩储层的岩性与测井介质之间的关系；

基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定所述目标区块内火成岩储层的实际顶界位置和实际底界位置；

确定每类火成岩储层的至少一个产量等级；

基于所述实际顶界位置和所述实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围；

基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对所述目标区块内不同类型的火成岩储层的边界进行刻画。

可选地，所述确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，包括：

获取所述目标区块内多口油气井中每口油气井的岩心观察资料、薄片鉴定资料、录井资料和测井资料；

基于每口油气井的岩心观察资料和薄片鉴定资料确定每口油气井沿深度方向的火成岩储层的类型；

基于每口油气井的录井资料和测井资料确定每口油气井沿深度方向的测井响应特征；

基于每口油气井沿深度方向的火成岩储层的类型和测井响应特征，确定所述目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。

可选地，所述基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定所述目标区块内火成岩储层的实际顶界位置和实际底界位置，包括：

基于每口油气井的岩心观察资料和薄片鉴定资料，确定所述目标区块的至少一个地质标志层；

提取每口油气井对应的地震子波，并从存储的所述目标区块的地震道数据中获取每口油气井的井旁地震道数据；

基于提取的地震子波和每口油气井的井旁地震道数据，在每个地质标志层的上层和/或下层对所述目标区块内的火成岩储层进行标定；

基于标定的火成岩储层和所述目标区块的地震道数据确定所述目标区块内火成岩储层的参考顶界位置和参考底界位置；

基于每类火成岩储层的测井响应特征，通过拟声波测井约束反演技术，对所述参考顶界位置和所述参考底界位置进行修正，并将修正后的参考顶界位置和参考底界位置确定为所述实际顶界位置和所述实际底界位置。

可选地，所述基于所述实际顶界位置和所述实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，包括：

从存储的所述目标区块的地震道数据中获取每类火成岩储层的地震道数据；

在纵向上对每类火成岩储层的地震道数据进行频谱分析，以确定每类火成岩储层分别对应的频率属性范围；

确定所述目标区块内所述实际顶界位置和所述实际底界位置之间的时间域长度；

基于所述时间域长度，确定所述目标区块内平均瞬时频率的平面分布；

基于每类火成岩储层分别对应的频率属性范围和所述目标区块内平均瞬时频率的平面分布，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

可选地，所述基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对所述目标区块内不同类型的火成岩储层的边界进行刻画，包括：

基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，通过分频技术提取所述目标区块内的预设控制时窗内火成岩储层对应的至少一个分频切片；

基于所述至少一个分频切片，对所述目标区块内不同类型的火成岩储层的边界进行刻画。

第二方面，提供了一种火成岩储层边界的刻画装置，所述装置包括：

第一确定模块，用于确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，所述目标区块是指待研究开发的区块，所述测井响应特征是指火成岩储层的岩性与测井介质之间的关系；

第二确定模块，用于基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定所述目标区块内火成岩储层的实际顶界位置和实际底界位置；

第三确定模块，用于确定每类火成岩储层的至少一个产量等级；

第四确定模块，用于基于所述实际顶界位置和所述实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围；

刻画模块，用于基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对所述目标区块内不同类型的火成岩储层的边界进行刻画。

可选地，所述第一确定模块包括：

第一获取单元，用于获取所述目标区块内多口油气井中每口油气井的岩心观察资料、薄片鉴定资料、录井资料和测井资料；

第一确定单元，用于基于每口油气井的岩心观察资料和薄片鉴定资料确定每口油气井沿深度方向的火成岩储层的类型；

第二确定单元，用于基于每口油气井的录井资料和测井资料确定每口油气井沿深度方向的测井响应特征；

第三确定单元，用于基于每口油气井沿深度方向的火成岩储层的类型和测井响应特征，确定所述目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。

可选地，所述第二确定模块包括：

第四确定单元，用于基于每口油气井的岩心观察资料和薄片鉴定资料，确定所述目标区块的至少一个地质标志层；

第二获取单元，用于提取每口油气井对应的地震子波，并从存储的所述目标区块的地震道数据中获取每口油气井的井旁地震道数据；

标定单元，用于基于提取的地震子波和每口油气井的井旁地震道数据，在每个地质标志层的上层和/或下层对所述目标区块内的火成岩储层进行标定；

第五确定单元，用于基于标定的火成岩储层和所述目标区块的地震道数据确定所述目标区块内火成岩储层的参考顶界位置和参考底界位置；

第六确定单元，用于基于每类火成岩储层的测井响应特征，通过拟声波测井约束反演技术，对所述参考顶界位置和所述参考底界位置进行修正，并将修正后的参考顶界位置和参考底界位置确定为所述实际顶界位置和所述实际底界位置。

可选地，所述第四确定模块包括：

第三获取单元，用于从存储的所述目标区块的地震道数据中获取每类火成岩储层的地震道数据；

第七确定单元，用于在纵向上对每类火成岩储层的地震道数据进行频谱分析，以确定每类火成岩储层分别对应的频率属性范围；

第八确定单元，用于确定所述目标区块内所述实际顶界位置和所述实际底界位置之间的时间域长度；

第九确定单元，用于基于所述时间域长度，确定所述目标区块内平均瞬时频率的平面分布；

第十确定单元，用于基于每类火成岩储层分别对应的频率属性范围和所述目标区块内平均瞬时频率的平面分布，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

可选地，所述刻画模块主要用于：

基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，通过分频技术提取所述目标区块内的预设控制时窗内火成岩储层对应的至少一个分频切片；

基于所述至少一个分频切片，对所述目标区块内不同类型的火成岩储层的边界进行刻画。

第三方面，提供了一种计算机可读存储介质，所述存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现上述第一方面提供的任一所述的方法。

本发明提供的技术方案带来的有益效果是：本发明通过确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，由于每类火成岩储层的岩性不同，造成每类火成岩储层的测井响应特征也不相同，因此，可以基于确定的每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置。之后，由于每类火成岩储层所在的位置不同时，该类火成岩储层的产量也可能不同，因此，为了明确目标区块内可以重点勘探开发的火成岩储层类型，可以预先确定每类火成岩储层的产量等级。进而可以基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，并基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，实现对目标区块内不同类型的火成岩储层的边界进行刻画，从而明确每类火成岩储层的边界，有效指导火成岩储层的井位部署，实现火成岩储层更为高效的勘探开发。

附图说明

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

图1是本发明实施例提供的一种火成岩储层边界的刻画方法的流程示意图；

图2是本发明实施例提供的一种火成岩储层边界的刻画方法的流程示意图；

图3a是本发明实施例提供的一种火成岩储层边界的刻画装置的结构示意图；

图3b是本发明实施例提供的一种火成岩储层边界的刻画装置的结构示意图；

图3c是本发明实施例提供的一种火成岩储层边界的刻画装置的结构示意图；

图3d是本发明实施例提供的一种火成岩储层边界的刻画装置的结构示意图；

图4是本发明实施例提供的一种火成岩储层边界的刻画装置的结构示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。

图1是本发明实施例提供的一种火成岩储层边界的刻画方法的流程示意图。参见图1，该方法包括如下步骤。

步骤101：确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，目标区块是指待研究开发的区块，该测井响应特征是指火成岩储层的岩性与测井介质之间的关系；

步骤102：基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置；

步骤103：确定每类火成岩储层的至少一个产量等级；

步骤104：基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围；

步骤105：基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对目标区块内不同类型的火成岩储层的边界进行刻画。

本发明实施例中，通过确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，由于每类火成岩储层的岩性不同，造成每类火成岩储层的测井响应特征也不相同，因此，可以基于确定的每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置。之后，由于每类火成岩储层所在的位置不同时，该类火成岩储层的产量也可能不同，因此，为了明确目标区块内可以重点勘探开发的火成岩储层类型，可以预先确定每类火成岩储层的产量等级。进而可以基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，并基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，实现对目标区块内不同类型的火成岩储层的边界进行刻画，从而明确每类火成岩储层的边界，有效指导火成岩储层的井位部署，实现火成岩储层更为高效的勘探开发。

可选地，确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，包括：

获取目标区块内多口油气井中每口油气井的岩心观察资料、薄片鉴定资料、录井资料和测井资料；

基于每口油气井的岩心观察资料和薄片鉴定资料确定每口油气井沿深度方向的火成岩储层的类型；

基于每口油气井的录井资料和测井资料确定每口油气井沿深度方向的测井响应特征；

基于每口油气井沿深度方向的火成岩储层的类型和测井响应特征，确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。

可选地，基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置，包括：

基于每口油气井的岩心观察资料和薄片鉴定资料，确定目标区块的至少一个地质标志层；

提取每口油气井对应的地震子波，并从存储的目标区块的地震道数据中获取每口油气井的井旁地震道数据；

基于提取的地震子波和每口油气井的井旁地震道数据，在每个地质标志层的上层和/或下层对目标区块内的火成岩储层进行标定；

基于标定的火成岩储层和目标区块的地震道数据确定目标区块内火成岩储层的参考顶界位置和参考底界位置；

基于每类火成岩储层的测井响应特征，通过拟声波测井约束反演技术，对参考顶界位置和参考底界位置进行修正，并将修正后的参考顶界位置和参考底界位置确定为实际顶界位置和实际底界位置。

可选地，基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，包括：

从存储的目标区块的地震道数据中获取每类火成岩储层的地震道数据；

在纵向上对每类火成岩储层的地震道数据进行频谱分析，以确定每类火成岩储层分别对应的频率属性范围；

确定目标区块内实际顶界位置和实际底界位置之间的时间域长度；

基于该时间域长度，确定目标区块内平均瞬时频率的平面分布；

基于每类火成岩储层分别对应的频率属性范围和目标区块内平均瞬时频率的平面分布，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

可选地，基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对目标区块内不同类型的火成岩储层的边界进行刻画，包括：

基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，通过分频技术提取目标区块内的预设控制时窗内火成岩储层对应的至少一个分频切片；

基于该至少一个分频切片，对目标区块内不同类型的火成岩储层的边界进行刻画。

上述所有可选技术方案，均可按照任意结合形成本发明的可选实施例，本发明实施例对此不再一一赘述。

图2是本发明实施例提供的一种火成岩储层边界的刻画方法的流程示意图。参见图2，该方法包括如下步骤。

步骤201：确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。

其中，目标区块是指待研究开发的区块，测井响应特征是指火成岩储层的岩性与测井介质之间的关系，测井响应特征可以包括自然伽马响应特征、补偿密度响应特征、电阻率响应特征、声波响应特征和自然电位响应特征等。

在对目标区块内的火成岩储层进行勘探开发时，由于不同类型的火成岩储层的含油性或含气性之间存在差异，因此，可以预先确定目标区块内火成岩储层的类型，以选择含油量或含气量较大的火成岩储层，实现对目标区块内火成岩储层的高效开发。在确定目标区块内火成岩储层的类型后，还可以确定每类火成岩储层的测井响应特征。具体地，可以通过如下步骤(1)-(4)实现。

(1)、获取目标区块内多口油气井中每口油气井的岩心观察资料、薄片鉴定资料、录井资料和测井资料。

为了方便对位于地下的火成岩储层的研究，可以通过获取目标区块内多口油气井中每口油气井的油气井资料来了解目标区块内的地下情况，从而实现对目标区块内火成岩储层的研究，进而确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。其中，油气井资料包括岩心观察资料、薄片鉴定资料、录井资料和测井资料等。

在获取目标区块内每口油气井的油气井资料时，由于每口油气井的钻井时间可能不同，造成每口油气井的油气井资料可能不全。比如，对于目标区块内的不同时期开发的油气井，可能由于重点研究的地层不同，导致部分油气井的地层未进行岩心观察和薄片鉴定，从而造成该部分油气井缺少岩心观察资料和薄片鉴定资料，也可能导致部分油气井未进行录井解释，从而造成该部分油气井缺少录井资料。因此，在获取目标区块内每口油气井的油气井资料时，还可以对目标区块内未进行岩心观察和薄片鉴定的油气井补充岩心观察资料和薄片鉴定资料，对目标区块内未进行录井解释的油气井补充录井解释资料和测井解释资料，以保证目标区块内该多口油气井的油气井资料更为完整，提高目标区块内火成岩储层的类型的准确度，以及每类火成岩储层的测井响应特征的准确度。

(2)、基于每口油气井的岩心观察资料和薄片鉴定资料确定每口油气井沿深度方向的火成岩储层的类型。

由于岩心观察资料和薄片鉴定资料中都可以包括油气井沿深度方向的地层的岩性，而地层的岩性可以反映地层的岩石类型，也即是，可以通过地层的岩性确定地层是否为火成岩储层，以及火成岩储层中的哪一类，因此，可以从每口油气井的岩心观察资料和薄片鉴定资料中确定目标区块内每口油气井沿深度方向的火成岩储层的类型。

其中，每口油气井沿深度方向可以包括一类火成岩储层，也可以包括多类火成岩储层，本发明实施例对此不做限定。

(3)、基于每口油气井的录井资料和测井资料确定每口油气井沿深度方向的测井响应特征。

由于录井资料和测井资料中都可以包括油气井沿深度方向的测井曲线，测井曲线可以包括自然伽马测井曲线、补偿密度测井曲线、电阻率测井曲线、声波测井曲线和自然电位测井曲线等。因此，对于每口油气井的录井资料和测井资料包括的测井曲线，当该测井曲线为自然伽马测井曲线时，由于自然伽马测井曲线用于反映每口油气井在不同深度时自然伽马射线的强度，因此，可以将自然伽马测井曲线确定为用于反映自然伽马射线强度的自然伽马响应特征。当该测井曲线为补偿密度测井曲线时，由于补偿密度测井曲线用于反映每口油气井在不同深度时地层的岩石体积密度，因此，可以将补偿密度测井曲线确定为用于反映岩石体积密度的补偿声波响应特征。当该测井曲线为电阻率测井曲线时，由于电阻率测井曲线用于反映每口油气井在不同深度时岩石介质中的电场分布情况，因此，可以将电阻率测井曲线确定为用于反映岩石介质中电场分布的电阻率响应特征。当该测井曲线为声波测井曲线时，由于声波测井曲线用于反映每口油气井在不同深度时的声波传播速度，因此，可以将声波测井曲线确定为用于反映声波传播速度的声波响应特征。当该测井曲线为自然电位测井曲线时，由于自然电位测井曲线用于反映每口油气井在不同深度时在电化学作用下产生的电位，因此，可以将自然电位测井曲线确定为用于反映在电化学作用下产生的电位的自然电位响应特征。

(4)、基于每口油气井沿深度方向的火成岩储层的类型和测井响应特征，确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。

在通过上述步骤(2)确定了每口油气井沿深度方向的火成岩储层的类型后，由于每口油气井是在目标区块内进行布置的，相应地，每口油气井沿深度方向的火成岩储层也位于目标区块的地层内，因此，可以将每口油气井沿深度方向的火成岩储层的类型确定为目标区块内火成岩储层的类型。比如，第一口油气井的火成岩储层的类型为玄武岩储层，第二口油气井的火成岩储层的类型为浅灰绿色辉绿岩储层和安山岩储层，则可以将玄武岩储层、浅灰绿色辉绿岩储层和安山岩储层确定为目标区块内火成岩储层的类型。

由于每口油气井在同一深度的火成岩储层为同一类型，也即是，每口油气井在同一深度只有一种类型的火成岩储层，并且由于不同油气井中同一类火成岩储层的岩性可能不相同，且同一油气井中同一类火成岩储层在不同深度的岩性也可能不同，因此，对于目标区块内任一类型的火成岩储层来说，为了便于说明，将该任一类型的火成岩储层称为目标类型的火成岩储层，这样，可以将目标区块内所有油气井中目标类型的火成岩储层的深度范围内的测井响应特征确定为目标区块内目标类型的火成岩储层的测井响应特征。

比如，第一口油气井的a类火成岩储层位于第一深度至第二深度的深度范围内，且随着深度的不同a类火成岩储层的岩性也不相同，第二口油气井的a类火成岩储层位于第三深度至第四深度的深度范围内，且随着深度的不同a类火成岩储层的岩性不相同，另外，位于第一深度至第二深度的深度范围内的a类火成岩储层的岩性与位于第三深度至第四深度的a类火成岩储层的岩性也不同。因此，将目标区块内的a类火成岩储层称为目标类型的火成岩储层，这样，可以将第一口油气井中第一深度至第二深度的深度范围内的测井响应特征和第二口油气井中第三深度至第四深度的深度范围内的测井响应特征确定为目标区块内目标类型的火成岩储层的测井响应特征。

在确定了目标区块内每类火成岩储层的测井响应特征后，可以基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置。具体地，可以通过如下步骤202-步骤206实现。其中，地震-测井联合的方法是指目标区块的地震资料与目标区块内每口油气井的测井资料相互结合、相互印证的一种方法。其中，地震资料包括用于反映目标区块内地质构造的地震道数据等。

步骤202：基于每口油气井的岩心观察资料和薄片鉴定资料，确定目标区块的至少一个地质标志层。

由于每口油气井的岩心观察资料和薄片鉴定资料均包括地层岩石的岩性，而目标区块的地质标志层均具有所含化石和岩性特征明显、层位稳定、分布范围广、易于鉴别的特点，因此，可以通过对每口油气井的岩心观察资料和薄片鉴定资料包括的沿深度方向的地层的岩性进行分析，将地层的岩性出现明显变化的位置所处的地层确定为每口油气井对应的地质标志层，从而得到目标区块的至少一个地质标志层。当然，还可以通过其他方法确定每口油气井对应的地质标志层，比如，可以将每口油气井的岩心观察资料和薄片鉴定资料包括的沿深度方向的地层的岩性和预先存储的地质标志层的岩性进行对比，当岩心观察资料和薄片鉴定资料包括的沿深度方向的地层的岩性和预先存储的地质标志层的岩性的匹配度大于预设数值时，将该岩性所处的地层确定为地质标志层。其中，预设数值可以为0.7、0.8或0.9等。

比如，目标区块的地质标志层可以为馆陶组地层、沙一段地层和孔二段地层等。

步骤203：提取每口油气井对应的地震子波，并从存储的目标区块的地震道数据中获取每口油气井的井旁地震道数据。

其中，在提取每口油气井对应的地震子波时，可以通过确定性子波提取方法和统计性子波提取方法进行提取，具体的提取过程可以参见现有技术，本发明实施例在此不再赘述。然而，由于提取的地震子波受地质变化的影响较大，且地震子波是一段具有确定的起始时间和一定延续长度的信号，因此，可以通过上述方法多次多反复提取每口油气井对应的地震子波。

步骤204：基于提取的地震子波和每口油气井的井旁地震道数据，在每个地质标志层的上层和/或下层对目标区块内的火成岩储层进行标定。

具体地，基于提取的地震子波确定每口油气井的合成地震记录，进而基于每口油气井的合成地震记录和每口油气井的井旁地震道数据，在每个地质标志层的上层和/或下层对目标区块内的火成岩储层进行标定。

其中，对于每口油气井的合成地震记录，由于提取的地震子波受地质变化的影响较大，且地震子波是一段具有确定的起始时间和一定延续长度的信号，因此，可以基于上述步骤203多次提取的地震子波确定每口油气井的最优合成地震记录。

在确定每口油气井的最优合成地震记录时，可以通过每口油气井的声波测井资料和密度测井资料确定每口油气井的反射系数序列，再将确定的反射系数序列与提取的地震子波通过反褶积算法得到每口油气井的多个合成地震记录，进而从多个地震合成记录中，查找与每口油气井的井旁地震道数据反映的地震剖面层位吻合的现象最多的一个合成地震记录，并将层位吻合现象最多的合成地震记录确定为最优合成地震记录。其中，通过声波测井资料和密封测井资料确定反射系数序列的具体计算方法可以参考现有技术，本发明实施例在此不做详细赘述。当然，本发明实施例还可以通过其他方法确定每口油气井的合成地震记录，本发明实施例对此不做限定。

其中，在每个地质标志层的上层和/或下层对目标区块内的火成岩储层进行标定时，可以先通过预先存储的目标区块内的时深关系表将每口油气井在深度域的每类火成岩储层转换成时间域的每类火成岩储层，进而基于每类火成岩储层所处时间域的时长，通过每口油气井的合成地震记录与井旁地震道数据，对每口油气井的每类火成岩储层进行标定，从而实现对目标区块内每类火成岩储层的标定。其中，时深关系表是指地震波传输至一定深度的地层并反射至地表所需要的时间与该地层的深度之间的对应关系表。基于每类火成岩储层所处时间域的时长，通过每口油气井的合成地震记录与井旁的地震道数据，对每类火成岩储层的标定的具体方法可以参考现有技术，本发明实施例在此不做详细赘述。

步骤205：基于标定的火成岩储层和目标区块的地震道数据确定目标区块内火成岩储层的参考顶界位置和参考底界位置。

在对目标区块内的火成岩储层进行标定后，可以明确目标区块内火成岩储层所处的位置，进而基于目标区块内火成岩储层所处的位置，以及目标区块的地震道数据对标定的火成岩储层的顶界和底界的延伸方向进行追踪，以确定目标区块内火成岩储层的参考顶界位置和参考底界位置。也即是，将目标区块内所有类型的火成岩储层作为一个整体，该参考顶界位置是指这个整体的最顶端的位置，该参考底界位置是指这个整体的最底端的位置。

其中，基于目标区块内火成岩储层所处的位置，以及目标区块的地震道数据对标定的火成岩储层的顶界和底界的延伸方向进行追踪的方法可以参考现有技术，本发明实施例在此不做详细赘述。

步骤206：基于每类火成岩储层的测井响应特征，通过拟声波测井约束反演技术，对参考顶界位置和参考底界位置进行修正，并将修正后的参考顶界位置和参考底界位置确定为实际顶界位置和实际底界位置。

目标区块的地震道数据主要是通过向地层内发射地震波，并在地面接收地层反射的地震波的方式获取的，由于地震波在沿目标区块的深度方向的分辨率较低，可能造成目标区块内火成岩储层的参考顶界位置和参考底界位置不太精确，因此，为了提高目标区块内每类火成岩储层边界的精确度，可以对参考顶界位置和参考地界位置进行修正。

具体地，基于每类火成岩储层的测井响应特征，确定每类火成岩储层的拟声波测井曲线，进而基于拟声波测井曲线通过拟声波测井约束反演技术对参考顶界位置和参考底界位置进行修正，并将修正后的参考顶界位置和参考底界位置确定为目标区块内火成岩储层的实际顶界位置和实际底界位置。

其中，在确定每类火成岩储层的拟声波测井曲线时，可以通过整流和滤波处理方法提取自然伽玛测井曲线的高频分量，通过滤波处理方法提取声波测井曲线的低频分量，进而将提取的高频分量和低频分量进行组合，得到既能反映地层速度和波阻抗的变化，又能反映岩性的细微差异的拟声波测井曲线；也可以通过整流和滤波处理方法提取电阻率测井曲线的高频分量，通过滤波处理方法提取声波测井曲线的低频分量，进而将提取的高频分量和低频分量进行组合，得到既能反映电阻率的高频特征，又能反映储层特征的拟声波测井曲线，本发明实施例对此不做限定。

其中，通过拟声波测井约束反演技术对参考顶界位置和参考底界位置的修正可以参考现有技术，本发明实施例在此不再赘述。

步骤207：确定每类火成岩储层的至少一个产量等级。

对于目标区块内的每类火成岩储层，由于火成岩储层的分布位置不同，可能造成火成岩储层的产量也不相同，因此，为了明确目标区块内产量较高的每类火成岩储层的分布位置，可以预先基于目标区块内每口油气井的试油试采产量，以及预先设置的产量等级划分标准，对每类火成岩储层进行划分，从而得到每类火成岩储层的至少一个产量等级。同时，对于任一口油气井，在通过该口油气井的试油试采产量对该口油气井中的一类火成岩储层进行产量等级划分时，可以将该口油气井所在的位置区域确定为该类火成岩储层的产量等级所分布的区域，从而可以得到目标区块内每类火成岩储层的至少一个产量等级的平面分布。

其中，对每类火成岩储层进行划分的具体操作可以为：基于预先设置的产量等级划分标准对每口油气井的试油试采产量进行划分，以确定每口油气井的产量等级，进而将每口油气井的产量等级确定为每口油气井中的火成岩储层的产量等级，从而实现对目标区块内每类火成岩储层的划分。比如，第一口油气井和第二口油气井的火成岩储层均为a类火成岩储层，第一口油气井的试油试采产量为18立方米每千米每天，第二口油气井的试油试采产量为3立方米每千米每天，如果18立方米每千米每天属于产量等级中的高产等级，则将高产等级确定为第一口油气井中a类火成岩储层的产量等级，如果3立方米每千米每天属于产量等级中的低产等级，则将低产等级确定为第二口油气井中a类火成岩储层的产量等级，从而实现对目标区块内a类火成岩储层的产量等级的划分。

其中，预先设置的产量等级划分标准可以预先根据油气井每千米井深的稳定产量进行设置。比如，产量等级划分标准可以设置为第一等级、第二等级、第三等级和第四等级，第一等级为高产等级，千米井深的稳定产量大于或等于15立方米每千米每天；第二等级为中产等级，千米井深的稳定产量大于或等于5立方米每千米每天且小于15立方米每千米每天；第三等级为低产等级，千米井深的稳定产量大于或等于1立方米每千米每天且小于5立方米每千米每天；第四等级为特低产等级，千米井深的稳定产量小于1立方米每千米每天。

步骤208：基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

具体地，从存储的目标区块的地震道数据中获取每类火成岩储层的地震道数据，在纵向上对每类火成岩储层的地震道数据进行频谱分析，以确定每类火成岩储层分别对应的频率属性范围。之后，确定目标区块内实际顶界位置和实际底界位置之间的时间域长度，基于该时间域长度，确定目标区块内平均瞬时频率的平面分布。进而基于每类火成岩储层分别对应的频率属性范围和目标区块内平均瞬时频率的平面分布，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

其中，可以通过频谱分析仪对每类火成岩储层的地震道数据进行频谱分析，以确定每类火成岩储层分别对应的频率属性范围，具体分析方法可以参考相关技术，本发明实施例在此不做赘述。

其中，对于目标区块内实际顶界位置和实际底界位置之间的时间域长度，可以将实际顶界位置对应的深度值和实际底界位置对应的深度值，通过时深关系表分别转换为实际顶界位置对应的时间域数值和实际底界位置对应的时间域数值，进而将实际顶界位置对应的时间域数值和实际底界位置对应的时间域数值之间的差值确定为实际顶界位置和实际底界位置之间的时间域长度。在确定了实际顶界位置和实际底界位置之间的时间域长度，以及每类火成岩储层分别对应的频率属性范围后，在火成岩储层的实际顶界位置和实际底界位置之间的时间域长度内，获取目标区块内每一地震道数据沿时间轴的瞬时频率，进而确定该时间域长度内每一地震道数据的所有瞬时频率的平均瞬时频率，从而得到目标区块内每一地震道数据的平均瞬时频率的平面分布，也即是目标区块内平均瞬时频率的平面分布。

在确定了目标区块内平均瞬时频率的平面分布后，可以通过每类火成岩储层分别对应的频率属性范围，在目标区块内的平均瞬时频率的平面分布上查找位于每类火成岩储层分别对应的频率属性范围的至少一个区域，进而将该至少一个区域的分布确定为目标区块内每类火成岩储层的平面分布。之后，通过每类火成岩储层对应的至少一个产量等级的平面分布，从每类火成岩储层的平面分布中，确定每类火成岩储层的至少一个产量等级中每个产量等级的平面分布范围，进而将每个产量等级所在的平面分布范围内的平均瞬时频率所在的范围，确定为每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

继续上述举例，玄武岩储层的中产等级对应的频率属性范围为18赫兹到22赫兹，浅灰绿色辉绿岩储层的特低产等级对应的频率属性范围为25赫兹到30赫兹。

进一步地，对于目标区块内无法采用平均瞬时频率进行区分的每类火成岩储层，还可以基于实际顶界位置和实际底界位置之间的时间域长度，提取每类火成岩储层的振幅属性范围或相位属性范围等，同时，确定目标区块内平面上的平均振幅属性的分布或平均相位属性的分布等，从而确定每类火成岩储层的至少一个产量等级中每个产量等级对应的振幅属性范围或相位属性范围等，具体地实现过程可以参考上述确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围的确定方法，本发明实施例在此不做详细赘述。

步骤209：基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对目标区块内不同类型的火成岩储层的边界进行刻画。

具体地，基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，通过分频技术提取目标区块内的预设控制时窗内火成岩储层对应的至少一个分频切片，基于该至少一个分频切片，对目标区块内不同类型的火成岩储层的边界进行刻画。

其中，通过分频技术提取分频切片可以参考相关技术，本发明实施例在此不做详细赘述。预设控制时窗可以基于目标区块内每类火成岩储层所处的位置进行预先设置，比如，目标区块内火成岩储层的顶部储层裂缝相对更发育，因此，可以将火成岩储层的实际顶界位置以及向下20毫秒确定为控制时窗，本发明实施例对此不做限定。

其中，分频切片用于指示目标区块每类火成岩储层的至少一个产量等级在平面上的分布，当分频切片的分频频率属于某一个产量等级对应的频率属性范围时，该产量等级的边界在该分频切片中显示的越明显。因此，以每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围包括的至少一个频率作为基础提取的至少一个分频切片，可以很明显地显示目标区块内每类火成岩储层的不同产量等级的分布边界，从而实现对目标区块内不同类型的火成岩储层的边界的刻画。

继续上述举例，以玄武岩储层的至少一个产量等级中每个产量等级对应的频率属性范围18～22赫兹包括的18赫兹、20赫兹和22赫兹提取分频切片时，可以从18赫兹、20赫兹和22赫兹分别对应的三个分频切片中选择一个玄武岩储层的中产等级的边界最为明显的分频切片，进而沿边界最明显的分频切片的边界对玄武岩储层的中产等级的边界进行刻画。

需要说明的是，在通过分频技术提取目标区块内的预设控制时窗内火成岩储层对应的至少一个分频切片时，还可以直接以每类火成岩储层的每个产量等级对应的频率属性范围作为基础，提取预设控制时窗内至少一个不同频率范围的分频切片，本发明实施例对此不做限定。

进一步地，对于目标区块内每类火成岩储层的不同的产量等级，如果无法通过提取的至少一个分频切片对不同的产量等级的分布边界进行区分时，还可以基于每个产量等级对应的振幅属性范围或相位属性范围等，通过分频切片上的振幅属性或相位属性等对不同的产量等级的边界进行区分，从而实现对目标区块内不同类型的火成岩储层的边界的刻画。

本发明实施例中，由于油气井资料中包括油气井沿深度方向的岩性特征和测井响应特征，因此，可以通过获取目标区块内每口油气井的油气井资料，确定目标区块内火成岩储层的类型和每类火成岩储层的测井响应特征。之后，通过每口油气井包括的能够反映延伸度方向的岩性特征的油气井资料确定目标区块的至少一个地质标志层，进而在至少一个地质标志层的上层和/或下层，通过每口油气井的合成地震记录和井旁的地震道数据标定火成岩储层的参考顶界位置和参考底界位置，以提高标定效率。为了避免参考顶界位置和参考底界位置存在偏差，可以通过拟声波测井约束反演技术对参考顶界位置和参考底界位置进行修正。之后，通过纵向上对地震道数据的频谱分析，以及平面上平均瞬时频率的提取，明确目标区块内平面上每类火成岩储层的分布，由于每类火成岩储层的位置不同时含油或含气量也可能不相同，因此，可以预先对每类火成岩储层进行产量分类，进而确定每类火成岩储层的不同产量等级对应的频率属性范围。进而基于每类火成岩储层的不同产量等级对应的频率属性范围，通过分频技术得到可以反映每类火成岩储层的不同产量等级的边界，从而实现对不同类型的火成岩储层的边界的刻画，以实现后期对每类火成岩储层的高产量等级的勘探开发，提高火成岩储层的开发效率。

图3a是本发明实施例提供的一种火成岩储层边界的刻画装置的结构示意图。参见图3a，该装置包括：

第一确定模块301，用于确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，目标区块是指待研究开发的区块，该测井响应特征是指火成岩储层的岩性与测井介质之间的关系；

第二确定模块302，用于基于每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置；

第三确定模块303，用于确定每类火成岩储层的至少一个产量等级；

第四确定模块304，用于基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围；

刻画模块305，用于基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，对目标区块内不同类型的火成岩储层的边界进行刻画。

可选地，参见图3b，第一确定模块301包括：

第一获取单元3011，用于获取目标区块内多口油气井中每口油气井的岩心观察资料、薄片鉴定资料、录井资料和测井资料；

第一确定单元3012，用于基于每口油气井的岩心观察资料和薄片鉴定资料确定每口油气井沿深度方向的火成岩储层的类型；

第二确定单元3013，用于基于每口油气井的录井资料和测井资料确定每口油气井沿深度方向的测井响应特征；

第三确定单元3014，用于基于每口油气井沿深度方向的火成岩储层的类型和测井响应特征，确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征。

可选地，参见图3c，第二确定模块302包括：

第四确定单元3021，用于基于每口油气井的岩心观察资料和薄片鉴定资料，确定目标区块的至少一个地质标志层；

第二获取单元3022，用于提取每口油气井对应的地震子波，并从存储的目标区块的地震道数据中获取每口油气井的井旁地震道数据；

标定单元3023，用于基于提取的地震子波和每口油气井的井旁地震道数据，在每个地质标志层的上层和/或下层对目标区块内的火成岩储层进行标定；

第五确定单元3024，用于基于标定的火成岩储层和目标区块的地震道数据确定目标区块内火成岩储层的参考顶界位置和参考底界位置；

第六确定单元3025，用于基于每类火成岩储层的测井响应特征，通过拟声波测井约束反演技术，对参考顶界位置和参考底界位置进行修正，并将修正后的参考顶界位置和参考底界位置确定为实际顶界位置和实际底界位置。

可选地，参见图3d，第四确定模块304包括：

第三获取单元3041，用于从存储的目标区块的地震道数据中获取每类火成岩储层的地震道数据；

第七确定单元3042，用于在纵向上对每类火成岩储层的地震道数据进行频谱分析，以确定每类火成岩储层分别对应的频率属性范围；

第八确定单元3043，用于确定目标区块内实际顶界位置和实际底界位置之间的时间域长度；

第九确定单元3044，用于基于该时间域长度，确定目标区块内平均瞬时频率的平面分布；

第十确定单元3045，用于基于每类火成岩储层分别对应的频率属性范围和目标区块内平均瞬时频率的平面分布，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围。

可选地，刻画模块305主要用于：

基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，通过分频技术提取目标区块内的预设控制时窗内火成岩储层对应的至少一个分频切片；

基于该至少一个分频切片，对目标区块内不同类型的火成岩储层的边界进行刻画。

本发明实施例中，通过确定目标区块内火成岩储层的类型，以及每类火成岩储层的测井响应特征，由于每类火成岩储层的岩性不同，造成每类火成岩储层的测井响应特征也不相同，因此，可以基于确定的每类火成岩储层的测井响应特征，通过地震-测井联合的方法，确定目标区块内火成岩储层的实际顶界位置和实际底界位置。之后，由于每类火成岩储层所在的位置不同时，该类火成岩储层的产量也可能不同，因此，为了明确目标区块内可以重点勘探开发的火成岩储层类型，可以预先确定每类火成岩储层的产量等级。进而可以基于实际顶界位置和实际底界位置，确定每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，并基于每类火成岩储层的至少一个产量等级中每个产量等级对应的频率属性范围，实现对目标区块内不同类型的火成岩储层的边界进行刻画，从而明确每类火成岩储层的边界，有效指导火成岩储层的井位部署，实现火成岩储层更为高效的勘探开发。

需要说明的是：上述实施例提供的火成岩储层边界的刻画装置在对火成岩储层的边界进行刻画时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的火成岩储层边界的刻画装置与火成岩储层边界的刻画方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。

图4示出了本发明一个示例性实施例提供的终端400的结构框图。参见图4，该终端400可以是：智能手机、平板电脑、笔记本电脑或台式电脑。终端400还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。参见图4，终端400可以包括处理器401和存储器402。

处理器401可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器401可以采用dsp(digitalsignalprocessing，数字信号处理)、fpga(field－programmablegatearray，现场可编程门阵列)、pla(programmablelogicarray，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器401也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称cpu(centralprocessingunit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器401可以在集成有gpu(graphicsprocessingunit，图像处理器)，gpu用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器401还可以包括ai(artificialintelligence，人工智能)处理器，该ai处理器用于处理有关机器学习的计算操作。

存储器402可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器402还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器402中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器401所执行以实现本申请中方法实施例提供的一种火成岩储层边界的刻画方法。

在一些实施例中，终端400还可选包括有：外围设备接口403和至少一个外围设备。处理器401、存储器402和外围设备接口403之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口403相连。具体地，外围设备包括：射频电路404、显示屏405、定位组件406和电源407中的至少一种。

外围设备接口403可被用于将i/o(input/output，输入/输出)相关的至少一个外围设备连接到处理器401和存储器402。在一些实施例中，处理器401、存储器402和外围设备接口403被集成在同一芯片或电路板上；在一些其他实施例中，处理器401、存储器402和外围设备接口403中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。

射频电路404用于接收和发射rf(radiofrequency，射频)信号，也称电磁信号。射频电路404通过电磁信号与通信网络以及其他通信设备进行通信。射频电路404将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路404包括：天线系统、rf收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路404可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：万维网、城域网、内联网、各代移动通信网络(2g、3g、4g及5g)、无线局域网和/或wifi(wirelessfidelity，无线保真)网络。在一些实施例中，射频电路404还可以包括nfc(nearfieldcommunication，近距离无线通信)有关的电路，本申请对此不加以限定。

显示屏405用于显示ui(userinterface，用户界面)。该ui可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏405是显示屏时，显示屏405还具有采集在显示屏405的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器401进行处理。此时，显示屏405还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏405可以为一个，设置终端400的前面板；在另一些实施例中，显示屏405可以为至少两个，分别设置在终端400的不同表面或呈折叠设计；在再一些实施例中，显示屏405可以是柔性显示屏，设置在终端400的弯曲表面上或折叠面上。甚至，显示屏405还可以设置成非矩形的不规则图形，也即异形屏。显示屏405可以采用lcd(liquidcrystaldisplay，液晶显示屏)、oled(organiclight-emittingdiode,有机发光二极管)等材质制备。

定位组件406用于定位终端400的当前地理位置，以实现导航或lbs(locationbasedservice，基于位置的服务)。定位组件406可以是基于美国的gps(globalpositioningsystem，全球定位系统)、中国的北斗系统或俄罗斯的伽利略系统的定位组件。

电源407用于为终端400中的各个组件进行供电。电源407可以是交流电、直流电、一次性电池或可充电电池。当电源407包括可充电电池时，该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池，无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。

本领域技术人员可以理解，图4中示出的结构并不构成对终端400的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。

在上述实施例中，还提供了一种包括指令的非暂态的计算机可读存储介质，该计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器所执行以实现上述图1或图2所示实施例提供的方法。

本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。