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.     

Post‐CO2 injection alteration of the pore network and intrinsic permeability tensor for a Permo‐Triassic sandstone

The aim of this study was to determine the process–structure–property relationships between the pre‐ and post‐CO₂ injection pore network geometry and the intrinsic permeability tensor for samples of core from low‐permeability Lower Triassic Sherwood Sandstone, UK. Samples were characterised using SEM‐EDS, XRD, MIP, XRCT and a triaxial permeability cell both before and after a three‐month continuous‐flow experiment using acidic CO₂‐rich saline fluid. The change in flow properties was compared to those predicted by pore‐scale numerical modelling using an implicit finite volume solution to the Navier–Stokes equations. Mass loss and increased secondary porosity appeared to occur primarily due to dissolution of intergranular cements and K‐feldspar grains, with some associated loss of clay, carbonate and mudstone clasts. This resulted in a bulk porosity increase from 18 to 25% and caused a reduction in mean diameter of mineral grains with an increase in apparent pore wall roughness, where the fractal dimension, D_f, increased from 1.68 to 1.84. All significant dissolution mass loss occurred in pores above c. 100 μm mean diameter. Relative dilation of post‐treatment pore area appeared to increase in relation to initial pore area, suggesting that the rate of dissolution mass loss had a positive relationship with fluid flow velocity; that is, critical flow pathways are preferentially widened. Variation in packing density within sedimentary planes (occurring at cm‐scale along the ‐z plane) caused the intrinsic permeability tensor to vary by more than a factor of ten. The bulk permeability tensor is anisotropic having almost equal value in ‐z and ‐y planes but with a 68% higher value in the ‐x plane (parallel to sedimentary bedding planes) for the pretreated sample, reducing to only 30% higher for the post‐treated sample. The intrinsic permeability of the post‐treatment sample increased by one order of magnitude and showed very close agreement between the modelled and experimental results.     

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.     

Modeling CO2 generation,migration, and titration in sedimentary basins

High mole fraction CO₂ gases pose a significant risk to hydrocarbon exploration in some areas. The generation and movement of CO₂ are also of scientific interest, particularly because CO₂ is an important greenhouse gas. We have developed a model of CO₂ generation, migration, and titration in basins in which a high mole fraction CO₂ gas is generated by the breakdown of siderite (FeCO₃) and magnesite (MgCO₃) where parts of the basin are being heated above approximately 330°C. The CO₂ reacts with Fe‐, Mg‐, and Ca‐silicates as it migrates upward and away from the generation zone (CO₂‐kitchen). Near the kitchen, where the Fe‐, Mg‐, and Ca‐silicates have been titrated and destroyed by previous packets of migrating CO₂, gas moves upward without lowering its CO₂ mole fraction. Further on, where Fe‐ and Mg‐silicates are still present but Ca‐silicates are absent in the sediments, the partial pressure of CO₂ is constrained to 0.1–30 bars and reservoirs contain a few mole percent CO₂ as described by Smith & Ehrenberg (1989) . Still further from the source, where Ca‐silicates have not been titrated, partial pressure of CO₂ in migrating methane gas are orders of magnitude lower. A 2D numerical model of CO₂ generation, migration, and titration quantifies these buffer relations and makes predictions of CO₂ risk in the South China Sea that are compatible with exploration experience. Reactive CO₂ transport models of the kind described could prove useful in determining how gases migrate in faulted sedimentary basins.     

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.     

Groundwater discharges in the Baltic Sea: survey and quantification using a schlieren technique application

Groundwater seeps are known to occur in Eckernförde Bay, Baltic Sea. Their discharge rate and dispersion were investigated with a new schlieren technique application, which is able to visualize heterogeneous water parcels with density anomalies down to Δσ_t = 0.049 on the scale of millimeters. With the use of an inverted funnel, discharged fluids can be captured and the outflow velocity can be determined. Overall, 46 stations could be categorized by three different cases: active vent sites, seep‐influenced sites, and non‐seep sites. New seep locations were discovered, even at shallow near‐shore sites, lacking prominent sediment depression, which indicate submarine springs. The detection of numerous seeps was possible and the groundwater‐influenced area was defined to be approximately 6.3 km². Flow rates of between 0.05 and 0.71 l m^?2 min^?1 were measured. A single focused fluid plume, which was not disturbed by the funnel was recorded and revealed a flux of 59.6 ± 20 ml cm^?2 min^?1 and it was calculated that this single focused plume would be strong enough to produce a flow rate through the funnel of 1.32 ± 0.44 l m^?2 min^?1. The effect of different seep‐meter funnel sizes is discussed.     

Bleached mudstone,iron concretions,and calcite veins: a natural analogue for the effects of reducing CO2‐bearing fluids on migration and mineralization of iron,sealing properties,and composition of mudstone cap rocks

This study investigates the Wangfu Depression of the Songliao Basin, China, as a natural analogue site for Fe migration (bleaching) and mineralization (formation of iron concretions) caused by reducing CO₂‐bearing fluids that leak along fractures after carbon capture, utilization, and storage. We also examined the origin of fracture‐filling calcite veins, the properties of self‐sealing fluids, the influence of fluids on the compositions of mudstone and established a bleaching model for the study area. Our results show that iron concretions are the oxidative products of precursor minerals (pyrite and siderite) during uplift and are linked to H₂S and CO₂ present in early stage fluids. The precipitation of calcite veins is the result of CO₂ degassing and is related to CO₂, CH₄, and minor heavy hydrocarbons in the main bleaching fluids. In our model, fluids preferentially enter high‐permeability fracture systems and result in the bleaching of surrounding rocks and precipitation of calcite veins. The infilling of calcite veins significantly decreases the permeability of fractures and forces the fluids to slowly enter and bleach the mudstone rocks. The Fe²⁺ released during bleaching migrates to elsewhere with the solutions or is reprecipitated in the calcite veins and iron concretions. The formation of calcite veins reduces the fracture space and effectively prevents fluid flow. The fluids have an insignificant effect on minerals within the mudstone. In terms of the chemistry of the mudstone, only the contents of Fe₂O₃, U, and Mo change significantly, with the content of U increasing in the mudstone and the contents of Fe₂O₃ and Mo decreasing during bleaching.     

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.     

Decreasing CO2 partial pressure triggered Mg–Fe–Ca carbonate formation in ancient Martian crust preserved in the ALH84001 Meteorite

We retrace hydrogeochemical processes leading to the formation of Mg–Fe–Ca carbonate concretions (first distinct carbonate population, FDCP) in Martian meteorite ALH84001 by generic hydrogeochemical equilibrium and mass transfer modeling. Our simple conceptual models assume isochemical equilibration of orthopyroxenite minerals with pure water at varying water‐to‐rock ratios, temperatures and CO₂ partial pressures. Modeled scenarios include CO₂ partial pressures ranging from 10.1325 to 0.0001 MPa at water‐to‐rock ratios between 4380 and 43.8 mol mol^?1 and different temperatures (278, 303 and 348 K) and enable the precipitation of Mg–Fe–Ca solid solution carbonate. Modeled range and trend of carbonate compositional variation from magnesio‐siderite (core) to magnesite (rim), and the precipitation of amorphous SiO₂ and magnetite coupled to magnesite‐rich carbonate are similar to measured compositional variation. The results of this study suggest that the early Martian subsurface had been exposed to a dynamic gas pressure regime with decreasing CO₂ partial pressure at low temperatures (approximately 1.0133 to 0.0001 MPa at 278 K or 6 to 0.0001 MPa at 303 K). Moderate water‐to‐rock ratios of ca. 438 mol mol^?1 and isochemical weathering of orthopyroxenite are additional key prerequisites for the formation of secondary phase assemblages similar to ALH84001’s ‘FDCP’. Outbursts of water and CO_2(g) from confined ground water in fractured orthopyroxenite rocks below an unstable CO₂ hydrate‐containing cryosphere provide adequate environments on the early Martian surface.     

The geochemical evolution of very dilute CO2‐rich water in Chungcheong Province,Korea: processes and pathways

A geochemical study was carried out on the CO₂‐rich water occurring in granite areas of Chungcheong Province, Korea. In this area, very dilute and acidic CO₂‐rich waters [62–242 mg l^?1 in total dissolved solid (TDS), 4.0–5.3 in pH; group I) occur together with normal CO₂‐rich waters (317–988 mg l^?1 in TDS, 5.5–6.0 in pH; group II). The concentration levels and ages of group I water are similar to those of recently recharged and low‐mineralized groundwater (group III). Calculation of reaction pathways suggests that group I waters are produced by direct influx of CO₂ gas into group III type waters. When the groundwater is injected with CO₂, it develops the capacity to accept dissolved solids and it can evolve into water with very high solute concentrations. Whether the water is open or closed to the CO₂ gases becomes less important in controlling the reaction pathway of the CO₂‐rich groundwater when the initial pco ₂ is high. Our data show that most of the solutes are dissolved in the CO₂‐rich groundwater at pH > 5 where the weathering rates of silicates are very slow or independent of pH. Thus, groundwater age is likely more important in developing high solute concentrations in the CO₂‐rich groundwaters than accelerated weathering kinetics because of acidic pH caused by high pco ₂.     

Chemical and isotopic signatures of Na/HCO3/CO2‐rich geofluids,North Portugal

Geochemical and isotopic studies have been undertaken to assess the origin of CO₂‐rich waters issuing in the northern part of Portugal. These solutions are hot (76°C) to cold (17°C) Na–HCO₃ mineral waters. The δ²H and δ¹⁸O signatures of the mineral waters reflect the influence of altitude on meteoric recharge. The lack of an ¹⁸O‐shift indicates there has been no high temperature water–rock interaction at depth, corroborating the results of several chemical geothermometers (reservoir temperature of about 120°C). The low ¹⁴C activity (up to 9.9 pmC) measured in some of the cold CO₂‐rich mineral waters (total dissolved inorganic carbon) is incompatible with the presence of ³H (from 1.7 to 4.1 TU) in those waters, which indicates relatively short subsurface circulation times. The δ¹³C values of CO₂ gas and dissolved inorganic carbon range between ?6‰ and ?1‰ versus Vienna‐Peedee Belemnite, indicating that the total carbon in the recharge waters is being diluted by larger quantities of CO₂ (¹⁴C‐free) introduced from deep‐seated (upper mantle) sources, masking the ¹⁴C‐dating values. The differences in the ⁸⁷Sr/⁸⁶Sr ratios of the studied thermal and mineral waters seem to be caused by water–rock interaction with different granitic rocks. Chlorine isotope signatures (?0.4‰ < δ³⁷Cl < +0.4‰ versus standard mean ocean chloride) indicate that Cl in these waters could be derived from mixing of a small amount of igneous Cl from leaching of granitic rocks.     

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     

Geodynamically induced variations in the emission of CO2 gas at San Faustino (Central Apennines,Italy)

The Central Apennines are affected by frequent earthquakes of moderate magnitude that occur mainly within the upper part of the crust at depths of <15 km. A large number of cold gas emissions that are rich in CO₂ are also found in the region. One particular vent with a high rate of degassing was equipped with a sensor to measure flow rates, which were recorded for a number of different periods between 2005 and 2010. Factors that could affect potentially CO₂ flow rates include barometric pressure, atmospheric temperature, precipitation and local seismicity. Our analysis indicates that the periods of anomalous flow rate were related not to the environmental factors but probably to the deformative processes of the crust associated with the local seismicity. Local seismic events as expression of geodynamic processes occurred always before and during these anomalous gas flow periods. This correlation exists only for events that occurred eastwards of the gas emission site close to the Martana fault zone. We herein consider this correlation as indication for a continuous interaction between the field of static strain and the deep fluid pressure. An approximation of the fluid pressure transmission towards the gas emission site gives reasonable values of 1–10 m² sec^?1. To make comparisons with the long‐term effects of the static strain, we also recorded the short‐term effects of the dynamic release of strain induced by the series of strong earthquakes that took place in L’Aquila in 2009. We detected a significant anomalous flow rate that occurred at the same time as this seismic sequence, during which widespread degassing was induced around the focal zone.     

Gas breakthrough experiments on fine-grained sedimentary rocks

The capillary sealing efficiency of fine‐grained sedimentary rocks has been investigated by gas breakthrough experiments on fully water saturated claystones and siltstones (Boom Clay from Belgium, Opalinus Clay from Switzerland and Tertiary mudstone from offshore Norway) of different lithological compositions. Sand contents of the samples were consistently below 12%, major clay minerals were illite and smectite. Porosities determined by mercury injection lay between 10 and 30% while specific surface areas determined by nitrogen adsorption (BET method) ranged from 20 to 48 m² g ^{? 1}. Total organic carbon contents were below 2%. Prior to the gas breakthrough experiments the absolute (single phase) permeability (k_abs) of the samples was determined by steady state flow tests with water or NaCl brine. The k_abs values ranged between 3 and 550 nDarcy (3 × 10^?21 and 5.5 × 10^?19 m²). The maximum effective permeability to the gas‐phase (k_eff) measured after gas breakthrough on initially water‐saturated samples extended from 0.01 nDarcy (1 × 10^?23 m²) up to 1100 nDarcy (1.1 × 10^?18 m²). The residual differential pressures after re‐imbibition of the water phase, referred to as the ‘minimum capillary displacement pressures’ (P_d), ranged from 0.06 to 6.7 MPa. During the re‐imbibition process the effective permeability to the gas phase decreases with decreasing differential pressure. The recorded permeability/pressure data were used to derive the pore size distribution (mostly between 8 and 60 nm) and the transport porosity of the conducting pore system (10^‐5–10^‐2%). Correlations could be established between (i) absolute permeability coefficients and the maximum effective permeability coefficients and (ii) effective or absolute permeability coefficients and capillary sealing efficiency. No correlation was found between the capillary displacement pressures determined from gas breakthrough experiments and those derived theoretically by mercury injection.     

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.     

Characteristics of CO2‐driven cold‐water geyser,Crystal Geyser in Utah: experimental observation and mechanism analyses

Geologic carbon capture and storage (CCS) is an option for reducing CO₂ emissions, but leakage to the surface is a risk factor. Natural CO₂ reservoirs that erupt from abandoned oil and gas holes leak to the surface as spectacular cold geysers in the Colorado Plateau, United States. A better understanding of the mechanisms of CO₂‐driven cold‐water geysers will provide valuable insight about the potential modes of leakage from engineered CCS sites. A notable example of a CO₂‐driven cold‐water geyser is Crystal Geyser in central Utah. We investigated the fluid mechanics of this regularly erupting geyser by instrumenting its conduit with sensors and measuring pressure and temperature every 20 sec over a period of 17 days. Analyses of these measurements suggest that the timescale of a single‐eruption cycle is composed of four successive eruption types with two recharge periods ranging from 30 to 40 h. Current eruption patterns exhibit a bimodal distribution, but these patterns evolved during past 80 years. The field observation suggests that the geyser's eruptions are regular and predictable and reflect pressure and temperature changes resulting from Joule–Thomson cooling and endothermic CO₂ exsolution. The eruption interval between multiple small‐scale eruptions is a direct indicator of the subsequent large‐scale eruption.     

Dry CO2–CO fluid as an important potential deep Earth solvent

Transport properties of reduced carbonic fluid have been studied experimentally at P = 2 kbar and T = 700–1000°C in internally heated pressure vessel (IHPV). Synthetic FeCO₃ and natural siderite were used to generate fluid during experiments using a platinum double‐capsule technique. A natural CaTiSiO₅ aggregate was placed into the inner capsule as an additional source of trace elements. The outer capsule was loaded with albite glass. No water was introduced to the system and oxygen fugacity was established near to graphite–oxygen (CCO) buffer due to transformation of FeCO₃ into a magnetite aggregate during decarbonation to yield CO and CO₂. The carbonates decomposed during initial heating of the experiments, causing their some constituent components to be dissolved in and transferred by the fluid to the pore space of the albite glass matrix. After temperature reached 1000°C glass, the shards annealed and then melted, as evidenced by a vesiculated glass in the quench products. Micro‐Raman investigation of the fluid in bubbles in the albite glass in experiments with decomposition of natural siderite yielded CO–CO₂ mixture where CO mole fraction was 0.15–0.16. We observe significant concentrations of Pt, Mn, P, and REE in the albite glass; in contrast, no Fe or Mg transfer was detected. LA‐ICP‐MS analysis of the albite glass product yielded the average Pt content of 2 ppm. Such high Pt signal came from Pt particles (100–500 nm in size), which were observed on the walls of the bubbles embedded in the glass. Olivines and aluminous spinel were observed in the Fe‐oxide aggregate, demonstrating transfer of SiO₂ and Al₂O₃ from the albite melt by the reduced carbonic fluid from the albite glass (large capsule). Our results demonstrate that dry CO–CO2 fluid can be important agents of dissolution and transport, especially for Pt and other metals. The data imply that metals are chiefly dissolved as carbonyl complexes.     

Permeability evolution in sorbing media: analogies between organic‐rich shale and coal

Shale gas reservoirs like coalbed methane (CBM) reservoirs are promising targets for geological sequestration of carbon dioxide (CO₂). However, the evolution of permeability in shale reservoirs on injection of CO₂ is poorly understood unlike CBM reservoirs. In this study, we report measurements of permeability evolution in shales infiltrated separately by nonsorbing (He) and sorbing (CO₂) gases under varying gas pressures and confining stresses. Experiments are completed on Pennsylvanian shales containing both natural and artificial fractures under nonpropped and propped conditions. We use the models for permeability evolution in coal (Journal of Petroleum Science and Engineering, Under Revision) to codify the permeability evolution observed in the shale samples. It is observed that for a naturally fractured shale, the He permeability increases by approximately 15% as effective stress is reduced by increasing the gas pressure from 1 MPa to 6 MPa at constant confining stress of 10 MPa. Conversely, the CO₂ permeability reduces by a factor of two under similar conditions. A second core is split with a fine saw to create a smooth artificial fracture and the permeabilities are measured for both nonpropped and propped fractures. The He permeability of a propped artificial fracture is approximately 2‐ to 3fold that of the nonpropped fracture. The He permeability increases with gas pressure under constant confining stress for both nonpropped and propped cases. However, the CO₂ permeability of the propped fracture decreases by between one‐half to one‐third as the gas pressure increases from 1 to 4 MPa at constant confining stress. Interestingly, the CO₂ permeability of nonpropped fracture increases with gas pressure at constant confining stress. The permeability evolution of nonpropped and propped artificial fractures in shale is found to be similar to those observed in coals but the extent of permeability reduction by swelling is much lower in shale due to its lower organic content. Optical profilometry is used to quantify the surface roughness. The changes in surface roughness indicate significant influence of proppant indentation on fracture surface in the shale sample. The trends of permeability evolution on injection of CO₂ in coals and shales are found analogous; therefore, the permeability evolution models previously developed for coals are adopted to explain the permeability evolution in shales.     

Geometry‐coupled reactive fluid transport at the fracture scale: application to CO2 geologic storage

Water acidification follows CO₂ injection and leads to reactive fluid transport through pores and rock fractures, with potential implications to reservoirs and wells in CO₂ geologic storage and enhanced oil recovery. Kinetic rate laws for dissolution reactions in calcite and anorthite are combined with the Navier‐Stokes law and advection–diffusion transport to perform geometry‐coupled numerical simulations in order to study the evolution of chemical reactions, species concentration, and fracture morphology. Results are summarized as a function of two dimensionless parameters: the Damköhler number Da which is the ratio between advection and reaction times, and the transverse Peclet number Pe defined as the ratio between the time for diffusion across the fracture and the time for advection along the fracture. Reactant species are readily consumed near the inlet in a carbonate reservoir when the flow velocity is low (low transverse Peclet number and Da > 10^?1). At high flow velocities, diffusion fails to homogenize the concentration field across the fracture (high transverse Peclet number Pe > 10^?1). When the reaction rate is low as in anorthite reservoirs (Da < 10^?1), reactant species are more readily transported toward the outlet. At a given Peclet number, a lower Damköhler number causes the flow channel to experience a more uniform aperture enlargement along the length of the fracture. When the length‐to‐aperture ratio is sufficiently large, say l/d > 30, the system response resembles the solution for 1D reactive fluid transport. A decreased length‐to‐aperture ratio slows the diffusive transport of reactant species to the mineral fracture surface, and analyses of fracture networks must take into consideration both the length and slenderness of individual fractures in addition to Pe and Da numbers.     

Basin‐scale fluid movement patterns revealed by veins: Wessex Basin,UK

Calcite veins at outcrop in the Mesozoic, oil‐bearing Wessex Basin, UK, have been studied using field characterization, petrography, fluid inclusions and stable isotopes to help address the extent, timing and spatial and stratigraphic variability of basin‐scale fluid flow. The absence of quartz shows that veins formed at low temperature without an influence of hydrothermal fluids. Carbon isotopes suggest that the majority of vein calcite was derived locally from the host rock but up to one quarter of the carbon in the vein calcite came from CO₂ from petroleum source rocks. Veins become progressively enriched in source‐rock‐derived CO₂ from the outer margin towards the middle, indicating a growing influence of external CO₂. The carbon isotope data suggest large‐scale migration of substantial amounts of CO₂ around the whole basin. Fluid inclusion salinity data and interpreted water‐δ¹⁸O data show that meteoric water penetrated deep into the western part of the basin after interacting with halite‐rich evaporites in the Triassic section before entering fractured Lower and Middle Jurassic rocks. This large‐scale meteoric invasion of the basin probably happened during early Cenozoic uplift. A similar approach was used to reveal that, in the eastern part of the basin close to the area that underwent most uplift, uppermost Jurassic and Cretaceous rocks underwent vein formation in the presence of marine connate water suggesting a closed system. Stratigraphically underlying Upper Jurassic mudstone and Lower Cretaceous sandstone, in the most uplifted part of the basin, contain veins that resulted from intermediate behaviour with input from saline meteoric water and marine connate waters. Thus, while source‐rock‐derived CO₂ seems to have permeated the entire section, water movement has been more restricted. Oil‐filled inclusions in vein calcite have been found within dominant E‐W trending normal faults, suggesting that these may have facilitated oil migration.