Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

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.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

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.

Porosity and density of the AND-1B sediment core

A study of density and porosity is presented for the 1285-m-long AND-1B core recovered from a flexural moat in the McMurdo Sound (Antarctica) in order to interpret sediment consolidation in an ice-proximal location on the Antarctic shelf. Various lithologies imply environmental changes from open marine to subglacial, and are numerically expressed in high-resolution whole-core wet-bulk density (WBD). Grain density data interpolated from discrete samples range from 2.14 to 3.85 g/cm3 and are used to calculate porosity from WBD in order to avoid the 5%-15% overestimation and underestimation of porosities obtained by standard methods. The trend of porosity extends from 0.5 near the top (Pleistocene) to 0.2 at the bottom (Miocene). Porosity fluctuations in different lithologies are superimposed with 0.2-0.3 in sequences younger than ca. 1 Ma and 0.5-0.8 in Pliocene diatomites. The AND-1B porosities and void ratios of Pliocene diatomites and Pleistocene mudstones exhibit a large negative offset compared to modern lithological analogs and their consolidation trends. This offset cannot be explained in terms of the effective stress at the AND-1B site. The effective stress ranges from 0 to 4000 kPa in the upper 600 m, and reaches 13,000 kPa at the base of the AND-1B hole. We suggest an excess of effective overburden stress of ~1700 and ~6000 kPa to explain porosities in Pliocene diatomites and Pleistocene mudstones, respectively. This is interpreted as glacial preconsolidation by subsequently grounded ice sheets under subpolar to polar, followed by colder polar types of glaciations. Information on Miocene consolidation is sparse due to alteration by diagenesis.

Susceptibility, conductivity and porosity along a GEM-SHARK profile

The coastal systems, are often subjected to high anthropogenic pressure, which makes it necessary to develop new techniques to assess the environmental impacts caused by such human activity. This paper presents the first results obtained during the development and implementation of a new equipment of submarine geophysics survey oriented to integrated coastal zone management (ICZM). It is based on the drag of a submarine in contact with the sea-bottom. The submarine is equipped with an electromagnetic sensor which allows the measurement of the magnetic susceptibility and electrical conductivity of the surface sediments continuously and to a depth of sediment of 40 cm. This system, once improved, will allow us to obtain valuable information for monitoring the environmental quality of coastal areas.

Porosity contributions of radiolarians in ODP Site 110-671 and Hole 110-672A sediments

This study quantitatively addresses the significance of porosity within radiolarian tests in the décollement zone at the toe of the northern Barbados accretionary prism. Quantification was accomplished using scanning electron microscope images of core samples taken from Ocean Drilling Program (ODP) Sites 671 and 672, representing the décollement and proto-décollement, respectively. The décollement is localized to a radiolarian claystone, and its depth correlates with a low-density anomaly that has been attributed to high porosity at all relevant ODP drilling sites in the area (Moore, Klaus, et al., 1998, doi:10.2973/odp.proc.ir.171A.1998; Shipley, Ogawa, Blum, et al., 1995, doi:10.2973/odp.proc.ir.153.1995; Mascle, Moore, et al., 1988, doi:10.2973/odp.proc.ir.110.1988). Porosity in the décollement zone is presumably lost between Sites 672 and 671 because of shear enhanced consolidation (Moore et al., 1998, doi:10.1130/0091-7613(1998)026<0811:CPDIAE>2.3.CO;2).

(Table 1) Bulk and grain density, porosity, and p- and s-wave velocity characteristics of ODP Hole 163-990A sediments

Seismic velocities in rocks are influenced by the properties of the solid, the pore fluid, and the pore space. Cracks dramatically affect seismic velocities in rocks; their influence on the effective elastic moduli of rocks depends on their shape and concentration. Thin cracks (or fractures) substantially lower the moduli of a rock relative to the effect of spherical voids (or vesicles), and lower moduli are reflected by lower P- and S-wave velocities. The objective of this research is to determine the types and concentrations of cracks and their influence on the seismic properties of subaerially erupted basalts drilled from Hole 990A on the Southeast Greenland margin during Ocean Drilling Program Leg 163. Ellipsoidal cracks are used to model the voids in the rocks. The elastic moduli of the solid (grains) are also free parameters in the inverse modeling procedure. The apparent grain moduli reflect a weighted average of the moduli of the constituent minerals (e.g., plagioclase, augite, and clay minerals). The results indicate that (1) there is a strong relationship between P-wave velocity and porosity, suggesting a similarity of pore shape distributions, (2) the distribution of crack types within the massive, central region of aa flows from Hole 990A is independent of total porosity, (3) thin cracks are the first to be effectively sealed by alteration products, and (4) grain densities (an alteration index) and apparent grain moduli of the basalt samples are directly related.

Sediment description, CaCO3, density and porosity at DSDP Site 89-585

Siliceous sediments and sedimentary rocks occur as chert and silicified chalk, limestone, and claystone in Site 585 lower Miocene to Campanian sediments, with one older occurrence of chert near the Cenomanian/Turonian boundary. The recovered drill breccia in the Miocene to middle Eocene interval is dominated by bright red, orange, yellow, and brown chips and fragments of chert. In early Eocene and older sediments gray silicified limestone and yellowish brown chert fragments predominate. Recovery is poor in cores with chert because chert tends to fracture into smaller pieces that escape the drill and because the hard chert fragments grind away other sediments during rotary drilling. Thin-section and hand-sample studies show complex diagenetic histories of silicification (silica pore infill) and chertification (silica replacement of host rock). Multiple events of silicification can occur in the same rocks, producing chert from silicified limestone. Despite some prior silicification, silicified limestone is porous enough to provide conduits for dissolved silica-charged pore waters. Silicification and chert are more abundant in the coarser parts of the sedimentary section. These factors reflect the importance of porosity and permeability as well as chemical and lithologic controls in the process of silica diagenesis.