第八章蒙特卡罗模拟

1、1第八章第八章 蒙特卡罗模拟法蒙特卡罗模拟法1 蒙特卡罗法概述2 随机数及其产生方法3 随机变量的抽样4 油气资源量的估算5 地质风险分析6 应用简例21 蒙特卡罗法概述蒙特卡罗法概述 蒙特卡罗法蒙特卡罗法(Monte Carle):以数值解不确定问题为以数值解不确定问题为对象，对计算模型中的各变量进行随机抽样对象，对计算模型中的各变量进行随机抽样(随机随机试验试验)，进而求问题概率解的一种统计学方法。因，进而求问题概率解的一种统计学方法。因此，蒙特卡罗法又称为统计试验法。此，蒙特卡罗法又称为统计试验法。 蒙特卡罗法是以概率论与数理统计理论为指导蒙特卡罗法是以概率论与数理统计理论为指导的、有着

2、广泛应用领域的通用性统计学方法，其的、有着广泛应用领域的通用性统计学方法，其核心是对随机变量的抽样。常用于油气资源量预核心是对随机变量的抽样。常用于油气资源量预测。测。 美国第二次全美石油资源评价的主要方法就是蒙美国第二次全美石油资源评价的主要方法就是蒙特卡罗法。目前各主要产油国及各大石油公司都特卡罗法。目前各主要产油国及各大石油公司都把该法作为油气资源评价的重要方法之一。把该法作为油气资源评价的重要方法之一。3 如何利用蒙特卡罗法进行油气资源量的计算？如何利用蒙特卡罗法进行油气资源量的计算？在概念模型清楚的基础上，首先要构造表示所研在概念模型清楚的基础上，首先要构造表示所研究问题概率解的数学

3、模型究问题概率解的数学模型(计算模型计算模型)。 我国我国1979年把蒙特卡罗法用于油气资源评价年把蒙特卡罗法用于油气资源评价,该该方法目前仍是我国预测油气资源量的常用方法之方法目前仍是我国预测油气资源量的常用方法之一一,也是以统计预测为主的应用软件评价系统的核也是以统计预测为主的应用软件评价系统的核心算法之一。心算法之一。 对一个实际地区的油气资源总量而言，它是各对一个实际地区的油气资源总量而言，它是各局部含油气地质单元油气资源量之和。而局部含局部含油气地质单元油气资源量之和。而局部含油气地质单元既可以是生油凹陷中的一个生油层油气地质单元既可以是生油凹陷中的一个生油层系，也可以是次一级构造单

4、元中的生油层系，还系，也可以是次一级构造单元中的生油层系，还可以是局部构造等。可以是局部构造等。4(8-1)nljljvinljlijXcXDQ111油气总资源量：油气总资源量：mjjQQ1 Qj -第第j个局部地质单元油气资源量个局部地质单元油气资源量( j=1,2, m); Xjl -第第j个局部地质单元第个局部地质单元第l个地质变量个地质变量; cj -第第j个局部地质单元中个局部地质单元中v个地质常数个地质常数Di之积。之积。 虽然虽然“局部地质单元局部地质单元”含义不同，但任何一个含义不同，但任何一个局部含油气地质单元的油气资源量都可归纳为与局部含油气地质单元的油气资源量都可归纳为与

5、油气资源量相关的地质常数和变量的乘积，即油气资源量相关的地质常数和变量的乘积，即: 因此，求油气资源量的问题，就归结为求上述因此，求油气资源量的问题，就归结为求上述两个计算模型的概率解问题。两个计算模型的概率解问题。(8-2)5一、蒙特卡罗法的基本思想一、蒙特卡罗法的基本思想 蒙特卡罗的基本思想可概括为：为求研究问题蒙特卡罗的基本思想可概括为：为求研究问题的概率解的概率解, 构造一个表示所研究问题概率解的数学构造一个表示所研究问题概率解的数学模型模型(计算模型计算模型)，记为，记为:),(21nXXXY(8-3) 依据计算模型中各随机变量依据计算模型中各随机变量Xi所服从的分布进行所服从的分布

6、进行随机抽样，并按计算模型计算随机抽样，并按计算模型计算Y的多个估计值，最的多个估计值，最终用终用频率统计法频率统计法求出求出Y的概率解。的概率解。 方法的核心是随机抽样，而随机抽样的关键在方法的核心是随机抽样，而随机抽样的关键在于产生于产生0,1区间上均匀分布的随机数。因为服从其区间上均匀分布的随机数。因为服从其它分布的随机数一般可通过该随机数变换得到。它分布的随机数一般可通过该随机数变换得到。6二、蒙特卡罗法求解的基本过程二、蒙特卡罗法求解的基本过程由基本思想，求解过程大致分为四步：由基本思想，求解过程大致分为四步： 分析并拟定给定问题中的随机变量，构造表分析并拟定给定问题中的随机变量，构

7、造表示给定问题概率解的数学模型；示给定问题概率解的数学模型； 对模型中的随机变量对模型中的随机变量X1, X2, Xn各进行各进行L次随次随机抽样，获得机抽样，获得L组抽样值：组抽样值：x1k , x2k , xnk (k=1,2, ,L) 把把L组抽样值代入计算模型，求出随机变量组抽样值代入计算模型，求出随机变量Y的的L个估计值个估计值Y1 , Y2 , , YL； 利用频率统计法，由利用频率统计法，由Y的估计值的估计值Y1 , Y2 , , YL求求出描述出描述Y分布特征的分布曲线分布特征的分布曲线(概率解概率解) ，如图，如图8-1。 7图图8-1 随机变量随机变量Y的分布曲线的分布曲线

8、yP(Yy)10.50P=0.582 随机数及其产生方法随机数及其产生方法一、随机数一、随机数 随机数是随机变量的观测值，由其构成的数据随机数是随机变量的观测值，由其构成的数据序列叫做随机数序列，它是一个无周期的数据序序列叫做随机数序列，它是一个无周期的数据序列。列。 蒙特卡罗法需要对变量进行数以千计、万计、蒙特卡罗法需要对变量进行数以千计、万计、甚至是百万计的抽样，这在实现过程中几乎是不甚至是百万计的抽样，这在实现过程中几乎是不可能的。因此考虑用计算机模仿实际抽样过程，可能的。因此考虑用计算机模仿实际抽样过程，形成一个有周期的抽样数据序列。称这种数据序形成一个有周期的抽样数据序列。称这种数据

9、序列为列为“伪随机数伪随机数”序列，其中的元素叫做序列，其中的元素叫做“伪随伪随机数机数”。 伪随机数显然不是真正意义下的随机数，即这伪随机数显然不是真正意义下的随机数，即这种抽样值并非是随机变量的真实观测值。种抽样值并非是随机变量的真实观测值。9 尽管如此，只要对伪随机数序列进行一系列严尽管如此，只要对伪随机数序列进行一系列严格的统计检验，证明它可以满足统计的要求，则格的统计检验，证明它可以满足统计的要求，则伪随机数就可以做为真随机数使用。伪随机数就可以做为真随机数使用。 为了满足抽样问题的需要，在计算机上产生的为了满足抽样问题的需要，在计算机上产生的伪随机数序列不仅要有伪随机数序列不仅要有

10、足够长的周期足够长的周期，而且应当，而且应当具有符合要求的具有符合要求的概率统计性质概率统计性质。 从理论上讲，只要有一种连续分布的随机数，从理论上讲，只要有一种连续分布的随机数，就可以采用数学变换的方法产生其它分布的随机就可以采用数学变换的方法产生其它分布的随机数。数。0,1区间上均匀分布的随机变量的抽样值是最区间上均匀分布的随机变量的抽样值是最简单、最基本的一种连续分布的随机数，其它分简单、最基本的一种连续分布的随机数，其它分布的随机数都可以借助它来产生，所以说：布的随机数都可以借助它来产生，所以说：0,1均均匀分布随机数技术是实现随机抽样的最基本工具。匀分布随机数技术是实现随机抽样的最基

11、本工具。 10 xn , xn+1第第n次和第次和第n+1次产生的伪随机数；次产生的伪随机数； -乘子系数；乘子系数；M- 模；模； rn+1- 0,1区间上的伪随机数。区间上的伪随机数。 xn+1xn (mod M)叫做以叫做以 M 为模的同余式，表示为模的同余式，表示xn+1取值为：取值为：与与xn的积除以的积除以 M 的余数部分。的余数部分。 该方法产生伪随机数序列的递推同余式为：该方法产生伪随机数序列的递推同余式为：MxrMxxnnnn/)(mod111 二、伪随机数的产生方法二、伪随机数的产生方法 1. 乘同余法乘同余法112. 混合同余法混合同余法 该方法产生伪随机数序列的递推同余

12、式如下：该方法产生伪随机数序列的递推同余式如下：MxrMxxnnnn/)(mod111 混合同余法比乘同余法仅是增加了一个增量混合同余法比乘同余法仅是增加了一个增量,其其它含义与乘同余法相同。它含义与乘同余法相同。 如：如：M =219 = 524288 , = 55 = 3125 时，时，x0= 23 , 11 ,19 , 37;= 3 , 7 , 11 , 17分分4套配合使用，混合同余法可产套配合使用，混合同余法可产生周期为生周期为524288伪随机数序列。伪随机数序列。 注注:所获得的伪随机数是否能代表真正意义上的所获得的伪随机数是否能代表真正意义上的随机数随机数,还需进行检验还需进行

13、检验(检验方法略检验方法略)。12 一、随机变量的经验分布函数一、随机变量的经验分布函数 在概率论中，随机变在概率论中，随机变量量X的分布函数是随机变的分布函数是随机变量量X的取值不大于实数的取值不大于实数 x的概率。通常记为：的概率。通常记为： F(x)=P(Xx)图图8-2 随机变量分布曲线随机变量分布曲线1F(x)x0如图如图8-2所示。所示。3 随机变量的抽样随机变量的抽样 对随机变量的抽样，有经验分布函数抽样法、对随机变量的抽样，有经验分布函数抽样法、直接抽样法和变换抽样法等。在此仅介绍经验分直接抽样法和变换抽样法等。在此仅介绍经验分布函数抽样法布函数抽样法。13 经验分布函数是由经

14、验分布函数是由X的的 n 个观测值个观测值x1, x2, xn ,用统计方法得到的分布函数，记为用统计方法得到的分布函数，记为Fn(x) 。如图如图8-3。 但是，在油气资源评但是，在油气资源评价中价中,人们总希望得到资人们总希望得到资源量不小于某个实数源量不小于某个实数 x 的概率的概率P(Xx), 显然显然:P(Xx)=1-P(Xx)=1-F(x)图图8-3 分布函数曲线分布函数曲线P(Xx)01x14)(1)()(xFxXPxAFn(8-4)以统计所得的以统计所得的Fn(x)代替代替F (x), 并记：并记：AF(x)就就是油气资源评价中的经验分布函数。是油气资源评价中的经验分布函数。二

15、、随机变量经验分布函数的构造方法二、随机变量经验分布函数的构造方法 频率统计法构造经验分布函数的条件及步骤频率统计法构造经验分布函数的条件及步骤: (1)使用条件：使用条件： 随机变量必须有观测值随机变量必须有观测值, 且个数足够多且个数足够多(如如30)。 (2)构造分布函数的步骤构造分布函数的步骤 确定频率统计区间数确定频率统计区间数 要将要将n个观测值所在的大区间划分为个观测值所在的大区间划分为k个小区间个小区间 ,一般一般n/k3 ，且，且k最好为奇数。最好为奇数。15 计算区间端值计算区间端值 将将(xmin , xmax) k等分等分, 各小区间的各小区间的k+1个端值为：个端值为

16、：) 1, 2, 1() 1(minmaxminkiikxxxxi其中其中xmin , xmax是观测值的最大和最小值。是观测值的最大和最小值。 经验分布函数经验分布函数 记记ni 为观测值落入区间为观测值落入区间(x i , xi+1)内的频数内的频数(个数个数), fi 为累加频率，则为累加频率，则 kijjikinnf),2, 1(1显然显然 f1=1.0 , fk=nk/n , f1 f2fk n=n1 +n2 +nk16kjjnnf111kjjnnf221AF(x)01xx1x2x3xkxk+1nnfkk图图 8-4 累计频率直方图累计频率直方图17图图8-5 经验分布函数曲线经验分

17、布函数曲线AF(x)01x.,0;,;,;,1)(1132221xxxxxfxxxfxxxxAFkkkk经验分布函数表示为：经验分布函数表示为：18 将将坐标原点坐标原点设设为为(xmin ,0)，若已知，若已知0,1区间上均匀区间上均匀分布的随机数分布的随机数ri ，则在图，则在图8-6纵轴上可确定点纵轴上可确定点(xmin , ri ), 过该点作横轴的平行线交分布曲线于点过该点作横轴的平行线交分布曲线于点(xi , ri ), xi 则则是对应于随机数是对应于随机数 ri的一次随机抽样值。的一次随机抽样值。 三、经验分布函数的抽样三、经验分布函数的抽样 由此可得随机变由此可得随机变量的一

18、系列抽样值。量的一系列抽样值。称该抽样方法为经称该抽样方法为经验分布函数抽样法。验分布函数抽样法。图图8-6 分布函数抽样过程示意图分布函数抽样过程示意图(xmin , 0)x i (出口值出口值)(xi , r i) x min ,r i (入口值入口值)xAF (x)(xmin , 1)19 如前所述如前所述,局部含油气地质单元是估算资源量的局部含油气地质单元是估算资源量的基本地质体。对第基本地质体。对第j个局部地质单元来说，其油气个局部地质单元来说，其油气资源量的通式为资源量的通式为:4 油气资源量的估算油气资源量的估算一、局部地质单元资源量概率模型一、局部地质单元资源量概率模型niji

19、jjXcQ1式中式中: Qj -第第j个局部地质单元的油气资源量个局部地质单元的油气资源量; Xji -第第j个局部地质单元第个局部地质单元第i个地质随机变量个地质随机变量; cj -第第j个局部地质单元中所有常数之积。个局部地质单元中所有常数之积。 20 二、资源量的计算二、资源量的计算 1. 最大最小可能资源量及累计频率区间间隔值最大最小可能资源量及累计频率区间间隔值 为了求资源量的分布曲线，在资源量计算之前，为了求资源量的分布曲线，在资源量计算之前，先求出资源量的最大、最小可能值和累积频率小先求出资源量的最大、最小可能值和累积频率小区间的端值区间的端值:),2,1(1maxmaxmjxc

20、Qnijijj),2, 1(1minminmjxcQnijijj) 1, 2,1() 1(minmaxminkhhkQQQQjjjjh最小可能值最小可能值: 式中式中k是小区间个数。是小区间个数。最大可能值最大可能值:区区 间间 端端 值值: 21 对对n个随机变量个随机变量Xji (i=1,2,n)各进行一次抽样各进行一次抽样,得第一次抽样值得第一次抽样值xji1 (i=1,2,n),由第一次抽样值的由第一次抽样值的积得到资源量的第积得到资源量的第1个估计值，即：个估计值，即：nijijjxkQ111 2. 资源量的计算资源量的计算 对随机变量对随机变量 Xji (i=1,2, n)再各进行

21、一次抽样，得再各进行一次抽样，得第二次抽样值第二次抽样值xji2 (i=1,2,n),由第二次抽样值的积得由第二次抽样值的积得到资源量的第到资源量的第2个估计值，即：个估计值，即：nijijjxkQ12222jNjjjjQQQQQ,4321按上述做法，可得按上述做法，可得: 第第 j 个局部地质单元的资源量的分布函数个局部地质单元的资源量的分布函数 根据区间间隔值，用频率统计的方法求出根据区间间隔值，用频率统计的方法求出Qj的分的分布函数布函数AF(Qj)，如图，如图8-7。图图 8-7 资源量分布函数资源量分布函数QjAF(Qj )Qjmin23图图8-8 油气资源量的概率加示意图油气资源量

22、的概率加示意图三、总资源量三、总资源量 总资源量是局部地质单元资源量的概率加总资源量是局部地质单元资源量的概率加(对各对各资源量分布函数的随机抽样加资源量分布函数的随机抽样加)。资源量概率加的大致过程如下图：资源量概率加的大致过程如下图：Q1Q2QmQ=+总资源量的分布函数总资源量的分布函数24四、估算资源量的基本过程四、估算资源量的基本过程v对于一个局部地质单元：对于一个局部地质单元：选择预测选择预测(估算估算)方法方法(建立概率估算模型建立概率估算模型)；确定参数中的随机变量；确定参数中的随机变量；构造随机变量的分布函数；构造随机变量的分布函数； 对各随机变量的分布函数进行重复抽样，计算对

23、各随机变量的分布函数进行重复抽样，计算出资源量的多个估计值；出资源量的多个估计值；由资源量的多个估计值求资源量分布函数。由资源量的多个估计值求资源量分布函数。v对于对于m(m1)个局部地质单元：个局部地质单元： 需要概率加求总的资源量需要概率加求总的资源量(局部地质单元资源量局部地质单元资源量分布函数是求总资源量的基础分布函数是求总资源量的基础)。25mQQQQ21niiXcQ1111niiXcQ1222nimimmXcQ1构造各随机变量经构造各随机变量经验分布函数验分布函数构造各随机变量经构造各随机变量经验分布函数验分布函数构造各随机变量经构造各随机变量经验分布函数验分布函数Q1经验分布函数

24、经验分布函数Q2经验分布函数经验分布函数Qm经验分布函数经验分布函数Q的经验分布函数的经验分布函数(概率解概率解)各随机变量抽样各随机变量抽样各随机变量抽样各随机变量抽样各随机变量抽样各随机变量抽样26 由于石油勘探的未来成效具有不确定性，因而由于石油勘探的未来成效具有不确定性，因而需要对估算的石油资源量进行风险分析。石油勘需要对估算的石油资源量进行风险分析。石油勘探中的风险是多种多样的，如勘探区是否具备形探中的风险是多种多样的，如勘探区是否具备形成油气藏的地质条件的风险；在具备形成油气藏成油气藏的地质条件的风险；在具备形成油气藏地质条件的含油气区内，能否找到一定规模油气地质条件的含油气区内，

25、能否找到一定规模油气藏的勘探风险；勘探后发现的油气藏是否具备开藏的勘探风险；勘探后发现的油气藏是否具备开采价值的经济风险；勘探过程中人与设备是否安采价值的经济风险；勘探过程中人与设备是否安全的风险；对于勘探地区，特别是大陆架地区是全的风险；对于勘探地区，特别是大陆架地区是否有国际争议的政治风险等。否有国际争议的政治风险等。5 地质风险分析地质风险分析一、风险分析一、风险分析27 油气资源评价为油气勘探提供决策依据，而油油气资源评价为油气勘探提供决策依据，而油气资源估算是资源评价的组成部分。因此，要对气资源估算是资源评价的组成部分。因此，要对上述风险下估算的油气资源量做好地质风险分析。上述风险下

26、估算的油气资源量做好地质风险分析。 地质风险分析是指对模型中的参数做进一步的地质风险分析是指对模型中的参数做进一步的深入研究，对各参数存在的可能性做出概率估计，深入研究，对各参数存在的可能性做出概率估计，在此基础上对预测的资源量进行可靠性分析。在此基础上对预测的资源量进行可靠性分析。 实际工作中，地质风险分析可以在不同的层次实际工作中，地质风险分析可以在不同的层次进行，如单一地质圈闭的风险分析，一组地质圈进行，如单一地质圈闭的风险分析，一组地质圈闭的风险分析闭的风险分析(国外对地质条件相似的一组地质圈国外对地质条件相似的一组地质圈闭称作一个勘探层闭称作一个勘探层)，一个油气聚集带的风险分析，一

27、个油气聚集带的风险分析，整个含油气盆地的风险分析等。整个含油气盆地的风险分析等。 地质风险分析大多从圈闭开始进行，其风险的地质风险分析大多从圈闭开始进行，其风险的计算公式为：计算公式为：28 式中式中K是圈闭的风险值是圈闭的风险值, ki是第是第i个因素的风险值。个因素的风险值。例如，若用容积法估算圈闭的石油储量，则有：例如，若用容积法估算圈闭的石油储量，则有：WDHSQ 对储量进行风险分析，要由熟悉工区的地质人员对储量进行风险分析，要由熟悉工区的地质人员对上述对上述5项参数的风险逐个分析论证。通常来说：项参数的风险逐个分析论证。通常来说：niikK111)( Q-石油储量石油储量 S-含油面

28、积含油面积 H-储层厚度储层厚度-储层孔隙度储层孔隙度 D- 石油充满系数石油充满系数 W- 采收率采收率29 S(含油面积含油面积)的风险决定于地质调查或地震勘探的风险决定于地质调查或地震勘探资料的可靠性；资料的可靠性； 的风险决定于储集层孔隙是否有次生改造或后的风险决定于储集层孔隙是否有次生改造或后期充填的影响；期充填的影响； H(储层厚度储层厚度)的风险受岩性岩相变化的影响；的风险受岩性岩相变化的影响； D(充满系数充满系数)的风险受生油岩的成熟度和油气运的风险受生油岩的成熟度和油气运移通道的制约移通道的制约; W(采收率采收率)的风险与原油性质及驱动类型有关。的风险与原油性质及驱动类型

29、有关。 对每个参数的风险值对每个参数的风险值k(0k1)，目前尚无完善的，目前尚无完善的方法确定。一种方法是由地质人员凭经验指定；另方法确定。一种方法是由地质人员凭经验指定；另一种方法是借用地质条件类似的邻区风险值。如某一种方法是借用地质条件类似的邻区风险值。如某地质圈闭经风险分析论证后给出如下风险值地质圈闭经风险分析论证后给出如下风险值:30地地 质质 参参 数数 风险值风险值/k保险值保险值/(1-k) 含油面积含油面积(S) 储层厚度储层厚度(H) 储层孔隙度储层孔隙度() 充满系数充满系数(D) 采收率采收率(W)00.500.3010.510.71 根据表内数据，计算该圈闭的风险系数

30、为根据表内数据，计算该圈闭的风险系数为0.65，而保险值为而保险值为0.35，即经过风险分析后该地质圈闭的，即经过风险分析后该地质圈闭的石油储量的可能性仅是风险分析前的石油储量的可能性仅是风险分析前的35%。表表8-1 某个地质圈闭风险值某个地质圈闭风险值31图图8-9 风险分析前后石油储量关系风险分析前后石油储量关系风险分析前风险分析前风险分析后风险分析后0.51.00.0F(q)q32 若第若第j个局部单元油气资源量分布函数为个局部单元油气资源量分布函数为AF(qj)，地质风险值为地质风险值为kj ，保险值为，保险值为1-kj，则需重新估算风，则需重新估算风险分析后的油气资源量。既在原分布

31、函数基础上险分析后的油气资源量。既在原分布函数基础上再次抽样，入口值区间由再次抽样，入口值区间由0,1变为变为0 ,1-kj，如图。如图。图图8-10 风险分析前后石油储量关系风险分析前后石油储量关系二、风险分析后资源总量的估算二、风险分析后资源总量的估算风险分析前风险分析前风险分析后风险分析后0.51.00.0F(q)q1-kj33 将将0,1上分布的随机数上分布的随机数rs改造为改造为0,1-kj上分布的随上分布的随机数机数rs*的变换公式为的变换公式为rs*=(1-kj )rs。然后进行随机抽。然后进行随机抽样，求出风险分析后局部单元油气资源量的值样，求出风险分析后局部单元油气资源量的值

32、: 求得求得qj(=1,2, N )后，可用频率统计法求出各局后，可用频率统计法求出各局部单元的资源量分布函数部单元的资源量分布函数AF(qj)( j=1,2, m)，最后，最后再利用再利用“抽样加抽样加”和和“频率统计法频率统计法”求出地质风求出地质风险分析后资源总量的分布函数险分析后资源总量的分布函数AF*(q) 。AF*(q)=AF(q2)+AF(q2)+AF(qm)qj(=1,2, N )346 应用简例应用简例21kkADHSQjjjj 某沉积盆地中的一个地质凹陷有三套生油层系。某沉积盆地中的一个地质凹陷有三套生油层系。用氯仿沥青法估算该凹陷的远景石油资源量，各用氯仿沥青法估算该凹陷

33、的远景石油资源量，各层系的石油资源量估算公式如下：层系的石油资源量估算公式如下： Qj - 第第j( j=1,2,3)套生油层系的石油资源量套生油层系的石油资源量; Sj - 生油岩的分布面积生油岩的分布面积; Hj - 生油岩的厚度生油岩的厚度; D - 生油岩密度生油岩密度; Aj - 生油岩氯仿沥青含量生油岩氯仿沥青含量 k1 - 排烃系数排烃系数; k2 - 聚集系数聚集系数35 凹陷石油资源量为三套生油层系资源量的和。三凹陷石油资源量为三套生油层系资源量的和。三套生油层系的地质参数见表套生油层系的地质参数见表8-2。 层系层系 地质参数地质参数层系层系1层系层系2层系层系3生油岩分布

34、面积生油岩分布面积1400070003000岩石密度岩石密度232323排烃系数排烃系数0.440.480.43聚集系数聚集系数0.1110.1110.111生油岩厚度生油岩厚度数据个数数据个数1407030取值范围取值范围0.1-1.00.1-1.00.1-0.5氯仿沥青含量氯仿沥青含量数据个数数据个数373721取值范围取值范围0.03-1.740.02-2.080.03-1.70表表8-2 三套生油层系的地质参数三套生油层系的地质参数 36第一套层系氯仿沥第一套层系氯仿沥青分布密度青分布密度第一套层系氯仿第一套层系氯仿沥青分布函数沥青分布函数第一套层系厚第一套层系厚度分布函数度分布函数第

35、一套层系厚第一套层系厚度分布密度度分布密度第一套层系石油资第一套层系石油资源量分布函数源量分布函数0.00.20.40.60.81.0F(a)a0.061.970.00.20.40.60.81.0f(a)a0.061.970.00.20.40.60.81.0F(h)h0.101.000.00.20.40.60.81.0f(h)h0.101.000.00.20.40.60.81.0F(q)q2.2458.0396.01037第二套层系氯仿第二套层系氯仿沥青分布函数沥青分布函数第二套层系氯仿沥第二套层系氯仿沥青分布密度青分布密度第二套层系厚第二套层系厚度分布函数度分布函数第二套层系厚第二套层系厚度

36、分布密度度分布密度第二套层系石油资第二套层系石油资源量分布函数源量分布函数0.00.20.40.60.81.0F(a)a0.051.540.00.20.40.60.81.0f(a)a0.051.540.00.20.40.60.81.0f(h)h0.101.000.00.20.40.60.81.0F(h)h0.101.000.00.20.40.60.81.0F(q)q3.8169.4897.48338第三套层系氯仿第三套层系氯仿沥青分布函数沥青分布函数第三套层系氯仿沥第三套层系氯仿沥青分布密度青分布密度第三套层系厚第三套层系厚度分布函数度分布函数第三套层系厚第三套层系厚度分布密度度分布密度第三套层系石油资第三套层系石油资源量分布函数源量分布函数0.00.20.40.60.81.0F(a)a0.021.740.00.20.40.60.81.0f(a)a0.021.740.00.20.40.60.81.0f(h)h0.101.000.00.20.40.60.81.0F(h)h0.101.000.00.20.40.60.81.0F(q)q1.9538.6546.75839第一套层系石油资第一套层系石油资源量分布函数源量分布函数0.00.20.40.60.81.0F(q)q2.2458.0396.010第二套层系石油资第二套层系石油资源量分布函数源量分布函数0.00.20.40.60