石油地球物理勘探  2021, Vol. 56 Issue (6): 1229-1235  DOI: 10.13810/j.cnki.issn.1000-7210.2021.06.004
0
文章快速检索     高级检索

引用本文 

林娟, 蒋立, 潘龙, 冷雪梅, 张欣吉, 游伟. 复杂近地表多次波分析及压制方法探讨——准噶尔盆地腹部地震资料处理实例. 石油地球物理勘探, 2021, 56(6): 1229-1235. DOI: 10.13810/j.cnki.issn.1000-7210.2021.06.004.
LIN Juan, JIANG Li, PAN Long, LENG Xuemei, ZHANG Xinji, YOU Wei. Analysis of complex near-surface multiple waves and discussion on suppression methods: a case of seismic data in the hinterland of the Junggar Basin. Oil Geophysical Prospecting, 2021, 56(6): 1229-1235. DOI: 10.13810/j.cnki.issn.1000-7210.2021.06.004.

作者简介

林娟   高级工程师, 1982年生; 2005年获成都理工大学信息工程专业工学学士学位, 2008年获该校应用地球物理专业工学硕士学位; 现就职于中国石油新疆油田公司勘探开发研究院石油地球物理研究所, 主要从事地震勘探数据处理及其方法研究

林娟, 新疆维吾尔自治区乌鲁木齐市北京北路397号新疆油田公司勘探开发研究院地球物理所, 830013。Email: linj0103@petrochina.com.cn

文章历史

本文于2020年12月27日收到,最终修改稿于2021年9月22日收到
复杂近地表多次波分析及压制方法探讨——准噶尔盆地腹部地震资料处理实例
林娟 , 蒋立 , 潘龙 , 冷雪梅 , 张欣吉 , 游伟     
中国石油新疆油田分公司勘探开发研究院地球物理研究所, 新疆乌鲁木齐 830013
摘要:准噶尔盆地腹部地区地表主要是沙漠区和公益林区,从东西向剖面分析得知这两种地表对应的波组特征、频率及相位特征差异明显,造成了地震剖面上明显的“阴阳脸”现象。为此,针对不同地表的单炮记录,分析初至频谱存在“陷波”特征的原因;结合表层结构及正演模型论证分析,认为复杂近地表多次波是产生“阴阳脸”现象的主要因素;明确了陷波强弱及产生的具体位置,进而提出基于旅行时差异及自相关第一波谷时间差异的空变预测步长反褶积技术,对该区近地表多次波压制效果明显,改善和消除了“阴阳脸”现象。该方法在准噶尔盆地腹部其他区块中推广应用,均见到良好效果。
关键词复杂近地表    多次波    陷波    变步长反褶积    准噶尔盆地    
Analysis of complex near-surface multiple waves and discussion on suppression methods: a case of seismic data in the hinterland of the Junggar Basin
LIN Juan , JIANG Li , PAN Long , LENG Xuemei , ZHANG Xinji , YOU Wei     
Institute of Geophysics, Research Institute of Exploration and Development, Xinjiang Oilfield Company, PetroChina, Urumqi, Xinjiang 830013, China
Abstract: The surface in the hinterland of the Junggar Basin is mainly desert and ecological forest. From the analysis of the east-west section, it is known that the wave group characteristics and the frequency and phase characteristics corresponding to the two types of surface are greatly different, which results in the obvious "Yin-yang face" phenomenon on the seismic section. Therefore, considering the single shot records on different types of surface, we analyzed the reasons for the "trapping" characteristics in the first-break frequency spectrum. Combined with the demonstration and analysis of the surface structure and forward mo-del, we believe that complex near-surface multiple waves are the main factor causing the "Yin-yang face" phenomenon. Moreover, the strength and specific location of trapping are clarified, and thus this paper proposes a predictive variant-step deconvolution technology based on the difference of tra-vel time and of the time of the first trough in autocorrelation, which has an obvious suppression effect on near-surface multiple waves in this area and improves and eliminates the phenomenon of "Yin-yang face". This method has been widely used and applied in other blocks in the hinterland of the Junggar Basin, and good results have been achieved.
Keywords: complex near-surface    multiple waves    "trapping"    variant-step deconvolution    Junggar Basin    
1 研究背景

中国新疆准噶尔盆地腹部前哨地区的地震剖面存在公益林区与沙漠区波组差异大的特征(图 1),通过原始单炮初至的频谱分析,可看出该区频谱普遍存在“陷波”特征。针对此现象,从原始单炮对应的不同地表条件的激发及接收因素、表层结构、正演子波模型分析等多角度进行分析,发现此“陷波”现象是由近地表多次波引起,当多次波与有效反射波重叠时,地震资料频谱就会出现明显的“两凸一凹”,即“陷波”现象。

图 1 公益林及沙漠区初叠剖面

多次波压制是地震数据处理的一个难题,现今主要应用抛物线拉东变换、SRME等多次波去噪方法[4-6],同时由于多次波在时间上呈现一定规律性,因此人们就利用预测反褶积去除多次波。预测反褶积步长可通过自相关得到,同时运行速度快,目前已经成为压制多次波的一种常用方法[7]

针对近地表多次波引起的“陷波”,结合子波模型正演模拟分析,本文研究了基于旅行时的空变预测步长反褶积方法对近地表引起的多次波进行压制;结合压制后剖面自相关分析结果,串联应用自相关第一波谷时间空变步长反褶积,更彻底地压制残存的近地表多次波[7],消除或改善了“陷波”特征及剖面上“阴阳脸”现象,提升了资料品质。

2 方法原理 2.1 近地表多次波产生机理 2.1.1 不同震源类型分析

准噶尔盆地腹部前哨1井三维及相邻的前哨2井三维地震均使用可控震源激发,其单炮记录频谱均存在“两凸一凹”的“陷波”特征。为了进一步分析此现象与震源类型是否相关,在该区相同位置分别选取一条2007年井炮采集的二维测线及2017年可控震源采集的二维测线。从单炮的相同初至范围的频谱分析(图 2)可见,井炮激发记录中看不出明显的“陷波”特征。究其原因是井炮的激发点处于高速层以下,地震波在传播过程中打穿高速层,因此即使存在近地表多次波也不明显。

图 2 井炮(蓝色)和可控震源(红色)激发单炮的频谱
2.1.2 激发、接收因素分析

由于该区主要为公益林和沙漠两种类型的地表,因此选取激发位置分别在沙漠、公益林和过渡带等三类典型地区的多个单炮记录,进而结合接收位置分别处于公益林和沙漠区的初至波频谱分析,发现所有位置均存在“陷波”特征(图 3~图 5)。由此可知:无论激发位置和接收位置处于公益林或沙漠区,均存在“陷波”现象;对频谱细致对比可看出,公益林区的“陷波”特征较沙漠区更加明显。

图 3 公益林区激发典型单炮及不同位置点(红框)的频谱

图 4 沙漠区激发典型单炮及不同位置点(红框)的频谱

图 5 过渡区激发典型单炮及不同位置点(红框)的频谱

结合本区Inline方向(图 6a)表层模型(图 6b)和不同位置微测井时距曲线分析(图 7),可见目标区自公益林区向沙漠区(从西向东)方向表层厚度逐渐增厚,公益林区表层厚度约6~8m,沙漠区表层厚度近40m。根据沙层吸收衰减一般规律,随着沙漠区表层厚度增大,吸收作用增强,这也是沙漠区接收的波场“陷波”现象较公益林区弱的原因之一。

图 6 目标区卫星图(a)及表层模型(b)

图 7 公益林区(a)、过渡区(b)及沙漠区(c)三种典型地表的微测井时距曲线 下标“0”代表低速层,下标“1、2、3”对应不同降速层;速度超过1800m/s即认为是高速层;VaH分别为地表层的平均速度和总厚度
2.2 理论子波模型分析

为了进一步分析“陷波”产生的原因及具体位置,首先设计了7个不同反射系数模型[8-10],有效波子波振幅为1.0,多次波反射系数分别为0.1、0.2、0.3、0.4、0.5、0.6、0.7,时间间隔为20ms,从这7个反射系数模型的频谱(图 8)可知:多次波反射系数越大,其“陷波”作用越强。又针对0.5反射系数分别采用10、20和30ms间隔构建三个反射系数模型,从该三模型频谱(图 9)可见,时间间隔越小,“陷波”频带越宽。

图 8 不同反射系数滤波频谱

图 9 10ms(a)、20ms(b)、30ms(c)旅行时对应的子波频谱

图 10a是用30Hz雷克子波与时间间隔为30ms的反射系数模型褶积产生的合成记录,该模型上、下层反射系数分别是1.0和-0.5,并得到相应的频谱分析结果(图 10c);图 10b是对该合成记录应用30ms时间步长的预测反褶积处理后记录,可见多次波旁瓣得到较好的压制,同时从其频谱分析图(图 10d)中也可看出“陷波”现象基本消除。

图 10 10ms(a)、20ms(b)、30ms(c)旅行时对应的子波频谱 预测反褶积前、后的子波(a、b)及频谱(c、d)
2.3 近地表多次波压制方法

前已述及,旅行时及反射系数对“陷波”位置均存在影响。因此,本文提出两种方法求取每个炮检点的预测步长Δtij。第一种是通过一次波及多次波在近地表传输过程中旅行时时差求取预测步长,第二种是通过炮检点自相关函数,结合一次波与多次波的时间差确定每个炮检点的预测步长,然后基于最小平方理论做预测反褶积处理[11-13]

$ \varepsilon = \sum\limits_{\tau = 1}^L {{{\left[ {x\left( {t - \tau } \right)*h\left( \tau \right) - x\left( {t + \Delta t} \right)} \right]}^2}} $

式中:x(t)为地震信号;h(τ)为反褶积算子;Δt为预测步长;ε为误差;L表示窗口内样点(τ)总数。

2.3.1 基于旅行时时差的变步长求取方法

根据一次波(图 11a)、多次波(图 11b)在地下介质的传播射线路径,假设存在近地表多次波,那么炮点在地表激发后,接收点既能接收到一次波能量,也能接收到多次波能量,此处仅考虑能量最大的一阶多次波影响[9-10],较复杂的高阶多次波暂不讨论。

图 11 一次波(a)和多次波(b)的传播路径(过程)示意图

设定地表激发能量是1.0,地表与高速层顶界之间传播的一阶多次波传播时间为双程低速带旅行时,假设高速层底界反射系数为0.5,地表反射系数为-1.0,正极性的一阶多次波经过高速层底界反射回0.5倍能量,再经过地表反射转换为负极性。另外,在水平层状介质假设下,零炮检距处有效波与一阶多次波的时差为低速带时深的两倍,这也是预测反褶积应选取的预测步长。

2.3.2 基于自相关剖面的步长求取方法

与旅行时求取步长方法类似,根据波传播规律进行模拟,当存在近地表多次波时,从单道自相关剖面示意图(图 12)可看出,地震波向下传播过程中,设定反射系数为1.0的波峰对应的是一次波,而反射系数为-0.5位置对应的即是多次波。依据该结果,结合前文试验的预测步长与反射系数时差的关系,本文将一次波到多次波的时间差作为预测步长,故此次从自相关剖面中提取第一波谷的时间T1作为预测步长进行反褶积处理,进而消除目标区近地表多次波的影响[14-15]

图 12 单道自相关剖面示意图
3 实际数据应用

按照上述两种变步长求取方法分别得到每一个炮检点对应的预测步长,针对前哨2井区三维实际资料进行两种方法结合的变步长预测反褶积处理,对比处理前、后的单炮(图 13)及频谱(图 14)可知,近地表引起的多次波得到较好压制,频谱图上能看出“陷波”现象也得到较好改善,同时结合变步长预测反褶积前、后自相关剖面(图 15)可见旁瓣压制效果较好。

图 13 变步长反褶积前(a)、后(b)单炮

图 14 变步长反褶积前(a)、后(b)频谱分析图

图 15 变步长反褶积前(a)、后(b)的自相关

图 16是经过两次变预测步长反褶积消除近地表多次波前后,应用相同切除、速度和静校正条件下的叠加剖面。从处理后的结果来看,剖面上公益林区与沙漠区近地表多次波引起的波组差异大的现象大体上被消除,剖面整体的一致性趋好,公益林区频谱中的“陷波”现象也得到明显改善(图 17)。

图 16 变步长反褶积前(a)、后(b)叠加剖面对比

图 17 变步长反褶积前、后不同地表区对应频谱 (a)反褶积前公益林区;(b)反褶积前沙漠区;(c)反褶积后公益林区;(d)反褶积后沙漠区
4 结论与认识

准噶尔盆地腹部前哨地区地震剖面上公益林区与沙漠区波组特征差异大的现象一直普遍存在。本文首先系统分析了波组特征不一致现象的产生原因,认为主要缘于近地表结构中高速层顶界面是一较强的波阻抗界面,且该界面产生近地表多次波;进而结合理论模型试验,提出基于旅行时及自相关求取步长的变预测步长反褶积方法,压制该近地表多次波;通过理论与实际地震资料的处理,结果表明所提方法能有效压制这种由近地表引起的多次波,使地震剖面的一致性明显改善。

通过上述研究,笔者初步认为准噶尔盆地除了前哨地区,在其他地区(尤其是地表存在稳定的高速层顶界面且下伏地层的速度是上覆地层的2倍以上)也有可能产生近地表多次波,只是有些地区这种现象明显,有些地区不明显,这与低降速层的厚度、速度、反射系数有关。

参考文献
[1]
凌云, 高军, 张汝杰. 地表一致性变预测步长反褶积与沙丘鸣震的压制[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
[2]
吕公河. 地震勘探虚反射界面的测定及其利用[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.
[3]
刘斌, 宋智强, 王超, 等. 顺北沙漠区地震资料陷波现象原因分析与消除[J]. 石油地球物理勘探, 2019, 54(6): 1237-1245.
LIU Bin, SONG Zhiqiang, WANG Chao, et al. Trapping-wave phenomenon elimination on seismic data in Shunbei desert area[J]. Oil Geophysical Prospecting, 2019, 54(6): 1237-1245.
[4]
刘振, 张军华, 韩双, 等. 与界面有关的层间多次波预测及压制[J]. 石油地球物理勘探, 2014, 49(1): 61-67.
LIU Zhen, ZHANG Junhua, HAN Shuang, et al. La-yer related interbed multiple elimination[J]. Oil Geophysical Prospecting, 2014, 49(1): 61-67.
[5]
叶月明, 姚根顺, 赵昌垒, 等. 利用地震干涉法衰减海底相关层间多次波[J]. 石油地球物理勘探, 2015, 50(2): 225-231.
YE Yueming, YAO Genshun, ZHAO Changlei, et al. Sea bottom peg-leg multiple suppression based on seismic interferometry[J]. Oil Geophysical Prospecting, 2015, 50(2): 225-231.
[6]
谭军, 宋鹏, 李金山, 等. 基于同相轴追踪的三维地震资料多次波压制方法[J]. 石油地球物理勘探, 2017, 52(5): 894-905.
TAN Jun, SONG Peng, LI Jinshan, et al. 3D multiple suppression based on event tracking[J]. Oil Geophy-sical Prospecting, 2017, 52(5): 894-905.
[7]
张军华, 缪彦舒, 郑旭刚, 等. 预测反褶积去多次波几个理论问题探讨[J]. 物探化探计算技术, 2009, 31(1): 6-10.
ZHANG Junhua, MIAO Yanshu, ZHANG Xugang, et al. Discussion of several theoretical questions to remove seismic multiples using predictive deconvolution[J]. Computing Technologies for Geophysical and Geochemical Exploration, 2009, 31(1): 6-10. DOI:10.3969/j.issn.1001-1749.2009.01.003
[8]
黄饶, 陈小宏, 李景叶. 层间多次波模拟与影响分析[J]. 地球物理学进展, 2009, 24(3): 981-987.
HUANG Rao, CHEN Xiaohong, LI Jingye. Simulation and effect analysis of interbed multiples[J]. Progress in Geophysics, 2009, 24(3): 981-987. DOI:10.3969/j.issn.1004-2903.2009.03.022
[9]
吴静, 吴志强, 胡天跃, 等. 基于构建虚同向轴压制地震层间多次波[J]. 地球物理学报, 2013, 56(3): 985-994.
WU Jing, WU Zhiqiang, HU Tianyue, et al. Seismic internal multiple attenuation based on constructing virtual events[J]. Chinese Journal of Geophysics, 2013, 56(3): 985-994.
[10]
崔永福, 刘嘉辉, 陈猛, 等. 虚同向轴方法及其在陆上地震层间多次波压制中的应用[J]. 石油地球物理勘探, 2019, 54(6): 1228-1236.
CUI Yongfu, LIU Jiahui, CHEN Meng, et al. Land seismic peg-leg multiple attenuation with the virtual event method[J]. Oil Geophysical Prospecting, 2019, 54(6): 1228-1236.
[11]
Weglein A B. Multiple attenuation: an overview of recent advances and the road ahead[J]. The Leading Edge, 1999, 18(1): 40-44. DOI:10.1190/1.1438150
[12]
Liu F Q, Sen M K, Stoffa P L. Dip selective 2-D multiple attenuation in the plane wave domain[J]. Geophysics, 2000, 65(1): 264-274. DOI:10.1190/1.1444717
[13]
宋家文, Verschuur D J, 陈小宏. 多次波压制的研究现状与进展[J]. 地球物理学进展, 2014, 29(1): 240-247.
SONG Jiawen, Verschuur D J, CHEN Xiaohong. Research status and progess in multiple elimination[J]. Progress in Geophysics, 2014, 29(1): 240-247.
[14]
Kelamis P G, Verschuur D J. Surface-related multiple elimination on land seismic data: Strategies via case studies[J]. Geophysics, 2000, 65(3): 719-734. DOI:10.1190/1.1444771
[15]
Herrmann F J, Wang D, Verschuur D J. Adaptive curvelet-domain primary-multiple separation[J]. Geophysics, 2008, 73(3): A17-A21. DOI:10.1190/1.2904986