资源描述：
《应用序贯指示模拟方法模拟沉积微相》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第23卷第2期西南石油学院学报Vol.23No.22001年4月JournalofSouthwestPetroleumInstituteApr2001文章编号:1000-2634(2001)02-0001-04X应用序贯指示模拟方法模拟沉积微相1123冯国庆,李允,林作华,王玉明(1.西南石油学院,四川南充637001;2.华北油田测井公司;3.吐哈油田鄯善采油厂)摘要:沉积微相的正确划分对储层非均质性研究以及储层参数场的分布有重要的指导意义。在以往的沉积微相研究中,微相边界的确定大多是由地质人员根据专
2、家知识推断,具有主观性,难以适应高精度地质模型的要求。引入地质统计学中的序贯指示模拟方法,对沉积微相的平面分布特征进行随机模拟。通过在X油田S2-2小层的沉积微相研究中的应用,可产生多个反映现有微相平面分布的等概率模型,克服了因资料不全等原因而引起的模型的不确定性,结合地质人员的经验,可确定出满意的模型,应用效果良好。关键词:沉积微相;储集层边界;随机模拟;地质统计学中图分类号:TE122.2文献标识码:A可确定出满意的模型。引言储层沉积相的研究是储层描述和评价的最基础1序贯指示模拟原理[2]的工作,沉
3、积微相的正确划分对储层非均质性研究以及储层参数场的分布有重要的指导意义。因此,设{Z(x),x∈A}是区域化变量所对应的随机它是现代油藏描述中不可缺少的环节之一。在油田函数,其中A表示某个区域。对区域化变量的一个阈开发的中晚期阶段,油田主要进入三次采油和寻找值z0,则定义二维指示随机变量I(x;z0)为:剩余油阶段,对储层地质模型的精细程度要求越来1Z(x)≤z0越高,这同时也需要展开储层沉积微相精细研究。I(x;z0)=0Z(x)>z0在以往沉积微相研究中,微相边界的确定大多是由如果随机函数Z(x)有
4、一组观测值{Z(xα),α地质人员根据专家知识推断,具有主观性,难以适应∈(n)},则Z(x)基于这组观测值的条件概率分布[1]高精度地质模型的要求。等于I(x;z0)的条件期望,即近年来出现的随机模拟技术是解决上述问题的E{I(x;z0)
5、z(xα),α∈(n)}=有利工具,它已经发展成为地质统计学的主要内容p{Z(x)≤z0
6、z(xα),α∈(n)}之一。随机模拟是以地质统计学为基础,综合地质从而,可由指示条件期望值E{I(x;z0)
7、z(xα)α,学、沉积学等学科的现有知识,对沉积相单元、岩相、∈
8、(n)}的估计来估计条件概率等组合或具体的流动单元的空间分布以及物性参数p{Z(x)≤z0
9、z(xα)α,∈(n)}在空间的变化进行模拟,从而产生一系列等概率的由线性无偏估计克里格法,条件概率的估计为:储层一维或多维成像或称实现。这些实现表达了储F3(x;z(xα),α∈(n))=0
10、z层各种尺度的变化特征和内部结构,是细致的、分辨p3{Z(x)≤z0
11、z(xα),α∈(n)}=率高的、定量的储层表征方式,而且易于在计算机上n重复产生、展示、编辑和更新。∑λα(x;z0)I(xα;z0)α=1本文将随机
12、模拟中的序贯指示模拟方法应用于其中,上标星号表示估计值;I(x;z0)为取样点沉积微相的研究,结果表明,该方法可产生多个反映在阈值z0的指示值;λα(xα;z0)为对应的指示克里现有参数分布的等概率模型,结合地质人员的经验,格权重,它由如下的克里格方程组确定:X收稿日期:2000-09-22作者简介:冯国庆(1974-),男(汉族),山东荷泽人,现在西南石油学院攻读博士学位,研究方向为油藏描述。2西南石油学院学报2001年n个指示值;∑λβ(x;z0)CI(xβ-xα;z0)+μ(x;z0)=β=1(f
13、)将模拟结果归入条件指示数据集中。CI(x-xα;z0)Pα=1,2,⋯,n(6)重复(5)直到所有的结点都被模拟。n∑λβ(x;z0)=1β=13应用实例式中,CI(x;z0)为协方差函数;(x;z0)为拉格朗日乘数。这样,对于给定的阈值z0就需要一个指示协使用上述模拟方法对西部X油田S2—2小层的方差函数CI(x,z0),并产生一个指示克里格系统沉积微相进行了模拟。研究区S2—2小层主要发育来估计分布函数。4种沉积微相,分别为分流河道、分流间区、心滩及这样,在给定K个阈值下,对任一待估位置,每漫堤,
14、它们以编码1、2、3、4表示。在该区内的149个阈值都对应一个方程组。在变量Z(x)的变化范口井中,四种微相所占比例分别为44.3%、30.2%、围内,用K个阈值对该范围离散化,通过求解每个14.1%和11.4%(见表1)。阈值对应的克里格方程组,可以得到每个阈值下对根据不同微相类型定义了4种指示变量,若某应的累积分布函数F{Zk,x
15、(n)},然后可采用Monto处存在相应的沉积微相,则其指示值为1,否则为0:Carlo法求得随机函数Z(
