Seismic impedance inversion using approximation zero-norm minimization with SQP method

Inversion for seismic impedance is an inherently complicated problem. It is ill-posed and band-limited. Thus the inversion results are non-unique and the inversion process is unstable. To take the sparsity of the seismic reflectivity into consideration, the inversion problem can be transformed into a zero-norm constraint optimization. The resulting minimization problem is proposed to be efficiently solved by the trust region based Sequential Quadratic Programming (SQP) method in this paper, which possesses global convergence and local quadratic convergence rate under suitable assumptions. Null space technique is considered to reduce the computational complexity and greatly accelerate computation speed. Numerical experiments for synthetic simulations are given. The results illustrate that the proposed approach has the advantages of high precision, robustness to noise and fast calculation.