Normalized shaping regularization for robust separation of blended data

Simultaneous source acquisition has attracted more and more attention from geophysicists because of its cost savings, whereas it also brings some challenges that have never been addressed before. Deblending of simultaneous source data is usually considered as an underdetermined inverse problem, which can be effectively solved with a least-squares (LS) iterative procedure between data consistency ([Formula: see text]-norm) and regularization ([Formula: see text]-norm or [Formula: see text]-norm). However, when it comes to abnormal noise that follows non-Gaussian distribution and possesses high-amplitude features (e.g., erratic noise, swell noise, and power line noise), the [Formula: see text]-norm is a nonrobust statistic that can easily lead to suboptimal deblended results. Although abnormal noise can be attenuated in the common source domain at first, it is still challenging to apply a coherency-based filter due to the sparse receiver or crossline sampling, e.g., that commonly found in ocean bottom node (OBN) acquisition. To address this problem, we have developed a normalized shaping regularization to make the inversion-based deblending approach robust for the separation of blended data when abnormal noise exists. Its robustness comes from the normalized shaping operator defined by the confidence interval of normal distribution, which minimizes the abnormal risk to a normal level to satisfy the assumption of LS shaping regularization. In special cases, the proposed approach will revert to the classic LS shaping regularization once the normalized coefficient is large enough. Experimental results on synthetic and field data indicate that the proposed method can effectively restore the separated records from blended data at essentially the same convergence rate as the LS shaping regularization for the abnormal noise-free scenario, but it can obtain better deblending performance and less energy leakage when abnormal noise exists.