Simultaneous seismic interpolation and denoising based on sparse inversion with a 3D low redundancy curvelet transform

The simultaneous seismic interpolation and denoising problem can be solved as a sparse inversion problem by using the sparseness of seismic data in a transformed domain as the a priori information, where the properties of the sparse transform will significantly influence the numerical results and computational efficiency. Curvelet transform has nearly optimal sparse expression for seismic data, thus seismic signal processing based on this transform tends to result in preferable results. However, the redundancy of this transform can be 24–32 for three dimensional data which is not computationally cost efficient. This paper introduces a low redundancy curvelet transform to simultaneously interpolate and denoise. The redundancy of the proposed transform can be reduced to 10 for three dimensional data, and this property will improve the computational efficiency of the curvelet transform-based signal processing. The iterative soft thresholding method was chosen to solve the sparse inversion problem. Some practical principles on how to choose the regularisation parameters are discussed, due to the crucial nature of the regularisation parameter for simultaneous interpolation and denoising. Numerical experiments on synthetic and field data demonstrate that the low redundancy transform can provide reliable results while at the same time improving the computational efficiency.