In this research, we present a seismic trace interpolation method which uses seismic data with surface-related multiples. It is different from conventional seismic data interpolation using information transformation or extrapolation of adjacent channels for reconstruction of missing seismic data. In this method there are two steps, first, we construct pseudo-primaries by cross-correlation of surface multiple data to extract the missing near- offset information in multiples, which are not displayed in the acquired seismic record. Second, we correct the pseudo-primaries by applying a Least-squares Matching Filter （LMF） and RMS amplitude correction method in time and space sliding windows. Then the corrected pseudo-primaries can be used to fill the data gaps. The method is easy to implement, without the need to separate multiples and primaries. It extracts the seismic information contained by multiples for filling missing traces. The method is suitable for seismic data with surfacerelated multiples.
Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be one of highly accurate methods. In this paper, we propose a new non-reflecting recursive algorithm for reverse-time migration by introducing the wave impedance function into the acoustic wave equation and the algorithm for the surface relief case is derived from the coordinate transformation principle. Using the exploding reflector principle and the zero-time imaging condition of poststack reverse- time migration, poststack numerical simulation and reverse-time migration with complex conditions can be realized. The results of synthetic and real data calculations show that the method effectively suppresses unwanted internal reflections and also deals with the seismic imaging problems resulting from surface relief. So, we prove that this method has strong adaptability and practicality.
