Les géostructures énergétiques sont des ouvrages de génie civil qui intègrent un circuit de fluide caloporteur. Dans cette étude l’installation de circuits géothermiques dans des remblais est envisagé afin d’injecter ou d’extraire de la chaleur. Lors de la construction de ces remblais, deux paramètres sont à contrôler : la teneur en eau w (%) et la masse volumique sèche rd (Mg/m3). Ces deux paramètres ainsi que la nature du sol, le chargement appliqué et la variation cyclique de la température pourraient influencer la capacité de stockage et la stabilité du système. Après une partie synthétisant l’état de l’art sur ces thématiques, cette thèse comprend une partie expérimentale puis une partie numérique appliquée à un cas test de stockage. Dans la partie expérimentale, l’influence de la variation de température induite par l’échange thermique sur les propriétés thermiques et mécaniques du sol est étudiée. Pour cela des essais sont réalisés en laboratoire sur un limon du bassin parisien compacté et soumis à une gamme de variation de température de 20° à 50°C. Les effets de ces variations sur les paramètres mécaniques et les paramètres thermiques sont mesurés. Les résultats obtenus montrent une évolution des paramètres thermiques du matériau du côté sec de la courbe de compactage, particulièrement visible pour les températures importantes. La réversibilité de l’effet du chauffage est obtenue après plusieurs cycles de chauffage-refroidissement. Les propriétés mécaniques mesurées sont la résistance à la compression simple et les paramètres pressiométriques. Pour cela, des essais en laboratoire, à l’échelle métrique, sont réalisés. Un ramollissement du sol sous l’effet du chauffage et une rigidification sous l’effet des cycles thermiques sont constatés. Pour optimiser le système du stockage, une modélisation thermo-hydrique d’un remblai modèle réalisé avec le limon de Plaisir compacté légèrement du côté humide de l’optimum Proctor est effectuée. La comparaison des cas tests a permis de déterminer l’entraxe qui minimise l’interaction entre les sondes thermiques, ainsi que la distance optimale entre les dernières sondes et le pied de remblai. Les scenarii d’injection/extraction étudiés ont montré que le stockage de chaleur dans le remblai est faisable en choisissant un programme thermique adéquat pour l’exploitation du stockage. La mise en place de plusieurs lits de sondes thermiques et d’une isolation thermique sur la partie supérieure du remblai de stockage permet d’augmenter la capacité du remblai à accumuler de la chaleur dans le noyau à la fin de la période de relaxation. La simulation de ce scénario sur 10 ans avec la prise en compte des cycles de température, montre un réchauffement du remblai sur plusieurs années jusqu’à atteindre une évolution équilibrée au bout de la 7 ème année. Arrivé à cet état d’équilibre, l’évolution des températures est reproductible d’une année sur l’autre / Energy geostructures are civil engineering structures equipped with energy exchanger elements in order to store heat seasonally. The aim of this study is the use of compacted soil to store energy through installation of horizontal exchangers in an embankment. During the soil compaction, two parameters were controlled: the water content w and the dry density ρd. This parameters as well as the soil nature, the mechanical load path and the cyclic temperature variation may affect the storage capacity and the system stability. This work includes an experimental part and a modeling part. In the experimental part, the impact of the temperature variation is studied on the thermal and mechanical properties of the compacted soil. The Plaisir loam (PL) extracted from the Paris region was investigated through laboratory tests at a temperatures range of 20 to 50°C. The results showed that the thermal properties of the compacted soils increased on the dry side of the compaction curve. This evolution was clearly confirmed for higher temperatures. The application of cyclic temperature variations showed reversible evolutions in the thermal properties after one cycle. The results of unconfined compressive tests and pressuremeter tests showed that heating induces a softening of the material whereas several temperature cycles induce a stiffening of the material. Thereafter, a coupled thermo-hydraulic modeling of an embankment made of the compacted Plaisir loam is performed in order to optimize the storage system. The comparison of different modelling results, fixed the temperature sensors spacing that minimizes the interactions between them, and the optimal distance between the last sensors and the bottom of the slope. It is shown that if an appropriate thermal program is chosen, the heat storage in the embankment could be possible. A better efficiency of the storage capacity can be reached by introducing 3 temperature sensors rows in the storage and by covering the storage with a thermal insulation. The simulation of this scenario over 10 years including temperature cycles shows a heating of the embankment for several years, until to reach an equilibrium state after the 7th year

Cette thèse cible l'étude de la structure thermique de la croûte supérieure (<10km) dans les arcs magmatiques continentaux, et son influence sur l'enregistrement thermochronologique de leur exhumation et de leur évolution topographique. Nous portons notre regard sur deux chaînes de montagne appartenant aux Cordillères Américaines : Les Cascades Nord (USA) et la zone de faille Motagua (Guatemala). L'approche utilisée est axée sur l'utilisation de la thermochronologie (U-Th-Sm)/He sur apatite et zircon, couplée avec la modélisation numérique de la structure thermique de la croûte. Nous mettons en évidence la variabilité à la fois spatiale et temporelle du gradient géothermique, et attirons l'attention du lecteur sur l'importance de prendre en compte la multitude des processus géologiques perturbant la structure thermique dans les chaînes de type cordillère, c'est à dire formées lors de la subduction océanique sous un continent. / This thesis focuses on the influence of the dynamic thermal structure of the upper crust (<10km) on the thermochronologic record of the exhumational and topographic history of magmatic continental arcs. Two mountain belts from the American Cordillera are studied: the North Cascades (USA) and the Motagua fault zone (Guatemala). I use a combined approach coupling apatite and zircon (U-Th-Sm)/He thermochronology and thermo-kinematic numerical modelling. This study highlights the temporal and spatial variability of the geothermal gradient and the importance to take into account the different geological processes that perturb the thermal structure of Cordilleran-type mountain belts (i.e. mountain belts related to oceanic subduction underneath a continent).

Ce travail s'inscrit dans le cadre de la recherche des solutions d’économie d'énergie du bâtiment tout en utilisant des sources naturelles et renouvelables (Energie solaire pour le chauffage et la géothermie pour le rafraîchissement). Il est nul besoin de rappeler que l'Algérie dispose d'un potentiel énergétique hélio géothermique important. Dans ce travail de thèse on s'intéresse particulièrement à l'étude du confort thermique (hiver et été) dans l'habitat alimenté par un plancher hydraulique réversible. Pour mener à bien cette étude, nous disposons d'un gisement solaire important d'une part et d'autre part d'un dispositif expérimental à échelle réelle. Une pièce munie d'un plancher hydraulique réversible (chauffant ou rafraichissant) est instrumentée. Une citerne de stockage enfuie à quelques mètres de la surface du sol afin de bénéficier du rafraichissement géothermique. Un service d'asservissement permettant la régulation du système en fonctionnement mode chauffage ou mode rafraichissement. Plusieurs sondes de mesures reliées à une station d'acquisition qui est reliée à un ordinateur permettent le suivi des évolutions de températures. La modélisation de la structure de l'enveloppe de la cellule et l'évolution de la température de l'air intérieur et celle des parois sont réalisées sous le logiciel TRNSYS. A l'aide des résultats obtenus par TRNSYS, logiciel FLUENT nous a permis de modéliser la tache solaire et son influence sur le plancher chauffant sous les conditions climatiques de la ville d'Oran.Après validation, la simulation numérique est utilisée pour étudier le comportement thermique de la cellule, les performances énergétiques du plancher réversible et le calcul des économies d'énergie que l'on pourrait réaliser avec de tels systèmes. / This work is part of the search for energy saving solutions in the building industry while using natural and renewable sources, such as solar energy for heating and geothermal energy for refreshment. There is no need to recall that Algeria has a very large geothermal gravitational energy potential in view of its geographical position.In this thesis work, we are particularly interested in the study of thermal comfort in the case of a housing powered with a reversible hydraulic floor (heating and cooling).To carry out this study, as we can see Algeria have an important solar field and on the other hand we use an experimental system representing a real scale local. To do this, a room with a reversible hydraulic floor (heated or refreshing) is instrumented. A storage tank buried in the ground at few meters from the ground surface is used for thegeothermal refreshment during the warm periods. A service system allows us to regulate the system in heating or cooling mode. Several measuring probes used are connected to an acquisition station which is connected to a computer for monitoring of temperature évolutions.The modeling of the structure of the cell envelope is carried out under the TRNSYS software. With this, we have access to evolutions of the temperatures of the indoor air and to that of the walls. These results obtained by TRNSYS are used in a second step as input data for the FLUENT software. This allows us to model the solar spot and its influence on the heating floor under the climatic conditions of the city of Oran.After validation, numerical simulation is used to study the thermal behavior of the cell, the energy performance of the reversible floor and the calculation of the energy savings that could be achieved with such systems.

La caldera du volcan Rinjani contient un lac d’un volume de 1 km³ qui est probablement le plus grand lac volcanique au monde présentant une anomalie thermique nette. Ce lac présente une composition neutre chlorure sulfate bicarbonate inhabituelle pour les lacs volcaniques. Sa TDS (2600 mg/l) et conductivité (3500µs/cm) élevées indiquent un apport de fluides hydrothermaux très important. Enfin, son alcalinité élevée (520 mg/l), indique un apport important de dioxyde de carbone dans le lac.<p>Les sources thermales situées autour du Gunung Baru (cône volcanique situé dans la caldera) ont une composition chimique en éléments majeurs et une composition isotopique proche de celles du lac volcanique indiquant qu’elles sont essentiellement le résultat du recyclage du lac par le système hydrothermal. Les variations de compositions entre les différentes sources ont permis de montrer que leurs compositions est le résultat du mélange entre un fluide hydrothermal profond de composition neutre chlorure, dont la température a été estimée à 270°C, et d’un fluide plus superficiel riche en magnésium et en sulfate.<p>Le flux de dioxyde de carbone à la surface du lac a été estimé à l’aide de la méthode de la chambre d’accumulation et par calcul à environ 2300 t/j, ce qui représente un apport significatif de gaz. Cependant, comme le lac présente une structure polymictique, le risque d’accumulation de dioxyde de carbone en profondeur et donc d’éruption limnique peut être exclus.<p>Pour la première fois dans cette thèse, le modèle d’estimation des flux thermiques émis par les lacs volcaniques mis au point par Stevenson (1992) a été contraint par des mesures des paramètres météorologiques mesurés en continu, ce qui a permis de valider le modèle. De plus, nous avons pu montrer que l’essentiel des variations de températures des lacs volcaniques est dû à des variations météorologiques. En utilisant le flux thermique plutôt que la température, il est dès lors possible d’avoir accès à des variations de l’activité volcanique.<p>Le flux thermique estimé pour le lac du Rinjani est de 1700 MW, ce qui représente le flux le plus élevé jamais mesuré sur un lac volcanique aérien. Ce flux thermique est aussi plus élevé que le flux thermique mesuré sur des lacs de lave à 800°C. Ce paradoxe apparent s’explique par la plus grande dimension des lacs volcaniques, la capacité calorifique de l’eau quatre fois plus importante que celle du magma et la viscosité de l’eau 1 million de fois inférieure, ce qui fait de l’eau un excellent fluide caloporteur pour transporter les calories vers la surface.<p> / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Cette thèse s’est déroulée dans le cadre du projet FONGEOSEC, qui vise à développer la filière de la géothermie profonde en France avec la conception d’un démonstrateur d’une centrale de production d’énergie géothermique dans le bassin Rhénan. Ce projet est piloté par Fonroche Géothermie, qui gère un consortium de plus de dix acteurs du milieu académique et industriel. Le financement du projet est réalisé avec participation de l’Agence de l’Environnement et de la Maîtrise de l’Énergie (ADEME). Ainsi, les travaux exposés dans ce document se sont intéressés à la caractérisation thermodynamique des fluides géothermaux (saumures chaudes contenant des gaz dissous) de la région ciblée par le projet. Il est nécessaire de déterminer la solubilité de chacun des gaz dissous dans ces saumures aux conditions de pression, température et salinité de l’exploitation géothermique. Des modèles thermodynamiques de prévision des équilibres entre phases liquide et vapeur peuvent être utilisés pour estimer ces solubilités. Néanmoins, en absence des mesures expérimentales dans les conditions de pression, température et salinité d’intérêt, pour effectuer la régression de paramètres de ces modèles, il sera difficile d’obtenir avec précision ces solubilités à partir de simulations. Ainsi, cette thèse est centrée sur l’étude expérimentale des solubilités des gaz dans des saumures représentatives des fluides du bassin Rhénan. La gamme de pression de FONGEOSEC va de 6.0 MPa à 40.0 MPa pour des températures de 333.15 K et 453.15 K. Le dispositif expérimental utilisé dans cette thèse fonctionne dans ces conditions. Les gaz dissous dans les saumures visées par le projet sont constitués essentiellement de dioxyde de carbone (CO2), puis d’azote (N2) et enfin de méthane (CH4) en plus faibles quantités. Les sels dissous dans ces fluides sont surtout du chlorure de sodium et du chlorure de calcium, à molalité de 1.2 mol NaCl-0.2 mol CaCl2.Kg H2O-1. Dans cette thèse, nous avons effectué la détermination expérimentale de la solubilité du dioxyde de carbone dans des saumures typiques du bassin Rhénan aux conditions de pression et de température du projet FONGEOSEC. Des réflexions sont proposées quant à une méthodologie d’analyse de solubilité du méthane et de l’azote dans des phases aqueuses. Nous observons aussi que dans les conditions de pression et température de fond du puits, la solubilité du dioxyde de carbone dans les saumures typiques du bassin Rhénan est la plus élevée parmi toutes les conditions caractérisées. Une étude du sating-out effect dans ces saumures est également proposée dans cette thèse. Enfin, il est remarqué que le modèle de Pitzer (Pitzer.dat sur PhreeqC) semble prédire de façon correcte nos mesures expérimentales à 333.15 K, mais il perd son efficacité à 453.15 K. Dans cette condition, le modèle E-NRTL (Simulis®) semblerait être plus approprié. / This thesis was part of the FONGEOSEC project, which aims to develop the deep geothermal energy sector in France through the the design of a geothermal power pilot plant on the Upper Rhine Graben. This project is controled by Fonroche Géothermie, which manages a consortitium of more than ten academic and industrial partners. The French Environment & Energy Management Agency (ADEME) participates at the fundings of the project.Therefore, the work exposed in this document concerns the thermodynamic characterisation of geothermal fluids (hot brines containing dissolved gases) from the target region of this project. It is thus necessary to determine the solubility of each gas dissolved in these brines at the pressure, temperature and salinity conditions of geothermal energy exploitation.Thermodynamic models that predict liquid-vapour phase equilibrium can be used to estimate these solubilities. However, if there is a lack of experimental measures on the pressure, temperature and salinity conditions of interest, it will not be possible to regress these models interaction parameters and, therefore, it will be difficult to have precise solubility results from these thermodynamic simulations. Thus, this thesis has focused on the experimental study of gas solubilities in brines representing the Upper Rhine Graben fluids. The pressure range of the FONGEOSEC project goes from 6.0 MPa to 40.0 MPa for temperatures of 333.15 K and 453.15 K. The experimental setup used on this thesis can operate at these conditions. Dissolved gases in the brines concerned by this project are mainly composed by carbon dioxyde (CO2), and then by nitrgen (N2) and methane (CH4) at lower quantitites. Dissolved salts in these fluids are basically chloride sodium and chloride calcium, at molalities of 1.2 mol NaCl-0.2 mol CaCl2.Kg H2O-1.On the scope of this thesis, we have performed the experimental determination of carbon dioxyde solubility in Upper Rhine Graben-type brines at the pressure and temperature conditions of the FONGEOSEC project. We propose a discussion about an analysis methodology for measuring nitrogen and methane solubility in aqueous phases. We also observed that at the pressure and temperature conditions found at the bottom of the production well, carbon dioxyde solubility in the Upper Rhine Graben-type brines reaches its highest value among all the conditions studied in this thesis. A salting-out effect study in these brines is also proposed in this document. Finally, it is noticed that the Pitzer model (Pitzer.dat at PhreeqC) seems to predict properly our experimental data at 333.15 K, but it is less efficient at 453.15 K. In this condition, the E-NRTL model (Simulis®) seems to be more appropriate.

Évaluer le potentiel géothermique d’une mine inondée est une tâche complexe qui nécessite la réalisation d’un test hydraulique et la modélisation de l’écoulement de l’eau souterraine à travers un réservoir non conventionnel. Cette tâche peut être efficacement complétée en réutilisant les anciennes infrastructures minières. Lors de cette étude, un puits de ventilation débouchant dans les galeries souterraines a été converti en ouvrage de captage d’eau profond. L’ancien puits de ventilation 1100 des Mines Gaspé à Murdochville au Canada a été utilisé pour effectuer un essai de pompage durant 3 semaines à un débit moyen de 0.062 m3/s. L’objectif principal des travaux encourus était d’évaluer le potentiel géothermique du site minier. Lors de l’essai, moins de 3,65 m de rabattement ont été observés et la température moyenne de l’eau pompée a été de 6,7 °C. Le comportement de la nappe souterraine durant l’essai a été reproduit à l’aide d’un modèle d’éléments finis tridimensionnel simulant l’écoulement de l’eau à travers les galeries. Les prédictions du modèle ainsi qu’un bilan énergétique simplifié suggèrent que le taux d’extraction d’énergie durable est atteint à un débit de pompage de 0.049 m3/s, ce qui indique un potentiel géothermique de 765 kW. Cette énergie pourrait être extraite à l’aide de pompes à chaleur géothermiques afin de chauffer les bâtiments du parc industriel de Murdochville. / Assessing the low-temperature geothermal potential of a flooded mine site is a complex task involving hydraulic testing and modelling of an unusual man-made reservoir. It can be achieved efficiently taking advantages of the former mine infrastructures such as mining shaft that provided here a deep well directly connected to the mine workings. The former mining shaft P1100 of the flooded Gaspé Mines near Murdochville, Canada, was used to perform a pumping test during a study with the objective of assessing the geothermal potential of the mine site. Water was pumped during 3 weeks at a rate averaging 0.062 m3/s with a mean recovery temperature equal to 6.7 °C and less than 3.65 m drawdown was observed. The hydraulic response of the pumping test was reproduced with a threedimensional finite element model that simulates groundwater flow through the mine workings. Model predictions and a simplified energy balance calculation suggested that the sustainable energy extraction rate is attained at a pumping rate of 0.049 m3/s which yield a geothermal potential of 765 kW. This energy could be extracted with geothermal heat pumps used for space heating at Murdochville industrial park.

Dans le cadre de la présente étude, un modèle géothermique tridimensionnel du pergélisol a été développé pour la toute première fois à l’échelle locale d’un village nordique afin d’évaluer les impacts des changements climatiques sur le pergélisol. Ce village est celui de Salluit au Québec nordique dans la zone de pergélisol continu. La majeure partie des infrastructures de ce village sont fondées sur un pergélisol salin riche en glace vulnérable aux changements climatiques. Le code numérique d’HydroGeoSphere a été modifié pour simuler la transmission de la chaleur par conduction avec changement de phase et tenir compte de la variation des propriétés thermiques du sol avec la température. Différents scénarios de changements climatiques lors du prochain siècle ont été utilisés pour alimenter le modèle géothermique. La variabilité spatiale du régime thermique du pergélisol simulée à l’échelle d’un village nordique est majeure lorsque l’hétérogénéité des dépôts, leur contenu en glace et les différentes conditions de surface observées sur le terrain sont considérés dans une modélisation numérique tridimensionnelle. L’influence des flux de chaleur latéraux dans cette modélisation produit des résultats qui sont dans certains cas significativement différents de ceux simulés au moyen d’un modèle unidimensionnel, particulièrement pour les sols sous d’importantes congères. Selon les scénarios de changements climatiques jugés pessimistes pour le maintien du pergélisol, un talik permanent se formera à partir de la seconde moitié du 21ième siècle sous les congères, en bordure d’un remblai routier où la neige et l’eau de fonte s’accumule à la surface d’un sol mal drainé et dans l’argile saline du secteur nord-est du village de Salluit. L’épaisseur maximale du talik à la fin du 21ième siècle varie entre 1,3 et 19,3 m respectivement dans l’argile et le roc sous les congères, alors qu’elle est de 2,8 m dans l’argile saline du secteur nord-est. Seuls les tills et le roc recouverts par d’importantes congères seront affectés par la présence d’un talik à la toute fin du 21ième siècle advenant un réchauffement climatique optimiste. Une dégradation plus importante du pergélisol sous les remblais routiers que pour un même sol en terrain naturel sous un couvert nival mince est à prévoir selon les scénarios climatiques pessimistes. Les infrastructures contruites sur des remblais seront potentiellement affectées par des tassements différentiels.

Pour la climatisation de bâtiments ou d'installations industrielles, les nappes d'eaux superficielles représentent une source de frigories très convoitée. Leurs exploitations intensives depuis plusieurs dizaines d'années conjuguées au redéploiement de la filière géothermique ces dernières années, soulèvent des préoccupations quant à la préservation des ressources en eau. Dans ce contexte, la présente étude vise à évaluer l'impact de variations locales de température sur la qualité physico-chimique et microbiologique des eaux souterraines sur la base (i) de suivis in-situ au niveau de 3 installations réelles exploitant les nappes d'eaux superficielles et, (ii) d'expérimentations sur un pilote (BIOTHERMEX) permettant de reproduire, en conditions parfaitement maitrisées, l'effet de la propagation d'un panache thermique dans un modèle réduit d'aquifère. Dans la gamme de température relevée sur site, les principaux résultats obtenus montrent que les impacts thermiques sont circonscrits au voisinage immédiat de l'installation, pouvant altérer jusqu'à plus d'une dizaine de degrés la sténothermie des nappes. Le suivi des paramètres physico-chimiques n'ont pas fait apparaitre de perturbations significatives sur la période de surveillance, constat étayé par des modélisations hydrogéochimiques. En revanche, une influence significative a été relevée au niveau des principaux descripteurs microbiologiques (activité, diversité de la microflore totale). Enfin, les expériences menées à l'échelle du laboratoire ont permis d'appréhender finement le comportement réactionnel du système et de définir une température de réinjection critique, au-delà de laquelle des désordres potentiels sont attendus.

La compréhension des écoulements souterrains est importante pour de nombreuses applications comme l’énergie ou le stockage des déchets nucléaires. Cette thèse, effectuée en collaboration avec le Bureau de Recherches Géologiques et Minières (BRGM), est dédiée à la simulation des écoulements diphasiques compositionnels thermiques en milieux poreux et ses applications à la géothermie haute énergie et plus particulièrement au champ géothermique de Bouillante (Guadeloupe). Tout d’abord, deux formulations à variables persistantes sont comparées en termes d’implémentation et de convergence numérique. Dans ces deux formulations, les fractions molaires d’une phase absente sont étendues par celles à l’équilibre thermodynamique avec la phase présente. Il en résulte que l’ensemble des variables principales et des équations ne dépend pas de l’ensemble de phases présentes. De plus, l’équilibre thermodynamique est exprimé par une contrainte de complémentarité pour chacune des phases, ce qui permet l’utilisation de méthodes de type semi-smooth Newton pour résoudre les systèmes non-linéaires. D’autre part, cette thèse présente une nouvelle méthodologie combinant des discrétisations centrées aux noeuds (le schéma Vertex Approximate Gradient - VAG) et aux faces (le schéma Hybrid Finite Volume - HFV) sur une partition arbitraire des ensembles de mailles ou de faces, dans le but d’adapter le choix du schéma aux différentes parties du maillage. En effet, les maillages hybrides composés de différents types de mailles sont plus adaptés à la discrétisation de la géologie et de la géométrie des différents domaines d’un système géothermique. Ainsi le schéma peut être choisi localement en fonction de la géométrie de la maille et des propriétés pétrophysiques. L’analyse de convergence est effectuée dans le cadre des discrétisations Gradient pour des problèmes de diffusion du second ordre et la convergence est confirmée numériquement sur différents types de maillages hybrides 3D. Ensuite la discrétisation VAG-HFV est étendue au cas des écoulements de Darcy diphasiques non-isothermes compositionnels et est appliquée au cas test 2D représentant le plan de faille vertical du réservoir géothermique de Bouillante. Un autre aspect important de la modélisation des flux géothermiques consiste à prendre en compte les interactions entre le flux dans le milieu poreux et l’atmosphère. Puisque le couplage entre le modèle poreux et un modèle 2D surfacique ou 3D atmosphérique n’est pas réaliste en terme de coût de calcul aux échelles spatiale et temporelle géologiques, l’interaction sol-atmosphère est modélisée grâce à une condition limite prenant en compte l’équilibre de matière et d’énergie à l’interface. Ce modèle considère une couche limite atmosphérique avec transfert convectif molaire et thermique (en supposant l’évaporation de la phase liquide), une condition de débordement liquide aux surfaces d’infiltration, ainsi que le rayonnement thermique et la recharge en eau douce due aux précipitations. Cette condition limite est évaluée à l’aide d’une solution de référence couplant les écoulements non-isothermes liquide-gaz en milieu poreux et le gaz dans le milieu libre. Elle est ensuite étudiée numériquement en terme de convergence et de solution sur des cas tests géothermiques, dont le plan de faille vertical du réservoir géothermique de Bouillante. En complément est présenté le travail issu d’une collaboration lors de l’école d’été du CEMRACS 2016. Le projet consistait à ajouter un modèle de puits multi-branche thermique au code ComPASS, un nouveau simulateur géothermique parallèle basé sur des maillages non-structurés avec la possibilité de représenter des fractures. / The study of the subsurface flows is important for various applications such as energy or nuclear waste storage. This thesis, performed in collaboration with the French Geological Survey (BRGM), is dedicated to the simulation of non-isothermal compositional two-phase flows in porous media and its applications to high-energy geothermal fields and more precisely to the Bouillante field (Guadeloupe, French West Indies). First of all, two persistent variable formulations are compared in terms of implementation and numerical convergence. In these two formulations, the choice of the principal variables is based on with the extension of the phase molar fractions by the one at thermodynamic equilibrium with the present phase. It results that the set of principal variables and equations does not depend on the set of present phases. It also has the advantage to express the thermodynamic equilibrium as complementarity constraints, which allows the use of semi-smooth Newton methods to solve the non-linear systems. Moreover, this thesis presents a new methodology to combine a node-centered discretization (the Vertex Approximate Gradient scheme - VAG) and a face-centered discretization (the Hybrid Finite Volume scheme - HFV) on arbitrary subsets of cells or faces in order to choose the best-suited scheme in different parts of the mesh. Indeed, hybrid meshes composed of different types of cells are best suited to discretize the geology and geometry of the different parts of the geothermal system. Then, the scheme is adapted locally to the type of mesh/ cells and to petrophysical properties. The convergence analysis is performed in the gradient discretization framework over second order diffusion problems and the convergence is checked numerically on various types of hybrid three-dimensional meshes. Then, the VAG-HFV discretization is extended to non-isothermal compositional liquid-gas Darcy flows and is applied on the two dimensional cross-section of the Bouillante high temperature geothermal reservoir. Another important aspect of the geothermal flows modelling consists in considering the interactions between the porous medium and the atmosphere. Since the coupling between the porous medium and the 2D surface of 3D atmospheric flows is not computationally realistic at the space and time scales of a geothermal flow, the soil-atmosphere interaction is modelled using an advanced boundary condition accounting for the matter (mole) and energy balance at the interface. The model considers an atmospheric boundary layer with convective molar and energy transfers (assuming the vaporization of the liquid phase in the atmosphere), a liquid outflow condition at seepage surfaces, as well as the heat radiation and the precipitation influx. This boundary condition is assessed using a reference solution coupling the Darcy flow to a full-dimensional gas free flow. Then, it is studied numerically in terms of solution and convergence of the Newton-min non-linear solvers on several geothermal test cases including two-dimensional simulations of the Bouillante geothermal field. In addition is presented the collaborative project which took place during the CEMRACS summer school 2016. The project consisted in adding a multibranch thermal well model into the ComPASS code, a new geothermal simulator based on unstructured meshes and adapted to parallel distributed architectures with the ability to represent fractures.

Les doublets géothermiques sur aquifères superficiels jouent un rôle important en France pour le chauffage, le rafraichissement et la production d’eau chaude sanitaire des bâtiments. La pérennité de ces installations est cependant conditionnée par la possibilité de pouvoir garantir dans le temps la production puis la réinjection de l’eau souterraine. Le colmatage de la boucle géothermale est un problème majeur qui affecte de nombreuses installations sur l’ensemble du territoire français et remet en cause leur viabilité technico-économique. La compréhension et la prédiction de ce phénomène nécessite de nouveaux moyens de caractérisation et de modélisation des processus biogéochimiques couplés au fonctionnement des boucles géothermales dans des environnements de subsurface hétérogènes. Cette thèse a ainsi fait l’objet d’une nouvelle collaboration entre le laboratoire Géosciences Rennes et l’entreprise Antea group pour identifier les mécanismes à l’origine de ces phénomènes de colmatage et d’en caractériser les paramètres influents dans le but de construire un outil d’analyse des risques permettant l’anticipation de ces processus de colmatage. Le premier volet de la thèse expose la synthèse des typologies de colmatage (biogéochimique, chimique, biologique et physique). Cette synthèse est issue de l’étude de la littérature sur la problématique de colmatage des forages d’eau et des retours d’expérience concernant les doublets géothermiques recensés et étudiés durant cette thèse. Le second volet présente les inventaires régionaux des problématiques de colmatage recensées en France dans différents contextes hydrogéologiques. Ces inventaires ont permis de délivrer un état des lieux à grande échelle de la problématique de colmatage et d’étudier les contextes associés à chaque type de colmatage. Dans le troisième volet, de nouvelles méthodologies de caractérisation in-situ des phénomènes de colmatage biogéochimique liés à l’oxydation du manganèse et du fer sont présentées. Elles ont été développées dans le cadre de ces travaux de thèse à partir des investigations menées sur différents sites confrontés à des problèmes d’exploitation. Ces méthodologies d’investigations pluridisciplinaires couplent la mesure des propriétés hydrauliques, des concentrations en éléments chimiques et de la diversité microbienne afin d’identifier la problématique impactant le fonctionnement de l’installation et d’en définir in fine les causes. Le quatrième volet expose les résultats d’une campagne de terrain effectuée sur un doublet géothermique impacté par un processus de colmatage biogéochimique. Cette campagne a permis d’identifier les composantes clés de la réactivité biochimique impliquée dans le colmatage : distribution du flux, hétérogénéité chimique et diversité microbiologique. A l’aide du code de calcul PHREEQC, un modèle géochimique simulant les cinétiques de précipitation observées a été développé. Il permet l’étude quantitative des mécanismes biogéochimiques favorisant l’apparition rapide du colmatage. Enfin, les retours d’expériences de l’exploitation des doublets sur nappes superficielles ont démontré le besoin d’une méthodologie fiable d’analyse des risques, permettant d’anticiper l’apparition des processus de colmatage à chaque étape de la vie du projet. Ainsi, à partir de la synthèse de la littérature scientifique et technique et des conclusions des études menées durant cette thèse, les facteurs de risques d’apparition des phénomènes de colmatage ont été déterminés. L’analyse des risques de colmatage intégrant ces différents facteurs a été implémentée à travers l’élaboration de méthodes développées sous Python 3. L’outil ARCADE (Analyse des Risques de Colmatage et Aide à la Décision) a été conçu d’une part pour évaluer le risque et d’autre part pour informer l’utilisateur averti des bonnes pratiques et moyens préventifs. Ces bonnes pratiques d’analyse et de gestion préventive sont présentées dans le dernier volet de cette thèse. / The sustainability of geothermal systems using shallow aquifers for heating, cooling and hot water production depends on the possibility to ensure, over long time-scales, the production and the reinjection of groundwater in the aquifer. Clogging of the geothermal loop is a major issue affecting the technical and economic viabilities of numerous operations in France. The understanding and prediction of this phenomenon requires new methods of characterization and modelling of biogeochemical processes coupled to the operation of geothermal loops in heterogeneous subsurface environments. This thesis is thus the result of a new collaboration between the Géosciences Rennes lab and Antea group to identify the mechanisms at the origin of clogging phenomena and characterize their controlling parameters, in order to establish a risk assessment tool allowing the anticipation of clogging processes. The first part of the thesis describes the main clogging processes (biogeochemical, chemical, biological and physical). This synthesis is the result of (1) the study of the literature dealing with the clogging of water wells and (2) our feedbacks on the geothermal doublets identified and studied during this thesis. The second part presents the regional inventories of shallow geothermal systems impacted by clogging problems identified in different hydrogeological contexts in France. These inventories provide a large-scale perspective of clogging phenomena and allow to study the contexts associated with each type of clogging processes. In the third part, new methodologies for the in-situ characterization of biogeochemical clogging phenomena linked to manganese and iron oxidation are presented through case studies of sites affected by clogging issues. These interdisciplinary studies couple the measurement of hydraulic properties, chemical element concentrations and bacterial diversity, to identify the specific issue impacting the operation and to define its causes. The fourth part presents the results of an interdisciplinary field campaign carried out on a geothermal doublet impacted by a biogeochemical clogging process. This campaign documented the key components involved in mixing induced biogeochemical reactivity: flow distribution, chemical heterogeneity and microbiological diversity. Using PHREEQC, a geochemical model simulating observed kinetics of precipitation was developed in order to quantitatively explore the biogeochemical mechanisms favoring rapid clogging. Feedback from shallow geothermal systems operation has demonstrated the need for a reliable risk analysis methodology that allowed to anticipate the apparition of clogging processes at each stage of the project life (part 5). From the synthesis of the scientific & technical literature and the conclusions of the studies carried out during this thesis, the risk factors for the appearance of clogging phenomena were determined. A clogging risk analysis integrating these factors was implemented through the development of methods developed under Python 3. The methodology of the ARCADE tool (Analyse des Risques de Colmatage et Aide à la Décision) is designed to assess the risk and to inform users of good practices and preventive methods. These good practices for analysis and preventive methods are presented in the last part of this thesis.