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.

WAVELET METHODS FOR THE SIDEWAYS HEAT EQUATION

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

An accurate regularization technique for the backward heat conduction problem with time-dependent thermal diffusivity factor

A quasi-reversibility method for solving nonhomogeneous sideways heat equation

A quasi-boundary value regularization method for the spherically symmetric backward heat conduction problem

Wavelet-Based Regularization of Dynamic Data Reconciliation

A variational regularization method for solving the non-characteristic Cauchy problem in multiple dimensions

The Symmetric Regularized-Long-Wave Equation: Ill-posedness and Long Period Limit

Sharp Global Existence for Semilinear Wave Equation with Small Data

Solving two kinds of inverse source problems for the heat equations by a mollification regularization method with Dirichlet kernel

Two Regularization Methods for Identifying the Spatial Source Term Problem for a Space-Time Fractional Diffusion-Wave Equation

Multiscale-Wavelet Regularization Method for the Inverse Problem of Electrocardiography

Recovering the temperature distribution for multi-term time-fractional sideways diffusion equations

Regularized solutions for some backward nonlinear parabolic equations with statistical data

Regularizing a two-dimensional time-fractional inverse heat conduction problem by a fractional Landweber iteration method

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

Simultaneous identification of the initial data in a degenerate and singular hyperbolic wave equation

A wavelet method for solving a class of nonlinear boundary value problems

On a degenerate boundary value problem to the two-dimensional self-similar nonlinear wave system

Least squares formulation for ill-posed inverse problems and applications

Method of fundamental solutions with an optimal regularization technique for the Cauchy problem of the modified Helmholtz equation