A time-domain multisource Bayesian/Markov chain Monte Carlo formulation of time-lapse seismic waveform inversion

Abstract

Seismic waveform inversion has been shown to be well-suited for identification and characterization of time-lapse changes in a reservoir. However, the subtle medium variations associated with enhanced oil recovery and/or [Formula: see text] storage problems give uncertainty quantification within such inverse approaches a heightened importance. To analyze both of these features of the time-lapse inverse problem, we formulate a Bayesian full-waveform inversion (FWI) procedure, based on a Markov chain Monte Carlo (MCMC) algorithm. The formulation uses several existing strategies, such as the use of a double-difference time-lapse FWI (DDFWI), incorporation of time-domain multisource data, and application of a local-updating target-oriented inversion. However, it incorporates these within a stochastic framework, involving computation of model covariance with an adaptive Metropolis algorithm, and a method to estimate data error statistics based on the features of time-lapse difference data that is incorporated. A random walk Metropolis-Hastings MCMC is adopted for optimization. In conventional, i.e., deterministic, DDFWI, inversions are carried out for the baseline and monitoring models; in the MCMC approach, a deterministic FWI procedure is carried out for the baseline model, and the MCMC algorithm is applied in the monitoring inversion stage; and the final time-lapse model is the difference between these. A feasibility study is carried out using synthetic 2D acoustic models and data, including time-lapse model estimation and uncertainty quantification. We compare the MCMC approach with conventional deterministic optimization DDFWI and remark on benefits derived, which appear to justify the expanded complexity and cost of a global approach. In addition to the availability of posterior distributions, which are critical for the assessment of the estimations, we observe that the MCMC approach tends to produce monitoring images with clearer edges and fewer coherent errors.