CN107807407B - 一种油气区带有效性评价方法和装置 - Google Patents

Abstract

本发明提供的一种油气区带有效性评价方法和装置，所述油气区带有效性评价方法包括：获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；根据所述控因参数分布获取每个网格坐标点的成藏参数值；根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值获取每个所述网格坐标点的评价参数处理值；根据所述网格坐标点的评价参数处理值获取所述目的层段的区带评价值。本发明提供一种能提高评价符合率的油气区带有效性评价方法和装置。

Description

一种油气区带有效性评价方法和装置

技术领域

本发明涉及石油天然气勘探的区带评价方法技术领域，尤其涉及一种油气区带有效性评价方法和装置。

背景技术

油气区带有效性评价就是评价区带是否存在工业价值油气藏。油气区带的有效性评价是决定钻探成功率高低，并决定油气发现和储量增长速度的关键。

现有技术中的油气区带有效性评价方法包括：R因子主成分分析方法、多因素加权评价方法、多图叠合区带评价方法、区带半定量评价方法和基于资源量及经济性的评价方法。

R因子主成分分析方法无法做到对单个区带的评价，区带条件相差较大，评价结果符合率低。多因素加权评价方法的权系数人为确定，不确定性较大，没有考虑控因参数的决定作用。多图叠合区带评价方法无法做到定量评价，对于勘探程度较低的地区，无法准确获取被评价区的单因素图件，评价符合率较低。区带半定量评价方法中不同油气区带的控制因素参数值相差较大，制定的评价参数标准适用范围局限，同时对于勘探程度较低的地区无法获取准确的参数值，评价符合率低。基于资源量及经济性的评价方法资源量大小无法确定区带是否能够获得工业价值油气藏，区带在未钻探前很难评价其经济性，评价符合率不高。

现有油气区带评价技术，主要基于地质条件类比、数学模型计算、资源量评价、多因素叠合等评价方法为主，受人为经验和掌握资料程度影响很大，无法制定统一的区带有效性评价标准，即无法做到对区带有效性的统一评价，区带评价符合率较低。

发明内容

本发明的目的是提供一种能提高评价符合率的油气区带有效性评价方法和装置。

本发明的上述目的可采用下列技术方案来实现：一种油气区带有效性评价方法，包括：获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；根据所述控因参数分布计算每个网格坐标点的成藏参数值；根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值计算每个所述网格坐标点的评价参数处理值；根据所述网格坐标点的评价参数处理值计算所述目的层段的区带评价值；所述根据所述控因参数分布计算每个网格坐标点的成藏参数值，具体包括：根据所述控因参数获取每个所述网格坐标点的无量纲标志值；根据每个所述网格坐标点的所述无量纲标志值获取每个所述网格坐标点的成藏参数值；

根据以下公式获取每个所述网格坐标点的成藏参数值：

其中，AR—成藏参数值，当AR≥1时为有效成藏区，当AR<1时为非有效成藏区，无量纲；Para_i—所述无量纲参数值Sth_S,Ts_S,Da_S,Rd_S,F_S；n—Para_i的个数；Sth_S—盖层厚度标志值，无量纲；Ts_S—输导体系标志值，无量纲；Da_S—地层倾角标志值，无量纲；Rd_S—油气运移距离标志值，无量纲；F_S—断层标志值；

根据以下公式确定所述盖层厚度标志值：

其中，Sth_S—盖层厚度标志值，无量纲；Sth—目的层段的盖层厚度，m；Sth_limt—工业价值油气藏的盖层厚度下限值，m；

根据以下公式确定所述输导体系标志值：

其中，Ts_S—输导体系标志值，无量纲；H—目的层段的储层厚度，m；H_limt—工业价值油气藏的储层厚度下限值，m；UF—与烃源岩有效生烃范围内直接接触的不整合和/或断层；当存在与烃源岩有效生烃范围内直接接触的不整合和/或断层时，为1；当与烃源岩有效生烃范围内直接接触的不整合或断层均不存在时，为0；

根据以下公式确定所述地层倾角标志值：

Da_S＝sin(α)，

其中，Da_S—地层倾角标志值，无量纲；α—地层倾角，度；

根据以下公式确定所述油气运移距离标志值：

其中，Rd_S—油气运移距离标志值，无量纲；Rd—有效生烃范围等效半径，km；L_hg—有效生烃范围外距有效生烃范围边界的距离，km，当处于有效生烃范围之内时，为0；a—经验系数，当目的层段为常规油气时，取3，当目的层段为致密油气时，取1.2，当目的层段为页岩油气时，取1；

当断穿目的层段的盖层的断层为开启断层时，所述断层标志值为-1；当断穿目的层段的盖层的断层为封闭断层时，所述断层标志值为1；

根据以下公式计算每个所述网格坐标点的评价参数处理值：

其中，P_i—第i个所述网格坐标点的评价参数处理值；S_i—第i个所述评价参数值；S_{i_limt}—第i个所述评价参数下限值；n—所述评价参数的个数；

根据以下公式计算所述目的层段的区带评价值：

其中，V_Play—区带评价值；-1～1，当V_Play≥0时，为有效区带分布区，当V_Play<0时，为无效区带分布区；P_i—第i个所述网格坐标点的评价参数处理值；n—所述网格坐标点的个数。

作为一种优选的实施方式，根据以下步骤确定所述烃源岩有效生烃范围：获取工业价值油气藏的生油丰度下限值、生气丰度下限值；根据所述生油丰度下限值和所述生气丰度下限值确定所述烃源岩有效生烃范围。

作为一种优选的实施方式，根据以下公式确定所述f(P_i)：

其中，P_i—第i个所述网格坐标点的评价参数处理值。

作为一种优选的实施方式，所述获取所述目的层段的地层倾角分布，具体包括：根据目的层段的地球物理资料获取所述目的层段的构造；根据所述目的层段的构造获取所述目的层段的地层倾角分布。

作为一种优选的实施方式，根据以下步骤确定所述目的层段的地层倾角：获取烃源岩有效生烃中心与所述目的层段的顶面构造等值线的垂直线段；获取所述垂直线段的倾斜角；根据所述垂直线段与水平线之间的位置关系获取所述倾斜角的倾斜方向；根据所述倾斜角和所述倾斜方向获取所述地层倾角。

作为一种优选的实施方式，所述断层包括：断穿所述目的层段的盖层的第一断层、沟通所述目的层段的烃源岩有效生烃范围内烃源岩与储层的第二断层和沟通所述目的层段的不整合与储层的第三断层。

作为一种优选的实施方式，还包括：根据所述目的层段的测井资料和烃源岩岩心分析的有机碳含量，获取所述目的层段的所述有机碳含量分布；根据所述目的层段的烃源岩岩心分析的镜质体反射率，获取所述镜质体反射率分布；对所述有机碳含量分布和所述镜质体反射率分布进行盆地模拟运算，获取所目的层段的生油丰度分布和生气丰度分布。

作为一种优选的实施方式，所述获取工业价值油气藏的储层厚度的下限值，具体包括：根据工业价值油气藏的储层厚度和试油气资料，获取储层厚度的下限值，其中，所述工业价值油气藏的储层类型与所述目的层段的储层类型相同。

作为一种优选的实施方式，所述获取工业价值油气藏的孔隙度的下限值，具体包括：根据工业价值油气藏的孔隙度分析资料和试气资料，获取孔隙度的下限值。

作为一种优选的实施方式，所述获取工业价值油气藏的盖层厚度的下限值，具体包括：根据工业价值油气藏的盖层厚度，获取盖层厚度的下限值，其中，所述工业价值油气藏的盖层岩性与所述目的层段的盖层岩性一致。

一种油气区带有效性评价装置，包括：控因参数获取模块，用于获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；评价参数下限值获取模块，用于获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；成藏参数值计算模块，用于根据所述控因参数分布计算每个网格坐标点的成藏参数值；评价参数处理值计算模块，用于根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值计算每个所述网格坐标点的评价参数处理值；区带评价值计算模块，用于根据所述网格坐标点的评价参数处理值计算所述目的层段的区带评价值；所述成藏参数值计算模块用于根据所述控因参数获取每个所述网格坐标点的无量纲标志值；并根据每个所述网格坐标点的所述无量纲标志值获取每个所述网格坐标点的成藏参数值；

根据以下公式获取每个所述网格坐标点的成藏参数值：

其中，AR—成藏参数值，当AR≥1时为有效成藏区，当AR<1时为非有效成藏区，无量纲；Para_i—所述无量纲参数值Sth_S,Ts_S,Da_S,Rd_S,F_S；n—Para_i的个数；Sth_S—盖层厚度标志值，无量纲；Ts_S—输导体系标志值，无量纲；Da_S—地层倾角标志值，无量纲；Rd_S—油气运移距离标志值，无量纲；F_S—断层标志值；

根据以下公式确定所述盖层厚度标志值：

其中，Sth_S—盖层厚度标志值，无量纲；Sth—目的层段的盖层厚度，m；Sth_limt—工业价值油气藏的盖层厚度下限值，m；

根据以下公式确定所述输导体系标志值：

其中，Ts_S—输导体系标志值，无量纲；H—目的层段的储层厚度，m；H_limt—工业价值油气藏的储层厚度下限值，m；UF—与烃源岩有效生烃范围内直接接触的不整合和/或断层；当存在与烃源岩有效生烃范围内直接接触的不整合和/或断层时，为1；当与烃源岩有效生烃范围内直接接触的不整合或断层均不存在时，为0；

根据以下公式确定所述地层倾角标志值：

Da_S＝sin(α)，

其中，Da_S—地层倾角标志值，无量纲；α—地层倾角，度；

根据以下公式确定所述油气运移距离标志值：

其中，Rd_S—油气运移距离标志值，无量纲；Rd—有效生烃范围等效半径，km；L_hg—有效生烃范围外距有效生烃范围边界的距离，km，当处于有效生烃范围之内时，为0；a—经验系数，当目的层段为常规油气时，取3，当目的层段为致密油气时，取1.2，当目的层段为页岩油气时，取1；

当断穿目的层段的盖层的断层为开启断层时，所述断层标志值为-1；当断穿目的层段的盖层的断层为封闭断层时，所述断层标志值为1；

所述评价参数处理值计算模块用于根据以下公式计算每个所述网格坐标点的评价参数处理值：

其中，P_i—第i个所述网格坐标点的评价参数处理值；S_i—第i个所述评价参数值；S_{i_limt}—第i个所述评价参数下限值；n—所述评价参数的个数；

所述区带评价值计算模块用于根据以下公式计算所述目的层段的区带评价值：

其中，V_Play—区带评价值；-1～1，当V_Play≥0时，为有效区带分布区，当V_Play<0时，为无效区带分布区；P_i—第i个所述网格坐标点的评价参数处理值；n—所述网格坐标点的个数。

本申请提供的一种油气区带有效性评价方法的有益效果是：本申请实施方式所述的油气区带有效性评价方法根据网格坐标点的评价参数、评价参数下限值和成藏参数值计算每个网格坐标点的评价参数处理值，从而确定目的层段的成藏范围。然后根据网格坐标点的评价参数处理值计算目的层段的区带评价值，如此将油气区带有效性评价定量处理。从而相比于现有技术，本申请所述的油气区带有效性评价方法能制定统一的油气区带有效性评价标准。从而能排除人为因素的影响。又由于采用已获得工业价值油气藏的评价参数下限值，从而提高了油气区带有效性的识别率，从而提高了油气区带有效性评价的符合率，如此能提高钻探成功率及发现有效区带的速度，降低油气勘探成本。

附图说明

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

图1本发明一个实施方式提供的油气区带有效性评价方法的流程图；

图2本发明一个实施方式提供的油气区带有效性评价方法的另一流程图；

图3本发明一个实施方式提供的油气区带有效性评价方法的另一流程图；

图4本发明一个实施方式提供的油气区带有效性评价方法的另一流程图；

图5本发明一个实施方式提供的油气区带有效性评价方法的另一流程图；

图6本发明一个实施方式提供的油气区带有效性评价方法的另一流程图；

图7本发明一个实施方式提供的目的层段的储层厚度分布特征图；

图8本发明一个实施方式提供的已获工业价值天然气区带的储层厚度下限值特征图；

图9本发明一个实施方式提供的目的层段的孔隙度分布特征图；

图10本发明一个实施方式提供的已获工业价值天然气区带的孔隙度下限值特征图；

图11本发明一个实施方式提供的目的层段的盖层厚度分布特征图；

图12本发明一个实施方式提供的已获工业价值天然气区带的盖层厚度下限值特征图；

图13本发明一个实施方式提供的目的层段的烃源岩镜质体反射率(Ro)分布特征图；

图14本发明一个实施方式提供的目的层段的烃源岩生气丰度分布特征图；

图15本发明一个实施方式提供的目的层段的烃源岩生油丰度分布特征图；

图16本发明一个实施方式提供的目的层段的区带评价结果特征图；

图17本发明一个实施方式提供的油气区带有效性评价装置的流程图。

具体实施方式

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

请参阅图1。本申请一种实施方式提供的油气区带有效性评价方法，其可以包括：S1：获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；S3：获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；S5：根据所述控因参数分布计算每个网格坐标点的成藏参数值；S7：根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值计算每个所述网格坐标点的评价参数处理值；S9：根据所述网格坐标点的评价参数处理值计算所述目的层段的区带评价值。

从以上技术方案可以看出：本申请实施方式所述的油气区带有效性评价方法根据网格坐标点的评价参数、评价参数下限值和成藏参数值获取每个网格坐标点的评价参数处理值，从而确定目的层段的成藏范围。然后根据网格坐标点的评价参数处理值获取目的层段的区带评价值，如此将油气区带有效性评价定量处理。从而相比于现有技术，本申请所述的油气区带有效性评价方法能制定统一的油气区带有效性评价标准。从而能排除人为因素的影响。又由于采用已获得工业价值油气藏的评价参数下限值，从而提高了油气区带有效性的识别率，从而提高了油气区带有效性评价的符合率，如此能提高钻探成功率及发现有效区带的速度，降低油气勘探成本。

如图1所示，在本实施方式中，S1：获取目的层段的控因参数分布，其中控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；

如图2所示，在一个具体的实施方式中，获取目的层段的地层倾角分布，具体包括：

S101：根据目的层段的地球物理资料获取目的层段的构造；

S103：根据目的层段的构造获取目的层段的地层倾角分布。

如图3所示，进一步地，根据以下步骤确定目的层段的地层倾角：

S1010：获取烃源岩有效生烃中心与所述目的层段的顶面构造等值线的垂直线段；

S1030：获取所述垂直线段的倾斜角；

S1050：根据所述垂直线段与水平线之间的位置关系获取所述倾斜角的倾斜方向；

具体地，当烃源岩有效生烃中心与目的层段的顶面构造等值线的垂直线，位于烃源岩有效生烃中心的水平线之上时，地层倾角为正；当烃源岩有效生烃中心与目的层段的顶面构造等值线的垂直线，位于烃源岩有效生烃中心的水平线之下时，地层倾角为负；当烃源岩有效生烃中心与目的层段的顶面构造等值线的垂直线，与位于烃源岩有效生烃中心的水平线重合时，地层倾角为0。

S1070：根据所述倾斜角和所述倾斜方向获取所述地层倾角。

在一个具体的实施方式中，断层包括：断穿目的层段的盖层的第一断层、沟通目的层段的烃源岩有效生烃范围内烃源岩与储层的第二断层和沟通目的层段的不整合与储层的第三断层。

在一个具体的实施方式中，与烃源岩有效生烃范围有效的不整合可以是烃源岩层与盖层之间的不整合。

在一个具体的实施方式中，获取目的层段的储层厚度分布可以通过采集目的层段的地球物理资料，获取目的层段的储层厚度分布。

在一个具体的实施方式中，获取目的层段的孔隙度分布可以通过采集目的层段的孔隙度分析资料和地球物理资料，获取目的层段的孔隙度分布。

在一个具体的实施方式中，获取目的层段的盖层厚度分布可以通过采集目的层段的地球物理资料，获取目的层段的盖层厚度分布。

如图4所示，在一个具体的实施方式中，获取目的层段的生油丰度分布和生气丰度分布，具体包括：

S105：采集目的层段内的测井资料和烃源岩岩心分析的有机碳含量，获取目的层段的有机碳含量分布；

S107：采集目的层段的烃源岩岩心分析的镜质体反射率，获取镜质体反射率分布；

S109：对有机碳含量分布和镜质体反射率分布进行盆地模拟运算，获取生油丰度分布和生气丰度分布。

如图1所示，在本实施方式中，S3：获取工业价值油气藏的评价参数下限值，其中，评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度。

在一个具体的实施方式中，获取工业价值油气藏的储层厚度的下限值具体包括：

根据工业价值油气藏的储层厚度和试油气资料，获取储层厚度的下限值，其中，工业价值油气藏的储层类型与目的层段的储层类型相同。

在一个具体的实施方式中，获取工业价值油气藏的孔隙度的下限值，具体包括：根据工业价值油气藏的孔隙度分析资料和试气资料，获取孔隙度的下限值。

在一个具体的实施方式中，获取工业价值油气藏的盖层厚度的下限值，具体包括：根据工业价值油气藏的盖层厚度，获取盖层厚度的下限值，其中，工业价值油气藏的盖层岩性与目的层段的盖层岩性一致。

在一个具体的实施方式中，工业价值油气藏的生油丰度下限值、生气丰度下限值用于确定烃源岩有效生烃范围。

利用本申请所述的油气区带有效性评价方法，对目的层段塔里木盆地库车山前地区的油气区带有效性进行评价。目的层段为三叠系、侏罗系烃源岩与白垩系巴什基奇克组碎屑岩储层组成的油气区带，目的层段圈闭以聚集天然气为主。

采集目的层段的地球物理资料，获取目的层段的构造、地层倾角分布、断层以及与烃源岩有效生烃范围有效的不整合。

如图7所示，采集目的层段的地球物理资料，获取目的层段的储层厚度分布。目的层段主要发育气藏为主，储层为裂缝型砂岩。

如图8所示，采集裂缝型砂岩储层已获得工业价值气藏的储层厚度、试油气资料，获取裂缝型砂岩储层工业价值气藏的储层厚度下限值为7米。

如图9所示，采集目的层段的孔隙度分析资料、地球物理资料，获取目的层段的孔隙度分布。

如图10所示，采集裂缝型砂岩储层已获工业价值气藏的孔隙度分析资料、试气资料、获取裂缝型砂岩储层工业价值气藏的孔隙度下限值为3.5％。

如图11所示，采集目的层段的地球物理资料，获取目的层段的盖层厚度分布。目的层段的盖层岩性为膏盐岩。

如图12所示，采集已获工业价值气藏的盖层厚度，获取工业价值气藏的膏盐岩盖层厚度下限值为8米。

采集目的层段的测井资料、烃源岩岩心分析的有机碳含量(TOC)，获取目的层段的有机碳含量分布。

如图13所示，采集目的层段的烃源岩岩心分析的镜质体反射率(Ro)，获取Ro分布。

如图14、图15所示，通过盆地模拟，获取生气丰度、生油丰度分布图。

采集目的层段的烃源岩岩心分析的有机碳含量(TOC)、岩心分析的镜质体反射率(Ro)资料，确定目的层段的天然气来源于侏罗系和三叠系。采集目的层段的测井资料、地震资料，通过地震解释确定烃源岩厚度，通过测井解释确定TOC分布。通过盆地模拟确定目的层段的生气丰度分布(图14)、生油丰度分布(图15)。根据对已获工业价值油气藏的含油气盆地解剖，当生气丰度大于10×10⁸m³/km²，生油丰度大于400×10⁴t/km²时，能够形成工业价值油气藏，从而确定有效生烃范围。

如图1所示，在本实施方式中，S5：根据控因参数分布获取每个网格坐标点的成藏参数值。

如图5所示，在一个具体的实施方式中，步骤S5：根据参考点的控因参数获取每个参考点的成藏参数值，具体包括：

S51：根据控因参数获取每个参考点的无量纲标志值；

S53：根据无量纲标志值获取每个参考点的成藏参数值。

在本实施方式中，根据以下公式获取每个参考点的成藏参数值：

其中，AR—成藏参数值，当AR≥1时为有效成藏区，当AR<1时为非有效成藏区，无量纲；Para_i—无量纲参数值Sth_S,Ts_S,Da_S,Rd_S,F_S；n—Para_i的个数；Sth_S—盖层厚度标志值，无量纲；Ts_S—输导体系标志值，无量纲；Da_S—地层倾角标志值，无量纲；Rd_S—油气运移距离标志值，无量纲；F_S—断层标志值。

在本实施方式中，根据以下公式确定盖层厚度标志值：

其中，Sth_S—盖层厚度标志值，无量纲；Sth—目的层段的盖层厚度，m；Sth_limt—工业价值油气藏的盖层厚度下限值，m。

在本实施方式中，根据以下公式确定输导体系标志值：

其中，Ts_S—输导体系标志值，无量纲；H—目的层段的储层厚度，m；H_limt—工业价值油气藏的储层厚度下限值，m；UF—与烃源岩有效生烃范围内直接接触的不整合和/或断层；当存在与烃源岩有效生烃范围内直接接触的不整合和/或断层时，为1；当与烃源岩有效生烃范围内直接接触的不整合或断层均不存在时，为0。

如图6所示，在本实施方式中，根据以下步骤确定烃源岩有效生烃范围：

S55：采集工业价值油气藏的生油丰度下限值、生气丰度下限值；

S57：根据生油丰度下限值和生气丰度下限值确定烃源岩有效生烃范围。

在本实施方式中，根据以下公式确定地层倾角标志值：

Da_S＝sin(α)，

其中，Da_S—地层倾角标志值，无量纲；α—地层倾角，度。

在本实施方式中，根据以下公式确定油气运移距离标志值：

其中，Rd_S—油气运移距离标志值，无量纲；Rd—有效生烃范围等效半径，km；L_hg—有效生烃范围外距有效生烃范围边界的距离，km，当处于有效生烃范围之内时，为0；a—经验系数，当目的层段为常规油气时，取3，当目的层段为致密油气时，取1.2，当目的层段为页岩油气时，取1。

在本实施方式中，当断穿目的层段的盖层的断层为开启断层时，断层标志值为-1；当断穿目的层段的盖层的断层为封闭断层时，断层标志值为1。

如图1所示，在本实施方式中，S7：根据网格坐标点的评价参数、评价参数下限值和成藏参数值获取每个网格坐标点的评价参数处理值。

在本实施方式中，根据以下公式确定每个网格坐标点的评价参数处理值：

其中，P_i—第i个网格坐标点的评价参数处理值；S_i—第i个评价参数值；S_{i_limt}—第i个评价参数下限值；n—评价参数的个数。

如图1所示，在本实施方式中，S9：根据网格坐标点的评价参数处理值获取目的层段的区带评价值。具体地，可以根据以下公式确定目的层段的区带评价值：

其中，V_Play—区带评价值；-1～1，当V_Play≥0时，为有效区带分布区，当V_Play<0时，为无效区带分布区；P_i—第i个网格坐标点的评价参数处理值；n—网格坐标点的个数。

在本实施方式中，根据以下公式确定f(P_i)：

其中，P_i—第i个网格坐标点的评价参数处理值。

图16为研究区目的层的区带评价结果。从已经钻探井的油气发现情况看：区带评价值V_Play≥0时的区域为有效区带分布区，其范围内的有效圈闭可以获得工业价值油气藏；区带评价值V_Play<0时的区域为无效区带分布区，其范围内的区域，所钻井均未获得工业价值油气藏。

如图17所示，本申请实施方式提供一种油气区带有效性评价装置，包括：控因参数获取模块11，用于获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；评价参数下限值获取模块13，用于获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；成藏参数值计算模块15，用于根据所述控因参数分布计算每个网格坐标点的成藏参数值；评价参数处理值计算模块17，用于根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值计算每个所述网格坐标点的评价参数处理值；区带评价值计算模块19，用于根据所述网格坐标点的评价参数处理值计算所述目的层段的区带评价值。

本申请实施方式所述的油气区带有效性评价装置根据网格坐标点的评价参数、评价参数下限值和成藏参数值计算每个网格坐标点的评价参数处理值，从而确定目的层段的成藏范围。然后根据网格坐标点的评价参数处理值计算目的层段的区带评价值，如此将油气区带有效性评价定量处理。从而相比于现有技术，本申请所述的油气区带有效性评价装置能制定统一的油气区带有效性评价标准。从而能排除人为因素的影响。又由于采用已获得工业价值油气藏的评价参数下限值，从而提高了油气区带有效性的识别率，从而提高了油气区带有效性评价的符合率，如此能提高钻探成功率及发现有效区带的速度，降低油气勘探成本。

在本实施方式中，控因参数获取模块11，用于获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度。

在一个具体的实施方式中，控因参数获取模块可以根据以下步骤获取目的层段的地层倾角分布：

1)根据目的层段的地球物理资料获取目的层段的构造；

2)根据目的层段的构造获取目的层段的地层倾角分布。

进一步地，根据以下步骤确定目的层段的地层倾角：

1)获取烃源岩有效生烃中心与所述目的层段的顶面构造等值线的垂直线段；

2)获取所述垂直线段的倾斜角；

3)根据所述垂直线段与水平线之间的位置关系获取所述倾斜角的倾斜方向；

具体地，当烃源岩有效生烃中心与目的层段的顶面构造等值线的垂直线，位于烃源岩有效生烃中心的水平线之上时，地层倾角为正；当烃源岩有效生烃中心与目的层段的顶面构造等值线的垂直线，位于烃源岩有效生烃中心的水平线之下时，地层倾角为负；当烃源岩有效生烃中心与目的层段的顶面构造等值线的垂直线，与位于烃源岩有效生烃中心的水平线重合时，地层倾角为0。

4)根据所述倾斜角和所述倾斜方向获取所述地层倾角。

在一个具体的实施方式中，断层包括：断穿目的层段的盖层的第一断层、沟通目的层段的烃源岩有效生烃范围内烃源岩与储层的第二断层和沟通目的层段的不整合与储层的第三断层。

在一个具体的实施方式中，与烃源岩有效生烃范围有效的不整合可以是烃源岩层与盖层之间的不整合。

在一个具体的实施方式中，控因参数获取模块11获取目的层段的储层厚度分布可以通过采集目的层段的地球物理资料，获取目的层段的储层厚度分布。

在一个具体的实施方式中，控因参数获取模块11获取目的层段的孔隙度分布可以通过采集目的层段的孔隙度分析资料和地球物理资料，获取目的层段的孔隙度分布。

在一个具体的实施方式中，控因参数获取模块11获取目的层段的盖层厚度分布可以通过采集目的层段的地球物理资料，获取目的层段的盖层厚度分布。

在一个具体的实施方式中，控因参数获取模块11获取目的层段的生油丰度分布和生气丰度分布，具体包括：

1)采集目的层段内的测井资料和烃源岩岩心分析的有机碳含量，获取目的层段的有机碳含量分布；

2)采集目的层段的烃源岩岩心分析的镜质体反射率，获取镜质体反射率分布；

3)对有机碳含量分布和镜质体反射率分布进行盆地模拟运算，获取生油丰度分布和生气丰度分布。

在一个具体的实施方式中，评价参数下限值获取模块13可以根据工业价值油气藏的储层厚度和试油气资料，获取储层厚度的下限值，其中，工业价值油气藏的储层类型与目的层段的储层类型相同。

在一个具体的实施方式中，评价参数下限值获取模块13可以根据工业价值油气藏的孔隙度分析资料和试气资料，获取孔隙度的下限值。

在一个具体的实施方式中，评价参数下限值获取模块13可以根据工业价值油气藏的盖层厚度，获取盖层厚度的下限值，其中，工业价值油气藏的盖层岩性与目的层段的盖层岩性一致。

在一个具体的实施方式中，工业价值油气藏的生油丰度下限值、生气丰度下限值用于确定烃源岩有效生烃范围。

利用本申请所述的油气区带有效性评价装置，对目的层段塔里木盆地库车山前地区的油气区带有效性进行评价。目的层段为三叠系、侏罗系烃源岩与白垩系巴什基奇克组碎屑岩储层组成的油气区带，目的层段圈闭以聚集天然气为主。

采集目的层段的地球物理资料，获取目的层段的构造、地层倾角分布、断层以及与烃源岩有效生烃范围有效的不整合。

如图7所示，采集目的层段的地球物理资料，获取目的层段的储层厚度分布。目的层段主要发育气藏为主，储层为裂缝型砂岩。

如图8所示，采集裂缝型砂岩储层已获得工业价值气藏的储层厚度、试油气资料，获取裂缝型砂岩储层工业价值气藏的储层厚度下限值为7米。

如图9所示，采集目的层段的孔隙度分析资料、地球物理资料，获取目的层段的孔隙度分布。

如图10所示，采集裂缝型砂岩储层已获工业价值气藏的孔隙度分析资料、试气资料、获取裂缝型砂岩储层工业价值气藏的孔隙度下限值为3.5％。

如图11所示，采集目的层段的地球物理资料，获取目的层段的盖层厚度分布。目的层段的盖层岩性为膏盐岩。

如图12所示，采集已获工业价值气藏的盖层厚度，获取工业价值气藏的膏盐岩盖层厚度下限值为8米。

采集目的层段的测井资料、烃源岩岩心分析的有机碳含量(TOC)，获取目的层段的有机碳含量分布。

如图13所示，采集目的层段的烃源岩岩心分析的镜质体反射率(Ro)，获取Ro分布。

如图14、图15所示，通过盆地模拟，获取生气丰度、生油丰度分布图。

采集目的层段的烃源岩岩心分析有机碳(TOC)、镜质体反射率(Ro)资料，确定目的层段的天然气来源于侏罗系和三叠系。采集目的层段的测井资料、地震资料，通过地震解释确定烃源岩厚度，通过测井解释确定TOC分布。通过盆地模拟确定目的层段的生气丰度分布(图14)、生油丰度分布(图15)。根据对已获工业价值油气藏的含油气盆地解剖，当生气丰度大于10×10⁸m³/km²，生油丰度大于400×10⁴t/km²时，能够形成工业价值油气藏，从而确定有效生烃范围。

成藏参数值计算模块15可以根据以下步骤计算每个参考点的成藏参数值：

1)根据控因参数获取每个参考点的无量纲标志值；

2)根据无量纲标志值获取每个参考点的成藏参数值。

在本实施方式中，根据以下公式获取每个参考点的成藏参数值：

其中，AR—成藏参数值，当AR≥1时为有效成藏区，当AR<1时为非有效成藏区，无量纲；Para_i—无量纲参数值Sth_S,Ts_S,Da_S,Rd_S,F_S；n—Para_i的个数；Sth_S—盖层厚度标志值，无量纲；Ts_S—输导体系标志值，无量纲；Da_S—地层倾角标志值，无量纲；Rd_S—油气运移距离标志值，无量纲；F_S—断层标志值。

在本实施方式中，根据以下公式确定盖层厚度标志值：

其中，Sth_S—盖层厚度标志值，无量纲；Sth—目的层段的盖层厚度，m；Sth_limt—工业价值油气藏的盖层厚度下限值，m。

在本实施方式中，根据以下公式确定输导体系标志值：

其中，Ts_S—输导体系标志值，无量纲；H—目的层段的储层厚度，m；H_limt—工业价值油气藏的储层厚度下限值，m；UF—与烃源岩有效生烃范围内直接接触的不整合和/或断层；当存在与烃源岩有效生烃范围内直接接触的不整合和/或断层时，为1；当与烃源岩有效生烃范围内直接接触的不整合或断层均不存在时，为0。

如图6所示，在本实施方式中，根据以下步骤确定烃源岩有效生烃范围：

1)采集工业价值油气藏的生油丰度下限值、生气丰度下限值；

2)根据生油丰度下限值和生气丰度下限值确定烃源岩有效生烃范围。

在本实施方式中，根据以下公式确定地层倾角标志值：

Da_S＝sin(α)，

其中，Da_S—地层倾角标志值，无量纲；α—地层倾角，度。

在本实施方式中，根据以下公式确定油气运移距离标志值：

其中，Rd_S—油气运移距离标志值，无量纲；Rd—有效生烃范围等效半径，km；L_hg—有效生烃范围外距有效生烃范围边界的距离，km，当处于有效生烃范围之内时，为0；a—经验系数，当目的层段为常规油气时，取3，当目的层段为致密油气时，取1.2，当目的层段为页岩油气时，取1。

在本实施方式中，当断穿目的层段的盖层的断层为开启断层时，断层标志值为-1；当断穿目的层段的盖层的断层为封闭断层时，断层标志值为1。

如图1所示，在本实施方式中，评价参数处理值计算模块17可以根据以下公式确定每个网格坐标点的评价参数处理值：

其中，P_i—第i个网格坐标点的评价参数处理值；S_i—第i个评价参数值；S_{i_limt}—第i个评价参数下限值；n—评价参数的个数。

在本实施方式中，区带评价值计算模块19可以根据以下公式确定目的层段的区带评价值：

其中，V_Play—区带评价值；-1～1，当V_Play≥0时，为有效区带分布区，当V_Play<0时，为无效区带分布区；P_i—第i个网格坐标点的评价参数处理值；n—网格坐标点的个数。

在本实施方式中，根据以下公式确定f(P_i)：

其中，P_i—第i个网格坐标点的评价参数处理值。

图16为研究区目的层的区带评价结果。从已经钻探井的油气发现情况看：区带评价值V_Play≥0时的区域为有效区带分布区，其范围内的有效圈闭可以获得工业价值油气藏；区带评价值V_Play<0时的区域为无效区带分布区，其范围内的区域，所钻井均未获得工业价值油气藏。

以上仅为本发明的几个实施例，本领域技术人员依据申请文件公开的内容可以对本发明实施例进行各种改动或变型而不脱离本发明的精神和范围。

Claims (11)

1.一种油气区带有效性评价方法，其特征在于，包括：

获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；

获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；

根据所述控因参数分布计算每个网格坐标点的成藏参数值；

根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值计算每个所述网格坐标点的评价参数处理值；

根据所述网格坐标点的评价参数处理值计算所述目的层段的区带评价值；

所述根据所述控因参数分布计算每个网格坐标点的成藏参数值，具体包括：

根据所述控因参数获取每个所述网格坐标点的无量纲标志值；

根据每个所述网格坐标点的所述无量纲标志值获取每个所述网格坐标点的成藏参数值；

根据以下公式获取每个所述网格坐标点的成藏参数值：

其中，AR—成藏参数值，当AR≥1时为有效成藏区，当AR<1时为非有效成藏区，无量纲；Para_i—无量纲参数值Sth_S,Ts_S,Da_S,Rd_S,F_S；n—Para_i的个数；Sth_S—盖层厚度标志值，无量纲；Ts_S—输导体系标志值，无量纲；Da_S—地层倾角标志值，无量纲；Rd_S—油气运移距离标志值，无量纲；F_S—断层标志值；

根据以下公式确定所述盖层厚度标志值：

其中，Sth_S—盖层厚度标志值，无量纲；Sth—目的层段的盖层厚度，m；Sth_limt—工业价值油气藏的盖层厚度下限值，m；

根据以下公式确定所述输导体系标志值：

其中，Ts_S—输导体系标志值，无量纲；H—目的层段的储层厚度，m；H_limt—工业价值油气藏的储层厚度下限值，m；UF—与烃源岩有效生烃范围内直接接触的不整合和/或断层；当存在与烃源岩有效生烃范围内直接接触的不整合和/或断层时，为1；当与烃源岩有效生烃范围内直接接触的不整合或断层均不存在时，为0；

根据以下公式确定所述地层倾角标志值：

Da_S＝sin(α)，

其中，Da_S—地层倾角标志值，无量纲；α—地层倾角，度；

根据以下公式确定所述油气运移距离标志值：

其中，Rd_S—油气运移距离标志值，无量纲；Rd—有效生烃范围等效半径，km；L_hg—有效生烃范围外距有效生烃范围边界的距离，km，当处于有效生烃范围之内时，为0；a—经验系数，当目的层段为常规油气时，取3，当目的层段为致密油气时，取1.2，当目的层段为页岩油气时，取1；

当断穿目的层段的盖层的断层为开启断层时，所述断层标志值为-1；当断穿目的层段的盖层的断层为封闭断层时，所述断层标志值为1；

根据以下公式计算每个所述网格坐标点的评价参数处理值：

其中，P_i—第i个所述网格坐标点的评价参数处理值；S_i—第i个所述评价参数值；S_{i_limt}—第i个所述评价参数下限值；n—所述评价参数的个数；

根据以下公式计算所述目的层段的区带评价值：

其中，V_Play—区带评价值；-1～1，当V_Play≥0时，为有效区带分布区，当V_Play<0时，为无效区带分布区；P_i—第i个所述网格坐标点的评价参数处理值；n—所述网格坐标点的个数。

2.根据权利要求1所述的油气区带有效性评价方法，其特征在于：根据以下步骤确定所述烃源岩有效生烃范围：

获取工业价值油气藏的生油丰度下限值、生气丰度下限值；

根据所述生油丰度下限值和所述生气丰度下限值确定所述烃源岩有效生烃范围。

3.根据权利要求1所述的油气区带有效性评价方法，其特征在于：根据以下公式确定所述f(P_i)：

其中，P_i—第i个所述网格坐标点的评价参数处理值。

4.根据权利要求1所述的油气区带有效性评价方法，其特征在于：所述获取所述目的层段的地层倾角分布，具体包括：

根据目的层段的地球物理资料获取所述目的层段的构造；

根据所述目的层段的构造获取所述目的层段的地层倾角分布。

5.根据权利要求4所述的油气区带有效性评价方法，其特征在于：根据以下步骤确定所述目的层段的地层倾角：

获取烃源岩有效生烃中心与所述目的层段的顶面构造等值线的垂直线段；

获取所述垂直线段的倾斜角；

根据所述垂直线段与水平线之间的位置关系获取所述倾斜角的倾斜方向；

根据所述倾斜角和所述倾斜方向获取所述地层倾角。

6.根据权利要求1所述的油气区带有效性评价方法，其特征在于：所述断层包括：断穿所述目的层段的盖层的第一断层、沟通所述目的层段的烃源岩有效生烃范围内烃源岩与储层的第二断层和沟通所述目的层段的不整合与储层的第三断层。

7.根据权利要求1所述的油气区带有效性评价方法，其特征在于，还包括：

根据所述目的层段的测井资料和烃源岩岩心分析的有机碳含量，获取所述目的层段的所述有机碳含量分布；

根据所述目的层段的烃源岩岩心分析的镜质体反射率，获取所述镜质体反射率分布；

对所述有机碳含量分布和所述镜质体反射率分布进行盆地模拟运算，获取所目的层段的生油丰度分布和生气丰度分布。

8.根据权利要求1所述的油气区带有效性评价方法，其特征在于：所述获取工业价值油气藏的储层厚度的下限值，具体包括：

根据工业价值油气藏的储层厚度和试油气资料，获取储层厚度的下限值，其中，所述工业价值油气藏的储层类型与所述目的层段的储层类型相同。

9.根据权利要求1所述的油气区带有效性评价方法，其特征在于：所述获取工业价值油气藏的孔隙度的下限值，具体包括：

根据工业价值油气藏的孔隙度分析资料和试气资料，获取孔隙度的下限值。

10.根据权利要求1所述的油气区带有效性评价方法，其特征在于：所述获取工业价值油气藏的盖层厚度的下限值，具体包括：

根据工业价值油气藏的盖层厚度，获取盖层厚度的下限值，其中，所述工业价值油气藏的盖层岩性与所述目的层段的盖层岩性一致。

11.一种油气区带有效性评价装置，其特征在于，包括：

控因参数获取模块，用于获取目的层段的控因参数分布，其中所述控因参数包括：地层倾角、断层、与烃源岩有效生烃范围有效的不整合、储层厚度、盖层厚度、生油丰度、生气丰度；

评价参数下限值获取模块，用于获取工业价值油气藏的评价参数下限值，其中，所述评价参数包括：储层厚度、孔隙度、盖层厚度、生油丰度、生气丰度；

成藏参数值计算模块，用于根据所述控因参数分布计算每个网格坐标点的成藏参数值；

评价参数处理值计算模块，用于根据所述网格坐标点的所述评价参数、所述评价参数下限值和所述成藏参数值计算每个所述网格坐标点的评价参数处理值；

区带评价值计算模块，用于根据所述网格坐标点的评价参数处理值计算所述目的层段的区带评价值；

所述成藏参数值计算模块用于根据所述控因参数获取每个所述网格坐标点的无量纲标志值；并根据每个所述网格坐标点的所述无量纲标志值获取每个所述网格坐标点的成藏参数值；

根据以下公式获取每个所述网格坐标点的成藏参数值：

其中，AR—成藏参数值，当AR≥1时为有效成藏区，当AR<1时为非有效成藏区，无量纲；Para_i—无量纲参数值Sth_S,Ts_S,Da_S,Rd_S,F_S；n—Para_i的个数；Sth_S—盖层厚度标志值，无量纲；Ts_S—输导体系标志值，无量纲；Da_S—地层倾角标志值，无量纲；Rd_S—油气运移距离标志值，无量纲；F_S—断层标志值；

根据以下公式确定所述盖层厚度标志值：

其中，Sth_S—盖层厚度标志值，无量纲；Sth—目的层段的盖层厚度，m；Sth_limt—工业价值油气藏的盖层厚度下限值，m；

根据以下公式确定所述输导体系标志值：

其中，Ts_S—输导体系标志值，无量纲；H—目的层段的储层厚度，m；H_limt—工业价值油气藏的储层厚度下限值，m；UF—与烃源岩有效生烃范围内直接接触的不整合和/或断层；当存在与烃源岩有效生烃范围内直接接触的不整合和/或断层时，为1；当与烃源岩有效生烃范围内直接接触的不整合或断层均不存在时，为0；

根据以下公式确定所述地层倾角标志值：

Da_S＝sin(α)，

其中，Da_S—地层倾角标志值，无量纲；α—地层倾角，度；

根据以下公式确定所述油气运移距离标志值：

其中，Rd_S—油气运移距离标志值，无量纲；Rd—有效生烃范围等效半径，km；L_hg—有效生烃范围外距有效生烃范围边界的距离，km，当处于有效生烃范围之内时，为0；a—经验系数，当目的层段为常规油气时，取3，当目的层段为致密油气时，取1.2，当目的层段为页岩油气时，取1；

当断穿目的层段的盖层的断层为开启断层时，所述断层标志值为-1；当断穿目的层段的盖层的断层为封闭断层时，所述断层标志值为1；

所述评价参数处理值计算模块用于根据以下公式计算每个所述网格坐标点的评价参数处理值：

其中，P_i—第i个所述网格坐标点的评价参数处理值；S_i—第i个所述评价参数值；S_{i_limt}—第i个所述评价参数下限值；n—所述评价参数的个数；

所述区带评价值计算模块用于根据以下公式计算所述目的层段的区带评价值：

其中，V_Play—区带评价值；-1～1，当V_Play≥0时，为有效区带分布区，当V_Play<0时，为无效区带分布区；P_i—第i个所述网格坐标点的评价参数处理值；n—所述网格坐标点的个数。