landmark分频处理工作流程

1、分频 处理工作流程1 本章主要讲述分频处理的主要工作流程,以及解释和分析分频数据。我们相信经过本章工作流程的描述，会使你更清晰的了解储层的精细模型。 在你获得清楚的认识之前，你需要明白分频处理中的一些参数问题。对于一个新手来说，我们推荐你从第二步（二维模块分析）开始。第二步从模型数据开始，演示各种参数对输出效果的影响。这将有助于我们刚接触的新手对分频处理理论进一步的理解。 注意：流程中使用的一些数据在CD中，请注意版本说明如何去加载一般的数据。“real-data” 是Landmark训练类数据体的命名样式。 通过模型的操作，初学者会看到叠加道和多样化地形地层是如何干扰影响分频中的离散F变的。

2、在所有的流程中，你会知道如何去拿捏这些参数，诸如一个单形体的时窗的选择，都会大大的影响分频处理的结果。 当你对模型数据的处理感到很惬意的时候，你就有能力进行你的实际勘探数据的参数化，并很快得跑完整个流程。 具体的流程如下： 数据准备初测调谐二维模块三维流程时窗目标调谐图离散频率体储层厚度估算 2数据准备（Data Preparation）3 输入的地震数据体须经精细的处理，去建立正确的子波相位和真实的道间振幅信息。一个定态的子波对分频处理是至关重要的。地表一致性处理和二维的滤波去噪，会加强你的分频效果。输入数据的带通滤波器也许会消除一些地质信息，而这部分信息有的却是你准备去分解的。如果可能尽量

3、用没有未滤波的偏移数据作为输入去处理。虽然分频处理是在振幅和相位域操作的，但是最好选择16道32位的数据，因为8位的数据不能做细节分析，虽然可以做初步的调谐工作。 在具体处理之前，你需作到所选择的层必须经过内插而且要贯穿所有断层。频谱分析通常能使断层更加明了，这使你能重新去认识这些断层，在频谱成像上。如果你的层位具有相对高的信噪比和稳定的地震波同相轴， 还得考虑对数据作剩余静校正。 这一步将清除许多由于工作站上自动追踪的同相轴所引起的跳跃现象。你也可以通过简单的圆滑地层来处理。频谱分析前的这步会敏锐地提高成像质量。这种提高能使你恰好得到地质信息的一个重要部分。 输入数据须适当的偏移，分频处理假

4、定已经校正倾斜地层和菲涅耳带的影响。简单的入射模型能够决定你的砂体是否发育在目的层段。假设某个砂体可以成像，然而，不能说她就是有形的。对于各种地震技术都显示的砂体，她必然有一个突出的波阻抗差异。从声波曲线的分析到合成地震记录形成，你可以测定你的砂体是不是有形的。 这一步回答了这个问题：你的数据准备好了吗？4初测调谐（Reconnaissance Tuning Cube）5 在单位厚度、不连续叠加的层和断层上，最优参数化的调谐体拥有最大的几率显示出露变化。细微的道间变化无论在厚度或者叠层上，都会使目的层段的数据元素调协成不同的分立的频率成分。高频部分显示细微的、低分辨率的断层；低频部分则趋向于去

5、噪和展示宏观结构和断层趋势。 这步流程给予一个目的层的调谐的过程的浏览，她将回答如下问题： 你看到目的层段了吗？ 流程主要分三部分： 一、Parameterization（参数化） 使用一个完全的长时窗（100200ms），中心在目的层段层上，有10到20%的余弦递变。在确定可以获得频率从DC到Nyquists 后，使用一个宽泛的频率范围让你去限定在解释种有用的信号带宽。 二、Animation（制作） 通过启动制作调谐体，观察并最终决定带宽；也就决定这相干的地质信息能否被观测到。 三、Selection（选择） 选择一对相对独立的频率体成像显示当前的目的层地质地形。 如果在这个流程中没有相干

6、的地质信息被观测到，那么接着，SpecDecomp工具就不能被用来加深你对目的层层位的解释/理解。如果这个流程KO了，你就能够最优化你的调谐参数，作进一步的详细分析。 6快速浏览流程一 1.在启动OpenWorks中选择Application(模块)SpecDecomp(分频)Tuning Cube (调谐体) 2.选择输入的项目名称如tns_train。 3.选择输入的数据体到seimic.cmp(CMP是由EJB容器自动完成对数据库的操作，你所有做的，就是在实体bean重写入SetXXX或getXXX方法，然后在ejb-jar.xml中定义cmp-field 4.设置分析时窗基准层到实际层

7、位。 5.改变cosine（余弦）梯度有20%左右的递变。 6.给予数据体前缀名。 注意输出的幅度体被命名位Workflow_1_0_DFT_AMP.bri. DFT(离散F变)表征着变换的方法；Amp 表示体属性为振幅体。SpecDecomp输出的地震数据体都被写成程序块格式。77.默认其他参数并选择需计算数据体。计算会给予你5分多钟的休息时间，具体在于你机器的性能。如果你的机器不是那么快，你可以去小憩一会。88.确定计算完成计算每道DFT过滤、标准、归类频率切片制作： 现在你预备去启动地震波剖面，选择一些分立的频率对目的层段地质地形成像。1.运行SeisWorks 3D,选择tns_tra

8、in 项目。2.打开SD_wf_1-0-1.t.w3s 时域。3.选择解释人员（程序），不要井和断层数据。4.从剖面上看，选择 TaskSeismic 打开一个新地震剖面显示窗口。95.改变输入地震数据设置为Workflow_1_0_DFT_Amp.bri,并重新选择时间频率至0到250Hz 。频率谐和体剖面图实际工作中上图是很难被确切解释的，然而你能注意到在35Hz左右有很多大的幅度差异现象。注意：在（seismic view）地震波剖面图上，35Hz显示的是35ms同时可以看到，70Hz以下的数据都显得像噪声。106.从新的地震剖面窗口，选择 SeismicSelect from List

9、。并选择T0的时间切片。这实际上是0Hz的频率切片。7.改变（color bar）色标到blkwht（自行调节）。8.选择ViewFrame Control,在seismic view工具条中。9.单击Frame Control forward 按纽图标“”启动制作调谐体。1110. 制作停止在30Hz你可以看到一个清楚的河道轨迹。当你在这种类型的显示图上看到了非地质或者是很线性化的信息时，你的原始数据可能有点问题。那样的话，你必须检测你的输入层位是否完全内差并覆盖目的层段；当然你还需要去确认你的地震数据是不是符合数据准备流程的要求。注意：Motion （引导线）是一个非常强有力的工具，解释人

10、员可以用来在地震勘探中提取数值。你可以启动制作整个数据体来探测细微的变化，这些变化在静态地震剖面图上是很难被检测到的。在频谱分析中，你可以发挥motion的作用去解释储层各向异性和厚度的复杂变化。1211. 向前并复查从0到折叠频率的数据体制作。如果你想去控制制作速度，你需要在EarthCube 或者 OpenVision中去设置。如果数据的通频带宽在60Hz左右，你就基本上得不到60Hz以外的信息。当你操作你自己的未滤波的数据时，你会发现有用信息的频率是从0到折叠频率的。13选取从上面的分析中，你会意识到频率从10到40Hz包含了有价值的地质信息，也就是你想进一步作研究的这部分。14为了充分

11、理解分析这些调谐体的成像，你必须明白砂体越厚，振幅谱的V形切痕（凹口）就越多。注意： 振幅谱上，楔形模型的最厚部分有12个切口，最薄部分仅就一个。如果你通过楔形模型（从低频到高频）来制作调谐体，你会发现幅度等值线是从厚到薄。在下章中将具体分析楔形模型。你能注意到等值线从河道中心漂移吗，Page34 当你在从10到40Hz制作的时候？ 15二维模块（2D Scope-out Analysis）16 一旦你已经通过初测调谐，你可以最优化频谱分析参数，这步解决这个问题： 时窗和（taper）梯度的关系？ 最优参数化流程的目的是为你提供一些参数，这些参数是用作覆盖整个区块的三维时窗分析流程。这步流程还

12、用来作训练流程，为了大家更好地了解频谱分析的细枝末节。这一步一般只能作横截面上一个独立单元的分析。二维断面测试方法效果不错，因为时窗的移动范围小，并且DFTs的迭代算法可以很快，还可以在原地震剖面上看到是否有对地质现象突出的频率响应。 流程分二步： .选择一个横截面作分析 .使用一个递变范围和分析时窗的长度来生成和观测频谱分析数据 尽管流程是分为二步的，我们将指导您通过一个楔形模型分形的详细分析来加强你对频谱分析的理解。17数据质量控制 先提供一个模型数据的基本情况，楔形模型在下面相应地作出分析。 瞬时厚度1.从一个地震剖面窗口，改变输入地震数据设置为wedge 3dv。 2.从 地震剖面窗口

13、，选择SeismicSelect from List. 并选择线号1131，重新选择时窗 900到1150ms。 2.你现在可以看到这个楔形的模型经过滤波整形 ，注意到瞬时厚度怎么没有楔形模型上那么尖呢，你也许想改变一下显示的比例。 183.高亮显示the Top_Wedge_Interpreted 和 base_Wedge_Interpreted 地层。 注意为什么楔形解释的精度不能和实际的模型相匹配呢？地下地形因薄层调谐而发生扭曲。如果你把没有带限的地震楔型变换到频率域，你就能看到如下所示的图象。注意振幅谱上切口会随着厚度的增加而增多。这是由于地层厚度大余1/4时引起的干扰和破坏，高频成分

14、短波长。也要看到没有频率能使0厚度的地层成像。（不能分辨尖灭现象）19上面图中的白色条带代表最主要的频率。在下一节中你会建立一个相似的振幅谱分析图，第三步在带限的地震剖面上。 初始调谐分析本节指导我们在带限楔型上去贯通一个调谐体的过程。 1.选择Applications SpecDecompTuning Cube 从 Openworks 平台上。 2.选择tsp_train 项目。 3.设置输入体为 wedge 3dv。 4.把开始结束线号调到1131。 5.改变分析时窗间断型到Constant Window(连续型). 6.设置基准时间为1025ms。 7.选择最大时窗长度，使用MB3 在时

15、窗长度输入栏中。注意你可以使用MB3在输入各项，可以得到可用输入值的列表。 8.最大频率改为250. 9.将余弦梯度给定为20%。 20 10.为Workflow_2_1设置数据体前缀. 11.选择Output Horizons. 12.选择Autocreate (自动生成) 如下层位。 13.完成后关闭输出设置。2114.默认其他参数选择Compute Volume。计算数据体15.关掉Computations Completed 信息，这步一分钟就搞定了。16.显示 1131线 在地震剖面图上Workflow_2_1_DFT_Amp.bri.重新设置列表作新的brick文件。2217.你需

16、要再选择一次时窗范围在地震剖面图上（从0到250Hz）注意：为什么通频带下面的所有数据都成了噪声。也还要看到噪声在72Hz左右每10Hz成一条带，你能想到是什么引起的吗？也许是归一化引起这些条带，因为这个问题出现在模型的1560到1570道之间。2318.地层演示（1025_ConstantTime_DFT_FisrstPeak_Freq）凭直觉你能看到已调谐地层的高振幅部分，得注意在此图上色标是反转的。是不是FirstPeak_Frep 层看来就像是沿着首峰的轨迹而来的？移动光标沿着并观看显示图底部振幅值的输出，你会看到一些振幅的峰值在计算出的FirstPeak_Frep 层之上。为什么呢？

17、首峰的计算是正确的，问题在于你的显示是在归一化数据上的。因此说明FirstPeak_Frep 层是在归一化前计算出来的。归一化后，所有的频率切片被提高到同一个幅度水平。归一化的细节具体讨论在P68，首峰计算在P136。24为了分析数据，研究楔型体的瞬时厚度(t) 是如何决定振幅谱伴随频率的削减周期(Pf )就很重要。同样的，你需要知道频率成分（f）的数值是如何决定振幅谱伴随地层厚度的削减周期(Pt )。在你作任何有关厚度估算流程前，必须充分理解这些关系。25调谐厚度的数学表达式如下：举例来说，一个10Hz的子波可能有50ms的调谐厚度：20Hz的子波的调谐厚度为25ms：30Hz的子波的调谐厚

18、度为17ms：注意： Widess(1973)和Wood(1982)论述了一个求取调谐厚度的方法表达式：2610，20和30Hz的调谐厚度（前面计算的），其带限振幅谱和频率域图示图下。你需要确定你明白这2条振幅和频率示意曲线的相互关系。如30Hz时，0振幅在2条曲线上的交叉点在哪儿？（34ms）27你也许会得到这样的结论：如果我简单的反演最主要的频率，就会得到相对的调谐厚度。那就是说，主频越高，地层越薄。对于初步的分析，你可以使用地震层位估算并反演xxx_DFT_FirstPeak_Freq 层，就可以获得一个相对的调谐厚度图。流程6.1指导你进行这个计算的细节。然而，你必须保持这样的观点：2

19、0ms的瞬时厚度，相对于更新世沉积的砂体，在8000f/s图示为80ft，当20ms相对泥盆纪石灰岩时，在18000f/s图示为180ft。就像我们所知道的，地质学会变得很容易，假使我们所打的井是时间意义上的而不是深度的时候。另一个分析得出的结论是：你花在主频分析上的时间越长，你的解就越向惟一靠拢。28二维时窗分析在区块输出参数测试中，我们推荐你从观测二维线的DFT剖面开始。你首先要通过使用一个窄的窗口（或者，你可以使用梯度减小有当前的时窗）聚焦一个小的时窗。缩短时窗或拉长梯度对DFTs 的输出效果是相似的。其效果就是DFTs 在垂向上（在频率的方向上）明显得圆滑了。圆滑后，能很容易地清除缺切

20、并高度了联合振幅谱分析，因而，一些目的层段的地质地形在35Hz下的，可以显示在100ms窗口中，当在160ms时窗中观测时。（P2-50）因此，我们建议“window carpentry” 窗口工序化，在选择任何分立的频率作进一步研究之前。1.在地震工作中，显示Workflow_2_1_DFT_Top, Workflow_2_1_DFT_Middle, Workflow_2_1_DFT_Bottom分析窗口在wedge01.3dv seismic view 上。29302.在SpecDecomp-Runing Cube中，改变窗口大小从218到100ms。 这是小时窗的最好输出。如果小时窗不能

21、显示任何地质信息，则逐渐增大其大小直到你看到相干的同相轴为止。接着，采准这个时窗直到很多冗余信息的出现。 如果你有时间，你可以从30ms时窗开始，按每次增加10ms一直增加到200ms。3.关闭并归一化切片。4.默认其他参数开始计算，按Compute Volume.如果你选择输出的数据体是已经存在的，你需准备是否删除或者覆盖原来的体。5.选择覆盖/合成文件。当你测试/分析你的调谐体的时候,这个覆盖和删除按纽使你能够使用规范化易于管理的输出数据设置。警告： 虽然覆盖和删除按纽是很有用途的,你使用的时候还是必须很小心.尤其是,当你不能够再生成特殊的切片(你正在地震工作中观测)的时候，删除足以引起问

22、题。请参考日志文件去解决已有的问题。316.演示/刷新 加载在地震工作中的数据体。你会发现默认的显示比例对没有归一化的数据是不能正确显示的。7.从地震演示上，选择Seismic Parameters 变化Dynamic Scale&Clip （动态比例按纽）从No Scaling 到 Sample*15/1。这些与218ms分析的相比会怎么样呢？它会不会更圆滑？首峰地层看起来象不象确切地计算出来的？地层噪声条带时够已经被清楚？你能想到在左边的演示中当楔型mod变得非常薄的时候会发生什么变化吗？328.重复操作生成一个50ms的分析窗口。这些与218和100ms的时窗分析比较又会怎么样呢？请把注

23、意力集中到第一条谱带密集区的凹口处，当凹口特别明显的时候，你已经对你的时窗分析参数做到了最优化。9.使用这个程序去测试不同的递变项。10.继续测试直到你对所设置的参数满意为止。33二维教学模型 下面的实例在SpecDecomp的安装盘上。这些实物是用来让你测试不同的参数如何反应不同的地质条件。请查看发布的日志文件来加载这些模型。 wedge01.3dv 是在前面练习中使用的楔形模型。这个模型涵盖了40条线（1111到1150）和460道（1110到1561）。时窗范围是900到1150ms，而频率范围是从10-20Hz（1111线）到10-100Hz（1150线）。wedge01.3dv对于w

24、edge01.3dv，尝试生成并阐明下面的展示图。34 chanv201.3dv 是从1101到1229线的一个河道模型，时窗2800到3200ms。这个模型从薄到厚随着线号的增加。chanv201.3dv chanv401.3dv 是一个河道叠加模型，就一条线1195。是chanv201.3dv的3个拷贝的混成，时窗从2900到3100ms。所有道都有不同的振幅。chanv401.3dv35 thinch01.3dv 是一个抽道集河道模型从1101道1229线，时窗是2800到3200ms。它和chanv201.3dv的道集模型厚度递减是一样，所以一般看来都是薄的。thinch01.3dv

25、chx2v101.3dv 是一个由chanv201.3dv和thinch01.3dv叠加的河道模型。从1101道1229线，时窗是2800到3200ms。chx2v101.3dv36 chx2v201.3dv 是两个浅河道模型的集合，线号为1196，时窗2900到3100ms。 chx2v201.3dv thinv201.3dv thinv201.3dv是浅河道的振幅变化模型，但没有薄层！线道号是1175和1420，时窗2800到3200ms。37 thinv301.3dv 浅河道模型，从thinch01.3dv，上面有一个强同相轴。线道号是1175和1420，时窗2800到3200ms。 t

26、hinmd01.3dv 是尖峰模型用来形成浅河道模型，线道号是1175和1420，时窗2800到3200ms。38 fubar401.3dv 是.wedge01.3dv 1130线的静态频移，相位旋转，振幅衰弱的变型。这个模型线号1130，时窗为900到1150ms。 除fubar401.3dv 和wedge01.3dv 是1ms采样，其他的.3dv都是2ms采样。39三维时窗（3D Running-window Analysis）40 对实际数据来说，与模型数据是有很多不同的，大多数勘探需要在可变分析时窗上进行3D分析。这步流程主要解决以下问题： “对多重叠加的产层怎么处理？”“如何与地质现

27、象同步？”“能观测到泛滥平原的分层面吗？” “信道被时间约束的地方在哪儿？”“能否看到其他的地震相界面？”“时窗分析的最佳位置点在什么？” 我们在开眼界之前，就可以看到流程带来的地形地貌。此流程有3个主要步骤。 Selecting(选取) 选择一个区带并在研究区上下适度地延伸一段距离。如果你的工区位于一个倾斜的地层，在生成运行时窗分析前需要把这个层拉平。 Generating(生成) 生成运行时窗分析，频率成分从初测调谐体中选取，覆盖全区。与初测调谐体相比，这里对每个时间点的局部分析是很重要的。一个80ms的时窗使用高斯分级可以提供所需的局部化。 Animating(演化) 通过得到的数据体演

28、化测定： 区带间的信号渗滤（例如在已知现象作用前，你移动窗口的时距不再显示在图象中。）41工区是由实际的层组成的，层是需要拉平的，为了消除分析中倾斜干扰的可能性。在这个流程中，你仅分析那些显现最高的信噪比的频率。视具体情况使用10，20和30Hz，和初测调谐体定义的一样（流程一）。选择ToolsVolume Recon 从调谐体菜单栏。这个项目会被自动设置为当前的项目。如果必要的话，设置输入数据体为real_data.cmp。在时窗范围选项中，改变时窗范围为Flatten To Horizon。对于层位，点击List. 按纽，选择实际层位。默认输出采样频率为2ms。改变中心窗口定义层的时间长度

29、为200ms。生成初测调谐体在程序选项中，设置时窗为100ms。设置变换函数为Gaussian。高斯变换集中于窗口中心，P120 有各种变换的对比。设置体前缀为WF_3_1_Gaussian_100。体前缀中“100”标志着时窗长度。“Gaussian”标志变换类型。42设置频率列表：10，20，30。默认其他参数并选择Compute Volume。选择OK到Computations Create Warning. 计算将花费几个小时，就看你工作站的CPU，内存和硬盘传输速度。本工区在800MHz Sun-Blade-1000 工作站上花费了40分钟 。清除完成信息。43演化初测调谐体加载 W

30、F_3_1_Gaussian_100_Amp_030_Hz_DFT.bri 数据体到地震时间切片观测图上。你需要重新选择地震时间，调谐体才会被显示出来。 点击“”演化整个分立的能量体。 分立的能量体在30Hz时的窗口。就像你在地质时期上下旅行的演化一样，0ms表示输出体实际层的新位置点。尝试用各种色标；blkwht_color_max,Sopacity2, 和Sspecdecomp。你看到有趣的地方了吗？请注意边滩（1260线，1170道附近，它什么时候出现和消失的呢？你又怎么解释这种现象呢？44演化10到20Hz的调谐体。当你演化整个“运行时窗分析能量体”时，你须要能够分选出在什么样的级数上

31、最能涵盖你的目的层段。你变更的时窗在什么时候不能使你看到主要目标的效果？在任何变更范围内底图和切片的调节会降低分辨率吗？你演化的时候，你可以去阐明其地质时期发生的变化。例如，你可以看到由主河道发展的分叉河道的现象吗？你可以看到边滩在地质年代上的发展吗？当你通过地质时期去演化时，你也能看到河道充填在其上的变化？当你演化一个河道的时候，你首先看到一个暗黑色条带在河道下面的中心上。然后，暗黑色条带相称性地向河道两边转动。河道沉积充填简单模式检索Geology 101会发现大多数的河道种砂体/沉积物水流较缓的边缘。因此，砂体一般都是中间薄由于河流的快速流动。能量体的演化于调谐体不一样，你看到同样厚度的

32、河道转化为不同的频率切片。要紧记调谐体的演化在时间上是冻结的（窗口一直不动）。运行窗口种分立的能量体演化冻结的只是在单一的频率上，而且它的谐和图可以在截面上上下移动。更重要的是，当观测复杂储层的成像时，地质学家，地球物理学家和储量工作者以及现场工程师可以看到在各自独立的学科上的直接效应。地震勘探队伍就能够尽快而协同地影响到全局。45当你有时间的时候，可以重复整个流程，测试不同的时窗的效果。当你有时间的时候，可以重复整个流程，测试不同的变换梯度。在本流程的最后，你也许对高保真度三维频谱成像（通过分频处理）已经有了一个粗略的概念，三维高保真度来自于直接对反射系数的成像（如子波的影响被清除），而且抽

33、样的反射系数包是独立于倾角的。这个流程输出的一些参数是为时窗分析和Taper设置的。下面列举一些合理的参数：.Analysis Window Reference Horizon real_data horizon.Reference Horizon Time Shift 0ms.Window Length 100 ms.Taper 20% Cosine Taper 时窗分析 参考层位 实际层位参考层位时窗变更0ms时窗长度100ms变换20%余弦递变未拉平当进行分析的时候，你没有必要一定拉平层数据。然而，拉平有两个主要的优点： 可以快速生成小数据体。 可以在拉平的地质年代上演化切片。 不拉平时

34、跑完整个流程，你首先要运行Volume Recon 从100ms到730ms。因为输入层的范围是200到730ms而你还需要一个上下100ms的延伸。然后你可以在Seismic View种选择Horizons Flatten.交互式得拉平你输出的横截面。不幸的是，地震工作不允许在时间切片上层拉平。因此，你需要在演化你的时间切片前，使用叠后资料拉平输出的数据体。在这个流程中归一化是不一样的，相邻未拉平道集会影响归一化的数量。（P65）46交互式拉平使用拉平数据的一个问题是：它使拉平数据和未拉平数据对比起来非常困难。如何去克服，一个办法是在地震转换中使用持续拉平，当然这种转换要求能同时显示拉平和未

35、拉平的地震数据体。这个流程发表在July/August 2002 UserNet article title:Spectral Decomposition Persistent Flattening Workflow for Volume Recon 你可以从两方面获得。这篇文章的另一个Adpbe Acrobat版本包含在$SPECDHOME/docs 文件夹里叫SpecD_PersistentFlattening_Workflow.pdf.你可以打开Adpbe Acrobat查看这个文件。47目标调谐图（Targeted Tuning Map Analysis ）48目标调谐图分析 这步也叫

36、做目标调谐体。此流程假定频率范围，时窗分析和变化都已经定义好了。这个流程能回答如下问题：“哪儿才是钻井的亮点？”流程输出包括：. 调谐体和子波叠加分离.用叠加补偿性频谱化的调谐体制作全带宽的振幅图（在子波被分离之后）信号带宽综合部分经常可以生成更好的概括性的断层及其边界图。不同的断层和边界用不同的频率成分（一般高频时比较清晰）来成像。总和（叠加）调谐体（在清除叠加的子波后）允许每个频率成分给予一个权重，因而当强化边界和断层的时候，由于厚度的可变性，会导致叠加出的调谐效应。然而，这存在一个冒险，这种叠加中你可能会抹掉一些频谱异常。如果你要成图的地形恰好是这样的，调入的大于这一半，而调出的大于另一

37、半，叠加就会严重的消除这些异常。在这种情况下，单独的窄带，调谐图可能会提供优秀的图象。.分析时窗层位的上、中和下.（First）频谱最大层.由工区中波长和多边形形成的瞬时厚度图.主要切片“调谐图” 保存为地震属性层位 49基础：清除子波叠加 调谐体由3部分组成：薄层干涉，子波叠加和噪声。既然地质响应是解释人员最感兴趣的部分，所以人们都很谨慎在没有降低地质信息的情况下来抵消子波的幅度。这一步将调谐体缩减为薄层干涉，噪音。50通用频谱均衡技术依靠稀疏不变量稳态统计学。那就是，先为所有点假定一个连续的子波。如果在拉平的层上呈现出地质调谐多样化，然后你可以通过每个频率切片的中间振幅值作均衡补偿来抵消子

38、波波谱。归一谱白化以后，子波叠加效应缩到最小程度，调谐体保留两个部分：薄层干涉，噪音。51 在调谐图上（频率切片）地层中，薄层干涉是作为相干振幅的变量。随机噪声斑点化干涉的样式就像低质量的雪花点电视带来的感觉一样。在主频上，相对高信噪比（SNR）带来了相对清晰的薄层调谐的图片。在有用带宽外的频率上，低信噪比则带来了雪花电视。生成目标性调谐体 (流程4.0) 这步和初测调谐体的生成(流程一)是相似的，那是确定你的研究工区是否具有有利的地质信息。然而，本流程目标调谐体在运行时，是有最优的参数输入的。如果你对3D时窗分析进程的分析流程三所依助的参数比流程一中的快速浏览差异较多的话，你就必须要再运行流

39、程1.0，使用更新过的参数。 生成目标性调谐图 (流程4.1) 一旦生成了一个正常的调谐体，就伴随着一些细节分析，应用于各个频率切片上。储层厚度估计就是这么一个流程，可以在频率切片上显示。然而，这个流程，需要把精力集中在调谐图对地震解释层位的转换以及精确的频率图计算上。 选择Application SpecDecomp Tuning Mapper 从蓝玛上启动。 选择spt_train 项目。 设置输入的数据体为Workflow_1_0_DFT_Amp.bri。 设置切片范围在Range List上为10，20，30，40。 在Merge Slices 下，打开Sum并生成一个输出层位叫作WF

40、_1.0_Merged_Frequencies. 设置输出层位前缀未“WF_1.0” “WF_1.0”表示输入体的来源。52点击Compute Horizons from Slices。点击OK,对Horizon Create Warning。清除 the computations completed 信息。53在地震解释地图上显示WF_1.0_030Hz层。 注意看清晰的河道特征。 54调整色标，只显示河道。 55调谐成图解释 (流程4.2) 作为一个解释人员，一个主要的任务就是生成图件并计算圈闭面积。蓝玛的Auto-Polygon Generator 工具就是用来成图和面积计算的，而且是更

41、加有效更为精确的。 在SeisWorks Map View上，选择Horizons Auto-Polygon Generator. 这个工具将：计算所有polygon的面积；生成适合在EarthCube 或 OpenVision 中观测的3D polygons；生成适合在Z-MAPPlus中成图的区域对照（控制）图。 参数化的Auto-Polygon Generator 显示如右图：56点击Process Horizon 后，你会等到如下（左图）的Completed信息。 观察生成的Polygons的数量，以及最大的20个的面积。 再运行一次Auto-Polygon Generator 作业，

42、将 Minimum Polygon Area (acres) to 100 。新的作业将只计算出7个Polygons。57在SeisWorks Map view 中显示这些Polygons及其面积。确认并关闭指示层。 58向前，关闭Auto-Polygon Generator 对话框，在地震横截面，你可以开始审阅河道部分如何变化的情况。在Seismic View 显示real_data.bri的1175线。在Seismic View 只打开Polygons层。这将显示河道的全局范围。如果观看Polygons层对你来说有困难，请选择View Contents 并选上 Horizon Inters

43、ection Circles。现在显示1475线。注意：在3D中查看Polygons ，你需要进入EarthCube 和OpenVision。流程6.1 在EarthCube2003.3 Help Workflows 下描述了一个功能不大的3D分析选项。59在Map View 和 Seismic View 之间使用光标追踪你所看到河道的地方。60如果你还有时间，重复流程4.2的步骤，作出10，20，30和40Hz的Polygons 图。10，20，30和40Hz的Polygons 图能显示河道随频率的变化吗？这或许是一个暗示，暗示着河道厚度的不协调。调谐图合成 (流程4.3) WF_1.0_M

44、erged_Frequencies 演示图注意：为什么频率合成在一起后，移除演示图上一些河道信息。回到 Tuning Mapper，不选Extract Slices as Horizons 并改变切片范围从Range 0 到 80 。确认其他参数的匹配性（参考流程4.1）61保持其他的参数不变并点击Compute Horizons from Slices。当完成调谐成图工作后，刷新Map View ，演示修改后的 WF_1.0_Merged_Frequencies 图。这个演示显示地震切片到分立的频率切片的强大的分解能力。通过分解后，你可以避免来自未调谐频率成分的破坏性的干扰。62RGB超览看

45、图器 (流程4.4) 超览允许你推算你的河道的RGB演示。在一张图上，你可以把RGB演示当成3个不同的估值。 如果须要打开 PowerView，Tools PowerView 从Tuning Mapper上。 通过选择 File Load Data 或者点击Data Loader 释放 Data Loader。如下所示：点击层位栏，选择WF_1.0_020Hz， WF_1.0_030Hz， WF_1.0_040Hz 层。点击“ transfer arrow” 移动所选择的层位到当前列表中。63从主窗口，选择Tools Image Composer 弹开如下的对话框。64点击RGB Blendi

46、ng 栏在下面的边框中，再点Preview 栏。你的混合图将显示在PowerView 主窗口中。演示图中20，30和40Hz分别对应红，绿和蓝。在图上显示为白色的地方，所有的频率都有一个相似的较高的地质响应；显示为单色的地方（红，绿和蓝），仅其对应的频率在区域内已调谐。这样的话， RGB混合的演示图把效果调谐到单一的图上。你可以叠加RGB演示，形成一个和分立频率能量体相似的动画。你可以通过更改最大最小值或者他们的滑动条，修改各种颜色演示的振幅范围。当你对演示图满意的时候，输入图名并点击Create Composite 按纽。65离散频率体（Discrete Frequency Volume I

47、nterpretation ）66 离散数据体解释、检测，在目的层分析之后进行。调谐图把问题限定在了一个局部的目的层段的范围上。而表征大型地震数据体时，则需要另外不同的方法。对于超出单反射数据包或目的层段的分频处理，我们推荐你使用“Discrete frequency energy cubes”（离散频率能量调谐体）。这些能量的数据体，本质来说是整个地震数据体的一个三维时窗分析，用连续的时窗。 离散频率体是由一个单一的地震数据体计算到多元的离散频率、振幅和相位的数据体。 通过一个滚动谱分析时窗,为每个模块计算出其振幅或者相位。这些谱分量会被归类成相同频率成分的模块。这种方法通常是再分析及参数化

48、目的层后进行。 67 本流程解决如下问题: 还有其他新的用来钻井的井位吗? 如何提高断层的解释精度? 流程将产生在单一频率下的振幅信息的数据体。对体的解释流程5.0 此流程与三维时窗分析一样，不同的是一个频率需要计算一个数据体。举个例子，你可以围绕产层（产油层）通过休整输入的初始数据，进行分析，并运行“Volume Recon ”仅选30Hz。68储层厚度估算（Thickness Estimation）69 频谱分析提供一个独立相位的稳定的方法来估算地震上的储层厚度。借助分频处理，即使是一个相对低的频率成分，如10Hz，都可以定量分析薄层的变化情况。 在油田上，任何定量分析的关键都在于有井控制

49、。正象地震反演技术一样，分频处理厚度分析需要井控来保证最终结果的合理性。 从块状模型薄层分析的理论来看，振幅谱的周期于频率Pf 存在V型切口的关系,而频率Pf与薄层的时间厚度t有关,如下:t = (2 * z) / Vt = 1 / Pf 其中z等于薄层厚度, V等于速度.可求得时间厚度z:z = V / (2 * Pf) 直达波的频率,f0等于1/2 Pf ,将这个关系式代入厚度等式,你会得到:z (厚度) = V / (4 * fo) t (双程旅行时) = 1 / (2 * f0) 很多时候,基于频谱分析的薄层厚度分析,可以确保薄层厚度与有效砂岩厚差在1米的范围内.而当传统的厚度分析产生歧义或者不能解决问题的时候,分频处理的厚度分析就突显了它的作用.显然地,伴随着反射系数模型复杂性和综合性的不断增加,它必将可以完善振幅谱的解释,并结合其相互的联系.这项研究正在进行中. 本流程解决如下问题: 储层的相对厚度是多少?70储层相对厚度图 (流程6.1) SpecDecomp的(R2003.0)版, 主要集中于定性分析.分析结果将显示地层是厚是薄,期望的有利井区等.但它不能使你获得精确的储集岩厚度. 再次运