Affordable Access

Publisher Website

A Markov chain Monte Carlo inverse procedure for radioactive contaminant transport in fractured rock: 2. Inverse model

Authors
Journal
Annals of Nuclear Energy
0306-4549
Publisher
Elsevier
Publication Date
Volume
27
Issue
14
Identifiers
DOI: 10.1016/s0306-4549(00)00030-x
Disciplines
  • Computer Science
  • Mathematics

Abstract

Abstract Model calibration for radioactive contaminant transport in fractured rock has been formulated as an inverse problem and a solution proposed using a Bayesian approach. In particular, Markov random field (MRF) models have been proposed as a priori distribution functions for fracture apertures. The fracture reconstruction is based on the a posteriori distribution function simulation using a Markov chain Monte Carlo procedure. A two-level estimation-reconstruction algorithm has been developed. On the first level, the fracture apertures are reconstructed using simulated annealing (SA), with the Metropolis-Hastings dynamics. The MRF parameters are assumed known. On the second level, after several sweeps of the entire field using SA, the MRF parameters are updated using the method of pseudolikelihood. The competing models are finally compared using an asymptotic form of the Bayes factor.

There are no comments yet on this publication. Be the first to share your thoughts.