Local uncertainty quantification for 3D time-domain full waveform inversion with ensemble Kalman filters: application to a North sea OBC dataset

Alexandre Hoffmann,Romain Brossier,Ludovic Métivier,Alizia Tarayoun
DOI: https://doi.org/10.1093/gji/ggae114
IF: 3.352
2024-03-21
Geophysical Journal International
Abstract:Summary Full waveform inversion has emerged as the state-of-the art high resolution seismic imaging technique, both in seismology for global and regional scale imaging and in the industry for exploration purposes. While gaining in popularity, full waveform inversion, at an operational level, remains a heavy computational process involving the repeated solution of large-scale 3D wave propagation problems. For this reason it is a common practice to focus the interpretation of the results on the final estimated model. This is forgetting full waveform inversion is an ill-posed inverse problem in a high dimensional space for which the solution is intrinsically non-unique. This is the reason why being able to qualify and quantify the uncertainty attached to a model estimated by full waveform inversion is key. To this end, we propose to extend at an operational level the concepts introduced in a previous study related to the coupling between ensemble Kalman filters and full waveform inversion. These concepts had been developed for 2D frequency-domain full waveform inversion. We extend it here to the case of 3D time-domain full waveform inversion, relying on a source subsampling strategy to assimilate progressively the data within the Kalman filter. We apply our strategy to an ocean bottom cable field dataset from the North Sea to illustrate its feasibility. We explore the convergence of the filter in terms of number of elements, and extract variance and covariance information showing which part of the model are well constrained and which are not. Analyzing the variance helps to gain insight on how well the final estimated model is constrained by the whole full waveform inversion workflow. The variance maps appears as the superposition of a smooth trend related to the geometrical spreading and a high resolution trend related to reflectors. Mapping lines of the covariance (or correlation matrix) to the model space helps to gain insight on the local resolution. Through a wave propagation analysis, we are also able to relate variance peaks in the model space to variance peaks in the data space. Compared to other posterior-covariance approximation scheme, our combination between Ensemble Kalman filter and full waveform inversion is intrinsically scalable, making it a good candidate for exploiting the recent exascale high performance computing machines.
geochemistry & geophysics
What problem does this paper attempt to address?