Closed-loop SRME based on 3D L1-norm sparse inversion

In many situations, the quality of seismic imaging is largely determined by a proper multiple attenuation as preprocessing step. Despite the widespread application of surface-related multiple elimination (SRME) and estimation of primaries by sparse inversion (EPSI) for the removal of multiples, there still exist some limitations in the process of prediction and subtraction (SRME) or inversion (EPSI), which make the efficiency of multiple attenuation less satisfactory. To solve these problems, a new fully data-driven method called closed-loop SRME was proposed, which combines the robustness of SRME and the multi-dimensional inversion strategy of EPSI. Due to the selection of inversion approach and constraint, primary estimation by closed-loop SRME may fall into a local optimum during the solving process, which lowers the accuracy of deep information and weakens the continuity of seismic events. To avoid these shortcomings, we first modified the solving method for closed-loop SRME to an L1 norm-based bi-convex optimization method, which stabilizes the solution. Meanwhile, in the L1 norm constraint-based optimization process, the 3D sparsifying transform, being a 2D Curvelet-1D wavelet transform, is brought in as a 3D sparse constraint. In the 3D sparsifying domain, the data become sparser, thus making the result of optimization more accurate, the information of seismic events more continuous and the resolution higher. Examples on both synthetic and field data demonstrate that the method proposed in this paper, compared with the traditional SRME and closed-loop SRME, have an excellent effect on primary estimation and suppress multiples effectively.