where is the well-sampled seismic data, denotes the sampling operator, and denotes the observed data Chen et al. (2014a). Because of the missing traces, the operator is highly singular and thus make equation 1 highly under-determined.

There have existed many different algorithms for solving equation 1 by adding different constraints. The POCS algorithm Abma and Kabir (2006) is one of the most widely used methods to interpolate seismic data with irregularly missing traces, and has the following iterative expression:

where denotes the reconstructed data after th iteration, and are the forward and inverse sparsity-promoting transforms, is an identity matrix and is a thresholding operator. However, for regularly missing traces, POCS can not obtain satisfying results because of the strong aliasing noise in the sparse domain

2020-02-21