地质油藏建模流程

主要模块介绍一、数据准备本实例中的数据整理如下：wellhead井位坐标文件jinghaoXYkbtopdepthbottomdepthX21-233973816364714261433.0821502195X21-243974070364716291433.082156.12193.1X21-253974257364718491433.082154.42190.4X21-263974480364720961436.52154.82189.8X22-193972535364705161407.562120.32152.3X22-203972803364707951417.462139.12165.1X22-213973010364710401379.72102.62135.6welltop分层文件XYhbwellpointsurfacejinghao397381636471426-716.92Horizonc811X21-23397381636471426-724.92Horizonc8121X21-23397381636471426-735.92Horizonc8122X21-23397381636471426-755.92Horizonc813X21-23397381636471426-761.92Horizonc821X21-23397407036471629-723.02Horizonc811X21-24397407036471629-731.02Horizonc8121X21-24397407036471629-742.02Horizonc8122X21-24397407036471629-754.02Horizonc813X21-24397407036471629-760.02Horizonc821X21-24测井文件准备DEPTHPERM_KPOR_KSW_KVSH_KNTG2140.1250.005901002140.250.005901012140.3750.005901002140.50.00590010二、数据输入1输入WellHeader（井位坐标文件）右键点击输入WellHeader:文件类型里选：wellheads(*.*)2输入WellTops（分层文件）:右键点击WellTops文件夹并选择Import(onSelection)；文件类型里选：PetrelWellTops(ASCII)3输入输入WellLogs右键点击Wells文件夹，选择Import(onSelection)；文件类型：welllogs(ASCII)inputDatalogsspecifylogstobeload加载per,perm,swvash,ntg等数据。设置templatesettings-9999左击OK然后下图设置全选左击OKForAll。在流程窗口左键双击出现下图，在窗口中输入建模名称点apply,再点cancel三、PillarGridding建一个网格边界（工区范围）边界标定了3D网格的侧向延伸。仅在边界内形成3D网格，因此在边界外不会进行储量计算，也不存在构造层面和属性单元。可用创建边界工具，在2D窗口中创建一个边界。同时再这里设置最终模型的横向网格大小。步骤：左击下拉菜单中的并把资源管理器中的wells打勾出现建模的2D图形：把pillargridding单击点亮：左击在2D里浏览全图，出现下图：左击，然后用鼠标左键画边界：边界将要封闭前双击左键，就封闭了.然后左键双击左击apply，出现下图左击是（Y）。Pillar网格化的过程就是一个空间结构生成的过程。在I、J方向上定义网格单元的大小。生成的骨架网格（也叫作pillar网格）定义出了空间结构。创建出的骨架网格不代表任何表面，而是代表了pillar顶部、中部和底部的位置。在下一个进程中（创建地层层面）地层层面会被插入（makehorizon），并连接到pillar上，Z方向上的网格单元也将被定义。Pillar网格化进程完成后，首先会生成一个3DGRID网格。网格化的目的就是要创建均匀分布的矩形网格单元。对生成的网格结果感觉满意后，点击OK以开始构建顶部和底部的骨架网格。在弹出窗口（询问是否将开始构建顶部和底部骨架网格）中点击Yes。四、MakeHorizon在3D骨架网格中加入层面左键双击：窗口中出现下图：左击5次，建立5个层面。出现下图：左击1,然后依次左击2、3，4把层面加上去。出现下图。并把光滑窗口数字0改为2。1234双击窗口中出现下图：直接左击。再cancel.五、Laying细分层仅仅是网格精细化的过程，不是所有输入资料都用在了这个进程中。用户可通过设置单元的厚度、单元的个数或用比例数，来定义网格垂向的分辨率。给定单元厚度时，zone的划分既可以跟随顶部也可以按底部。小层本应该根据将要建立的属性模型来定义。通常，小层的厚度应该是模拟的最薄相的厚度。但是，有很重要的一点应当记住，小层厚度减少时，单元数目会增加，所以不应该插入太多的细节。步骤：下面创建细分层：左键双击出现下图:把分层数改成如下图，两个主力层各分成10个小层。左击apply,再cancel.六、建立几何建模几何属性是通过先前定义好的方法，如网格高度（CellHeight）、总体积（BulkVolume）、创建出的3D属性。每个网格都将赋予一个与所选方法相对应的数值。在进行储量计算和岩石物理属性间的数学运算（如生成含水饱和度属性）时可能会用到。几何属性建模进程允许用户建立几何属性模型，另外还可进行简单的建模操作，如接触面之上的计算，它是计算用户定义的接触面之上的网格单元的高度。选择“单元体积”方法；用总体积作为属性模板，然后点击应用即可生成；左键双击出现下图，左击apply,再左击cancel.打开input资源管理器（出现下图）对右击出现下图：左击设置如下图：把Vsh测井曲线处理成离散的值0（代表砂岩段）和1（代表泥质段）七、离散化测井曲线离散化进程就是给井曲线穿过的网格单元赋值。因为每个网格单元仅能得到一个值，那就要求测井曲线要均匀分布，即离散化。其目的就是要在属性建模时能把井的信息作为输入，即控制井间的属性分布。有一点要明确，离散化之后得到的网格单元将作为属性的一部分，而不是独立出的一项。沿井轨迹的网格单元内分布的值与整个3D离散化之后得到的属性分布是一致的。左键双击定义离散化设置。算法选平均法，以线数据处理测井曲线，使用Neighborcell方法。.注意：*对Vshale设置成如下方式，然后左击*对于连续数据Perm离散化时应该选用Hamonic方法，其它设置完后左击*每次设置完一道测线参数并Apply后，都要把点上，以免被冲掉（那就白做了）。*对于por,sw,ntg选用的方法一样，都用Arithmetic八、对Vsh数据分析数据分析的结果可以直接被相建模和属性建模的模块调用，数据分析分为两类：对离散数据的分析和连续数据的分析。VSH离散数据的分析:1.打开分析窗口,界面如下:2.选择分析的对象,是经过离散化的井点,还是未离散化的测井曲线,或整个模型（默认）.3.是否使用滤波功能，很少使用（默认不用）。4.左击标签,进行相应的分析：选中vsh[U]的Zone1依次左击1、2、3、45、Vsh变差函数的分析对Vshale的4个层（zone1,zone2,zone3,zone4）分别进行变差函数的分析：以zone1层的泥shale为例,其它2、3、4层的操作相同左击标签,进行相应的变差分析：出现下图：1432首先设置主方向的分析参数,包括带宽,搜索半径,步长等,然后再设置次方向和垂向上的参数,这些参数的设置需要用户对本地区数据的大概了解的基础上,否则分析的结果的可信度大大降低.在该例中的分析参数和结果如下图.在分析变差之前,首先大概了解数据的分布情况,然后再调整这些分析参数,这样才能达到比较好的分析效果。1）先把MajorDirection（主方向值）左击，（出现下图）然后设置其它参数，手动调参完毕后左击Apply,手动调整下图的1和2的幅度，使得蓝线尽量与小黑点重合，调好后块金值Nugget设为0。注意：Majorrange值比Minorrange值一般要大些。Type值均设置成spherical，Nolags一般根据经验设置，,本例设为30,这样条柱较多，容易调参。2）MajorDirection手动调参完毕后，再左击次方向（出现下图）,开始手动调参，调完后，左击Apply3）MajorDirection调参完毕后，左击VerticalDirection（出12现下图）垂向厚度一般小于10米，但有时主力层厚要超过10米。块金值Nugget设为0至此，Vshalede的Zone1层的泥质操作完毕。然后是Zone1层的砂层是同样操作。本例中调好的图形如下：主方向调试如下：次方向调试如下：垂向调试如下：然后调试Vsh[U]的Zone2,zone3,zone4，重复前面操作过程，4个层都依次变差分析完成后，对Vsh的数据分析结束。九、相建模Petrel中有几种方法可以用来生成相模型：在这里常使用序贯指示模拟法的相模型随机计算（随机建模）。用序贯指示模拟法(SIS)，建立一个基本的相模型的进程1、双击下的。2、打开相建模流程3、选择Useexistingproperty,属性选择4、选择Zonesettings，选择Zone1(顶层)，左击，把打勾，并设为2，并把LeaveZoneUnchanged按钮按下。以zone层为例：5、选择zone1。取消选择LeaveZoneUnchanged按钮以改变设置并选择序贯指示模拟方法。、左击、、、、各一次，并左击、两次、出现下图：选中zone1打勾左击Apply运行，然后依次对zone2,zone3,zone4层做和上述相同步骤。（本例做出的2、3、4层图形如下）zone2层相建模如下：zone3层相建模如下：23451zone4层相建模如下：这样相模型就做完了。十、对连续数据进行分析：双击流程图中的出现下图：分别对perm,por,sw等连续数据进行如下操作：用鼠标选中中的Perm，对Perm进行操作。21３操作：分别对此窗口中的、、左击，然后依次点击加入到右侧窗口。出现下图：先用鼠标选中（反白）进行如下操作：依次左击1、2、3、4、5、6、7对不做任何操作（图形如下）再用鼠标选中（反白）进行如下操作：依次左击1、2、3、4、5、61234567再用鼠标左击Variograms,出现下图：进行perm变差函数的分析：Perm的主方向调变程图如下（前面已经对Vsh进行了分析，perm变差分析与Vsh变差分析的步骤相似），分别调好主方向，次方向，和垂向的变差。（下面选主方向变差图，次方向和垂向省略）再选中permzone1的砂体如下图进行操作（步骤重复以上）：142365解释1.选择分析对象,是经过离散化的井点,还是未离散化的测井曲线,或整个模型（默认即可）2.是否使用滤波功能;是否用相约束（打勾）3.分别按以下2个标签,进行相应的分析4.打开每个标签后,按键,刷新显示5.在进行transformation分析时,可以能够多种数据的处理,包括,输入截断，对数变换,正态分布变换等。6.例如要进行输入截断和正态分布转换,处理的对话框如下:十一、属性建模在建立油藏模型时，PETREL提供了确定性模拟和随机模拟两种算法：输入数据为粗化的井模型和趋势数据以及各种设置参数。当建立属性模型时，所有的网格都赋于数值.井数据和趋势数据分布在3D网格中.在建立模型之前,用户必须进行详细的数据分析,定义变差函数。随机性建模Petrel可以根据序贯高斯模拟建立随机岩石物理模型.这是用于产生多变量高斯域实现的直接算法。该算法可以产生局部变量.意思是如果基于相同的输入定义100个实现(使用不同的种子点)，可以得到100个不同的结果，他们都能和输入相匹配，但既然通过分布来输入，每个网格的数值的值会根据这个分布赋值.如果运行50-10