Abstract：Reverse time migration (RTM) is a pre-stack depth migration technique commonly used for imaging complex structures.In conventional RTM,high-order finite difference and post-imaging filter-ing algorithms are usually used，but spatial dispersion will be caused by large computing girds and effective low-frequency information will be lost when using post-imaging filtering to remove low-frequency noises.A GPU-based Fourier finite-difference RTM method is developed.It uses the Fourier finite-difference algorithm to calculate wavefield propagation,and explicit wavefield decomposition and cross-correlation imaging are directly performed in time-wavenumber domain du-ring wavefield extrapolation.In addition,the GPU-based algorithm can significantly improve the calculation speed of Fourier finite difference RTM.The application in field data shows that the method can protect and utilize the high and low frequency information effectively，so the imaging resolution can be improved apparently.
Yoon K,Marfurt K J,Starr W.Challenges in reverse-time migration[C].SEG Technical Program Expanded Abstracts,2004,23:1057-1060.
[2]
Symes W W.Reverse time migration with optimal check pointing[J].Geophysics,2007,72(5):SM213-SM221.
[3]
Zhang Y,Sun J,Gray S.Reverse-time migration:amplitude and implementation issues[C].SEG Technical Program Expanded Abstracts,2007,26:2145-2149.
[4]
Zhang Y,Sun J.Practical issues in reverse time migration:true amplitude gathers,noise removal and harmonic source encoding[J].First Break,2009,26(1):29-35.
[5]
刘定进,杨勤勇,方伍宝,等.叠前逆时深度偏移成像的实现与应用[J].石油物探,2011,50(6):545-549.LIU Dingjin,YANG Qinyong,FANG Wubao,et al.Realization and practices of prestack reverse-time depth migration[J].Geophysical Prospecting for Petroleum,2011,50(6):545-549.
[6]
王延光,匡斌.起伏地表叠前逆时深度偏移与并行实现[J].石油地球物理勘探,2012,47(2):266-271.WANG Yanguang,KUANG Bin.Prestack reverse-time depth migration on rugged topography and parallel computation realization[J].Oil Geophysical Prospecting,2012,47(2):266-271.
[7]
孔祥宁,张慧宇,刘守伟,等.海量地震数据叠前逆时偏移的多GPU联合并行计算策略[J].石油物探,2013,52(5):288-293.KONG Xiangning,ZHANG Huiyu,LIU Shouwei,et al.Multi-GPUs parallel computational stratery of prestack reverse-time migration for mass seismic data[J].Geophysical Prospecting for Petroleum,2013,52(5):288-293.
[8]
Guitton A,Kaelin B,Biondi B.Least-squares attenuation of reverse-time-migration artifacts[J].Geophy-sics,2007,72(1):S19-S23.
[9]
陈康,吴国忱.逆时偏移拉普拉斯算子滤波改进算法[J].石油地球物理勘探,2012,47(2):249-255.CHEN Kang,WU Guochen.An improved filter algorithm based on Laplace operator in reverse-time migration[J].Oil Geophysical Prospecting,2012,47(2):249-255.
[10]
闫贫,何樵登,王晓春.波动方程差分法正演的频散与分辨率[J].长春地质学院学报,1992,22(2):223-228.YAN Pin,HE Qiaodeng,WANG Xiaochun.Frequency-dispersion and resolution in wave motion forward modelling by differencing method[J].Journal of Changchun University of Earth Sciences,1992,22(2):223-228.
[11]
董良国,李培明.地震波传播数值模拟中的频散问题[J].天然气工业,2004,24(6):53-56.DONG Liangguo,LI Peiming.Dispersive problem in seismic wave propagation numerical modeling[J].Natural Gas Industry,2004,24(6):53-56.
[12]
吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[J].地球物理学进展,2005,20(1):58-65.WU Guochen,WANG Huazhong.Analysis of numerical dispersion in wave-field simulation[J].Progress in Geophysics,2005,20(1):58-65.
[13]
WuW,LinesLR,LuH.Analysisofhigher-order,finite-difference schemes in 3D reverse-time migration[J].Geophysics,1996,61(3):845-856.
[14]
Liu Y,Sen M K.A new time-space domain high-order finite difference method for the acoustic wave equation[J].Journal of Computational Physics,2009,228(23):8779-8806.
[15]
Gazdag J.Modeling of the acoustic wave equation with transform methods[J].Geophysics,1981,46(6):854-859.
[16]
Kosloff D,Baysal E.Forward modeling by the Fourier method[J].Geophysics,1982,47(7):1402-1412.
[17]
Fornberg B.The pseudo spectral method:Compari-sons with finite difference for the elastic wave equation[J].Geophysics,1987,52(4):483-501.
[18]
Fornberg B.The pseudo spectral method:Accurate representation of interfaces in elastic wave calculations[J].Geophysics,1988,53(5):625-637.
[19]
Chen J B.Modeling the scalar wave equation with Nyström methods[J].Geophysics,2006,71(5):T151-T158.
[20]
Carcione J M.A generalization of the Fourier pseudo spectral method[J].Geophysics,2010,75(6):A53-A56.
[21]
Etgen J T,Brandsberg-Dahl S.The pseudo-analytical method:Application of pseudo-Laplacians to acoustic and acoustic anisotropic wave propagation[C].SEG Technical Program Expanded Abstracts,2009,28:2552-2556.
[22]
Song X L,Fomel S.Fourier finite-difference wave propa-gation[J].Geophysics,2011,76(5):T123-T129.
[23]
Fei T W,Luo Y,Schuster G T.De-blending reverse-time migration[C].SEG Technical Program Expanded Abstracts,2010,29:3130-3134.
[24]
Fei T W,Luo Y and Qin F.An endemic problem in reverse-time migration[C].SEG Technical Program Expanded Abstracts,2014,33:3811-3815.
[25]
Liu F Q,Zhang G Q,Morton S A,et al.An effective imaging condition for reverse-time migration using wavefield decomposition[J].Geophysics,2011,76(1):S29-S39.
[26]
胡江涛,王华忠.基于解析时间波场外推与波场分解的逆时偏移方法研究[J].地球物理学报,2015,58(8):2886-2895.HU Jiangtao,WANG Huazhong.Reverse time migration using analytical time wavefield extrapolation and separation[J].Chinese Journal of Geophysics,2015,58(8):2886-2895.
[27]
Shen P,Albertin U.Up-down separation using Hilbert transformed source for causal imaging condition[C].SEG Technical Program Expanded Abstracts,2015,34:4175-4179.
[28]
Soubaras R,Zhang Y.Two-step explicit marching method for reverse time migration[C].SEG Technical Program Expanded Abstracts,2008,27:2272-2276.
[29]
Zhang Y,Zhang G Q.One-step extrapolation method for reverse time migration[J].Geophysics,2009,74(4):A29-A33.
[30]
刘文卿,王宇超,雍学善,等.基于GPU/CPU叠前逆时偏移研究及应用[J].石油地球物理勘探,2012,47(5):712-716.LIU Wenqing,WANG Yuchao,YONG Xueshan,et al.Prestack reverse time migration on GPU/CPU co-parallel computation[J].Oil Geophysical Prospecting,2012,47(5):712-716.
[31]
赵虎,武泗海,尹成,等.基于OpenACC编程模型的逆时偏移多级并行的设计与优化[J].石油地球物理勘探,2018,53(6):1307-1313.ZHAO Hu,WU Sihai,YIN Cheng,et al.Multi-level parallel design and optimization for reverse time migration based on OpenACC programming model[J].Oil Geophysical Prospecting,2018,53(6):1307-1313.
[32]
吴吉忠,左虎.叠前衰减补偿时间偏移及GPU实现[J].石油地球物理勘探,2019,54(1):84-92.WU Jizhong,ZUO Hu.Attenuation compensation in prestack time migration and its GPU implementation[J].Oil Geophysical Prospecting,2019,54(1):84-92.
