The main objective of this research is to present a robust numerical framework based upon Isogeometric analysis (IGA) for simulation of thermo-hydro-mechanical (THM) processes in variably saturated soils. The proposed platform employs the Bézier extraction operator to connect IGA to the conventional finite element analysis (FEA), allowing to take advantage of features offered by the two methods. In the first part, the formulation and numerical implementation for fully coupled numerical simulation of THM problems in saturated porous media are presented. The results are compared against analytical solutions and experimental tests available in the literature. In the second part, the proposed method is used to study the temperature effect on the hydro-mechanical response of sd supporting hydrocarbon pipelines, an aspect that has been overlooked in the majority of previous studies. The results highlight the need for considering nonisothermal behavior in different analysis and design stages of sd-buried pipelines. In the third part, the proposed IGA-FEA framework is extended to evaluate the nonisothermal elasto-plastic behavior of unsaturated soils. Drucker-Prager yield surface is used as criterion to limit the modified effective stress where the model follows small strain, quasi-static loading conditions. The framework is used to simulate strain localization of unsaturated dense sand subjected to undrained compression loading. In comparison with FEA, the present method smoothly distributes plastic strain over the adjacent elements. The parametric study highlights the importance of considering temperature effects in elasto-plastic analysis of unsaturated soils.

Esta Tesis está dirigida al análisis teórico y experimental de problemas acoplados Termo-Hidro Mecánico (THM) que se desarrollan en formaciones geológicas profundas destinadas al almacenamiento de residuos radiactivos de alta actividad. En las últimas décadas, han sido estudiadas las formaciones arcillosas para ser utilizadas como barreras geológicas debido a su reducida conductividad hidráulica. La degradación de las rocas arcillosas producida por efectos de temperatura y por efectos de variación en el grado de saturación, es un factor de fundamental importancia, que es actualmente investigado en ensayos in situ a gran escala, como así también en ensayos de laboratorio. En ésta tesis, la roca Opalinus Clay ha sido ampliamente caracterizada mediante ensayos de laboratorios. Desde un punto de vista macro-estructural se ha obtenido la curva de retención de agua, conductividad hidráulica, resistencia y deformación. El análisis micro-estructural está enfocado a la caracterización mineralógica obtenida por difracción de rayos X, la distribución del tamaño de los poros determinada por porosimetría de mercurio (MIP) y microscopía electrónica (SEM). La tesis describe también un ensayo in situ de calentamiento diseñado para analizar la interacción entre la barrera de ingeniería (bloques de bentonita compactada) y la barrera geológica (Opalinus clay). Esta interacción ha sido analizada a través de simulaciones numéricas realizadas con el código de elementos finitos CODE_BRIGHT. Una célula termo-hidráulica fue especialmente diseñada para observar el comportamiento THM de la roca en condición drenada y no drenada, a través de pulsos de calor. Parámetros térmicos e hidráulicos de la roca fueron determinados por retro análisis a través de simulaciones numéricas realizadas con CODE_BRIGHT. Desde el punto de vista mecánico, un modelo constitutivo ha sido formulado en 3D e implementado en CODE_BRIGHT con el objetivo de reproducir el comportamiento mecánico anisótropo y rotura frágil de las rocas arcillosas. El modelo es formulado en un marco viscoplástico y considera la resistencia y deformabilidad de la matriz y de las juntas. El criterio de falla de la matriz y de las juntas es definido por superficies de fluencias hiperbólicas en el espacio de tensiones p-J y τ−σ, respectivamente. El comportamiento frágil de las rocas arcillosas es simulado por un reblandecimiento isótropo y cinemático definido en términos de trabajo de deformación plástico. El modelo constitutivo ha sido calibrado mediante ensayos triaxiales de laboratorio realizados en especimenes con diferentes ángulos de buzamiento. El modeloconstitutivo anisótropo ha sido aplicado a la simulación numérica en 3D de un ensayo de calentamiento in-situ. Una simulación numérica en 3D de un ensayo de ventilación in-situ realizado en un micro-túnel sin recubrimiento ha sido realizada para reproducir el brusco cambio de permeabilidad por efectos de secado de la roca. En este caso, un modelo hidráulico que considera la apertura de las juntas por efectos de secado ha sido implementado para reproducir los cambios de permeabilidad en excavaciones subterráneas. / This thesis deals with the theoretical and experimental analysis of the coupled Thermo- Hydro-Mechanical (THM) processes developed in geological formations suitable for the repository of radioactive waste of high activity. In the last decades, the argillaceous formations have been studied to be used as geological barriers, due to its reduced hydraulic conductivity. The degradation of clay shales induced by temperature and saturation effects is an important factor which is currently being investigated in large scale in situ tests as well as in laboratory studies. In this thesis, the Opalinus clay rock has been widely characterized by means of laboratory tests. From a macro-structural point of view, the water retention curve, hydraulic conductivity, strength and deformability parameters have been determined. The micro-structural analysis is focused to the mineralogical characterization obtained by means of X ray diffraction, pore size distribution (PSD) determined by means of mercury intrusion porosimetry (MIP) and scanning electronic microscopy (SEM). The thesis describes also a large scale heating in situ test designed to analyze the interaction between the engineer barrier (compacted bentonite blocks) and by the geological barrier, (Opalinus clay). This interaction has been analyzed by means of numerical simulations performed with the finite element code CODE_BRIGHT. A thermo hydraulic cell was specially designed to observe the coupled THM behaviour of the clay shale rock under drained and undrained conditions by means of heat pulses. Thermal and hydraulic parameters of rock were determined by means of back-analysis performed with the help of CODE_BRIGHT. In order to reproduce the anisotropic and brittle behaviour of the clay shale, a 3D mechanical constitutive model has been formulated and implemented in CODE_BRIGHT. The constitutive model is formulated in a viscoplastic framework and it considers the strength and deformability of both matrix and discontinuities (joints). The failure criterion of the matrix and the joints is defined by means of hyperbolic yield surfaces in the p-J and τ-σ stress space, respectively. The brittle behaviour of clay shale is simulated by means of isotropic and kinematic softening defined in terms of a workhardening criterion. The anisotropic constitutive model has been calibrated against triaxial laboratory tests performed on specimens with a main family of discontinuities having different dip angles. The constitutive model has been applied to a 3D numerical simulation of an "in-situ" heating test. A 3D numerical simulation of a ventilation test performed in an unlined micro tunnel was also performed in order to reproduce the changes of the rock permeability by drying effects. In this case, a hydraulic model able to consider the changes in joint thickness by drying effects has been developed to reproduce the changes of permeability in underground excavations.

Glaciation is considered one of the main natural processes that can have a significant impact on the long term performance of DGRs. The northern part of the American continent has been subjected to a series of strong glaciation and deglaciation events over the past million years. Glacial cycles cause loading and unloading, temperature changes and hydraulic head changes at the ground surface. These changes can be classified as transient boundary conditions. It is widely accepted that the periodic pattern of past glacial cycles during the Late Quaternary period are resultant of the Earth’s orbital geometry changes that is expected to continue in the future. Therefore, from the safety perspective of DGRs, such probable events need to be taken into account. The objective of this thesis is to develop a numerical model to investigate the thermo-hydro-mechanical-chemical (THMC) coupled processes that have resulted from long term past and future climate changes and glaciation cycles on a proposed DGR in sedimentary rocks in southern Ontario. The first application is done on a large geological cross section that includes the entire Michigan basin by using a hydro-mechanical (HM) coupled process. The results are compared with field data of anomalous pore water pressures from deep boreholes in sedimentary rocks of southern Ontario. In this work. The modeling results seem to support the hypothesis that at least the underpressures in the Ordovician formation could be partially attributed to past glaciation. The second application is made on site conditions by using the THMC model. The results for the pore water pressure, tracer profiles, permafrost depth and effective stress profile are compared with the available field data, the results show that the solute transport in the natural limestone and shale barrier formations is controlled by diffusion, which provide evidence that the main mechanism of transport at depth is diffusion-dominant. The third application is made on site conditions to determine the effect of underground changes in DGRs due to DGR construction. The results show that future glaciation loads will induce larger increases in effective stresses on the shaft. Furthermore, it is found that hypothetical nuclide transport in a failed shaft can be controlled by diffusion and advection. The simulation results show that the solute transported in a failed shaft can reach the shallow bedrock groundwater zone. These results might imply that a failed shaft will substantially lose its effectiveness as a barrier. The fourth application is proposed to investigate the geochemical evolution of sedimentary host rock in a near field scale. In this part, a new thermo-hydro-mechanical-geochemical simulator (COMSOL-PHREEQC) is developed. It is anticipated that there will be a geochemical reaction within the host rock that results from interaction with the water enriched with the CO2 generated by nuclear waste.

Le traitement à la chaux est une technique qui améliore considérablement la maniabilité et le comportement mécanique des sols à problèmes. Cependant, la durabilité de ce traitement dans les ouvrages en terre sur le long terme représente un enjeu important pour leur stabilité. En outre, la procédure de mise en place, et par conséquent, la taille des agrégats qui en résulte, est un paramètre essentiel qui peut influencer le comportement des sols traités à la chaux utilisés dans le domaine de la construction d’ouvrages. Ce travail de thèse vise à étudier le comportement thermo-hydro-mécanique des sols traités à la chaux, et plus particulièrement les effets du temps de cure et de la taille des agrégats. Des échantillons de sols limoneux et argileux traités à la chaux ont été préparés avec des agrégats de différentes tailles puis soumis à des temps de cures plus ou moins longs. Ces matériaux ont ensuite été étudiés à travers des observations de la microstructure, des analyses minéralogiques, des mesures de la conductivité thermique, de la perméabilité à l’air et de la capacité de rétention d'eau, complétées par de la détermination de la compressibilité et des mesures des modules de cisaillement en petites déformations. Les résultats montrent que le traitement à la chaux modifie de manière significative le comportement thermo-hydro-mécanique des sols. De plus, le comportement des sols traités est fortement influencé par la taille des agrégats. Plus celle-ci est grande, plus la conductivité thermique et la perméabilité à l'air est importante. En revanche, la capacité de rétention en eau est diminuée de même que la compressibilité et la rigidité du sol / Lime treatment is a technique which greatly improves the workability and the mechanical behaviour of problematic soils. However, the sustainability of this treatment in the earthworks for the long term is an important issue for their stability. Besides, the aggregate size resulting from the construction procedure is an essential parameter that may influence the behaviour of treated soils in field construction. The present work deals with the thermo-hydro-mechanical properties of lime-treated soils, with an emphasis put on the curing time and the aggregate size effects. Lime-treated soil samples (both silt and clay) were prepared with different sizes of aggregates and cured during different periods. Afterwards, these soils were studied through microstructural observations, mineralogical analyses, thermal conductivity, air permeability and water retention capacity measurements, as well as the determinations of compressibility and small strain shear modulus. The results show that significant changes of thermo-hydro-mechanical behaviour of soils are induced by lime treatment after curing. Moreover, the aggregate size also plays an essential role in the behaviour of treated soils. Samples prepared with the large aggregates present higher thermal conductivity and air permeability, but with lower water retention capacity, poorer compression behaviour and smaller stiffness

The shear strength of soil-structure interfaces is relevant to the stability of energy piles. The thermo-hydro-mechanical processes can have a strong effect on the behaviour of interfaces between unsaturated soils and piles. Temperature changes lead to water movement in the soil. The moisture loss or gain in the soil causes drying or wetting. In addition, water movement influences the heat transfer properties of the soil. Temperature and moisture content changes affect the magnitude of soil suction in unsaturated soils. Changes in soil suction alter the strength and deformation characteristics of the soil mass and soil-structure interfaces. Similar to the effects of temperature changes, the mechanical loading and the changes in hydraulic conditions in the ground would cause changes in the void ratio, degree of saturation, suction, strength and deformation characteristics of soil. The interface behaviour under varying thermo-hydro-mechanical (THM) conditions is classified as a coupled problem and this is the subject of the present research. In the present investigation, laboratory studies and numerical analyses are carried out to evaluate the THM effect on the behaviour of interfaces between an energy pile material and an unsaturated soil. A 3D interface apparatus (Fakharian and Evgin 1996) has been modified (Fu et al. 2013) to allow the behaviour of an interface to be studied under thermo-mechanical loading conditions. In the present study, the experiments are conducted on soil samples with low degree of saturation and high degree of saturation. It is found that in interface tests using soil samples with low degree of saturation, the adhesion increased due to a positive effect of suction on strength than the negative effect of increasing temperatures. However, in interface tests on soil samples with high degree of saturation, the adhesion decreased with increasing temperatures while the positive effect of suction was not large enough to overcome the negative effect of increasing temperatures. This is a new finding that has not been reported anywhere in the literature. The friction angle for both soil samples (with different degrees of saturation) changed slightly with temperature change. Coupled finite element analyses conducted in the present study provide the following geotechnical information that would be useful for the design of energy piles: (a) Bearing capacity of the pile with and without the effect of temperature, (b) The effect of degree of saturation (or suction) on the strength and deformation characteristics of both the soil and the soil-structure interface, (c) Temperature effects on the amount of pile head movements (up or down), (d) Temperature induced stresses in the pile, (f) Amount of heat that can be stored or extracted from the ground as a function of time. At the initial stages of this study, THM effects on the behaviour of energy piles under saturated and unsaturated conditions are analyzed by using SIGMA/W and VADOSE/W finite element codes of GeoStudio 2012. Although these codes are not multi-physics FE codes, it is possible to use them sequentially to obtain results that will show the trends in pile behaviour. This numerical approach is used first to analyze the behaviour of an energy pile installed partially in unsaturated soil. The moisture content and temperature distributions around a 10 m long, bored pile are calculated using transient analyses. Changes taking place in the stress state along the pile shaft and the bearing capacity of the pile at different temperatures are calculated. In the second part of the numerical analysis of the present study, THM effects on the behaviour of energy piles under saturated and unsaturated conditions are analyzed by using PLAXIS 2D finite element code. PLAXIS is a fully couples finite element code. In order to enhance present understanding of the behaviour of energy piles and do the analysis correctly, a fully coupled analysis involving thermo-hydro-mechanical processes was carried out. Three simulations (mechanical loading only, thermo-mechanical coupling and thermo-hydro-mechanical coupling) are conducted using case studies that are available in the literature. In addition, the behaviour of a generic energy pile, which is installed in a kaolin-sand mixture, is studied by taking into consideration of thermo–hydro-mechanical processes. The coupled analysis provided the distributions of temperature, degree of saturation, suction and heat flux in the analysis domain. Numerical results of the fully-coupled method are compared with the results of sequential method of analysis.

Glaciation is considered one of the main natural processes that can have a significant impact on the long term performance of DGRs. The northern part of the American continent has been subjected to a series of strong glaciation and deglaciation events over the past million years. Glacial cycles cause loading and unloading, temperature changes and hydraulic head changes at the ground surface. These changes can be classified as transient boundary conditions. It is widely accepted that the periodic pattern of past glacial cycles during the Late Quaternary period are resultant of the Earth’s orbital geometry changes that is expected to continue in the future. Therefore, from the safety perspective of DGRs, such probable events need to be taken into account. The objective of this thesis is to develop a numerical model to investigate the thermo-hydro-mechanical-chemical (THMC) coupled processes that have resulted from long term past and future climate changes and glaciation cycles on a proposed DGR in sedimentary rocks in southern Ontario. The first application is done on a large geological cross section that includes the entire Michigan basin by using a hydro-mechanical (HM) coupled process. The results are compared with field data of anomalous pore water pressures from deep boreholes in sedimentary rocks of southern Ontario. In this work. The modeling results seem to support the hypothesis that at least the underpressures in the Ordovician formation could be partially attributed to past glaciation. The second application is made on site conditions by using the THMC model. The results for the pore water pressure, tracer profiles, permafrost depth and effective stress profile are compared with the available field data, the results show that the solute transport in the natural limestone and shale barrier formations is controlled by diffusion, which provide evidence that the main mechanism of transport at depth is diffusion-dominant. The third application is made on site conditions to determine the effect of underground changes in DGRs due to DGR construction. The results show that future glaciation loads will induce larger increases in effective stresses on the shaft. Furthermore, it is found that hypothetical nuclide transport in a failed shaft can be controlled by diffusion and advection. The simulation results show that the solute transported in a failed shaft can reach the shallow bedrock groundwater zone. These results might imply that a failed shaft will substantially lose its effectiveness as a barrier. The fourth application is proposed to investigate the geochemical evolution of sedimentary host rock in a near field scale. In this part, a new thermo-hydro-mechanical-geochemical simulator (COMSOL-PHREEQC) is developed. It is anticipated that there will be a geochemical reaction within the host rock that results from interaction with the water enriched with the CO2 generated by nuclear waste.

With the widening and increasing use of nuclear energy, it is very important to design and build long-term deep geological repositories (DGRs) to manage radioactive waste. The disposal of nuclear waste in deep rock formations is currently being investigated in several countries (e.g., Canada, China, France, Germany, India, Japan and Switzerland). In Canada, a repository for low and intermediate level radioactive waste is being proposed in Ontario’s sedimentary rock formations. During the post-closure phase of a repository, significant quantities of gas will be generated from several processes, such as corrosion of metal containers or microbial degradation of organic waste. The gas pressure could influence the engineered barrier system and host rock and might disturb the pressure-head gradients and groundwater flows near the repository. An increasing gas pressure could also cause damage to the host rock by inducing the development of micro-/macro-cracks. This will further cause perturbation to the hydrogeological properties of the host rock such as desiccation of the porous media, change in degree of saturation and hydraulic conductivity. In this regard, gas generation and migration may affect the stability or integrity of the integrate barriers and threaten the biosphere through the transmitting gaseous radionuclides as long-term contaminants. Thus, from the safety perspective of DGRs, gas generation and migration should be considered in their design and construction. The understanding and modelling of gas migration within the host rock (natural barrier) and the associated potential impacts on the integrity of the natural barrier are important for the safety assessment of a DGR. Therefore, the key objectives of this Ph.D. study include (i) the development of a simulator for coupled modelling of gas migration in the host rock of a DGR for nuclear waste; and (ii) the numerical investigation of gas migration in the host rock of a DGR for nuclear waste in Ontario by using the developed simulator. Firstly, a new thermo-hydro-mechanical-chemical (THMC) simulator (TOUGHREACT-COMSOL) has been developed to address these objectives. This simulator results from the coupling of the well-established numerical codes, TOUGHREACT and COMSOL. A series of mathematical models, which include an elastoplastic-damage model have been developed and then implemented into the simulator. Then, the predictive ability of the simulator is validated against laboratory and field tests on gas migration in host rocks. The validation results have shown that the developed simulator can predict well the gas migration in host rocks. This agreement between the predicted results and the experimental data indicates that the developed simulator can reasonably predict gas migration in DGR systems. The new simulator is used to predict gas migration and its effects in a potential DGR site in Ontario. Valuable results regarding gas migration in a potential DGR located in Ontario have been obtained. The research conducted in this Ph.D. study will provide a useful tool and information for the understanding and prediction of gas migration and its effect in a DGR, particularly in Ontario.

Better understanding and control of crack growth direction during hydraulic fracturing are essential for enhancing productivity of geothermal and petroleum reservoirs. Structural analysis of fracture propagation and impact on fluid flow is a challenging issue because of the complexity of rock properties and physical aspects of rock failure and fracture growth. Realistic interpretation of the complex interactions between rock deformation, fluid flow, heat transfer, and fracture propagation induced by fluid injection is important for fracture network design. In this work, numerical models are developed to simulate rock failure and hydraulic fracture propagation. The influences of rock deformation, fluid flow, and heat transfer on fracturing processes are studied using a coupled thermo-hydro-mechanical (THM) analysis. The models are used to simulate microscopic and macroscopic fracture behaviors of laboratory-scale uniaxial and triaxial experiments on rock using an elastic/brittle damage model considering a stochastic heterogeneity distribution. The constitutive modeling by the energy release rate-based damage evolution allows characterizing brittle rock failure and strength degradation. This approach is then used to simulate the sequential process of heterogeneous rock failures from the initiation of microcracks to the growth of macrocracks. The hydraulic fracturing path, especially for fractures emanating from inclined wellbores and closed natural fractures, often involves mixed mode fracture propagation. Especially, when the fracture is inclined in a 3D stress field, the propagation cannot be modeled using 2D fracture models. Hence, 2D/3D mixed-modes fracture growth from an initially embedded circular crack is studied using the damage mechanics approach implemented in a finite element method. As a practical problem, hydraulic fracturing stimulation often involves fluid pressure change caused by injected fracturing fluid, fluid leakoff, and fracture propagation with brittle rock behavior and stress heterogeneities. In this dissertation, hydraulic fracture propagation is simulated using a coupled fluid flow/diffusion and rock deformation analysis. Later THM analysis is also carried out. The hydraulic forces in extended fractures are solved using a lubrication equation. Using a new moving-boundary element partition methodology (EPM), fracture propagation through heterogeneous media is predicted simply and efficiently. The method allows coupling fluid flow and rock deformation, and fracture propagation using the lubrication equation to solve for the fluid pressure through newly propagating crack paths. Using the proposed model, the 2D/3D hydraulic fracturing simulations are performed to investigate the role of material and rock heterogeneity. Furthermore, in geothermal and petroleum reservoir design, engineers can take advantage of thermal fracturing that occurs when heat transfers between injected flow and the rock matrix to create reservoir permeability. These thermal stresses are calculated using coupled THM analysis and their influence on crack propagation during reservoir stimulation are investigated using damage mechanics and thermal loading algorithms for newly fractured surfaces.

<p>One of many applications of geohydraulic modelling is assessing the suitability of a site to host a nuclear waste repository. This modelling task is complicated by scale-dependent heterogeneity and coupled thermo-hydro-mechanical (THM) processes. The objective here was to develop methods for (i) upscaling flow and transport in fractured media from detailed-scale data and (ii) accounting for THM-induced effects on regional-scale transport. An example field data set was used for demonstration.</p><p>A systematic framework was developed where equivalent properties of flow, transport, and stress-effects were estimated with discrete fracture network (DFN) modelling, at some block scale, and then transferred to a regional-scale stochastic continuum (SC) model. The selected block scale allowed a continuum approximation of flow, but not of transport. Instead, block-scale transport was quantified by transit time distributions and modelled with a particle random walk method at the regional scale.</p><p>An enhanced SC-upscaling approach was developed to reproduce the DFN flow results more simply. This required: (i) weighting of the input well-test data by their conductivity-dependent test volumes and (ii) conductivity-dependent correlation structure. Interestingly, the best-fitting correlation structure resembled the density function of DFN transmissivities. </p><p>Channelized transport, over distances exceeding the block scale, was modelled with a transport persistence length. A linear relationship was found between this persistence length and the macroscale dispersion coefficient, with a slope equal to a representative mean block-scale dispersion coefficient.</p><p>A method was also developed to combine well-test data and rock-mechanical data in estimating fracture transmissivities, and its application was demonstrated.</p><p>Finally, an overall sequential THM analysis was introduced allowing the estimation of the significance of waste-related thermo-mechanical (TM) effects on regional transport; here TM effects are calculated separately and their impact on fracture transmissivities were incorporated into the hybrid framework. For the particular case, their effects on regional-scale transport were small.</p>

In this thesis, fractured rock masses are treated asequivalent continua for large-scale analyses of rockengineering projects. Systematic developments are made for thedetermination of equivalent mechanical and hydraulic propertiesof fractured rock masses using a hybrid discrete fracturenetwork - distinct element method (DFN-DEM) approach. Thedetermined equivalent properties are then used for a far-fieldfinite element analysis of the thermo-mechanical impacts on thestress, deformation and permeability of fractured rockssurrounding a hypothetical geological repository of nuclearwaste. The geological data were extracted from the results ofan extensive site investigation programme at Sellafield, UK,conducted by Nirex UK Ltd. The scale dependencies of the hydraulic and mechanicalproperties were investigated by using multiple realizations ofthe fracture system geometry with increasing model sizes untilproperly defined hydraulic and mechanical representativeelementary volumes (REVs) were reached. The validity of thesecond order permeability tensor and the fourth-ordermechanical compliance tensor were tested for continuum analysesat larger scales. The REV was determined to be around 5 m formechanical and hydraulic data in this study. Analysis of the stress-dependent mechanical and hydraulicproperties shows that the effect of rock stresses is crucial.The elastic moduli increase significantly with the increase ofstress and an empirical equation of stress-dependent elasticmodulus is suggested based on results of numerical experiments.Calculations of the Poisson's ratios suggest greater valuesthan are normally assumed in practice. Depending on the stateof stress, permeability decreases or increases with increasingcompressive stress. Stress-induced flow channeling effect iscaptured by numerical modeling for the first time and detailedmechanisms of shear dilation of fractures are provided. Basedon the numerical experiments, a set of empirical equations wassuggested for the stress-dependent permeability, consideringboth normal deformation and shear dilation of fractures. Thermo-mechanical impact on the performance of ahypothetical repository at a far-field scale (5 km by 1 km) wasinvestigated with the stress-dependent equivalent propertiesdetermined at the REV scale. This analysis shows thatmechanical responses vary significantly depending on how themechanical properties were determined. The change ofpermeability due to the thermal loading is, however, notsignificant in this particular case. The thesis provides a framework for systematic analysis oflarge-scale engineering applications in fractured rock masses,such as geological repositories of nuclear wastes. Keyword:Fractured rock masses, Equivalent Continuum,Discrete Fracture Network (DFN), Distinct Element Method (DEM),Finite Element Method (FEM), Nuclear Waste Disposal, CoupledThermo-Hydro-Mechanical Processes