Multilevel Strategies Improve History Matching of Complex Reservoir Models
Authors: Carpenter Chris 
Published: March 2020

This article, written by JPT Technology Editor Chris Carpenter, contains highlights of paper SPE 193895, “Multilevel Strategies and Geological Parameterizations for History Matching Complex Reservoir Models,” by Yimin Liu, SPE, and Louis J. Durlofsky, SPE, Stanford University, prepared for the 2019 SPE Reservoir Simulation Conference, Galveston, Texas, 10–11 April. The paper has not been peer reviewed.

The complete paper explores the use of multilevel derivative-free optimization for history matching, with model properties described using principal component analysis (PCA) -based parameterization techniques. The parameterizations applied are optimization-based PCA (O-PCA) and convolutional-neural-network-based PCA (CNN-PCA). Mesh adaptive direct search (MADS), a pattern search method that parallelizes naturally, is used for the optimizations required to generate posterior models. The use of PCA-based parameterization reduces considerably the number of variables that must be determined during history matching, but the optimization problem can still be computationally demanding. The multilevel strategy addresses this issue by reducing the number of simulations that must be performed at each MADS iteration. History-matching results demonstrate that substantial uncertainty reduction is achieved in all cases considered and that the multilevel strategy is effective in reducing the number of simulations required.

Overview

Parameterization has two key advantages in history matching:

Far fewer parameters must be determined, which simplifies greatly the minimization in some cases.

Posterior (history-matched) geomodels will be consistent with the geological scenario or training image for which the parameterization was constructed.

PCA is the foundation for many geological-dimension-reduction parameterization methods. PCA is based on the eigen-decomposition of the prior geomodel covariance matrix and, thus, honors only two-point spatial statistics.

History-matching algorithms largely can be classified as either optimization-based or ensemble-based procedures. Optimization-based algorithms generate one history-matched model at a time, while ensemble-based algorithms update a large set of models simultaneously. The resulting ensemble of posterior models is used to quantify uncertainty. In this study, only optimization-based history matching is considered, though the parameterizations applied also can be used with ensemble-based methods.

Optimization-based history-matching algorithms can be classified further into gradient-based and derivative-free optimization (DFO) methods. Gradient-based algorithms often are combined with adjoint procedures that provide the required gradient efficiently. Unlike adjoint-based procedures, DFO algorithms (which are considered in this work) are nonintrusive. They treat the reservoir simulator as a so-called “black box,” and, as a result, they are much easier to implement and can be applied with any simulator.

In this study, MADS, a pattern-search DFO, is applied for history matching. In pattern-search (or stencil-based) algorithms, the number of simulations required at each iteration scales with the number of parameters to be determined. Because the latter typically is much smaller than the number of grid-blocks in the model, geological parameterization leads to large computational savings when a pattern search, or other DFO, is used for history matching. The performance of DFO for history matching with non-Gaussian models, a primary topic in the complete paper, has not been studied extensively.

The authors devise a multilevel strategy, applicable for PCA-based parameterization methods, to further reduce the computational cost of history matching using DFO. A multilevel strategy lends itself naturally to PCA-based parameterizations because the PCA representation can be ordered from the most-important to the least-important principal component. Thus, PCA parameters can be determined sequentially, level by level, in groups, rather than all at once.