  Full waveform inversion (FWI) enables us to obtain high-resolution velocity
models of the subsurface. However, estimating the associated uncertainties in
the process is not trivial. Commonly, uncertainty estimation is performed
within the Bayesian framework through sampling algorithms to estimate the
posterior distribution and identify the associated uncertainty. Nevertheless,
such an approach has to deal with complex posterior structures (e.g.,
multimodality), high-dimensional model parameters, and large-scale datasets,
which lead to high computational demands and time-consuming procedures. As a
result, uncertainty analysis is rarely performed, especially at the industrial
scale, and thus, it drives practitioners away from utilizing it for
decision-making. This work proposes a frugal approach to estimate uncertainty
in FWI through the Stein Variational Gradient Descent (SVGD) algorithm by
utilizing a relatively small number of velocity model particles. We warm-start
the SVGD algorithm by perturbing the optimized velocity model obtained from a
deterministic FWI procedure with random field-based perturbations. Such
perturbations cover the scattering (i.e., high wavenumber) and the transmission
(i.e., low wavenumber) components of FWI and, thus, represent the uncertainty
of the FWI holistically. We demonstrate the proposed approach on the Marmousi
model; we have learned that by utilizing a relatively small number of
particles, the uncertainty map presents qualitatively reliable information that
honours the physics of wave propagation at a reasonable cost, allowing for the
potential for industrial-scale applications. Nevertheless, given that
uncertainties are underestimated, we must be careful when incorporating them
into downstream tasks of seismic-driven geological and reservoir modelling.
