La géothermie est considérée comme une source d'énergie propre et inépuisable à échelle humaine. Actuellement, le rendement des centrales géothermiques est limité à l'exploitation de fluides de températures inférieures à 300°C. L'association de l'activité tectonique et volcanique aux dorsales océaniques fait de l'Islande un lieu où l'extraction de fluides supercritiques (T°>400°C) peut être envisagée. Cette exploitation pourrait multiplier par dix la puissance électrique délivrée par les puits géothermiques. Ces fluides peuvent-ils circuler dans la croûte océanique ? Ce travail propose de contraindre les observations géophysiques et de prédire le fonctionnement des réservoirs hydro-géothermaux de très haute température par l'étude des propriétés physiques des basaltes. La première approche est focalisée sur l'étude de roches ayant accueilli une circulation hydrothermale par le passé. L'étude de ces roches au site ODP 1256, montre que leur porosité est associée à la présence de minéraux d'altération hydrothermale du faciès amphibolite (T°>500°C ). La seconde approche a consisté à recréer en laboratoire les conditions des systèmes hydrothermaux à très haute température afin de prédire les propriétés mécaniques et électriques des basaltes dans ces conditions. Les résultats mécaniques indiquent que la transition fragile/ductile, souvent associée à une forte décroissance de perméabilité, intervient à une température d'environ 550°C. La mise en place d'une cellule de mesure de la conductivité électrique de haute températures a fournit les premiers résultats utiles à l'analyse des données géophysiques. Appliqués aux conditions de la croûte basaltique Islandaise, ces résultats indiquent que des fluides hydrothermaux pourraient circuler au moins transitoirement à l'état supercritique jusqu'à ~ 5 km de profondeur. / Geothermal energy is considered as a green and infinite energy source at human scale. Currently, the yield of geothermal power plants is limited to temperatures of the operating fluid which 300 °C. From tectonic and volcanic activity at mid-ocean ridges, Iceland is a locuswhere supercritical fluid extraction (T > 400° C) can considered for the near future. Exploiting such fluids could theoretically multiply by a factor of ten the electrical power delivered by geothermal wells. Can such fluids circulate at the base of brittle oceanic crust? This work investigates the petrophysical properties of basalts in order to constrain geophysical observations in Iceland and predict the behavior of very high temperature hydro-geothermal reservoirs. The first approach consisted in studying the physical properties of rocks that have hosted deep hydrothermal circulations at oceanic ridges. The study of these rocks at ODP Site 1256 shows that the porosity mea sured both in the field and in the lab is associated with amphibolite facies alteration minerals (T > 500° C). The second approach was to recreate in the laboratory the conditions of pressure, temperature and pore fluid pressure of high temperature to supercritical hydrothermal systems to predict the mechanical and electrical properties of basalts under these conditions. The mechanical results indicate that the brittle/ductile transition occurs at a temperature of about 550° C, where a strong permeability decrease is expected. The implementation and calibration of a new cell for measuring electrical conductivity at high temperature provide the first results for the interpretation of geophysical data. When applied to basaltic crustal conditions in Iceland, these results indicate that hydrothermal fluids could circulate, at least temporarily, in a supercritical state up to 5 km depth.

Le développement et l'exploitation de réservoirs géothermiques profonds s'accompagnent généralement d'une sismicité induite - un effet secondaire indésirable. Cette recherche est axées sur l'utilisation de simulations numériques pour étudier la propagation des fractures hydrauliques et la réactivation de failles préexistantes lors de la stimulation hydraulique des réservoirs afin de mieux comprendre le comportement du réservoir fracturé et de réduire le risque potentiel de sismicité induite.La sismicité induite est d'abord étudiée du point de vue de l'utilisation de la loi de conservation de l'énergie afin d'expliquer le mécanisme de génération d'ondes élastiques à partir d'une rupture de roche. Ensuite, une approche approximative est proposée pour calculer les accélérations de pointe (PGA) induites par le glissement de faille. Les PGA calculés à la surface du sol servent à évaluer la perception humaine des ondes sismiques et le potentiel de dégradation des structures. / The development and the exploitation of deep geothermal reservoirs are usually accompanied with induced seismicity - an unwanted side effect. This research is focused on using numerical simulations to investigate the propagation of hydraulic fractures and the reactivation of pre-existing faults during the hydraulic stimulation of the reservoirs in an effort to better understand the fractured reservoir behavior and to reduce the potential risk of induced seismicity.The induced seismicity is studied first from the standpoint of using the law of energy conservation in order to explain the mechanism of generating elastic waves from rock failure. Then an approximate approach is proposed to calculate the peak ground accelerations (PGAs) that are induced by the fault slip. The computed PGAs on ground surface are used to assess the human perception of the seismic waves and the damage potential to structures.

Cette thèse de doctorat a pour objet l'étude d'échangeurs géothermiques hélicoïdaux implantés dans le proche sous-sol. Ces objets font partie de la famille des échangeurs géothermiques compacts, et sont une alternative aux sondes verticales et aux nappes horizontales plus couramment utilisées pour chauffer et rafraîchir des bâtiments (géothermie sèche, très basse enthalpie). L'étude à la fois théorique et expérimentale de ces échangeurs mène à différents modèles, numériques comme analytiques, permettant de prévoir le comportement en température du sous-sol comme du fluide caloporteur. Des résultats de simulation intégrant le couplage aux bâtiments sont présentés, et des outils logiciels ont été créés pour aider au dimensionnement d'installations. De surcroît, un "test de réponse thermique" adapté à ces échangeurs a été développé, afin de permettre d'estimer les propriétés thermiques moyennes du sous-sol et la résistance thermique liée à l'installation de l'échangeur. Enfin, des conseils pratiques sont donnés, dont le but est d'obtenir les températures d'utilisation les mieux adaptées aux pompes à chaleur, et en conséquence d'améliorer les coefficients de performance du système. / This thesis consists in a study of geothermal helical heat exchangers buried in the close underground. These objects belong to the compact exchangers kind and may be used instead of borehole heat exchanger or geothermal horizontal exchangers to cope with the needs in heating and cooling buildings with a heat pump. A theoretical and experimental study of these exchangers leads to numerical and analytical models which can be used to estimate the temperatures in the underground and in the exchanger over years. Results of simulations including a building are given, and softwares were developed to help with the sizing of installations. Moreover, a specific "thermal response test" for helical heat exchangers was developed, in order to estimate the mean thermal properties of the underground and the thermal resistance between the heat transfer fluid and the underground. At least, practical advices are given, which goal is to get more valuable temperatures at the heat pump and as a consequence better coefficients of performance of the geothermal system.

Dans de nombreux réservoirs géothermiques profonds, l'échangeur thermique est un massif de roches chaudes fracturées. Comment la morphologie des fractures influence-t-elle le flux hydraulique et le champ de température lorsque de l'eau froide est injectée dans ce milieu ? Cette question est abordée à l'échelle de la fracture, de manière numérique, à l'aide de deux méthodes : par différences finies, en 2D, et par méthode de Boltzmann sur réseau, en 3D. Les calculs en différences finies montrent qu'une ouverture auto-affine induit une variabilité spatiale du flux hydraulique et du champ de température. Dans certains cas, il se crée une chenalisation du fluide, qui réduit l'efficacité du transfert thermique entre la roche et le fluide. Ces caractéristiques hydro-thermiques sont essentiellement engendrées par les plus grandes longueurs d'ondes de ces ouvertures auto-affines. Une première application à grande échelle, dimensionnée suivant les caractéristiques du site géothermique de Soultz-sous-Forêts (Alsace, France), est proposée. Une deuxième application consiste à évaluer la perméabilité du sous-sol de Draix (Alpes, France). Pour cela, la géométrie de fractures présentes dans le sous-sol est caractérisée précisément à partir de carottes géologiques. La topographie des fractures, mesurée avec un laser, est utilisée pour reconstituer la morphologie de l'ouverture des fractures. Ces données sont utilisées pour quantifier la perméabilité, ici élevée -- de l'ordre de 10^-9-10^-8 m^2. À l'aide de méthodes de Boltzmann sur réseau, la modélisation du comportement hydro-thermique dans des fractures présentant une topographie avec de fortes pentes est abordée. Au voisinage d'une aspérité de forme simple, le comportement hydro-thermique modélisé est très différent (recirculation) de celui dans le reste de la fracture. L'ensemble de ces modélisations suggère qu'il est nécessaire de prendre en compte la morphologie des fractures pour estimer le comportement hydro-thermique en géothermie.

Cette thèse s’inscrit dans le cadre du projet ANR « Uncertain flow optimization » (UFO) consacré au développement et à l’application de méthodes efficaces de quantification d’incertitudes pour l’analyse et l’optimisation d’écoulements. Dans ce contexte, ces méthodes sont appliquées à des gazosiphons de grande échelle utilisés comme pompe. Plus particulièrement, on s’intéresse à de tels gazosiphons choisis pour constituer l’organe central d’un système géothermique innovant,  compatible avec un environnement urbain. On souhaite en quantifier le potentiel énergétique par voie numérique avec la recherche d’un compromis entre justesse des résultats et efficacité optimale. La simulation de l’écoulement diphasique produit dans le gazosiphon est fondée sur un modèle quasi-1D à flux de dérive et s’appuie sur une démarche de résolution implicite. Les résultats sont validés sur les études expérimentales les plus pertinents de la littérature, dont aucune toutefois n’atteint les longueurs requises de l'ordre du kilomètre. Le code de simulation du gazosiphon fait ensuite l’objet d’une démarche de prise en compte d’incertitudes physiques et de modélisation, précédée par une analyse de deux méthodes de quantification d’incertitude : une méthode non-intrusive de type chaos polynomial, et une méthode plus récente dite semi-intrusive qui fut développée en amont du projet UFO. Cet outil est intégré dans une modélisation simplifiée du système géothermique urbain dans son ensemble impliquant les composants en surface, notamment le compresseur d'air. Il en résulte une optimisation énergétique robuste préliminaire de deux variantes du système géothermique urbain proposé, respectivement de récupération de chaleur et de production d’électricité. / This PhD thesis is part of the ANR project « Uncertain Flow Optimization » (UFO). The project is devoted to the development and application of efficient uncertainty quantification methods for flow analysis and optimization. In this framework, these methods are applied to the study of a large-scale airlift pump. The airlift pump is selected to be part of an innovative geothermal system, which can be exploited within an urban environment. We wish to quantify and optimize the energy potential of this new system with numerical tools. They provide both good accuracy and efficiency properties. The airlift two-phase flow simulation is based on a quasi one-dimensional drift flux model, which is implicitly solved. The solver is validated by comparison with relevant experimental airlift studies from the literature. However, these studies remain below the kilometric-targeted pipe length. Thanks to the analysis of two uncertainty quantification methods, a non-intrusive approach relying on polynomial chaos expansion and a new semi-intrusive method developed ahead of the UFO project, we perform airlift pump simulations taking into account physical and modelling uncertainties. This numerical tool is inserted into a simplified model of the complete urban geothermal system that involves surface devices, such as an air compressor. Finally, a robust preliminary optimization process is performed for two versions of the proposed geothermal urban system. They are designed respectively for heat recovery and electricity production.

En France métropolitaine, les projets de géothermie haute température pour la production d’électricité sont principalement localisés dans le socle des fossés d’effondrement liés à la mise en place du Rift Ouest Européen. Le socle de ces fossés a été étudié sur deux analogues à l’affleurement sur les épaules du rift : les Vosges du Nord pour le fossé Rhénan et la bordure Est du Massif central pour la fosse de Valence. Cette étude a permis de montrer que le réseau de failles s’organise selon trois ordres de grandeurs de longueurs et d’espacements caractéristiques qui individualisent des blocs structuraux. Les orientations et l’espacement des failles formant ces blocs et la présence ou l’absence de certains ordres de grandeurs sont le résultat de l’héritage anté-rift du socle, ainsi que du mécanisme d’ouverture du bassin. Le potentiel réservoir des formations de socle et de la couverture surincombante a été analysé au regard de ces zones de failles et de l’altération supergène qui affecte le toit du socle. Ainsi, les lithologies potentiellement rencontrées en base des fossés ont pu être classées en fonction du potentiel de développement de porosité et de perméabilité matricielle dans les cœurs de failles, les zones endommagées et le réseau pervasif de fractures dans le protolithe. L’évolution de la fracturation dans les zones de failles a également pu être appréhendé, et une méthodologie de modélisation double milieu a été élaborée pour caractériser la porosité et la perméabilité de fractures et modéliser le fonctionnement d’un doublet géothermique dans une faille synthéthique équivalente aux cibles des projets géothermiques / High temperature geothermal projects for electricity production are in France mostly localized in the basement of basins linked to the West European Rifting event. The basement of theses basins have been studied on two outcrop analogues at the shoulders of the rift: the Northern Vosges mountains for the Upper Rhine Graben and the Eastern border of the Massif central for the Valence Graben. This study has shown that the fault network is organized in three orders of size with characteristic length and spacing, and that form characteristic structural blocks. The orientation and spacing of these faults and the presence or absence of some size orders are the result of structural inheritance of the basement and the mechanism of the basin opening. The reservoir potential of the basement rocks and the surrounding sedimentary cover has been analysed in light of the fault zones structure and the weathered layer at the top of the basement. Thus, the basement rocks of the basins has been classified in light of the potential of matrix porosity and permeability development in the fault core, the damaged zone and the fractured protolith. The evolution of the fracture network in the fault zone has been studied too, with the development of a double-porosity model in order to characterize the fracture porosity and permeability, and therefore simulate the working of a geothermal doublet in a synthethic fault zone analogue of the geothermal drilling targets

Le Buntsandstein, localisé dans le graben du Rhin supérieur, apparait comme une cible intéressante pour la géothermie, associant une formation argilo-gréseuse à l’anomalie thermique régionale. Cette étude vise à caractériser les propriétés pétrophysiques de ces grès ainsi que la fracturation les affectant, dans le but de fournir un modèle conceptuel de la formation qui servira de guide pour son exploitation futur. Les faciès sédimentaires sont composés par cinq faciès pétrographiques (grès propres, grès à enrobage argileux, grès à matrice argileuse, grès silicifiés et grès à ciment carbonaté), qui se répartissent dans des proportions variables et contrôlent une partie des propriétés pétrophysiques mesurées à l’échelle matricielle. La comparaison des données pétrophysiques, des données macroscopiques issues d’une analyse de gradients de température, des données de modélisation et de la fracturation permet de construire un modèle de circulation dans le réservoir. Ces analyses mettent en avant le rôle de la zone endommagée des zones de faille pour le transfert de fluides à grande échelle, mais aussi celui de deux faciès sédimentaires : les grès déposés dans un environnement de Playa Lake et fluvio-éolien. L’analyse de différents affleurements montre que la fracturation évolue en fonction de la situation dans la pile sédimentaire et en fonction de la situation par rapport aux accidents tectoniques majeures. / The Buntsandstein, located in the Upper Rhine Graben, appears to be an easy target forgeothermal exploitation, linking sandstone and clay with the regional thermal anomaly. This study aims at characterizing petrophysical characteristics of these sandstones as well as the fracturation affecting them, with the intention of providing a conceptual model of the formation which will act as guide for future exploitation. The sedimentary facies are composed by five petrographical facies (clean sandstones, sandstones with clayed coating, clay matrix sandstones, silicified sandstones and carbonated matrix sandstones) which split with variable proportions and control a part of petrophysical properties measured at matrix scale. The comparison between petrophysical data, macroscopic data from temperature gradient analysis, modelling data and fracturing, allows the building of a Buntsandstein Sandstones fluids circulation conceptual model. This analysis points the role of the damage zone of fault zones for fluids transfer at large scale, but also that of two sedimentary facies: marginal erg and Playa Lake. The analysis of different outcrops shows that the fracturation evolves according to the situation in the sedimentary pile and according to the situation in comparison with major tectonic accidents.

Les milieux cristallins fracturés constituent un potentiel géothermique non négligeable. Il est essentiel d'améliorer son exploitation, pour la géothermie basse et haute énergie, afin de répondre au processus de transition énergétique. Pour cela, la compréhension des mécanismes de transport thermique dans les milieux fracturés est fondamentale. Le transport de chaleur est fortement influencé par l'hétérogénéité hydrodynamique des milieux fracturés et par la géométrie des fractures et des blocs matriciels. A travers des travaux basés sur des développements analytiques et numériques ainsi que des expériences sur site, l'objectif de cette thèse est ainsi de mieux évaluer l'impact de la géométrie des fractures, que ce soit à l'échelle d'un réseau de fractures, ou à l’échelle d’une fracture, sur le transport et le stockage d’énergie thermique dans les milieux cristallins fracturés. Des simulations numériques du transport de chaleur dans un réseau simple de fractures planes et bien connectées ont permis de caractériser l'impact de la géométrie du système de fractures sur le stockage thermique. Deux régimes sont mis en évidence. Tout d'abord, à court terme, la densité de fractures ou de chemins préférentiels, caractérisant la surface d'échange, contrôle l'échange thermique. Puis, à long terme, c'est le volume de roche total entre les fractures qui contrôle le stockage thermique. Ce modèle ne prend toutefois pas en compte la variabilité des ouvertures à l'échelle de la fracture qui est particulièrement présente dans les réseaux de fractures naturels. Des tests de traçage thermique et de soluté ont ainsi été réalisés pour caractériser le transport de chaleur dans un milieu fracturé sur le site de Ploemeur (SNO H+). Pour interpréter les traçages, les expressions analytiques du retard et de l'amplitude du pic de la courbe de restitution thermique ont été développées pour différentes géométries de fractures : fractures planes et chenaux. Ces expressions constituent un outil puissant et novateur pour caractériser la géométrie des fractures lors de tests de traçage thermique mais aussi pour prédire le déplacement du front thermique et la durée de vie des systèmes géothermiques à partir de tests de traçage de soluté. La comparaison de ces expressions avec les résultats expérimentaux permet de mettre en évidence l'importante chenalisation des flux, induisant l'arrivée anticipée du traceur thermique. / Fractured crystalline media has a significant geothermal potential. Its exploitation, for low and high enthalpy geothermal power generation, could be enhanced to satisfy the energy transition process. For this, understanding thermal transport processes in fractured media is fundamental. Heat transport is strongly influenced by hydrodynamics heterogeneity of fractured media and by fracture and matrix block geometry. Through analytical and numerical modelling and field site experiments, the aim of this thesis is thus to better assess the impact of fracture geometry on thermal transport and storage in fractured crystalline rock, at fracture and fracture network scale. Numerical simulations of heat transport in a simple network of well connected plane fractures allowed us to characterize the impact of the fracture system geometry on thermal storage. Two regimes are highlighted. First, at short term, the density of fractures, or preferential paths, controls heat exchanges. Then at long term, the total rock volume between the fractures controls thermal storage. This model does not take into account the aperture variability at fracture scale, which is particularly present in natural fracture networks. Thus, thermal and solute tracer tests have been achieved to characterize heat transport in a fractured media at Ploemeur field site (SNO H +). To interpret the tracer tests, analytical expressions of thermal breakthrough peak retardation and amplitude have been developed for different fracture geometries : parallel plate fractures and channels. Those expressions are a powerful and innovative tool to characterize fracture geometries from thermal tracer tests, and also to predict thermal front transit time and lifetime of geothermal systems from solute tracer tests. Confrontation of those expressions to experimental results shows that observed differences between thermal and solute breakthrough can be explained only by channeling flow inducing low thermal transit times.

Les géostructures énergétiques consistent à établir un échange thermique direct avec le sol grâce à des systèmes intégrés dans les fondations ou les structures géotechniques. L’incorporation des échangeurs de chaleur aux géostructures provoque une variation cyclique de la température du sol adjacent. Des questions se posent sur l'impact de ces variations thermiques sur les paramètres géotechniques des sols en général, et en particulier des sols sensibles argileux. L’objectif de cette thèse est d’améliorer la compréhension et la quantification de l’impact de la variation de la température sur la capacité portante des pieux géothermiques. Actuellement, le dimensionnement des capacités portantes des fondations profondes est basé sur les résultats d’essais pénétrométriques ou pressiométriques. Des méthodes expérimentales ont été développées afin de permettre la réalisation de ces essais dans les conditions du laboratoire. Des essais mini-pénétrométriques sont réalisés sur des éprouvettes compactées à différents états initiaux et soumises à des températures variant de 1 à 70 °C. Les résultats montrent une évolution sensible des paramètres étudiés, la résistance en pointe (qc) et le frottement latéral (fs), pour un matériau illitique, lorsqu’il est compacté du côté sec de l’optimum. Les essais mini-pressiométriques, réalisés sur des massifs de sol illitique compactés en modèle réduit d’échelle métrique dans une cuve thermo-régulée, ont montré une diminution de la pression de fluage (pf) et de la pression limite (pl) avec l’augmentation de la température, tandis que la variation du module pressiométrique (EM) est moins marquée. Les résultats montrent une quasi-réversibilité des effets d’un cycle de chauffage dans la gamme de température testée alors que l’effet d’un cycle de refroidissement n’est que partiellement réversible. Pour les essais soumis à plusieurs cycles thermiques, le premier cycle induit des variations de paramètres toujours plus importantes que les cycles suivants. Une analyse approfondie de l’évolution des propriétés thermiques (la conductivité thermique (λ), la capacité thermique volumique (Cv) et la diffusivité thermique (D)) des sols en fonction de la teneur en eau, de la masse volumique sèche et de la température montre une augmentation de ces paramètres avec l’augmentation de w et ρd et une augmentation de λ des éprouvettes illitiques du côté sec de l’optimum avec l’augmentation de la température de 1 à 70 °C. En résumé, pour les pieux énergétiques, les résultats obtenus en laboratoire montrent une modification de la capacité portante due à la variation des paramètres du sol illitique sous l’effet des cycles thermiques / Energy geostructures involve providing a direct heat exchange with the ground through integrated systems in the foundations or geotechnical structures. The incorporation of heat exchangers in geostructures produces a cyclic variation of the temperature in the adjacent soil. Therefore, there are important scientific questions about the effect of temperature variations on hydro-mechanical soil parameters in general, and particularly for sensitive clay soils. The main objective of this thesis is to improve the understanding and the quantification of the impact of temperature variation on the bearing capacity of geothermal piles. Currently, the design of the bearing capacity of deep foundations is based on the results of in situ penetrometer and pressuremeter tests. Herein, experimental methods are developed to carry out these tests in laboratory conditions. Mini-penetrometer tests were carried out on samples compacted at different initial states and subjected to temperature variations ranging from 1 to 70 °C. The results showed a significant change in the studied parameters: the cone resistance (qc) and the friction sleeve resistance (fs) for an illitic material compacted on the dry side of the compaction curve. Mini-pressuremeter tests performed on the same illitic compacted soil in a thermo-regulated metric scale container, showed a decrease in creep pressure (pf) and limit pressure (pl) with increasing temperature, while the variation of pressuremeter modulus (EM) is less pronounced. The results showed a quasi-reversibility of the effect of a heating cycle through the temperature range tested, while the effect of a cooling cycle was only partially reversible. In the case of several thermal cycles, the first cycle induced more important parameter variations than the subsequent cycles, and at the end of the experimentation. Further analysis of the evolution of the thermal properties (thermal conductivity (λ), heat capacity (Cv) and thermal diffusivity (D)) within heating and cooling process as a function of soil water content and dry density showed an increase of these parameters with the increase of initial values of w and ρd, and an increase of λ in the dry side of the compaction curve with increasing temperature from 1 to 70 °C. In summary, for the energy piles driven in the clay soils, some modifications in the bearing capacity have to be taken into account due to the variation of the hydro-mechanical parameters of the soil induced by thermal cycles

L’objectif de notre étude est de modéliser un réservoir géothermique. Si nous supposons que le réservoir géothermique n’est composé que d’eau pure le transfert de matière et d’énergie est classiquement décrit par deux équations de conservation : la conservation de la matière et la conservation de l’énergie. À ces deux équations vient s’ajouter la vitesse du fluide classiquement donnée par la loi de Darcy tandis que les propriétés thermodynamiques, obtenues grâce à des équations théoriques ou empiriques (les équations d’état), ferment le modèle mathématique. Dès lors, ce modèle fermé, il existe différents schémas de résolution. Le premier est de résoudre en pression et température puis de procéder à un changement de variables lors du passage de monophasique à diphasique ou de diphasique à monophasique. TOUGH2 utilise le couple pression-saturation dans la zone diphasique.La seconde approche est de résoudre en pression et enthalpie afin d’accroître la stabilité lors de la transition entre l’état monophasique et l’état diphasique (voir Hydrotherm). Nous avons adopté la seconde option, résoudre en pression et enthalpie. De plus la résolution spatiale est faite avec les volumes finis.La modélisation d’un réservoir géothermique fait intervenir des équations fortement dépendantes l’une de l’autre. Cependant nous avons fait le choix de découpler la résolution afin de se libérer de la complexité de la résolution du système couplé. En effet, cette méthode possède l’avantage d’être moins consommatrice de mémoire puisque nous travaillons toujours avec le même nombre de données, mais dans une matrice deux fois moins importante. Nous montrerons que cette méthode demeure suffisamment précise pour une utilisation aussi bien dans le domaine industriel que dans celui de la recherche.Nous offrons à l’utilisateur une grande liberté grâce à l’implémentation de plusieurs méthodes : Euler implicite, explicite, Runge-Kutta ou BDF2 pour les solveurs temporels ou GMRES et BICGSTAB pour les solveurs linéaires. Nous pouvons gérer des conditions aux limites très variées telles que des flux nuls (décrivant une frontière qui n’échange pas de matière avec l’extérieur) ou une condition mixte (un Dirichlet sur la pression et un Dirichlet ou condition « sortie libre » sur la température… Cette dernière situation décrit une zone de recharge ou de décharge. Nous avons développé un outil multilangage : Python,Fortran et C++ (une implémentation de l’IAPWS provenant du projet freesteam incluant la zone supercritique). Tous ces langages sont orientés objet. L’IAPWS est l’outil permettant de calculer les propriétés physiques inconnues et par conséquent il ferme le système.Enfin nous avons appliqué le modèle sur le bassin parisien, France, sur plusieurs systèmes 1D et un autre système 2D réalisés par Coumou avec la plateforme CSMP++. Le bassin parisien est un réservoir exploité pour produire de la chaleur par le biais du pompage d’une eau à 70 _C et réinjecté à 40 _C. Les simulations 1D permettent de visualiser le déplacement d’un front de chaleur en haute enthalpie. La simulation 2D montre la convection naturelle de l’eau dans une faille. Chaque simulation a été comparée aux résultats obtenus avec un autre code (CSMP++, HYDROTHERM ou TOUGH2) et les résultats sont en accord. / The purpose of our study is to model a geothermal reservoir. When geothermal reservoirs are assumed to be composed of pure water, the transfer of mass and energy is classically described by two balance equations: the mass balance equation and the energy balance equation. In addition to those equations, fluid velocity is classically given by the Darcy law while thermodynamic properties, inferred from theoretical or empirical equations of state, are used to close the mathematical system. Once this system is closed, there exist different solutions. The first one is to solve for pressure and temperature with a variable switch to saturation in the two-phase region (e.g. TOUGH2). The second one is to solve for pressure and enthalpy to increase the stability of phase transition between single and two-phase states (e.g. Hydrotherm). We adopted the second option. We solve the system by using a splitting method — to get rid of the complexity of coupling equations — and a finite volume method for the spatial discretization. We offer some freedom to users thanks to the implementation of several methods like explicit or implicit Euler, Runge-Kutta or BDF2 for time solvers or GMRES and BICGSTAB for the linear solver. We can handle several boundary conditions like no-flow — describing a boundary which cannot exchange matter withthe exterior — or like a mixed-therm condition — a Dirichlet condition to the pressure and a Dirichlet or an outflow condition to the temperature in order to describe a recharge or a discharge zone — …Selecting object-oriented languages, we developed a multi-language framework, combining Python, Fortran and a C++ implementation of IAPWS (from the freesteam project) including the supercritical equations. To close the system physical propertiesare determined by the IAPWS- IF97 thermodynamic formulation. We’ve applied this simulation model to the dogger in Paris, France, to several onedimensional systems and a two-dimensional one made by Coumou with the CSMP++platform. The dogger is a reservoir exploited to produce heat by pumping water at 70 _C and reinjecting it in the reservoir at 40 _C. In the one-dimensional systems we wanted to observe the process of heat transfer from a higher temperature boundary to a smaller one in a high-energy domain. The last simulation shows the natural convection of water in a fault. For every simulation we compared the solutions we found with another code (TOUGH2 or CSMP++) and they agreed.