RadEx Pro处理电火花数据操作步骤(上)主要讲述RadEx Pro读取电火花数据,查看数据的质量,以及简单的滤波和振幅纠正。
6、海底拾取,建立流程060 seafloor pick
Trace Input加载stack数据集
Trace Header Math,设置海底起始ms数
正在上传…重新上传取消
First Break Picking,拾取海底。
选择Global(全局),Threshold设置为50,选择Min或者Max。
Derivative(派生),Window to calculate derivative设置为50
效果如下:
如果你觉得海底高了或者低了,可以调整海底的高度。
Trace Header Math,向上移动 1.3ms。
我们只是流程里修改了 FBPICK道头值并观察了其结果,我们需要将结果保存,可以Trace Output,但鉴于数据本身并没有在这修改,我们最好是只保存道头而不用另外建立一个数据集。我们可以用Header<->Dataset Transfer,选择 stack数据集用来保存道头,设置 CDP为 相匹配的原始记录,FBPICK为赋值记录。
7、除鬼影,建立一个流程070 deghosting。
Trace Input选择 stack dataset, sort为 CDP。
Trace Editing模块实现顶切,在 Trace Editing Parameters选项卡中设置Muting类型为 Top muting,Taper windowlength设置为 1ms。
Horizon选项卡指定Trace header为FBPICK。
顶切结果如下图。海底之上的振幅为0。
放大海底某处反射,估计海底真值和鬼波之间的延迟,大约为1ms。
Editing Trace 前插入 Predictive Deconvolution(预测反褶积)模块 。设置 Prediction Gap 为 1.0ms,White Noise Level 为0.10
结果如下图。结果较原来有一些噪声。
在 Predictive Deconvolution和Trace Editing之间插入带通滤波模块。
Trace Output模块保存结果,数据集命名为stack_dgh。
使用Deghosting替换 Predictive Deconvolution,为了便于对比,我将未做反褶积、Predictive Deconvolution和Deghosting的图像做个gif图,结果如下图。
8、建立流程080 postprocessing。
最后,我们可能希望最终叠加频谱来稍微提高一下分辨率并让结果成图更好看,有更多的 “关注点”。有几个模块可以用来扩展和白化数据频谱,包括谱白化、谱整形和一些典型的反褶积。
F-K Amplitude Power,可提高 F-X域(频率-空间域)或 F-K域(频率-波数域)的振幅频谱到任意能量范围并进行傅里叶逆变换( 2D or 1D)到原始 )到原始 T-X域(时间-空间域)。取决于指数和应用的不同,算法结果从压制噪音到谱白化都会不同:大体上,谱白化的获得指数值<1,当指数值 >1时频谱窄带会抑制噪声。选择 F-K和 F-X域是很灵活的。 (这一段完全看不懂,照抄。)
谱白化(Spectral Whitening)是地震勘探数据处理的一种常用方法,可以展宽信号的频谱提高信号分辨能力,是一种“纯振幅”的滤波过程。实现方法:谱白化处理既可以在频率域中实现,也可以在时间域中实现。频率域谱白化方法的基本方法是将信号进行傅里叶变换,由时间域变到频率域;再在有效的频率范围内进行分频率;对分频信号进行傅里叶反变换,由频率域变到时间域;对各时间域信号进行时变增益;最后将时变增益信号重新合成,生成白化后的时间域信号。
Exponent 设置为0.5。
Trace Output模块保存结果,数据集命名为stack_final。
结果如下图。
9、打印处理结果,建立一个新的流程090 plotting。
RadEx Pro通过Pick来画地层,地层画好以后,先设置Pick的Header参数,即位置,然后导出Pick即可。
RadEx Pro提供了一个视频“Demultiple of single-channel boomer data in RadExPro 2016”,利用Zero-Offset Demultiple模块压制单道的二次波,下面是前后对比图。
从上面的对比图可看出,效果真好,但是悲催的是,我们完全按照它的设置,二次波怎么也不能压制。
结果见下图,二次波清晰可见,这是什么道理?完全搞不懂。
后来我们找到了中科院海洋所的钱老师,在他的指导和启发下,修改了一些参数,并查看了运行效果以后,总结了Zero-Offset Demultiple模块各参数的影响,见下图。
将Number of traces设置成50,High Frequency设置为1500,结果如下图。效果勉强凑合,不如上面视频的效果。
小结
从上面的图像来看,这个海洋高分辨率多道地震数据处理同样不能有效压制二次波,其处理结果与电火花自带软件GeoSuite相当。如果使用RadEx Pro专门单道地震的二次压制模块Zero-Offset Demultiple,效果勉强凑合。RadEx Pro主要是用来处理多道地震的,相对于GeoSuite而言,其信号处理和反褶积方面的工具更为丰富,但是,在使用的便利性上,不如GeoSuite。