Une méthode mixte multi-échelles pour un simulateur de réservoir biphasé

A multiscale hybrid mixed finite element method is presented in this paper to solve two-phase flow equations on heterogeneous media under the effect of gravitational segregation. It is designed to cope with the complex geometry and inherent multiscale nature of the rocks, leading to stable and accurate multi-physics reservoir simulations. This multiscale approach makes use of coarse scale fluxes between subregions (macro domains) that allow to reduce substantially the dominant computational costs associated with the flux/pressure kernel embedded in the numerical model. As such, larger scale problems can be approximated in a reasonable computational time. Dividing the problems into macro domains leads to a hierarchy of meshes and approximation spaces, allowing the efficient use of static condensation and parallel computation strategies. The method documented in this work utilizes discretizations based on a general domain partition formed by poly-hedral subregions. The normal flux between these subregions is associated with a finite dimensional trace space. The global system to be solved for the fluxes and pressures is expressed only in terms of the trace variables and of a piecewise constant pressure associated with each subregion. The fine scale features are resolved by mixed finite element approximations using fine flux and pressure representations inside each subregion, and the trace variable (i.e. normal flux) as Neumann boundary conditions. This property implies that the flux approximation is globally H(div)-conforming, and, as in classical mixed formulations, local mass conservation is observed at the micro-scale elements inside the subregions, an essential property for flows in heterogeneous media

A novel block non-symmetric preconditioner for mixed-hybrid finite-element-based flow simulations

In this work we propose a novel block preconditioner, labelled Explicit Decoupling Factor Approximation (EDFA), to accelerate the convergence of Krylov subspace solvers used to address the sequence of non-symmetric systems of linear equations originating from flow simulations in porous media. The flow model is discretized blending the Mixed Hybrid Finite Element (MHFE) method for Darcy's equation with the Finite Volume (FV) scheme for the mass conservation. The EDFA preconditioner is characterized by two features: the exploitation of the system matrix decoupling factors to recast the Schur complement and their inexact fully-parallel computation by means of restriction operators. We introduce two adaptive techniques aimed at building the restriction operators according to the properties of the system at hand. The proposed block preconditioner has been tested through an extensive experimentation on both synthetic and real-case applications, pointing out its robustness and computational efficiency