An efficient local operator-based Q-compensated reverse time migration algorithm with multistage optimization

Abstract

Most existing [Formula: see text]-compensated reverse time migration ([Formula: see text]-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 [Formula: see text]-RTM algorithm — the multistage optimized [Formula: see text]-RTM method. This [Formula: see text]-RTM algorithm uses a finite-difference method to compensate the amplitude and the phase simultaneously by uniquely combining two techniques: (1) a negative [Formula: see text] 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 [Formula: see text]-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 [Formula: see text]-RTM methods, this [Formula: see text]-RTM approach offers nearly identical imaging quality. Based on local numerical differential operators, this [Formula: see text]-RTM method is very suitable for parallel computing and graphic processing unit implementation, an important feature for large 3D seismic surveys.