Hydrogeology of the Krafla geothermal system,northeast Iceland

The Krafla geothermal system is located in Iceland's northeastern neovolcanic zone, within the Krafla central volcanic complex. Geothermal fluids are superheated steam closest to the magma heat source, two‐phase at higher depths, and sub‐boiling at the shallowest depths. Hydrogen isotope ratios of geothermal fluids range from ?87‰, equivalent to local meteoric water, to ?94‰. These fluids are enriched in ¹⁸O relative to the global meteoric line by +0.5–3.2‰. Calculated vapor fractions of the fluids are 0.0–0.5 wt% (~0–16% by volume) in the northwestern portion of the geothermal system and increase towards the southeast, up to 5.4 wt% (~57% by volume). Hydrothermal epidote sampled from 900 to 2500 m depth has δD values from ?127 to ?108‰, and δ¹⁸O from ?13.0 to ?9.6‰. Fluids in equilibrium with epidote have isotope compositions similar to those calculated for the vapor phase of two‐phase aquifer fluids. We interpret the large range in δD_EPIDOTE and δ¹⁸O_EPIDOTE across the system and within individual wells (up to 7‰ and 3.3‰, respectively) to result from variable mixing of shallow sub‐boiling groundwater with condensates of vapor rising from a deeper two‐phase reservoir. The data suggest that meteoric waters derived from a single source in the northwest are separated into the shallow sub‐boiling reservoir, and deeper two‐phase reservoir. Interaction between these reservoirs occurs by channelized vertical flow of vapor along fractures, and input of magmatic volatiles further alters fluid chemistry in some wells. Isotopic compositions of hydrothermal epidote reflect local equilibrium with fluids formed by mixtures of shallow water, deep vapor condensates, and magmatic volatiles, whose ionic strength is subsequently derived from dissolution of basalt host rock. This study illustrates the benefits of combining phase segregation effects in two‐phase systems during analysis of wellhead fluid data with stable isotope values of hydrous alteration minerals when evaluating the complex hydrogeology of volcano‐hosted geothermal systems.     

3D numerical modeling of hydrothermal processes during the lifetime of a deep geothermal reservoir

Understanding hydrothermal processes during production is critical to optimal geothermal reservoir management and sustainable utilization. This study addresses the hydrothermal (HT) processes in a geothermal research doublet consisting of the injection well E GrSk3/90 and production well Gt GrSk4/05 at the deep geothermal reservoir of Groß Schönebeck (north of Berlin, Germany) during geothermal power production. The reservoir is located between ?4050 to ?4250 m depth in the Lower Permian of the Northeast German Basin. Operational activities such as hydraulic stimulation, production (T = 150°C; Q = ?75 m³ h^?1; C = 265 g l^?1) and injection (T = 70°C; Q = 75 m³ h^?1; C = 265 g l^?1) change the HT conditions of the geothermal reservoir. The most significant changes affect temperature, mass concentration and pore pressure. These changes influence fluid density and viscosity as well as rock properties such as porosity, permeability, thermal conductivity and heat capacity. In addition, the geometry and hydraulic properties of hydraulically induced fractures vary during the lifetime of the reservoir. A three‐dimensional reservoir model was developed based on a structural geological model to simulate and understand the complex interaction of such processes. This model includes a full HT coupling of various petrophysical parameters. Specifically, temperature‐dependent thermal conductivity and heat capacity as well as the pressure‐, temperature‐ and mass concentration‐dependent fluid density and viscosity are considered. These parameters were determined by laboratory and field experiments. The effective pressure dependence of matrix permeability is less than 2.3% at our reservoir conditions and therefore can be neglected. The results of a three‐dimensional thermohaline finite‐element simulation of the life cycle performance of this geothermal well doublet indicate the beginning of thermal breakthrough after 3.6 years of utilization. This result is crucial for optimizing reservoir management. Geofluids (2010) 10 , 406–421     

Sampling and analysis of geothermal fluids

Sampling of geothermal fluids presents some problems not encountered when sampling surface and nonthermal ground waters. Specific collection techniques are required to obtain representative samples because of the elevated temperature and boiling of these fluids, the effect of exposing them to the atmosphere and cooling of the samples. Sample treatment during collection depends on the analytical method to be used. When sampling wet‐steam wells, both the liquid and the vapour fractions should be collected at the same fluid separation pressure. When sampling fumarole steam, maximum information is obtained if the total discharge is collected into a single container without separating the gas and the steam condensate fractions. Silica polymerization affects the solution pH. The only way to obtain reliable pH measurement of a water sample supersaturated with respect to amorphous silica is to measure it on site, before the onset of polymerization. This paper provides an outline of the geothermal sampling techniques and analytical methods currently in use in Iceland. Sampling of hot‐water and wet‐steam wells is described, as is sampling of hot springs, fumaroles and gas bubbling through hot‐spring waters. Detailed procedures are given for the analysis of total carbonate carbon and total sulphide sulphur in geothermal water and steam condensate samples.     

Hydrothermal,multiphase convection of H2O‐NaCl fluids from ambient to magmatic temperatures: a new numerical scheme and benchmarks for code comparison

Thermohaline convection of subsurface fluids strongly influences heat and mass fluxes within the Earth's crust. The most effective hydrothermal systems develop in the vicinity of magmatic activity and can be important for geothermal energy production and ore formation. As most parts of these systems are inaccessible to direct observations, numerical simulations are necessary to understand and characterize fluid flow. Here, we present a new numerical scheme for thermohaline convection based on the control volume finite element method (CVFEM), allowing for unstructured meshes, the representation of sharp thermal and solute fronts in advection‐dominated systems and phase separation of variably miscible, compressible fluids. The model is an implementation of the Complex Systems Modelling Platform CSMP++ and includes an accurate thermodynamic representation of strongly nonlinear fluid properties of salt water for magmatic‐hydrothermal conditions (up to 1000°C, 500 MPa and 100 wt% NaCl). The method ensures that all fluid properties are taken as calculated on the respective node using a fully upstream‐weighted approach, which greatly increases the stability of the numerical scheme. We compare results from our model with two well‐established codes, HYDROTHERM and TOUGH2, by conducting benchmarks of different complexity and find good to excellent agreement in the temporal and spatial evolution of the hydrothermal systems. In a simulation with high‐temperature, high‐salinity conditions currently outside of the range of both HYDROTHERM and TOUGH2, we show the significance of the formation of a solid halite phase, which introduces heterogeneity. Results suggest that salt added by magmatic degassing is not easily vented or accommodated within the crust and can result in dynamic, complex hydrologies.     

Thermal springs and heat flow in North America

Thermal springs are commonly thought to be an indicator of geothermal resource potential. However, there have been few analyses of the relationship between thermal springs and the underlying thermal regime. An examination of temperature and discharge rates for a large database of thermal springs in North America demonstrates that there is not a simple relationship between these measurements made at the surface and subsurface heat flow. Hydrogeological factors appear to exert strong controls on the temperature and discharge at these springs and should be carefully considered in geothermal resource assessments.     

The stress regime of the Western Canadian Sedimentary Basin

All available information relevant to in situ stress orientations and magnitudes in the Western Canadian Sedimentary Basin (WCSB) were examined to provide a better understanding of how regional stress fields may affect geothermal development. The smallest principal stress is horizontal over most of the Western Canadian Sedimentary Basin, and it varies in magnitude across the region. Horizontal stress trajectories show that S_Hmax axes are generally aligned SW–NE. A total of 1643 measurements of microfracture and minifracture closure pressures, leak‐off pressures and fracture breakdown pressures have been harnessed to map S_Hmin gradients across the basin at depths of 156–500, 500–1000, 1000–4185 and 2000–4185 m. Vertical stress magnitudes, calculated in 91 wells, showed that at constant depth, S_V increases towards the Canadian Rocky Mountains. Resultant regional stress maps show consistent trends in orientation of stress axes. As a result, predictions can be made that propagation axes of subsurface hydraulic fractures will be dominantly SW–NE, except over the Peace River Arch area, where they will trend more towards SSW–NNE. Engineered geothermal systems in the WCSB can be optimised by drilling horizontal wells parallel to S_Hmin.     

Mineralogy and fluid inclusion gas chemistry of production well mineral scale deposits at the Dixie Valley geothermal field,USA

At the Dixie Valley geothermal field, Nevada, USA, fluid boiling triggered the precipitation of carbonate scale minerals in concentric bands around tubing inserted into production well 28–33. When the tubing was removed, this mineral scale was sampled at 44 depth intervals between the wellhead and 1227 m depth. These samples provide a unique opportunity to evaluate the effects of fluid boiling on the scale mineralogy and geochemistry of the vapor and liquid phase. In this study, the mineralogy of the scale deposits and the composition of the fluid inclusion gases trapped in the mineral scales were analyzed. The scale consists mainly of calcite from 670–1112 m depth and aragonite from 1125 to 1227 m depth, with traces of quartz and Mg‐smectite. Mineral textures, including hopper growth, twinning, and fibrous growth in the aragonite and banded deposits of fine grained calcite crystals, are the result of progressive boiling. The fluid inclusion noncondensable gas was dominated by CO₂. However, significant variations in He relative to N₂ and Ar provide evidence that the geothermal reservoir consists of mixed source deeply circulating reservoir water and shallow, air saturated meteoric water. Gas analyses for many inclusions also showed higher CH₄ and H₂ relative to CO₂ than measured in gas sampled from this well, other production wells, and fumaroles. These inclusions are interpreted to have trapped CH₄‐ and H₂‐enriched gas resulting from early stages of boiling.     

Thermal waters of ‘tectonic origin’: the alkaline,Na‐HCO3 waters of the Rio Valdez geothermal area (Isla Grande de Tierra del Fuego,Argentina)

The geothermal area of Rio Valdez is located in the central portion of the Isla Grande de Tierra del Fuego (South Argentina), ten kilometers south of the southeastern sector of the Fagnano Lake. It consists of a series of thermal springs with low discharge rates (≤1 L/s) and temperatures in the range of 20–33°C distributed in an area of <1 km². The thermal springs are characterized by alkaline, Na‐HCO₃ waters with low salinity (0.53÷0.58 g/L), but relatively high fluoride contents (up to 19.4 mg/L). Their composition is the result of a slow circulation at depth, possibly through deep tectonic discontinuities connected with the Magallanes‐Fagnano Fault (MFF) system. According to geothermometric calculations, thermal waters reach temperatures in the range of 100–150°C and an almost complete chemical equilibrium with the alkali‐feldspars in the metavolcanic country rocks. The relatively high fluorine contents can be explained by the slow ascent and cooling of deep groundwaters followed by a progressive re‐equilibration with F‐bearing, hydrated Mg‐silicates, such as chlorite, which has been recognized as an abundant mineral in the metavolcanics of the Lemaire Formation and metapelites and metagraywackes of the Yahgán Formation. Finally, the isotopic composition of the investigated samples is consistent with the infiltration from local snow melting at altitudes in the range of 610–770 m asl. The comparison of our data with those collected in 1991 seems to suggest a possible progressive decline of the bulk thermal output in the near future. This possibility should be seriously considered before planning a potentially onerous exploitation of the resource. Presently, the only ways to exploit this geothermal resource by the population scattered in the area are the direct use of thermal waters and/or spa structures.     

Permeability of the continental crust: dynamic variations inferred from seismicity and metamorphism

The variation of permeability with depth can be probed indirectly by various means, including hydrologic models that use geothermal data as constraints and the progress of metamorphic reactions driven by fluid flow. Geothermal and metamorphic data combine to indicate that mean permeability ( k ) of tectonically active continental crust decreases with depth ( z ) according to log  k  ≈ −14–3.2 log  z , where k is in m² and z in km. Other independently derived, crustal-scale k – z relations are generally similar to this power-law curve. Yet there is also substantial evidence for local-to-regional-scale, transient, permeability-generation events that entail permeabilities much higher than these mean k – z relations would suggest. Compilation of such data yields a fit to these elevated, transient values of log  k  ≈ −11.5–3.2 log  z , suggesting a functional form similar to that of tectonically active crust, but shifted to higher permeability at a given depth. In addition, it seems possible that, in the absence of active prograde metamorphism, permeability in the deeper crust will decay toward values below the mean k – z curves. Several lines of evidence suggest geologically rapid (years to 10³ years) decay of high-permeability transients toward background values. Crustal-scale k – z curves may reflect a dynamic competition between permeability creation by processes such as fluid sourcing and rock failure, and permeability destruction by processes such as compaction, hydrothermal alteration, and retrograde metamorphism.     

Stable–unstable flow of geothermal fluids in fractured rock

Density-driven geothermal flow in 3-D fractured rock is investigated and compared with density-driven haline flow. For typical matrix and fracture hydraulic conductivities, haline flow tends to be unstable (convecting) while geothermal flow is stable (non-convecting). Thermal diffusivity is generally three orders of magnitude larger than haline diffusivity and, as a result, large heat conduction diminishes growth of geothermal instabilities while low mass diffusion enables formation of unstable haline 'fingering' within fractures. A series of thermal flow simulations is presented to identify stable and unstable conditions for a wide range of hydraulic conductivities for matrix and fractures. The classic Rayleigh stability criterion can be applied to classify these simulations when fracture aperture is very small. However, the Rayleigh criterion is not applicable when the porous matrix hydraulic conductivity is very small, because stabilizing fracture–matrix heat conduction is independent of matrix hydraulic conductivity. In that case, the numerically estimated critical fracture conductivity is nine orders of magnitude larger than the theoretically calculated critical fracture conductivity based on Rayleigh theory. The numerical stability analysis presented here may be used as a guideline to predict if a geothermal system in 3-D fractured rock is stable or unstable.