English

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

Numerical Analysis 2022-06-13 v1 Numerical Analysis

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 0\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 0\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 0\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 1\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 1\ell_1 norm in terms of the SNR values and visual quality of the restored synthetic seismograms.

Keywords

Cite

@article{arxiv.2206.05129,
  title  = {Inverting Incomplete Fourier Transforms by a Sparse Regularization Model and Applications in Seismic Wavefield Modeling},
  author = {Tingting Wu and Yuesheng Xu},
  journal= {arXiv preprint arXiv:2206.05129},
  year   = {2022}
}
R2 v1 2026-06-24T11:46:38.881Z