Exploring the potential linkages between oil‐field brine reinjection,crystalline basement permeability,and triggered seismicity for the Dagger Draw Oil field,southeastern New Mexico,USA, using hydrologic modeling

We used hydrologic models to explore the potential linkages between oil‐field brine reinjection and increases in earthquake frequency (up to M_d 3.26) in southeastern New Mexico and to assess different injection management scenarios aimed at reducing the risk of triggered seismicity. Our analysis focuses on saline water reinjection into the basal Ellenburger Group beneath the Dagger Draw Oil field, Permian Basin. Increased seismic frequency (>M_d 2) began in 2001, 5 years after peak injection, at an average depth of 11 km within the basement 15 km to the west of the reinjection wells. We considered several scenarios including assigning an effective or bulk permeability value to the crystalline basement, including a conductive fault zone surrounded by tighter crystalline basement rocks, and allowing permeability to decay with depth. We initially adopted a 7 m (0.07 MPa) head increase as the threshold for triggered seismicity. Only two scenarios produced excess heads of 7m five years after peak injection. In the first, a hydraulic diffusivity of 0.1 m² s^?1 was assigned to the crystalline basement. In the second, a hydraulic diffusivity of 0.3 m² s^?1 was assigned to a conductive fault zone. If we had considered a wider range of threshold excess heads to be between 1 and 60 m, then the range of acceptable hydraulic diffusivities would have increased (between 0.1–0.01 m² s^?1 and 1–0.1 m² s^?1 for the bulk and fault zone scenarios, respectively). A permeability–depth decay model would have also satisfied the 5‐year time lag criterion. We also tested several injection management scenarios including redistributing injection volumes between various wells and lowering the total volume of injected fluids. Scenarios that reduced computed excess heads by over 50% within the crystalline basement resulted from reducing the total volume of reinjected fluids by a factor of 2 or more.     

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.     

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.     

Interpretation of observed fluid potential patterns in a deep sedimentary basin under tectonic compression: Hungarian Great Plain, Pannonian Basin

The ≈ 40 000 km² Hungarian Great Plain portion of the Pannonian Basin consists of a basin fill of 100 m to more than 7000 m thick semi‐ to unconsolidated marine, deltaic, lacustrine and fluviatile clastic sediments of Neogene age, resting on a strongly tectonized Pre‐Neogene basement of horst‐and‐graben topography of a relief in excess of 5000 m. The basement is built of a great variety of brittle rocks, including flysch, carbonates and metamorphics. The relatively continuous Endr?d Aquitard, with a permeability of less than 1 md (10^?15 m²) and a depth varying between 500 and 5000 m, divides the basin's rock framework into upper and lower sequences of highly permeable rock units, whose permeabilities range from a few tens to several thousands of millidarcy. Subsurface fluid potential and flow fields were inferred from 16 192 water level and pore pressure measurements using three methods of representation: pressure–elevation profiles; hydraulic head maps; and hydraulic cross‐sections. Pressure–elevation profiles were constructed for eight areas. Typically, they start from the surface with a straight‐line segment of a hydrostatic gradient (γ_st = 9.8067 MPa km^?1) and extend to depths of 1400–2500 m. At high surface elevations, the gradient is slightly smaller than hydrostatic, while at low elevations it is slightly greater. At greater depths, both the pressures and their vertical gradients are uniformly superhydrostatic. The transition to the overpressured depths may be gradual, with a gradient of γ_dyn = 10–15 MPa km^?1 over a vertical distance of 400–1000 m, or abrupt, with a pressure jump of up to 10 MPa km^?1 over less than 100 m and a gradient of γ_dyn > 20 MPa km^?1. According to the hydraulic head maps for 13 100–500 m thick horizontal slices of the rock framework, the fluid potential in the near‐surface domains declines with depth beneath positive topographic features, but it increases beneath depressions. The approximate boundary between these hydraulically contrasting regions is the 100 m elevation contour line in the Duna–Tisza interfluve, and the 100–110 m contours in the Nyírség uplands. Below depths of ≈ 600 m, islets of superhydrostatic heads develop which grow in number, areal extent and height as the depth increases; hydraulic heads may exceed 3000 m locally. A hydraulic head ‘escarpment’ appears gradually in the elevation range of ? 1000 to ? 2800 m along an arcuate line which tracks a major regional fault zone striking NE–SW: heads drop stepwise by several hundred metres, at places 2000 m, from its north and west sides to the south and east. The escarpment forms a ‘fluid potential bank’ between a ‘fluid potential highland’ (500–2500 m) to the north and west, and a ‘fluid potential basin’ (100–500 m) to the south and east. A ‘potential island’ rises 1000 m high above this basin further south. According to four vertical hydraulic sections, groundwater flow is controlled by the topography in the upper 200–1700 m of the basin; the driving force is orientated downwards beneath the highlands and upwards beneath the lowlands. However, it is directed uniformly upwards at greater depths. The transition between the two regimes may be gradual or abrupt, as indicated by wide or dense spacing of the hydraulic head contours, respectively. Pressure ‘plumes’ or ‘ridges’ may protrude to shallow depths along faults originating in the basement. The basement horsts appear to be overpressured relative to the intervening grabens. The principal thesis of this paper is that the two main driving forces of fluid flow in the basin are gravitation, due to elevation differences of the topographic relief, and tectonic compression. The flow field is unconfined in the gravitational regime, whereas it is confined in the compressional regime. The nature and geometry of the fluid potential field between the two regimes are controlled by the sedimentary and structural features of the rock units in that domain, characterized by highly permeable and localized sedimentary windows, conductive faults and fracture zones. The transition between the two potential fields can be gradual or abrupt in the vertical, and island‐like or ridge‐like in plan view. The depth of the boundary zone can vary between 400 and 2000 m. Recharge to the gravitational regime is inferred to occur from infiltrating precipitation water, whereas that to the confined regime is from pore volume reduction due to the basement's tectonic compression.     

Injection‐induced seismicity in Carbon and Emery Counties,central Utah

Utah is one of the top producers of oil and natural gas in the United States. Over the past 18 years, more than 4.2 billion gallons of wastewater from the petroleum industry has been injected into the Navajo Sandstone, Kayenta Formation, and Wingate Sandstone in Carbon and Emery Counties, central Utah, where seismicity has increased during the same period. Previous studies have attributed this seismicity to coal mining. Here, we present evidence for wastewater injection being a major cause of the increased seismicity. We show that, in the coal mining area, seismicity rate increased significantly 1–5 years following the wastewater injection, and the earthquakes, mostly with magnitudes <3.0, are concentrated in areas seismically active prior to the injection. Using simple analytical and numerical models, we show that the injection in central Utah can sufficiently raise pore pressure to trigger seismicity within 10–20 km of the injection wells, and the time needed for the diffusion of pore pressure may explain the observed lag of seismicity increase behind the commencement of injection. The b‐value of these earthquakes increased following the wastewater injection, which is consistent with these events being injection‐induced. We conclude that the marked increase in seismicity rate in central Utah is induced by both mining activity and wastewater injection, which raised pore pressure along preexisting faults.     

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.     

Modeling of hydrogen production by serpentinization in ultramafic‐hosted hydrothermal systems: application to the Rainbow field

The production of hydrogen by serpentinization in ultramafic‐hosted hydrothermal systems is simulated by coupling thermodynamic and dynamic modeling in the framework of a thermo‐hydraulic single‐pass model where a high‐temperature hydrothermal fluid moves preferentially through a main canal of high permeability. The alteration of ultramafic rocks is modeled with a first‐order kinetic formulation, wherein the serpentinization rate coefficient, K_r, takes the form: K_r = A exp(?α(T ? T₀)²). In this formulation, α determines the temperature range of the reaction and T₀ is the temperature at which the serpentinization rate reaches its maximum. This model is applied to the Rainbow hydrothermal system, which is situated on the Mid‐Atlantic Ridge, and characterized by a high temperature, a high mass flux, and a very high hydrogen concentration. The results show that a first‐order kinetic law gives a useful representation of the kinetics of serpentinization. The estimated value for the parameter A in the temperature‐dependent formulation of the serpentinization rate coefficient lies in the range (1–5) × 10^?11 s^?1. This effective parameter is several orders of magnitude lower than the values obtained from small grain‐size experiments, but in agreement with other published modeling studies of natural systems. Numerical simulations show that the venting site is able to produce the observed high concentration of hydrogen during the whole continuous lifetime of the Rainbow site.     

Large‐scale chemical stratification of fluids in the crust: hydraulic and chemical data from the geothermal research site Urach,Germany

The Urach 3 research borehole in SW Germany has been drilled through a sedimentary cover sequence and reached gneisses of the Variscan crystalline basement at 1604 m below surface. An additional 2840 m has been drilled through fractured basement rocks. The borehole has been used for hydraulic tests in the context of a ‘hot dry rock’ (HDR) project. The sedimentary cover ranges from the Carboniferous to the Middle Jurassic (Dogger) in age and comprises mostly clastic sediments in the Paleozoic and limestone and shale in the Mesozoic. Water composition data from 10 different depths include samples from all major lithological units. The total dissolved solids (TDS) increases from the surface to about 650 m where it reaches 4.1 g l^?1 in Triassic limestone. In lower Triassic sandstones, TDS increases very sharply to 28.5 g l^?1 and the water is saturated with pure CO₂ gas. With increasing depth, TDS does not change much in the clastic sediments of the Permian and Carboniferous. The crystalline basement is marked by a very sharp increase in TDS to 55.5 g l^?1 at about 1770 m depth. TDS increases within the basement to more than 78.5 g l^?1 at about 3500 m depth. The data suggest that there is limited vertical chemical communication over long periods of time. The CO₂ gas cap in the lower Triassic sandstones requires a gastight cover. The chemical stratification of the fluids relates to the permeability structure of the crust at the Urach site and fits well with hydraulic and thermal data from the site.     

Effects of episodic fluid flow on hydrocarbon migration in the Newport‐Inglewood Fault Zone,Southern California

Fault permeability may vary through time due to tectonic deformations, transients in pore pressure and effective stress, and mineralization associated with water‐rock reactions. Time‐varying permeability will affect subsurface fluid migration rates and patterns of petroleum accumulation in densely faulted sedimentary basins such as those associated with the borderland basins of Southern California. This study explores the petroleum fluid dynamics of this migration. As a multiphase flow and petroleum migration case study on the role of faults, computational models for both episodic and continuous hydrocarbon migration are constructed to investigate large‐scale fluid flow and petroleum accumulation along a northern section of the Newport‐Inglewood fault zone in the Los Angeles basin, Southern California. The numerical code solves the governing equations for oil, water, and heat transport in heterogeneous and anisotropic geologic cross sections but neglects flow in the third dimension for practical applications. Our numerical results suggest that fault permeability and fluid pressure fluctuations are crucial factors for distributing hydrocarbon accumulations associated with fault zones, and they also play important roles in controlling the geologic timing for reservoir filling. Episodic flow appears to enhance hydrocarbon accumulation more strongly by enabling stepwise build‐up in oil saturation in adjacent sedimentary formations due to temporally high pore pressure and high permeability caused by periodic fault rupture. Under assumptions that fault permeability fluctuate within the range of 1–1000 millidarcys (10^?15–10^?12 m²) and fault pressures fluctuate within 10–80% of overpressure ratio, the estimated oil volume in the Inglewood oil field (approximately 450 million barrels oil equivalent) can be accumulated in about 24 000 years, assuming a seismically induced fluid flow event occurs every 2000 years. This episodic petroleum migration model could be more geologically important than a continuous‐flow model, when considering the observed patterns of hydrocarbons and seismically active tectonic setting of the Los Angeles basin.     

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

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

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.     

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.     

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.     

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.     

Pore pressure evolution and fluid flow during visco‐elastic single‐layer buckle folding

Pore pressure and fluid flow during the deformational history of geologic structures are directly influenced by tectonic deformation events. In this contribution, 2D plane strain finite element analysis is used to study the influence of different permeability distributions on the pore pressure field and associated flow regimes during the evolution of visco‐elastic single‐layer buckle folds. The buckling‐induced fluid flow regimes indicate that flow directions and, to a lesser degree, their magnitudes vary significantly throughout the deformation and as a function of the stratigraphic permeability distribution. The modelling results suggest that the volumetric strain and the permeability distribution significantly affect the resulting flow regime at different stages of fold development. For homogeneous permeability models (k > 10^?21 m²), low strain results in a mostly pervasive fluid flow regime and is in agreement with previous studies. For larger strain conditions, fluid focusing occurs in the buckling layer towards the top of the fold hinge. For low permeabilities (<10^?21 m²), local focused flow regimes inside the buckling layer emerge throughout the deformation history. For models featuring a low‐permeability layer embedded in a high‐permeability matrix or sandwiched between high‐permeability layers, focused flow regimes inside the folded layer result throughout the deformation history, but with significant differences in the flow vectors of the surrounding layers. Fluid flow vectors induced by the fold can result in different, even reversed, directions depending on the amount of strain. In summary, fluid flow regimes during single‐layer buckling can change from pervasive to focused and fluid flow vectors can be opposite at different strain levels, that is the flow vectors change significantly through time. Thus, a complete understanding of fluid flow regimes associated with single‐layer buckle folds requires consideration of the complete deformation history of the fold.     

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.     

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.     

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.     

Fluid inclusion and stable isotope constraints on the genesis of the Jian copper deposit,Sanandaj–Sirjan metamorphic zone,Iran

The Jian copper deposit, located on the eastern edge of the Sanandaj–Sirjan metamorphic zone, southwest of Iran, is contained within the Surian Permo‐Triassic volcano‐sedimentary complex. Retrograde metamorphism resulted in three stages of mineralization (quartz ± sulfide veins) during exhumation of the Surian metamorphic complex (Middle Jurassic time; 159–167 Ma), and after the peak of the metamorphism (Middle to Late Triassic time; approximately 187 Ma). The early stage of mineralization (stage 1) is related to a homogeneous H₂O–CO₂ (XCO₂ > 0.1) fluid characterized by moderate salinity (<10 wt.% NaCl equivalent) at high temperature and pressure (>370°C, >3 kbar). Early quartz was followed by small amounts of disseminated fine‐grained pyrite and chalcopyrite. Most of the main‐ore‐stage (stage 2) minerals, including chalcopyrite, pyrite and minor sphalerite, pyrrhotite, and galena, precipitated from an aqueous‐carbonic fluid (8–18 wt.% NaCl equivalent) at temperatures ranging between 241 and 388°C during fluid unmixing process (CO₂ effervescence). Fluid unmixing in the primary carbonaceous fluid at pressures of 1.5–3 kbar produced a high XCO₂ (>0.05) and a low XCO₂ (<0.01) aqueous fluid in ore‐bearing quartz veins. Oxygen and hydrogen isotope compositions suggest mineralization by fluids derived from metamorphic dehydration (δ¹⁸O_fluid = +7.6 to +10.7‰ and δD = ?33.1 to ?38.5‰) during stage 2. The late stage (stage 3) is related to a distinct low salinity (1.5–8 wt.% NaCl equivalent) and temperatures of (120–230°C) aqueous fluid at pressures below 1.5 kbar and the deposition of post‐ore barren quartz veins. These fluids probably derived from meteoric waters, which circulated through the metamorphic pile at sufficiently high temperatures and acquire the characteristics of metamorphic fluids (δ¹⁸O_fluid = +4.7 to +5.1‰ and δD = ?52.3 to ?53.9‰) during waning stages of the postearly Cimmerian orogeny in Surian complex. The sulfide‐bearing quartz veins are interpreted as a small‐scale example of redistribution of mineral deposits by metamorphic fluids. This study suggests that mineralization at the Jian deposit is metamorphogenic in style, probably related to a deep‐seated mesothermal system.     

Experimental determination of CePO4 and YPO4 solubilities in H2O–NaF at 800°C and 1 GPa: implications for rare earth element transport in high‐grade metamorphic fluids

Monazite (CePO₄) and xenotime (YPO₄) are important accessory minerals in metasediments. They host significant rare earth elements (REE) and are useful for geochronology and geothermometry, so it is essential to understand their behavior during the metasomatic processes that attend high‐grade metamorphism. It has been proposed that F‐bearing fluids enhance solubility and mobility of REE and Y during high‐grade metamorphism. We assessed this possibility by determining the solubility of synthetic CePO₄ and YPO₄ crystals in H₂O–NaF fluids at 800°C and 1 GPa. Experiments used hydrothermal piston‐cylinder and weight‐loss methods. Compared to the low solubilities of CePO₄ and YPO₄ in pure H₂O (0.04 ± 0.04 and 0.25 ± 0.04 millimolal, respectively), our results indicate an enormous increase in the solubility of both phosphates with increasing NaF concentration in H₂O: CePO₄ solubility reaches 0.97 molal in 20 mol.% NaF, and YPO₄ shows an even stronger solubility enhancement to 0.45 molal in only 10 mol.% NaF. The greatest relative solubility increases occur at the lowest NaF concentration. The solubilities of CePO₄ and YPO₄ show similar quadratic dependence on NaF, consistent with possible dissolution reactions of: CePO₄ + 2NaF =  CeF₂⁺ + Na₂PO₄^? and YPO₄ + 2NaF = YF₂⁺ + Na₂PO₄^?. Solubilities of both REE phosphates are significantly greater in NaF than in NaCl at equivalent salt concentration. A fluid with 10 mol.% NaCl and multiply saturated with fluorite, CePO₄, and YPO₄ would contain 1.7 millimolal Ce and 3.3 millimolal Y, values that are respectively 2.1–2.4 times greater than in NaCl‐H₂O alone. The results indicate that Y, and by extension heavy rare earth elements (HREE), can be fractionated from LREE in fluorine‐bearing saline brines which may accompany granulite‐facies metamorphism. The new data support previous indications that REE/Y mobility at these conditions is enhanced by complexing with F in the aqueous phase.