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.     

Salt tectonics and shallow subseafloor fluid convection: models of coupled fluid‐heat‐salt transport

Thermohaline convection associated with salt domes has the potential to drive significant fluid flow and mass and heat transport in continental margins, but previous studies of fluid flow associated with salt structures have focused on continental settings or deep flow systems of importance to petroleum exploration. Motivated by recent geophysical and geochemical observations that suggest a convective pattern to near‐seafloor pore fluid flow in the northern Gulf of Mexico (GoMex), we devise numerical models that fully couple thermal and chemical processes to quantify the effects of salt geometry and seafloor relief on fluid flow beneath the seafloor. Steady‐state models that ignore halite dissolution demonstrate that seafloor relief plays an important role in the evolution of shallow geothermal convection cells and that salt at depth can contribute a thermal component to this convection. The inclusion of faults causes significant, but highly localized, increases in flow rates at seafloor discharge zones. Transient models that include halite dissolution show the evolution of flow during brine formation from early salt‐driven convection to later geothermal convection, characteristics of which are controlled by the interplay of seafloor relief and salt geometry. Predicted flow rates are on the order of a few millimeters per year or less for homogeneous sediments with a permeability of 10^?15 m², comparable to compaction‐driven flow rates. Sediment permeabilities likely fall below 10^?15 m² at depth in the GoMex basin, but such thermohaline convection can drive pervasive mass transport across the seafloor, affecting sediment diagenesis in shallow sediments. In more permeable settings, such flow could affect methane hydrate stability, seafloor chemosynthetic communities, and the longevity of fluid seeps.     

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.     

Analysis of strain‐induced ground‐water fluctuations at Devils Hole,Nevada

The pool in Devils Hole is a sensitive indicator of crustal strain and fluctuates in response to changes in atmospheric pressure, earth tides, earthquakes, large‐scale tectonic activity and ground‐water development. Short‐term and cyclic water‐level fluctuations caused by atmospheric pressure and earth tides were found to be on the order of millimeters to centimeters. The 1992 Landers/Little Skull Mountain earthquake sequence and the 1999 Hector Mine earthquake induced water‐level offsets of greater than ?12 and ?3.6 cm, respectively. The results of a dislocation model used to compute volumetric strain for each earthquake indicates that the coseismic water‐level offsets are consistent in magnitude and sense with poroelastic responses to earthquake‐induced strain. Theoretical postseismic fluid‐flow modeling indicates that the diffusivity of the system is on the order of 0.03 m² sec^?1, and identified areas of anomalous water‐level fluctuations. Interpretation of model results suggests that while the persistent post‐Landers rise in water‐level can be attributed to deformation‐induced channeling of fluid to the Devils Hole fault zone, the cause of the pre‐Hector Mine water‐level rise may be related to postseismic excess fluid pressures or preseismic strain accumulation.     

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.     

Modelling the Lost City hydrothermal field: influence of topography and permeability structure

The Lost City hydrothermal field (LCHF) is hosted in serpentinite at the crest of the Atlantis Massif, an oceanic core complex close to the mid‐Atlantic Ridge. It is remarkable for its longevity and for venting low‐temperature (40–91°C) alkaline fluids rich in hydrogen and methane. IODP Hole U1309D, 5 km north of the LCHF, penetrated 1415 m of gabbroic rocks and contains a near‐conductive thermal gradient close to 100°C km^?1. This is remarkable so close to an active hydrothermal field. We present hydrothermal modelling using a topographic profile through the vent field and IODP site U1309. Long‐lived circulation with vent temperatures similar to the LCHF can be sustained at moderate permeabilities of 10^?14 to 10^?15 m² with a basal heatflow of 0.22 W m^?2. Seafloor topography is an important control, with vents tending to form and remain in higher topography. Models with a uniform permeability throughout the Massif cannot simultaneously maintain circulation at the LCHF and the near‐conductive gradient in the borehole, where permeabilities <10^?16 m² are required. A steeply dipping permeability discontinuity between the LCHF and the drill hole is required to stabilize venting at the summit of the massif by creating a lateral conductive boundary layer. The discontinuity needs to be close to the vent site, supporting previous inferences that high permeability is most likely produced by faulting related to the transform fault. Rapid increases in modelled fluid temperatures with depth beneath the vent agree with previous estimates of reaction temperature based on geochemical modelling.     

Evolution of deformation and fault‐related fluid flow within an ancient methane seep system (Eocene,Varna, Bulgaria)

Faults are often important in fuelling methane seep systems; however, little is known on how different components in fault zones control subsurface fluid circulation paths and how they evolve through time. This study provides insight into fault‐related fluid flow systems that operated in the shallow subsurface of an ancient methane seep system. The Pobiti Kamani area (NE Bulgaria) encloses a well‐exposed, fault‐related seep system in unconsolidated Lower Eocene sandy deposits of the Dikilitash Formation. The Beloslav quarry and Beloslav N faults displace the Dikilitash Formation and are typified by broad, up to 80 m wide, preferentially lithified hanging wall damage zones, crosscut by deformation bands and deformation band zones, smaller slip planes and fault‐related joints. The formation of a shallow plumbing system and chimney‐like concretions in the Dikilitash Formation was followed by at least two phases of fault‐related methane fluid migration. Widespread fluid circulation through the Dikilitash sands caused massive cementation of the entire damage zones in the fault hanging walls. During this phase, paths of ascending methane fluids were locally obstructed by decimetre‐thick, continuous deformation band zones that developed in the partly lithified sands upon the onset of deformation. Once the entire damage zone was pervasively cemented, deformation proceeded through the formation of slip planes and joints. This created a new network of more localized conduits in close vicinity to the main fault plane and around through‐going slip planes. ¹³C‐depleted crustiform calcite cements in several joints record the last phase of focused methane fluid ascent. Their formation predated Neogene uplift and later meteoric water infiltration along the joint network. This illustrates how fault‐related fluid pathways evolved, over time, from ‘plumes’ in unconsolidated sediments above damage zones, leading to chimney fields, over widespread fluid paths, deflected by early deformation structures, to localized paths along fracture networks near the main fault.     

Stress‐dependent transport properties of fractured arkosic sandstone

We measure the fluid transport properties of microfractures and macrofractures in low‐porosity polyphase sandstone and investigate the controls of in situ stress state on fluid flow conduits in fractured rock. For this study, the permeability and porosity of the Punchbowl Formation sandstone, a hydrothermally altered arkosic sandstone, were measured and mapped in stress space under intact, microfractured, and macrofractured deformation states. In contrast to crystalline and other sedimentary rocks, the distributed intragranular and grain‐boundary microfracturing that precedes macroscopic fracture formation has little effect on the fluid transport properties. The permeability and porosity of microfractured and intact sandstone depend strongly on mean stress and are relatively insensitive to differential stress and proximity to the frictional sliding envelope. Porosity variations occur by elastic pore closure with intergranular sliding and pore collapse caused by microfracturing along weakly cemented grain contacts. The macroscopic fractured samples are best described as a two‐component system consisting (i) a tabular fracture with a 0.5‐mm‐thick gouge zone bounded by 1 mm thick zones of concentrated transgranular and intragranular microfractures and (ii) damaged sandstone. Using bulk porosity and permeability measurements and finite element methods models, we show that the tabular fracture is at least two orders of magnitude more permeable than the host rock at mean stresses up to 90 MPa. Further, we show that the tabular fracture zone dilates as the stress state approaches the friction envelope resulting in up to a three order of magnitude increase in fracture permeability. These results indicate that the enhanced and stress‐sensitive permeability in fault damage zones and sedimentary basins composed of arkosic sandstones will be controlled by the distribution of macroscopic fractures rather than microfractures.     

Inferring intermediate‐scale fluid flow in a heterogeneous metasedimentary multilayer sequence during progressive deformation: evidence from the Monts d’Arrée slate belt (Brittany,France)

Quartz veins in the early Variscan Monts d’Arrée slate belt (Central Armorican Terrane, Western France), have been used to determine fluid‐flow characteristics. A combination of a detailed structural analysis, fluid inclusion microthermometry and stable isotope analyses provides insights in the scale of fluid flow and the water–rock interactions. This research suggests that fluids were expelled during progressive deformation and underwent an evolution in fluid chemistry because of changing redox conditions. Seven quartz‐vein generations were identified in the metasedimentary multilayer sequence of the Upper Silurian to Lower Devonian Plougastel Formation, and placed within the time frame of the deformation history. Fluid inclusion data of primary inclusions in syn‐ to post‐tectonic vein generations indicate a gradual increase in methane content of the aqueous–gaseous H₂O–CO₂–NaCl–CH₄–N₂ fluid during similar P–T conditions (350–400°C and 2–3.5 kbar). The heterogeneous centimetre‐ to metre‐scale multilayer sequence of quartzites and phyllites has a range of oxygen‐isotope values (8.0–14.1‰ Vienna Standard Mean Ocean Water), which is comparable with the range in the crosscutting quartz veins (10.5–14.7‰ V‐SMOW). Significant differences between oxygen‐isotope values of veins and adjacent host rock (Δ = ?2.8‰ to +4.9‰ V‐SMOW) suggest an absence of host‐rock buffering on a centimetre scale, but based on the similar range of isotope values in the Plougastel Formation, an intraformational buffering and an intermediate‐scale fluid‐flow system could be inferred. The abundance of veins, their well‐distributed and isolated occurrence, and their direct relationship with the progressive deformation suggests that the intermediate‐scale fluid‐flow system primarily occurred in a dynamically generated network of temporarily open fractures.     

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.     

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.     

The relevance of vein texture in understanding the past hydraulic behaviour of a crystalline rock mass: reconstruction of the palaeohydrology of the Mecsekalja Zone,south Hungary

This study reconstructs the palaeohydrogeologic evolution of the shallow‐to‐moderate Mesozoic subsidence history for the Mecsekalja Zone (MZ), a narrow metamorphic belt in the eastern Mecsek Mountains, Hungary. Brittle deformation of the MZ produced a vein system with a cement history consisting of five sequential carbonate generations and one quartz phase. Vein textures suggest different fluid‐flow mechanisms for the parent fluids of subsequent cement generations. Combined microthermometric and stable‐isotope measurements permit reconstruction of the character of subsequent fluid generations with different flow types, as defined by vein textures, yielding new information regarding the hydraulic behaviour of a metamorphic crystalline complex. Textural observations and geochemical data suggest that fracture‐controlled flow pathways and externally derived fluids were typical of some flow events, while percolation through the rock matrix and the relationship to the Cretaceous volcanism and dyke emplacement were typical of others. The difference in the mode of calcite deposition from pervasive fluids (i.e. pervasive carbonatisation along grain boundaries versus deposition in antitaxial veins) between two calcite generations related to the volcanism inspired a stress‐dependent model of antitaxial vein growth. Textural and isotope variations in a vein generation produced by the same parent fluid indicate rock‐dependent hydraulic behaviour for different rock types, distinct action of the contemporaneous fracture systems and different extents of fluid–rock interaction. Cathodoluminescence microscopy and fluid‐inclusion microthermometry shed light on the possible role of hydraulic fracturing in the formation of massive calcite. The time of formation was estimated from the isotope composition of the oldest calcite generation and its presumptive relationship with the sedimentary sequences to the north, whereas microthermometry permitted conciliation of the reconstructed flow sequence with the Mesozoic subsidence history of the Mórágy Block (including the MZ).     

Topography‐driven flow versus buoyancy‐driven flow in the U.S. midcontinent: implications for the residence time of brines

Topography‐driven flow is normally considered to be the dominant groundwater flow system in uplifted sedimentary basins. In the U.S. midcontinent region east of the Rocky Mountains, the presence of brines derived from dissolution of halite suggests that significant topography‐driven flushing has occurred to remove older brines that presumably formed concurrently with Permian evaporites in the basin. However, the presence of evaporites and brines in the modern basin suggests that buoyancy‐driven flow could limit topography‐driven flushing significantly. Here we used numerical models of variable‐density fluid flow, halite dissolution, solute transport, and heat transport to quantify flow patterns and brine migration. Results indicate the coexistence of large‐scale topography‐ and buoyancy‐driven flow. Buoyancy‐driven flow and low permeability evaporites act to isolate brines, and the residence time of the brines was found to be quite long, at least 50 Myr. The modern distribution of salinity appears to reflect near‐steady‐state conditions. Results suggest that flushing of original evaporatively‐concentrated brines occurred tens of millions of years ago, possibly concurrent with maximum uplift ca. 60 Ma. Simulations also suggest that buoyancy‐driven convection could drive chemical exchange with crystalline basement rocks, which could supply significant Ca²⁺, Sr²⁺, and metals to brines.     

Numerical model of pore‐pressure diffusion associated with the initiation of the 2010–2011 Guy–Greenbrier,Arkansas earthquakes

We model pore‐pressure diffusion caused by pressurized waste‐fluid injection at two nearby wells and then compare the buildup of pressure with the observed initiation and migration of earthquakes during the early part of the 2010–2011 Guy–Greenbrier earthquake swarm. Pore‐pressure diffusion is calculated using MODFLOW 2005 that allows the actual injection histories (volume/day) at the two wells to diffuse through a fractured and faulted 3D aquifer system representing the eastern Arkoma basin. The aquifer system is calibrated using the observed water‐level recovery following well shut‐in at three wells. We estimate that the hydraulic conductivities of the Boone Formation and Arbuckle Group are 2.2 × 10^?2 and 2.03 × 10^?3 m day^?1, respectively, with a hydraulic conductivity of 1.92 × 10^?2 m day^?1 in the Hunton Group when considering 1.72 × 10^?3 m day^?1 in the Chattanooga Shale. Based on the simulated pressure field, injection near the relatively conductive Enders and Guy–Greenbrier faults (that hydraulically connect the Arbuckle Group with the underlying basement) permits pressure diffusion into the crystalline basement, but the effective radius of influence is limited in depth by the vertical anisotropy of the hydraulic diffusivity. Comparing spatial/temporal changes in the simulated pore‐pressure field to the observed seismicity suggests that minimum pore‐pressure changes of approximately 0.009 and 0.035 MPa are sufficient to initiate seismic activity within the basement and sedimentary sections of the Guy–Greenbrier fault, respectively. Further, the migration of a second front of seismicity appears to follow the approximately 0.012 MPa and 0.055 MPa pore‐pressure fronts within the basement and sedimentary sections, respectively.     

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.     

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.     

Fracture spacing during hydro‐fracturing of cap‐rocks

Layered low permeability rock units, like shales, represent seals or ‘cap‐rocks’ in a variety of geological settings. A continuous increase in the fluid pressure gradients across a virtually impermeable rock layer will ultimately lead to hydro‐fracturing. Depending on the boundary conditions, such fracturing may lead to the formation of a set of sub‐parallel cracks oriented more or less perpendicular to the cap‐rock layer. In this article, we propose a new numerical model that describes interactions between multiple cross‐cutting fractures in an elastic low permeability rock layer. The width of each fracture and the spacing between them are modeled as a force balance between the fluid pressure and the elastic forces in the cap‐rock and between each fracture. The model indicates that the system of fractures evolves toward a spatially periodic steady‐state distribution with a fixed fracture spacing and aperture. The results are similar for incompressible and compressible fluids. The steady‐state conditions depend on only two dimensionless parameters, and the fracture spacing is only weakly dependent on the cap‐rock thickness. This is in contrast to fracturing produced by simple extension of an elastic rock layer beyond the fracture strength, in which case fracture spacing is proportional to layer thickness.     

Impacts of Pleistocene glaciation on large‐scale groundwater flow and salinity in the Michigan Basin

Pleistocene melting of kilometer‐thick continental ice sheets significantly impacted regional‐scale groundwater flow in the low‐lying stable interiors of the North American and Eurasian cratons. Glacial meltwaters penetrated hundreds of meters into the underlying sedimentary basins and fractured crystalline bedrock, disrupting relatively stagnant saline fluids and creating a strong disequilibrium pattern in fluid salinity. To constrain the impact of continental glaciation on variable density fluid flow, heat and solute transport in the Michigan Basin, we constructed a transient two‐dimensional finite‐element model of the northern half of the basin and imposed modern versus Pleistocene recharge conditions. The sag‐type basin contains up to approximately 5 km of Paleozoic strata (carbonates, siliciclastics, and bedded evaporites) overlain by a thick veneer (up to 300 m) of glacial deposits. Formation water salinity increases exponentially from <0.5 g l^?1 total dissolved solids (TDS) near the surface to >350 g l^?1 TDS at over 800 m depth. Model simulations show that modern groundwater flow is primarily restricted to shallow glacial drift aquifers with discharge to the Great Lakes. During the Pleistocene, however, high hydraulic heads from melting of the Laurentide Ice Sheet reversed regional flow patterns and focused recharge into Paleozoic carbonate and siliciclastic aquifers. Dilute waters (<20 g l^?1 TDS) migrated approximately 110 km laterally into the Devonian carbonate aquifers, significantly depressing the freshwater‐saline water mixing zones. These results are consistent with ¹⁴C ages and oxygen isotope values of confined groundwaters in Devonian carbonates along the basin margin, which reflect past recharge beneath the Laurentide Ice Sheet (14–50 ka). Constraining the paleohydrology of glaciated sedimentary basins, such as the Michigan Basin, is important for determining the source and residence times of groundwater resources, in addition to resolving geologic forces responsible for basinal‐scale fluid and solute migration.     

Effects of pore fluid flow and chemistry on compaction creep of calcite by pressure solution at 150°C

Uni‐axial compaction creep experiments were performed on crushed limestone and analytical grade calcite powders at 150°C, a pore fluid pressure of 20 MPa, and effective axial stresses of 30 and 40 MPa. Previous experiments have shown that compaction under these conditions is dominated by intergranular pressure solution (IPS). The aim of the present tests was to determine the inter‐relationship between pore fluid chemistry, compaction rate and the rate‐controlling process of IPS. Intermittent flow‐through runs conducted using CaCO₃ solution showed no effect on creep rate at low strains (<4–5%) but a major acceleration at high strains (5–10%). Measurements of the Ca concentration present in fluid samples revealed the build‐up of a high super‐saturation of CaCO₃ during compaction under zero flow conditions, especially at high strains. Active flow‐through led to a drop in Ca concentration, which corresponded with creep acceleration. Addition of NaCl to the pore fluid, at a concentration of 0.5 m , increased the creep rate of the analytical grade calcite samples roughly in proportion to the enhancement of calcite solubility. Addition of Mg²⁺ and to the pore fluid, in concentrations of 0.05 and 0.001 m, respectively, caused major retardation of compaction creep. Integrating our mechanical, flow‐through and chemical data points strongly to diffusion‐controlled IPS being the dominant deformation mechanism in the calcite‐water system under closed‐system (zero flow) conditions at low strains (<4–5%), giving way to precipitation control at higher strains. Our fluid composition data suggest that this transition is because of accumulation of impurities in the pore fluid. As Mg²⁺ and phosphate ions are common in natural pore fluids, it is likely that retarded precipitation will be the rate‐limiting step of IPS in carbonates in nature. To quantify diagenetic compaction and porosity‐permeability reduction rates by IPS in carbonates needs to account for this.     

Fracture‐focused fluid flow in an acid and redox‐influenced system: diagenetic controls on cement mineralogy and geomorphology in the Navajo Sandstone

Differential cement mineralogy is influenced by depositional textures, structural deformation, pore fluid chemistry, and ultimately influences landscape evolution by introducing heterogeneities in erodibility. In Southern Utah, the region West of the Kaibab uplift known as Mollies Nipple (Mollies) in Grand Staircase‐Escalante National Monument exhibits a complex history of fluid–sediment interactions, which has resulted in a localized zone of anomalous diagenetic iron sulfate (jarosite) mineralogy in a well‐exposed dune–interdune deposit within the Navajo Sandstone. Mineralogy and geochemistry of cements within this region are examined using reflectance and imaging spectroscopy, field investigations, microscopy, and whole‐rock geochemical analyses. These data show that the in‐situ jarosite cement is localized to a plane along the highest ridge of the butte, providing an armor along with other secondary cements, which controls the butte's geomorphic evolution. The jarosite cement is associated with other mineralogies suggesting that the sulfate was one of the latest fluid‐related precipitates in the paragenetic sequence. It was preceded by a regional bleaching event, precipitation of clay cements, some localized concretionary iron oxide precipitation, and formation of deformation bands. At least one generation of dense iron oxide mineralization is associated with cataclastic brittle deformation predating the sulfate precipitation. Trace element geochemistry of cements shows certain metal oxide populations associated with extremely high (>2000 ppm) arsenic values. We interpret the combination of spatial mineral distribution, observed paragenetic sequence, and trace element geochemistry to suggest this region experienced acid sulfate diagenesis along fracture‐controlled fluid conduits related to weathering of proximal, unidentified, sulfides, or H₂S associated with deep source beds. Jarosite is highly soluble, and its presence suggests that abundant fluid flow has not occurred in this region since its formation. This terminal cement‐forming event must have occurred prior to sandstone exhumation and erosion to form the current extreme landscape at Mollies. This site exhibits the influence that fluid geochemistry, sedimentary mineralogy, and structural fabric have on geomorphic evolution.