Eliminating the redundant source effects from the cross correlation reverse time migration using a modified stabilized division
Cross correlation reverse time migration is the kernel of two way wave equation migration and inversion. However, it more or less tapers the spectrum of receiver data due to a redundant overlay of the source wavelet, whose amplitude spectrum is usually bandlimited and non flat. To circumvent this issue, there are two optional strategies: whitening the source directly, or preconditioning the seismic traces by division with the amplitude spectrum of the source in the frequency domain. In this paper, we choose the latter one because the source signature is crucial to illumination compensation and seismic inversion. To avoid division by zero, a modified stabilized division algorithm based on the Taylor expansion is developed. The modified division is easy to complete with computers and can be extend to any order. Moreover, when simulating 2 D source wavefield, the half integral effect is also considered. We will demonstrate our proposed scheme using the Sigsbee2b synthetic data and a real field data. Conventional RTM suffers from the redundant source effect due to its non flat spectrum.We precondition the seismic traces by division with source the amplitude spectrum in the frequency domain.No artificial source is employed.A modified stabilized division is considered.The half integral effect of 2 D source wavefield is considered.
Comput. Geosci.