A one-way propagator is proposed for more accurately modeling wide-angle wavefields in the presence of severe lateral variations of the velocity. The method adds a higher-order correction to improve the split-step Fourier method by directly designing a cascaded operator that matches the exact phase-shift operator of a varying velocity. Using an optimization scheme, the coefficients in the cascaded operator are determined according to the local velocity distribution and the prescribed angular range of wavefield propagation. The proposed algorithm is implemented alternately in spatial and wavenumber domains using fast Fourier transforms, as in the split-step Fourier and generalized-screen methods. This algorithm can achieve higher accuracy than the generalized-screen method for wide-angle wavefields, although the same numerical scheme is used with comparable computational cost. No extra error arises for the proposed algorithm when used for 3D wave propagation, in contrast to methods that introduce an implicit finite–difference higher-order correction to the split-step Fourier method, such as the Fourier finite difference (FFD) and wide-angle screen methods. A detailed comparison of the proposed one-way propagator with the split-step Fourier, generalized-screen, and FFD methods is presented. The 2D Marmousi and 3D SEG/EAEG overthrust data sets are used to test the prestack depth-migration schemes developed based on the proposed one-way propagators.

You do not have access to this content, please speak to your institutional administrator if you feel you should have access.