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.

Modified Fourier-Hankel Method Based On Analysis Of Errors In Abel Inversion Using Fourier Transform Techniques

A New Method for Numerical Abel-Inversion

Tomography of cylindrical objects: comparison of noise property and accuracy of Abel inversion techniques with and without noise filtering

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

High-precision Fourier forward modeling of potential fields

Dealing with Asymmetrical Abel Inversion Using Fourier-Bessel Expansions

Abel inversion using fast Fourier transforms.

Interpolant Method in Analyzing the Laser-Induced Plasma Spectroscopy

Fast Calculation of Laplace-Beltrami Eigenproblems Via Subdivision Linear Subspace

Computation Of The Laplace Inverse Transform By Application Of The Wavelet Theory

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

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

DEALING WITH ABEL INVERSION BY USING LENGENDRE EXPANSIONS

A Numerical Inversion of the Laplace Transform By Use of the Scaling Function Transform of Wavelet Theory

NUMERICAL INVERSION OF LAPLACE TRANSFORS IN VISCOELASTIC PROBLEMS BY FOURIER SERIES EXPANSION

Inexact line search method in full waveform inversion

Improved Smoothing Frequency Shifting and Filtering Algorithm for Harmonic Analysis With Systematic Error Compensation

On The Aliasing Artifacts In Direct Fourier Methods

Spurious eigenvalue-free algorithms of the method of fundamental solutions for solving the Helmholtz equation in bounded multiply connected domains

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory