Flux-corrected transport for full-waveform inversion

Type
Article

Authors
Kalita, Mahesh
Alkhalifah, Tariq Ali

KAUST Department
Earth Science and Engineering Program
Physical Science and Engineering (PSE) Division
Seismic Wave Analysis Group

Online Publication Date
2019-02-25

Print Publication Date
2019-06-01

Date
2019-02-25

Abstract
Conventional full-waveform inversion (FWI) often fails to retrieve the unknown model parameters from noisy seismic data. A successful FWI implementation usually requires to follow a multistage recovery approach, starting from the retrieval of the lower model wavenumbers (tomography) to those with the higher resolutions (migration). Here, we propose a new method based on the flux-corrected transport (FCT) technique often used in computational fluid dynamics for the removal of instabilities in a shock profile. FCT involves three finite-difference steps: a transport, a diffusion and an antidiffusion process. This third step involves non-linear operators such as maximum and minimum, which are non-differentiable in a classic sense. However, since the seismic source wavelet and the corresponding wavefield are relatively smooth and continuous in nature without any strong ripples like shock waves, we exclude the non-linear step from FCT, which allows us to evaluate the novel FWI gradient efficiently. As a result, we achieve a converging FWI model by gradually reducing the diffusive flux-correction. We demonstrate the versatility of FCT-based FWI on a noisy synthetic data set from the Marmousi II model and a marine field data set from offshore Australia.

Citation
Kalita, M., & Alkhalifah, T. (2019). Flux-corrected transport for full-waveform inversion. Geophysical Journal International, 217(3), 2147–2164. doi:10.1093/gji/ggz105

Acknowledgements
The research reported in this publication was supported by funding from King Abdullah University of Science and Technology (KAUST). We would like to thank KAUST for its support, as well as, all members of the Seismic Wave Analysis Group (SWAG) for the fruitful discussions. We also thank Dr. David Ketcheson from Applied Mathematics and Computational Science at KAUST for his valuable suggestions. For computer time, this research used the resources of the Supercomputing Laboratory and IT Research Computing at KAUST. In addition, this research used the resources of the Core Labs of KAUST. We acknowledge Thomas Theussl from Visualization Core Laboratory at KAUST, for his assistance in displaying some of the results reported in this work. The real data shown in this study are proprietary and provided by courtesy of CGG.

Publisher
Oxford University Press (OUP)

Journal
Geophysical Journal International

DOI
10.1093/gji/ggz105

Additional Links
https://academic.oup.com/gji/article/217/3/2147/5364552

Permanent link to this record