Gelfill (GF) is made of tailings, water, binder and chemical additives (Fillset, sodium silicate gel). The components of GF are combined and mixed on the surface and transported (by gravity and/or pumping) to the underground mine workings, where the GF can be used for both underground mine support and tailings storage. Thermal (T), hydraulic (H), and mechanical (M) properties are important performance criteria of GF. The understanding of these engineering properties and their evolution with time are still limited due to the fact that GF is a new cemented backfill material. In this thesis, the evolution of the thermal, hydraulic, mechanical, and microstructural properties of small GF samples are determined. Various binder contents of Portland cement type I (PCI) are used. The GF is cured for 3, 7, 28, 90, and 120 days. It is found that the thermal, hydraulic and mechanical properties are time-dependent or affected by the degree of binder hydration index. Furthermore, a relationship is found between the compressive strength and the saturated hydraulic conductivity of the GF samples. The unsaturated hydraulic properties of GF samples have also been investigated. The outcomes show that unsaturated hydraulic conductivity is influenced by the degree of binder hydration index and binder content, especially at low suction ranges. Simple functions are proposed to predict the evolution of air-entry values (AEVs), residual water content, and fitting parameters from the van Genuchten model with the degree of hydration index (α). Furthermore, two columns are built to simulate the coupled thermo-hydro-mechanical (THM) behaviour of GF under drained and undrained conditions. The obtained results from the GF columns are compared with the small samples. It is observed that the mechanical properties, hydraulic properties (suction and water content), and temperature development are strongly coupled. The magnitude of these THM coupling factors is affected by the size of the GF. The findings also show that the mechanical, hydraulic and thermal properties of the GF columns are different from samples cured in plastic moulds.

Gelfill (GF) is made of tailings, water, binder and chemical additives (Fillset, sodium silicate gel). The components of GF are combined and mixed on the surface and transported (by gravity and/or pumping) to the underground mine workings, where the GF can be used for both underground mine support and tailings storage. Thermal (T), hydraulic (H), and mechanical (M) properties are important performance criteria of GF. The understanding of these engineering properties and their evolution with time are still limited due to the fact that GF is a new cemented backfill material. In this thesis, the evolution of the thermal, hydraulic, mechanical, and microstructural properties of small GF samples are determined. Various binder contents of Portland cement type I (PCI) are used. The GF is cured for 3, 7, 28, 90, and 120 days. It is found that the thermal, hydraulic and mechanical properties are time-dependent or affected by the degree of binder hydration index. Furthermore, a relationship is found between the compressive strength and the saturated hydraulic conductivity of the GF samples. The unsaturated hydraulic properties of GF samples have also been investigated. The outcomes show that unsaturated hydraulic conductivity is influenced by the degree of binder hydration index and binder content, especially at low suction ranges. Simple functions are proposed to predict the evolution of air-entry values (AEVs), residual water content, and fitting parameters from the van Genuchten model with the degree of hydration index (α). Furthermore, two columns are built to simulate the coupled thermo-hydro-mechanical (THM) behaviour of GF under drained and undrained conditions. The obtained results from the GF columns are compared with the small samples. It is observed that the mechanical properties, hydraulic properties (suction and water content), and temperature development are strongly coupled. The magnitude of these THM coupling factors is affected by the size of the GF. The findings also show that the mechanical, hydraulic and thermal properties of the GF columns are different from samples cured in plastic moulds.

Gelfill (GF) is made of tailings, water, binder and chemical additives (Fillset, sodium silicate gel). The components of GF are combined and mixed on the surface and transported (by gravity and/or pumping) to the underground mine workings, where the GF can be used for both underground mine support and tailings storage. Thermal (T), hydraulic (H), and mechanical (M) properties are important performance criteria of GF. The understanding of these engineering properties and their evolution with time are still limited due to the fact that GF is a new cemented backfill material. In this thesis, the evolution of the thermal, hydraulic, mechanical, and microstructural properties of small GF samples are determined. Various binder contents of Portland cement type I (PCI) are used. The GF is cured for 3, 7, 28, 90, and 120 days. It is found that the thermal, hydraulic and mechanical properties are time-dependent or affected by the degree of binder hydration index. Furthermore, a relationship is found between the compressive strength and the saturated hydraulic conductivity of the GF samples. The unsaturated hydraulic properties of GF samples have also been investigated. The outcomes show that unsaturated hydraulic conductivity is influenced by the degree of binder hydration index and binder content, especially at low suction ranges. Simple functions are proposed to predict the evolution of air-entry values (AEVs), residual water content, and fitting parameters from the van Genuchten model with the degree of hydration index (α). Furthermore, two columns are built to simulate the coupled thermo-hydro-mechanical (THM) behaviour of GF under drained and undrained conditions. The obtained results from the GF columns are compared with the small samples. It is observed that the mechanical properties, hydraulic properties (suction and water content), and temperature development are strongly coupled. The magnitude of these THM coupling factors is affected by the size of the GF. The findings also show that the mechanical, hydraulic and thermal properties of the GF columns are different from samples cured in plastic moulds.

Les tunnels profonds sont souvent construits dans les roches sédimentaires et métamorphiques stratifiées qui présentent habituellement des propriétés anisotropes en raison de leur structure et des propriétés des constituants. Le présent travail vise à étudier les tunnels profonds dans un massif rocheux anisotrope élastique en portant une attention particulière sur les effets des couplages hydromécaniques par des approches analytiques et numériques. Une solution analytique pour un tunnel creusé dans un massif rocheux anisotrope saturé est développée en tenant compte du couplage hydro-mécanique dans le régime permanent. Cette solution analytique est utilisée pour réaliser une série d’études paramétriques afin d'évaluer les effets des différents paramètres du matériau anisotrope sur le comportement du tunnel. Dans un deuxième temps la solution analytique est élargie pour décrire le comportement du tunnel pendant la phase transitoire hydraulique. Afin de compléter ces études analytiques qui prennent en compte seulement un couplage unilatéral (dans le sens que seul le comportement hydraulique influence le comportement mécanique et pas l’inverse) de l’analyse numérique avec un couplage complet, ont été réalisés. Une application de la solution analytique sur la méthode de convergence-confinement est aussi bien abordée qui peut prendre en compte l'influence du front de taille du tunnel sur le travail du soutènement ainsi que sur le massif. La solution obtenue peut servir comme un outil de dimensionnement rapide des tunnels en milieux poreux en le combinant avec des approches de dimensionnement comme celle de convergence confinement. / Deep tunnels are often built in the sedimentary and metamorphic foliated rocks which exhibits usually the anisotropic properties due to the presence of the discontinuity. The analysis of rock and liner stresses due to tunnel construction with the assumption of homogeneous and isotropic rock would be incorrect. The present thesis aims to deal with the deep tunnel in anisotropic rock with a particular emphasis on the effects of hydraulic phenomenon on the mechanical responses or reciprocal effects of hydraulic and mechanical phenomena by combining analytical and numerical approach. On that point of view, a closed-formed solution for tunnel excavated in saturated anisotropic ground is developed taking into account the hydromechanical coupling in steady-state. Based on the analytical solution obtained, parametric studies are conducted to evaluate the effects of different parameters of the anisotropic material on the tunnel behavior. The thesis considers also to extend the analytical solution with a time-dependent behavior which takes into account the impact of the pore pressure distribution on mechanical response over time, i.e., one way coupling modeling. In addition, some numerical analysis based on fully-coupled modeling, i.e., two ways coupling, are conducted which are considered as the complete solution for the analytical solution. An application of the closed-form solution on convergence-confinement method is as well referred which can take into account the influence of the tunnel face on the work of the support as well as the massif. The obtained solution could be used as a quick tool to calibrate tunnels in porous media by combining with design approaches such as convergence-confinement method.

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. / We present two different approaches to describe the hydromechanical behaviour of geomaterials. In the first approach the porous media is studied through an equivalent continuum media where the interaction between the fluide and solid phases caracterize the coupling behaviour at the macroscale.We take into account this approach to model experimental tests performed over a hollow cylinder sample of clay rock (Boom Clay), considered for nuclear waste storage. The experimental results clearly show that the mechanical behaviour of the material is strongly anisotropic. For this reason we chose an elasto-plastic model based on Drucker-Prager criterion where the elastic part is characterized by cross anisotropy.The numerical results of boundary value problem clearly show localised strains around the inner hollow section. In order to regularize the numerical problem we consider a second gradient local continuum media with an enriched kinematic where an internal lenght can be introduced making the results mesh independent. The uniqueness study is carried out showing that changing the temporal discretization of the problem leads to different solutions.In the second approach we study the hydromechanical behaviour of a porous media that it is characterised by the microstructure of the material. The microstructure taken into account is composed by elastic grains, cohesives interfaces and a network of fluid channels. Using a periodic media a numerical homogenization (square finite element method) is considered to compute mass flux, stress and density of the mixture. In this way a pure numerical constitutive law is built from the microstructure of the media. This method has been implemented into a finite element code (Lagamine, Université de Liège) to obtain results at the macroscale. A validation of this implementation is performed for a pure mechanical boundary value problem and for a hydromechanical one.

Este trabalho caracterizou os siltitos do Grupo Itararé, rocha branda proveniente da amostragem de pedreira localizada nas proximidades de futuro empreendimento brasileiro, o Trem de Alta Velocidade. Objetivou entender as causas da pronunciada e rápida alteração da rocha intacta, e simular, em condições controladas de laboratório, a degradação hidromecânica causada pela simples variação da umidade relativa do ar. Para tal, estabeleceu-se técnica de amostragem de testemunhos por meio da extração via seca com ar-comprimido, e processo de intemperismo artificial, composto por 4 ciclos de 2 etapas (secagem e umedecimento). O processo se deu em dessecadores, sendo a variação da umidade ocasionada pela técnica do vapor e sílica gel. O acompanhamento da degradação dos corpos de prova efetuou-se com duas técnicas não-destrutivas de propagação de ondas: o pulso ultrassônico e a excitação por impulso (sistema Sonelastic). Em paralelo, realizou-se a caracterização geotécnica dos siltitos em umidade natural, com ensaios mineralógicos, físicos, físico-mecânicos, físico-químicos e de alteração laboratorial. Por meio desta metodologia, foi possível estimar as minorações do módulo de deformabilidade e da resistência, além das causas físico-químicas da degradação. / This work characterized the siltstones of the Itararé Group, a soft rock sampled in a quarry located near a future Brazilian venture, the High Speed Train. It aims at understanding the causes of the pronounced and fast changes of the intact rock, and simulate in controlled laboratory conditions the hydro degradation process caused by simple variation of relative air humidity. Dry sampling technique with compressed air was used to extract samples from blocks, and the process of artificial weathering consisted on 4 cycles of 2 steps (drying and wetting). The process occurred in desiccators, where moisture variation was caused by the vapor technique and silica gel. The monitoring of the specimens degradation was performed with two non-destructive wave propagation techniques: the ultrasonic pulse and the impulse excitation (Sonelastic system). A geotechnical characterization of the siltstones in natural moisture condition was made in parallel, consisting on mineralogical, physical, physical-mechanical, physical-chemical and durability tests. By this methodology, it was possible to estimate the mitigation of the deformability modulus and strength besides the physical- chemical causes of degradation.

The deep geological repository concept is based on “engineered barriers systems (EBS)” that are constructed in the repository and “natural barriers” provided by the surrounding geological environment. The EBS comprises of variety of sub-systems or components, such as the waste form, canister, buffer, backfill, seals, and plugs. Geological disposal is based on the concept of multiple barriers that work together to provide containment. The buffer is made up of densely compacted bentonite or bentonite-sand mix. Bentonite has both mechanical and physico-chemical functions, to fulfill as a barrier material in DGR. The bentonite buffer should hold the containers in place and prevent collapse of the excavation. A plastic deformability of the bentonite is desired to redistribute the stresses that can result from creep in the rock, and prevent transfer of excessive stresses to the canisters. The bentonite buffer must create an impermeable zone around the containers to ensure that the radionuclide released from the vitrified waste is limited by diffusive transport rather than advective transport in groundwater. Another important property of the highly compacted bentonite is its swelling potential. Its swelling potential should be as high as possible, to guarantee the sealing of any cracks occurring in the buffer material or in the storage gallery and thus ensure good imperviousness. Besides its mechanical function, bentonite buffer must sorb escaping radionuclides and thus retard their migration to the geo-environment. The bentonite buffer must retain its mechanical and physico-chemical functions over a span of several hundred thousand years to fulfill its role as a containment barrier in DGR. The bentonite buffer should maintain its physico-chemical and hydro-mechanical integrity on exposure to groundwater. Nuclear power agencies of several countries have identified suitable bentonites for use as buffer in DGR through laboratory experiments and large scale underground testing facilities. Japan has identified Kunigel VI bentonite, South Korea-Kyungju bentonite, China-GMZ bentonite, Belgium-FoCa clay, Sweden-MX-80 bentonite, Spain-FEBEX bentonite and Canada-Avonseal bentonite as candidate bentonite buffer for deep geological repository program. Bentonite from Barmer (Rajasthan State) was identified as suitable buffer for use in Indian deep geological repositories. The influence of moisture and dissolved salt migration on the physico-chemical and hydro-mechanical properties of Barmer bentonite has not been examined. The study is important to understand the clay’s behaviour under deep geological repository conditions, where, the bentonite buffer would come in contact with groundwater. Infiltration of groundwater with variable chemical composition could alter the physico-chemical and hydro-mechanical properties of the clay. The objectives of the thesis are as follows: Examine the influence of permeation of distilled water (DW) and synthetic ground water (SGW) solutions under constant volume condition on suction, physico-chemical and moisture content/dry density characteristics of compacted Barmer bentonite specimens as function of permeation period (maximum permeation period– 30 days). Examine the influence of variation in dry density and gravimetric water content as consequence of DW and SGW solution permeation on swell pressure and unconfined compression strength of Barmer bentonite specimens. Compare experimental swell pressures of re-constituted bentonite specimens with swell pressures predicted by diffuse double layer models. Examine the influence of total dissolved solids (TDS) concentration of permeating solution on the unsaturated permeability of compacted Barmer bentonite specimens. Organization of thesis: After the first introductory chapter, a detailed review of literature is performed in Chapter 2 to review the physicochemical, mineralogical and hydro-mechanical properties of bentonites identified as buffer materials for deep geological repositories of various countries. Based on current understanding and need to perform similar studies with Barmer clay, the chapter develops the scope and objectives of the study. Chapter 3 presents a detailed experimental program of the study. Chapter 4 examines the influence of permeation of distilled water (DW) and synthetic groundwater (SGW) solutions (under constant volume conditions) on the total suction of compacted bentonite specimens at two locations in the clay. The influence of variation in dry density on the moisture migration-suction inter-relations of compacted bentonite specimens is also examined. The associated changes of DW and SGW solution migration under constant volume conditions on the physico-chemical properties, water content and dry density of compacted Barmer bentonite specimens are also examined. The experimental results brought out that matric suction mainly contributed (75 to 92 %) to total suction of the permeated specimens; the permeated specimens experienced reduction in matric suction with increase in gravimetric water content from increase in degree of saturation. Osmotic suction contributed to 10 to 25 % of the total suction of the permeated specimens and was observed to increase with gravimetric water content due to solubilization of salts contained in the voids of the compacted bentonite specimens. The total suction of compacted Barmer bentonite specimen was responsive to the total dissolved solids concentration of the permeating solutions as the specimen permeated with more saline solution (higher TDS value) exhibited lesser total suction. Upon permeation with DW and SGW solutions, the CEC of bentonite was unaltered, while, pH and TDS values were affected. Softening of the bentonite clay occurred from increase in water content and existence of compression zones (material used to seal 1mm gap in relative humidity probe aperture) that in turn facilitated dissipation of swelling stress leading to reduction in dry density values. Chapter 5 examines influence of reduction in dry density and increase in water content on the swell pressure and compression strength characteristics of compacted Barmer bentonite specimens upon DW and SGW solution migration as the results could provide insight into possible deviations from the design properties upon wetting of bentonite buffer by groundwater under deep geological repository conditions. The experimental swelling pressures are also compared with those predicted by Gouy-Chapman diffuse double layer theory. The dry density of 1.6 Mg/m specimens permeated with DW and SGW solutions reduced to 1.59 to 1.36 Mg/m and water contents increased to 18.9 to 27 % on permeation with distilled water and SGW solutions for 30 days. The reductions in dry density and increase in water content caused 30 to 70 % reductions in swell pressures and 31 to 74 % decrease in unconfined compression strength values. Specimens initially compacted to dry density of 1.8 Mg/m, experienced reduction in dry density ranging from 1.79 to 1.52 Mg/m and increase in water content from 18.6 to 24.2 % on permeation of DW and SGW solutions for 30 days. These reductions in dry density and increase in water caused the swell pressures to reduce from 4 to 55 % and unconfined compressive strengths to reduce by 31 to 67 %. Comparison of swell pressures gave -8 to 127 % variations between theoretical (from DDL theory) and experimental values due to errors associated with estimation of surface area and dissolved salt concentrations in pore water. Chapter 6 examines the influence of salinity of permeating solution on the unsaturated permeability of compacted Barmer bentonite specimens. The salinity of permeants was varied by permeating distilled water (DW) and synthetic ground water solutions under constant volume conditions over maximum period of 30 days. Experimental results showed that the saturated permeability coefficients (ksat) of specimens compacted to 1.6 Mg/m, responded to variations in TDS of the permeant. Comparatively, the ksat values of specimens compacted to 1.8 Mg/mwere unaffected by variation in TDS of the permeant. Permeation of DW and SGW solutions decreased the ksat values with time from cation hydration and growth of diffuse ion layers for both, 1.6 and 1.8 Mg/mseries specimens. Increase in gravimetric water content from DW and SGW permeation increased the kunsat values of 1.6 Mg/m specimens from reduction in total suction. Re-orientation of soil structure mobilized larger kunsat values for specimens permeated with SGW solutions than DW at similar total suction. Permeation of DW and SGW solutions had lesser impact on kunsat values of the 1.8 Mg/m specimens in comparison to the 1.6 Mg/m series specimens. Further at both densities, the influence of permeation was more evident at location closer to hydration surface. Chapter 7 summarizes the main findings of this study.

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.

Gelfill (GF) is made of tailings, water, binder and chemical additives (Fillset, sodium silicate gel). The components of GF are combined and mixed on the surface and transported (by gravity and/or pumping) to the underground mine workings, where the GF can be used for both underground mine support and tailings storage. Thermal (T), hydraulic (H), and mechanical (M) properties are important performance criteria of GF. The understanding of these engineering properties and their evolution with time are still limited due to the fact that GF is a new cemented backfill material. In this thesis, the evolution of the thermal, hydraulic, mechanical, and microstructural properties of small GF samples are determined. Various binder contents of Portland cement type I (PCI) are used. The GF is cured for 3, 7, 28, 90, and 120 days. It is found that the thermal, hydraulic and mechanical properties are time-dependent or affected by the degree of binder hydration index. Furthermore, a relationship is found between the compressive strength and the saturated hydraulic conductivity of the GF samples. The unsaturated hydraulic properties of GF samples have also been investigated. The outcomes show that unsaturated hydraulic conductivity is influenced by the degree of binder hydration index and binder content, especially at low suction ranges. Simple functions are proposed to predict the evolution of air-entry values (AEVs), residual water content, and fitting parameters from the van Genuchten model with the degree of hydration index (α). Furthermore, two columns are built to simulate the coupled thermo-hydro-mechanical (THM) behaviour of GF under drained and undrained conditions. The obtained results from the GF columns are compared with the small samples. It is observed that the mechanical properties, hydraulic properties (suction and water content), and temperature development are strongly coupled. The magnitude of these THM coupling factors is affected by the size of the GF. The findings also show that the mechanical, hydraulic and thermal properties of the GF columns are different from samples cured in plastic moulds.

Le stockage géologique du CO2 dans des réservoirs aquifères de type calcaire et grès, du charbon non exploité est une des solutions envisagées pour réduire les émissions de gaz à effet de serre dans l’atmosphère. Cependant, l’injection de CO2 peut perturber les propriétés pétrophysiques (porosité et perméabilité), minéralogiques (transformations) et mécaniques (déformations, résistance à la rupture) des roches réservoir (calcaire, grès, charbon). Dans le cas du charbon, l’injection de CO2 peut également se traduire par des phénomènes de gonflement de la matrice liés au processus d’adsorption. L’objectif de ce travail de thèse est de traduire en termes de modèles phénoménologiques les comportements et les couplages chimio-poromécaniques des roches réservoir de type charbon. Dans ce travail, nous nous sommes focalisés en particulier sur l’étude de l’injection de CO2 dans le charbon. Pour cela, deux modèles homogénéisés de porosité du charbon ont été développés avec la prise en compte du phénomène d’adsorption, connu pour être le principal mécanisme de production ou de séquestration de CO2 dans de nombreux réservoir de charbon. Le premier modèle permet d’étudier le comportement poro-élastique du charbon pour une injection simple de CO2 et le second permet d’étudier le comportement poro-élastique du charbon pour une injection de CO2 avec une récupération assistée de méthane CH4. Le processus d’adsorption est classiquement modélisé à l’aide de l’isotherme d’adsorption de Langmuir (pour un gaz dans le premier modèle et pour deux gaz dans le second modèle). L’implantation de ces modèles dans le Code_Aster (code d’analyse de calcul de structures entièrement couplé THM, développé par Electricité De France - EDF) nous a permis de faire des simulations numériques de stockage de CO2 dans le charbon. Pour une injection simple du CO2 dans le charbon (premier modèle), la matrice du charbon s’est comportée de deux façons différentes : elle gonfle (ce qui induit une diminution de la porosité du charbon) avec la prise en compte du phénomène d’adsorption et se contracte (ce qui induit une augmentation de la porosité du charbon) dans le cas contraire. Etant en bon accord avec les résultats de la littérature spécialisée, cela montre la capacité du modèle à prédire le comportement poro-élastique du charbon durant l’injection de CO2. Toujours avec le premier modèle, nous avons en particulier étudié l’influence des propriétés hydro-mécaniques du charbon (coefficient de Biot, module de Young/module d’incompressibilité), les paramètres d’adsorption de Langmuir et la pression initiale du liquide interstitiel dans le charbon, sur la réponse du charbon à l’injection du CO2. Dans le cas d’une récupération assistée du méthane CH4 (le second modèle), un couplage du Code_Aster et un code de transport réactif HYTEC (HYdrological Transport coupled with Equilibrium Chemistry, développé par MINES Paris Tech) était nécessaire pour gérer surtout le calcul des pressions partielles des deux gaz (CO2 et CH4) à chaque pas de temps. Un travail de développement numérique sur les deux codes de calcul était alors nécessaire. Ce travail de thèse a proposé une méthode de couplage entre les deux codes (Code_Aster et HYTEC) dont les techniques sont largement décrites dans le manuscrit. / The geological storage of CO2 in aquifers reservoirs such as limestone and sandstone, coal is a possible way to reduce greenhouse gas emission into the atmosphere. However, the injection of CO2 may modify petrophysical (porosity and permeability), mineralogical (transformations) and mechanical (deformations, strength) properties of reservoir rocks (limestone, sandstone, coal). In the case of coal, the injection of CO2 can also induce matrix swelling due to adsorption processes. The focus of this thesis is to translate in terms of phenomenological models, the behaviors and chemo-poromechanical coupling of reservoir rocks of coal type. In this work, we focused particularly on the study of CO2 injection into coal. For this, two models of homogenized coal porosity have been developed by taking into account the adsorption phenomenon, known to be the main mechanism of production or sequestration of CO2 in many coal reservoirs. The first model allows the study of the poroelastic behavior of coal in the case of a single injection of CO2, and the second model allows the study of the poroelastic behavior of coal in the case of an injection of CO2 with methane CH4 recovery. The adsorption process is classically modelled using Langmuir’s isotherm (for one gas in the first model and for two gases in the second model). The implementation of these models in Code_Aster (a fully coupled Thermo-Hydro-Mechanical analysis code for structures calculations, developed by Electricity of France - EDF) allowed us to make numerical simulations of CO2 storage in coal. For a single injection of CO2 into coal (first model), the coal matrix behaved in two different ways: it swells (resulting in the decrease of coal porosity) when the adsorption phenomenon is taken into account and shrinks (resulting in the increase of coal porosity) otherwise. Being in good agreement with the results in specialized literature in this field, it shows the ability of the model to predict the poroelastic behaviour of coal to CO2 injection. Also with the first model, we studied particularly through numerical simulations the influence of coal’s hydro-mechanical properties (Biot’s coefficient, bulk modulus), Langmuir’s adsorption parameters and the initial liquid pressure in rock mass during CO2 injection in coal. In the case of methane recovery (second model), a coupling of Code_Aster and a reactive transport code, HYTEC (Hydrological Transport coupled with Equilibrium Chemistry, developed by Mines Paris Tech) was needed to handle the above calculation of partial pressures of the two gases (CO2 and CH4) at each time step. Digital development work on the two computers codes (Code_Aster and HYTEC) was then necessary. This thesis proposed a method of coupling between the two codes whose techniques are widely described in the manuscript.