ABSTRACT
For vertical seismic profile (VSP), free-surface multiples can provide much wider subsurface illumination when compared with primaries. However, migration of multiple reflections generates not only the desired image of reflection interfaces but also many crosstalk artifacts. Therefore, the least-squares reverse time migration method is used to image the VSP downgoing free-surface multiples (receiver-side ghosts) and iteratively suppress crosstalks, in which full downgoing data (including direct waves) and downgoing multiples are used as sources and observed data, respectively. To reduce the computational cost, we have developed the simultaneous imaging of different common-receiver gathers that are dynamically blended together with iterations through the altered realizations of the phase-encoding function. Relative to the popular encoding function with a combination of random time delays and polarities, only the random polarities can be applied for further increasing the computational efficiency. Synthetic experiments on Sigsbee2B and Pluto1.5 models indicate that the proposed method can effectively eliminate crosstalk artifacts and improve imaging resolution while calculated even more efficiently than reverse time migration of VSP ghosts.