A numerical investigation is conducted on the impacts of the thermal loading history on the evolution of mechanical response and permeability field of a fractured rock mass containing a hypothetical nuclear waste repository. The geological data are extracted from the site investigation results at Sellafield, England. A combined methodology of discrete and continuum approaches is presented. The results of a series of simulations based on the DFN-DEM (Discrete Fracture Network-Distinct Element Method) approach provide the mechanical and hydraulic properties of fractured rock masses, and their stress-dependencies. These properties are calculated on a representative scale that depends on fracture network characteristics and constitutive models of intact rock and fractures. In the present study, data indicate that the large scale domain can be divided into four regions with different property sets corresponding to the depth. The results derived by the DFN-DEM approach are then passed on to a large-scale analysis of the far-field problem for the equivalent continuum analysis. The large-scale far-field analysis is conducted using a FEM code, ROCMAS for coupled thermo-mechanical process. The results show that the thermal stresses of fractured rock masses vary significantly with mechanical properties determined at the representative scale. Vertical heaving and horizontal tensile displacement are observed above the repository. Observed stress and displacement fields also shows significant dependency on how the mechanical properties are characterized. The permeability changes induced by the thermal loading show that it generally decreases close to the repository. However, change of permeability is small, i.e., a factor of two, and thermally induced dilation of fracture was not observed. Note that the repository excavation effects were not considered in the study. The work presented in this paper is the result of efforts on a benchmark test (BMT2) within the international co-operative projects DECOVALEX III and BENCHPAR.

Initiated in 1992, the DECOVALEX project is an international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical (THM) processes in geologic systems. The project has made important scientific achievements through three stages and is progressing in its fourth stage. It has played a key role in the development of mathematical modeling and in situ testing of coupled THM processes in fractured rock and buffer/backfill materials, a subject of importance for performance assessment of radioactive waste geologic repositories. This paper summarizes studies under the most recent stage of the project, DECOVALEX III (2000-2003). These studies include those of two major field experiments: (a) the FEBEX experiment at Grimsel, Switzerland, investigating coupled THM processes in a crystalline rock-bentonite system, and (b) the Drift Scale Test (DST) experiment at Yucca Mountain, Nevada, investigating coupled THM processes in unsaturated tuff. These are two of the largest multiyear heater tests undertaken to date for the study of coupled THM processes in geological systems. In addition, three so-called benchmark tests are also studied to evaluate the impact of coupled THM processes under different scenarios and geometries. Within the DECOVALEX project, multiple research teams participated in each of the studies, using different approaches and computer codes. Comparisons of results have provided insight into coupled THM processes, which in turn has stimulated further development of our modeling capabilities. Lessons learned from these studies are discussed. The scientific advances and enhanced insight gained through this kind of international cooperation illustrate the effectiveness of the DECOVALEX project.

We investigate the stress-dependent permeability issue in fractured rock masses considering the effects of nonlinear normal deformation and shear dilation of fractures using a two-dimensional distinct element method program, UDEC, based on a realistic discrete fracture network realization. A series of "numerical" experiments were conducted to calculate changes in the permeability of simulated fractured rock masses under various loading conditions. Numerical experiments were conducted in two ways: (1) increasing the overall stresses with a fixed ratio of horizontal to vertical stresses components; and (2) increasing the differential stresses (i.e., the difference between the horizontal and vertical stresses) while keeping the magnitude of vertical stress constant. These numerical experiments show that the permeability of fractured rocks decreases with increased stress magnitudes when the stress ratio is not large enough to cause shear dilation of fractures, whereas permeability increases with increased stress when the stress ratio is large enough. Permeability changes at low stress levels are more sensitive than at high stress levels due to the nonlinear fracture normal stress-displacement relation. Significant stress-induced channeling is observed as the shear dilation causes the concentration of fluid flow along connected shear fractures. Anisotropy of permeability emerges with the increase of differential stresses, and this anisotropy can become more prominent with the influence of shear dilation and localized flow paths. A set of empirical equations in closed-form, accounting for both normal closure and shear dilation of the fractures, is proposed to model the stress-dependent permeability. These equations prove to be in good agreement with the results obtained from our numerical experiments.