Efficient characterization of uncertain model parameters with a reduced-order ensemble Kalman filter

Abstract

Spatially variable model parameters are often highly uncertain and di fficult to observe. This has prompted the widespread use of Bayesian characterization methods that can infer parameter values from measurements of related variables, while explicitly accounting for uncertainty. Ensemble versions of Bayesian characterization are particularly convenient when uncertain variables have complex spatial structures that do not conform to Gaussian descriptions. However, ensemble methods can be time-consuming for high-dimensional problems. This paper describes a reduced-order approach to ensemble characterization that is particularly well-suited for subsurface flow and transport problems. It uses a truncated discrete cosine transform (DCT) to reduce the dimensionality of spatially variable time-invariant model parameters and a nonlinear extension of principle orthogonal decomposition (POD) to reduce the dimensionality of dynamic model states. The resulting nonlinear reduced-order model can be included in the forecast step of a reduced-order ensemble Kalman fi lter. These concepts are illustrated in a subsurface solute transport problem using ensembles produced by full and reduced-order order models. These ensembles are very similar when there are no measurement updates. When the forecast ensemble is recursively updated with measurements the reduced-order Kalman fi lter does at least as well as the full-order fi lter in characterizing a dynamic solute plume, even though its augmented state dimension is only 2% of the dimension of the full-order state. This substantial increase in effi ciency implies that a reduced-order fi lter with the same ensemble size as its full-order counterpart can give comparable performance for orders of magnitude less computational e ffort or can use a much larger ensemble for the same computational e ffort. The possibility of substantial increases in ensemble size could lead to performance improvements through reductions in sampling error and in the rank of the ensemble null space. Also, a reduced-order model similar to the one described here could be used in ensemble real-time control applications, where it can decrease the eff ort required for both characterization and control.Shell Oil Compan