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)

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.

Ce mémoire présente l'étude expérimentale et la modélisation à l'échelle microscopique du comportement hydromécanique des argilites, roche hôte potentielle pour le stockage souterrain des déchets radioactifs. Le champ de déformation est mesuré par microscopie électronique à balayage environnementale et corrélation d'images numériques. En premier lieu, on étudie le cas du chargement hydrique pur. Le champ de déformation obtenu est très hétérogène, et montre une anisotropie. La non-linéarité de déformation pour HR élevée est le résultat combiné d'une fissuration et d'un gonflement non-linéaire de la phase argileuse dû à des mécanismes différents selon humidité relative (HR). On constate une déformation irréversible lors d'un cycle hydrique, ainsi qu'un réseau de microfissures localisées dans la phase argileuse ou aux interfaces grain-matrice. Ensuite, on étudie le cas du chargement combiné hydrique et mécanique dans le MEBE. Trois types de bandes de déformation apparaissent au cours du chargement mécanique : horizontales (compaction), verticales (fissuration), et inclinées (cisaillement). Les bandes de cisaillement apparaissent plus tôt à HR plus élevée. Finalement, le matériau sous chargement hydrique est modélisé comme un composite constitué par de inclusions non gonflantes au sein d'une matrice gonflante. On calcule d'abord le champ de contrainte interne dû aux interactions inclusion-matrice ainsi qu'au gradient d'humidité, et ensuite la déformation globale. On en dérive un modèle micromécanique du type du problème d'Eshelby. De plus, des modélisations 2D aux éléments finis sont effectuées.

Ce rapport rassemble des travaux portant sur l'endommagement et la rupture de roches à des échelles allant de celle de l'échantillon de laboratoire à la celle de la croûte terrestre en passant par celle des massifs rocheux (falaises et mines souterraines). On y associe l'analyse de données acquises au laboratoire et in situ , leur modélisation statistique et conceptuelle et la simulation numérique. Cette approche permet d'appréhender la complexité du comportement des roches selon différents aspects. L'accent est d'abord mis sur les propriétés statistiques et les invariances d'échelles des structures d'endommagement et des événements sismiques observés pendant le processus de déformation et de rupture. Un modèle numérique, basé sur une approche mésoscopique de l'endommagement, est présenté pour la simulation des processus de rupture à court et long terme. Il permet de reproduire de nombreuses observations expérimentales en leur donnant un cadre mécanique et en permettant ainsi de préciser leurs conditions d'apparition. Cette approche permet de considérer le processus de déformation et de rupture des roches comme un système complexe dont les propriétés macroscopiques émergent de l'interaction entre éléments de plus petite échelle ayant des propriétés simples. Le rôle des fluides dans ce processus est abordé par l'observation de la sismicité induite dans un massif rocheux fracturé et la simulation numérique du couplage roche fluides. Enfin un projet de recherche est développé sur le thème de l'endommagement à long terme et du couplage hydromécanique.

Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique.

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

La perméabilité des roches est fortement influencée par la présence de fractures car ces dernières constituent un chemin préférentiel pour l’écoulement des fluides. Ainsi la présence de fractures naturelles est un facteur déterminant pour la productivité d’un réservoir. Dans le cas de roches à faible conductivité, des techniques de stimulation telle que la fracturation hydraulique sont utilisées pour en augmenter la perméabilité et rendre le réservoir exploitable d’un point de vue économique. A l’inverse, dans le cas du stockage géologique, la présence de fractures dans la roche représente un danger dans la mesure où elle facilite le transport et la migration des espèces disséminées dans la roche. Pour le stockage de CO2, les fuites par les fractures présentes dans les couvertures du réservoir et la réactivation des failles constituent un risque majeur. Et en ce qui concerne le stockage géologique de déchets radioactifs, la circulation de fluide dans des réseaux de fractures nouvelles ou réactivées au voisinage de la zone de stockage peut aboutir à la migration de matériaux nocifs. Il est donc important de prévoir les effets de la présence de fractures dans un réservoir. Le but de cette thèse est le développement d’un outil numérique pour la simulation d’un réseau de fractures et de son évolution sous sollicitation hydro-mécanique. Grâce à sa commodité, la méthode des éléments finis étendue (XFEM) sera retenue et associée à un modèle de zone cohésive. La méthode XFEM permet en effet l’introduction de fissures dans le modèle sans nécessairement remailler en cas de propagation des fissures. L’écoulement du fluide dans la fissure et les échanges de fluide entre la fissure et le milieu poreux seront pris en compte via un couplage hydro-mécanique. Le modèle est validé avec une solution analytique asymptotique pour la propagation d’une fracture hydraulique plane dans un milieu poroélastique en 2D comme en 3D. Puis, nous étudions la propagation de fractures hydrauliques sur trajets inconnus. Les fissures sont initialement introduites comme des surfaces de fissuration potentielles étendues. Le modèle de zone cohésive sépare naturellement les domaines adhérents et ouverts. Les surfaces potentielles de fissuration sont alors actualisées de manière implicite par un post-traitement de l’état cohésif. Divers exemples de réorientation de fissures hydrauliques et de compétition entre fissures voisines sont analysés. Enfin, nous présentons l’extension du modèle aux jonctions de fractures hydrauliques / The permeability of rocks is widely affected by the presence of fractures as it establishes prevailing paths for the fluid flow. Natural cracks are then a critical factor for a reservoir productiveness. For low permeability rocks, stimulation techniques such as hydrofracturing have been experienced to enhance the permeability, so that the reservoir becomes profitable. In the opposite, when it comes to geological storage, the presence of cracks constitutes a major issue since it encourages the leak and migration of the material spread in the rock. In the case of CO2 storage, the scenario of leakage across the reservoir seal through cracks or revived faults is a matter of great concern. And as for nuclear waste storage, the fluid circulation in a fracture network around the storage cavity can obviously lead to the migration of toxic materials. It is then crucial to predict the effects of the presence of cracks in a reservoir. The main purpose of this work is the design of a numerical tool to simulate a crack network and its evolution under hydromechanical loading. To achieve this goal we chose the eXtended Finite Element Method (XFEM) for its convenience, and a cohesive zone model to handle the crack tip area. The XFEM is a meshfree method that allows us to introduce cracks in the model without necessarily remeshing in case of crack propagation. The fluid flow in the crack as well as the exchanges between the porous rock and the crack are accounted for through an hydro-mechanical coupling. The model is validated with an analytical asymptotic solution for the propagation of a plane hydraulic fracture in a poroelastic media, in 2D as well as in 3D. Then we study the propagation of hydraulic fractures on non predefined paths. The cracks are initially introduced as large potential crack surfaces so that the cohesive law will naturally separate adherent and debonding zones. The potential crack surfaces are then updated based on a directional criterion appealing to cohesive integrals only. Several examples of crack reorientation and competition between nearby cracks are presented. Finally, we extend our model to account for the presence of fracture junctions

Le présent travail montre la modélisation constitutive d'un géomatériau composé d'une matrice poreuse saturée et déformable contenant une distribution périodique de fissures évolutives remplies de fluide. La méthode d'homogénéisation des développements asymptotiques est utilisée afin de déduire un modèle capable de décrire le couplage hydro-mécanique macroscopique. Prenant en considération l'évolution de fissures et sans faire des hypothèses phénoménologiques, une analyse énergétique mésoscopique couplé avec un schéma d'homogénéisation a été développée et elle fournit une loi d'évolution d'endommagement macroscopique. De cette façon, un lien direct entre les phénomènes de rupture de la structure mésoscopique et l'endommagement macroscopique correspondant est établie. Finalement, on présente une étude numérique du comportement macroscopique d'endommagement hydro-mécanique. / This work presents the constitutive modeling of a geomaterial consisting of a deformableand saturated porous matrix including a periodic distribution of evolving fluid-filledcavities. The homogenization method based on two-scale asymptotic developments isused in order to deduce a model able to describe the macroscopic hydro-mechanicalcoupling. By taking into account the cavity growth and without any phenomenologicalassumption, it is proposed a mesoscopic energy analysis coupled with the homogenizationscheme which provides a damage evolution law. In this way, a direct link between themeso-structural fracture phenomena and the corresponding macroscopic damage isestablished. Lastly, a numerical study of the local macroscopic hydro-mechanical damage behaviour is presented. / In questa tesi si presenta la modellazione costitutiva di un geomateriale composto da unamatrice porosa satura e deformabile contenente una distribuzione periodica di cavitàriempite da fluido che si propagano. Il metodo di omogeneizzazione basato sugli sviluppiasintotici a doppia scala viene utilizzato con l'obiettivo di dedurre un modello capace didescrivere l'accoppiamento idro-meccanico macroscopico. Prendendo in considerazione lapropagazione delle cavità e senza nessuna ipotesi fenomenologica, si propone un'analisienergetica mesoscopica accoppiata ad uno schema di omogeneizzazione che fornisce unalegge di evoluzione del danno.In questo modo, una relazione diretta tra i fenomeni difrattura meso-strutturali ed il corrispondente danno macroscopico viene stabilita. Infine,uno studio numerico del comportamento macroscopico locale di danno idro-meccanico viene presentato.

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