Accelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet

jingjie cao,yanfei wang,benfeng wang

DOI: https://doi.org/10.1071/EG14016

2015-01-01

Exploration Geophysics

Abstract:Seismic interpolation, as an efficient strategy of providing reliable wavefields, belongs to large-scale computing problems. The rapid increase of data volume in high dimensional interpolation requires highly efficient methods to relieve computational burden. Most methods adopt the L-1 norm as a sparsity constraint of solutions in some transformed domain; however, the L-1 norm is non-differentiable and gradient-type methods cannot be applied directly. On the other hand, methods for unconstrained L-1 norm optimisation always depend on the regularisation parameter which needs to be chosen carefully. In this paper, a fast gradient projection method for the smooth L-1 problem is proposed based on the tight frame property of the curvelet transform that can overcome these shortcomings. Some smooth L-1 norm functions are discussed and their properties are analysed, then the Huber function is chosen to replace the L-1 norm. The novelty of the proposed method is that the tight frame property of the curvelet transform is utilised to improve the computational efficiency. Numerical experiments on synthetic and real data demonstrate the validity of the proposed method which can be used in large-scale computing.