Generating Azimuth–Reflection Angle Gathers From Reverse Time Migration Using the High-Dimensional Local Phase Space Approximation of Seismic Wavefields

Amplitude-preserving angle gathers are ideal inputs for seismic prestack inversion. However, due to the limitation of computational efficiency, generating subsurface azimuth–reflection angle gathers from 3-D seismic imaging is still a very difficult task. In this article, we propose a new method to generate azimuth–reflection angle gathers from 3-D reverse time migration (RTM). The proposed method approximately reconstructs the source wavefield using high-dimensional wavelets and the excitation information. After using directional vectors to calculate the subsurface observation angles and applying the cross correlation imaging condition, we can generate azimuth–reflection angle gathers by angle binning. Without storing source wavefields or reconstructing source wavefields using boundary conditions, the proposed method has high computational efficiency. Numerical experiments on a synthetic model and a real marine seismic dataset demonstrate that compared with the excitation amplitude imaging condition, the proposed method can generate azimuth–reflection angle gathers with continuous complete events and high signal-to-noise ratio. The image quality and resolution of angle gathers are significantly improved. At the same time, the computational complexity does not increase much.