Le transport des composants chimiques dans le sol par l'eau est une cause principale de la pollution du sol. Ce transport a lieu autour des remblais et de la zone de stockage du déchet nucléaire, des rejets et des pertes de mine, etc. Un grand nombre de ces sites est non saturé de l'eau et dans certains cas la chaleur peut changer le destin des espèces chimiques, qui nous mènent à un problème couplé dans le milieu multiphasique. Dans ce mémoire, les simulations numériques avec un modèle thermo-hydro-mécanique existant et les modélisations théoriques et les simulations numériques du transport et les interactions d'une espèce chimique dans les milieux multiphasiques sont présentées. Les formulations THM intégrées dans le Code_Aster sont présentés. L'excavation, la mise en place de barrière ouvragée et le chargement thermique de puits de stockage des déchets nucléaires sont modélisés. La vérification du modèle est presentée à travers de ces simulations numériques. Une modélisation théorique du comportement thermo-hydro-mécanique couplée aux phénomènes chimiques est ensuite présentée dans le cadre d'une approche entièrement couplée dans laquelle l'eau, le gaz, l'espèce chimique et le squelette de sol ont été considérés comme constituants; les inconnues correspondantes sont la température, la pression de l'eau, la pression de gaz, la concentration chimique et les déplacements. Pour chaque constituant, l'équation de bilan de masse et la loi de comportement appropriée sont écrites et résolues simultanément. Les résultats de ce modèle ont été comparés à des résultats théoriques et expérimentaux existant dans la littérature. Les résultats de quelques applications du modèle sont également présentées. Enfin, quelques perspectives de recherche sont proposées. En particulier, il y a un besoin d'effectuer des expériences pour obtenir les paramètres nécessaires du sol pour permettre de modéliser plus précisément le transport de la chaleur et du soluté dans les sols non saturés et déformables.

L'étude du comportement du béton en température est un problème majeur qui vise in fine à évaluer le niveau de sécurité des structures sous des sollicitations thermiques sévères, lors d'incendies par exemple. Pour cela, de nombreux modèles sont développés dans un cadre de couplage thermo-hydro-mécanique (THM), pour prendre en compte les différents processus physico-chimiques et mécaniques mobilisés par ces sollicitations et conditionnant la stabilité du matériau béton.Cependant, ces modélisations ne prennent souvent pas en compte explicitement la nature hétérogène du matériau. En effet, le béton est un matériau constitué d'inclusions noyées dans une matrice cimentaire pouvant être schématisée comme un milieu poreux ouvert partiellement saturé en eau. Les inclusions sont caractérisées par leurs natures minéralogiques, leurs morphologies et leurs tailles classées en fuseaux granulaires. Cette hétérogénéité introduit une complexité supplémentaire : la nécessité de prendre en compte la microstructure pour quantifier l'effet de l'incompatibilité (thermique, hydrique et mécanique) inclusion-matrice sur le comportement THM du béton. Ce travail constitue une première avancée dans ce sens. A ce titre, une modélisation élément fini multi-échelle tridimensionnelle (3D) est développée permettant d'affecter des comportements spécifiques à la matrice et aux inclusions. Pour la matrice, siège des transports de masse dans son réseau poreux connecté, une approche thermo-hydrique à trois fluides (eau, vapeur et air sec) est adoptée et est couplée à une loi de comportement poro-mécanique endommageable. Les inclusions (granulats) sont considérées hydriquement inertes, une approche thermo-mécanique avec endommagement est alors retenue.Une analyse, par simulations numériques, des effets de la nature minéralogique des granulats (calcaires ou silico-calcaires) de leurs distributions et de leurs morphologies a été menée sur des configurations 2D et 3D. Les effets étudiés ont notamment concerné l'influence de ces paramètres sur les fluctuations locales des champs de température, de pression de gaz et d'endommagement au regard de la dispersion des mesures expérimentales. L'analyse est limitée à l'échelle mésoscopique, celle où les granulats de taille caractéristique supérieure à 5 mm sont considérés, le reste (stable et pâte de ciment) étant une matrice homogène. Enfin, cette analyse a mis en évidence le besoin de mettre en œuvre une approche expérimentale cohérente avec une analyse mutli-échelle, à la fois pour la caractérisation des propriétés (thermiques, hydriques et mécaniques) de chaque constituant et pour l'étude des évolutions des champs lors des changements d'échelles. Un protocole expérimental a été définit et des premiers résultats d'essais sont présentés et analysés au travers de résultats obtenus dans la partie modélisation

Les matériaux soumis à de grandes déformations présentes pour la plupart l’apparition de déformations inélastiques. Ce phénomène est souvent accompagné d’une localisation des déformations dans une zone étroite, précurseur de la rupture. Un cas particulier, mais très fréquent, est les bandes de cisaillement qui apparaissent pour beaucoup de géomatériaux. Ces bandes peuvent être rencontrées à des échelles allant de l’échelle kilométrique pour les zones de subduction à l’échelle micrométrique à l’intérieur des zones de faille. Etudier et modéliser la création de ces zones d’instabilité est fondamental pour décrire la rupture des géomatériaux et des phénomènes associés comme les glissements sismiques dans les zones de faille mature de la lithosphère. Les conditions de pression, de température, l’interaction de l’eau interstitielle avec un matériau finement fracturé conduisent à l’apparition de multiples processus physiques impliqués dans les glissements sismiques. Dans ce travail, nous nous attachons à modéliser la création de bandes de cisaillement à l’intérieur des gouges de faille en prenant en compte l’effet de la microstructure par l’intermédiaire des milieux continus de Cosserat, ainsi que les couplages thermo-hydro-mécanique. L’utilisation de la théorie de Cosserat permet non seulement de régulariser le problème de localisation des déformations par l’introduction d’une longueur interne dans les lois constitutives, mais en même temps de prendre en compte l’effet de la microstructure. Deux approches sont employées pour étudier le système d’équations couplées aux dérivées partielles non linéaires : L’analyse de stabilité linéaire et la méthode des éléments finis. L’analyse de stabilité linéaire permet d’examiner les conditions d’apparitions d’instabilités pour un système mécanique avec des couplages multi-physiques. Par ailleurs, des considérations sur les perturbations appliquées au système permettent aussi de déterminer l’épaisseur de la zone de cisaillement, un paramètre clé pour la compréhension du mécanisme mécanique des failles. Ces estimations sont confirmées par l’intégration numérique pour des déformations restant dans une gamme donnée. Elles sont confrontées aux observations expérimentales et in situ et présentent une bonne corrélation. D’autre part, les simulations numériques permettent d’obtenir la réponse mécanique de la gouge de faille et de donner des informations sur l’influence des différents couplages dans le budget énergétique d’un tremblement de terre / When materials are subjected to large deformations, most of them experience inelastic deformations. It is often accompanied by a localization of these deformations into a narrow zone leading to failure. One particular case of strain localization is the formation of shear bands which are the most common patterns observed in geomaterials. In geological structures, they appear at very different scales, from kilometer scale for subduction zones, to micrometric scale inside fault cores. Studying their occurrence and evolution is of key importance to describe the failure of geomaterials and model seismic slip for mature crustal faults. The pressure and temperature conditions in these faults and the interaction with the pore water inside a highly fractured materials highlight the importance of different physical processes involved in the nucleation of earthquakes. In this thesis, we study the occurrence and evolution of shear bands inside fault gouges taking into account the material microstructure by resorting to elastoplastic Cosserat continua and also the effect of thermo-hydro mechanical couplings. The use of Cosserat theory introduces information about the gouge microstructure, namely the grain size, and permits to regularize the mathematical problem of in the post-localization regime by introducing an internal length into the constitutive equations. Two approaches are used to study the coupled non-linear partial differential set of equations: linear stability analysis and finite element simulations. Linear stability analysis allows to study the occurrence of localized deformation in a mechanical system with multi-physical couplings. Considerations on the dominant wave length of the perturbations permit also to determine the width of the localized zone. This shear band thickness is confirmed by numerical integration in the post-localization regime for a certain range of deformation. The obtained widths of the localized zone are key parameters for understanding fault behavior, are in agreement with experimental and field observations. Moreover, numerical finite element computations enable to model the mechanical response of a fault gouge during seismic slip and give insights into the influence of various physical couplings on the energy budget

Les roches argileuses profondes de très faible perméabilité (10-20 m2), telles que l’argilite du Callovo-Oxfordien (COx) en France ou l’Argile à Opalinus en Suisse, sont des roches hôtes potentielles pour le stockage géologique des déchets radioactifs. Lors des différentes phases du stockage, ces roches seront soumises à des sollicitations thermo-hydro-mécaniques (THM) couplées. La détermination de leurs caractéristiques THM reste à compléter pour une meilleure compréhension de la réponse du champ proche des galeries de stockage.L’étude expérimentale des matériaux de faible perméabilité est délicate et plusieurs questions sur leur comportement THM restent posées. Dans ce contexte, deux systèmes expérimentaux originaux avec mesures précises des déformations locales et chemin de drainage (H) réduit ont été développés et utilisés pour la réalisation d’essais saturés drainés : une cellule de compression isotrope (H = 10 mm) et une cellule triaxiale standard dont le chemin de drainage a été réduit à l’aide d’un géotextile placé autour de l’échantillon (H = 19 mm) .Une caractérisation poroélastique détaillée de l’argilite du COx en cellule isotrope a permis, par des approches directes et indirectes compatibles, une détermination fiable des paramètres poroélastiques isotropes transverses du matériau, identifiés dans un cadre théorique permettant de déterminer les composantes du tenseur de Biot (coefficients b1 et b2, respectivement perpendiculaire et parallèle au litage). Un paramètre clé peu documenté à ce jour est le module de compression des grains solides, déterminé à l’aide d’un essai de compression sans membrane (Ks = 21.7 GPa).Un contrôle de température et un étalonnage précis de ces effets parasites ont permis l’étude de la réponse volumique thermique de l’argilite du COx lors d’un essai de chauffage drainé sous confinement isotrope constant proche de l’état in-situ. Une expansion thermoélastique suivie d’une contraction thermoplastique a été observée, avec une température de transition à 48°C, proche de la plus grande température supportée pendant l’histoire géologique de l’argilite. On confirme ainsi le fait que les argilites gardent en mémoire la température maximale supportée. La précision des mesures de déformations a aussi permis d’identifier un fluage volumique, dont l’amplitude est accentuée à 80°C.Des calculs poroélastiques avec les conditions aux limites du système de drainage amélioré ont permis de déterminer l’ordre de grandeur des taux de déformations axiales permettant un bon drainage lors du cisaillement triaxial (6.6×10-8 s-1). Les essais drainés sur l’argilite du COx présentent une bonne compatibilité avec un critère déjà publié. Des résultats cohérents ont aussi été obtenus sur le critère de rupture de l’Argile à Opalinus, par rapport auquel ont été comparés ceux d’essais de compression simple à différentes vitesses sur des échantillons équilibrés à 94% d’humidité relative qui a mis en évidence une dépendance de la résistance au pic vis-à-vis du taux de déformation axial.L’ensemble de ces résultats permet de réduire l’incertitude concernant les propriétés THM des argilites et devrait permettre une meilleure estimation de la réponse du champ proche des galeries au cours des différentes phases de leur période service.Mot clés : argilites, perméabilité, drainage, saturation, poroélasticité, coefficients de Biot, isotropie transverse, cisaillement triaxial, chauffage, fluage / Deep low permeability claystones (10-20 m2), such as Callovo-Oxfordian claystone (COx) in France or the Opalinus Clay in Switzerland, are potential host rocks for deep geological radioactive waste disposal. During the various phases of the storage, these rocks will be subjected to thermo-hydro-mechanical (THM) coupled effects. The determination of their THM parameters remains to be completed for a better understanding of the near-field response of the storage galleries.The experimental study of low permeability geomaterials is difficult and several questions about their THM behavior still remain. In this context, two original experimental systems with high precision local strain measurements and reduced drainage lengths (H) were developed and used for saturated drained tests: an isotropic compression cell (H = 10 mm) and a standard triaxial cell with a reduced drainage length using a geotextile placed around the sample (H = 19 mm).A detailed poroelastic characterization of the COx argillite in the isotropic cell provided a set of compatible transverse isotropic poroelastic parameters of the material, identified in a theoretical framework allowing to determine the Biot tensor components (coefficients b1 and b2, perpendicular and parallel to bedding plans, respectively). A key parameter not well documented to date is the unjacketed modulus that was determined by means of an unjacketed compression test (Ks = 21.7 GPa).A temperature control and an accurate calibration of thermal parasite effects allowed the investigation of the thermal volumetric response of the COx argillite during a drainage test under constant isotropic confining stress close in-situ state conditions. A thermoelastic expansion followed by a thermoplastic contraction was observed, with a transition at a temperature of 48 °C, close to the highest temperature supported during the geological history of the claystone. It is thus confirmed that such claystones keep in memory the maximum supported geological temperature. The precision of the deformation measurements also made it possible to identify a volumetric creep that is enhanced at 80 °C.Poroelastic calculations with the boundary conditions of the improved drainage system allowed to determine the magnitude of axial strain rates allowing good drainage during triaxial drained shearing (6.6×10-8 s-1). The drained tests carried out on the COx claystone showed a good compatibility with a criterion already published. Coherent data were also obtained on the Opalinus Clay failure criterion, that were compared to those of uniaxial compression tests at different speeds on samples equilibrated at 94% relative humidity, that exhibited a dependence of the peak strength on the shear rate.These results make it possible to reduce the uncertainties concerning the THM properties of claystones and should allow a better estimation of the response of the near field close to the galleries during the different phases of their service period.Key words: claystone, permeability, drainage, saturation, poroelasticity, Biot coefficient, transverse isotropy, triaxial testing, heating, creep

Prolonger la durée de vie d’un parc nucléaire impose de garantir la pérennité des enceintes de confinement, troisième et dernière barrière de protection contre les potentielles émissions radioactives. Ces structures sont attentivement auscultées depuis la construction grâce à un dispositif de mesures spécifiques incluant des extensomètres noyés dans le béton. De conception intrinsèquement fiable et bien que largement éprouvés, ces capteurs de déformation commencent à présenter des taux de défaillance significatifs, nécessitant, de fait, le recours à des chaînes de mesure de substitution. En ne pouvant être mises en œuvre qu’au voisinage de la surface du béton, ces solutions ne permettent qu’un accès indirect à la déformation initialement observée. Des doutes peuvent dès lors être émis sur la représentativité des mesures et par conséquent sur la qualité de la surveillance.Pour analyser l'influence de la localisation de la mesure, le développement et l'exploitation de deux modèles numériques éléments finis complémentaires sont proposés. Il s’agira principalement, pour l’un, de quantifier les possibles écarts de comportement sur le long terme et en épreuve liés aux gradients de température et d’humidité ; et pour l’autre, la variabilité liée au caractère hétérogène et particulier de la géométrie des enceintes. Suite à la calibration des modèles de comportement thermo-hydro-mécanique, s’appuyant sur des résultats d'essais matériaux en laboratoire, une prévision à l'échelle de la structure, utilisant des chargements sur ouvrage reconstruits à partir de données historiques ou de synthèses fondées sur le retour d’expérience du parc nucléaire français exploité par EDF, est menée. L'historique de déformation est ainsi reconstruit en tout point de l'enceinte pour en déduire un modèle analytique des variations spatiales dans l'enceinte des déformations afin de faciliter son exploitation.En mettant en évidence de faibles écarts sur le long terme entre position à cœur et en surface, il est d'abord démontré que des chaînes de mesure de substitution sont utilisables avec un bon niveau de confiance en dehors des zones singulières. Par ailleurs, l’usage de fonctions de correction basées sur les transferts physiques régnant dans la structure permet de minimiser l’influence des changements de position de mesure dans l’épaisseur. / Extending the life of nuclear plant needs to ensure the sustainability of containment, third and last barrier against potential radioactive emissions. These structures are closely monitored from construction through a specific measuring system including strain gauges embedded in the concrete. Inherently reliable, though widely proven technology, these strain sensors start presenting significant failure rates, requiring the use of substitution measure chains. These solutions, which can only be placed near surface of concrete, allow only an indirect access to the initially observed strain. Doubts may therefore be issued on the representativeness of measurements and therefore on the quality of monitoring.To analyze the influence of the location of the measurement, the development and exploitation of two finite elements numerical models are proposed. This will be mainly for one to quantify the possible misbehavior in the long term and during tests related to temperature and moisture gradients; and for the other, the variability related to heterogeneous and particular nature of the geometry of nuclear power plant. As a result of the calibration of the thermo-hydro-mechanical behavior models, based on laboratory materials test, a forecast throughout the structure is conducted, using reconstructed loads from historical data or syntheses based on feedback from the french nuclear power plant operated by EDF. The history of deformation is thus reconstructed in any enclosure point to settle an analytical model of the spatial variations of the enclosure deformation to facilitate its operations.Highlighting small long-term differences between embedded and surface position, it is first shown that alternative measuring systems can be used with a high confident level outside singular position. Furthermore, the use of correction functions based on physical transfers prevailing in the structure makes it possible to minimize the influence of position measurement of changes in thickness.

Le présent travail est consacré à l'étude du comportement thermo-hydro-mécanique linéaire et non linéaire des roches poreuses de type argilites par approche de changement d'échelle. A partir des observations microstructurales de ces matériaux, un modèle conceptuel a été proposé. Dans ce modèle, le volume élémentaire représentatif du milieu hétérogène est composé d'une phase matricielle argileuse contenant des inclusions sphériques de minéraux de quartz et de calcite et des inclusions ellipsoïdales aplaties représentant l'espace poreux. Dans un premier temps, le procédé de la modélisation a été exploité par la détermination des propriétés effectives isotropes et isotrope transverses des argilites : la conductivité thermique et les propriétés thermo-hydro-mécanique croisées. En outre, de nombreuses études numériques ont mis en évidence l'influence de la morphologie de l'espace poreux, de la minéralogie et des schémas d'estimation sur les résultats prédictifs. Dans un deuxième temps, nous avons modélisé le comportement mécanique non linéaire (élasto-plastique, élasto-viscoplastique) des roches argileuses. La comparaison entre les simulations numériques et les résultats expérimentaux disponibles (essai de compression triaxiale, essai de fluage) a confirmé la validation du modèle développé / The present work deals with the linear and non-linear thermo-hydro-mechanical behaviour of porous rocks such as the argillite by the multiscale modelling approach. Based on microstructure observations, a conceptual model was proposed. In this model, the representative elementary volume of a heterogeneous medium is composed of an argillaceous matrix containing spherical inclusions of minerals quartz and calcite and ellipsoidal inclusions representing the pore space. In a first step, the process of modelling has been exploited by determining the isotropic and transversely isotropic effective properties of the argillite: thermal conductivity and thermo-hydro-mechanical properties. Furthermore, many numerical studies have highlighted the influence of the morphology of the pore space, of the mineralogy and of the estimate schemes to the predictive results. In a second step, we modelled the non linear mechanical behaviour (elasto-plastic, elasto-viscoplastic) of argillaceous rocks. The comparison between numerical simulations and available experimental results (triaxial compression test, creep test) confirmed the validation of the model developed

Après une dernière phase de forage et de stimulation hydraulique, le site expérimental de Soultz-sous-Forêts comporte désormais trois puits connectés à un réservoir géothermique profond (5 km). A l'avenir, du fait de la mise en place d'une unité de production électrique, seule la modélisation numérique renseignera sur l'état du réservoir. Une compréhension poussée des propriétés géométriques de ce dernier est donc incontournable. C'est ce que propose ce travail, notamment en étudiant la prise en compte de la zone endommagée dans les modèles de faille.<br />Au cours de cette étude, les propriétés physiques de la zone endommagée (porosité, perméabilité, surface spécifique, conductivité thermique) ont été mesurées à différentes échelles et sur des échantillons présentant des structures et des faciès d'altération variés. La combinaison de ces mesures, et leur intégration dans différents modèles numériques, a notamment permis de mettre en évidence l'impact de la zone endommagée sur les transferts de chaleur et de matière, ainsi qu'une évolution particulière de la géométrie du réseau poreux avec l'altération.

Le mouvement de Séchilienne est situé à une trentaine de kilomètre au sud-est de l'agglomération grenobloise sur la bordure Sud du Massif de Belledonne. Il affecte le versant rocheux en rive droite de la vallée alpine de la Romanche façonné par les nombreux cycles de glaciation et déglaciation et qui lors de la dernière phase glaciaire a été complètement recouvert par le glacier. Le raidissement des pentes et l'approfondissement du pied du versant associés au phénomène de décompression affectant le versant et correspondant au relâchement des contraintes liées à la disparition progressive du glacier seraient donc les facteurs prédominants dans l'initiation et l'évolution du mouvement de Séchilienne. Les récentes contraintes temporelles obtenues sur l'escarpement sommital ont cependant montré une initiation différée du mouvement par rapport à la déglaciation et suggèrent une initiation lors de l'Optimum climatique de l'Holocène suite à une modification brusque du climat. L'initiation et l'évolution de grands mouvements de versants rocheux en domaine montagneux peuvent être le résultat de la contribution de différents facteurs géologiques (la fracturation, la lithologie ou le régime tectonique), géomorphologiques (relief, pente) mais également de facteurs de forçage externe tels que les séismes, le pergélisol, les variations hydrologiques, le retrait du glacier ou encore les effets climatiques. Ces facteurs étant tous potentiellement actifs sur le mouvement de Séchilienne, l'identification des facteurs contrôlant l'initiation et l'évolution du mouvement et pouvant expliquer le délai observé entre la déglaciation et l'initiation du mouvement est donc délicate.Nous avons donc pour cela mené une étude multidisciplinaire permettant à partir (1) des données géophysiques et d'essais en laboratoire de déterminer en profondeur les limites du mouvement et d'évaluer le volume déstabilisé ainsi que son état de déstructuration, (2) d'une analyse de la structure détaillée du versant en surface et en profondeur de préciser l'origine des fractures observées mais également leur rôle de contrôle sur le mouvement (3) d'une analyse géomorphologique des objets présents à la surface du versant, d'identifier de discrétiser et de sélectionner les structures glaciaires et gravitaires sur lesquels nous avons par la suite appliqué la méthode de datation permettant ainsi d'obtenir le modèle cinématique du retrait du glacier de la Romanche le long du versant de Séchilienne et de l'initiation et de l'évolution du mouvement (4) d'une modélisation thermique du versant basée sur la courbe d'évolution des températures reconstruite depuis la dernière phase glaciaire quaternaire de caractériser la géométrie et la persistance du pergélisol depuis la dernière phase glaciaire jusqu'à la période actuelle.La confrontation de l'ensemble de ces résultats et leur replacement dans le contexte climatique régnant depuis la dernière phase glaciaire il ya 21000ans jusqu'à l'actuelle permet ainsi de proposer un modèle conceptuel de l'initiation et de l'évolution du mouvement de Séchilienne et du rôle des différents facteurs étudiés. Le mouvement se serait donc initié suite à la décompression du versant en réponse au retrait du glacier. La mise en place et le développement du pergélisol au sein du versant aurait cependant limité la propagation du mouvement en renforçant mécaniquement le versant. Son action ne c'est cependant pas limitée à renforcer le versant mais aurait également sur une échelle de temps plus importante favoriser l'endommagement et l'affaiblissement mécanique de la zone affecté imposant ainsi au mouvement sa géométrie actuelle. / The Sechilienne landslide is located thirty kilometers at the Southeast of Grenoble on the Southern edge of the massif of Belledonne. It affects the rock slope on the right bank of the Romanche alpine valley shaped by many cycles of glaciation and deglaciation and which was completely covered by the glacier during the last glacial phase. The oversteepening of the slopes and the deepening of the foot of the slope associated with the phenomenon of debutressing affecting the slope and corresponding to the lateral stress release resulting from ice melting would be so predominant factors in the initiation and development of the Sechilienne landslide. However the recent temporal constraints obtained on the scarp showed a delayed motion relative to the initiation and deglaciation suggest initiation during the Holocene climatic Optimum after a sudden climate change. The initiation and evolution of large movements of rock slopes in mountainous area may be the result of the contribution of different geological factors (fracturing , lithology or tectonic regime) , geomorphology (relief, slope) but also external forcing factors such as earthquakes , permafrost , hydrological changes , the glacier disappearance or climatic effects. These factors are all potentially active in the Sechilienne landslide, so the identification of factors controlling the initiation and development of the movement and which can explain the delay observed between deglaciation and the initiation of movement is tricky.Therefore, we conducted a multidisciplinary study from (1) geophysical data and laboratory tests to determine the depth limits of the mass movement and to assess the volume and the breakdown state, (2) an analysis of the detailed structures of the slope at the surface and in depth to identify the origin of the observed fractures and also their role in the control of the movement (3) a geomorphological analysis of objects on the surface of the slope to discretize , to identify and to select the glacial and gravitational structures on which then we applied the datation method to obtain the kinematic model of the Romanche glacier retreat along the slope of Séchilienne and of the initiation and the evolution of the movement (4) a thermal modeling of the slope based on the temperature curve proposed since the last quaternary glacial stage to characterize the geometry and the persistence of permafrost since the last glacial phase to the current period.The comparison of all these results and their replacement in the climatic conditions prevailing since the last glacial stage 21000 years ago to the present, allows to propose a conceptual model for the initiation and the evolution of Séchilienne landslide and the role of the various studied factors. So, the landslide would be initiated by slope debutressing in response to the glacier retreat. However the establishment and the development of the permafrost within the slope would limit the propagation of the movement by reinforcing mechanically the slope.

Le cadre général de la thèse concerne la valorisation sous forme de production d'énergie électrique de la chaleur présente à quelques kilomètres de profondeur (3 à 5 km), en général dans des milieux peu perméables et fracturés. Notre objectif principal est d'étudier le phénomène des microséismes induits relativement au refroidissement, en nous basant sur une expérience de terrain de longue durée, menée sur le site de Rosemanowes (Cornwall, UK). Pour cela nous avons procédé à la mise en place d'un outil de calcul, FRACAS, capable de simuler ce phénomène en introduisant une approche à double milieu thermique pour mieux simuler le refroidissement du réservoir dû à l'injection de fluide à long terme, responsable des nouveaux mécanismes de ruptures dus à la traction de la roche. Dans ce contexte nous avons introduit un nouvel algorithme pour prendre en compte les manifestations d'instabilités, un mécanisme de « stick-slip » avec prise en compte d'une friction statique et d'une friction dynamique. La possibilité d'induire des microséismes est ensuite étudiée à partir des données issues d'un site particulier, avec deux modèles 3D proposant des approches géométriques différentes, un modèle déterministe et un modèle stochastique, dont les propriétés géométriques et physiques ont été tirées des observations et travaux antérieurs effectués sur ce site de Rosemanowes. La simulation thermo-hydro-mécanique (THM) du modèle déterministe nous a permis de modéliser les échanges thermiques en régime transitoire dans le réservoir formé par le système de forages RH12/RH15 et d'estimer un ordre de grandeur des tractions d'origine thermique. Pour mieux étudier l'effet induit par la contraction des blocs de roche dans le temps nous utilisons le modèle 3D stochastique dont l'objectif principal est de simuler de façon plus réaliste la progression dans l'espace les ruptures en cisaillement. Avec ce modèle nous avons constaté l'apparition différée d'une activité et l'effet d'un cycle de pulses de pression, ce qui suggère un moyen d'atténuer les fortes magnitudes potentielles des ruptures en cisaillement dues au refroidissement. / The general framework of our research deals with the development of geothermal energy for electricity production using the heat stored in geological formations at depths ranging in 3 to 5 km, Generally the environment is poorly permeable and fractured. Our main objective is to study the phenomenon of induced micro-earthquakes in relation to the cooling of the rock. The work is based on field experiences including long duration tests, conducted on the Rosemanowes site (Cornwall, UK). For this, we proceeded to the development of a calculation tool, FRACAS, able to simulate this phenomenon by introducing a dual thermal approach to better simulate the cooling of the reservoir due to long term fluid injections, which might be responsible for new failure mechanisms due to the induced tractions. In this context, we introduced a new algorithm to describe shear in stabilities, a mechanism of "stick-slip" type with the consideration of static/dynamic friction coefficients. The possibility of inducing micro-seismicity is then studied using the in situ data base, with two 3D models offering different geometric approaches, a deterministic model and a stochastic model whose geometrical and physical properties were obtained from observations and previous work on this Rosemanowes site. The Thermo-Hydro-Mechanical (THM) simulation using the deterministic model has allowed us to calibrate the transient heat transfer in the reservoir formed by the drilling system RH12/RH15 and to give an estimate of tensile stress of thermal origin. To better study the effect induced by the contraction of the rock during time, we use the stochastic 3D model whose main objective is to simulate a more realistic spatial migration of shear ruptures. With this model we found a delayed onset of shear activity and discuss the effect of pressure step tests. The results suggest a way to mitigate the potential impact of shear ruptures due to cooling.

Durant les différentes phases du stockage profond des déchets radio-actifs exothermiques (excavation, exploitation) jusqu'à la fermeture définitive, la roche hôte sera soumise à des sollicitations mécaniques, hydriques et thermiques couplées. Afin de connaître et de modéliser le comportement à court et long terme des dispositifs de stockage, une investigation approfondie du comportement de la roche est nécessaire afin de compléter les données existantes. C'est dans ce but que cette étude du comportement thermo-hydro-mécanique de la formation argileuse du Callovo-Oxfordien (COx) considérée par l'ANDRA comme roche hôte potentielle, a été développée. Dans un premier temps le comportement en compression–gonflement de l'argilite du COx a été étudié par la réalisation d'un programme d'essais oedométriques haute pression. Les résultats, interprétés en termes de couplage endommagement-gonflement, ont montré que l'ampleur du gonflement était lié à la densité de fissuration engendrée lors de la compression. Dans un second temps, le comportement hydromécanique et thermo-hydro-mécanique de l'argilite saturée sous une contrainte moyenne proche de l'in-situ a été étudié à l'aide de cellules à faible chemin de drainage (10 mm), dont une cellule isotrope et une cellule triaxiale à cylindre creux avec mesures locales d'un type nouveau. Ces appareils ont permis de résoudre deux problèmes difficiles typiques des argilites de très faible perméabilité : i) une saturation préalable correcte, attestée par de bonnes valeurs du coefficient de Skempton et ii) de bonnes conditions de drainage. Les paramètres caractéristiques du comportement de l'argilite à température ambiante (coefficients de Skempton et de Biot, compressibilité drainée et non-drainée) ont été déterminés à partir d'essais de compression isotrope qui ont également confirmé l'isotropie transverse du matériau. La cohérence des paramètres obtenus a été vérifiée dans un cadre poro-élastique saturé. Deux aspects du comportement thermo-hydro-mécanique de l'argilite du COx ont été étudiés à l'aide d'essais de chauffage et de compression volumique en température (80°C) : les effets de la température sur le comportement intrinsèque de l'argilite et le phénomène de pressurisation thermique. Un essai de chauffage drainé sous contrainte in-situ a mis en évidence, apparemment pour la première fois, un comportement plastique contractant de l'argilite (comme les argiles normalement consolidées), non pris en compte dans les modélisations thermo-élastiques actuelles des systèmes de stockage profond. Un autre élément nouveau et important observé est l'augmentation de la compressibilité avec la température, à la différence des argiles. L'étude de la pressurisation thermique (engendrée par la faible perméabilité de l'argilite et la forte différence entre les coefficients de dilatation thermique de l'eau et de la phase solide), a été réalisée à l'aide d'essais de chauffage non drainés, suite à une analyse détaillée des effets perturbateurs du système de mesure lors de variations de pression et de température (peut-être à considérer également dans les mesures in-situ). Le coefficient de pressurisation thermique s'est révélé être très sensible aux variations de température et de contrainte, il diminue de 0,14 à 0,1 MPa/°C entre 25 et 80°C. La nature des différentes réponses hydro-mécaniques et thermo-hydro-mécaniques obtenues au cours de ce travail permettront une interprétation et une modélisation plus précises du comportement du massif d'argilite autour des galeries, dans des zones qui sont pour la plupart saturées, sauf à proximité immédiate des galeries (quelques décimètres) / During the different phases of the exothermic radioactive waste deep disposal (excavation, operation) and after permanent closure, the host rock is submitted to various coupled mechanical, hydraulic and thermal phenomena. Hence, a thorough investigation of the thermo-hydro-mechanical behaviour of the rock is necessary to complete existing data and to better understand and model the short and long term behaviour of the Callovo-Oxfordian (COx) clay formation in Bure (Meuse/Haute-Marne - M/HM), considered by ANDRA as a potential host rock in France.In this work, the compression – swelling behaviour of the COx claystone was first investigated by carrying out a series of high-pressure oedometric tests. The results, interpreted in terms of coupling between damage and swelling, showed that the magnitude of swelling was linked to the density of the fissures created during compression. In a second step, the hydromechanical and thermo-hydro-mechanical behaviour of the saturated claystone under a mean stress close to the in situ one were investigated by using two devices with short drainage path (10 mm), namely a isotropic cell and a newly designed hollow cylinder triaxial cell with local displacement measurements. These devices helped to solve two majors problems related to testing very low permeability materials : i) a satisfactory previous sample saturation (indicated by good Skempton values) and ii) satisfactory drainage conditions. Some typical constitutive parameters (Skempton and Biot's coefficients, drained and undrained compressibility coefficients) have been determined at ambient temperature through isotropic compression tests that also confirmed the transverse isotropy of the claystone. The consistency of the obtained parameters has been checked in a saturated poroelastic framework. Two aspects of the thermo-hydro-mechanical behaviour of the COx claystone have then been investigated through different heating tests and through drained and undrained isotropic compression tests at elevated temperature (80°C) : the effects of temperature on the behaviour of claystone and thermal pressurization. A drained heating test under in-situ stress conditions evidenced, probably for first time, a plastic contractant response of the claystone (like normally consolidated clays), a feature not considered in the presently conducted numerical modelling of deep disposal systems. Another new important observed feature is the increase in drained compressibility of the COx claystone with temperature, not observed in clays. The investigation of thermal pressurization (caused by the low claystone permeability and by the significant difference in thermal expansion between water and the solid phase) has been carried out by means of undrained heating tests, after a detailed analysis of the major effects of the measurement system (which should perhaps be also analyzed when performing in-situ measurements). The thermal pressurization coefficient appeared to be quite sensitive to changes in temperature and stress, it decreased between 0.14 and 0.1 MPa/°C between 25 and 80°C. It is believed that the different thermo-hydro-mechanical volumetric responses obtained here allow a better interpretation and modelling of the behaviour of the claystone formation around the galleries in areas that are mostly saturated, except close to the galleries (a few decimetres)