含油气性评价数学模型

如题所述

一般使用聚类分析等综合评判和专家系统对一个圈闭乃至一个区域的含油气性进行评价。

6.2.1 聚类分析

聚类分析又称群分析,簇分析。聚类分析的方法较多,主要是根据一批研究对象(样品或变量)在性质、特征、数量等方面的相似程度进行分类。油气资源评价中聚类分析通常思路是:选择部分勘探和研究程度高、含油气性已明确的圈闭、区块、盆地(包括含油气的和不含油气的)等作为标型单元,与勘探和研究程度相对较低,含油气性不确定的评价对象一起,依据这些样品(标型单元和评价对象)的多个变量指标(包括性质、特征等)应用多种方法确定距离系数或相似系数,然后编制分类谱系图,直观显示样品间相似程度,最后在不同相似性尺度上进行粗细不等的分类。

距离系数确定如下:对于几个样品,若我们选定m个变量指标进行分析,则可以把每个样品看成是m维空间上由其m个变量指标确定的一个点。两两样品(即点)之间的距离可由下式(即Minkowski距离公式)计算:

油气资源评价方法与实践

式中q为1时称为绝对(值)距离,为2时称为欧几里得距离。研究中通常应用欧几里得距离。

上述距离确定方法中存在两个问题:①距离与各变量的量纲有关,也就是说与变量的数值大小(量纲变化,数值会相应变化)有关,解决此难题的一个办法即是数据标准化,如标准差标准化、极差标准化、极差正规化等;②在计算距离时没考虑变量之间的相关性,解决这一问题的办法主要是要对变量进行分析,筛选主要变量,剔除次要变量和非独立变量。

表6-1 资源评价方法分类及其主要方法适用范围列表

相似系数通常指夹角余弦和相似系数。

夹角余弦定义:若把每个样品看成m维变量空间中的一个向量,那么样品i和样品j之间的相似性就可用两个向量间夹角余弦cosθij来表示,其计算公式为

油气资源评价方法与实践

相似系数定义:它是经过标准化后的夹角余弦,一般用rij表示,计算公式为

油气资源评价方法与实践

图6-1即是夹角余弦绘制的聚类分析谱系图,可见,该图非常明确地显示了17个样品间相似程度,当我们以cosθ=-0.2为划分界线时,可将17个样品分为2个组,但以cosθ=0为划分界线时,则分为3个组,若以cosθ=0.4为划分界线则又分为5个组,以此类推。

图6-1 聚类分析谱系图

6.2.2 判别分析

油气资源评价中,常常会遇到判断某个目标(如圈闭、区块、盆地)是否含油的问题,这就可用判别分析,即由一组已知样品,建立样品含油气性(或其他性质)与地质变量间的定量关系,即建立判别方程,并确定出类型归属的界线值。然后对于一个新样品(即未知样品,通常是评价目标),用该判别方程计算出判别值,再把此判别值与界线值比较,从而确定新样品的类型归属。

判别分析通常分为两组判别、多组判别和逐步判别。

两组判别分析是指把样品归类为两种,如一个勘探目标(如圈闭、区块)要么含油,要么不含油。当然,对于地质目标而言,用某个单项指标(即单一变量)进行判断是比较困难的,故一般都用多项指标(即多个变量)建立判别方程。为简化问题,判别方程一般是用线性方程,对于每个样品,当用m个可供判别分类的变量建立线性判别方程时,一般采用待定系数法。即设线性判别方程为

油气资源评价方法与实践

式中c1,c2,…,cm为待定系数。应用n个样品(其中n1个为1类,n2个为2类)的m个变量的观测值,按照两类样品间差别最大,而同类样品间差别最小的原则,通过解偏微分方程等一系列算法求取待定系数。

所谓多组判别就是指样品归属的类型多于两个,如某一勘探目标可能是大型油田,也可能是中型油田,也可能是小型油田甚至还可能不含油。多组判别方程建立的方法与两组判别方程建立方法相同,只是两组判别只一个判别函数,把空间分成两个域;三组判别有3个判别函数,把空间分成3个域;而4组判别需要6个判别函数,以此类推。因此,对于多组判别分析,计算组间的判别函数很不方便。

逐步判别就是对变量按其对判别分类的重要性,在计算过程中有进有出,保留那些对判别类型起主要作用的变量,剔除那些不起作用或作用不大的变量。

6.2.3 特征分析

在油气资源评价中,特征分析是通过已知圈闭、区块或盆地(统称为标型区)含油气性评价的主要控油地质变量,经最优化模型计算后,确定各地质变量的权系数,建立关联值计算公式,作为各地质变量具有的综合特征,依此进行未知盆地、区块、圈闭(统称为评价区)含油性评价。因此,特征分析更能反映事物的本质。特征分析法的基本思路是:通过标型区建立最优模型公式,并计算标型区和评价区的关联值,然后用下式计算标型区与评价区的相对偏差,即:

油气资源评价方法与实践

并定义为:|相对偏差|≤0.2,完全相似,地质类比系数为1;|相对偏差|=0.2~0.4,比较相似,地质类比系数为1.2或0.8;|相对偏差|=0.4~0.6,基本相似,地质类比系数1.4或0.6;|相对偏差|>0.6,基本不相似,不能类比。

特征分析法一般评价步骤是:

(1)依据各项地质变量含油气性判别的准则,确定评价的阈限值,将各项地质参数转换成1、0、-1三元逻辑表达式。

(2)用乘积矩阵的平方和法、主分量法及匹配概率矩阵主分量法等对标型区进行最优化模型特征的定量化,确立各项地质变量的权系数,建立关联值计算公式:

油气资源评价方法与实践

式中:y为关联值,

xi项地质变量的三元逻辑值,

ai为xi的权系数。

(3)计算标型区与待判区的关联值,并进行类比,确定未知盆地、区块、圈闭的含油气性。类比准则是:凡与已知油气田(区块、盆地)关联值可比的列入Ⅰ级含油气圈闭(区块、盆地),与已知含油气圈闭(区块、盆地)关联值可比的列入Ⅱ级含油气圈闭(区块、盆地);关联值介于含油与非含油圈闭(区块、盆地)之间的,列入Ⅲ级含油气圈闭(区块、盆地);关联值在非含油气圈闭(区块、盆地)范围内的,则圈闭(区块、盆地)不含油气。

6.2.4 逻辑信息法

逻辑信息法属于“数量化”理论的范畴,是数理逻辑、组合分析及数理统计为基础的综合评价方法。主要是通过已知不同级别含油气性的圈闭(如已知油气田、含油气圈闭、有油气显示圈闭和干圈闭),按一定序列建立圈闭含油气性评价的变异序列,以此作为评价模型。依据“位移帕斯卡三角形”,求得各地质参数的相对标志权,挑选评价的主要地质变量,计算各地质变量的标志信息权、标志分权及各圈闭的对象权,然后与已知圈闭对象权分布范围进行比较,作出圈闭含油气性评价。具体步骤如下:

(1)选择已知圈闭为标型,对各控油地质变量进行0,1二态变量转换,构成m组变异序列。

(2)依据“位移帕斯卡三角形”计算各地质变量取值为1的相对标志权Pk*(1)及取值为0的相对标志权Pk*(0)(k为地质变量,k=1,2,…,p),然后将m个变异序列中相同级别的圈闭进行相互调换,分别计算各地质变量的Pk*(1)、Pk*(0)。若在所有变异序列中几乎都是Pk*(1)≥-Pk*(O)或Pk*(0)≥Pk*(1),则地质变量k与圈闭含油性关系密切,挑选该项变量,反之则剔除。

(3)对已选地质变量,划分最小区分标志组合,计算各标志信息权Pk:Pk=包含k标志的最小区分标志组合数/总的最小区分标志组合数信息权愈大,反映该标志对圈闭含油气性评价意义愈大。

(4)依据各地质变量相对标志权及信息权计算标志分权及标志分权的绝对差。标志分权Pk(1)及Pk(0)按下式求得

Pk(1)=Pk(1)·Pk Pk(0)=Pk(0)·Pk

标志分权的绝对差Rk=|Pk(1)-Pk(0)|

若Pk(1)>Pk(0),则将分权差Rk赋予1,而将0值赋予0;

若Pk(1)<Pk(0),则将分权差Rk赋予0,而将0值赋予1。

(5)据各圈闭每个标志的取值,应用下式求得各圈闭对象权:

油气资源评价方法与实践

式中:Ii为对象权(i=1,2,…,m)

Rk(1,0):据第i个圈闭k标志的取值(0或1)而赋予Rk(0)或Rk(1)(k=1,…,n)。

(6)将各圈闭对象权与标型圈闭的对象权分布范围进行类比,按其相应范围,评价圈闭含油气性。

6.2.5 油气资源评价专家系统

油气资源评价经历了几十年的发展,而今已形成了以统计学为基础的综合评价系统和以数字模拟为基础的盆地模拟系统,但此二系统适用范围局限,当定量资料欠缺时无法评价或评价结果可靠性很低,为充分利用各种定量和定性资料,必须引入专家系统技术。

所谓专家系统即是应用专家知识,按照专家思路解决某一特定范围内的专业问题,得出与专家分析类似的结论的计算机软件系统。

在我国,油气资源评价专家系统的研制还是近十几年的事,已研制的较成熟油气资源评价专家系统有两个:一个是海洋石油勘探开发研究中心与吉林大学“七五”期间合作研制的油气资源评价专家系统(Petroleum Resources Evaluation System,简称PRES),包括6个推理及程序模块、9个知识库和3个数据库,可在PC机和VAX11/78O机上运行,对凹陷与圈闭进行资源评价。另一个是地矿部石油地质研究所等多家合作研制的油气资源评价专家系统,其评价对象主要是圈闭,但也可对凹陷进行评价。

严格地讲,不管国内还是国外,目前油气资源评价专家系统尚属发展趋向,处在开发试验阶段,成功者很少,上述二系统即是成功者之一。自1994年下半年以来,地矿部研制的油气资源评价专家系统分别在东海陆架盆地西湖凹陷、松辽盆地南部地区、四川盆地川东北地区、鄂尔多斯盆地伊陕斜坡、南沙海域万安盆地进行了推广应用,特别是在新星公司圈闭管理中得到了广泛应用,取到了令人满意的结果,系统的先进性、科学性和合理性正越来越得到专家们的承认。下面即对该专家系统进行介绍。

6.2.5.1 油气资源评价专家系统简介

油气资源评价专家系统由知识库、知识库维护模块、综合数据库、推理机、解释模块组成(图6-2)。其中知识库、推理机和解释模块是核心,实现评价模型的具体算法。

图6-2 专家系统结构示意图

知识库主要是规则库。规则库之基本构件是节点(地质事实)和边(地质事实间相互关系),边依产生式原则和框架式结构将节点连接起来,边与边之间的关系或用特性表示,或用“与”、“或”、“异或”、“加权”、“新与”和“独立”六种关系表示,并加上边信度(即规则信度),形成有机的网络,其中最底层节点没有前提,叫叶节点;顶层节点没有结论,叫根节点;其余叫中间节点(如图6-3所示)。该专家系统就是以此网络实现评价模型的。建立评价模型,构造知识库特别是规则库,是建立专家系统的中心任务,也是难点所在,即所谓“瓶颈”问题。由于各地区气藏成藏模式不同,故知识库采取模块化结构,推理时,根据具体研究区实际成藏模式,调入相应的模块生成一个适合本区的临时知识库,以提高推理效率。现在系统已依据松辽盆地南部东南隆起区泉头组—登楼库组、东海陆架盆地西湖凹陷平湖组—花港组、鄂尔多斯盆地伊陕斜坡奥陶系顶部风化壳、四川盆地川东北地区二叠系—中三叠统和川西坳陷上三叠统5种不同的成藏模式“调入生成”了相应含油气性评价的临时知识库。

图6-3 节点关系示意图

该专家系统推理机推理为不确定推理,推理方式主要为向前推理,类似人类专家思维那样,依据叶节点或中间节点(前提)按一定规则推理、判断另一些中间节点或根节点(结论)成立的可能性,其实质是进行大量的信度运算,以此实现评价过程。该系统的解释模块以人机会话方式,在推理的一些重要环节允许用户询问并向用户解答诸如“本中间结论是怎样得到的?”、“为什么需要在这里提供信息?”等问题,增加了推理透明度,使用户能更好地参与推理。

6.2.5.2 评价模型在专家系统里的实现

6.2.5.2.1 规则库建立

按产生式规则和框架式结构,应用特性表和“与”、“或”、“异或”、“加权”、“新与”、“独立”6种关系将具体研究区利用典型油气藏解剖、油气藏形成模式及控制因素分析、油气成藏条件及油气藏分布规律研究等成果而建立的评价模型中各种规则有机连接起来,并输入专家系统知识库内,即建成一个规则库。在规则库建立中,采取了模块化结构,分为总模块、含油气标志子模块、油气源条件子模块、储集条件子模块、圈闭条件子模块、聚集条件子模块和保存条件子模块共7个模块,各模块结构如图6-4~6-10所示。

6.2.5.2.2 推理方式及6种推理算法

专家系统评价的实质就是根据规则库里一系列规则及其相互关系(即评价模型)进行推理。本专家系统主要推理方式是不确定推理,即结论信度可在(-1,1)范围内连续取值,符合地质现象间关系复杂,多为不确定性的这一特点。按照推理方向以向前推理为主,辅以反向推理的设计思路,系统选用了大量的计算方法模型,其中主要的是推理算法,以下仅就6种逻辑关系的推理算法作一简介:

图6-4 总模块结构图

图6-5A 含油气标志子模块结构图

图6-5B 含油气标志子模块结构图

(1)只有一个前提的情形

若H是结论,E是前提,并假定E的信度cf(E)和从E到H的边信度a,b已知,则:

图6-5C 含油气标志子模块结构图

1)若cf(E)>0,则cf(H)=a·cf(E);

2)若cf(E)<0,则cf(H)=b·cf(E)。

由于a,b取值都可在[-1,1]之间变化,因此可以反映E对H的多种影响,例如E是H的充分条件(a=1,b=0),必要条件(a=0,b=1),充要条件(a=b=1)等等。

(2)具有多个前提的情形

设H是结论,且,E1,E2,…,En是前提,cfi(H)表示由单一前提Ei所求出的H的信度,i=1,2,…,n,则H的最终信度cf(H)计算方法如下:

1)若诸Ei之间的关系为“与”,则令

油气资源评价方法与实践

2)若诸Ei之间的关系为“或”,则令

油气资源评价方法与实践

3)若诸Ei之间的关系为“加权”,则令

油气资源评价方法与实践

图6-6 油气源条件子模块结构图

诸Wi是权值。在特殊情形下,规定cf(H)》ΣWi·cfi(H),以表示“突变”关系。

4)若诸Ei之间的关系为“异或”,则诸Ei中必定只有且恰有一个满足cf(H)>t>0(t是阈值),因此,可令

油气资源评价方法与实践

图6-7A 储集条件子模块结构图

图6-7B 储集条件子模块结构图

图6-8 圈闭条件子模块结构图

5)若诸Ei之间的关系为“独立”,为了简单起见,先设n=2,则令

油气资源评价方法与实践

当n>2时,先将同号的诸cfi(H)分别按上式合并,再求其代数和得到cf(H)。

6)若诸Ei之间的关系是“新与”,且Ei为主要节点(起决定性作用),则令:

油气资源评价方法与实践

图6-9 聚集条件子模块结构图

如果一个结论H的诸前提E1,E2,…,En之间存在着多种关系,则通过设置运算节点(虚节点)来反映各种关系之间的优先级。

如果对信度计算有特定要求,例如需经查表、计算(调用子程序)、类比或询问来确定结论H的信度值,则知识库中标明之。另外,对于边信度可能需要修改的情形,某个前提可以单独肯定或否定结论的情形等,也在规则库中注明,以便通过有关元知识及时处理。

图6-10A 保存条件子模块结构图

6.2.5.2.3 解释功能

当系统推出某个重要中间结果时,系统显示这一结论并暂停推理,等待用户提问并予以回答;另一方面,当推理到某一步时,系统需要用户提供某种信息,也允许用户在输入信息前询问并进行回答,然后等候用户输入此信息;最后,系统还能通过列表方式显示最终结论或中间结论是如何推得的。

图6-10B 保存条件子模块结构图

温馨提示:答案为网友推荐,仅供参考
相似回答