Fluid effect on hydraulic fracture propagation behavior: a comparison between water and supercritical CO2‐like fluid

The initiation of hydraulic fractures during fluid injection in deep formations can be either engineered or induced unintentionally. Upon injection of CO₂, the pore fluids in deep formations can be changed from oil/saline water to CO₂ or CO₂ dominated. The type of fluid is important not only because the fluid must fracture the rock, but also because rocks saturated with different pore fluids behave differently. We investigated the influence of fluid properties on fracture propagation behavior by using the cohesive zone model in conjunction with a poroelasticity model. Simulation results indicate that the pore pressure fields are very different for different pore fluids even when the initial field conditions and injection schemes (rate and time) are kept the same. Low viscosity fluids with properties of supercritical CO₂ will create relatively thin and much shorter fractures in comparison with fluids exhibiting properties of water under similar injection schemes. Two significant times are recognized during fracture propagation: the time at which a crack ceases opening and the later time point at which a crack ceases propagating. These times are very different for different fluids. Both fluid compressibility and viscosity influence fracture propagation, with viscosity being the more important property. Viscosity can greatly affect hydraulic conductivity and the leak‐off coefficient. This analysis assumes the in‐situ pore fluid and injected fluid are the same and the pore space is 100% saturated by that fluid at the beginning of the simulation.     

Fluid environment for preservation of pore spaces in a deep dolomite reservoir

Well TS1 reveals many uncemented pores and vugs at depths of more than 8000 m in a deep Cambrian dolomite reservoir in the Tarim Basin, northwestern China. The fluid environment and mechanism required for the preservation of reservoir spaces have yet not been well constrained. Carbon, oxygen, and strontium isotope compositions and fluid inclusion data suggest two types of fluids, meteoric water and hydrothermal fluid, affecting the Lower Paleozoic carbonate reservoirs in the Tarim Basin. Based on simulation using a thermodynamic model for H₂O‐CO₂‐NaCl‐CaCO₃ system, meteoric water has the ability to continuously dissolve carbonate minerals during downward migration from the surface to deep strata until it reaches a transition depth, below which it will begin to precipitate carbonate minerals to fill preexisting pore spaces. In contrast, hydrothermal fluid has the ability to dissolve carbonate in deep strata and precipitate carbonate in shallow strata during upward migration. Based on the dissolution–precipitation characteristics of the two types of fluids, the ideal fluid environment for the preservation of preexisting reservoir spaces occurs when carbonate reservoir is neither in the CaCO₃ precipitation domain of meteoric water nor in the CaCO₃ precipitation domain of hydrothermal fluid. Taking the Lower Paleozoic carbonate reservoirs in the north uplift area as an example, the spaces in the deep Cambrian dolomite reservoir near well TS1 were seldom filled because thick Ordovician deposits blocked meteoric water from migrating downward into the Cambrian dolomite reservoir and because the Cambrian dolomite reservoir has been in the domain of hydrothermal dissolution since the Permian. The deep carbonate layers in basins elsewhere with a similar fluid environment may have high uncemented porosity and consequently have good hydrocarbon exploration potential.     

Main and trace elements in KTB-VB fluid: composition and hints to its origin

The pilot hole (VB) of the German Continental Deep Drilling Program (KTB) was drilled to a depth of 4000 m, where large amounts of free fluids were met. The KTB‐VB 4000 m fluid can be related to either Mesozoic seawater or formation water from Permo‐Carboniferous sedimentary rocks of the Weiden embayment. During the Upper Cretaceous uplift of the Bohemian Massif both fluids could have passed organic‐rich Triassic to Carboniferous formations of the Weiden embayment before invading the uplifted and fractured basement rocks of Devonian amphibolites and metagabbros, where the chemical composition of the fluids was changed by albitization, adularization, and chloritization. Results of chemical mass balances for both sources are presented. In order to concentrate the formation water from the Weiden embayment significant amphibolitization has to be assumed. During a 1‐year pumping test the chemical composition of the 4000 m fluids remained constant. The accuracy of chemical analyses is critically reviewed. An improved preconcentration method of rare earth elements and yttrium in high‐Ca‐bearing saline fluids is described.     

Infrared radiofluorescence (IR-RF) dating of middle pleistocene fluvial archives of the Heidelberg Basin (Southwest Germany)

The infrared radiofluorescence (IR-RF) dating technique was applied to eight fluvial samples that were collected from two sediment cores at the Heidelberg Basin located near Viernheim and Ludwigshafen in southwest Germany. Based on the IR-RF derived ages of the samples it was possible to establish a chronological framework for the Mid-Pleistocene fluvial deposits of the Heidelberg Basin. The results allow us to distinguish between four main periods of aggradation. The lowermost sample taken from 100 m core depth lead to an IR-RF age of 643 ± 28 ka pointing to a Cromerian period of aggradation (OIS 17–16). For the Elsterian it is now possible to distinguish between two aggradation periods, one occurring during the Lower Elsterian period (OIS 15) and a second during the Upper Elsterian period (OIS 12–11). For the so called Upper interlayer (or “Oberer Zwischenhorizont” — a layer of organic-rich and finer-grained deposits), the IR-RF results point to a deposition age of around 300 ka, with samples taken directly on top and out of this layer yielding IR-RF ages of 288 ± 19 ka and 302 ± 19 ka, respectively. Hence, the measured IR-RF ages clearly point to a deposition during the Lower Saalian period (OIS 9–8) whereas earlier studies assumed a Cromerian age for the sediments of the Upper Interlayer based on pollen records and also mollusc fauna. The new IR-RF dataset indicates that significant hiatuses are present within the fluvial sediment successions. In particular the Eemian and Upper Saalian deposits are missing in this part of the northern Upper Rhine Graben, as the 300 ka deposits are directly overlain by Weichselian fluvial sediments. It is obvious that time periods of increased fluvial aggradation were interrupted by time periods of almost no aggradation or erosion which should have been mainly triggered by phases of increased and decreased subsidence of the Heidelberg Basin.     

Impacts of Pleistocene glacial loading on abnormal pore‐water pressure in the eastern Michigan Basin

The hydromechanical effects of Pleistocene glacial loading on the Michigan Basin are assessed using numerical analysis based on coupled stress and pore‐water pressure. The two‐dimensional model domain included the Basin cross section and extended 10 km into the Precambrian. In the analysis, we considered the effects of the number of glacial loading cycles, the presence and connectedness of a deep Cambrian aquifer, the direction of glacial advance, the effect of a wet versus dry glacier/soil interface, topographic effects, density‐driven flow effects, and lithosphere flexure on the development of anomalous pressures. The modeling results were compared with data collected from highly instrumented wells completed in the eastern margin of the Basin. The present‐day results define regions of significant underpressure in the upper Ordovician and lower Silurian formations characterized by very low hydraulic conductivity and regions of overpressure where hydraulic conductivity is higher. To achieve the degree of underpressure observed in the instrumented wells using the model, a specific loading cycle applied over 100 000 years was repeated four times. As the number of loading cycles increased, the Paleozoic formations reached a state where the underpressures remain constant. Our results also illustrate the difference in poroelastic modeling between the application of mechanical loads on the land surface and the application of an equivalent hydraulic head, where the latter developed overpressures rather than the observed underpressures. The modeling also shows that the overpressures observed in the Cambrian formations are most likely to be the result of density‐driven flow. Finally, the simulations show that the effects of lithosphere flexure in the hydromechanical models results in the development of lateral stresses that generate overpressures rather than underpressures in the southern half of the domain. As there are no suitable observation points, these results remain unconfirmed, and further study is warranted.     

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.     

Hydrogeochemical dynamics affecting steam‐heated pools at El Chichón Crater (Chiapas – Mexico)

El Chichón is an active volcano located in the north‐western Chiapas, southern Mexico. The crater hosts a lake, a spring, named Soap Pool, emerging from the underlying volcanic aquifer and several mud pools/hot springs on the internal flanks of the crater which strongly interact with the current fumarolic system (steam‐heated pools). Some of these pools, the crater lake and a cold spring emerging from the 1982 pumice deposits, have been sampled and analysed. Water–volcanic gas interactions determine the heating (43–99°C) and acidification (pH 2–4) of the springs, mainly by H₂S oxidation. Significantly, in the study area, a significant NH₃ partial pressure has been also detected. Such a geochemically aggressive environment enhances alteration of the rock in situ and strongly increases the mineralization of the waters (and therefore their electrical conductivity). Two different mineralization systems were detected for the crater waters: the soap pool‐lake (Na⁺/Cl^? = 0.4, Na/Mg>10) and the crater mud pools (Na⁺/Cl^? > 10, Na/Mg < 4). A deep boiling, Na⁺‐K⁺‐Cl^?‐rich water reservoir generally influences the Soap Pool‐lake, while the mud pool is mainly dominated by water‐gas–rock interactions. In the latter case, conductivity of sampled water is directly proportional to the presence of reactive gases in solution. Therefore, chemical evolution proceeds through neutralization due to both rock alteration and bacterial oxidation of ammonium to nitrate. The chemical compositions show that El Chichón aqueous fluids, within the crater, interact with gases fed by a geothermal reservoir, without clear additions of deep magmatic fluids. This new geochemical dataset, together with previously published data, can be used as a base line with which to follow‐up the activity of this deadly volcano.     

Fluid mapping in deeply buried Ordovician paleokarst reservoirs in the Tarim Basin,western China

The storage spaces within deeply buried Ordovician paleokarst reservoirs in the Tarim Basin are mostly secondary and characterized by strong heterogeneity and some degree of anisotropy. The types of fluids that fill the spaces within these reservoirs are of great importance for hydrocarbon exploration and exploitation. However, fluid identification from seismic data is often controversial in this area because the seismic velocity for this particular reservoir could be significantly influenced by many factors, including pore shapes, porosity, fluid types, and mineral contents. In this study, we employ the differential effective medium‐Gassmann rock physics model to interpret and discuss the characteristics of conventional karstic carbonate reservoirs in the Tarim Basin that are filled with different fluids (oil, gas, and water) using logging data and thus objectively build corresponding fluid identification criteria. These criteria are subsequently evaluated by amplitude versus offset (AVO) forward analysis based on typical logging data and further applied to ascertain the reservoir fluid types in two different areas in the Tarim Basin based on prestack inversion results. For conventional carbonate reservoirs, gas can be distinguished from heavy oil and water, but heavy oil and water are broadly similar on seismic data. For condensate carbonate reservoirs, water can be differentiated from light oil (i.e., condensates) and gas, but light oil and gas demonstrate substantial similarities in terms of their seismic responses. The predicted fluid results are in good agreement with the results of drilling and oil testing. In particular, modeling the seismically resolvable reservoirs in the carbonate strata in the Tarim Basin, which have needle‐ and sphere‐shaped storage spaces (pore aspect ratio > 0.3) and clay content that is lower than 5%, indicates that fluid properties could be properly evaluated if the porosity is larger than 5% for conventional carbonate reservoirs and >7% for condensate carbonate reservoirs.     

From deep to shallow fluid reservoirs: evolution of fluid sources during exhumation of the Sierra Almagrera,Betic Cordillera,Spain

The combination of structural, geochemical and palaeotopographic data proves to be an efficient tool to understand fluid transfers in the crust. This study discriminates shallow and deep fluid reservoirs on both sides of the brittle–ductile transition under an extensional regime and points out the role of major transcurrent fault activity in this palaeohydrogeological setting. Palaeofluids trapped in quartz and siderite–barite veins record the transfer of fluids and metal solute species during the Neogene exhumation of the Sierra Almagrera metamorphic belt. Ductile then brittle–ductile extensional quartz veins formed from a deep fluid reservoir, trapping metamorphic secondary brines containing low‐density volatile phases derived from the dissolution of Triassic evaporites. During exhumation, low‐salinity fluids percolated within the brittle domain, as shown by transgranular fluid inclusion planes affecting previous veins. These observations indicate the opening of the system during Serravalian to early Tortonian times and provide evidence for the penetration of surficial fluids of meteoric or basinal origin into the upper part of the brittle–ductile transition. During exhumation, synsedimentary transcurrent tectonic processes occurred from late Tortonian times onwards, while marine conditions prevailed at the Earth's surface. At depth in the brittle domain, quartz veins associated with haematite record a return to high‐salinity fluid circulation suggesting an upward transfer fed from a lower reservoir. During the Messinian, ongoing activity of the trans‐Alboran tectono‐volcanic trend led to the formation of ore deposits. Reducing fluids caused the formation of siderite and pyrite ores. The subsequent formation of galena and barite may be related to an increase of temperature. The high salinity and Cl/Br ratio of the fluids suggest another source of secondary brine derived from dissolved Messinian evaporites, as corroborated by the δ³⁴S signature of barite. These evaporites preceded the main sea‐level drop related to the peak of the salinity crisis (5.60–5.46 Ma).     

Downward hydrocarbon migration predicted from numerical modeling of fluid overpressure in the Paleozoic Anticosti Basin,eastern Canada

The Anticosti Basin is a large Paleozoic basin in eastern Canada where potential source and reservoir rocks have been identified but no economic hydrocarbon reservoirs have been found. Potential source rocks of the Upper Ordovician Macasty Formation overlie carbonates of the Middle Ordovician Mingan Formation, which are underlain by dolostones of the Lower Ordovician Romaine Formation. These carbonates have been subjected to dissolution and dolomitization and are potential hydrocarbon reservoirs. Numerical simulations of fluid‐overpressure development related to sediment compaction and hydrocarbon generation were carried out to investigate whether hydrocarbons generated in the Macasty Formation could migrate downward into the underlying Mingan and Romaine formations. The modeling results indicate that, in the central part of the basin, maximum fluid overpressures developed above the Macasty Formation due to rapid sedimentation. This overpressured core dissipated gradually with time, but the overpressure pattern (i.e. maximum overpressure above source rock) was maintained during the generation of oil and gas. The downward impelling force associated with fluid‐overpressure gradients in the central part of the basin was stronger than the buoyancy force for oil, whereas the buoyancy force for gas and for oil generated in the later stage of the basin is stronger than the overpressure‐related force. Based on these results, it is proposed that oil generated from the Macasty Formation in the central part of the basin first moved downward into the Mingan and Romaine formations, and then migrated laterally up‐dip toward the basin margin, whereas gas throughout the basin and oil generated in the northern part of the basin generally moved upward. Consequently, gas reservoirs are predicted to occur in the upper part of the basin, whereas oil reservoirs are more likely to be found in the strata below the source rocks. Geofluids (2010) 10 , 334–350     

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     

On the geochemistry of gases and noble gas isotopes (including 222Rn) in deep crustal fluids: the 4000 m KTB-pilot hole fluid production test 2002–03

The concentrations of H₂, O₂, CO₂, and concentrations and isotopic composition of the noble gases (including ²²²Rn), N₂, CH₄, and higher hydrocarbons dissolved in 4000 m deep‐seated fluids from a 12‐month fluid production test in the KTB pilot hole were analyzed. This determination of the gas geochemistry during the test in combination with the knowledge of the hydraulic data provides relevant information about the fluid hydraulics of the deep system. All gas concentrations and isotopic signatures, except for ²²²Rn, showed constancy during the course of the test. This, in combination with large fluid flow rates at a moderate water table drawdown, imply an almost infinite fluid reservoir in 4000 m depth. From the change in ²²²Rn‐activity as a function of pump rate, the contribution of smaller and wider pores to the overall fluid flow in an aquifer can be deduced. This ²²²Rn‐activity monitoring proved therefore to be a valuable instrument for the qualitative observation of the scavenging of pore and fracture surfaces, a hydraulic feature invisible to standard hydraulic testing tools. The observance of this scavenging effect is due to (i) the continuous on‐line geochemical monitoring, (ii) the durability of the test, (iii) a change in pump rate during the course of the test, and (iv) due to the short half‐life of ²²²Rn. The fluids have a 5.9% mantle He component, and a δ²¹Ne excess of 14%, and a noble gas model age of about (5.5–6.2) ± 2.0 Myr. The mean N₂/Ar‐ratio of 516 and δ¹⁵N‐data of about +1.5‰ indicates sedimentary or metamorphic origin of N₂. The hydrocarbons, amounting to 33 vol.% in the gas phase, are derived from thermal decomposition of marine organic matter of low maturity. But a key question, the identification of the potential source region of the fluids and the migration pathway, is still unidentified.     

Water for the City,Ruins for the Country: Archaeology of the New York City Watershed

Two New York City institutions, the Board of Water Supply (BWS) and Department of Environmental Protection (DEP), have shaped rural regions far outside city boundaries. The BWS depopulated places selected for reservoir construction. Residents were evicted and towns were demolished then submerged. Those who remained struggled to reorganize their lives amid the landscape clearance. Once the reservoirs were complete, the DEP replaced the BWS as the institution in charge of ensuring the city’s water supply. The DEP Police patrols around the reservoirs and enforces land-use regulations. Archaeological survey of city-owned watershed lands reveals a scarred landscape of ongoing colonial conflict.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

Residence times and source ages of deep crustal fluids: interpretation of 129I and 36Cl results from the KTB‐VB drill site,Germany

We present here ¹²⁹I/I and ³⁶Cl/Cl ratios, together with halogen concentrations in crustal fluids from the continental deep drill site (KTB‐VB) in Germany, where fluids were collected from 4000 m depth during a pump test carried out in 2002 and 2003. Compared with seawater, the fluids are enriched by factors of 2, 8 and 40 for Cl^?, Br and I, respectively, and show little variation over the test period. The ¹²⁹I/I ratios are between 1700 and 4100 × 10^?15; the ³⁶Cl/Cl ratios are below 10 × 10^?15. Co‐variation between ¹²⁹I and ³⁶Cl concentrations in the fluids indicates that anthropogenic components are absent and that the ratios reflect an addition from crustal sources. The results suggest residence times of 10 Ma or more for the fluids in formations with uranium concentrations of 1 ppm. A minimum age of 30 Ma for the iodine source was derived from the correlation between ¹²⁹I and ³⁶Cl concentrations in the fluids. The results demonstrate that the halogen characteristics of the KTB fluids are very similar to those of other deep crustal fluids and that the combination of ¹²⁹I and ³⁶Cl systematics allows determination of residence times and source ages of such fluids.     

Spatial variations in the salinity of pore waters in northern deep water Gulf of Mexico sediments: implications for pathways and mechanisms of solute transport

Spatial variations in the salinity of pore waters in sedimentary basins can provide important insight into basin-scale hydrogeologic processes. Although there have been numerous studies of brine seeps in the deep water Gulf of Mexico, much less is known about porewater salinities in the vast areas between seeps. A study has been made of spatial variation in pore water salinities in sediments in an approximately 500 km by 200 km area of the northern deep water (water depth >500 m) Gulf of Mexico sedimentary basin (GOM) to provide insight into pathways and mechanisms of solute transport in this portion of the basin. A second objective was to document salinities in the upper 500 m of the sedimentary section, the approximate depth to which methane hydrates, a potential future energy resource, may be stable. Elevated salinities would reduce the P – T stability range of hydrates. Salinities were calculated from borehole logs using a dual electrical conductivity model. Even though much of the northern GOM is underlain by allochthonous salt most of the undisturbed shallow sedimentary section has not been permeated by hypersaline waters. These waters are limited to areas near brine seeps. Hypersaline waters having salinities in excess of 100 g l⁻¹ become more common at subseafloor depths of 2 km and greater. A field study at Green Canyon 65 and published numerical simulations of fluid flow above tabular salt bodies suggest that brines produced by salt dissolution migrate laterally and pond above salt and/or within minibasins and that the dominant mechanism of vertical solute transport is a combination of compaction-driven advection and diffusion, not large-scale thermohaline overturn. Superimposed on this diffuse upward flux of dissolved salt is the more focused and localized expulsion of saline fluids up faults.     

The upper continental crust, an aquifer and its fluid: hydaulic and chemical data from 4 km depth in fractured crystalline basement rocks at the KTB test site

Detailed information on the hydrogeologic and hydraulic properties of the deeper parts of the upper continental crust is scarce. The pilot hole of the deep research drillhole (KTB) in crystalline basement of central Germany provided access to the crust for an exceptional pumping experiment of 1‐year duration. The hydraulic properties of fractured crystalline rocks at 4 km depth were derived from the well test and a total of 23100 m³ of saline fluid was pumped from the crustal reservoir. The experiment shows that the water‐saturated fracture pore space of the brittle upper crust is highly connected, hence, the continental upper crust is an aquifer. The pressure–time data from the well tests showed three distinct flow periods: the first period relates to wellbore storage and skin effects, the second flow period shows the typical characteristics of the homogeneous isotropic basement rock aquifer and the third flow period relates to the influence of a distant hydraulic border, probably an effect of the Franconian lineament, a steep dipping major thrust fault known from surface geology. The data analysis provided a transmissivity of the pumped aquifer T = 6.1 × 10^?6 m² sec^?1, the corresponding hydraulic conductivity (permeability) is K = 4.07 × 10^?8 m sec^?1 and the computed storage coefficient (storativity) of the aquifer of about S = 5 × 10^?6. This unexpected high permeability of the continental upper crust is well within the conditions of possible advective flow. The average flow porosity of the fractured basement aquifer is 0.6–0.7% and this range can be taken as a representative and characteristic values for the continental upper crust in general. The chemical composition of the pumped fluid was nearly constant during the 1‐year test. The total of dissolved solids amounts to 62 g l^?1 and comprise mainly a mixture of CaCl₂ and NaCl; all other dissolved components amount to about 2 g l^?1. The cation proportions of the fluid (X_Ca approximately 0.6) reflects the mineralogical composition of the reservoir rock and the high salinity results from desiccation (H₂O‐loss) due to the formation of abundant hydrate minerals during water–rock interaction. The constant fluid composition suggests that the fluid has been pumped from a rather homogeneous reservoir lithology dominated by metagabbros and amphibolites containing abundant Ca‐rich plagioclase.     

Long‐term chemical evolution and modification of continental basement brines – a field study from the Schwarzwald,SW Germany

Highly saline, deep‐seated basement brines are of major importance for ore‐forming processes, but their genesis is controversial. Based on studies of fluid inclusions from hydrothermal veins of various ages, we reconstruct the temporal evolution of continental basement fluids from the Variscan Schwarzwald (Germany). During the Carboniferous (vein type i), quartz–tourmaline veins precipitated from low‐salinity (<4.5wt% NaCl + CaCl₂), high‐temperature (≤390°C) H₂O‐NaCl‐(CO₂‐CH₄) fluids with Cl/Br mass ratios = 50–146. In the Permian (vein type ii), cooling of H₂O‐NaCl‐(KCl‐CaCl₂) metamorphic fluids (T ≤ 310°C, 2–4.5wt% NaCl + CaCl₂, Cl/Br mass ratios = 90) leads to the precipitation of quartz‐Sb‐Au veins. Around the Triassic–Jurassic boundary (vein type iii), quartz–haematite veins formed from two distinct fluids: a low‐salinity fluid (similar to (ii)) and a high‐salinity fluid (T = 100–320°C, >20wt% NaCl + CaCl₂, Cl/Br mass ratios = 60–110). Both fluids types were present during vein formation but did not mix with each other (because of hydrogeological reasons). Jurassic–Cretaceous veins (vein type iv) record fluid mixing between an older bittern brine (Cl/Br mass ratios ~80) and a younger halite dissolution brine (Cl/Br mass ratios >1000) of similar salinity, resulting in a mixed H₂O‐NaCl‐CaCl₂ brine (50–140°C, 23–26wt% NaCl + CaCl₂, Cl/Br mass ratios = 80–520). During post‐Cretaceous times (vein type v), the opening of the Upper Rhine Graben and the concomitant juxtaposition of various aquifers, which enabled mixing of high‐ and low‐salinity fluids and resulted in vein formation (multicomponent fluid H₂O‐NaCl‐CaCl₂‐(SO₄‐HCO₃), 70–190°C, 5–25wt% NaCl‐CaCl₂ and Cl/Br mass ratios = 2–140). The first occurrence of highly saline brines is recorded in veins that formed shortly after deposition of halite in the Muschelkalk Ocean above the basement, suggesting an external source of the brine's salinity. Hence, today's brines in the European basement probably developed from inherited evaporitic bittern brines. These were afterwards extensively modified by fluid–rock interaction on their migration paths through the crystalline basement and later by mixing with younger meteoric fluids and halite dissolution brines.     

Hydrothermal dolomitization,mixing corrosion and deep burial porosity formation: numerical results from 1‐D reactive transport models

Major corrosion has been found at depth in carbonate hydrocarbon reservoirs from different geologic provinces. Fluid inclusion microthermometry and stable isotopic compositions of carbonate cements, predating major corrosion, constrain the interpretation of the evolution of parental fluids during progressive burial and prior to the major corrosion event. Post‐major corrosion mineral paragenesis includes pyrite (‐marcasite), anhydrite, kaolinite (dickite) and fluorite. Although the post‐corrosion mineral paragenesis represents minor volumes of rock, it may provide valuable insights into the post‐corrosion brine chemistry. Using reactive transport numerical models, the roles of cooling and/or mixing of brines on corrosion have been evaluated as controls for dolomitization, deep burial corrosion and precipitation of the post‐corrosion mineral paragenesis. Modelling results show that cooling of deep‐seated fluids moving upward along a fracture may cause minor calcite dissolution and porosity generation. Significant dolomitization along a fracture zone and nearby host‐rock only occurs when deep‐seated fluids have high salinities (4 mol Cl kg^?1 of solution) and Darcian flow rates are relatively high (1 m³ m^?2 year^?1). Only minor volumes of quartz and fluorite precipitate in the newly formed porosity. Moreover, modelling results cannot reproduce the authigenic precipitation of kaolinite (dickite at high temperatures) by cooling. As an alternative to cooling as a cause of corrosion, mixing between two brines of different compositions and salinities is represented by two main cases. One case consists of the flow up along a fracture of deep‐seated fluids with higher salinities than the fluid in the wall rock. Dolomite does not precipitate at a fracture zone. Nevertheless, minor volumes of dolomite are formed away from the fracture. The post‐corrosion mineral paragenesis can be partly reproduced, and the results are comparable to those obtained from cooling calculations. Minor volumes of quartz and fluorite are formed, and kaolinite‐dickite does not precipitate. The major outputs of this scenario are calcite dissolution and slight net increase in porosity. A second case corresponds to the mixing of low salinity deep‐seated fluids, flowing up along fractures, with high salinity brines within the wall rock. Calculations predict major dissolution of calcite and precipitation of dolomite. The post‐corrosion mineral paragenesis can be reproduced. High volumes of quartz, fluorite and kaolinite‐dickite precipitate and may even completely occlude newly formed porosity.     

Deep geothermal groundwater flow in the Seferihisar–Balçova area,Turkey: results from transient numerical simulations of coupled fluid flow and heat transport processes

The Seferihisar–Balçova Geothermal system (SBG) is characterized by complex temperature and hydrochemical anomalies. Previous geophysical and hydrochemical investigations suggest that hydrothermal convection in the faulted areas of the SBG and recharge flow from the Horst may be responsible for the observed patterns. A numerical model of coupled fluid flow and heat transport processes has been built in order to study the possible fluid dynamics of deep geothermal groundwater flow in the SBG. The results support the hypothesis derived from interpreted data. The simulated scenarios provide a better understanding of the geophysical conditions under which the different fluid dynamics develop. When recharge processes are weak, the convective patterns in the faults can expand to surrounding reservoir units or below the seafloor. These fault‐induced drag forces can cause natural seawater intrusion. In the Melange of the Seferihisar Horst, the regional flow is modified by buoyant‐driven flow focused in the series of vertical faults. As a result, the main groundwater divide can shift. Sealing caprocks prevent fault‐induced cells from being overwhelmed by vigorous regional flow. In this case, over‐pressured, blind geothermal reservoirs form below the caprocks. Transient results showed that the front of rising hot waters in faults is unstable: the tip of the hydrothermal plumes can split and lead to periodical temperature oscillations. This phenomenon known as Taylor–Saffman fingering has been described in mid‐ocean ridge hydrothermal systems. Our findings suggest that this type of thermal pulsing can also develop in active, faulted geothermal systems. To some extent, the role of an impervious fault core on the flow patterns has been investigated. Although it is not possible to reproduce basin‐scale transport processes, this first attempt to model deep groundwater geothermal flow in the SBG qualitatively supported the interpreted data and described the different fluid dynamics of the basin. Geofluids (2010) 10 , 388–405