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.

Weak-norm Posterior Contraction Rate of the 4DVAR Method for Linear Severely Ill-Posed Problems

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

Numerical Solution for Nonlinear 4D Variational Data Assimilation (4D-Var) via ADMM

A scalable space-time domain decomposition approach for solving large-scale nonlinear regularized inverse ill-posed problems in 4D variational data assimilation

Saddle point preconditioners for weak-constraint 4D-Var

On the use of the saddle formulation in weakly-constrained 4D-VAR data assimilation

A Time-parallel Approach to Strong-constraint Four-dimensional Variational Data Assimilation

Strong posterior contraction rates via Wasserstein dynamics

Analysis of 4D-Variational Data Assimilation Problems in Low Regularity Spaces

Posterior contraction rates in a sparse non-linear mixed-effects model

A posterior contraction for Bayesian inverse problems in Banach spaces

Weak convergence and invariant measure of a full discretization for non-globally Lipschitz parabolic SPDE

Strong and weak convergence rates of finite element method for stochastic partial differential equation with non-sided Lipschitz coefficient

On strong posterior contraction rates for Besov-Laplace priors in the white noise model

A 4D-Var Method with Flow-Dependent Background Covariances for the Shallow-Water Equations

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

A posteriori error analysis of an ultra-weak local discontinuous Galerkin method for nonlinear fourth-order boundary-value problems

l 1–2 -norm regularized basis pursuit seismic inversion based on exact Zoeppritz equation

Oracle-type Posterior Contraction Rates in Bayesian Inverse Problems

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