Hollow Cylinder Tests on Boom Clay: Modelling of Strain Localization in the Anisotropic Excavation Damaged Zone

Boom Clay is extensively studied as a potential candidate to host underground nuclear waste disposal in Belgium. To guarantee the safety of such a disposal, the mechanical behaviour of the clay during gallery excavation must be properly predicted. In that purpose, a hollow cylinder experiment on Boom Clay has been designed to reproduce, in a small-scale test, the Excavation Damaged Zone (EDZ) as experienced during the excavation of a disposal gallery in the underground. In this article, the focus is made on the hydro-mechanical constitutive interpretation of the displacement (experimentally obtained by medium resolution X-ray tomography scanning). The coupled hydro-mechanical response of Boom Clay in this experiment is addressed through finite element computations with a constitutive model including strain hardening/softening, elastic and plastic cross-anisotropy and a regularization method for the modelling of strain localization processes. The obtained results evidence the directional dependency of the mechanical response of the clay. The softening behaviour induces transient strain localization processes, addressed through a hydro-mechanical second grade model. The shape of the obtained damaged zone is clearly affected by the anisotropy of the materials, evidencing an eye-shaped EDZ. The modelling results agree with experiments not only qualitatively (in terms of the shape of the induced damaged zone), but also quantitatively (for the obtained displacement in three particular radial directions).     

Analysis of Fracture Mechanics Tests on Opalinus Clay

Many studies have recently been conducted to evaluate various mechanical characteristics of the Opalinus Clay (OPA) formation in view of its potential use as the hosting rock for the Swiss nuclear waste repositories. Its sedimentary bedding makes OPA a transversely isotropic rock and its directional mechanical properties need to be measured. This paper reports on an experimental and computational approach that was adopted to define the parallel-to-bedding fracture mechanics (FM) parameters of OPA in Mode-I. OPA cores from Mont Terri Underground Research Laboratory (URL) were submitted to laboratory tests on notched semi-circular specimens under three-point bending (SCB). In these tests, crack propagation is forced along the notch direction. However, the 45° bedding inclination of the specimen axis frequently deviated the crack from the expected direction. An analysis of the SCB tests was performed by means of non-linear FM techniques and the pertinent Mode-I parameters along the bedding were estimated.     

Estimating Abrasivity of Rock by Laboratory and In Situ Tests

The degree to which a rock abrades another rock is called its “abrasivity”. Laboratory tests of abrasivity can be broadly divided into four kinds: drilling, rubbing, turning-operation and tumbling tests. The present study was initiated 30 years ago with the objective of investigating and developing methods for measuring rock abrasivity, and making some contribution towards understanding the relationships between the above test methods. Within the range of tests conducted, the turning-operation test turned out to be superior to the drilling test, albeit slightly, in terms of practicality. We have also conducted in situ tests using rock drills for the last 20 years. The results of those tests have been investigated and compared with the results of laboratory tests. There is a large degree of scatter in the data on gauge loss in button bits, which has obscured any correlations with laboratory data. Some correlations were found between height loss in button bits and laboratory findings.     

Hollow Cylinder Simulation Experiments of Galleries in Boom Clay Formation

In the context of nuclear waste disposal in clay formations, laboratory experiments were performed to study at reduced scale the excavation damaged zone (EDZ) induced by the construction of galleries in the Boom clay formation. For this purpose, thick-walled hollow cylindrical samples were subjected (after recovery of in situ stress conditions) to a decrease in the inner confining pressure aiming at mimicking a gallery excavation. X-ray computed tomography (XRCT) scans of the specimens were carried out through the testing cell before and after the mechanical unloading and allowed to quantify the displacements undergone by the clay as a result of the mechanical unloading. The deformation of the hollow cylinders and the inferred extent of the damaged zone around the central hole are found to depend on the orientation of the specimen with respect to the bedding planes and show a great similarity with in situ observations around galleries and boreholes at Mol URL in the Boom clay formation. In the experiments performed on samples cored parallel to the bedding, the damaged zone is not symmetrical with respect to the hole axis and extends more in the direction parallel to the bedding. It is the same for the radial convergence of the hole walls which is larger in the direction parallel to bedding than in the perpendicular one. In contrast, a test on a sample cored perpendicularly to the bedding did not show any ovalisation of the central hole after the mechanical unloading. These observations confirm the significance of the pre-existing planes of weakness (bedding planes) in Boom clay and the need for a correct consideration of the related mechanical anisotropy.     

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.     

Temperature and Damage Impact on the Permeability of Opalinus Clay

The effects of temperature on the water transport properties of intact and damaged Opalinus (OPA) clay are investigated by using a recently developed hollow cylinder triaxial cell [Monfared (Int J Rock Mech Min Sci 48:637–649, 2011b)] that allows full saturation and drainage conditions in low-permeability clays and shales. The volumetric response of saturated OPA clay sample during a drained heating test shows an irreversible contraction after a temperature threshold. The permeability tests which are performed before and after the heating test show that the induced irreversible sample contraction by thermal loading reduces the permeability of OPA clay sample. In order to study the effect of temperature on the permeability of a damaged sample of OPA clay, the permeability tests are performed on a saturated sample previously sheared by a standard drained triaxial loading. The test results show no significant effect of shear-type damage on the permeability of the sample at 25 °C and 80 °C. The experimental results presented in this paper show the crucial role of the thermally induced strains on sample permeability. Thermo-elastic dilation leads to a slight increase of the permeability, whereas thermoplastic contraction leads to a reduction.     

In Situ Agrobiodiversity Conservation in the Swiss Inner Alpine Zone

This paper considers Swiss experiences of in situ agrobiodiversity conservation. The Swiss recognise the multifunctionality of agriculture, particularly within the marginal, mountainous areas of the country and are spending public and private funds to maintain crop landraces, agroecosystems, landscapes, agricultural practices, and rural communities. State programs in marginal agricultural regions tie economic assistance to biodiversity conservation practices, through the use of direct payment systems. With the removal of many import restrictions and state-led marketing, alternative innovative programs linked to local values of diversity provide opportunities for conservation. In situ agrobiodiversity conservation ideas being pursued include awareness-raising schemes, scientific initiatives, the establishment of diversity gardens and community-based programs. Any extension of the in situ conservation programs could have substantial implications for regional development. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

Consequences of the Thermal Transient on the Evolution of the Damaged Zone Around a Repository for Heat-Emitting High-Level Radioactive Waste in a Clay Formation: a Performance Assessment Perspective

A proper evaluation of the perturbations of the host rock induced by the excavation and the emplacement of exothermic wastes is essential for the assessment of the long-term safety of high-level radioactive waste disposals in clay formations. The impact of the thermal transient on the evolution of the damaged zone (DZ) has been explored in the European Commission project TIMODAZ (thermal impact on the damaged zone around a radioactive waste disposal in clay host rocks, 2006–2010). This paper integrates the scientific results of the TIMODAZ project from a performance assessment (PA) point of view, showing how these results support and justify key PA assumptions and the values of PA model parameters. This paper also contextualises the significance of the thermal impact on the DZ from a safety case perspective, highlighting how the project outcomes result into an improved understanding of the thermo–hydro–mechanical behaviour of the clay host rocks. The results obtained in the TIMODAZ project strengthen the assessment basis of the safety evaluation of the current repository designs. There was no evidence throughout the TIMODAZ experimental observations of a temperature-induced additional opening of fractures nor of a significant permeability increase of the DZ. Instead, thermally induced plasticity, swelling and creep seem to be beneficial to the sealing of fractures and to the recovery of a very low permeability in the DZ, close to that of an undisturbed clay host rock. Results from the TIMODAZ project indicate that the favourable properties of the clay host rock, which guarantee the effectiveness of the safety functions of the repository system, are expected to be maintained after the heating–cooling cycle. Hence, the basic assumptions usually made in PA calculations so far are expected to remain valid, and the performance of the system should not be affected in a negative way by the thermal evolution of the DZ around a radioactive waste repository in clay host rock.     

Modeling the Effects of Heterogeneity and Anisotropy on the Excavation Damaged/Disturbed Zone (EDZ)

Summary  When modeling the mechanical behavior of underground excavations, it is necessary to include the influence of the rock mass characteristics on the Excavation Damaged/Disturbed Zone (EDZ). In this paper, the Realistic Failure Process Analysis code, RFPA, is used to model the extent of the EDZ. The inhomogeneous characteristics of rock at the mesoscopic level are included by assuming that the material properties of the constituent elements conform to a Weibull distribution; the anisotropy is incorporated as a transversely isotropic medium; the non-elastic characteristic is simulated via an elastic damage-based constitutive law. A finite element program is adopted as the basic stress analysis tool. In this study, a notable feature is that no a priori assumptions need to be made about where and how fracture and failure will occur – cracking can take place spontaneously and can exhibit a variety of mechanisms when certain local stress conditions are met. The deformation and failure process of anisotropic rock around excavations of different geometries is analyzed, and compared to experimental tests, showing similar fracture patterns. Additionally, the effect of confining stress and of different material layers is modeled and discussed. It is found that the model clearly illustrates that fracturing, both initiation and propagation, occurs as a combination of the stress concentrations and weakness planes introduced via the transverse anisotropy – which could represent either foliations or ubiquitous joint sets. Correspondence: Dr. Shuhong Wang, Box 265, School of Resource and Civil Engineering, Northeastern University, Shenyang 110004, P.R. China     

Rates of microbial hydrogen oxidation and sulfate reduction in Opalinus Clay rock

Hydrogen gas (H₂) may be produced by the anoxic corrosion of steel components in underground structures, such as geological repositories for radioactive waste. In such environments, hydrogen was shown to serve as an electron donor for autotrophic bacteria. High gas overpressures are to be avoided in radioactive waste repositories and, thus, microbial consumption of H₂ is generally viewed as beneficial. However, to fully consider this biological process in models of repository evolution over time, it is crucial to determine the in situ rates of microbial hydrogen oxidation and sulfate reduction. These rates were estimated through two distinct in situ experiments, using several measurement and calculation methods. Volumetric consumption rates were calculated to be between 1.13 and 1.93 μmol cm⁻³ day⁻¹ for H₂, and 0.14 and 0.20 μmol cm⁻³ day⁻¹ for sulfate. Based on the stoichiometry of the reaction, there is an excess of H₂ consumed, suggesting that it serves as an electron donor to reduce electron acceptors other than sulfate, and/or that some H₂ is lost via diffusion. These rate estimates are critical to evaluate whether biological H₂ consumption can negate H₂ production in repositories, and to determine whether sulfate reduction can consume sulfate faster than it is replenished by diffusion, which could lead to methanogenic conditions.     

Drained triaxial testing of shales: insight from the Opalinus Clay

Acta Geotechnica - The investigation of the mechanical behaviour of swelling claystones and shales is challenging because of their very low permeability and of their high sensitivity to changes in...     

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.     

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.     

Characterizing Excavation Damaged Zone and Stability of Pressurized Lined Rock Caverns for Underground Compressed Air Energy Storage

In this paper, we investigate the influence of the excavation damaged zone (EDZ) on the geomechanical performance of compressed air energy storage (CAES) in lined rock caverns. We conducted a detailed characterization of the EDZ in rock caverns that have been excavated for a Korean pilot test program on CAES in (concrete) lined rock caverns at shallow depth. The EDZ was characterized by measurements of P- and S-wave velocities and permeability across the EDZ and into undisturbed host rock. Moreover, we constructed an in situ concrete lining model and conducted permeability measurements in boreholes penetrating the concrete, through the EDZ and into the undisturbed host rock. Using the site-specific conditions and the results of the EDZ characterization, we carried out a model simulation to investigate the influence of the EDZ on the CAES performance, in particular related to geomechanical responses and stability. We used a modeling approach including coupled thermodynamic multiphase flow and geomechanics, which was proven to be useful in previous generic CAES studies. Our modeling results showed that the potential for inducing tensile fractures and air leakage through the concrete lining could be substantially reduced if the EDZ around the cavern could be minimized. Moreover, the results showed that the most favorable design for reducing the potential for tensile failure in the lining would be a relatively compliant concrete lining with a tight inner seal, and a relatively stiff (uncompliant) host rock with a minimized EDZ. Because EDZ compliance depends on its compressibility (or modulus) and thickness, care should be taken during drill and blast operations to minimize the damage to the cavern walls.     

In Situ Pore Water Pressure Dissipation Testing of Marine Clay under Reclamation Fills

In situ dissipation tests provide a means of evaluating the in situ coefficient of horizontal consolidation and horizontal hydraulic conductivity of soft clays. Dissipation tests by means of piezocone (CPTU), dilatometer (DMT), self-boring pressuremeter (SBPT) and BAT permeameter (BAT) were utilized in the characterization of the coefficient of horizontal consolidation and horizontal hydraulic conductivity of Singapore marine clay at Changi in a land reclamation project. Dissipation tests were carried out prior to reclamation as well as after ground improvement with vertical drains to compare the changes in the coefficient of horizontal consolidation and horizontal hydraulic conductivity prior to and after ground improvement. Tests were carried out in a vertical drain area as well as in an adjacent untreated control area after 23 months of surcharge loading, for comparison purposes. The purpose of this research is to determine the horizontal consolidation parameters of Singapore marine clay prior to reclamation as well as after 23 months of surcharge loading with and without vertical drains by means of in situ dissipation tests.     

Characterization of dissolved organic matter in anoxic rock extracts and in situ pore water of the Opalinus Clay

Dissolved organic matter (DOM) from the Opalinus Clay, a potential host rock for the disposal of radioactive waste, was isolated under strictly anoxic conditions from ground rock material and compared with DOM of in situ pore water samples. For the extractions, deionized water, synthetic pore water (SPW, water containing all major ions at pore water concentrations but no organic matter) and 0.1 M NaOH were used. The influence of the solid-to-liquid ratio, extraction time, acid-pretreatment and O₂ exposure of the rock material on the isolated DOM were investigated. Liquid chromatography coupled with a total organic C detector (LC-OCD) and reverse-phase ion chromatography were used to characterize the DOM size distributions and to determine the low molecular weight organic acid (LMWOA) contents in the pore water samples and the rock extracts.The results revealed that only a small portion of the total organic C of the rock material (<0.38%) was extractable, even after removal of carbonates by acid-pretreatment. The concentrations of dissolved organic C (DOC) were found to range from 3.9 ± 0.4 to 8.0 ± 0.8 mg/L in the anoxic extracts. The pore waters exhibited similar DOC concentrations ranging from 1.2 to 15.8 ± 0.5 mg/L. The analysis by LC-OCD showed that the DOM extracted under anoxic conditions and the pore water DOM mainly consisted of hydrophilic compounds of less than 500 Da. The DOM extracted with SPW was most similar in size to the pore water DOM. Grinding the rock under oxic conditions increased the DOC yields and shifted the size distribution toward higher molecular weight compounds compared to the strictly anoxic treatment. Acetate, lactate and formate were identified in all extracts and in the pore water. In total, LMWOA accounted for 36% of the total DOC in both pore water and SPW extracts. The results imply that controlled anoxic conditions and the use of SPW as an extractant are required to isolate DOM from Opalinus Clay rocks which most resembles the in situ pore water DOM with respect to its size distribution and the LMWOA contents.     

Thermal Conductivity of Argillaceous Rocks: Determination Methodology Using In Situ Heating Tests

This study focuses on the characterisation of thermal conductivity for three potential host rocks for radioactive waste disposal. First, the heat conduction process is reviewed on the basis of an analytical solution and key aspects related to anisotropic conduction are discussed. Then the existing information on the three rocks is summarised and a broad uncertainty range of thermal conductivity is estimated based on the mineralogical composition. Procedures to backanalyse the thermal conductivity on the basis of in situ heating tests are assessed and a methodology is put forward. Finally, this methodology is used to estimate the impact of experimental uncertainties and applied to the four in situ heating tests. In the three potential host rocks, a clear influence of the bedding planes was identified and anisotropic heat conduction was shown to be necessary to interpret the observed temperature field. Experimental uncertainties were also shown to induce a larger uncertainty on the anisotropy ratio than on the equivalent thermal conductivity defined as the geometric mean of the thermal conductivity in the three principal directions.