Acoustic full waveform inversion using Discrete Cosine Transform (DCT)

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