In Belgium, the Boom Clay is currently studied as the reference host formation for the disposal of high-level and long-lived radioactive waste. In case of direct disposal of spent fuel, uranium isotopes are important contributors along with their daughters to the dose rate at very long term. Therefore, it is essential to study the migration of uranium in the host formation. The present work contributes to improve the knowledge of uranium speciation in the Boom Clay, U(IV) versus U(VI), and of the mechanisms controlling the uranium mobility such as solubility, sorption and complexation by organic matter. The information necessary to interpret the migration behaviour is derived from the study of naturally occurring uranium in the rock and from laboratory experiments conducted under conditions representative for the Boom Clay. Uranium naturally present in the Boom Clay is concentrated in detrital heavy minerals and in authigenic iron(II)-bearing minerals such as siderite and glauconite. Despite its reducing capacity, pyrite is surprisingly depleted in uranium relative to the mean content. Furthermore, uranium is also associated with the surfaces of clay minerals. The clayey fraction contains about 4 ppm uranium and is the main contributor to the total uranium content of the Boom Clay since it constitutes up to 60 wt. % of the rock. The correlation observed on the field between uranium and organic matter suggests that uranium is reduced, likely during the early diagenesis process of bacterially-mediated sulphate reduction. If hexavalent oxidation state of uranium predominates as predicted by geochemical calculations based on the most recent thermodynamic data of the Nuclear Energy Agency (NEA), less than 5% of uranium is complexed by humic acids in the Boom Clay pore water. The U(VI) speciation is dominated by the inorganic carbonate complexes, merely UO2(CO3)34-. The conditional constant determined for the complexation of U(VI) by humic acids under in situ Boom Clay conditions is log exp = 12.4. However, experimental studies show that UO2(CO3)34- is reduced by interaction with pyrite, the main reducing mineral present in the rock, and precipitates as a mixed oxide of U(IV)/U(VI), i.e. UO2+x. Moreover, electromigration experiments suggest that U(VI) is not stable in the Boom Clay: U(VI) is reduced and precipitates as U(IV) oxy-hydroxides. The experimentally measured solubility of U(IV) amorphous oxide, UO2(am), in Boom Clay pore water is about 10 8 mol•l 1. This solubility value is not increased by complexation of U(IV) with dissolved organic matter. The dominant effect of organic matter on the dissolution of UO2(am) is the stabilisation of U(IV) real colloids which increase the uranium concentration by three orders of magnitude. However, the mobility of these colloids is expected to be very limited because of the compaction level of the Boom Clay and its ultra-filtrating feature. The diffusive transport of dissolved uranium is furthermore retarded by significant sorption onto clay minerals. In conclusion, the presence of organic matter in the Boom Clay has no negative effect on the uranium retention which is dominated by the solubility and sorption of U(IV) species

Le stockage des déchets radioactifs dans les formations géologiques profondes peu perméables comme les argilites et les argiles plastiques est envisagée comme une solution possible et fait l'objet de nombreuses études depuis une trentaine d'années. Dans le cadre du projet européen TIMODAZ, l'accent a été mis sur l'étude des effets d'une augmentation de la température engendrée par les déchets exothermiques sur la zone endommagée autour d'une galerie souterraine de stockage. Dans le cadre de ce projet, une étude expérimentale sur le comportement thermique de l'argile de Boom et de l'argile à Opaline a été réalisée. Afin de surmonter les difficultés reliées à l'étude expérimentale des matériaux peu perméables en laboratoire, une nouvelle cellule triaxiale à court chemin de drainage a été mise en œuvre. Les essais ainsi qu'une modélisation numérique montrent que la re-saturation des échantillons désaturés par le processus d'excavation, transport, stockage et préparation peut être réalisée beaucoup plus rapidement par ce dispositif. Les essais de chargement mécanique et thermique en condition drainée (c'est-à-dire avec une surpression interstitielle engendrée négligeable) peuvent être réalisés également dans cette cellule avec des vitesses de chargement plus élevée comparée aux cellules triaxiales classiques. La possibilité de réactivation d'une bande de cisaillement par pressurisation thermique du fluide interstitiel dans un échantillon de l'argile de Boom est mise en évidence. On observe qu'un plan de rupture préexistant dans l'échantillon agit comme un plan de faiblesse pouvant être réactivé de façon préférentielle au moment de la rupture. La résistance au cisaillement obtenue sur le plan de rupture est inférieure à celle de matériau intact pour l'argile de Boom. Le comportement thermique de l'argile à Opaline a été étudié à partir d'essais de chauffage en condition drainée et non drainée sur des échantillons saturés. L'essai de chauffage drainé montre un comportement thermo-elasto-plastique avec limite expansion/contraction à 65°C. Ce comportement est similaire au comportement des argiles faiblement surconsolidées. L'analyse des résultats de l'essai de chauffage non drainé met en évidence que l'eau interstitielle dans l'argile à Opaline a un coefficient de dilation thermique plus important comparé à celui de l'eau libre. Dans la gamme de températures étudiées (25°C-80°C), les mesures de perméabilité sur les échantillons endommagés par un chargement déviatorique ne montrent aucun effet de l'endommagement sur la perméabilité, ce qui prouve la bonne capacité de scellement de l'argile de Boom et l'argile à Opaline saturées / Storage of exothermic radioactive waste in deep low permeability geological formations such as clayey rocks and plastic clays is a solution considered for long term repositories. However the excavation of underground galleries creates a damaged zone (EDZ). The effect of the damage zone on the transport properties of the geological barrier has been widely studied. Within the framework of the TIMODAZ European project, emphasis has been put on the effect of temperature. As a partner of this project, the current work is performed to investigate the coupling effect between temperature, damage and permeability on Boom clay and Opalinus clay through an experimental study. View to the experimental difficulties related to the low permeability materials, a new hollow cylinder triaxial cell with short drainage path specifically designed to study the thermo-hydro-mechanical behaviour of very low permeable materials is developed during this work. The tests and the numerical analysis show that the short sample drainage path reduces significantly the time needed to resaturate an initially unsaturated sample and it also permits to achieve drained conditions (i.e. negligible excess pore pressure during testing) with a higher loading rate. For Boom clay, the effect of the pore water thermal pressurisation on a sample with a pre-existing shear band is investigated. The undrained heating under shear stress decreases the effective stress on the sample which leads to its failure. An existing failure plane in the sample behaves like a preferential weakness plane which can be reactivated by pore water thermal pressurisation. The estimated shearing resistance along the sheared plane is smaller than that of the intact material. For the Opalinus claystone, drained heating on a saturated sample shows that this claystone behaves like a slightly overconsolidated material (thermo-elasto-plastic behaviour) with transition from expansion to contraction at 65°C. The decrease of the permeability of the sample before and after the heating-cooling cycles proves the irreversible volumetric compaction of the sample. The undrained heating test on the Opalinus claystone induces an excess pore pressure in the sample which cannot be explained by the difference between the free water thermal expansion coefficient and that of the solid matrix. The back analysis of the results shows a higher value for the water thermal expansion coefficient as compared to that of free water. Permeability measurements at 25°C and 80°C on samples previously damaged by deviatoric loading show no significant effect of damage on the permeability of the samples. These results confirm the good sealing capacity of both clays at various temperatures

La présente étude vise à comprendre le comportement chimico-hydro-mécanique des argiles raides à travers deux formations géologiques, l'argile de Boom et les argiles yprésiennes, qui sont censées être des formations hôtes potentielles pour le stockage de déchets radioactifs en Belgique. Le comportement volumique a été étudié tant à l'état intact qu'à l'état reconstitué, et sous différentes conditions : K0 et isotrope, sous des boucles de chargement - déchargement. Les résultats obtenus montrent que le comportement volumique de ces argiles est gouverné par la compétition entre l'effet physico-chimique et l'effet mécanique, caractérisée par une contrainte seuil qui correspond à la contrainte de gonflement en termes de changements de structure. Une loi de comportement volumique a été ainsi développée afin de décrire cet aspect. La perméabilité a été déterminée, comparée avec les résultats dans la littérature et corrélée avec les paramètres comme l'indice des vides. La variation de la perméabilité avec la profondeur ont mis en évidence le rôle déterminant des macro-pores dans le transfert des fluides. Le comportement volumique et la perméabilité des argiles de Boom et yprésiennes intactes sont aussi influencés par la variation de la composition chimique de l'eau de pore, qui modifie la double couche diffuse et favorise l'agrégation des particules argileuses. Les caractéristiques élastiques, la surface de charge et l'enveloppe de rupture ont été identifiées pour le comportement déviatorique des argiles de Boom et yprésiennes. Un modèle élasto-plastique conceptuel a été développé permettant de tenir compte des effets du gonflement et de la compétition entre l'effet mécanique et l'effet physico-chimique / The present research aims to understand the chemo-hydro-mechanical behavior of stiff clays through two geological formations, the Boom Clay and the ypresian clays which are considered as possible host formations for the radioactive wastes disposal in Belgium. The volume change behavior was studied in both intact and reconstituted states, and under different conditions: under K0 and isotropic loading, under loading/unloading loops. The results show that the volume change behavior is governed by the competition between the physico-chemical effect and the mechanical effect, characterized by a threshold stress which corresponds to the swelling stress in terms of structure changes. A constitutive law was developed to capture this aspect. The permeability was determined, compared with the results in literature and correlated with the parameters as void ratio. The permeability variation with depth shows the important role of macro-pores in fluids' transfer. The volume change behavior and permeability of intact Boom Clay and ypresian clays are also influenced by pore water chemical composition changes which modify the diffuse double layer and give rise to the aggregation of clay particles. The elastic parameters, yield curve and failure envelope of Boom Clay and ypresian clays were identified. A conceptual elasto-plastic model was developed, accounting for the swelling effects and the competition between the physico-chemical effect and the mechanical effect

Dans le cadre de la déposition géologique profonde des déchets radioactifs, l'argile de Boom est choisie comme une des formations hôtes potentielles dans le programme belge. Par sa formation géologique, cette argile est considérée comme un matériau isotrope transversal. En effet, l'anisotropie de ses propriétés hydraulique et thermique a été mise en évidence dans plusieurs études. Il apparait maintenant nécessaire de mener une étude complète sur son comportement anisotrope. Sur le plan expérimental, l'anisotropie des propriétés thermo-hydro-mécaniques de l'argile de Boom a été mise en évidence à l'aide des mesures de la conductivité thermique, de la conductivité hydraulique, et du module de cisaillement. Grâce à ces mesures sur des carottes fraichement forées, l'endommagement dû à l'excavation de la galerie a été identifié (dans la zone près de la galerie « Connecting » à Mol, Belgique) et modélisé par un modèle empirique dont la variable d'endommagement est basée sur le volume des macro-pores. Ensuite, l'anisotropie du comportement en condition K0 a été étudiée en analysant le comportement de trois types d'éprouvettes différents (0°, 45° et 90° au plan du litage) pendant et après la re-saturation dans des cellules oedométriques et triaxiales à plusieurs états de contraintes différents. Enfin, le comportement anisotrope en condition triaxiale a été étudié en réalisant des essais triaxiaux sur trois types d'éprouvettes. Au cours de ces essais, les mesures des vitesses d'ondes Vs et Vp ont été également effectuées à l'aide des éléments piézoélectriques encastrés dans deux embases supérieure et inférieure d'une cellule triaxiale. Ces résultats ont permis de déterminer les paramètres de l'élasticité anisotrope. Sur le plan de modélisation, un modèle anisotrope avec huit paramètres anisotropes (cinq paramètres élastiques et trois paramètres plastiques) a été développé, en utilisant la théorie de Bohler (Boehler et Sawczuk, 1977). Ce modèle est basé sur un modèle élasto-plastique isotrope à deux surfaces de charge validé précédemment pour l'argile de Boom (Hong, 2013). La validation de ce modèle anisotrope a été réalisée grâce aux résultats expérimentaux des essais oedométriques et triaxiaux obtenus dans cette étude / In the program of deep geological radioactive waste disposal in Belgium, Boom Clay has been chosen as one of the potential host rocks. Due to the geological stratification, this stiff clay has been regarded as a transverse isotropic material. The anisotropy of its hydraulic and thermal properties was shown in several studies. It seems necessary now to conduct a more in-depth study on the anisotropic behaviour of Boom Clay. In terms of experimental works, the anisotropy of the thermo-hydro-mechanical properties of Boom Clay was evidenced using measurements of thermal conductivity, hydraulic conductivity and small-strain shear modulus. Through these measurements on freshly cored Boom Clay samples, the damage due to excavation of the gallery was identified (in the zone near the Connecting gallery at Mol, Belgium) and modelled using an empirical model in which the damage variable is based on the volume of macro-pores. Afterwards, the anisotropy behaviour under K0 condition was studied by analysing the behaviour of three types of specimens (0°, 45°, and 90° to the bedding plane) during and after the re-saturation in the triaxial and oedometer cells under different stress states. Finally, the anisotropic behaviour under triaxial condition was investigated through several triaxial tests on three types of specimens. During these tests, the velocity measurements of seismic waves Vs et Vp were performed thanks to the bender elements installed in the upper and bottom bases of a triaxial cell. These results were used to determine the parameters of anisotropic elasticity. In terms of modelling works, an anisotropic model with eight anisotropic parameters (fives elastic parameters and three plastic parameters) was developed using the theory of Boehler (Boehler et Sawczuk, 1977). This model is based on an elasto-plastic isotropic model with two yield surfaces elaborated previously for Boom Clay (Hong, 2013). The validation of this anisotropic model was made based on the results obtained from oedometer and triaxial tests performed in this study

Les Argiles de Boom ont été identifiées par le SCK-CEN comme un éventuel site de stockage de déchets nucléaires en couche géologique profonde : elles font l’objet d’études dans le laboratoire souterrain de Mol (Belgique). Dans ce contexte, il est important d’évaluer les conséquences du creusement de galeries sur les propriétés de ces Argiles. Ce travail de thèse cible plus particulièrement les effets d’oxydation à l’air et de biodégradation anaérobie sur la MO fossile. Les expériences d’oxydation ont combiné des suivis expérimentaux (série artificielle) et des études sur des échantillons altérés in situ (série naturelle) dans les galeries du laboratoire. Elles ont ciblé le kérogène, la MOE et MOD. La confrontation des données des deux séries révèle que nos simulations expérimentales sont représentatives des processus ayant lieu dans les galeries. Ces travaux démontrent aussi que l’oxydation induit (i) une augmentation importante de la quantité de MOD et (ii) une modification majeure de la chimie de la MOE et de la MOD caractérisée par la formation de molécules oxygénées de faible poids moléculaire. Par ailleurs, l’étude comparative des eaux issues des échantillons altérés avec celles prélevées dans les piézomètres du site démontre que ces dernières ne sont pas affectées par des processus d’oxydation et sont comparables aux eaux issues des échantillons sains. Des expériences de biodégradation menées sur des argiles saines et préalablement oxydées artificiellement n’ont montré aucune évolution significative de la MO fossile (MOE et MOD); la biodégradation anaérobie n’est donc pas un processus dominant dans le contexte des perturbations induites par les excavations / The Boom Clay was focused because it is identified by SCK-CEN as a possible radioactive waste storage in the geological disposal site and in situ experiments are performed in the underground laboratory of Mol (Belgium). In this context, it is important to assess the consequences of galleries excavation on the properties of the Boom Clay. The particular focus of this study is the effects of air oxidation and anaerobic biodegradation on the OM. The experiments dealing with the effects of air oxidation have combined studies of artificial oxidized samples (artificial series) and samples altered in the gallery of the underground laboratory (natural series). These experiments focus on the Kerogen, the EOM and the DOM. The comparison of data from artificial and natural series shows firstly that our experimental simulations are the representative of processes taking place in the galleries. These studies show that air oxidation induced (i) a significant increase in the amount of DOM and (ii) a major change in the chemistry of the EOM and DOM characterized by the formation of low molecular weight oxygenated molecules. Moreover, comparison between water extracted from altered samples and piezometers shows that the water of the site is not affected by oxidation processes. The piezometer water samples are similar to that extracted from non-altered samples. Biodegradation experiments conducted on non altered and artificially oxidized clay did not show significant changes of fossil and dissolved organic matter. It seems that anaerobic biodegradation is not a major process in the context of disturbances induced by the excavation

Cette étude est consacrée à la modélisation du comportement thermo-mécanique des argiles raides saturées et au développement d'un algorithme d'intégration efficace de contrainte correspondant. Le comportement mécanique de l'argile de Boom naturelle dans des conditions isothermes a été caractérisé. Le modèle Cam Clay modifié (MCC) a été ensuite appliquée pour simuler le comportement de l'argile de Boom naturel. Il a été constaté que le MCC donne des prédictions de mauvaise qualité pour le comportement de l'argile de Boom naturel. Ainsi, un modèle Cam Clay (ACC-2) adapté a été développé en introduisant une nouvelle surface de charge et un nouveau potentiel plastique ainsi que d'un mécanisme plastique de Deux surfaces. Ce modèle permet la description satisfaisante des caractéristiques principales du comportement mécanique de l'argile de Boom naturelle. De plus, les équations de ce modèle peuvent être formulées mathématiquement comme dans un modèle élasto-plastique classique. L'algorithme d'intégration de contrainte classique peut donc être appliqué. Les effets thermiques ont été examinés par l'évaluation de la pertinence de trois lois thermomécaniques avancées (Cui et al, 2000; Abuel-Naga et al, 2007; Laloui et François, 2008; 2009). Il apparaît que tous les trois modèles peuvent décrire les caractéristiques principales du comportement thermo-mécanique des argiles saturées. Cependant, chaque modèle a ses limites ou des points peu clairs du point de vue théorique. L'algorithme d'intégration de contrainte du modèle thermo-mécanique de Cui et al. (2000) au point de contrainte a également été développé spécifiquement en utilisant une méthode adaptive du pas de temps. Le temps de calcul nécessaire pour obtenir une précision donnée est ainsi largement réduit pour des chemins de chargements thermiques et mécaniques. Un modèle thermo-mécanique à Deux surfaces (modèle TEAM) a été développé en se basant sur le mécanisme plastique de Deux surfaces. Le modèle proposé a étendu le modèle de Cui et al. (2000) à une formulation de Deux surfaces considérant le couplage entre les déformations plastiques des chemins de chargements thermiques et mécaniques. La simulation des essais drainés montre que ce modèle peut décrire les caractéristiques principales thermo-mécaniques de l'argile de Boom naturelle le long de différents chemins de chargements. Le modèle TEAM a finalement été étendu à des conditions non drainées. Après la clarification du concept des contraintes effectives et la définition d'une condition de déformation volumique, le processus d'échauffement non drainé est analysé. La validité des équations thermo-hydro-mécaniques de ce modèle a été examinée en se basant sur des résultats d'essais typiques / This study is devoted to the thermo-mechanical constitutive modeling for saturated stiff clays and the development of a corresponding efficient stress integration algorithm. The mechanical behavior of natural Boom Clay in isothermal conditions was first characterized. The Modified Cam Clay model (MCC) was then applied to simulate the natural Boom Clay behavior. It has been found that the MCC gives poor-quality predictions of the natural Boom Clay behavior. Thereby, an adapted Cam Clay model (ACC-2) was developed by introducing a new yield surface and a new plastic potential as well as a Two-surface plastic mechanism. This model allows satisfactory prediction of the main features of the mechanical behavior of natural Boom Clay. Moreover, the constitutive equations of this model can be formulated mathematically as in a classic elasto-plastic model. Thus, the classic stress integration algorithm can be applied. The thermal effects were considered by assessing the performance of some advanced thermo-mechanical models (Cui et al., 2000; Abuel-Naga et al., 2007; Laloui and François, 2008; 2009). It appears that all the three models can capture the main features of the thermo-mechanical behavior of saturated clays. However, each constitutive model has its own limitations or unclear points from the theoretical point of view. The stress integration algorithm of the thermo-mechanical model proposed by Cui et al. (2000) at the stress point level was also developed using a specifically designed adaptive time-stepping scheme. The computation time required to achieve a given accuracy is largely reduced with the adaptive sub-stepping considered for both mechanical and thermal loadings. A Two-surface thermo-mechanical model (TEAM model) was developed based on the Two-surface plastic mechanism. The proposed model extends the model of Cui et al. (2000) to a Two-surface formulation, considering the plastic strain coupling between the thermal and the mechanical loading paths. The simulation of drained tests shows that this model can capture the main thermo-mechanical features of natural Boom Clay along different loading paths. The TEAM model was finally extended to undrained conditions. After setting up an appropriate effective stress principle and defining a volumetric strain condition, the undrained heating process was analyzed. The validity of the thermo-hydro-mechanical constitutive equations was examined based on the data from typical tests

En el caso del almacenamiento de los residuos radioactivos los flujos osmóticos pueden ser relevantes y requieren un análisis en detalle. El residuo nuclear bituminizado (BW) será almacenado mediante contenedores en cavidades excavadas en la Boom Clay, que es una arcilla marina que presenta propiedades favorables para limitar y retrasar la migración de los contaminantes radioactivos. La interacción entre los dos materiales es un proceso acoplado químico-hidro-mecánico y depende de la respuesta hidromecánica de la Boom Clay y del BW. En condiciones de almacenamiento, el contacto del BW, que contienen cantidades importantes de NaNO3, con el agua subterránea induce la hidratación por gradientes osmóticos y el consiguiente hinchamiento, además de la difusión de la sal disuelta hacia la Boom Clay. Se pueden distinguir dos tipos de afecciones: la perturbación geomecánica causada por el hinchamiento del BW y el aumento de presión en el BW y cambio de las distribución de tensiones en la roca, y la perturbación físico química por la migración de grandes cantidades de sales. El objetivo de esta tesis es: (i) Mejorar la comprensión de los procesos que controlan la absorción de agua y el consecuente hinchamiento del BW que contengan sales (NaNO3), y (ii) Investigar los posibles efectos de la concentración de fluidos de los poros sobre el hinchamiento, la compresibilidad y comportamiento de corte de la Boom Clay. En primer lugar, se ha desarrollado una formulación para el análisis de la deformación inducida por la disolución de sales en medio poroso con contacto con agua. Las ecuaciones planteadas incluyen los flujos acoplados de agua y soluto. Se presenta también un trabajo teórico que ayuda a la comprensión del comportamiento mecánico del BW. Se considera este material como una mezcla de bitumen y cristales de NaNO3. Se ha desarrollado un modelo elasto-viscoplástico que describe el comportamiento de fluencia del BW considerando el comportamiento de fluencia de sus constituyentes. El modelo constitutivo elasto-viscoplástico ha sido implementado en el programa CODE_BRIGHT. Los resultados se han comparado con observaciones experimentales. Se ha estudiado el comportamiento a largo plazo del BW en contacto con agua al simular ensayos de hinchamiento por absorción de agua bajo condiciones confinadas. El análisis numérico ha demostrado ser capaz de proporcionar una representación satisfactoria de los principales patrones observados en su comportamiento. En lo que respecta al segundo objetivo de la tesis, se ha propuesto una formulación para el análisis de las deformaciones inducidas por procesos osmóticos en un medio poroso de doble estructura. Esta formulación distingue dentro del material un nivel micro-estructural y otro macro-estructural con cambios químicos que tienen un efecto significativo en la micro-estructura. Se han obtenido las ecuaciones básicas que describen los flujos acoplados de agua y solutos y el transporte de sus componentes a través de los macroporos así como las ecuaciones de balance de masa para agua y soluto en los macroporos y microporos. La formulación propuesta ha sido aplicada particularmente para analizar cualitativamente el efecto de la succión osmótica sobre el hinchamiento de los suelos arcillosos. Se han analizado los efectos a corto y largo plazo. Se ha investigado también la influencia del aumento de la concentración del fluido en los poros sobre las propiedades geotécnicas y el comportamiento de la Boom Clay no saturada. Se ha llevado a cabo un programa sistemático de investigación experimental, con control de succión osmótica y matricial, con el fin de investigar el efecto del incremento de la concentración del fluido de poros sobre la resistencia de corte y el cambio volumétrico bajo condiciones edométricas. Se ha observado, que bajo condiciones parcialmente saturadas, un cambio en la salinidad provoca una disminución en la compresibilidad y en la resistencia de corte del material. / For deep storage of high-level nuclear waste osmotic flows can be significant and so require a careful analysis. In Belgium, The bituminized nuclear waste (BW) named Eurobitum contained in metallic drums will be placed inside a tunnel or a shaft excavated in the Boom Clay, which is 100 m thick marine clay presenting favourable properties to limit and delay the migration of the leached radionuclides over extended periods of time. In Geological disposal conditions, contact of the bituminized radioactive waste which contains high amounts of highly soluble salt (NaNO3) with groundwater will result in water uptake and swelling of the waste and in subsequent diffusion of the dissolved salt through the host clay formation. Basically, two types of disturbance can be distinguished: A geo-mechanical perturbation, caused by the swelling of the waste and the increase of the pressure in and around the waste and a physico-chemical perturbation by the release of large amounts of NaNO3 and other soluble salts. In this context the aim of this thesis is: (i) to improve the understanding of the processes controlling the water uptake and the subsequent swelling of bituminized waste containing soluble salts (NaNO3), and (ii) to investigate of the possible effects of the increase of pore fluid concentration on swelling, compressibility and shear behaviour of Boom Clay. A formulation has been proposed for the analysis of deformation induced by dissolution of salts in porous media in contact with water. The equations include the effect of coupled transport phenomena and the formulation has been included as an extension in the coupled THM program CODE_BRIGHT. A theoretical and experimental work aiming at understanding the mechanical behaviour of the Bituminized Waste has been presented.This material is considered for this purpose as a mixture of bitumen and crystals of NaNO3. An elasto-viscoplastic model has been developed that describes the creep behaviour of BW considering the constituents' creep behaviour. The elasto-viscoplastic constitutive model has been implemented into CODE_BRIGHT. The modelling results have been compared with the experimental data. The impact of osmotic forces on the swelling of the material has been investigated by simulating water uptake swelling tests under confined conditions and comparing the predictions with experimental results. The numerical analysis has proven to be able to furnish a satisfactory representation of the main observed patterns of the behaviour. In regard to the second objective of this thesis, a formulation has been proposed for the analysis of deformations induced by osmotic processes in double structure porous media. The formulation is based on the distinction within the material of a microstructural and a macrostructural levels with chemical changes having a significant effect on the microstructure. A macroscopic description of the system is provided. Then the basic equations describing coupled flows of water and solutes and the transport of its components through macropores and mass balance equations for water and solute in macro and micro pores have been obtained. The proposed formulation has been particularly applied to analyze qualitatively the effect of osmotic suction on swelling of clayey soils. Transient and long term effects have been analyzed. The influence of pore fluid concentration on the geotechnical properties and behavior of Boom Clay under partially saturated conditions has been investigated. A systematic experimental research program involving osmotic suction and matric suction controlled experiments has been carried to investigate the effect of the increase of pore fluid concentration on shear strength and on the volume change behaviour under odometer stress state conditions. It has been observed that under partially saturated conditions a change in salinity causes a decrease in compressibility and shear strength.

This study is devoted to the thermo-mechanical constitutive modeling for saturated stiff clays and the development of a corresponding efficient stress integration algorithm. The mechanical behavior of natural Boom Clay in isothermal conditions was first characterized. The Modified Cam Clay model (MCC) was then applied to simulate the natural Boom Clay behavior. It has been found that the MCC gives poor-quality predictions of the natural Boom Clay behavior. Thereby, an adapted Cam Clay model (ACC-2) was developed by introducing a new yield surface and a new plastic potential as well as a Two-surface plastic mechanism. This model allows satisfactory prediction of the main features of the mechanical behavior of natural Boom Clay. Moreover, the constitutive equations of this model can be formulated mathematically as in a classic elasto-plastic model. Thus, the classic stress integration algorithm can be applied. The thermal effects were considered by assessing the performance of some advanced thermo-mechanical models (Cui et al., 2000; Abuel-Naga et al., 2007; Laloui and François, 2008; 2009). It appears that all the three models can capture the main features of the thermo-mechanical behavior of saturated clays. However, each constitutive model has its own limitations or unclear points from the theoretical point of view. The stress integration algorithm of the thermo-mechanical model proposed by Cui et al. (2000) at the stress point level was also developed using a specifically designed adaptive time-stepping scheme. The computation time required to achieve a given accuracy is largely reduced with the adaptive sub-stepping considered for both mechanical and thermal loadings. A Two-surface thermo-mechanical model (TEAM model) was developed based on the Two-surface plastic mechanism. The proposed model extends the model of Cui et al. (2000) to a Two-surface formulation, considering the plastic strain coupling between the thermal and the mechanical loading paths. The simulation of drained tests shows that this model can capture the main thermo-mechanical features of natural Boom Clay along different loading paths. The TEAM model was finally extended to undrained conditions. After setting up an appropriate effective stress principle and defining a volumetric strain condition, the undrained heating process was analyzed. The validity of the thermo-hydro-mechanical constitutive equations was examined based on the data from typical tests