资源描述
Geolog软件技术手册
Full Sonic Wave Processing -SWB
帕拉代姆公司北京代表处
2023年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.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模块的使用教程,通过典型工作流程和程序模块的使用,使用户熟悉阵列声波解决过程的基本工具的使用。
l 熟悉阵列声波波形
l 查看/创建一个阵列声波工具模版
l 创建一个给定波形的其他属性
l 解装测井记录
l 预解决及波形数据分析
l 使用平均深度、平均时间和频率滤波器去噪
l 去除数据采集中振幅增益和标准化
l 解决偶极和单极声波测井仪下波形
l 自动和交互式的拾取波至
l 频散修正
l 传播时间覆盖
l 查看结果的相关性图形
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
函数:
图像:
井: polaris_03, polaris_04
2、阵列声波全波形
典型的阵列声波工具在接受剖面上有8个接受器和3个不同的发射器。也要注意的是,不同服务公司出产的阵列工具也不相同。
在数据采集上,基于不同的调查研究目的采用不同的工作模式。在Geolog里,“WF”是惯用的波形记录名字后跟两位数字,首位数代表操作模式索引,第二位数表达接受器索引。例如,WF21代表波形记录为模式2的第一记录道(基于偶极模式)。
每一个接受器的声波波形以数组记录存载。用Geolog layout(模版)可以在图像道和阵列道显示一个声波波形。在这两个道里,水平轴表达信号的旅行时,垂直轴为参考深度。当选定一个图像道后,用户可以打开一个声波阵列查看窗口,对数据做进一步分析。
图1:原始声波波形显示
图1为原始声波波形显示图像道(第一道和第三道),阵列道(第二道和最后一道)。最初的两道显示的是第一接受器的偶极波形记录,最后两道是第一接受器的单极波形记录。
2.1数据准备
通过本节可以熟悉声波全波形数据准备
l 查看/创建一个列阵声波工具模版
l 创建给定波形的其他属性
l 手动解装波形
l 进行深度转换
2.1.1查看/创建一个声波列阵工具模版
在波形解决前先了解工具的相关信息,它是以文献夹形式存储。斯伦贝谢公司的DSI测得的井资料在Geolog中已经存储在Site目录下,作为四个缺省文献(dsi_m1 to 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
图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的存储位置以及值
提醒:不管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所示。
原始数据 原始数据频率谱 原始数据瞬时相位 原始数据瞬时振幅
图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
图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
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 解决
l 通过本节可以使用户熟悉全波形声波解决流程
l 在声波数组查看窗口操作和解释数据
l 使用平均时间和平均深度以及Butterworth频率滤波器去除波形噪音
l 设计一个滤波器
l 去除数据采集中的振幅增益和标准化
2.2.1 数据分析
为了保证解决结果质量,在数据解决前要检查数据保证无误。波形数据通常包含噪声。通常需要选择适当的工具来尽也许多的去噪。
1、打开模板swp_raw_sonic_waveforms.layout。
2、选择WF21的图像道,选择View > Array Sonic View。
数组声波查看窗口显示波形WF21以及它当前深度频率域的位置。(见图9)
) 要打开阵列声波查看窗口,需要选择打开带数据的图像道。
阵列声波查看窗口
图9:图像道和声波数组窗口
阵列声波窗口用来具体的分析波形,它有两个显示区域。显示区在当前深度位置显示信号波形,可以估算出时间间隔(接受器的时间延时)。频率域显示区显示每个波形相应的的频率范围。
3、点击声波数组记录- Log Select 图标打开可选择多条曲线。
提醒:挑选出来的记录必须来自同一个集合和模型工具(例如:有相同的模型索引-字符名字后的第一位数相同)在图像道上显示。
4、选择记录SONIC.WF21_1 to SONIC.WF28_1,点击OK。
多重曲线以及其相应的频率域在数组声波的窗口显示。
5、在模版上第一图像道(SONIC.WF21)点击任何位置。声波数组显示窗口的信号在每次选择深度点的时候会刷新。
) 另一个方式是通过深度区域,键入范围值或者使用上下按键。
6、设立深度点和时间范围如下所示:
在显示区放大或者改变时间范围
7、将光标定位在信号显示区任意位置,按住鼠标左键拖曳定义一个矩形区域,放开鼠标左键。通过显示的时间范围值来反映被选择的区域。
) 将时间范围值的起始和终点值输入与用鼠标改变信号显示区域的效果是同样的。
8、再次点击鼠标左键使区域变为本来大小。
在波形上估算时间间隔:
提醒:这个功能只在多重波形在当前活动的时候起作用。
)当时间间隔固定期,功能不可选时一方面将信号显示区放大。
9、固定期间间隔,在慢度剖面显示。
10、显示区,在第一个接受器(WF21)的信号,按住鼠标左键,拖至到最后一个接受器的信号处(WF28),放开鼠标左键。如图10
粗淡红色的陡线显示第一个和最后一个接受器初至时间间隔。
淡红色显示估算时间间隔,波至慢度及接受间隔
图10:时间间隔估算
8个接受器的信号以及频率域在当前选择的深度显示。
11、不选Moveout复选栏,不显示时间间隔数据。
12、关闭阵列声波显示窗口。
2.2.2去噪
波形图总是会包含噪声,太多的噪声会不能辨认可检测的波至。预解决是提高数据质量保证初至被辨认的重要环节。
Geolog 有一系列基本的预解决工具可以满足一个波形的任意去噪规定。
l 去除时间均值,在时间方向使用平均时滤波器
l 去除深度均值,在深度方向使用平均深度滤波器
l Butterworth 滤波器,是一个频率滤波器
根据具体的需求,用户可决定针对某个波形使用相应的滤波器
(1)去除时间均值噪声
这个模块提供沿时间方向去除噪声滤波器。先在每个窗口计算出平均值,然后从原始波形移除噪声。
) 这是用于移除低频噪声背景的滤波器。
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。
有低频噪声的原始图像
运营滤波器后的结果图
图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_IN WF41
图13:设立平均深度滤波器模型值
5、点击Start运营模型
6、打开绘图模版swp_depth_average.layout。
套管和仪器形成的信号可以很强,使得重要的波至难以辨认。图14左边道显示的波形是一个这样的例子,直直的垂线占重要部分。右边的道显示的是不需要信息被移除后的清楚的波至。
运营滤波器后的结果图
有低频噪声的原始图像
图14:用深度滤波器滤波前后对比图
(3)使用频率滤波器
Butterworth 滤波器是一个频率滤波器,可以通过给定一个高的和低的截止值以及创建滤过顺序来实现。Butterworth 滤波器选择的顺序(从1到8)决定了从通频带到克制带的过渡陡度。单位设为千赫兹而不是赫兹,是为了便于进行声波波形的解决。假如数据是来自于标准的DSI工具采集,可以忽略采样率,由于Geolog可以自动地读取。
有关滤波器的调整适应更多信息参考菜单选项下的Design Filter。
1、保存对Polaris_04的改动,打开井Polaris_03。
2、选择Petrophysics > Full Waveform Sonic > Pre-Processing > Frequency Filter.
3、对打开窗口做如下改变(如图15所示)
— Input Set SONIC_DEPTH_SHIFT
— Output Set 改为FILTERED03
— HIGH_END 设为3HZ(移除所有高于这个值之外的频率信号)
— SAMPLE_RATE 10
— 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:WF滤波后图形
8、关闭阵列声波查看窗口。
9、选择tp_array_sonic_filter模块,反复第二步和第三步改变频率范围为0-5赫兹,输出名为FILTERED05。
10、点击Start运营模块
11、打开swp_frequency05.layout
12、打开声波数组查看窗口。
大于5Hz的频率信号都被移除,如图17所示,在当前深度位置显示的频率谱。
图17:滤除高于5Hz波后的WF21波形图
13、关闭阵列声波窗口。
2.2.3 设计滤波器
Butterworth滤波器将作为频率滤波器用在下面的例子中。Butterworth 是一个很好的带通滤波器,在损失了部分频率域后有较为合理的时间衰减。高频和低频截止值用于带通和带止。对于低通滤波器,高频波被截止,反之亦然。Butterworth 滤波器的排列控制着在通频带之外信号变弱的快慢。
设计一个合适的滤波器或是子波,必须要注意采样率和频率值。设计一个Butterworth子波,采样率(SR)要使用0.001秒,这样可以得到Nyquist frequency = 1 / (2 * SR) = 500 Hz。输入频率必须要小于尼奎斯特频率(Nyquist frequency)频率。也就是说,假如提供的值大于尼奎斯特频率(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:运营创建滤波器模块后的结果
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
— 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阵列声波解决
通过本节,用户可以:
l 解决偶极波形
l 解决单极波形
l 使用自动交互式拾取波至
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:偶极波形解决模型窗口
5、点击Start运营模型。
当模型运营结束,结果以相关性图件显示,如图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
— 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:阵列声波数组窗口查看图件
13、关闭相关性图件,不保存改变。
14、打开模版swp_projection.layout。
上面例子中的偶极和单极子映射记录如图26所示。
偶极波形解决后图形 单极子波形解决后图形 水平轴显示慢度(START)到(END)
图26:深度相关函数映射
15、在WELL下点击SAVE.
2.3.4 拾取标志波至
映射记录在慢度-深度平面内代表了一个或多个相关性。它们与某些波至有较好的一致性。模型在输入的映射记录中自动或者交互式的寻找相关最大值位置,从慢度投影记录拾取相应一组波至。搜索从一个起始值开始(初始慢度值在映射记录里)来标定某一初至,搜寻范围以从映射记录估计的起始值为中心。估算输入起始值最简便方法是将投影记录插入到图像道,打开阵列生波可视窗口,移动鼠标标定出信号显示区的目的位置,光标位置X的值会在位置对话框的状态栏显示出来。这个值就是从图上读取的慢度值,它可以用来作为搜索起始值。
对于交互式拾取,测井曲线一方面是自动被拾取(参见Curve Insert联机帮助),编辑,最后作为备用初始记录被保存。
重要有:
1、自动拾取一个偶极测井记录
2、自动拾取一个单极测井记录
3、对于第二步中不合适点做交互式拾取。
★自动拾取一个偶极测井记录
1、选择Petrophysics > Full Waveform Sonic > Processing > Picking来打开the tp_array_sonic_label窗口。
2、选取输入起始值(SLOWNESS_S1)。
—让swp_projection为活动窗口,选择偶极投影道来打开阵列声波可视窗口。
—在信号显示区移动鼠标到拾取位置,读取显示出的X值。这个值即是从图上读取的慢度值,可作为起始值(SLOWNESS_S1),如例子中图27示。
—关闭阵列声波可视窗口。
图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_FROM 30
— SLOWNESS_TO 298.7
— SLOWNESS_STEP 1.98
— SLOWNESS_S1 60
— SLOWNESS_S2 85
图31:交互式拾取两个井记录模型界面
11、运营模型。
12、选择模板swp_projection作为查看模板,如图32所示。
图32:自动拾取井记录覆盖在单极图像道上
13、菜单栏WELL下点击SAVE
★Interactive 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 OPTION SINGLE
— SLOWNESS_S1 使用DT_SEED作为初始记录
— DT_PICK1 使用DTS_PICK
图33:单井交互式拾取模型界面
20、运营模块
21、在曲线模板中查看结果,使用图34作为参考。假如效果不够好,调整搜索范围(例如:SLOWNESS_WINDOW和SLOWNESS_S1值),反复先前的环节。
图34:偶极和单极的投影记录的自动式和交互式拾取显示
在单极投影区,绿色的曲线表达作为初始记录的自动拾取曲线(DT_PICK2)。蓝色曲线是基于初始子记录的交互式拾取(DTS_PICK)。
提醒:解决的慢度记录需要进行深度转换,这涉及到井的参考深度。在Geolog中不会自动转换,由于对于参考深度工具,没有足够的信息可供选择。
2.4后期解决
通过本节,用户将熟悉后期解决流程
l 计算并应用频散校正,保证每个井记录在同一参考系下
l 应用频散校正
l 导出旅行时
l 查看相关性图形显示
l 阵列声波数组重解决,计算平均慢度
2.4.1综述
后期解决以及后期质量控制是声波解决的重要部分。从解决过的记录中导出的旅行时可以作为质量控制来描述声波波形波至。频散以及井眼情况等将影响最终解决结果。
频散校正模型计算出频散曲线,在每个深度位置使用校正值来修正频散影响。波形将按照一定间隔被重解决(在后期解决中)来增强垂直分辨率或者补偿井眼较大的直径差值。
2.4.2频散校正
这个模块输出修正后的慢度记录。它可以通过覆盖阵列道随意输出频散和频谱记录,也可以输出坏道来标记有问题区域。
1、菜单栏WELL下点击SAVE。
2、选择Petrophysics > Full Waveform Sonic > Post-Processing > Dispersion Correction。
提醒:离散修正是一个很费时的解决过程,本教程中为了实现这个操作只是用了少量资料来解决。
3、如图35所示,设立解决参数。
图35:设立解决范围
4、输入参数如下:如图36所示
— Input Set SONIC_DEPTH_SHIFT
— Output Set RAW
— TOOL_SPEC dsi_m2
— START 0
— STOP 10000
— CALI WIRE.HCAL
— BS WIRE.BS
— RHOB_F 1.24 G/C3(1240 K/M3)
— RHOB_S WIRE.RHOZ
— DT_F 245.8 US/F (75 US/M)
— DT_P MONOPOLE.DT_PICK1
— DT_S DIPOLE.DT_PICK1
—LOG_IN 从 SONIC_DEPTH_SHIFT 集合, 选择8个波形WF21_1 to WF28_1
图36:频散校正窗口
5、运营模块
) 解决速度取决于计算机的性能和数据量的大小。
6、反复第三和第四步,改变值。
— Input /Output set FILTERED03
— DT_S FILTERED03.DT_PICK1
7、反复第三和第四步,改变值
— Input /Output set FILTERED05
— DT_S FILTERED05.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 Set DIPOLE
— DT DT_PICK1_1
— CALI WIRE.HCAL_1
— DT_F 203
图38:建立旅行时模块
4、运营模块
5、打开曲线模板swp_traveltime.layout。
图39:第一个接受器计算的旅行时(从解决得到的慢度中导出)
2.4.4 相关性显示
相关性图形是由慢度时间相关分析(STC)得到,假如STC解决有差错,可通过这个图件做检测分析。
1、打开swp_semblance.lay
展开阅读全文