PETREL操作流程(详细)

PETREL操作流程1.前期数据准备地震数据体，断层线FAULT LINS OR 断层棍FAULT STICKS，FAULTPOL YGONS，数字化的等值线。

工区内各井的坐标，顶深，海拔，底深（完钻井深），东西偏移，方位角，倾角，砂岩分层数据，砂层等厚图，测井曲线（公制单位），单井相，各层沉积相图，砂岩顶面构造图，单井岩性划分，测井解释成果表，含油面积图。

（在编辑数据的过程中，命名文件时最好数据文件名都和井名一致）2.数据加载①加载井口数据（WELL HEADERS）WELL_NAME X Y KB TOP BOTTOM SYMBOL井名X坐标Y坐标海拔顶深底深（完钻）井的类型②加载井斜数据（WELL PATH）第一种数据格式MD TVD DX DY AZIM INCL斜深垂深东西偏移南北偏移方位角倾角第二种数据格式MD INCL AZIM第三种数据格式TVD DX DY（单井用WELL LOGS，多井加井斜可用PRODUCTION LOGS）③加载分层数据（WELL TOPS）（包括断点数据）MD WELLPOINT 层名WELL NAME-1500 HORIZON Nm31 NP1-1600 FAULT Nm32 NP1以WELL TOPS加载之后删除系统的缺省项，新建4项，对应输入数据的列，名称进行编辑，Sub-sea Z values must be negative!(低于海平面的Z值都为负)，该选项在编辑时不要选中④加载测井曲线（WELL LOGS）LAS格式文件MD RESIS AC SP GR曲线采用0.125m的点数据（1m8个点数据），注意有的曲线单位要由英制转换为公制，如：AC 英制单位μs/in要换成工制单位μs/m，再用转换程序转换为LAS格式文件进行输入，以提高数据的加载速度。

如果有孔渗饱数据，按相同格式依次排列即可。

在/INPUT DATA中设置数据的排列顺序，曲线内容较多，系统缺省项只有MD，所以要用SPECIFY TO BE LOADED定义新的曲线，对应加载数据的列数，名称和属性进行编辑。

如：DEPT 1 MEASURED DEPTH DEPTSP 4 SPONTANEOUS POTEN SP在/UNITS中可设置输入数据的单位在/SETTING 中设置参考高度：矫正海拔OR平均海平面OR其他⑤加载单井相（WELL LOGS）PRN格式文件MD（斜深）FACE岩性代码在曲线上划分不同层所属的相，划相时读取每层的顶值，不同的相取不同的代码，如：河道1，天然堤2，决口扇 3 …..以WELL LOGS加载之后，使用默认的MD，使用类似曲线的添加方式，属性取岩性LITHOLOGIES。

备注点沙坝BAR POINT天然堤NATURAL LEVEE河道CHANNEL决口扇CREV ASSE SPLAY泛滥平原BACKGROUND FLOODPLAIN河道沙坝CHANNEL BAR在加载完数据之后，GLOBAL WELL LOGS 中会出现FACE项，在/SETTING中设置每个相代码所代表的名称，颜色和图案。

⑥加载测井解释成果表（WELL LOGS）PRN格式文件TOP BOTTOM EXPLAIN WELL NAME顶深底深解释结果井名加载方式和加相一样，将油层，差油层，气层，水层等进行编码，在加载完数据之后，GLOBAL WELL LOGS 中会出现EXPLAIN项，在/SETTING 中设置每个代码所代表的名称，颜色和图案。

※使用甲方给定的最后一次测井解释成果表⑦加载构造图X坐标Y坐标类型Z值加载之后可选择以点或线输入，然后设置列数⑧加载断层POLYGON（闭合）数据X Y 类型（Z）⑨加载地震数据3.构造模型STRACTURAL MODELINGMAKE/EDIT SURFACE 成面选择INPUT 数据中加载进来的每个层的等值线或等厚线或POLYGONS，用户自己根据对应的层名命名，如果输入的有边界，就用边界，没输入边界的就点击选择自动生成边界。

在GEOMETRY中，几何网格模型选择自动模式，设置GRID INCREMENT，如果输入的有边界，就不选BOUNDARY。

在ALGORITHM中选择成面所用的算法，不同的面用不同的方法，常用的有构造面，等厚图，相图，在对所研究区域的地质情况有充分了解的时候可以在不同的算法里进行不同的设置。

在做相图的时候可选择协克里金方法或者模拟算法。

DEFINE MODEL新建一个模型，用户自己命名。

FAULT MODELING断层模型在3D窗口显示断层POLYGON，用户输入的断层没有Z值，要给断层赋Z值：双击断层POLYGON，在/CALCULATION下用公式A=SURFACE，SURFACE选择该断层所对应的面。

编辑断层POLYGON数据，例如，在3D窗口显示编辑好的1号断层POLYGON，该断层在每个层位的POLYGON都显示出来，然后CREATE NEW FAULT FROM POLYGONS，然后根据断层的掉向和组合方式对新建的断层进行编辑。

新生成的断层可根据用户需要，自己定义编号，以便管理。

PILLAR GRIDING生成三维框架编辑完断层之后，系统自动弹出2D窗口，显示出等值线和断层线，用户自己绘制一个边界，（如要设定断层走向要加趋势线），点击CREATE NEW 创建新3D GRID，设定网格的平面步长INCREMENT，并且选择断层处的网格编辑方式。

其他设置按用需要进行设置。

APPIY，观察网格分布，不符合要求再次调整，调整好后，按OK完成。

MAKE HORIZONS在HORIZON中根据用户分层添加层数，添加对应层位的分层，和SURFACE，在FAULT中设置断层的距离，一般情况下，不同的层有不同的断层，在列表中选择断层，某一层中没有的断层，就不激活。

设置完所有的层位，在 /SETTING 中完成相关设置，点击完成。

MAKE ZONES对每个层进行编辑，在隔层存在的情况下，每层有顶底两个SURFACE，再加上该层的分层，如有其他需要，在 SETTING和WELL ADJUSTMENT 中再做设置，如不整合设置和光滑设置。

LAYERING在LAYERING中，设置每个层在纵向上的分层，根据甲方要求划分纵向分层，隔层用1层就行了。

4.相建模㈠SCALE UP WELL LOGS 离散化测井曲线选取工区内有单井相数据的井，对单井相数据进行离散化。

㈡DATA ANALYSIS 数据分析①相数据分析PROPORTION（比例分析）原井数据会形成直方图，点击FIT ACTIVE/ALL CURVE（S）TO HISTOGRAME，观察在比例分布里是否有异常值，如有，则要进行更改THICKNESS（厚度分析）观察各相的厚度比例。

PROBABILITY（概率分析）V ARIOGRAMS（变差函数分析）根据数据分布的相关性，调节主，次和垂向变程，调节物源方向，选择分布模型。

河道中，曲流河的主变程一般范围在1200-1500m，辫状河主变程一般在900-1200m左右。

（经验值）一般情况下，曲流河的变程比辨状河的稍大，次变程一般去主变程的1/3左右垂向上的变程一般取河道沉积的单砂体厚度？数据分析中的地震属性体是起约束作用？变差函数是区域化变量空间变异性的一种度量，反映了空间变异程度随距离而变化的特征。

强调三维空间上的数据构形，从而可定量描述区域化变量的空间相关性，使克里金技术以及随机模拟的一个重要工具。

变程(Range) ：指区域化变量在空间上具有相关性的范围。

在变程范围之内，数据具有相关性；而在变程之外，数据之间互不相关，即在变程以外的观测值不对估计结果产生影响。

块金值(Nugget) ：变差函数如果在原点间断，在地质统计学中称为“块金效应”，表现为在很短的距离内有较大的空间变异性，无论h多小，两个随机变量都不相关。

它可以由测量误差引起，也可以来自矿化现象的微观变异性。

在数学上，块金值c0相当于变量纯随机性的部分。

基台值(Sill)：代表变量在空间上的总变异性大小。

即为变差函数在h大于变程时的值，为块金值c0和拱高cc之和。

拱高为在取得有效数据的尺度上，可观测得到的变异性幅度大小。

当块金值等于0时，基台值即为拱高。

已知井位数据，为了计算某个方向的实验变差函数，通常计算该方向上的若干不同距离的实验变差函数值。

此时选取某个距离为计算实验变差函数的基本距离L，称之为步长，分别计算L, 2L,3L,…..mL距离的实验变差函数值，此时就可以得到m个实验变差函数点，m称之为步长个数。

单个步长可选为指定方位上的平均井距，步长数目15--25个。

实际情况下是不可能在精确的某个方向上，或精确的某个步长上能获得需要的点对数目来计算实验变差函数值，此时。

一个变通的方法就是给方向一个容许的范围，称之为角度容限，同时，可以给每个步长一个容许范围，称之为步长容限，只要点对记录落入到该容限内就可以认为该点对参与计算。

这样就可以计算出某个方向相应特定步长的实验变差函数。

步长容限可选为1/2该方位上的井距，角度容限可选为π/8，值得注意的是，随着步长增加，虽然有时点对符合方向容限与步长容限，但是偏差将会增大，为此需要用偏离方向主线的一个固定宽度的带子来限制，使超出该范围的点对不参与计算，这样有利于得到合理实验变差函数计算值，带宽可选为2倍井距。

在拟合过程中，对于短距离的实验变差函数值应予以重视，是拟合模型形状尽量照顾到这些点。

因为这些点所提供的信息对空间变异型及预测结果有重要影响。

另外注意大多数情况下，由于单个步长一般不会太小，从而对选取块金常数不能提供精确的信息。

通常的做法是，若研究目标为区域上的物性参数变化情况，那么小的块金常数会告诉我们该物性参数具有很好的连续性。

因此多数情况下，可置块金常数为0。

变程的选取依据是适当步长大到何值时，实验变差函数值基本稳定在某个值范围邻近摆动，这个值就是拱高，而相应的步长就是变程。

数据分析中提供3种分布模型球状模型接近原点处，变差函数呈线性形状，在变程处达到基台值。

原点处变差函数的切线在变程的2/3处与基台值相交。

指数模型变差函数渐近地逼近基台值。

在实际变程处，变差函数为0.95c。

模型在原点处为直线高斯模型变差函数渐近地逼近基台值。

在实际变程处，变差函数为0.95c。

模型在原点处为抛物线。

㈢FACIES MODELING 相建模选择数据分析的相，分ZONES，选择算法，一般用序贯指示模拟算法，利用数据分析的结果，并加上每层的沉积相分布图加以约束，相关系数取值？几种不同的算法：克里金方法是一种实用的、有效的插值方法，它优于三角剖分法、距离反比加权法，它不仅考虑到被估点位置与已知数据位置的相互关系，而且还考虑到已知点位置之间的相互关系，因此估值精度相对较高。