Abstract:We propose a sparse regularization model for inversion of incomplete Fourier transforms and apply it to seismic wavefield modeling. The objective function of the proposed model employs the Moreau envelope of the $\ell_0$ norm under a tight framelet system as a regularization to promote sparsity. This model leads to a non-smooth, non-convex optimization problem for which traditional iteration schemes are inefficient or even divergent. By exploiting special structures of the $\ell_0$ norm, we identify a local minimizer of the proposed non-convex optimization problem with a global minimizer of a convex optimization problem, which provides us insights for the development of efficient and convergence guaranteed algorithms to solve it. We characterize the solution of the regularization model in terms of a fixed-point of a map defined by the proximity operator of the $\ell_0$ norm and develop a fixed-point iteration algorithm to solve it. By connecting the map with an $\alpha$-averaged nonexpansive operator, we prove that the sequence generated by the proposed fixed-point proximity algorithm converges to a local minimizer of the proposed model. Our numerical examples confirm that the proposed model outperforms significantly the existing model based on the $\ell_1$-norm. The seismic wavefield modeling in the frequency domain requires solving a series of the Helmholtz equation with large wave numbers, which is a computationally intensive task. Applying the proposed sparse regularization model to the seismic wavefield modeling requires data of only a few low frequencies, avoiding solving the Helmholtz equation with large wave numbers. Numerical results show that the proposed method performs better than the existing method based on the $\ell_1$ norm in terms of the SNR values and visual quality of the restored synthetic seismograms.

Analysis of the SQP Method for Hyperbolic PDE-Constrained Optimization in Acoustic Full Waveform Inversion

Numerical Analysis for a Hyperbolic PDE-Constrained Optimization Problem in Acoustic Full Waveform Inversion

Convergence of the SQP method for quasilinear parabolic optimal control problems

A stabilized SQP method: superlinear convergence

A Globally Convergent Stabilized SQP Method

On the Performance of SQP Methods for Nonlinear Optimization

The reduced space Sequential Quadratic Programming (SQP) method for calculating the worst resonance response of nonlinear systems

Inverting Incomplete Fourier Transforms by a Sparse Regularization Model and Applications in Seismic Wavefield Modeling

An SQP method for mathematical programs with vanishing constraints with strong convergence properties

Superlinear convergence of a stabilized SQP-type method for nonlinear semidefinite programming

Global and local convergence of a nonmonotone SQP method for constrained nonlinear optimization

A sequential quadratic programming method for nonsmooth stochastic optimization with upper-C^2 objective

Sequential quadratic programming methods for parametric nonlinear optimization

Modified Symplectic Schemes with Nearly-Analytic Discrete Operators for Acoustic Wave Simulations

A New SQP Algorithm for Large-Scale Nonlinear Programming

Analysis of Sequential Quadratic Programming through the Lens of Riemannian Optimization

A method to guarantee local convergence for sequential quadratic programming with poor Hessian approximation

A stabilized filter SQP algorithm for nonlinear programming

Symplectic Interior Penalty Discontinuous Galerkin Method for Solving the Seismic Scalar Wave Equation

Quasi-Periodic Solutions of a Damped Nonlinear Quasi-Periodic Mathieu Equation by the Incremental Harmonic Balance Method With Two Time Scales