石油地球物理勘探  2019, Vol. 54 Issue (6): 1237-1245  DOI: 10.13810/j.cnki.issn.1000-7210.2019.06.007
0
文章快速检索     高级检索

引用本文 

刘斌, 宋智强, 王超, 陈洪芝. 顺北沙漠区地震资料陷波现象原因分析及消除. 石油地球物理勘探, 2019, 54(6): 1237-1245. DOI: 10.13810/j.cnki.issn.1000-7210.2019.06.007.
LIU Bin, SONG Zhiqiang, WANG Chao, CHEN Hongzhi. Trapping-wave phenomenon elimination on seismic data in Shunbei desert area. Oil Geophysical Prospecting, 2019, 54(6): 1237-1245. DOI: 10.13810/j.cnki.issn.1000-7210.2019.06.007.

本项研究受中石化石油工程公司项目"顺北沙漠区高密度空间采样先导试验"(SG18-84X)资助

作者简介

刘斌  高级工程师, 1969年生; 1992年获西南石油大学勘查地球物理专业学士学位; 现在中国石化石油工程地球物理公司胜利分公司从事地震数据采集方法研究和管理工作, 在国内外专业期刊上发表相关论文30余篇

宋智强, 山东省东营市东营区牛庄镇中国石化石油工程地球物理公司胜利分公司研究所, 257086。Email:15577209@qq.com

文章历史

本文于2019年1月18日收到,最终修改稿于同年8月28日收到
顺北沙漠区地震资料陷波现象原因分析及消除
刘斌 , 宋智强 , 王超 , 陈洪芝     
① 中国石化石油工程地球物理公司胜利分公司研究所, 山东东营 257086;
② 中国石化西北油田分公司, 新疆轮台 841600;
③ 中国石油华北油田公司勘探开发研究院, 河北任丘 062552
摘要:位于塔克拉玛干沙漠腹地的顺北沙漠区,单炮地震资料中普遍存在陷波现象,陷波频率约为15Hz,正好位于奥陶系目的层有效频带内,对地震成像影响较大。针对此现象进行了理论剖析,发现是由于单炮记录中存在能量较强的浅层多次折射波,与地下有效反射波混叠在一起,在频谱中表现为与地下反射波主频之间形成"两高夹一低"的形态,形成了视觉上的"陷波"现象。采用分离波场的方法,分离出浅层多次折射波,可消除"陷波"现象的影响,提高地震资料的品质。
关键词顺北沙漠区    陷波    频谱    近地表    浅层多次折射波    信噪比    压制方法    
Trapping-wave phenomenon elimination on seismic data in Shunbei desert area
LIU Bin , SONG Zhiqiang , WANG Chao , CHEN Hongzhi     
① Shengli Branch, Geophysical Corporation, SINOPEC, Dongying, Shandong 257086, China;
② Northwest Oilfield Branch Co., SINOPEC, Luntai, Xinjiang 841600, China;
③ Exploration and Development Research Institute, PetroChina Huabei Oilfield Company, Renqiu, Hebei 062552, China
Abstract: Trapping-wave phenomena are quite common on seismic shot records in Shunbei desert area, Taklimakan Desert hinterland.The frequency of trapping wave is around 15Hz just within the band of exploration targets, which affects greatly the seismic imaging.Our theoretical analysis demonstrates that these phenomena are caused by strong-energy shallow multi-refractions.These multi-refractions are mixed with underground reflections.In the spectrum, it shows a form of "two high sandwiched one low" with the dominant frequency of underground reflections, seems like a visual "trapping-wave" phenomenon.With the wavefield separation, shallow multi-refractions are separated and eliminated, and the seismic imaging is greatly improved.
Keywords: Shunbei area    trapping wave    spectrum    near-surface    shallow multi-refractions    signal-to-noise ratio (SNR)    suppression method    
1 问题的提出

顺北沙漠区是中国石化西部重要的油气勘探区,奥陶系目的层较深,地表为沙漠,因此地震资料有效频带范围较窄。而该地区地震单炮中普遍存在“陷波”现象,“陷波”点位于12~15Hz内,正好位于地震资料有效频带范围内(图 1)。

图 1 顺北沙漠区地震资料单炮频谱分析 (a)不同炸药量试验单炮频谱对比;(b)不同井深试验单炮频谱对比

以往认为陷波现象可能与地震波的虚反射有关[1-3]。首先根据微测井数据近地表调查数据,建立该区的近地表速度模型,再通过虚反射原理计算陷波点的位置,可以分析这种陷波现象是否与地震波虚反射有关。

图 2为微测井解释结果,位置与与图 1试验炮相同。可以看出,该处近地表为两层结构,低速层厚度为2.48m,速度为344m/s,降速层厚度为3.87m,速度为515m/s,高速层顶界面深度为6.35m,高速层速度为1678m/s。图 3为该处近地表速度模型。虚反射是由于炸药在高速层顶面以下激发时,地震波向上传播,遇到强波阻抗界面(高速层顶界面)向下反射,与向下传播的地震波叠合,由于两个地震波存在时差,对部分频率进行了压制[4-5]。采用虚反射原理计算陷波点的位置[6](图 4),第一个陷波点位于60Hz,与实际资料中15Hz陷波点不符。另外根据该工区表层厚度为21.2m的微测井数据,计算的陷波点位置是28Hz,也与实际资料陷波点位置不符。因此,该陷波现象并非由虚反射引起。

图 2 微测井解释结果

图 3 建立的近地表模型

图 4 理论计算的虚反射陷波现象
2 陷波现象成因分析 2.1 地震记录不同区段的频谱分析

进一步对单炮地震记录沿着反射波同相轴开较小的时窗(图 5左方框所示)进行频谱分析,并没有出现陷波现象。

图 5 单炮记录反射波同相轴(左图中黄色区域)的频谱(右)

另外,从单炮记录中可以明显看到浅层多次折射波(图 6左),多次折射波速度约为1700m/s,与高速层速度相近并且与初至波平行,是在低降速带内经多次折射滑行的多次波,其主频为10Hz(图 6右)。从图 7所示单炮记录可见,浅层多次折射波不仅存在于炮记录的浅层,而且贯穿于整个记录,在深层与有效反射波的波场混叠在一起(图 7方框所示)。

图 6 浅层多次折射波(左)及频谱分析(右)

图 7 单炮中的浅层多次折射波

根据图 5~图 7,可能是分析时窗包含了反射波和浅层多次折射波,频谱中出现了两个主频峰值、导致了峰值间低凹处出现了视觉陷波现象。

将浅层多次折射波和有效反射波频谱叠合显示(图 8),可以明显看到两个频谱峰值,一个是浅层多次折射波的主频,而另一个是有效反射波的主频,该叠合显示的频谱与实际资料中放大时窗所显示的频谱特征相同(图 9),陷波点的频率都为15Hz,说明由于近地表中的浅层多次折射能量太强,与地下有效反射波的混叠,使频谱分析中产生了“两高夹一低”的现象,形成了视觉上的陷波现象。

图 8 两种波场的频谱叠合显示

图 9 时窗包括两种波场的频谱
2.2 浅层折射波的正演模拟

为了进一步说明问题,通过微测井数据(图 10a)采用VSP速度解释方法,即拾取初至时间并计算井中相邻两炮间的速度曲线[7-11](图 10b),建立精细的近地表模型(图 10c)。将近地表模型与地下地质模型结合建立该区的综合模型(图 11),运用波动方程进行正演模拟。从模拟记录(图 12a)可以看到明显的浅层多次折射波,对该模拟单炮进行频谱分析,可以看到频谱中明显存在两个峰值频率,在这两个峰值频率之间存在着视觉上的陷波现象(图 12b)。

图 10 精细近地表模型 (a)微测井数据;(b)近地表速度解释结果;(c)近地表模型

图 11 建立的综合模型及其正演波场快照

图 12 模拟单炮记录(a)及其频谱(b)
2.3 不同近地表厚度浅层多次折射波对地震资料的影响

对SZ工区资料中的浅层多次折射波调查发现,并不是所有地方都有明显的浅层多次折射波,且强度也有差异,这主要由近地表结构决定。利用SZ工区中的单炮记录及近地表结构图分析不同近地表厚度与单炮中浅层多次折射波现象之间的对应关系(图 13)。图中E、F方框所示分别为无和有明显浅层多次折射影响的两个时窗,图 14为这两个时窗对应的频谱。可以看到,E处表层较薄,为4m,单炮中浅层多次折射波不发育,频谱中没有伪陷波现象。而F处表层较厚,大于20m,单炮中表现为强烈的浅层多次折射波,频谱中可以见到明显的伪陷波现象。因此低降速带厚度与浅层多次折射波的发育及伪陷波现象有关。

图 13 不同近地表厚度与浅层多次折射波的对应关系 上图为近地表厚度图,中图为近地表结构图,下图为对应的单炮记录

图 14 图 13中E(a)、F(b)两个窗口对应的频谱

为分析整个工区近地表结构对浅层多次折射波的影响,按照近地表厚度将全区分为四类:小于5m;5~15m;15~20m;大于20m。在图 13中,A区近地表厚度为17m,B区近地表厚度为16m,C区近地表厚度为9m,D区近地表厚度为22m。从单炮记录面貌看,D区浅层多次折射波能量较强,贯穿于记录的始终。图 15是A区2.5、3.5、4.0、4.5和5.0s处时窗(图 13炮集黄框所示)的频谱,可以看出,小于4.5s都存在明显的“伪陷波”现象,而大于4.5s“伪陷波”现象不再明显。因此近地表厚度为15~20m时浅层多次折射波能量较强,能影响到4.5s,也就是奥陶系目的层反射。图 16是C区1.7、2.5和3.3s处时窗(图 13炮集蓝框所示)的频谱,小于2.5s存在明显的“伪陷波”现象,大于2.5s“伪陷波”现象不明显。因此近地表厚度为5~15m时浅层多次折射波最大影响到2.5s左右。

图 15 图 13A区2.5s(a)、3.5s(b)、4.0s(c)、4.5s(d)和5.0s(e)处时窗(图 13炮集黄框所示)的频谱

图 16 图 13C区1.7s(a)、2.5s(b)和3.5s(c)处时窗(图 13炮集蓝框所示)的频谱

图 13~图 16的分析可知,近地表厚度小于5m时,单炮中的浅层多次折射波不明显;近地表厚度大于5m时,单炮中有明显的浅层多次折射波;近地表厚度大于20m时,浅层多次折射波能量最强,干扰范围最大。因此近地表厚度越大,浅层多次折射波能量越强,影响范围越大,伪陷波现象越明显。

2.4 不同近地表厚度的波场模拟

为什么近地表厚度越大,浅层多次折射波能量越强?为此,建立不同近地表厚度模型(图 17a),通过波动方程正演模拟,对模拟数据进行分析。

图 17 不同近地表厚度模型(a)及模拟的单炮记录(b) 从左到右近地表厚度依次为4、10、20、30、40和50m

由模拟炮记录可见,随着近地表模型厚度的增加,浅层多次折射波能量逐渐增强(图 17b)。图 18是近地表厚度为50m时地震波在近地表内反复震荡波场快照。首先震源激发,地震波从低降速带之下传入低降速带内,当遇到地表时,地震波被反射回来向下传播,当遇到高速层顶界面,地震波再次被反射回传,如此地震波在地表和高速层顶界面间来回震荡,形成浅层多次折射波[12-14]

图 18 近地表厚度为50m时的波场模拟快照

图 19是近地表厚度为4m时的波场快照,地震波在近地表内波长是10m,在高速层内是50m,当近地表厚度较小时,即小于1/2波长,地震波从高速层传入近地表时,地震波直接透射过去,不会在近地表内震荡,从而不会形成浅层多次折射波。

图 19 近地表厚度为4m时的波场模拟快照

当近地表厚度较大时,即大于1/2波长,地震波从高速层传入近地表时,地震波在近地表内形成多次震荡,从而形成浅层多次折射波。且近地表厚度越大,浅层多次折射波能量越强,影响范围越大。

3 陷波现象压制方法

由以上分析可知,顺北沙漠区地震资料中存在的陷波现象是由近地表浅层多次折射波能量较强导致的,严重影响了地震资料的信噪比。浅层多次折射波具有固定的视速度(约为1700m/s)和频率(主频为10Hz),在近地表内沿高速层顶界面进行传播,是典型的线性干扰[15-17],可通过频率和速度特征进行约束剔除[18-20]。线性干扰压制技术较FK滤波具有更好的压制效果,不受空间假频影响。图 20为浅层多次折射压制前、后结果对比,可见浅层多次折射波分布于整个单炮中,剔除浅层多次折射波后单炮信噪比有了较大提高,反射同相轴更加清晰连续,单炮频谱中的伪陷波现象已不明显。

图 20 多次折射波压制前、后效果对比 (a)原始记录;(b)剔除的浅层多次折射波;(c)压制后记录;(d)原始记录频谱;(e)压制后记录频谱
4 结论

顺北沙漠区地震单炮记录中的“陷波”现象,是由于近地表结构中高速层顶界面是一个较强的波阻抗界面,地震波在近地表内传播时反复震荡,在单炮记录中形成能量较强的浅层多次折射波,速度为高速层顶界面速度,主频为10Hz左右。该多次折射波贯穿于整个单炮记录中,与地下有效反射波混叠在一起,影响了深层反射波的信噪比,在频谱中表现为与地下反射波主频之间形成“两高夹一低”的形态,形成了视觉上的“陷波”现象。当近地表厚度大于5m时,单炮中有明显的浅层多次折射波,并且当近地表厚度大于20m时,浅层多次折射波能量更强,所产生的陷波现象影响到了目的层反射。由于浅层多次折射波具有固定的速度和频率,因此可以在处理中加以压制,达到消除陷波现象、提高资料信噪比的目的。

参考文献
[1]
陈学强, 段孟川, 钟海, 等. 塔里木盆地碳酸盐岩提高分辨率地震采集技术[J]. 石油地球物理勘探, 2004, 39(西部会议专刊): 50-54.
CHEN Xueqiang, DUAN Mengchuan, ZHONG Hai, et al. High-resolution seismic acquisition technology for carbonate rocks in Tarim basin[J]. Oil Geophysical Prospecting, 2004, 39(Western Conference Special Issue): 50-54.
[2]
杨平, 孙赞东, 李海银, 等. 影响岩溶缝洞体地震反射特征的关键因素分析[J]. 石油地球物理勘探, 2015, 50(3): 523-529.
YANG Ping, SUN Zandong, LI Haiyin, et al. Analysis of key factors affecting seismic reflection characteristics of karst fractures and caves[J]. Oil Geophysical Prospecting, 2015, 50(3): 523-529.
[3]
蔡希玲, 张俊桥. 复杂地表区噪声分析技术与压噪方法应用[J]. 石油地球物理勘探, 2002, 37(增刊1): 1-4.
CAI Xiling, ZHANG Junqiao. Noise analysis techno-logy and application of noise suppression method in complex surface area[J]. Oil Geophysical Prospecting, 2002, 37(S1): 1-4.
[4]
李天树, 陈宝德. 双井微测井技术在表层调查结构中的应用[J]. 石油物探, 2004, 43(5): 471-474.
LI Tianshu, CHEN Baode. Application of dual well micrologging in surface investigation structure[J]. Geophysical Prospecting for Petroleum, 2004, 43(5): 471-474. DOI:10.3969/j.issn.1000-1441.2004.05.013
[5]
张庆淮, 李玉兴. 沙漠区宽频高能噪声对地震记录的影响[J]. 石油地球物理勘探, 2008, 43(2): 142-145.
ZHANG Qinghuai, LI Yuxing. Effects of broadband high-energy noise on seismic records in desert areas[J]. Oil Geophysical Prospecting, 2008, 43(2): 142-145. DOI:10.3321/j.issn:1000-7210.2008.02.004
[6]
吕公河. 地震勘探虚反射界面的测定及其利用[J]. 石油地球物理勘探, 2002, 37(3): 296-300.
LYU Gonghe. Determination and utilization of ghost reflection interface in seismic exploration[J]. Oil Geophysical Prospecting, 2002, 37(3): 296-300.
[7]
李合群, 孟小红, 赵波, 等. 塔里木沙漠区地震数据品质与沙层Q吸收[J]. 石油地球物理勘探, 2010, 45(1): 28-34.
LI Hequn, MENG Xiaohong, ZHAO Bo, et al. Seismic data quality and sand Q absorption in Tarim desert[J]. Oil Geophysical Prospecting, 2010, 45(1): 28-34.
[8]
谭绍泉, 徐锦玺, 张庆怀, 等. 虚反射和岩性对激发井深选择的影响[J]. 油气地球物理, 2003, 1(3): 11-14.
TAN Shaoquan, XU Jinxi, ZHANG Qinghuai, et al. Influence of ghost reflection and lithology on selection of excitation well depth[J]. Oil and Gas Geophysics, 2003, 1(3): 11-14.
[9]
凌云, 高军, 张汝杰. 地表一致性变预测步长反褶积与沙丘鸣震的压制[J]. 石油地球物理勘探, 1998, 33(1): 23-31.
LING Yun, GAO Jun, ZHANG Rujie. Surface consis-tent deconvolution using variable prediction step and suppression of sand dune ringing[J]. Oil Geophysical Prospecting, 1998, 33(1): 23-31. DOI:10.3321/j.issn:1000-7210.1998.01.004
[10]
陈学强, 白文杰, 钟海. 沙漠区高信噪比地震采集方法[J]. 石油地球物理勘探, 2010, 45(3): 321-325.
CHEN Xueqiang, BAI Wenjie, ZHONG Hai. High signal-to-noise ratio seismic acquisition method in desert area[J]. Oil Geophysical Prospecting, 2010, 45(3): 321-325.
[11]
李庆忠. 走向精确勘探道路[M]. 北京: 石油工业出版社, 1994: 34-38.
[12]
龚明平, 张军华, 王延光, 等. 分方位地震勘探研究现状及进展[J]. 石油地球物理勘探, 2018, 53(3): 642-658.
GONG Mingping, ZHANG Junhua, WANG Yanguang, et al. Research status and progress of azimuth seismic exploration[J]. Oil Geophysical Prospecting, 2018, 53(3): 642-658.
[13]
赵玲芝, 谷跃民, 张建中, 等. 多信息融合的近地表速度建模技术及应用[J]. 石油地球物理勘探, 2017, 52(1): 34-41.
ZHAO Lingzhi, GU Yuemin, ZHANG Jianzhong, et al. Multi-information fusion near-surface velocity modeling technology and its application[J]. Oil Geophysical Prospecting, 2017, 52(1): 34-41.
[14]
王晓涛, 谭佳, 毛海波, 等. 强衰减条件下能量比法的修正及其应用[J]. 石油地球物理勘探, 2017, 52(1): 42-47.
WANG Xiaotao, TAN Jia, MAO Haibo, et al. Modification and application of energy ratio method under strong attenuation conditions[J]. Oil Geophysical Prospecting, 2017, 52(1): 42-47.
[15]
张华, 陈小宏, 李红星, 等. 曲波变换三维地震数据去噪技术[J]. 石油地球物理勘探, 2017, 52(2): 226-232.
ZHANG Hua, CHEN Xiaohong, LI Hongxing, et al. Seismic data denoising technology based on Curvelet transform[J]. Oil Geophysical Prospecting, 2017, 52(2): 226-232.
[16]
白英哲, 孙赞东, 周新源. Q偏移在塔里木盆地碳酸盐岩储层成像中的应用[J]. 石油地球物理勘探, 2011, 46(增刊1): 7-11.
BAI Yingzhe, SUN Zandong, ZHOU Xinyuan. Application of Q migration in carbonate reservoir imaging in Tarim basin[J]. Oil Geophysical Prospecting, 2011, 46(S1): 7-11.
[17]
何振香. 改进南疆地震勘探工作的几点意见[J]. 石油地球物理勘探, 1980, 15(增刊1): 92-94.
HE Zhenxiang. Suggestions on improving seismic exploration in southern Xinjiang[J]. Oil Geophysical Prospecting, 1980, 15(S1): 92-94.
[18]
周翼, 陈学强, 江民, 等. 塔克拉玛干沙漠区浅表层对地震波的吸收衰减作用[J]. 石油地球物理勘探, 2016, 51(2): 218-223.
ZHOU Yi, CHEN Xueqiang, JIANG Min, et al. Absorption and attenuation of seismic waves in the shal-low surface of Taklamakan desert[J]. Oil Geophysical Prospecting, 2016, 51(2): 218-223.
[19]
吴群苍. 塔克拉玛干沙漠区静校正数据的采集与整理[J]. 石油地球物理勘探, 1985, 20(2): 146-152.
WU Quncang. Collection and collation of statics data in Taklamakan desert[J]. Oil Geophysical Prospecting, 1985, 20(2): 146-152.
[20]
于承业, 周志才. 利用双井微测井资料估算近地表Q值[J]. 石油地球物理勘探, 2011, 46(1): 89-92.
YU Chengye, ZHOU Zhicai. Using dual-well micro-logging data to estimate Q value near the surface[J]. Oil Geophysical Prospecting, 2011, 46(1): 89-92.