目录一、认识 Petrel 界面1、快捷键设置2、菜单栏3、系统设置1基础效果2地震数据设置4、项目设置1坐标与单位2自动保存3平滑处理二、导入数据三、二维/三维窗口常用功能四、构造建模-角点网格1、创建边界2、构造面 - 层位1基础信息2平滑处理3矫正先粗化后细化4检查重叠3、构造面 - 断层去除异常点4、地质网格1网格边界2平面网格划分3层位网格划分5、海拔属性建模五、构造建模-结构框架模型1、设置断层2、边界设置3、模型构建1导入层位2模型细化4、结构网格六、属性建模1、单井单属性粗化2、沉积相定量筛选3、数据分析4、插值建模一、认识 Petrel 界面1、快捷键设置左上角可自定义分别表示保存、撤回、下一步、窗口、调整颜色表、检查器、播放器、质量报告等2、菜单栏Home 下可以选择视角不同的视角对应不同的菜单栏功能3、系统设置涉及到一些个人使用的习惯1基础效果常用鼠标直接重命名节点和禁止旋转动画惯性跟随效果各功能如下所示机翻2地震数据设置后面几个不常用4、项目设置1坐标与单位选择坐标参考系统 CRS搜索关键字选择这一步在新创建工区的时候也会弹出来单位设置国内一般就是米制右边的customize是自定义2自动保存在其他设置中可以设置自动保存时间也可以选择自动静默保存不推荐3平滑处理Global surface filter全局曲面滤波器控制层面的平滑和形态约束能够对生成的曲面进行平滑处理Minimum curvature algorithm最小曲率算法曲面拟合的数学方法原理生成一个曲面使其曲率变化最小最大误差容限允许曲面偏离原始数据点的最大距离值越小曲面越贴近原始数据但可能保留噪点最大迭代次数迭代越多曲面越精细二、导入数据详见【Petrel】基础教程1-数据导入三、二维/三维窗口常用功能以三维窗口 3D windows 为例如图标注出的功能从左到右1、veiw all一键展示出包含所有数据信息的视角2、view all右边的是view from...表示一键展示从某个方向观察到的视角3、toggle background切换背景颜色只有黑/白4、show/hide display elements显示/隐藏元素常用 Axis 轴展示数据在坐标轴中的位置另外当我们处在某个窗口下点击右上角红色的 Window 时菜单栏会展示出该窗口更详细的功能以3D windows为例详细菜单栏信息如下所示与窗口中的快捷功能键对应四、构造建模-角点网格菜单栏对应为Structual Modeling1、创建边界1选中层位数据在二维窗口中展示2菜单栏中找到Polygon editing多边形编辑3点击Add points to polygon通过添加点来创建多边形最后一个点双击封闭4点击Polygon editing可以进行调整2、构造面 - 层位在菜单栏中找到Make surface另一种打开方法先在Home中找到Panes → Processes点击后可以打开左下角的Peocesses窗口包含了建模各流程对应的功能Make surface在Utilities下1基础信息包含层位数据、边界、几何设置等网格尺寸和位置、大小等根据实际情况自己设置如果要批量创建则选中设置好后点击Apply创建2平滑处理选中一个层位点击设置Setting找到操作中的平滑操作Smooth运行Run若生成的面与层位相差太远说明井震数据不匹配需要进行矫正3矫正先粗化后细化① 可以通过Calculator进行平移② 或者在Setting中找到算术运算Arithmetic operation效果与Calculator相同③ 上述调整都是粗调接下来需要通过Make surface细调矫正把该层位重新作为输入设置好边界扩展在Well adjustment中输入井数据对应的层位后 Apply4检查重叠① 插入交平面Insert intersection plane② 通过Snap 2 points用两个点选定交平面方向和位置Manipulate plane可以移动交叉面③ 选中交平面创建交集窗口Create Intersection Window④ 选中层位就可以从交平面截面观察层位之间是否重叠还可以进一步添加播放器观察⑤ 在快捷键中找到Players添加一个交叉面播放器就可以更严谨地观察 层位之间是否有动态重叠3、构造面 - 断层逻辑与层位一致只介绍不同的地方本文介绍的断层数据是以点的形式展示的去除异常点① 菜单栏中找到Surface editing在功能中选中Point editing编辑点选择移除点功能就可以去除异常点了4、地质网格新建一个角点模型corner point grids在左下角模型窗口中显示1网格边界将之前创建的层位边界移到该模型中2平面网格划分用Pillar gridding进行网格划分可以手动设置增量默认划分方向是坐标轴方向如果要调整切割方向在Edit fault model中设置通过画出两点添加新的I方向再次在Pillar gridding中 Apply就可以重新划分出沿该方向的网格在二维窗口中完成网格划分后会在模型中生成一个三维网格3层位网格划分在菜单栏中找到Horizons用第一个快捷键Append添加层位数对应的列数再用multiple drop加入对应层位效果如下在模型中双击某个单个层段可以看到对应的统计信息面板包括坐标参数、网格参数等信息其中Average Zinc是埋深统计参数表示该层段的平均埋深得到数据后接下来可以进行层段垂向细分也就是把大层段在垂向上细分成若干给更薄的模拟小层可以根据层段的重要程度考虑好垂向分辨率也就是目标层的单层厚度厚度越小分层越密然后用平均埋深厚度➗目标层的单层厚度得到分层总数在Layering的Number of layers中设置Apply后在各层生成网格格架5、海拔属性建模在菜单栏Property Modeling 属性建模中找到Geometrical 几何建模在方法中选择Absolute or relative depth 绝对或相对深度属性中选择Real 实数和Positive 正数工具栏中点击Clear display清除当前显示的内容左下角打开刚才创建的海拔属性显示双击Elevation general在Colors中设置颜色显示打开自动图例Auto legend右键模型可以打开网格显示五、构造建模-结构框架模型关闭前视图新建一个结构框架模型创建结果依然可以在左下角看到1、设置断层自动导入断层注意顺序断层中的属性设置创建好后可以在Fault relationships中看到断层间的关系2、边界设置导入之前创建的边界Polygons3、模型构建1导入层位与上文同理快速导入层位2模型细化4、结构网格创建结构网格Structural gridding导入刚才创建的结构模型Structural model设置好各层的细分数量Number of layers参考创建角点网格时的划分依据六、属性建模1、单井单属性粗化用Well log upscaling单属性粗化的本质是测井曲线的垂直尺度升级测井曲线是高分辨率采样而我们使用的网格厚度远大于采样间距因此需要把各井各属性按网格的垂向分层边界把区域内的测井采样点合并成一个值赋值给对应网格算术平均方法在单个网格对应的井段范围内将所有测井采样点的属性值相加再除以采样点总数得到的平均值就是该网格的粗化属性值插值实现属性建模的核心数学方法储层参数实测数据只来自有限的钻井井与井之间的地下参数是未知的。插值就是通过空间统计算法根据井点的离散数据结合地质规律比如沉积相约束、变差函数推算出整个三维空间里每一个网格的属性值把离散的点数据变成连续的三维属性体。2、沉积相定量筛选通过Well Calculator筛选以自然伽马 GR作为筛选属性GR 是最核心、最通用的指标数值高低直接反映岩石的泥质含量是区分砂体与泥岩、判断水动力强弱的黄金标准几乎所有沉积相划分都会用到生成沉积相曲线数据可以自定义编码对应的岩相信息创建沉积相属性3、数据分析用Data analysisZones按层分析Facies按相沉积相沉积环境形成的、地层性质接近的储层类型分析以计算生成的 facies 为例可以分析岩相比例和随垂向深度的变化趋势上述窗口是岩相比例部分除此以外还有厚度、概率等内容可以分析Thickness厚度统计统计每种相的单层厚度分布规律最终用来给三维相建模提供砂体厚度的统计约束上图标准差远大于平均值说明厚度数据波动极大离散性很强绝大多数相段集中在左侧说明厚度越小的相段出现频率高正偏态分布Probability概率统计相的空间接触与出现概率分析次级变脸和砂岩出现概率的对应关系上图展示的是不同深度下砂岩出现的概率Variograms变差函数是随机建模最核心的参数描述相在空间上的相关性量化储层属性在空间上的相似程度随距离的变化规律核心作用是定义岩体的几何形态和空间连续性核心参数①Nugget块金值距离趋近于 0 时的微小差异理想情况下接近 0②Sill基台值曲线走平时的高度即最大差异值理想情况下接近 1③Range变程曲线到达基台时对应的距离变差函数方向①Major direction主方向岩体延伸最长的方向变程最大②Minor direction次方向垂直于主方向、岩体延伸最短的平面方向变程比主方向小③Vertical direction垂向方向地层上下方向对应单岩体厚度变程通常最小常用理论模型拟合散点的数学公式①Spherical球状模型上升到变程后立刻平稳对应边界清晰的岩体②Exponential指数模型缓慢上升、渐进平稳永远不会完全走平对应渐变型、边界模糊的储层③Gaussian高斯模型上升最平缓对应极度均值、渐变很强的储层Experimental variogram computation 实验变差函数计算表格参数含义列名含义Direction方向名称垂向 / 主方向 / 次方向Azimuth平面方位角正北为 0°顺时针旋转Dip方向倾角垂向 90°平面 0°Number of lags分多少个距离组Lag distance每个组的距离宽度Search radius最大搜索距离超过就不再统计Band width方向带的宽度增加样本量Tolerance angle角度宽松范围Lag tolerance距离的宽松范围Thickness垂向厚度容差Declustering去聚类修正井分布不均匀导致的统计偏差4、插值建模用 Facies 进行相建模