Parallel Geothermal Numerical Model with Fractures and Multi-Branch Wells

¹ Université Côte d’Azur, Inria, CNRS, LJAD, UMR 7351 CNRS, team Coffee, Parc Valrose 06108 Nice Cedex 02, France,
laurence.beaude@unice.fr
² CEA Saclay, DEN/DANS/DM2S/STMF/LMEC
thibaud.beltzung@cea.fr
³ Université Côte d’Azur, Inria, CNRS, LJAD, UMR 7351 CNRS, team Coffee, Parc Valrose 06108 Nice Cedex 02, Romania
thibaud.beltzung@cea.fr
⁴ BRGM, 3 avenue Claude-Guillemin, BP 36009, 45060 Orléans Cedex 2, France,
s.lopez@brgm.fr
⁵ Université Côte d’Azur, Inria, CNRS, LJAD, UMR 7351 CNRS, team Coffee, Parc Valrose 06108 Nice Cedex 02, Romania
roland.masson@unice.fr
⁶ BRGM, 3 avenue Claude-Guillemin, BP 36009, 45060 Orléans Cedex 2, France,
f.smai@brgm.fr
⁷ Storengy
jean-frederic.thebault@storengy.com
⁸ Université Côte d’Azur, Inria, CNRS, LJAD, UMR 7351 CNRS, team Coffee, Parc Valrose 06108 Nice Cedex 02, BRGM,
feng.xing@outlook.com

Abstract

To answer the need for an efficient and robust geothermal simulation tool going beyond existing code capabilities in terms of geological and physical complexity, we have started to develop a parallel geothermal simulator based on unstructured meshes. The model takes into account complex geology including fault and fracture networks acting as major heat and mass transfer corridors and complex physics coupling the mass and energy conservations to the thermodynamic equilibrium between the gas and liquid phases. The objective of this Cemracs project was to focus on well modeling which is a key missing ingredient in our current simulator in order to perform realistic geothermal studies both in terms of monitoring and in terms of history matching. The well is discretized by a set of edges of the mesh in order to represent efficiently slanted or multi-branch wells on unstructured meshes. The connection with the 3D matrix and the 2D fracture network at each node of the well is accounted for using Peaceman’s approach. The non-isothermal flow model inside the well is based on the usual single unknown approach assuming the hydrostatic and thermodynamical equilibrium inside the well. The parallelization of the well model is implemented in such a way that the assembly of the Jacobian at each Newton step and the computation of the pressure drops inside the well can be done locally on each process without MPI communications.

Résumé

Afin de dépasser les limites des codes actuels de simulation des systemes géothermiques en matière de complexité géologique et physique, nous avons initié le développement d’un nouveau simulateur d’écoulements géothermiques parallèle à base de maillages non structurés. Le modèle prenden compte une géologie complexe incorporant notamment les réseaux de failles qui jouent un r ôle majeur dans le transport de masse et d’énergie, ainsi qu’une physique complexe couplant les conservations de la masse et de l’énergie à l’équilibre thermodynamique entre les phases liquide et gazeuse. L’objectif de ce projet Cemracs était d’y incorporer un modèle de puits qui constitue un ingrédient essentiel pour réaliser des études géothermiques réalistes à la fois pour la surveillance du réservoir et la reproduction des historiques de production. Le puits est discrétisé par un sous ensemble d’arêtes du maillage de façon a pouvoir représenter efficacement des puits déviés ou multi-branches. La connection avec la matrice 3D et le réseau de failles 2D repose sur des indices de productivité de type Peaceman. Le transport de masse et d’énergie dans le puits se base sur un modèle classique en simulation de réservoir a une inconnue par puits qui suppose l’équilibre hydrostatique et thermodynamique dans le puits. La parallélisation du modèle de puits est réalisée de façon à pouvoir assembler la Jacobienne et à calculer les pertes de charge dans le puits localement sur chaque processus sans nécessiter de communications MPI.