油藏地质建模技术

油藏地质建模技术第一页，共一百零四页，编辑于2023年，星期日油藏地质模型类型和建模技术一、油藏地质模型类型二、储层地质建模方法三、储层地质建模工作流程第二页，共一百零四页，编辑于2023年，星期日

油藏地质模型是指反映油气藏分布的基本特征和空间分布规律的地质实体。

目标：对油藏描述研究的最终成果的概括，它是对油藏的形态、储层性质、规模大小及分布、流体性质及空间展布等的高度概括。是能够如实反映实际地质体特征的模型。类型：从不同角度或强调某一特征，如岩石分类、成因结构等，所以不同学者提出各种差异较大的分类方案。一、油藏地质模型类型第三页，共一百零四页，编辑于2023年，星期日1、按地质的研究内容划分

（S.R.Jackson(1989),美国BellGreek油田）

（1）地质模型构造子模型沉积子模型成岩子模型地球化学子模型

（2）渗透层模型

沉积模型+岩石物性数据（3）流动单元模型

最终的模型应该是由许多流动单元块镶嵌叠砌组成的。各单元块的界线应与构造断层的位置、岩性、岩相带以及成岩胶结物类型的分布相对应。（4）定量的流体动态模型一、油藏地质模型类型油藏结构模型流体模型储层参数模型第四页，共一百零四页，编辑于2023年，星期日构造子模型：油藏几何形态及分布，包括断层、褶曲、地应力等构造要素。沉积子模型：描述储层结构特征，这是流体流动模型的基础，它标征沉积环境、沉积相类型、成积和剥蚀过程、沉积事件序列、单描述储集体形状、大小及内部结构。成岩子模型：储层性质，包括物性、孔隙结构、储层非均质性、渗透性分布等。表征成岩序列和成岩事件，储层质量的演变。地球化学模型：描述地层流体类型、分布及流动机制（单相、多相）、流动单元体。渗流层模型：主要依据沉积子模型，把岩石物性数据加进去，使模型定量化，确定哪些构造相带属于高渗透层、低渗透层。模型是三维的，反映不同相带垂直渗透率和水平渗透率的分布。流动单元模型：流动单元是一个横向上、垂向上连续的储集相带，在该单元体内各部位岩性相似，影响流体流动的岩石物质也相似。定量的流体动态模型：可作为油藏的代表，以此为依据计算油藏的产量、产能、分析注水前缘和推进速度以及注水后残余油的空间展布。一、油藏地质模型类型1、按研究内容划分第五页，共一百零四页，编辑于2023年，星期日（粗—细，定性—定量的原则，强调储层的作用分类）（1）概念地质模型（ConceptualModel）:

针对某一种沉积类型或成因储层，把它代表性的储层（非均质性、连续性等）抽象结果加以典型化和概念化，建立一个对这类储层在研究区内具有普通意义的储层地质模型。这种模型不是一具体的储层地质模型而是某一地区（油田）某一类储层的基本面貌模式图。

这个阶段油田仅有少数大井距的探井和评价井的岩心、测井及测试资料以及二维和三维地震资料，因而不能详细地描述储层细致的非均质特征，只能依据少量的信息、借鉴理论上的沉积模式、成岩模式建立工区储层概念模型一、油藏地质模型类型2、按研究程度划分（裘怿楠，1990）第六页，共一百零四页，编辑于2023年，星期日（粗—细，定性—定量的原则，强调储层的作用分类）（2）静态地质模拟型（StaticModel）：

针对某一具体的油藏或开发区一个（或一套）储层，将某储层特征在三维空间的变化和分布如实地加以描述而建立的地质模型。真实反映地下油藏特征，如小层分布的“五等图”

这一模型主要为编制开发方案和调整方案服务，如确定注采井别、射孔方案、作业施工、配产配注及油田开发动态分析等。一、油藏地质模型类型2、按研究程度划分（裘怿楠，1990）第七页，共一百零四页，编辑于2023年，星期日（粗—细，定性—定量的原则，强调储层的作用分类）（3）预测地质模型（PreductualModle）:

在静态地质模型基础上，建立描述井网数十米级或甚至数米级规模的储层变化的地质模型。针对多井控制井点网的储量参数作一定的内插和外推预测图。这一模型主要为在二次采油之后对地下剩余油需进行开发调整、井网加密或进行三次采油一、油藏地质模型类型2、按研究程度划分（裘怿楠，1990）Ng53层含油饱和度立体图第八页，共一百零四页，编辑于2023年，星期日3、按储层结构形式分

（K.J.Weber和L.C.VanGeuns,1989）一、油藏地质模型类型

储层结构指的是储集砂体的几何形态及其三维空间的分布。这一模型是储层地质模型的骨架，也是决定油藏数值模拟中模拟网块大小和数量的重要依据。储层结构模型的核心是沉积模型。壳牌石油公司将不同沉积相形成的储层结构类型归纳为三类，即千层饼状储层结构、拼合板状储层结构和迷宫状储层结构第九页，共一百零四页，编辑于2023年，星期日3、按储层结构形式分

（K.J.Weber和L.C.VanGeuns,1989）一、油藏地质模型类型第十页，共一百零四页，编辑于2023年，星期日Weber(1986)的分类Weber不仅考虑了储层非均质性的规模，同时还考虑了非均质性对流体渗流的影响。第十一页，共一百零四页，编辑于2023年，星期日（1）（HelgeqiH.Haldorson(1986),美国分四级）巨大规模：构造+岩相——油气藏组合大型规模：相带+储层——油气藏宏观规模：储层+物性——内部结构微观规模：孔隙结构+流体——流动单元4、按模型组成规模一、油藏地质模型类型第十二页，共一百零四页，编辑于2023年，星期日（3）储层非均质模型：

(Pettijion(1973)的分类)

4、按模型组成规模一、油藏地质模型类型第十三页，共一百零四页，编辑于2023年，星期日（4）（信荃麟（1992），中国分五级）4、按模型组成规模一、油藏地质模型类型油田规模的地质模型小层规模的地质模型砂体规模的地质模型岩心规模的地质模型孔隙规模的地质模型油田规模油层规模砂体内韵律性和沉积结构岩心规模孔隙规模韩大匡（1995年）提出河流相五级非均质模型：第十四页，共一百零四页，编辑于2023年，星期日（1）地面原型模型：露头观测（2）地下原型模型：密井网水平井井间地震1、确定性建模2、随机建模(StochasticModeling)（1）定量的储层地质知识库推理法（2）条件模拟(ConditionalSimulation)（3）各种地质统计模型(Geostatistic)3、确定性+随机建模：约束建模二、储层地质建模方法第十五页，共一百零四页，编辑于2023年，星期日1、确定性建模（1）露头原型模型建模地质家提出重返露头的口号。选择各种沉积上有代表性的露头，进行详细的研究，把各种技术手段放到露头上去进行试验，一方面从露头上得到很多储层的定量地质知识，建立知识库，作为地下储层建模的依据；二是试验各种技术手段在储层描述中应如何发挥作用。最有名的就是由英国石油公司(BP)和美国能源部投资进行的Gypcy露头工作，已形成了一个综合各种技术手段的油藏描述现场实验室。二、储层地质建模方法第十六页，共一百零四页，编辑于2023年，星期日（1）露头原型模型建模①露头储层研究方法和手段方法：高分辨率层序地层学、储层沉积学和沉积动力学（岩石相、砂体成因单元、沉积体系）、层次结构分析手段：露头实测：航拍，照片镶嵌法取样：通常10x10x8cm大剖面写实二、储层地质建模方法1、确定性建模第十七页，共一百零四页，编辑于2023年，星期日第十八页，共一百零四页，编辑于2023年，星期日第十九页，共一百零四页，编辑于2023年，星期日第二十页，共一百零四页，编辑于2023年，星期日第二十一页，共一百零四页，编辑于2023年，星期日露头区砂体几何形态定量研究PhilipLowry等研究了三角洲前亚相中三种主要储集砂岩体的大小（分支河道、分支河口坝和远砂坝），在野外露头上采集砂岩体大小有意义的参数是：最大长度、最大宽度和最小厚度，统计了每种砂体的宽度与厚度、厚度与长度、宽度与长度的相关关系图，并给出了每种河道砂体几何特征极限值：②露头区储层原型模型内容二、储层地质建模方法1、确定性建模第二十二页，共一百零四页，编辑于2023年，星期日第二十三页，共一百零四页，编辑于2023年，星期日第二十四页，共一百零四页，编辑于2023年，星期日第二十五页，共一百零四页，编辑于2023年，星期日第二十六页，共一百零四页，编辑于2023年，星期日第二十七页，共一百零四页，编辑于2023年，星期日第二十八页，共一百零四页，编辑于2023年，星期日90%的分支河道的宽度小于850m，50%的小于325m，20%的小于75m。90%的分支河道的厚度小于18m，50%的小于9m，20%的小于3.5m。90%的分支河道的宽厚比小于130：1，50%的小于40：1，20%的小于10：1。90%的分支河口坝的厚度小于30m，50%的小于18m，20%的小于7m。90%的分支河口坝的长度（走向）小于9Km，50%的小于6.5Km，20%的小于3.5Km。90%的分支河口坝砂岩的宽度（横向）小于6.5Km，50%的小于3.25Km，20%的小于1.2Km。90%的远砂坝的厚度小于0.55m，70%的小于0.10m。OlavWalderhang等研究了英国约克郡侏罗系河流体系的砂体形态，该区河流沉积物主要是两种类型的砂体：河道砂岩和决口扇砂岩。萨尔特维克组河道砂岩宽100-700m，厚4-19m，主要是垂向加积的有85%的河道砂体的宽/厚比在15-40之间。斯卡尔贝组河道砂体大部分的宽度约100m，厚度4-12m，其中主要由垂向加积形成的带状砂体，厚4-12m，宽/厚比小于15；侧向加积的河道砂体，厚4-6m，宽/厚比小于20-25，由几个带状砂体复合而成的河道砂体，总宽度有几百米，宽/厚比35。决口扇砂体厚2.5m，宽2.5公里，最大长度3-4公里，宽厚比和长厚比分别为100-2000和500-2000。所测得的河道砂体和决口扇砂体的宽厚比大小可以使我们能够以井为依据预测砂体的几何形态，增加连井对比的可靠性。第二十九页，共一百零四页，编辑于2023年，星期日露头沉积体的非均质性C.Ravenne等研究了英国Yorkshire地区悬崖上的中侏罗统露头，悬崖长10公里，高30—150米，通过直升飞机对悬崖处全部中侏罗统露头连续拍摄了照片将近1000张，通过摄影地质研究，区分了不同沉积单元的几何特征，砂岩的内部构造、砂体与周围页岩的接触特征以及非渗透层的连续性等。

用覆盖悬崖全高的16条剖面确定出不同岩相和构造的岩石学特征及其区别，在600米长、30米高的范围对其中的Scalby组地层进行了详细的摄影地质研究，通过照片解释和纵剖面资料区别出一些沉积体，对摄影地质资料进行简化，经区域倾角校正后，又作了数字化，整个悬崖被简化成174行3163列的网格。

用数字化计算了各种主要岩性的纵横向百分含量②露头区储层原型模型内容二、储层地质建模方法1、确定性建模第三十页，共一百零四页，编辑于2023年，星期日第三十一页，共一百零四页，编辑于2023年，星期日第三十二页，共一百零四页，编辑于2023年，星期日第三十三页，共一百零四页，编辑于2023年，星期日在纵向上，下部砂岩百分比很高，约为40—50%，同时泥质砂岩的百分比比较高，在第100行左右是一个泥质为主区（95%），之上有一个狭窄的厚砂带，最上部是一个完全的泥质区。在横向上指示了不同岩性的相对位置除了上述二维研究之外，在悬崖附近进行了三维研究，首先进行了钻井取心，井网为正方形，井距为300米，为了研究小规模的非均质性，井距从300逐渐减小到10米，在每口取心井中，都进行了SP、GR、DEN、含氧指数、倾角测井等，声波测井只在两口井中，每隔25厘米远提取了一块直径23毫米、长4厘米的岩心塞，用于ф、K测量。露头调查研究河流三角洲储集层中沉积体的非均质性第三十四页，共一百零四页，编辑于2023年，星期日露头储层中的渗透率分布方式

TorginJacobsen，HansRendall对英格兰东北部约克郡河流相砂岩露头剖面进行了研究，通过详细测量露头中渗透率的变化，提供各种砂体类型骨渗透率变化的资料。该文主要研究了两种类型的河道砂体（透镜层状砂体、侧向加积砂体）和决口扇砂体。②露头区储层原型模型内容二、储层地质建模方法1、确定性建模第三十五页，共一百零四页，编辑于2023年，星期日第三十六页，共一百零四页，编辑于2023年，星期日第三十七页，共一百零四页，编辑于2023年，星期日第三十八页，共一百零四页，编辑于2023年，星期日第三十九页，共一百零四页，编辑于2023年，星期日第四十页，共一百零四页，编辑于2023年，星期日露头采集约克郡中侏罗流斯卡尔贝组河流－三角洲相沉积的二维剖面，悬崖沿海岸长10公里，高40—50米，过去已有人用这些露头的二维剖面进行地质研究（描述储层的非均质性）作为河流—三角洲储层几何形状的条件模拟基础。微渗透率计是一个可携式机械装置，它能测量空气注入岩石表面时的流动速度和注入压力，用达西定律来确定渗透率，50个岩心塞样品在实验室测量了孔隙度、渗透率，并对5块样品作X—衍射全岩矿物分析。利用露头调查研究储层中的渗透率分布方式在野外用微渗透率计测量了500个测点的渗透率，并取到50个测点的岩心塞，以网格状方式对每种类型的砂体测量渗透率，透镜层状砂体取样间隔横向为12米，垂向为0.25米，侧向加积层状砂体的取样间隔横向为2米，垂向为0.15米，决口扇砂体的横向为15米，垂向为0.10米。透镜状砂层达3米厚25米宽，具大型交错层理，是辫状河道侧向迁移形成的多层结构中的河道砂层，侧向加积砂岩层是小型曲流河的点坝沉积，横向上广泛分布薄层砂岩（0.25—0.5米厚）是同一个曲流河体系的决口扇沉积。第四十一页，共一百零四页，编辑于2023年，星期日分析结果①透镜层状砂岩（辫状河体系）中测量的渗透率值最高，这个砂体纵向上分为四个单元，绝对渗透率的平均值分别为3125md、2152md、1100md、246md。侧向加积层状砂岩（曲流河体系）中测量的渗透率为中等（314md），横向上广泛分布的薄砂层（决口扇沉积）的渗透率值最低（75md）。②沉积相标志（结构、构造）与渗透率之间没有明显的正相关性，因此，认为所观察到的渗透率分布规律主要受成岩作用的控制。③自生石英是主要的成岩矿物相，似乎主要是纯化学反应的结果，压溶作用提供的物质是次要的。石英加大似乎受沉积相类型的控制，仅在很小程度上受沉积构造和结构的控制。④岩心塞资料与野外测量值的对比，野外测量的渗透率值与在测量点所取的岩心塞实验室测量的渗透率的平均相关系数为0.7，在500~2000md之间相关性最好，当渗透率很大时，微渗透率计测出的渗透率偏大，这是由于高渗透率情况下存在湍流，而低渗透率情况下气体泄漏引起的。⑤渗透率的纵横向变化规律不明显，研究表明，成岩作用对渗透率有重要的控制作用，定量的成岩作用研究表明，这套地层的成岩作用受砂体类型的控制，而沉积构造和结构对其影响较小，换句话说，成岩作用受砂岩类型、体积、几何形状和周围沉积的控制，而不受粒度和沉积构造的控制。利用露头调查研究储层中的渗透率分布方式第四十二页，共一百零四页，编辑于2023年，星期日（3）储层地震模型

以高分辩率的三维地震为基础，利用其覆盖率高的优势，直接追踪井间砂体和求取储层参数。目前遇到的关键问题是分辩率还满足不了油田开发研究单砂体的要求。但对其前景大家都寄以很大的厚望。二、储层地质建模方法1、确定性建模三维多道反演结果第四十三页，共一百零四页，编辑于2023年，星期日这是国际上近20年来兴起，发展很快的一项热门技术。其主要思路是：选择储层砂体在地面出露的露头，进行详细测量和描述，取样密度达到几十厘米的网络(1英尺×1英尺)，把这类砂体的储层物性(如渗透率)的空间分布，原原本本的揭示出来，以此作为原型模型。从中利用地质统计技术寻找其物性空间分布的统计规律，以此统计规律就可以去预测井下各类储层的物性分布。国外一些主要石油科研机构、院校和油公司都在开展这一工作。二、储层地质建模方法2、随机建模方法第四十四页，共一百零四页，编辑于2023年，星期日（1）随机建模的必要性地下储层本身是确定的储层的随机性质是指那些在现有资料不完善的条件下，需要通过猜测确定的储层性质随机模拟结果比拟合结果更贴近“真实”二、储层地质建模方法2、随机建模方法

挪威NorskHydroA/S公司HelgeH.Haldorsen和EivindDamsleth（1990）对随机模拟技术作了较全面的介绍。他们认为，在油藏工程上应用随机模拟有6个原因，即

1）有关储层空间大小、内部几何结构和岩石性质变化的资料不完整；2）储层建块或相空间配置复杂；3）岩石性质随空间位置和方向的变化难以掌握；4）在平均化时特征值和岩石体积之间的关系不清楚；5）静态油藏资料（沿井的点值，如K、H、φ、Sw和地震数据）相对于动态油藏资料（与时间有关的影响因素、岩石结构影响采收率等）的丰度；6）方便和快速。

第四十五页，共一百零四页，编辑于2023年，星期日EstimationversusSimulationEstimation:locallyaccurateandsmooth,goodforvisualizingtrends,noassessmentofglobaluncertaintySimulation:honorsspatialvariability,goodforflowsimulation,assessmentofuncertaintywithalternativerealizations插值与模拟第四十六页，共一百零四页，编辑于2023年，星期日PermeabilityCrossSection:KrigingGRIDSTATGriddingSize-1000-1550-1500-1450-1400PermeabilityDistance-10001000300050007000SubseaDepth-1550-1500-1450-1400.051581270S2S3S4GRIDSTAT5-16-97PermeabilityCrossSection:SimulationPermeabilityCrossSection:SimulationGRIDSTATGriddingSizeDistance-10001000300050007000-1550-1500-1450-1400Distance-10001000300050007000SubseaDepth-1550-1500-1450-1400.051581270S2S3S4GRIDSTAT5-16-97第四十七页，共一百零四页，编辑于2023年，星期日StochasticReservoirModelingRealityDistributionofRock/FluidPropertiessingletruedistributionRecoveryProcessactualprocessimplementedFieldResponse

ModelDistributionoftheRock/FluidPropertiesmultiplestochasticmodelsRecoveryProcessnumericalmodelofprocessFieldResponse岩石和流体属性分布—确定性的油田开采过程—单个真实分布油田响应—真实实现 岩石和流体属性分布—模拟的油田开采过程—多个随机模型油田响应—生产可能性分布第四十八页，共一百零四页，编辑于2023年，星期日（2）储层地质知识库建立二、储层地质建模方法2、随机建模方法储层地质知识库是根据储层地质概念模型，结合静、动态资料对单井模型进行统计分析的基础上建立的能够定量表征储层特征参数统计分布规律的地质知识，它们是储层随机建模的控制条件。①露头调查②地震预测与解释③密井网统计分布概率、平均值、中值、标准偏差、最大值、最小值、变异系数及各参数间的相关系数等。第四十九页，共一百零四页，编辑于2023年，星期日利用计算机研究一些随机建模算法。针对不同沉积类型储层本身的物性参数空间分布特点，形成一些模拟其规律的算法。至今已推出一些商业性软件。如法国石油研究院和地质统计中心等合作推出的一个"HereSim”软件，曾声称最适合河流——三角洲砂体。美国斯坦福大学已有一些通用算法公开出售。挪威一个叫ODEN公司的推出一个STORM软件在市场上推销。美国雪菲龙公司曾大力研究“分形”在储层建模中的应用，等等。这些随机建模方法，总的追求：利用其技术建立的地质模型，交给数值模拟进行历史拟合，一次运算就在大多数井点拟合得较好，这就说明模型基本符合地下实际。或者利用建立的模型布井，选择成功机率最大的井位去实施。我国东部各主力油田缺乏储层的直接露头可供研究，在西部进行了一些工作，主要针对河流砂体储层。大庆的三次加密试验区，井网已达到100口井／km2，利用这一难得的资料，进行精细描述，可供其他地区和油田借鉴。随机建模方法，一些院校都有一些软件研制出来，在西部油田早期评价阶段进行油藏描述建立地质模型上已取得很好效果，在老油田中为研究剩余油饱和度还没有实践经验。

（3）地质统计学方法和技术二、储层地质建模方法2、随机建模方法第五十页，共一百零四页，编辑于2023年，星期日①变差函数分析

储层参数空间结构分析就是对储层参数进行一系列统计分析，研究储层的空间结构及各储层参数的空间分布的连续性和各向异性。目前，研究数据变量空间连续性或各向异性的方法和手段有多种，如协方差函数、相关函数等，但是这些方法只能概括地质体某一特征的整体概况而无法反映其局部变化特征，地质统计学中引入变差函数。

二、储层地质建模方法（3）地质统计学方法和技术①变差函数分析②克立格③离散变量指示型模拟④高斯型模拟⑤渗透率指示模拟⑥模拟退火第五十一页，共一百零四页，编辑于2023年，星期日变差函数

假设空间点只沿某一方向上变化，那么我们把区域化变量在两点处的值之差的方差之半定义为在该特定方向上的变差函数，它表示区域化变量的两值Z(x)、Z(x+h)的变异程度。记为：

在实际工作中，实验变差函数通常用计算算术平均值的方法来获得，即：①变差函数分析

二、储层地质建模方法（3）地质统计学方法和技术第五十二页，共一百零四页，编辑于2023年，星期日CalculatingExperimentalVariograms...Example:StartingWithOneLag(i.e.#6)Startatanode,andcomparevaluetoallnodeswhichfallinthelagandangletolerance.①变差函数分析

第五十三页，共一百零四页，编辑于2023年，星期日(1)半变异函数：(2)交叉半变异函数：(3)指示半变异函数：先对变量进行指示化，再针对指示化变量求取半变异函数。

(4)协方差半变异函数：式中xi，yi为地质区域化变量，它既具有随机性，又具有空间结构性，h为方向矢量，N(h)为数据点对数，实验变异函数求取第五十四页，共一百零四页，编辑于2023年，星期日(5)对数半变异函数:(6)相对半变异函数：(7)绝对半变异函数：(8)平方根半变异函数：式中xi，yi为地质区域化变量，它既具有随机性，又具有空间结构性，h为方向矢量，N(h)为数据点对数，第五十五页，共一百零四页，编辑于2023年，星期日变异函数模型的拟合

对实验变异函数做最优拟合，从而确定变异函数的块金常数Co、变程a和基台值C(0)。选择并建立合适的变异函数模型是地质统计学的基础。常用的变异函数模型有四种：

(1)球状模型：(2)指数模型：

第五十六页，共一百零四页，编辑于2023年，星期日(3)高斯模型：(4)幂模型：第五十七页，共一百零四页，编辑于2023年，星期日②克立格方法WeightedLinearEstimators（加权线性估计）DerivationoftheKrigingEquations（克里格方程的导出）DifferentTypesofKriging（不同类型的克立格）二、储层地质建模方法（3）地质统计学方法和技术第五十八页，共一百零四页，编辑于2023年，星期日WeightedLinearEstimatorsTheBasicideaistoestimatetheattributevalue(say,porosity)atalocationwherewedonotknowthetruevaluewhereureferstoalocation,Z*(u)isanestimateatlocationu,therearendatavaluesZ(ui),

i=1,...,n,andirefertoweights.Whatfactorscouldbeconsideredinassigningtheweights?closenesstothelocationbeingestimatedredundancybetweenthedatavaluesanisotropiccontinuity(preferentialdirection)magnitudeofcontinuity/variability②克立格第五十九页，共一百零四页，编辑于2023年，星期日WeightedLinearEstimatorsAssignalloftheweighttothenearestdata(polygonal-typeestimate)Assigntheweightsinverselyproportionaltothedistancefromthelocationbeingestimated(inversedistanceschemes) wherediisthedistancebetweendataiandthelocationbeing estimated,cisasmallconstant,andisapower(usuallybetween1to3).Howaboutusingthevariogram?thatiskriging②克立格第六十页，共一百零四页，编辑于2023年，星期日KrigingMinimumerrorvarianceestimateoftheunknown:Simplekriging(SK):noconstraint,themeanmisknown.Ordinarykriging(OK):thesumoftheweightsn=1(u)

constrainedtoequal1.Krigingwithatrendmodel(KT)considersthatmisunknownandthatithasamorecomplextrendofknownshapebutunknownparameters.

wherem(u)isthelocalmean,a1,l=0...,Lareunknowncoefficientsofthetrendmodel,andfl(u)arelowordermonomialsofthecoordinates.GSLIBallowsninedifferentmonomialsuptoquadraticfunctionsofthecoordinates(x,y,z,xx,yy,zz,xy,xz,and,yz).KrigingwithanexternaldriftisanextensionofKT.Considersasingletrendfunctionf1(u)definedateachlocationfromsomeexternal(secondary)variable.Cokrigingisforkrigingwithdifferenttypesofdata.IndicatorKriginghasadifferentgoaltobuildaconditionalcdfratherthanderiveanestimate.②克立格第六十一页，共一百零四页，编辑于2023年，星期日③离散变量指示型模拟

（IndicatorSimulationforCategoricalVariables）二、储层地质建模方法④高斯型模拟

（GaussianSimulationforPorosityModeling）⑤渗透率的指示型模拟

（IndicatorSimulationforPermeabilityModeling）⑥模拟退火

（SimulatedAnnealingTechniques）（3）地质统计学方法和技术第六十二页，共一百零四页，编辑于2023年，星期日IndicatorCodingofDatalocalhardindicatordatai(u;z)originatingfromlocalharddataz(u):

i(u;z)=1ifz(u)z,=0ifnotlocalhardindicatordataj(u;z)originatingfromancillaryinformationthatprovideshardinequalityconstraintsonthelocalvaluez(u).Ifz(u)(a,b],then:

0ifza j(u;z)=

undefined(missing)ifz(a,b]

1ifz>blocalsoftindicatordatay(u;z)originatingfromancillaryinformationprovidingprior(pre-posterior)probabilitiesaboutthevaluez(u):

y(u;z)=Prob

{Z(u)z|localinformation} [0,1],andF(z):globalpriorasdefinedhereafterglobalpriorinformationcommontoalllocationsuwithinthestationaryareaA:

F(z)=Prob{Z(u)z},uA第六十三页，共一百零四页，编辑于2023年，星期日IndicatorCodingofDatalocalhardindicatordatai(u;z)originatingfromlocalharddataz(u):

i(u;z)=1ifz(u)z,=0ifnotlocalhardindicatordataj(u;z)originatingfromancillaryinformationthatprovideshardinequalityconstraintsonthelocalvaluez(u).Ifz(u)(a,b],then:

0ifza j(u;z)=

undefined(missing)ifz(a,b]

1ifz>blocalsoftindicatordatay(u;z)originatingfromancillaryinformationprovidingprior(pre-posterior)probabilitiesaboutthevaluez(u):

y(u;z)=Prob

{Z(u)z|localinformation} [0,1],andF(z):globalpriorasdefinedhereafterglobalpriorinformationcommontoalllocationsuwithinthestationaryareaA:

F(z)=Prob{Z(u)z},uA第六十四页，共一百零四页，编辑于2023年，星期日模型种类技术/方法算法离基条示性点过程法散于件马尔科夫随机域法型目模截断高斯法模标拟两点直方图法型体指标模拟法非条件模拟布尔法连基条件模拟模拟退火模拟法续于序贯指示模拟法型象马尔科夫随机域法模素LU分解法型非条件模拟转向带法二、储层地质建模方法（3）随机模拟方法分类第六十五页，共一百零四页，编辑于2023年，星期日离散模型是为了描述具有不连续性质的地质特征而开发的，如河流相地层中砂体的位置和几何分布，砂岩中页岩夹层的分布和规模，裂缝和断层的分布、方向和长度，以及岩相模拟等。在以上各种情况下，空间中的一个点属于有限分类数中的一个，而且仅有一个，随机模型控制着在每一点该分类的数值如何交互影响。例如，该模型可以控制一个砂体可以怎样侵入到另一个砂体中；裂缝是否横切，怎样横切，不同的岩相怎样相互牵制和排斥等等。对于重要层位，如可注水砂岩段，可以不经过实际流动模拟而计算出来。建模期间地质指导是十分重要的。二、储层地质建模方法（3）随机模拟方法分类①离散模型第六十六页，共一百零四页，编辑于2023年，星期日Object-based3-DmodelingofwelldefinedgeometricobjectsPseudogeneticallysimulatedepositionalhistoryStatisticalcontrols:fromcores,logs,seismicandoutcrops基于目标体的模拟（Object-BasedModeling）地质几何体的三维模拟模拟沉积历史统计数据来自:岩芯、测井、地震、和露头二、储层地质建模方法（3）随机模拟方法分类①离散模型第六十七页，共一百零四页，编辑于2023年，星期日

离散模型用于描述离散性质的地质特征，如河流沉积环境的砂体（河道、决口）位置和大小，砂岩中悬浮的页岩的分布和大小，裂缝和断层的分布、方向和长度，相模拟等。随机模拟可以控制一个砂体如何侵蚀另一个砂体，裂缝如何走向，不同相如何相对应等。不用实际流动模拟就可以计算出可注水砂的比例。在建模型时,地质思路是很关键的。常见的离散型模型有：示性点过程、马尔科夫随机场、截断高斯、两点直方图。基于目标体的模拟（Object-BasedModeling）二、储层地质建模方法（3）随机模拟方法分类①离散模型第六十八页，共一百零四页，编辑于2023年，星期日对于每一变量，储层空间内的每一个点都有其不同的值。随机模型可以描述这些变量的中值情况，或可能的横向、纵向变化趋势；中值的可变性；与邻点之间的对比可靠程度。如果不仅是一个变量，则还可描述所研究的各变量之间的协变性。除了中值的变化趋势，大多数连续模型假定储层中各参数具有某种程度的稳定性。即储层的统计特征在空间中并不改变。这一假设并不总是正确的。对所研究储层的地质认识和经验，对于建立连续模型也是十分重要的，但是这种方法似乎更为机械。连续模型用来描述连续变化的地质现象，如岩石特征、地震速度和量纲参数(如储层顶面、OWC等)。二、储层地质建模方法（3）随机模拟方法分类②连续模型第六十九页，共一百零四页，编辑于2023年，星期日Pixel-based3-Dmodelingoflithology,porosityandpermeabilityPreservesmajorheterogeneitiesandstatisticalfeaturesStatisticalcontrols:fromcores,logs,seismicandoutcropsLithologyPorosityPermeability基于象素的模拟（Pixel-BasedModeling）岩性孔隙度渗透率岩性、孔隙度、渗透率的三维模拟表征主要非均质性和统计特征统计数据来自:岩芯、测井、地震、和露头二、储层地质建模方法（3）随机模拟方法分类②连续模型第七十页，共一百零四页，编辑于2023年，星期日

连续型模型用于描述储层参数连续变化的特征，如渗透率、孔隙度和残余饱和度等、地震速度、油藏顶点和油/水界面。模拟退火模拟法顺序指标模拟法分形随机函数法马尔科夫随机域法LU分解法转带法二、储层地质建模方法（3）随机模拟方法分类②连续模型第七十一页，共一百零四页，编辑于2023年，星期日

离散模型更接近于地质师对储层的解释，也更适合于模拟大规模的非均质性和储层不连续性。连续模型则适合于描述岩石特征的空间分布，但却或多或少采取了假设为稳定的形式。

人们自然会想到将这两种模型结合形，成一种混合模型(Alabert和Massonnat.1990；Damsleth等．1992)，这种设想最近已引起普遍关注。大多数离散模型和连续模型可以以一种混合的方法进行结合。

第一步，用—个离散模型描述模型中的大型非均质性。第二步，用不同的连续模型描述离散模型中各单元内的岩石物理参数的空间变化。

二、储层地质建模方法（3）随机模拟方法分类③混合模型第七十二页，共一百零四页，编辑于2023年，星期日二、储层地质建模方法3、确定性+随机建模：约束建模（1）断控定量建模：分断块构造建模（2）相控定量建模：分相带储层建模（3）井约束地震建模：（4）已知模型类比约束建模：输入已知模型第七十三页，共一百零四页，编辑于2023年，星期日离散模型用于描述离散性质的地质特征，如河流沉积环境的砂体（河道、决口）位置和大小，砂岩中悬浮的页岩的分布和大小，裂缝和断层的分布、方向和长度，相模拟等。随机模拟可以控制一个砂体如何侵蚀另一个砂体，裂缝如何走向，不同相如何相对应等。不用实际流动模拟就可以计算出可注水砂的比例。在建模型时,地质思路是很关键的。常见的离散型模型有：示性点过程、马尔科夫随机场、截断高斯、两点直方图。

二、储层地质建模方法4、储层建模方法的适用性（1）相和砂体分布定量模拟第七十四页，共一百零四页，编辑于2023年，星期日相带模拟SimulationResults:岩相分布横剖面P17I27I31P47P01I9I38P08R42I20P36NFaciesSimulation第七十五页，共一百零四页，编辑于2023年，星期日垂直岩性比例矩阵VPM井的岩性资料层的地震属性岩相模拟第七十六页，共一百零四页，编辑于2023年，星期日岩相模拟垂直岩性比例曲线岩性单元井中岩相第七十七页，共一百零四页，编辑于2023年，星期日储层随机模拟第七十八页，共一百零四页，编辑于2023年，星期日连续型模型用于描述储层参数连续变化的特征，如渗透率、孔隙度和残余饱和度等、地震速度、油藏顶点和油/水界面。

模拟退火模拟法顺序指标模拟法分形随机函数法马尔科夫随机域法LU分解法转带法

二、储层地质建模方法4、储层建模及适用方法（2）储层物性定量模拟第七十九页，共一百零四页，编辑于2023年，星期日从岩相到岩石物性P17I27I31P47P01I9I38P08R42I20P36NWESimulationResults:孔隙度分布横剖面第八十页，共一百零四页，编辑于2023年，星期日物性三维方向上的变化第八十一页，共一百零四页，编辑于2023年，星期日某门槛值的显示及体积计算第八十二页，共一百零四页，编辑于2023年，星期日三维地质模型立体图第八十三页，共一百零四页，编辑于2023年，星期日栅状剖面图第八十四页，共一百零四页，编辑于2023年，星期日等深切片图第八十五页，共一百零四页，编辑于2023年，星期日交叉剖面图第八十六页，共一百零四页，编辑于2023年，星期日孔隙度各种方法比较第八十七页，共一百零四页，编辑于2023年，星期日渗透率各种方法比较第八十八页，共一