1、 现代油藏描述以建立定量三维油藏地质模型为最现代油藏描述以建立定量三维油藏地质模型为最终目标。这是计算机技术在油藏描述中广泛应用的结终目标。这是计算机技术在油藏描述中广泛应用的结果,也是提高油藏模拟和开采动态预测精度的要求。果,也是提高油藏模拟和开采动态预测精度的要求。由于计算机技术的发展,地质和数学更进一步的结合,由于计算机技术的发展,地质和数学更进一步的结合,以及地质工作本身向定量化的深入发展,使过去只能以及地质工作本身向定量化的深入发展,使过去只能以各种二维图件来表现油藏地质面貌的传统地质工作以各种二维图件来表现油藏地质面貌的传统地质工作方法已逐步被应用计算机技术建立和显示三维的、定方法
2、已逐步被应用计算机技术建立和显示三维的、定量的地质模型所代替,各种建模技术和计算机软件、量的地质模型所代替,各种建模技术和计算机软件、不断地问世,成为近十几年来油藏描述向油藏表征推不断地问世,成为近十几年来油藏描述向油藏表征推进的主要标志。进的主要标志。 一个完整的油藏地质模型应包括:一个完整的油藏地质模型应包括:构造模型构造模型 储层模型储层模型 流体模型流体模型根据油田不同开发阶段的任务,对油藏地质模根据油田不同开发阶段的任务,对油藏地质模型的精细程度要求不同,依此通常可以把油藏型的精细程度要求不同,依此通常可以把油藏地质模型分为三类地质模型分为三类: :概念模型概念模型 静态模型静态模型
3、 预测模型预测模型 依据油藏描述规模的地质模型分类。为配合油藏模拟进依据油藏描述规模的地质模型分类。为配合油藏模拟进行不同开发问题的研究,实际工作经常需要建立不同规行不同开发问题的研究,实际工作经常需要建立不同规模的地质模型,常用的有:模的地质模型,常用的有: 一维单井地质模型一维单井地质模型 二维砂体剖面模型二维砂体剖面模型 二维砂体平面模型二维砂体平面模型 三维砂体模型三维砂体模型 二维层系剖面模型二维层系剖面模型 三维井组模型三维井组模型 三维油藏整体摸型三维油藏整体摸型 二维层内隔层模型二维层内隔层模型 三维层内隔层模型三维层内隔层模型 油藏地质模型建模技术中的关键点,是如何根据已油藏
4、地质模型建模技术中的关键点,是如何根据已知的控制点资料内插、外推资料点间及以外的油藏特性。知的控制点资料内插、外推资料点间及以外的油藏特性。根据这一特点,建立油藏地质模型方法可分两大类:根据这一特点,建立油藏地质模型方法可分两大类: 确定性模型确定性模型 随机性模型随机性模型目前通行的软件建模方法目前通行的软件建模方法 一般是把整个油藏网块化。先建立井模型,把各井一般是把整个油藏网块化。先建立井模型,把各井同层位网块等时对比相连建立层模型,以同层位网块高同层位网块等时对比相连建立层模型,以同层位网块高程表征油藏构造特征,以非储层网块分隔的储层网块表程表征油藏构造特征,以非储层网块分隔的储层网块
5、表征储层的格架,以储层网块中记入各种储层属性的量值征储层的格架,以储层网块中记入各种储层属性的量值表征这些参数空间的分布和非均质面貌。网块尺寸的大表征这些参数空间的分布和非均质面貌。网块尺寸的大小反映模型的粗细程度;属性量值的精度,特别是无资小反映模型的粗细程度;属性量值的精度,特别是无资料控制点处的内插外推值的精度则反映模型的精度。料控制点处的内插外推值的精度则反映模型的精度。 (一一) 确定性建模方法认为资料控制点间的插确定性建模方法认为资料控制点间的插值是唯一解,确定性的。传统地质工作方法值是唯一解,确定性的。传统地质工作方法的内插编图,就属于这一类。克里金作图和的内插编图,就属于这一类
6、。克里金作图和一些数学地质方法作图也属这一类建模方法。一些数学地质方法作图也属这一类建模方法。开发地震的储层解释成果和水平井沿层直接开发地震的储层解释成果和水平井沿层直接取得的数据或测井解释成果,都是确定性建取得的数据或测井解释成果,都是确定性建模的重要依据。模的重要依据。设研究区域中有设研究区域中有n n个已知数据点,采用线性组个已知数据点,采用线性组合的方式,可以获得影响范围内任意点的估计合的方式,可以获得影响范围内任意点的估计值:值: 式中,式中,ii是与已知数据点是与已知数据点Z(Z(xixi) )有关的加有关的加权系数,它表示各个已知数据点对待估计权系数,它表示各个已知数据点对待估计
7、点的数值的贡献。点的数值的贡献。1.无偏性条件 若要使*VZ为VZ的无偏估计,即要求11Nii 因为 mdxxZEVZEVV)(1 又因 NiiiNiiiNiiVmZEZEZE111*)( 所以要使*VVZEZE,就有11Nii。 2.最优性条件(即估计方差最小条件) 估计方差为 ninjiiiiniEyxCVxCVVC11j1i2),(),(2),( 在无偏条件11Nii下,要求出2E达到极小的权系数i(i=1,2,n),这是个求条件极值的问题,要用拉格朗日乘数法。 令FEiin2121,求 F 对i及的偏导,并整理得: niinjijijniVxCxxC111, 2 , 1),(),( 上
8、式为上式为n+1n+1个方程的普通克里格方程组。个方程的普通克里格方程组。该方程组有该方程组有n1个未知数和个未知数和n1个方程组,因此是有解的。根据克里金方法求得各网格点的估个方程组,因此是有解的。根据克里金方法求得各网格点的估计值后即可以用图形函数库进行编程实现来绘制三维图。计值后即可以用图形函数库进行编程实现来绘制三维图。 储层三维地质建模步骤储层三维地质建模步骤数据库流动单元物性参数模型微构造沉积微相井模型断裂系统模型单砂体模型单砂体属性模型多砂体模型多砂体属性模型构造模型模型应用油田开发油田勘探数模储量工业绘图定向井构造图油藏剖面图砂岩等厚图砂岩顶面形态等孔图等渗图等饱和度图砂体连通
9、图栅状图 数据准备数据准备数据来源:岩心、测井、地震、试井、开发动态。数据来源:岩心、测井、地震、试井、开发动态。从建模内容来看,基本数据类型包括以下四类:坐标数据;分层数据;从建模内容来看,基本数据类型包括以下四类:坐标数据;分层数据;断层数据;储层参数数据。断层数据;储层参数数据。储层数据又分为以下三种:储层数据又分为以下三种: 井眼储层数据;岩心分析和测井解释井眼储层数据;岩心分析和测井解释硬数据:包括井内相、砂体、硬数据:包括井内相、砂体、隔夹层、孔隙度、渗透率、含油饱和度等数据,即井模型。隔夹层、孔隙度、渗透率、含油饱和度等数据,即井模型。 地震储层数据:主要为速度、波阻抗、频率等,
10、为储层建模的软数据。地震储层数据:主要为速度、波阻抗、频率等,为储层建模的软数据。 试井(包括地层测试)储层数据:其一为储层连通性信息,可作为储试井(包括地层测试)储层数据:其一为储层连通性信息,可作为储层建模的硬数据;其二为储层参数数据,因它为井筒周围一定范围内的层建模的硬数据;其二为储层参数数据,因它为井筒周围一定范围内的渗透率平均值,精度相对较低,一般做为储层建模的软数据。渗透率平均值,精度相对较低,一般做为储层建模的软数据。 储层三维地质建模步骤储层三维地质建模步骤 数据集成及质量检查数据集成及质量检查数据集成是多学科综合一体化储层表征和建模的重要前提。集成各种不数据集成是多学科综合一
11、体化储层表征和建模的重要前提。集成各种不同比例尺、不同来源的数据(井数据、地震数据、试井数据、二维图形同比例尺、不同来源的数据(井数据、地震数据、试井数据、二维图形数据等),形成统一的储层建模数据库,以便于综合利用各种资料对储数据等),形成统一的储层建模数据库,以便于综合利用各种资料对储层进行一体化分析和建模。层进行一体化分析和建模。对不同数据来源的数据进行质量检查也是储层建模的十分重要的环节。对不同数据来源的数据进行质量检查也是储层建模的十分重要的环节。为了提高储层建模精度,必须尽量保证用于建模的原始数据特别是硬数为了提高储层建模精度,必须尽量保证用于建模的原始数据特别是硬数据的准确可靠性,
12、而应用错误的原始数据进行建模不可能得到符合地质据的准确可靠性,而应用错误的原始数据进行建模不可能得到符合地质实际的储层模型。实际的储层模型。构造模型的建立构造模型的建立构造模型反映储层的空间格架。因此,在建立储层属性的空间分布之前,构造模型反映储层的空间格架。因此,在建立储层属性的空间分布之前,应进行构造建模。构造模型由断层模型和层面模型组成。应进行构造建模。构造模型由断层模型和层面模型组成。储层三维地质建模步骤储层三维地质建模步骤 储层参数模型建立储层参数模型建立在构造模型基础上,建立储层属性的三维分布数值模型。在构造模型基础上,建立储层属性的三维分布数值模型。在构造模型的基础上,利用井数据
13、和(或)地震数据,按照一定的插值(或模在构造模型的基础上,利用井数据和(或)地震数据,按照一定的插值(或模拟)方法对每个三维网块进行赋值,建立储层属性(离散和连续属性)的三维拟)方法对每个三维网块进行赋值,建立储层属性(离散和连续属性)的三维数据体,即储层数值模型。模型网块尺寸越小,标志着模型越细;每个网块上数据体,即储层数值模型。模型网块尺寸越小,标志着模型越细;每个网块上参数值与实际误差值越小,标志着模型的精度越高。参数值与实际误差值越小,标志着模型的精度越高。 模型精度及可信度分析模型精度及可信度分析资料丰富程度及解释精度:资料丰富程度不同,所建模型精度亦不同。对于给资料丰富程度及解释精
14、度:资料丰富程度不同,所建模型精度亦不同。对于给定的工区及给定的赋值方法,可用的资料越丰富,所建模型精度越高。另一方定的工区及给定的赋值方法,可用的资料越丰富,所建模型精度越高。另一方面,对于已有的原始资料,其解释的精度亦严重影响储层模型的精度。如沉积面,对于已有的原始资料,其解释的精度亦严重影响储层模型的精度。如沉积相类型的确定、测井资料的解释精度,等等。相类型的确定、测井资料的解释精度,等等。赋值方法:赋值方法很多,就井间插值(或模拟)而言,有传统的插值方法赋值方法:赋值方法很多,就井间插值(或模拟)而言,有传统的插值方法(如中值法、反距离平方法等)、各种克里金方法等。不同的赋值方法将产生
15、(如中值法、反距离平方法等)、各种克里金方法等。不同的赋值方法将产生不同精度的储层模型。因而,建模方法的选择是储层建模的关键。不同精度的储层模型。因而,建模方法的选择是储层建模的关键。此外,建模人员的技术水平,包括储层地质理论水平及对工区地质的掌握程度、此外,建模人员的技术水平,包括储层地质理论水平及对工区地质的掌握程度、计算机应用水平及对建模软件的掌握程度等,也是影响储层模型精度的因素。计算机应用水平及对建模软件的掌握程度等,也是影响储层模型精度的因素。储层三维地质建模步骤储层三维地质建模步骤 建立数值模型即三维数据体图形显示建立数值模型即三维数据体图形显示主要包括三维图形显示、任意旋转、不
16、同方向切片、从不同角度显示储主要包括三维图形显示、任意旋转、不同方向切片、从不同角度显示储层的外部形态及其内部特点,地质人员和油藏管理人员可据此三维图件层的外部形态及其内部特点,地质人员和油藏管理人员可据此三维图件进行三维储层非均质分析和进行油藏开发管理。进行三维储层非均质分析和进行油藏开发管理。 据三维储层模型进行油气储量计算据三维储层模型进行油气储量计算层总体积;储层总体积以及不同相(或流动单元)的体积;储层孔隙体层总体积;储层总体积以及不同相(或流动单元)的体积;储层孔隙体积及含烃孔隙体积;油气体积及油气储量;连通体积(连通的储层岩石积及含烃孔隙体积;油气体积及油气储量;连通体积(连通的
17、储层岩石体积、孔隙体积及油气储量);可采储量。体积、孔隙体积及油气储量);可采储量。 储层数值模型输出应用与油藏数值模拟储层数值模型输出应用与油藏数值模拟一般需要对储层数值模型进行模型粗化,使细网格的精细地质模型一般需要对储层数值模型进行模型粗化,使细网格的精细地质模型“转转化化”为粗网格模型,使等效粗网格模型能反映原模型的地质特征及流动为粗网格模型,使等效粗网格模型能反映原模型的地质特征及流动响应。响应。储层三维地质建模步骤储层三维地质建模步骤 (二二) 随机建模(随机建模(Stochastic ModelingStochastic Modeling):就是以地质统计学):就是以地质统计学为
18、基础,综合地质学、沉积学等学科的现有知识,根据岩为基础,综合地质学、沉积学等学科的现有知识,根据岩心分析、测井解释、地震勘探、生产动态以及露头观察等心分析、测井解释、地震勘探、生产动态以及露头观察等多种来源的已知数据,对沉积相单元、岩相、砂体、断层、多种来源的已知数据,对沉积相单元、岩相、砂体、断层、裂缝或具体的流动单元的空间分布以及物性参数在空间的裂缝或具体的流动单元的空间分布以及物性参数在空间的变化性进行模拟,从而产生一系列等概率的储层一维或多变化性进行模拟,从而产生一系列等概率的储层一维或多维图象或实现。这些实现表达了储层各种尺度的变化特征维图象或实现。这些实现表达了储层各种尺度的变化特
19、征和内部结构,是高分辨率的、数字化的、定量的储层表征和内部结构,是高分辨率的、数字化的、定量的储层表征方式,而且易于在计算机上重复产生多个这样的实现。每方式,而且易于在计算机上重复产生多个这样的实现。每个实现都是对现实的合理抽样,实现之间的差别反映了由个实现都是对现实的合理抽样,实现之间的差别反映了由于资料缺乏等原因引起的不确定性。于资料缺乏等原因引起的不确定性。 随机建模就是对于一个非均质场中变量随机建模就是对于一个非均质场中变量Z(u) 的分的分布,人工合成反映布,人工合成反映Z(u) 空间分布等概率的模型过程。如空间分布等概率的模型过程。如果模拟中,每个实现在它的已知点位置处的值与原来的
20、果模拟中,每个实现在它的已知点位置处的值与原来的样品值一致,则称之为条件模拟。所谓等概率是指模拟样品值一致,则称之为条件模拟。所谓等概率是指模拟的各个实现,其总体的统计量符合样品或理论的统计量,的各个实现,其总体的统计量符合样品或理论的统计量,这里的统计量是指直方图、累积频率图、变异函数等。这里的统计量是指直方图、累积频率图、变异函数等。由于对应每个模拟点都有一个分布,所以,对预测值不由于对应每个模拟点都有一个分布,所以,对预测值不确定性就有一个定量的描述,可以指出预测值在某一区确定性就有一个定量的描述,可以指出预测值在某一区间的概率。随机建模方法承认地质参数的分布有一定的间的概率。随机建模方
21、法承认地质参数的分布有一定的随机性,而人们对它的认识总会存在一些不确定的因素,随机性,而人们对它的认识总会存在一些不确定的因素,因此建立地质模型时考虑这些随机性引起的多种可能出因此建立地质模型时考虑这些随机性引起的多种可能出现的实现,供地质人员选择。现的实现,供地质人员选择。1 1、储层随机建模原理、储层随机建模原理2 2、储层随机建模方法、储层随机建模方法 离散型模型离散型模型用来描述离散性的地质特征,如砂体的分布,用来描述离散性的地质特征,如砂体的分布,隔层的分布,岩石类型的分布等。隔层的分布,岩石类型的分布等。 连续型模型连续型模型用来描述储层参数连续变化的特性,如孔隙用来描述储层参数连
22、续变化的特性,如孔隙度、渗透率、饱和度的空间分布。度、渗透率、饱和度的空间分布。常用随机建模技术表 技 术 模 型 种类 方法 算 法 示点性过程法 马尔科夫随机域法 截断高斯法 两点直方图法 条 件 模 拟 指标模拟法 离 散 型 模 型 以 目 标 物 体 基 础 非条件模拟 布尔法 模拟退火模拟法 顺序指标模拟法 分形随机函数法 马尔科夫随机域法 条 件 模 拟 LU 分解法 连 续 型 模 型 以 象 元 为 基 础 非条件模拟 转带法 常见随机模拟方法的比较表 随机模拟方法 分类 名称 变量 类型 适用条件 评述 布尔方法 离散型 可以重复而易描述的形状, 主要用于勘探早期砂体和泥岩
23、夹层描述 原理简单,计算量小,易将沉积学知识溶入模拟中;难于条件化 以目标物体为单元 示性点过程 离散型 可以重复而易描述的形状, 如河道等 模拟的结果直观上更容易接受,符合地质规律;难以完全条件化,数学模型复杂 顺序高斯 连续型 变量必须是正态或多元正态分布,要计算变差函数 计算速度快,数学上具有一致性;很难考虑间接信息,要求变量服从正态分布 截断高斯 离散型 变量必须是正态或多元正态分布 适合解决具有排序分布的相组合 顺序指示模拟 连续和离散型 主要用于渗透率和微相的分布, 要知道各指示类型的变差函数 能综合多种信息, 适合解决极值分布问题;计算量大,需要推断很多协方差函数 模拟退火 连续
24、和离散 要构造目标函数, 通常都包括变差函数在内 能综合多种信息,是最灵活的随机模拟方法;计算量大,不易收敛 以象元为模拟单元 分形随机模拟 连续型 变量具分形特征, 如渗透率、 裂缝的分布 快速和经验性强;难考虑间接信息 3 3、储层随机建模步骤、储层随机建模步骤 原始数据库 定性地质概念模型 构造建摸 地层坐标变换 确定参数统计特征值 骨架模型 属性模型 分组实现 随机模型优选 网格粗化 油藏模拟输入 原始模型 野外露头 现代沉积 储层评价等 开发成熟油田 密井网区 1.1.原始数据库的建立原始数据库的建立 原始数据库,即基础地质数据库,如坐标数据、分层数据、原始数据库,即基础地质数据库,
25、如坐标数据、分层数据、断层数据和储层数据等。主要用于断层数据和储层数据等。主要用于建立定性的地质概念模建立定性的地质概念模型,以指导随机建模的过程;型,以指导随机建模的过程;用作模拟的条件限制;用作模拟的条件限制;模模拟参数(统计特征值)的确定。拟参数(统计特征值)的确定。建立模型的构造格架。建立模型的构造格架。 2.2.建立定性的地质概念模型建立定性的地质概念模型 根据原始数据库及其它基础地质资料,建立定性的地质概根据原始数据库及其它基础地质资料,建立定性的地质概念模型,如沉积相分布模型、砂体连续性模型等。定性的地念模型,如沉积相分布模型、砂体连续性模型等。定性的地质概念模型的建立主要是为随
26、机建模选择模拟方法、确定模质概念模型的建立主要是为随机建模选择模拟方法、确定模拟参数、指导实现的优选等服务。拟参数、指导实现的优选等服务。3 3、储层随机建模步骤、储层随机建模步骤 3.3.构造建模构造建模主要根据地震等资料,建立研究区的构造模型。应用高主要根据地震等资料,建立研究区的构造模型。应用高分辨层序地层学建立地层的等时间界面,为后续的骨架分辨层序地层学建立地层的等时间界面,为后续的骨架模型和物性参数模型的建立做准备。模型和物性参数模型的建立做准备。4.4.地层坐标转换地层坐标转换储层形成以后,由于受到构造、压实、剥蚀等一系列地储层形成以后,由于受到构造、压实、剥蚀等一系列地质作用,使
27、得储层变的厚薄不均。而储层特征的分布及质作用,使得储层变的厚薄不均。而储层特征的分布及其持续性和有利方向是沿地层坐标的,并不是笛卡尔坐其持续性和有利方向是沿地层坐标的,并不是笛卡尔坐标。同一地层它的厚度在横向上有变化,但它代表的时标。同一地层它的厚度在横向上有变化,但它代表的时间间隔是一样的。因此,在进行模拟之前,要进行坐标间间隔是一样的。因此,在进行模拟之前,要进行坐标转换。模拟前把地层厚度转换成时间间隔,模拟时垂向转换。模拟前把地层厚度转换成时间间隔,模拟时垂向上按时间间隔划分网格,模拟后再转换成地层厚度。上按时间间隔划分网格,模拟后再转换成地层厚度。3 3、储层随机建模步骤、储层随机建模
28、步骤 5. 5.确定模拟输入的统计特征参数确定模拟输入的统计特征参数主要包括两类:统计特征参数;条件限制参数(原始数据)。主要包括两类:统计特征参数;条件限制参数(原始数据)。 6.6.建立储层骨架模型建立储层骨架模型在构造模型和地层坐标变换的基础上,首先建立能够表征储层较大规在构造模型和地层坐标变换的基础上,首先建立能够表征储层较大规模非均质性的骨架模型。这种非均质性主要是由不同地质体或不同沉模非均质性的骨架模型。这种非均质性主要是由不同地质体或不同沉积相的空间分布引起的。根据地质概念模型、研究目的及现有的技术积相的空间分布引起的。根据地质概念模型、研究目的及现有的技术条件选择合适的随机模拟
29、方法。条件选择合适的随机模拟方法。 7.7.建立储层物性参数模型建立储层物性参数模型在骨架模型的基础上,对不同沉积相带内各种物性参数分别建模。这在骨架模型的基础上,对不同沉积相带内各种物性参数分别建模。这些模型主要用来表征储层各地质体或沉积相内部岩石性质小范围的变些模型主要用来表征储层各地质体或沉积相内部岩石性质小范围的变化。化。3 3、储层随机建模步骤、储层随机建模步骤 3 3、储层随机建模步骤、储层随机建模步骤4 4、相控随机建模原理及方法、相控随机建模原理及方法 目前,储层随机建模技术已经得到充分发展,并形成一些目前,储层随机建模技术已经得到充分发展,并形成一些较成熟、合理的建模策略,在
30、油田应用中取得了良好的效果。较成熟、合理的建模策略,在油田应用中取得了良好的效果。相控建模策略就是其中发展较成熟,应用较广的建模技术。相控建模策略就是其中发展较成熟,应用较广的建模技术。 越来越多的学者认识到沉积微相对储层物性参数有控制作越来越多的学者认识到沉积微相对储层物性参数有控制作用,沉积物的沉积属性用,沉积物的沉积属性(成分、粒度等成分、粒度等)是影响储层孔隙度、渗透是影响储层孔隙度、渗透率等物性参数的最主要因素。在建立储层物性模型时加入相的率等物性参数的最主要因素。在建立储层物性模型时加入相的控制,必然使模拟结果更具合理性和准确性。控制,必然使模拟结果更具合理性和准确性。Damsle
31、sh阐述了阐述了相控建模的基本思路。首先应用离散随机模拟方法建立三维沉相控建模的基本思路。首先应用离散随机模拟方法建立三维沉积相储层结构模型,然后根据不同沉积相积相储层结构模型,然后根据不同沉积相(砂体类型砂体类型)的储层参数的储层参数定量分布规律,分相定量分布规律,分相(或砂体类型或砂体类型)进行连续变量的随机模拟,建进行连续变量的随机模拟,建立三维储层参数分布模型。立三维储层参数分布模型。 4 4、相控随机建模原理及方法、相控随机建模原理及方法具体步骤如下:具体步骤如下:1、利用能够获得的有限钻井的定量资料。研究工区是否存在沉、利用能够获得的有限钻井的定量资料。研究工区是否存在沉积微相对物
32、性参数的明显控制。在单并相分析基础上,采用统积微相对物性参数的明显控制。在单并相分析基础上,采用统计方法,统计各微相内物性参数的平均值、均方差列成表格后计方法,统计各微相内物性参数的平均值、均方差列成表格后比较。如果各微相的上述统计参数表现出明显的差异,表明沉比较。如果各微相的上述统计参数表现出明显的差异,表明沉积微相控制了物性参数展布,对物性建模时应该加入相的控制。积微相控制了物性参数展布,对物性建模时应该加入相的控制。2、利用序贯指示模拟技术建立微相的骨架模型。、利用序贯指示模拟技术建立微相的骨架模型。3 3、建立各沉积微相的相控模型,使用能够获得的有限井的岩心、建立各沉积微相的相控模型,
33、使用能够获得的有限井的岩心分析或测井解释定量资料,获得不同微相物性参数分布规律。分析或测井解释定量资料,获得不同微相物性参数分布规律。对不同微相采用高斯模拟技术建立物性空间分布模型。对不同微相采用高斯模拟技术建立物性空间分布模型。4 4、相控随机建模原理及方法、相控随机建模原理及方法 顺序指示模拟方法可以在模拟时对不同的变量采用顺序指示模拟方法可以在模拟时对不同的变量采用不同的变差函数,从而在模拟过程中同时考虑不同变量不同的变差函数,从而在模拟过程中同时考虑不同变量的各自特点,所以一般选用该方法对微相进行模拟。对的各自特点,所以一般选用该方法对微相进行模拟。对于物性参数,由于顺序高斯模拟是应用
34、广泛的连续性变于物性参数,由于顺序高斯模拟是应用广泛的连续性变量的随机模拟方法,该方法被认为是模拟连续型变量的量的随机模拟方法,该方法被认为是模拟连续型变量的首选方法。因此,针对不同微相采用顺序高斯模拟分别首选方法。因此,针对不同微相采用顺序高斯模拟分别进行模拟。下面对这两种方法分别加以介绍。进行模拟。下面对这两种方法分别加以介绍。 4 4、相控随机建模原理及方法、相控随机建模原理及方法顺序模拟的基本原理顺序模拟的基本原理顺序模拟的基本思路是沿着给定的随机路径顺序地求取各网格点的顺序模拟的基本思路是沿着给定的随机路径顺序地求取各网格点的局部条件概率分布(局部条件概率分布(lcdflcdf),然
35、后从),然后从lcdflcdf中随机抽取模拟值。基本中随机抽取模拟值。基本算法如图所示算法如图所示 随机地选择一个还没有被模拟的网格点。随机地选择一个还没有被模拟的网格点。 在该处估计局部条件概率分布(在该处估计局部条件概率分布(LCPD)。)。 从从LCPD中随机地抽取一个值。中随机地抽取一个值。 把新模拟的值包括到条件数据集中。把新模拟的值包括到条件数据集中。 重复步骤重复步骤,直到所有的网格点都被模拟,直到所有的网格点都被模拟4 4、相控随机建模原理及方法、相控随机建模原理及方法顺序指示模拟顺序指示模拟顺序指示模拟既可用于离散变量,又可用于连续变量的顺序指示模拟既可用于离散变量,又可用于
36、连续变量的随机模拟。该方法不需要对原始条件数据分布的参数形随机模拟。该方法不需要对原始条件数据分布的参数形式作任何假设,而是在现有资料的基础上,通过一系列式作任何假设,而是在现有资料的基础上,通过一系列的门槛值把条件数据转化成指示数据。根据各离散变量的门槛值把条件数据转化成指示数据。根据各离散变量的指示变差函数,采用指示克里金法对每个网格点处的的指示变差函数,采用指示克里金法对每个网格点处的局部条件概率分布(局部条件概率分布(LCPD)进行估计。其主要特点是)进行估计。其主要特点是变量的指示变换,指示克里金和顺序模拟算法。变量的指示变换,指示克里金和顺序模拟算法。 4 4、相控随机建模原理及方
37、法、相控随机建模原理及方法在进行模拟计算之前,首先要进行指示变换,即根据不同的门槛值把原始数据编码成0或1的过程。设Z(X)为X处的参数值对于门槛值为Z0的指示变换可写成 00001),(ZZZZZZI 假定对变量Z进行观测时外界条件不变,对变量Z进行n次观测,得到Zi (i=12n)。当n值较大时,可以用ZiZ0的个数与n的比值来表示变量ZZ0的概率。 在实际地质研究中,在某一点上对变量一般只取一个观测值。通常我们已经假定随机过程是二阶平稳的。 因此, 当样本容量n较大时, 我们可以利用Zi (第i个样本的变量值) Zo的个数与样本容量n的比值来表示变量ZZo的概率,即 niiZZInnZZ
38、PnZZF1000),(1)( |)( |, 对于类型变量,同样可进行指示变换,对于模拟目标区内的每一类相,当它出现于某一位置时,指示变量为1,否则为0(即出现其它相类型时,该相的指示值为0)。 指示值也可以是地质家的解释和推断,因此,可把不同种类和精度的信息都转化成1和0的数据,从而可以进行数据综合。 4 4、相控随机建模原理及方法、相控随机建模原理及方法顺序指示模拟采用指示克里金来估计LCPD,指示克里金不同于其他克里金方法,它主要用作指示预测。该方法通常取待估样品周围一定范围内的样品来进行估计,根据样品的相对间位置及承载的大小而赋予不同的权值。此时,上式可以写成 ),(),(),()(
39、|,1iniixzIxzaxzInxzF 其中I z x( , ) 为预测的量,a z xi( , )为权值,它可通过解下列方程组求得 njxxzCxzxxzCxzajIjiInii 2 , 1,1 1,1xzanii 就某一位置来说,对于每一个门槛值都对应着一个方程组。在变量Z的变化范围内,我们可以用K个门槛值对该范围离散化,因此要求解K个方程组才能求出离散的累计分布函数F Zx nk, |( ),对于ZZkk,1之间的累计分布函数值可以用线性插值等插值方法求得,这样就求出了待估处的局部条件概率分布。 4 4、相控随机建模原理及方法、相控随机建模原理及方法应用顺序指示模拟对离散变量进行条件模
40、拟应用顺序指示模拟对离散变量进行条件模拟), 2 , 1(KkSk 将离散变量kS进行指示变换,变换成指示变量。 设)(uik是kS的指示值,当kSu 时)(uik为1,否则为0。所有的样品均保证K个离散变量是相互排斥的,即保证下列关系式成立。 kkuiuikk, 0)()( Kkkui11)( 计算每种指示变量)(uik的指示变差函数,如果原始数据有丛聚效应时应先进行数据去丛聚效应处理。 4 4、相控随机建模原理及方法、相控随机建模原理及方法应用应用顺序顺序指示指示模拟模拟对离对离散变散变量进量进行条行条件模件模拟拟), 2 , 1(KkSk 进行顺序模拟,模拟的主要步骤如下: a.确定随机
41、访问每个网格节点路径。指定估计网格点的邻域条件数据(包括原始y数据和先前模拟的网格节点的y值)的个数(最大值和最小值)。 b.对指示变量)(uik应用指示克里金来估计该节点处的变量类型属于kS的概率。例如当采用简单指示克里金时,kS在节点u的概率为: nkkkkpIpnIob1*)u()( |1)u(Pr 式中 1 , 0)(uIEpkk是通过推断得到的kS指示协方差类型简单克里金的边缘频率。权系数由使用具有kS指示协方差的简单克里金方程组给出。 c.确定k个离散变量kS的一个顺序(如1,2, ,K) ,这个顺序定义了k个离散变量kS在概率区间0,1上的分布顺序。 d.在0,1上随机产生一个随
42、机数,确定该随机数对应的离散变量的类型,即为该节点处的变量类型。 e.用模拟值更新所有k个指示数据集,并沿随机路径处理下一个网格节点,直到每个节点都被模拟,就可得到一个实现。 顺序指示模拟的输入参数主要为各变量的所占的比例、 各变量的指示变差函数模型及指示代码、条件数据、网格的划分等。 4 4、相控随机建模原理及方法、相控随机建模原理及方法顺序高斯模拟顺序高斯模拟 ), 2 , 1(KkSk 高斯随机域是最经典的随机函数。该模型的最大特征是随高斯随机域是最经典的随机函数。该模型的最大特征是随机变量符合高斯分布(正态分布)。该方法主要用于连续变量机变量符合高斯分布(正态分布)。该方法主要用于连续
43、变量(如孔隙度、厚度)的随机模拟。顺序高斯模拟是一种应用高(如孔隙度、厚度)的随机模拟。顺序高斯模拟是一种应用高斯概率理论和顺序模拟算法产生连续变量空间分布的随机模拟斯概率理论和顺序模拟算法产生连续变量空间分布的随机模拟方法。模拟过程是从一个象元到另一个象元顺序进行的,用于方法。模拟过程是从一个象元到另一个象元顺序进行的,用于计算某象元计算某象元LCPDLCPD的条件数据除包括在给定有效范围内的原始数的条件数据除包括在给定有效范围内的原始数据和已被模拟的网格数据。据和已被模拟的网格数据。 4 4、相控随机建模原理及方法、相控随机建模原理及方法顺序高斯模拟顺序高斯模拟 ), 2 , 1(KkSk
44、顺序高斯模拟变量的步骤如下:顺序高斯模拟变量的步骤如下: 确定代表整个研究区的单变量分布函数(确定代表整个研究区的单变量分布函数(cdf)。如)。如果果Z数据分布不均,则应先对其进行去丛聚效应分析。数据分布不均,则应先对其进行去丛聚效应分析。 利用变量的分布函数,对利用变量的分布函数,对Z数据进行正态得分变换转数据进行正态得分变换转换成换成y数据,使之具有标准正态分布的分布函数。数据,使之具有标准正态分布的分布函数。 检验检验y数据的二元正态性。如果符合则可使用该方数据的二元正态性。如果符合则可使用该方法,否则应考虑其它随机模型法,否则应考虑其它随机模型 4 4、相控随机建模原理及方法、相控随
45、机建模原理及方法顺序高斯模拟顺序高斯模拟 ), 2 , 1(KkSk 如果多变量高斯模型适用于如果多变量高斯模型适用于y变量,则可按下列步骤进行顺序模拟:变量,则可按下列步骤进行顺序模拟: a.确定随机访问每个网格节点路径。指定估计网格点的邻域条件数据(包确定随机访问每个网格节点路径。指定估计网格点的邻域条件数据(包括原始括原始y数据和先前模拟的网格节点的数据和先前模拟的网格节点的y值)的个数(最大值和最小值)。值)的个数(最大值和最小值)。b.应用简单克里金来确定该节点处随机函数的条件分布函数(应用简单克里金来确定该节点处随机函数的条件分布函数(ccdf)的参数)的参数(均值和方差)。(均值
46、和方差)。c.从从ccdf随机地抽取模拟值。随机地抽取模拟值。d.将模拟值加入已有的条件数据集。将模拟值加入已有的条件数据集。e.沿随机路径处理下一个网格节点,直到每个节点都被模拟,就可得到一沿随机路径处理下一个网格节点,直到每个节点都被模拟,就可得到一个实现。个实现。 把模拟的正态值经过逆变换变回到原始变量的模拟值。在逆变换过程中把模拟的正态值经过逆变换变回到原始变量的模拟值。在逆变换过程中可能需要进行数据的内插和外推。可能需要进行数据的内插和外推。4 4、相控随机建模原理及方法、相控随机建模原理及方法), 2 , 1(KkSk 如果多变量高斯模型适用于如果多变量高斯模型适用于y变量,则可按
47、下列步骤进行顺序模拟:变量,则可按下列步骤进行顺序模拟: a.确定随机访问每个网格节点路径。指定估计网格点的邻域条件数据(包确定随机访问每个网格节点路径。指定估计网格点的邻域条件数据(包括原始括原始y数据和先前模拟的网格节点的数据和先前模拟的网格节点的y值)的个数(最大值和最小值)。值)的个数(最大值和最小值)。b.应用简单克里金来确定该节点处随机函数的条件分布函数(应用简单克里金来确定该节点处随机函数的条件分布函数(ccdf)的参数)的参数(均值和方差)。(均值和方差)。c.从从ccdf随机地抽取模拟值。随机地抽取模拟值。d.将模拟值加入已有的条件数据集。将模拟值加入已有的条件数据集。e.沿随机路径处理下一个网格节点,直到每个节点都被模拟,就可得到一沿随机路径处理下一个网格节点,直到每个节点都被模拟,就可得到一个实现。个实现。 把模拟的正态值经过逆变换变回到原始变量的模拟值。在逆变换过程中把模拟的正态值经过逆变换变回到原始变量的模拟值。在逆变换过程中可能需要进行数据的内插和外推。可能需要进行数据的内插和外推。整个顺序模拟过程可以按一条新的随机路径重复以上步骤,以获取一个新整个顺序模拟过程可以按一条新的随机路径重复以上步骤,以获取一个新的实现,通常的做法是改变用于产生随机路径的随机种子数。的实现,通常的做法是改变用于产生随机路径的随机种子数。