本发明涉及油气勘探，尤其涉及一种新近系岩性圈闭潜力快速优选方法。
背景技术：
1、在渤海湾含油气盆地中，渤海海域新近系油气藏是勘探评价的主要目标。勘探早期，从油藏类型上看，它们均以构造油气藏为主，随着渤海海域油气勘探程度逐渐提高，新近系的构造圈闭越来越少，圈闭的规模也越来越小，所发现油气藏的规模也将越来越小。因此，渤海海域下一步的勘探思路必须从单一构造圈闭转到岩性圈闭领域，岩性油气藏逐步成为最具潜力、最现实的油气勘探领域。
2、近几年，渤海新近系岩性勘探进程日渐加深，尤其新近系河道型复合砂体贡献了近两年大部分的储量发现。在勘探评价的过程中，搜索浅层岩性潜力目标区十分关键，在以往一直是通过人为逐层搜索，对全区潜力区缺乏整体把控，不能够高效锁定多套含油气砂体叠置的优势潜力区，同时，浅层岩性勘探目标的搜索，是在制作严格的地层切片后，先分层寻找规模性岩性砂体进行追踪解释，再根据追踪完毕之后的砂体叠置关系进行井位设计，选择最优靶点，这样缺乏先期砂体发育或油气富集的区域概念，导致前期工作目标性差。因此有必要寻找一种快速优选潜力区的方法，为浅层岩性圈闭的高效勘探评价提供支持。
技术实现思路
1、有鉴于此，本发明旨在提出一种新近系岩性圈闭潜力快速优选方法，其能够快速应用于勘探评价过程中有利岩性目标的搜索，为前期靶区优选效率的提高发挥了作用，为浅层岩性圈闭的高效勘探评价提供支持。
2、为达到上述目的，本发明的技术方案是这样实现的：一种新近系岩性圈闭潜力快速优选方法，包括以下步骤：
3、第一步，将原始地震数据和顶底层位数据转成标准格式；
4、第二步，通过顶、底两个层位制作严格的地层切片；
5、第三步，提取每一层的表征储层属性或含油气属性；
6、第四步，通过读取表征目的层储层属性或含油气属性，自动搜索各小层内有利区域；
7、第五步，针对多个小层有利区域进行叠置分析，自动寻找叠置较好、有利于后续井位设计的潜力区，精准锁定目标砂体，为后续岩性砂体的追踪解释提供支持。
8、进一步的，第一步包括，对地震数据进行预处理，对于segy地震数据，隐藏道头信息，将其转化为时间、联络线数、主测线数的三维数据体，同时通过显示地震数据多个剖面检查转换后的地震数据是否正确；对于软件中导出的层位数据是根据精细的井震标定进行解释的目的层的顶界和底界，解释密度为12.5m×25m，并要进行插值和平滑，将其导出后转化为以inline号为横坐标，cdp号为纵坐标，其内的每一个点表示对应的时间值，层位数据仅显示顶、底两个大层即可。
9、进一步的，第二步具体为，在第一步预处理之后的地震数据基础上，基于一定的沉积模式，用顶底界面作为约束层，在顶、底两个大层之间通过等比例的内插得到位于两个界面之间的切片，以顶底界面的时间厚度除以顶、底界面之间的间隔数作为两个切片之间的厚度，如公式1所示，随后每个切片所处的时间位置可以通过公式2计算得到，该地层切片能够真实地反映地层产状；
10、
11、s(j,k,i0＝t顶(j,k)+i×δt(j,k) (2)
12、其中，s(j,k,i)为每个切片所处的时间位置，j表示地震主测线号，k表示地震联络线号，δt表示两个切片之间的厚度，(j,k)表示地震数据上第j主测线与第k联络线交点位置处的采样点，i表示第i个切片，n表示切片的个数，t顶表示顶界面，t底表示底界面。
13、进一步的，所述第三步中，由于地震振幅属性切片与地质体外形以及含油气属性息息相关，因此提取每一层的表征储层属性或含油气属性时，以地层切片上的采样点作为中心点进行纵向开窗，并提取时窗内的最小振幅属性。
14、进一步的，所述第四步中，自动搜索各小层内有利区域是指通过输入最大井斜以及属性阈值作为约束，首先筛选出属性阈值之上的区域，在满足井轨迹合理的前提下，尽可能多地经过振幅大的区域并求取振幅和作为最终结果；在此步骤中，可以选择单个切片结果，也可以选择包含两两、三三或更多切片组合下的结果，并在下一步中将可选择结果进行归一化或不归一化，并进行排序。
15、进一步的，所述第五步中，自动寻找叠置较好、有利于后续井位设计的潜力区是指在上一步得到的归一化或者未归一化结果的基础上，通过阈值点数和阈值百分数两个参数进行约束，阈值点数是指优选出潜力区范围内所有切片属性阈值后的下限，阈值百分数是指切片上符合条件的面积占时窗总面积的百分比；这两个参数决定了两种确定优势潜力区的方式，参数的具体选取要根据实际工区来定，在优选时，可选择全部切片参与优选潜力区计算，也可根据目标区的实际情况选择主要切片参与计算，通过设置搜索步长以及时窗大小最终得到优势潜力区的三点坐标并在切片平面属性图上投影；具体地，通过定义最小振幅属性阈值统计一个面元内大于该阈值的采样点个数，随后定义采样点阈值对所有切片进行筛选，最后按照两种定义方式确定优势潜力区并对潜力区进行排序，两种定义方式如下所示：
16、定义方式一：所有切片上符合阈值下的点取平均，并求取平均值与面元上所有采样点的比值，当大于某一值时就定义为需要的潜力区：
17、a＝sort_number/n (3)
18、p1＝a/sum_sample (4)
19、其中，sort_number为所有切片上符合阈值下的点；n为切片的个数；sum_sample为某一切片上一个面元内所有采样点数；a为所有切片上符合阈值下的点的平均值；p1为所有切片上符合阈值下的点的平均值与面元上所有点的比值；
20、定义方式二：求取所有切片上符合阈值下的采样点所占面积求和与面元总面积的比值，当大于某一值时就定义为需要的潜力区：
21、area＝sort_number×bin1×bin2 (5)
22、p2＝area/sum_area (6)
23、其中，sort_number为所有切片上符合阈值下的点；bin1为cdp方向长度；bin2为inline方向长度；sum_area为面元总面积；area为所有切片上符合阈值下的采样点所占的面积；p2为所有切片上符合阈值下的采样点所占的面积与面元总面积的比值。
24、相对于现有技术，本发明所述的一种新近系岩性圈闭潜力快速优选方法具有以下优势：
25、(1)本发明所述方法克服了常规人为逐层搜索方式确定潜力区的不足，通过岩性勘探优势潜力区方法可集成为一体化流程，提高了计算效率，而且有效避免了人为搜索时产生的误差，能够快速应用于勘探评价过程中有利岩性目标的搜索；
26、(2)本发明所述方法自动寻找叠置较好、有利于后续井位设计的潜力区，尤其对于多层叠置砂体目标区具有较好的应用效果，同时最终输出三点坐标以及图形化显示可以更直观地展示有利目标区。
1.一种新近系岩性圈闭潜力快速优选方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种新近系岩性圈闭潜力快速优选方法，其特征在于：第一步包括，对地震数据进行预处理，对于segy地震数据，隐藏道头信息，将其转化为时间、联络线数、主测线数的三维数据体，同时通过显示地震数据多个剖面检查转换后的地震数据是否正确；对于软件中导出的层位数据是根据精细的井震标定进行解释的目的层的顶界和底界，解释密度为12.5m×25m，并要进行插值和平滑，将其导出后转化为以inline号为横坐标，cdp号为纵坐标，其内的每一个点表示对应的时间值，层位数据仅显示顶、底两个大层即可。
3.根据权利要求2所述的一种新近系岩性圈闭潜力快速优选方法，其特征在于：第二步具体为，在第一步预处理之后的地震数据基础上，基于一定的沉积模式，用顶底界面作为约束层，在顶、底两个大层之间通过等比例的内插得到位于两个界面之间的切片，以顶底界面的时间厚度除以顶、底界面之间的间隔数作为两个切片之间的厚度，如公式1所示，随后每个切片所处的时间位置可以通过公式2计算得到，该地层切片能够真实地反映地层产状；
4.根据权利要求1所述的一种新近系岩性圈闭潜力快速优选方法，其特征在于：所述第三步中，由于地震振幅属性切片与地质体外形以及含油气属性息息相关，因此提取每一层的表征储层属性或含油气属性时，以地层切片上的采样点作为中心点进行纵向开窗，并提取时窗内的最小振幅属性。
5.根据权利要求1所述的一种新近系岩性圈闭潜力快速优选方法，其特征在于：所述第四步中，自动搜索各小层内有利区域是指通过输入最大井斜以及属性阈值作为约束，首先筛选出属性阈值之上的区域，在满足井轨迹合理的前提下，尽可能多地经过振幅大的区域并求取振幅和作为最终结果；在此步骤中，可以选择单个切片结果，也可以选择包含两两、三三或更多切片组合下的结果，并在下一步中将可选择结果进行归一化或不归一化，并进行排序。
6.根据权利要求5所述的一种新近系岩性圈闭潜力快速优选方法，其特征在于：所述第五步中，自动寻找叠置较好、有利于后续井位设计的潜力区是指在上一步得到的归一化或者未归一化结果的基础上，通过阈值点数和阈值百分数两个参数进行约束，阈值点数是指优选出潜力区范围内所有切片属性阈值后的下限，阈值百分数是指切片上符合条件的面积占时窗总面积的百分比；这两个参数决定了两种确定优势潜力区的方式，参数的具体选取要根据实际工区来定，在优选时，可选择全部切片参与优选潜力区计算，也可根据目标区的实际情况选择主要切片参与计算，通过设置搜索步长以及时窗大小最终得到优势潜力区的三点坐标并在切片平面属性图上投影；具体地，通过定义最小振幅属性阈值统计一个面元内大于该阈值的采样点个数，随后定义采样点阈值对所有切片进行筛选，最后按照两种定义方式确定优势潜力区并对潜力区进行排序，两种定义方式如下所示：
