Efficient full waveform inversion using the excitation representation of the source wavefield

Kalita, Mahesh; Alkhalifah, Tariq Ali

Efficient full waveform inversion using the excitation representation of the source wavefield

Files

document.pdf (5.2 MB)

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

2017-05-17

Print Publication Date

2017-09

Date

2017-05-17

Abstract

Full waveform inversion (FWI) is an iterative method of data-fitting, aiming at high-resolution recovery of the unknown model parameters. However, its conventional implementation is a cumbersome process, requiring a long computational time and large memory space/disk storage. One of the reasons for this computational limitation is the gradient calculation step. Based on the adjoint state method, it involves the temporal cross-correlation of the forward propagated sourcewavefield with the backward propagated residuals, inwhichwe usually need to store the source wavefield, or include an extra extrapolation step to propagate the source wavefield from its storage at the boundary. We propose, alternatively, an amplitude excitation gradient calculation based on the excitation imaging condition concept that represents the source wavefield history by a single, specifically the most energetic arrival. An excitation based Born modelling allows us to derive the adjoint operation. In this case, the source wavelet is injected by a cross-correlation step applied to the data residual directly. Representing the source wavefield through the excitation amplitude and time, we reduce the large requirements for both storage and the computational time. We demonstrate the application of this approach on a two-layer model with an anomaly, the Marmousi II model and a marine data set acquired by CGG.

Citation

Kalita M, Alkhalifah T (2017) Efficient full waveform inversion using the excitation representation of the source wavefield. Geophysical Journal International 210: 1581–1594. Available: http://dx.doi.org/10.1093/gji/ggx214.

Acknowledgements

For computer time, this research used the resources of the Super-computing Laboratory and IT Research Computing at King Abdullah-University of Science&Technology (KAUST) in Thuwal, Saudi Arabia. The real data shown in this study are proprietary to and provided courtesy of CGG. The well-log information is provided by Geoscience Australia. We are grateful to KAUST for financial support and all the members of seismic wave analysis group (SWAG), especially Christos Tzivanakis, Zedong Wu, Juwon Oh, Yunseok Choi and Vladimir Kazei for their fruitful discussions. In addition, we are especially thankful to Zedong Wu for his assistance in the real data example.