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

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 设计滤波器 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 后期处理 34 2 4 1 综述 34 2 4 2 频散校正 35 2 4 3 传播时间叠加 37 2 4 4 相关性显示 38 2 4 5 阵列声波重处理 40 3 机械性质 机械性质 44 3 1 综述 44 3 2 计算动力学弹性性质 44 附录附录 I 快速运行快速运行 45 附录附录 II 频散校正讨论频散校正讨论 46 1 综述 综述 欢迎阅读Geolog软件SWB指导教程 这是帕拉代姆公司Geolog软件SWB模块的使用教程 通过典型工作流程和程序模块的使 用 使用户熟悉阵列声波处理过程的基本工具的使用 熟悉阵列声波波形 查看 创建一个阵列声波工具模版 创建一个给定波形的其他属性 解装测井记录 预处理及波形数据分析 使用平均深度 平均时间和频率滤波器去噪 去除数据采集中振幅增益和标准化 处理偶极和单极声波测井仪下波形 自动和交互式的拾取波至 频散修正 传播时间覆盖 查看结果的相关性图形 1 1 预备知识预备知识 Geolog基本测井记录中文档 所有的输入文件为Bold Courier New 输出为Bold Courier New 不加粗 1 2 数据数据 教程中使用了下面的附加文件 文件并非由软件提供 数据 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 04 2 阵列声波全波形 阵列声波全波形 典型的阵列声波工具在接收剖面上有8个接收器和3个不同的发射器 也要注意的是 不 同服务公司出产的阵列工具也不相同 在数据采集上 基于不同的调查研究目的采用不同的工作模式 在Geolog里 WF 是惯用的波形记录名字后跟两位数字 首位数代表操作模式索引 第二位数表示接收器索引 例如 WF21代表波形记录为模式2的第一记录道 基于偶极模式 每一个接收器的声波波形以数组记录存载 用Geolog layout 模版 可以在图像道和阵列 道显示一个声波波形 在这两个道里 水平轴表示信号的旅行时 垂直轴为参考深度 当选 定一个图像道后 用户可以打开一个声波阵列查看窗口 对数据做进一步分析 图图1 原始声波波形显示 原始声波波形显示 图1为原始声波波形显示图像道 第一道和第三道 阵列道 第二道和最后一道 最初的两道显示的是第一接收器的偶极波形记录 最后两道是第一接收器的单极波形记录 3 2 1 数据准备数据准备 通过本节可以熟悉声波全波形数据准备 查看 创建一个列阵声波工具模版 创建给定波形的其他属性 手动解装波形 进行深度转换 2 1 1 查看查看 创建一个声波列阵工具模版创建一个声波列阵工具模版 在波形处理前先了解工具的相关信息 它是以文件夹形式存储 斯伦贝谢公司的DSI测 得的井资料在Geolog中已经存储在Site目录下 作为四个缺省文件 dsi m1to dsi m4 当 有改动或者初次运行此工具 相关信息会保存在本工具文件夹下的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 4 图2 array sonic create toolspec 中中dsi m2的缺省值 Spec文件中的参数如下表所示 参数描述 TOOL NAME声波采集仪器的名称如 DSI XMAC TOOL SPEC声波工具模版 以模版文件命名 TOOL MODE声波工具的工作模式 通常在Text Constants表格中定义 TOOL R声波仪器的直径 SAMPLE RATE声波采样率 通常以微次生法测量 TR DISTANCE发射器与第一接收器的距离 X INTERV两个邻近接收器之间的距离 DEPTH SHIFT工具深度的平差不同于工具的参考深度和井的参考深度 负的值表示工具参考深度在 井的参考深度之下 5 如果需要的话可以对缺省值做调整 6 点击Start运行模块 一个新的模块文件在本工区的SPECS目录下产生了 7 选择Launcher Close来关闭阵列声波工具模块 8 检查新模块dsi m2的值 如图3所示 图图3 dsi m2的存储位置以及值的存储位置以及值 5 提示 不管Geolog中使用的量度制 米制的或是特定的 声波阵列工具模版总是会保存为米制的 当 要运行这个工具 首先请确认工具模版的值 以确保得到正确的结果 2 1 2 练习指导练习指导 2 创建其他波形属性创建其他波形属性 声波波形可以用反射振幅正常显示 频率谱 瞬时相位 和瞬时振幅属性可以从解释角 度很好的量化突出这些特性 原始声波频率谱显示的是不同到达时的主要频率范围 这样可 以在频率域查看各到达模式的属性特征 瞬时相使弱信号加强 但是也同样使噪音增强 瞬 时振幅是包络的振幅 包含更多侧向变化 但是降低了垂向分辨率 tp array sonic attribute模块用于计算一个给定波形的这些属性 三个测井记录可以从这个 模块输出 频率谱 瞬时相位以及瞬时振幅 可以在图像道或者阵列道中查看这些属性 1 选择Petrophysics Full Waveform Sonic Data Preparation Create Other Attributes来显示tp array sonic attribute模块 2 如图4所示 对窗口做如下值改变 Input Set改为SONIC Output Set改为SONIC ADJ 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 波形分解波形分解 来自每个接收器的声波波形有时被叠加记录到一起 这个模块用于把波形分解为每个接 收器的单个记录 通常 每个模块有自有的测深记录 这个记录要作为输出的参考深度 在 这个模块里 需要提供分解波形的名字前缀 要遵守Geolog的命名约定 为了做示范 本例 中的未分解波形都刻意不进行自动分解 1 选择Petrophysics Full Waveform Sonic Data Preparation Log Unpacking来显示tp array sonic unpack模块 2 如图4所示 对窗口做如下值改变 Input Set改为SONIC PACKED PREFIX类型为WG DEPTH选择PWD2 LOG WF选择PWG2 7 图图6 改改变变分分解解波波形形模模型型参参数数 3 点击Start运行 4 将PREFIX值变为WN LOG WF值变为PWN2 再运行模型 5 打开Text 检查未分解记录值 如下图7所示 图图7 分解波形输出记录 分解波形输出记录 2 1 4 深度转换深度转换 DSI工具数据的参考深度与井记录的参考深度不同 差值从测量工具的底部到接收器剖 面的中段进行修正 如上节的depth shift 1 选择曲线道模版swp other attributes 2 选择Tools Depth Shift Apply Logs 8 3 如图8所示 对窗口值做如下改变 Input Set 改为SONIC Output Set 改为SONIC DEPTH SHIFT DEPTH OFFSET 输入参数为 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 sonic waveforms layout 2 选择WF21的图像道 选择View Array Sonic View 9 数组声波查看窗口显示波形WF21以及它当前深度频率域的位置 见图9 要打开阵列声波查看窗口 需要选择打开带数据的图像道 阵列声波查看窗口 图图9 图图像像道道和和声声波波数数组组窗窗口口 阵列声波窗口用来详细的分析波形 它有两个显示区域 显示区在当前深度位置显示信 号波形 可以估算出时间间隔 接收器的时间延时 频率域显示区显示每个波形相应的的 频率范围 3 点击声波数组记录 Log Select 图标打开可选择多条曲线 10 提示 提示 挑选出来的记录必须来自同一个集合和模型工具 例如 有相同的模型索引 字符名字后的 第一位数相同 在图像道上显示 4 选择记录 SONIC WF21 1 to SONIC WF28 1 点击 OK 多重曲线以及其对应的频率域在数组声波的窗口显示 5 在模版上第一图像道 SONIC WF21 点击任何位置 声波数组显示窗口的信号在每 次选择深度点的时候会刷新 另一个方式是通过深度区域 键入范围值或者使用上下按键 6 设置深度点和时间范围如下所示 在显示区放大或者改变时间范围 7 将光标定位在信号显示区任意位置 按住鼠标左键拖曳定义一个矩形区域 放开鼠 标左键 通过显示的时间范围值来反映被选择的区域 将时间范围值的起始和终点值输入与用鼠标改变信号显示区域的效果是一样的 8 再次点击鼠标左键使区域变为原来大小 在波形上估算时间间隔 提示 这个功能只在多重波形在当前活动的时候起作用 当时间间隔固定时 功能不可选时首先将信号显示区放大 9 固定时间间隔 在慢度剖面显示 10 显示区 在第一个接收器 WF21WF21 的信号 按住鼠标左键 拖至到最后一个接收器 的信号处 WF28WF28 放开鼠标左键 如图10 11 图图10 时时间间间间隔隔估估算算 8个接收器的信号以及频率域在当前选择的深度显示 11 不选Moveout复选栏 不显示时间间隔数据 12 关闭阵列声波显示窗口 2 2 2 去噪去噪 波形图总是会包含噪声 太多的噪声会不能识别可检测的波至 预处理是提高数据质量 保证初至被识别的重要步骤 Geolog 有一系列基本的预处理工具可以满足一个波形的任意去噪要求 去除时间均值 在时间方向使用平均时滤波器 去除深度均值 在深度方向使用平均深度滤波器 Butterworth 滤波器 是一个频率滤波器 根据具体的需求 用户可决定针对某个波形使用对应的滤波器 1 去除时间均值噪声 这个模块提供沿时间方向去除噪声滤波 淡红色显示 估算时间间 隔 波至慢 度及接收间 隔 粗淡红色的 陡线显示第 一个和最后 一个接收器 初至时间间 隔 12 器 先在每个窗口计算出平均值 然后从原始波形移除噪声 这是用于移除低频噪声背景的滤波器 1 选择 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 打开绘图模版swp array sonic comparison layout 有低频噪声的原始 图像 运行滤波器后的结 果图 13 图图12 用时间平均滤波器滤波前后对比图 用时间平均滤波器滤波前后对比图 2 去除深度均值噪声 这个模块提供沿深度方向去除噪声滤波器 先在每个活动窗口计算出平均值 然后从原 始波形移除噪声 这个滤波器可去除典型的诸如套管信号或者仪器带来的噪声 1 打开井Polaris 04 保存对先前井所作的改变 2 关闭所有打开窗口 3 选择Petrophysics Full Waveform Sonic Pre Processing RemoveDepth Average 4 在tp array sonic depth average窗口对值做如下改变 如图13 Input Set改为SONIC WINDOW设置滤波器长度为21 LOG INWF41 图图13 设设置置平平均均深深度度滤滤波波器器模模型型值值 5 点击Start运行模型 6 打开绘图模版swp depth average layout 套管和仪器形成的信号可以很强 使得重要的波至难以识别 图14左边道显示的波形是 一个这样的例子 直直的垂线占主要部分 右边的道显示的是不需要信息被移除后的清晰的 波至 14 图图14 用深度滤波器滤波前后对比图 用深度滤波器滤波前后对比图 3 使用频率滤波器 Butterworth 滤波器是一个频率滤波器 可以通过给定一个高的和低的截止值以及创建 滤过顺序来实现 Butterworth 滤波器选择的顺序 从1到8 决定了从通频带到抑制带的过 渡陡度 单位设为千赫兹而不是赫兹 是为了便于进行声波波形的处理 如果数据是来自于 标准的DSI工具采集 可以忽略采样率 因为Geolog可以自动地读取 有关滤波器的调整适应更多信息参考菜单选项下的Design Filter 1 保存对Polaris 04的改动 打开井Polaris 03 2 选择Petrophysics Full Waveform Sonic Pre Processing Frequency 有低频噪声的 原始图像 运行滤波器后的 结果图 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 layout 6 选择View Array Sonic View来查看滤波后图形 7 设置当前深度到如图16所示 在3HZ以上的频率已经在当前深度位置被移除 16 图图16 WF滤波后图形滤波后图形 8 关闭阵列声波查看窗口 9 选择tp array sonic filter模块 重复第二步和第三步改变频率范围为0 5赫兹 输 出名为FILTERED05 10 点击Start运行模块 11 打开swp frequency05 layout 12 打开声波数组查看窗口 大于5Hz的频率信号都被移除 如图17所示 在当前深度位置显示的频率谱 17 图图17 滤除高于 滤除高于5Hz波后的波后的WF21波形图波形图 13 关闭阵列声波窗口 2 2 3 设计滤波器设计滤波器 Butterworth滤波器将作为频率滤波器用在下面的例子中 Butterworth 是一个很好的带 通滤波器 在损失了部分频率域后有较为合理的时间衰减 高频和低频截止值用于带通和带 止 对于低通滤波器 高频波被截止 反之亦然 Butterworth 滤波器的排列控制着在通频 带之外信号变弱的快慢 设计一个合适的滤波器或是子波 必须要注意采样率和频率值 设计一个Butterworth子 波 采样率 SR 要使用0 001秒 这样可以得到Nyquist frequency 1 2 SR 500 Hz 输入频率必须要小于尼奎斯特频率 Nyquist frequency 频率 也就是说 如果提供 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 查看创建的滤波器 如图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 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 阵列声波处理阵列声波处理 通过本节 用户可以 处理偶极波形 处理单极波形 21 使用自动交互式拾取波至 2 3 1 处理模块简介处理模块简介 声波波形处理基于慢度时间相关法 这个模块在活动窗口中通过时间和慢度方向可以查 看多个声波 默认的 时间范围可以包括所有的记录时间 慢度的范围至少包括目的区 例 如 慢度范围为100 800 us m 1200 10000 m s 可以包含大多数区域的纵波波至 建议用 较大的慢度范围这样可以尽可能包含所有波至 但这也带来更多的计算量 首先选择少量数据来决定时间和慢度方向的最优选择 要注意慢度步长的选择 过大的步长会降低最终声波分辨率 在每个深度点计算出的最大值作为与深度相关的函数映射到慢度轴上 这是最终输出 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 TOOL 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 plot 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 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 阵列声波数组窗口查看图件 阵列声波数组窗口查看图件 13 关闭相关性图件 不保存改变 14 打开模版swp projection layout 上面例子中的偶极和单极子映射记录如图26所示 26 偶极波形处理后图形 单极子波形处理后图形 水平轴显示慢度 START 到 END 图图26 深度相关函数映射 深度相关函数映射 15 在WELL下点击SAVE 2 3 4 拾取标志波至拾取标志波至 映射记录在慢度 深度平面内代表了一个或多个相关性 它们与某些波至有较好的一致 性 模型在输入的映射记录中自动或者交互式的寻找相关最大值位置 从慢度投影记录拾取 相应一组波至 搜索从一个起始值开始 初始慢度值在映射记录里 来标定某一初至 搜寻 范围以从映射记录估计的起始值为中心 估算输入起始值最简便方法是将投影记录插入到图 像道 打开阵列生波可视窗口 移动鼠标标定出信号显示区的目的位置 光标位置X的值会 在位置对话框的状态栏显示出来 这个值就是从图上读取的慢度值 它可以用来作为搜索起 始值 对于交互式拾取 测井曲线首先是自动被拾取 参见Curve Insert联机帮助 编辑 最后作为备用初始记录被保存 主要有 1 自动拾取一个偶极测井记录 2 自动拾取一个单极测井记录 27 3 对于第二步中不合适点做交互式拾取 自动拾取一个偶极测井记录 1 选择Petrophysics Full Waveform Sonic Processing Picking来打开来打开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 因为仅有一个波至 对 于偶极使用一个大点的范围是可靠的 窗口长度必须包括波至的最大波 动幅度 但是也不能过大要注意排除邻区波至 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 输入模型中 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 24 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 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 OPTION SINGLE SLOWNESS S1 使用DT SEED作为初始记录 DT PICK1 使用DTS PICK 图图3333 单井交互式拾取模型界面单井交互式拾取模型界面 20 运行模块 21 在曲线模板中查看结果 使用图34作为参考 如果效果不够好 调整 搜索范围 例如 SLOWNESS WINDOW和和SLOWNESS S1值 重复先前 的步骤 图图34 偶极和单极的投影记录的自动式和交互式拾取显示偶极和单极的投影记录的自动式和交互式拾取显示 在单极投影区 绿色的曲线表示作为初始记录的自动拾取曲线 DT PICK2 蓝色曲线是基于初始子记录的交互式拾取 DTS PICK 提示 处理的慢度记录需要进行深度转换 这涉及到井的参考深度 在Geolog中不会 自动转换 因为对于参考深度工具 没有足够的信息可供选择 2 4 后期处理后期处理 通过本节 用户将熟悉后期处理流程 计算并应用频散校正 确保每个井记录在同一参考系下 应用频散校正 导出旅行时 查看相关性图形显示 阵列声波数组重处理 计算平均慢度 2 4 1 综述综述 后期处理以及后期质量控制是声波处理的重要部分 从处理过的记录中导 出的旅行时可以作为质量控制来描述声波波形波至 频散以及井眼情况等将影 响最终处理结果 频散校正模型计算出频散曲线 在每个深度位置使用校正值来修正频散影 响 波形将按照一定间隔被重处理 在后期处理中 来增强垂直分辨率或者补 偿井眼较大的直径差值 2 4 2 频散校正频散校正 这个模块输出修正后的慢度记录 它可以通过覆盖阵列道随意输出频散和 频谱记录 也可以输出坏道来标记有问题区域 1 菜单栏WELL下点击SAVE 2 选择Petrophysics 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 F245 8 US F 75 US M DT PMONOPOLE DT PICK1 DT SDIPOLE DT PICK1 LOG IN从 SONIC DEPTH SHIFT 集合 选择8 个波形WF21 1 to WF28 1 图图36 频散校正窗口 频散校正窗口 5 运行模块 处理速度取决于计算机的性能和数据量的大小 6 重复第三和第四步 改变值 Input Output setFILTERED03 DT SFILTERED03 DT PICK1 7 重复第三和第四步 改变值 Input Output setFILTERED05 DT SFILTERED05 DT PICK1 8 打开模板swp 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 传播时间叠加传播时间叠加 在前面步骤中得到的慢度记录 可以从中导出波形的旅行时 导出的旅行 时可以作为质量控制参数 稍后可以将这个旅行时覆盖在波形上来形象化的检 查处理结果的质量 导出旅行时的模块需要处理出来的慢度记录以及钻孔直径记录图 这些资 料都必须与参考深度匹配 1 关闭所有打开的窗口 保存所做的改动 2 选择Petrophysics Full Waveform Sonic Post Processing Traveltime打开旅行时模块 3 如图37所示 输入以下参数 如图38所示 Input Output SetDIPOLE DTDT PICK1 1 CALIWIRE HCAL 1 DT F203 图图3838 建立旅行时模块 建立旅行时模块 4 运行模块 5 打开曲线模板swp traveltime layout 图图39 第一个接收器计算的旅行时 从处理得到的慢度中导出 第一个接收器计算的旅行时 从处理得到的慢度中导出 2 4 4 相关性显示相关性显示 相关性图形是由慢度时间相关分析 STC 得到 如果STC处理有差错 可通过这个图件做检测分析 1 打开swp semblance layout 2 打开阵列声波窗口改变下列值 从FILTERED03集合增加声波阵列记录 WF21 WF28 在显示的Display Waveform s Plus里 选择Semblance 并且选择记 录为FILTERED03 SEMBLANCE 选择the FILTERED03 SEMBLANCE log 结果如图40所示 图图40 在阵列声波数组中显示的相关性图 在阵列声波数组中显示的相关性图 3 菜单栏WELL下点击SAVE 4 关闭所有打开的窗口 2 4 5 阵列声波重处理阵列声波重处理 波形将按照一定间隔被重处理 在后期处理中 来增强垂直分辨率或者补 偿井眼的直径差值 重处理是基于multiple shot processing 技术 它包含了 由于与工具位置点重合而带来的冗余信息 选择若干子阵列 这样子阵列的大 小可以覆盖住调查区的地层厚度 有两个方法 configurations 来进行处理 RECEIV和TRANSM RECEIV使用的是共源极配置 可以用来提高垂向分辨 率 TRANSM使用的是共接收配置 由这两个配置得来的平均慢度用于补偿因 井眼大小改变而带来的影响 可使用模块Calculate Mean完成 提示 当信号波长非常长时 这个方法是无效的 厚度需要大于波长的1 4 这个技 术也易于带来噪声 较少的接收器删除的噪声也少 RE PROC可以开启和关闭 一旦打开重处理 子波列接收器NUM R在 subarray 即定义的地层厚度也可以输入 但是图像和相关性输出是不可用的 1 选择Petrophysics Full Waveform Sonic Processing Waveform Processing打开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 Yes 1 Hsu K and Cheng S 1987 X M Tang and A Cheng 2004 METHOD RECEIV NUM 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 LOGIN2 DIPOLE TRANSM4 DT PI