Pre‐Stack Reverse Time Migration for Rugged Topography and GPU Acceleration Technology

Pre-stack reverse time migration (RTM) for rugged topography is a very useful tool for seismic imaging in the areas with strong relief on the surface and complex subsurface structures. In this paper, we illustrate the implementation process of RTM for rugged topography. To deal with the difficulty of dealing with the rugged free boundary condition with the finite difference method, we employ a simplified boundary condition which could avoid the abundant logical judgment. On this basis, we use the Graphic Processing Unit (GPU) architecture to accelerate RTM and get an order of magnitude higher speedup ratio compared to the traditional CPU architecture. The tests on synthetic data examples and comparison with the pre-stack one-way wave method for rugged topography prove that RTM does not have the imaging dip limit, and the imaging results for the areas with rugged surface topography and subsurface steep structures are significantly improved. The problems of imaging noise removal and massy memory demand of RTM have been stated in our previous paper and will not be discussed here.

[1]  John R. Berryhill,et al.  Wave-equation datuming , 1979 .

[2]  J. R. Berryhill.  Submarine canyons: Velocity replacement by wave-equation datuming before stack , 1986 .

[3]  Huazhong Wang,et al.  Wave Equation Datuming From Irregular Surfaces Using Finite Difference Scheme , 1999 .

[4]  G. Margrave,et al.  Prestack depth imaging from topography with a Fourier method , 2001 .

[5]  Craig J. Beasley,et al.  The zero‐velocity layer: Migration from irregular surfaces , 1992 .

[6]  Xiong Xiao-jun,et al.  One-way wave equation seismic prestack forward modeling with irregular surfaces , 2006 .

[7]  J. Faye,et al.  Wave field datuming for land prestack migration , 1987 .

[8]  William A. Schneider,et al.  Wave-equation velocity replacement of the low-velocity layer for overthrust-belt data , 1995 .

[9]  J. R. Berryhill.  Wave‐equation datuming before stack , 1984 .

[10]  G. Margrave,et al.  Downward Extrapolation From Topography With a Laterally Variable Depth Step , 2000 .

[11]  D. Bevc.  Flooding the topography: Wave-equation datuming of land data with rugged acquisition topography , 1997 .

[12]  Paul L. Stoffa,et al.  Split-Step Fourier Migration , 1990 .

[13]  M. Reshef.  Depth migration from irregular surfaces with depth extrapolation methods , 1991 .

[14]  Peyman P. Moghaddam,et al.  Industrial-Scale Reverse Time Migration On GPU Hardware , 2009 .

[15]  Oz Yilmaz,et al.  Prestack layer replacement , 1986 .

[16]  H. D. Geiger,et al.  Equivalent Offset Prestack Migration For Rugged Topography , 1996 .

[17]  J. W. Wiggins.  Kirchhoff integral extrapolation and migration of nonplanar data , 1984 .

[18]  Nick Ellis,et al.  Recursive Implementation of Redatuming, Imaging, And Layer Replacement For Irregular Topography , 1989 .

[19]  Paulius Micikevicius,et al.  3D finite difference computation on GPUs using CUDA , 2009, GPGPU-2.

[20]  Paul Sava,et al.  Wave-equation Migration From Topography , 2005 .

[21]  George A. McMechan,et al.  Prestack processing of land data with complex topography , 1995 .

[22]  Hong Liu,et al.  A Method of Using GPU to Accelerate Seismic Pre‐Stack Time Migration , 2009 .

[23]  Overthrust imaging with tomo‐datuming: A case study , 1998 .