High-precision elastic reverse-time migration (ERTM) imaging has always been one of the trends in the development of geophysics. However, current wavefield simulations using time-domain finite-difference (FD) approaches in ERTM have second-order temporal accuracy, resulting in travel time changes and waveform distortion in wavefield propagation with large time steps, i.e., temporal dispersion. Errors caused by the temporal dispersion can lead to erroneous imaging locations and out-of-focus diffraction events. A new ERTM and its workflow are established here using temporal and spatial high-order FD accuracy wavefields and the vector-based imaging condition. Our method computes elastic vector-based wavefields by solving a P- and S-wave decomposition form of a quasi-stress–velocity equation. An advanced finite-difference scheme is employed in the wavefield solution to achieve simulation with temporal fourth-order accuracy and spatial arbitrary even-order accuracy. The normalized dot-product imaging condition of the source and receiver P/S wavefields is then applied to generate high-quality images. The elastic wavefield simulation and ERTM numerical examples presented here reveal that the anti-dispersion workflow can improve modeling and imaging accuracy. In addition, the field data application shows that our method can achieve reasonable and reliable ERTM images. This method can integrate the most advanced imaging techniques into this computational framework to improve imaging accuracy.

Three-dimensional (3D) elastic reverse-time migration (ERTM) can image the subsurface 3D seismic structures, and it is an important tool for the Earth’s interior imaging. A common simulation kernel used in 3D ERTM is the current staggered-grid finite-difference (SGFD) method of the first-order elastic wave equation. However, the mere second-order accuracy in time of the current SGFD method can bring non-negligible time dispersion, which reduces the simulation accuracy and further leads to the distortion of the imaging results. This paper proposes a vector-based 3D ERTM using the high-order accuracy SGFD method in time to obtain high-accuracy images. This approach is a new high-resolution ERTM workflow that improves the imaging accuracy of conventional ERTM from numerical simulation. The proposed ERTM workflow is established on a quasi-stress–velocity wave equation and its vector wavefield decomposition form. Advanced SGFD schemes and their corresponding coefficients with fourth-order temporal accuracy solve the quasi-linear wave equation system. The normalized dot product imaging condition produces high-quality images using high-accuracy vector wavefields solved using the SGFD method. Through the numerical examples, we test the simulation efficiency and analyze how temporal accuracy in numerical simulations affects migration imaging quality. We include that the proposed method obtains highly accurate images.