Le but de cette thèse est de mettre en place un nouveau modèle de comportement hydromécanique permettant de prendre en compte les anisotropies initiale et induite et l'impact de l'endommagement mécanique sur la perméabilité. Afin de construire ce modèle, une nouvelle approche de modélisation permettant d'étendre les lois de comportement mécaniques des matériaux isotropes aux matériaux anisotropes est développée. Cette approche, employée dans le cadre des milieux continus à variables internes, est utilisée pour construire une loi de comportement elasto-viscoplastique qui distingue les régimes de sollicitation en compression et en traction. Un tenseur de second ordre est introduit pour décrire l'anisotropie induite suite à des sollicitations de traction et une variable interne scalaire est utilisée pour traduire le durcissement/adoucissement du matériau suite à des sollicitations de compression. Sous des sollicitations complexes, ces deux mécanismes sont couplés et l'effet de fermeture/réouverture des fissures est traité. Le couplage endommagement-perméabilité est ensuite modélisé par l'introduction d'une loi phénoménologique reliant la perméabilité intrinsèque du matériau aux variables internes de la mécanique.Ce modèle a été appliqué dans le cas des ouvrages souterrains du site de Bure afin de comprendre les mécanismes d'altération des propriétés hydromécaniques autour des galeries et des alvéoles de stockage causée non seulement par les opérations de creusement mais également par les surpressions dues à la production d'hydrogène gazeux suite à la corrosion des parties métalliques des modules de déchets. / This thesis aims to introduce a new hydromechanical constitutive model taking into account both initial and induced anisotropies and the impact of the mechanical damage on the permeability. To build this model, a new modeling approach is developed allowing the extension of mechanical behavior laws from isotropic materials to transversely isotropic materials. This approach is used, within the framework of continuous media with internal variables, to propose an elasto-viscoplastic behavior law that distinguishes between compressive and tensile loading regimes. A second order tensor is introduced to describe the induced anisotropy due to tensile loadings, and a scalar internal variable is employed to account for hardening and softeningof the material due to compressive loadings. Under complex loadings, these two mechanisms are coupled, and the effect of cracks closing/reopening is taken into consideration. The damage-permeability coupling is modeled by the introduction of a phenomenological law linking the material intrinsic permeability to the mechanical internal variables.The developed model is applied to the case of the underground drifts of Bure site in order to better understand the mechanisms of hydromechanical properties alteration, around drifts and storing cells. Not only the impact of the excavation operations is considered but also the consequences of the overpressures caused by the produced hydrogen due to the corrosion of the metallic parts of nuclear waste containers.

Ce travail vise à comprendre et quantifier le comportement de U(VI) dans les argilites du Callovo-Oxfordien(COₓ), prévues pour accueillir les déchets nucléaires de haute activité en France. L’effet de la température sur ce comportement est particulièrement étudié. La première partie du travail s'est concentrée sur les propriétés thermodynamiques des complexes ternaires Ca-U(VI)-CO₃ qui contrôlent la spéciation de U(VI) en solution. Ces dernières ont été mesurées par une méthode de compétition en présence d’une résine sous pCO₂ et pH contrôlés. Les résultats indiquent que la température ne favorise pas la formation de CaUO₂(CO₃)₃²⁻ (log₁₀ βº₁₁₃ = 27,3 ± 0,3 ; ΔrHº = -27,4 ± 8 kJ/mol) et n'affecte pas la formation de Ca₂UO₂(CO₃)₃(aq)(log₁₀ βº₂₁₃ = 29,7 ± 0,3 ; ΔrHº = 0 ± 2 kJ/mol). Une approche « bottom-up » avec le modèle « 2SPNE SC /CE » publié dans la littérature a été suivie pour décrire les phénomènes de sorption, en considérant que la fraction argileuse du COₓ (Illite et I/S) gouverne la sorption de U(VI). Ce modèle a été appliqué avec succès pour reproduire une multitude de données expérimentales obtenues avec l'illite, la fraction argileuse du COₓ et les argilites du COₓ en fonction de paramètres clés (pH, pCO₂, [U(VI)], [Ca]) à 20 °C. Le phénomène de rétention dans les conditions in-situ est principalement régi par la sorption des complexes U(VI)-carbonate et une nouvelle constante d´équilibre de réaction de complexation de surface est proposée pour l'illite. Une augmentation de la température à 80 °C conduit à une augmentation de la rétention de U(VI) sur le COₓ. Cette augmentation est accompagnée d'un changement de pCO₂et de pH. Le modèle de rétention testé à 20 °C combiné avec les paramètres thermodynamiques décrivant le comportement de U(VI) en solution expliquent cette augmentation sans pour autant obtenir un accord satisfaisant avec l’expérience. Le modèle est amélioré en intégrant des valeurs de ΔrHº obtenues pour les réactions de complexation de surface à partir du système U(VI)/illite. / The aim of this study was to understand and quantify the behaviour of U(VI) on the Callovo-Oxfordian(COx) clay which is envisioned to host high-level radioactive waste in France. The temperature effect up to 80°C on this behaviour was particularly studied. The first part of the work focussed on the thermodynamic properties of the calcium uranyl carbonate aqueous complexes which govern U(VI) speciation in solution. They were measured indirectly by sorption-based methodologies under controlled pCO₂ and pH. The results indicate that the temperature does not favour the formation of CaUO₂(CO₃)₃²⁻ (log₁₀ βº₁₁₃ = 27,3 ± 0,3 ; ΔrHº = -27,4 ± 8 kJ/mol) and does not affect the formation of Ca₂UO₂(CO₃)₃(aq)(log₁₀ βº₂₁₃ = 29,7 ± 0,3 ; ΔrHº = 0 ± 2 kJ/mol). A bottom-up approach with the published “2SPNE SC/CE”model was used for describing the sorption processes, with the assumption that the clay fraction of the COx (Illite, andI/S) governs U(VI) sorption.The model was successfully applied to reproduce a wealth of experimental data obtained with illite, the COₓ clay fractionand the COₓ clay rock as a function of key parameters (pH, pCO2, [U(VI)], [Ca]) at 20 °C. The sorption on COₓ conditions is mainly governed by the sorption of U(VI)-CO3 complexes and a new sorption constant is proposed for illite. An increase in temperature to 80 °C leads to an in-crease of U(VI) retention on COx. This increase is ac-companiedby a change of both pCO₂ and pH. The sorption model developed at 20 °C, together with the thermodynamic parameters describing U(VI) speciation in solution, can explain this increase but without obtaining a good agreement with the experiment. The model is improved by considering ΔrHº values for sur-face complexation reactions obtained for the U(VI))/illite system.

L'Agence nationale pour la gestion des déchets radioactifs (Andra) a commencé en 2000 la construction du Laboratoire Souterrain de Meuse / Haute-Marne (LS-M/HM) avec l'objectif principal de démontrer la faisabilité d’un stockage géologique dans l’argilite du Callovo-Oxfordien. Un réseau de galeries expérimentales a été excavé, principalement en suivant les directions des contraintes horizontales (majeure et mineure), avec des variations sur : la méthode d'excavation, la géométrie de la structure et le soutènement. Chaque galerie a été instrumentée en différentes sections pour suivre le comportement hydromécanique de la roche face à l’excavation. Le suivi de la zone autour des galeries excavées au niveau principal (-490 m) a révélé le développement d'une zone fracturée (fractures en extension et en cisaillement) induite par l'excavation. La distribution de la zone fracturée dépend à la fois de l'orientation de la galerie et du champ de contraintes in-situ et a une influence importante sur la déformation des galeries. En effet, les mesures de convergence ont montré une fermeture anisotrope de la section de la galerie. De plus, il a été observé un champ de distribution anisotrope de la pression de pores ainsi que des surpressions autour des galeries.Afin d’analyser la réponse anisotrope du massif pendant l’excavation et après celle-ci, les travaux effectués dans le cadre de la thèse sont axés principalement sur une étude directe des mesures de convergence in-situ. Cette analyse s’effectue à l’aide de la loi semi-empirique proposée par Sulem et al. (1987) [Int J Rock Mech Min Sci Geomech Abstr 24: 145–154]. A cet égard, différentes galeries excavées dans le LS-M/HM ont été étudiées. Ces galeries présentent certaines différences dans leurs orientations et l’état initial des contraintes, dans la méthode et la vitesse d’excavation ainsi que dans les diamètres de la section et les types de soutènements installées. Cette analyse permet d’obtenir des prédictions fiables de la convergence à long-terme, ce qui peut servir pour le dimensionnement et la prévision de la performance du soutènement à long-terme.En outre, nous avons étudié la réponse anisotrope du champ de pression interstitielle observée in-situ. Cette analyse est basée sur une approche poroélastique anisotrope. L’objectif principal est de reproduire qualitativement l’évolution de la pression des pores autour des galeries avec une approche simple qui prend en compte l’anisotropie intrinsèque du matériau. Enfin, une analyse de l’apparition de la rupture montre le rôle clé que joue le couplage hydromécanique dans l’extension de la zone fracturée / The French National Radioactive Waste Management Agency (Andra) began in 2000 the construction of an Underground Research Laboratory (URL) with the main goal of demonstrating the feasibility of a geological repository in Callovo-Oxfordian claystone. Several research programs have taken place to improve the knowledge of the rock properties and its response to the excavation progress. A network of experimental drifts has been constructed with variations on: excavation method, structure geometry, supports system and orientations with respect to principal stresses’ directions. In each drift different sections have been instrumented to monitor the hydro-mechanical behavior of the rock mass formation. Continuous monitoring of the excavated zone around the drifts in the main level (-490 m) revealed the development of a fractured zone (extensional and shear fractures) induced by the excavation. The extent of this fractured zone depends on the drift orientation regarding the in-situ stress field. Accordingly, the convergence measurements showed an anisotropic closure which depends also on the drifts’ orientations. Moreover, marked overpressures and an anisotropic pore pressure field around the drifts have been also observed.The approach proposed in this work is mainly based on a direct analysis of the convergence measurements, for studying the anisotropic response of the rock formation during and after excavation. The convergence evolution is analyzed on the basis of the semi-empirical law proposed by Sulem et al. (1987) [Int J Rock Mech Min Sci Geomech Abstr 24: 145–154]. The monitoring and analysis of convergence data can provide a reliable approach of the interaction between rock mass and support. Therefore, the anisotropy and the variability of the closure are analyzed taking into account different field cases: drifts excavated in two different orientations (i.e. influence of the initial stress state), different methods, sizes and rates of excavation and different supports systems with different conditions of installation. This broad range of cases permits to refine the analysis for reliable predictions of the convergence evolution in the long term. This approach can thus be used for the design of various types of support and the evaluation of its performance in the long term.On the other hand, the pore pressure evolution induced by excavation of drifts as recorded in situ has been analyzed. The anisotropic response observed in-situ suggests that the intrinsic anisotropy of the material plays a key role in the response of the rock formation. To understand these phenomena, an anisotropic poroelastic analysis of the pore pressure evolution induced by the drift excavation is performed. The main goal is to simulate the main trends of the pore pressure evolution with a simple model taking into account the inherent anisotropy of the material. Finally, an analysis of the onset of failure shows the key role of the hydro-mechanical coupling on the extension of the failed zone around the drifts