Le stockage de CO2 dans les formations géologiques profondes est considéré comme un moyen envisageable de réduire les émissions de gaz à effet de serre dans l'atmosphère. L'état des roches constituant le réservoir est un paramètre clé dont dépendent aussi bien la sûreté que l'efficacité du stockage. L'objectif de cette thèse est de caractériser les risques générés par une modification possible des propriétés mécaniques et de transfert des matériaux du sous-sol suite à une injection de CO2. Des simulations à grande échelle visant à représenter le processus d'injection de CO2 à l'état supercritique dans ut réservoir souterrain ont été réalisées. Une analyse des champs de contraintes obtenus met en évidence la possibilité de rencontrer pour d'importants débits d'injection divers régimes de détérioration des matériaux. Le travail est consacré à l'étude de l'apparition de fissures ouvertes. A la suite d'une étude analytique simplifiée de la formation et de ta propagation de fissures ouvertes dans le milieu considéré basée sur un modèle probabiliste on montre que la formation d'un réseau de fissures est possible. On s'attache alors à développer dans te Code_Aster un outil numérique permettant de simuler la formation de réseaux de fissures. Un modèle non local basé sur une régularisation des contraintes est proposé. Un critère portant sur le facteur d'intensité des contraintes est utilisé pour décrire la propagation des fissures. L'amorçage de nouvelles fissures est décrit par un modèle de Poisson-Weibull. Les paramètres utilisés sont identifiés par une campagne expérimentale menée sur des échantillons prélevés sur un site de stockage géologique de CO2. Le modèle développé est ensuite validé sur des cas numériques, puis par rapport à la campagne expérimentale réalisée.

Dans le cadre d’un programme de recherche mené par le GISOS (Groupement d’Intérêt Scientifique de Recherche sur l’Impact et la Sécurité des Ouvrages Souterrains), le site pilote de Cerville-Buissoncourt (Lorraine, France) a fait l’objet d’une importante instrumentation géophysique et géotechnique pour assurer la surveillance d’une cavité saline à 200 m de profondeur, depuis son état stationnaire jusqu’à l’effondrement des terrains du recouvrement. Les objectifs principaux de cette thèse consistaient à : 1) valider la technique de surveillance basée sur l’écoute microsismique dans un contexte salin, et 2) modéliser numériquement le comportement mécanique complexe du recouvrement, particulièrement l’initiation des microfissures et leur propagation. L’analyse de la sismicité induite enregistrée a permis de caractériser l’état initial de la cavité en terme d’activité microsismique. Deux types d’événements ont été identifiés : les événements isolés correspondant aux ruptures localisées, et les événements en rafale, d’une dizaine de secondes de durée. D’après les résultats de localisation d’hypocentres, la totalité de la sismicité enregistrée est générée au niveau de la cavité dans le gisement de sel, ou bien dans les faciès marneux qui composent le toit immédiat de la cavité actuelle. Les déclenchements en rafale seraient liés à des phénomènes de délitement puis de décrochement de blocs de marne, suivis des chutes de blocs dans la cavité remplie de saumure. Le travail de modélisation numérique a été focalisé sur la possibilité de rendre compte de l’endommagement dans les couches fragiles du recouvrement. Nous avons mis en oeuvre un modèle géomécanique hybride à l’échelle du site pilote qui intègre les différentes formations géologiques présentes dans le recouvrement, ainsi que l’initiation, la propagation et la coalescence des microfissures dans le banc raide, à l’aide des logiciels FLAC et PFC2D. La calibration du modèle discret PFC2D pour reproduire le comportement en traction du banc raide a été vérifiée numériquement à l’échelle du site pilote. Cette vérification a été basée sur la comparaison, en termes de la réponse élastique et d’apparition des ruptures dans le banc raide, entre l’approche hybride FLACPFC 2D et la modélisation purement continue avec FLAC. Le modèle hybride ainsi défini pourra être utilisé dans le cadre d’une retro-analyse une fois que les mesures in-situ, notamment les enregistrements microsismiques et les données de déformation, seront disponibles à Cerville-Buissoncourt / Within the framework of a research program carried out by the GISOS (Scientific Grouping of Research Interest on the Impact and Safety of Underground Works), the pilot site of Cerville-Buissoncourt (Lorraine, France) was the subject of a large geophysical and geotechnical instrumentation to ensure the monitoring of a salt cavity at 200 m depth, from its stationary state to the final overburden collapse. The main objectives of this work consisted on : 1) the validation of the microseismic monitoring technique in a salt mine environment, and 2) the numerical modelling of the mechanical behavior of the overburden, particularly the initiation and the propagation of microcracks. The analysis of the recorded induced seismicity allowed to characterize the initial state of the cavity in terms of microseismic activity. Two types of events were identified : isolated events corresponding to localized ruptures, and swarms of events, of tens of seconds of duration. According to hypocenter location results, the totality of the recorded seismicity is generated either in the cavity surroundings within the salt layer, or in the marly facies of the current cavity roof. Swarms would be related to delamination of clayley marls in the immediate roof, followed by rock debris falling in the brine filled cavity. The numerical modelling was focused on the possibility of accounting for the damage in the fragile layers of the overburden. We implemented a hybrid geomechanical model of the pilot site which integrates the various geological formations present in the overburden, as well as the initiation, the propagation and the coalescence of microcracks in the stiff layer, using FLAC and PFC2D softwares. The calibration of the discrete PFC2D model to reproduce the tensile behaviour of the stiff layer was numerically checked on the site scale. The validation was based on the comparison, in terms of the elastic response and the damage onset in the stiff layer, between the hybrid approach FLAC-PFC2D and the purely continuous modelling using FLAC. The hybrid model thus defined would be used for back-analysis studies once in-situ measurements, in particular microseismic recordings and deformation data, will be available at Cerville-Buissoncourt