[go: up one dir, main page]
More Web Proxy on the site http://driver.im/ Skip to main content

Advertisement

Springer Nature Link
Log in

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

  • Published:
Journal of Scientific Computing Aims and scope Submit manuscript

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. This makes the proposed model particularly suitable for the seismic wavefield (SW) modeling. 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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Subscribe and save

Springer+ Basic
£29.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price includes VAT (United Kingdom)

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8

Similar content being viewed by others

Data Availability

The datasets generated and analysed during the current study are not publicly available but are available from the corresponding author on reasonable request.

References

  1. Babuška, I., Sauter, S.A.: Is the pollution effect of the FEM avoidable for the Helmholtz equation considering high wave numbers? SIAM Rev. 42, 451–484 (2000)

    MathSciNet  MATH  Google Scholar 

  2. Borwein, J.M., Li, G., Tam, M.K.: Convergence rate analysis for averaged fixed point iterations in common fixed point problems. SIAM J. Optim. 27, 1–33 (2017)

    Article  MathSciNet  Google Scholar 

  3. Brigham, E.: The fast Fourier transform and its application. Prentice-Hall Inc., NJ (1988)

    Google Scholar 

  4. Cai, J., Chan, R., Shen, L., Shen, Z.: Convergence analysis of tight framelet approach for missing data recovery. Adv. Comput. Math. 31, 87–113 (2009)

    Article  MathSciNet  Google Scholar 

  5. Candès, E.J., Romberg, J., Tao, T.: Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 52, 489–509 (2006)

    Article  MathSciNet  Google Scholar 

  6. Candès, E.J., Romberg, J., Tao, T.: Stable signal recovery from incomplete and inaccurate measurements. Commun. Pure Appl. Math. 59, 1207–1223 (2006)

    Article  MathSciNet  Google Scholar 

  7. Chen, F., Shen, L., Xu, Y., Zeng, X.: The Moreau envelope approach for the L1/TV image denoising model. Inverse Problems and Imaging 8, 53–77 (2014)

    Article  MathSciNet  Google Scholar 

  8. Chen, Z., Cheng, D., Feng, W., Wu, T.: An optimal 9-point finite difference scheme for the Helmholtz equation with PML. Int. J. Numer. Anal. Model. 10, 389–410 (2013)

    MathSciNet  MATH  Google Scholar 

  9. Combettes, P.L., Yamada, I.: Compositions and convex combinations of averaged nonexpansive operators. J. Math. Anal. Appl. 425, 55–70 (2015)

    Article  MathSciNet  Google Scholar 

  10. Donoho, D.L.: Compressed sensing. IEEE Trans. Inf. Theory 52, 1289–1306 (2006)

    Article  MathSciNet  Google Scholar 

  11. Fan, J., Li, R.: Variable selection via non-concave penalized likelihood and its oracle properties. J. Am. Stat. Assoc. 96(456), 1348–1360 (2001)

    Article  Google Scholar 

  12. Goldstein, T., Osher, S.: The split Bregman method for \(l^1\) regularization problems. SIAM J. Imag. Sci. 2, 323–343 (2009)

    Article  Google Scholar 

  13. Horn, R.A., Johnson, C.R.: Matrix Analysis. Cambridge University Press, Cambridge (1985)

    Book  Google Scholar 

  14. Hustedt, B., Operto, S., Virieux, J.: Mixed-grid and staggered-grid finite-difference methods for frequency domain acoustic wave modelling. Geophys. J. Int. 157, 1269–1296 (2004)

    Article  Google Scholar 

  15. Jo, C.-H., Shin, C., Suh, J.H.: An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator. Geophysics 61, 529–537 (1996)

    Article  Google Scholar 

  16. Krol, A., Li, S., Shen, L., Xu, Y.: Preconditioned alternating projection algorithms for maximum a posteriori ECT reconstruction. Inverse Prob. 28(11), 115005 (2012)

    Article  MathSciNet  Google Scholar 

  17. Lalush, D., Tsui, B.: Simulation evaluation of Gibbs prior distributions for use in maximum a posteriori SPECT reconstructions. IEEE Trans. Med. Imaging 11, 267–275 (1992)

    Article  Google Scholar 

  18. Lebed, E., Herrmann, F.J.: A hitchhiker’s guide to the galaxy of transform-domain sparsification. In SEG Technical Program Expanded Abstracts, SEG, 27 (2008)

  19. Li, Q., Shen, L., Xu, Y., Zhang, N.: Multi-step fixed-point proximity algorithms for solving a class of optimization problems arising from image processing. Adv. Comput. Math. 41(2), 387–422 (2015)

    Article  MathSciNet  Google Scholar 

  20. Lin, T.T.Y., Herrmann, F.J.: Compressed wavefield extrapolation. Geophysics, 72, SM77-SM93 (2007)

  21. Lin, T.T.Y., Lebed, E., Erlangga, Y., Herrmann, F.J.: Interpolating solutions of the Helmholtz equation with compressed sensing. In SEG Technical Program Expanded Abstracts, SEG 27, 2122–2126 (2008)

    Google Scholar 

  22. Lustig, M., Donoho, D., Pauly, J.M.: Sparse MRI: The application of compressed sensing for rapid MR imaging, Magnetic Resonance in Medicine: An Official Journal of the International Society for. Magn. Reson. Med. 58, 1182–1195 (2007)

    Article  Google Scholar 

  23. Lysmer, J., Drake, L.A.: A finite element method for seismology. Methods Comput. Phys. 11, 181–216 (1972)

    Google Scholar 

  24. Marfurt, K.J.: Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations. Geophysics 49, 533–549 (1984)

    Article  Google Scholar 

  25. Micchelli, C.A., Shen, L., Xu, Y.: Proximity algorithms for image models: Denoising. Inverse Prob. 27, 45009–45038 (2011)

    Article  MathSciNet  Google Scholar 

  26. Micchelli, C.A., Shen, L., Xu, Y., Zeng, X.: Proximity algorithms for the \(l_1\)/TV image denosing models. Adv. Comput. Math. 38, 401–426 (2013)

    Article  MathSciNet  Google Scholar 

  27. Moreau, J.-J.: Fonctions convexes duales et points proximaux dans un espace hilbertien. C.R. Acad. Sci. Paris Sér. A Math. 255, 1897–2899 (1962)

    MathSciNet  MATH  Google Scholar 

  28. Riyanti, C.D., Kononov, A., Erlangga, Y.A., Vuik, C., Oosterlee, C.W., Plessix, R.E., Mulder, W.A.: A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation. J. Comput. Phys. 224, 431–448 (2007)

    Article  MathSciNet  Google Scholar 

  29. Shen, L., Xu, Y., Zeng, X.: Wavelet inpainting with the \(l_0\) sparse regularization. Appl. Comput. Harmon. Anal. 41, 26–53 (2016)

    Article  MathSciNet  Google Scholar 

  30. Stuart, A.M.: Inverse problems: A Bayesian perspective. Acta Numer 19, 451–559 (2010)

    Article  MathSciNet  Google Scholar 

  31. Wu, T., Shen, L., Xu, Y.: Fixed-point proximity algorithms solving an incomplete Fourier transform model for seismic wavefield modeling. J. Comput. Appl. Math. 385, 113208 (2021)

    Article  MathSciNet  Google Scholar 

  32. Xu, Y.: Sparse regularization with the \(\ell _0\) norm. Analysis and Applications, accepted (2022)

  33. Zeng, X., Shen, L., Xu, Y.: A convergent fixed-point proximity algorithm accelerated by FISTA for the \(\ell _0\) sparse recovery problem. In: Imaging, Vision and Learning Based on Optimization and PDEs, X.-C. Tai et al. (eds.), Springer, 27-45 (2018)

  34. Zheng, W., Li, S., Krol, A., Schmidtlein, C.R., Zeng, X., Xu, Y.: Sparsity promoting regularization for effective noise suppression in SPECT image reconstruction. Inverse Prob. 35, 115011 (2019)

    Article  MathSciNet  Google Scholar 

Download references

Funding

T. Wu was supported in part by the Natural Science Foundation of Shandong Province of China under grants ZR2021MA049, ZR2020MA031, and Shandong Province Higher Educational Science and Technology Program of China under grant J18KA221. Y. Xu was supported in part by the US National Science Foundation under grant DMS-1912958 and by the US National Institutes of Health under grant R21CA263876.

Author information

Authors and Affiliations

  1. School of Mathematics and Statistics, Shandong Normal University, Jinan, 250358, P.R. China

    Tingting Wu

  2. Department of Mathematics and Statistics, Old Dominion University, Norfolk, VA, 23529, USA

    Yuesheng Xu

Authors
  1. Tingting Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yuesheng Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Yuesheng Xu.

Ethics declarations

Conflict of Interest

The authors have not disclosed any competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, T., Xu, Y. Inverting Incomplete Fourier Transforms by a Sparse Regularization Model and Applications in Seismic Wavefield Modeling. J Sci Comput 92, 48 (2022). https://doi.org/10.1007/s10915-022-01906-8

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s10915-022-01906-8

Keywords

Mathematics Subject Classification

Search

Navigation