A simple algorithm to generate small geostatistical ensembles for subsurface flow simulation

This note describes the theory behind a simple Matlab program to generate ensembles of geological realizations of relatively small-scale reservoir models (permeability fields) for subsurface flow simulation. The algorithm makes use of Principal Component Analysis (PCA) to parameterize the spatial covariance of a reference image, and a Bayesian...