Abstract: High-velocity salt domes and steeply dipping structures in subsurface media pose challenges to conventional imaging methods.In conventional seismic exploration,multiples are usually treated as noise to be suppressed or eliminated.However,due to their long propagation paths and smaller reflection angles within subsurface media,multiples show a wider imaging range and richer information on subsurface structures.Since seismic interferometry enables accurate prediction of multiples,the reverse time migration(RTM) of multiples based on seismic interferometry is crucial for imaging complex subsurface structures.Based on the analysis of the fundamental principles of seismic interferometry,this study derived the corresponding equations and algorithms.Furthermore,this study proposed a multiple prediction technique based on seismic interferometry.By circumventing conventional multiple processing methods,this technique enhances imaging efficiency, providing a research basis for subsequent RTM of multiples.To address the significant crosstalk noise in conventional RTM,this study innovatively utilized seismic interferometry to explore the RTM of order-divided multiples.Finally,this study proposed a RTM method for order-divided multiples based on seismic interferometry,effectively mitigating the problem of crosstalk noise in multiple imaging.Numerical experiments show that compared to the RTM of primary waves,the RTM of multiples exhibited significant improvements in imaging pre-salt structures.The RTM of order-divided multiples based on seismic interferometry could further suppress the crosstalk noise in conventional RTM,enabling accurate imaging of complex pre-salt structures and improving imaging resolution.