Considering the 3D propagation characteristics of seismic waves, theoretically, 3D surface-related multiples elimination (3D SRME) can suppress multiples with high accuracy. However, 3D SRME has strict requirements for acquisition geometry, which makes it difficult to be implemented in practice. In the process of 3D SRME, the multiple contribution gather (MCG) is a collection of wavefields with different propagation paths. The accuracy of the multiple propagation paths in the MCGs can be directly characterized by the inclination of the wavefields, which can achieve the weighted superposition of the wavefields. The direct summation of the sparse MCGs in the crossline direction produces serious spatial aliasing, which can easily cause the contamination of primaries. Based on the kinematic characteristics of multiple propagation, MCGs can be considered as a set of hyperbolas with temporal and spatial characteristics. Then, the direct summation of the sparse MCGs can be transformed into a process of superposition along the hyperbolic integration paths. However, as the stable phase points of the events, the apexes of the hyperbola have different spatial distributions in complex geological structures. Such hyperbolic stacking paths are difficult to be controlled by conventional Radon transform or constrained inversion. In this paper, we modify the apex-shifted hyperbolic Radon transform (ASHRT) to implement the summation of crossline MCGs with variable stable phase points along the hyperbolic integration paths. Improved ASHRT uses local similarity to locate the position of stable phase points, which can improve the stability of the algorithm and the efficiency of the computation. The proposed method is demonstrated on a 3D synthetic data set, as well as on a 3D marine data set, effectively avoiding the spatial aliasing caused by sparse crossline MCGs and improving the accuracy of multiple suppression.

Passive seismic has recently attracted a great deal of attention because non-artificial source is used in subsurface imaging. The utilization of passive source is low cost compared with artificial-source exploration. In general, constructing virtual shot gathers by using cross-correlation is a preliminary step in passive seismic data processing, which provides the basis for applying conventional seismic processing methods. However, the subsurface structure is not uniformly illuminated by passive sources, which leads to that the ray path of passive seismic does not fit the hyperbolic hypothesis. Thereby, travel time is incorrect in the virtual shot gathers. Besides, the cross-correlation results are contaminated by incoherent noise since the passive sources are always natural. Such noise is kinematically similar to seismic events and challenging to be attenuated, which will inevitably reduce the accuracy in the subsequent process. Although primary estimation for transient-source seismic data has already been proposed, it is not feasible to noise-source seismic data due to the incoherent noise. To overcome the above problems, we proposed to combine focal transform and local similarity into a highly integrated operator and then added it into the closed-loop surface-related multiple elimination based on the 3D L1-norm sparse inversion framework. Results proved that the method was capable of reliably estimating noise-free primaries and correcting travel time at far offsets for a foresaid virtual shot gathers in a simultaneous closed-loop inversion manner.