Effects of drilling and stress release on hydraulic properties and porewater chemistry of crystalline rocks

The experimental verification of matrix diffusion in crystalline rocks largely relies on indirect methods performed in the laboratory. Such methods are prone to perturbations of the rock samples by collection and preparation and therefore the laboratory-derived transport properties and fluid composition might not represent in situ conditions. We investigated the effects induced by the drilling process and natural rock stress release by mass balance considerations and sensitivity analysis of analytical out-diffusion data obtained from originally saturated, large-sized drillcore material from two locations drilled using traced drilling fluid. For in situ stress-released drillcores of quartz-monzodiorite composition from the Aspo HRL, Sweden, tracer mass balance considerations and 1D and 2D diffusion modelling consistently indicated a contamination of <1% of the original pore water. This chemically disturbed zone extends to a maximum of 0.1 mm into the drillcore (61.8 mm x 180.1 mm) corresponding to about 0.66% of the total pore volume (0.77 vol.%). In contrast, the combined effects of stress release and the drilling process, which have influenced granodioritic drillcore material from 560 m below surface at Forsmark. Sweden, resulted in a maximum contamination of the derived porewater Cl(-) concentration of about 8%. The mechanically disturbed zone with modified diffusion properties covers the outermost similar to 6 mm of the drillcore (50 mm x 189 mm), whereas the chemically disturbed zone extends to a maximum of 0.3 mm based on mass balance considerations, and to 0.15 mm to 0.2 mm into the drillcore based on fitting the observed tracer data. This corresponds to a maximum of 2.4% of the total pore volume (0.62 vol.%) being affected by the drilling-fluid contamination. The proportion of rock volume affected initially by drilling fluid or subsequently with experiment water during the laboratory diffusion and re-saturation experiments depends on the size of the drillcore material and will become larger the smaller the sample used for the experiment. The results are further in support of matrix diffusion taking place in the undisturbed matrix of crystalline rocks at least in the cm range.

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Matrix Porewater In Crystalline Rocks: Extraction and Analysis

The chemical and isotopic characterization of porewater residing in the inter- and intragranular pore space of the low-permeability rock matrix is an important component with respect to the site characterization and safety assessment of potential host rocks for a radioactive waste disposal. The chemical and isotopic composition of porewater in such low permeability rocks has to be derived by indirect extraction techniques applied to naturally saturated rock material. In most of such indirect extraction techniques – especially in case of rocks of a porosity below about 2 vol.% – the original porewater concentrations are diluted and need to be back-calculated to in-situ concentrations. This requires a well-defined value for the connected porosity – accessible to different solutes under in-situ conditions. The derivation of such porosity values, as well as solute concentrations, is subject to various perturbations during drilling, core sampling, storage and experiments in the laboratory. The present study aims to demonstrate the feasibility of a variety of these techniques to charac-terize porewater and solute transport in crystalline rocks. The methods, which have been de-veloped during multiple porewater studies in crystalline environments, were applied on four core samples from the deep borehole DH-GAP04, drilled in the Kangerlussuaq area, Southwest Greenland, as part of the joint NWMO–Posiva–SKB Greenland Analogue Project (GAP). Potential artefacts that can influence the estimation of in situ porewater chemistry and isotopes, as well as their controls, are described in detail in this report, using specific examples from borehole DH-GAP04

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

Structural, metamorphic and geochronological relations between the Zanskar Shear Zone and the Miyar Shear Zone (NW Indian Himalaya): Evidence for two distinct tectonic structures and implications for the evolution of the High Himalayan Crystalline of Zanskar

The geologic structures and metamorphic zonation of the northwestern Indian Himalaya contrast significantly with those in the central and eastern parts of the range, where the high-grade metamorphic rocks of the High Himalayan Crystalline (HHC) thrust southward over the weakly metamorphosed sediments of the Lesser Himalaya along the Main Central Thrust (MCT). Indeed, the hanging wall of the MCT in the NW Himalaya mainly consists of the greenschist facies metasediments of the Chamba zone, whereas HHC high-grade rocks are exposed more internally in the range as a large-scale dome called the Gianbul dome. This Gianbul dome is bounded by two oppositely directed shear zones, the NE-dipping Zanskar Shear Zone (ZSZ) on the northern flank and the SW-dipping Miyar Shear Zone (MSZ) on the southern limb. Current models for the emplacement of the HHC in NW India as a dome structure differ mainly in terms of the roles played by both the ZSZ and the MSZ during the tectonothermal evolution of the HHC. In both the channel flow model and wedge extrusion model, the ZSZ acts as a backstop normal fault along which the high-grade metamorphic rocks of the HHC of Zanskar are exhumed. In contrast, the recently proposed tectonic wedging model argues that the ZSZ and the MSZ correspond to one single detachment system that operates as a subhorizontal backthrust off of the MCT. Thus, the kinematic evolution of the two shear zones, the ZSZ and the MSZ, and their structural, metamorphic and chronological relations appear to be diagnostic features for discriminating the different models. In this paper, structural, metamorphic and geochronological data demonstrate that the MSZ and the ZSZ experienced two distinct kinematic evolutions. As such, the data presented in this paper rule out the hypothesis that the MSZ and the ZSZ constitute one single detachment system, as postulated by the tectonic wedging model. Structural, metamorphic and geochronological data are used to present an alternative tectonic model for the large-scale doming in the NW Indian Himalaya involving early NE-directed tectonics, weakness in the upper crust, reduced erosion at the orogenic front and rapid exhumation along both the ZSZ and the MSZ.

39Ar-40Ar geochronology of mono- and polymetamorphic basement rocks

The 39Ar-40Ar technique is often used to date the metamorphic evolution of basement rocks. The present review article examines systematic aspects of the K-Ar decay system in different mineral chronometers frequently found in mono- and polymetamorphic basements (amphibole, biotite, muscovite/phengite, K-feldspar). A key observation is that the measured dissolution rate of silicates in aqueous fluids is many orders of magnitude faster, and has a much lower activation energy, than the rate of Fickian diffusion of Ar. The effects of this inequality are patchy age zonations, very much like those observed in many U-Pb chronometers, unaccompanied by intra-crystalline bell¬shaped Ar loss profiles. Recognizing the importance of the respective rate constants in field situations leads to re-evaluating the ages and the interpretive paradigms in classic examples such as the Central Alpine "Lepontine" amphibolite event and the Western Alpine eclogitic event.