petrel 技巧

Petrel中Polygon编辑技巧1）首先在2D 窗口中显示出想要编辑的Polygon，在其上右键，选择Polygon编辑工具。

2）如要编辑Polygon中的点，可以使用Polygon editing 工具条中的第一个工具-Polygon editing，快捷键为E。

一旦激活此工具，Polygon会显示出其上的编辑点。

3）将鼠标点击2D 窗口中的Polygon上的点，此时激活的点显示为黄色，可以对要编辑的点进行拖拽。

4）如要对Polygon进行切断，则选择工具条的最后一个工具，Cut Polygon，此时鼠标变为类似刀锋的模式，用刀锋尖点击要剪断处，此时点击出会添加一个点，用上一步的工具去拖动此点，此时可以发现Polygon已断开。

5）如要在切断的Polygon上继续编辑或者画此Polygon，需要先激活画Polygon的工具-Add point to polygon工具，快捷键为A。

6）激活Add point 工具后，在所要编辑的Polygon上直接点击，查看切断的Polygon的首尾点，一般此时点击Petrel默认是从Polygon的末尾点继续来画，直接点击鼠标左键即可编辑。

7）如要从断开的Polygon的首点开始画，则需要将首点选中为黄色（可以使用第一步中的工具来选中），选中Add Point工具，按住Shift键，在靠近选中的黄色的点附近点击鼠标左键，即可继续以首点来画。

（注意此处，根据测试，按住Shift后点击鼠标左键时距离Polygon上点的位置越近，其越能编辑其附近的点。

）8）最后编辑完，如要闭合，则可以直接右键点击，在弹出的工具条上选择Close工具即可。

9）对于两条Polygon需要进行连接的，如果两条Polygon没有在同一个文件中，首先需要将两条Polygon合并到一个Polygon文件中。

此时可以使用Polygon Setting中的Operation中Polygon Operation中的Append工具来完成。

主要模块介绍一、数据准备本实例中的数据整理如下：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-23397381636471426-724.92Horizon c8121X21-23397381636471426-735.92Horizon c8122X21-23397381636471426-755.92Horizon c813X21-23397381636471426-761.92Horizon c821X21-23397407036471629-723.02Horizon c811X21-24397407036471629-731.02Horizon c8121X21-24397407036471629-742.02Horizon c8122X21-24397407036471629-754.02Horizon c813X21-24397407036471629-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 输入WellHeader（井位坐标文件）右键点击输入Well Header:文件类型里选：well heads(*.*)2 输入Well Tops（分层文件）:右键点击Well Tops文件夹并选择Import (on Selection)；文件类型里选：Petrel Well Tops (ASCII)3 输入输入Well Logs右键点击Wells文件夹，选择Import (on Selection)；文件类型：well logs(ASCII)input Data logs specify logs to be load加载per,perm,sw vash,ntg 等数据。

Petrel中统计特定深度测井成果的平均值在进行平面对比工作中，研究人员通常需要统计不同深度段的测井成果，这个深度不仅只在某个特定层位中，也可能是某特定深度段。

那么，下面就介绍如何在Petrel中统计特定深度段测井成果的平均值。

具体操作如下：1、在原有Project的基础上，新建well tops，名为well tops 1。

2、建立两个分层面，例如一个为深度-1900起始，一个深度为-2200止。

用make/edit surface来做，给定边界，geometry选automatic，剩下设置如图所示。

3、将两个做好的surface变成新的分层进入到well tops 1中进行管理。

• 激活新的well tops1• 双击wells，进入report标签下，将做好的两个层位按顺序给入，其他设置如图• 点击run，新的well tops1中就有这两个分层4、在well tops 1 的attributes中添加新的属性，算法选平均，以Por、Perm为例，这样就可以得到深度为1900-2200m这段中整个孔渗的平均值，取名avg 1900-2200。

5、如果想计算NTG。

右键测井曲线calculator，键入相应的门槛值即可。

6、还可以计算在1900-2200m深度段中NTG的储层厚度，在well tops 1 的attributes中添加新的属性储层厚度”thickness of NTG“，以计算Por>0.15&Perm>3 code值为1的储层厚度为例。

7、那么如果想计算NTG为1的储层中孔渗等属性的平均值，那么用孔渗NTG做交会图，建立符合门槛值的filter，然后右键filter转换为boolean log，并存放至global well logs下。

8、再次在well tops 1 的attributes中添加新的属性Por和SW. 算法选平均，我以Por Perm为例，但是要勾选boolean log作为过滤器，这样就可以得到1900-2200m深度段中且在NTG储层段中的孔渗段平均值，我取名1900-2200。

Petrel中如何精确定位图片坐标日常工作中，通常需要把一些图片格式的研究成果加到petrel中以方便后续的一系列项目研究。

要想把图片的信息精确的利用起来，准确定位图片是一个关键步骤。

现在介绍一种比较实用的在petrel中精确定位图片方法—通过图片井位信息定位图片位置的方式，下面是介绍的详细定位过程：第一步：鼠标放在图片上查看图片的像素为945*715，如下图1，然后在input标签下鼠标右键菜单选import file，找到相应路径下选上该图片，文件类型选bitmap image格式，open加载进来。

在input标签的最下面，可以看到加进来的图片，双击该图片到settings>settings下，按照图片像素提供的大小设置图片相应位置的三点坐标，如下图2坐标一(0,0),坐标二(0,715),坐标三(945,0),然后ok：打开一个2D 窗口，显示该图片，打开point editing工具，用add points to pointset工具在工区的井位M5、M2、M1上依次点三个点，如下图3：在input下新产生一个新文件points2，在points2文件上点右键选择菜单上的spreadsheet 打开，分别拷贝X，Y列的数据，粘贴到excel表中，如下图4，同时把图片上想要定位的井在petrel工区的wells的well manager里拷贝然后粘贴到excel里，分别拟合X坐标组、Y坐标组的数据点，可以看到都具有一个线性关系。

如下图5、图6。

按照图片需要放大的要求，需要拟合一个斜率和截距，首先拟合一个斜率k=(B3-B2)/(A3-A2)公式求出线性关系式的斜率k=13.7，然后按照线性方程X的横截距A=B3-k*A3=227156.8143，同理求纵截距B=B9-k*A9=8557049.8,然后基于上面已知的斜率和截距，用公式X0=A14*k+A, Y0=B14*k+B就得到X0、Y0的大地坐标，同理用相同的公式计算出（X0、Y1），（X1、Y0），如下表1：然后在petrel里双击该图片到settings>settings下，设置三点坐标如下图7：然后打开一个map窗口，投上该图片，再投上wells下定位的三口井坐标，可以看到图片上的井和实际井坐标是完全吻合的，说明图片已经按照一定的比例放大定位到正确的工区位置上了。

如何切割属性模型在Petrel中，有时为了数模的需要，常常不需要对全区模型都进行数模，只需要对部分模型进行数模，因此常常需要在Petrel中输出部分模型，下面为输出部分模型的几种方法:第一种:我们可以通过创建boundary来协助对属性模型进行切割。

Boundary用做控制切割模型的边界.一共可分成四个步骤。

1, 可先选取一个解释层位,最好在二维窗口操作,在显示窗口显示上一个模型层位,激活Process下的Utilities, 激活Make/edit polygons,对模型里感兴趣的区块用polygon圈出来并闭合polygon(如图1)。

2, 使用polygon来裁剪显示的构造层位(如图2)。

针对显示的某一构造层位鼠标右键双击打开Setting,在Operation下,找到Eliminate where文件夹,选择Eliminate outside, 然后在蓝色按钮处,将构造层位选进去。

3, 裁剪属性体先打开你想要裁剪的模型的某一属性的setting,然后在Operation键下打开Property operations文件夹,找到Property=surface(x,y)并选择之。

并将上一步裁剪的层位用蓝色箭头选择进去(如图3)。

4, 为裁剪出的小模型赋属性值。

点击剪切出来的小模型右键选择calculator。

在三维窗口显示小模型,将大模型的值通过公式赋值给小模型(如图4)。

第二种:1，准备属性模型。

2，第一种方法是简单选取部分网格。

鼠标左键双击或者右击打开要选取部分模型的属性模型的Settings,如图5所示。

在Output tab下，选中Use segment filter, Use zone filter, Use cell index limits，此三个过滤器可以任意组合使用，当选中Use cell index limits后，可以在下面的IJ选取栏内选取要显示的IJ方向的网格的范围。

Petrel中蚂蚁体的运算技巧
现在的Petrel用户中对蚂蚁体属性的使越来多，大部分都能取得较好的效果，下面是目前Petrel中常用的蚂蚁体的研究流程：
该流程基本涵盖了蚂蚁体生成的流程，但有时虽然各环节做的都很好了，但最终结果仍令人不太满意，实际上通过下面的办法可能对最终结果有进一步提高，用户在下面研究时可以试一下。

它的指导思想是，在现有蚂蚁体基础上再多次运算蚂蚁体，以期得到效果上的改变。

例子1：单次运算结果
1. 单次运算被动模式的蚂蚁体（Passive）得到不好的结果：信号非常差，非常差的连续性。

2. 单次运算主动式的蚂蚁体(aggressive)得到不好的结果：信号差，不好的联系性，很难理解断裂的分布。

例子2：两次运行结果
1. 两次运行：先被动后主动方式，获得较强的信号和比较好的连续性。

2. 两次运行：先主动后被动方式，获得更强的信号，更好的连续性，很多重要的特征显现出来。

例子3：连续3次运行结果
1. 3次运行：被动+主动+被动（极大地呈现了主断裂系统特征，清晰的图像，比较容易进行断面的提取）
2. 3次运行：主动+被动+主动方式（最强的结果，可能不太容易进行提取，但很好展示整个项目的断裂的分布）
如果能够遵循下面原则，蚂蚁体会得到较好的结果：
1. 数据预处理：通过Chaos或Variance属性来减少噪音和不连续性增强
2. 断层增强：通过Ant Tracking计算，把上面的重要的不连续性分解为片。

3. 断层提取：验证、编辑、把蚂蚁体结果合并为面。

4. 建立最终三维断层模型。

构造与岩性共同控制的油藏中如何设置多个油水界
面
对于构造和岩性双重控制且油水关系复杂的油藏，在petrel中往往比常规单一的油藏做油水界面要复杂难处理一些。

碰到类似的情况，可以采取下面的工作流来解决。

1、在2D窗口下，显示构造层面，选择合适的构造等值线间距，如果有确定的油水边界文件，可以输入到petrel中，如果没有这类数据，可以激活make/edit polygons 进程，按照分析的确定性成果，在构造图上圈取油水边界线及岩性控制的边界线(图1)。

双击make/edit surface进程，弹出界面下输入多边形文件，选取Artificial algorithms，赋对应的Z值，生成油水界面surface文件。

对于构造控制的断块，surface对应一个油水边界的海拔深度值；对于岩性控制的断块，surface对应岩性边界所在的最低构造线值，最后形成了多个油水界面surface文件，双击一个surface文件，弹出settings界面下<operations<surface-surface operations，用A="" union="" B,keep="" B操作，联合其它的surface，合并成一个surface 文件（图2）。

</operations
双击make contacts进程，弹出界面下输入此surface文件，就可以求取相应层位的油水界面了（图3）。

使用地质体雕刻和人工神经网络来描述盐丘下面将阐述如何利用Petrel的geobody interpretation 和neural network这两种方法来描述盐丘。

具体实施步骤如下：1．生成多种地震属性选择你认为对刻画盐丘这种属性有帮助的一些地震属性，比如：混沌体（Chaos），Envelop, RMS振幅，方差，PCA-方位角，主频等等。

然后分别查看这些地震属性体（图1）。

在图1中用到的是曲率体(vertical radius=14;inline/xlineradius=1;method=strike curvature;),局部构造倾角（local structural dip）,边缘探测体（edge detecting cube）混沌体和方差体来做。

2．地质体雕刻第一阶段使用Petrel的geobody interpretation这个功能，根据雕刻体需求个数设置相同的blob然后转换成散点（图2，图3）。

然后选取你想要的一个体来进行下面的分类过程。

3. 人工神经网络训练完成第二步以后，如果需要一个作为背景的点集，可以通过手工拾取一些点，或者通过拾取多边形或断层然后再转换成散点都可以。

然后给这些散点附上属性（图4）。

再定义一个分类（图5）：一旦完成你想定义的分类，你需要把所有的点集合并起来。

在setting>common operations>append point with attributes)。

这个过程需要把点集放在一个文件夹里。

（图6）结果可以通过三维窗口和Spreedsheet进行显示，如图7。

4. 分类在上述训练属性计算好之后，然后可以进入到人工神经网络分类，如图8所示。

可以使用督导的分类方式，并选用合适的模板，可以调整迭代次数，误差范围等来调整分类结果（如图9）。

5.结果显示通过上述步骤，在正确的操作步骤下会得到下面的结果，如图10所示：。

[Petrel使用技巧] Petrel蚂蚁体的问题拿一个实际数据体说说各个参数的影响力。

原始数据的分析：一般作构造解释的地震数据并不特别需要做保幅处理，所以有些适度的AGC是比较好的作AGC前后的剖面对比如下：Petrel工作流程中推荐做“Structural smoothing”，注意这个属性的参数范围：效果对比：不同的平滑参数的对比效果：原始数据、AGC数据、Structural smoothing数据的方差体结果的对比：同数据源，不同的蚂蚁步长“Ant step size”对比效果：时间切片的效果对比，顺序跟上图不一致：雕刻一下：提取Fault paches的参数当然同样意义非常模糊，但是也同样重要，其实原问题问的就是这个部分了：对于面积、倾角等作过滤非常重要：合并Patches：适度平滑：转解释数据：其实有些Patches被错误合并到一起也是导致解释成果转化错误的一个原因。

蚂蚁体不好做出效果的原因，我们把事情说穿就好——工作流程过长，参数过多互相牵制。

原始数据体——强化后处理——方差体-——蚂蚁体——断层提取——Patches编辑——生成解释数据。

这个长流程中任何一个环节其实都可能增加或者强化出我们不希望出现的东西。

而且蚂蚁体的目标就是强化方差体，尤其要解决方差体因为搜索时窗在大断层带来的阶梯效果和有些小断层可以被重新识别出来。

这两个主要目标本身就是矛盾的，你希望模糊大断层上的时窗效果势必需要比较大的蚂蚁步来过滤（Filter）细节，而要保留甚至强化小断层效果则肯定要保证蚂蚁的敏感度，更适合较小的蚂蚁步。

这样的思考矛盾也融合在长流程之中了。

说一个不恰当的比喻，从婚姻为目的来看过长时间的恋爱。

太短的恋爱就结婚固然是比较冒险的行为，但是那种爱情长跑也非常容易变成爱情杀手，因为增加了中间过程的变数概率。

而蚂蚁体和断层提取两个步骤的思考过程中那些参数，就好比你对女友有若干互相矛盾的要求，女友上得厅堂可能就下不了厨房。