Model-Order Reduction Applied To Coupled Flow and Geomechanics Conference Paper uri icon


  • 2018 European Association of Geoscientists and Engineers EAGE. All rights reserved. Coupled flow and geomechanics computations are highly complex and require solving large nonlinear systems. Such simulations are intense for both runtime and memory standpoint, which strongly hints at employing Model-Order Reduction (MOR) techniques to speed them up. Different types of Reduced-Order Models (ROM) have been proposed to alleviate this computational burden. MOR approaches rely on projection. We first execute a computationally expensive "offline" stage, during which we carefully study the full-order model (FOM). Upon creating a ROM basis, we then perform the cheap "online" stage. Our reduction strategy estimates a ROM using Proper Orthogonal Decomposition (POD). We determine a family of solutions of the problem, for a suitable sample of input conditions, where every single realization is so-called a "snapshot." We then ensemble all snapshots to determine a compressed subspace that spans the solution. Usually, POD employs a fixed reduced subspace of global basis vectors. The usage of global basis is not convenient to tackle problems characterized by different physical regimes, parameter changes, or high-frequency features. Having many snapshots to capture all these variations is unfeasible, which suggests seeking adaptive approaches based on the closest regional basis. We thus develop such strategy based on local POD basis to reduce one-way coupled flow and geomechanics computations. We focus on the mechanics and consider factors such as challenging variable boundary conditions and the role of the heterogeneity. We also assess how to tackle different degrees of freedom, such as the displacements (intercalated and coupled) and pressure, with MOR. It seems to be convenient to treat DOF separately since this reduces the rank of the snapshot matrix. Preliminary 3-D results show significant compression ratios up to 90% for the mechanics part. We formally compare FOM and ROM and provide time data to demonstrate the speedup of the procedure. Examples focus on linear and nonlinear poroelasticity. We employ continuous Galerkin finite elements for most of the discretizations.

name of conference

  • ECMOR XVI - 16th European Conference on the Mathematics of Oil Recovery

published proceedings

  • Conference Proceedings

author list (cited authors)

  • Florez, H., & Gildin, E.

citation count

  • 8

complete list of authors

  • Florez, H||Gildin, E

publication date

  • September 2018