An optimised one-way wave migration method for complexmedia with arbitrarily varying velocity based on eigen-decomposition

The classical one-way generalised screen propagator (GSP) and Fourier finite-difference (FFD) schemes have limitations in imaging large angles in complex media with substantial lateral variations in wave velocity. Some improvements to the classical one-way wave scheme have been proposed with optimised methods. However, the performance of these methods in imaging complex media remains unsatisfying. To overcome this issue, a new strategy for wavefield extrapolation based on the eigenvalue and eigenvector decomposition of the Helmholtz operator is presented herein. In this method, the square root operator is calculated after the decomposition of the Helmholtz operator at the product of the eigenvalues and eigenvectors. Then, Euler transformation is applied using the best polynomial approximation of the trigonometric function based on the infinite norm, and the propagator for one-way wave migration is calculated. According to this scheme, a one-way operator can be computed more accurately with a lower-order expansion. The imaging performance of this scheme was compared with that of the classical GSP, FFD and the recently developed full-wave-equation depth migration (FWDM) schemes. The impulse responses in media with arbitrary velocity inhomogeneity demonstrate that the proposed migration scheme performs better at large angles than the classical GSP scheme. The wavefronts calculated in the dipping and salt dome models illustrate that this scheme can provide a precise wavefield calculation. The applications of theMarmousi model further demonstrate that the proposed approach can achieve better-migrated results in imaging small-scale and complex structures, especially in media with steep-dipping faults.