储层随机建模方法研究进展

1、第 31卷 第 5期物探化探计算技术 2009年 9月 基金项目 :中国石油天然气股份有限公司工程技术重点攻关项目 (2006-2008 收稿日期 :2008-08-02 改回日期 :2009-02-24文章编号 :1001 1749(2009 05 0454 06储层随机建模方法研究进展李 霞 1, 王铜山 1, 王建新2(1. 中国石油勘探开发研究院 , 北京 100083;2. 东方公司东部事业部 长庆经理部 , 宁夏 银川 751000摘 要 :油气储层随机建模方法是近年来新发展起来的一门油气预测技术 , 它在综合了地质 、 地球物理信息的基础上 , 充分利用已知点的储层信息 , 根据

2、所要预测的储层属性来选择不同的随机 模型 , 从而实现了对油气储层的定量表征及对各种尺度非均质性的定量刻画 , 其最终模拟结果可 以有任意多个以供油藏描述选择 。 这里简要介绍了目前比较常用的储层随机建模方法 、 原理 , 各 方法的优缺点 , 以及不同随机模型的地质适用性 , 并针对我国储层的特点 , 术的发展方向 。关键词 :储层随机建模 ; 随机模型 ; ; 中图分类号 :TE 13:A0 前言储层随机建模技术的产生 , 与地质统计学及油 气生产的需要密切相关 。地质统计学创建于二十 世纪六十年代初 , 由法国著名学者 G . Mather on 教 授提出 , 他将传统统计学理论与区域

3、化变量的概念 相结合 , 发展出一套以变差函数为工具 , 来研究矿 产矿化特征区域分布的数学技术 。在其后的发展 中 , 地质统计学先后出现了克里格技术 1、 2和随机模拟技术3。 在 1985年之前 , 随机模拟技术以非条件模拟为主 , 如转向带法 、 傅立叶谱估计法和 LU 矩阵分解法等 。 到九十年代 , 克里格技术不但被用 作插值方法 , 而且越来越多地用于建立数据的条件 累积分布函数 (CCDF , 使随机建模得到了飞速发 展 。 该技术逐渐被用来解决储层表征中的一些问 题 , 如建立储层物性非均质模型 , 储层内部非渗透 性隔夹层模型 , 以及储层空间连续性模型等 , 为油 藏早期

4、评价及开发阶段制定方案服务47。近几年来 , 人们应用计算机技术 , 发展出一套利用计算机存储和显示的三维储层建模方法 。 运用此方法 ,在油藏评价 , 乃至油田开发的不同阶段 , 均可建立 三维储层地质模型 。 目前 , 建立储层地质模型有确 定性建模和随机建模二种方法 , 而后者是近年来国 内 、 外学者研究的热点 。1 储层随机建模概述建立储层地质模型 , 是现代油藏描述的核心与 关键 。 建立储层地质模型主要有二种途径 , 即确定 性建模与随机建模 。确定性建模是以确定性资料 为基础 , 推测井间确定的 、 唯一的 、 真实的储层参 数 。 然而 , 地下储层是复杂的 , 它是许多复杂

5、地质 过程综合作用的结果 , 具有复杂的储层结构空间配 置及储层参数的空间变化 。以非常有限的资料描 述地下储层属性的分布是困难的 , 必然存在着不确 定性与随机性 。 因此 , 储层随机建模技术就应运而 生 , 并得到了越来越广泛的应用 。 随机建模是以已 知的信息为基础 , 以随机函数为理论 , 应用随机模 拟方法产生多个可选的 、 等概率的储层模型方法 。 虽然预测结果可能有任意多个 , 但是其结果可以真 实地反映储层属性的非均质性 。 当然 , 已知点的属 性是不变的 , 而未知点的属性会有不同的结果出 现 。 多种结果也表现出了储层属性的不确定性 , 其 中有 “ 悲观 ” 的结果

6、, 也有 “ 乐观 ” 的结果 , 可以给地 质人员提供不同参考 。 Haldorsen 8提出了应用随 机建模技术的理由 : 不同尺度下地质结构 , 岩石 物性的变化不确定 ; 相分布空间变化的复杂性 ; 储层动态和静态关系问题 ; 岩石特性和构造位 置空间不确定 ; 基于对建模速度和效率方面的考 虑 。 实践表明 , 储层随机建模在描述储层非均质 性 , 综合各类信息评价不确定性等方面 , 具有明显 优势 。2 储层随机建模技术2. 1 随机模型的分类随机模型的分类有很多种 , Haldorsen 等 8根 据研究现象的随机特征 ,型 、 连续模型和混合模型 。 9单元的特征 ,: 按数据

7、分布特征 :; 按模 拟的数据条件 :条件模拟和非条件模拟 ; 按使用 变量的个数 :单变量模拟和多变量协同模拟 。 目前 比较常用的方法是混合模拟 , 亦称二步模拟 , 即先 应用离散模型描述储层大规模非均质性 , (如沉积 相等 , 然后应用连续模型描述各个沉积相内部的 岩石物理参数的变化特征 , 这样会使预测结果更合 理 。2. 2 随机建模方法2. 2. 1 基于目标的随机建模方法(1 布尔模拟 。该方法是随机建模技术中最 简单的一种方法 , 是一个“ 逐步逼近过程 ” 。布尔 模拟认为在结构模拟中 , 应从具有成因意义的目标 创建 , 不应一个个节点来模拟 。该方法的基本思 路 ,

8、是根据一定的概率定律 , 按照空间中几何物体 的分布统计规律 , 产生这些物体中心点的空间分 布 , 并通过多个随机函数的联合分布 , 确定中心点 处的几何物体形状 、 大小和方向 。 布尔模拟目前主 要用于非条件模拟 , 用来建立离散模型 , 适用于建 立砂体格架模型或隔 (夹 层分布模型 , 如砂体格 架平面 、 剖面或者三维空间分布模型 。 以泊松点过 程为基础的模拟算法 , 适合模拟砂岩背景上存在小 尺度泥岩隔层这类现象 ; 以吉布斯点过程为基础的 模拟算法 , 适合模拟河道砂岩带内各河道砂体相互 镶嵌的现象 , 如模拟河流或河流 三角洲河道及相 关的沉积相带 。布尔模拟方法的优点在于

9、 : 很容易用于二维 和三维建模 ; 所用的参数较少 ; 计算速度快 。 其缺点在于 : 井数据对模拟结果的限制作用小 , 可能会导致人为低估砂体的连续性 ; 目标形状是 高度简化后的结果 , 往往与实际情况相差甚远 。 近 年来 , 为了使布尔方法适合于多井情况 (条件数据 较多 , 许多学者对布尔方法进行了改进 , 以忠实 于条件数据 10。(2 示性点过程模拟 。该模拟也称为标点法 , 其模型为布尔模型 , 即储层模型根据具有某些成因 意义的目标来建立 , 而不是暂时建立一个基本的节 点或象元 。 它主要模拟离散性质的地质特征 , 如沉 积相分布 、 岩相分布 、 。从地质统计学 角度来

10、讲 , 亦即将目标体投放于背景相 , 这种方法适合于具有背景相的目标 (物 体或相 模拟 , 如冲积体系的河道和决口扇 (其背 景相为泛滥平原 , 三角洲分流河道和河口坝 (其 背景相为河道间和湖相泥岩 , 浊积扇中的浊积水 道 (其背景相为深水泥岩 , 滨浅海障壁砂坝 , 潮汐 水道 (其背景相为泻湖或浅海泥岩 等 。另外 , 砂 体中的非渗透泥岩夹层 、 非渗透胶结带 、 断层 、 裂缝 等 , 均可利用此方法来模拟 。示性点过程法的优点是运算速度快 , 方法简 单 , 模拟结果服从几何体形状的空间分布 。但是 , 在其应用中 , 要有很强的先验地质知识 , 如各相的 体积含量 , 各相几

11、何形态 (长 、 宽 、 厚等 。该方法 主要适用于油田勘探阶段 , 以及开发早期井间砂体 和非渗透隔夹层的描述 。 当钻井增加时 , 示性点过 程法存在的一个很大问题 , 是难以条件化到已知的 实验数据 11。 另外 , 由于没有考虑沉积过程中 , 先 沉积的砂体对后面的沉积有重要的控制作用 , 标点 过程法模拟的结果往往过于乐观 , 砂体的连通性过 好 12。2. 2. 2 基于象元的随机建模方法(1 序贯高斯模拟 。该模拟可用于模拟连续 的地质现象 , 如孔隙度 、 地震反射界面的构造图等 。 其模拟过程是从一个象元到另一个象元序贯进行 的 。 在序贯高斯模拟中 , 条件数据首先要转换为

12、标 准高斯值 , 以确定转换后数据的变差函数 。 用简单 克里格得到条件分布的估计值 。 在每个节点处 , 根 5545期 李霞等 :储层随机建模方法研究进展 据模拟值和条件数据得到的克里格估计值 , 以及相 关的克里格方差 , 被用来确定条件高斯分布 ; 然后 从分布中随机抽取 , 得到一系列实现 ; 最后 , 将高斯 模拟结果反转化为原始数据空间 。序贯高斯模拟的主要优点在于 : 数据的条件 化是模拟的一个整体部份 , 无需作为一个单独的步 骤进行处理 ; 可自动地处理各向异性问题 ; 适 合于任意类型的协方差函数 ; 在运行过程中 , 只 需要一个有效的克里格算法 。主要缺点是要求变 量

13、分布服从高斯分布 。(2 序贯指示模拟 。该模拟既可用于模拟连 续的变量 , 又可用于模拟离散的特征 9、 13。它不要 求对条件分布的参数形式作任何假设 , 在现有数据 和其它相关资料的基础上 , 根据门槛值将条件数据 转化为指示数据 。 根据相应的指示变差函数 , 应用 简单克里格给出每个网格节点处条件分布的估计 量 。 该方法多用于模拟复杂的岩相或岩相组合 、 沉 积相的空间分布结构 , 如模拟复杂的河 湖相沉 积 、 河流 三角洲沉积 、成因的碳酸盐储层的指示变差函数 ,的储层空间分布 。序贯指示模拟最大的优点 11, 是可以模拟复 杂各向异性的地质现象及连续性分布的极值 。对 具有不

14、同连续性分布的变量 (如沉积微相 , 可利 用不同的变差函数进行表征 , 从而建立各向异性 的模拟结果 。 另外 , 指示模拟除忠实硬数据 (如井 数据 外 , 还可方便地把软数据 (如地震 、 试井数 据 通过编码加入到模拟中 。其缺点是 : 模拟结 果有时并不能很好地再现输入的变差函数 ; 在条 件数据点较少且模拟目标各向异性较强时 , 难以计 算各类型变量的变差函数 ; 不能很好地恢复指定 模拟目标的几何形态 (尤其是相边界 , 一些类型 变量是以一个或多个象元为单元零星地分布 。 (3 截断高斯模拟 。该模拟属于离散随机模 型 , 用于分析离散型或类型变量 。 模拟过程是通过 一系列门

15、槛值及截断规则网格中的三维连续变量 , 建立离散物体的三维分布 14。由于离散物体的分 布取决于一系列门槛值对连续变量的截断 , 因此 , 模拟实现中的相分布是排序的 。这一方法适合于 相带呈排序分布的沉积相模拟 , 如三角洲 (平原 、 前缘和前三角洲 , 呈同心分布的湖相 (滨湖 、 浅 湖 、 深湖 , 以及滨面相 (上滨 、 中滨 、 下滨 的随机 模拟 。 该方法的优点是 : 易于实现 、 速度快 ; 可 在模拟中考虑地质因素 ; 可以对模拟结果进行条 件限制 , 使之与条件数据相吻合 。 缺点是多次截断 不能控制具体的指示协方差 , 特别是当相与相之间 各项异性彼此互不相同时 。胡

16、雪涛等 15通过比较 截断高斯模拟和序贯高斯模拟方法在沉积微相 (共八种 上的模拟 , 认为截断高斯模拟不适用于 非均质性复杂的低渗储层微相展布的描述 。这说 明该方法并非一种普适性方法 , 对相序关系非常简 单微相展布 , 可能得到较好的模拟结果 ; 但对微相 类型较多 , 展布关系较复杂的微相分布 , 模拟结果 的可靠性较低 。而序贯指示模拟适用于非均质性 复杂的低渗储层微相展布的描述 , 结果具有可信性 和预测性 。(4 马尔柯夫随机场模拟 。该模拟既可用于 离散物体 , 也可用于离散化连续变量类别的随机模 拟 。 , 某类型变量条件概率 , 条件 。 (如 Metr opolis -H

17、astings , , 然后 逐步进行迭代 , 直到满足指定的条件概率分布为 止 。 九十年代初又提出了用于相和岩性模拟的半 马尔柯夫随机域 (Tjel m eland, and Holen, 1993 模 型 。 在该模型中 , 象元岩性不是取决于局部邻区 , 而是取决于较大的区域 。输入参数包括类型变量 的分布范围 、 边界关系 、 各类型变量的含量等 。 马尔柯夫随机域及半马尔柯夫随机域模拟法 的优点 , 在于再现每一种状态复杂的非均质性能力 较强 , 适合于镶嵌状分布的相 (或岩性 的随机模 拟 , 以及单一类型的相或岩性分布 (如砂体内钙质 胶结层的分布 。 其缺点是 : 条件概率的

18、确定相 当复杂 , 特别是在条件数据有限时更为困难 ; 难 以很好地恢复相的几何形态 ; 难以应用软数据 (虽然很容易忠实硬数据 ; 模拟收敛很慢 。目 前这类模型应用很少 , 且主要限于二维空间 。 (5 二点直方图法的随机模拟 。该模拟主要 用于类型变量的随机模拟 , 它属于二点统计学的范 畴 。 其主要特征是在空间范围内 , 二个相距一定距 离的象元 , 分属于不同类型变量的转换概率分布 , 在特定偏移距所有二元类型变量的转移概率 , 即构 成二点直方图模型 9, 17, 它主要应用优化算法 (如 模拟退火 进行随机模拟 。由于转换概率的计算 , 在条件数据有限时很难进行 , 所以在实际

19、应用中 , 要有一个与待模拟地区地质条件相近的 , 数据密度 较高的原型模型 。二点直方图适用于镶嵌状分布654物探化探计算技术 31卷 的沉积相 (或岩性 的随机模拟 , 也可用于只有二 个相的沉积相随机模拟 。 在实际应用中 , 二点直方 图常应用于模拟退火中作为其它随机实现的后处 理 , 但也不能很好地恢复相几何形态 。目前 , 国内 的相关研究较少 。(6 模拟退火法 。该法适于模拟连续或离散 的特征 , 它的算法来自于金属淬火后冷却过程中分 子有序排列的过程 , 但是最低温度的求取比较困 难 , 在温度参数确定方面有不少人做了探索 1820, 其储层模型是通过迭代试算法来建立的 。一

20、般要 构建一个能量函数 , 使其为网格的统计量和目标值 之间的差 , 模拟的目的就是使能量函数 E 的值最 小 。 模拟退火法的优点是 : 综合利用各种来源信 息能量强 , 如野外露头 、 岩心 、 测井 、 地震 , 以及井间 示踪剂试井资料等 ; 模拟重现实验样品非均质性 的效果好 ; 能够灵活的忠实于所给出的任何地质 参数 。 缺点是计算量大 。该方法可以很好地保持 数据所反映的空间结构 , 即实验变差函数 。 在复杂 的碳酸盐岩条件下也可使用 ,来源的数据 21。(7 分形模拟的相似性 ,性 , 即在任一规模上 ,量的方差成正比 , 其比率取决于分形维数 (或间断 指数 。 He we

21、tt 22将分形随机域引入岩石物理参 数的随机模拟 , 其基本理论为分数布朗运动和分数 高斯噪音 。 Painter 23在分形随机模拟中引入了 Levy 稳态概率分布 , 避免了高斯分布的假设 , 可适 用于成层性很强的地层条件下随机变量的分形模 拟 。 分形模拟一般采用误差模拟算法 , 其模拟实现 为克里格估值加上随机“ 噪音 ” , 分形随机域的自 相似性是它最大的优点 。在确定变量符合分形特 征后 , 便可根据自相似性原理 , 应用少量数据 , 预测 整个模拟目标区的变量分布 。 然而 , 在分形模拟的 应用中 , 一定要检验待模拟变量是否具有分形特 征 。 由于地质情况的复杂性 ,

22、不同规模的地质特征 受控于不同的地质控制因素 , 而很多地质变量并不 一定符合分形特征 。另外 , 还要检验垂向与平面 上 , 分形特征的差别 。在很多分形模拟的应用中 , 由于数据点比较稀少 , 往往用垂向分形维数代替平 面分形维数 。 虽然很多学者证明 , 垂向与平面上分 形特征的相似性 , 但当模拟目标区纵 、 横向相变不 符合沃尔特相序时 , 垂向和平面上的分形特征便不 相似了 。 分形模拟可用于模拟连续的特征 , 也可用 于模拟天然裂缝的分布模式 , 多用于模拟孔 、 渗参 数的空间变化 。2. 3 随机建模方法及模型的选择在实际应用中 , 应从以下几方面考虑所选取的 建模方法是否合

23、理 : 应考虑变量的随机模型 , 不 同的随机模型有不同的模拟方法 ; 应考虑模拟的 对象 , 如对于沉积相等类型的变量和孔隙度等 , 具 有连续性的变量 , 模拟的侧重点不同 , 选取的模拟 方法也不同 ; 应考虑条件信息的种类 , 确定出是 来自单变量的信息 , 还是多变量的信息 ; 是软信息 , 还是这些信息的综合 。 不同的信息量 , 对应不同的 模拟方法 ; 要考虑区域地质情况 。 最后要选取不 同的优化算法 , 对模拟结果进行优化 。随机模型理论及应用研究表明 , 不同随机模型 对不同的地质条件 (如不同的沉积相 有一定的实 用性 。 A lbert 24曾讨论过多种随机模型的优缺

24、点 , 。从随 , 可将不同随机模型 :(和岩性 随机建模的随机模 , 示性点过程 , 截断高斯 , 序贯指 示 , 马尔柯夫随机场和二点直方图法 。 在待模拟目 标区存在多种沉积相 (或岩性 的情况下 , 布尔模 拟用于模拟砂体在空间的形态 、 大小 、 位置和排列 方式 ; 示性点过程适用于具有背景相的沉积相 (或 岩性 的随机模拟 ; 截断高斯适用于具有排序规律 的沉积相 (或岩性 的随机模拟 ; 序贯指示模拟 , 马 尔柯夫随机场和二点直方图适用于具镶嵌结构的 沉积相 (或岩性 的随机模拟 。在待模拟目标区 , 只有二种相 (或岩性 的情况下 , 上述各种模型均 可使用 。 但是 ,

25、由于所有基于象元的随机模型均不 能很好地恢复相 (或岩性 的几何形态 (尤其是相 边界 , 因此在这种情况下 , 应尽量使用基于目标 的随机模型 (示性点过程 。(2 用于岩石物理参数随机建模的随机模型 主要有序贯高斯 , 序贯指示 , 分形随机域 , 模拟退火 和马尔柯夫随机域 。序贯高斯模拟适用于各向异 性不强的条件下连续变量的随机模拟 ; 序贯指示模 拟适用于复杂各向异性的 、 具奇异值分布的连续变 量的随机模拟 ; 分形随机域适用于在数据点很少 , 且随机变量具有统计自相似性条件下连续变量的 随机模拟 ; 模拟退火适于模拟连续或离散的特征 , 可以很好地保持数据所反映的空间结构 , 在

26、复杂的 碳酸盐岩条件下也可使用 , 特别适用于综合多种来 源的数据 ; 马尔柯夫随机域可用于复杂各向异性条 件下连续变量的随机模拟 , 但由于其统计推断和参 7545期 李霞等 :储层随机建模方法研究进展 数求取十分复杂 (要求有训练图像 , 因此目前应 用很少 。3 储层建模的原则我国含油气盆地类型多 , 储层以陆相碎屑岩及 海相碳酸盐岩为主 , 储层成因复杂 , 非均质性严重 。 如河流 、 三角洲及冲积扇等环境形成的储层 , 在纵 、 横向上相变快 , 不同规模的非均质性严重 。因此 , 对这类储层进行勘探与开发 , 将面临储层非均质性 的问题 。 为了建立尽量符合地质实际情况的储层 模

27、型 , 针对我国储层的特点 , 制定如下建模原则 。 3. 1 确定性建模与随机建模相结合的原则确定性建模是根据确定性资料 , 推测出井间确 定的 、 唯一的储层特征分布 。 而随机建模是对井间 未知区 , 应用随机模拟方法建立可选的 , 等概率的 储层地质模型 。 应用随机建模方法 , 可建立一簇等 概率的储层三维模型 , 因而可评价储层的不确定 性 , 进一步把握井间储层的变化 。 在实际建模的过 程中 ,3. 2 等时建模原则沉积地质体是在不同的时间段形成的 。为了 提高建模精度 , 在建模过程中 , 应进行等时地质约 束 , 即应用高分辨率层序地层学原理来确定等时界 面 , 并利用等时

28、界面将沉积体划分为若干等时层 。 在建模时 , 按层建模 , 然后再将其组合为统一的三 维沉积模型 。同时 , 针对不同的等时层 , 输入反映 各自地质特征的不同的建模参数 , 这样可使所建模 型能更客观地反映地质实际 。3. 3 相控储层建模原则相控建模 , 即首先建立沉积相 , 储层结构或流 动单元模型 , 然后根据不同沉积相 (砂体类型或流 动单元 的储层参数定量分布规律 , 分相 (砂体类 型或流动单元 进行井间插值或随机模拟 , 进而建 立储层参数分布模型 。4 储层随机建模技术发展方向不同随机模拟方法适用于不同类型的油气藏 以及不同的阶段 , 考虑到中国油气储层的特殊性以 及储层随

29、机建模方法本身的一些局限性 , 需要在引 进国外技术的基础上 , 结合我国油气储层的特点 , 形成一系列特有的储层随机建模技术 。国内以张 昌民教授 25为首的 3RG 研究小组 , 经过十多年的 研究 , 不断发展创新 , 形成了一系列具有鲜明特色 的储层随机建模技术 。 其主要研究成果有 : 变差 函数球状模型的自动拟合 ; 储层随机建模技术与 GI S 相结合 ; 采用混合建模方法建立微相模型 ; 多物源条件下储层随机建模的方法 ; 布尔模拟 的改进 ; 基于储层骨架的多点地质统计学随机建 模方法 。 储层随机建模虽然经过不断完善 , 不断 推陈出新 , 然而应该看到 , 地质现象是非常

30、复杂 的 , 储层随机建模在许多方面仍然需要努力 , 以 满足油田生产开发的需求 10。4. 1 目标之间相互影响与函数的确定目标之间相互影响具有成因上的意义 。对于 河道沉积 , 先期沉积的河道影响后期河道的分布 。 一般来说 , 先期沉积河道区由于地形的变化 , 往往 成为后期河道沉积时的背景相区 。而同期沉积河 道之间 , 与沉积时构造 、 地形 、 坡度 、 物源等密切相 关 ,化 ,够 , 其沉积具有成因 、 决口扇等沉积微 、 分布等 , 对河流系统的建模 , 应该考虑这 些因素 。 而基于目标的标点过程法 , 由于是分目标 的模拟 , 很难考虑到目标间这种成因关系 。 如何在 模

31、拟中考虑目标之间的相互影响 , 并用数学模型 来表征 , 应该是储层建模以后研究的一个重要方 向 。4. 2 获取反应实际储层结构特征的参数随机建模需要输入大量参数 :在基于目标的建 模中 , 需要输入目标的几何形态参数 , 以及目标相 互之间影响参数 ; 在二点统计学中 , 需要输入变差 函数特征参数 ; 在多点地质统计学中 , 需要输入反 映储层结构特征的训练图像 。 这些参数的获取 , 需 要地质学家仔细分析 , 其获取是相当的困难 。 在基 于目标的建模中 , 目标的几何形态参数化往往难以 进行 , 因为地下储层形态很难用简单的长宽高来表 述 , 甚至难以用数学函数来反映 。在二点统计

32、学 中 , 变差函数仅反映了二点相关的关系 , 其准确推 断和拟合是模拟的基础 。 在多点统计学中 , 训练图 像决定着多点概率的求取 , 不同的训练图像将会产 生不同的模拟结果 。如何获取反映地下储层结构 特征的训练图像 , 一直是研究的热点和难点 。 4. 3 先验概率的求取在储层随机建模中 , 先验概率的获得往往被忽 视 , 而先验概率也是影响模型不确定性的重要因 素 。 在勘探阶段 , 先验概率的提取相当困难 , 这是 因为此时数据量过少 , 不足以推断出目标参数准确 分布 , 导致模拟结果带有不确定性 。在开发阶段 ,854物探化探计算技术 31卷 5期 李霞等 : 储层随机建模方法

33、研究进展 459 数据量增加 ,使得先验概率的推断成为可能 。尽管 如此 ,这种先验概率的获得仍然依赖于人的主观思 考 ,因而其中必然存在不确定性 , 因此对先验概率 的求取 ,也应该引起相当重视 。 4. 4 随机模拟结果评价 对随机模拟结果的评价 ,一直受到储层建模家 的关注 。为了对随机模拟结果进行合理评价 ,不少 学者给出了不同的评判标准 。例如通过统计多个 等概率随机模拟结果中每个网格点内目标出现的 概率 ,来获取最优的模拟结果 ; 通过抽稀检验评价 模拟结果 ; 通过与实际储层展布模式对比评价模拟 结果等 。由于随机建模结果将会输入到数值模拟 中 ,因此利用数值模拟结果进行检验 ,

34、 将是对随机 模拟结果评价的一个方向 。 10 尹艳树 ,吴胜和 . 储层随机建模研究进展 J . 天然 气地球科学 , 2006, 17 ( 2 : 210. 11 李少华 ,张昌民 ,尹艳树 . 河流相储层随机建模的几 种方法 J . 西安 石油学 院学 报 (自 然 科 学 版 , 2003, 18 (5 : 10. 12 HENR I UEZ A , TYLER K J. , HURST A , et al . Q . si ulation modeling C . SPE16753, 1987: 591. m Charcterizat ion of fluvial sedim ent

35、ology for reservoir 13 JOURNEL A G Non - parametric estim ation of spatial . distributions J . M athematical Geology, 1983, ( 15 : 445. 14 MATHERON G, BEUCHER H , DE FOUQUET C, et al taic reservoirs J . SPE16753, 1987: 591. 5 结语 储层随机建模已成为油气勘探与开发的核心 技术 ,其发展前景十分广阔 。我国含油气盆地类型 多 ,储集层以陆相碎屑岩及海相碳酸盐岩为主 , 储

36、层成因复杂 ,各种尺度的非均质性严重 。如何针对 国内不同类型 ,不同勘探与开发阶段的陆相高度非 均质复杂油气藏 ,采用不同的建模方法实现对油气 藏的精细描述 , 是当前面临的一项艰巨任务 。因 此 ,研究和使用各种随机建模方法 , 既具有理论意 义 ,又具有巨大的潜在经济价值 。随着计算机计算 能力的提高 ,以及新理论 、 新方法的引入 ,特别是多 点地质统计学的出现 ,随机模拟技术将在各个领域 显示其广阔的应用前景 。 参考文献 : 1 王仁铎 ,胡光道 . 线性地质统计学 M . 武汉 : 武汉 15 胡雪涛 ,李允 . 储层沉积微相的随机模拟及其对比研 究 J . 西南石油学院学报 ,

37、 2000 , 22 (1 : 16. 16 BESAGE J E. Spatial interaction and statistical analysis ries B , 1974, 36: 192. 17 FARM ER C L. The generation of stochastic fields of 18 姚姚 . 地球物理非线性反演模拟退火法的改进 J . 地球物理学报 , 1995, 38 ( 5 : 643. 19 井西利 ,杨长春 ,李丽 . 地球物理问题的自适应模拟 退火方法求解 J . 地球物理学进展 , 2003, 18 (2 : 327. 20 夏红敏 ,黄旭日

38、 ,王尚旭 ,等 . 区域特性约束下的油藏 物性模 拟 J . 地 球 物 理 学 进 展 , 2005, 20 ( 3 : 769. 21 QUENES A. App lication of sim ulated annealing and other global op tim ization methods to reservoir descrip tion M . M yths and Realities SPE28415, 1994. . 地质学院 , 1984. 2 侯景儒 ,尹镇南 ,李维明 , 等 . 实用地质统计学 (空间 信息统计学 M . 北京 : 地质出版社 , 199

39、8. 3 JOURNEL A G, HU IJBEREGTS C J . M ining Geosta2 tisticsM . New York : Academ ic Press , 1978. 4 张一伟 ,熊琦华 ,王志章 ,等 . 陆相油藏描述 M . 北 京 : 石油工业出版社 , 1997. 版社 , 1996. 5 裘怿楠 ,陈子琪 . 油藏描述 M . 北京 : 石油工业出 6 王家华 ,张团峰 . 油气储层随机建模 M . 北京 : 石 油工业出版社 , 2001: 93. 7 贾爱林 ,肖敬修 . 油藏评价阶段建立地质模型的技术 与方法 M . 北京 : 石油工业出版社 ,

40、 2001. ing J . JPT , 1990: 404. 8 HALDORSEN H H , DAM SLET H E. Stochastic model2 9 DEUTSCH C V , JOURNEL A G GSL I : Geostatistical . B 作者简介 : 李霞 ( 1982 - , 女 , 博士 ,主要从事油气 储层地球物理综合研究工作 。 22 HEWW TT T A. Fractal distribution of reservoir hetero2 geneity and their influence on fluid transport J . SPE

41、15386, 1986: 1. 23 SCOTT PA I TER. Numerical method for conditional N ology, 1995, 30 ( 2 : 163. 24 ALBER T F , MODO T V. Stochastic models of reser2 voir heterogeneities: I pact on connectivity and average m permeabilities J . SPE24893, 1992: 355. M . 北京 : 石油工业出版社 , 2007. 25 李少华 , 尹艳树 , 张昌民 . 储层随机建

42、模系列技术 si ulation of levy random fields J . M athematical Ge2 m . Conditional sim ulation of the geometry of fluvio - del2 of lattice system s J . Journal Royal statistical Soc, Se2 bution A . M athematics in O il Production C . Ox2 ford, Clarendon Press, 1988: 235. reservoir parameters w ith specifi

43、ed geostatistical distri2 Softw are L ibrary and U sers Guide M . Oxford Univer2 sity Press, 1996. Sep t 2009 COMPUTING TECHN I UES FOR GEOPHYSICAL AND GEOCHEM ICAL EXPLORATI N . Q O 3 Environmental Science, Zhengzhou 450011, China . COM PU T I N G TECHN IQU ES FOR GEO PHYS ICAL AND GEOCHEM ICAL EX

44、PLORA T ION , 2009, 31 ( 5 : 442 Karst topography effective investigation is the s important aspects to guarantee the p roject construction security and the quality in karst area. B y thinking the app lication condition, the ground penetrating radar method is more suitable for the shallow layer comp

45、 lex karst topography investigation compared w ith other ge2 ophysics method. This paper introduced ground pene2 trating radar p rincip le in karst topography investiga2 s tion, analyzes radar wave dissem ination characteristic and the curve characteristic in the karst area, and p resented the p roj

46、ect examp les to discuss ground pen2 etrating radar app lication in karst topography investi2 s gation. The results show s that this technology can ef2 fectively comp lete the karst area engineering geology investigate work, and the methods have the good p ros2 pects and the app lication value. Key

47、words: ground penetrating radar; karst; frac2 ture; sinkholes; curve characteristic; data interp reta2 tion The Huanglong for ation of carboniferous is m main gas reservoir of the W ubaiti gas field in eastern Sichuan p rovince. B ased on the thoroughly exp lorato2 ry developm ent, it seem s very st

48、ringency that system 2 atically studying the inside building of the depositional system and inter2space allocation, building the model of Sedim entary facie, figuring out the rule of accumu2 late reservoirs distribution, reservoir and penetration units gas reservoir Based on the geology, logging, .

49、thin section, etc. , we make a particular study on the Huanglong formation petrophysical association and s depositional environment by using the analysis tech2 niques of deposit system. W e identify out three kinds of sub 2facies, eight kinds of m icro facies, deal w ith the sedim entary characteris

50、tics and conditions and de2 fine the depositional model of the Huanglong forma2 tion of carboniferous as a saltiness lagoon carbonate depositional model w ith lim netic p lane going up and THE SED IM ENTARY M ICRO FAC IES CHAR2 ACTER IST I O F CARBO N IFERO US RESER2 C VO IR, AND ITS RESERVO IR AND

51、PENETRA 2 T I N UN ITS AT W UBA IT I, EASTERN S I2 O CHUAN PRO V I NCE XU Guo 2sheng, ZHU J ian 2 in, ZHOU Cun 2jian, et m al ( State Key Laboratory of O il and Gas Reservoir . Geology and Exp loration, Chengdu University of Tech2 nology, Chengdu 610059, China . COM PU T I N G TECHN IQU ES FOR GEO P

52、HYS ICAL AND GEO 2 CHEM ICAL EX PLORA T ION , 2009, 31 ( 5 : 447 down. Moreover, since the carboniferous reservoir s p roperties are generally bad and aeolotrop is is pow 2 m erful, it per s meability become a important element of controlling the gas well deliverability capacity, s therefore, to div

53、ide the reservoir and penetration units into m iddle 2tall penetration division and low penetra2 tion division of gas reservoir The study show s that the . reserve and penetration pattern disp lay as two high and three low in gas reservoir range which is con2 firmed by well drilling, which p rovide a science foun2 dation for the formulation of development regulatory excogitation in intermediary and later stage of gas field. Key words: carbonate; sedi entary m icrofacies; m reservoir and penetration uni