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.

Using learned priors to regularize the Helmholtz equation least-squares method

Fast Bayesian hierarchical inference via sparsity enforcing a priori for aeroacoustic source imaging

Robust Bayesian Super-Resolution Approach Via Sparsity Enforcing a Priori for Near-Field Aeroacoustic Source Imaging

A Variational Bayesian Approximation Approach Via A Sparsity Enforcing Prior In Acoustic Imaging

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

An Efficient Variational Bayesian Inference Approach Via Studient's-t Priors for Acoustic Imaging in Colored Noises

A Hierarchical Variational Bayesian Approximation Approach in Acoustic Imaging

Proceedings of Meetings on Acoustics

STUDY ON RECONSTRUCTION OF COHERENT ACOUSTIC FIELD WITH MULTI-SOURCE BASED ON HELMHOLTZ EQUATION AND LEAST SQUARE METHOD

Statistical reduced order modelling for the parametric Helmholtz equation

A Novel Least Squares Method for Helmholtz Equations with Large Wave Numbers

Sound field reconstruction using sparse Bayesian learning equivalent source method with hyperparametric-coupled prior

Theorical and Experimental Research on the Hels Method in the Cyclostationary Acoustic Field

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

Half-Space Sound Field Reconstruction Based on the Combination of the Helmholtz Equation Least-Squares Method and Equivalent Source Method

An Investigation on Regularization Methods for Sound Source Reconstruction of Magnetoacoustic Tomography (Mat)

Learning based numerical methods for acoustic frequency-domain simulation with high frequency

Eigenvector Model Descriptors for Solving an Inverse Problem of Helmholtz Equation: Extended Materials

An iterative method for Helmholtz boundary value problems arising in wave propagation

Machine learning aided near-field acoustic holography based on equivalent source method

A Discontinuous Least Squares Finite Element Method for the Helmholtz Equation