第十章蒙特卡罗模拟介绍.ppt

若第j个局部单元油气资源量分布函数为AF qj ，地质风险值为kj ，保险值为1-kj，则需重新估算风险分析后的油气资源量。既在原分布函数基础上再次抽样，入口值区间由[0,1]变为[0 ,1-kj]，如图。 图8-10 风险分析前后石油储量关系 二、风险分析后资源总量的估算 风险分析前 风险分析后 0.5 1.0 0.0 F q q 1-kj 将[0,1]上分布的随机数rs改造为[0,1-kj]上分布的随机数rs*的变换公式为rs* 1-kj rs。然后进行随机抽样，求出风险分析后局部单元油气资源量的值: 求得qjα α 1,2,…, N 后，可用频率统计法求出各局部单元的资源量分布函数AF qj j 1,2,…, m ，最后再利用“抽样加”和“频率统计法”求出地质风险分析后资源总量的分布函数AF* q 。 AF* q AF q2 +AF q2 +…+AF qm qjα α 1,2,…, N §6 应用简例 某沉积盆地中的一个地质凹陷有三套生油层系。用氯仿沥青法估算该凹陷的远景石油资源量，各层系的石油资源量估算公式如下： Qj - 第j j 1,2,3 套生油层系的石油资源量; Sj - 生油岩的分布面积; Hj - 生油岩的厚度; D - 生油岩密度; Aj - 生油岩氯仿沥青含量 k1 - 排烃系数; k2 - 聚集系数 凹陷石油资源量为三套生油层系资源量的和。三套生油层系的地质参数见表8-2。 0.03-1.70 0.02-2.08 0.03-1.74 取值范围 21 37 37 数据个数 氯仿沥青含量 0.1-0.5 0.1-1.0 0.1-1.0 取值范围 30 70 140 数据个数 生油岩厚度 0.111 0.111 0.111 聚集系数 0.43 0.48 0.44 排烃系数 23 23 23 岩石密度 3000 7000 14000 生油岩分布面积 层系3 层系2 层系1 层系 地质参数 表8-2 三套生油层系的地质参数 第一套层系氯仿沥青分布密度 第一套层系氯仿沥青分布函数 第一套层系厚度分布函数 第一套层系厚度分布密度 第一套层系石油资源量分布函数 0.0 0.2 0.4 0.6 0.8 1.0 F a a 0.06 1.97 0.0 0.2 0.4 0.6 0.8 1.0 f a a 0.06 1.97 0.0 0.2 0.4 0.6 0.8 1.0 F h h 0.10 1.00 0.0 0.2 0.4 0.6 0.8 1.0 f h h 0.10 1.00 0.0 0.2 0.4 0.6 0.8 1.0 F q q 2.245 8.039 6.010 第八章 蒙特卡罗模拟法 §1 蒙特卡罗法概述 §2 随机数及其产生方法 §3 随机变量的抽样 §4 油气资源量的估算 §5 地质风险分析 §6 应用简例 §1 蒙特卡罗法概述 蒙特卡罗法 Monte Carle :以数值解不确定问题为对象，对计算模型中的各变量进行随机抽样 随机试验 ，进而求问题概率解的一种统计学方法。因此，蒙特卡罗法又称为统计试验法。 蒙特卡罗法是以概率论与数理统计理论为指导的、有着广泛应用领域的通用性统计学方法，其核心是对随机变量的抽样。常用于油气资源量预测。 美国第二次全美石油资源评价的主要方法就是蒙特卡罗法。目前各主要产油国及各大石油公司都把该法作为油气资源评价的重要方法之一。 如何利用蒙特卡罗法进行油气资源量的计算？在概念模型清楚的基础上，首先要构造表示所研究问题概率解的数学模型 计算模型 。 我国1979年把蒙特卡罗法用于油气资源评价,该方法目前仍是我国预测油气资源量的常用方法之一,也是以统计预测为主的应用软件评价系统的核心算法之一。 对一个实际地区的油气资源总量而言，它是各局部含油气地质单元油气资源量之和。而局部含油气地质单元既可以是生油凹陷中的一个生油层系，也可以是次一级构造单元中的生油层系，还可以是局部构造等。 8-1 油气总资源量： Qj -第j个局部地质单元油气资源量 j 1,2,…, m ; Xjl -第j个局部地质单元第l个地质变量; cj -第j个局部地质单元中v个地质常数Di之积。 虽然“局部地质单元”含义不同，但任何一个局部含油气地质单元的油气资源量都可归纳为与油气资源量相关的地质常数和变量的乘积，即: 因此，求油气资源量的问题，就归结为求上述两个计算模型的概率解问题。 8-2 一、蒙特卡罗法的基本思想 蒙特卡罗的基本思想可概括为：为求研究问题的概率解, 构造一个表示所研究问题概率解的数学模型 计算模型 ，记为: 8-3 依据计算模型中各随机变量Xi所服从的分布进行随机抽样，并按计算模型计算Y的多个估计值，最终用频率