摘要
Pre-stack reverse time migration (RTM) for rugged topography is a very useful tool for seismic imaging in the areas with rugged surface topography andcomplex subsurface structures. In this paper, we illustrate the implementation process of RTM for rugged topography. For the difficulty in dealing with the rugged free boundary condition with the finite difference method, we employ a simplified boundary condition which can avoid the abundant logical judgment. On this basis, we use the Graphic Processing Unit (GPU) algorithm to accelerate RTM and get an order of magnitude higher speedup ratio compared to the traditional CPU algoritm. The tests onsynthetic data and the 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 near surface structures and subsurface steep structures are significantly improved. The problems of imaging noise removal and massy memory demand of RTM have been stated previously and will not be discussed in this paper.