High energy geothermal simulations in faulted reservoirs with code ComPASS
- Team: this project started in 2014 in collaboration between
- BRGM: Simon Lopez and Farid Smai
- Joint INRIA-LJAD team Coffee: Feng Xing, Laurence Beaude, Konstantin Brenner and Roland Masson
- The ANR project CHARMS starting in october 2016 for 4 years will strenghen the project and enlarge
it to the additional partners Storengy (Delphine Patriarche, Vincent Seignole, Jean-Frédéric Thébault), La Maison de la Simulation (Michel Kern, Julien Bigot, Pascal Tremblin, Thomas Dufaut) and the Laboratoire Jacques Louis Lions (Cindy Guichard, Robert Eymard).
Objectives: our objectives are to develop efficient parallel numerical methods to simulate gas liquid compositional thermal
flow in high energy geothermal reservoirs.
The spatial discretization is adapted to unstructured polyhedral meshes and takes into account discrete fracture
or fault networks represented as 2D surfaces connected to the surrounding 3D heterogeneous anisotropic
porous medium (the matrix).
The difficulties result from the highly contrasted spatial and temporal scales
and petrophysical properties between the 3D matrix and the fault network.
The second difficulty is to account for the strong couplings and high nonlinearities induced by the large range of pressure and temperature and by the phase appearance and disapearance.
These numerical methods and physical models are implemented in code ComPASS (Computing Architecture to Speed up Simulation).
- Brief history of the code ComPASS
- Current version of ComPASS: it has the following features
- Computer science
- Version management of the code using Inria Forge and Git
- SPMD Paradigm
- Fortran 2003 + C/C++ + MPI
- Parallelism based on a partitioning of the mesh with Metis
- One layer of ghost cells
- Connected to the solver libraries Petsc + Hypre + Trilinos
- Visualization ouputs using parallel vtk format
- Checkpointing using HDF5
- Numerical methods
- Multiphase compositional thermal Darcy flow model
- 2D discrete fracture or fault network coupled with the surrounding 3D matrix