Low‐temperature catalytic CO2 hydrogenation with geological quantities of ruthenium: a possible abiotic CH4 source in chromitite‐rich serpentinized rocks

Metal‐catalysed CO₂ hydrogenation is considered a source of methane in serpentinized (hydrated) igneous rocks and a fundamental abiotic process germane to the origin of life. Iron, nickel, chromium and cobalt are the catalysts typically employed in hydrothermal simulation experiments to obtain methane at temperatures >200°C. However, land‐based present‐day serpentinization and abiotic gas apparently develop below 100°C, down to approximately 40–50°C. Here, we document considerable methane production in thirteen CO₂ hydrogenation experiments performed in a closed dry system, from 20 to 90°C and atmospheric pressure, over 0.9–122 days, using concentrations of non‐pretreated ruthenium equivalent to those occurring in chromitites in ophiolites or igneous complexes (from 0.4 to 76 mg of Ru, equivalent to the amount occurring approximately in 0.4–760 kg of chromitite). Methane production increased with time and temperature, reaching approximately 87 mg CH₄ per gram of Ru after 30 days (2.9 mgCH₄/g_ru/day) at 90°C. At room temperature, CH₄ production rate was approximately three orders of magnitude lower (0.003 mgCH₄/g_ru/day). We report the first stable carbon and hydrogen isotope ratios of abiotic CH₄ generated below 100°C. Using initial δ¹³C_CO2 of ‐40‰, we obtained room temperature δ¹³C_CH4 values as ¹³C depleted as ?142‰. With time and temperature, the C‐isotope separation between CO₂ and CH₄ decreased significantly and the final δ¹³C_CH4 values approached that of initial δ¹³C_CO2. The presence of minor amounts of C₂‐C₆ hydrocarbons is consistent with observations in natural settings. Comparative experiments at the same temperatures with iron and nichel catalysts did not generate CH₄. Ru‐enriched chromitites could potentially generate methane at low temperatures on Earth and on other planets.     

Constraining the rate and extent of mantle serpentinization from seismic and petrological data: implications for chemosynthesis and tectonic processes

We used seismic velocity as a proxy for serpentinization of the mantle, which occurred beneath thinned but laterally continuous continental crust during continental break up, prior to opening of the Atlantic Ocean. The serpentinized sub‐continental mantle is now exhumed, beneath the Iberia Abyssal Plain and was accessed by scientific drilling on Ocean Drilling Program legs 149 and 173. Chromatographic modelling of kinetically limited transport of the serpentinization front yields a front displacement of 2197 ± 89 m, a time‐integrated fluid flux of 1098 ± 45 m³ m^?2 and a Damköhler number of 6.0 ± 0.2. Whether either surface reaction or chemical transport limit the rate of reaction, we calculate timescales for serpentinization of approximately 10⁵–10⁶ years. This yields time‐average fluid flux rates for H₂O, entering and reacting with the mantle, of 60–600 mol m^?2 a^?1 and for CH₄, produced as a by‐product of oxidation of Fe⁺⁺ to magnetite and exiting the mantle, of 0.55–5.5 mol m^?2 a^?1. This equates to a CH₄‐flux of 0.18–1.8 Tg a^?1 for coeval serpentinization of the mantle that was exhumed west of Iberia. This represents 0.03–0.3% of the present‐day annual CH₄‐flux from all sources and a higher fraction of pre‐anthropogenic (lower) CH₄ levels. CH₄ released by serpentinization at or beneath the seafloor could provide substrate for biological chemosynthesis and/or promote gas‐hydrate formation. Finally, noting its volumetric extent and rapidity (<10⁶ years), we interpret serpentinization to be a reckonable component of tectonic processes, contributing both diapiric and expansional forces and helping to ‘lubricate’ extensional processes. Given its anisotropic permeability, actively deforming serpentinite might impede melt migration which may be of interest, given the apparent lack of melt in some rifted margins.     

Migration of metamorphic CO2 into a coal seam: a natural analog study to assess the long‐term fate of CO2 in Coal Bed Carbon Capture,Utilization, and Storage Projects

This study assesses the displacement of coalbed methane by CO₂ migration along a fault into the coal seam in the Yaojie coalfield. Coal and gas samples were collected continuously at various distances in NO.2 coal seam from F19 fault. Vitrinite reflectance, maceral, and pore distributions and proximate analysis of fourteen coal samples were performed. Gas components, concentrations, carbon isotopes of 28 gas samples were determined. We examined the coal–gas trace characteristics of coalbed methane displaced away from the fault by CO₂ injection after geological ages. From east to west, away from the F19 fault, the CO₂ concentration decreased, whereas the CH₄ concentration increased gradually. The δ¹³C values for CO₂ varied between ?9.94‰ and 1.12‰, suggesting a metamorphic origin. A wider range of values (from ?9.94‰ to 20‰) was associated with the mixing of microbial carbon dioxide, isotopic fractionation during CO₂ migration through the microporous structures of coals, and/or carbon isotope fractionation during gas–water exchange and dissolution of CO₂. Away from the F19 fault, the volumes of micropores, mesopores and macropores decrease gradually. The Dubinin–Radushkevich (DR) micropore volume decreased from 0.0059 to 0.0037 cm³ g^‐1, and the mesopore and macropore volumes decreased from 0.066 to 0.026 cm³ g^‐1. The CO₂ injection can mobilize aromatic hydrocarbons and mineral matter from coal matrix, resulting in the decrease in the absorption peak intensity for coal samples after supercritical CO₂ treatment, which indicates that chemical reactions occur between coal and CO₂, not only physical adsorption.     

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     

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.     

The effect of deformation on permeability development in anhydrite and implications for caprock integrity during geological storage of CO2

Geological storage of CO₂ in depleted oil and gas reservoirs is one of the most promising options to reduce atmospheric CO₂ concentrations. Of great importance to CO₂ mitigation strategies is maintaining caprock integrity. Worldwide many current injection sites and potential storage sites are overlain by anhydrite‐bearing seal formations. However, little is known about the magnitude of the permeability change accompanying dilatation and failure of anhydrite under reservoir conditions. To this extent, we have performed triaxial compression experiments together with argon gas permeability measurements on Zechstein anhydrite, which caps many potential CO₂ storage sites in the Netherlands. Our experiments were performed at room temperature at confining pressures of 3.5–25 MPa. We observed a transition from brittle to semi‐brittle behaviour over the experimental range, and peak strength could be described by a Mogi‐type failure envelope. Dynamic permeability measurements showed a change from ‘impermeable’ (<10^?21 m²) to permeable (10^?16 to 10^?19 m²) as a result of mechanical damage. The onset of measurable permeability was associated with an increase in the rate of dilatation at low pressures (3.5–5 MPa), and with the turning point from compaction to dilatation in the volumetric versus axial strain curve at higher pressures (10–25 MPa). Sample permeability was largely controlled by the permeability of the shear faults developed. Static, postfailure permeability decreased with increasing effective mean stress. Our results demonstrated that caprock integrity will not be compromised by mechanical damage and permeability development. Geofluids (2010) 10 , 369–387     

Permeability of the Mercia Mudstone: suitability as caprock to carbon capture and storage sites

The Upper Triassic Mercia Mudstone is the caprock to potential carbon capture and storage (CCS) sites in porous and permeable Lower Triassic Sherwood Sandstone reservoirs and aquifers in the UK (primarily offshore). This study presents direct measurements of vertical (k_v) and horizontal (k_h) permeability of core samples from the Mercia Mudstone across a range of effective stress conditions to test their caprock quality and to assess how they will respond to changing effective stress conditions that may occur during CO₂ injection and storage. The Mercia samples analysed were either clay‐rich (muddy) siltstones or relatively clean siltstones cemented by carbonate and gypsum. Porosity is fairly uniform (between 7.4 and 10.7%). Porosity is low either due to abundant depositional illite or abundant diagenetic carbonate and gypsum cements. Permeability values are as low as 10^?20 m² (10nD), and therefore, the Mercia has high sealing capacity. These rocks have similar horizontal and vertical permeabilities with the highest k_h/k_v ratio of 2.03 but an upscaled k_h/k_v ratio is 39, using the arithmetic mean of k_h and the harmonic mean of k_v. Permeability is inversely related to the illite clay content; the most clay‐rich (illite‐rich) samples represent very good caprock quality; the cleaner Mercia Mudstone samples, with pore‐filling carbonate and gypsum cements, represent fair to good caprock quality. Pressure sensitivity of permeability increases with increasing clay mineral content. As pore pressure increases during CO₂ injection, the permeability of the most clay‐rich rocks will increase more than carbonate‐ and gypsum‐rich rocks, thus decreasing permeability heterogeneity. The best quality Mercia Mudstone caprock is probably not geochemically sensitive to CO₂ injection as illite, the cause of the lowest permeability, is relatively stable in the presence of CO₂–water mixtures.     

Gas breakthrough experiments on pelitic rocks: comparative study with N2, CO2 and CH4

The capillary‐sealing efficiency of intermediate‐ to low‐permeable sedimentary rocks has been investigated by N₂, CO₂ and CH₄ breakthrough experiments on initially fully water‐saturated rocks of different lithological compositions. Differential gas pressures up to 20 MPa were imposed across samples of 10–20 mm thickness, and the decline of the differential pressures was monitored over time. Absolute (single‐phase) permeability coefficients (k_abs), determined by steady‐state fluid flow tests, ranged between 10^?22 and 10^?15 m². Maximum effective permeabilities to the gas phase k_eff(max), measured after gas breakthrough at maximum gas saturation, extended from 10^?26 to 10^?18 m². Because of re‐imbibition of water into the interconnected gas‐conducting pore system, the effective permeability to the gas phase decreases with decreasing differential (capillary) pressure. At the end of the breakthrough experiments, a residual pressure difference persists, indicating the shut‐off of the gas‐conducting pore system. These pressures, referred to as the ‘minimum capillary displacement pressures’ (P_d), ranged from 0.1 up to 6.7 MPa. Correlations were established between (i) absolute and effective permeability coefficients and (ii) effective or absolute permeability and capillary displacement pressure. Results indicate systematic differences in gas breakthrough behaviour of N₂, CO₂ and CH₄, reflecting differences in wettability and interfacial tension. Additionally, a simple dynamic model for gas leakage through a capillary seal is presented, taking into account the variation of effective permeability as a function of buoyancy pressure exerted by a gas column underneath the seal.     

Low‐temperature dunite hydration: evaluating CH4 and H2 production from H2O and CO2

Abiotic methane (CH₄) and hydrogen (H₂) produced after hydration of mafic/ultramafic rocks represent energy sources for microbes that may thrive in the deep subsurface regions of Earth and possibly on other planets. While H₂ is a direct product of serpentinization, CH₄ can form via Fischer–Tropsch Type (FTT) reactions (carbon reduction) that, due to potential H₂ migration, can be spatially and temporally detached from serpentinization. We tested an alternative process hypothesized by some scholars, in which CO₂ can be reduced through dunite hydration without initially added H₂, implying that CH₄ can form in the same serpentinized fluid–rock system. The experiment used natural dunite sand (Forsterite 92), CO₂ with δ¹³C ~ ?25‰ (VPDB), and a 1 mm dissolved SiO₂ solution mixed in 30 glass bottles (118 mL) stored for up to 8 months at low temperature (50°C) to simulate land‐based serpentinization systems. In addition, 30 control bottles without olivine were used as blanks. Trivial amounts of CH₄ (orders of 0.2–0.9 ppmv) were detected in both samples and blanks, likely representing analytical noise; essentially, no significant amount of CH₄ formed under the experimental conditions used in this work. Low amounts of H₂ (~2.55 ± 1.39 ppmv) were generated, with production yields that were one order of magnitude lower than in previously published experiments. Moderate concentrations of SiO₂ appeared to hinder low‐temperature H₂ production. Our experiment confirms that the low‐temperature reduction of CO₂ into CH₄ through direct olivine hydration, without initial H₂, is sluggish and not straightforward, which is consistent with previous studies. The presence of substantial amounts of H₂, as well as suitable metal catalysts, appears to be essential in the low‐temperature production of abiotic CH₄, as observed in published FTT experiments.     

Geochemistry and origin of natural gases dissolved in brines from gas fields in southwest Japan

Previous geochemical studies indicated that most natural gases dissolved in brines in Japan are of microbial origin, consisting of methane produced via carbonate reduction. However, some of those from gas fields in southwest Japan contain methane relatively enriched in ¹³C, whose origin remains to be clarified. To address this issue, chemical and isotopic analyses were performed on natural gases and brines from the gas fields in Miyazaki and Shizuoka prefectures, southwest Japan. Methane isotopic signatures (δ¹³C ≈ ?68‰ to ?34‰ VPDB; δ²H ≈ ?183‰ to ?149‰ VSMOW) suggest that these gases are of microbial (formed via carbonate reduction) or of mixed microbial and thermogenic origin. The relatively high δ²H‐CH₄ values and their relationship with the δ²H‐H₂O values argue against the possibility of their formation via acetate fermentation. The δ¹³C‐CO₂ values (≈?5‰), together with the slope of the correlation between δ²H‐CH₄ and δ¹³C‐CH₄ (Δδ²H‐CH₄/Δδ¹³C‐CH₄ ≈ 1), contradict the possibility of their formation via carbonate reduction followed by partial oxidation by methanotrophs. The ³He/⁴He ratios of the gases from Miyazaki (≈0.11–1.3 Ra) and their low correlation with δ¹³C‐CH₄ values do not support an abiogenic origin. It is inferred therefore that the high δ¹³C‐CH₄ values of natural gases dissolved in brines from gas fields in southwest Japan are indications of the contribution of thermogenic hydrocarbons, although whether abiogenic hydrocarbons contribute significantly to the gases from Shizuoka requires further investigation. This study has clarified that, for the future exploration of natural gases in southwest Japan, we should adopt the strategies for conventional thermogenic gas accumulations, such as checking the content, type and maturity of organic matter in the underlying sedimentary rocks.     

Carbon dioxide controlled earthquake distribution pattern in the NW Bohemian swarm earthquake region,western Eger Rift,Czech Republic – gas migration in the crystalline basement

The ascent of magmatic carbon dioxide in the western Eger (Oh?e) Rift is interlinked with the fault systems of the Variscian basement. In the Cheb Basin, the minimum CO₂ flux is about 160 m³ h^?1, with a diminishing trend towards the north and ceasing in the main epicentral area of the Northwest Bohemian swarm earthquakes. The ascending CO₂ forms Ca‐Mg‐HCO₃ type waters by leaching of cations from the fault planes and creates clay minerals, such as kaolinite, as alteration products on affected fault planes. These mineral reactions result in fault weakness and in hydraulically interconnected fault network. This leads to a decrease in the friction coefficient of the Coulomb failure stress (CFS) and to fault creep as stress build‐up cannot occur in the weak segments. At the transition zone in the north of the Cheb Basin, between areas of weak, fluid conductive faults and areas of locked faults with frictional strength, fluid pressure can increase resulting in stress build‐up. This can trigger strike‐slip swarm earthquakes. Fault creep or movements in weak segments may support a stress build‐up in the transition area by transmitting fluid pressure pulses. Additionally to fluid‐driven triggering models, it is important to consider that fluids ascending along faults are CO₂‐supersaturated thus intensifying the effect of fluid flow. The enforced flow of CO₂‐supersaturated fluids in the transitional zone from high to low permeability segments through narrowings triggers gas exsolution and may generate pressure fluctuations. Phase separation starts according to the phase behaviour of CO₂‐H₂O systems in the seismically active depths of NW Bohemia and may explain the vertical distribution of the seismicity. Changes in the size of the fluid transport channels in the fault systems caused, or superimposed, by fault movements, can produce fluid pressure increases or pulses, which are the precondition for triggering fluid‐induced swarm earthquakes.     

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.     

The Design and Development of a Closed Chamber for the in‐situ Quantification of Dryland Soil Carbon Dioxide Fluxes

This paper describes the design features and capabilities of a portable automated in‐situ closed chamber (ISCC) for the quantification of CO₂ fluxes in dryland soils where both photosynthetic and respiratory components may be associated with a cyanobacterial crust. The processes of CO₂ flux in dryland soils are briefly described in order to clarify the conditions that make quantification of these fluxes problematic. The instrumentation currently available for in‐situ soil CO₂ flux measurements is then reviewed demonstrating their inadequacies for the dryland environment. The ISCC described here is a member of the closed or enrichment class of soil respiration chambers. The ISCC, however, features an optical window possessing high (>90%) transmission in the photosynthetic active region (PAR) of the solar irradiance spectrum, permitting observations of photosynthesis. The ISCC possesses automatic venting and purging so that gaseous concentrations inside the chamber do not change from ambient sufficiently to significantly affect diffusion. The ISCC features both active and passive cooling employing internal solid‐state Peltier coolers and external aluminised Mylar respectively. This avoids severe disturbance of the microclimate within the chamber due to admission of high fluxes of PAR and permits in‐situ operation under a wide range of ambient field temperatures (~ ?5 to 40°C). Sensors internal to the chamber monitor temperature, relative humidity, irradiance and pressure. In this implementation the ISCC is coupled to a portable gas chromatograph (Agilent GC‐3000) to sample the chamber atmosphere. Indicative data for Kalahari Sand soils of Botswana are presented as an illustration of the general performance characteristics.     

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.     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

Study of coal gas wettability for CO2 storage and CH4 recovery

To quantify and rank gas wettability of coal as a key parameter affecting the extent of CO₂ sequestration in coal and CH₄ recovery from coal, we developed a contact angle measuring system based on a captive gas bubble technique. We used this system to study the gas wetting properties of an Australian coal from the Sydney Basin. Gas bubbles were generated and captivated beneath a coal sample within a distilled water‐filled (pH 5.7) pressurised cell. Because of the use of distilled water, and the continuous dissolution and shrinkage of the gas bubble in water during measurement, the contact angles measured correspond to a ‘transient receding’ contact angle. To take into account the mixed‐gas nature (CO₂, CH₄, and to a lesser extent N₂) of coal seam gas in the basin, we evaluated the relative wettability of coal by CH₄, CO₂ and N₂ gases in the presence of water. Measurements were taken at various pressures of up to 15 MPa for CH₄ and N₂, and up to 6 MPa for CO₂ at a constant temperature of 22°C. Overall, our results show that CO₂ wets coal more extensively than CH₄, which in turn wets coal slightly more than N₂. Moreover, the contact angle reduces as the pressure increases, and becomes < 90° at various pressures depending on the gas type. In other words, all three gases wet coal better than water under sufficiently high pressure.     

Hydrogeochemical modelling of CO2 equilibria and mass transfer induced by organic–inorganic interactions in siliciclastic petroleum reservoirs

Different feldspar types control complex hydrogeochemical processes in hydrocarbon‐bearing siliciclastic reservoirs, which have undergone different degrees of degradation. To test such processes generically, carbon dioxide equilibria and mass transfers induced by organic–inorganic interactions have been modelled for different hydrogeochemical scenarios. The approach is based on and compared with data from the Norwegian continental shelf ( Smith & Ehrenberg 1989 ) and assumes local thermodynamic equilibrium among solids and fluids. Equilibrating mineral assemblages (different feldspar types, quartz, kaolinite, calcite) are based on the primary reservoir composition. Equilibration and coupled mass transfer were triggered by the addition and reaction of different amounts of CO₂, CH₄ and H₂ (plus acetic acid) at temperatures between 50 and 95°C (323 and 368 K). These components occur in oil fields as products of anaerobic bacterial degradation, hydrolytic disproportionation of hydrocarbons and/or thermal maturation of kerogen. We apply two different computer codes and two different thermodynamic data bases to calculate the results. Reaction of 0.32–0.6 mol CO₂, 0.16–0.3 mol CH₄ and 0.8–1.5 mol H₂ with K‐feldspar, quartz, kaolinite and calcite in 1 l of pore water results in modelled values of 0.3–2.3 mol% CO₂ in a multicomponent gas phase that resembles measured data (0.2–1.5 mol%). Similar CO₂ contents result from acetic acid addition (CO₂, CH₄, H₂ + 0.016 mol CH₃COOH). Equilibration with albite or anorthite reduces the release of CO₂ into the multicomponent gas phase dramatically, by 1 or 4 orders of magnitude compared with the equilibration with K‐feldspar. Minor differences in the modelled CO₂ content (0.1–0.2 mol%) result from calculations with different computer codes if the same thermodynamic data base is applied. Relevant differences (up to 1.9 mol% CO₂) result from calculations using different thermodynamic data bases.     

Regional groundwater flow and interactions with deep fluids in western Apennine: the case of Narni‐Amelia chain (Central Italy)

The elemental fluxes and heat flow associated with large aquifer systems can be significant both at local and at regional scales. In fact, large amounts of heat transported by regional groundwater flow can affect the subsurface thermal regime, and the amount of matter discharged towards the surface by large spring systems can be significant relative to the elemental fluxes of surface waters. The Narni‐Amelia regional aquifer system (Central Italy) discharges more than 13 m³ sec^?1 of groundwater characterised by a slight thermal anomaly, high salinity and high pCO₂. During circulation in the regional aquifer, groundwater reacts with the host rocks (dolostones, limestones and evaporites) and mixes with deep CO₂‐rich fluids of mantle origin. These processes transfer large amounts of dissolved substances, in particular carbon dioxide, and a considerable amount of heat towards the surface. Because practically all the water circulating in the Narni‐Amelia system is discharged by few large springs (Stifone‐Montoro), the mass and energy balance of these springs can give a good estimation of the mass and heat transported from the entire system towards the surface. By means of a detailed mass and balance of the aquifer and considering the soil CO₂ fluxes measured from the main gas emission of the region, we computed a total CO₂ discharge of about 7.8 × 10⁹ mol a^?1 for the whole Narni‐Amelia system. Finally, considering the enthalpy difference between infiltrating water and water discharged by the springs, we computed an advective heat transfer related to groundwater flow of 410 ± 50 MW.     

FT-IR measurements of petroleum fluid inclusions: methane, n-alkanes and carbon dioxide quantitative analysis

A recent advancement in petroleum geochemistry is to model fossil oil composition using microthermometric and volumetric data acquired from individual fluid inclusion analysis. Fourier transform infrared (FT‐IR) microspectroscopy can record compositional information related to gas (CH₄ and CO₂) and alkane contents of petroleum inclusions. In this study, a quantitative procedure for FT‐IR microspectrometry has been developed to obtain, from individual fluid inclusions, mol percentage concentrations of methane, alkanes and carbon dioxide as constraints to thermodynamic modelling. A petroleum inclusion in a sample from the Québec City Promontory nappe area was used as standard to record a reference spectrum of methane. The analytical procedure is based on the measurement of CH₄/alkane and CH₄/CO₂ band area ratios. CH₄/alkane infrared band area ratio is obtained after spectral subtraction of the reference methane spectrum. This area ratio, affected by absolute absorption intensities of methane, methyl and methylene, provides a molar CH₄/alkane ratio. Methyl/methylene ratio (CH₂/CH₃) ratio is obtained following procedures established in previous work. CO₂/CH₄ concentration ratio is estimated from relative absolute absorption intensities. Application to natural inclusions from different environments shows good correlation between FT‐IR quantification and PIT (petroleum inclusion thermodynamic) modelling.     

Petroleum type determination through homogenization temperature and vapour volume fraction measurements in fluid inclusions

Physical parameters of petroleum‐bearing fluid inclusions such as bulk density (ρ), molar volume (V_m), vapour volume fraction (?_vap) and homogenization temperature (T_h) are essential information to model petroleum composition (x) in inclusions and to reconstruct palaeotemperature and palaeopressure of trapping. For the main petroleum types contained in a fluid inclusion, we can follow how ?_vap and T_h are simultaneously influenced by a change of bulk density in a ?_vap versus T_h projection. We have correlated T_h and ?_vap for different petroleum compositions for a large range of bulk density values. However, postentrapment events under new pressure (P) and temperature (T) conditions can greatly modify the initial fingerprints of physical conditions and chemical composition of fluid inclusions. Re‐equilibration is frequent, especially in the case of fragile minerals. Stretching and leakage phenomenon have been simulated using the Petroleum Inclusion Thermodynamics (pit ) software, from virtual petroleum inclusions with known hydrocarbon composition. The aim of these simulations is to understand how ?_vap and T_h evolve with these re‐equilibration phenomena, with respect to the oil composition. Results of stretching simulations show a characteristic increase of T_h and ?_vap along correlation curves, with the curve shape dependent on petroleum composition. Leakage simulations show an increase of T_h and a smaller increase or even a decrease in ?_vap. Consequently, the better preserved inclusions in a given population can be presumed to be those that have the lowest T_h. Applications of T_h and ?_vap measurements of natural inclusions in calcite and in quartz showed that the fragility of the host mineral is a key factor allowing the recording of post‐entrapment events. Inclusions that have stretched or leaked are identified and the best preserved inclusions selected for evaluation of P–T–x trapping conditions. Moreover, petroleum types trapped in inclusions can be identified from ?_vap and T_h measurements without compositional modelling.