Abstract:Seismic interpolation can provide a complete wavefield from incomplete wavefield information and is crucial for some processings such as multiple elimination, wave equation migration and denoising. Currently, most interpolation methods are suitable for high signal-to-noise ratio data. However, filed seismic data often contain different kinds of noise, especially the random noise. The existence of noise can deteriorate the result of interpolation. The weighted projection onto convex sets (POCS) method is a method for random noise contained data interpolation. However, there is no specific explanation of this method except the iteration formula. This paper analyzed the algorithm theory of the weighted POCS method, and proved that it is a method for unconstrained optimization problem, the weighting factor can be regarded as the coefficient of the fitting error. Meanwhile, we also proposed an improved POCS method for simultaneous interpolation and denoising. This algorithm requires no additional computational effort compared with the original POCS method. Seismic interpolation especially the high dimensional interpolation belongs to large-scale computational problems, thus fast algorithms should be chosen to improve the efficiency of seismic interpolation. Sparse transform methods change the interpolation problem into a sparse optimization based on the assumption that seismic data can be expressed sparsely in some transform domain. POCS method is such a method with high efficiency and robust results. However, this method can not get acceptable results for random noise contained data interpolation. The weighted POCS method is a method for random noise contained data interpolation on the basis of the POCS method. But there is no specific explanation of this method up to now. This paper deduced the weighted POCS method from the iterative thresholding method, we proved that the weighted POCS method is a method for unconstrained optimization, the weighted factor can be seen as the coefficient of the fitting error. However, the original FOGS method is aiming at an equality constrained optimization. This is the essential difference between them. Meanwhile, we also proposed an improved FOGS method with the same computational complexity as the original POCS method. The improved POCS method can realize interpolation and denoising simultaneously by adjustment of the least thresholds. In order to test the interpolation and denoising capability of the improved POCS method, four data sets were used to make numerical experiments. The first experiment is a synthetic noiseless data interpolation. The dimension of the synthetic data is 256 X 256 with trace distance 25 m and time sampling rate 4 ms. The original POCS method is chosen as the benchmark. Both the computational efficiency and numerical result are comparable for them. The second data is field data with 115 traces. The trace distance is 25 m and the time sampling rate is 4 ms. This experiment also proved that both of them can get comparable results for noiseless data interpolation. The third example is a noise contained data interpolation. The original FOGS method can not get acceptable result no matter how to choose the least thresholds, while the improved POCS method can provide good result by choosing the least thresholds carefully. The fourth example also proved that the original FOGS method failed in the noise contained field data interpolation, and the improved FOGS method is suitable for noise contained data interpolation. Numerical experiment on first-interpolation-then-denoising proved the rationality and robustness of simultaneous interpolation and denoising. We deduced the weighted FOGS method from the iterative thresholding method, and Proved that it is a method for unconstrained optimization, the weighting factor is the coefficient of the fitting error. However, the original FOGS method is a method for equality constrained optimization. This is the essential difference between them. We also proposed an improved FOGS method to realize simultaneous interpolation and denoising. This method just need to adjust the least thresholds to realize noiseless data and random noise contained data interpolation. The computational effort of the improved FOGS method is the same as the original FOGS method but with more broad applications.

A NOVEL METHOD FOR SIMULTANEOUS SEISMIC DATA INTERPOLATION AND NOISE REMOVAL BASED ON THE L-0 CONSTRAINT

An Improved Weighted Projection Onto Convex Sets Method for Seismic Data Interpolation and Denoising

AN IMPROVED PROJECTION ONTO CONVEX SETS METHOD FOR SIMULTANEOUS INTERPOLATION AND DENOISING

Comparison of the Projection Onto Convex Sets and Iterative Hard Thresholding Methods for Seismic Data Interpolation and Denoising

Simultaneous Seismic Data Interpolation and Denoising with a New Adaptive Method Based on Dreamlet Transform

A fast-convergence POCS seismic denoising and reconstruction method

A Novel Method for Simultaneous Interpolation and Denoising Based on L0 Norm Constraint

Comparison of the POCS and IHT Methods for Seismic Data Interpolation and Denoising

Separation Method for Multi-Source Blended Seismic Data

Highly efficient methods of seismic data acquisition based on sparse constraints and blended-source excitation

Simultaneous Seismic Data Interpolation and Denoising Based on Nonsubsampled Contourlet Transform Integrating With Two-Step Iterative Log Thresholding Algorithm

An Anti-Aliasing POCS Interpolation Method for Regularly Undersampled Seismic Data Using Curvelet Transform

Curvelet-based POCS Interpolation of Nonuniformly Sampled Seismic Records

A Projection-Onto-convex-sets Network for 3D Seismic Data Interpolation

An Efficient POCS Interpolation Method in the Frequency-Space Domain

Seismic Random Noise Attenuation Using Nonlocal-Means Via Smooth Patch Ordering

On Analysis-Based Two-Step Interpolation Methods for Randomly Sampled Seismic Data

Simultaneous seismic data de‐aliasing and denoising with a fast adaptive method based on hybrid wavelet transform

Dreamlet-based interpolation using POCS method

Accurate and Efficient Seismic Data Interpolation in the Principal Frequency Wavenumber Domain

An iterative half thresholding method for seismic data interpolation