CN109870719B

CN109870719B - 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统

Info

Publication number: CN109870719B
Application number: CN201910072027.4A
Authority: CN
Inventors: 何争光; 张亚娟; 李金付
Original assignee: China National Petroleum Corp; BGP Inc
Current assignee: China National Petroleum Corp; BGP Inc
Priority date: 2019-01-25
Filing date: 2019-01-25
Publication date: 2021-01-29
Anticipated expiration: 2039-01-25
Also published as: CN109870719A

Abstract

本说明书实施例公开了一种碳酸盐岩致密薄储层的井位布设方法、装置及系统，所述方法包括三方面主要内容：今构造高有利于成藏；古构造高控制沉积相，沉积相控制储层分布；低泊松比反演与流体活动性检测预测含气性，通过三方面内容综合确定碳酸盐岩有效储层分布。今构造高利用二维十字线控制或三维地震控制下，利用大比例尺小网格变速成图方法确定目标工区的今构造高点；对所述今构造高点对应目标层段进行层位拉平恢复古构造，并在恢复的古构造的基础上进行地震相分析处理，获得沉积相带分布结果；根据所述沉积相带分布结果确定所述目标工区的储层分布结果。在此基础上，通过泊松比反演、流体活动性检测技术预测碳酸盐岩气层分布有利区。

Description

一种碳酸盐岩致密薄储层的井位布设方法、装置及系统

技术领域

本发明涉及石油天然气勘探技术领域，特别地适用于大型潜台式碳酸盐岩薄储层地震预测。

背景技术

世界油气产量的60％来自于碳酸盐岩储层，但由于碳酸盐岩致密薄储层的非均质性，地震横向预测一直是地球物理勘探的一个难点，国内外地球物理界虽经多年攻关，但迄今未有较理想的解决方法。

碳酸盐岩储层的储集空间主要为溶孔、溶洞和微细裂缝，储层单层有效厚度小(小于5m)，非均质性强烈等特点。通过对大量已知井地质规律进行分析归纳，主力储层的厚度、空间展布形态等情况均受侵蚀地貌控制。即在侵蚀量大的潜沟中，因主力储层缺失而不含气，相反，地层保存较全的潜台区含气好。但碳酸盐岩储层波阻抗差异较小、地震反射弱，岩性横向变化大，白云岩储存及有效储层均为薄储层，含气性预测困难，目标层段成藏主控因素也认识不清。因此，目前亟需一种可以有效预测识别碳酸盐岩储层有利区的方法。

发明内容

本说明书实施例的目的在于提供一种碳酸盐岩致密薄储层的井位布设方法、装置及系统，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

本说明书提供一种碳酸盐岩致密薄储层分布确定方法、装置及系统是包括如下方式实现的：

一种碳酸盐岩致密薄储层分布确定方法，包括：

基于预先构建的速度场，利用大比例尺小网格变速成图方法确定目标工区的今构造高点；

对所述今构造高点对应目标层段进行层位拉平恢复古构造，并在恢复的古构造的基础上进行地震相分析处理，获得沉积相带分布结果；

根据所述沉积相带分布结果确定所述目标工区的储层分布结果。

本说明书提供的所述方法的另一个实施例中，所述方法还包括：

根据模型层析法构建速度场；

相应的，所述确定目标工区的今构造高点包括基于所述速度场利用大比例尺小网格变速成图方法确定目标工区的今构造高点。

本说明书提供的所述方法的另一个实施例中，所述确定目标工区的今构造高点，包括：

基于合成地震记录进行层位标定，确定标志层反射及目标层段小层反射特征；

根据所述标志层反射及目标层段小层反射特征进行地震层位小层段解释，获得层位解释结果；

根据所述层位解释结果，利用大比例尺小网格变速成图方法确定目标工区的今构造高点。

本说明书提供的所述方法的另一个实施例中，所述在恢复的古构造的基础上进行地震相分析处理，获得沉积相带分布结果，包括：

在恢复的古构造的基础上进行地震相分析处理，确定古构造高点以及地震相异常特征；

根据所述古构造高点以及地震相异常特征确定有利沉积相带。

本说明书提供的所述方法的另一个实施例中，所述确定所述目标工区的储层分布结果，包括：

获取目标工区的泊松比数据，根据所述泊松比数据对所述目标工区进行含气性预测；

根据所述沉积相带分布结果以及含气性预测结果确定所述目标工区的储层分布结果。

本说明书提供的所述方法的另一个实施例中，所述获取目标工区的泊松比数据，包括：

对目标工区的纵波或者转换波的叠前地震数据进行叠前地震反演处理，获得所述目标工区的泊松比数据。

本说明书提供的所述方法的另一个实施例中，所述对所述目标工区进行含气性预测，包括：

对目标工区的地震波能量随频率的变化特征进行分析处理，确定地震波的低频能量异常反射特征；

根据所述泊松比数据以及地震波的低频能量异常反射特征对所述目标工区进行含气性预测。

另一方面，本说明书实施例还提供一种碳酸盐岩致密薄储层分布确定装置，所述装置包括：

构造构建模块，用于基于预先构建的速度场，利用大比例尺小网格变速成图方法确定目标工区的今构造高点；

沉积相带确定模块，用于对所述今构造高点对应目标层段进行层位拉平恢复古构造，并在恢复的古构造的基础上进行地震相分析处理，获得沉积相带分布结果；

储层分布确定模块，用于根据所述沉积相带分布结果确定所述目标工区的储层分布结果。

另一方面，本说明书实施例还提供一种碳酸盐岩致密薄储层分布确定设备，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：

另一方面，本说明书实施例还提供一种碳酸盐岩致密薄储层分布确定系统，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个实施例所述方法的步骤。

本说明书一个或多个实施例提供的一种碳酸盐岩致密薄储层分布确定方法、装置及系统，可以基于预先构建的速度场，利用变速大比例尺小网格变速成图技术准确确定现今微幅度构造高点及相应的构造形态。并可以进一步对今构造高点对应的目标层段进行地震层拉平恢复古构造，利用古地貌构造及地震相综合分析准确确定有利沉积相带分布，为寻找储层高产富集区提供依据。利用本说明书各实施例，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

附图说明

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

图1为本说明书提供的一种碳酸盐岩致密薄储层分布确定方法实施例的流程示意图；

图2为本说明书提供的另一种碳酸盐岩致密薄储层分布确定方法实施例的流程示意图；

图3为本说明书提供的一个实施例中的采用叠前同时反演获得泊松比属性剖面图；

图4为本说明书提供的一个实例中的地震-地质层位标定面示意图；

图5为本说明书提供的一个实例中的桃*井过井剖面及奥陶系顶部构造图；

图6为本说明书提供的一个实例中的桃*过井剖面及Tc2层平均方根振幅切片示意图；

图7为本说明书提供的一个实例中的三维Line测线偏移、泊松比反演、流体活动性属性剖面图；

图8是本说明书提供的一种碳酸盐岩致密薄储层分布确定装置实施例的模块结构示意图；

图9是本说明书提供的另一种碳酸盐岩致密薄储层分布确定装置实施例的模块结构示意图。

具体实施方式

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

碳酸盐岩储层的储集空间主要为溶孔、溶洞和微细裂缝，储层单层有效厚度小(小于5m)，非均质性强烈等特点。通过对大量已知井地质规律进行分析归纳，主力储层的厚度、空间展布形态等情况均受侵蚀地貌控制，即在侵蚀量大的潜沟中，因主力储层缺失而不含气；相反，地层保存较全的潜台区含气好。但碳酸盐岩储层波阻抗差异较小、地震反射弱，岩性横向变化大，白云岩储存及有效储层均为薄储层，含气性预测困难，目标层段成藏主控因素认识不清。

相应的，本说明书实施例提供了一种碳酸盐岩致密薄储层分布确定方法，可以预先构建速度场，并利用变速大比例尺小网格变速成图技术准确确定现今微幅度构造高点以及相应的构造形态。并可以进一步对今构造高点对应的目标层段进行地震层拉平恢复古构造，并基于古构造高点及地震相分析准确确定有利沉积相带分布，为寻找储层高产富集区提供依据。利用本说明书各实施例，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

图1是本说明书提供的所述一种碳酸盐岩致密薄储层分布确定方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。

具体的一个实施例如图1所示，本说明书提供的碳酸盐岩致密薄储层分布确定方法的一个实施例中，所述方法可以包括：

S2：基于预先构建的速度场，利用大比例尺小网格变速成图方法确定目标工区的今构造高点。

可以先合成地震记录，并进行精细地震地质层位标定，确定主要标志层反射及目标层段小层反射，以进行目标层小层精细解释。在层位精细解释的基础上，可以进一步通过二维十字线控制、三维地震资料获得相对地震解释T₀图(等时间构造图)，在时间域T₀图基础上构建的速度场变速成图获得深度域构造图。然后，可以在深度域构造图上落实现今构造高点，确定今构造高点位置。

再通过大比例尺小网格变速成图技术准确预测该今构造高点的幅度、面积等构造形态。所述大比例尺小网格变速成图技术可以是指在大比例尺、小网格尺寸的设定下，对时间域的地震解释图形进行变速成图，获得深度域构造图。

利用二维十字线控制、三维精细小层目标解释，以及大比例尺小网格变速成图技术精细构造成图，相比传统方法确定的构造图，构造成图精度更高，对微幅度构造高点特征反映更清楚。构造等值线可以由原来25米间隔变为5米间隔，对于5-20米的微幅度刻画更加精确，有利于准确圈定下古生界碳酸盐岩储层薄储层微幅度构造圈闭及面积。

所述合成地震记录可以包括用声波测井或垂直地震剖面资料经过人工合成转换成的地震记录(地震道)。合成地震记录从井出发，可以准确地标定出主要地质界面在地震剖面上所对应的反射同相轴(波峰或波谷)。通过层位标定可以将地震资料与钻井资料进行相互关联，使二者之间建立一个准确的对应关系，使得解释的目的层更明确。一些实施方式中，可以根据声波和密度测井曲线计算得到反射系数，将反射系数与提取的地震子波进行褶积得到合成地震记录。

本说明书的一个实施例中，可以基于合成地震记录进行层位标定，确定标志层反射及目标层段小层反射特征，根据所述标志层反射及目标层段小层反射特征进行地震层位小层段解释，获得层位解释结果。进一步根据所述层位解释结果，利用大比例尺小网格变速成图方法确定今构造高点。

可以对已完钻井合成记录进行精细地震地质层位标定，确定主要标志层反射及目标层段小层反射，指导目标层小层精细解释。由于碳酸盐岩储层非均质性强，层系多、储层薄、横向变化大，成藏主要受岩性、构造控制。但由于碳酸盐岩之间阻抗差异较小，在地震剖面上反射较弱，局部空白或杂乱反射，横向上也难以连续追踪。对主要目标层进行地震层位小层段精细解释，有利于明确目标层段的平面展布特征，为进一步精细落实构造特征奠定了坚实的基础。

本说明书的另一个实施例中，可以根据模型层析法构建速度场。传统的速度场一般采用测井层速度或者地震数据迭加速度DIX公式转换为平均速度，本说明书一些实施方式中，可以进一步采用模型层析法构建速度场，以提高构建的速度场的准确性，进而提高反射层构造形态及埋深解释的准确性。所述根据模型层析法构建速度场可以包括：在原有DIX速度转换的基础上，考虑到地下地质体横向差异，分别从浅、中、深三个不同的层段进行层位横向约束与控制，构建速度场。

一些实施方式中，可以通过钻井、地震相结合，通过钻井层位标定，完成地震资料浅、中、深层段地震层位解释，将目标工区纵向待分析区域划分成浅、中、深三个层段。然后，对浅、中、深三个层段分别进行DIX速度转换，构建速度场。地层纵向上地层速度变化较为明显，浅、中、深层段速度差异较为明显。通过预先划分浅、中、深层段，然后在不同层段中分别构建速度场，并进行层位约束，可以使得构建的速度场更符合地质体空间变化，进而提高获得的构造图的精准性。

S4：对所述今构造高点对应目标层段进行层位拉平恢复古构造，并在恢复的古构造的基础上进行地震相分析处理，获得沉积相带分布结果。

确定所述今构造高点对应目标层段，然后，可以利用研究区内二维、三维地震资料，利用层位拉平技术恢复古构造。并可以在古地貌构造分析的基础上，对目标层段的地震属性、地震波形等地震相变化进行分析，获得沉积相带分布结果。

一些实施方式中，可以综合地层厚度、测井相特征明确目标层段沉积特征及反射特征，并结合上述步骤中进行的层位解释结果确定目标层段附近的地震标志层，然后，可以基于所述地震标志层进行层位拉平恢复古构造。然后，可以在古构造基础上，分析地震波形特征及瞬时振幅等地震属性横向变化差异，划分研究区沉积相带分布，确定有利沉积相带。经过分析可知，有利地震相可以表现为以下两方面特征：

①沿目标层段顶界地震层拉平后表现为相对的构造位置高，即为古构造高。一般碳酸盐岩台地古构造高幅度较小(5米-20米)，古构造一方面本身有利于高能碳酸盐岩沉积形成高孔高渗储层，另一方面古构造高有利于后期风化淋滤改造使储层物性变好。

②拉平层后目标层段地震波形特征横向具有一定变化，振幅增强，地震外形包络线时差变大，地震属性瞬时振幅等地震属性表现异常，一般呈现丘状反射。

相应的，本说明书的一个实施例中，可以根据古构造高点以及地震属性异常特征确定有利地震相带，从而确定有利沉积相带。基于古地貌分析以及地震波形特征及瞬时振幅等地震属性横向变化差异特征分析，可以更加准确的确定研究区内的有利沉积相带分布。

S6：根据所述沉积相带分布结果确定所述目标工区的储层分布结果。

可以进一步根据上述步骤S4中确定的有利沉积相带确定出目标工区的有效储层分布区。碳酸盐岩储层分布往往受到沉积相的控制，利用上述实施例准确确定有利沉积相带分布区，可以提高储层高产富集区确定的准确性。当然，具体实施时，也获取已有储层分布分析结果，在已有储层分布结果上进一步结合有利沉积相分布结果，确定最终的有效储层分布区，为寻找高产富集区奠定了良好基础。

图2是本说明书提供的另一种碳酸盐岩致密薄储层分布确定方法实施例流程示意图。如图2所示，本说明书的另一个实施例中，所述方法还可以包括：

S8：基于泊松比数据对所述目标工区进行含气性预测；

S10：根据所述沉积相带分布结果以及含气性预测结果确定所述目标工区的储层分布结果。

可以基于泊松比数据对所述目标工区进行含气性预测，然后，综合步骤S4中的沉积相带分布结果以及含气性检测结果确定所述目标工区的储层分布结果。

具体实施时，可以在步骤S4中确定的有效储层分布区基础上，对所述有效储层分布区进行含气性预测。可以先利用纵波阻抗和横波阻抗识别岩性，然后，利用泊松比识别相应岩性储层的含气性。一些实施方式中，可以利用叠前弹性反演技术，获得弹性波阻抗，纵、横波阻抗，密度，泊松比等多种与岩性有关的参数体，从而反映储层及含气储层特征。所述叠前弹性反演技术可以包括利用纵波或转换波的叠前地震资料，按不同的入射角进行反演获得弹性波阻抗，纵、横波阻抗，密度，泊松比等多种与岩性有关的参数体。

经过对实测井的全波列测井曲线井岩石物理分析可知，白云岩与灰岩相比具有低纵波阻抗、高横波阻抗特征，气层具有低泊松比特征。通过叠前弹性参数反演获得泊松比数据识别含气性，可以提高相应岩性储层含气性预测的准确性。如图3所示，图3表示采用叠前同时反演获得泊松比属性剖面图，可以准确方便的利用图3中的泊松比属性剖面图检测白云岩储层的含气性。

因碳酸盐岩储层非均质性强，含气之后地震反射特征多变，目前对碳酸盐岩储层含气性预测尚未形成统一的、准确的认识，世界范围内碳酸盐岩储层含气性预测技术方法仍处于探索阶段。本说明书实施例，通过利用泊松比属性剖面图进行储层含气性预测，为储层含气性预测提供了一种准确高效的方法。

本说明书的另一个实施例中，可以对目标工区的地震波能量随频率的变化特征进行分析处理，确定地震波的低频能量异常反射特征，根据所述泊松比数据以及地震波的低频能量异常反射特征对所述目标工区进行含气性预测。通过进一步结合地震波的低频能量异常反射特征进行碳酸盐岩储层含气性预测，可以更加准确的确定高产富集区内的井位布设位置，以取得良好的勘探效果。

固体与流体之间相对位移的一个重要作用是引起双相介质中地震波能量的再分配，即同单相介质情况相比，地震波各个频率成分所占的比重将发生明显的变化。主要表现为地震波能量向低频方向移动，低频段地震振幅能量增强。因此，利用上述实施例的方案，通过采用泊松比反演与流体活动性属性预测结果相互约束验证，对所述有效储层分布区进行含气性预测，可以明显提高含气性检测的准确性，降低多解性。

相应的，本说明书还提供一个应用上述实施例方案的具体实例。具体的，以下桃*井在下古生界马家沟组马五4段白云岩储层为例，详细说明本说明书上述实施例提供的方案。

鄂尔多斯盆地古隆起东侧奥陶系盐下作为奥陶系深层碳酸盐岩的天然气成藏新领域，成藏地质条件有利，具备一定规模的勘探潜力，目前，盐下勘探已经取得重要进展。但碳酸盐岩储层波阻抗差异较小、地震反射弱，岩性横向变化大，白云岩储存及有效储层薄(储层累计厚度5.0-15.5m，单层2.0-9.5m，气层1.0-5.0m)，致使储层及含气性预测困难，目标层段成藏主控因素认识不清。

1、精细构造解释，寻找有利成藏的今构造高点。

由于碳酸盐岩储层非均质性强，阻抗差异较小，在地震剖面上反射较弱，局部空白或杂乱反射，横向上也难以连续追踪。因此，需要进行精细的层位标定，确定各主要目标层小层空间分布，全区共标定主要地震层位有Tc2、TC、To1m5⁵，通过精细解释主要目的层地震层位。

图4表示研究区内地震-地质层位标定面示意图。如图4所示，该研究区内古生界主要存在4个主要地震反射界面。

T_C2：盆地主要的区域标志层之一，为太原组与本溪组顶部厚煤层之间的反射。该反射同相轴为强振幅波峰反射，全区稳定，可以连续追踪对比，代表地层界面为太原组与本溪组地层分界面。

T_O1m₅₅：相当于奥陶系马家沟组第五段五亚段(马五5)底部附近反射，能量中～强。鄂尔多斯盆地下古生界奥陶系马家沟组顶部风化壳主要标志层，是马五5段灰岩与马五6段顶部含泥、膏、盐白云岩之间的反射，主要表现为波峰反射，局部由于岩性变化会形成波谷反射。

T_O1m₄：相当于马家沟组马四段底部附近反射，能量中～强，是马四段底部白云岩与马三段顶部盐岩之间的反射，其反射同相轴为强波峰。

T_O1m₂：相当于马家沟组底部附近反射，能量中～强，是马二段底部白云岩与马一段顶部盐岩之间的反射，由于马一段地层厚度薄，横向变化大，部分地区缺失，故其反射定位马家沟组底部反射，同相轴为强波峰。

在层位精细解释的基础上，可以进一步通过二维十字线控制、三维地震资料通过变速成图技术落实微幅度构造，如图5所示。图5表示桃*井过井剖面及奥陶系顶部构造图。分析图5可知，桃*井马五4段今构造高，奥陶系顶部反射层(TC)在井点处存在一个6ms构造幅度(图5左所示)。进一步通过变速大比例尺小网格变速成图技术预测该构造形态，构造幅度为12米，构造面积为8.3km²(图5右所示)。钻前地震预测Tc海拔为-2081米，实钻为-2080.7米，预测与实钻结果高度吻合。

2、沉积相结合岩溶古地貌分析，预测古构造高点，预测有效储层分布。

通过选取目标层段附近地震标志层，该区Tc2与马五5段为鄂尔多斯盆地内主要标志层，分别代表了上古生界太原组底界与下古生界奥陶系马家沟组马五5段反射。图6表示桃*过井剖面及Tc2层(1850ms)平均方根振幅切片示意图。如图6所示，通过地震目标层段顶(Tc2反射层)或底(To1m5⁵反射层)层层拉平之后目标层段地震反射表现为丘状结构，横向变化存在一定变化，井点及其附近多出一套地震相位，存在明显地震相位加宽，为典型的古构造高(图6左所示)。

古构造高有利于高能碳酸盐岩沉积，又有利于后期风化淋滤改造使储层物性。幅度为5-25米，面积5.5km²(图6右所示)。钻前预测该井点奥陶系顶部层位为马五3³，实钻奥陶系顶部层位为马五3³；预测马五4储层厚度发育，实钻马五4段储层发育，储层厚度为12.1米，高孔渗储(孔隙度大于8％，渗透率大于1.1MD)层共发育4段累计厚度为9.2米。

因此，通过岩溶古地貌结合地震相预测有效储层分布，为寻找高产富集区奠定了良好基础。

3、基于泊松比反演控制下流体活动性属性预测碳酸盐岩储层含气性。

由上述实施例可知，泊松比反演能有较好识别碳酸盐岩储层含气性，而流体活动性利用双相介质中地震波场的传播理论进行烃类检测可以取得良好效果。因此，采用泊松比反演与流体活动性属性预测结果相互约束验证，可以提高碳酸盐岩储层含气性预测结果。

在实际泊松比反演中，可以将多个分偏移距叠加资料(至少3个)，在具有横波测井的数据和地质层位的约束下，通过全波动方程反演出储层的纵波阻抗、横波阻抗、密度、与泊松比等弹性参数。图7表示三维Line测线偏移、泊松比反演、流体活动性属性剖面图。如图7所示，图7最上边的图表示三维Line测线偏移确定现今构造高点的示意图，图7中间的图为过桃*井钻前测线的叠前同时反演泊松比剖面图，泊松比剖面上马五4段表现为低值异常，反映此段储层含气性较好。实钻马五4段储层12.1m，测井解释气层3.9m，含气层4.0m，为Ⅰ类井。

目前，国内外多数地震属性预测主要是基于振幅属性预测。难以有效利用地震频率域信息进行含气性检测。本次采用自主研发GeoEast软件进行频率域含气性属性预测。GeoEast软件具有强大的多域(振幅类、频率类、统计类)属性提取和分析子系统。图7中的最下边的图为桃*井钻前流体活动属性含气性检测剖面，马五4段目的层段在流体活动性表现为频谱低频段高异常，反映储层流体活动性较好。实钻马五4段储层12.1m，测井解释气层3.9m，含气层4.0m，马五4段试气获日产百万方的高产工业气流。

因此，采用泊松比及流体活动性属性共同控制下含气性检测可以提高碳酸盐岩储层含气性预测的准确性，降低多解性，取得良好的勘探效果。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。

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

本说明书一个或多个实施例提供的一种碳酸盐岩致密薄储层分布确定方法，可以基于预先构建的速度场，利用变速大比例尺小网格变速成图技术准确确定现今微幅度构造高点以及相应的构造形态。并可以进一步对今构造高点对应的目标层段进行地震层拉平恢复古构造，利用古地貌构造及地震相综合分析准确确定有利沉积相带分布，为寻找储层高产富集区提供依据。利用本说明书各实施例，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

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

构造构建模块102，可以用于基于预先构建的速度场，利用大比例尺小网格变速成图方法确定目标工区的今构造高点；

沉积相带确定模块104，可以用于对所述今构造高点对应目标层段进行层位拉平恢复古构造，并在恢复的古构造的基础上进行地震相分析处理，获得沉积相带分布结果；

储层分布确定模块106，可以用于根据所述沉积相带分布结果确定所述目标工区的储层分布结果。

利用上述实施例中的方案，可以提高有效储层分布区确定的准确性，为寻找高产富集区奠定了良好基础。

本说明书的另一个实施例中，所述构造构建模块102可以包括：

速度场构建单元，可以用于根据模型层析法构建速度场。

利用上述实施例中的方案，使得最终建立的平均速度场更加精确，进一步提高有效储层分布区确定的准确性。

层位标定单元，可以用于基于合成地震记录进行层位标定，确定标志层反射及目标层段小层反射特征；

层位解释单元，可以用于根据所述标志层反射及目标层段小层反射特征进行地震层位小层段解释，获得层位解释结果；

构造构建单元，可以用于根据所述层位解释结果，利用大比例尺小网格变速成图方法确定目标工区的今构造高点。

利用上述实施例中的方案，可以进一步提高层位精细解释结果，进而提高构造图对地下构造刻画的精确性。

本说明书的另一个实施例中，所述沉积相带确定模块104可以包括：

地震相确定单元，可以用于在恢复的古构造的基础上进行地震相分析处理，确定古构造高点以及地震相异常特征；

沉积相带确定单元，可以用于根据所述古构造高点以及地震相异常特征确定有利沉积相带。

利用上述实施例中的方案，利用古地貌结合地震相预测有效储层分布，可以进一步提高有效储层分布区确定的准确性。

图9表示说明书提供的一种碳酸盐岩致密薄储层分布确定装置实施例的模块结构示意图。如图9所示，本说明书的另一个实施例中，所述装置还可以包括含气性检测模块105，其中，

所述含气性检测模块105可以用于获取目标工区的泊松比数据，根据所述泊松比数据对所述目标工区进行含气性预测；

相应的，所述储层分布确定模块106可以用于根据所述沉积相带分布结果以及含气性预测结果确定所述目标工区的储层分布结果。

利用上述实施例中的方案，利用气层泊松比特征进行储层含气性检测，可以提高储层含气性检测的准确性，进而提高有效储层分布区确定的准确性。

本说明书的另一个实施例中，所述含气性检测模块105可以包括：

泊松比数据确定单元，可以用于对目标工区的纵波或者转换波的叠前地震数据进行叠前地震反演处理，获得所述目标工区的泊松比数据。

本说明书的另一个实施例中，所述含气性检测模块105还可以包括：

地震波能量特征确定单元，可以用于对目标工区的地震波能量随频率的变化特征进行分析处理，确定地震波的低频能量异常反射特征；

含气性检测模块单元，用于根据所述泊松比数据以及地震波的低频能量异常反射特征对所述目标工区进行含气性预测。

利用上述实施例中的方案，可以明显提高含气性检测的准确性，降低多解性，取得良好的勘探效果。

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

本说明书一个或多个实施例提供的一种碳酸盐岩致密薄储层分布确定装置，可以基于预先构建的速度场，利用变速大比例尺小网格变速成图技术准确确定现今微幅度构造高点以及相应的构造形态。并可以进一步对今构造高点对应的目标层段进行地震层拉平恢复古构造，利用古地貌构造及地震相综合分析准确确定有利沉积相带分布，为寻找储层高产富集区提供依据。利用本说明书各实施例，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

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

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

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

上述实施例所述的一种碳酸盐岩致密薄储层分布确定设备，可以基于预先构建的速度场，利用变速大比例尺小网格变速成图技术准确确定现今微幅度构造高点以及相应的构造形态。并可以进一步对今构造高点对应的目标层段进行地震层拉平恢复古构造，利用古地貌构造及地震相综合分析准确确定有利沉积相带分布，为寻找储层高产富集区提供依据。利用本说明书各实施例，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

本说明书还提供一种碳酸盐岩致密薄储层分布确定系统，所述系统可以为单独的碳酸盐岩致密薄储层分布确定系统，也可以应用在多种石油勘探系统中。所述的系统可以为单独的服务器，也可以包括使用了本说明书的一个或多个所述方法或一个或多个实施例装置的服务器集群、系统(包括分布式系统)、软件(应用)、实际操作装置、逻辑门电路装置、量子计算机等并结合必要的实施硬件的终端装置。所述碳酸盐岩致密薄储层分布确定系统可以包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个或者多个实施例中所述方法的步骤。

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

上述实施例所述的一种碳酸盐岩致密薄储层分布确定系统，可以基于预先构建的速度场，利用变速大比例尺小网格变速成图技术准确确定现今微幅度构造高点以及相应的构造形态。并可以进一步对今构造高点对应的目标层段进行地震层拉平恢复古构造，利用古地貌构造及地震相综合分析准确确定有利沉积相带分布，为寻找储层高产富集区提供依据。利用本说明书各实施例，可以进一步提高碳酸盐岩有效储层分布区确定的准确性。

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

尽管本说明书实施例内容中提到的层位解释、地震属性相分析等获取、定义、交互、计算、判断等操作和数据描述，但是，本说明书实施例并不局限于必须是符合标准数据模型/模板或本说明书实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书的可选实施方案范围之内。

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

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

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

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

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

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

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

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

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

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

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

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

Claims

1.一种碳酸盐岩致密薄储层的井位布设方法，其特征在于，包括：

基于合成地震记录对目标工区进行层位标定，确定标志层反射及目标层段小层反射特征；

基于所述层位解释结果，对目标工区纵向划分层段，对划分的层段分别进行DIX速度转换，构建速度场；

基于构建的速度场，利用大比例尺小网格变速成图方法确定目标工区的今构造高点；其中，大比例尺小网格变速成图方法是指在大比例尺、小网格尺寸的设定下，对时间域的地震解释图形进行变速成图得到深度域构造图的方法；

根据所述今构造高点对应目标层段的沉积特征、反射特征以及所述层位解释结果确定目标层段附近的地震标志层；基于所述地震标志层对所述今构造高点对应目标层段进行层位拉平恢复古构造；

在恢复的古构造的基础上进行地震相分析处理，确定古构造高点以及地震相异常特征；根据古构造高点以及地震属性异常特征确定有利沉积相带；根据确定的所述有利沉积相带确定所述目标工区的有效储层分布区；

对所述有效储层分布区的纵波或者转换波的叠前地震数据进行叠前地震反演处理，获得所述有效储层分布区的泊松比属性剖面图；

对有效储层分布区的地震波能量随频率的变化特征进行分析处理，确定地震波的低频能量异常反射特征；

根据所述泊松比属性剖面图以及地震波的低频能量异常反射特征对所述有效储层分布区进行含气性预测；

根据所述含气性预测结果确定所述有效储层分布区的井位布设位置。

2.一种碳酸盐岩致密薄储层的井位布设装置，其特征在于，所述装置包括：

构造构建模块，用于基于合成地震记录对目标工区进行层位标定，确定标志层反射及目标层段小层反射特征；根据所述标志层反射及目标层段小层反射特征进行地震层位小层段解释，获得层位解释结果；基于所述层位解释结果，对目标工区纵向划分层段，对划分的层段分别进行DIX速度转换，构建速度场；基于构建的速度场，利用大比例尺小网格变速成图方法确定目标工区的今构造高点；其中，大比例尺小网格变速成图方法是指在大比例尺、小网格尺寸的设定下，对时间域的地震解释图形进行变速成图得到深度域构造图的方法；

沉积相带确定模块，用于根据所述今构造高点对应目标层段的沉积特征、反射特征以及所述层位解释结果确定目标层段附近的地震标志层；基于所述地震标志层对所述今构造高点对应目标层段进行层位拉平恢复古构造；在恢复的古构造的基础上进行地震相分析处理，确定古构造高点以及地震相异常特征；根据古构造高点以及地震属性异常特征确定有利沉积相带；

储层分布确定模块，用于根据确定的所述有利沉积相带确定所述目标工区的有效储层分布区；以及，对所述有效储层分布区的纵波或者转换波的叠前地震数据进行叠前地震反演处理，获得所述有效储层分布区的泊松比属性剖面图；对有效储层分布区的地震波能量随频率的变化特征进行分析处理，确定地震波的低频能量异常反射特征；根据所述泊松比属性剖面图以及地震波的低频能量异常反射特征对所述有效储层分布区进行含气性预测；根据所述含气性预测结果确定所述有效储层分布区的井位布设位置。

3.一种碳酸盐岩致密薄储层的井位布设系统，其特征在于，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现权利要求1所述方法的步骤。