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.

A Regularization Homotopy Method for the Inverse Problem of 2-D Wave Equation and Well Log Constraint Inversion

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

Regularization with two differential operators and its application to inverse problems

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

Regularity of an inverse problem in wave propagation

Well-posedness and Tikhonov regularization of an inverse source problem for a parabolic equation with an integral constraint

Acoustic- and elastic-waveform inversion with total generalized p-variation regularization

A Case Study of the Regularization Methods for Ill-Posed Inverse Problems

An asymptotical regularization with convex constraints for inverse problems

An iteration regularization method with general convex penalty for nonlinear inverse problems in Banach spaces

General Optimization Framework for Robust and Regularized 3D Full Waveform Inversion

A Homotopy Method for the Constrained Inverse Problem in the Multiphase Porous Media Flow

Positive Jacobian constraints for elliptic boundary value problems with piecewise-regular coefficients arising from multi-wave inverse problems

Weighted inhomogeneous regularization for inverse problems with indirect and incomplete measurement data

Least squares formulation for ill-posed inverse problems and applications

An Analysis of Tikhonov Regularization for Nonlinear Ill-Posed Problems under a General Smoothness Assumption

Tikhonov regularization for simultaneous inversion of initial value and source term of a time‐fractional Black‐Scholes equation

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory

A Quasi Tikhonov Regularization for a Two-Dimensional Backward Heat Problem by a Fundamental Solution

Study on Global Optimization Methods for Wave Equation Inversion Problems