Geolog-全波列声波测井中文手册-

1、Geolog 软件技术手册Full Sonic Wave Processing -SWB帕拉代姆公司北京代表处2006 年 12 月1 、综述 - 1 -1.1 预备知识 -. .1. .-1.2 数据 -. .1. .-2、阵列声波全波形 - 2 -2.1 数据准备 -. .3. .-2.1.1 查看 /创建一个声波列阵工具模版 -. .3. .- 2.1.2 练习指导 2- 创建其他波形属性 -. .5. .- 2.1.3 波形分解 -. .6. .-2.1.4 深度转换 -. .7. .-2.2 处理 -. .8. .-2.2.1 数据分析 -. .8. .-2.2.2 去噪 - 11

2、-2.2.3 设计滤波器 - 17 -2.2.4 振幅恢复 - 19 -2.3 阵列声波处理 - 20 -2.3.1 处理模块简介 - 20 -2.3.2 偶极波形处理 - 21 -2.3.3 单极波形处理 - 23 -2.3.4 拾取标志波至 - 26 -2.4 后期处理 322.4.1 综述 322.4.2 频散校正 332.4.3 传播时间叠加 362.4.4 相关性显示 382.4.5 阵列声波重处理 393 、机械性质 443.1 综述 443.2 计算动力学弹性性质 44附录 I- 快速运行 46附录 II- 频散校正讨论 471、综述欢迎阅读 Geolog软件SWB指导教程。这是

3、帕拉代姆公司 Geolog软件 SWB模块的使用教程，通过典型工作流程和程序模块的使 用，使用户熟悉阵列声波处理过程的基本工具的使用。熟悉阵列声波波形 查看/创建一个阵列声波工具模版 创建一个给定波形的其他属性 解装测井记录 预处理及波形数据分析 使用平均深度、平均时间和频率滤波器去噪 去除数据采集中振幅增益和标准化 处理偶极和单极声波测井仪下波形 自动和交互式的拾取波至频散修正传播时间覆盖查看结果的相关性图形1.1 预备知识Geolog 基本测井记录中文档，所有的输入文件为 Bold Courier New ，输出为 Bold Courier New，不加粗。1.2 数据教程中使用了下面的附

4、加文件（文件并非由软件提供）数据： geophysics_master.unl曲线模板： (copy from layouts ) swp_array_sonic_comparison swp_dispersion swp_depth_average swp_frequency03 swp_frequency05 swp_other_attributes swp_projection swp_raw_sonic_waveforms swp_reprocessing swp semblance swp traveltime函数：图像：井：polaris 03, polaris 042、阵列声波全

5、波形典型的阵列声波工具在接收剖面上有 8个接收器和 3个不同的发射器。也要注意的是，不 同服务公司出产的阵列工具也不相同。在数据采集上，基于不同的调查研究目的采用不同的工作模式。在Geolog 里，“WF”是惯用的波形记录名字后跟两位数字，首位数代表操作模式索引，第二位数表示接收器索引。例 如， WF21代表波形记录为模式 2的第一记录道（基于偶极模式）。每一个接收器的声波波形以数组记录存载。 用Geolog layout（ 模版）可以在图像道和阵列道 显示一个声波波形。在这两个道里，水平轴表示信号的旅行时，垂直轴为参考深度。当选定 一个图像道后，用户可以打开一个声波阵列查看窗口，对数据做进一

6、步分析。图1：原始声波波形显示图1为原始声波波形显示图像道（第一道和第三道），阵列道（第二道和最后一道）。最 初的两道显示的是第一接收器的偶极波形记录，最后两道是第一接收器的单极波形记录2.1 数据准备通过本节可以熟悉声波全波形数据准备 查看/创建一个列阵声波工具模版 创建给定波形的其他属性 手动解装波形 进行深度转换2.1.1 查看/ 创建一个声波列阵工具模版在波形处理前先了解工具的相关信息， 它是以文件夹形式存储。 斯伦贝谢公司的 DSI测得 的井资料在 Geolog中已经存储在 Site目录下，作为四个缺省文件（ dsi_m1 to dsi_m4）。当有 改动或者初次运行此工具，相关信息

7、会保存在本工具文件夹下的 specs路径中。工具模版的数据信息可以在 Text中的 Constants 表格中看到1、启动 Geolog软件打开 STARS工区2、点开 Well 3、点击 well>open ，打开 Polaris_03 4 、选择 Petrophysics > Full Waveform Sonic > Data Preparation > Create Tool Specification模块显示斯伦贝谢公司的 spec文件dsi_m2 ，给定的为缺省值，如下图 2图 2： array_sonic_create_toolspec中 dsi_m2 的

8、缺省值Spec文件中的参数如下表所示：参数描述TOOL NAME声波采集仪器的名称如： DSI, XMACTOOL SPEC声波工具模版，以模版文件命名TOOL MODE声波工具的工作模式，通常在 Text> Constants 表格中定义TOOL R声波仪器的直径SAMPLE RATE声波采样率，通常以微次生法测量TR DISTANCE发射器与第一接收器的距离X INTERV两个邻近接收器之间的距离DEPTH_SHIFT工具深度的平差不同于工具的参考深度和井的参考深度。负的值表示工具参考深度在井 的参考深度之下5、如果需要的话可以对缺省值做调整6、点击 Start 运行模块 一个新的模

9、块文件在本工区的 SPECS目录下产生了7、选择 Launcher > Close 来关闭阵列声波工具模块8 、检查新模块 dsi_m2 的值 ，如图 3 所示图3： dsi_m2 的存储位置以及值提示：不管Geolog 中使用的量度制，米制的或是特定的，声波阵列工具模版总是会保存为米制的。当要 运行这个工具，首先请确认工具模版的值，以确保得到正确的结果 。2.1.2 练习指导 2- 创建其他波形属性声波波形可以用反射振幅正常显示。频率谱、瞬时相位、和瞬时振幅属性可以从解释角 度很好的量化突出这些特性。原始声波频率谱显示的是不同到达时的主要频率围，这样可以 在频率域查看各到达模式的属性特

10、征。瞬时相使弱信号加强，但是也同样使噪音增强。瞬时 振幅是包络的振幅，包含更多侧向变化，但是降低了垂向分辨率。 tp_array_sonic_attribute 模 块用于计算一个给定波形的这些属性。 三个测井记录可以从这个模块输出 - 频率谱、瞬时相位 以及瞬时振幅，可以在图像道或者阵列道中查看这些属性。1、选择 Petrophysics > Full Waveform Sonic > Data Preparation > CreateOther Attributes 来显示 tp_array_sonic_attribute 模块。2、如图 4所示，对窗口做如下值改变Inp

11、ut Set改为SONICOutput Set改为SONIC_ADJLOG_IN类型为 WF21FREQ_OUT类型为 WF21FREQINST_PHASE类型为 WF21PHASEINST_AMP类型为 WF21AMP图 4: 改变属性模型参数3、点击 Start 运行4、打开曲线模板 swp_other_attributes.layout 来查看结果 ，如图 5所示。原始数据 原始数据频率谱 原始数据瞬时相位 原始数据瞬时振幅 图 5 ：波形属性显示2.1.3 波形分解来自每个接收器的声波波形有时被叠加记录到一起。这个模块用于把波形分解为每个接 收器的单个记录。通常，每个模块有自有的测深记

12、录，这个记录要作为输出的参考深度。在 这个模块里，需要提供分解波形的名字前缀。要遵守 Geolog 的命名约定。为了做示，本例中 的未分解波形都刻意不进行自动分解。1、选择 Petrophysics > Full Waveform Sonic > Data Preparation > LogUnpacking 来显示 tp_array_sonic_unpack 模块2、如图 4所示，对窗口做如下值改变 Input Set改为SONIC_PACKED PREFIX类型为WG DEPTH选择 PWD2 LOG_WF选择 PWG2图 6: 改变分解波形模型参数3、点击 Start

13、运行4、将PREFIX值变为 WN ，LOG_WF值变为 PWN2 ，再运行模型5、打开 Text ，检查未分解记录值。如下图 7所示图 7 ：分解波形输出记录2.1.4 深度转换DSI工具数据的参考深度与井记录的参考深度不同， 差值从测量工具的底部到接收器剖面 的中段进行修正，如上节的 depth-shift 。1、选择曲线道模版 swp_other_attributes2、选择 Tools > Depth Shift > Apply Logs3、如图 8所示，对窗口值做如下改变 Input Set 改为 SONICOutput Set改为SONIC_DEPTH_SHIFTDEP

14、TH_OFFSET输入参数为 -7.9248LOG_IN选择波形 WF21_1为WF28_1 WF41_1为WF48_1CORRECTED选择DEPTH图 8： 深度转换模型4、点击 Start 运行5、关闭所有窗口2.2 处理通过本节可以使用户熟悉全波形声波处理流程 在声波数组查看窗口操作和解释数据 使用平均时间和平均深度以及 Butterworth 频率滤波器去除波形噪音 设计一个滤波器去除数据采集中的振幅增益和标准化2.2.1 数据分析 为了保证处理结果质量，在数据处理前要检查数据确保无误。波形数据通常包含噪声。通 常需要选择适当的工具来尽可能多的去噪。1、打开模板 swp_raw_so

15、nic_waveforms.layout 。2、选择 WF21 的图像道，选择 View > Array Sonic View 。 数组声波查看窗口显示波形 WF21 以及它当前深度频率域的位置。（见图 9）要打开阵列声波查看窗口，需要选择打开带数据的图像道阵列声波查看窗口图 9 ：图像道和声波数组窗口阵列声波窗口用来详细的分析波形，它有两个显示区域。显示区在当前深度位置显示信 号波形，可以估算出时间间隔（接收器的时间延时）。频率域显示区显示每个波形相应的的 频率围。3、点击声波数组记录 - Log Select 图标打开可选择多条曲线。 提示：挑选出来的记录必须来自同一个集合和模型工具

16、（例如：有相同的模型索引- 字符名字后的第一位数相同）在图像道上显示。4、选择记录 SONIC.WF21_1 to SONIC.WF28_1 ，点击 OK。 多重曲线以及其对应的频率域在数组声波的窗口显示。5、在模版上第一图像道 （SONIC.WF21） 点击任何位置。声波数组显示窗口的信号在每次 选择深度点的时候会刷新。另一个方式是通过深度区域，键入围值或者使用上下按键。6、设置深度点和时间围如下所示：在显示区放大或者改变时间围7、将光标定位在信号显示区任意位置， 按住鼠标左键拖曳定义一个矩形区域，放开鼠标 左键。通过显示的时间围值来反映被选择的区域。将时间围值的起始和终点值输入与用鼠标改变

17、信号显示区域的效果是 一样的。8、再次点击鼠标左键使区域变为原来大小。 在波形上估算时间间隔： 提示： 这个功能只在多重波形在当前活动的时候起作用。当时间间隔固定时，功能不可选时首先将信号显示区放大9、固定时间间隔，在慢度剖面显示。10、显示区，在第一个接收器 （WF21）的信号，按住鼠标左键，拖至到最后一个接收器 的信号处 （WF28 ），放开鼠标左键。如图 10淡红色显示 估算时间间 隔，波至慢度 及接收间隔粗淡红色的 陡线显示第 一 个和 最后 一个接收器 初至时间间图 10 ：时间间隔估算 8个接收器的信号以及频率域在当前选择的深度显示。11、不选 Moveout 复选栏，不显示时间间

18、隔数据。12、关闭阵列声波显示窗口。2.2.2 去噪 波形图总是会包含噪声，太多的噪声会不能识别可检测的波至。预处理是提高数据质量 保证初至被识别的重要步骤。Geolog 有一系列基本的预处理工具可以满足一个波形的任意去噪要求。 去除时间均值，在时间方向使用平均时滤波器 去除深度均值，在深度方向使用平均深度滤波器Butterworth 滤波器，是一个频率滤波器 根据具体的需求，用户可决定针对某个波形使用对应的滤波器 （1）去除时间均值噪声 这个模块提供沿时间方向去除噪声滤波器。先在每个窗口计算出平均值，然后从原始波形移除噪声。) 这是用于移除低频噪声背景的滤波器 。1 、选择 Petrophy

19、sics > Full Waveform Sonic > Pre-Processing > RemoveTime-Average2、如图 10，对打开窗口中值做如下修改。改为SONIC_DEPTH_SHIFT Input Set UTIME WINDOW LOG_IN例如，设置滤波窗口值为 520 in micro-seconds (usecs)忽略采样场。这个模块以取得合适的结果。设为WF21 (可以同时导入多个文件)图11 ：时间平均滤波器参数设置是在Geolog不能识别采样率时需要采样场)改变窗口的取值，多次运行3、点击 Start 运行模块4、打开绘图模版 swp_a

20、rray_sonic_comparison.layout 。图13 ：设置平均深度滤波器模型值图12：用时间平均滤波器滤波前后对比图（2）去除深度均值噪声 这个模块提供沿深度方向去除噪声滤波器。先在每个活动窗口计算出平均值，然后从原 始波形移除噪声。这个滤波器可去除典型的诸如套管信号或者仪器带来的噪声。1、打开井 Polaris_04 ，保存对先前井所作的改变。 2、关闭所有打开窗口。3、选择 Petrophysics>Full Waveform Sonic>Pre-Processing> RemoveDepth-Average.4、在 tp_array_sonic_dept

21、h_average 窗口对值做如下改变（如图 13 ） Input Set改为SONIC WINDOW设置滤波器长度为 21 LOG_INWF415、点击 Start 运行模型6、打开绘图模版 swp_depth_average.layout 。套管和仪器形成的信号可以很强，使得重要的波至难以识别。图 14 左边道显示的波形是一个这样的例子，直直的垂线占主要部分。右边的道显示的是不需要信息被移除后的清晰的 波至。运行滤波器后的结果图图14：用深度滤波器滤波前后对比图3）使用频率滤波器Butterworth 滤波器是一个频率滤波器， 可以通过给定一个高的和低的截止值以及创建滤 过顺序来实现。 B

22、utterworth 滤波器选择的顺序 （从 1到8）决定了从通频带到抑制带的过渡陡度。单位设为千赫兹而不是赫兹，是为了便于进行声波波形的处理。如果数据是来自于标准 的DSI工具采集，可以忽略采样率，因为 Geolog 可以自动地读取。有关滤波器的调整适应更多信息参考菜单选项下的 Design Filter 。1、保存对 Polaris_04的改动，打开井 Polaris_03 。2、选择 Petrophysics > Full Waveform Sonic > Pre-Processing > Frequency Filter3、对打开窗口做如下改变（如图 15所示） In

23、put SetSONIC_DEPTH_SHIFT Output Set HIGH_END改为FILTERED03 SAMPLE_RATE 10 LOG_IN从输入集合选择8个波形： SONIC_DEPTH_SHIFT.WF21_1 to设为 3HZ（移除所有高于这个值之外的频率信号 ）WF28_1图 15 ：改变 butterworth 滤波器的频率围和大小4、点击 Start 运行这个模型5、打开 swp_frequency03.layout6、选择 View > Array Sonic View 来查看滤波后图形7、设置当前深度到如图 16所示。在 3HZ以上的频率已经在当前深度位置

24、被移除图16：WF滤波后图形8、关闭阵列声波查看窗口。9、选择 tp_array_sonic_filter 模块，重复第二步和第三步改变频率围为 0-5 赫兹，输出名 为FILTERED05。10、点击 Start 运行模块11、打开 swp_frequency05.layout12、打开声波数组查看窗口。大于5Hz 的频率信号都被移除，如图 17所示，在当前深度位置显示的频率谱。图17：滤除高于 5Hz波后的 WF21 波形图13、关闭阵列声波窗口2.2.3 设计滤波器Butterworth 滤波器将作为频率滤波器用在下面的例子中。 Butterworth 是一个很好的带通 滤波器， 在损失

25、了部分频率域后有较为合理的时间衰减。 高频和低频截止值用于带通和带止。 对于低通滤波器， 高频波被截止， 反之亦然。 Butterworth 滤波器的排列控制着在通频带之外 信号变弱的快慢。设计一个合适的滤波器或是子波，必须要注意采样率和频率值。设计一个Butterworth 子波，采样率( SR)要使用 0.001秒，这样可以得到 Nyquist frequency = 1 / (2 * SR) = 500 Hz 。 输入频率必须要小于尼奎斯特频率( Nyquist frequency )频率。也就是说，如果提供的值大 于尼奎斯特频率( Nyquist frequency )，采样率也要更好

26、。1、选择 Petrophysics > Full Waveform Sonic > Pre-Processing > DesignFilter2、在 filter_create 窗口对值做如下改变（如图 18 ）FILTER_TYPE改为ButterworthFILTER_NAMEbw1040LOW_START1000HIGH_END3000PLOT_SR0.00004DISP_SCALEAMPLITUDE图18：滤波器创建模块3、点击 Start 运行这模块。4、查看创建的滤波器（如图 19）。 Input Set NORM_OPT GAIN_OPT REFERENCE选

27、择UNPACK.DEPTH图19 ：运行创建滤波器模块后的结果2.2.4 振幅恢复在数据采集中，振幅增益以及标准化是为了增强波形信号的质量。可基于有效的增益和标准化因素（归一化因素），使波形回到原来形状。默认的，去除标准化因素选择inverse ；去除增益选择为 multiply 。1、选择Petrophysics > Full Waveform Sonic > Pre-Processing > Amplitude Recovery.2、在 tp_array_sonic_tar 窗口对值做改变（如图 20所示）改为SONIC_DEPTH_SHIFT设为INVERSE设为MUL

28、TIPLY LOG_WF选择 WF21_1至WF28_1波形 LOG_NORM 选择井 UNPACK.WN21_1 到 UNPACK.WN28_1 LOG_GAIN选择井 UNPACK.WG21_1 到UNPACK.WG28_1图20 ：振幅恢复模块3、点击 Start 运行模型。4、打开text 窗口在原来的集合中查看井数据(如图 21所示)。图21 ：Text 窗口查看去除振幅增益后的波形记录2.3 阵列声波处理通过本节，用户可以：处理偶极波形处理单极波形使用自动交互式拾取波至2.3.1 处理模块简介声波波形处理基于慢度时间相关法。这个模块在活动窗口过时间和慢度方向可以查看多个声波。默认的

29、，时间围可以包括所有的记录时间；慢度的围至少包括目的区。例如，慢度围为100800 us/m (1200-10000 m/s) 可以包含大多数区域的纵波波至。建议用较大的慢度围这样可以尽可能包含所有波至，但这也带来更多的计算量首先选择少量数据来决定时间和慢度方向的最优选择要注意慢度步长的选择，过大的步长会降低最终声波分辨率。在每个深度点计算出的最大值作为与深度相关的函数映射到慢度轴上；这是最终输出Projection log 。这个模块有很多控制参数，详细参数解释可以查看模块帮助。2.3.2 偶极波形处理在本节，被处理的波形filtered03 和 filtered05 将用于下文的离散。1、

30、保存井文件 Polaris_03。2、关闭所有查看窗口。3、选择 Petrophysics> Full Waveform Sonic > Processing > Waveform Processing来打开 tp_array_sonic_process 窗口。4、设置输入参数（如图 22所示） Input SetSONIC_DEPTH_SHIFT Output SetDIPOLE TOOL_SPECdsi_m2 START0 STOP10000SLOWNESS_FROM 35 SLOWNESS_TO500 SLOWNESS_STEP4 LOG_IN由SONIC_DEPTH_

31、SHIFT集合，选择8个波形： WF21_1 to WF28_1图22：偶极波形处理模型窗口5、点击 Start 运行模型。当模型运行结束，结果以相关性图件显示，如图 23所示。数据在当前深度位置显示（当深度未规定时是在井记录的顶部）图23 ：相关性图 - WF21 - WF28 波形处理结果图6、保存相关性图为 semblance_plot_01 。将相关性图以唯一文件名保存，便于以后查看，否则图件会被覆盖7 、改变输入和输出集合的名字为FILTERED03，重复 4-6 步，将相关性图形保存为semblance_plot_03 。8、改变输入和输出集合的名字为FILTERED05，重复4-

32、5 步，将相关性图形为semblance_plot_05 。9、关闭相关性图。10、保存井文件。2.3.3 单极波形处理11、在 tp_array_sonic_process窗口，设置输入参数，如图 24所示： Input SetSONIC_DEPTH_SHIFT Output SetMONOPOLE TOOL_SPECdsi_m4 START0 STOP5120 UTIME_WINDOW400 UTIME_STEP50 SLOWNESS_FROM30 SLOWNESS_TO298.7 SLOWNESS_STEP1.98 LOG_IN从SONIC_DEPTH_SHIFT集合，选择 8 个波形

33、WF41_1to WF48_1图24：单极波形处理模型窗口2D 相关12、点击 Start 运行模块。当模型运行结束输出一个映射和一个相关性图件，图件在每个深度点存储的是 性图像， 可以通过阵列声波数组查看窗口进行质量控制。如下图所示图25 ：阵列声波数组窗口查看图件13、关闭相关性图件，不保存改变。14、打开模版 swp_projection.layout 。上面例子中的偶极和单极子映射记录如图 26所示偶极波形处理后图形 单极子波形处理后图形 水平轴显示慢度 (START)到 (END) 图26：深度相关函数映射15、在WELL下点击 SAVE.2.3.4 拾取标志波至映射记录在慢度 -

34、深度平面代表了一个或多个相关性。它们与某些波至有较好的一致性。 模型在输入的映射记录中自动或者交互式的寻找相关最大值位置，从慢度投影记录拾取相应 一组波至。搜索从一个起始值开始(初始慢度值在映射记录里)来标定某一初至，搜寻围以 从映射记录估计的起始值为中心。估算输入起始值最简便方法是将投影记录插入到图像道， 打开阵列生波可视窗口， 移动鼠标标定出信号显示区的目的位置， 光标位置 X的值会在位置对 话框的状态栏显示出来。这个值就是从图上读取的慢度值，它可以用来作为搜索起始值。对于交互式拾取，测井曲线首先是自动被拾取(参见 Curve Insert 联机帮助)，编辑，最 后作为备用初始记录被保存。

35、主要有：1、自动拾取一个偶极测井记录2、自动拾取一个单极测井记录3、对于第二步中不合适点做交互式拾取自动拾取一个偶极测井记录1 、选择 Petrophysics > Full Waveform Sonic > Processing > Picking 来打开 the tp_array_sonic_label 窗口。2、选取输入起始值 (SLOWNESS_S1)。让swp_projection 为活动窗口，选择偶极投影道来打开阵列声波可视窗口。在信号显示区移动鼠标到拾取位置， 读取显示出的 X值。这个值即是从图上读取的慢度 值，可作为起始值 (SLOWNESS_S1)，如例子中

36、图 27示。关闭阵列声波可视窗口。图27 ：确定搜索起始值3、选择模型窗口，设置值如下图 28所示 Input/Output Set选择 DIPOLE，这是处理后的输出集合图29：自动拾取的起始记录覆盖在偶极映射图像上 SLOWNESS_WINDOW 本例中设为 121 us/f （因为仅有一个波至，对于 偶极使用一个大点的围是可靠的）窗口长度必须包括波至的最大波动幅 度，但是也不能过大要注意排除邻区波至。 SLOWNESS_S1 在前一步操作中， 位置对话框里确定的 X 值，必须要接近值121图28偶极自动拾取模型设置4、运行模型 5、重复地3-4 步二次，改变输入输出集合名字分别为FILT

37、ERED03，FILTERED05。6、设定 swp_projection 为活动窗口，查看曲线道（如图 29示） PICK_OPTIONBOTH7、打开属性对话框，注意 DIPOLE.DT_PICK1井的左右边界为处理的慢度起 始和终止值(35 and 500 US/F )，输入模型中。8、关闭对话框。单极投影记录自动拾取9、确定输入 tp_array_sonic_label 模型所需要的的起始值 (SLOWNESS_S1 and SLOWNESS_S2)选择单极图像道打开阵列声波可视窗口。在信号显示区使用鼠标定位， 在定位对话框确定两个初始值， 如图 30所示 为确定第二波至的例子。图30

38、 ：确定 SLOWNESS_S2的输入初始值10、使 tp_array_sonic_label 为活动窗口，对模型值做如下设置(如图 31所示)： Input / Output Set选择MONOPOLE集合，它是处理中的输出集合 Slowness_Window设为15.24 us/f 为搜索窗口大小 SLOWNESS_FROM 30 SLOWNESS_TO 298.7 SLOWNESS_STEP 1.98 SLOWNESS_S1 60 SLOWNESS_S2 85图31 ：交互式拾取两个井记录模型界面11、运行模型 12、选择模板 swp_projection 作为查看模板，如图 32所示。

39、图32：自动拾取井记录覆盖在单极图像道上13、菜单栏 WELL下点击 SAVEInteractive Picking 交互式拾取14 、 用 swp_projection 模 板 ， 在 单 极 投 影 道 图 像 道 上 选 择 拾 取 记 录DT_PICK2。15、选择 Tools > Curve Insert 。16、手动的重新拾取不合适的横波波至。17、打开文本查看重命名被编辑的井记录 DT_PICK2 version 2 为DT_SEED。18、菜单栏 WELL下点击 SAVE。19、选择 tp_array_sonic_label 模型，对值做如下改变（如图 33）： PICK

40、 OPTION SINGLE SLOWNESS_S1使用DT_SEED作为初始记录 DT_PICK1使用DTS_PICK图33：单井交互式拾取模型界面20、运行模块21、在曲线模板中查看结果，使用图 34作为参考。 如果效果不够好， 调整搜 索围（例如： SLOWNESS_WINDOW 和SLOWNESS_S1值），重复先前的步骤。图34：偶极和单极的投影记录的自动式和交互式拾取显示在单极投影区，绿色的曲线表示作为初始记录的自动拾取曲线 (DT_PICK2) 蓝色曲线是基于初始子记录的交互式拾取 (DTS PICK)。提示：处理的慢度记录需要进行深度转换，这涉及到井的参考深度。在Geolog

41、中不会自动转换，因为对于参考深度工具，没有足够的信息可供选择。2.4 后期处理通过本节，用户将熟悉后期处理流程 计算并应用频散校正，确保每个井记录在同一参考系下 应用频散校正导出旅行时查看相关性图形显示 阵列声波数组重处理，计算平均慢度2.4.1 综述后期处理以及后期质量控制是声波处理的重要部分。 从处理过的记录中导出 的旅行时可以作为质量控制来描述声波波形波至。 频散以及井眼情况等将影响最终处理结果频散校正模型计算出频散曲线，在每个深度位置使用校正值来修正频散影 响。波形将按照一定间隔被重处理 （在后期处理中） 来增强垂直分辨率或者补偿 井眼较大的直径差值。2.4.2 频散校正这个模块输出修

42、正后的慢度记录。 它可以通过覆盖阵列道随意输出频散和频 谱记录，也可以输出坏道来标记有问题区域。1、菜单栏 WELL下点击 SAVE。2、选 择Petrophysics > Full Waveform Sonic > Post-Processing >Dispersion Correction 。提示：离散修正是一个很费时的处理过程，本教程中为了实现这个操作只是用了少量 资料来处理。3、如图 35所示，设置处理参数 。图35：设置处理围4、输入参数如下：如图 36所示 Input SetSONIC_DEPTH_SHIFT Output SetRAWTOOL_SPECdsi_m

43、2START0STOP10000CALIWIRE.HCALBSWIRE.BSRHOB_F1.24 G/C3(1240 K/M3)RHOB_SWIRE.RHOZDT_F245.8 US/F (75 US/M )DT_PMONOPOLE.DT_PICK1DT_SDIPOLE.DT_PICK1LOG_IN从 SONIC_DEPTH_SHIFT集合, 选择8个波形WF21_1 to WF28_ 1图36：频散校正窗口5、运行模块处理速度取决于计算机的性能和数据量的大小6、重复第三和第四步，改变值 Input /Output setFILTERED03 DT_S7、重复第三和第四步，改变值FILTERE

44、D03.DT_PICK1 Input /Output setFILTERED05FILTERED05.DT_PICK1 DT_S8、打开模板 swp_dispersion.layout 。提示： 离散修正结果如图 37 所示，处理围为 13150 - 13450 ft.图37 ：离散结果对比实例（数据为 13150 - 13450 ft ）1,2,3道分别表示覆盖在叠加频率域的波形离散曲线， 4道表示低频滤波后的 波形， 5道波形大于5 Khz的频率被滤过， 6道大于3 Khz频率被滤过。对原始波形做处理时， 发现在做频散校正后， 有不稳定的波动现象。 如上例， 围如箭头所示，通过0-5 KH

45、z带通滤波处理，修正的剪切波慢度稳定， 同时1-3 KHz 滤波器给出的离散都很清晰。9、菜单栏 WELL下点击 SAVE。2.4.3 传播时间叠加在前面步骤中得到的慢度记录，可以从中导出波形的旅行时。导出的旅行 时可以作为质量控制参数， 稍后可以将这个旅行时覆盖在波形上来形象化的检查 处理结果的质量。导出旅行时的模块需要处理出来的慢度记录以及钻孔直径记录图， 这些资料 都必须与参考深度匹配。1、关闭所有打开的窗口，保存所做的改动。2、选择 Petrophysics > Full Waveform Sonic > Post-Processing > Traveltime打开旅

46、行时模块3、如图 37所示，输入以下参数（如图 38所示） Input/Output SetDIPOLE5、打开曲线模板 swp_traveltime.layout DT CALIDT_PICK1_1203WIRE.HCAL_1 DT_F图38 ：建立旅行时模块4、运行模块图39：第一个接收器计算的旅行时(从处理得到的慢度中导出)2.4.4 相关性显示相关性图形是由慢度时间相关分析( STC)得到，如果 STC处理有差错，可 通过这个图件做检测分析。1、打开 swp_semblance.layout .2、打开阵列声波窗口改变下列值：从FILTERED03集合增加声波阵列记录 (WF21 -

47、WF28)在显示的 Display Waveform(s) Plus 里，选择 Semblance ，并且选择记录 为FILTERED03.SEMBLANCE。选择the FILTERED03.SEMBLANCE log结果如图 40所示。图40 ：在阵列声波数组中显示的相关性图3、菜单栏 WELL下点击 SAVE。4、关闭所有打开的窗口2.4.5 阵列声波重处理TRANSM波形将按照一定间隔被重处理 (在后期处理中) 来增强垂直分辨率或者补偿 井眼的直径差值。 重处理是基于 multiple-shot processing 技术，它包含了由于 与工具位置点重合而带来的冗余信息。 选择若干子阵列， 这样子阵列的大小可以 覆盖住调查区的地层厚度。有两个方法 (configurations) 来进行处理， RECEIV和TRANSM ，RECEIV使用的是共源极配置，可以用来提高垂向分辨率；使用的是共接收配置。 由这两个配置得来的平均慢度用于补偿因井眼大小改变而 带来的影响。可使用模块 Calculate Mean 完成。提示： 当信号波长非常长时，这个方法是无效的