An exhumed palaeo-hydrocarbon migration fairway in a faulted carrier system, Entrada Sandstone of SE Utah, USA

The Moab Anticline, east‐central Utah, is an exhumed hydrocarbon palaeo‐reservoir which was supplied by hydrocarbons that migrated from the Moab Fault up‐dip towards the crest of the structure beneath the regional seal of the Tidwell mudstone. Iron oxide reduction in porous, high permeability aeolian sandstones records the secondary migration of hydrocarbons, filling of traps against small sealing faults and spill pathways through the Middle Jurassic Entrada Sandstone. Hydrocarbons entered the Entrada Sandstone carrier system from bends and other leak points on the Moab Fault producing discrete zones of reduction that extend for up to 400 m from these leak points. They then migrated in focused stringers, 2–5 m in height, to produce accumulations on the crest of the anticline. Normal faults on the anticline were transient permeability barriers to hydrocarbon migration producing a series of small compartmentalized accumulations. Exsolution of CO₂ as local fault seals were breached resulted in calcite cementation on the up‐dip side of faults. Field observations on the distribution of iron oxide reduction and calcite cements within the anticline indicate that the advancing reduction fronts were affected neither by individual slip bands in damage zones around faults nor by small faults with sand: sand juxtapositions. Faults with larger throws produced either sand: mudstone juxtapositions or sand: sand contacts and fault zones with shale smears. Shale‐smeared fault zones provided seals to the reducing fluid which filled the structural traps to spill points.     

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.     

Communication between overpressured compartments

Overpressure build‐up in compartments, and communication between overpressured compartments across faults are studied with simple analytical and numerical models. It is shown that the excess pressure in a (vertical) one‐dimensional, one‐compartment model can be written as the sum of the excess pressure generated in the seal above the compartment, and a second part, which is due to the expulsion of fluid from the compartment and the rocks below. The one‐compartment model is generalized to a two‐compartment model, which accounts for the fluid communication between the compartments through a fault zone. The volume rates of flow through the seals and the fault zone are shown to be the weighted mean of the volume rates of the one‐dimensional, one‐compartment model. The normalized weights are given by dimensionless numbers, called fault–seal numbers, which control the communication between the compartments. A fault–seal number much less than unity implies that the fault is a stronger barrier for the fluid flow than the seal. A fault–seal number larger than unity implies the opposite: that the seal is a stronger barrier than the fault. The conditions for isolated compartments and other regimes are identified in terms of the fault–seal numbers. It is discussed how the compartment fault–seal numbers can be computed when the permeability is given in the fault zone. The results given by the analytical compartment models are demonstrated and validated with two‐dimensional numerical (finite element) simulations.     

Modeling fluid flow in a heterogeneous, fault-controlled hydrothermal system

Previous studies have shown that most hydrothermal systems discharging at the land surface are associated with faulting, and that the location, temperature and rate of discharge of these systems are controlled by the geometry and style of the controlling fault(s). Unfortunately, the transport of heat and fluid in fault-controlled hydrothermal systems is difficult to model realistically; although heterogeneity and anisotropy are assumed to place important controls on flow in faults, few data or observations are available to constrain the distribution of hydraulic properties within active faults. Here, analytical and numerical models are combined with geostatistical models of spatially varying hydraulic properties to model the flow of heat and fluid in the Borax Lake fault of south-east Oregon, USA. A geometric mean permeability within the fault of 7 × 10⁻¹⁴ m² with 2× vertical/horizontal anisotropy in correlation length scale is shown to give the closest match to field observations. Furthermore, the simulations demonstrate that continuity of flow paths is an important factor in reproducing the observed behavior. In addition to providing some insight into possible spatial distributions of hydraulic properties at the Borax Lake site, the study highlights one potential avenue for integrating field observations with simulation results in order to gain greater understanding of fluid flow in faults and fault-controlled hydrothermal and petroleum reservoirs.     

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.     

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.     

Capillary seal capacity of faults under hydrodynamic conditions

Many fault bound traps are underfilled despite the top seal capacity being secure. The hydrocarbon sealing performance of faults themselves can be compromised either by mechanical or capillary process. Capillary process can be important either due to juxtaposition or to fine‐grained clay or cataclastic material within the fault zone itself. There is debate about how important each of these mechanisms is over geological timescales of hydrocarbon trapping. Recent work has provided insights into fine‐tuning capillary‐related fault seal calibration methodologies. Over the last 15 years, vigorous scientific debate with multiple published laboratory experiments and modelling studies has led some researchers and industry technologists to theorise that for water‐wet conventional hydrocarbon reservoirs, the relative water permeability in the reservoir (towards the top of the hydrocarbon column) may become very small, but in practice never reach zero. While not advocating for either side in this debate, the importance of accounting for hydrodynamic conditions regardless of the capillary sealing mechanism is demonstrated. Additionally, it is noted that nonzero relative water permeability has implications on how a seal's capillary threshold pressure for the nonwetting hydrocarbon phase is estimated from field data. In the particular case where there are pressure differences between unproduced hydrocarbon reservoirs on either side of a fault, then the hydrocarbon saturation must be discontinuous across the fault. For hydrocarbon leakage to occur across the entire thickness of the fault zone, the hydrocarbon pressure must exceed the threshold pressure on the side of the fault zone with the highest formation water hydraulic head. This approach to estimating across‐fault pressure difference will result in an improved calibration data set used for predrill estimation of capillary fault seal capacity.     

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.     

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.     

Influence of fault map resolution on pore pressure distribution and secondary hydrocarbon migration; Tune area,North Sea

Pressure and hydrocarbon migration modelling was carried out in the Tune Field area, Viking Graben, offshore Norway. The pressures are considered to be controlled by compartments bounded by mapped faults. Two different interpreted fault maps at the top reservoir level (Brent Group) are used as input to the modelling. First, a low‐resolution fault map is used, with only the large faults interpreted, and next, both large and small faults are included. The simulations show high overpressures generated in the western area, in the deeper part of the Viking Graben, and hydrostatic in the eastern areas. A sharp transition zone results from using the low‐resolution fault map in the simulations. Small N–S striking faults situated in between the wells have to have higher sealing capacity than expected from juxtaposition analysis alone, to be able to match the overpressures measured in well 30/5‐2 and 30/8‐1S in the Tune Field, and well 30/8‐3 east of Tune. The intermediate pressure in the western part is probably related to flow in the deeper parts of the sedimentary column in the compartment, where well 30/8‐3 is situated. The secondary oil migration models show that overpressures have major effects on the migration pathways of hydrocarbons. The level of detail in the fault interpretation is important for simulation results, both for pressure distribution and for hydrocarbon migration.     

Development of Neural Network Based Hysteretic Models for Steel Beam-Column Connections Through Self-Learning Simulation

Beam-column connections are zones of highly complex actions and deformations interaction that often lead to failure under the effect of earthquake ground motion. Modeling of the beam-column connections is important both in understanding the behavior and in design. In this article, a framework for developing a neural network (NN) based steel beam-column connection model through structural testing is proposed. Neural network based inelastic hysteretic model for beam-column connections is combined with a new component based model under self-learning simulation framework. Self-learning simulation has the unique advantage in that it can use structural response to extract material models. Self-learning simulation is based on auto-progressive algorithm that employs the principles of equilibrium and compatibility, and the self-organizing nature of artificial neural network material models. The component based model is an assemblage of rigid body elements and spring elements which represent smeared constitutive behaviors of components; either nonlinear elastic or nonlinear inelastic behavior of components. The component based model is verified by a 3-D finite element analysis. The proposed methodology is illustrated through a self-learning simulation for a welded steel beam-column connection. In addition to presenting the first application of self-learning simulation to steel beam-column connections, a framework is outlined for applying the proposed methodology to other types of connections.     

Hydrodynamics and membrane seal capacity

The impact of hydrodynamic groundwater movement on the capacity of seals is currently in debate. There is an extensive record of publication on seals analysis and a similar history on petroleum hydrodynamics yet little work addresses the links between the two. Understanding and quantifying the effects of hydrodynamic flow has important implications for calibrating commonly used seal capacity estimation techniques. These are often based on measurements such as shale gouge, clay smear or mercury porosimitry where membrane sealing is thought to occur. For standard membrane seal analysis, seal capacity is estimated by quantifying capillary pressure‐related measurements and calibrating them with a large observational database of hydrocarbon column heights and measured buoyancy pressures. The seal capacity estimation process has historically been adjusted to account for a number of different generic trapping geometries. We define the characteristics of these geometries from a hydrodynamics viewpoint in order to fine‐tune the seal capacity calibration process. From theoretical analyses of several simplified trapping geometries, it can be concluded that generally, the high pressure side of the seal should be used as the water pressure gradient with which to calculate buoyancy pressure. Secondly, trap geometries where hydrocarbon is reservoired on both sides of a fault are not useful for estimating across fault seal capacity.     

Deformation bands and their impact on fluid flow in sandstone reservoirs: the role of natural thickness variations

Cataclastic deformation bands, which are common in sandstone reservoirs and which may negatively affect fluid flow, are generally associated with notable thickness variations. It has been suggested previously that such thickness variations represent an important control on how deformation bands affect fluid flow. The effects of such thickness variations are tested in this study though statistical analysis and fluid flow simulation of an array of cataclastic deformation bands in Cretaceous sandstones in in the Bassin de Sud‐Est in Provence, France. Spatial outcrop data are statistically analyzed for incorporation in flow simulation models, and numerical simulations are used to investigate the effects of notable thickness variations on how the deformation bands influence effective permeability and flow dynamics. A suite of simulations is performed using a combination of fine‐scale and coarse‐scale grids, revealing that the effective permeability of the simulated reservoir is reduced by a factor of 15–25. More interestingly, the simulations further demonstrated that, as compared to the overall effect of the deformation band array on fluid flow, thickness variations along the bands proved to have negligible effects only. Thus, our simulations indicate that the configuration and connectivity of the deformation bands, together with the permeability contrast between the bands and the host rock and the mean band thickness, are the most important controls on the effective permeability. Our findings represent new insight into the influence of deformation bands on fluid flow in subsurface aquifers and reservoirs, indicating that thickness variations of individual deformation bands are of less significance than previously thought.     

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.     

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.     

Upper crustal fluid flow in the outboard region of the Southern Alps,New Zealand

The currently active fluid regime within the outboard region of the Southern Alps, New Zealand was investigated using a combination of field observations, carbon‐ and oxygen‐stable isotopes from fault‐hosted calcites and interpretation of magnetotelluric (MT) data. Active faulting in the region is dominated by NE striking and N striking, oppositely dipping thrust fault pairs. Stable isotopic analyses of calcites hosted within these fault zones range from 10 to 25‰δ¹⁸O and from ?33 to 0‰δ¹³C. These values reflect mixing of three parent fluids: meteoric water, carbon‐exchanged groundwater and minor deeper rock‐exchanged fluids, at temperatures of 10–90°C in the upper 3–4 km of the crust. A broad, ‘U‐shaped’ zone of high electrical conductivity (maximum depth c. 28 km) underlies the central Southern Alps. In the ductile region of the crust, the high‐conductivity zone is subhorizontal. Near‐vertical zones of high conductivity extend upward to the surface on both sides of the conductive zone. On the outboard side of the orogen, the conductive zone reaches the surface coincident with the trace of the active Forest Creek Faults. Near‐surface flow is shown to dominate the outboard region. Topographically driven meteoric water interacts, on a kilometre scale, with either carbon‐exchanged groundwater or directly with organic material within Pliocene gravels, resulting in a distinctive low ¹³C signal within fault‐hosted calcites of the outboard region. The high‐strain zone in the lower crust focuses the migration of deeply sourced fluids upward to the base of the brittle–ductile transition. Interconnected fluid is imaged as a narrow vertical zone of high conductivity in the upper crust, implying continuous permeability and possibly buoyancy‐driven flow of deeply sourced fluids to higher levels of the crust where they are detected by the isotopic analysis of the fault‐hosted calcites.     

Linked thermal convection of the basement and basinal fluids in formation of the unconformity‐related uranium deposits in the Athabasca Basin,Saskatchewan, Canada

World‐class unconformity‐related U deposits in the Athabasca Basin (Saskatchewan, Canada) are generally located within or near fault zones that intersect the unconformity between the Athabasca Group sedimentary basin rocks and underlying metamorphic basement rocks. Two distinct subtypes of unconformity‐related uranium deposits have been identified: those hosted primarily in the Athabasca Group sandstones (sediment‐hosted) and those hosted primarily in the underlying basement rocks (basement‐hosted). Although significant research on these deposits has been carried out, certain aspects of their formation are still under discussion, one of the main issues being the fluid flow mechanisms responsible for uranium mineralization. The intriguing feature of this problem is that sediment‐hosted and basement‐hosted deposits are characterized by oppositely directed vectors of fluid flow via associated fault zones. Sediment‐hosted deposits formed via upward flow of basement fluids, basement‐hosted deposits via downward flow of basinal fluids. We have hypothesized that such flow patterns are indicative of the fluid flow self‐organization in fault‐bounded thermal convection (Transport in Porous Media, 110, 2015, 25). To explore this hypothesis, we constructed a simplified hydrogeologic model with fault‐bounded thermal convection of fluids in the faulted basement linked with fluid circulation in the overlying fault‐free sandstone horizon. Based on this model, a series of numerical experiments was carried out to simulate the hypothesized fluid flow patterns. The results obtained are in reasonable agreement with the concept of fault‐bounded convection cells as an explanation of focused upflow and downflow across the basement/sandstone unconformity. We then discuss application of the model to another debated problem, the uranium source for the ore‐forming basinal brines.     

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.     

Evaluation of recoverable energy potential from enhanced geothermal systems: a sensitivity analysis in a poro‐thermo‐elastic framework

This study presents application of an efficient approach to simulate fluid flow and heat transfer in naturally fractured geothermal reservoirs. Fluid flow is simulated by combining single continuum and discrete fracture approaches. The local thermal nonequilibrium approach is used to simulate heat transfer by conduction in the rock matrix and convection (including conduction) in the fluid. Fluid flow and heat transfer models are integrated within a coupled poro‐thermo‐elastic framework. The developed model is used to evaluate the long‐term response of a geothermal reservoir with specific boundary conditions and injection/production schedule. A comparative study and a sensitivity analysis are carried out to evaluate the capability of the integrated approach and understand the degree by which different reservoir parameters affect thermal depletion of Soultz geothermal reservoir, respectively. Also observed, there exists an optimum fracture permeability after which the reservoir stimulation does not change the recovery factor significantly. Estimation of fluid temperature by the assumption of local thermal nonequilibrium heat transfer between the fracture fluid and the rock matrix gives fluid temperature of about 3°C less than that of estimated by thermal equilibrium heat transfer at early stage of hot water production.     

Structural controls on hydrothermal flow in a segmented rift system, Taupo Volcanic Zone, New Zealand

Vigorous hydrothermal convection transfers 10 times the average continental heat flow through the central Taupo Volcanic Zone (TVZ), a region of active extension (approximately 8 mm year^?1) and productive rhyolitic volcanism. Over 20 high‐temperature (>250°C) geothermal fields occur within Quaternary pyroclastic basins, with convective circulation to depths of 7–8 km presumably extending through basement rocks. Parallel‐striking normal faults, fractures and dikes dissect the convective regime, interacting with fluids to either enhance or restrict flow according to the relative permeability of structure and host rock. In the basement, high bulk permeability is maintained by focussed flow through faults and associated fractures well oriented for reactivation in the prevailing stress field. In contrast, distributed flow through fault‐bounded compartments prevails within Quaternary basins, masking any signal of deeper structural control. Exceptions occur where more competent rocks are exposed at the surface. As in narrow magmatic rifts elsewhere, the extensional fabric is partitioned into discrete rift segments linked along strike by accommodation zones. Eighty per cent of TVZ geothermal fields correlate spatially with rift architecture, with 60% located in accommodation zones. We suggest that segmented rift fabrics generate bulk permeability anisotropy that is to some extent predictable, with rift segments characterized by enhanced axial flow, and accommodation zones characterized by locally enhanced vertical permeability that is tectonically maintained. This provides a plausible explanation for the common occurrence of geothermal fields within accommodation zones and their notable absence within densely faulted rift segments. Maintenance of structural permeability in zones of active hydrothermal precipitation necessarily requires repeated brittle failure. Geothermal plumes therefore exploit tectonically maintained permeability within accommodation zones, with rift segments functioning mostly as drawdown regions. The influence of rift architecture on flow paths has important implications for geothermal extraction and epithermal mineral exploration within the TVZ and other structurally segmented hydrothermal systems, both active and extinct.