Enhanced linearized reduced-order models for subsurface flow simulation

Authors:
J. He;J. Sætrom;L. J. Durlofsky
Affiliations:
Department of Energy Resources Engineering, Stanford University, Stanford, CA 94305, USA;Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim NO-7034, Norway1Now at Statoil.1;Department of Energy Resources Engineering, Stanford University, Stanford, CA 94305, USA
Venue:
Journal of Computational Physics
Year:
2011

Citing 12
Cited 1

A multiscale finite element method for elliptic problems in composite materials and porous media

Journal of Computational Physics
Accurate multiscale finite element methods for two-phase flow simulations

Journal of Computational Physics
Goal-oriented, model-constrained optimization for reduction of large-scale systems

Journal of Computational Physics
Stabilizing schemes for piecewise-linear reduced order models via projection and weighting functions

Proceedings of the 2007 IEEE/ACM international conference on Computer-aided design
Hybrid Simulations of Reaction-Diffusion Systems in Porous Media

SIAM Journal on Scientific Computing
Guaranteed stable projection-based model reduction for indefinite and unstable linear systems

Proceedings of the 2008 IEEE/ACM International Conference on Computer-Aided Design
Linearized reduced-order models for subsurface flow simulation

Journal of Computational Physics
Stable reduced models for nonlinear descriptor systems through piecewise-linear approximation and projection

IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems
A trajectory piecewise-linear approach to model order reduction and fast simulation of nonlinear circuits and micromachined devices

IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems
Macromodel Generation for BioMEMS Components Using a Stabilized Balanced Truncation Plus Trajectory Piecewise-Linear Approach

IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems
A Piecewise-Linear Moment-Matching Approach to Parameterized Model-Order Reduction for Highly Nonlinear Systems

IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems
Efficient linear circuit analysis by Pade approximation via the Lanczos process

IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems

Reduced-order flow modeling and geological parameterization for ensemble-based data assimilation

Computers & Geosciences

Quantified Score

Hi-index	31.45

Visualization

Abstract

Trajectory piecewise linearization (TPWL) represents a promising approach for constructing reduced-order models. Using TPWL, new solutions are represented in terms of expansions around previously simulated (and saved) solutions. High degrees of efficiency are achieved when the representation is projected into a low-dimensional space using a basis constructed by proper orthogonal decomposition of snapshots generated in a training run. In recent work, a TPWL procedure applicable for two-phase subsurface flow problems was presented. The method was shown to perform well for many cases, such as those with no density differences between phases, though accuracy and robustness were found to degrade in other cases. In this work, these limitations are shown to be related to model accuracy at key locations and model stability. Enhancements addressing both of these issues are introduced. A new TPWL procedure, referred to as local resolution TPWL, enables key grid blocks (such as those containing injection or production wells) to be represented at full resolution; i.e., these blocks are not projected into the low-dimensional space. This leads to high accuracy at selected locations, and will be shown to improve the accuracy of important simulation quantities such as injection and production rates. Next, two techniques for enhancing the stability of the TPWL model are presented. The first approach involves a basis optimization procedure in which the number of columns in the basis matrix is determined to minimize the spectral radius of an appropriately defined amplification matrix. The second procedure incorporates a basis matrix constructed using snapshots from a simulation with equal phase densities. Both approaches are compatible with the local resolution procedure. Results for a series of test cases demonstrate the accuracy and stability provided by the new treatments. Finally, the TPWL model is used as a surrogate in a direct search optimization algorithm, and comparison with results using the full-order model demonstrates the efficacy of the enhanced TPWL procedures for this application.