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.     

Changes in chemical composition caused by water–rock interactions across a strike‐slip fault zone: case study of the Atera Fault,Central Japan

In some cases, water–rock interactions in fault zones can affect radionuclide migration. Here, we analyzed the chemical compositions of well‐exposed fault rocks from the strike‐slip Atera Fault, Central Japan, in order to understand the variability and behavior of major and selected trace elements. The fault zone has a 1.2‐m‐wide, smectite‐rich fault core and paired damage zones that developed within welded tuff on one side of the core and within granite on the other side. The 30‐cm‐wide, kaolinite‐rich fault gouge is developed in granite cataclasite, and it shows indications of the latest fault activity, while the 1.2‐m‐wide fault core appears to be older. Hydrogen and oxygen isotope ratios in the clay‐rich fault gouges, and carbon and oxygen isotope ratios in carbonates indicate that the two major clay‐rich zones formed in bedrock near the surface, consistent with observed deformation structures. Based on chemical analyses, we identified (1) a slight depletion in SiO₂, Na₂O, K₂O, and light rare earth elements at the edges of the 1.2‐m‐wide fault core, (2) a clear depletion in SiO₂, Na₂O, K₂O, and all rare earth elements except Eu in the 30‐cm‐wide fault gouge, and (3) an increase in CaO, MnO, and heavy rare earth elements across the entire 1.2‐m‐wide fault core. Findings (1) and (2) reflect water–rock interactions in the 1.2‐m‐wide fault core and in the 30‐cm‐wide fault gouge that resulted in the formation of smectite and kaolinite. Finding (3) reflects carbonate precipitation caused by the addition of basalt fragments from a nearby site to the 1.2‐m‐wide fault core during faulting, and subsequent sorption reactions of heavy rare earth elements via processes such as complexation with the carbonates.     

Evolution of porosity, permeability and fluid pressure in dilatant faults post-failure: implications for fluid flow and mineralization

Mineralization of brittle fault zones is associated with sudden dilation, and the corresponding changes in porosity, permeability and fluid pressure, that occur during fault slip events. The resulting fluid pressure gradients cause fluid to flow into and along the fault until it is sealed. The volume of fluid that can pass through the deforming region depends on the degree of dilation, the porosity and permeability of the fault and wall rocks, and the rate of fault sealing. A numerical model representing a steep fault cutting through a horizontal seal is used to investigate patterns of fluid flow following a dilatant fault slip event. The model is initialized with porosity, permeability and fluid pressure representing the static mechanical state of the system immediately after such an event. Fault sealing is represented by a specified evolution of porosity, coupled to changes in permeability and fluid pressure, with the rate of porosity reduction being constrained by independent estimates of the rate of fault sealing by pressure solution. The general pattern of fluid flow predicted by the model is of initial flow into the fault from all directions, followed by upward flow driven by overpressure beneath the seal. The integrated fluid flux through the fault after a single failure event is insufficient to account for observed mineralization in faults; mineralization would require multiple fault slip events. Downward flow is predicted if the wall rocks below the seal are less permeable than those above. This phenomenon could at least partially explain the occurrence of uranium deposits in reactivated basement faults that cross an unconformity between relatively impermeable basement and overlying sedimentary rocks.     

Tectonic evolution of a Paleozoic thrust fault influences the hydrogeology of a fractured rock aquifer,northeastern Appalachian foreland

In polyorogenic regions, the superposition of structures during a protracted tectonic history produces complex fractured bedrock aquifers. Thrust‐faulted regions, in particular, have complicated permeability patterns that affect groundwater flow paths, quantity, and quality. In the Appalachian foreland of northwestern Vermont, numerous bedrock wells that are spatially related to the Paleozoic Hinesburg thrust have elevated naturally occurring radioactivity and/or low yields. The association of groundwater quality and quantity issues with this thrust was a unique opportunity to investigate its structural and hydrogeologic framework. The Hinesburg thrust juxtaposed metamorphic rocks of the hanging wall with sedimentary rocks of the footwall during the Ordovician. It was then deformed by two orthogonal Devonian fold sets and was fractured during the Cretaceous. Median well yields in the hanging wall aquifer are significantly lower than those of the footwall aquifer, consistent with the respective permeability contrast between metamorphic and carbonate rocks. For wells drilled through the Hinesburg thrust, those completed closest (vertically) to the thrust have the highest median yields, whereas others completed farther below have yields in the footwall range. The geochemical signature of the hanging wall and footwall aquifers correlates with their whole‐rock geochemistry. The hanging wall aquifer is enriched in alpha radiation, Na+K‐Cl, Ba, and Sr, whereas the footwall aquifer is enriched in Ca‐Mg‐HCO₃ and alkalinity. Wells that penetrated the Hinesburg thrust generally have hanging wall geochemical signatures. A simple hydrogeologic model for the permeability evolution of the Hinesburg thrust involves the ductile emplacement of a low‐K hanging wall onto a high‐K footwall, with subsequent modification by fractures.     

Weakening processes in thrust faults: insights from the Monte Perdido thrust fault (southern Pyrenees,Spain)

The Monte Perdido thrust fault (southern Pyrenees) consists of a 6‐m‐thick interval of intensely deformed clay‐bearing rocks. The fault zone is affected by a pervasive pressure solution seam and numerous shear surfaces. Calcite extensional‐shear veins are present along the shear surfaces. The angular relationships between the two structures indicate that shear surfaces developed at a high angle (70°) to the local principal maximum stress axis σ₁. Two main stages of deformation are present. The first stage corresponds to the development of calcite shear veins by a combination of shear surface reactivation and extensional mode I rupture. The second stage of deformation corresponds to chlorite precipitation along the previously reactivated shear surfaces. The pore fluid factor λ computed for the two deformation episodes indicates high fluid pressures during the Monte Perdido thrust activity. During the first stage of deformation, the reactivation of the shear surface was facilitated by a suprahydrostatic fluid pressure with a pore fluid factor λ equal to 0.89. For the second stage, the fluid pressure remained still high (with a λ value ranging between 0.77 and 0.84) even with the presence of weak chlorite along the shear surfaces. Furthermore, evidence of hydrostatic fluid pressure during calcite cement precipitation supports that incremental shear surface reactivations are correlated with cyclic fluid pressure fluctuations consistent with a fault‐valve model.     

Isotopic and petrological evidence of fluid–rock interaction at a Tethyan ocean–continent transition in the Alps: implications for tectonic processes and carbon transfer during early ocean formation

We report overprinting stable isotope evidence of fluid–rock interaction below two detachment faults along which mantle rocks were exhumed to the seafloor, between the respective landward and seaward limits of oceanic and continental crust, at a Tethyan ocean–continent transition (OCT). This OCT, which is presently exposed in the Tasna nappe (south‐eastern Switzerland) is considered an on‐land analogue of the well‐studied Iberian OCT. We compare our results with the fault architecture (fault core–damage zone–protolith) described by Caine et al. [Geology (1996) Vol. 24, pp. 1025–1028]. We confirm the existence of a sharp boundary between the fault core and damage zone based on isotopic data, but the boundary between the damage zone and protolith is gradational. We identify evidence for: (1) pervasive isotopic modification to 8.4 ± 0.1‰ which accompanied or post‐dated serpentinization of these mantle rocks at an estimated temperature of 67–109°C, (2) either (i) partial isolation of some highly strained regions [fault core(s) and mylonite] from this pervasive isotopic modification, because of permeability reduction (Caine et al.) or (ii) subsequent isotopic modification caused by structurally channelled flow of warm fluids within these highly strained regions, because of permeability enhancement, and (3) isotopic modification, which is associated with extensive calcification at T = 54–100°C, primarily beneath the younger of the two detachment faults and post‐dating initial serpentinization. By comparing the volumetric extent of calcification with an experimentally verified model for calcite precipitation in veins, we conclude that calcification could have occurred in response to seawater infiltration, with a calculated flux rate of 0.1–0.2 m year^?1 and a minimum duration of 0.2–4.0 × 10⁴ years. The associated time‐averaged uptake flux of carbon during this period was 8–120 mol m^?2 year^?1. By comparison with the estimated area of exhumed mantle rocks at the Iberian OCT, we calculate a maximum annual uptake flux for carbon of 2–30 Tg year^?1. This is an order of magnitude greater than that for carbon exchange at the mid‐ocean ridges and 0.1–1.4% of the global oceanic uptake flux for carbon.     

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.     

Analytical and numerical models of hydrothermal fluid flow at fault intersections

Fault intersections are the locus of hot spring activity and Carlin‐type gold mineralization within the Basin and Range, USA. Analytical and numerical solutions to Stokes equation suggest that peak fluid velocities at fault intersections increase between 20% and 47% when fracture apertures have identical widths but increase by only about 1% and 8% when aperture widths vary by a factor of 2. This suggests that fault zone intersections must have enlarged apertures. Three‐dimensional finite element models that consider intersecting 10‐ to 20‐m wide fault planes resulted in hot spring activity being preferentially located at fault zone intersections when fault zones were assigned identical permeabilities. We found that the onset of convection at the intersections of the fault zones occurred in our hydrothermal model over a narrow permeability range between 5 × 10^?13 and 7 × 10^?13 m². Relatively high vertical fluid velocities (0.3–3 m year^?1) extended away from the fault intersections for about 0.5–1.5 km. For the boundary conditions and fault plane dimensions used, peak discharge temperatures of 112°C at the water table occurred with an intermediate fault zone permeability of 5 × 10^?13 m². When fault plane permeability differed by a factor of 2 or more, the locus of hot spring activity shifted away from the intersections. However, increasing the permeability at the core of the fault plane intersection by 40% shifted the discharge back to the intersections. When aquifer units were assigned a permeability value equal to those of the fault planes, convective rolls developed that extend about 3 km laterally along the fault plane and into the adjacent aquifer.     

Hydrothermal control on organic matter alteration and illite precipitation, Mt Isa Basin, Australia

Abundant illite precipitation in Proterozoic rocks from Northern Lawn Hill Platform, Mt Isa Basin, Australia, occurred in organic matter‐rich black shales rather than in sandstones, siltstones and organic matter‐poor shales. Sandstones and siltstones acted as impermeable rocks, as early diagenetic quartz and carbonate minerals reduced the porosity–permeability. Scanning and transmission electron microscopy (SEM and TEM) studies indicate a relation between creation of microporosity–permeability and organic matter alteration, suitable for subsequent mineral precipitation. K–Ar data indicate that organic matter alteration and the subsequent illite precipitation within the organic matter occurred during the regional hydrothermal event at 1172 ± 50 (2σ) Ma. Hot circulating fluids are considered to be responsible for organic matter alteration, migration and removal of volatile hydrocarbon, and consequently porosity–permeability creation. Those rocks lacking sufficient porosity–permeability, such as sandstones, siltstones and organic matter poor shales, may not have been affected by fluid movement. In hydrothermal systems, shales and mudstones may not be impermeable as usually assumed because of hydrocarbons being rapidly removed by fluid, even with relatively low total organic carbon.     

The influence of poorly interconnected fault zone flow paths on spring geochemistry

Thermal springs commonly occur along faults because of the enhanced vertical permeability afforded by fracture zones. Field and laboratory studies of fault zone materials document substantial heterogeneities in fracture permeabilities. Modeling and field studies of springs suggest that spatial variations in permeability strongly influence spring locations, discharge rates and temperatures. The impact of heterogeneous permeability on spring geochemistry, however, is poorly documented. We present stable isotope and water chemistry data from a series of closely spaced thermal springs associated with the Hayward Fault, California. We suggest that substantial spatial variations observed in δ¹⁸O and chloride values reflect subsurface fluid transport through a poorly connected fracture network in which mixing of subsurface waters remains limited. Our measurements provide insight into the effect of fracture zone heterogeneities on spring geochemistry, offer an additional tool to intuit the nature of tectonically induced changes in fault zone plumbing, and highlight the need to consider local variations when characterizing fracture zone fluid geochemistry from spring systems with multiple discharge sites.     

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.     

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.     

Fluid evolution along multistage composite fault systems at the southern margin of the Lower Palaeozoic Anglo-Brabant fold belt, Belgium

Mineralised vein systems have been investigated at nine localities at the southern margin of the Anglo‐Brabant fold belt in Belgium. During the late Silurian to early Middle Devonian Caledonian orogeny, shear zones formed, inferred to be associated with granitoid basement blocks in the subsurface. The circulation of a metamorphic fluid, possibly originating in the Cambrian core of the fold belt, along these shear zones resulted in the formation of mesozonal orogenic mineralisation at the southern margin of the Anglo‐Brabant fold belt. The fluid had a composition dominated by H₂O–CO₂–X–NaCl–KCl. The shear zones form part of a greater fault zone, the Nieuwpoort–Asquempont fault zone, which is characterised by normal faulting that started before the Givetian and by the reactivation of the shear zones. Two fluid generations are associated with this normal faulting. First, a low salinity H₂O–NaCl(–KCl) fluid migrated through the Palaeozoic rocks after the Silurian. Based on the isotopic composition, this fluid could be a late‐metamorphic Caledonian fluid or a younger fluid that originated from the Rhenohercynian basin and interacted with Lower Devonian rocks along its migration path. Second, a high salinity H₂O–NaCl–CaCl₂ fluid was identified in the fault systems. Similar fluids have been found in southern and eastern Belgium, where they produced Mississippi Valley‐type Zn–Pb deposits. These fluids are interpreted as evaporative brines that infiltrated the Lower Palaeozoic basement, from where they were expelled during extensional tectonism in the Mesozoic.     

Hydrothermal convection in and around mineralized fault zones: insights from two- and three-dimensional numerical modeling applied to the Ashanti belt, Ghana

Among hydrogeological processes, free convection in faults has been cited as a possible cause of gold mineralization along major fault zones. Here, we investigate the effects of free convection to determine whether it can cause giant orogenic gold deposits and their regular spatial distribution along major fault/shear zones. The approach comprises: (i) coupled two- and three-dimensional numerical heat- and fluid-flow simulations of simplified geological models; and (ii) calculation of the rock alteration index (RAI) to delineate regions where precipitation/dissolution can occur. Then, comparing the deduced alteration patterns with temperature distribution, potential areas of gold mineralization, defined by T  > 200°C and RAI < 0, are predicted. The models are based on the orogenic Paleoproterozoic ore deposits of the Ashanti belt in western Africa. These deposits occur in the most permeable parts of the fault zone, where the lateral permeability contrast is the highest. For a simple geometry, with a fault zone adjacent to a sedimentary basin half as permeable, we note a transition from three-dimensional circulation within the fault to a two-dimensional convective pattern in the basin far from the fault. Moreover, whereas two-dimensional undulated isotherms dominate in the basin, three-dimensional corrugated isotherms result from the preferred convective pattern within the fault, thus enhancing a periodic distribution of thermal highs and lows. In our most elaborate three-dimensional model with an imposed lateral permeability gradient, the RAI distribution indicates that fluid circulation in fault zones gives rise to a spatial periodicity of alteration patterns consistent with field data.     

Polyphasic hydrothermal and meteoric fluid regimes during the growth of a segmented fault involving crystalline and carbonate rocks (Barcelona Plain,NE Spain)

A polyphasic tectonic‐fluid system of a fault that involves crystalline and carbonate rocks (Hospital fault, Barcelona Plain) has been inferred from regional to thin section scale observations combined with geochemical analyses. Cathodoluminescence, microprobe analyses and stable isotopy in fracture‐related cements record the circulation of successive alternations of hydrothermal and low‐temperature meteoric fluids linked with three main regional tectonic events. The first event corresponds to the Mesozoic extension, which had two rifting stages, and it is characterized by the independent tectonic activity of two fault segments, namely southern and northern Hospital fault segments. During the Late Permian‐Middle Jurassic rifting, these segments controlled the thickness and distribution of the Triassic sediments. Also, dolomitization was produced in an early stage by Triassic seawater at shallow conditions. During increasing burial, formation of fractures and their dolomite‐related cements took place. Fault activity during the Middle Jurassic–Late Cretaceous rifting was localized in the southern segment, and it was characterized by hydrothermal brines, with temperatures over 180°C, which ascended through this fault segment precipitating quartz, chlorite, and calcite. The second event corresponds to the Paleogene compression (Chattian), which produced exhumation, folding and erosion, favouring the percolation of low‐temperature meteoric fluids which produced the calcitization of the dolostones and of the dolomite cements. The third event is linked with the Neogene extension, where three stages have been identified. During the syn‐rift stage, the southern segment of the Hospital fault grew by tip propagation. In the relay zone, hydrothermal brines with temperature around 140°C upflowed. During the late postrift, the Hospital fault acted as a unique segment and deformation occurred at shallow conditions and under a low‐temperature meteoric regime. Finally, and possibly during the Messinian compression, NW‐SE strike‐slip faults offset the Hospital fault to its current configuration.     

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.     

Numerical modeling of the origin of calcite mineralization in the Refugio‐Carneros fault,Santa Barbara Basin,California

Many faults in active and exhumed hydrocarbon‐generating basins are characterized by thick deposits of carbonate fault cement of limited vertical and horizontal extent. Based on fluid inclusion and stable isotope characteristics, these deposits have been attributed to upward flow of formation water and hydrocarbons. The present study sought to test this hypothesis by using numerical reactive transport modeling to investigate the origin of calcite cements in the Refugio‐Carneros fault located on the northern flank of the Santa Barbara Basin of southern California. Previous research has shown this calcite to have low δ¹³C values of about ?40 to ?30‰PDB, suggesting that methane‐rich fluids ascended the fault and contributed carbon for the mineralization. Fluid inclusion homogenization temperatures of 80–125°C in the calcite indicate that the fluids also transported significant quantities of heat. Fluid inclusion salinities ranging from fresh water to seawater values and the proximity of the Refugio‐Carneros fault to a zone of groundwater recharge in the Santa Ynez Mountains suggest that calcite precipitation in the fault may have been induced by the oxidation of methane‐rich basinal fluids by infiltrating meteoric fluids descending steeply dipping sedimentary layers on the northern basin flank. This oxidation could have occurred via at least two different mixing scenarios. In the first, overpressures in the central part of the basin may have driven methane‐rich formation waters derived from the Monterey Formation northward toward the basin flanks where they mixed with meteoric water descending from the Santa Ynez Mountains and diverted upward through the Refugio‐Carneros fault. In the second scenario, methane‐rich fluids sourced from deeper Paleogene sediments would have been driven upward by overpressures generated in the fault zones because of deformation, pressure solution, and flow, and released during fault rupture, ultimately mixing with meteoric water at shallow depth. The models in the present study were designed to test this second scenario, and show that in order for the observed fluid inclusion temperatures to be reached within 200 m of the surface, moderate overpressures and high permeabilities were required in the fault zone. Sudden release of overpressure may have been triggered by earthquakes and led to transient pulses of accelerated fluid flow and heat transport along faults, most likely on the order of tens to hundreds of years in duration. While the models also showed that methane‐rich fluids ascending the Refugio‐Carneros fault could be oxidized by meteoric water traversing the Vaqueros Sandstone to form calcite, they raised doubts about whether the length of time and the number of fault pulses needed for mineralization by the fault overpressuring mechanism were too high given existing geologic constraints.     

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.     

Precipitation of fracture fillings and cements in the Buntsandstein (NW Germany)

The relationship between fracturing and fracture filling in opening‐mode fractures in the Triassic Buntsandstein in the Lower Saxony Basin (LSB; NW Germany) has been studied by an integration of petrographic and structural analysis of core samples, strontium isotope analysis and microthermometry on fluid inclusions. This revealed the relationship between the timing of the fracturing and the precipitation of different mineral phases in the fractures by constraining the precipitation conditions and considering the possible fluid transport mechanisms. The core was studied from four different boreholes, located in different structural settings across the LSB. In the core samples from the four boreholes, fractures filled with calcite, quartz and anhydrite were found, in addition to pore‐filling calcite cementation. In boreholes 2 and 3, calcite‐filled fractures have a fibrous microstructure whereas in borehole 1, fractures are filled with elongate‐blocky calcite crystals. Anhydrite‐filled fractures have, in all samples, a blocky to elongate‐blocky microstructure. Fractures that are filled with quartz are observed in borehole 2 only where the quartz crystals are ‘stretched’ with an elongated habit. Fluid inclusion microthermometry of fracturing‐filling quartz crystals showed that quartz precipitation took place at temperatures of at least 140°C, from a fluid with NaCl–CaCl₂–H₂O composition. Melting phases are meta‐stable and suggest growth from high salinity formation water. Strontium isotopes, measured in leached host rock, indicate that, in boreholes 2 and 3, the fluid which precipitated the calcite cements and calcite‐filled fractures is most likely locally derived whereas in borehole 1, the ⁸⁷Sr/⁸⁶Sr ratios from the pore‐filling cements and in the elongate‐blocky calcite‐filled fracture can only be explained by mixing with externally derived fluids. The elongate‐blocky anhydrite‐filled fractures, present in boreholes 1, 3 and 4, precipitated from a mixture of locally derived pore fluids and a significant quantity of fluid with a lower, less radiogenic, ⁸⁷Sr/⁸⁶Sr ratio. Taking into account the structural evolution of the basin and accompanying salt tectonics, it is likely that the underlying Zechstein is a source for the less radiogenic fluids. Based on the samples in the LSB, it is probable that fibrous fracture fillings in sedimentary rocks most likely developed from locally derived pore fluids whereas elongate‐blocky fracture fillings with smooth walls developed from externally derived pore fluids.     

Transient fluid flow in and around a fault

We present the results of simple numerical experiments in which we study the evolution with time of fluid flow around and within a permeable fault embedded in a less permeable porous medium. Fluid movement is driven by an imposed vertical pressure gradient. The results show that fluid flow is controlled by two timescales: τ_f = Sl²/κ_F and τ_F = Sl²/κ_M, where S is the specific storage of the porous material, l the length of the fault, and κ_M and κ_F are the hydraulic conductivities of the porous material and the fault, respectively. Fluid flow and the associated fluid pressure field evolve through three temporal stages: an early phase [t < τ_f] during which the initial fluid pressure gradient within the fault is relaxed; a second transient stage [τ_f < t < τ_F] when fluid is rapidly expelled at one end of the fault and extracted from the surrounding rocks at the other end leading to a reduction in the pressure gradient in the intact rock; a third phase [t < τ_F] characterized by a steady‐state flow. From the numerical experiments we derive an expression for the steady‐state maximum fluid velocity in the fault and the values of the two timescales, τ_f and τ_F. A comparison indicates excellent agreement of our results with existing asymptotic solutions. For km‐scale faults, the model results suggest that steady‐state is unlikely to be reached over geological timescales. Thus, the current use of parameters such as the focusing ratio defined under the assumption of steady‐state conditions should be reconsidered.