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.

Exponential relaxation data analysis by parametrized regularization of severely ill-posed Fredholm integral equations of the first kind

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

An adaptive RKHS regularization for Fredholm integral equations

Regularization with two differential operators and its application to inverse problems

Regularization of linear ill-posed problems involving multiplication operators

Influence of the order between discretization and regularization in solving ill-posed problems

A relaxed iterated Tikhonov regularization for linear ill-posed inverse problems

Regularized solutions for some backward nonlinear parabolic equations with statistical data

A fast iterative regularization method for ill-posed problems

Adaptive optimal regularization of the linear ill posed integral equations

A sparse optimization approach to infinite infimal convolution regularization

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

Regularization of the discrete source problem in the nonlinear diffusive-logistic equation

Error Estimates for Data-driven Weakly Convex Frame-based Image Regularization

On convergence and convergence rates for Ivanov and Morozov regularization and application to some parameter identification problems in elliptic PDEs

A relaxation approach to modeling properties of hyperbolic-parabolic type models

On regularization methods based on dynamic programming techniques

Discrepancy Curves for Multi-Parameter Regularization

Projected Regression Methods for Inverting Fredholm Integrals: Formalism and Application to Analytical Continuation

On Learning the Optimal Regularization Parameter in Inverse Problems

Refined Exponential Filter with Applications to Image Restoration and Interpolation