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.

METHODS FOR THE MODEL OF l_1-NORM MINIMIZATION IN SEISMIC INVERSE

Geophysical Data Sparse Reconstruction Based on L0-Norm Minimization

L0+L1+L2 mixed optimization: a geometric approach to seismic imaging and inversion using concepts in topology and semigroup

L-1/2 NORM REGULARIZATION FOR 3D SEISMIC DATA INTERPOLATION

Riemannian Linearized Proximal Algorithms for Nonnegative Inverse Eigenvalue Problem.

Linear Convergence of Inexact Descent Method and Inexact Proximal Gradient Algorithms for Lower-Order Regularization Problems

L1−2 Minimization for Exact and Stable Seismic Attenuation Compensation

Adaptive mixed-norm seismic inversion for non-Gaussian errors

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

Developing Lagrangian-based Methods for Nonsmooth Nonconvex Optimization

A Proximal Iteratively Regularized Gauss–Newton Method for Nonlinear Inverse Problems

Model Misfit Minimization

REGULARIZING ITERATIVE METHOD FOR SOLVING INVERSE PROBLEM OF ONE-DIMENSIONAL WAVE EQUATION

Seismic Acoustic Impedance Inversion Using Reweighted L1-Norm Sparse Constraint

Smoothed Proximal Lagrangian Method for Nonlinear Constrained Programs

Pre-stack seismic inversion based onl1-2-norm regularized logarithmic absolute misfit function

On the Convergence of an Active-Set Method for ℓ1 Minimization

On Partly Smoothness, Activity Identification and Faster Algorithms of L 1 over L 2 Minimization

A NEW INVERSION METHOD OF TIME-LAPSE SEISMIC

The variant of the iterative shrinkage-thresholding algorithm for minimization of the l 1 over l ∞ norms

Seismic Impedance Inversion Using Approximation Zero-norm Minimization with SQP Method