Turning wave false image removal based on wavelength-dependent smoothing operator in reverse time migration
Tang Yongjie1,2, Song Guiqiao3, Liu Shaoyong1,2, Gu Hanming1,2, Yan Zhe1,2
1. Institute of Geophysics and Geomatics, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China; 2. Hubei Subsurface Multi-scale Imaging Key Laboratory, Wuhan, Hubei 430074, China; 3. Exploration & Development Department, SINOPEC, Beijing 100728, China
Abstract:The conventional reverse time migration (RTM) based on cross-correlation imaging conditions not only produces low-frequency artifacts,but also generates the false images of turning wave in the strong velocity gradient regions,and the false images cannot be removed by a simple Laplacian filter.In this paper,we introduce the imaging condition by separating up-going and down-going wavefields to analyze the causes of the false images in RTM and it does indicate that the false images of turning wave are generated by the cross-correlation of the up-going source and down-going receiver wavefields in the strong velocity gradient areas.In order to eliminate the false images and take into account the computational efficiency,we propose a false image removal method based on the wavelength-dependent smoothing operator.This method includes the wavelength effect which makes the propagation characteristics of wavefields in the strong velocity gradient areas more close to the propagation characteristics in the true velocity.Numerical experiments demonstrate that the proposed wavelength-dependent smoothing algorithm can suppress the turning wave false images and perform RTM efficiently under the classical imaging framework.
Yoon K,Marfurt K J,Starr W.Challenges in reverse-time migration.SEG Technical Program Expanded Abstracts,2004,23:1057-1060.
[2]
Zhu J,Lines L R.Comparison of Kirchhoff and reverse time migration methods with application to prestack depth imaging of complex structures.Geophysics,1998,63(4):1166-1176.
[3]
Sava P,Hill S J.Overview and classification of wavefield seismic imaging methods.The Leading Edge,2009,28(2):170-183.
[4]
杨午阳,张厚柱,撒利明等. 逆时偏移关键问题探讨. 石油地球物理勘探,2016,51(6):1251-1262.Yang Wuyang,Zhang Houzhu,Sa Liming et al.Some issues about reverse time migration.OGP,2016,51(6):1251-1262.
[5]
王鹏飞,何兵寿.基于行波分离的三维弹性波矢量场点积互相关成像条件.石油地球物理勘探,2017,52(3):477-483.Wang Pengfei,He Bingshou.Vector field dot product cross-correlation imaging based on 3D elastic wave separation.OGP,2017,52(3):477-483.
[6]
叶月明,郭庆新,孙小东等.多次波对逆时偏移构造成像的影响.石油地球物理勘探,2014,49(3):523-529.Ye Yueming,Guo Qingxin,Sun Xiaodong et al.Multiples influence on structure imaging in reverse time migration.OGP,2014,49(3):523-529.
[7]
Deeks J,Lumley D.Prism waves in seafloor canyons and their effects on seismic imaging.Geophysics,2015,80(6):S213-S222.
[8]
Weglein A B.Multiples:Signal or noise? Geophysics,2016,81(4):V283-V302.
[9]
Whitmore D N.Iterative depth migration by backward time propagation.SEG Technical Program Expanded Abstracts,1983,2:382-385.
[10]
Baysal E,Kosloff D D,Sherwood J W C.Reverse time migration.Geophysics,1983,48(11):1514-1524.
[11]
McMechan G A.Migration by extrapolation of time-dependent boundary values.Geophysical Prospecting,1983,31(3):413-420.
[12]
Fricke J R.Reverse-time migration in parallel:A tutorial.Geophysics,1988,53(9):1143-1150.
[13]
Harris C E,McMechan G A.Using downward continu-ation to reduce memory requirements in reverse-time migration.Geophysics,1992,57(6):848-853.
[14]
Symes W.Reverse time migration with optimal check-pointing.Geophysics,2007,72(5):SM213-SM221.
[15]
刘红伟,李博,刘洪等.地震叠前逆时偏移高阶有限差分算法及GPU实现.地球物理学报,2010,53(7):1725-1733.Liu Hongwei,Li Bo,Liu Hong et al.The algorithm of high order finite difference prestack reverse time migration and GPU implementation.Chinese Journal of Geophysics,2010,53(7):1725-1733.
[16]
Anderson J E,Tan L,Wang D.Time-reversal check pointing methods for RTM and FWI.Geophysics,2012,77(4):S93-S103.
[17]
张慧,蔡其新,秦广胜等.基于GPU并行加速的逆时偏移成像方法.石油地球物理勘探,2013,48(5):707-710.Zhang Hui,Cai Qixin,Qin Gunagsheng et al.GPU-accelerated reverse time migration and its application.OGP,2013,48(5):707-710.
[18]
唐祥功,匡斌,杜继修等.多GPU协同三维叠前逆时偏移方法研究与应用.石油地球物理勘探,2013,48(6):910-914.Tang Xianggong,Kuang Bin,Du Jixiu et al.3D pre-stack reverse time migration based on multiple-GPU collaboration.OGP,2013,48(6):910-914.
[19]
Clapp R G.Seismic processing and the computer revolution.SEG Technical Program Expanded Abstracts,2015,34:4832-4837.
[20]
Yang P,Brossier R,Virieux J.Wavefield reconstruction by interpolating significantly decimated boundaries.Geophysics,2016,81(5):T197-T209.
[21]
Loewenthal D,Stoffa P L,Faria E L.Suppressing the unwanted reflections of the full wave equation.Geophysics,1987,52(2):1007-1012.
[22]
Mulder W,Plessix R.One-way and two-way wave-equation migration.SEG Technical Program Expanded Abstracts,2003,22:881-884.
[23]
Valenciano A,Biondi B.2D deconvolution imaging condition for shot-profile migration.SEG Technical Program Expanded Abstracts,2003,22:1059-1062.
Yoon K,Marfurt K J.Reverse-time migration using the poynting vector.Exploration Geophysics,2006,37(1):102-107.
[26]
Guitton A,Kaelin B,Biondi B.Least-square attenuation of reverse-time migration artifacts.SEG Technical Program Expanded Abstracts,2006,25:2348-2352.
[27]
Liu F,Zhang G Q,Morton S A.Reverse-time migration using one-way wavefield imaging condition.SEG Technical Program Expanded Abstracts,2007,26:2170-2174.
[28]
Zhang Y,Sun J.Practical issues of reverse time mi-gration:true amplitude gathers,noise removal and harmonic-source encoding.ASEG Extended Abstracts,2009,1-5.
[29]
Liu F,Zhang G Q,Morton S A et al.An effective imaging condition for reverse-time migration using wavefield decomposition.Geophysics,2011,76(1):S29-S39.
[30]
Rocha D,Tanushev N,Sava P.Acoustic wavefield imaging using the energy norm.Geophysics,2016,81(4):S151-S163.
[31]
Fei T W,Luo Y,Yang J et al.Removing false images in reverse time migration:The concept of de-primary.Geophysics,2015,80(6):S237-S244.