landmark分频处理流程PPT课件.ppt

分频处理工作流程 Landmark 1 本章主要讲述分频处理的主要工作流程 以及解释和分析分频数据 我们相信经过本章工作流程的描述 会使你更清晰的了解储层的精细模型 在你获得清楚的认识之前 你需要明白分频处理中的一些参数问题 对于一个新手来说 我们推荐你从第二步 二维模块分析 开始 第二步从模型数据开始 演示各种参数对输出效果的影响 这将有助于我们刚接触的新手对分频处理理论进一步的理解 注意 流程中使用的一些数据在CD中 请注意版本说明如何去加载一般的数据 real data 是Landmark训练类数据体的命名样式 通过模型的操作 初学者会看到叠加道和多样化地形地层是如何干扰影响分频中的离散F变的 在所有的流程中 你会知道如何去拿捏这些参数 诸如一个单形体的时窗的选择 都会大大的影响分频处理的结果 当你对模型数据的处理感到很惬意的时候 你就有能力进行你的实际勘探数据的参数化 并很快得跑完整个流程 具体的流程如下 数据准备 初测调谐 二维模块 三维流程时窗 目标调谐图 离散频率体 储层厚度估算 2 数据准备 DataPreparation 3 输入的地震数据体须经精细的处理 去建立正确的子波相位和真实的道间振幅信息 一个定态的子波对分频处理是至关重要的 地表一致性处理和二维的滤波去噪 会加强你的分频效果 输入数据的带通滤波器也许会消除一些地质信息 而这部分信息有的却是你准备去分解的 如果可能尽量用没有未滤波的偏移数据作为输入去处理 虽然分频处理是在振幅和相位域操作的 但是最好选择16道32位的数据 因为8位的数据不能做细节分析 虽然可以做初步的调谐工作 在具体处理之前 你需作到所选择的层必须经过内插而且要贯穿所有断层 频谱分析通常能使断层更加明了 这使你能重新去认识这些断层 在频谱成像上 如果你的层位具有相对高的信噪比和稳定的地震波同相轴 还得考虑对数据作剩余静校正 这一步将清除许多由于工作站上自动追踪的同相轴所引起的跳跃现象 你也可以通过简单的圆滑地层来处理 频谱分析前的这步会敏锐地提高成像质量 这种提高能使你恰好得到地质信息的一个重要部分 输入数据须适当的偏移 分频处理假定已经校正倾斜地层和菲涅耳带的影响 简单的入射模型能够决定你的砂体是否发育在目的层段 假设某个砂体可以成像 然而 不能说她就是有形的 对于各种地震技术都显示的砂体 她必然有一个突出的波阻抗差异 从声波曲线的分析到合成地震记录形成 你可以测定你的砂体是不是有形的 这一步回答了这个问题 你的数据准备好了吗 4 初测调谐 ReconnaissanceTuningCube 5 在单位厚度 不连续叠加的层和断层上 最优参数化的调谐体拥有最大的几率显示出露变化 细微的道间变化无论在厚度或者叠层上 都会使目的层段的数据元素调协成不同的分立的频率成分 高频部分显示细微的 低分辨率的断层 低频部分则趋向于去噪和展示宏观结构和断层趋势 这步流程给予一个目的层的调谐的过程的浏览 她将回答如下问题 你看到目的层段了吗 流程主要分三部分 一 Parameterization 参数化 使用一个完全的长时窗 100 200ms 中心在目的层段层上 有10到20 的余弦递变 在确定可以获得频率从DC到Nyquists后 使用一个宽泛的频率范围让你去限定在解释种有用的信号带宽 二 Animation 制作 通过启动制作调谐体 观察并最终决定带宽 也就决定这相干的地质信息能否被观测到 三 Selection 选择 选择一对相对独立的频率体成像显示当前的目的层地质地形 如果在这个流程中没有相干的地质信息被观测到 那么接着 SpecDecomp工具就不能被用来加深你对目的层层位的解释 理解 如果这个流程KO了 你就能够最优化你的调谐参数 作进一步的详细分析 6 快速浏览流程一 1 在启动OpenWorks中选择Application 模块 SpecDecomp 分频 TuningCube 调谐体 2 选择输入的项目名称如tns train 3 选择输入的数据体到seimic cmp CMP是由EJB容器自动完成对数据库的操作 你所有做的 就是在实体bean重写入SetXXX或getXXX方法 然后在ejb jar xml中定义cmp field4 设置分析时窗 基准层到实际层位 5 改变cosine 余弦 梯度有20 左右的递变 6 给予数据体前缀名 注意输出的幅度体被命名位Workflow 1 0 DFT AMP bri DFT 离散F变 表征着变换的方法 Amp表示体属性为振幅体 SpecDecomp输出的地震数据体都被写成程序块格式 7 7 默认其他参数并选择需计算数据体 计算会给予你5分多钟的休息时间 具体在于你机器的性能 如果你的机器不是那么快 你可以去小憩一会 8 8 确定计算完成 计算每道DFT 过滤 标准 归类频率切片 制作 现在你预备去启动地震波剖面 选择一些分立的频率对目的层段地质地形成像 1 运行SeisWorks3D 选择tns train项目 2 打开SD wf 1 0 1 t w3s时域 3 选择解释人员 程序 不要井和断层数据 4 从剖面上看 选择File NewTask Seismic打开一个新地震剖面显示窗口 9 5 改变输入地震数据设置为Workflow 1 0 DFT Amp bri 并重新选择时间频率至0到250Hz 频率谐和体剖面图 实际工作中上图是很难被确切解释的 然而你能注意到在35Hz左右有很多大的幅度差异现象 注意 在 seismicview 地震波剖面图上 35Hz显示的是35ms 同时可以看到 70Hz以下的数据都显得像噪声 10 6 从新的地震剖面窗口 选择Seismic SelectfromList 并选择T0的时间切片 这实际上是0Hz的频率切片 7 改变 colorbar 色标到blkwht 自行调节 8 选择View FrameControl 在seismicview工具条中 9 单击FrameControlforward按纽图标 启动制作调谐体 11 10 制作停止在30Hz 你可以看到一个清楚的河道轨迹 当你在这种类型的显示图上看到了非地质或者是很线性化的信息时 你的原始数据可能有点问题 那样的话 你必须检测你的输入层位是否完全内差并覆盖目的层段 当然你还需要去确认你的地震数据是不是符合数据准备流程的要求 注意 Motion 引导线 是一个非常强有力的工具 解释人员可以用来在地震勘探中提取数值 你可以启动制作整个数据体来探测细微的变化 这些变化在静态地震剖面图上是很难被检测到的 在频谱分析中 你可以发挥motion的作用去解释储层各向异性和厚度的复杂变化 12 11 向前并复查从0到折叠频率的数据体制作 如果你想去控制制作速度 你需要在EarthCube或者OpenVision中去设置 如果数据的通频带宽在60Hz左右 你就基本上得不到60Hz以外的信息 当你操作你自己的未滤波的数据时 你会发现有用信息的频率是从0到折叠频率的 13 选取 从上面的分析中 你会意识到频率从10到40Hz包含了有价值的地质信息 也就是你想进一步作研究的这部分 14 为了充分理解分析这些调谐体的成像 你必须明白砂体越厚 振幅谱的V形切痕 凹口 就越多 注意 振幅谱上 楔形模型的最厚部分有12个切口 最薄部分仅就一个 如果你通过楔形模型 从低频到高频 来制作调谐体 你会发现幅度等值线是从厚到薄 在下章中将具体分析楔形模型 你能注意到等值线从河道中心漂移吗 Page34当你在从10到40Hz制作的时候 15 二维模块 2DScope outAnalysis 16 一旦你已经通过初测调谐 你可以最优化频谱分析参数 这步解决这个问题 时窗和 taper 梯度的关系 最优参数化流程的目的是为你提供一些参数 这些参数是用作覆盖整个区块的三维时窗分析流程 这步流程还用来作训练流程 为了大家更好地了解频谱分析的细枝末节 这一步一般只能作横截面上一个独立单元的分析 二维断面测试方法效果不错 因为时窗的移动范围小 并且DFTs的迭代算法可以很快 还可以在原地震剖面上看到是否有对地质现象突出的频率响应 流程分二步 选择一个横截面作分析 使用一个递变范围和分析时窗的长度来生成和观测频谱分析数据尽管流程是分为二步的 我们将指导您通过一个楔形模型分形的详细分析来加强你对频谱分析的理解 17 数据质量控制 先提供一个模型数据的基本情况 楔形模型在下面相应地作出分析 瞬时厚度 1 从一个地震剖面窗口 改变输入地震数据设置为wedge3dv 2 从地震剖面窗口 选择Seismic SelectfromList 并选择线号1131 重新选择时窗900到1150ms 2 你现在可以看到这个楔形的模型经过滤波整形 注意到瞬时厚度怎么没有楔形模型上那么尖呢 你也许想改变一下显示的比例 18 3 高亮显示theTop Wedge Interpreted和base Wedge Interpreted地层 注意为什么楔形解释的精度不能和实际的模型相匹配呢 地下地形因薄层调谐而发生扭曲 如果你把没有带限的地震楔型变换到频率域 你就能看到如下所示的图象 注意振幅谱上切口会随着厚度的增加而增多 这是由于地层厚度大余1 4 时引起的干扰和破坏 高频成分短波长 也要看到没有频率能使0厚度的地层成像 不能分辨尖灭现象 19 上面图中的白色条带代表最主要的频率 在下一节中你会建立一个相似的振幅谱分析图 第三步在带限的地震剖面上 初始调谐分析 本节指导我们在带限楔型上去贯通一个调谐体的过程 1 选择Applications SpecDecomp TuningCube从Openworks平台上 2 选择tsp train项目 3 设置输入体为wedge3dv 4 把开始结束线号调到1131 5 改变分析时窗 间断型到ConstantWindow 连续型 6 设置基准时间为1025ms 7 选择最大时窗长度 使用MB3在时窗长度输入栏中 注意你可以使用MB3在输入各项 可以得到可用输入值的列表 8 最大频率改为250 9 将余弦梯度给定为20 20 10 为Workflow 2 1设置数据体前缀 11 选择OutputHorizons 12 选择Autocreate 自动生成 如下层位 13 完成后关闭输出设置 21 14 默认其他参数选择ComputeVolume 计算数据体 15 关掉ComputationsCompleted信息 这步一分钟就搞定了 16 显示1131线在地震剖面图上Workflow 2 1 DFT Amp bri 重新设置列表作新的brick文件 22 17 你需要再选择一次时窗范围在地震剖面图上 从0到250Hz 注意 为什么通频带下面的所有数据都成了噪声 也还要看到噪声在72Hz左右每10Hz成一条带 你能想到是什么引起的吗 也许是归一化引起这些条带 因为这个问题出现在模型的1560到1570道之间 23 18 地层演示 1025 ConstantTime DFT FisrstPeak Freq 凭直觉你能看到已调谐地层的高振幅部分 得注意在此图上色标是反转的 是不是FirstPeak Frep层看来就像是沿着首峰的轨迹而来的 移动光标沿着并观看显示图底部振幅值的输出 你会看到一些振幅的峰值在计算出的FirstPeak Frep层之上 为什么呢 首峰的计算是正确的 问题在于你的显示是在归一化数据上的 因此说明FirstPeak Frep层是在归一化前计算出来的 归一化后 所有的频率切片被提高到同一个幅度水平 归一化的细节具体讨论在P68 首峰计算在P136 24 为了分析数据 研究楔型体的瞬时厚度 t 是如何决定振幅谱伴随频率的削减周期 Pf 就很重要 同样的 你需要知道频率成分 f 的数值是如何决定振幅谱伴随地层厚度的削减周期 Pt 在你作任何有关厚度估算流程前 必须充分理解这些关系 25 调谐厚度的数学表达式如下 举例来说 一个10Hz的子波可能有50ms的调谐厚度 20Hz的子波的调谐厚度为25ms 30Hz的子波的调谐厚度为17ms 注意 Widess 1973 和Wood 1982 论述了一个求取调谐厚度的方法表达式 26 10 20和30Hz的调谐厚度 前面计算的 其带限振幅谱和频率域图示图下 你需要确定你明白这2条振幅和频率示意曲线的相互关系 如30Hz时 0振幅在2条曲线上的交叉点在哪儿 34ms 27 你也许会得到这样的结论 如果我简单的反演最主要的频率 就会得到相对的调谐厚度 那就是说 主频越高 地层越薄 对于初步的分析 你可以使用地震层位估算并反演xxx DFT FirstPeak Freq层 就可以获得一个相对的调谐厚度图 流程6 1指导你进行这个计算的细节 然而 你必须保持这样的观点 20ms的瞬时厚度 相对于更新世沉积的砂体 在8000f s图示为80ft 当20ms相对泥盆纪石灰岩时 在18000f s图示为180ft 就像我们所知道的 地质学会变得很容易 假使我们所打的井是时间意义上的而不是深度的时候 另一个分析得出的结论是 你花在主频分析上的时间越长 你的解就越向惟一靠拢 28 二维时窗分析 在区块输出参数测试中 我们推荐你从观测二维线的DFT剖面开始 你首先要通过使用一个窄的窗口 或者 你可以使用梯度减小有当前的时窗 聚焦一个小的时窗 缩短时窗或拉长梯度对DFTs的输出效果是相似的 其效果就是DFTs在垂向上 在频率的方向上 明显得圆滑了 圆滑后 能很容易地清除缺切并高度了联合振幅谱分析 因而 一些目的层段的地质地形在35Hz下的 可以显示在100ms窗口中 当在160ms时窗中观测时 P2 50 因此 我们建议 windowcarpentry 窗口工序化 在选择任何分立的频率作进一步研究之前 1 在地震工作中 显示Workflow 2 1 DFT Top Workflow 2 1 DFT Middle Workflow 2 1 DFT Bottom分析窗口在wedge01 3dvseismicview上 29 30 2 在SpecDecomp RuningCube中 改变窗口大小从218到100ms 这是小时窗的最好输出 如果小时窗不能显示任何地质信息 则逐渐增大其大小直到你看到相干的同相轴为止 接着 采准这个时窗直到很多冗余信息的出现 如果你有时间 你可以从30ms时窗开始 按每次增加10ms一直增加到200ms 3 关闭并归一化切片 4 默认其他参数开始计算 按ComputeVolume 如果你选择输出的数据体是已经存在的 你需准备是否删除或者覆盖原来的体 5 选择覆盖 合成文件 当你测试 分析你的调谐体的时候 这个覆盖和删除按纽使你能够使用规范化易于管理的输出数据设置 警告 虽然覆盖和删除按纽是很有用途的 你使用的时候还是必须很小心 尤其是 当你不能够再生成特殊的切片 你正在地震工作中观测 的时候 删除足以引起问题 请参考日志文件去解决已有的问题 31 6 演示 刷新加载在地震工作中的数据体 你会发现默认的显示比例对没有归一化的数据是不能正确显示的 7 从地震演示上 选择Seismic Parameters 变化DynamicScale Clip 动态比例按纽 从NoScaling到Sample 15 1 这些与218ms分析的相比会怎么样呢 它会不会更圆滑 首峰地层看起来象不象确切地计算出来的 地层噪声条带时够已经被清楚 你能想到在左边的演示中当楔型mod变得非常薄的时候会发生什么变化吗 32 8 重复操作生成一个50ms的分析窗口 这些与218和100ms的时窗分析比较又会怎么样呢 请把注意力集中到第一条谱带密集区的凹口处 当凹口特别明显的时候 你已经对你的时窗分析参数做到了最优化 9 使用这个程序去测试不同的递变项 10 继续测试直到你对所设置的参数满意为止 33 二维教学模型 下面的实例在SpecDecomp的安装盘上 这些实物是用来让你测试不同的参数如何反应不同的地质条件 请查看发布的日志文件来加载这些模型 wedge01 3dv是在前面练习中使用的楔形模型 这个模型涵盖了40条线 1111到1150 和460道 1110到1561 时窗范围是900到1150ms 而频率范围是从10 20Hz 1111线 到10 100Hz 1150线 wedge01 3dv 对于wedge01 3dv 尝试生成并阐明下面的展示图 34 chanv201 3dv是从1101到1229线的一个河道模型 时窗2800到3200ms 这个模型从薄到厚随着线号的增加 chanv201 3dv chanv401 3dv是一个河道叠加模型 就一条线1195 是chanv201 3dv的3个拷贝的混成 时窗从2900到3100ms 所有道都有不同的振幅 chanv401 3dv 35 thinch01 3dv是一个抽道集河道模型从1101道1229线 时窗是2800到3200ms 它和chanv201 3dv的道集模型厚度递减是一样 所以一般看来都是薄的 thinch01 3dv chx2v101 3dv是一个由chanv201 3dv和thinch01 3dv叠加的河道模型 从1101道1229线 时窗是2800到3200ms chx2v101 3dv 36 chx2v201 3dv是两个浅河道模型的集合 线号为1196 时窗2900到3100ms chx2v201 3dv thinv201 3dv thinv201 3dv是浅河道的振幅变化模型 但没有薄层 线道号是1175和1420 时窗2800到3200ms 37 thinv301 3dv浅河道模型 从thinch01 3dv 上面有一个强同相轴 线道号是1175和1420 时窗2800到3200ms thinmd01 3dv是尖峰模型用来形成浅河道模型 线道号是1175和1420 时窗2800到3200ms 38 fubar401 3dv是 wedge01 3dv1130线的静态频移 相位旋转 振幅衰弱的变型 这个模型线号1130 时窗为900到1150ms 除fubar401 3dv和wedge01 3dv是1ms采样 其他的 3dv都是2ms采样 39 三维时窗 3DRunning windowAnalysis 40 对实际数据来说 与模型数据是有很多不同的 大多数勘探需要在可变分析时窗上进行3D分析 这步流程主要解决以下问题 对多重叠加的产层怎么处理 如何与地质现象同步 能观测到泛滥平原的分层面吗 信道被时间约束的地方在哪儿 能否看到其他的地震相界面 时窗分析的最佳位置点在什么 我们在开眼界之前 就可以看到流程带来的地形地貌 此流程有3个主要步骤 Selecting 选取 选择一个区带并在研究区上下适度地延伸一段距离 如果你的工区位于一个倾斜的地层 在生成运行时窗分析前需要把这个层拉平 Generating 生成 生成运行时窗分析 频率成分从初测调谐体中选取 覆盖全区 与初测调谐体相比 这里对每个时间点的局部分析是很重要的 一个80ms的时窗使用高斯分级可以提供所需的局部化 Animating 演化 通过得到的数据体演化测定 区带间的信号渗滤 例如在已知现象作用前 你移动窗口的时距不再显示在图象中 41 工区是由实际的层组成的 层是需要拉平的 为了消除分析中倾斜干扰的可能性 在这个流程中 你仅分析那些显现最高的信噪比的频率 视具体情况使用10 20和30Hz 和初测调谐体定义的一样 流程一 选择Tools VolumeRecon从调谐体菜单栏 这个项目会被自动设置为当前的项目 如果必要的话 设置输入数据体为real data cmp 在时窗范围选项中 改变时窗范围为FlattenToHorizon 对于层位 点击List 按纽 选择实际层位 默认输出采样频率为2ms 改变中心窗口定义层的时间长度为200ms 生成初测调谐体 在程序选项中 设置时窗为100ms 设置变换函数为Gaussian 高斯变换集中于窗口中心 P120有各种变换的对比 设置体前缀为WF 3 1 Gaussian 100 体前缀中 100 标志着时窗长度 Gaussian 标志变换类型 42 设置频率列表 10 20 30 默认其他参数并选择ComputeVolume 选择OK到ComputationsCreateWarning 计算将花费几个小时 就看你工作站的CPU 内存和硬盘传输速度 本工区在800MHzSun Blade 1000工作站上花费了40分钟 清除完成信息 43 演化初测调谐体 加载WF 3 1 Gaussian 100 Amp 030 Hz DFT bri数据体到地震时间切片观测图上 你需要重新选择地震时间 调谐体才会被显示出来 点击 演化整个分立的能量体 分立的能量体在30Hz时的窗口 就像你在地质时期上下旅行的演化一样 0ms表示输出体实际层的新位置点 尝试用各种色标 blkwht color max Sopacity2 和Sspecdecomp 你看到有趣的地方了吗 请注意边滩 1260线 1170道附近 它什么时候出现和消失的呢 你又怎么解释这种现象呢 44 演化10到20Hz的调谐体 当你演化整个 运行时窗分析能量体 时 你须要能够分选出在什么样的级数上最能涵盖你的目的层段 你变更的时窗在什么时候不能使你看到主要目标的效果 在任何变更范围内底图和切片的调节会降低分辨率吗 你演化的时候 你可以去阐明其地质时期发生的变化 例如 你可以看到由主河道发展的分叉河道的现象吗 你可以看到边滩在地质年代上的发展吗 当你通过地质时期去演化时 你也能看到河道充填在其上的变化 当你演化一个河道的时候 你首先看到一个暗黑色条带在河道下面的中心上 然后 暗黑色条带相称性地向河道两边转动 河道沉积充填简单模式 检索Geology101会发现大多数的河道种砂体 沉积物水流较缓的边缘 因此 砂体一般都是中间薄由于河流的快速流动 能量体的演化于调谐体不一样 你看到同样厚度的河道转化为不同的频率切片 要紧记调谐体的演化在时间上是冻结的 窗口一直不动 运行窗口种分立的能量体演化冻结的只是在单一的频率上 而且它的谐和图可以在截面上上下移动 更重要的是 当观测复杂储层的成像时 地质学家 地球物理学家和储量工作者以及现场工程师可以看到在各自独立的学科上的直接效应 地震勘探队伍就能够尽快而协同地影响到全局 45 当你有时间的时候 可以重复整个流程 测试不同的时窗的效果 当你有时间的时候 可以重复整个流程 测试不同的变换梯度 在本流程的最后 你也许对高保真度三维频谱成像 通过分频处理 已经有了一个粗略的概念 三维高保真度来自于直接对反射系数的成像 如子波的影响被清除 而且抽样的反射系数包是独立于倾角的 这个流程输出的一些参数是为时窗分析和Taper设置的 下面列举一些合理的参数 AnalysisWindow ReferenceHorizon real datahorizon ReferenceHorizonTimeShift 0ms WindowLength 100ms Taper 20 CosineTaper 时窗分析 参考层位 实际层位 参考层位时窗变更 0ms 时窗长度 100ms 变换 20 余弦递变 未拉平 当进行分析的时候 你没有必要一定拉平层数据 然而 拉平有两个主要的优点 可以快速生成小数据体 可以在拉平的地质年代上演化切片 不拉平时跑完整个流程 你首先要运行VolumeRecon从100ms到730ms 因为输入层的范围是200到730ms而你还需要一个上下100ms的延伸 然后你可以在SeismicView种选择Horizons Flatten 交互式得拉平你输出的横截面 不幸的是 地震工作不允许在时间切片上层拉平 因此 你需要在演化你的时间切片前 使用叠后资料拉平输出的数据体 在这个流程中归一化是不一样的 相邻未拉平道集会影响归一化的数量 P65 46 交互式拉平 使用拉平数据的一个问题是 它使拉平数据和未拉平数据对比起来非常困难 如何去克服 一个办法是在地震转换中使用持续拉平 当然这种转换要求能同时显示拉平和未拉平的地震数据体 这个流程发表在July August2002UserNetarticletitle SpectralDecompositionPersistentFlatteningWorkflowforVolumeRecon你可以从两方面获得 这篇文章的另一个AdpbeAcrobat版本包含在 SPECDHOME docs文件夹里叫SpecD PersistentFlattening Workflow pdf 你可以打开AdpbeAcrobat查看这个文件 47 目标调谐图 TargetedTuningMapAnalysis 48 目标调谐图分析 这步也叫做目标调谐体 此流程假定频率范围 时窗分析和变化都已经定义好了 这个流程能回答如下问题 哪儿才是钻井的亮点 流程输出包括 调谐体和子波叠加分离 用叠加补偿性频谱化的调谐体制作全带宽的振幅图 在子波被分离之后 信号带宽综合部分经常可以生成更好的概括性的断层及其边界图 不同的断层和边界用不同的频率成分 一般高频时比较清晰 来成像 总和 叠加 调谐体 在清除叠加的子波后 允许每个频率成分给予一个权重 因而当强化边界和断层的时候 由于厚度的可变性 会导致叠加出的调谐效应 然而 这存在一个冒险 这种叠加中你可能会抹掉一些频谱异常 如果你要成图的地形恰好是这样的 调入的大于这一半 而调出的大于另一半 叠加就会严重的消除这些异常 在这种情况下 单独的窄带 调谐图可能会提供优秀的图象 分析时窗层位的上 中和下 First 频谱最大层 由工区中波长和多边形形成的瞬时厚度图 主要切片 调谐图 保存为地震属性层位 49 基础 清除子波叠加 调谐体由3部分组成 薄层干涉 子波叠加和噪声 既然地质响应是解释人员最感兴趣的部分 所以人们都很谨慎在没有降低地质信息的情况下来抵消子波的幅度 这一步将调谐体缩减为薄层干涉 噪音 50 通用频谱均衡技术依靠稀疏不变量稳态统计学 那就是 先为所有点假定一个连续的子波 如果在拉平的层上呈现出地质调谐多样化 然后你可以通过每个频率切片的中间振幅值作均衡补偿来抵消子波波谱 归一谱白化以后 子波叠加效应缩到最小程度 调谐体保留两个部分 薄层干涉 噪音 51 在调谐图上 频率切片 地层中 薄层干涉是作为相干振幅的变量 随机噪声斑点化干涉的样式就像低质量的雪花点电视带来的感觉一样 在主频上 相对高信噪比 SNR 带来了相对清晰的薄层调谐的图片 在有用带宽外的频率上 低信噪比则带来了雪花电视 生成目标性调谐体 流程4 0 这步和初测调谐体的生成 流程一 是相似的 那是确定你的研究工区是否具有有利的地质信息 然而 本流程目标调谐体在运行时 是有最优的参数输入的 如果你对3D时窗分析进程的分析 流程三所依助的参数比流程一中的快速浏览差异较多的话 你就必须要再运行流程1 0 使用更新过的参数 生成目标性调谐图 流程4 1 一旦生成了一个正常的调谐体 就伴随着一些细节分析 应用于各个频率切片上 储层厚度估计就是这么一个流程 可以在频率切片上显示 然而 这个流程 需要把精力集中在调谐图对地震解释层位的转换以及精确的频率图计算上 选择Application SpecDecomp TuningMapper从蓝玛上启动 选择spt train项目 设置输入的数据体为Workflow 1 0 DFT Amp bri 设置切片范围在RangeList上为10 20 30 40 在MergeSlices下 打开Sum并生成一个输出层位叫作WF 1 0 Merged Frequencies 设置输出层位前缀未 WF 1 0 WF 1 0 表示输入体的来源 52 点击ComputeHorizonsfromSlices 点击OK 对HorizonCreateWarning 清除thecomputationscompleted信息 53 在地震解释地图上显示WF 1 0 030Hz层 注意看清晰的河道特征 54 调整色标 只显示河道 55 调谐成图解释 流程4 2 作为一个解释人员 一个主要的任务就是生成图件并计算圈闭面积 蓝玛的Auto PolygonGenerator工具就是用来成图和面积计算的 而且是更加有效更为精确的 在SeisWorksMapView上 选择Horizons Auto PolygonGenerator 这个工具将 计算所有polygon的面积 生成适合在EarthCube或OpenVision中观测的3Dpolygons 生成适合在Z MAPPlus中成图的区域对照 控制 图 参数化的Auto PolygonGenerator显示如右图 56 点击ProcessHorizon后 你会等到如下 左图 的Completed信息 观察生成的Polygons的数量 以及最大的20个的面积 再运行一次Auto PolygonGenerator作业 将MinimumPolygonArea acres to100 新的作业将只计算出7个Polygons 57 在SeisWorksMapview中显示这些Polygons及其面积 确认并关闭指示层 58 向前 关闭Auto PolygonGenerator对话框 在地震横截面 你可以开始审阅河道部分如何变化的情况 在SeismicView显示real data bri的1175线 在SeismicView只打开Polygons层 这将显示河道的全局范围 如果观看Polygons层对你来说有困难 请选择View Contents并选上HorizonIntersectionCircles 现在显示1475线 注意 在3D中查看Polygons 你需要进入EarthCube和OpenVision 流程6 1在EarthCube2003 3 Help Workflows下描述了一个功能不大的3D分析选项 59 在MapView和SeismicView之间使用光标追踪你所看到河道的地方 60 如果你还有时间 重复流程4 2的步骤 作出10 20 30和40Hz的Polygons图 10 20 30和40Hz的Polygons图能显示河道随频率的变化吗 这或许是一个暗示 暗示着河道厚度的不协调 调谐图合成 流程4 3 WF 1 0 Merged Frequencies演示图 注意 为什么频率合成在一起后 移除演示图上一些河道信息 回到TuningMapper 不选ExtractSlicesasHorizons并改变切片范围从Range0到80 确认其他参数的匹配性 参考流程4 1 61 保持其他的参数不变并点击ComputeHorizonsfromSlices 当完成调谐成图工作后 刷新MapView 演示修改后的WF 1 0 Merged Frequencies图 这个演示显示地震切片到分立的频率切片的强大的分解能力 通过分解后 你可以避免来自未调谐频率成分的破坏性的干扰 62 RGB超览看图器 流程4 4 超览允许你推算你的河道的RGB演示 在一张图上 你可以把RGB演示当成3个不同的估值 如果须要打开PowerView Tools PowerView从TuningMapper上 通过选择File LoadData或者点击DataLoader释放DataLoader 如下所示 点击层位栏 选择WF 1 0 020Hz WF 1 0 030Hz WF 1 0 040Hz层 点击 transferarrow 移动所选择的层位到当前列表中 63 从主窗口 选择Tools ImageComposer弹开如下的对话框 64 点击RGBBlending栏在下面的边框中 再点Preview栏 你的混合图将显示在PowerView主窗口中 演示图中20 30和40Hz分别对应红 绿和蓝 在图上显示为白色的地方 所有的频率都有一个相似的较高的地质响应 显示为单色的地方 红 绿和蓝 仅其对应的频率在区域内已调谐 这样的话 RGB混合的演示图把效果调谐到单一的图上 你可以叠加RGB演示 形成一个和分立频率能量体相似的动画 你可以通过更改最大最小值或者他们的滑动条 修改各种颜色演示的振幅范围 当你对演示图满意的时候 输入图名并点击CreateComposite按纽 65 离散频率体 DiscreteFrequencyVolumeInterpretation 66 离散数据体解释 检测 在目的层分析之后进行 调谐图把问题限定在了一个局部的目的层段的范围上 而表征大型地震数据体时 则需要另外不同的方法 对于超出单反射数据包或目的层段的分频处理 我们推荐你使用 Discretefrequencyenergycubes 离散频率能量调谐体 这些能量的数据体 本质来说是整个地震数据体的一个三维时窗分析 用连续的时窗 离散频率体是由一个单一的地震数据体计算到多元的离散频率 振幅和相位的数据体 通过一个滚动谱分析时窗 为每个模块计算出其振幅或者相位 这些谱分量会被归类成相同频率成分的模块 这种方法通常是再分析及参数化目的层后进行 67 本流程解决如下问题 还有其他新的用来钻井的井位吗 如何提高断层的解释精度 流程将产生在单一频率下的振幅信息的数据体 对体的解释 流程5 0 此流程与三维时窗分析一样 不同的是一个频率需要计算一个数据体 举个例子 你可以围绕产层 产油层 通过休整输入的初始数据 进行分析 并运行 VolumeRecon 仅选30Hz 68 储层厚度估算 ThicknessEstimation 69 频谱分析提供一个独立相位的稳定的方法来估算地震上的储层厚度 借助分频处理 即使是一个相对低的频率成分 如10Hz 都可以定量分析薄层的变化情况 在油田上 任何定量分析的关键都在于有井控制 正象地震反演技术一样 分频处理厚度分析需要井控来保证最终结果的合理性 从块状模型薄层分析的理论来看 振幅谱的周期于频率Pf存在V型切口的关系 而频率Pf与薄层的时间厚度t有关 如下 t 2 z V t 1 Pf 其中z等于薄层厚度 V等于速度 可求得时间厚度z z V 2 Pf 直达波的频率 f0等于1 2Pf 将这个关系式代入厚度等式 你会得到 z 厚度 V 4 fo t 双程旅行时 1 2 f0 很多时候 基于频谱分析的薄层厚度分析 可以确保薄层厚度与有效砂岩厚差在1米的范围内 而当传统的厚度分析产生歧义或者不能解决问题的时候 分频处理的厚度分析就突显了它的作用 显然地 伴随着反射系数模型复杂性和综合性的不断增加 它必将可以完善振幅谱的解释 并结合其相互的联系 这项研究正在进行中 本流程解决如下问题 储层的相对厚度是多少 70 储层相对厚度图 流程6 1 SpecDecomp的 R2003 0 版 主要集中于定性分析 分析结果将显示地层是厚是薄 期望的有利井区等 但