Abstract

We consider numerical solution methods for seismic impedance inversion problems in this paper. The inversion process is ill-posed. To tackle the ill-posedness of the problem and take the sparsity of the reflectivity function into consideration, an l ₁ norm regularization model is established. In computation, a nonmonotone gradient descent method based on Rayleigh quotient for solving the minimization model is developed. Theoretical simulations and field data applications are performed to verify the feasibility of our methods.