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

1、Geolog 软件技术手册 Full Sonic Wave Processing -SWB帕拉代姆公司北京代表处帕拉代姆公司北京代表处20062006 年年 1212 月月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.3

2、 设计滤波器.- 17 -2.2.4 振幅恢复.- 19 -2.3 阵列声波处理.- 20 -2.3.1 处理模块简介.- 21 -2.3.2 偶极波形处理.- 21 -2.3.3 单极波形处理.- 23 -2.3.4 拾取标志波至.- 26 -2.4 后期处理.342.4.1 综述 .342.4.2 频散校正 .352.4.3 传播时间叠加.372.4.4 相关性显示 .382.4.5 阵列声波重处理.403、机械性质、机械性质.443.1 综述.443.2 计算动力学弹性性质.44附录附录 I-快速运行快速运行.45附录附录 II-频散校正讨论频散校正讨论.461、综述、综述欢迎阅读Geo

3、log软件SWB指导教程。这是帕拉代姆公司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- 2 -图像：函数：井：polaris_03, polaris_0

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

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

7、i_m4）。当有改动或者初次运行此工具，相关信息会保存在本工具文件夹下的specs路径中。) 工具模版的数据信息可以在Text中的Constants 表格中看到。1、启动Geolog软件打开STARS工区2、点开 Well 3、点击wellopen，打开Polaris_034、选择Petrophysics Full Waveform Sonic Data Preparation Create Tool Specification.模块显示斯伦贝谢公司的spec文件dsi_m2，给定的为缺省值，如下图2- 4 -图2：array_sonic_create_toolspec 中中dsi_m2的缺省

8、值Spec文件中的参数如下表所示：参数描述TOOL_NAME声波采集仪器的名称如：DSI, XMACTOOL_SPEC声波工具模版，以模版文件命名TOOL_MODE声波工具的工作模式，通常在TextConstants表格中定义TOOL_R声波仪器的直径SAMPLE_RATE声波采样率，通常以微次生法测量TR_DISTANCE发射器与第一接收器的距离X_INTERV两个邻近接收器之间的距离DEPTH_SHIFT工具深度的平差不同于工具的参考深度和井的参考深度。负的值表示工具参考深度在井的参考深度之下5、如果需要的话可以对缺省值做调整6、点击Start运行模块一个新的模块文件在本工区的SPECS目

9、录下产生了。7、选择Launcher Close来关闭阵列声波工具模块8、检查新模块dsi_m2的值，如图3所示图图3：dsi_m2的存储位置以及值的存储位置以及值- 5 -提示：不管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所示，对窗口做如下值改变 Input Set改为SONIC Output Set改为SONIC_ADJ

11、 LOG_IN类型为WF21 FREQ_OUT类型为WF21FREQ INST_PHASE类型为WF21PHASE INST_AMP类型为WF21AMP图图4: 改改变变属属性性模模型型参参数数3、点击Start运行4、打开曲线模板swp_other_attributes.layout来查看结果，如图5所示。- 6 -原始数据原始数据 原始数据频率谱原始数据频率谱 原始数据瞬时相位原始数据瞬时相位 原始数据瞬时振幅原始数据瞬时振幅图图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- 7 -图图6: 改改变变分分解解波波形形模模型型参参数数3、点击Start运行4、

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

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

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

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

17、值输入与用鼠标改变信号显示区域的效果是一样的。8、再次点击鼠标左键使区域变为原来大小。在波形上估算时间间隔：提示：这个功能只在多重波形在当前活动的时候起作用。)当时间间隔固定时，功能不可选时首先将信号显示区放大。9、固定时间间隔，在慢度剖面显示。10、显示区，在第一个接收器（WF21WF21）的信号，按住鼠标左键，拖至到最后一个接收器的信号处（WF28WF28），放开鼠标左键。如图10- 11 -图图10：时时间间间间隔隔估估算算8个接收器的信号以及频率域在当前选择的深度显示。11、不选Moveout复选栏，不显示时间间隔数据。12、关闭阵列声波显示窗口。2.2.2 去噪去噪波形图总是会包含噪

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

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

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

21、rray_sonic_depth_average窗口对值做如下改变（如图13） Input Set改为SONIC WINDOW设置滤波器长度为21 LOG_INWF41图图13：设设置置平平均均深深度度滤滤波波器器模模型型值值5、点击Start运行模型6、打开绘图模版swp_depth_average.layout。套管和仪器形成的信号可以很强，使得重要的波至难以识别。图14左边道显示的波形是一个这样的例子，直直的垂线占主要部分。右边的道显示的是不需要信息被移除后的清晰的波至。- 14 -图图14：用深度滤波器滤波前后对比图：用深度滤波器滤波前后对比图（3）使用频率滤波器Butterworth

22、 滤波器是一个频率滤波器，可以通过给定一个高的和低的截止值以及创建滤过顺序来实现。Butterworth 滤波器选择的顺序（从1到8）决定了从通频带到抑制带的过渡陡度。单位设为千赫兹而不是赫兹，是为了便于进行声波波形的处理。如果数据是来自于标准的DSI工具采集，可以忽略采样率，因为Geolog可以自动地读取。有关滤波器的调整适应更多信息参考菜单选项下的Design Filter。1、保存对Polaris_04的改动，打开井Polaris_03。2、选择Petrophysics Full Waveform Sonic Pre-Processing Frequency 有低频噪声的原始图像运行滤波

23、器后的结果图- 15 -Filter.3、对打开窗口做如下改变（如图15所示） Input Set SONIC_DEPTH_SHIFT Output Set 改为FILTERED03 HIGH_END 设为3HZ(移除所有高于这个值之外的频率信号) SAMPLE_RATE10 LOG_IN 从输入集合选择8个波形：SONIC_DEPTH_SHIFT.WF21_1 to WF28_1图图15 ：改改变变butterworth滤波器的频率范围和大小滤波器的频率范围和大小4、点击Start运行这个模型5、打开swp_frequency03.layout6、选择View Array Sonic Vie

24、w来查看滤波后图形7、设置当前深度到如图16所示。在3HZ以上的频率已经在当前深度位置被移除。- 16 -图图16：WF滤波后图形滤波后图形8、关闭阵列声波查看窗口。9、选择tp_array_sonic_filter模块，重复第二步和第三步改变频率范围为0-5赫兹，输出名为FILTERED05。10、点击Start运行模块11、打开swp_frequency05.layout12、打开声波数组查看窗口。大于5Hz的频率信号都被移除，如图17所示，在当前深度位置显示的频率谱。- 17 -图图17：滤除高于：滤除高于5Hz波后的波后的WF21波形图波形图13、关闭阵列声波窗口。2.2.3 设计滤波

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

26、。也就是说，如果提供- 18 -的值大于尼奎斯特频率（Nyquist frequency），采样率也要更好。1、选择Petrophysics Full Waveform Sonic Pre-Processing DesignFilter.2、在filter_create窗口对值做如下改变（如图18） FILTER_TYPE 改为Butterworth FILTER_NAME bw1040 LOW_START 1000 HIGH_END 3000 PLOT_SR 0.00004 DISP_SCALE AMPLITUDE图图18：滤波器创建模块：滤波器创建模块3、点击Start运行这模块。4、查看

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

28、 NORM_OPT 设为INVERSE GAIN_OPT 设为MULTIPLY- 20 - REFERENCE 选择UNPACK.DEPTH 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 阵列声波处理阵列声波处理通过

29、本节，用户可以：处理偶极波形处理单极波形- 21 -使用自动交互式拾取波至2.3.1 处理模块简介处理模块简介声波波形处理基于慢度时间相关法。这个模块在活动窗口中通过时间和慢度方向可以查看多个声波。默认的，时间范围可以包括所有的记录时间；慢度的范围至少包括目的区。例如，慢度范围为100800 us/m (1200-10000 m/s)可以包含大多数区域的纵波波至。建议用较大的慢度范围这样可以尽可能包含所有波至，但这也带来更多的计算量。) 首先选择少量数据来决定时间和慢度方向的最优选择。要注意慢度步长的选择，过大的步长会降低最终声波分辨率。在每个深度点计算出的最大值作为与深度相关的函数映射到慢度

30、轴上；这是最终输出Projection log。这个模块有很多控制参数，详细参数解释可以查看模块帮助。2.3.2 偶极波形处理偶极波形处理在本节，被处理的波形filtered03 和filtered05将用于下文的离散。1、保存井文件Polaris_03。2、关闭所有查看窗口。3、选择Petrophysics Full Waveform Sonic Processing Waveform Processing来打开tp_array_sonic_process窗口。4、设置输入参数（如图22所示） Input Set SONIC_DEPTH_SHIFT Output Set DIPOLE TOO

31、L_SPEC dsi_m2 START 0 STOP 10000 SLOWNESS_FROM 35 SLOWNESS_TO 500 SLOWNESS_STEP 4 LOG_IN 由SONIC_DEPTH_SHIFT集合，选择8个波形： WF21_1 to WF28_1- 22 -图图22：偶极波形处理模型窗口：偶极波形处理模型窗口5、点击Start运行模型。当模型运行结束，结果以相关性图件显示，如图23所示。数据在当前深度位置显示（当深度未规定时是在井记录的顶部）- 23 -图图23：相相关关性性图图- WF21 - WF28波波形形处处理理结结果果图图6、保存相关性图为semblance_p

32、lot_01。)将相关性图以唯一文件名保存，便于以后查看，否则图件会被覆盖。7、改变输入和输出集合的名字为FILTERED03，重复4-6步，将相关性图形保存为semblance_plot_03。8、改变输入和输出集合的名字为FILTERED05，重复4-5步，将相关性图形为semblance_plot_05。9、关闭相关性图。10、保存井文件。2.3.3 单极波形处理单极波形处理11、在tp_array_sonic_process窗口，设置输入参数，如图24所示： Input Set SONIC_DEPTH_SHIFT Output Set MONOPOLE- 24 - TOOL_SPEC

33、dsi_m4 START 0 STOP 5120 UTIME_WINDOW 400 UTIME_STEP 50 SLOWNESS_FROM 30 SLOWNESS_TO 298.7 SLOWNESS_STEP 1.98 LOG_IN 从SONIC_DEPTH_SHIFT集合，选择8个波形WF41_1to WF48_1图图24：单极波形处理模型窗口单极波形处理模型窗口12、点击Start运行模块。当模型运行结束输出一个映射和一个相关性图件，图件在每个深度点存储的是2D 相关性图像， 可以通过阵列声波数组查看窗口进行质量控制。如下图所示。- 25 -图图25：阵列声波数组窗口查看图件：阵列声波数组

34、窗口查看图件13、关闭相关性图件，不保存改变。14、打开模版swp_projection.layout。上面例子中的偶极和单极子映射记录如图26所示。- 26 -偶极波形处理后图形 单极子波形处理后图形 水平轴显示慢度(START)到(END)图图26：深度相关函数映射：深度相关函数映射15、在WELL下点击SAVE.2.3.4 拾取标志波至拾取标志波至映射记录在慢度-深度平面内代表了一个或多个相关性。它们与某些波至有较好的一致性。模型在输入的映射记录中自动或者交互式的寻找相关最大值位置，从慢度投影记录拾取相应一组波至。搜索从一个起始值开始（初始慢度值在映射记录里）来标定某一初至，搜寻范围以从

35、映射记录估计的起始值为中心。估算输入起始值最简便方法是将投影记录插入到图像道，打开阵列生波可视窗口，移动鼠标标定出信号显示区的目的位置，光标位置X的值会在位置对话框的状态栏显示出来。这个值就是从图上读取的慢度值，它可以用来作为搜索起始值。对于交互式拾取，测井曲线首先是自动被拾取（参见Curve Insert联机帮助），编辑，最后作为备用初始记录被保存。主要有：1、自动拾取一个偶极测井记录2、自动拾取一个单极测井记录- 27 -3、对于第二步中不合适点做交互式拾取。自动拾取一个偶极测井记录1、选择Petrophysics Full Waveform Sonic Processing Pickin

36、g来打开来打开the tp_array_sonic_label窗口。2、选取输入起始值(SLOWNESS_S1)。让swp_projection为活动窗口，选择偶极投影道来打开阵列声波可视窗口。在信号显示区移动鼠标到拾取位置，读取显示出的X值。这个值即是从图上读取的慢度值，可作为起始值（SLOWNESS_S1），如例子中图27示。关闭阵列声波可视窗口。- 28 -图图27：确定搜索起始值：确定搜索起始值3、选择模型窗口，设置值如下图28所示 Input/Output Set选择DIPOLE，这是处理后的输出集合 SLOWNESS_WINDOW 本例中设为121 us/f（因为仅有一个波至，对于

37、偶极使用一个大点的范围是可靠的）窗口长度必须包括波至的最大波动幅度，但是也不能过大要注意排除邻区波至。 SLOWNESS_S1在前一步操作中，位置对话框里确定的X值，必须要接近值121。图28偶极自动拾取模型设置4、运行模型5、重复地3-4步二次，改变输入输出集合名字分别为FILTERED03，FILTERED05。6、设定swp_projection为活动窗口，查看曲线道（如图29示）图图29：自动拾取的：自动拾取的起始起始记录覆盖在偶极映射图像上记录覆盖在偶极映射图像上7、打开属性对话框，注意DIPOLE.DT_PICK1井的左右边界为处理的慢度起始和终止值（35 and 500 US/F

38、），输入模型中。8、关闭对话框。单极投影记录自动拾取9、确定输入tp_array_sonic_label模型所需要的的起始值(SLOWNESS_S1 and SLOWNESS_S2)选择单极图像道打开阵列声波可视窗口。在信号显示区使用鼠标定位，在定位对话框确定两个初始值，如图30所示为确定第二波至的例子。图图30：确定：确定SLOWNESS_S2的输入初始值的输入初始值10、使tp_array_sonic_label为活动窗口，对模型值做如下设置（如图31所示）： Input / Output Set 选择MONOPOLE集合，它是处理中的输出集合 Slowness_Window 设为15.2

39、4 us/f为搜索窗口大小 PICK_OPTION BOTH SLOWNESS_FROM30 SLOWNESS_TO 298.7 SLOWNESS_STEP 1.98 SLOWNESS_S1 60 SLOWNESS_S2 85图图3131：交互式拾取两个井记录模型界面：交互式拾取两个井记录模型界面11、运行模型。12、选择模板swp_projection作为查看模板，如图32所示。图图32：自动拾取井记录覆盖在单极图像道上：自动拾取井记录覆盖在单极图像道上13、菜单栏WELL下点击SAVE InteractiveInteractive PickingPicking交互式拾取14、用swp_pr

40、ojection模板，在单极投影道图像道上选择拾取记录DT_PICK2。15、选择Tools Curve Insert。16、手动的重新拾取不合适的横波波至。17、打开文本查看重命名被编辑的井记录DT_PICK2 version 2为DT_SEED。18、菜单栏WELL下点击SAVE。19、选择tp_array_sonic_label模型，对值做如下改变（如图33）： PICK OPTION SINGLE SLOWNESS_S1 使用DT_SEED作为初始记录 DT_PICK1 使用DTS_PICK图图3333：单井交互式拾取模型界面单井交互式拾取模型界面20、运行模块21、在曲线模板中查看结

41、果，使用图34作为参考。如果效果不够好，调整搜索范围（例如：SLOWNESS_WINDOW和和SLOWNESS_S1值），重复先前的步骤。图图34：偶极和单极的投影记录的自动式和交互式拾取显示偶极和单极的投影记录的自动式和交互式拾取显示在单极投影区，绿色的曲线表示作为初始记录的自动拾取曲线(DT_PICK2)。蓝色曲线是基于初始子记录的交互式拾取(DTS_PICK)。提示：处理的慢度记录需要进行深度转换，这涉及到井的参考深度。在Geolog中不会自动转换，因为对于参考深度工具，没有足够的信息可供选择。2.4 后期处理后期处理通过本节，用户将熟悉后期处理流程计算并应用频散校正，确保每个井记录在同

42、一参考系下应用频散校正导出旅行时查看相关性图形显示阵列声波数组重处理，计算平均慢度2.4.1 综述综述后期处理以及后期质量控制是声波处理的重要部分。从处理过的记录中导出的旅行时可以作为质量控制来描述声波波形波至。频散以及井眼情况等将影响最终处理结果。频散校正模型计算出频散曲线，在每个深度位置使用校正值来修正频散影响。波形将按照一定间隔被重处理（在后期处理中）来增强垂直分辨率或者补偿井眼较大的直径差值。2.4.2 频散校正频散校正这个模块输出修正后的慢度记录。它可以通过覆盖阵列道随意输出频散和频谱记录，也可以输出坏道来标记有问题区域。1、菜单栏WELL下点击SAVE。2、选择Petrophysi

43、cs Full Waveform Sonic Post-Processing Dispersion Correction。提示：离散修正是一个很费时的处理过程，本教程中为了实现这个操作只是用了少量资料来处理。3、如图35所示，设置处理参数。图图35：设置处理范围：设置处理范围4、输入参数如下：如图36所示 Input SetSONIC_DEPTH_SHIFT Output SetRAW TOOL_SPECdsi_m2 START0 STOP10000 CALIWIRE.HCAL BSWIRE.BS RHOB_F1.24 G/C3(1240 K/M3) RHOB_SWIRE.RHOZ DT_F2

44、45.8 US/F (75 US/M) DT_PMONOPOLE.DT_PICK1 DT_SDIPOLE.DT_PICK1LOG_IN从 SONIC_DEPTH_SHIFT 集合, 选择8个波形WF21_1 to WF28_1图图36：频散校正窗口：频散校正窗口5、运行模块) 处理速度取决于计算机的性能和数据量的大小。6、重复第三和第四步，改变值。 Input /Output setFILTERED03 DT_SFILTERED03.DT_PICK17、重复第三和第四步，改变值 Input /Output setFILTERED05 DT_SFILTERED05.DT_PICK18、打开模板s

45、wp_dispersion.layout。提示：离散修正结果如图37所示，处理范围为13150 - 13450 ft.图图37：离散结果对比实例（数据为：离散结果对比实例（数据为13150 - 13450 ft）1,2,3道分别表示覆盖在叠加频率域的波形离散曲线，4道表示低频滤波后的波形，5道波形大于5 Khz的频率被滤过，6道大于3 Khz频率被滤过。对原始波形做处理时，发现在做频散校正后，有不稳定的波动现象。如上例，范围如箭头所示，通过0-5 KHz带通滤波处理，修正的剪切波慢度稳定，同时1-3 KHz滤波器给出的离散都很清晰。9、菜单栏WELL下点击SAVE。2.4.3 传播时间叠加传播

46、时间叠加在前面步骤中得到的慢度记录，可以从中导出波形的旅行时。导出的旅行时可以作为质量控制参数，稍后可以将这个旅行时覆盖在波形上来形象化的检查处理结果的质量。导出旅行时的模块需要处理出来的慢度记录以及钻孔直径记录图，这些资料都必须与参考深度匹配。1、关闭所有打开的窗口，保存所做的改动。2、选择Petrophysics Full Waveform Sonic Post-Processing Traveltime打开旅行时模块。3、如图37所示，输入以下参数（如图38所示） Input/Output SetDIPOLE DTDT_PICK1_1 CALIWIRE.HCAL_1 DT_F203图图3

47、838：建立旅行时模块：建立旅行时模块4、运行模块5、打开曲线模板swp_traveltime.layout。图图39：第一个接收器计算的旅行时（从处理得到的慢度中导出）：第一个接收器计算的旅行时（从处理得到的慢度中导出）2.4.4 相关性显示相关性显示相关性图形是由慢度时间相关分析（STC）得到，如果STC处理有差错，可通过这个图件做检测分析。1、打开swp_semblance.layout.2、打开阵列声波窗口改变下列值：从FILTERED03集合增加声波阵列记录(WF21 - WF28)。在显示的Display Waveform(s) Plus里，选择Semblance，并且选择记录为F

48、ILTERED03.SEMBLANCE。选择the FILTERED03.SEMBLANCE log结果如图40所示。图图40：在阵列声波数组中显示的相关性图：在阵列声波数组中显示的相关性图3、菜单栏WELL下点击SAVE。4、关闭所有打开的窗口。2.4.5 阵列声波重处理阵列声波重处理波形将按照一定间隔被重处理（在后期处理中）来增强垂直分辨率或者补偿井眼的直径差值。重处理是基于multiple-shot processing 技术，它包含了由于与工具位置点重合而带来的冗余信息。选择若干子阵列，这样子阵列的大小可以覆盖住调查区的地层厚度。有两个方法(configurations)来进行处理，R

49、ECEIV和TRANSM，RECEIV使用的是共源极配置，可以用来提高垂向分辨率；TRANSM使用的是共接收配置。由这两个配置得来的平均慢度用于补偿因井眼大小改变而带来的影响。可使用模块Calculate Mean完成。提示：当信号波长非常长时，这个方法是无效的，厚度需要大于波长的1/4。这个技术也易于带来噪声，较少的接收器删除的噪声也少。RE_PROC可以开启和关闭，一旦打开重处理，子波列接收器NUM_R在subarray，即定义的地层厚度也可以输入，但是图像和相关性输出是不可用的。1、选择Petrophysics Full Waveform Sonic Processing Wavefor

50、mProcessing打开tp_array_sonic_process窗口。2、选择Launcher Datum设置处理参数Depth 13123.35 - Depth 13189.00 ft (4000 - 4020m)3、键入输入参数如下图41所示： Input Set SONIC_DEPTH_SHIFT Output Set DIPOLE_RECEIV4 TOOL_SPECdsi_m2 START 0 STOP 8000 RE_PROC Yes1.Hsu, K. and Cheng, S., 1987; X.M.Tang and A.Cheng, 2004 METHOD RECEIV N

51、UM_R 4 LOG_IN 选择波形SONIC_DEPTH_SHIFT WF21_1 to WF28_1图图41：使用：使用RECEIV方法重处理声波波形方法重处理声波波形4、运行模型。5、改变输出集合为DIPOLE_TRANSM4，方法改为TRANSM，重复以上步骤。6、关闭相关性图，不保存。7、对创建的投影记录执行自动拾取。选择Petrophysics Full Waveform Sonic Processing Picking来打开tp_array_sonic_label窗口。改变输入/输出集合为DIPOLE_RECEIV4使用上步的慢度值(SLOWNESS_S1) 121运行模型改变输入/输出集合为DIPOLE_TRANSM4运行模型8、打开layout swp_reprocessing，如下图42图图42：RECEIV和和 TRANSM方法得到的慢度记录对比图方法得到的慢度记录对比图9、选择Petrophysics Full Waveform Sonic Post-Processing Calculate Mean来打开tp_array_sonic mean模块。10、键入输入参数，如图43所示 Input Set DIPOLE_RECEIV4 Output Set DIPOLE_MEAN LOGIN1 DT_PICK1 LOG