Le stockage géologique de CO2 est une des solutions techniques envisagées pour faire face aux émissions de gaz à effet de serre dans l'atmosphère. Déjà mis en oeuvre aujourd'hui dans quelques sites pilotes, son acceptabilité à grande échelle requiert la démonstration de la tenue à long terme des propriétés de confinement des couvertures sus-jacentes aux réservoirs, et des ouvrages en ciment, particulièrement sensibles, qui constituent notamment les scellements de fond de puits. Pendant la phase d'injection, une bulle de CO2 supercritique se forme, migre vers le sommet du réservoir sous l'effet de la poussée d'Archimède et se solubilise lentement (plusieurs dizaines d'années pour que la solubilisation soit complète). Alors que les réactions de dissolution et de carbonatation en présence d'eau liquide sont largement étudiées et commencent à être bien connues, la réactivité de la bulle de CO2 en phase supercritique ou vapeur est encore mal appréhendée. C'est face à ces incertitudes que nous avons mené une série d'expériences visant à décrire, à comprendre et à quantifier la réactivité d'un matériel minéral en présence d'une phase vapeur riche en CO2. Un minéral a été particulièrement étudié : la portlandite. Il s'agit d'un constituant clef des ciments riches en calcium, et à ce titre, peut, en première approche, être considéré comme un analogue des ciments de scellements de fond de puits. La cinétique de carbonatation de la portlandite a été déterminée à 80, 120 et 200°C, sous une pression de 160 bar et au contact de CO2 initialement anhydre, ou d'un mélange CO2-H2O. Un protocole original a été élaboré : il associe mesures de l'avancement de la réaction (principalement, analyses par DRX du taux de carbonatation de la fraction solide et enregistrement de la consommation de CO2 au cours de la réaction) et modélisation de l'évolution des équilibres de phases et des propriétés volumétriques du fluide réactionnel. D'un point de vue phénoménologique, nous avons montré que l'état de phase du système (vapeur ou biphasique liquide-vapeur) est déterminant. Il contrôle en grande partie le mode de réactivité et l'évolution cinétique de la réaction. Un modèle géochimique prenant en compte les spécificités de notre système réactionnel pauvre en eau a alors été développé pour simuler la carbonatation de la portlandite et à terme celle d'un ciment complet de type "fond de puits".

L'effet de serre lié à l'émission de CO2 a conduit à des projets de stockage de ce gaz dans des formations réservoirs. Ces formations peuvent être traversées de failles et l'examen de la sûreté du stockage nécessite alors l'étude du risque de réactivation et de propagation de ces failles. Cette étude passe par des investigations approfondies portant sur des conditions de propagation des fractures sous sollicitations hydromécaniques. Cette thèse a pour objectif l'étude théorique et numérique de ces conditions ainsi que la simulation numérique de la propagation. La modélisation numérique des processus thermo-hydro-mécaniques dans les milieux poreux fracturés par la méthode des éléments finis (MEF) permet de simuler des phénomènes complexes et non linéaires. Les difficultés liées à l'intégration des équations d'échanges de fluide entre la fracture et la matrice environnante avec la MEF ont été résolues dans des travaux récents et nos simulations numériques ont pu être basées sur cette méthode. Dans un premier temps, nous avons modélisé l'écoulement transitoire dans et au voisinage d'une fracture soumise à une injection de fluide et nous avons étudié le facteur d'intensité des contraintes (FIC) à l'extrémité de la fracture dans le cadre de la théorie de la poroélasticité linéaire. Si les conditions d'injection sont maintenues constantes et la fracture n'évolue pas, l'écoulement tend vers un état stationnaire. Le FIC évolue au cours de la phase transitoire pour atteindre une valeur limite dans l'état stationnaire. La modélisation de l'écoulement transitoire est très coûteuse en temps de calcul et il est intéressant de trouver un moyen d'exploiter au mieux les résultats d'un calcul en état stationnaire. L'analyse théorique et les résultats des simulations numériques montrent en effet que le FIC calculé à l'état stationnaire peut fournir certaines bornes pour la propagation des fractures sous l'écoulement transitoire. Dans le cadre de la poroélasticité linéaire et de l'écoulement de Poiseuille dans les fractures, des expressions semi-analytiques pour le FIC à l'état d'écoulement stationnaire ont pu être dérivées. Pour des géométries simples, ces formules approximatives se révèlent efficaces pour discuter des conditions de propagation des fractures pour des cas typiques et simples de géométrie de la fracture et des conditions d'injection de fluide. Dans un deuxième temps, un Modèle de Fracture Cohésive (MFC) a été utilisé pour modéliser la propagation de fracture sur la base de l'endommagent. Ce modèle, basé sur un critère de rupture de Mohr–Coulomb modifié, permet de simuler l'endommagement de l'interface à la fois sous sollicitations en mode I et II. Une relation d'équivalence entre les paramètres de ce modèle et du modèle de Mécanique Linéaire de la Rupture (MLR) a été établie sur la base de la longueur de propagation de fracture sous des charges similaires. Cette relation permet l'extension de l'équivalence théorique entre MLR et MFC établie pour les matériaux fragiles et sur la base de critères énergétiques, à des matériaux quasi-fragiles et ductiles. On a d'ailleurs montré que le MFC permet de simuler certains phénomènes spécifiques tels qu'instabilités de propagation en mode I et II et le branchement de la fracture en mode II. Enfin, la prise en compte de la pression de fluide dans la fracture a permis d'obtenir un modèle de MFC couplé avec l'hydraulique qui a été implémenté dans un code numérique aux éléments finis en vue d'étudier la propagation des fractures sous sollicitations hydromécaniques. Des simulations numériques ont été réalisées afin d'étudier le risque de réactivation et de propagation des failles dans le contexte de stockage du CO2 en particulier dans une configuration de formation réservoir du type Bassin de Paris / Global warming effect related to CO2 emission has led to sequestration projects of this gas in reservoir formations. These formations can be crossed by faults and safety issue of storage requires the study of fault reactivation and propagation risk. This study goes through in-depth investigations of fracture propagation conditions under hydromechanical solicitations. This thesis aims at theoretical and numerical studies of these conditions and the numerical simulation of fracture propagation. Numerical modelling of thermo-hydro-mechanical processes in fractured porous media using Finite Element Method (FEM) allows the simulation of complex and nonlinear phenomena. Difficulties in integrating fluid mass exchange between fracture and surrounding matrix in the equations with FEM have been solved in recent works and our numerical simulations have been based on this progress. In a first step, we modelled transient flow subjected to a fluid injection and we studied the Stress Intensity Factor (SIF) at fracture tip in the framework of linear poroelasticity theory. If injection conditions are kept constant and the fracture does not evolve, the flow tends to a steady state. The SIF develops during transient phase to reach a limit value in the steady state. Modelling of transient flow is very time consuming and it is interesting to find a method to exploit the results of a calculation in steady state. Theoretical analysis and results of numerical simulations show that the SIF calculated at steady state can provide some bounds for fracture propagation under transient flow. In the framework of linear poroelasticity and Poiseuille flow in fractures, some semi-analytical expressions of SIF at steady state could be derived. For simple geometries, these approximate formulations are efficient to discuss fracture propagation conditions for typical and simple cases of fracture geometry and fluid injection conditions. In a second step, a Cohesive Zone/Fracture Model (CFM) was used to model fracture propagation on the basis of damage. This model, based on a modified Mohr-Coulomb failure criterion, simulates interface damage under both mode I and II loads. An equivalence relation between parameters of CFM and Linear Elastic Fracture Mechanics model (LEFM) was established on the basis of fracture propagation length under similar loads. This relationship allows the extension of theoretical equivalence between LEFM and CFM established for brittle materials and on the basis of energy criteria, for quasi-brittle and ductile materials. It has also shown that CFM can simulate specific phenomena such as propagation instabilities for mode I and II and fracture kinking under mode II. Finally, taking into account the fluid pressure in the fracture permitted to obtain a CFM coupled with hydraulic processes which has been implemented in a numerical finite element code to study fracture propagation under hydromechanical solicitations. Numerical simulations were performed to study the risk of fault reactivation and propagation in the context of CO2 injection in Paris Basin reservoir formation