Geochemical fingerprinting of hydraulic fracturing fluids from Qusaiba Hot Shale and formation water from Paleozoic petroleum systems,Saudi Arabia

Provenance studies of produced water are essential to trace flow dynamics and reservoir compartmentalization in petroleum systems and to quantify fluid recovery rates from unconventional fracturing. Produced water from a hydraulically fractured well in the Qusaiba Hot Shale in the Northern Exploration Area, Saudi Arabia, was daily monitored and analyzed for water chemistry, and environmental (δ²H, δ¹³C, δ¹⁸O_H2O, δ¹⁸O_SO4, δ³⁴S_SO4, δ³⁷Cl, ⁸⁷Sr/⁸⁶Sr) and cosmogenic isotopes (³H, ¹⁴C, ³⁶Cl), to differentiate from reference fluids of supply water, fracturing fluids, and formation water from adjacent Paleozoic units. Initially, recovered water is composed of fracturing fluids and subsequently replaced by a homogeneous cut of pristine formation water. Formation water is composed of dominant meteoric water (approximately 84 vol%) and minor fossil evaporated seawater. The young ¹⁴C‐apparent age between 6000 and 6700 years BP and depleted δ¹⁸O/δ²H values for the meteoric component confirm the infiltration of surface water into the Qusaiba Hot Shale interval or adjacent units during the Early Holocene Pluvial Period under cooler and wetter climatic conditions than present, which suggest the presence of a very recent, dynamic hydraulic flow system. ³⁶Cl/Cl ratios between 102 × 10^?15 and 31 × 10^?15 are ambiguous and can be attributed to atmospheric recharge close to the coast, mixing of ³⁶Cl‐enriched Quaternary meteoric recharge with ³⁶Cl‐depleted fossil seawater, and/or hypogene production by U‐Th‐enriched host rock. Produced waters from Qusaiba Hot Shale are within the compositional range of Na‐Cl‐type formation water from Paleozoic reservoir units in northern Saudi Arabia with salinities from 30 000 to 130 000 mg l^?1. As a novel technological approach for exploration wells in Northern Saudi Arabia, multi‐isotopic methods were successfully implemented to quantify flowback volumes from hydraulic fracturing, and to fingerprint pristine formation water or pore water in Paleozoic systems on their provenance, residence time, migration pathways, and secondary alteration processes.     

Hydraulic observations from a 1 year fluid production test in the 4000 m deep KTB pilot borehole

A long‐term pump test was conducted in the KTB pilot borehole (KTB‐VB), located in the Oberpfalz area, Germany. It produced 22 300 m³ of formation fluid. Initially, fluid production rate was 29 l min^?1 for 4 months, but was then raised to an average of 57 l min^?1 for eight more months. The aim of this study was to examine the fluid parameters and hydraulic properties of fractured, crystalline crusts as part of the new KTB programme ‘Energy and Fluid Transport in Continental Fault Systems’. KTB‐VB has an open‐hole section from 3850 to 4000 m depth that is in hydraulic contact with a prominent continental fault system in the area, called SE2. Salinity and temperature of the fluid inside the borehole, and consequently hydrostatic pressure, changed significantly throughout the test. Influence of these quantities on variations in fluid density had to be taken into account for interpretation of the pump test. Modelling of the pressure response related to the pumping was achieved assuming the validity of linear Darcy flow and permeability to be independent of the flow rate. Following the principle ‘minimum in model dimension’, we first examined whether the pressure response can be explained by an equivalent model where rock properties around the borehole are axially symmetric. Calculations show that the observed pressure data in KTB‐VB can in fact be reproduced through such a configuration. For the period of high pumping rate (57 l min^?1) and the following recovery phase, the resulting parameters are 2.4 × 10^?13 m³ in hydraulic transmissivity and 3.7 × 10^?9 m Pa^?1 in storativity for radial distances up to 187 m, and 4.7 × 10^?14 m³ and 6.0 × 10^?9 m Pa^?1, respectively, for radial distances between 187 and 1200 m. The former pair of values mainly reflect the hydraulic properties of the fault zone SE2. For a more realistic hydraulic study on a greater scale, program FEFLOW was used. Parameter values were obtained by matching the calculated induced pressure signal to fluid‐level variations observed in the KTB main hole (KTB‐HB) located at 200 m radial distance from KTB‐VB. KTB‐HB is uncased from 9031 to 9100 m and shows indications of leakage in the casing at depths 5200–5600 m. Analysis of the pressure record and hydraulic modelling suggest the existence of a weak hydraulic communication between the two boreholes, probably at depths around the leakage. Hydraulic modelling of a major slug‐test in KTB‐HB that was run during the pumping in KTB‐VB reveals the effective transmissivity of the connected formation to be 1 to 2 orders of magnitude lower than the one determined for the SE2 fault zone.     

Origin of deep saline groundwaters in the Vienne granitic rocks (France): constraints inferred from boron and strontium isotopes

As part of a preliminary geological characterization programme to assess the feasibility of an underground laboratory in granitic rock, a series of 17 deep boreholes (maximum depth, 900 m) was drilled by ANDRA in the Vienne district, France. A salinity gradient was demonstrated in the granitic waters with concentrations varying from approximately 1 g L^?1 at 150 m depth at the top of the basement (beneath the sedimentary cover) to 10 g L^?1 in the deeper part (from 400 to 600 m depth). Sr and B isotope ratios were measured in order to better understand the origin of the salinity and to evaluate the degree of water–rock interaction in the system. The results obtained were compared to those of mineral spring waters emerging from the granitic basement in the Massif Central. Evidence in support of a significant marine contribution include: (i) the Cl–Br investigations agree with a marine origin for the saline groundwaters without evolution from seawater; (ii) the ⁸⁷Sr/⁸⁶Sr ratio of the Vienne deep groundwaters (0.7078–0.7084) is in agreement with a palaeo‐seawater isotopic signature; (iii) measured δ¹¹B values for the deepest brine samples are enriched in ¹¹B (up to 36.1‰) relative to the granitic springs. The combined use of δ¹¹B, Cl, B, Br, Sr contents and ⁸⁷Sr/⁸⁶Sr ratios makes it possible to define and quantify a mixing model between marine and crustal end‐members in order to explain the origin of the deep saline groundwaters in the Vienne granitic rocks.     

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.     

Br/Cl signature of hydrothermal fluids: liquid–vapour fractionation of bromine revisited

Br/Cl ratios of hydrothermal fluids are widely used as geochemical tracers in marine hydrothermal systems to prove fluid phase separation processes. However, previous results of the liquid–vapour fractionation of bromine are ambiguous. Here we report new experimental results of the liquid–vapour fractionation of bromine in the system H₂O–NaCl–NaBr at 380–450°C and 22.9–41.7 MPa. Our data indicate that bromine is generally more enriched than chlorine in the liquid phase. Calculated exchange coefficients K_D(Br‐Cl)^{liquid‐vapour} for the reaction Br^vapour + Cl^liquid = Br^liquid + Cl^vapour are between 0.94 ± 0.08 and 1.66 ± 0.14 within the investigated P–T range. They correlate positively with D_Cl^{liquid‐vapour} and suggest increasing bromine–chlorine fractionation with increasing opening of the liquid–vapour solvus, i.e. increasing distance to the critical curve in the H₂O–NaCl system. An empirical fit of the form K_D(Br‐Cl)^{liquid‐vapour} = a*ln[b*(D_Cl^{liquid‐vapour}?1) + e^1/a] yields a = 0.349 and b = 1.697. Based on this empirical fit and the well‐constrained phase relations in the H₂O–NaCl system we calculated the effect of fluid phase separation on the Br/Cl signature of a hydrothermal fluid with initial seawater composition for closed and open adiabatic ascents along the 4.5 and 4.8 J g^?1 K^?1 isentropes. The calculations indicate that fluid phase separation can significantly alter the Br/Cl ratio in hydrothermal fluids. The predicted Br/Cl evolutions are in accord with the Br/Cl signatures in low‐salinity vent fluids from the 9 to 10°N East Pacific Rise.     

Temperature‐dependent Li isotope ratios in Appalachian Plateau and Gulf Coast Sedimentary Basin saline water

Lithium (Li) concentrations of produced water from unconventional (horizontally drilled and hydraulically fractured shale) and conventional gas wells in Devonian reservoirs in the Appalachian Plateau region of western Pennsylvania range from 0.6 to 17 mmol kg^?1, and Li isotope ratios, expressed as in δ⁷Li, range from +8.2 to +15‰. Li concentrations are as high as 40 mmol kg^?1 in produced waters from Plio‐Pleistocene through Jurassic‐aged reservoirs in the Gulf Coast Sedimentary Basin analyzed for this study, and δ⁷Li values range from about +4.2 to +16.6‰. Because of charge‐balance constraints and rock buffering, Li concentrations in saline waters from sedimentary basins throughout the world (including this study) are generally positively correlated with chloride (Cl), the dominant anion in these fluids. Li concentrations also vary with depth, although the extent of depth dependence differs among sedimentary basins. In general, Li concentrations are higher than expected from seawater or evaporation of seawater and therefore require water–mineral reactions that remove lithium from the minerals. Li isotope ratios in these produced waters vary inversely with temperature. However, calculations of temperature‐dependent fractionation of δ⁷Li between average shale δ⁷Li (?0.7‰) and water result in δ⁷Li_water that is more positive than that of most produced waters. This suggests that aqueous δ⁷Li may reflect transport of water from depth and/or reaction with rocks having δ⁷Li lighter than average shale.     

Solubility of corundum + kyanite in H2O at 700°C and 10 kbar: evidence for Al‐Si complexing at high pressure and temperature

The solubility of the assemblage corundum + kyanite in H₂O was determined at 700°C and 10 kbar, using a piston‐cylinder apparatus and rapid‐quench/fluid‐extraction techniques. Weighted mean concentrations of total Al and Si were 5.80 ± 0.03 mmol kg^?1 H₂O and 0.308 ± 0.003 mol kg^?1 H₂O, respectively (1σ errors). The Al concentration is nearly five times higher than that of corundum solubility in pure H₂O. This difference is interpreted to arise from complexing between Si and Al to form HAlSiO_4,aq species. Charged or more polymerized species are also possible, but their abundance cannot be constrained based on these experiments. Assumption of a single aqueous aluminosilicate complex permits calculation of the thermodynamic consequences of Al‐Si interaction in high‐pressure fluids, as well as phase diagrams for the system Al‐Si‐O‐H. Formation of Al‐Si complexes leads to a large increase in dissolved Al with increasing Si in solution, such that Al concentration in equilibrium with kyanite + quartz is predicted to be 7.1 mmolal, higher than with kyanite + corundum. The elevated concentration of Si in deep‐crustal and mantle aqueous fluids suggests that Al must readily be dissolved and transported by Al‐Si complexing in high‐pressure metamorphic and metasomatic environments. The results provide a simple explanation for the common observation of kyanite + quartz segregations in eclogites and Barrovian metamorphic rocks.     

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.     

Constraining groundwater flow in the glacial drift and saginaw aquifers in the Michigan Basin through helium concentrations and isotopic ratios

³He and ⁴He concentrations in excess of those in water in solubility equilibrium with the atmosphere by up to two and three orders of magnitude are observed in the shallow Glacial Drift and Saginaw aquifers in the Michigan Basin. A simplified He transport model shows that in situ production is negligible and that most He excesses have a source external to the aquifer. Simulated results show that ³He and ⁴He fluxes entering the bottom of the Saginaw aquifer are 7.5 × 10^?14 and 6.1 × 10^?7 cm³STPcm^?2 yr^?1, both of which are lower than fluxes entering the underlying Marshall aquifer, 1.0 × 10^?13 and 1.6 × 10^?6 cm³STPcm^?2 yr^?1 for ³He and ⁴He, respectively. In contrast, He fluxes entering the Saginaw aquifer are higher than fluxes entering the overlying Glacial Drift aquifer of 5.2 × 10^?14 and 1.5 × 10^?7 cm³STPcm^?2 yr^?1 for ³He and ⁴He, respectively. The unusually high He fluxes and their decreasing values from the lower Marshall to the upper Glacial Drift aquifer strongly suggest the presence of an upward cross‐formational flow, with increasing He dilution toward the surface by recharge water. These fluxes are either comparable to or far greater than He fluxes in deeper aquifers around the world. Model simulations also suggest an exponential decrease in the horizontal groundwater velocity with recharge distance. Horizontal velocities vary from 13 to 2 myr^?1 for the Saginaw aquifer and from 18 to 6 myr^?1 for the Marshall aquifer. The highly permeable Glacial Drift aquifer displays a greater velocity range, from 250 to 5 myr^?1. While Saginaw ⁴He ages estimated based on the simulated velocity field display an overall agreement with ¹⁴C ages, ¹⁴C and ⁴He ages in the Glacial Drift and Marshall aquifers deviate significantly, possibly due to simplifications introduced in the He transport model leading to calculation of first‐order approximation He ages and high uncertainties in Glacial Drift ¹⁴C ages.     

Fluid flow and the Heart Mountain fault: a stable isotopic, fluid inclusion, and geochronologic study

Numerous studies have proposed that movement along the 3400‐km² low‐angle (<3°) Heart Mountain detachment fault of north‐western Wyoming and south‐western Montana was facilitated by the presence of lubricating fluids. A recent stable isotope study suggested that the fluids along the Heart Mountain fault originated from large hydrothermal systems associated with Eocene intrusive centers. Herein, we present results from a combined stable isotope, fluid inclusion, and ⁴⁰Ar/³⁹Ar geochronology study of the relationship between shallow crustal fluids from Eocene intrusive centers in the New World Mining district and fluids associated with the Heart Mountain fault. Our results suggest that: (i) Eocene intrusive bodies within the New World Mining District focused hydrothermal fluids along specific units and structures from the ore deposit centers along the Heart Mountain fault detachment surface; (ii) hydrothermal waters were focused along the Heart Mountain fault from the breakaway region near Silvergate, Montana toward Heart Mountain in Cody, Wyoming; and (iii) hydrothermal activity along the Heart Mountain fault occurred at the same time as emplacement of the Eocene intrusives, between 50.1 and 48.1 Ma. These data, taken together, suggest that the hydrothermal activity generated by intrusion of Eocene rhyodacites and dacites provided the source of fluids found along the Heart Mountain fault plane.     

Sources and composition of fluids associated with fluorite deposits of Asturias (N Spain)

The fluorite deposits of Asturias (northern Iberian Peninsula) are hosted by rocks of Permo‐Triassic and Palaeozoic age. Fluid inclusions in ore and gangue minerals show homogenization temperatures from 80 to 170°C and the presence of two types of fluids: an H₂O–NaCl low‐salinity fluid (<8 eq. wt% NaCl) and an H₂O–NaCl–CaCl₂ fluid (7–13 wt% NaCl and 11–14 wt% CaCl₂). The low salinity and the Cl/Br and Na/Br ratios (Cl/Br_molar 100–700 and Na/Br_molar 20–700) are consistent with an evaporated sea water origin of this fluid. The other end‐member of the mixture was highly saline brine with high Cl/Br and Na/Br ratios (Cl/Br_molar 700–13 000 and Na/Br_molar 700–11 000) generated after dissolution of Triassic age evaporites. LA‐ICP‐MS analyses of fluid inclusions in fluorite reveal higher Zn, Pb and Ba contents in the high‐salinity fluids (160–500, 90–170, 320–480 p.p.m. respectively) than in the low‐salinity fluid (75–230, 25–150 and 100–300 p.p.m. respectively). The metal content of the fluids appears to decrease from E to W, from Berbes to La Collada and to Villabona. The source of F is probably related to leaching of volcanic rocks of Permian age. Brines circulated along faults into the Palaeozoic basement. Evaporated sea water was present in permeable rocks and faults along or above the unconformity between the Permo‐Triassic sediments and the Palaeozoic basement. Mineralization formed when the deep brines mixed with the surficial fluids in carbonates, breccias and fractures resulting in the formation of veins and stratabound bodies of fluorite, barite, calcite, dolomite and quartz and minor amounts of sulphides. Fluid movement and mineralization occurred between Late Triassic and Late Jurassic times, probably associated with rifting events related to the opening of the Atlantic Ocean. This model is also consistent with the geodynamic setting of other fluorite‐rich districts in Europe.     

Experimental study of aqueous fluid infiltration into quartzite: implications for the kinetics of fluid redistribution and grain growth driven by interfacial energy reduction

To investigate the kinetics of interfacial energy‐driven fluid infiltration, experiments were carried out in a quartzite–water system at 621–925°C and 0.8 GPa. Infiltration couples were made by juxtaposing presynthesized dry quartzite cylinders and fluid reservoirs. The infiltration process was confirmed by the presence of pores at the quartzite grain edges. As predicted from theoretical considerations and previous experiments, wetting fluids such as pure water and NaCl aqueous solution infiltrated into quartzite, whereas nonwetting CO₂‐rich fluids did not. Newly precipitated quartz layers at the surfaces of the infiltrated sample proved that infiltration took place by a dissolution–precipitation mechanism. The enhancement of grain growth by fluid infiltration was observed over the entire range of experimental temperatures. The fluid fraction, gauged by the porosity of the run products, increases at the infiltration front and then decreases towards the fluid reservoir to form a high‐porosity zone with a maximum porosity of 2.3–2.9%. As infiltration proceeds, the high‐porosity zone advances like a travelling wave. This porosity wave is probably caused by a grain curvature gradient resulting from preferential grain growth in the infiltrated part of the quartzite, perhaps combined with other factors. The infiltration kinetics were modelled with a steady‐state diffusion model over the high‐porosity zone. The solubility difference between dissolving and precipitating grains was deduced to be 2 × 10^?2?3 × 10^?1 wt %. The experimentally obtained infiltration rate of aqueous fluid in the steady‐state diffusion regime (2 ± 0.5 × 10^?8 m sec^?1 at 823°C) is much faster than the estimated metamorphic fluid flux rates, so that interfacial energy‐driven fluid redistribution in quartz‐rich layers could significantly contribute to the fluid flux in high‐grade metamorphism, at least over a short distance. Cathodoluminescence observations of the run products revealed that the grain growth of quartzite in the presence of fluid proceeds extensively, which would promote the chemical equilibration between fluid and rock more effectively than would volume diffusion in quartz crystals.     

Molecular transport of methane, ethane and nitrogen and the influence of diffusion on the chemical and isotopic composition of natural gas accumulations

Laboratory experiments have been performed to determine diffusion coefficients of natural gas components (methane, ethane and nitrogen) and isotope fractionation effects under simulated in situ pressure (up to 45 MPa effective stress) and temperature conditions (50–200°C) in water‐saturated pelitic and coarse‐grained rocks. Effective diffusion coefficients of molecular nitrogen (0.39 × 10^?11 to 21.6 × 10^?11 m² sec^?1 at 90°C) are higher than those for methane (0.18 × 10^?11 to 18.2 × 10^?11 m² sec^?1 at 90°C). Diffusive flux rates expressed in mass units are generally higher for N₂ than for CH₄. Both methane and (to a lesser extent) nitrogen diffusion coefficients decrease with increasing total organic carbon (TOC) content of the rock samples because of sorption processes on the organic matter. This effect decreases with increasing temperature. Effective diffusion coefficients increase upon a temperature increase from 50 to 200°C by a factor of four. Effective diffusion coefficients and steady‐state diffusive flux decrease with effective stress. Stationary diffusive fluxes drop by 50–70% for methane and 45–62% for nitrogen while effective diffusion coefficients are reduced by 38% (CH₄) and 32–48% (N₂), respectively. Isotope fractionation coefficients of diffusive transport are higher for methane (?1.56 and ?2.77‰) than for ethane (?0.84 and ?1.62‰). Application of the experimental results to geological systems show that diffusive transport has only a low transport efficiency. Significant depletion of natural gas reservoirs by molecular diffusion is only expected in cases of very poor caprock qualities (in terms of thickness and/or porosity) and over extended periods of geological time. Under these circumstances, the chemical and isotopic composition of a gas reservoir will change and maturity estimates based on these parameters may be deceptive. To account for these potential effects, nomograms have been developed to estimate diffusive losses and apply maturity corrections.     

Composition and evolution of fluids during skarn development in the Monte Capanne thermal aureole,Elba Island,central Italy

We describe the chemistry of the fluids circulating during skarn formation by focusing on fluids trapped in calcsilicate minerals of the inner thermal aureole of the Late Miocene Monte Capanne intrusion of western Elba Island (central Italy). Primary, CH₄‐dominant, C‐O‐H‐S‐salt fluid inclusions formed during prograde growth of the main skarn‐forming mineral phases: grossular/andradite and vesuvianite. The variable phase ratios attest to heterogeneous entrapment of fluid, with co‐entrapment of an immiscible hydrocarbon–brine mixture. Chemical elements driving skarn metasomatism such as Na, K, Ca, S and Cl, Fe and Mn were dominantly partitioned into the circulating fluid phase. The high salinity (apparent salinity between 58 and 70 wt% NaCl eq.) and the C‐component of the fluids are interpreted as evidence for a composite origin of the skarn‐forming fluids that involves both fluids derived from the crystallizing intrusion and contributions from metamorphic devolatilization. Oxidation of a Fe‐rich brine in an environment dominated by fluctuation in pressure from lithostatic to hydrostatic conditions (maintained by active crack‐sealing) contributed to skarn development. Fluid infiltration conformed to a geothermal gradient of about 100°C km^?1, embracing the transition from high‐temperature contact metamorphism and fluid‐assisted skarn formation (at ca 600°C) to a barren hydrothermal stage (at ca 200°C).     

Solubility of quartz in crustal fluids: experiments and general equations for salt solutions and H2O–CO2 mixtures at 400–800°C and 0.1–0.9 GPa

The solubility of quartz has been measured in a wide range of salt solutions at 800°C and 0.5 GPa, and in NaCl, CaCl₂ and CsCl solutions and H₂O–CO₂ fluids at six additional P–T conditions ranging from 400°C at 0.1 GPa to 800°C at 0.9 GPa. The experiments cover a wide range of compositions along each binary. At P–T conditions where the density of pure water is low (0.43 g cm^?3), addition of most salts produces an enhancement of quartz solubility at low to moderate salt concentrations (salt‐in effect), although quartz solubility falls with further decrease in X_H2O. At higher fluid densities (0.7 g cm^?3 and greater), the salt‐in effect is generally absent, although this depends on both the cation present and the actual P–T conditions. The salt‐in effect is most readily produced by chloride salts of large monovalent cations, while CaCl₂ only produced a salt‐in effect at the most extreme conditions of high‐T and low‐P investigated (800°C at 0.2 GPa). Under most crustal conditions, the addition of common salts to aqueous fluids results in a lowering of quartz solubility relative to that in pure water (salt‐out effect). Comparing quartz solubility in different fluids by calculating X_H2O on the basis that all salts are fully associated under all conditions yields higher quartz solubility in solutions of monovalent salts than in solutions of divalent salts, absolute values are also influenced by cation radius. Quartz solubility measurements have been fitted to a Setchenow‐type equation, modified to take account of the separate effects of both the lowering of X_H2O and the specific effects of different salts, which are treated as arising through distinct patterns of non‐ideal behaviour, rather than the explicit formation of additional silica complexes with salt components. Quartz solubility in H₂O–CO₂ fluids can be treated as ideal, if the solvation number of aqueous silica is taken as 3.5. For this system the solubility (molality) of quartz in the binary fluid, S is related to its solubility in pure water at the same P–T conditions, S^o, by: Quartz solubility in binary salt systems (H₂O–RCl_n) can be fitted to the relationship: where salt concentration mRCl_n is expressed as molality and the exponent b has a value of 1 except under conditions where salting‐in is observed at low salt concentrations, in which case it is <1. Under most crustal conditions, the solubility of quartz in NaCl solutions is given to a good approximation by: We propose that quartz solubility in multicomponent fluids can be estimated from an extended expression, calculating X_H2O based on the total fluid composition (including dissolved gasses), and adding terms for each major salt present. Our experimental results on H₂O–NaCl–CO₂ fluids are satisfactorily predicted on this basis. An important implication of the results presented here is that there are circumstances where the migration of a fluid from one quartz‐bearing host into another, if it is accompanied by re‐equilibration through cation exchange, may lead to dissolution or precipitation of quartz even at constant P and T, with concomitant modification of the permeability structure of the deep crust.     

Geological setting, nature of ore fluids and sulphur isotope geochemistry of the Fu Ning Carlin-type gold deposits, Yunnan Province, China

Carlin‐type gold deposits in southern China are present in Palaeozoic to Mesozoic siliciclastic and carbonate rocks. The border region of Yunnan, Guizhou and Guangxi Provinces contains gold deposits on the south‐western margin of the Pre‐Cambrian South China Craton in south‐eastern Yunnan Province. The Fu Ning gold deposits host epigenetic, micron‐sized disseminated gold in: (i) Middle Devonian (D₁p) black carbonaceous mudstone at the Kuzhubao gold deposit and (ii) fault breccia zones at the contact between Triassic gabbro (β ) and the Devonian mudstone (D₁p) at the Bashishan gold deposit. The deposits are associated with zones of intense deformation with enhanced permeability and porosity that focused hydrothermal fluid flow, especially where low‐angle N‐S striking thrust faults are cut by NW striking strike‐slip and/or NE striking normal faults. Major sulphide ore minerals in the Fu Ning gold deposits are pyrite, arsenopyrite, arsenic‐rich pyrite, stibnite and minor iron‐poor sphalerite. Gangue minerals are quartz, sericite, calcite, ankerite and chlorite. Hypogene ore grades range from 1 to 7 g t^?1 Au and up to 18 g t^?1 Au at the Kuzhubao gold deposit and are generally less than 3 g t^?1 Au at the Bashishan gold deposit. Sub‐microscopic gold mineralization is associated with finely disseminated arsenic‐rich pyrite in the Stage III mineral assemblage. Two types of primary fluid inclusions have been recorded: Type I liquid–vapour inclusions with moderate‐to‐high liquid/vapour ratios, and Type II inclusions containing moderate liquid/vapour ratios with CO₂ as determined from laser Raman analysis. Temperature of homogenization (T_h) data collected from these primary fluid inclusions in gold‐ore Stage III quartz ranged from 180 to 275°C at the Kuzhubao gold deposit and 210 to 330°C at the Bashishan gold deposit. Salinity results indicate that there were possibly two fluids present during gold deposition, including: (i) an early fluid with 0.8–6.5 wt.% NaCl equivalent, similar to salinity in shear‐zone‐hosted gold deposits with metamorphic derived fluids; and (ii) a late fluid with 11.8–13.4 wt.% NaCl equivalent, indicating possible derivation from connate waters and/or brine sources. CO₂ and trace CH₄ were only detected by laser Raman spectrometry in gold‐ore‐stage primary fluid inclusions. Results of sulphur isotope studies showed that δ³⁴S values for pyrite and arsenopyrite associated with gold‐ore mineralization during Stage III at the Kuzhubao and Bashishan gold deposits are isotopically similar and moderately heavy with a range from +9 to +15 per mil, and also fall into the range of δ³⁴S values reported for Carlin‐type gold deposits. Sulphur isotopes suggest that the Fu Ning gold deposits were formed from connate waters and/or basinal brines. Fluid geochemistry data from the Fu Ning gold deposits suggest a Carlin‐type genetic model, involving fluid mixing between: (i) deep CO₂‐rich metamorphic fluids, (ii) moderately saline, reduced connate waters and/or basinal brines; and (iii) evolved meteoric waters.     

Strontium isotopic characterization of the Palmottu hydrosystem (Finland): water–rock interaction and geochemistry of groundwaters

The Palmottu hydrosystem is located in a granitic host rock in southern Finland. Along well‐defined pathways in the fractured crystalline rock, strontium isotopes are used to trace the degree of water–rock interaction (WRI) and mixing processes in groundwaters. The ⁸⁷Sr/⁸⁶Sr ratios range between 0.716910 and 0.735606 in the surface waters and between 0.719991 and 0.750787 in the groundwaters, but are between 0.720 and 0.735 in most of the samples. Moreover, the results show a lack of correlation between the water chemistries determining the classification into different water‐types (Na–Cl, Na–SO₄, etc.) and the results of the strontium (Sr) contents and Sr isotopic ratios. From a WRI standpoint, this implies that the Sr behaviour is independent of the water chemistry; the occurrence of large ⁸⁷Sr/⁸⁶Sr variations is site specific and mainly dependent on the lithology. A model to determine the ⁸⁷Sr/⁸⁶Sr ratio of water after interaction with granite was developed. This model is based on the assumption that Sr was derived from three minerals: plagioclase, K‐feldspar and biotite. The results of the calculation indicate that around half of the water analysed within the Palmottu hydrosystem can be explained by the weathering of the granites. However, clearly lower ⁸⁷Sr/⁸⁶Sr are observed in waters when compared to the calculated ⁸⁷Sr/⁸⁶Sr and other sources of Sr, with low ⁸⁷Sr/⁸⁶Sr, rather than the calculated granite–water interaction, which may be suspected. When comparing the ⁸⁷Sr/⁸⁶Sr and ion ratios (Ca/Na, Mg/Na, Sr/Na, Cl/Na), the scattering of the data can be explained by the presence of four end‐members: a brine component (low ⁸⁷Sr/⁸⁶Sr and Ca/Na ratios…), a deep granitic component (high ⁸⁷Sr/⁸⁶Sr ratios and low Ca/Na ratios…), a subsurface component (intermediate ⁸⁷Sr/⁸⁶Sr ratios associated with high Ca/Na ratios…) and a surface end‐member:snow and river drainage (low ⁸⁷Sr/⁸⁶Sr and low Ca/Na ratios…). These extreme end‐members define a series of WRI‐mixing line within a rather complex hydrosystem.     

Effective gas permeability of Tight Gas Sandstones as a function of capillary pressure – a non‐steady‐state approach

Single‐ and two‐phase (gas/water) fluid transport in tight sandstones has been studied in a series of permeability tests on core plugs of nine tight sandstones of the southern North Sea. Absolute (Klinkenberg‐corrected) gas permeability coefficients (k_gas__inf) ranged between 3.8 × 10^?16 and 6.2 × 10^?19 m² and decreased with increasing confining pressure (10–30 MPa) by a factor 3–5. Klinkenberg‐corrected (intrinsic) gas permeability coefficients were consistently higher by factors from 1.4 to 10 than permeability coefficients determined with water. Non‐steady‐state two‐phase (He/water) flow experiments conducted up to differential pressures of 10 MPa document the dynamically changing conductivity for the gas phase, which is primarily capillary‐controlled (drainage and imbibition). Effective gas permeability coefficients in the two‐phase flow tests ranged between 1.1 × 10^?17 and 2.5 × 10^?22 m², corresponding to relative gas permeabilities of 0.03% and 10%. In the early phase of the nonstationary flow regime (before establishment of steady‐state conditions), they may be substantially (>50%) lower. Effective gas permeability measurements are affected by the following factors: (i) Capillary‐controlled drainage/imbibition, (ii) viscous–dynamic effects (iii) and slip flow.     

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.