ARTICLE

Acoustic full waveform inversion using Discrete Cosine Transform (DCT)

SUMIN KIM1 WOOKEEN CHUNG1,2 JONGHYUN LEE3
Show Less
1 Department of Convergence Study on the Ocean Science and Technology, Ocean Science and Technology (OST) School, Korea Maritime and Ocean University, Busan, South Korea.,
2 Department of Energy and Resource Engineering, Korea Maritime and Ocean University, Busan, South Korea.,
3 Department of Civil and Environmental Engineering, and Water Resources Research Center, University of Hawaii at Manoa, Honolulu, HI 96822, U.S.A.,
JSE 2021, 30(4), 365–380;
Submitted: 9 June 2025 | Revised: 9 June 2025 | Accepted: 9 June 2025 | Published: 9 June 2025
© 2025 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract

Kim, S., Chung, W.K. and Lee, J.H., 2021. Acoustic full waveform inversion using Discrete Cosine Transform (DCT). Journal of Seismic Exploration, 30: 365-380. Full waveform inversion (FWI) has been implemented widely to reconstruct high-quality velocity model in the subsurface. However, due to the advance in geophysics data acquisition with increasing demand for high dimension velocity model estimation, computational costs become prohibitive if not impossible. To alleviate computational burdens, we incorporate discrete cosine transform, one of the most widely used compression techniques in image processing, into FWI while estimation results are still kept comparable to those from the full-model based FWI. The unknown velocity fields are transformed in the DCT domain and only a small number of DCT coefficients are included in FWI to describe common velocity model features without losing much reconstruction accuracy. Generally, DCT coefficients can be chosen specific window (e.g., square window). However, there can be more dominant DCT coefficients out of this specific window. To take more dominant DCT coefficients, we sort the absolute value of DCT coefficients in descending order and determine DCT coefficients with compression ratio. Through the comparison reconstructed velocity models, our proposed method generate more accurate velocity model than case of using square window. We investigate the applicability of our DCT-based FWI method to two numerical examples. It is shown that the proposed method can reduce the computational cost significantly and produce satisfactory results. Through these results, we expect that our FWI method can contribute to enhance computational efficiency for FWI with enormous amount of unknown parameters.

Keywords
acoustic
seismic
full waveform inversion (FWI)
Discrete Cosine Transform (DCT)
References
  1. Baraniuk, R.G. and Steeghs, P., 2017. Compressive sensing: A new approach to seismic
  2. data acquisition. The Leading Edge, 36: 642-645.
  3. Ben-Hadj-Ali, H., Operto, S. and Virieux, J., 2011. An efficient frequency-domain full
  4. waveform inversion method using simultaneous encoded sources. Geophysics,
  5. 76: R109-R124.
  6. Chen, Y., 2017. Fast dictionary learning for noise attenuation of multidimensional
  7. seismic data. Geophys. J. Internat., 209: 21-31.
  8. Dalmau, F.R., Hanzich, M., de la Puente, J. and Gutiérrez, N., 2014. Lossy data
  9. compression with DCT transforms. EAGE Workshop on High Performance
  10. Computing for Upstream, HPC30.
  11. Gao, F., Atle, A. and Williamson, P., 2010. Full waveform inversion using deterministic
  12. source encoding. Expanded Abstr., 80th Ann. Internat. SEG Mtg., Denver:
  13. 1013-1017.
  14. Habashy, T.M., Abubakar, A., Pan, G. and Belani, A., 2011. Source-receiver
  15. compression scheme for full-waveform seismic inversion. Geophysics, 76:
  16. R95-R108.
  17. Jafarpour, B. and McLaughlin, D.B., 2008. History matching with an ensemble Kalman
  18. filter and discrete cosine parameterization. Computat. Geosci., 12: 227-244.
  19. Jain, A.K., 1989. Fundamentals of Digital Image Processing. Prentice Hall, Englewood
  20. Cliffs, NJ.
  21. Krebs, J.R., Anderson, J.E., Hinkley, D., Neelamani, R., Lee, S., Baumstein, A. and
  22. Lacasse, M.D., 2009. Fast full-wavefield seismic inversion using encoded
  23. sources. Geophysics, 74: WCC177-WCC188.
  24. Lailly, P., 1983. The seismic inverse problem as a sequence of before stack migrations:
  25. Conference on Inverse Scattering, Theory and Application, Society for
  26. Industrial and Applied Mathematics, Expanded Abstr.: 206-220.
  27. Lam, E.Y. and Goodman, J.W., 2000. A mathematical analysis of the DCT coefficient
  28. distributions for images. IEEE Transact. Image Process., 9: 1661-1666.
  29. Lee, J., Yoon, H., Kitanidis, P.K., Werth, C.J. and Valocchi, A.J., 2016. Scalable
  30. subsurface inverse modeling of huge data sets with an application to tracer
  31. concentration breakthrough data from magnetic resonance imaging. Water
  32. Resour. Res., 52: 5213-5231.
  33. Li, X., Aravkin, A.Y., van Leeuwen, T. and Herrmann, F.J., 2012. Fast randomized full-
  34. waveform inversion with compressive sensing. Geophysics, 77: A13-A17.
  35. Li, X., Esser, E. and Herrmann, F.J., 2016. Modified Gauss-Newton full-waveform
  36. inversion explained - Why sparsity-promoting updates do matter. Geophysics,
  37. 81: R125-R138.
  38. Lin, Y., Abubakar, A., and Habashy, T. M., 2012. Seismic full-waveform inversion using
  39. truncated wavelet representations. Expanded Abstr., 82nd Ann. Internat. SEG
  40. Mtg., Las Vegas: 1-6..
  41. Martin, G. S., Wiley, R., and Marfurt, K. J., 2006. Marmousi2: An elastic upgrade for
  42. Marmousi. The Leading Edge, 25: 156-166.
  43. Nocedal, J., and Wright, S., 2006. Numerical optimization. Springer Science & Business
  44. Media.
  45. Pratt, R.G., Shin, C. and Hick, G.J., 1998. Gauss-Newton and full Newton methods in
  46. frequency-space seismic waveform inversion. Geophys. J. Internat., 133: 341-
  47. Son, W., Pyun, S., Jang, D., Park, Y. and Shin, C., 2012. A new algorithm adapting
  48. encoded simultaneous-source full waveform inversion to the marine-streamer
  49. data. Expanded Abstr., 82nd Ann. Internat. SEG Mtg., Las Vegas: 1-5.
  50. Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation.
  51. Geophysics, 49: 1259-1266.
  52. Virieux, J. and Operto, S., 2009. An overview of full-waveform inversion in exploration
  53. geophysics. Geophysics, 74: WCCI-WCC26.
  54. Zhou, H. and Li, Q., 2013. Increasing waveform inversion efficiency of GPR data using
  55. compression during reconstruction. J. Appl. Geophys., 99: 109-113.
  56. Zhu, L., Liu, E. and McClellan, J-H., 2015. Seismic data denoising through multiscale
  57. and sparsity-promoting dictionary learning. Geophysics, 80: WD45-WD57.
  58. Zhu, L., Liu, E. and McClellan, J.H., 2017. Sparse-promoting full-waveform inversion
  59. based on online orthonormal dictionary learning. Geophysics, 82: R87-R107.
Share
Back to top
Journal of Seismic Exploration, Electronic ISSN: 0963-0651 Print ISSN: 0963-0651, Published by AccScience Publishing