Petrel地震地质解释和建模使用技巧.pdf

Petrel地震地质解释和建模使用技巧 Petrel 合成记录工作流 制作合成地震记录 进行层位标定和确定时深关系是地震解释工作中非常重要的环节 从Petel2009 1 1 开始Petrel里有两个制 作合成记录的模块 一个叫Synthetics 一个叫Seismic Well tie 这里介绍如何使用Synthetics模块制作合成地震记录 从Petrel 2007开始Synthetics模块有了很大改进 最重要的变化是其结果可在Global well logs下有相应的synthetic目录 其相 应时深关系可在数据表中显示 对同一口井可产生多个合成记录 如图1 1 1 2所示 Synthetics模块制作合成记录工作流主要分为两大步骤 按照已有数据产生合成记录 通过welltop 进行时深关系调整 bulkshift或sqeeze stretch 一 生成合成记录 1 双击synthetic模块 打开合成记录主界面 如下图 选择create new folder 从界面中well 到well seismic 四个界面对合 成记录中所需数据进行选择或创建 如图2所示 Well 选择要做合成记录的井 可多选 但每口井必须有相应的数据 DT和子波 Sonic and time 确定原始输入数据及时深关系 根据实际数据品质 如果有checkshot 可用来做DT曲线校正 所有井上时深关系以工区井目录 以及每口井的Settings界面里 Time界面下设置为准 Synthetics界面里的Overwrite global time log项不启用 Create synthetic seismogram 创建合成记录 选择创建合成记录所需数据 Density Acoustic Impedence Reflectiotion coefficients和Wavelet 如果这些数据都不存 在 或者希望修改参数重新创建 则点击黄色星状按钮创建新数据 如果要使用已有数据 则点击三角下拉菜单选择即可 如图 3所示 不同的数据组合产生不同的合成记录 同一合成记录文件夹中可保存多个合成记录 Create well Seismic 产生井旁道地震数据 进行合成记录的对比 在所有数据选择好后 点击apply 进行合成记录计算 其结果有相应的目录保存 在Input数据的global well logs下产生Synthetic1文件夹 其中保存了所有数据 AI synthetic seismic 鼠标右键点击Synthetic1文件夹下 选择菜单上的spreadsheet选项 打开数据表格 查看该次合成地震记录分析得到的时深 表 且该时深关系自动作为井当前时深关系 其Time log 也会按该时深表调整 如图4 1 4 2所示 2 打开well section 窗口 进行合成记录显示 同时显示井和Synthetic1目录数据 如图5所示 如果合成记录与地震数据有偏差 可以进行适当调节 二 时深关系的调节 1 打开合成记录主界面 在Sonic and Time下 打开overwrite global time logs 选择DT曲线作为井上时深关系 不启用 manual adjustment 如图6所示 2 从Spreadsheet 中清空well top 表中TWT picked列 见图7 3 激活well top edit edit well top time按钮 如图8 通过well top进行深度调节 此时well top 显示出调节前后虚实两个位 置 见图9 Apply后所有的调节最终生效 在时深关系调节中 通常有两种可能 整体移动 Bulkshift 或拉伸压缩 Squeeze stretch 这里详细介绍具体操作 1 整体移动 Bulkshift 启用manual adjustment apply 在 well top Spreadsheet 表中出现一个TWT picked值 在合成记录主界面重新调整 welltop 直到合适的位置 此时时深关系会整体移动 Apply 后产生更新合成记录及时深关系 见图10 图11 注意 manual adjustment启用后同一welltop不会再出现两条线 一旦welltop spreadsheet 中出现TWT picked值 该welltop在 wellsection 界面上显示为实线 2 拉伸压缩 Squeeze stretch 启用manual adjustment 在well section中调节多个welltop 在well top Spreadsheet中出现多个TWT picked值 此时的 welltop将整井段分成几段 调节某个welltop 受其影响井段相对拉伸或压缩 见图12 Apply后更新合成记录 结论 在Petrel 的Synthetics模块中制作合成记录 调节时间是通过调整welltop实现的 如果工区中没有welltop 可以生成虚拟的 welltop 从而实现时深关系的调节 一旦创建合成记录 所对应的时深关系会自动保存为合成记录文件夹下的时深表 且作为当前时深关系 时深关系调整生效 后 其相应的时深关系表自动调整 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 构造面上叠合显示其它属性的等值线 有时需要将物性 例如孔隙度等 与构造信息叠合在一起显示 这样可以直观地发现构造和储层物性都比较有利的部位 这里介 绍一种方法 实现这种不同属性的叠合显示 一 创建物性的等值线数据 这里以孔隙度属性为例 假设已经有了某一层段里的孔隙度层面数据 如图1所示 打开该层面的Settings界面 在Style面板下 点击Create Contours按钮 创建该层面独立的等值线数据 如图2所示 会创建出独立的等值线数据 存放在Input窗口里 属于Polygon的一种 具体类型为Contours 如图3所示 二 给等值线赋Z值 将深度层面的Z值赋给孔隙度等值线 如图4所示 这样该等值线就可以沿着构造面显示 而等值线形态可以表现物性变化 三 成果展示 显示构造面 在该面的Settings界面里隐藏等值线显示 同时显示分离出的孔隙度等值线 颜色设为黑色 该等值线可以在深度 面的上下两侧看到 就如同常规显示的等值线效果 如图5所示 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 动态计算显示沿层不同时间段的地震属性 如果能够动态显示某一个层段范围内不同时间 深度的沿层地震属性 将对我们分析目的层范围内的地震属性变化 推测岩性物 性变化很有帮助 这一过程可以通过Petrel的自动计算即Workflow模块 调用沿层地震属性提取和层面计算功能 自动实现 这一这里介绍详细的操作步骤 一 复制原始层面数据 首先将原始层面数据复制保存以防损坏数据 需要的话可以复制后的层面改名 给出一个易于区别的名字 以备后续计算使用 并将其显示在3D窗口 二 创建Workflow 创建一个Worklfow 自动运算和显示数据 见图1 语句详细解释如下 第一句 添加一个循环 Loop 并指定循环的次数 第二句 插入层位计算Z Z Constant 并给出Constant的值 这里的Z是原始层面的时间 深度值 Constant是要动态改变的 层面的时间 深度增量步长 通过增量和循环控制要研究的目的层范围 所以 Loop次数和Constant增量要根据自己研究的目标 范围决定 第三句 调用Surface attribute功能 提取感兴趣的地震属性 具体参数见图2 第四句 调用Pause Duration功能插入一定显示的时间间隔 以方便用户浏览数据 第五句 Endloop 结束循环 整个流程通过动态变换目标层的位置 并参考目标层提取地震属性 从而实现动态扫描感兴趣层段内某一地震属性的变化 三 运行Workflow 点击图1中的Run按钮 运行该工作流 结果会显示在3D窗口 并按照所设定的时间间隔 Pause Duration 不断更新 结果 如图3所示 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 自适应河道建模方法介绍 从Petrel2007开始引入了一个新的基于目标体的河道建模方法 自适应河道建模 即Adaptive channel 自适应河道建模方法 的主要目的是使河道相的模拟与井数据更加吻合 更符合河道的实际分布 这里详细介绍一下该方法的主要思路 一 工作流程 1 解释索引体曲线 首先 在Well section window里 创建一条新的相曲线 该曲线的主要目的是指示哪些井上的河流相是属于一个河道的 从而 做为河道建模的控制 指导数据 激活well correlation process 创建一条新的离散曲线 即点击Create discrete curve的按钮 选择Bodies做为新曲线的模 板 如图1所示 则在Well目录的Global well logs目录下产生一条新的离散相曲线 索引体曲线 如图2所示 将该曲线显示在Well section窗口里 解释该数据 操作同普通沉积相解释 解释依据是对照原始河流相解释曲线 经过对比分 析 将所有井上属于同一个河道的相对应解释为同一个body 如图3所示 将索引体曲线和河流相曲线都显示在Well section窗口里 图中每口井上左边的一道是要解释的索引体曲线 右边 一道是河流相曲线 图中将Well20和Well22上的最上边的一套Channel sand相对应的索引体曲线都解释为Body1 那么意味着 后期如果以河道相数据建模时 这两个井上的浅层的Channel sand模拟出的河道将会连接在一起 是一条河道 一次类推 可 以控制不同井上 不同期次发育的河道之间的连接关系 指导河道建模 2 测井曲线粗化 启动Scale up well logs模块 分别对索引体曲线和河道相曲线进行粗化 3 自适应河道建模 启动Facies modeling进程 基于粗化后的河道相曲线进行相建模 选择Object modeling做为建模算法 在Facies bodies面板下点击 Add a new adaptive channel 按钮 添加一条自适应河道 第二个按钮是原来的河道建模方 法 并设置参数 需要设置的参数为 1 打开Use body选项 2 Facies项下选择对应的河道相 例如Channel sand 3 选择 Use Output Body Property as input for connecting channels 选项 4 同普通河道相建模一样去Layout Section和Trends项中去设置参数 如图5所示 打开Backgroud面板 设置背景相 打开Other output面板 启动Body property选项 点击三角箭头 从下拉菜单上选择上 一步粗化后的索引体曲线 如图6所示 该数据会呼应Use body选项 从而使所选择数据参与河道相建模的控制 4 结果展示 自适应河道建模算法得到两个结果 一个是河道相的分布 一个是索引体的分布 二 算法拟合步骤及数据影响 自适应河道建模算法中数据拟合顺序 1 已有数据 例如井上的河道相和索引体数据 2 河道几何形状参数 3 相分布比例 即Facies proportion参数 因此 使用该算法 启用体索引数据做为河道连接关系的指示 将增加已有数据 拟合结果与不使用体索引数据大不相同 图8 展示的是使用体索引数据进行控制和不进行控制的结果比较 左为没有启动 Use output body property as input for channel connection 选项的 右为开启该选项的结果 如果使用了体索引数据进行控制 那么修改体索引数据的解释结果会直接影响自适应河道建模的结果 请看图8 1 8 2 使用两条体索引曲线制河道相建模 得到的结果也不一样 图9中左边是使用图8 1的体索引曲线得到的结果 可以看到 红色的 body同时穿过Well20和Well2 右边是使用图8 2的体索引曲线得到的结果 在红色body发育的地方 同时又出现了浅蓝色的 body体 三 Adaptive channel与fluvial channel对比 1 相同点 二者都使用随机算法 可以产生不同的建模实现 realization 河道的分布都受原始数据 例如井上的河流相解释曲线 和河 道几何形状等参数控制 2 不同点 Adaptive channel加入 relative sinuosity 控制河道的相对变形度 使得河道的形状更符合自然界河道的真实情况 调整该参数 可以模拟出不同形状特点的河道 例如辫状河 曲流河或废弃河道 Adaptive channel可以通过体索引数据 body index 定义河道间的连接关系 使结果更容易控制 符合地质观点 吻合井 上的河流相解释结果 四 两种方法使用条件 1 Adaptive channel 适合数据丰富 地质情况明了 且有河道连接信息需要模拟的情况 2 Fluvial channel 适合数据稀少 没有河道连接信息 河道形状规则 变化较小的情况 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 属性建模的层面控制技巧 在做属性建模时 有时用户不一定掌握充分的数据 或者没有足够好的井上 地震数据来控制属性的分布 这样可能无法分析出 很好的变差函数变化规律 或者模拟出的结果不能很好地体现地质专家的分析观点 那么在Petrophysical Modeling中使用协 克里金面数据控制就是解决此问题的一个有效尝试了 这一方法的主要思路是创建一个多边形 即Polygon 给该多边形里的每条线赋不同的值 该值是通过各种地质分析得到的 还 句话说是地质家希望某处某属性 例如孔隙度应该出现的值 使用Make Edit surface功能将Polygon数据做成属性面 1 创建多边形 使用Make Edit polygons模块在需要做控制的区域画出任意线 注意一个Polygon名字下可以有多条Polygon线 2 给每条任意线赋值 将鼠标设置为选择状态 左键点击一条polygon激活 然后点击按钮 激活赋值选项在Z 之后的空 格里输入数字 然后点击左边蓝色的 按钮为Polygon赋值 如图1所示 给Polygon赋值 其Z值设定为0 05 粉色 0 15 蓝色 以及0 25 绿色 用以控制孔隙度属性最大值不会超 过0 3 3 创建面 启动Make Edit surface进程 以多边形作为主输入 Main input 来做Surface 此多边形的边界作为油藏控制边界使用 如 图2所示 检查计算完的Surface的值 在该数据的Settings statistics界面里 查看数据范围 如果有值在0 0 0 3范围以外 可以通过 Calculator用如下公式截除 最终结果如图3所示 3 使用创建好的Surface参与Petrophysical Modeling 层面值与要模拟的数据之间的相关系数可以通过Coefficient参数进行调整 点击Estimate按钮软件将会自动估算二者的相关 性 建模结果见图5 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 Petrel2008 2009按Zone进行连续属性统计的workflow 摘要 此 workflow 可帮助你很容易的对于每一个连续性属性 Continoues property 按Zone进行统计 并将结果显示到一个表单 中 总体上有三个步骤 1 创建一个新的连续性属性 continuous property 2 给定 workflow的输入 3 运行 workflow 并检查结果 过程细节 Statistics to output sheet 在 workflow 编辑器中没有设计到对应于任何的过滤 Filters 这使得你必需手动检查每一个 Settings面板的 Statistics tab 才得以获得Zone内统计信息 此workflow同样是利用这语句实现按Zone 批量统计的效果 此 workflow的核心思路是构建一个暂时性的属性 temporary property 确保其中只有一个Zone的数据 然后就可以利 用 Statistics to output sheet 语句来获得统计结果了 其它一些语句主要是设计把统计结果显示到一个报表内 三个步骤如下 1 创建一个新的连续性属性 continuous property 有很多方式可以创建 例如 a 计算器 Property calculator 可以使用计算器计算一个新属性 例如General Statistics 语句为 General statistics U 如图1 b Petrophysical modeling 也可以通过Petrophysical modeling 用赋值的建模方法创建新属性 如图2所示 c Geometrical modeling 也可以使用Geometrical modeling创建新属性 如图3所示 2 给定输入 首先要把离散性属性 discrete properties 移至另一个文件夹 比如说叫做 Facies 在程序的开始如图4所示设置四个输 入 完整的工作流见附录的图5 1 5 2 3 运行 运行Workflow 得到类似下表的结果 附录 完整Workflow示意图见后面的图5 1 5 2 由于语句比较多 所以分两张图片显示 合在一起是完整工作流 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 Petrel 中实现地震反演的方法和技巧 介绍 地震反演的目标 在井控范围以外使用地震数据来描述岩性和流体属性 所需井曲线 声波和密度 确定性反演的步骤包括 使用合成记录对井曲线进行标定和校正 精确的时间 深度标定 确切的子波提取或者描述以及井曲线 和地震数据的对应关系 叠前和叠后的比较 叠后的结果是AI 声波阻抗 信噪比比较好 不能解决入射角相关的属性问题 叠前的结果是泊松比 拉梅系数等 针对入射角相关属性 可提取剪切模量信息 Petrel遗传算法反演 遗传反演是一种结合神经网络算法和遗传算法的反演技术 主要用于通过数据间非线性关系解析获取声波阻抗信息 其所需的输 入数据为 地震振幅体和测井曲线 其优势在于 运算快 误差小 可以应用于多种地球物理属性 人工神经网络 ANN 是一种模拟人类神经细胞行为的算法 它的运行基于复杂的神经元关联网络 遗传算法 基于DNA链的运行原则 基因可以从父本向子本作转移和交叉互换 数学运行框图如图5所示 工作流程示意如图6所示 基本工作流程框图见图7 实例分析 步骤一 从GeoFrame中输出Petrel格式数据包 Basemap send export to Petrel Geology Office tools export well data to Petrel 步骤二 安装 GeoFrameDataConnectorForPetrel2009 1 msi 插件对应功能显示如下 步骤三 输入压缩包 步骤四 检查数据 步骤五 合成记录标定 步骤六 时间 深度转换 步骤七 反演 步骤八 结果评价与分析 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 Petrel中计算某个层段的砂泥比或砂层厚度 在做项目研究时 有时关注目标层段的砂泥比或者想统计砂层厚度 下面介绍如何实现这一目的 1 准备岩性曲线 外部加载 手工解释 Calculator计算或者神经网络判别 准备好岩性曲线 2 计算地质层段的砂泥比属性 鼠标右键点击某个地质分层目录 即Well tops目录 Insert new attribute 如图1所示 在弹出的界面里 类型设置为 continuous ok 如图2所示 在弹出的界面里Info面板里如图所示 给属性一个名字 template可以选择为Percent或者Fraction apply 勿关 闭此界面 如图3所示 继续到Attribute operations面板里 计算范围设置为To the zones at level 1 选择第一步计算出来的相曲线 计 算公式选择成Percentage或者Fraction 选择按照MD或者SSTVD作为计算深度的参考 选择您关心的相 点击Run进行计算 点击OK关闭此界面 3 计算地质层段的砂层厚度属性 重复以上的步骤 插入一个新的仍然是连续型的属性 此次算法Average method设置成Thickness 将相选择成砂层 则会统 计出砂层的厚度 类似地 用户还可以根据自己的需要统计其他关心的属性 4 属性的应用 通过以上操作可以获得井点处关心的地质层段里的砂泥比或者砂层厚度 有时需要根据井点值做平面预测 此时可以使用 Make edit surface模块进行计算 如图4所示 Main Input选择某个关心的地质层段 在Welltops目录的Stratigraphy目录里可以找到 点击Attributes旁的 下拉三角 找到需要的属性 设置合适的参数进行计算即可 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 对某个地质层段内的测井曲线做交汇和直方图分析 在做测井曲线分析时 有时需要对某个关注的地质层段内的数据进行分析 例如仅对某两个地质分层之间的测井曲线进行分析 下面介绍如何实现这一功能 一 插入地质层段曲线Zone log 鼠标右键点击某个地质分层目录 即Welltops目录 从菜单上选择Insert update zone log 如图1所示 结果如图1右侧图所示 会在Input面板的Wells目录的Global well logs目录里生成一条zone曲线 该曲线会提示是与那个 welltops目录关联的 注意 一个Petrel工区里 可以有多个welltops目录 即多套分层方案 每个welltops目录对应只能有一条zone曲线 第一次执 行以上操作时会为某分层方案生成zone log曲线 当分层方案做调整时 再执行该操作会对其zone log进行更新 二 创建过滤器 如图2所示 在Input面板里 找到叫Filter folder的目录 右键点击Filter folder Create 1D filters 在弹出的过滤器Settings界面里 在Info面板下 给过滤器起个名字 在Definition面板里 Object to filter处 鼠标左键点击 Input面板里的wells目录 然后点击蓝色箭头将Wells目录发送进来 点击下拉箭头选择第一步创建好的Zone曲线 在列表里可 以看到该Zone曲线包括的地质层段 接下来需要对哪个或哪几个zone进行分析 就从列表上将其选中 成蓝色高亮显示 三 使用过滤器 1 直方图窗口分析 如图3所示 打开一个直方图窗口 显示某口井的某条测井曲线 左图 此时是全井段数据的统计分析 然后选择刚才创建的 过滤器 则此时高亮显示的 是该井中 该地质层段内曲线的分布 2 交汇图窗口分析 打开一个Function窗口 交汇某两条测井曲线 类似直方图窗口 选择应用创建好的过滤器 双击过滤器 在Settings界面的 Style面板下 可以设置只显示当前zone里曲线的交汇结果 图4 1 或者是设置成将当前zone里曲线交汇结果高亮显示 图4 2 此时若对交汇的数据拟合相关性方程 也是仅对过滤后的数据计算的方程 四 修改过滤器条件 如图5所示 打开刚才创建好的过滤器的Settings界面 到Definition面板下 选择其他的zone 然后点击Apply或OK 这样 就修改了过滤器的条件 可以对其他感兴趣的zone做分析 当然用户愿意的话 也可以根据步骤二 对每个关心的zone创建一个过滤器 对不同的窗口使用不同的过滤器进行数据分析 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 使用地震剖面照片进行二维地震解释 如果手头只有地震剖面截图或照片 而并非实际的地震数据 但是需要进行二维地震解释的时候 根据不同的解释数据用途 可以在 Petrel中有三种方式来实现 根据解释数据的不同用途 1 如果解释数据用于归档数据库 那么断层可以直接在地震剖面截图上解 释 而层位解释必须基于真是的地震数据 所以需要使用插件Blueback 来将图片转换为解释数据 2 解释数据用来在Petrel中创建 三维构造模型 断层可以直接在图像上解释 或解释为多边形 层位可以解释为多边形 3 如果解释数据是用来为其他软件生成输 入数据 如IGEOSS Dynel 3D Dynel 3D需要地震解释数据作为输入数据或者使用构造三角网格 如 ts文件 此处介绍如何在 petrel里不用插件进行解释 解释步骤分解如下 第一步 输入Bitmap图片 1 使用Bitmap格式加载图片 如图1 2 设置选项右上角Independent edges保持depth Z 垂直 如图2 第二步 插入general intersection进行解释 1 在任意文件夹右键插入general intersection 2 选择三点确定general intersection的位置 第三步 使用Make edit polygon开始解释层位 1 创建一个 pseudo 的interpretation filter来区分不同的polygon的图像来源 如图3 2 使用append polygons的功能将属于同一层位的polygon合并成一个 如图4 第四步 移动Intersection面板 1 一旦移动了intersection 设置正确的视角 此处设置为west如图5 2 将前一图像解释的层位当作 neighbor 也显示在解释窗口协助 如图6 3 解释该图像 如图7 第五步 进行断层解释 1 因为断层解释独立与地震道 可直接激活Seismic interpretation进程进行解释 图8 2 重新命名 第六步 使用Make edit polygon进行断层解释 与第五步互为替换 1 创建 pseudo 的interpretation filter来区分不同断层的图像来源 图9 2 重新命名 第七步 使用Make edit surface生成断面 1 使用第五步生成的polygon作为输入数据生成断面 图10 11 2 使用Convergent interpolation作为插值算法 第八步 生成3D Grid 仅供选择 1 在Fault Modeling下直接用fault polygons或者fault interpretation生成断层模型 2 在Make Horizon进程下用解释的horizon polygon 或surface 转换 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 在Petrel中合并多个地震数据体 有时由于受到地震体数据大小的局限或是从工作站中连续切割出线道号连续的地震体需要在加入到Petrel之后进行合并 在 Petrel当中可通过简单几步完成连续线道号不同地震体的合并 第一步 在seismic main folder下创建一个survey 第二步 先加载第一个单个segy文件 并编辑survey的geometry 既是线道号的起始位置 如图1 第三步 分别加入其它剩余segy文件 如图2 第四步 加载完成之后 选中右键合并 如图3 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 如何切割属性模型 在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方向的网格的范围 注意不要点击Get from filter 点击Copy global grid 此时会在Models面板里看到 Part of 模型的名字 这时候你打开3Dwindow可以查看你选取的部分 模型 如图6所示 3 第二种方法是用边界选取部分网格 首先在3D或者2D window中画出你想要选取部分网格的polygon边界 如图7所示 应用Geometrical Modeling进程 建立一个新的模型 建模型的时候选择Assign between surfaces and polygons 在下面 的polygon处放入刚才编辑好的边界polygon 如图8所示 点击Apply 这时在已有模型中的Properties下会看到你刚做的Elevation depth模型 选中显示在3D窗口 可以发现其为Polygon边界内的 模型 如图9所示 然后设置过滤器 双击Properties 打开Settings for Properties 在Filter tab下 选中value filter 在窗口最下面的value filter中 将刚才用polygon边界建立的Elevation depth选中 右侧选中Use filter 此时左侧的Elevation depth模型显示为紫 色 说明已经应用过滤器 如图10所示 点击Apply OK 关闭此窗口 鼠标左键双击或者右击打开要选取部分模型的属性模型的Settings 如图10所示 在Output tab下 选中Use segment filter Use zone filter Use cell index limits 此三个过滤器可以任意组合使用 当选中Use cell index limits后 在此不需要在下面 的IJ选取栏内选取要显示的IJ方向的网格的范围 直接点击Get from filter 然后再点击下面的Copy global grid 此时会在 Models面板里看到 Part of 模型的名字 这时候你打开3Dwindow可以查看你选取的部分模型 如图11所示 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 Petrel中如何在Map window下显示不同layer的属性 在Petrel中 在map window下如何查看不同layer的属性 1 准备属性模型 2 新建一个map window 鼠标左键点击选中要显示的属性模型 比如说bulkvolume 此时mapwindow中显示为Bulk volume这个属性模型 Petrel默 认将此属性模型的第一个layer的属性显示在map window中 有时第一层网格比较少时在map window中不容易看到 如图2所示 在map window中只显示第一层的网格属性 如图3所示 双击或者右击Properties 打开Settings for Properties 在Style tab的最下端 可以看到Property layer to show in mapping选项 在K后面输入要显示的layer号 Apply 即可在Map window中看到要显示的layer的属性 3 在Map window中检查输入layer的属性值 如图4所示 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 在Stereonet中显示断层 有时我们需要查看断层的倾角和方位角 我们可以用下面的流程把断层显示在stereonet window里 1 把断层解释数据转到深度域 如图1 2 把深度域断层转成点 如图2 3 Create fracture network 以转成点的断层为输入 选择deterministic方法 如图3 4 在Stereonet Window里显示相应结果 如图4 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 地震解释窗口中Ghost的应用 应用一 叠合显示多种属性 1 点击Ghost图标 画矩形圈定范围 确认在Windows栏Ghost被激活 如图1 2 在Input勾选另一属性体 注意不是Inline或Xline 如图1显示的为Envelope属性 此时如果移动窗口或改变窗口大小 可以点击右键更新对应数据 3 Windows栏双击Ghost 在settings里可以设置Ghost窗口的透明度 右键菜单可以删除Ghost 应用二 辅助地震解释 Ghost还可以用来进行断层上下盘同相轴的对比 与邻近剖面的对比 帮助确定解释的层位 此时 Ctrl Shift可以旋转Ghost窗 口 Shift RMB恢复垂直状态 如图2 应用三 在Ghost窗口里显示井曲线 通过压缩和拉伸Ghost窗口 可以与地震剖面之间进行初步的对比 如图3所示 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 断层建模中如何较好的处理复杂接触关系的断层 在做构造模型的过程中 遇到断层发育的油藏 断层交接关系又很复杂 即使勉强处理的断层接触关系 在网格化之后 网格的 质量很难保证 往往容易形成畸形的网格 给数模的网格质量也很不理想 在处理这种问题的时候 petrel2012给我们提供了一 个好的解决办法 避免了上述的诸多问题 工作流程如下 1 在Structural framework下 双击Geometry definition进程 图1 创建一个新的几何体 从地震数据体得到工区范围 选择对应的深度域或时间域 点击ok 图2 2 双击Fault framework modeling进程 弹出窗口如下 图3 从input中输入所有的解释断层 设置相应的参数 点击 ok 则自动生成含有所有断层相互连接关系的断层格架 图4 在3D窗口中查看断层之间的交接关系 对接触关系复杂又处理 不好的断层 记下它们的断层名字 3 展开Corner point gridding文件夹 双击Fault model from structural framework进程 图5 在这个进程界面下 对第二 步中记下的断层Closing Fault East Closing Fault South 在include标签下对应的断层名字处 不打勾 让它们不参与计 算 设置相应的参数 点击ok 图6 即对处理不好的复杂断层先不让它参与形成断层模型 以保证整体断层模型的网格质 量 然后按照常规的构造建模流程 进行pillar gridding horizons 一系列的小层构造模型建完之后 开始添加复杂的断层进入到成型的构造框架中 采用阶梯状网格算法 保证复杂断层附近的网 格不发生扭曲变形 双击Stair step faulting进程 图7 弹出界面窗口中即为没有参与断层模型建立的复杂断层 图8 在 Stair step标签下勾选 点击ok 针对复杂接触关系的断层 通过建立阶梯状断层 就把矛盾极小化了 难以处理的断层通过这种方式插入到已经形成的高质量小 层构造模型中 按照断层的位置插入到网格中 并兼顾网格的分布 保证了网格的质量 图9 图10 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 构造与岩性共同控制的油藏中如何设置多个油水界面 对于构造和岩性双重控制且油水关系复杂的油藏 在petrel中往往比常规单一的油藏做油水界面要复杂难处理一些 碰到类似的 情况 可以采取下面的工作流来解决 1 在2D窗口下 显示构造层面 选择合适的构造等值线间距 如果有确定的油水边界文件 可以输入到petrel中 如果没有这 类数据 可以激活make edit polygons进程 按照分析的确定性成果 在构造图上圈取油水边界线及岩性控制的边界线 图1 双击make edit surface进程 弹出界面下输入多边形文件 选取Artificial algorithms 赋对应的Z值 生成油水界面surface文 件 对于构造控制的断块 surface对应一个油水边界的海拔深度值 对于岩性控制的断块 surface对应岩性边界所在的最低构 造线值 最后形成了多个油水界面surface文件 双击一个surface文件 弹出settings界面下 双击make contacts进程 弹出界面下输入此surface文件 就可以求取相应层位的油水界面了 图3 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 Workflow制作等比例沿层时间切片 在三维地震体中 可以任意插入Time slice 在三维空间内对任意水平时间切片进行查看 通过Surface attribute也可以提取沿 层的地震属性 然而对于某一特定地层段 获知等比例地层切片需借助Workflow进行实现 从而更好的帮助用户分析与地层等 时沉积有关地震属性的变化特征 1 数据准备 在Input面板下 依据地震体 解释地震层位 创建新的文件夹 并命名为 surface attribute result 用于存放最终生成等比例 时间切片 2 原理 先计算目标层面的厚度差 按一定的比例划分 求得单个小层的厚度面 将目的层顶面构造按照等比例层厚下移 利用surface attribute 逐次提取不同深度相应的地震属性值 最终获得等比例沿层地震属性切片 3 创建Worklfow解释 Step1 将地震解释层面结果插值生成surface 变量A B分别代表顶 底的插值结果 Step2 Copy 顶面surface 将其中一个命名为Thickness between two surface Step3 Thickness between two surface 变量C 由顶面和地面surface计算求得两层面厚度 Copy 该厚度面得到新的层面 变量D Step4 将Step3生成的厚度进行n等分 命名等分后的厚度面为 Thickness Pro Step5 依据等比例划分的个数n 建立n个循环计算步骤 逐次将地层顶部构造面等比例下移一定距离 获得新的surface 变 量E 并复制产生一个新的层面 变量F Step6 加载显示的地震体 以每次下移后的构造面为参考地层 选择Extract value 提取邻近地震道的数值 振幅值或其他属 性值 选择合适的波形位置 生成结果赋值给已存在的面 变量F Step7 依次更改新生成的面为surface attribute 1 3 n 将最终结果移动到surface attribute result的文件保存 结束循环 最后 点击Test和Run按钮 运行该工作流 在3D Window等窗口下 查看surface attribute result中产生的等比例不同位置 沿层切片 返回首页 返回成功案例 Petrel地震地质解释和建模使用技巧 基于地震数据体的属性建模方法 在井资料稀少的地区 合理的利用地震资料是对井间岩性 物性进行解释的有效方法 利用Geobody模块 通过交汇基于遗传 反演算法得到波阻抗体与主频体 定义不同岩性的值区间 最终得到由能够区别出各种岩性的离散型地震数据体 Train Estimation Medol主要基于神经网络算法 能够模拟人脑 识别属性曲线与地震体之间的规律 最终完成对孔隙度 渗 透率等属性的运算过程 算法优势 1 不依赖三维网格模型的相建模 属性建模技术 2 神经网络算法保证计算结果具有更高的准确性 3 无网格引入和粗化过程 4 快速 高效的计算过程 1 提取地震属性 依据地震数据体和测井曲线 提取多个与岩性 物性有关的属性体 如遗传反演AI体 主频体 构造平滑体 Envelope 瞬时 频率 甜点等属性体 图1 2 Geobody体解释 激活地震体 Gen inversion 点击Insert box probe 插入三维Probe体 通过Manipulate Probe选择合适的体积大小 参 考Grid Model 目的层位 在Input 面板下的Geobody Interpretation Probes文件夹里自动产生一个Probe体 右键 Setting 在弹出的对话框里选择Volumes 在2nd Cube里导入第二地震属性体 如 Dominant Freq 体 点击对话框下 方的Apply键 切换至Opacity子界面 由默认的Histogram切换到Crossplot 在右侧区域选择某一特定颜色 点击Tools下的第 一排按钮 在两种属性的交汇图上 勾画特定的区域范围 例如 砂岩通常具有较高的波阻抗 呈现低频特征 观察3D Windows地震体的变化 按住Ctrl键 切换另外一种颜色 选择泥岩等其它岩性 最终 将连续的地震数据体转变成离散的岩相 数据体 图2 3 将Probe体转换为地震体 右键Probe box体 点击convert to seimic cube 将解释的Probe体转换成可以为Grid Model采样的标准地震数据体 离散数 据 图3 右键转化的岩相地震体 点击Create well seismic 在Well Section界面 对比解释结果与原始测井曲线 图4 在Geometrical Modeling里 利用Seismic Resampling 将地震体重采样到激活的网格模型里 调整为相模板显示 图5 4 利用Train estimation Model生成孔隙度 渗透率模型 在Input Wells Global well logs Checkshot Attribute TWT picked 右键双击TWT picked 选择Convert to log 在 Input面板Global Well Logs下 生成一个TWT 曲线 右键Global Well Logs 图6 点击Calculator 键入 Seismic log name S