Accelerated parabolic Radon domain 2D adaptive multiple subtraction with fast iterative shrinkage thresholding algorithm and its application in parabolic Radon domain hybrid demultiple method

Adaptive multiple subtraction is an important step for successfully conducting surface-related multiple elimination in marine seismic exploration. 2D adaptive multiple subtraction conducted in the parabolic Radon domain has been proposed to better separate primaries and multiples than 2D adaptive multiple subtraction conducted in the time-offset domain. Additionally, the parabolic Radon domain hybrid demultiple method combining parabolic Radon filtering and parabolic Radon domain 2D adaptive multiple subtraction can better remove multiples than the cascaded demultiple method using time-offset domain 2D adaptive multiple subtraction and the parabolic Radon transform method sequentially. To solve the matching filter in the optimization problem with L1 norm minimization constraint of primaries, traditional parabolic Radon domain 2D adaptive multiple subtraction uses the iterative reweighted least squares (IRLS) algorithm, which is computationally expensive for solving a weighted LS inversion in each iteration. In this paper we introduce the fast iterative shrinkage thresholding algorithm (FISTA) as a faster alternative to the IRLS algorithm for parabolic Radon domain 2D adaptive multiple subtraction. FISTA uses the shrinkage-thresholding operator to promote the sparsity of estimated primaries and solves the 2D matching filter with iterative steps. FISTA based parabolic Radon domain 2D adaptive multiple subtraction reduces the computation time effectively while achieving similar accuracy compared with IRLS based parabolic Radon domain 2D adaptive multiple subtraction. Additionally, the provided examples show that FISTA based parabolic Radon domain 2D adaptive multiple subtraction can better separate primaries and multiples than FISTA based time-offset domain 2D adaptive multiple subtraction. Furthermore, we introduce FISTA based parabolic Radon domain 2D adaptive multiple subtraction into the parabolic Radon domain hybrid demultiple method to improve its computation efficiency while maintaining its accuracy. Synthetic and field data examples demonstrate the effectiveness of the proposed FISTA based parabolic Radon domain 2D adaptive multiple subtraction and its application in the parabolic Radon domain hybrid demultiple method.