CN104933282A - 多层圈闭含油气概率的获取方法和系统 - Google Patents
多层圈闭含油气概率的获取方法和系统 Download PDFInfo
- Publication number
- CN104933282A CN104933282A CN201410109198.7A CN201410109198A CN104933282A CN 104933282 A CN104933282 A CN 104933282A CN 201410109198 A CN201410109198 A CN 201410109198A CN 104933282 A CN104933282 A CN 104933282A
- Authority
- CN
- China
- Prior art keywords
- trap
- probability
- layer
- oily
- multilayer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
本发明公开了一种多层圈闭含油气概率的获取方法和系统,该方法包括:根据多层圈闭中各层圈闭的每个成藏条件之间的关系来确定各层圈闭的共性成藏条件和独立成藏条件,并获取各层圈闭的每个成藏条件的概率值;基于各个成藏条件的概率值计算多层圈闭中各层圈闭的边际概率和条件概率,进而得到各层圈闭的含油气概率;根据各层圈闭的边际概率值范围确定各层圈闭之间的决定关系,进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,作为多层圈闭的含油气概率。本发明方法在计算多层圈闭的含油气概率时,充分考虑各层圈闭的成藏要素之间的相互依赖关系和各层圈闭含油气概率的各种组合关系,使得计算出的多层圈闭含油气概率更为合理和准确。
Description
技术领域
本发明涉及油气勘探领域,尤其涉及一种多层圈闭含油气概率的获取方法和系统。
背景技术
早期的油气资源分析一般会忽视地质风险分析,这样会导致油气资源的分析不准确。因此,在后期引进了统计学与概率理论的分析方法,采用地质风险概率法,根据圈闭的各项成藏条件发生的概率对圈闭含油气概率进行综合分析。
对于地质风险概率法(详见《油气勘探风险分析与实物期权法经济评价》.贾承造,杨树锋,张永锋,等.2004.北京:石油工业出版社)而言,由于地下地质的复杂性与不确定性,油气资源分析所反映出的把握性是不确定的。可以采用不同概率条件下的把握性大小分析来科学地反映这种不确定性。
根据含油气系统理论,把烃源条件、储层条件、圈闭条件及保存条件等成藏要素看成是相互独立的,如果圈闭含油气,则上述几项条件缺一不可。因此,圈闭含油气的概率可以表示为:
其中,P表示含油气概率(0≤P≤1),Pi表示单项地质条件发生的概率(0≤Pi≤1)。
地质风险概率法可用于计算单个油气聚集单元(次级圈闭)的含油气概率,但是无法综合评价含有多个油气聚集单元(次级圈闭)的评价目标的含油气概率。
普通概率树模型(详见《概率树在全概率公式中的应用》.李晓红.高等数学研究.2008,Jul.Vol.11(4))为常见的组合概率计算技术手段之一,其具有典型的二叉树数据结构。同一节点的两个分支可被赋予当前油气聚集单元(次级圈闭)含有和不含有油气的概率值,N个油气聚集单元(次级圈闭)就等同进行N次含油气性试验,会形成2n个叶节点,每一叶节点代表不同油气聚集单元(次级圈闭)含油气状态的组合,该组合的概率值就等于从根节点到此叶节点路径上,所有分支概率值的乘积。
但是,普通概率树假定N次试验是相互独立,在不同油气聚集单元(次级圈闭)衔接时,并不能够很好地反映不同油气聚集单元(次级圈闭)间可能存在的各种类型的含油气依赖关系,也不能够确定含油气聚合单元的可能含油气组合。因此,亟需一种解决方案,以能够确定多层圈闭的含油气概率。
发明内容
本发明所要解决的技术问题之一是需要提供一种多层圈闭含油气概率的获取方法,该方法能够根据各层圈闭之间的依赖关系,获取不同含油气组合的多层圈闭的含油气概率。另外,还提供了一种多层圈闭含油气概率的获取系统。
为了解决上述技术问题,本发明提供了一种多层圈闭含油气概率的获取方法,包括:步骤一,根据多层圈闭中各层圈闭的每个成藏条件之间的关系来确定各层圈闭的共性成藏条件和独立成藏条件,并获取各层圈闭的每个成藏条件的概率值;步骤二,基于所述各个成藏条件的概率值计算多层圈闭中各层圈闭的边际概率和条件概率,进而得到各层圈闭的含油气概率;步骤三,根据各层圈闭的边际概率值范围确定各层圈闭之间的决定关系,进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,作为多层圈闭的含油气概率。
在一个实施例中,所述步骤三进一步包括:采用融合有多层圈闭的边际概率和各层圈闭的条件概率的概率树,以获取各种层圈闭组合下的多层圈闭的含油气概率进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,其中,层圈闭组合基于决定关系和各层圈闭的含油气概率而得到。
在一个实施例中,所述决定关系包括:完全决定关系、部分决定关系以及完全独立关系,在完全决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率等于各层圈闭中含油气概率的最大值,则对应含油气概率最大值的那层圈闭的含油气有无直接决定其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为完全相关的各种层圈闭组合下的含油气概率的加和;在部分决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率大于各层圈闭中含油气概率的最大值并小于1,则对应含油气概率最大值的那层圈闭的含油气的有无间接反映其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为部分相关的各种层圈闭组合下的含油气概率的加和;在完全独立关系中,各层圈闭间不存在共性成藏条件,该多层圈闭的边际概率为1,每一层圈闭的含油气有无仅取决于自身的成藏条件,所述多层圈闭的含油气概率为完全独立的各种层圈闭组合下的含油气概率的加和。
在一个实施例中,在所述概率树中,以多层圈闭的边际概率为对象,作为根节点,各层圈闭分别作为其他不同级别节点,且各级别节点的状态对应各层圈闭的条件概率。
在一个实施例中,通过如下表达式计算多层圈闭中至少有一层发现油气的圈闭组合概率P含油气概率:
其中,P边际概率表示多层圈闭的边际概率,P第i个条件概率表示第i层圈闭的条件概率,n表示圈闭的层数。
在一个实施例中,在步骤二中,采用贝叶斯概率模型计算多层圈闭中各层圈闭的边际概率和条件概率;基于各层圈闭的边际概率和条件概率得到各层圈闭的含油气概率,其中,所述边际概率为各个共性成藏条件的概率值乘积,所述条件概率为各个独立成藏条件的概率值乘积,所述各层圈闭的含油气概率为非条件概率,其为所述条件概率和所述边际概率的乘积。
根据本发明的另一方面,还提供了一种多层圈闭含油气概率的获取系统,包括:第一单元,其根据多层圈闭中各层圈闭的每个成藏条件之间的关系来确定各层圈闭的共性成藏条件和独立成藏条件,并获取各层圈闭的每个成藏条件的概率值;第二单元,其基于所述各个成藏条件的概率值计算多层圈闭中各层圈闭的边际概率和条件概率,进而得到各层圈闭的含油气概率;第三单元,其根据各层圈闭的边际概率值范围确定各层圈闭之间的决定关系,进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,作为多层圈闭的含油气概率。
在一个实施例中,所述第三单元进一步用于:采用融合有多层圈闭的边际概率和各层圈闭的条件概率的概率树,以获取各种层圈闭组合下的多层圈闭的含油气概率进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,其中,层圈闭组合基于决定关系和各层圈闭的含油气概率而得到。
在一个实施例中,所述决定关系包括:完全决定关系、部分决定关系以及完全独立关系,在完全决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率等于各层圈闭中含油气概率的最大值,则对应含油气概率最大值的那层圈闭的含油气有无直接决定其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为完全相关的各种层圈闭组合下的含油气概率的加和;在部分决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率大于各层圈闭中含油气概率的最大值并小于1,则对应含油气概率最大值的那层圈闭的含油气的有无间接反映其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为部分相关的各种层圈闭组合下的含油气概率的加和;在完全独立关系中,各层圈闭间不存在共性成藏条件,该多层圈闭的边际概率为1,每一层圈闭的含油气有无仅取决于自身的成藏条件,所述多层圈闭的含油气概率为完全独立的各种层圈闭组合下的含油气概率的加和。
在一个实施例中,在所述概率树中,以多层圈闭的边际概率为对象,作为根节点,各层圈闭分别作为其他不同级别节点,且各级别节点的状态对应各层圈闭的条件概率。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:
针对现有技术的不足,本发明的方法在计算含有多个次级油气聚集单元的圈闭(勘探目标)的含油气概率时,充分考虑各次级油气聚集单元的成藏要素之间的相互依赖关系和各油气聚集单元含油气概率的各种组合关系,使得计算出的多层圈闭含油气概率更为合理和准确。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1是根据本发明一实施例的多层圈闭含油气概率的获取方法的流程图;
图2是根据本发明一示例的多层圈闭含油气概率的获取方法的流程图;
图3是根据本发明一示例的多层圈闭(含有A、B两个层圈闭)的含油气概率的概率树模型的示意图;
图4是根据本发明另一示例的多层圈闭(具有四次级圈闭)的含油气概率组合的概率树模型的示意图;
图5是根据本发明另一实施例的多层圈闭含油气概率的获取系统5的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下结合附图对本发明作进一步地详细说明。
需要说明的是,多层圈闭是指在一个勘探目标(二级、三级或四级圈闭)中,包含多个单一的层圈闭(最小油气聚集单元)。多层圈闭的含油气概率就是指多个层圈闭中至少有一个层圈闭勘探成功,即有油气发现的概率。
第一实施例
图1是根据本发明一实施例的多层圈闭含油气概率的获取方法的流程图。下面结合图1对本方法进行说明。
步骤S110,根据多层圈闭中各层圈闭的每个成藏条件之间的关系来确定各层圈闭的共性(全局)风险因子和独立(局部)风险因子,并获取各层圈闭的每个成藏条件的概率值。
一般,成藏条件包括圈闭条件、储层条件、充注条件以及保存条件。具体地,圈闭条件主要包括侧网密度、剖面质量、解释可信度、圈闭类型、圈闭面积和幅度以及圈闭闭合形式。储层条件包括岩性、沉积相、储层厚度、储层物性以及成岩作用等。充注条件包括烃源岩厚度、有机质丰度、有机质类型、热演化程度、运移通道、运移距离以及流线形式等。保存条件包括盖层条件、断层、构造运动以及水动力等。而所谓的共性(全局)风险因子和独立(局部)风险因子,也可以分别被称为共性(全局)成藏条件和独立(局部)成藏条件。
具体地,通过分析目标圈闭的构造可以判定出各个成藏条件的关系,根据这种关系,进一步地将成藏条件归类为共性成藏条件和独立成藏条件。可以理解,共性成藏条件为各层圈闭(各个油气聚集单元)的某一成藏条件是相同或相似的,而独立成藏条件为各个油气聚集单元所拥有的不同的成藏条件。例如,若某一目标圈闭的目的层均来自于相同的沉积环境,且各层圈闭的油气充注也是以断层作为主要的运移通道,则这两个因素为共性成藏条件,即为共性风险因子。另外两个成藏条件,即圈闭条件和保存条件为独立成藏条件,即为独立风险因子。
然后,通过建立评价单元,选择评价模板来获取每个油气聚集单元的各个成藏条件的概率值(也称为评价值)。此步骤为现有技术,在此不再赘述。
步骤S120,基于上述各个成藏条件的概率值计算多层圈闭中各层圈闭的边际概率、条件概率,进而得到各层圈闭的含油气概率(非条件概率)。
具体地,采用贝叶斯概率模型来计算每层圈闭的边际概率、条件概率以及非条件概率。其中,边际概率为单层圈闭共性风险因子的概率值乘积(即相同或相似成藏条件的成功概率乘积),条件概率为单层圈闭独立风险因子的概率值乘积(即独立成藏条件的成功概率乘积),非条件概率为条件概率和边际概率的乘积,亦为单层圈闭的含油气概率。
举例而言,具体单层圈闭地质因子风险概率设置表格可如下表1所示:
表1
步骤S130,根据各层圈闭的边际概率值范围确定各层圈闭之间的决定关系(关系计算模型),进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,作为多层圈闭的含油气概率。
在一个优选示例中,可以采用融合有多层圈闭的边际概率和各层圈闭的条件概率的概率树,以获取各种层圈闭组合下的多层圈闭的含油气概率进而得到不同圈闭间影响下的多层圈闭的含油气概率,其中,层圈闭组合基于关系计算模型和各层圈闭的含油气概率得到的。
需要说明的是,多层圈闭含油气概率的计算与各层圈闭之间的相互关系密切相关,而现有技术中没有考虑这一因素。因此利用现有技术计算得到的圈闭含油气概率不太准确合理。
在本步骤中,关系计算模型包括:完全决定、部分决定以及完全独立这三种关系模型:
完全决定关系模型:各层圈闭的成藏条件具有完全依赖关系,含油气概率较大的层圈闭决定含油气概率较小的层圈闭。也可以说,各层圈闭间存在至少一个共性成藏条件,共性成藏条件相乘所得的概率(边际概率)等于各层圈闭中含油气概率最大值,意味着某一层圈闭(含油气概率最大者)的含油气有无直接决定了其它层圈闭的含油气性有无,多层圈闭的含油气概率为完全相关的各种层圈闭组合下的含油气概率的加和。
部分决定关系模型:各层圈闭的成藏条件存在部分依赖关系,含油气概率较大的层圈闭决定含油气概率较小的层圈闭。也可以说,各层圈闭间存在至少一个共性成藏条件,共性成藏条件相乘所得的概率(边际概率)大于各层圈闭中含油气概率最大值,且小于1.0,意味着某一层圈闭(含油气概率最大者)的含油气有无间接反映了其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为部分相关的各种层圈闭组合下的含油气概率的加和。
完全独立关系模型:各层圈闭的成藏条件彼此完全独立,即圈闭条件、充注条件、储层条件和保存条件等成藏要素完全独立,任何层圈闭的任意成藏要素与其他层圈闭相应的成藏要素没有依赖关系。也就是说,各层圈闭间不存在共性成藏条件,也就是意味着边际概率值1.0,意味着每一层圈闭的含油气有无仅仅取决于自生成藏条件,互相之间没有任何关系,所述多层圈闭的含油气概率为完全独立的各种层圈闭组合下的含油气概率的加和。
需要说明的是,本实施例提出的概率树改进方法,很好地融合了不同层圈闭的边际概率和条件概率,进而能够完整地进行“完全独立”、“部分决定”和“完全决定”三种类型的组合概率计算。
该方法将多层圈闭的“边际概率”作为“层圈闭”对象,作为根节点。其它不同级别节点对应不同层圈闭,各级别节点的状态对应各层圈闭的条件概率值,这就意味着:①当“边际概率”节点的值等于1时,不同叶节点计算结果等同普通概率树的组合概率计算结果,体现的是“完全独立”情况下的含油气概率;②当“边际概率”节点的概率值在所有层圈闭概率最大值到1.0之间时,不同叶节点计算结果等同“部分决定”情况下的组合概率;③当“边际概率”的值等于所有层圈闭中含油气概率最大值,且将含油气概率最大节点的两个分支调整为1.0和0.0,此时不同叶节点计算结果就表示在某一层圈闭“完全决定”后续层圈闭情况下的组合概率。
下面结合图2以及表2-表6对多层圈闭含油气概率的计算的过程稍作说明。
如图2所示,首先对有关目标圈闭(多层圈闭)的数据进行整理,建立评价单元和选择评价模板。然后针对不同的层圈闭设置共性因子。输入风险参数概率进而对每一层圈闭的含油气性进行评价。再后,设置风险相关性,即各层圈闭之间的关系。最终,基于改进后的概率树计算得到不同情况下的多层含油气概率值,并输出计算结果。
在表2到表4中,各个成藏条件的概率值是根据现有的评价模板得到的。
如表2所示,其表示“部分决定”地质背景下单层圈闭共性因子设置及含油气概率计算,其中包括边际概率、条件概率和非条件概率。
表2
从以上表2可以看出,该层圈闭的共性因子为圈闭条件和充注条件,独立因子为储层条件和保存条件。根据上面的计算方法,可知该层圈闭的边际概率为共性因子概率的乘积,即0.8×0.7=0.56,条件概率为独立因子概率的乘积,即0.8×0.5=0.4,非条件概率为边际概率和条件概率的乘积,为0.224。
表3表示“完全独立”地质背景下单层圈闭共性因子设置及含油气概率计算,其中包括边际概率、条件概率和非条件概率。
表3
该层圈闭的边际概率、条件概率和非条件概率的计算方法与表2的计算方法相同,在此不再赘述。
表4表示“完全决定”地质背景下单层圈闭共性因子设置及含油气概率计算,其中包括边际概率、条件概率和非条件概率。
表4
有关表4的层圈闭的边际概率、条件概率和非条件概率的计算方法亦与表2的计算方法相同,在此不再赘述。
在工程实践中,通过对地质条件分析可以判定各个成藏条件之间的关系,进而通过各个成藏条件的概率值可以获得三种关系模型中单层圈闭的含油气概率。
图3是多层圈闭(含有A、B两个层圈闭)的含油气概率的概率树模型的示意图。其中实心圆和T代表成功,空心圆和F代表失败,数字代表发生的概率。
结合图3而言,对应“多层圈闭”的节点为整个圈闭的边际概率(为共性因子概率的乘积)为0.48,且为圆圈(表示失败)。A层圈闭所对应的节点状态中的条件概率0.56,B层圈闭所对应的节点状态中的条件概率为0.72。而样本空间包括四种:TT、TF、FT和FF。通过将每条分支上所有节点上的数值进行乘积运算得到不同样本空间发生的可能概率。
如表5表示获取的为至少一个次级圈闭含油气的概率,表6表示获取的为至少两个次级圈闭含油气的概率。
表5
表6
下面结合图3,以结点状态TT为示例来说明如何计算两个次级圈闭均包含油气的条件概率和非条件概率。PA,B(条件概率)=0.56×0.72=0.4032;两个次级圈闭均包含油气的非条件概率(即含油气概率)为PA,B(非条件概率)=P(边际概率)×PA,B(条件概率)=0.1935。
而,至少一个次级圈闭含油气的次级圈闭组合对应的节点状态包括:TT、TF以及FT,因此,至少一个次级圈闭含油气的概率为这三种组合的概率值加和。至少一个次级圈闭含油气的概率含油气概率为:0.1935+0.0753+0.1521=0.4209(如表5所示)。至少两个次级圈闭含油气的概率组合只有一种情况,即次级圈闭A和B均含有油气,从表6中可看出该概率为:0.1935。
在另一个优选示例中,对不同层圈闭间含油气影响的量化表达,即考虑了不同层圈闭间含油气影响下的多层圈闭成功概率(多层圈闭中至少一层次级圈闭含油气的圈闭组合概率)的计算方法,可利用如下公式计算:
其中,P含油气表示多层圈闭含油气概率,P群风险表示由于共性成藏条件好坏而造成的所有层圈闭都不含有油气的概率,P所有个体风险表示由于局部成藏条件好坏而造成的单个层圈闭不含有油气的概率,P边际概率表示多层圈闭(整个圈闭)的边际概率,P第i个条件概率表示第i层圈闭的条件概率,n表示圈闭的层数。
具体地,在各层圈闭之间完全独立时,各层圈闭的地质风险因子成功概率不受其它层圈闭地质风险因子概率的影响,此时多层圈闭的边际概率为1,多层圈闭的含油气概率为完全独立的各个圈闭组合下的多层圈闭含油气概率的加和。
在多层圈闭之间部分决定时,各层圈闭的部分地质风险因子成功概率受其它层圈闭地质风险因子概率的影响,此时多层圈闭含油气边际概率为全局地质因子(共性成藏条件)的概率值的乘积,多层圈闭的含油气概率为部分相关的层圈闭组合下的含油气概率的加和。
当某一层圈闭是否有油气发现的情况完全决定其它层圈闭是否有油气发现时,该层圈闭的所有地质风险因子成功概率都影响到其它层圈闭地质风险因子概率,此时所有地质风险因子均为全局因子,多层圈闭含油气的边际概率为所有风险因子概率的乘积,多层圈闭的含油气概率为完全相关的层圈闭组合下多层圈闭含油气概率的加和。
示例
为了进一步理解本发明,下面对一示例进行说明。
以某断块圈闭为例,该断块包含四个次级圈闭(即油气聚集单元)。从地质条件分析,该断块的烃源来自构造西北低洼部位的烃源岩,目的层位为碎屑岩,盖层为上覆厚层泥岩。该圈闭目的层来自于相同的沉积环境,圈闭的油气充注也是以断层作为主要的运移通道,此二因素为全局成藏条件,为共性条件;圈闭条件和保存条件为局部成藏条件,为个性条件。该圈闭中各次级圈闭成藏要素评价结果见表7。
表7
根据以上地质认识可对T1断块各次级圈闭和整个圈闭群的含油气概率进行综合分析。
图4为四层圈闭的油气概率组合的概率树模型的示意图。图4中最左边的根节点数字表示整个圈闭的边际概率,第一级分支上的两个节点表示第一个次级圈闭(这里指Ⅰ号圈闭)的勘探预测结果(实心圆为成功,圆圈为失败,数字为成功或失败的条件概率);第二级分支上的节点表示在Ⅰ号圈闭预测结果的基础上第二个次级圈闭(Ⅱ号圈闭)的勘探预测结果;同理,第三级分支上的节点表示在Ⅰ、Ⅱ号圈闭预测结果的基础上第三个次级圈闭(Ⅲ号圈闭)的勘探预测结果;第四级分支的节点表示在Ⅰ、Ⅱ、Ⅲ号圈闭预测结果的基础上第四个次级圈闭(Ⅳ号圈闭)的勘探预测结果;样本空间表示该分支下各次级圈闭的预测结果(T为有油气发现,F为无油气发现)。在图4中,每一分支最后的概率表示该分支发现油气的可能概率,即非条件概率,非条件概率数值为该分支所有节点上数值的乘积。
具体的计算方法和计算过程如下:
整个圈闭含油气的边际概率为:
P(边际概率)=P(储层)×P(充注)=0.8×0.85=0.68;
各次级圈闭含油气的条件概率为:
PⅠ(条件概率)=P(圈闭)×P(保存)=0.8×0.9=0.72;
PⅡ(条件概率)=P(圈闭)×P(保存)=0.8×0.85=0.68;
PⅢ(条件概率)=P(圈闭)×P(保存)=0.75×0.9=0.675;
PⅣ(条件概率)=P(圈闭)×P(保存)=0.75×0.8=0.6;
各次级圈闭含油气概率(即非条件概率):
PⅠ(含油气)=P(边际概率)×PⅠ(条件概率)=0.68×0.72=0.489;
PⅡ(含油气)=P(边际概率)×PⅡ(条件概率)=0.68×0.68=0.462;
PⅢ(含油气)=P(边际概率)×PⅢ(条件概率)=0.68×0.675=0.459;
PⅣ(含油气)=P(边际概率)×PⅣ(条件概率)=0.68×0.6=0.408;
如前所述,该断块圈闭各次级圈闭中,储层条件和充注条件为全局成藏条件(共性条件),同时各次级圈闭又有局部成藏条件(个性条件),因此各次级圈闭之间具有部分依赖性,即含油气概率较高的次级圈闭(Ⅰ号圈闭)部分决定含油气概率较低的次级圈闭(Ⅱ、Ⅲ、Ⅳ号圈闭),这种情况下,各次级圈闭发现油气的可能组合情况及各自发生的概率如下:
(1)4个次级圈闭均获得勘探成功(即均发现工业油气流),其概率为:
P4=(PⅠ(条件概率)×PⅡ(条件概率)×PⅢ(条件概率)×PⅣ(条件概率))×P(边际概率)=(0.72×0.68×0.675×0.60)×0.68≈0.135。
(2)只有3个次级圈闭获得勘探成功,包括4种情况(Ⅰ+Ⅱ+Ⅲ、Ⅰ+Ⅱ+Ⅳ、Ⅰ+Ⅲ+Ⅳ和Ⅱ+Ⅲ+Ⅳ),所以,有3个次级圈闭获得勘探成功的概率为这4种组合的概率加和,即:
P3=P123+P124+P134+P234=0.052+0.090+0.065+0.064≈0.27。
(3)只有2个次级圈闭获得勘探成功,包括6种情况(Ⅰ+Ⅱ、Ⅰ+Ⅲ、Ⅰ+Ⅳ、Ⅱ+Ⅲ、Ⅱ+Ⅳ和Ⅲ+Ⅳ),因此,有2个次级圈闭获得勘探成功的概率为这6种组合的概率加和,即:P2=P12+P13+P14+P23+P24+P34=0.043+0.042+0.031+0.035+0.025+0.025=0.201。
(4)只有1个次级圈闭获得勘探成功的概率包括这4种组合(即Ⅰ、Ⅱ、Ⅲ和Ⅳ分别单独成功)的概率加和,即:P1=P1-Ⅰ+P1-Ⅱ+P1-Ⅲ+P1-Ⅳ=0.021+0.017+0.017+0.012=0.067。
(5)所有4个次级圈闭均勘探失败,其概率为:
P0=(1-PⅠ(条件概率))×(1-PⅡ(条件概率))×(1-PⅢ(条件概率))×(1-PⅣ(条件概率))×P(边际概率)=(1-0.72)×(1-0.68)×(1-0.675)×(1-0.6)×0.68=0.012×0.68=0.008
由以上分析可知,T1断块圈闭获得勘探成功的概率就等于其4个次级圈闭至少有一个获得勘探成功的概率,即上述(1)-(4)四种情况的概率加和:
P=P4+P3+P2+P1=0.135+0.27+0.201+0.066=0.672
计算结果见表8:
表8
下面通过如下表达式对上面所求出的至少一层次级圈闭含油气的概率进行验证:
在本示例中,应用该公式的具体计算过程为:
P含油气概率=0.68×[1-(1-0.72)×(1-0.68)×(1-0.675)×(1-0.60)]≈0.68×(1-0.0115)=0.68×0.9885≈0.672
由此可知,通过本方法获取的圈闭的含油气概率不仅非常准确,并且还可以获取含油气的次级圈闭的各自可能的组合。
综上所述,本实施例的多层圈闭含油气概率的获取方法在获取多含油气概率时,充分考虑了各个油气聚集单元的成藏要素之间的互相依赖关系,并且获取各油气聚集单元可能的含油气组合以及不含油气组合的概率,在获取合理的圈闭的综合含油气概率同时,获得了各油气聚集单元含油气的组合,在工程勘探中具有实际的指导意义。
第二实施例
图5是根据另一实施例的多层圈闭含油气概率的获取系统的结构示意图。
如图5所示,该系统包括第一单元50、第二单元51和第三单元53。这三个单元分别对应执行第一实施例的步骤S110、S120和S130。在此不再赘述。
以上所述,仅为本发明的具体实施案例,本发明的保护范围并不局限于此,任何熟悉本技术的技术人员在本发明所述的技术规范内,对本发明的修改或替换,都应在本发明的保护范围之内。
Claims (10)
1.一种多层圈闭含油气概率的获取方法,包括:
步骤一,根据多层圈闭中各层圈闭的每个成藏条件之间的关系来确定各层圈闭的共性成藏条件和独立成藏条件,并获取各层圈闭的每个成藏条件的概率值;
步骤二,基于所述各个成藏条件的概率值计算多层圈闭中各层圈闭的边际概率和条件概率,进而得到各层圈闭的含油气概率;
步骤三,根据各层圈闭的边际概率值范围确定各层圈闭之间的决定关系,进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,作为多层圈闭的含油气概率。
2.根据权利要求1所述的获取方法,其特征在于,所述步骤三进一步包括:
采用融合有多层圈闭的边际概率和各层圈闭的条件概率的概率树,以获取各种层圈闭组合下的多层圈闭的含油气概率进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,其中,层圈闭组合基于决定关系和各层圈闭的含油气概率而得到。
3.根据权利要求2所述的获取方法,其特征在于,
所述决定关系包括:完全决定关系、部分决定关系以及完全独立关系,
在完全决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率等于各层圈闭中含油气概率的最大值,则对应含油气概率最大值的那层圈闭的含油气有无直接决定其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为完全相关的各种层圈闭组合下的含油气概率的加和;
在部分决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率大于各层圈闭中含油气概率的最大值并小于1,则对应含油气概率最大值的那层圈闭的含油气的有无间接反映其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为部分相关的各种层圈闭组合下的含油气概率的加和;
在完全独立关系中,各层圈闭间不存在共性成藏条件,该多层圈闭的边际概率为1,每一层圈闭的含油气有无仅取决于自身的成藏条件,所述多层圈闭的含油气概率为完全独立的各种层圈闭组合下的含油气概率的加和。
4.根据权利要求2所述的获取方法,其特征在于,
在所述概率树中,以多层圈闭的边际概率为对象,作为根节点,各层圈闭分别作为其他不同级别节点,且各级别节点的状态对应各层圈闭的条件概率。
5.根据权利要求1所述的获取方法,其特征在于,通过如下表达式计算多层圈闭中至少有一层发现油气的圈闭组合概率P含油气概率:
其中,P边际概率表示多层圈闭的边际概率,P第i个条件概率表示第i层圈闭的条件概率,n表示圈闭的层数。
6.根据权利要求2所述的获取方法,其特征在于,在步骤二中,
采用贝叶斯概率模型计算多层圈闭中各层圈闭的边际概率和条件概率;
基于各层圈闭的边际概率和条件概率得到各层圈闭的含油气概率,
其中,所述边际概率为各个共性成藏条件的概率值乘积,所述条件概率为各个独立成藏条件的概率值乘积,所述各层圈闭的含油气概率为非条件概率,其为所述条件概率和所述边际概率的乘积。
7.一种多层圈闭含油气概率的获取系统,包括:
第一单元,其根据多层圈闭中各层圈闭的每个成藏条件之间的关系来确定各层圈闭的共性成藏条件和独立成藏条件,并获取各层圈闭的每个成藏条件的概率值;
第二单元,其基于所述各个成藏条件的概率值计算多层圈闭中各层圈闭的边际概率和条件概率,进而得到各层圈闭的含油气概率;
第三单元,其根据各层圈闭的边际概率值范围确定各层圈闭之间的决定关系,进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,作为多层圈闭的含油气概率。
8.根据权利要求7所述的获取系统,其特征在于,所述第三单元进一步用于:
采用融合有多层圈闭的边际概率和各层圈闭的条件概率的概率树,以获取各种层圈闭组合下的多层圈闭的含油气概率进而得到多层圈闭中至少有一层发现油气的圈闭组合概率,其中,层圈闭组合基于决定关系和各层圈闭的含油气概率而得到。
9.根据权利要求8所述的获取系统,其特征在于,
所述决定关系包括:完全决定关系、部分决定关系以及完全独立关系,
在完全决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率等于各层圈闭中含油气概率的最大值,则对应含油气概率最大值的那层圈闭的含油气有无直接决定其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为完全相关的各种层圈闭组合下的含油气概率的加和;
在部分决定关系中,各层圈闭间存在至少一个共性成藏条件,且共性成藏条件相乘所得的概率大于各层圈闭中含油气概率的最大值并小于1,则对应含油气概率最大值的那层圈闭的含油气的有无间接反映其它层圈闭的含油气性有无,所述多层圈闭的含油气概率为部分相关的各种层圈闭组合下的含油气概率的加和;
在完全独立关系中,各层圈闭间不存在共性成藏条件,该多层圈闭的边际概率为1,每一层圈闭的含油气有无仅取决于自身的成藏条件,所述多层圈闭的含油气概率为完全独立的各种层圈闭组合下的含油气概率的加和。
10.根据权利要求8所述的获取系统,其特征在于,
在所述概率树中,以多层圈闭的边际概率为对象,作为根节点,各层圈闭分别作为其他不同级别节点,且各级别节点的状态对应各层圈闭的条件概率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410109198.7A CN104933282B (zh) | 2014-03-21 | 2014-03-21 | 多层圈闭含油气概率的获取方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410109198.7A CN104933282B (zh) | 2014-03-21 | 2014-03-21 | 多层圈闭含油气概率的获取方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104933282A true CN104933282A (zh) | 2015-09-23 |
CN104933282B CN104933282B (zh) | 2017-12-12 |
Family
ID=54120448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410109198.7A Active CN104933282B (zh) | 2014-03-21 | 2014-03-21 | 多层圈闭含油气概率的获取方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104933282B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107180293A (zh) * | 2016-03-11 | 2017-09-19 | 中国石油化工股份有限公司 | 一种面向勘探目标的地质评价水平测定方法 |
CN107292411A (zh) * | 2016-03-31 | 2017-10-24 | 中国石油化工股份有限公司 | 成藏组合技术可采储量规模预测方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1164651A (zh) * | 1996-10-29 | 1997-11-12 | 中国石油天然气总公司西北地质研究所 | 地震-测井多信息储层参数反演方法 |
-
2014
- 2014-03-21 CN CN201410109198.7A patent/CN104933282B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1164651A (zh) * | 1996-10-29 | 1997-11-12 | 中国石油天然气总公司西北地质研究所 | 地震-测井多信息储层参数反演方法 |
Non-Patent Citations (3)
Title |
---|
JAN-ERIKK ET AL.: "The CCOP guidelines for risk assessment of petroleum prospects", 《PROFEESSIONAL PAPER OF COORDINATING COMMITTEE FOR OFFSHORE PROSPECTING IN ASIA》 * |
盛秀杰: "成臧体系油气资源评价中的统计方法体系", 《石油与天然气地质》 * |
蔡利学 等: "新区预探圈闭优选的含油气概率与战略价值评价方法", 《大庆石油学院学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107180293A (zh) * | 2016-03-11 | 2017-09-19 | 中国石油化工股份有限公司 | 一种面向勘探目标的地质评价水平测定方法 |
CN107292411A (zh) * | 2016-03-31 | 2017-10-24 | 中国石油化工股份有限公司 | 成藏组合技术可采储量规模预测方法及装置 |
CN107292411B (zh) * | 2016-03-31 | 2020-11-13 | 中国石油化工股份有限公司 | 成藏组合技术可采储量规模预测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104933282B (zh) | 2017-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109478208B (zh) | 用于石油勘探和生产评估的综合数据和过程集成的迭代且可重复的工作流程 | |
US10767471B2 (en) | Resource density screening tool | |
US11693140B2 (en) | Identifying hydrocarbon reserves of a subterranean region using a reservoir earth model that models characteristics of the region | |
Rudolph et al. | Benchmarking exploration predictions and performance using 20+ yr of drilling results: One company’s experience | |
US7467044B2 (en) | Method and system for assessing exploration prospect risk and uncertainty | |
RU2496972C2 (ru) | Устройство, способ и система стохастического изучения пласта при нефтепромысловых операциях | |
US11486230B2 (en) | Allocating resources for implementing a well-planning process | |
EA011054B1 (ru) | Байесова сеть в применении к геологии и геофизике | |
US11815650B2 (en) | Optimization of well-planning process for identifying hydrocarbon reserves using an integrated multi-dimensional geological model | |
Bertolini et al. | A methodology to evaluate and reduce reservoir uncertainties using multivariate distribution | |
Yu | Hydraulic unit classification of un-cored intervals/wells and its influence on the productivity performance | |
CN104933282A (zh) | 多层圈闭含油气概率的获取方法和系统 | |
EP3948365A1 (en) | Automatic calibration of forward depositional models | |
Samimi et al. | SENSITIVITY & UNCERTAINTY ANALYSIS OF ORIGINAL OIL-IN-PLACE IN CARBONATE RESERVOIR MODELING, A CASE STUDY. | |
Meddaugh et al. | Reservoir forecast optimism–Impact of geostatistics, reservoir modeling, heterogeneity, and uncertainty | |
Popova | Analogy in the World of Geological Uncertainties, or How Reservoir Analogs May Refine Your Probabilistic Geomodel | |
Gorbovskaia | Probabilistic Geological Uncertainties Assessment: Lessons Learned | |
Singh et al. | Optimal Subsurface Appraisal: A Key Link to the Success of Development Projects-Few Examples | |
Tømmerås et al. | Prewell and postwell predictions of oil and gas columns using an iterative Monte Carlo technique with three-dimensional petroleum systems modeling | |
Becerra et al. | Petroleum reservoir uncertainty mitigation through the integration with production history matching | |
CN107038327B (zh) | 多目的层圈闭地质资源量确定方法及装置 | |
Abdul Majeed et al. | Volumetric equation approaches in Mishrif formation/Amara oil field: As a case study | |
Sokhin et al. | An Application of Direct Sequential Simulation and Co-simulation for Evaluating the Resource and Uncertainty of the Ncondezi Coal Deposit, Mozambique | |
Okotie et al. | Resources and Reserves | |
Lichang et al. | Application research of data mining on reservoir characterization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |