# Elastic Reverse Time Migration for Weakly Illuminated Structure
## 1 Introduction  
Elastic reverse time migration (ERTM) possesses the ability to reposition the multi-component seismic data into underground structure information, thus providing a scientific basis for further seismic interpretation and hydrocarbon development [1; 2]. Particularly, subsalt has become a significant hydrocarbon exploration target due to the appearance of salt beds as regional overburden rock [3]. However, owing to the shielding effect generated by the complex overburden rock on seismic signals, the imaging results of the subsurface geological structure are often unsatisfactorily [4; 5]. Therefore, it is essential to apply the elastic migration technique for insufficient illumination areas. The wide-azimuth measurement, illumination compensation and acquisition aperture correction approaches have been developed to enhance the imaging quality of shadow regions [6; 7; 8]. In addition, the staining algorithm can highlight the wavefields associated with weakly illuminated structures among the full wavefields [9; 10; 11], which is a powerful means to realize target-oriented imaging within the ERTM framework.  
ERTM images produced by multi-component seismic data always involve crosstalks due to the reason that P- and S-waves are concurrently excited when the waves propagate to a reflecting interface [12]. To suppress the crosstalk artifacts, the elastic wavefield should be adequately decoupled before imaging. A common approach is to apply the curl and divergence operators to the coupled wavefield [13]. However, the decoupled S-wave presenting the polarity reversal phenomenon ultimately destroys the continuity of thereflection lineups of multi-shot stacked profiles. [PERSON] et al. [14] argued that the polarity of the decoupled S-wave is reversed on both sides of the vertical incident direction of the P-wave. With the assistance of the Poynting vector, a sign factor can be estimated to correct the polarity of PS imaging [15]. However, the divergence and curl operators will change the amplitude and phase of wavefields. As a consequence, it is difficult for ERTM to obtain a definite physical meaning from the converted wave imaging profile. To address this problem, the pure P- and S-wave can be effortlessly decoupled by solving the vector-decoupled wave equation [16].  
When applying the correlation imaging condition, the noises inevitably accompanying valuable information will debase the quality of migration profiles, especially on the interface where the velocity varies dramatically. Although the high-pass filter can suppress migration noises, it alters the waveform of the original image [17]. Another image denoising approach whereby the high angle information of the common imaging gathers in the angle domain is eliminated has advantages in amplitude-preserving [18; 19]. [PERSON] et al. [20] first suggested that the migration noises are related to the wave propagation path through the scattering point, and finally proposed an imaging condition with wavefield decomposition to remove the noises. Unfortunately, the initial wavefield decomposition method is realized in the wavenumber-frequency domain, which may result in a huge burden in computation, storage and input/output. [PERSON] et al. [21] imported a cost-effective scheme to decompose the wavefield based on the Hilbert transform. Further, [PERSON] et al. [22] decomposed the wavefield into several components according to the propagating directions. Particularly, the Poynting vector can indicate the direction of energy transfer, thereby flexibly realizing wavefield decomposition [23].  
The paper is arranged as follows: after the introduction section, the vector-decoupled wave equation is reviewed. On this basis, a novel equation is derived to extract the wavefield associated with insufficient illumination. Next, an inner product imaging condition with filters is described. Finally, some synthetic examples are employed to demonstrate the effectiveness of our approach.