第八章蒙特卡罗模拟ppt课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,局部含油气地质单元是估算油气资源量的基本地质体,。它的含义可以不同，既可以是生油凹陷中的一个生油层系，又可以是次一级构造单元中的生油层系，还可以是局部构造等。关于资源量估算方法，基于不同的找油理论，估算方法也不一样，但是,任何一个局部含油气地质单元的油气资源量却都是与油气资源量相关的地质常数和变量的连乘积，即,第八章 蒙特卡罗模拟法,一、油气资源量估算通式及其关键问题,1.,油气资源量估算通式,(8-1),局部含油气地质单元是估算油气资源量的基本地质,1,2. 关键问题,油气资源量是个随机变量，对其预测的实质是估计它的数学期望，也就是与油气资源量相关的地质随机变量连乘积的数学期望。由式,(8-1),可知，油气资源量的数,学期望和方差为：,式中：,Q,局部地质单元油气资源量；,k ,与油气资源量有关的常数或经验系数的积 ；,与油气资源量有关的地质变量的积。,（8-2）,由上式可得以下认识：,估算局部含油气地质单元油气资源量的,实质是求,2. 关键问题式中： Q局部地质单元油,2,油气资源量的数学期望。,它的可靠性依赖于参与估算的随机变量。,若,k,中含有面,积,(常数),参,数,，,那么,预测资源量的关键则是求单位面积的资源量,。局部含油气地质单元油气资源量的期望值与单位面积资源量、,k,成正比，而方差与单位面积的资源量、,成正比。由此表明，,单位面积资源量极大地影响着资源量的期望值,。但是，它却是一个难以准确估算的问题。当单位面积资源量的偏差和,k,的数量级较大时，预测结果将明显偏离油气资源量的期望值，并导致资源量估计值的分布区间急剧变大。,由上可知，,油气资源量估算不仅是个典型的数值解不,油气资源量的数学期望。它的可靠性依赖于参与估算的随机变量。,3,确定问题，而且也不能考查传统方法每次估算结果的可靠程度。为此，在油气资源量估算中引入了蒙特卡罗法，它,是对数值解不确定问题求概率解的一种统计试验法。利用这种方法不仅可以求出问题的近似解，而且还可以给出与之对应的概率。,二、特卡罗法的基本思想及其求解的基本过程,1. 蒙特卡罗法的基本思想,蒙特卡罗法又称统计试验法，其基本思想可概括为：,欲求研究问题的概率解，就构造一个表征所研究问题概率,解的数学模型，记为：,（8-3）,确定问题，而且也不能考查传统方法每次估算结果的可靠程度。为此,4,使所求的概率解是式（8 -3）的某个数字特征（如数学期望、方差等），而且这个数字特征又可用统计的方法求得其估计值，那么就把该估计值作为研究问题概率解。,2. 求解的基本过程,由方法的基本思想可知，求解过程大致分为四步：,(1)分,析并拟定给定问题中的随机变量，构造表征,给定问题概率解的数学模型。,(2)对模型中的随机变量,进行,m,次随机抽样，获得 m,组抽样(观测)值：,(,k=,1,2,m,),使所求的概率解是式（8 -3）的某个数字特征（如数学期望、方,5,(3),把,m,组抽样值代入概率解模型，求出随机变量,Y,的,m,个估计值,y,1,y,2,y,m,；,y,F,(,y,),1,0.5,0,P=0.5,图8-1 随机变量,Y,的分布曲线,(4),利用频率统计法，由,Y,的估计值,求出,描述,Y,分布特征的分布曲线，如图,8-1,。,04,(3)把 m 组抽样值代入概率解模型,6,由上可知,蒙特卡罗法并不是估算油气资源量的专用方法，而是以概率论与数理统计理论为指导的有着广泛应用领域的通用性统计学方法，,其核心是对随机变量的抽样。,美国于1975年完成的第二次全美石油资源评价的主要方法就是蒙特卡罗法。目前各主要产油国及各大石油公司都把该法作为油气资源评价的重要方法之一。,1979年我国把蒙特卡罗法用于油气资源评价，是我国以统计预测为主的应用软件评价系统的核心算法之一。,由上可知,蒙特卡罗法并不是估算油气资源量的专用方法，,7,随机数是随机变量的观测值，由其构成数列,x,1,x,2,x,n,是一个无周期的数据序列，叫做随机数序列,。我们不可能对随机变量进行无限次的观测获得随机数序列。目前多用数学的方法在计算机上产生用于随机抽样的随机数序列,，它是一个周期性的有限序列，称其为伪随机数序列，其中的元素叫做伪随机数。当伪随机数序列的周期足够长,时，就可以满足对地质随机变量随机抽样,的需要。,二、随机数及其产生方法,1.随机数,2. 伪随机数的产生方法,（1）乘同余法,随机数是随机变量的观测值，由其构成数列 x1,8,式中,x,n,x,n+,1,第,n,次和第,n,+1,次产生的伪随机数；,乘子系数；,M,模；,r,n,0,1,上的第,n,个伪随机数。,x,n,+1,x,n,(,mod,M,)叫做以,M,为模的同余式，表示,x,n,+1,为,与,x,n,的积除以,M,的余数部分。,该方法产生伪随机数序列的递推同余公式为：,用乘同余法产生伪随机数时，可取初值,而,(,l,=1,2, )的奇数，,k,为计算机字长。,式中 x n , xn+1第n次和第n+1次产生的伪,9,三、随机变量的抽样,(2)混合同余法,对随机变量的抽样有经验分布函数抽样法、直接抽,混合同余法产生伪随机数序列的递推同余式如下：,混合同余法比乘同余法仅是增加了一个增量,，其它含义与乘同余法相同。,在,M,=,2,19,=,524288,=,5,5,=,3125,时，,x,0,=,23,11,19,37；,=,3,7,11,17,分,4,套配合使用，混合同余法产生周期为,524288,伪随机数序列。,三、随机变量的抽样 (2)混合同余法,10,随机变量,X,的分布函数是随机变量,X,的取值不大于实数,x,的概率，通常记为：,它的图形如图,8-2。,样法和变换抽样法。在此介绍经验分布函数抽样法。,1. 随机变量的经验分布函数,1,图8-2 随机变量分布曲线,随机变量X的分布函数是随机变量X的取值不,11,经验分布函数是由,X,的,n,个观测值,，,用统计方法得到的分布函数，并记为,。,在油气资源评价中，人们总希望得到资源量不小于某个实数,x,的概率,P,(,X,x,),显然:,它的图形如图8-3。,),(,1,),(,1,),(,x,F,x,X,P,x,X,P,-,=,-,=,P,(,Xx,),图,8-3,分布函数曲线,0,1,经验分布函数是由X 的 n 个观测值,12, 计算区间间隔值,计算,m,+1个区间间隔值：,（8-4）,用,代替,F,(,x,)并记为：,上式是油气资源评价中的经验分布函数。,2. 经验分布函数的构造方法,用频率统计法构造经验分布函数的条件及步骤如下：,(1),使用条件：,随机变量观测值的个数不少于,30。,(2),构造分布函数的步骤, 确定频率统计区间数,设有,n,个观测值，区间数为,m,，按照,n,/,m,等于35的原则确定,m,的值，且,m,为奇数。, 计算区间间隔值（8-4）用,13, 经验分布函数,记,n,i,为观测值落入区间（,x,i,，,x,i+,1,）内的频数，,f,i,为累加频率，则：,经验分布函数是：,式中：,分别是观测值的最大和最小值。, 经验分布函数,14,坐标原点为,(,x,min,0），,r,i,是随机数序列中的第,i,个随机数，在图8 -4纵轴上确定点,(,x,min,r,i,),，过该点作横轴的平行线交分布曲线于点(,x,i,r,i,)，,x,i,则是对应于随机数,r,i,的一次随机抽样值。用上述方法可得到随机变量的一系列随机抽样值 ，称该抽样方法称为随机数技术。,x,min,0,x,i,(出口值),(,x,i,r,i,),x,min,r,i,(入口值),x,A,F,(,x,),图8-4 分布函数抽样过程示意图,3. 经验分布函数的抽样,坐标原点为(x min ,0），r i 是,15,如前所述，局部含油气地质单元是估算资源量的基本地质体。不同估算方法，其含义有所不同，例如：,容积法，地质体一个油藏或一个油层；单储系数法，地质体局部构造断块等；,氯仿沥青法：地质体生油凹陷等。,但对第,j,个局部地质单元来说，其油气资源量的通式为：,四、局部含油气地质单元油气资源量的估算,1. 资源量概率模型,如前所述，局部含油气地质单元是估算资源量的,16,式中:,Q,j, 第,j,个局部地质单元的油气资源量；,X,j,i,第,j,个局部地质单元中第,i,个地质随机变量；,k,j,第,j,个局部地质单元中,l,个地质常数与经验系数的积。,2. 资源量的计算,1) 最大最小可能资源量及累积频率区间间隔值,为了求资源量的分布曲线，在资源量计算之前，先求出资源量的最大、最小可能值和累积频率区间间隔值：,最大可能值,式中: Q j 第 j 个局部地质单元的油气资源量；,17,2) 资源量的计算, 对随机变量,X,j i,(,i,=1,2,n,)各进行一次抽样，得第一次抽样值,x,ji,1,(,i,=1,2,)，,由第一次抽样值的积得到资源量的第一个估计值，即：,最小可能值,区间间隔 ,式中 是区间间隔个数,2) 资源量的计算 最小可,18,对随机变量,X,j i,(,i,=1,2,n,) 再进行一次抽样，得第二次抽样值,x,ji,2,(,i,=1,2,n,)，由第二次抽样值的积得到资源量的第二个估计值，即：,重复上述做法，可以得到 :,对随机变量 X j i ( i =1,19, 第,j,个局部地质单元的资源量的分布函数,根据,1,),中给出的区间间隔值，用频率统计的方法求出,Q,j,的分布函数,AF,(,Q,j,)，它的图形如图8-5所示。,Q,j,AF,(,Q,j,),Q,j,min,图 8-5 资 源 量 分 布 图, 第 j 个局部地质单元的资源量的分布函,20,图8-6 油气资源量的概率和示意图,Q,1,Q,2,Q,m,Q,=,+,+,+,五、总资源量,这是个局部地质单元资源量的概率加（对各资源量分布函数的随机抽样加）。,资源量概率加的大致过程如图8-6所示：,图8-6 油气资源量的概率和示意图Q1Q2Q mQ=+,21,六、估算资源量的基本过程,对于一个局部地质单元：,选择预测(估算)方法(建立概率估算模型)；,确定参数中的随机变量；,构造随机变量的分布函数；,对各随机变量的分布函数进行重复抽样，计算资源量的估计值；,由资源量的估计值求资源量分布函数。,对于,m(m1),个局部地质单元：,还需要求总的资源量(局部地质单元资源量分布函数是求总资源量的基础）,六、估算资源量的基本过程对于一个局部地质单,22,图8-6 油气资源量的概率和示意图,Q,1,Q,2,Q,m,Q,=,+,+,+,图8-6 油气资源量的概率和示意图Q1Q2Q mQ=+,23,由于石油勘探的未来成效具有不确定性，因而需要对估算的石油资源量进行风险分析。石油勘探中的风险是多种多样的，如勘探区是否具备,形成油气藏的地质条件的风险,；在具备形成油气藏地质条件的含油气区内，能否,找到一定规模的油气藏的勘探风险,；勘探后发现的油气藏是否具备,开采价值的经济风险,；勘探过程中人与设备是否安全的,环境风险,；对于勘探地区，特别是大陆架地区是否有国际争议的,政治风险,等。,七、地质风险分析,由于石油勘探的未来成效具有不确定性，因而,24,油气资源评价为油气勘探提供决策依据，而油气资源估算是资源评价的组成部分。因此，地质人员要对上述风险下估算的油气资源量做好,地质风险分析,。,地质风险分析是指对模型中的参数做进一步的深入研究，,对参数共同存在的可能性做出概率估计,，在此基础上对预测的资源量进行可靠性分析。,实际工作中，地质风险分析可以在不同的层次进行，,如,单一地质圈闭,的风险分析，,一组地质圈闭,(国外对地质条件相似的一组地质圈闭称作一个勘探层)的风险分析，,一个油气聚集带,的风险分析，整个,含油气盆地,的风险分析等。,地质风险分析大多从圈闭做起，,其风险的计算公式为,油气资源评价为油气勘探提供决策依据，而油气资源估算是,25,式中,K,是圈闭的风险值， 是第,i,个因素的风险值。例如，若用容积法估算圈闭的石油储量，则有,石油储量,含油面积,储层厚度,储层孔隙度,石油充满系数,采收率,对储量进行风险分析，就要由熟悉含油气区地质资料的人员对上述5项参数逐个分析论证。通常来说：,式中K是圈闭的风险值， 是第i个因素的风险值。例如，若用容积,26,S,的风险决定于地质调查或地震勘探资料的可靠性；,的风险决定于储集层孔隙是否有次生改造或后期充填的影响；,H,的风险受岩性岩相变化的影响；,D,的风险受生油岩的成熟度和油气运移通道的制约;,W,的风险则与原油性质及驱动类型有关。,经,风险分析论证后,对每个参数给出风险值,。目前尚无完善的方法确定风险值，一种方法是由分析论,证人凭经验指定；另一种方法是借用地质条件类似的临区,的风险值。如，某个地质圈闭经风险分析论证后给出如下,风险值,S的风险决定于地质调查或地震勘探资料的可靠性；,27,地 质 参 数 保,风险值/,k,保险值/（1-,k,）,含油面积,储层厚度,储层孔隙度,充满系数,采收率,0,0.5,0,0.3,0,1,0.5,1,0.7,1,根据表内数据，计算该圈闭的风险系数为0.65，而保险值为0.35，即经过风险分析后该地质圈闭的石油储量仅是风险分析前的35%。,地 质 参 数 保风险值/k保险值/（1-k）含油面积01,28,风险分析前,风险分析后,图风险分析前后石油储量关系,风险分析前风险分析后图风险分析前后石油储量关系,29,若含油气区内第j个局部含油气地质单元油气资源量的分布函数为 ，地质风险值为 ，保险值为 ，那么在重新估算含油气区地质风险分析后的油气资源总量,时，随机抽样入口值的分布区间应由0,1变为 ，,风险分析前,风险分析后,图风险分析前后石油储量关系,八、风险分析后资源总量的估算,若含油气区内第j个局部含油气地质单元油气资源量的分布,30,因此，应将0,1上分布的随机数改造为0, 上分布的随机数 ，且 。在以上基础上，求出风险分析后局部地质单元油气资源量的估计值,再由 求得油气资源总量估计,值,，最后利用频率统计法求出地质风险分析后资源总量的分布函数,因此，应将0,1上分布的随机数改造为0, 上分,31,九、应用算例,某沉积盆地中的一个地质凹陷有三套生油层系。用氯仿沥青法估算该凹陷的远景石油资源量，各生油层系的石油资源量估算公式如下：,生油层系的石油资源量,生油岩的分布面积,生油岩的厚度,生油岩密度,生油岩氯仿沥青含量,排烃系数,聚集系数,九、应用算例 某沉积盆地中的一个地质凹陷有三套生油层系,32,凹陷石油资源量为三套生油层系资源量的和。三套生油层系的地质参数见表8-4。,表8-4 三套生油层系的地质参数,层系,层系1,层系2,层系3,生油岩分布面积,14000,7000,3000,岩石密度,23,23,23,排烃系数,0.44,0.48,0.43,聚集系数,0.111,0.111,0.111,生油岩厚度,数据个数,140,70,30,取值范围,0.1-1.0,0.1-1.0,0.1-0.5,氯仿沥青含量,数据个数,37,37,21,取值范围,0.03-1.74,0.02-2.08,0.03-1.70,地质参数,凹陷石油资源量为三套生油层系资源量的和。三,33,第一套生油岩氯仿沥青分,布密度,第一套生油岩氯仿沥青分布函数,第一套生油岩厚度分布函数,第一套生油岩厚度分布密度,第一套生油岩石油资源量分布函数,第一套生油岩氯仿沥青分第一套生油岩氯仿沥青分布函数第一套生油,34,第二套生油岩氯仿沥青分布函数,第二套生油岩氯仿沥青分布密度,第二套生油岩厚度分布函数,第二套生油岩厚度分布密度,第二套生油岩石油资源量分布函数,第二套生油岩氯仿沥青分布函数第二套生油岩氯仿沥青分布密度第二,35,第三套生油岩厚度分布函数,第三套生油岩厚度分布密度,第三套生油岩氯仿沥青分布函数,第三套生油岩氯仿沥青分布密度,第三套生油岩石油资源量分布函数,第三套生油岩厚度分布函数第三套生油岩厚度分布密度第三套生油岩,36,图8-7 凹陷生油岩生油参数及石油资源量分布函数,凹陷石油资源量分布函数,=,第一套生油岩石油资源量分布函数,+,第二套生油岩石油资源量分布函数,+,第三套生油岩石油资源量分布函数,图8-7 凹陷生油岩生油参数及石油资源量分布函数凹陷石油资源,37,