PETREL操作流程
1.前期数据准备
地震数据体，断层线FAULT LINS OR 断层棍FAULT STICKS，FAULT
POL 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 DEPT
SP 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 MODELING
MAKE/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，选择算法，一般用序贯指示模拟算法，利用数据分析的结果，并加上每层的沉积相分
布图加以约束，相关系数取值？
几种不同的算法：
克里金方法是一种实用的、有效的插值方法，它优于三角剖分法、距离反比加权法，它不仅考虑到被估点位置与已知数据位置的相互关系，而且还考虑到已知点位置之间的相互关系，因此估值精度相对较高。
指示克里金插值-----是对原始数据的指示变换值进行克里金估计。所谓指示变换是指按照一定的原则将原始数据变换为0或1的过程。属于非参数统计方法，不同于其他克里金方法（通过参数的均值和方差来估值），它是以概率形式考虑特异值得存在，在不舍弃特异值数据的条件下进行有效的空间估计。
序贯模拟的总体思路是沿着随机路径序贯求取各节点的累积条件分布函数ccdf，并从ccdf中提取模拟值。其中用于求取ccdf的条件数据不仅包括原始的样品点，还包括已模拟好的点。概算法的主要目的是充分利用更多的条件数据来恢复变量的空间相互性。
序贯指示模拟------ 无需假设原始样本服从正态分布，而是通过给出一系列的门槛值，估计某一类型变量低于某一门槛值的概率，以此确定随机变量的分布。实际上是应用指示克里金求取ccdf的序贯模拟方法，主要特点是变量的指示变换、指示克里金和序贯模拟算法。
5.属性建模
㈠SCALE UP WELL LOGS离散化
对属性体数据进行离散化，一般建模过程是
PORE—PERM—SOIL——NTG
㈡DATA ANALYSIS 数据分析
选择要分析的属性体数据，分层，分相分析
TRANSFORMATIONS （数据转换分析）
选择属性体后，会出现数据的分布图，观察数据的分布，有异常值，在输入削截和输出削截中调节数值范围控制数据的分布
输入削截的数据不参与运算，输出削截的数据参与运算，而不向外输出？
其他参数的调节
V ARIOGRAMS（变差函数分析）
和相的变差函数分析一样，调节主次，垂向变程，调节物源方向，选择分布模型，拟合分布曲线，便于下步操作。
㈢PETROPHYSICAL MODELING 属性建模
选择要进行建立模型的属性体数据，选择ZONES，在单ZONE 里分相处理，例如对于Ng3-1，存在泥岩，砂岩和差砂岩三个相，，用序贯指示模拟的方法计算，采用数据分析的结果。在协-克里金方法下设置相关属性体或相关平面。其他设置
6.储量计算
MAKE CONTACTS 建立油水界面，油气界面。
划分界面时，要根据不同的断块，不同的层进行划分，每个层的界面应该都不一样，每个断块的也是一样。
VOLUME CALCULATION 储量计算操作界面
选择所在油藏的含有相，以南堡为例，有油气水三相，因此在前期建立的就应有油水和油气界面，
GEN.PROPS 中选择相关属性体NTG/POR，并设置相关系数。
OIL.PROPS中设置含油饱和度，含气饱和度（三相时设置），并且设置油的体积系数和溶解气油比(solution gas_oil radio)，并且设置采收率系数？
GAS.PROPS 在有气相存在时，设置气的体积系数和溶解油气比(vapourized oil_gas radio)，设置采收率系数？
RESULT/OUTPUT 中进行输出内容设置，在REPORT SETTING 中可以设置输出报告的内容以及输出结果的单位。/FACIES是怎么应用的？/BOUNDARIES 中可设置工区边界？井边界？
几个相关公式：
BulkV olume
V N NET=BulkV olume * N/G
V POR=NetV olume * Porosity
HCPV o= porvolume * So
HCPVg= porvolume * Sg
STOILP=HCPV o/Bo + (HCPVg/Bg) * Rv
GILP=HCPVg/Bg + (HCPV o/Bo) * Rs
Recoverable oil=STOILP * RecFo(oil recoverable factor)
Recoverable gas = GILP * RecFg(gas recoverable factor)
经过储量计算之后，得出一分报告，在报告中，会给出软件名，作者，制作时间，项目名称，模型名称以及所含相；作者的Bo，Bg,Rs,Rv,RecFo,RecFg,So,Sw,Sg等参数设置都会列表显示出来。列表显示上面每个公式所得到的值，最后得到可采储量。而后，分层ZONES，分断块SEGMENTS，显示结果。
7.井设计及相关预测
WELL DESIGN
画一条POLYGAN。新成一个剖面，显示PORE模型，ADD NEW POINTS，在预计井所过的位置加点设计轨迹。
在完成设计的井轨迹之后，在INPUT/WELLS中多出PROPSED WELLS ，打开，双击选择MAKE LOGS TAB，就会生成预测曲线；选择REPORTTAB/ICONNIZE POINT AS ：HORIZON，选择网格体，按EXECUTE就可以整理出分层数据，同时出一份报告。
8.出图及相关报告所需资料
绘制平面图MAP WINDOW
绘制截面图INTERSECTION WINDOW
FUNCTION WINDOW
蒙太奇绘图NEW PLOT WINDOW
9.模型粗化及输出
. 模型完成后，为了达到数值模拟的数据要求，保证网格总量能在数模中运算，我们要对网格进行粗化后输出（GENERIC ECLIPSE STYLE(ASCII) GRID GEOMETRY）；输出网格中包含粗化后的孔隙度模型，渗透率模型（三向），饱和度模型（ECLIPSE中需要的是含水饱和度Sw），NTG模型；最后还要输出模型中的井轨迹数据（ECLIPSE WELL COMPLATION DATA），射孔数据（ECLIPSE WELL CONNECTION DATA）。
下面是粗化网格及属性体的过程。
PILLAR GRIDDING，创建一个新的3D GRID，设定要粗化到的网格步长（INCREMENT），例如由粗化到100*100，设置断层
处的网格形式。
SCALE UP ZONES/SELECT GRID，选择要粗化的已建立好的精细地质模型，设置断层重采样模式，设置粗化过程中的网格变化形态，设置LAYING 的建立模式，粗化纵向分层（一般情况下纵向的分层是不变的），如果要因为网格数太大而粗化，就把原来的分层数压缩到1/2，例如，原来Ng2-1分10层，粗化到5层。
问题：隔层的粗化
ZONE MAPPING，参照对比粗化前后的纵向分层。
SCALE UP PROPERTIES，选择要粗化的模型，首先是孔隙度模型，选择取样模式及粗化精度，选择粗化的平均方法。并且选择各属性体的相关属性，进行控制。要注意的是粗化PERM时，要点击TENSOR，选择三向粗化（粗化使用的方法？），其他相关设置要再摸索，现在不是很清楚。
粗化完网格和属性体后，保存3D GRID。网格输出时附加属性体。射孔和完井数据要单独输出。
属性建模： 一、相模型的建立： 1、测井曲线离散化 双击：Process ——Proerty modelding——Scall up well logs； 弹出对话框：
在Select里选择需要离散化的相曲线数据facies(input到wells的沉积相数据)，点击all可以对需要离散的井进行选择，剔除没有曲线或者曲线数据不正确的井）。 在相模型建立时：Average选择“most of”、method选择“Simple”。单击“Apply”或“OK”确定。完成沉积相数据的离散化，离散化后，沉积相数据赋给井轨迹所通过的网格。离散化后models里的properties里新增了沉积相属性“facies”，可在3D视图里进行查看。
2、沉积相模型建立; 双击：Process ——Proerty modelding——Facies modeling。 弹出对话框：
对话框右上角选择离散化后的沉积相数据，依次选择各小层（zone）进行属性控制；点击解锁进行编辑控制。 目前的沉积相建模算法很多；通常，纵向上细分网格后用序贯高斯的算法，纵向上未细分用经典算法（此处的“纵向细分“是指layering里把zone细分为不同个数的网格。 ⑴、序贯高斯的算法； “Method for zone /facie”选项单击下拉菜单， 选择序贯高斯算法：“Sequential indicator simula”，在左侧选择该小层所以相类型（可从 左侧出现的百分比统计中看出）单击箭头，相 类型移动到右侧。
下侧空白区域新增两个选项卡“Variogram”，“Fraction”，点击按钮，弹 出对话框：
Petrel建模常用术语 Petrel引入了一些新的术语和公式表达式，现简要地解释如下。 3D Grid –是一个用来描述三维地质模型的由水平线和垂直线组成的网格。Petrel中应用了角点三维网格技术。 Artificial method –用于make surface进程中，意思是在建surface 时不用任何输入数据。 Attribute map –是一张地震属性图。可以从地震体中通过提取穿过某一层面的属性值来获得（分两种：一种是从某一表面开始的一定偏移量内的平均属性；另一是两个面之间的平均属性）。 Automatic legend - 一个预先确定好的用于显示窗口中目标体色标的模板 Bitmap image - 输入的位图，例如BMP和JPG格式的位图文件，它们都可以在UTM（通用横轴墨卡托投影坐标系）中显示出来。 Bulk Volume - 总的岩石体积 Cell Volume –三维网格中单位网格的体积。 Connected Volume –在离散的3D属性中计算相连体积的进程，可用来查找相连的河道。 Contact Level –油水或油气界面，通常是一个固定深度值。Contact Set –由用户自己定义的一组接触界面，用作储量计算的输入值，也可用作显示使用。 Cropping –通过定义主线、联络线和时间范围，创建真实的地震体。Crossline intersection –垂直于主测线方向的垂向地震切面。
Cross plot –两个或两个以上的数据相互间形成的交会图（也叫做scatter plot（散点图））。 Datum –在测定海拔时用到的一个固定深度、时间值或是一个层面。Depth Contours –层面的等高线，描述相同的深度或时间值。Depth Conversion –将Z值在深度域和时间域间相互转换。 Depth panel –井上的垂向深度标尺。 Display Window –用于显示模型的窗口，分为二维、三维两种类型。Dongle –硬件加密锁（hardware key），也叫做软件防盗锁(software protection key)，它控制着软件模块的使用时间。 Drainage Area –流域，指的是可能产生烃的区域。 Erosion Line –剥蚀线，用于定义层面间的相互削截。 Fault Center Line –3D网格中用于连接断层Pillar中点的线。 Fault Modeling - 在三维空间骨架中建立断面的过程。其第一步就是建立Key Pillar（主要断层柱子）。 Fault Polygon –断层平面和层面间的交线。 Fault Stick (fault dip line) –描述断层的线，通常是贯穿顶部和底部。Fluid Constants （流体常量）–地层体积系数，油Bo，气Bg。GOR：气油比。严格讲采收率不是流体常量，但在Petrel中将其列入了储量计算的流体常量菜单中。 Formation Volume Factor –地层体积系数。地表情况下的烃体积与油藏中的体积之比（油和气的分别为Bo和Bg）。 Function Bar –在微软术语中叫作工具栏（toolbar）。不同的进程中，
目录 1.加载数据 (4) 1.1 井位数据 (4) 1.2 井斜数据 (4) 1.3 测井曲线加载 (5) 1.4 分层数据加载 (9) 1.5 测井解释成果加载 (13) 1.6 断层加载 (14) 1.7 地震数据加载 (15) 1.8 制作地震子体 (17) 1.9 地震解释 (23) 2.Make surface (32) 2.1 圈定边界 (32) 2.2 做面 (32) 3.调节断层 (37) 3.1 双击加载的断层.TXT文件 (37) 3.2 删掉断层一盘 (37) 3.3 将断层赋给一个面 (38) 4.断层模型 (39) 4.1 初步调整 (39) 4.2 pillar Giidding (45) 4.3 Make horizons (47) 4.4 Make zones (50)
4.5 调节断层上下盘 (51) 4.6 补缺口/horizon (52) 4.7 做垂向网格/layering (56) 5.砂孔建模 (58) 5.1砂体模型（确定性） (58) 5.2砂体模型（指示建模） (66) 5.3夹层模型 (66) 6.沉积相模型—确定性 (70) 6.1 创建沉积相模型 (70) 6.2 相图加载 (71) 6.3 数字化位图 (72) 6.4 生成相多边形曲面/对每个相做surface (74) 6.5 生成相分布曲面 (76) 6.6 相建模 (77) 7.沉积相建模—随机性 (79) 7.1 PPT--序贯指示 (79) 7.2 阳光石油相模型建立--序贯指示 (80) 7.3 沉积相模型建立—聚类分析方法 (86) 8.沉积相相控属性建模 (103) 8.1 孔隙度模 (103) 8.2 渗透率模拟 (112) 8.3 含油饱和度模拟 (118) 9.计算储量 (126) 10.模型粗化 (134) 11 离散化测井曲线 (138)
一个综合Barnett页岩气藏的工作流程的建模与仿真 C. Du, SPE, X. Zhang, SPE, B. Melton, D. Fullilove, B. Suliman, SPE, S. Gowelly, SPE, D. Grant, SPE,J. Le Calvez, SPE, Schlumberger 这篇文章是准备在2009年5月31日至6月3号在哥伦比亚卡塔赫纳举行的拉丁美洲和加勒比石油工程会议上作为（会议）报告用的。 这篇文章根据作者所提出的包含在摘要中的信息被程序委员会选择出来作为一篇会议上的报告。石油工程师协会没有对本文的内容进行检查，需要作者自己进行校正。该文章不反映石油工程师协会、工作人员和会员的任何态度。电子复制品、分发品，没有经过石油工程师协会的书面同意，任何文件的一部分的存储都是禁止的。允许复制的（范围）限定在不超过300字的摘要，插图可能不能被复制。（被）复制印刷的摘要必须包含显眼的石油工程协会的版权信息。 摘要 密西西比Barnett页岩储层开辟了美国的天然气生产的新时代。做的许多油藏描述方面的努力和完成的一些实际生产，以帮助更加深刻的了解Barnett页岩储层。钻孔图像解译，钻井诱导产生的裂缝和连通的/闭合的裂缝，揭示（地层）应力方向，断层的形貌和方向等解释结果指导水平井设计，控制水力压裂方向和强度。常规测井和岩心分析已经用于对岩相的分类和评价油层物性和地球物理性质，以用于井的定位和储量计算。地震调查不仅用于水平层位和断层的解释，也用于3D物性的评价分析，如岩相分布，离散裂隙网络和应力场。在实际施工方面，多钻较长的水平井和进行大规模的多级、多层次水力压裂处理。大量的井的钻探和水力压裂都被广泛实施。微震（MS）对评价水力压裂所波及到的油藏的体积和压裂产生的断裂强度估算的起到重要作用。 尽管在这个方面巨大的努力和进展，但现有的文献中仍然缺乏一个系统
主要模块介绍 一、数据准备 本实例中的数据整理如下： wellhead井位坐标文件 jinghao X Y kb topdepth bottomdepth X21-233973816364714261433.0821502195 X21-243974070364716291433.082156.12193.1 X21-253974257364718491433.082154.42190.4 X21-263974480364720961436.52154.82189.8 X22-193972535364705161407.562120.32152.3 X22-203972803364707951417.462139.12165.1 X22-213973010364710401379.72102.62135.6 welltop分层文件 X Y hb wellpoint surface jinghao 397381636471426-716.92Horizon c811X21-23 397381636471426-724.92Horizon c8121X21-23 397381636471426-735.92Horizon c8122X21-23 397381636471426-755.92Horizon c813X21-23 397381636471426-761.92Horizon c821X21-23 397407036471629-723.02Horizon c811X21-24 397407036471629-731.02Horizon c8121X21-24 397407036471629-742.02Horizon c8122X21-24 397407036471629-754.02Horizon c813X21-24 397407036471629-760.02Horizon c821X21-24 测井文件准备 DEPTH PERM_K POR_K SW_K VSH_K NTG 2140.1250.00590100 2140.250.00590101 2140.3750.00590100 2140.50.00590010 二、数据输入 1 输入WellHeader（井位坐标文件）
Petrel建模流程 一、数据预备 二、数据输入 三、Pillar gridding 四、Make horizon 五、Laying 六、建立几何模型 七、离散化测井曲线 八、对Vsh数据进行分析 九、相建模 十、对连续数据进行分析 十一、属性建模 十二、网格粗化及属性粗化的操作 十三、储量运算 十四、产生STOIIP （烃体积密度分布图） 十五、输出数模所需要的文件
要紧模块介绍 一、数据预备 本实例中的数据整理如下： wellhead井位坐标文件 jinghao X Y kb topdepth bottomdepth X21-233973816364714261433.0821502195 X21-243974070364716291433.082156.12193.1 X21-253974257364718491433.082154.42190.4 X21-263974480364720961436.52154.82189.8 X22-193972535364705161407.562120.32152.3 X22-203972803364707951417.462139.12165.1 X22-213973010364710401379.72102.62135.6 welltop分层文件 X Y hb wellpoint surface jinghao 397381636471426-716.92Horizon c811X21-23 397381636471426-724.92Horizon c8121X21-23 397381636471426-735.92Horizon c8122X21-23 397381636471426-755.92Horizon c813X21-23 397381636471426-761.92Horizon c821X21-23 397407036471629-723.02Horizon c811X21-24 397407036471629-731.02Horizon c8121X21-24 397407036471629-742.02Horizon c8122X21-24 397407036471629-754.02Horizon c813X21-24 397407036471629-760.02Horizon c821X21-24 测井文件预备 DEPTH PERM_K POR_K SW_K VSH_K NTG 2140.1250.00590100 2140.250.0059010 1 2140.3750.00590100 2140.50.005900 1 0 二、数据输入 1 输入Well Header（井位坐标文件） 右键点击输入Well Header: 文件类型里选：Well heads (*.*)
储层建模的步骤 目前普遍的认识是，储层建模应分为油藏构造建模、沉积（微）相建模和油藏属性建模三步完成。构造模型反应储层的空间格架，在建立储层属性的空间分布之前，应进行构造建模。由于沉积相对储层物性有决定性的作用，油藏属性建模多采用相控建模，即先建立沉积微相模型，然后以此为基础进行油藏属性建模。 张天渠油田长2油藏的储层地质模型是以测井资料为基础资料，采用确定性建模的储层建模方法建立的。储层建模的整个过程包括4个主要环节，即数据准备、构造建模、油藏属性建模、模型的应用。 一、数据准备与预处理 1．数据准备 一般从数据来源看，建模数据包括岩心、测井、地震、试井、开发动态等方面的数据。从建模的内容来看，基本数据包括以下四类： ①坐标数据：包括井位坐标、地震测网坐标等； ②分层数据：各井的油组、砂组、小层、砂体划分对比数据；地震解释层面数据； ③断层数据：断层位置、断点、断距等； ④储层数据：储层数据是储层建模中最重要的数据。包括井眼储层数据、地震储层数据和试井数据。井眼数据为岩心和测井解释数据，包括井内相、砂体、隔夹层、孔隙度、渗透率、含油饱和度等数据，这是储层建模的硬数据。 对不同来源的数据进行质量检查是储层建模中十分重要的环节。为了提高储层建模的精度，必须尽量保证用于建模的原始数据特别是硬数据的准确性。因此，必须对数据进行全面的质量检查，如检查岩心分析的孔渗参数的奇异值是否符合地质实际，测井解释的孔渗饱是否正确等等。 建模过程中能被储层建模软件所采用的资料来源于这些基础资料，但它们有特殊的格式要求，需要转换成不同格式要求的文本文件才能以正确的格式导入到Petrel软件中。从文件类型上来看，它们包括井头文件（Well head）、井斜文件或井轨迹文件（Well deviation）和测井数据文件（Well log）。它们的格式和作用分别如下： ①井头文件： 文件内容包括井名、井位坐标（X、Y）、地面补心海拔（补心高与地面海拔之和）以及目标井段深度（井段顶部深度和测井段底部深度）。井头文件主要用来确定油藏中的井数、井位和各井的研究层段等井信息。
Contents 一、数据准备 二、数据输入 三、Pillar gridding 四、Make horizon 五、Laying 六、建立几何模型 七、离散化测井曲线 八、对Vsh数据进行分析 九、相建模 十、对连续数据进行分析 十一、属性建模 十二、网格粗化及属性粗化的操作 十三、储量计算 十四、产生STOIIP （烃体积密度分布图）十五、输出数模所需要的文件
主要模块介绍 一、数据准备 本实例中的数据整理如下： wellhead井位坐标文件 jinghao X Y kb topdepth bottomdepth X21-233973816364714261433.0821502195 X21-243974070364716291433.082156.12193.1 X21-253974257364718491433.082154.42190.4 X21-263974480364720961436.52154.82189.8 X22-193972535364705161407.562120.32152.3 X22-203972803364707951417.462139.12165.1 X22-213973010364710401379.72102.62135.6 welltop分层文件 X Y hb wellpoint surface jinghao 397381636471426-716.92Horizon c811X21-23 397381636471426-724.92Horizon c8121X21-23 397381636471426-735.92Horizon c8122X21-23 397381636471426-755.92Horizon c813X21-23 397381636471426-761.92Horizon c821X21-23 397407036471629-723.02Horizon c811X21-24 397407036471629-731.02Horizon c8121X21-24 397407036471629-742.02Horizon c8122X21-24 397407036471629-754.02Horizon c813X21-24 397407036471629-760.02Horizon c821X21-24 测井文件准备 DEPTH PERM_K POR_K SW_K VSH_K NTG 2140.1250.00590100 2140.250.0059010 1 2140.3750.00590100 2140.50.005900 1 0 二、数据输入 1 输入Well Header（井位坐标文件） 右键点击输入Well Header: 文件类型里选：Well heads (*.*)
P e t r e l中的属性建模 流程简介
属性建模： 一、相模型的建立： 1、测井曲线离散化 双击：Process ——Proerty modelding——Scall up well logs； 弹出对话框：
在Select里选择需要离散化的相曲线数据 facies(input到wells的沉积相数据)，点击all可以对需要离散的井进行选择，剔除没有曲线或者曲线数据不正确的井）。 在相模型建立时：Average选择“most of”、method选择“Simple”。单击“Apply”或“OK”确定。完成沉积相数据的离散化，离散化后，沉积相数据赋给井轨迹所通过的网格。离散化后models里的properties里新增了沉积相属性“facies”，可在3D视图里进行查看。 2、沉积相模型建立; 双击：Process ——Proerty modelding——Facies modeling。
弹出对话框： 对话框右上角选择离散化后的沉积相数据，依次选择各小层（zone）进行属性控制；点击解锁进行编辑控制。
目前的沉积相建模算法很多；通常，纵向上细分网格后用序贯高斯的算法，纵向上未细分用经典算法（此处的“纵向细分“是指layering里把zone细分为不同个数的网格。 ⑴、序贯高斯的算法； “Method for zone /facie”选项单击下拉菜单，选择序贯 高斯算法：“Sequential indicator simula”，在左侧选择 该小层所以相类型（可从左侧出现的百分比统计中看 出）单击箭头，相类型移动到右侧。
1 Fault Throw Background The offset along a fault can be measured by using the Measure Distance tool on faults that have been filtered for a certain horizon. Measuring this distance at several places along your faults will give you an approximation of the vertical displacement of a horizon. For more accurate data, you can use this workflow that calculates the exact throw at a regular interval along all your faults. Result The workflow calculates the throw along given faults for a specific horizon. This will be outputted as a point set with a Throw attribute. Intermediate results of this workflow that will be kept are a folder containing fault sticks for all faults in the model, and a polygon and point set that outline the intersection between the faults / fault sticks and the horizon. Prerequisites For this workflow one needs to have a Model with faults and horizons, for which you want to calculate the throw. These will be dropped into the first few lines of the workflow as input data. The project already contains some Reference Data that will be used during the workflow. These can be copied over to your project, when also transferring the Workflow. Alternatively, you can create your own. The workflow calculates the throw as distance along the pillars of the 3D grid. One can check the settings for the fault throw point set that is created to get more insight into the displacement along the faults for a given horizon. Under statistics the range and average values can be viewed, and in the histogram tab one can have a closer look at the distribution of the fault throw. The output point set can be used to make a surface with a throw attribute, allowing you to have a clear visual display of the displacement. Workflow Specifics Workflow Specifics Petrel Version 2009.1 Domain Geology Modules Geoscience Core or Combined Core Expert Level Intermediate
PETREL操作流程 1.前期数据准备 地震数据体，断层线FAULT LINS OR 断层棍FAULT STICKS，FAULT POL 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 DEPT SP 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项，在
Petrel 软件篇——裂缝建模Fracture Modeling 对于裂缝的认知可以帮助我们更加充分的了解和预测油藏特征。根据所建立的精确的裂缝模型，我们可以 充分了解相邻网格的空间相关性。在模型中，每一条裂缝都可以用一个面表示，Petrel 将以离散性数据形 式来描述裂缝，并建立“离散裂缝模型”。在Petrel2007 版本中，Schlumberger 与业内的领军者Golder 联手，共同为油藏裂缝建模打造完美工作流程。 裂缝建模需要多步实现，涉及到油藏描述和模拟的方方面面。其主要的宗旨是基于地质概念，充分利用基底解释、断层和成像测井的裂缝知识、通过类比野外露头建立的裂缝概念模型、可预测裂缝成因的地震属性等等，并将这些资料转换成裂缝强度等参数，建立三维的裂缝模型。 裂缝建模流程： 1．输入数据， a) 质量控制，并且显示由测井资料解释出的裂缝、由成像测井资料得到的 倾角和方位角； b) 井点资料，每个属性都可以描述裂缝的类型和质量； c) 产生蝌蚪图，用于显示倾角和方位角； d) 打开一个显示玫瑰图的窗口，将这些井点数据投放在玫瑰图上，估计裂 缝的类型； 2．如何生成蝌蚪图 a) 在Well Section上显示具有裂缝资料的井； b) 右键击一口井，选中Insert Points/tadpole panel； c) Select fracture values from points.选择裂缝值； d) 点击Apply，观察蝌蚪图； 3．数据分析 a) 产生新的点属性数据，通过Petrel计算器计算出裂缝相对于地层的倾角； b) 产生裂缝密度曲线； c) 利用玫瑰图中的提供的一些工具给裂缝单元赋值； 4．建立裂缝模型