  In inverse problems, uncertainty quantification (UQ) deals with a
probabilistic description of the solution nonuniqueness and data noise
sensitivity. Setting seismic imaging into a Bayesian framework allows for a
principled way of studying uncertainty by solving for the model posterior
distribution. Imaging, however, typically constitutes only the first stage of a
sequential workflow, and UQ becomes even more relevant when applied to
subsequent tasks that are highly sensitive to the inversion outcome. In this
paper, we focus on how UQ trickles down to horizon tracking for the
determination of stratigraphic models and investigate its sensitivity with
respect to the imaging result. As such, the main contribution of this work
consists in a data-guided approach to horizon tracking uncertainty analysis.
This work is fundamentally based on a special reparameterization of
reflectivity, known as "deep prior". Feasible models are restricted to the
output of a convolutional neural network with a fixed input, while weights and
biases are Gaussian random variables. Given a deep prior model, the network
parameters are sampled from the posterior distribution via a Markov chain Monte
Carlo method, from which the conditional mean and point-wise standard deviation
of the inferred reflectivities are approximated. For each sample of the
posterior distribution, a reflectivity is generated, and the horizons are
tracked automatically. In this way, uncertainty on model parameters naturally
translates to horizon tracking. As part of the validation for the proposed
approach, we verified that the estimated confidence intervals for the horizon
tracking coincide with geologically complex regions, such as faults.
