Petrel中文操作手册2010-（6~9章）

Petrel中文操作 手册 华为质量管理手册 下载焊接手册下载团建手册下载团建手册下载ld手册下载 2010-（6~9章） 第六章 相建模(Facies Modeling) 6.1 Petrel2010版本中相建模技术的大发展 Petrel相建模(Facies Modeling)现有方法主要包括:多点地质统计学相模拟、基于目标的河流相模拟，基于像元的序贯指示模拟、截断高斯模拟，带趋势的截断高斯模拟，指示克里金模拟、神经网络方法，用于详细表征相带分布特征的确定性和随机性相建模技术，而且可以交互使用。同时用户可以导入自己的算法和人工赋值的方法，建立沉积相模型。 Petrel2010在原有版本基础上对相建模方法做了较大的改进，主要体现在以下四个方面: 1)全新的MPS多点统计相模拟算法 在Petrel2010版本中，引进了多点地质统计学相模拟方法，该方法的引进改变了过去传统的两点统计地质学方法，而发展为多点地质学，解决了过去两点统计关系上变差函数的不足，特别是对储层非均质性描述上的不足，多点统计地质学能够充分描述复杂几何形状砂体的空间连续性和变异性。多点统计地质学是建立在多个点的相关关系上，它在解决描述空间变量的连续性和变异性方面得到越来越广泛的应用。斯坦福大学的Journel教授曾指出多点地质统计学是今后地质统计学发展的方向，它的优势已越来越显著。 2)基于快速傅立叶变换的高斯模拟算法 一种新的新的高斯模拟算法在Petrel 2010.1.中被引用，这种算法与GSLIB的序贯随机模拟方式不同。 A 它比SGS运算速度提高了很多 B 它不是序贯算法 C 它可以并行运算 D 它可以进行快速的协同模拟设定 如同上面提到的，这种高斯算法不同于序贯模拟的序贯算法，允许并行计算，采用的算法是傅立叶变换算法，这种算法具有快速、并行、在大的范围变程内优选最合理的变程等优点，这种算法的界面与序贯高斯模拟算法有些类似。 3)进一步改进克里金算法 在2010.1版本中引用了新的克里金算法，这是完全不同于标准GSLIB 克里金的一种 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 ，其搜索性能和并行运算都有很大改进。 克里金可以沿网格方向、也可以沿海平面进行插值。采用沿海平面设定的插值算法在某些状态下更加合理，例如对饱和度的插值模拟(属性插值)。 4)基于目标模型的条件模拟 在河道模拟流程中，不论是河流相模拟还是适应性模拟，在2010.1新版本中，都新增加了3D属性模拟约束条件。 1 在自适应河道模拟下，一种更强的叫做河道迁移的模拟方式也增加了进来，当选择这种运算方式时，河道就会按照概率属性体的趋势向上逐渐增长。 5)概率趋势模拟 一个新的“Trend modeling”流程被引入Petrel2010.1。在这个版本中，提供了一个估算以后用到相模拟中的概率趋势的方法。这可以被认为是垂向比例曲线的标准化。作为输入，算法在井处得到粗化的相数据，然后用块克里金算法估算每种相的平均概率。 其实，在Petrel2008版本中，就已经新增加和改进了一些算法及功能，主要体现在以下三个方面: 1)基于目标的河流相建模算法在原有方法的基础上增加了一个自适应河道建模技术，它不同于传统的河道模拟技术首先生成河道模型，然后将其放到模型中。新的算法运用了截断高斯模拟，充分考虑了所有井资料从而生成每条河道，这种方法的最大特征是能够以较快的速度完成模拟井密度较大的区块。 另外，可以在多口特定的井中通过指定同一条河道连同体的方式，对河道目标体的模拟加以限制，增强模拟的确定性信息。 2)在Petrel2008中增加了一个新的模拟方法就是截断高斯模拟，是适用于具有大批量数据的碳酸盐环境下的基于像元的模拟方法，同时用一种带有趋势约束的截断高斯模拟代替了原有的转相带模拟方法，以适用于模拟三角洲等具有分带特征的沉积环境。 3)Petrel2008的相建模最大改进就是增加了一个相控的功能，所以说相控建模不再是专指相控的物性建模了，对于离散化数据类型同样可以进行相控，这种模拟方法源于分级建模的思想，对于具有较多的相类型，尤其像大型三角洲等沉积环境，我们就可以采用从相,亚相,微相这种分布建模的方式，充分刻画出沉积相模型。 总之，通过这两年的创新及研发，Petrel软件的相建模方法不断得到了改进，为用户提供了更多的思路和好的实现方法。 但由于相建模涉及因素较多，要想建立较准确的相模型，还需要大量的练习，和丰富的工作经验，本手册只给出常规的操作步骤，还需要用户在以后的学习和工作中不断实践。 6.2 Petrel2010版本中多点统计地质学模拟方法 6.2.1技术发展 在Petrel2010版本中，相模拟算法的改进应该是其突出的技术特色，特别是多点地质统计学相模拟的更新。 传统的地质统计学是利用变差函数描述地质变量的相关性和变异性，通过建立在某个方向上两点之间地质变量的变化关系来描述空间变化特性。但是，建立在两点统计关系上的变差函数本身在描述储层非均质性上有很大的不足，它反映的仅仅是空间两点之间的相关性，不能充分描述复杂几何形状砂体的空间连续性和变异性(如河道砂体和冲积扇砂体)。当井资料较少时，用于计算变差函数的点对很少，它也就不能正确反应空 2 间两点之间的相关性。建立在两点基础上的变差函数在储层地质建模中存在一定的不足，而多点地质统计学是建立在多个点的相关关系上，所以它在解决描述空间变量的连续性和变异性方面得到越来越广泛的应用。斯坦福大学的Journel教授曾指出多点地质统计学是今后地质统计学发展的方向，它的优势已越来越显著。 因此，从目前的发展来看，建模发展具有以下两个优势: (1)多点地质统计学是今后地质统计学发展的主要方向，它可以联合反映空间多个位置点的几何形状和相互配位关系，在模拟具有复杂形状地质体分布时，它比序贯指示模拟算法具有更大的优势。 (2)利用多点地质统计学模拟算法具有快速灵活性，模拟的岩相展布图具有一定的真实性，它为储层参数的模拟奠定了基础。 基于上述发展趋势，在Petrel 2010.1版本中新增加了多点地质统计学相模拟算法，这种算法在分类上也是序贯模拟算法，但它比序贯指示模拟算法对相分布的模拟更加成熟。 这种更高级别算法的最大优势在于在模拟时能够比序贯指示模拟考虑到更多的地质方面的因素，在多点统计地质学模拟的模型中能够达到和超过基于目标模拟的效果，这种算法同基于目标模拟算法的优势在于算法的约束条件是基于图像的，它比布尔模型更容易来获得约束条件，这一点非常重要，特别是在井数非常多或者区域的地质情况已经可以通过特定的资料认识清楚地情况下，多点统计地质学模拟的优势就非常明显了。 举例来说，当Petrel的Geobodies技术从地震数据中提取了一些特定的地质体，并且想把这个从地震中局部提取的地质体做到模型中去，用户怎样在不改变模型其它部分(模拟好的相模型)的情况下而实现这一目的呢,或者模型新打了一口井，由于储层分布的变化，需要更新这个相模型，很显然，我们建模人员在具体实现的操作过程中显得非常困难。 多点统计地质学相模拟方法对这类问题的解决提供了可能，对于该方法的更多说明，可以参阅多点统计地质学相模拟用户手册，在第一次采用这种方法之前建议最好查阅手册说明。 图6-1 以训练图像作为多点统计地质学相模拟分布约束条件的模拟 多点统计地质学与变程相比的技术优势 3 变程应该是SIS(斯伦贝谢)随机储层预测的核心技术，变差函数是传统地质统计学中研究地质变量空间相关性的重要工具。然而，变差函数只能把握空间上两点之间的相关性，因而难于表征复杂的空间结构和再现复杂目标的几何形态(如弯曲河道)。 图6-2 变差函数与多点统计地质学对比 如上图，弯曲河道的3种不同的空间结构(图a，b，c)在横向上(东西方向，图d)和纵向上(南北方向，图e)的变差函数十分相似，这说明应用变差函数不能区分这3种不同的空间结构及几何形态，因此，基于变差函数的传统地质统计学插值和模拟方法难于精确表征具有复杂空间结构和几何形态的地质体。 在常规两点地质统计学中，变程起到了一个显著地作用，主要是在模拟时把变程用来作为控制条件模拟属性，而对于多点统计地质学相模拟，训练图像也起到了相应的显著作用，一个训练图像就是三维模型的基础，用来表明我们想要模拟模型的典型几何分布形态，在2010.1版本中，训练图像可以通过以下几个方式来得到，一是在用布尔模型在规则的网格上建立一个简单的无条件3D模型，二是通过手工画相，三是从Petrel软件外部导入模型(如3D数字露头模型)。 6.2.2多点地质统计学的原理 在阐述多点地质统计学之前，首先回顾一下变差函数的地质统计学方法是如何模拟储层岩相分布的，以序贯指示模拟算法为例进行说明。该方法的基本原理简述如下: 假设在模拟区域有k种岩相s1、s2、„„sk，对于模拟目标区域内的每一相，定义指示变量: 1 ，Z ( u) ?sk I ( u) = 0 ，Z ( u) | sk 对于任一待模拟点，其出现第k种相的概率为: 4 P( I = 1 | Z ( u) = s, α)，α为待估点所包括的条件区域，利用两点地质kαα 统计学方法计算该概率是采用克里格方法: 1α=P( I = 1 | Z ( u) = s,α) = E( I) +Σλ[1-E( I) ] kααknαk 其中，λα为克里格方法确定的权系数，它通过求解由变差函数或协方差函数建立的克里格方程组来确定。 多点地质统计学与两点地质统计学的主要区别在于上面的概率的确定方法不同，它首先引入一训练图像，通过在训练图像中寻找与待估点内条件数据分布完全相同的事件的个数来确定概率分布，因此它可以反映出多个位置的联合变异性。 例如，计算图1 (a) 中u 点的概率时，相应的条件数据场为d = { z ( u1) , z ( u2) , n z ( u3) , z ( u4) }，其基本方法是首先要在训练图像(b)中寻找与图(a)中数据分布完全相同的事件的个数，即要在训练图像中找出与图(a)几何完全相同的区域，同时在该区域中相同的位置处z ( u1) , z ( u2) , z ( u3) ,z ( u4) 的值完全相同。在训练图像中一共找到4个既能满足条件数据u1 , u2 , u3 , u4 数值，同时又能满足它们分空间几何形状的事件，在这4个事件中，3个事件的u点的值为0，只有1个事件中u点值为1，因此u 点岩相为1 的条件概率为P{ u = 1 | d} =1/ 4 ，而P{ u = 0 | d} = 3/ nn4 ，这样便可求出了u 点的条件概率。 图6-3 多点统计方法估计条件概率的原理 上述方法不仅考虑了区域内条件数据的值，而且也考虑了条件数据的几何形状。而两点地质统计学只是依靠z ( u1) , z ( u2) , z ( u3) , z ( u4) 的值及各点与u 点距离通过求解克里格方程组来确定u 点的概率，并没有考虑dn的几何形状和各条件数据的配位关系。 基于上述原理，Sebastien St rebelle 提出了snesim 模拟算法，利用该算法可以快速、灵活地模拟岩相分布。该方法的具体步骤为: 5 (1) 利用非条件模拟建立三维训练图像; (2) 定义通过所有待估结点的随机路径; (3) 对随机路径中的任意待估点l ( = 1 ,2 , ?,l) ?定义查找范围内的条件数据; ?保留邻区的数据点; ?在训练图像中寻找与该区域内条件数据完全相同的事件，计算该点岩相的分布概率;?由Monto Carlo 法得到位置处的一个模拟值;?将模拟结果归入条件指示数据集中。 (4) 重复上一步模拟，直到所有的点全被模拟。 训练图像既可以通过非条件模拟求出，也可以通过该地区的地质露头资料分析得出。对训练图像的非条件模拟可以选择非条件的布尔模拟方法。 6.2.3 多点地质统计学模拟方法操作 1)多点统计地质学相模拟的第一步是判别训练图像涵盖的信息，并转换为模拟算法直接识别的内容。这样的信息内容被保存在一种叫做多点统计地质学模式树中，尽管这种典型的模式树比较大，而且用了相对简单的数据格式，并且这种模式是不能看见的，内置到软件中的，但是非常重的一点是这种模式树能够非常简单的被产生和使用(这个图仅仅来演示模式，但具体的模式内容，用户是看不到的)(图7-4)。 A Pattern produced by a small training image. 图6-4 由小的训练图像产生的模式图 操作流程在Property Modeling下面的流程下，如下图: 6 图6-5 操作流程 这是根据训练图像产生模式的一种新方法(图7-6)，虚拟搜索面主要集中在要模拟的网格周围数据上，其结果主要用来模拟有条件的模型分布，多点统计地质学相模拟方法最初模拟是在一个粗的网格上，然后逐渐细化，这种粗的网格叫做多点网格，经常用3作为一个默认数字设定，对于这些参数的设定可以参考MPS手册。 图6-6 多点训练相模式设置 在此操作后，会在input的最下面产生一个新的多点训练相模式结果数据(图7-7)，这个数据可以直接应用到下一步的操作中去。 7 图6-7 多点训练相模式产生的数据 2)多点统计地质学相模拟的第二步就是模拟阶段，这个过程就是选择相模拟流程，并从下拉菜单中选择多点统计地质学相模拟方法，象可以选择模拟的小层一样，在这种模拟方法里还可以选择模拟的区域单元，这个模拟的区域单元可以通过生成的区域属性体来得到。 对不同的模拟区域单元可以此采用不同的训练图像模式(图7-8)，并且允许用户改变各相之间的相关系数，只有激活的区域才被模拟，各模拟区域单元之间可以单独或者按照顺序来模拟，后来模拟的区域单元会依赖于早先的区域单元模拟结果，其他系统的参数控制着模拟时相的分布比例。 图6-8 多点相模拟时训练模式加载 8 3)硬数据及软的第二变量数据在多点统计地质学模拟中，都可以作为约束条件添加到约束条件栏下面，作为第二变量的软数据体，可以形成一个包含可能概率的3D属性体，用它对每一个相进行约束。 同时，硬数据也可以作为约束条件被应用(图7-9)，这也是以3D属性体的形式来加入的，这个硬数据会作为部分约束条件来参与模拟运算，如上面提到的，这样可以将先前模拟的模型或者确定性建好的相数据体添加到新的模型中。 图6-9 多点相模拟时约束软数据加载 4) 最后，比例尺设定和旋转角度设定可以控制训练图像的压缩量和旋转(图7-10)，这样用户就可以方便的将训练图像调整到与模拟的模型几何形态和方向一致，这是由于训练图像的网格应该与模拟的网格一致，举例来说，如果建模工程师认为训练图像在I方向应该扩大一倍，那么他应该输入参数2、1、1作为比例参数，训练图像的旋转是用角度参数来设定的，这些设定可以是数值也可以是属性体。 图6-10 多点相模拟时图像旋转及调整 9 5)查看模拟结果 6.2.4 自定义目标体 用户定义目标体的生成 用户定义目标体的生成可以生成并使用包含一个或多个相代码的非参数目标体。当生成新的训练图像和训练模式的流程时这是非常有用的。 你还可以使用多边形定义的目标体和网格定义的目标体。第一种情况下，只能定义一种相，而第二种情况下，可以定义多种相。两种类型之间可以进行转换，尽管需要根据目标体的复杂情况而定，结果的质量却是千差万别。 1)运行User defined object creation 流程，弹出如下对话框; 2)选择Create new产生一个新目标体; 3)选择边界类型，如网格定义或多边形 定义(点击切换边界类型); 4)定义方向,包括Planar(IJ)平面， Longitudinal(JK)纵向，Transversal (IK)横截面三个方向定义目标体的空间 几何形态; 5)Resolution I、J、K，定义目标体的 分辨率; 6)OK 产生一个目标体，存放在Input文 件夹下。 10 6.2.5 训练图像和训练模式 随着Petrel2009中引入训练模式和多点地质统计学相模拟流程，更多的地质实现加入到模拟世界。新算法的基本输入-训练图像-是这个模拟的首先必备条件。得到训练图像的方法有很多种，许多公司应用露头来建立训练图像(以3D属性输入Petrel)。一种实际的解决方法是使用Petrel中已有的目标体建模的功能来生成一个训练图像。这一过程在2010.1版本中得到了很大的扩展。 在2010.1版本中，生成和使用训练图像以及已有的模式生成组合成一个新的流程“Training image and pattern creation”。操作不需要离开这个界面，训练图像生成后还可以自动进行检测。 打开对话框时，你可以通过定义网格数、网格大小和原点生成一个简单的箱状模型。给网格一个原点可以使你把它放在离工区中实 际数据近的地方，这样可以训练图像的网格可 以作为油藏数据被很容易的放在相同的窗口 下。 特殊训练图像网格会放在“Model”面板下 面。这种类型网格唯一的用途就是存储训练图 像。 ”对网格总体可以通过“Facies modeling 话框得到，也可以通过“Training image and pattern creation”流程中一个快捷键得到， 或者通过流程图直接得到。 除了已有的相建模功能，新的非条件的用 户定义的目标体可以作为“Objet modeling” 中的第四种选择被模拟。边界定义的目标体和 网格定义的目标体都可以被使用。用户定义的 目标体中插入了一个按钮，如果相建模流程打开训练图像的网格是激活的话，这个按钮只用于目标体建模(“Model”面板下黑体显示)。 11 根据多重网格的数目，为了帮助定义一个合适的搜索半径，“Make pattern”栏可以显示用网格数来显示I、J、K方向上有效的半径。 最后，流程中另外一个工具是对训练的图像进行测试，界面中有一个按钮是进行这个操作的，按了这个按钮后，Petrel在后台训练图像网格的基础上会运行多点相模拟，使用刚刚训练的相模式。理论上，使用训练图像测试的结果看起来应该与用来做训练模式的属性类似，这就意味着训练图像的质量与为训练图像设置的参数(搜索面、多重网格的数目和节点信息)的结合足够使训练模式捕捉相的几何形态之间的关系。 总之，现在通过用训练图像得到的概念模型来产生一个地质实现的结果会更容易。 在边界定义目标体的情况下，你选择的是线，而网格定义的目标体中，单个的网格是可以选择的。可以沿着IJ，IK或JK的平面进行编辑。对于网格定义的目标体，有一个附加的选项，可以在已有的属性中提取单个的目标体，提取后可以进一步编辑。这是一个非常有用的选项，可以用来再属性中复制地质体，编辑它们并且用它们来生成训练图像。另外，你可以使用3D属性中的刷子生成单一或多种相的目标体，并且为了生成训练图像把它们“拉”出来变成可以模拟的目标体。用户定义的目标体可以“推”进一个已有的或新生成的属性中来3D窗口检查用它们来模拟之前是什么样子。Petrel2010.1中用户定义目标体的主要目的是支持用更灵活的选项来新生成一个训练图像。在这个版本中，这个功能只能用在训练图像网格中的无条件模拟。 12 6.3概率趋势模拟 一个新的“Trend modeling”流程被引入Petrel2010.1。在这个版本中，提供了一个估算以后用到相模拟中的概率趋势的方法。这可以被认为是垂向比例曲线的标准化。作为输入，算法在井处得到粗化的相数据，然后用块克里金算法估算每种相的平均概率。对每种相来说，算法的结果是一组属性体的集合。对于一种特定相来说，X位置的相概率值等于以X位置为中心的一个计算窗口内该特定相的概率。计算窗口的大小在“Block averaging”一栏中设置。 练习一:如何建立相的概率趋势模型, 1)Upscaled well log data，相数据离 散化; 2)定义变差函数，也可以使用数据分 析功能进行变差函数分析; 3)定义克里金均值，设置2D或3D约 束(如下图)。 4)定义计算窗口大小 block size = (2*size_i+1, 2*size_j+1, 2*size_k+1). 13 使用第二变量来定义局部的克里金均值约束对概率的估计，第二变量可以是2D图(针对每种相)或3D体积(也是针对每种相)或地震属性体。后面的例子中，使用的是数据分析中的属性概率体(与指示模拟中的使用方法类似)。下面图中，使用了属性概率，因此，另外一种约束方法“由面/属性得到的局部克里金均值”被关掉。 在“Expert”一栏中，可以选择简单或普通克里金算法。举一个例子，下面井上识别出了4种相。红色对应的是河道砂，蓝色是决口坝，黄色是漫滩，灰色是泥。决口坝只出现在区域中的一部分和特定的层。我们可以看到河道和决口坝估算的概率体。 14 6.4 Petrel2010版本中基于目标模型的条件模拟 基于目标体的相建模对河流相的模拟起到了非常好的作用，它一般分为基于目标条件模拟和自适应河道模拟，在Petrel2010版本中新增加了3D属性约束条件及迁移河道模拟。 A 一般操作步骤 1) 激活模型New Model/3D Grid (DC)。 2) 双击。 Process diagram/ Property Modeling / Facies Modeling 3) 在建模设置标签下，在Existing Property Status下拉菜单中选择Fluvial Facies Modpel。 4) 从Zones下拉列表选择zone Tarbert 2并取消选择Leave Zone Unchanged 按钮创建一个 realization。 5) 从下拉菜单的方法选择的相体标签选择目标体建模(随机)并点击Add a new channel 。 6) 在设置标签下，选择适当的与Channel and Levee各自相适应的相代码。点击 兰色箭头允许用户从条件数据评估相比例 (图7-11)。 图6-11 基于目标模拟设置 7) 在面板的Channel and Levee标签下定义channel和levee几何体，默认设置可 以使用。 15 8) 进入Background标签并从固定的下拉列表选择泛滥平原背景。 9) 用相同的方法设置所有的zones，回到第4步。 10) 最后点击Apply生成模型。 B 自适应河道目标建模 自适应河道建模方式的最大特点是能够以较快的速度完成模拟井密度较大的区块，而且可以在多口特定的井中通过指定同一条河道连同体的方式，对河道目标体的模拟加以限制，增强模拟的确定性信息。 1)为同一条河道指定连同体，即将同一条河道定义为同一个body。具体做法是:首先为所有井产生一条离散属性的body曲线，然后在well section窗口为河道相设置body，例如将一条河道设置为body1，将另外一条河道设置为body2„„，如图。 图6-12 河道连同体 2)激活模型New Model/3D Grid (DC)。 3)双击Process diagram/ Property Modeling / Facies Modeling。 4)在建模设置标签下，在Existing Property Status下拉菜单中选择Fluvial Facies Modpel。 5)从Zones下拉列表选择zone Tarbert 2并取消选择Leave Zone Unchanged 按钮创建一个 realization。 6)从下拉菜单的方法选择的相体标签选择目标体建模(随机)并点击添加一个adaptive channel，如图6-13。 16 图6-13 自适应河道模拟设置 7)在Facies 栏选择要模拟的河道相channel，点击蓝色箭头自动评估相比例。 8)在other output 面板将Body property复选框选中如图7-13，然后选择一个要使用的body属性。 9)回到Settings面板将勾选，apply运行。 10)模拟结果如图7-14. 图6-14 自适应河道模拟结果 C 新增加功能介绍 在河道模拟流程中，不论是河流相模拟还是适应性模拟，在Petrel2010版本中，都新增加了3D属性模拟约束条件(图6-15)。 17 图6-15 自适应河道模拟2009新版本中增加到约束条件设置 在河道模拟时，作为河道体出现的概率约束条件可以在趋势菜单这一栏下，选择体积约束这一项，由于约束的属性体应该是概率体，因此数值应该在0和1之间。作为自适应河道模拟，选择框也在相应的趋势菜单栏下面。 在自适应河道模拟下，一种更强的叫做河道迁移的模拟方式也增加了进来(图6-16)，当选择这种运算方式时，河道就会按照概率属性体的趋势向上逐渐增长，也就是说，当选择这种方法时，随着河道的逐渐变化生长，它会在概率体的约束下，选择一个生长方向。 采用这种河道迁移算法，就会更加具有适应性，解决了简单的嵌入式点算法只能简单考虑河道初始点分布的问题，而能够使这种概率约束到河道的生长方向和偏移方向了。 18 图6-16 自适应河道模拟2009新版本中增加下切河道模拟设置 6.5 Petrel2010版本中的序贯指示模拟 序贯指示模拟(SIS)是基于象元的随机模拟方法，可用于离散化的连续变量的随机模拟，也可以用于离散型变量的随机模拟。该方法无须假设原始样本服从正态分布，而是通过给定出一系列的门槛值，估计某一类型变量或离散化连续变量低于某一门槛值的概率，以此确定随机变量的分布。是应用指示克里金求取ccdf(累积条件概率)的序贯模拟方法，其主要特点是变量的指示变换、指示克里金和序贯模拟算法。 序贯指示模拟是一种基于象元的随机算法，针对每种相分别设置各自的变差函数和所占比例，在沉积概念模型和地震属性模型约束下建立沉积相模型 A练习步骤 1) 在处理进程表里双击Facies Modeling打开设置对话框(如图7-17) 19 图6-17 序贯指示模拟设置 2) 先选择离散化属性数据Facies(测井曲线scale up或者用户自己解释的单井 相)。 3) 选择Zone settings，选择Zone 并点击Leave Zone Unchanged 按钮。使用 序贯指示模拟方法，并选择Undefined。 4) 由于不同的Zone不同算法的选择，要选择变差函数，概率曲线或者厚度 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 等 在数据分析模块中建立的。如果数据分析都完整作过，而且结果也已经保存， 可以使用这一系列按钮来保证每个Zone的分析结果被调用。 5) 然后进行数据分析(参阅前面数据分析)，例如变差分析，厚度数据分析，概 率曲线分析等，让后对每一个zone都使用同样的设置。 6) 点击Apply运行模型，然后对top zone的Tarbert 1通过选择leave zone unchanged按钮锁定进程 。 B 练习步骤 在序贯指示模拟时，可以考虑协同反演数据体或其它概率属性体来模拟，也可以考 虑沉积相平面图、砂体及储层厚度图、聚类出来的方程等来约束，从而降低单由变程引 起的不确定性(图6-18)。 20 图6-18 序贯指示模拟协同变量设置 6.6 多级相控建模 前面已经提到，Petrel2010的相建模增加了一个相控的功能(从2008版本中引进)，相控建模不再是专指相控的物性建模了，对于离散化数据类型同样可以进行相控，这种模拟方法源于分级建模的思想，对于具有较多的相类型，尤其像大型三角洲等沉积环境，我们就可以采用从相,亚相,微相这种分布建模的方式，充分刻画出沉积相模型。 以三角洲相分级建模为例: 首先用截断高斯模拟方法做一个三角洲亚相级别的模型如图6-6左，模型主要分4种相类型，0-background，5- reef，6-proximal，7-distal，然后用分级相控的方法在distal相中模拟了分流河道相，如图6-19右。 图6-19 多级相控细分前后对比图 21 具体操作步骤如下: 1)激活，打开相建模对话框，选择离散化属性数据Carbonate facies(测井曲线scale up或者用户自己解释的单井相)，如图7-20左; 2)选择模拟方法Truncated Gaussian simulation; 3)选择要模拟的相代码0，5，6，7，设置变差函数，图6-20右; 4)ok 执行，生成模型。 5)从新激活，打开相建模对话框，选择离散化曲线数据Carbonate Facies(Hierarchical); ，选择Carbonate facies相模型做控制，将相类型0,5,66)激活 的模拟方法选择Assign values，指定other property 为原相类型不变，如图6-8左; 7)将相类型7的模拟方法选定基于目标的模拟方法Objected modeling(stochastic)，选择要模拟的相类型1chanel和2levees，进行参数设置(图6-21); 8)ok 执行计算。。 图6-20 未采用分级相控的设置 22 图6-21 多级相控设置 6.7 截断高斯趋势相建模和截断高斯模拟 对于储层连续性性好的沉积以及碳酸盐岩的储层沉积模拟，可以考虑截断高斯趋势相建模和截断高斯模拟。 操作步骤 具体操作步骤如下: 1)激活Facies modeling，打开相建模对话框，选择离散化属性数据facies(测井曲线scale up或者用户自己解释的单井相)，如图6-22; 2)选择模拟方法Truncated Gaussian simulation with trends，如图6-22; 3)选择要模拟的相代码0，1，2，3，并选择下面的趋势为:Geometry trends; 4)设置模拟相的几何形态(Geometry)，图6-23; 5)设置变差函数，图6-24; 6)ok 执行计算，生成模型。 23 图6-22 图6-23 图6-24截断高斯趋势相建模变程设定 6.8 交互式相建模 地震体数据有时能够提供个别的相体图象。交互式相建模能画确定性相体。这个进程是不可返回的，因此在您开始画相之前，做一个模型备份。 操作步骤 1) 在Gullfaks项目下激活Gullfaks Final (DC) ; 2) 复制三维属性体(河流相模型)并重新命名为Facies Object Drawing。交互 式相建模是不能撤销操作的!!! 3) 在进程表中点击相建模进程，显示出相建模工具条。 4) 在3D窗口中显示选择的属性模型，使用Zone Filter显示一个层。 5) 点击刷子按钮。 24 6) 使用侧面按钮 选择一个宽度、高和刷子的形状。i 7) 点击Select Facies 按钮并选择相代码绘制。在3D窗口中显示的模型上可 以开始绘制。 8) 改变刷子类型并保持绘制(如图6-25)。 图6-25交互式相建模成果 25 6.9 其它相建模方法 在Petrel相建模操作中，还有几种建模方法，能够给建模带来很多的方便及快捷，这些 方法包括以下几种: (1)指示克里金插值模拟(图6-26) 快速插值运算。 图6-26指示克里金相建模设置 (2)赋值建模方法(图6-27) 赋值的可以是常数、属性体、面数据或者公式等。 26 图6-27赋值建模方法设置 27 (3)神经网络建模方法(图6-28) 图6-28神经网络建模方法设置 28 (4)用户自定义算法相模拟(图6-29、6-30) 图6-29用户自定义算法程序引进设置 图6-30用户自定义算法程序输入文件设置 29 6.10局部模型更新 这个新的功能可以使在一个多边形内的区域进行相模型和岩石物理属性模型的更新，同时，与相邻没有改变的网格之间保持连续性。当新的局部数据需要加入一个大的模型(比如一口新钻的井)，剩下的区域已经进行了历史匹配并且不想对以前的工作进行重做时，进行局部更新就是一个特别好的办法。 首先，在通过“Scale up well logs”来加入新的粗化数据的时候要保持先前模型的其它部分保持不变。其次，在相建模和岩石物理属性建模中，“Common”一栏中有一个新的选项，可以让用户拿定义好的侧向厚度作为附加的硬数据来对沿着封闭多边形的网格进行设置。当流程在运行时，只有多边形内的网格会被模拟，所有多边形外的网格保持不变。多边形是任意的并且可以把多个封闭的多边形组成一个多边形体。多数情况下是沿着新钻的井画一个多边形，用它新粗化的井曲线来进行属性模拟。 , 如果用的是目标体模拟，局部模型粗化不能对“zone”运行，如果在相模拟对 话框的“Common”栏中为局部模型更新设置了一个多边形，使用目标体模拟方 法的zone会被忽略掉。 , 多边形侧向厚度越厚，对多边形内区域进行模拟的时候受到以前模型数据的影 响就会越大。系统建议的缺省值是3个网格。缺省值的设置会适合大多数的算 法，但对“克里金插值”算法结果的表现会变差，任何情况下这种算法都保持 逆向兼容。所有新的应用都应该使用标准“克里金”算法，它可以用缺省值来 计算。 练习一:怎样运用局部模型更新, 这是一个局部模型更新的例子: 1. 在有大量井的区域进行了测井 曲线离散化后进行了相模拟。有一口新 钻的井，需要把它的相数据和岩石物理 属性数据加入模型。 2. 运行“Scale up well logs” 流程，编辑已有模型的相属性，或者， 先把属性复制一份再进行这个操作(参 见下面“数据分布和全局比例发生了什 么变化”部分来得到更多信息)。选择 “Replace and add new”测井数据并 检查选项“Leave other cells unchanged”。确保列表中包含这口新井 并被离散化。 3. 打开“Facies modeling”流程更 新相属性。在“Common”一栏中，“Local model update”中，使用蓝色的箭头把 闭合的多边形指入需要更新的模型。在 30 这个例子中，封闭的多边形是围绕着这口新井的。 4. 选择侧向网格厚度作为硬数据。 5. 对所有要更新的层设置模拟的参数并运行，所有多边形内的网格会得到一个新的值。 6. 重复步骤2-5对岩石物理属性模型进行更新，需要特别注意的是过程中要使用正确的相模型进行控制。 31 第七章 裂缝建模 裂缝建模需要多步实现，涉及到油藏描述和模拟的方方面面。其主要的宗旨是基于地质概念，充分利用基底解释、断层和成像测井的裂缝知识、通过类比野外露头建立的裂缝概念模型、可预测裂缝成因的地震属性等等，并将这些资料转换成裂缝强度等参数，建立三维的裂缝模型。 Petrel2010裂缝模块主要功能 Petrel2010裂缝模块主要有两部分功能，即建立离散型裂缝几何模型建立及和模型粗化裂缝等效物性模型建立。 图7-1 7.1 裂缝数据准备 1)裂缝数据源多样，数据格式通用、灵活 裂缝信息可以直接采用岩心观测统计数据，或根据测井资料得到的倾角和方位及裂 32 缝类型等解释成果等。数据有效格式为‘Point well data (ASCII)’，即(x，y，z，裂缝参数)，数据以点数据类型存储于Global Well Logs文件夹下。 2)数据观察及修正 单井裂缝信息可以通过右键点击选择显示电子数据表进行浏览，同时可以进行数据添加或编辑;裂缝属性也可以用计算器实现批量添加和浏览。 可以通过右击裂缝图标选择插入属性来添加属性。例如，计算裂缝相对地层层面的倾斜和方位: , 右击Global well logs文件夹下的裂缝点数据 , 在菜单中选择‘Insert new attribute’ , 在pop-up对话框中选择‘Derived’ 作为选项 , 在Derived设置对话框下，用 ‘Dip&azimuth’作为算法 , 选择要计算的层面 , 检查‘Relative’检验栏 , 点击Apply，新的地层倾角和地层倾 角方位被存储为裂缝点数据属性 7.2 裂缝性质分析 对于不同地区，不同的构造运动也会产生不同性质的裂缝类型(一些可能与构造变动有关，一些与沉积因素有关)，裂缝网络是代表裂缝的平面组。同类裂缝由同一个裂缝集产生，每个裂缝网络至少一个裂缝集合。因此对不同性质裂缝类型归类非常重要。 裂缝类型定义可以采用两种方法:一种是直接根据裂缝岩心观测结果(裂缝倾角、倾向、开度等)直接 记录 混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载 入数据库;第二种是通过计算器计算或立体图分析归类。 Petrel2010裂缝模块具有立体图分析功能，即将所有统计结果在同一方位坐标系统内显示，根据其本身性质加以滤波和定量归类。 33 图7-2 7.3 裂缝信息单井剖面建立 根据单井统计结果可以产生裂缝累加曲线、密度曲线和裂缝属性等参数，为井间预测奠定基础。 累积曲线:用来将油藏细分为力学区域的一种客观方法。他们决定于与强度计算密切相关的有效平滑形式，计算结果中带有一定主观性。 密度测井:对于模型密度是重要的元素，是裂缝网格的基本数据。测井曲线是把离散的数据转换成统计的数据来描述裂缝。对于密度测井构造窗的长度对于曲线的光滑度很重要，例如:排出个别的裂缝的作用，用平均的值画点，裂缝密度用累积对数的导数计算，与使用者定义的权一致或有偏差。 裂缝强度曲线是用来描述井上裂缝发育，它可以被粗化到模型中去，并进行模拟，得到的属性体可以作为基础数据用来产生裂缝网络。当产生裂缝强度时，用户需要提供窗口长度，这个窗口长度是用来平均平滑数值。 计算公式为: intensity(md) = (cumulative(md+w/2) -cumulative(md-w/2)) /w 建立裂缝密度测井曲线流程: , 在井剖面窗口显示裂缝数据点。 , 右击测井曲线文件夹下的裂缝数据，在下拉菜单选择插入全井密度曲线。 , 在起始密度点行处，选择密度测井曲线，作为测井曲线类型，核实离散钻孔 的准确性。 图7-3 裂缝密度曲线生成 34 图7-4 裂缝 累计密度曲线 密度曲线 7.4 裂缝参数确定及预测 Petrel 2010裂缝建模提供了确定性及随机模拟两种预测方法。 1)确定性建模方法 可以用蚂蚁追踪得到的断层数据(fault patches )，转换为planar features，直接用到裂缝模型中，或者使用已经存在的裂缝模型。 图7-5 蚂蚁体形成的裂缝片(已经删除断层及错误信息)及确定性建立DFN 2)随机建模方法 随机性建模可以用数学或者属性体计算，同时可以地震数据或其它二维成果图件加以约束模拟。 35 随机方法建模需要输入的参数主要包括:裂缝分布、裂缝几何形态、裂缝方向。 a.裂缝分布:定义模型哪个部分参与计算，垂向上是某个zone或者某个区域，还是整个网格。如果是整个网格，可以用密度属性体加以约束。 密度属性体定义每个单元体积内的裂缝数。可以表示为体密度、线密度和面密度。 裂缝面积/ 体积=裂缝表面积/模型体积，通常应用密度测井曲线 裂缝长度/ 体积=裂缝长度/模型体积 裂缝/体积= 裂缝数量/裂缝体积 使用密度曲线模拟时，允许用户通过定义裂缝长度，选择是否模拟最小的裂缝，因为对渗透率贡献小，因此，模拟结果可能与最初的密度属性值略有不同。 b.裂缝几何形态:定义裂缝外形及长度。 裂缝几何形态:在Petrel中裂缝被简化为面。 裂缝长度:裂缝长度决定了模型中裂缝长度的分配，并且可以用指数、幂数，常数和常量来表示最大和最小的切断。切断需要常态和幂数法分配(最大和最小的复选框变灰)，并且定义在裂缝网络明确模拟的并在第一标签中定义的裂缝强度。各参数可以是常量或者来自于属性。 需要指出的是指数和幂函数分布模型要求输入最大值和最小值，控制裂缝数字和裂缝尺度。密度属性通过这些限制考虑产生裂缝的比例。 c. 裂缝方向:由网格或者趋势面定义裂缝方向，裂缝倾角和方位角的平均数和集中度被输入到Fisher模型中。 Fisher模型的输入数据主要为平均倾角、方位及集中度。 这些值可以在空间中以实际方位来衡量，或者在3D模型或者平面模型中用坐标来计量。 Mean Dip平均倾角—0到90度 Mean dip azimuth平均方位值0-360度 Concentration集中度:数值从0开始到某个值(低值表示分散越大，高值则表示集中度越高)。 d. 裂缝发育长度:裂缝发育长度可以通过露头观察得到，，可以根据经验赋值或由蚂蚁追踪得到。 3)模型约束数据 储层裂缝是非常复杂的渗流系统，随机性强，比较复杂，在预测构成中各以根据地质经验，综合应用多学科手段加以约束。 裂缝发育密度通常与断层距离、油藏沉积类型和地层曲率有关，其中一些可以通过建立几何模型得到。另外地震属性也可以提供裂缝发育信息，外部应力分布图也可以作为趋势约束。 36 图7-6 裂缝密度及分布范围定义 图7-7 裂缝的几何形态定义 图7-8 裂缝的方位定义 37 7.5 裂缝几何模型建立 在以上工作基础上，可以直接生成裂缝几何模型。在Petrel2010中，可以根据实际需要，按整个油层组或分小层进行精细模拟，也可按照精细地质研究结果，按相带进行模拟。 , 在裂缝网络模型流程下打开‘Create discrete fracture network’ , 选择创建新的裂缝模型‘Create new fracture network’ , 在输入数据窗口中，选择自动断裂系统解释结果输入。 图7-9 7.6 等效介质(孔、渗)模型建立 用统计学或基于流体的方法粗化裂缝网格物性，生成可用于模拟运行的输出物性。 , 右击fracture，选择View Calculator 。用存在的Aperture and Transmissivity属性计算这些值。aperture / transmissivity关系确定主 要以来于裂缝渗流公式。 , 右击fracture，选择View Spreadsheet，这样可填充Aperture和 Transmissivity属性栏。 38 图7-10 , 粗化裂缝网格属性 计算裂缝网格属性:将裂缝模型转换为双孔、双渗模型中的裂缝属性。在Eclipse等模拟软件中，主要参数有裂缝孔隙度，裂缝渗透率和sigma 参数(描述基质和裂缝之间的连通性)。Petrel 2010中包含有两种渗透率离散方法，两种方法—Oda方法和基于流体的离散化方法，两种方法都是全张量模拟方法。。在运算过程中这两种方法都要直接使用Golder技术。 Oda方法是数据统计的计算方法，它以单个网格内裂缝的总面积及裂缝的不同参数为基准，估算每个单元渗透率和其它参数。该方法速度很快，但是忽略了裂缝间的连通性，因此，在裂缝密度较小的情况下，可能会低估裂缝的渗透率。 基于流体的离散化是一个严格的方法，它为每各网格都进行特别的限定，并在压力梯度下进行流动模拟，以计算每个方向的渗透率。该计算过程相比前者而言较为漫长，但是它考虑到了流动系统的全部可能几何态。 边界条件应用于每个网格单元。压力梯度用于单元时通常沿测量渗透率的方向放大，有四种选择方式: 不流动:流体不能从侧面流出，此时边界条件没有应用(边线)。 常压:选择界于压力最高值和最低值之间的压力为常压，适用于出口和入口界面。 线性压力:侧向压力将随单元一端(入口)的压力线性增加，低于另一端(出口)压力;线性压力将允许流体沿侧面网格流动。 周期:如果渗透率沿I方向计算，J+和 J-方向将具有同样的值，所有的流体从J+方向流出(或流进)将出J-流进(或流出)，W+/W-相同。 39 7.7 确定性裂缝建模 新类型的数据输入建立裂缝模型，点线面均可作为输入进行裂缝建模，这些数据类型也可以放入一个目录下作为一个输入。 1. 加入的反映裂缝分布的点线面 图8-11 2.点线面显示 图7-12 3.采用确定性裂缝建模方式，可以将加载的数据应用到裂缝建模中去。 图7-13 40 4.在DFN中建立的裂缝片 图7-14 41 第八章 油藏属性建模 这是一个将三维网格中的每个单元赋予属性值的过程，利用测井数据、钻井数据和各属性层面趋势图，采用序贯高斯模拟的算法进行工区内的确定性和随机性属性建模。随机建模可以采用岩相模型、地震属性模型等作为属性模拟的约束条件。同时Petrel特有的科学算法和强大的数据分析功能为合并已有的模型或计算新的模型提供了灵活的约束条件。 8.1 Moving average滑动平均模拟 PETREL包含几种确定性建模方法:克里格方法、函数法、滑动平均法、近点距离等方法， 确定性建模会产生光滑的效果. 克里格方法包含变差的信息，因此产生的非均质模型包含了3D空间各点的相互地质统计的相关性. 确定性方法不会产生局部变量，如果选择100次实现，各实现都相同，下面是滑动平均法的实现练习。 操作步骤 1) 激活Geomodel项目下的一个模型; 2) 打开Petrophysical Modeling 处理流程; 3) 选择 Existing Property 和 Porosity 属性进行建模; 4) 从下拉菜单中选择 一个Zone ，点击 Leave Zone Unchanged 按钮; 5) 选择 Moving average 方法， 其他设置为缺省设置; 6) 点击 OK 产生属性模型并在3D窗口中显示，使用 Zone Filter 显示某一层的属性 (如图9-1); 7) 从WINDOW窗口下打开 Histogram 窗口. 显示属性Porosity 属性使用Zone Filter 观察建立的层位. 比较井曲线，井模型和整个模型的直方图分布。 8) 在Properties 文件夹下双击 Porosity 属性模型，在Statistics 标签下检查统计 结果(如图9-2)。 42 图8-1 图8-2 注意:克里格、移动平均、Closet方法都是确定性方法，确定性方法不会产生局部变量，如果选择100次实现，各实现都相同。 8.2 序贯高斯模拟 Petrel 可以根据序贯高斯模拟建立随机岩石物理模型. 这是用于产生多变量高斯域实现的直接算法。该算法可以产生局部变量和重新产生输入的直方图. 意思是如果基于相同的输入定义100 个实现(使用不同的种子点)，可以得到100个不同的结果，他们 43 都能和输入相匹配，但既然通过分布来输入，每个网格的数值的值会根据这个分布赋值. 如果运行50-100个实现，各模型之间的差别反映了模型的不确定性. 注意, 如果模型本身是正确的,这种不确定性是合理的! 如果输入的数据是错误的，输出的结果必然是错误的。 操作步骤 1(产生一个模型: 1) 激活Geomodel项目下一个模型) 2) 打开Petrophysical Modeling 处理流程, 选择Use Existing Property 和 Porosity 模型. 3) 激活Zone层的Lock icon 按钮取消Zone1层的锁定状态. 4) 选择Sequential Gaussian Simulation方法. 5) 在Variogram选项卡, 选择 Exponential Variogram类型, 2000 =>Major Range, 1000 =>Minor Range, 3 => Vertical Range 和 ,44=》Azimuth.按 OK(如图 9-3) 图8-3 2(使用其他方法，观察其变化: 1) 从新打开Petrophysical Modeling 处理流程的Distribution选项卡 2) 在 Output data range 右边点击 Estimate 按钮 (将估计孔隙度的范围). 设定输出范围 (例如. 0.1-0.25). 图8-4 44 3) 按 OK，观察模型的变化 8.3 相控模拟 一般情况下岩石物理模型要受沉积相的约束， 这将约束岩石物理数值，很明显在泥岩背景下的地质沉积环境，如河流相沉积，主要是代表河道的网格拥有高的孔隙度值，并不是在整个层中平均化孔隙度值。 操作步骤 1) 运行 Petrophysical Modeling 处理流程. 2) 选择 Porosity 属性进行建模. 3) 打开 Zone Settings ，选择 Zone1 进行建模， 点击当前层的 Reset settings ，按 设置为缺省值. 4) 打开 Common settings 设置实现的数目10. 5) 选择 Sequential Gaussian Simulation 方法. 6) 单击 Facies 按钮选择 Fluvial Facies 模型，定义变差设置 图8-5 7) 定义其他相的变差设置. 8) 打开Petrophysical Modeling 处理流程的设置标签， 在Output data range 下点击 Estimate 定义 minimum 和 maximum 孔隙度值.对其他相作相同的操作。 图8-6 45 9) 按 OK 。图(9-7，9-8)是我们工作中实现的一个具体例子。 图8-7沉积相模型 图8-8孔隙度模型 10) 在 Properties 文件夹下, 右键单击，打开Calculator 标签. 选择所有的孔隙 度实现，并计算他们的平均值. Mean_por=(por01+por02+por03+por04+por05+por06+por07+por08+por09)/9 11) 显示平均Porosity 模型，和Fluvial Facies 模型进行对比. 使用 Zone filter 只 显示Zone1 层. 8.4 有趋势约束的模拟 1) 运行 Petrophysical Modeling (双击Petrophysical Modeling) 2) 选择 Porosity 属性进行建模. 3) 打开 Common settings 设置实现的数目10. 4) 打开 Zone Settings ，选择 Zone1 进行建模， 点击当前层的 Reset settings ，按 设置为缺省值. 5) 选择 Sequential Gaussian Simulation 方法. 6) 打开Petrophysical Modeling 处理流程的设置标签， 定义变差的变程，方向，在 Output data range 下点击 Estimate 定义 minimum 和 maximum 孔隙度值.对其 他相作相同的操作. 图8-9 7) 建立孔渗饱模型时，我们除了可以用相约束，同时也可以利用其他的趋势面进行约 束，具体的操作步骤为: a). 打开Trend标签(如图8-6) 46 b). 把描述该层的空隙度的平面趋势图输进来 c). 按APPLY 键观察模型的变化 图8-10 8) 在 Properties 文件夹下, 右键单击，打开Calculator 标签. 选择所有的孔隙度 实现，并计算他们的平均值. Mean_por=(por01+por02+por03+por04+por05+por06+por07+por08+por09)/9 9) 观察利用趋势图约束的效果和相约束效果的区别。 10)我们可以同时使用相约束和趋势约束(如图9-11，9-12是一个例子)。 图8-11用于约束的趋势面 图8-12约束后的孔隙度模型 8.5 高斯随机模拟 2010.1.中引用一种新的高斯模拟算法，这种算法与GSLIB的序贯随机模拟方式不同。 1) 它比SGS运算速度提高了很多 2) 它不是序贯算法 3) 它可以并行运算 4) 它可以进行快速的协同模拟设定 这种高斯算法不同于序贯模拟的序贯算法，允许并行计算，尽管很难给出不同机器的运算效率，在一个4核处理器上，运算的速度会提高2.5倍，这种运算是通过分解算法 47 来实现的，其分解公式可以描述为:条件模拟=克里金+非条件模拟，在这个基础上，运算就允许我们采用并行克里金算法的方式进行，利用正确的算法，非条件模拟速度会非常的快，这里用的算法是傅里叶转换算法，这种算法具有快速、并行、在大的范围变程内优选最合理的变程等优点，这种算法的界面与序贯高斯模拟算法有些类似。 新的协同高斯随机函数模拟通常比类似的序贯高斯模拟优越，主要原因如下: 1) 通常该算法比序贯模拟运算快 2) 快速改变主次变量之间的相关性参数设定 3) 在模拟的变量或者主次变量之间的相关性之间没有系统误差。 图8-13 操作步骤 1) 打开Petrophysical Modeling 处理流程, 选择Use Existing Property 和 Porosity 模型. 2) 激活Zone层的Lock icon 按钮取消Zone1层的锁定状态. 3) 选择Sequential Gaussian Simulation方法. 4) 在Variogram选项卡, 选择 Exponential Variogram类型, 2000 =>Major Range, 1000 =>Minor Range, 3 => Vertical Range 和 44=》Azimuth. 5) 打开Co-kriging页面，在Secondary variable下面勾选Property，然后在下来 菜单下选择约束储层的地震资料。 图8-14 48 第九章 模型粗化与输出输入 9.1 模型粗化 模型粗化是使细网格的精细地质模型“转化”为粗网格模型的过程。在这一过程中，用一系列等效的粗网格去“替代”精细模型中的细网格，并使该等效网格模型能反映原模型的地质特征及流动响应。Petrel2010引入了新的IJK fault 粗化算法，使得逆断层附近不再有洞。 Petrel里的粗化过程主要分成两个步骤: 1. Scale Up Structure: 构造粗化，包括重新划分zones，每个zones里边划分多 少个layers，建立起要做的数模网格(粗的，coarse)和地质模型网格(精细， fine)的对应关系。 2. Scale Up Properties: 应用适当的算法把精细网格中的属性逐个重新采样到 粗化网格。 9.1.1 创建粗化网格 1) 用Pillar Gridding建立一个粗化网格，名字3D Grid (SIMU)也可以自己命 名。 2) 定义粗化网格的网格大小，勾选Make zigzag type faults。 3) Apply。可能会报错，需要重新定义趋势线和断层方向使Pillar Gridding完 成， 4) Ok。 5) 修改新生成的3D Grid (SIMU)的Domain为Elevation Depth(如图10-1)。 图9-1 49 9.1.2 Scale Up Structure 1) 激活练习6-1创建的3D Grid (SIMU)。 2) 双击Process diagram/ Upscaling / Scale Up Structure:。 3) 确定勾选了 Keep locked nodes unchanged 。 4) 定义粗化模型中的Layer个数和基本分配 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 :删除Tarbert 2这个Zone，同 时Tarber 3更名为 Tarbert 3+2，Layer的划分视具体情况而定，可以将所 有的单元都划分1个layer(如图10-2)。 5) OK。这样3D Grid (SIMU)里边就有了Horizons 和Zones。 图9-2 9.1.3 Scale Up Properties 1) 双击Process diagram/ Upscaling /Scale Up Properties。 2) 在精细地质模型中选择要粗化的属性，比如净毛比模型，选中，然后在Scale Up Properties面板中点击按钮输入，如图10-3，同样可以把所有将要粗 化的属性都导入进来，然后分别设置所要采取的粗化方法。 图9,3 50 9.2 局部网格加密 Petrel中油藏地质模型粗化模块增加模拟网格设置技术，允许用户对完成的油藏模型进行局部网格加密。您可以对井周围或者是某个区域进行局部网格加密，下面是一个具体的练习例子是对一个区域的网格加密。 新的局部网格加密( LGR )方法: , 基于层面生成 ; , 在井射孔段周围生成 ; , 渐变的局部网格加密 (井, 层 面)。 操作步骤 1) 首先激活并显示您需要进行加 密操作的模型 图9,4局部网格加密 2) 使用make/edit polygon工具沿 着您需要加密的区域创建一个 封闭的Polygon。如图(10-5) 图9-5 3) 双击ProcessDiagram/structuralmodeling/make local grids 4) 弹出如图(10-6)的参数设置对话框，进行参数设置。 51 设置x，y，z要细化的网格数 选择加密区域所在选择要加密的zone 的Segment 输入上面步骤建立的 Polygong 图9-6 5) 参数设置好以后点击ok，然后观察结果。如图(10-7) 图9-7 9.3 模型输出 Petrel可以输出多种类型的网格数据，涵盖了绝大多数主流数值模拟软件的输入格 52 式。我们以ECLIPSE为例，说明粗化模型的输出操作。 1) 首先输出角点网格，在3D Grid (SIMU)上按右键选择Export„。弹出文件选 择对话框。 2) 文件名给Test.grdecl，文件类型选ECLIPSE Grid keywords (grdecl)(ASCII)(*.*)。该文件的例子如下: SPECGRID 25 25 7 1 F / COORDSYS 1 7 'INCOMP ' / -- Origin (upper left corner, UTM coords) -- (X,Y)=(0.000000,0.000000) -- COORD 494925.406 6695765.000 2816.368 494925.406 6695765.000 3117.095 495088.469 6695653.500 2810.917 495088.469 6695653.500 3102.346 ... ... / ZCORN 2816.37 2810.92 2810.92 2814.42 2814.42 2822.83 2822.83 2565.79 2565.79 2396.01 2396.01 3) 然后再输出断层网格，在3D Grid (SIMU)/Faults上按右键选择Export„。2290.90 弹出文件选择对话框。 ... 4) 文件名Test_fault.grid，文件类型ECLIPSE Fault Data(ASCII)(*.*)。 ... 该文件的例子如下: -- KEYWORD "FAULTS" HAS BEEN WRITTEN FROM PETREL. FAULTS --NAME IX1 IX2 IY1 IY2 IZ1 IZ2 FACE 'MainFaul' 22 22 90 90 1 40 'Y+' / 'MainFaul' 23 23 90 90 1 40 'Y+' / 'MainFaul' 23 23 90 90 1 40 'X+' / ... / 5) 输出Well Completion Data，在3D Grid (SIMU)上按右键选择Export„。弹 出文件选择对话框。 6) 文件名Test_WCompD.grid，文件类型选ECLIPSE Well Completion Data (ASCII)(*.*)。 53 该文件的例子如下: -- COMPLETION SPECIFICATION DATA -- -- WELL -LOCATION- -- NAME I J K1 K2 COMPDAT --------------------------- '34/10-A-21' 30 45 1 1/ '34/10-A-21' 30 46 1 1/ '34/10-A-21' 31 46 1 1/ ... 7) 输出Well Connection Data，在3D Grid (SIMU)上按右键选择Export„。弹/ 出文件选择对话框。 8) 文件名Test_WConnD.grid，文件类型选ECLIPSE Well Connection Data (ASCII)(*.*)。 该文件的例子如下: -- Well connection data for Schedule -- TRAJECTORY_COLUMN_ORDER MD_ENTRY GRID_I GRID_J GRID_K WELL_ENTRY ENTRY_FACE MD_EXIT WELL_EXIT EXIT_FACE PROP1 PROP2 WELLNAME 34/10-A-21 WELLHEAD_I 30 WELLHEAD_J 45 TRAJECTORY 1841.22 30 45 1 455549.38 6784225.50 1819.22 Z- 1820.13 455543.34 6784237.00 1830.52 X- -999.0 1222 1877.30 30 46 1 455543.34 6784237.00 1830.52 X+ 1841.22 455533.53 6784255.00 1848.91 Y- 333.0 1223 ... END_TRAJECTORY -2684.760995 455330.718750 6784624.500000 -2321.159912 54 9.4 模型输入 操作步骤 1) 激活MODEL标签，单击INSERT菜单中NEW MODEL选项，加入一个新的模型，在MODEL 标签下会产生一个名为NEW MODEL 的模型。(如图10-14) 图9-14 2) 鼠标右键单击该文件夹，选择IMPORT选项如图(10-15) 图9-15 55 3) 弹出对话框，选择默认选项，按OK(如图9-16) 图9-16 4) 弹出下面的对话框，按Y 图9-17 5) 弹出下面的对话框，按OK FOR ALL 图9-18 6) 输入的模型放在NEW MODEL文件夹下，打开某一模型浏览。 56