Data assimilation and uncertainty assessment for complex geological models using a new PCA-based parameterization

Abstract

In this paper, a recently developed parameterization procedure based on principal component analysis (PCA), which is referred to as optimization-based PCA (O-PCA), is generalized for use with a wide range of geological systems. In O-PCA, the mapping between the geological model in the full-order space and the low-dimensional subspace is framed as an optimization problem. The O-PCA optimization involves the use of regularization and bound constraints, which act to extend substantially the ability of PCA to model complex (non-Gaussian) systems. The basis matrix required by O-PCA is formed using a set of prior realizations generated by a geostatistical modeling package. We show that, by varying the form of the O-PCA regularization terms, different types of geological scenarios can be represented. Specific systems considered include binaryfacies, three-facies and bimodal channelized models, and bimodal deltaic fan models. The O-PCA parameterization can be applied to generate random realizations, though our focus here is on its use for data assimilation. For this application, O-PCA is combined with the randomized maximum likelihood (RML) method to provide a subspace RML procedure that can be applied to non-Gaussian models. This approach provides multiple history-matched models, which Electronic supplementary material The online version of this article (enables an estimate of prediction uncertainty. A gradient procedure based on adjoints is used for the minimization required by the subspace RML method. The gradient of the O-PCA mapping is determined analytically or semianalytically, depending on the form of the regularization terms. Results for two-dimensional oil-water systems, for several different geological scenarios, demonstrate that the use of O-PCA and RML enables the generation of posterior reservoir models that honor hard data, retain the large-scale connectivity features of the geological system, match historical production data, and provide an estimate of prediction uncertainty. MATLAB code for the O-PCA procedure, along with examples for three-facies and bimodal models, is included as online Supplementary Material. Mathematics Subject Classification (2010)15A04 · 49N45 · 60G60 · 78M34 · 90-08 · 94A08 · 90C30

A Deep-Learning-Based Geological Parameterization for History Matching Complex Models

“…In such systems, the value of m is either 0 or 1, where, e.g., 0 represents mud (shale) facies and 1 represents sand (channel) facies. The detailed formulation for O-PCA for bimodal systems is presented in Vo and Durlofsky (2015).…”

Big Data Management for Cloud-Enabled Geological Information Services

“…Effective data mining also could reduce manual intervention during information processing and make use of methods and tools of big data intelligent analysis [35,36]. Recently, there has been a growing interest in the geological big data mining through the use of some novel computational intelligent methods-for example, rough set [37] and fuzzy aggregation [38].…”

Flow-based dissimilarity measures for reservoir models: a spatial-temporal tensor approach

“…The numerical model of one realization has 3600 grid blocks of size 50 m × 50 m. A description of the physical parameters, wells configuration, and a link to the data files can be found in [39]. The sequential solvers of MRST, see [24], have been used to solve the pressure and saturation equations and the production has been simulated for a period of 15 years, time step of 5 days.…”