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)

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 / The investigation of the behavior of heated concrete is a major research topic which concerns the assessment of safety level of structures when exposed to high temperatures, for instance during a fire. For this purpose, several modeling approaches were developed within thermo-hydro-mechanical (THM) frameworks in order to take into account the involved physic-chemical and mechanical processes that affect stability of heated concrete. However, existing models often do note account explicitly for the heterogeneity of the material : concrete is composite material that may be schematized as an assembly of inclusions (aggregates) embedded in a cementitious matrix (cement paste). This latter may be described as a partially saturated open porous medium. The aggregates are characterized by their mineralogical nature together with their morphology and size distribution. The material heterogeneity bring an additional complexity : the need to take into account the microstructure in order to quantify the effect of matrix-inclusion thermal, hygral and mechanical incompatibilities on the THM behavior of concrete. This work is a first step in this direction. For this purpose, a three-dimensional (3D) multi-scale finite element model is developed. It allows affecting specific behaviors to matrix and inclusions. For the former, where mass transports occur within the connected porous network, a three-fluids approach (liquid water, vapor and dry air) is adopted and is coupled to a poro-mechanical damage based approach. For inclusions (aggregates) no hygral component arises a pure thermo-mechanical model is considered. The developed model is then used to investigate, either by 2D or 3D numerical simulations, effects of mineralogical nature, morphology and distribution of aggregates. Studied effects have mainly concerned the influence of these parameters on local fluctuations of simulated temperature, gas pressure and damage fields with regard to experimentally observed dispersion. The analysis is here limited to the mesoscale, at which only aggregates with a characteristic size above 5 mm are meshed while the remaining inclusions together with the cement paste are considered to be a homogeneous matrix. Finally, the numerical analysis carries out the need to perform an experimental campaign that is consistent with a multi-scale approach of the THM behavior of concrete : an experimental campaign that allows to identify thermal, hygral and mechanical properties of each concrete constituent and that permit to assess evolution of fields during upscalling. An experimental protocol is then elaborated for this purpose and some obtained results are presented and analyzed with regard to results obtained in the modeling part