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.

Chebyshev polynomial-based Fourier transformation and it's use in low pass filter of gravity data

High-resolution time-frequency hilbert transform using sparsity-aware weighting function

Sparsity-Enhanced Constrained Least-Squares Spectral Analysis with Greedy-FISTA

Time-Frequency Phase Retrieval for Deconvolutive Short-Time Fourier Transform

A Cascaded Synchrosqueezing Transform for Precise Analysis of Seismic Signal

Data-Driven Filter Design in FBP: Transforming CT Reconstruction with Trainable Fourier Series

Efficient anti-aliasing and anti-leakage fourier transform for high-dimensional seismic data regularization using cube removal and GPU

The Synchrosqueezed Method and Its Theory-Analysis-Based Novel Short-Time Fractional Fourier Transform for Chirp Signals

Application of Fractional Fourier Transform in Cepstrum Analysis

Chebyshev and The Fast Fourier Transform Methods for Signal Interpolation

Trend Filtering -- I. A Modern Statistical Tool for Time-Domain Astronomy and Astronomical Spectroscopy

Improved Reconstruction for Fourier-Sparse Signals

Real-Time Joint Filtering of Gravity and Gravity Gradient Data Based on Improved Kalman Filter

A robust method for instantaneous frequency estimation via inverse spectral decomposition

Optimized sparse fractional Fourier transform: Principle and performance analysis

Accurate Estimators of Correlation Functions in Fourier Space

Fourier-Based Approach To Interpolation In Single-Slice Helical Computed Tomography

Robust and high-resolution seismic complex trace analysis

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

A NEW TOOL FOR IMAGE ANALYSIS BASED ON CHEBYSHEV RATIONAL FUNCTIONS: CHEF FUNCTIONS