第8章蒙特卡罗模拟.pptVIP

第八章 蒙特卡罗模拟法;§1 蒙特卡罗法概述; 如何利用蒙特卡罗法进行油气资源量的计算？在概念模型清楚的基础上，首先要构造表示所研究问题概率解的数学模型(计算模型)。;(8-1);一、蒙特卡罗法的基本思想;二、蒙特卡罗法求解的基本过程;图8-1 随机变量Y的分布曲线;§2 随机数及其产生方法; 尽管如此，只要对伪随机数序列进行一系列严格的统计检验，证明它可以满足统计的要求，则伪随机数就可以做为真随机数使用。 ; xn , xn+1—第n次和第n+1次产生的伪随机数； α-乘子系数；M- 模； rn+1- [0,1]区间上的伪随机数。 xn+1≡αxn (mod M)叫做以 M 为模的同余式，表示xn+1取值为：α与xn的积除以 M 的余数部分。;2. 混合同余法; 一、随机变量的经验分布函数; 经验分布函数是由X的 n 个观测值x1, x2,…, xn ,用统计方法得到的分布函数，记为Fn(x) 。;(8-4); ② 计算区间端值 将(xmin , xmax) k等分, 各小区间的k+1个端值为：;AF(x);图8-5 经验分布函数曲线; 将坐标原点设为(xmin ,0)，若已知[0,1]区间上均匀分布的随机数ri ，则在图8-6纵轴上可确定点(xmin , ri ), 过该点作横轴的平行线交分布曲线于点(xi , ri ), xi 则是对应于随机数 ri的一次随机抽样值。; 如前所述,局部含油气地质单元是估算资源量的基本地质体。对第j个局部地质单元来说，其油气资源量的通式为:; 二、资源量的计算 1. 最大最小可能资源量及累计频率区间间隔值 为了求资源量的分布曲线，在资源量计算之前，先求出资源量的最大、最小可能值和累积频率小区间的端值:; ①对n个随机变量Xji (i=1,2,…,n)各进行一次抽样 ,得第一次抽样值xji1 (i=1,2,…,n),由第一次抽样值的积得到资源量的第1个估计值，即：;按上述做法，可得:;图8-8 油气资源量的概率加示意图;四、估算资源量的基本过程;······; 由于石油勘探的未来成效具有不确定性，因而需要对估算的石油资源量进行风险分析。石油勘探中的风险是多种多样的，如勘探区是否具备形成油气藏的地质条件的风险；在具备形成油气藏地质条件的含油气区内，能否找到一定规模油气藏的勘探风险；勘探后发现的油气藏是否具备开采价值的经济风险；勘探过程中人与设备是否安全的风险；对于勘探地区，特别是大陆架地区是否有国际争议的政治风险等。; 油气资源评价为油气勘探提供决策依据，而油气资源估算是资源评价的组成部分。因此，要对上述风险下估算的油气资源量做好地质风险分析。; 式中K是圈闭的风险值, ki是第i个因素的风险值。例如，若用容积法估算圈闭的石油储量，则有：; S(含油面积)的风险决定于地质调查或地震勘探资料的可靠性；;地 质 参 数 ;图8-9 风险分析前后石油储量关系; 若第j个局部单元油气资源量分布函数为AF(qj)，地质风险值为kj ，保险值为1-kj，则需重新估算风险分析后的油气资源量。既在原分布函数基础上再次抽样，入口值区间由[0,1]变为[0 ,1-kj]，如图。; 将[0,1]上分布的随机数rs改造为[0,1-kj]上分布的随机数rs*的变换公式为rs*=(1-kj )rs。然后进行随机抽样，求出风险分析后局部单元油气资源量的值:;§6 应用简例;;第一套层系氯仿沥青分布密度;第二套层系氯仿沥青分布函数;第三套层系氯仿沥青分布函数;第一套层系石油资源量分布函数; 概率/%; 本章复习要点 1. 蒙特卡罗法的概念及概率解的表达形式; 2. 形成[0,1]区间上伪随机数的两种方法 ; 3. 随机变量经验分布函数的形成; 4. 随机变量的抽样法种类及经验函数抽样法; 5. 地质风险分析及考虑风险后的经验分布函数; 6. 估算一个地区油气资源总量的一般步骤; 7. 实际应用。;思考与练习题 1. 简述蒙特卡罗法的基本思想及过程。 2. 对随机变量进行随机抽样时,为何要产生[0,1]上的随机数? 3. 如何构造地质研究中随机变量的经验分布函数?它与概率论中随机变量的分布函数有何不同? 4. 估算油气资源量时，如何实现对随机变量的随机抽样?常用的抽样方法有几种? 5.蒙特卡罗法估算油气资源量比传统估算方法有何优点? 6. 应用蒙特卡罗法求数值解不确定问题概率的基本过程是什么?;END