In situ behaviour of a stiff layered clay subject to thermal loading: Observations and interpretation

The paper presents an interpretation of an in situ heating test carried out on Opalinus clay in the Mont Terri underground laboratory. Opalinus clay is a stiff, strongly bedded, Mesozoic clay of marine origin. When subjected to thermal loading, saturated stiff clays exhibit a strong pore pressure response that significantly affects the hydraulic and mechanical behaviour of the material. The observations gathered in the in situ test have provided an opportunity to examine the integrated thermo-hydro-mechanical (THM) response of this sedimentary clay. Coupled THM numerical analyses have been carried out to provide a structured framework for interpretation, and to enhance understanding of THM clay behaviour. Numerical analyses have been based on a coupled theoretical formulation that incorporates a constitutive law especially developed for this type of material. The law includes degradation of bonding by damage. By performing three-dimensional computations, it has been possible to incorporate anisotropy of material parameters and of in situ stresses. The 3D simulation has proved able to furnish a satisfactory representation of the development of the in situ test and of the main observed patterns of behaviour. A sensitivity analysis has also been carried out to examine the potential effect of various key or uncertain parameters. The critical examination of test observations and the results of the numerical analyses have allowed the classification, by differing degrees of significance, of the various coupled phenomena present in the proble

RHEOLOGICAL CHARACTERIZATION OF A CLAY FORMATION FROM DRIFTS EXCAVATION : ELASTIC AND ELASTOPLASTIC APPROACH

International audienceAn extensive scientific programme has been carried out by Andra (French Agency in charge of radioactive waste management) for investigating feasibility of High Level Activity Waste disposal in deep geological formation. An Underground Research Laboratory (URL) is currently being constructed in North-eastern France to assess the adequacy of a hard-clay argillite layer (Callovo-Oxfordian formation) situated between 420 m and 550 m of depth. Geotechnical measurements have been carried out during the shafts and drifts excavation and particularly upon the main level of the laboratory (-490 m). The drifts are “horseshoe section” type with about 17 m² in area mainly supported by metallic ribs and rock bolts. The digging has been performed with classical pneumatic hammer. Measurement sections have been instrumented very close to the front face using convergencemeters and radial extensometers. This paper presents a comparison between in situ measurements and numerical modelling. Elastic calculations are not in agreement with the measured deformations. An elastoplastic constitutive model considering damage and using Hoek & Brown criteria has been developed and implemented in the FLAC3D numerical code. Mechanical parameters came from lab tests performed on core samples. For the first meters, model provides consistent displacements. Beyond 4 meters, a time dependent convergence takes place and has to be integrated in the model to take into account creep and/or hydromechanical behaviour

Hydromechanical response to a mine by test experiment in a deep claystone

International audienceIn order to demonstrate the feasibility of radioactive waste repository in deep geological formation, an underground research laboratory is being constructed by Andra (French national radioactive waste management agency) in eastern France, in a Callovo-Oxfordian claystone. 15 boreholes were drilled from a drift at -447 m to install sensors around the shaft (6 m diameter) at depth -460 m to -474 m in order to record the hydromechanical behaviour of the claystone during a shaft sinking (drill and blast method). This paper is devoted to the analysis of the mechanical and hydromechanical behaviour observed during the shaft sinking. Analytical approach used herein allows to realistically evaluate the undrained response of the shaft neighboring with agreement with the in situ measurements. For the transient phase, prediction is qualitatively in accordance with measurement. In addition, deformation and displacement measurements are successfully compared to a simple 3D elastic calculation performed with the real face advancing. This emphasizes the quality of the data set which would be compared in the Modex-Rep European project with complex numerical modelling (poro-elasto-plasic-damage models, creep behaviour,...).En France, l'Andra (Agence Nationale de gestion des déchets radioactifs) est en charge des études pour la faisabilité d'un stockage de déchet radioactif haute activité à vie longue, dans une formation géologique profonde comme des argiles. Pour cela, l'Andra a construit le laboratoire de recherche souterrain de Meuse Haute Marne dans le Nord-Est de la France (à 300 km environs de Paris) dans une formation d'argilite du Callovo-Oxfordien qui se trouve entre 420 m et 550 m au niveau du puits principal. La première expérimentation géomécanique réalisée dans le laboratoire est un " mine by test " autour du puits principal (diamètre 6 m). A partir d'une galerie se trouvant à -445 m, 15 forages ont été réalisés pour installer des capteurs et mesurer le comportement hydromécanique de l'argilite entre -460 et -470 m lors du creusement du puits principal. Les évolutions des mesures in situ sont présentées et comparées avec la solution analytique poro-élastique du creusement d'un puits infini. Les mesures de déplacement sont comparées avec les résultats d'un calcul élastique en 3D. Ces analyses simples montrent la cohérence et la qualité des différentes mesures in situ qui serviront de données de référence dans le projet Européen ModexRep, où le creusement du puits est simulé avec des modèles 3D complexes (poro-elasto-endo-plasticité, fluage,...)

In situ behaviour of a stiff layered clay subject to thermal loading: Observations and interpretation

George Stephenson Medal 2008, atorgada per la Institution of Civil Engineers del Regne UnitThe paper presents an interpretation of an in situ heating test carried out on Opalinus clay in the Mont Terri underground laboratory. Opalinus clay is a stiff, strongly bedded, Mesozoic clay of marine origin. When subjected to thermal loading, saturated stiff clays exhibit a strong pore pressure response that significantly affects the hydraulic and mechanical behaviour of the material. The observations gathered in the in situ test have provided an opportunity to examine the integrated thermo-hydro-mechanical (THM) response of this sedimentary clay. Coupled THM numerical analyses have been carried out to provide a structured framework for interpretation, and to enhance understanding of THM clay behaviour. Numerical analyses have been based on a coupled theoretical formulation that incorporates a constitutive law especially developed for this type of material. The law includes degradation of bonding by damage. By performing three-dimensional computations, it has been possible to incorporate anisotropy of material parameters and of in situ stresses. The 3D simulation has proved able to furnish a satisfactory representation of the development of the in situ test and of the main observed patterns of behaviour. A sensitivity analysis has also been carried out to examine the potential effect of various key or uncertain parameters. The critical examination of test observations and the results of the numerical analyses have allowed the classification, by differing degrees of significance, of the various coupled phenomena present in the problem.Peer ReviewedAward-winningPostprint (published version

Modeling the growth of stylolites in sedimentary rocks

[1] Stylolites are ubiquitous pressure solution seams found in sedimentary rocks. Their morphology is shown to follow two self-affine regimes. Analyzing the scaling properties of their height over their average direction shows that (1) at small scale, they are self-affine surfaces with a Hurst exponent around 1, and (2) at large scale, they follow another self-affine scaling with Hurst exponent around 0.5. In the present paper, we show theoretically the influence of the main principal stress and the local geometry of the stylolitic interface on the dissolution reaction rate. We compute how it is affected by the deviation between the principal stress axis and the local interface between the rock and the soft material in the stylolite. The free energy entering in the dissolution reaction kinetics is expressed from the surface energy term and via integration from the stress perturbations due to these local misalignments. The resulting model shows the interface evolution at different stress conditions. In the stylolitic case, i.e., when the main principal stress is normal to the interface, two different stabilizing terms dominate at small and large scales which are linked respectively to the surface energy and to the elastic interactions. Integrating the presence of small-scale heterogeneities related to the rock properties of the grains in the model leads to the formulation of a Langevin equation predicting the dynamic evolution of the surface. This equation leads to saturated surfaces obeying the two observed scaling laws. Analytical and numerical analysis of this surface evolution model shows that the crossover length separating both scaling regimes depends directly on the applied far-field stress magnitude. This method gives the basis for the development of a paleostress magnitude marker. We apply the computation of this marker, i.e., the morphological analysis, on a stylolite found in the Dogger limestone layer located in the neighborhood of the ANDRA Underground Research Laboratory at Bure (eastern France). The results are consistent with the two scaling regimes expected, and the practical determination of the major principal paleostress, from the estimation of a crossover length, is illustrated on this example

Hydromechanical modelling of shaft sealing for CO2 storage

The geological sequestration of CO2 in abandoned coal mines is a promising option to mitigate climate changes while providing sustainable use of the underground cavities. In order to certify the efficiency of the storage, it is essential to understand the behaviour of the shaft sealing system. The paper presents a numerical analysis of CO2 transfer mechanisms through a mine shaft and its sealing system. Different mechanisms for CO2 leakage are considered, namely multiphase flow through the different materials and flow along the interfaces between the lining and the host rock. The study focuses on the abandoned coal mine of Anderlues, Belgium, which was used for seasonal storage of natural gas. A two-dimensional hydromechanical modelling of the storage site is performed and CO2 injection into the coal mine is simulated. Model predictions for a period of 500 years are presented and discussed with attention. The role and influence of the interface between the host rock and the concrete lining are examined. In addition the impact of some uncertain model parameters on the overall performance of the sealing system is analysed through a sensitivity analysis

Gas transfer through clay barriers

Gas transport through clay-rocks can occur by different processes that can be basically subdivided into pressure-driven flow of a bulk gas phase and transport of dissolved gas either by molecular diffusion or advective water flow (Figure 1, Marschall et al., 2005). The relative importance of these transport mechanisms depends on the boundary conditions and the scale of the system. Pressure-driven volume flow (“Darcy flow”) of gas is the most efficient transport mechanism. It requires, however, pressure gradients that are sufficiently large to overcome capillary forces in the typically water-saturated rocks (purely gas-saturated argillaceous rocks are not considered in the present context). These pressure gradients may form as a consequence of the gravity field (buoyancy, compaction) or by gas generation processes (thermogenic, microbial, radiolytic). Dissolved gas may be transported by water flow along a hydraulic gradient. This process is not affected by capillary forces but constrained by the solubility of the gas. It has much lower transport efficiency than bulk gas phase flow. Molecular diffusion of dissolved gas, finally, is occurring essentially without constraints, ubiquitously and perpetually. Effective diffusion distances are, however, proportional to the square root of time, which limits the relevance of this transport process to the range of tens to hundreds of metres on a geological time scale (millions of years). 2 Process understanding and the quantification of the controlling parameters, like diffusion coefficients, capillary gas breakthrough pressures and effective gas permeability coefficients, is of great importance for up-scaling purposes in different research disciplines and applications. During the past decades, gas migration through fully water-saturated geological clay-rich barriers has been investigated extensively (Thomas et al., 1968, Pusch and Forsberg, 1983; Horseman et al., 1999; Galle, 2000; Hildenbrand et al., 2002; Marschall et al., 2005; Davy et al., 2009; Harrington et al., 2009, 2012a, 2014). All of these studies aimed at the analysis of experimental data determined for different materials (rocks of different lithotype, composition, compaction state) and pressure/temperature conditions. The clay-rocks investigated in these studies, ranged from unconsolidated to indurated clays and shales, all characterised by small pores (2-100 nm) and very low hydraulic conductivity (K < 10-12 m·s-1) or permeability coefficients (k < 10-19 m²). Studies concerning radioactive waste disposal include investigations of both the natural host rock formation and synthetic/engineered backfill material at a depth of a few hundred meters (IAEA, 2003, 2009). Within a geological disposal facility, hydrogen is generated by anaerobic corrosion of metals and through radiolysis of water (Rodwell et al., 1999; Yu and Weetjens, 2009). Additionally, methane and carbon dioxide are generated by microbial degradation of organic wastes (Rodwell et al., 1999; Ortiz et al., 2002; Johnson, 2006; Yu and Weetjens, 2009). The focus of carbon capture and storage (CCS) studies is on the analysis of the long-term sealing efficiency of lithologies above depleted reservoirs or saline aquifers, typically at larger depths (hundreds to thousands of meters). During the last decade, several studies were published on the sealing integrity of clay-rocks to carbon dioxide (Hildenbrand et al., 2004; Li et al., 2005; Hangx et al., 2009; Harrington et al., 2009; Skurtveit et al., 2012; Amann-Hildenbrand et al., 2013). In the context of petroleum system analysis, a significant volume of research has been undertaken regarding gas/oil expulsion mechanisms from sources rocks during burial history (Tissot & Pellet, 1971; Appold & Nunn, 2002), secondary migration (Luo et al., 2008) and the capillary sealing capacity of caprocks overlying natural gas accumulations (Berg, 1975; Schowalter, 1979; Krooss, 1992; Schlömer and Kross, 2004; Li et al., 2005; Berne et al., 2010). Recently, more attention has been paid to investigations of the transport efficiency of shales in the context of oil/gas shale production (Bustin et al., 2008; Eseme et al., 2012; Amann-Hildenbrand et al., 2012; Ghanizadeh et al., 2013, 2014). Analysis of the migration mechanisms within partly unlithified strata becomes important when explaining the 3 origin of overpressure zones, sub-seafloor gas domes and gas seepages (Hovland & Judd, 1988; Boudreau, 2012). The conduction of experiments and data evaluation/interpretation requires a profound process understanding and a high level of experience. The acquisition and preparation of adequate samples for laboratory experiments usually constitutes a major challenge and may have serious impact on the representativeness of the experimental results. Information on the success/failure rate of the sample preparation procedure should therefore be provided. Sample specimens “surviving” this procedure are subjected to various experimental protocols to derive information on their gas transport properties. The present overview first presents the theoretical background of gas diffusion and advective flow, each followed by a literature review (sections 2 and 3). Different experimental methods are described in sections 4.1 and 4.2. Details are provided on selected experiments performed at the Belgian Nuclear Research Centre (SCK-CEN, Belgium), Ecole Centrale de Lille (France), British Geological Survey (UK), and at RWTH-Aachen University (Germany) (section 4.3). Experimental data are discussed with respect to different petrophysical parameters outlined above: i) gas diffusion, ii) evolution of gas breakthrough, iii) dilation-controlled flow, and iv) effective gas permeability after breakthrough. These experiments were conducted under different pressure and temperature conditions, depending on sample type, burial depth and research focus (e.g. radioactive waste disposal, natural gas exploration, or carbon dioxide storage). The interpretation of the experimental results can be difficult and sometimes a clear discrimination between different mechanisms (and the controlling parameters) is not possible. This holds, for instance, for gas breakthrough experiments where the observed transport can be interpreted as intermittent, continuous, capillary- or dilation-controlled flow. Also, low gas flow rates through samples on the length-scale of centimetres can be equally explained by effective two-phase flow or diffusion of dissolved gas

Modelling of the hydromechanical response of a shaft sinking in a deep claystone

International audienceIn the framework of the feasibility study of a radioactive-waste repository in a geological formation, Andra (French radioactive waste management agency) has been built an underground research laboratory within a Callovo-Oxfordian argillite formation located in Eastern France. During the sinking of the laboratory's access shaft, the hydromechanical behaviour of the argillites was monitored through an in situ experimental program called REP. The experimental zone is located between 460 and 476m depth. From a drift located at -445 m, 15 instrumentated boreholes were drilled downwards and equipped with 120 mechanical and hydraulic sensors. A predictive modelling of the shaft sinking has been performed in the framework of the European Modex-Rep project, using a poro-elastoplastic model based on a generalized Hoek and Brown criterion. Results of blind prediction emphasize that the model reproduces the in situ phenomena, but is not able to reproduce the amplitude of the drop in pore pressure during the shaft sinking. This article presents the analysis of this discrepancy and discusses new approaches aimed at improving the previous model. Finally, permeability changes around the shaft based on the in situ measurements and/or damage-induced permeability changes implemented in the model allow to better reproduce in situ data