petrel建模学习笔记.pdf

1,在加分层时；注意输入模板%T 2)ell To ps:veil t o ps.pmCo lUBll*1234AttributeWellSurface|MDTypeAttribute nameWellSurfaceraTypeAttribute typeKIDTAGKIDTAGDOUBLEVOID IUnitmet ermet ermet ermi 11i seco|Connec t t o wel l t rac eSettings Filters V alues Filter 对 Bulk V olume进行滤波检查负的体积值.由于网格被三角切割，所以断块边界和边界周围的网格体积是精确的46,显示误差锥形图（Erro r Co ne）在井轨迹周围可以显示误差锥形图。该选项可以在Settings窗口中的Error Cone 标签中选择，误差锥形图围绕井轨迹，离起始点越远圆柱体越厚。（允许较大的 不确定性）.47,属性建模中粗化曲线时average met ho d的选择相模型中砂泥岩相的粗化方法用Most of all空隙度、饱和度曲线、及净毛比用Arithmetic mean渗透率用：（再仔细地看，好像翻译得有点问题）Harmonic mean（对低值敏感，不能有负值，每一个小层有持续的perm）Geometric mean（没有空间相关性及对数自然分布时，对低值敏感）RMS（Root Mean Square）提供偏低的平均值Most of all：仅适用于离散曲线。Median：仅适用于离散曲线，输入的数据按大小排序，然后取中间一个数。48,粗化曲线时其他选项（1）Treat log：有 As lines 和 As points 两种用Aspoints时，如果曲线没有数值点落到网格里，那么该网格不会 有这种属性，所以它不能用与离散曲线。用As lines：如果粗化算法中有体积加权，那么必然要用到As lines,其他几乎所有的曲线粗化都用它。As lines：如图所示，L表示曲线值的影响范围（作为线，在h内，任意点的值都相同且等于中间点的值），若选用的计算方法为算数平均,1.则计算公式为a=aX产 其中，ai为该点的曲线值，h为该点值的长度（取两个点的中kot al间值）由此，我们可以看出，如上图最右边的那个点，虽然点在网格外，但其值域有一部分在网格内，所以其仍然会参与网格内的值的计算，这个就称为weighted。（2）Method：有 simple、Through c ell 和 Neighbor c ell 三种。它指的是井轨迹 穿过时给哪些网格赋值。simple是给所有的井轨迹穿过的网格赋值。Through c ell是指只给井轨迹击穿的网格赋值（同时过网格顶底面或 两个相对的侧面）Neighbor c ell,若曲线穿过位与同一个layer的多个网格，那么这些 同一个layer的网格数值将会取平均值。一般情况下用simple和用Neighbor区别不大，但是如果遇到水平井的情况，那 么就不能使用Neighbor 了，这样会使很多网格都取一 个值，最好使用simple,不容易出错。49,mak e zones时如果使用的是小层的平均厚度thic k ness,那么在Build along 里应该使用TV T,这样才能统一50,off shore指软件按照真实情况给数据赋值，比如说海平面下的海拔值为负 值（注意：数据中多数值为负的情况下软件会将值转为负，但若反过来则保 持不变）On shore不管怎么输入，数据保持自己的输入值不变。但录像中为什么选on shore呢？因为什么原因？51,在well的setting中，Thic k ness标签中将厚度输入well tops是做什么用的？52,在输入测井曲线的时候，会弹出一个对话框，里面的setting中有一个other setting,跳过未定义曲线值是什么含义？为什么在录像里都要选这一选项并将值 改为-9999?增加曲线，但是NP132还是加不进去，所以这个和曲线加不进去 无关。它的含义是曲线里面的点如果值为-9999（自己根据曲线里面的值来设）那么就为无意义的点，系统自动跳过。53,用DC解释成果做成的Fluids模型（直接用LOG（*.*）的格式输入，选择 模板为Fluids）做的是一个油气水三相的一个模型，根据原始的油层，水层，气 层，干层等建立一个模型，可以直观的看见油水、油气界面。但是关于它的应用还待摸索，它和岩性模型，相模型在井剖面上对比,修改（还不知作用是什么）A：Fluids的颜色模板和其在连井剖面中的形态（最右边一个柱子）S3CodeBaseColoPatter0zm 1 1 11Oil Zone1 _ _ i _1 1 12Water Zone 1 1 13Undifferenciat1 1 1 1Cr，tac，ou Top Tarbert吩1Tamer.：一Gas Oil Contact-*一Tarbert l-OU Water Contact-加Nessl-a2125 TSSTVD1SSTVD2100-B：录像中的岩相模型（lith）的判断和他的颜色模Facies=lf(POR=0.151Facies)flu诙岫idS=-999犯恤$)l iWf|POR=0,1,2,0|lith=lf(POR=0.15,1Jith)jRunlfor l it hStyle|Info|Statistics|Histogram Colors|Operations|/尾 N目已口|南CodeNameParenti ColorPattern0shale口.一二一二-1Fine Sandt b -屏:sand幕录像q H号家月做岩相模型进行了利用POR进行了两次运算，第一次用Lith=if(POR=0.1,2,0)判断出砂子和泥，然后修改模版使之对应。第二次又从砂子中判断出好砂子利用左上图公式。C：录象中的Fanc ies模型的模版和三种资料的连井对比O Set t ing s for Fac ies0区Style InfoHistogram ColorsApply I U OK I X Cancel号-ils一 是怎么样来进行对比的:RMN 14.15RML 15.27 A11-2口3-SS7VD1:ine1shale将lith和fanc ies中的砂子、河道修改得和Fluids操作3,若lith和fanc ies中的相为泥，则消掉Fluids中对应的zm层69.乐Jll右-.操作4,若Fluids中为水层，而lith和fanc ies中为泥，则修改lith和fanc ies与Fluids 一致（水层对应砂层）总结：对比原则是以DC做的Fluids为主，增加lith和fanc ies中的砂岩，以lith 和fanc ies中的泥岩为标准，消去Fluids中的zm层。其他的情况未说明。54,mak e surfac e后可以如果看起来比较暗，则可以将它设为高亮状态Settingsoli-Material-Highlight55,用构造面做surfac e时，可以选高密度的点，也可以不选。不知道两者区别 是什么？做surfac e时候可以不选输入类型。56,连井剖面可以两边延长，在setting-plan setting中。57在 windows窗口下的well sec tion的曲线的设置中有一个Curve filling-Delete,它的作用是去掉曲线填充颜色。58,当两条曲线都在一个道里面的时候，可以将其中一条变为虚线，以便对比。曲线一setting-style-line type59,井设计的录像：，用等厚图来mak e zones,然后再将生成的horizons转成surfac e。另外再生成 一个网格，用转化而来的surfac e来mak e horizon（注意，不是用来mak e zones）,不做zone,选其中一个层做laying（只做一个laying）弄清楚到底是为什么？为什么不直接用第一个模型来做呢？是否牵扯到horizons和zones的区别？，做连井剖面的方法，将井和面放大，polygon点在井柱子上且离面要近，这 样做切面的时候才好一点。然后再切面的设置里面延伸切面的长度。，在一个horizons上做laying后只能再这个horizons上来做属性模型，做好属 性后，再在上面做连井剖面，点开粗化的por曲线，效果如下,在3D窗口中的切片中打开井、属性后，就可以开始井设计了，直接点 well design的图标，开始并设计。点击度图标（用polygons的那个工具来设计井轨迹也一样），然后在剖面中开 始设计井轨迹。在well中找到刚设计好的井，右键，打开表格，可以得到这口井的井斜等其他 各种参数。右键setting有一个mak e logs,可以做出新井的曲线，report可以得到C Settings for r Proposed 1Settings|Statistics|Time Report|Make logs|口 l Reset output sheet:Y11 intersectionI-Iconiz points第Ginto active wellHori zon弟 C Fa.xl“f断事却，Dip andUs simple dip Use aerial dip Dip influence Fine ridnd azimuth clo and azimuth calc井与层面交点等数据。Apply OK X C=clD*BcEGVeil report2Model name:zhneji3Grid name:3D-1 hid5 4Rell name:Proposed 161Exit ZoneEnter ZoneIIZXDSecient89EndpointZoneIIzDSeeient10Zone 3334165.04152576,5栩3200.25SeEinent 1D:，右键well sec tion-Layout panel By domain可以显示设计井做出来的 曲线。（但做出来的曲线不对，不知道是不是方法不对？）这个模版里面的功能 还有其他的专门找时间摸索。、v J vaobserved d*t*I*!r r r r r出厂mS G34-15B Propostd vllsEOF Prepesai 1三民 P Yll lo(s Y P(Synt 0p FOR Syrfouzto kl-l.pmhouda I1TJxJrrrppp 氐而国囱回囱国腮.国甲因由田由由World Cursor Li&t Sources 3D Vindow 1 2D Window 1 3D Window 2 3D Yindow 3 3D Window 4 tell Secti ob 1-最后效果图如下:，还有一种设计井的方法是在面上需要用来钻穿的目的层上用mak epolygons中的做散点的工具，来点上需要的点，然后再来点击well design输入点后自动生成所需要的井资料。60,尽量利用文件的复制，粘贴功能，能够减少工作量。例如：在设置数字化层面的polygon的宽度的时候，可以设置一个，然后再用复 制功能使文件夹中的其他polygon也改变设置。61,mak e surfac e的时候也可以用井点来控制，well a4justment中输入分层数据。但这么做有什么作用？其中算法等等也都可选，不同算法之间有些什么区别？62,在 surfac e 里的 operation 中有一个 eliminate where,圈一个 polygon 后可以去 掉圈外或圈内的层面。63,如果断层polygons相交，利用断层里面的T工具也可以生成断层。64,pillar girding 里面,用餐图标可以在2D窗口中移动断层的节点。65,在well sec tion中，按照某一个层面对齐，可以直接点击分层数据，右键即 可以对齐了（比进入窗口 setting设置方便）65,在断层卡片中Displacement：Min：|5 Allow hing9 Smooth：厂，关掉第二格，Allow hing,然后设置最小垂直断距，设置平滑次数。关闭Allowhing后，断层两盘的horizons会有垂直方向上的断距产生，设置最 小最大断距后，如果超过这个值那么断距按同一个相同的断距处理。66,；x Hint:Conformable means conformable to well tops Horizons with Steep slopesmak e zone中的一个选项，仅仅在从底或从顶建中使用。勾上了以后，在遇到了陡峭的坡面时，电脑会使用一种算法通过排除断层 来恢复原始地貌。67,做好油水界面以后，打开一个切片，然后生成一个打开一个new intersec tion 窗口，再打开油水界面和horizon,打开井，在切片的setting中调整井距，显示 名称和well tops就可以做一张图。Column:MD:1-MD decreases,GeoFrame styleXY68,在加载曲线的时候,如果曲线和一般的曲线不一样，MD是右大到小的，那么就选这一项。69,在基于目标的河流相模型中，除了依次设置外，还有一项值得注意的是 1囱Settings|数 Layout|Channel|Levee 纽 Trends口|Body insertion point probability:懿 Vertical:21flil人.7.勿 I 在 Trends有一个应用砂泥百分比的选项，注意一定不要漏选(选了此项后，在运算之前先在数据分析窗口中 统计砂泥百分比。/Distribution 园 Co-kriging|y Expert|9 HinF口|Hint:If modeling 的 Secondary Variable:*Propertyuse as secondary variableS I 0 POR Ur r3Horizontal surfaceVertical function屏幕录像专:Method:Local varying mean6 Collocated co-kriging蹲黑 Constant V(086049 1-1 to 1 egjnem/Lone5-1这个是segment/zone过滤器，关掉他，则显示所有的segment/zone,一般不关掉它。127,计算器From file:在里面打勾可以生成宏或者调用宏。经常使用系统变量生成总体积或者接触面上的高度等几何模型，他们在用计 算器计算的时候经常作为约束，例如，算SW的时候.If 镶嵌语句：newpro=if（a b,0,if（a c,1,2）如果 ab,贝U newpro=0,否则，若 a c,newpro=l,其他的=2。128,在接触面上右键operation,有一个属性计算，可以生成一个离散的油水界面属 性模型。通过过滤器可以只显示油或者只显示水等等。他也能够输入到属性计算 器中来参与运算。129,关于储量计算，体积高度图:生成一个基于多次实现的分布函数.例如:一个STOIIP图,反应了 对每一个x,y的stoiip的集中分布。在一般设置里指定了生成分布函数的属性后，在蒙特卡罗模拟界面里面设置个性说明，即可得到显示.，储量计算里pivot表格和标准表格是两种不同的显示计算结果的模式。，蒙特卡罗模拟（幻灯片3,第21页）：在实现标签里面选择实现的次数以及 种子数，然后在接触统计里面选择分布类型（同一、三角或者正态分布），制定 半宽和相关系数（具体说明见在线帮助），将结果复制到exc el里面生成饼状图可以更直观的看到储量的分布。129,使用注解功能（annotation func tion）能够在任何2D、3D和地图窗口显示你自 己的注解。你可以选择对petrel里的任何项目做注解（或写上你自己想要写的内 容），它可以旋转，可以抓着它移动到任何你需要它的地方。（具体在什么地方设 置还没找到）130,关于粗化网格，当使用来自于petrel之外的网格的时候，网格也可以被粗 化。这种粗化和使用正常流程粗化一样。（粗化的时候有可能会改变断层pillar 的几何形状（在一个标签里设置）131,关于粗化属性中的体积加权Property Weightingj Weight with the following properties:Upscaling:Weight with product of:Net Gross VolumePorosity Volume Net GrossSat.Oil Volume*Net Gross*PorositySat.Gas Volume*Net Gross*PorosityJ Volume is included in the product sum by toggling Use volume weighting.132,参数、语法汇总:在不确定流程中 用这个参数）$c on取两个数之间的随机值（比如设定油水边界的时候可以在流程窗口中；圆风neric Expression给一个变量命名，前面加上$符号，后面是具体的值或者一个表达式.例如$PI=3.14$random=RAN(l 00,5000$low=Normal(30,2),具体的语法为RAN(A,B)Normal=(mean,standard deviation):在A和B之间给一个随机的数，A和B为实数值。INT(RAN(A,B)在A和B之间给一个随机的整数值。0.40.6,0.8U ses the five values defined in the listbyselec ting the first number for the first run,the sec ond number for the sec ond run,etc.Number of runs(loops)must not exc eed number of list entries.其他语法:value of$Top is the name“Ek ofisk”.1 Clear all variablesDefines a string expression.Example:$Top=Ek ofiskTheall variables in the work flow(so that the variablenames c an be re-used with another meaning further down in the work flow,or you c an refer to another work flow/Read output sheetH|that uses the same variable names but other definitions).Reads numeric and string variables from a tab separated output sheet.The user must spec ify in whic h row and c olumn the variable is found.The text file must have beenread into memory by using the Load Output Sheet c ommand.计算器语法：i几何函数:True vertic al depth(only for the well log c alc ulator)Syntax:NewLog-TVD()Creates a new log with the values as the true vertic al depth along the well path(TV D is positive below sea level).Measured depth(only for the well log c alc ulator)Syntax:NewLog-MD()Creates a new log with the values as the measured depth along the well path(MD is always positive).Depth(only for the property c alc ulator)Syntax:NewProper ty=Dep th()Eac h c ell gets the value as the depth to the midpointof the c ell.AboveSyntax:NewProper ty=Abo ve(-2000)Eac h c ell gets the value as the height above the spec ified depth(-2000).Height(only for the property c alc ulator)Syntax:NewProperty=Height()Eac h c ell gets the average vertic al height of the c ell.V olume(only for the property c alc ulator)SyntaxNewProperty=Volume()Eac h c ell gets the value as the bulk volume of the c ell.I,J,K,X,Y,ZBy using these ic ons the user c an write expressions whic h inc lude information about the I,J,K or the X,Y and Z-value.The I,J and K are the c ell numbers in the grid and the X,Y and Z are the c oordinates.NewProperty-3*I+JNewProper ty-Proper ty.SumSyntax:ii变量属性SumSyntax:Generates a property with a c onstant value.E.g.Sum(Porosity)will summarize all porosity values within the 3D grid-and eac h c ell will have that sum value as its value.AverageSyntax:Ne wProper ty-Proper ty.MeanGenerates a property with a c onstant value.E.g.Mean(Porosity)to obtain average porosity value within the 3D grid-and eac h c ell will have that mean value as its value.MaximumSyntax:NewProper ty=Proper ty.MaxGenerates a property with a c onstant value.E.g.Max(Porosity)to obtain maximum value within the 3D grid-and eac h c ell will have that max value as its value.Standard deviationSyntax:NewProper ty=Proper ty.StdGenerates a property with a c onstant value.E.g.Std(Porosity)to c alc ulate the standard deviation of the 3D grid-and eac hc ell will have that Std valueas its value.MinimumSyntaxNewProper ty-Proper ty.MinGenerates a property with a c onstant value.E.g.Min(Porosity)to obtain minimum value within the 3D grid-and eac h c ell will have that min value as its value.Number of defined valuesSyntax:NewProper ty=Proper ty.NGenerates a property with a c onstant value.E.g.N(Porosity)to c alc ulate the number of defined values within the 3D grid-and eac h c ell will have that N value as its value.iii 一般函数Random valueSyntax:NewProper ty=Ran(0,1)Eac h c ell gets a random value between the spec ified numbers(0 and 1).Triangular distributionNewProperty=Triangular(0,1,2)Generates a triangular distributed property with a minimum value of 0,a mean of 1 and a maximum of 2.The values are randomly spread throughout the model.Normal distributionNewProper ty-Normal(lGenerates a normally distributed property with the mean value equal to 0 and the standard deviation equal to 1.The values are randomly spread throughout the model.Log normal distributionNewProper ty=lognormal(0,1)Generates a log normally distributed property with the mean value equal to 0 and the log of the standard deviation equal to 1.The values are randomly spread throughout the model.Trunc ated normal distributionSyntaxNewProper ty=TruncNormal(0,1,0.1,0.5)Generates a normally distributed property trunc ated at-0.1 and 0.5,with the mean value equal to 0 and the standard deviation equal to 1.The values are randomly spread throughout the model.Trunc ated log normal distributionNewProper ty=TruncLogNormal(0,1,0.1,0.5/Generates a log normally distributed property trunc ated at-0.1 and 0.5,with the mean value equal to 0 and the log of the standard deviation equal to 1.The values are randomly spread throughout the model.Inversec umulative normal distributionSyntax:NewProper ty=In vCumNormal(0,1,0.7)Returnsa spec ific value representing the requested probability on the c umulative distribution c urve for the normal distribution.In this example it will return the value at whic h 70%(0.7)of numbers pic k ed at random from the normal distribution will be below.iv逻辑综述OrSyntaxNewProperty=If(Property3)OR(Property5),Property,U)If the value in the c ell is between 2 and 3 or between4 and 5,then it will be unc hanged,otherwise it will be set as undefined.v，数学函数Round upwardsSyntaxNewProper ty=Round(Property)E.g.Round(3.76773);4.000000Absolute valueSyntax:NewProperty-Abs(Property)Exponential,exE.g.Abs(-2)Syntax:NewProperty=Exp(Property)Power,ySyntax:NewProperty=Pow(Property,4)Gives the values for Property4Square rootSyntaxNewProperty=Sqrt(Property)Integer valueSyntax:NewProperty-Int(Property)E.g.Int 76773)=3.000000Natural logarithmNewProperty-Ln(Property)Logarithm of base 10NewProperty-Log(Property)133,在做模型的时候，比如构造模型里面的Horizon和属性里面的孔隙度等等,可以在其中加入种子数，当设置了种子数运算后。在不确定流程里面算储量的时 候如果应用这个流程，那么在这个流程后面会自动带上随机种子函数。134,关于流程窗口的应用例子，见Release Note第12页。135,首要组分分析在聚类分析流程窗口里面，点击PC分析，即为首要组分分析，在进行聚类 分析之前优先做首要组分分析。他能够和相关系数表格里的数据进行对比，使曲线按要求被使用或者被拒绝。136,在Mak e/Edit polygons流程里面，有一个强 按钮，它允许用户直接在2D 或者3D窗口中生成点状的polygons,并且将其输入到井位设计里面来设计#（也 可以修改点的位置）137,在编辑井斜数据的时候，注意，若不是从。开始造斜，那么必须在第一行 写上0,0,0o否则会发生错误。井斜变化很大。138,两个断层做削截的时候，如果断层的公共pillar相互间相交，那么做出来的 层面有可能会卷，这时候重新来调整断层，拉开角度或者断层高度，或者重设公 共pillar,使之尽可能的少相交，这样才可能最大程度的避免地层打卷。139,在调整断层的时候把断层在构造图上的polygons输入到软件中，调整断层 的时候按照polygons范围来调整，这样做horizons的时候断层两边的开口会很准 确。140,在建NP1-1区的模型的时候，划分小层依据是根据沉积相来划分，不同的 沉积相所划分的标准都不一样，例如辩状河，典型的二元结构，沙砾下粗上细，所以划分小层的时候把泥岩来作为顶，砂子做底。141,在做孔隙度模型的时候，有时候出现no upsc aled c ell is found，问题是模型 用相来约束了，而粗化的相模型和粗化的孔隙度模型对应不上，所以会导致这个 问题。去掉相约束就可以计算，但是效果很差。142,在计算孔隙度模型的时候，用相做约束的时候，如果是泥岩或者其他没有 孔隙度的相直接给定0值的，如果粗化曲线上的值和相不是很对应，那么模型在 粗化曲线的那个网格效果会不好，就好像是发生突变一样。143,右键总曲线，有一个曲线编辑；在其中的算术操作中可以对曲线的值做加减 乘除和给定一个固定值。修剪操作（c lip）中可以给定两个值,使曲线值作一个高值、低值截断。阻隔（bloc k）可以给曲线一个固 定值，使用不同的平均方法。平滑（smooth）在给定的范围内 给曲线一个平均的值，使用不同的平均方法。削峰（Despik e）还有其他许多操作，需要自己去 摸索，上面写的只是我字面上的理解，错误很多。144,右键golbal well logs,可以看到里面可以生成很多其他的曲线，包括组合 曲线、校正声波曲线、波阻抗、反射系数等，需要仔细摸索。145,在well sec tion中更改曲线道的宽度，直接在窗口浏览器中，找到该井该曲 线，在设置中有一个宽度，打勾，然后就可以自己设置宽度了，同样，在其中还 可以设置背景、横向网格等其他的一些参数。146,在生成的油水界面上右键，有一个插入/更新流体曲线，点击以后在well sec tions中可以看到自动生成了一个离散曲线，它指出了哪处在油以下，哪些在 油以上。这个功能和在well tops上右键生成的一个zone曲线一样，在well sec tion 窗口中可以看见。147,在曲线拉齐的一个界面中，在他的下面，有一个井拉齐，选一口井，应用 以后，其他的井会放大或者缩小比例，使其他的井的分层和选定的井分层一致（全 部分层在剖面上都在一条线上）当很多井都用的一个曲线的时候，把很多井用这种方式拉齐，这样就可以 很容易看到分层时所用的依据，因为，这个时候主要突出分层，很容易看到分层 附近曲线的变化。148,在作图窗口中，右键horizone,在里面有一个等高线的设置，其中，start设 置的是等高线的间隔距离（两个等高线值之间的间隔，使很多等高线不至于互相 重合看不清），还有一个interver设置的是等高线的数值显示距离。149,P32页中关于克里金的计算方法应用，仔细阅读，然后试验。150,小工具的使用，第47页，仔细看151,幻灯片阅读心得（吴键）在储层建模的流程里面，完整的模型包括了沉积相模型和岩相模型，其中沉积 相模型是建立在沉积相研究的基础上，在这以后才是属性模型的建立在软件中怎么样来进行地震数据的时间深度区域校正？地震数据在进行时深 转换的时候不可能在每个井点都吻合得很好，所以地震数据和井上的数据难免有 些地方不会吻合得很好，这个时候就应该以井上的数据为准，如果相差太大再来考虑他们之间是否存在错误的问题。在油田的以往的工作中，小层界限的划分主要是针对单砂层的界面以及相互之 间的关系，并没有划分出含有泥岩在内的小层地质界限，即等地质时间的沉积界 面。在进行建模的小层划分的时候应该以泥岩段的地质层位归属为主，这需要在以往工作的基础上重新划分出最基 本的沉积单元或者沉积旋回。所以小层应该重新对比，对比的目的是确定小层的 沉积界面，在确定小层界限的时候应该主要参考地层的沉积韵律，小层应该是可对比、可分辨的最小沉积旋回或者 地层单元。目前国内的测井解释主要只针对有效储集层，薄层砂岩、泥质砂岩、致密砂岩、泥岩等非有效储层通常不做解释，全部以0值来代替，虽然这种解释方式对单井储