基于开源软件序贯指示沉积相建模后处理方法原理及改进.docVIP

基于开源软件序贯指示沉积相建模后处理方法原理及改进 　　摘要：文章介绍了一种基于开源软件进行沉积相建模的新方法，即采用序贯指示模拟（SISIM）和基于最大后验概率更新法（MAPS）相结合进行沉积相建模，并对传统的MAPS后处理方法提出了一种确实可行的、有效的改进方法。传统的MAPS方法使用5X5方形光滑窗口进行搜索，光滑效果有限，文中提出依据地质信息在软件参数文件中调整相应的光滑窗口的大小和方向，使用与主次变程相关的，沿物源方向的，大小合适区域地质经验的椭圆形光滑窗口，得到的模拟结果更加光滑、连续，更加符合实际地质认识。同时开源软件分析环境的搭建使得从事相建模及后处理的过程更透明，更可控，更好的指导矿产的勘探工作。 　　关键词：开源软件 序贯指示模拟（SISIM） 最大后验概率更新法（MAPS） 沉积相建模后处理 改进的MAPS 　　中图分类号：P618.13;TP391.4 文献标识码：A 文章编号：1007-9416（2014）09-0118-02 　　序贯指示模拟方法（SISIM）[1]是基于像元的随机模拟方法中适应性较强，能够综合多种信息并且使用灵活的一种方法，可以模拟具有复杂各向异性的地质体，对于具有不同连续性分布的变量（如沉积相），可利用不同的变差函数进行表征，从而可建立各向异性的模拟结果，是进行沉积相建模的经典方法，是目前公认的运用最广泛的相建模方法。该方法建立沉积相模型时，给出的模拟实现在局部会出现星星点点的零散分布，在边界部位尤为突出，微相显得不连续，这与地质体的实际形态差别很大，甚至不符合地质认识，非常有必要对输出的模型进行后处理以得到更加光滑连续的相模型。 　　最大后验概率更新法（MAPS）是针对序贯指示模拟法得到的模拟实现中微相边界不光滑的缺点而提出的一种后处理方法，使用它可以对基于像元的模拟实现进行光滑处理。文章将介绍一种利用开源软件（GSLIB程序库[2]）进行沉积相建模的方法和流程，即采用序贯指示模拟和基于最大后验概率更新法相结合进行沉积相建模。并对传统的MAPS后处理方法提出了一种确实可行的、有效的改进方法。 　　在众多的公开发布的程序中，应用最广泛的当属美国Stanford大学的开发的Gslib地质统计学程序库（可以在http：//免费下载）[2]。程序可以在任何操作系统下顺利运行。并且有源程序代码下载，便于感兴趣进行更深入的学习和应用。文章主要运用了GSLIB程序库中的SISIM程序和MAPS程序。 　　1 基于最大后验概率更新法相建模结果后处理及改进 　　1.1 相建模结果后处理的意义 　　在用序贯指示模拟（SISIM）建立沉积微相模型时，得到的模拟实现十分不连续，如图2.1-1所示。 　　假如我们用这些并不十分真实的模拟实现作为基础，用其控制储层属性建模，将导致得到的储层属性模型失真，使其失去意义，并给以储层属性模型为基础的油藏数值模拟以及储层剩余油空间分布的预测带来干扰甚至困难。因此，有必要对通过SISIM得到的沉积相模型进行光滑处理，从而真实再现储层沉积相的空间分布特征及规律[3]。如图2.3-2（b）是某研究区SISIM相建模结果运用传统MAPS法处理后的结果，传统的MAPS方法较好地对模拟实现进行了后处理，图像光滑程度高，砂岩连续性好，而且模拟结果更接近经验地质图形，但其光滑效果有限，还有很多零散的星点分布，对物源方向等重要地质信息反应不明确。 　　1.2 MAPS后处理方法的原理 　　传统MAPS后处理方法，通过5X5方形光滑窗口[如图2.3-1（a）]内各个岩相加权平均求取中心节点各个岩相出现的概率值，并直接将出现概率最高的岩相作为此节点的岩相类型[4]。其确定函数为： 　　（式2.2-1） 　　式（2.2-1）中为权重因子;为条件数据权重;为岩相k实现概率;为岩相k的边缘概率;为位置处初始实现值;S为MAPS口中方形光滑窗口中网格的总个数，如图2.3-1（a）所示光滑窗口中S值为25。 　　图2.3-1（a）中这是MAPS处理中经典的方形光滑窗口，我们运用公式（2.2-1）就可以计算出中心节点各种沉积相出现的概率，并直接将概率最高的岩相作为此节点的岩相类型，如图2.2-1（a）给出的SISIM模拟实现（1、2、3、4分别代表不同的沉积相类型），现根据MAPS原理，由公式2.2-1求出中心节点位置四种沉积相出现的比率分别为：沉积相1为0.4;沉积相2为0.2;沉积相3为0.15;沉积相4为0.25，如图2.2-2所示。 　　根据MAPS原理“直接将概率最高的岩相作为此节点的岩相类型”那么可以将中心待处理点的沉积相类型更新为沉积相1（概率最大为0.4），如图2.2-1（b）。这样就可以处理由于资料不足而导致的随机抽样误差以及概率极小的抽