Abstract: A method is described for seismic depth uncertainty analysis including receiving wavelet basis functions and cutoff thresholds and randomly perturbing wavelet coefficients in reduced wavelet space based on the wavelet basis functions and the cutoff thresholds to generate a plurality of random wavelet fields; receiving a reference model in a depth domain; transforming the plurality of random wavelet fields to the depth domain and combining them with the reference model to form candidate models; performing a hierarchical Bayesian modeling with Markov Chain Monte Carlo (MCMC) sampling methods using the candidate models as input to generate a plurality of realizations; and computing statistics of the plurality of realizations to estimate depth uncertainty. The method may be executed by a computer system.

Abstract: A method is described for estimating depth uncertainty including receiving seismic data, a reference model, and trial model realizations; generating realization gathers from the trial model realizations; generating reference gathers from the reference model; determining a reference data fit based on the reference gathers and a data fit for trial models based on the realization gathers; selecting refined models from the trial model realizations based on the reference data fit, the data fit for trial models, and a data fit tolerance criterion; and calculating depth uncertainty based on statistics of the refined models. The method may be executed by a computer system.

Abstract: A method is described for stochastic modeling of seismic velocity and anisotropic parameters, including receiving 3D bounds of normal moveout velocity (Vnmo) and anisotropic parameter ?; modeling 3D bounds for vertical velocity V and anisotropic parameter ? based on the 3D bounds of Vnmo and ?; generating 3D model realizations of V, ?, and ? within the 3D bounds; and testing detectability of each of the 3D model realizations to create a detectable subset of model realizations wherein the detectability identifies which 3D model realizations will produce images with flat migrated gathers. The method may be executed by a computer system.

Abstract: A method is described for seismic depth uncertainty analysis including receiving wavelet basis functions and cutoff thresholds and randomly perturbing wavelet coefficients in reduced wavelet space based on the wavelet basis functions and the cutoff thresholds to generate a plurality of random wavelet fields; receiving a reference model in a depth domain; transforming the plurality of random wavelet fields to the depth domain and combining them with the reference model to form candidate models; performing a hierarchical Bayesian modeling with Markov Chain Monte Carlo (MCMC) sampling methods using the candidate models as input to generate a plurality of realizations; and computing statistics of the plurality of realizations to estimate depth uncertainty. The method may be executed by a computer system.