Cordua, Knud Skou5; Hansen, Thomas Mejer2; Mosegaard, Klaus1
1 Center for Energy Resources Engineering, Center, Technical University of Denmark2 National Space Institute, Technical University of Denmark3 Geomagnetism, National Space Institute, Technical University of Denmark4 Mathematical and Computational Geoscience, National Space Institute, Technical University of Denmark5 Department of Applied Mathematics and Computer Science, Technical University of Denmark
Some multiple-point-based sampling algorithms, such as the snesim algorithm, rely on sequential simulation. The conditional probability distributions that are used for the simulation are based on statistics of multiple-point data events obtained from a training image. During the simulation, data events with zero probability in the training image statistics may occur. This is handled by pruning the set of conditioning data until an event with non-zero probability is found. The resulting probability distribution sampled by such algorithms is a pruned mixture model. The pruning strategy leads to a probability distribution that lacks some of the information provided by the multiple-point statistics from the training image, which reduces the reproducibility of the training image patterns in the outcome realizations. When pruned mixture models are used as prior models for inverse problems, local re-simulations are performed to obtain perturbed realizations. Consequently, these local re-simulations lead to additional pruning in the set of conditioning data, which further deteriorates the pattern reproduction. To mitigate this problem, it is here suggested to combine the pruned mixture model with a frequency matching model. The multiple-point statistics of outcome realizations from this combined model has improved degree of match with the statistics from the training image. An efficient algorithm that samples this combined model is suggested. Finally, a tomographic cross-borehole inverse problem with prior information expressed by the combined (prior) model is used to demonstrate the effect of pattern reproducibility on the resolution of an inverse problem.
Mathematical Geosciences, 2014, Vol 47, p. 317-343