The interplay of permeability and fluid properties as a first order control of heat transport, venting temperatures and venting salinities at mid-ocean ridge hydrothermal systems

While the fundamental influence of fluid properties on venting temperatures in mid-ocean ridge (MOR) hydrothermal systems is now well established, the potential interplay of fluid properties with permeability in controlling heat transfer, venting temperatures, and venting salinities has so far received little attention. A series of numerical simulations of fully transient fluid flow in a generic, across-axis model of a MOR with a heat input equivalent to magmatic supply at a spreading rate of 10 cm year⁻¹ shows a strong dependence of venting temperature and salinity on the system's permeability. At high permeability, venting temperatures are low because fluid fluxes are so high that the basal conductive heating cannot warm the large fluid masses rapidly enough. The highest venting temperature around 400°C as well as sub-seafloor fluid phase separation occur when the permeability is just high enough that the fluid flux can still accommodate all heat input for advection, or for lower permeabilities where advection is no longer capable to transfer all incoming magmatic heat. In the latter case, additional mechanisms such as eruptions of basaltic magma may become relevant in balancing heat flow in MOR settings. The results can quantitatively be explained by the 'fluxibility' hypothesis of Jupp & Schultz (Nature, 403 , 2000, 880), which is used to derive diagrams for the relations between basal heat input, permeability and venting temperatures. Its predictive capabilities are tested against additional simulations. Potential implications of this work are that permeability in high-temperature MOR hydrothermal systems may be lower than previously thought and that low-temperature systems at high permeability may be an efficient way of removing heat at MOR.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

Factors affecting fluid flow in strike–slip fault systems: coupled deformation and fluid flow modelling with application to the western Mount Isa Inlier, Australia

Deformation and focused fluid flow within a mineralized system are critical in the genesis of hydrothermal ore deposits. Dilation and integrated fluid flux due to coupled deformation and fluid flow in simple strike–slip fault geometries were examined using finite difference analysis in three dimensions. A series of generic fault bend and fault jog geometries consistent with those seen in the western Mount Isa Inlier were modelled in order to understand how fault geometry parameters influence the dilation and integrated fluid flux. Fault dip, fault width, bend/jog angle, and length were varied, and a cross-cutting fault and contrasting rock types were included. The results demonstrate that low fault dips, the presence of contrasts in rock type, and wide faults produce highest dilation and integrated fluid flux values. Increasing fault bend lengths and angles increases dilation and integrated fluid flux, but increasing fault jog length or angle has the opposite effect. There is minimal difference between the outputs from the releasing and restraining fault bend and jog geometries. Model characteristics producing greater fluid flows and/or gradients can be used in a predictive capacity in order to focus exploration on regions with more favorable fault geometries, provided that the mineralized rocks had Mohr–Coulomb rheologies similar to the ones used in the models.     

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.     

Factors controlling free thermal convection in faults in sedimentary basins: implications for the formation of zinc–lead mineral deposits

Stratiform sediment‐hosted Zn–Pb–Ag mineral deposits constitute about 40% of the Earth's zinc resources ( Allen 2001 ), and in most cases their genesis involves the discharge of basinal brines near or on the seafloor through syndepositional faults ( Sangster 2002 ). From the point of view of base metal exploration, it is therefore essential to identify all possible faults that formerly carried the upwelling ore‐forming solutions during mineralising events. This paper presents a numerical investigation of the relative importance of various physical parameters in controlling fluid discharge, recharge and heat transport in faults. A two‐dimensional, free convection of pure water, hydrogeological model is developed for the McArthur basin in northern Australia based on the surface geology, known stratigraphic and structural relationships and regional geophysical interpretations. Numerical experiments and sensitivity analyses reveal that faults with strong initial heat input, due to depth of penetration or magmatic activity, are the most likely candidates to carry discharge fluids to the sites of metal precipitation. Deeper, wider and more permeable faults are more likely to behave as the fluid discharge pathways, whereas shallow, narrow or less permeable faults act as marine water recharge pathways. Compared with these fault‐related factors, aquifer physical properties are less important in determining fluid flow patterns and the geothermal regime. These results are an important step in understanding hydrothermal fluid flow in sedimentary basins in order to develop effective exploration criteria for the location of stratiform Zn–Pb–Ag deposits.     

Deep geothermal groundwater flow in the Seferihisar–Balçova area,Turkey: results from transient numerical simulations of coupled fluid flow and heat transport processes

The Seferihisar–Balçova Geothermal system (SBG) is characterized by complex temperature and hydrochemical anomalies. Previous geophysical and hydrochemical investigations suggest that hydrothermal convection in the faulted areas of the SBG and recharge flow from the Horst may be responsible for the observed patterns. A numerical model of coupled fluid flow and heat transport processes has been built in order to study the possible fluid dynamics of deep geothermal groundwater flow in the SBG. The results support the hypothesis derived from interpreted data. The simulated scenarios provide a better understanding of the geophysical conditions under which the different fluid dynamics develop. When recharge processes are weak, the convective patterns in the faults can expand to surrounding reservoir units or below the seafloor. These fault‐induced drag forces can cause natural seawater intrusion. In the Melange of the Seferihisar Horst, the regional flow is modified by buoyant‐driven flow focused in the series of vertical faults. As a result, the main groundwater divide can shift. Sealing caprocks prevent fault‐induced cells from being overwhelmed by vigorous regional flow. In this case, over‐pressured, blind geothermal reservoirs form below the caprocks. Transient results showed that the front of rising hot waters in faults is unstable: the tip of the hydrothermal plumes can split and lead to periodical temperature oscillations. This phenomenon known as Taylor–Saffman fingering has been described in mid‐ocean ridge hydrothermal systems. Our findings suggest that this type of thermal pulsing can also develop in active, faulted geothermal systems. To some extent, the role of an impervious fault core on the flow patterns has been investigated. Although it is not possible to reproduce basin‐scale transport processes, this first attempt to model deep groundwater geothermal flow in the SBG qualitatively supported the interpreted data and described the different fluid dynamics of the basin. Geofluids (2010) 10 , 388–405     

Numerical modelling of 3D fluid flow and oxygen isotope exchange in fractured media: spatial distribution of isotope patterns

An understanding of fluid flow, mass transport and isotopic exchange in fractured rock is required to understand the origin of several geological processes including hydrothermal mineral deposits. The numerical model HydroGeoSphere simulates 3D advection, molecular diffusion, mechanical dispersion and isotopic exchange in a discretely fractured porous media, and can be used to better understand the processes of mass transport and isotopic exchange in fractured rocks. Study of ¹⁸O isopleth patterns for different types of fractures and fracture networks with a range of structural complexity and hydraulic properties shows that fracture properties and geometry control mass transport and isotopic exchange. The hydraulic properties, as well as the density, spacing, and connectivity of fractures determine the isotopic patterns. Asymmetries in the geometry of oxygen isotope patterns could be used to determine the direction of hydrothermal fluid flow.     

Rates of retrograde metamorphism and their implications for crustal rheology

The rate of reaction of a natural hornblende garnet granulite with water under a range of mid- to lower crustal conditions has been investigated experimentally. In runs of between 7 and 84 days small but measurable amounts of water were consumed, and sheet silicates (300°C, 300 MPa and 400°C, 400 MPa) and/or secondary actinolite (400°C, 400 MPa and 500°C, 500 MPa) were observed to have grown. When normalized to the surface area of the starting materials, hydration rates were in the range of 2–5 × 10⁻⁸ g m⁻² sec⁻¹. These reaction rates imply that a film of water that infiltrated a planar crack with a half width of 100 μm would be completely consumed within c. 100 years. These results imply that where water infiltrates the crust along faults or underlying shear zones in response to a deformation, it will remain as a free phase for only a finite period of time, which in some cases will be less than the repeat time for major earthquakes in the fault system. Thus, the rheology of fault zones and shear zones is likely to be cyclical, with the zone becoming stronger with time as water is consumed, and then weakened by infiltration of water after each rupture.     

The application of failure mode diagrams for exploring the roles of fluid pressure and stress states in controlling styles of fracture-controlled permeability enhancement in faults and shear zones

Permeability enhancement associated with deformation processes in faults and shear zones plays a key role in facilitating fluid redistribution between fluid reservoirs in the crust. Especially in high fluid flux hydrothermal systems, fracture-controlled permeability can be relatively short-lived, unless it is repeatedly regenerated by ongoing deformation. Failure mode diagrams in pore fluid factor and differential stress space, here termed λ–σ failure mode diagrams, provide a powerful tool for analysing how fluid pressure and stress states drive failure, associated permeability enhancement and vein styles during deformation in faults and shear zones. During fault-valve behaviour in the seismogenic regime, relative rates of recovery of pore fluid factor, differential stress and fault cohesive strength between rupture events impact on styles of veining and associated, fracture-controlled permeability enhancement in faults and shear zones. Examples of vein-rich fault zones are used to illustrate how constraints can be placed, not just on fluid pressure and stress states at failure, but also on the fluid pressurization and loading paths associated with failure and transitory permeability enhancement in faults and shear zones. This provides insights about when, during the fault-valve cycle, various types of veins can form. The use of failure mode diagrams also provides insights about the relative roles of optimally oriented faults and misoriented faults as hydraulically conductive structures. The analysis highlights the dynamics of competition between fluid pressures and loading rates in driving failure and repeated permeability regeneration in fracture-controlled, hydrothermal systems.     

Physical properties of carbonate fault rocks,fucino basin (Central Italy): implications for fault seal in platform carbonates

We documented the porosity, permeability, pore geometry, pore type, textural anisotropy, and capillary pressure of carbonate rock samples collected along basin‐bounding normal faults in central Italy. The study samples consist of one Mesozoic platform carbonate host rock with low porosity and permeability, four fractured host rocks of the damage zones, and four fault rocks of the fault cores. The four fractured samples have high secondary porosity, due to elongated, connected, soft pores that provide fluid pathways in the damage zone. We modeled this zone as an elastic cracked medium, and used the Budiansky–O'Connell correlation to compute its permeability from the measured elastic moduli. This correlation can be applied only to fractured rocks with large secondary porosity and high‐aspect ratio pores. The four fault rock samples are made up of survivor clasts embedded in fine carbonate matrices and cements with sub‐spherical, stiff pores. The low porosity and permeability of these rocks, and their high values of capillary pressure, are consistent with the fault core sealing as much as 77 and 140 m of gas and oil columns, respectively. We modeled the fault core as a granular medium, and used the Kozeny–Carmen correlation, assigning the value of 5 to the Kozeny constant, to compute its permeability from the measured porosities and pore radii. The permeability structure of the normal faults is composed of two main units with unique hydraulic characteristics: a granular fault core that acts as a seal to cross‐fault fluid flow, and an elastic cracked damage zone that surrounds the core and forms a conduit for fluid flow. Transient pathways for along‐fault fluid flow may form in the fault core during seismic faulting due to the formation of opening‐mode fractures within the cemented fault rocks.     

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

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

Hydrothermal,multiphase convection of H2O‐NaCl fluids from ambient to magmatic temperatures: a new numerical scheme and benchmarks for code comparison

Thermohaline convection of subsurface fluids strongly influences heat and mass fluxes within the Earth's crust. The most effective hydrothermal systems develop in the vicinity of magmatic activity and can be important for geothermal energy production and ore formation. As most parts of these systems are inaccessible to direct observations, numerical simulations are necessary to understand and characterize fluid flow. Here, we present a new numerical scheme for thermohaline convection based on the control volume finite element method (CVFEM), allowing for unstructured meshes, the representation of sharp thermal and solute fronts in advection‐dominated systems and phase separation of variably miscible, compressible fluids. The model is an implementation of the Complex Systems Modelling Platform CSMP++ and includes an accurate thermodynamic representation of strongly nonlinear fluid properties of salt water for magmatic‐hydrothermal conditions (up to 1000°C, 500 MPa and 100 wt% NaCl). The method ensures that all fluid properties are taken as calculated on the respective node using a fully upstream‐weighted approach, which greatly increases the stability of the numerical scheme. We compare results from our model with two well‐established codes, HYDROTHERM and TOUGH2, by conducting benchmarks of different complexity and find good to excellent agreement in the temporal and spatial evolution of the hydrothermal systems. In a simulation with high‐temperature, high‐salinity conditions currently outside of the range of both HYDROTHERM and TOUGH2, we show the significance of the formation of a solid halite phase, which introduces heterogeneity. Results suggest that salt added by magmatic degassing is not easily vented or accommodated within the crust and can result in dynamic, complex hydrologies.     

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.     

What do aqueous geothermometers really tell us?

The application of chemical geothermometry to shallow groundwaters or spring discharge assumes that there is minimal mixing or re-equilibration of water as it travels from depth to the surface. In this study, we examine the potential for mixing and re-equilibration by examining heat and fluid flow along crustal-scale faults in tectonic geothermal systems. Numerical modeling results indicate that maximum in situ temperatures could be under-predicted by up to 30% due to mixing of fluids that enter the fault at different depths. This, coupled with the depression of isotherms by downward groundwater flow in the hanging wall, could cause underestimates of maximum circulation depth of greater than 80% in extreme cases. Kinetics does not favor re-equilibration in the shallower portions of faults due to low temperatures and higher fluid velocities. However, in areas of deeper circulation or higher heat flow such reactions are possible.     

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.     

Migration of Mn-rich fluids through normal faults and fine-grained terrigenous sediments during early development of the Neogene Vallès-Penedès half-graben (NE Spain)

The Miocene siliciclastic sediments infilling the Vallès‐Penedès half‐graben are affected by two sets of structures developed during the extensional tectonics that created the basin. The first set, represented by extension fractures infilled with mud and sands, is attributed to seismically induced liquefaction. The second set, represented by normal faults, corresponds to a high‐permeability horsetail extensional fracture mesh developed near the surface in the hanging walls of normal faults. The incremental character of the vein‐fills indicates episodic changes in the tectonic stress state and fault zone permeability. Two episodes of fluid migration are recorded. The first episode occurred prior to consolidation and lithification when shallow burial conditions allowed oxidizing meteoric waters to flow horizontally through the more porous and permeable sandy layers. Development of clastic dikes allowed local upward flow and dewatering of the sandy beds. Liquefaction and expulsion of fluids were probably driven by seismic shaking. During the first episode of fluid migration there was no cementation of the sandstone or within the fractures, probably because little fluid was mobilized by the predominantly compaction‐driven flow regime. The second episode of fluid migration occurred synchronously with normal fault development, during which time the faults acted as fluid conduits. Fluids enriched in manganese, probably leached from local manganese oxyhydroxides soon after sedimentation, moved laterally and produced cementation in the sandstone layers, eventually arriving at the more porous and permeable fault pathways that connected compartments of different porosities and permeabilities. Carbonate probably precipitated in fractures saturated with meteoric water near the ground surface at a transitional redox potential. Once the faults became occluded by calcite cement, shortly after fault development, they became barriers to both vertical and horizontal fluid flow.     

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.     

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.     

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.     

Seismic evidence for fluid migration pathways from an overpressured system in the South China Sea

Overpressured systems and intense, anomalously hot fluid expulsion in the Yinggehai Basin of the South China Sea offer an opportunity to understand the history of fluid flow and the process of hydrocarbon accumulation in overpressured environments. Fluid migration pathways from overpressured compartments in the basin are largely controlled by the distribution of faults and fractures. Episodic opening of these faults are related to the dynamics of an overpressured system and tectonic movements during basin evolution. At the crests of diapiric structures, fluid expulsion is seismically imaged as chimney‐ or plume‐like features, low to middle seismic amplitudes, and intermittently chaotic and blank reflecting seismic facies. These fluid pathways are controlled by vertical faults, which commonly penetrate overpressured and overlying normally pressured zones. Fluid expulsion is also observed near the main faults, such as the No. 1 Fault at the north‐eastern margin of the basin. Investigation by sidescan sonar on onshore and offshore Hainan Island indicates that there are more than 100 gas seepages adjacent to the No. 1 Fault. Migration pathways in the diapiric structures are controlled by three types of fault and fracture. Penetrative faults formed by dextral strike‐slip movement of the Red River faults commonly occur in the centre of the diapirs, and may have been a triggering factor for the diapirism, and controlled their distribution. Hydrofractures occur in certain mud‐rich layers and may have been generated by hydraulic fracturing. Radial normal faults occur at the top of diapirs and were formed by the intrusive process. These fluid migration pathways played an important role in regional hydrocarbon accumulation.