【阿什卡精描专栏】一位精描“老专家”的系列文摘

ESSCA

sherry0923 发表于 2013-8-15 16:46
真的是受益匪浅。我原来是学油藏工程的，但是三年下来总是觉得自己不是很理直气壮，后来回过头来才发现是自 ...

ESSCA

问：02和04对比，在模型粗化方面的精度有哪些提高？

答：04新增的结构化模拟网格设计，可对主要目标区进行网格的细化，尽可能保留有效储层的地质特征；对于在走向和或倾向变化比较大的断层，04可对其进行断层锯齿化（包括走向或倾向），从而保留构造的精度，并提供网格的正交性。

问：04能够支持井事件的输出，并且能够实现不同时间的井措施，有什么实际意义？

答：如果没有井事件文件的话,那么Petrel只是输出地质模型，井与模型之间的连接以及射孔压裂酸化等措施在2004以前的版本里是没有这个功能。

问：利用04的新功能可以一步到位实现局部网格加密，节省数值模拟对模型的修改功能，把更多的时间和精力用在数值模拟研究预测中去，但是一步到位实现局部网格局部加密遵循什么原则和条件，合理性如何？

答：一般是进行增产措施改造的井如：压裂，酸化，*近断层的井一般需要进行局部加密,Petrel做的局部加密完全适用于数值模拟，其合理性与否主要取决于构造网格的质量.

问：使用PETREL进行地质建模的人员只需要对*近断层的井或构造形态出现突变的部位进行局部网格加密？至于相控方面需要做哪些工作？增产措施改造的井进行局部网格加密的工作是建模人员做的，应该是数模人员来做的吧？

答:Petrel建立的模型本身已经是非常精细的,在建立模型的时候已经考虑了相控因素. 局部网格加密是相对于粗化后的模型而言的,局部网格加密的目的是为了数值模拟工作能够更进一步的细化研究工作。如果模型不经过粗化而进行数值模拟在理论上是可行的，但是目前的硬件技术以及数学研究的计算方法而言还不能够完成的，这也就是为什么我们做的地质模型进行粗化后才能进行数值模拟计算，否则一个很小的区块就会因为计算迭代不收敛而失败的。而为了细致的研究某些方面，数值模拟人员就要对某些井附近进行局部网格加密工作来进行进一步的研究。局部网格加密原本是数值模拟人员的工作，现在Petrel软件能够实现此功能，能够从地质模型的工作中进行此工作。局部网格加密是指对粗化后的处理工作而言。

问：是不是认为在04中做流线模拟比在Eclipse中更快捷，修改更方便？

答：Petrel中如果要实现流线模拟功能的话则必须要用高级模块的。Petrel中添加流线模拟是为了验证建立的地质模型的精确度，从而让地质建模人员从一开始就对自己构建的模型有一个更深入的认识，对某些不确定因素引起的模型误差及时的进行修改。

问：Eclipse的模型计算结果能够十分方便的导入到Petrel里，动态部分能显示吗？

答：可以实现模型的动静态资料，只要是Eclipse中能够显示的Petrel都可以显示，04里Petrel和eclipse已经实现了无缝兼容。02能够支持输出ECL格式模型，但是部分文件（如well connection文件）需要进行重新设置, 04完全支持无缝连接ECL模型，不需要进行任何重新设置的工作。

问：如何在岩相模型中输出砂体等厚图？？？

答：在2004版是可以输出NET MAP的，只是输出位置发生了变化，在2002SE版中，你输 出哪一个层的，就显示那一层，然后如你提到的在面板中选择输出就可以了，在2004版可以不显示，而在输出中来选择即可。操作步骤是点击你要输出的属性体，选择SETTING，然后找到OPERATION那一项，里面有make map from property,但是需要注意的是不同的属性体，提示的做MAP的内容有一些差别，在相的类型里面是不让你来进行NET MAP的，我试验了一下，可以通过以下方式达到你所要的结果，首先用运算器算出一个属性体，这个属性体的类型选择为N/G的类型，那么N/G＝IF(FACES=1,1,0)，由于选择了净毛比的类型，就可以来进行砂岩厚度图的制作了，如果在你的FACE属性体上是不可以计算的，每个ZONE的选择，在下面的ZONE里面选择，也可以所有的ZONE 一起做。

问：属性建模中u与0的区别？

答：undefine是没有定义数值的值，0 是个有效数据，在运算中，0 将参与计算，而undefine 不会参与计算。

问：众所周知，目前勘探开发的大多数低渗透油气藏都具有相当的裂缝发育，另外，灰岩和火山岩油气藏的裂缝发育也是该类油藏描述的重要内容。那么，怎样来更好的预测裂缝的分布呢，裂缝的分布是否具有一定的方向性？或者根本就是随机分布的？是沿主应力方向还是别的什么方向？

答：对于裂缝建模，问题的关键是裂缝分布的不规则性，而传统的建模考虑更多的是基于高分辨率层序地层学下的等时沉积理论，也就是等时沉积的地层在平面上除了剥蚀以外具有比较好的连续性，因此，在进行裂缝预测的时候，不能过多的考虑地层的连续性，而应该把考虑的重点放在裂缝的方向和分布范围等问题上，首先在PETREL上，可以充分将井上的电测信息及生产信息充分考虑进去（在做这部分时，主要通过建立构造模型和曲线离散化来进行，但考虑到裂缝分布的不规则性，建议建立模型时，尽量用大的分层方案，分层不要过于细化，如果电测结果有裂缝分级最好，如果没有最好将井上的主力产油层认识清楚），在进行模拟的时候，需要处理好两个方面的内容，一是裂缝的方向性，一般来说，裂缝与该区的构造发育有着一定的联系，大的裂缝主要沿主应力方向，在进行模拟时，各方向的变程和方位角等是可以控制的，其二是必须充分运用地震上的信息，可以做多种地震属性体，通过优选和交互统计，找出与裂缝相关的属性体，用作协同模拟，这样预测出的效果会更好一些。

问：在petrel中，如何将三维图形的某一个方向，比如Z方向的比例放大，而水平方向的比例不变？

答：垂向比例的设定就在显示窗口的上沿，通常显示为5的。

Set Z-scale - Changing this parameter givesa higher/lower relief on the data set. It is unit sensitive, i.e. if XY are inmeters and Z is in feet, a distance of 100 units will be approximately a thirdin Z to what it is in X or Y.

问：如果我手头没有理想的断层数据，跳过fault modeling直接做pillar gridding，再通过general intersection来解释断层。那是否就意味着在完成断层建模以后还要重新进行pillar gridding ,再生成一个3D model，再在这个model中做后面的工作呢？因为我觉得第一个3D model是在没有建立断层模型之前产生的，断层模型应该是后面所有工作的基础。斑竹认为是这样吗？

答：你的想法是完全可以的，我曾经用过，pillar gridding不一定非得有断层，你可以把这个过程当作是一个中间过程，用完后将其删除就行。其实软件是死的，人是活的，我们不一定按部就班，无论通过什么方式，只要达到我们的目的就行！不过还有更简便的办法，就是直接make surface，生成面以后就可以直接利用intersection来进行解释断层了，可以省很多时间！

问：测井曲线数据文件一般都是多个层的，建模层位的深度范围往往小于测井范围，这样直接导入测井数据文件时，曲线会超出井轨迹，如何使得超出轨迹的曲线直接"去掉"？

答：一个办法可以试一试，Operationstab of Setting for “Wells” panel，set Top limit & Base limit。注意：这个过程是不可逆的！

问：在surface-surfaceoperation当中求A intersectB的操作过程中,箭头只能找到Surface B.请问Surface A是如何指定的？

答：其实在这个操作中你已经事先假定A的存在了，具体实现：分别找到你需要执行intersect的两个surface（A，B），选择其中一个（A）=〉右键=〉setting =〉surface-surface operation，然后用箭头输入另外那个surface（B），执行，那么原来的A就是你所需要的suface！

问：一个SEGY数据太大,我只用其中的一部分,应该怎么做?

答：先加载，然后作Crop。事实上加载SEGY并不会真地把数据放到工区内，Realize以后才是进入工区管理的数据。

问：该软件经常有一些窗口上的英文显示不全（同介绍和帮助上的相同窗口对比），是什么原因呀？

答：这是因为你装的系统是中文版的建议装英文版，就可以避免这个问题。

问：Petrel能否实现用对数形式表示不同曲线？

答：可以。例如我要用对数形式来表示孔隙度，先建立一个新曲线模板，设置settings，如图红色区域的设置参数，因为渗透率是用对数来表示，因此只要选择这个模板就能把孔隙度用对数来表示。要实现其他曲线的对数显示如法炮制就可以。

问：抽稀斜井建立模型，如何提取斜井的属性检验？那个软件可以？不需要自己写程序？

答：Petrel就可以阿，建立模型以后，再把抽稀的井加入到wells里面，对其make log就能沿井轨迹提取模型的各种属性了，这个道理跟井轨迹设计有点异曲同工！

问：如果在makehorizon之后的层面与输入的地震解释层面在深度上有细微的出入,是什么原因造成的?该如何校正呢?

参考解决方案（1）Petrel软件在Make Horizon时需要对原层面进行二次处理，处理的算法可能会产生一些系统误差！（2）Petrel软件在Make Horizon时在断层附近由于截断距离设置问题，也可能产生一些系统误差！

参考：（1）如果是第一种误差，而且误差超过允许范围，可以考虑输出horizon为surface，计算模型surface与实际surface的误差surface，在模型的horizon上加上或者减去这个误差surface。（2）如果是第二种误差，可以从Horizon生成断层线，通过编辑断层线来加以解决（该方法仅适用于2004以后版本）。

注：以上答案摘自顽石斑竹

问：在make zones时我需要用到Isochore面，但不知Isochore面如何得来呢？是不是必须自己把厚度数据输入呢？还是Petrel可以通过Well Top中的小层数据直接得到Isochore？

答：在输入well top的点上，先激活上一层的点，在下层的点上点右键 convert toisochore poins 然后，在 makesurface 形成面。

问：模型建好后,如何做连井栅状图???

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

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

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

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

问：如大家所知，渗透率是具有对数属性的参数，在做等值线图时应该用LOG间隔标线（如下图）， 中如何实现?

答：1. 将要作图的surface拷贝n份（按需要），放于同一文件夹中；

2. 每个surface用不同的方式作图，如第1个从10开始，间隔10，到90结束；第2个从100开始，间隔100，900结束，第3个从1000开始，间隔1000，9000结束；………………

3. 将它们同时显示在一张图上即可。

kofamazon

学习了。。。谢谢分享

ESSCA

kofamazon 发表于 2013-8-29 20:36
学习了。。。谢谢分享

不客气。欢迎和大家多多分享您的经验和知识。周末愉快！

ESSCA

问：Eclips和Petrel有什么区别和联系？

Petrel能否完全取代Eclips？Petrel能进行剩余油分布预测、水淹层解释、流体分布模型和提出开发调整建议?Petrel主要是针对储层，能否进行三维油藏数值模拟？

答： Eclips是数值模拟软件，Petrel是建模和油藏描述为主的一个软件，PETREL软件可以为Eclips建立精确的模型，Eclips的数模结果也可以很方便的加载到PETREL中进行研究。

Petrel在2004版中，将数值模拟的很多工作都进行了完善，PETREL可以进行流体模型的建立、相渗曲线、趋势控制、事件设定等数值模拟的设置，也可以进行生产及完井数据的加载及历史拟合，也可以进行完井设计，PETREL中可以完成流线模拟，因此可以说数值模拟的很多工作都可以在PETREL中完成。

如果你的机子已经安装了ECLIPSE介质，那么可以在PETREL中参数设置完后（包括构造模型）可以直接调用ECLIPSE运算器进行数值模拟。

问：现在，在孔、渗场的基础上通过J函数建立了一个饱和度场，然后利用测井解释的井点饱和度值来对这个饱和度场进行校正，应该怎么操作？

解释：你可以将由J函数建立的饱和度场作为第二变量，约束饱和度的模拟。J函数建立的饱和度与井点处的饱和度相关性如何？这个方案可以通过相关性控制来调节第二变量的可信度，最终获得井点数据校正模型原始数据的效果。

她给我们展示了一个反演的Gamma数据通过Gamma测井曲线粗化数据校正的结果与原始反演数据的比较，我们认为可以实现非常好的约束结果。所以使用第二变量适度调整相关系数值可以取得需要的最终效果。不过对于数据没有很好相关性的时候需要慎重。

问：构造如果上下层互相交*有什么好办法把它们拽开？

答：（jiayou版主）其实层的上下打架有一个很好的解决方案,打架的原因是由于井控制不住或者有的分层不正确,我建议大家这样来做,就会比较好的克服以上问题.

一般大家得到的构造是地震解释好的大层构造图,如果在MAKEHORIZON时控制层较少就会在上下层之间出现很多问题,那么大家可以考虑增加控制层位,即每隔几个小层在增加一个HORIZON来作为控制层.

这样的控制层怎么来做呢?可以采取以下方法,即趋势面约束做SURFACE的方法.

第二个值得大家注意的地方是:在MAKE HORIZON时,WELL井的分层调整作用,如果你用了上面的构造趋势面作为控制层,那么你如果再进行井的调整,它就会发生改变而起不到应有的作用了,所以加上分层数据以后,应该把调整的选项改变.

（火山版主）

避免打架的几种方法：

1、趋势面约束法：加友先生提到了，不赘述。

2、等厚图法：就是走Make Zone流程。

3、断层线约束法：利用拾取的断层与层面交线（Polygon）重新约束Make Surface。

4、人机交互修改法：（1）利用Petrel 2004版新功能－－断层线编辑和约束建模；（2）人机交互修改Surface；（3）人机交互修改构造模型网格。

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

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

1布尔模拟方法

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

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

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

2示性点模拟法

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

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

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

3序贯指示模拟方法

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

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

问：都说Petrel显示效果好，出图效果如何？

用了快两年了，一问才知道大家基本没出过大图（A0），为了让更多人了解，看到Petrel成果，介绍一下你的经验，谁出过漂亮的图让我们看看，不要屏幕截图。另外，3D图如何能打印大图而不发虚？

答：目前都实行无纸办公，我们在汇报时一个ppt 就包含了我们所有的汇报内容，报告附图也都比较小，通常是A4 的，很少出A0 以上的大图，所以大家好像都淡忘了出大图的事。

利用Petrel的地质绘图功能，我们可以高质量的打出任意比例尺的大图，如构造图、等值图、地质剖面图、油藏剖面图、地震剖面图等。

我曾在大庆物探公司出过1：10000的构造图，效果非常好。另在2004年中石油地质建模软件选型时，我们也利用petrel现场打出了各种地质图件（也只有petrel 这样做的）。

至于想将3D 模型打印出来，我建议将机器的分辨率调到最高，并将3D 窗口放大，再抓图，效果可能会好些。

问：用petrel建模后怎么网格都是不规则的形状，长方形的四个边都是扭曲的，请问是什么原因啊，有什么方法可以解决这个问题？

zhou：网格化的质量保证措施：

1、断层Pillar 间距要等距或均匀变化，Pillar 长度要相等或均匀变化；

2、所有的断层的Pillar 的高度要等高或均匀变化（尤其是对于斜度比较打的断层）；

3、在运行Pillar Gridding 前，要合理设置断层的网格方向，添加趋势线等；

4、针对出现奇异网格的地方，添加趋势线；

5、I或J方向的断层或趋势要平行。

问：petrel中处理zone中的单砂体（有有效厚度）时，直接将单砂体的顶底数据和做zone一样输入，还是在做好horizon和zone后单独解释出来做？？

zhou：如果你的单砂体的顶底界面可以区域对比出来，你可以按照zone 的方式建立单砂体的模拟单元。这是最理想的地层框架模型。

海纳百川1391

ESSCA 发表于 2013-8-15 14:04
据说不定期连载更那啥

好吧！我懂了！

ESSCA

海纳百川1391 发表于 2013-9-12 22:15
好吧！我懂了！

一笑nkk

很经典，仔细拜读中

ESSCA

大家节后工作顺利！