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.

Tangential Cone Condition for the Full Waveform Forward Operator in the Viscoelastic Regime: The Nonlocal Case

Full Waveform Inversion Based on Time-Integral-damping Wavefield

A Time-Domain Preconditioned Truncated Newton Approach to Visco-acoustic Multiparameter Full Waveform Inversion

Calderón problem for nonlocal viscous wave equations: Unique determination of linear and nonlinear perturbations

Well-posedness and inverse problems for semilinear nonlocal wave equations

Direct error in constitutive equation formulation for inverse problems in time harmonic elastography

Attenuation Sensitivity Kernel Analysis in Viscoelastic Full-Waveform Inversion Based on the Generalized Standard Linear Solid Rheology

Inverse Source Problems In Elastodynamics

Inversions of two wave-forward operators with variable coefficients

Full Waveform Inversion of Viscoelastic Media Based on Gradient Preconditioning

Multi-level, Multi-frequency, compressed Full Waveform Inversion

Immersed boundary parametrizations for full waveform inversion

Full waveform inversion in fractured media based on velocity–stress wave equations in the time domain

Numerical analysis of a deep learning formulation of elastic full waveform inversion with high order total variation regularization in different parameterization

Sensitivity Analysis of an Inverse Problem for the Wave Equation with Caustics

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

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

Efficient time-domain 3D elastic and viscoelastic full-waveform inversion using a spectral-element method on flexible Cartesian-based mesh

Iterative PDE-Constrained Optimization for Seismic Full-Waveform Inversion

3-D Time-Domain VTI Viscoacoustic Full Waveform Inversion: Application to Valhall Field Data

Inverse boundary value problem for the helmholtz equation with multi-frequency data