Open Access

ARTICLE

Multiscale Hybrid-Mixed Finite Element Method for Flow Simulation in Fractured Porous Media

Philippe Devloo¹, Wenchao Teng², Chen-Song Zhang^3,∗

1 FEC, Universidade Estadual de Campinas, Brazil.
2 Tracy Energy Technologies, Beijing, China.
3 NCMIS & LSEC, Academy of Mathematics and System Sciences, CAS, Beijing, China.
^∗Corresponding Author: Chensong Zhang. Email: zhangcs@lsec.cc.ac.cn.

(This article belongs to the Special Issue: Recent Developments of Immersed Methods for Fluid-structure Interactions)

Computer Modeling in Engineering & Sciences 2019, 119(1), 145-163. https://doi.org/10.32604/cmes.2019.04812

Abstract

The multiscale hybrid-mixed (MHM) method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures. The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements. The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux. A zero-dimensional pressure (point element) is used to express conservation of mass where fractures intersect. The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations. A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh. Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.

---

Keywords

---