petrel学习资料

答：在Petrel中可以通过和多种方法来实现过井剖面图和栅状图，下面用列举两种方法：1、用Polygon转换成intersection的方法来做联井井栅状图，具体做法是，在二维井位图上创建过井的Polygon，然后右键此Polygon=〉create vertical intersection

另外的一种方法：

步骤1：在Project Explorer 的 wells资料夹中点右键选Create Well Intersection Fence。此时wells资料夹的最下面出现Well Intersection fence

步骤2： 双击Well Intersection Fence，在设置属性页中加入井；

步骤3： 选定Well Intersection Fence，按下可视化按钮；

步骤4： 在属性模型中选择要显示的内容。

问：请问一般贵公司经常使用那些方法建立相模型

答：用得比较多的是岩相建模,从方法上来看序贯指示较多一些,但工区的情况不同,方法上也有一定的差异性.我对各种方法的区别上做了一下统计,可以作为参考.

1布尔模拟方法 

该方法主要适用于油田勘探阶段以及开发早期井间砂体和非渗透镉夹层的描述，也可以用于岩芯的描述。 

布尔模拟能够忠实某种离散参数的地质形态，如河道、沉积砂体等。该方法的主要优点是：（1）很容易用于二维和三维建模；（2）所用的参数较少；（3）非常灵活。 

它的主要缺点在于统计推导复杂且困难，模拟结果很难忠实于局部的数据，如钻井所遇到的岩性序列，这些缺点了这一方法更广泛的应用。 

2示性点模拟法 

这种方法的优点是难以条件化到已知的实验数据，优点是服从于几何形体的空间分布。 

对于辫状分流河道沉积相的模拟，可以建立相应的随机模型，并给出相应的结果。 

这种模拟方法非常优秀，能够模拟任意定义的几何形态。然而，麻烦就在于我们无法知道模拟对象（三类储层）的几何形态、长度、宽度、高度和垂向倾角，水平方向倒可以知道。 

3序贯指示模拟方法 

在序贯指示模拟方法中，不要求对条件分布的参数形式作任何假设，在现有数据和其它相关资料的基础上，根据门楷值将条件数据转化为指示数据。根据相应的指示变差函数，应用简单指示克里格给出每个网格节点处条件分布的估计值。 

序贯指示模拟最大的优点是可以模拟复杂各向异性的地质现象及连续分布的极值。对于具有不同连续性分布的类型变量（相），可指定不同的变差函数，从而可建立各向异性的模拟图像。其它优点在于：（1）数据的条件化是模拟的一个整体部分，无需作为一个单独的步骤进行处理；（2）变量的分布形式无需作任何假设；（3）协方差或变差函数不限定于某些特殊类型；（4）能综合软信息。主要缺点是：由用户定义的一些参数，参数如果选择不好，产生的结果会存在一些差异性，也就是随机性模拟的不确定性问题。

问：我看操作手册中介绍在计算储量时可以自动生成直方图的,是否在histogram  window下,在Result中的volumetrics下勾选相应的参数就可以了?我是这样试的,但什么也没显示,是否是要运算到一定的次数才可以生成直方图呢?请指教.

答：这种直方图是储量分布的概率累积分布图，需要多次储量运算结果才能得到这个图。

问：那具体需要多少次运算有要求吗?我按照储量计算不确定性工作流程,做了10个loop,还是得不到概率分布图

答：10个就已经可以绘制概率累积分布图了。

问：把断层面上的horizon lines 调整之后，重新 make horizon，为什么结果不符合我所做的调整呢？是因为断点的问题，还是我某些地方的选择不对？

答：

在Make Horizons 对话框里，需在Settings 表单下的Fault Re-sampling from the Fault Model项选中Lock all Re-sampled Horizon ，这时make horizon才能用到修改后的fault line数据。

如果fault line（horizon  line）修改幅度不是很大，可以用另一种快速方法调整断层线附近的层面：

在欲修改的层面上点击右键调用Settings菜单，选择Operation表项，在Resample to/from fault 下单选第一个按钮Fault model －> 3D grid后选择断层模型， 执行Resample，很快就可以修改此层面

问：welltops里面如何输出分层的TVD数据

Joseph：

wells里面的tvd（）＝well tops里面的Z＝SSTVD（Sub sea true vertical depth）

1、把Wellheader文件当作Production Log输入（KB为某种Production Log，其它忽略），

2、然后把KB曲线转换为Well Tops上的一个Attribute。

数据分析时带宽、搜索半径和步长应该怎样根据实际情况来定。

谢谢好心人回复，小弟第一次发帖，急切等待回复！

答：翻过了很多帖子，总结如下（所有权属于以前的作者）：数据分析时带宽——半个工区的宽度、搜索半径——半个工区的长度，和log distance ——1-1.5个井距。容差角需要考虑你的井网形状。我认为一般要大于你的井网的最小角。如果不对，欢迎指正啊。

有100口井在平面上一组裂缝线，井和井间知道动态链接的裂缝线，怎么能投影到模型里做裂缝模型，遵守裂缝线走向

答：Overview of Discrete Fracture Network (DFN) Approaches

The movement of hydrocarbons and other fluids in fractured reservoirs or in conventional reservoirs with significant fracture permeability often is not as expected or predicted. This behavior is seen in early water breakthroughs; reduced tertiary recovery efficiency due to channeling of injected gas or fluids; dynamic calculations of recoverable hydrocarbons that are much less than static mass balance ones due to reservoir compartmentalization; and dramatic production changes due to changes in reservoir pressure as fractures close down as conduits. These problems often lead to reduced ultimate recoveries or higher production costs. 

Production experience with fractured reservoirs has repeatedly shown how understanding and exploiting the fracture connectivity at the reservoir scale is an important factor for optimizing reservoir performance. Discrete fracture network (DFN) models portray fractures and their fracture connectivity very differently from other methods. For example, the conventional method to simulate fracture-dominated reservoirs is to represent the rock as a dual-porosity, dual-permeability continuum. The matrix is represented as blocks or slabs. The fractures are mathematically represented as another continuum spatially coincident with the block faces. The properties of each matrix block are represented as a symmetrical tensor, and the properties are continuous throughout the entire block. Smaller fractures within each matrix block are coupled to the rock by means of factors that are related to the fracture geometry, such as the Sigma factor. Values for the fractures are generally specified as transmissibility multipliers.

This type of numerical model cannot reproduce many commonly observed types of fractured reservoir behavior because it does not accurately reflect the geometry of fluid flow pathways. Conventional dual-porosity models do not explicitly model the geometry of discrete fractures, solution features, and bedding that control flow pathway geometry. This inaccurate model of discrete feature connectivity results in inaccurate flow predictions in areas of the reservoir where there is not good well control. 

DFN models afford an alternative approach. These models more realistically model the connectivity of the faults and joints that give rise to reservoir-scale and well-scale non-continuum flow behavior. In the DFN model, each conductive fracture is modeled explicitly as one or more 1D, 2D, or 3D element. Physical properties such as transmissivity or storage, and geometrical properties such as size, elongation and orientation are assigned to each polygon based upon measured data or geologically conditioned statistical distributions derived from measured values. 

A DFN model typically combines deterministic and stochastic discrete fractures. The deterministic fractures are those directly imaged through seismic or intersected in wells. Other, usually smaller-scale fractures may not have been detected through seismic, yet may be very important for reservoir performance. These fractures are generated stochastically. The geometrical and physical properties for these stochastic fractures are assigned through Monte Carlo sampling of relevant distributions, which may also be conditioned to both structural geology and depositional framework. Examples of fracture models conditioned to fault systems, folding, doming and three-dimensional stratigraphic reconstructions are shown. 

One of the many advantages of the DFN approach is that it makes consistent use of a wide variety of disparate geological, geophysical and production data, which conventional dual-porosity models cannot incorporate to the same extent. Data which can be used for constructing DFN models can be derived from lineament maps, outcrops, 2D and 3D seismic, well logs of various types, core, single well and multi-well production tests, flow logs, injectivity profiles, as well as structural or depositional conceptual models. Specialized tools have been developed to derive the necessary input data for DFN models from these sources, as illustrated in the presentation.

DFN models have been used for a wide variety of exploration and production purposes over the last decade. 

突然想到的，似乎可以实现，但是没有经过测试：

1、老套路，把Polygons做Surface，算法还是Artificial algorithms，参数已经给你说过。结果我们暂称SurfaceA

2、Make simple grid，做一个方盒子网格，垂向和横向要足够精细。结果我们暂称GridB

3、在GridB基础上构建Facies Model，使用Assign value的办法，输入Surface就是SurfaceA。结果我们暂称FaciesC

4、把FaciesC转换为Seismic（这需要一个地震数据作Template）。结果我们暂称SeismicD

5、把SeismicD改成Ant模板，然后使用SeismicD做Automatic fault extraction。其结果断片数据我们暂称FPatchesE

6、Create discrete fracture network，使用Deterministic算法，数据源使用FPatchesE

另一个可能性方案：

1、创建一个新模型

2、利用Polygons直接作垂直断层

3、把垂直断层转换为Fault Surfaces

4、然后从上面方案一所作步骤6，输入为这些Fault Surfaces

方案一的优点在于我们可以对于裂缝的倾角、方位角作一些变化或者朦胧化处理。缺点就是要应用到地震模块还有Ant Tracking许可的问题也比较突出。

3答：告诉你一个最笨，但也最简单、最可靠、最准确的方法：

做一个相模型，然后在模型上根据解释的裂缝一个K网格一个K网格的用刷子刷。裂缝的三维形态不管多复杂都可以准确的表达出来。

然后把刷好的相模型输出成地震数据体，以后你想怎么用都可以了。

试想一下，刷一个K网格最多需要一分钟。一个小时就可以刷60个K网格，一天坐着不动地方，完全可以把整个模型刷完。不一定比反复计算麻烦多少。

顺便透露一下，刷模型是我常用的一种方法。遇到依靠计算不好解决的问题，我一般都会直接上刷子。我将其称之为“模型解释”，归为地质综合研究的一部分。

我加载的射孔数据，最后不能输出来呢？

对射孔数据进行粗化，然后用Connection格式输出就可以

有谁知道辫状河流相沉积砂体主变程值

各位大侠？我在建的是渤海辫状河流砂体模型，(开发中后期，井距为300-400m）。由于没有经验，不知道该主变程值应该取多大合适，根据你们的经验，帮我看看下面的变差设置是否合适？有哪位曾建过辫状河流沉积的模型，辫状河流沉积的砂体通常主变程设为多少合适？是否可以给个大致区间范围？

这个变程具有不确定性，同样是辫状河流沉积，不同的研究区域，其变程也是有差异的，没有定数。所谓的变程，就是在变程范围内，井点数据具有相关性，属性预测的时候会起约束作用。建议你先做分析。

1、如果井较多，可能通过数据分析来获取变成，但最好还是结合地质认识，例如相图、小层平面图、地震相。。等

2、个人建议，密井网给小变程。结识实际情况

3、如果沉积相类型较多，形态较复杂，“两点地质统计”不好控制。建议使用“分级相控”（例如：基于目标模拟低弯曲河道，然后在背景相上用序贯高斯模拟其他相。比较灵活）或者“多点地质统计”来实现。

主变程与河流长度大致相当。次变程相当于河流宽度。不知道这种理解是否正确？ 

    1、可以参考河道几何性质给，但 河道的长宽是变程的“上限值”。还要结合井上的情况，可以把井上沉积相饼图投到相图上来观察。

   2、如果有相图，最好用相控。但变程不要给大，否则相图控制不住。

  3、如果用变差函数（两点统计）的话，那么对于低弯曲河道往往不好控制，因为两点统计只能给一个主变程方向，而河道一般不笔直的。要想做精准一些，就要把河道按不同方向分成不同region,然后不同region给不同的主变程方向。

4、应用基于目标的方法或者 多点地质统计

那位大侠能详细说说蚂蚁追踪技术？

我现在看的文献都是英文的，中文介绍太少了，那位大哥大姐能详细介绍蚂蚁追踪技术。

或给一些中文的资料， 谢谢！

选自之前在论坛上搜集的资料 发帖人：我给忘了

中油网消息：斯伦贝谢公司Ｐｅｔｒｅｌ自动构造解释模块荣获《世界石油》杂志２００５年"最佳勘探技术奖"，它利用先进的"蚂蚁追踪"算法，克服了解释主观性，有效提高了断层解释精度，大幅缩减了人工解释时间。

　　弄清断层体系断层面变化趋势及流体流动特征，是储层描述的最主要内容之一。虽然三维地震资料空间"立体"解释技术已经发展很多年了，但直到目前断层面解释仍然存在很大的主观性。 

　　伦贝谢公司最近推出的新型勘探技术---Ｐｅｔｒｅｌ自动构造解释模块，完全改变了这一状况。它利用先进的"蚂蚁追踪"算法，克服了解释工作中的主观性，有效提高了解释精度，大幅缩减了人工解释时间。

　　Ｐｅｔｒｅｌ自动构造解释模块通过产生三维地震体，清楚显示断层轮廓，并利用智能搜索功能和三维可视化技术，自动提取断层面。该技术重点研究构造地质，而不是常规的断层拾取，加深对地质及构造细节的认识，增加对储层特征的了解，最终更好地建立复杂储层模型。 

　　Ｐｅｔｒｅｌ自动构造解释模块在应用新的数据校正程序的基础上，采用先进的"蚂蚁追踪"算法，使地质专家以更宽的视野完成断层解释，增加构造解释的客观性、准确性及可重复性。

　　利用该技术的自动提取断层功能以及极坐标图和各种筛选程序，可抽提感兴趣的断层体系。"蚂蚁追踪"算法可根据工作流程需要，按任意比例自动提取断层。例如，在勘探阶段，可将工作重点集中在寻找跨盆地的大型构造断层体系以及确定它们对勘探前景的影响等方面；而在储层评价阶段以及开发和生产阶段，可采用同样的方法，将主要精力放在自动提取往往会影响油气最终采收率的那些局部的小型断层和断层体系上。 

　　"蚂蚁追踪"算法的工作流程分四步：增强边界特征，突出特殊的地层不连续性，预处理地震资料；生成蚂蚁追踪立方体，提取断层；确认、校验断层；创建最终断层解释模型。

　　流程的第一步包括利用边缘探测手段，增强地震资料中的空间不连续性，并通过噪声压制技术，随意预处理地震资料。

　　第二步建立蚂蚁追踪立方体。蚂蚁追踪算法遵循类似于蚂蚁在其巢穴和食物源之间，利用可吸引蚂蚁的信息素（一种化学物质）传达信息，以寻找最短路径的原理。在最短路径上，用更多的信息素做标记，使随后的蚂蚁更容易选择这一最短路径。 

　　该技术原理就是在地震体中设定大量这样的电子"蚂蚁"，并让每个"蚂蚁"沿着可能的断层面向前移动，同时发出"信息素"。沿断层前移的"蚂蚁"应该能够追踪断层面，若遇到预期的断层面将用"信息素"做出非常明显的标记。而对不可能是断层的那些面将不做标记或只做不太明显的标记。"蚂蚁追踪"算法建立了一种突出断层面特征的新型断层解释技术。通过该算法可自动提取断层组，或对地层不连续详细成图。 

　　流程的第三步需要人工交互操作。提取的断层必须要经过评价、校正和筛选，才能得到最后的解释结果。这一过程需要利用交互式立体网络和柱状图滤波器来完成。

　　流程的最后一步，确定的断层既可用于进一步地震解释，也可直接输入到断层模拟中。

　　利用Ｐｅｔｒｅｌ自动构造解释模块，三维地震资料解释人员不需要人工逐一勾绘断层面，从而可用更多的时间和精力研究断层面走向、分析自动提取断层之间的相互关系。 

　　该技术的优势主要包括以下几个方面：提高构造解释精度，改善地质细节描述；大大缩短枯燥的人工解释时间；提供客观、详尽、可重复的地层不连续性构造图；充分利用地质模拟技术，更好地建立复杂构造模型，优化应用三维地震资料。

从petrel英文说明里翻译的，鉴于我的英文水平，我鼓起很大的勇气才发的。打这些字的时候脸都烧烧得。呵呵

主要的目的是希望能对楼主有所帮助！

里面的很多专业术语是我自己翻译的，可能不太对。凑或看吧

在断层解释领域，这种介绍了新的范例独特的算法是一项创新的工作流程。包括四个步骤：第一步是通过减少信号中的噪音使地震数据处于良好状态，第二步是提高地震数据中的空间地震波速度突变面（断层特征产生，尖灭引起），第三步产生蚂蚁体，通过压制噪声和保持非断层同相轴来显著改善断层属性。这个被达到通过效仿生物界中蚂蚁的行动他们如何应用信息素标志路径为了最优化寻找食物。同样地，人造蚂蚁就像种子点在地震不连续体上寻找断层区域。虚的信息素被蚂蚁捕获的地震体中与断层有关的信息配置。结果是一个非常清楚详细地指示地震体的属性体。

不连续的三维断层片可以被蚂蚁体提取出来，断层可能会被断开在不同的非交叉子系统里。那些子系统减弱断层面的提取物在一个的步骤中。

蚂蚁属性的输入是一个边缘增强体，就像方差体或chaos（信号嘈杂程度）。

做蚂蚁体之前最好将地震体缩小，一般来说2G内存1.7G CPU的机子做100M的体需要20分钟。

连贯的追踪者基于“群体智能”为了寻找最理想的连通性for断层特征在一个边缘检测体中。结果是自动解释的主要成分并且允许方向滤波的独特性能。

作为蚂蚁体追踪算法中的第一步，每一个中介为在中介的地域中被识别的局部最大做一个初始的评估方向（蚂蚁边界参数）。这个方向评估将会指示行动者的追踪方向。在程序中，行动者被在最大背离常规评估方向15%，

准备了两个不同的预先确定的设置：消极的和积极的。

消极的蚂蚁设置群体行动是保守的，用这种方法提取大的整个地区的断层。

初始的蚂蚁边界

初始的蚂蚁边界参数定义了初始的行动者分类，每个蚂蚁的领土半径，被定义的体素数量，看下图。这个参数对算法的完成时间影响很大，由于它定义了在体中传播的行动者的总数量。

出于提取大的地区性断层的目的，这个分配可以粗糙一些（5-7个体素）。如果是出于精细工作的目的，并且想画出小断层和裂缝的图，那么就要设置一个更紧密的分配（3-4个体素）。一般来说，没有必要设置小于3个体素，由于行动者会跟踪相同的同相轴并且不会加新的信息到结果中。这个参数并不意味着行动者在体中是均匀分配的。这个参数意味着没有行动者在这个半径内取代另一个，由于每一个搜索for局部最大在数据中。

蚂蚁边界被一个半径定义，在体素中 蚂蚁行动者的分配。如果行动者不能鉴别局部最大或在这个半径内做一个方向评估，那么这个行动者将会被消除。

蚂蚁追踪背离（体素数量）

蚂蚁追踪背离控制追踪时的局部最大中最大的可以允许的背离。行动者假定一个平面形态，并且与初始方向只能背离15度。这个方法允许行动者接受局部最大的一个体素在任一个预定位置合法。如果最大值在蚂蚁追踪行列之外，那么追踪背离参数开始起作用。

如果一个区域最大背离在一个平面有大量，行动者不可能长久地追踪他们。这个参数将会允许行动者从区域最大中背离并且他们可能更长久地这个构造。

数值1将会允许行动者在任一方向从行动者位置到搜索到一个区域最大背离一个体素。如果最大值没有被找到，那么这将会被记做一个违法的步骤。如果找到一个最大值，那么据这个最大值最近的合法的位置将会被用作一个合法的步骤。看下图。

蚂蚁追踪背离参数控制可能的体素 行动者可以搜索到合法的位置 到发现局域最大值。如果一个最大值被发现，最近的合法位置被应用。

蚂蚁步长（体素的数量）

这个参数定义一个蚂蚁行动者在一个搜索步长中每一次增加提高的体素数量。增加这个值将会允许一个蚂蚁行动者搜索更远，但这将降低结果的分辨力。

蚂蚁步长参数决定一个蚂蚁行动者在一个的步长中可以在多远的距离搜索一个局域最大值。

违法的步长规定（体素的数量）

违法的步长规定定义一个行动者的最终可以持续多远当没有发现合法值的时候。

例如，如果行动者通常在一个有效的位置并且向没有探测到边缘的地方移动，这将构成一个违法步骤。如果这个参数被设为1，那么行动者将中止在这个方向的搜索。如果一个有效的体素被发现，那么行动者的被记录的路径将会包括这个有效的位置和违法的位置。违法的步骤仅在合法步骤已经被记录下来后被计算在内。

蚂蚁追踪违法步长参数决定了一个行动者在没有发现一个局域最大值的情况下可以持续多少步长.

要求的合法步长

要求的合法步长，控制如何连接一个检得边缘帮助辨别一个边缘从未导向噪声中。这个参数与违法步长规定结合应用。

蚂蚁追踪合法步长参数决定包括行动者持续的边缘值得步长的数量

停止标准

停止标准涉及到一个行动者生命中违法步长对被允许的全部的百分比。这个参数将会取得当行动者向前移动，发现一个违法步长，程序再次发现一个合法步长并且继续以这种方式，当逐渐地积聚违法步长知道他们表现出一个有意义的一分在行动者搜索区域。当这个值太大的时候，我们失去信心在这个我们追踪的即将合法的断层排列中，所以停止标准百分数允许我们停止这个追踪。

蚂蚁追踪停止标准参数可以结束一个行动者当在所有的路径中的违法步长的数量超过设定的百分数的时候。

蚂蚁建立了一个小的面 这个面在每一个边缘延伸。如果这个在个别边缘面上的体素的百分数不是一个局域最大值并且超过停止标准，在这个方向的延伸的部分将会被终止。

立体图标签

另一个参数对于应用者来说是可视的，一个方向滤波器。在这个图形设计下，使用者可以限定行动者的方位角和倾斜。这个有力的参数控制允许一个发生方向滤波体，这个滤波体可以将解释中的断层系统分开。

用这个模块提取断层区域得主要方向和子方向在数据中。

方向滤波器允许使用者搜索的方向。这个图已经被象限分割，这些象限可以被个别丢弃（未被行动者搜索的）。

这个有力的参数可以被用于滤掉与噪声相关的东西，像采集覆盖区，或在精细构造研究中隔离特种断层系统。

一些模糊区域将会在计算过程中被丢弃。这意味着没有断层信息被获取在这个方向和区域内，不管是否在这个区块内存在不存在。

How to use the stereonet

1.When generating an ant-track attribute, click on the Stereonet tab. 

点到立方图标签 

2.Set the Dip and Azimuth sectors, by either using the Cycle through the dip and azimuthicon or defining user specified values using the Azimuth increment and Dip increment . 

设置Dip和Azimuth，应用Cycle through the dip and azimuthicon按钮，或设置使用者设置值应用Azimuth increment和Dip increment 

3.Move the mouse pointer outside of the Stereonet. You can use the mouse pointer outside the Stereonet circle to drag and mark areas that cover dip ranges from 0 to 90 degrees. 

在立方图外部移动鼠标，你可以用鼠标点击立方图外部拖拽并且标记区域 

4.Move the mouse inside the Stereonet. You can use the mouse pointer inside the stereonet to point or drag to mark single rejection areas. 

在立方图内部移动鼠标，你可以在里面拖拽鼠标标记的拒绝区域 

5.When you click on a shaded (rejected) sector, the sector will turn white (accepted). 

当你点击阴暗区域，这个区域将会变白 

6.If you like to invert your selection, press the invert selection icon . 

如果你要点到选择，可以点击invert selection按钮 

7.Use the Reset Orientation to default icon to clear your selections. 

应用重置按钮清楚你的选择 

8.Continuous ReadoutA red cross with continuous readout values indicates the current dip and azimuth value of the position where the mouse pointer is located 

连续的读出器 一个红色十字在持续读出数值指示鼠标所在位置当前方向的角度和增量值