Large-scale tracer profiles in a deep claystone formation (Opalinus Clay at Mont Russelin,Switzerland): Implications for solute transport processes and transport properties of the rock

Natural tracers (Cl^? and stable water isotopes) in pore water of the Opalinus Clay and adjacent formations were studied in the motorway tunnel at Mont Russelin, Switzerland. The Opalinus Clay occurs in the core of an anticline which is cut by a complex system of thrust faults. Concentration profiles of natural tracers were taken from 17 boreholes along a 363 m long section. Pore waters of drillcore samples were analysed with indirect and direct methods. The Cl^? and stable water isotope distribution in the pore water shows a regular and well defined profile, with a conspicuous decrease towards the overlying Dogger limestone aquifer. The highest Cl^? values (approximately 23,000 mg/L) are found in the core of the anticline in Liassic claystones underlying the Opalinus Clay. To quantify the large-scale transport properties of the Opalinus Clay formation, a 2D transport model was constructed and used to reproduce the observed concentration profiles. The calculations indicate that the observed tracer distributions are consistent with diffusion as the dominant transport process. Groundwater flow in the overlying Dogger aquifer was initiated about 2–4 Ma ago, which is long after the folding of the Jura Mountains and probably coincides with the exposure of the aquifer to freshwater recharge following continued erosion of the anticline. The calculations suggest that tracer distributions are controlled by 1) the timing of freshwater recharge in the overlying limestone aquifer, 2) the shape of the anticline and 3) the magnitude and the anisotropy ratio of diffusion coefficients.     

Modelling the Mont Terri HE-D experiment for the Thermal–Hydraulic–Mechanical response of a bedded argillaceous formation to heating

Coupled thermal–hydrological–mechanical (THM) processes in the near field of deep geological repositories can influence several safety features of the engineered and geological barriers. Among those features are: the possibility of damage in the host rock, the time for re-saturation of the bentonite, and the perturbations in the hydraulic regime in both the rock and engineered seals. Within the international cooperative code-validation project DECOVALEX-2015, eight research teams developed models to simulate an in situ heater experiment, called HE-D, in Opalinus Clay at the Mont Terri Underground Research Laboratory in Switzerland. The models were developed from the theory of poroelasticity in order to simulate the coupled THM processes that prevailed during the experiment and thereby to characterize the in situ THM properties of Opalinus Clay. The modelling results for the evolution of temperature, pore water pressure, and deformation at different points are consistent among the research teams and compare favourably with the experimental data in terms of trends and absolute values. The models were able to reproduce the main physical processes of the experiment. In particular, most teams simulated temperature and thermally induced pore water pressure well, including spatial variations caused by inherent anisotropy due to bedding.     

A micro–macro approach of permeability evolution in rocks excavation damaged zones

Excavation damaged zone, with significant irreversible deformations and nonnegligible changes in flow and transport properties generally occurs in indurated clay around underground structures. The stress perturbation around the excavation could lead to a significant increase of the permeability physically due to diffuse and/or localized microcracks growth in the material. In the present study, we investigate microcracks-induced damage processes together with the subsequent modification in permeability. The proposed approach is based on a homogenization-based upper bound extended to the context of microcracked media in presence of initial stress. Application of this approach is done on a borehole excavation problem related to the Selfrac in situ experiments on Opalinus Clay. Although, the model fails to quantitatively account for the in situ permeability change (which may also originated from existing macrofractures), its prediction shows a significant evolution of the material permeability around the borehole. This is in qualitative agreement with available data.     

The use of two-part Hooke’s model (TPHM) to model the mine-by test at Mont Terri Site,Switzerland

The full-scale mine-by (MB) test conducted in 2008 in the Mont Terri underground rock laboratory (Switzerland) investigated the deformation and the coupled hydro-mechanical behavior of the Opalinus Clay in response to tunnel excavations. The Opalinus Clay is currently under investigation in Switzerland as a potential host rock for geologic disposal of high-level radioactive waste. To further improve the understanding and modeling of the coupled processes and their impact on the performance of a geologic repository in Opalinus Clay, a newly developed two-part Hooke’s model (TPHM) was implemented into a geomechanical simulator. A three-dimensional simulation model based on the TPHM was then developed to predict the deformation and pore pressure responses in the near field of the MB Niche 2 test at the Mont Terri Site. The usefulness and validity of the TPHM are demonstrated by the consistency between simulation results and field observations. Simulation results show that the pore pressure disturbance becomes visible at about 11 m ahead of the mine-by excavation advancing face (along the longitudinal direction of the MB Niche). The results also demonstrate that there exists a good correlation between the excavation damage zone (EDZ) and the pore-pressure evolution, which may have important practical implications for monitoring EDZ evolution with pore-pressure sensors. The simulation results, which are sensitive to the constitutive relationships used in the model, capture both the observed displacements and the size of the damage zone, whereas the approach based on the conventional Hooke’s law underestimates both. The comparison between simulated and observed results also indicates that laboratory-measured mechanical properties can be used to accurately predict field-scale mechanical deformations, as long as valid constitutive relationships are employed.     

Laboratory and In Situ Simulation Tests of the Excavation Damaged Zone Around Galleries in Opalinus Clay

In the context of nuclear waste disposal in clay formations, laboratory and in situ simulation experiments were performed to study at reduced scale the excavation damaged zone (EDZ) around tunnels in the indurated Opalinus Clay at Mont Terri, Switzerland. In the laboratory, thick-walled hollow cylindrical specimens were subjected to a mechanical unloading mimicking a gallery excavation. In samples cored parallel to bedding, cracks sub-parallel to the bedding planes open and lead to a buckling failure in two regions that extend from the borehole in the direction normal to bedding. The behaviour is clearly anisotropic. On the other hand, in experiments performed on specimens cored perpendicular to bedding, there is no indication of failure around the hole and the response of the hollow cylinder sample is mainly isotropic. The in situ experiment at Mont Terri which consisted in the overcoring of a resin-injected borehole that follows the bedding strike of the Opalinus Clay showed a striking similarity between the induced damaged zone and the fracture pattern observed in the hollow cylinder tests on samples cored parallel to bedding and such a bedding controlled “Excavation” Damaged Zone is as well consistent with the distinct fracture patterns observed at Mont Terri depending on the orientation of holes/galleries with respect to the bedding planes. Interestingly, the damaged zone observed in the hollow cylinder tests on samples cored parallel to bedding and in situ around URL galleries is found to develop in reverse directions in Boom Clay (Mol) and in Opalinus Clay (Mont Terri). This most probably results from different failure mechanisms, i.e. shear failure along conjugated planes in the plastic Boom Clay, but bedding plane splitting and buckling in the indurated Opalinus Clay.     

Hydrogeology of a fractured shale (Opalinus Clay): Implications for deep geological disposal of radioactive wastes

As part of the Swiss programme for high-level radioactive-waste disposal, a Jurassic shale (Opalinus Clay) is being investigated as a potential host rock. Observations in clay pits and the results of a German research programme focusing on hazardous waste disposal have demonstrated that, at depths of 10–30 m, the permeability of the Opalinus Clay decreases by several orders of magnitude. Hydraulic tests in deeper boreholes (test intervals below 300 m) yielded hydraulic conductivities <10^–12 m/s, even though joints and faults were included in some of the test intervals. These measurements are consistent with hydrogeological data from Opalinus Clay sections in ten tunnels in the Folded Jura of northern Switzerland. Despite extensive faulting, only a few indications of minor water inflow were encountered in more than 6,600 m of tunnel. All inflows were in tunnel sections where the overburden is less than 200 m. The hydraulic data are consistent with clay pore-water hydrochemical and isotopic data. The extensive hydrogeological data base – part of which derives from particularly unfavourable geological environments – provides arguments that advective transport through faults and joints is not a critical issue for the suitability of Opalinus Clay as a host rock for deep geological waste disposal. Electronic Publication     

Biogeochemical processes in a clay formation in situ experiment: Part E - Equilibrium controls on chemistry of pore water from the Opalinus Clay, Mont Terri Underground Research Laboratory, Switzerland

The chemistry of pore water (particularly pH and ionic strength) is an important property of clay rocks being considered as host rocks for long-term storage of radioactive waste. Pore waters in clay-rich rocks generally cannot be sampled directly. Instead, their chemistry must be found using laboratory-measured properties of core samples and geochemical modelling. Many such measurements have been made on samples from the Opalinus Clay from the Mont Terri Underground Research Laboratory (URL). Several boreholes in that URL yielded water samples against which pore water models have been calibrated. Following a first synthesis report published in 2003, this paper presents the evolution of the modelling approaches developed within Mont Terri URL scientific programs through the last decade (1997-2009). Models are compared to the composition of waters sampled during dedicated borehole experiments. Reanalysis of the models, parameters and database enabled the principal shortcomings of the previous modelling efforts to be overcome. The inability to model the K concentrations correctly with the measured cation exchange properties was found to be due to the use of an inappropriate selectivity coefficient for Na-K exchange; the inability to reproduce the measured carbonate chemistry and pH of the pore waters using mineral-water reactions alone was corrected by considering clay mineral equilibria. Re-examination of the measured Ca/Mg activity ratios and consideration of the mineralogical composition of the Opalinus Clay suggested that Ca/Mg cation exchange rather than dolomite saturation may control the ratio of these ions in solution. This re-examination also suggests that the Ca/Mg ratio decreases with increasing pore-water salinity. Several possible reasons for this are proposed. Moreover, it is demonstrated that feldspar equilibria must not be included in Opalinus Clay modelling because feldspars are present only in very small quantities in the formation and because Na/K ratios measured in pore water samples are inconsistent with feldspar saturation. The principal need to improve future modelling is additional or better data on rock properties, in particular: (i) a more detailed identification of phases in the Opalinus Clay that include redox-sensitive elements together with evaluation of their thermodynamic properties; (ii) an improved understanding of the distribution of celestite throughout the Opalinus Clay for Sr/SO₄ concentrations control; (iii) improvements in analytic and thermodynamic data for Ca-Mg rock cation exchange and mineral chemical properties and (iv) the measurement of composition and stability constants of clay minerals actually present in the formation.     

Characterization of drying-induced deformation behaviour of Opalinus Clay and tuff in no-stress regime

The drying-induced deformation behaviour of Opalinus Clay and tuff, which are being investigated under international and local collaborative projects for nuclear waste disposal in Switzerland and Japan, was investigated under a no-stress condition in the laboratory to evaluate their generic susceptibility to the formation of excavation damaged zone. The cylindrical core samples of Opalinus Clay and tuff were prepared to a one-dimensional drying condition and submitted to an uncontrolled laboratory environment. The strain evolution, evaporative water loss and environmental entities, such as temperature and relative humidity, were recorded simultaneously and quasi-continuously. It was observed that the drying phase induced significant strain magnitude and damage in Opalinus Clay samples, which was evidenced by the formation of hairy cracks on the surface parallel to the bedding. On the contrary, the strain occurrences in tuff samples were relatively insignificant, and no tendency of cracking was observed. In addition, the quasi-continuous availability of volumetric strains was further used in poroelastic relation for the estimation of capillary suction evolution. The calculated results were validated with pore size distributions obtained from mercury intrusion porosimetry.     

Thermal Impact on Damaged Boom Clay and Opalinus Clay: Permeameter and Isostatic Tests with μCT Scanning

Within the framework of the TIMODAZ project, permeameter tests and isostatic tests were performed on Boom Clay and Opalinus Clay in order to assess the impact of temperature, pore water composition, and confining stress on the sealing of damaged samples of Boom Clay and Opalinus Clay. A microfocus X-ray computed tomography technique was used to visualize the evolution of the sealing process. Compared to the fast sealing of Boom Clay, the sealing of Opalinus Clay was much slower. The heating showed a significant, favorable impact on the sealing behavior of Opalinus Clay under permeameter test conditions, while the sealing behavior of Boom Clay appeared to be unaffected. Tests performed under isostatic conditions did not reveal a significant influence of a heating–cooling cycle on the sealing behavior of these clays. The reappearance of the fractures or holes in the samples after dismantling confirms earlier observations which showed that after sealing, the original mechanical properties are not recovered. In other words, a heating cycle does not seem to induce healing.     

Constitutive analysis of the mechanical anisotropy of Opalinus Clay

This paper aims to analyse the anisotropic features of behaviour of Opalinus Clay using the theory of plastic multi-mechanisms. The results of triaxial tests conducted with different load levels and directions showed that the mechanical behaviour of this shale is cross-anisotropic. The stiffer samples are those in which the loading direction is parallel to the bedding plane. This indicates that the preconsolidation stress depends on the orientation of the load with respect to the fabric of Opalinus Clay. It is proposed to interpret the observed cross-anisotropy with an elastoplastic model based on four plastic strain mechanisms that may be successively mobilised depending on the loading direction. The predicted stress–strain responses vary according to the directions of the space as a result of the hardening process, depending on the number of plastic strain mechanisms that have been mobilised. The numerical predictions show overall good agreement with the experimental data in terms of deviatoric stress versus axial strain, demonstrating that multi-mechanism plasticity is a suitable constitutive tool for the interpretation of the mechanical anisotropy of this shale.     

Transport of U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of ²³⁸U, ²³⁴U and ²³⁰Th was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim was to assess the in situ, long-term migration behaviour of ²³⁴U in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both ²³⁸U and ²³⁰Th immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of ²³⁴U from the rock matrix to the pore water by the decay of ²³⁸U will be established. Diffusive redistribution separates ²³⁴U from its immobile parent ²³⁸U resulting in bulk rock ²³⁴U/²³⁸U activity disequilibria. These may provide a means of estimating the mobility of ²³⁴U in the rock if the diffusion rate of ²³⁴U is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by β-scanning, which shows that the (dominant) clay is richer in both elements.     

Numerical Modelling of the Anisotropic Mechanical Behaviour of Opalinus Clay at the Laboratory-Scale Using FEM/DEM

The Opalinus Clay (OPA) is an argillaceous rock formation selected to host a deep geologic repository for high-level nuclear waste in Switzerland. It has been shown that the excavation damaged zone (EDZ) in this formation is heavily affected by the anisotropic mechanical response of the material related to the presence of bedding planes. In this context, the purpose of this study is twofold: (i) to illustrate the new developments that have been introduced into the combined finite-discrete element method (FEM/DEM) to model layered materials and (ii) to demonstrate the effectiveness of this new modelling approach in simulating the short-term mechanical response of OPA at the laboratory-scale. A transversely isotropic elastic constitutive law is implemented to account for the anisotropic elastic modulus, while a procedure to incorporate a distribution of preferentially oriented defects is devised to capture the anisotropic strength. Laboratory results of indirect tensile tests and uniaxial compression tests are used to calibrate the numerical model. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in strong agreement with experimental observations. Subsequently, the calibrated model is validated by investigating the effect of confinement and the influence of the loading angle with respect to the specimen anisotropy. Simulated fracture patterns are discussed in the context of the theory of brittle rock failure and analyzed with reference to the EDZ formation mechanisms observed at the Mont Terri Underground Research Laboratory.     

The Mont Terri rock laboratory, a new international research project in a Mesozoic shale formation, in Switzerland

In many countries, argillaceous formations are being considered as potential host rocks for repositories of radioactive waste. Therefore, in 1995 several organisations decided to start an international research project in the reconnaissance gallery of the Mont Terri motorway tunnel, in north-western Switzerland, in a Mesozoic shale formation, the Opalinus Clay (Aalenian).

This project is under the patronage of the Swiss National Hydrological and Geological Survey. The following organisations are also partners in the project: ANDRA and IPSN (France), BGR (Germany), ENRESA (Spain), Nagra (Switzerland), JNC and Obayashi (Japan), SCKCEN (Belgium).

The aims of the project are to analyse the hydrogeological, geochemical and rock mechanical properties of an argillaceous formation, the changes of these properties induced by the excavation of galleries and to evaluate and improve appropriate investigation techniques. In January, 1996, eight niches were excavated for the project and 15 experiments were started in boreholes of no more than 30 m in length. In the winter of 1997–1998 a new research gallery was excavated to host further experiments. Different excavation, drilling and investigation techniques have been tested and improved. Significant results of the different experiments are already available. Some of those dealing with hydrochemistry, porewater pressure measurements and the evolution of the excavation disturbed zone (EDZ) are reported in this paper.     

Tunnelling and Hydrogeological Issues: A Short Review of the Current State of the Art

Since the 1960s, there has been an increasing interest in the understanding of the hydraulic flow inside a hard rock mass, since water inflow into deep tunnels constitute a hazard, in addition to being an important factor in controlling the advancement of excavation. The characterisation of fluid flow through hard rock masses is still one of the most challenging problems faced by geologists and engineers. A rock mass is characterised by networks of discrete and ubiquitous discontinuities that strongly affect its hydraulic properties, but detailed knowledge of the discontinuity properties allows for the evaluation of the hydraulic flow in the rock mass affected by the excavation of a tunnel. A geostructural field survey is fundamental in order to correctly define the discontinuity types, settings and networks. Numerous approaches have been proposed to estimate the water inflow based on empirical relations supported by field experience and case studies, as well as analytical solutions. Often, however, these approaches are not easily applicable in standard practice and in complex scenarios. The most appropriate approach to characterising the hydraulic flow of the rock mass and to predicting in the most effective way the expected water inflow during the excavation of a tunnel is based on a detailed geological model and geostructural analysis as described in this paper.     

A Laboratory Investigation on Thermal Properties of the Opalinus Claystone

Some aspects of the thermal behavior of the Opalinus claystone are investigated through laboratory tests conducted on a new hollow cylinder triaxial apparatus specially designed for studying the thermo-hydro-mechanical behavior of very low permeable materials. Two hollow cylinder samples are first resaturated under isotropic stress state equal to the mean effective in situ one in order to minimize swelling and induced damage during the resaturation phase. Two drained heating–cooling cycles are performed on the first sample of Opalinus claystone. During the first cycle, a thermo-elasto-plastic response similar to that of plastic clays with low overconsolidation ratio is obtained. The thermal hardening of the sample is demonstrated by the quasi-reversible behavior of the sample during the second heating–cooling cycle. An undrained heating test performed on the second sample of Opalinus claystone induces an excess pore pressure in this sample. This induced pore pressure is attributed to the higher thermal expansion coefficient of pore water compared to that of the solid phase. It is shown that the excess pore pressure generated in the sample by undrained heating cannot be modeled by considering the free water thermal expansion coefficient. The thermal expansion coefficient of the Opalinus claystone water is back-analyzed from the experimental results which show a higher value than free water.     

Biomarker transformations as constraints for the depositional environment and for maximum temperatures during burial of Opalinus Clay and Posidonia Shale in northern Switzerland

Samples from two argillaceous formations (Opalinus Clay and Posidonia Shale) of near-identical maturity from northern Switzerland were subjected to a geochemical characterisation of organic matter and to confined-system pyrolysis experiments. Throughout the study area, the characteristics of organic matter are similar, indicating a spatially homogeneous sedimentary facies. Posidonia Shale contains marine organic matter deposited in a reducing environment, while a predominantly terrigenous source and a more oxidising environment of deposition was identified for Opalinus Clay. In the western and central parts of the study area, organic maturity is close to the onset of oil generation. In the easternmost part, a higher maturity has been reached due to a deeper burial below thick Tertiary Molasse deposits.Isothermal pyrolysis experiments were conducted at temperatures between 250 and 390 °C over 24 h. Bitumen yields increase along similar pathways for both Opalinus Clay and Posidonia Shale, but the maximum values are displaced by 10–20 °C. Data pertaining to maturity were determined from GC–MS analyses of saturated hydrocarbons, and specific attention was given to C₂₉-sterane and C₃₂-hopane isomerisation ratios. The evolution of these parameters with rising temperature is slightly different in the two formations, which is attributed to the contrasting organic facies. The pyrolysis data, together with literature data from natural basins, were used to calculate kinetic parameters for C₂₉-sterane and C₃₂-hopane, assuming a single-step isomerisation scheme according to the Arrhenius law. The resulting values based on pyrolysis data alone are very similar to those based on the combination of pyrolysis and natural data. Activation energies are similar in both formations, while the frequency factors are up to one order of magnitude higher for Posidonia Shale when compared to Opalinus Clay. For the Benken site, maximum temperature during Cretaceous burial was calculated on the basis of the kinetic data, using the TTI approach. The resulting temperatures of 75–80 °C are 5–10 °C below those derived in the literature from apatite fission-track analysis, vitrinite reflectance and basin modelling.     

Seismic and mechanical properties of Opalinus Clay: comparison between sandy and shaly facies from Mont Terri (Switzerland)

For the safe disposal of high-level radioactive waste, different host rocks are currently being considered. The favorable properties of clay are low permeability, some retention capacity concerning radionuclides, and the ability to self-seal cracks and fissures, e.g. by swelling or time-dependent compaction creep. In Switzerland, the Jurassic Opalinus Clay is envisaged as a potential host rock which—at Mont Terri—is subdivided into the sandy, shaly, and carbonate-rich facies, the latter being less abundant. For long-term safety assessments, the understanding of the relations of properties (e.g. mineralogical composition and microstructure) and performance (e.g. mechanical behavior) of clays and claystones is essential. In the case of the sandy Opalinus Clay, the mechanical strength increases with increasing carbonate content, because carbonates form the matrix. The mineralogical investigation of a set of sandy facies samples proved a significantly larger carbonate content (20–40 mass %) when compared to the shaly facies (10–20 mass %). The carbonates of the shaly Opalinus Clay, on the other hand, are mostly localized fossils aligned parallel to the bedding, acting as predetermined breaking points. Image analysis of SEM images of polished sections proved the determined microstructural differences. In addition, carbonate particles of the sandy facies are mostly isometric, whereas carbonate particles of the shaly facies cover a greater range of shapes. The mechanical tests were accompanied by investigations of the p- and s-wave velocities, which revealed that the anisotropy of the sandy facies is less pronounced than sedimentological analyses would suggest. The mechanical strength, which, for the first time, presents results of real triaxial tests of the sandy facies. The samples of the sandy facies exhibit a failure strength of σ _eff,B, approximately twice as high as was found for the shaly facies considering the deformation axis parallel to the bedding. Similar values were obtained when measuring perpendicularly to the bedding.