Preserving the discontinuities in least-squares reverse time migration of simultaneous-source data

The simultaneous-source shooting technique can accelerate field acquisition and improve spatial sampling but it will cause strong interferences in the recorded data and artifacts in the final image. The previously proposed structural smoothing operator can effectively attenuate artifacts for relatively simple reflection structures during least-squares inversion, but it will cause damage to complicated reflection events such as discontinuities. To preserve discontinuities in a seismic image, we apply the singular spectrum analysis (SSA) operator to attenuate artifacts during least-squares inversion. Considering that global SSA cannot deal with over-complicated data very well, we use local SSA to remove noise and to better preserve the steeply dipping components. The local SSA operator corresponds to a local low-rank constraint applied in the inversion process. The migration operator used in the study is the reverse time migration (RTM) operator. Tests using the Marmousi model showed the superior performance of the proposed algorithm in preserving the discontinuities of seismic images.