Most existing -compensated reverse time migration (-RTM) algorithms are based on pseudospectral methods. Because of the global nature of pseudospectral operators, these methods are not ideal for efficient parallelization, implying that they may suffer from high computational cost and inefficient memory usage for large-scale industrial problems. In this work, we reported a novel -RTM algorithm — the multistage optimized -RTM method. This -RTM algorithm uses a finite-difference method to compensate the amplitude and the phase simultaneously by uniquely combining two techniques: (1) a negative method for amplitude compensation and (2) a multistage dispersion optimization technique for phase correction. To prevent high-frequency noise from growing exponentially and ruining the imaging results, we apply a finite impulse response low-pass filter using the Kaiser window. The theoretical analyses and numerical experiments demonstrate that this -RTM algorithm precisely recovers the decayed amplitude and corrects the distorted phase caused by seismic attenuation effects, and hence produces higher resolution subsurface images with the correct structural depth information. This new method performs best in the frequency range of 10–70 Hz. Compared with pseudospectral -RTM methods, this -RTM approach offers nearly identical imaging quality. Based on local numerical differential operators, this -RTM method is very suitable for parallel computing and graphic processing unit implementation, an important feature for large 3D seismic surveys.