Permeability contrasts of fault zones - from conceptual model to numerical simulation

2020 ◽  
Author(s):  
Ulrich Kelka ◽  
Thomas Poulet ◽  
Luk Peeters
Keyword(s):  
Fluid Flow ◽  
Fault Zone ◽  
Large Scale ◽  
Regional Scale ◽  
Damage Zone ◽  
Fluid Pressure ◽  
Fault Zones ◽  
Strong Impact ◽  
Fault Zone Architecture ◽  
The Impact

<p>Fault and fracture networks can govern fluid flow patterns in the subsurface and predicting fluid flow on a regional scale is of interest in a variety of fields like groundwater management, mining engineering, energy, and mineral resources. Especially the pore fluid pressure can have a strong impact on the strength of fault zones and might be one of the drivers for fault reactivation. Reliable simulations of the transient changes in fluid pressure need to account for the generic architecture of fault zones that comprises strong permeability contrast between the fault core and damage zone.</p><p>Particularly, the distribution and connectivity of large-scale fault zones can have a strong impact on the flow field. Yet, modelling numerically such features in their full complexity remains challenging. Often faults zones are conceptualized as forming exclusively either barriers or conduits to fluid flow. However, a generic architecture of fault zones often comprises a discrete fault core surrounded by a diffuse damage zone and conceptualizing large scale discontinuities simply as a barrier or conduit is unlikely to capture the regional scale fluid flow dynamics. It is known that if the fault zone is hosted in low-permeability strata, such as clays or crystalline rocks, a transversal flow barrier can form along the fault core whereas the fracture-rich fault damage zone represents a longitudinal conduit. In more permeable host-rocks (i.e. sandstones or carbonates) the reverse situation can occur, and the permeability distributions in the damage zones can be governed by the abundance of low-permeability deformation features. A reliable numerical model needs to account for the difference and strong contrasts in fluid flow properties of the core and the damage zone, both transversally and longitudinally, in order to make prediction about the regional fluid flow pattern.</p><p>Here, we present a numerical method that accounts for the generic fault zone architecture as lower dimensional interfaces in conforming meshes during fluid flow simulations in fault networks. With this method we aim to decipher the impact of fault zone architecture on subsurface flow pattern and fluid pressure evolution in fractures and faulted porous media. The method is implemented in a finite element framework for Multiphysics simulations. We demonstrate the impact of considering the more generic geological structure of individual faults on the flow field by conceptualizing discontinuities either as barriers, conduits or as a conduit-barrier system and show were these conceptualizations are applicable in natural systems. We further show that a reliable regional scale fluid flow simulation in faulted porous media needs to account for the generic fault zone architecture. The approach is finally used to evaluate the fluid flow response of statistically parameterised faulted media, in order to investigate the impact and sensitivity of each variable parameter.</p>

Partitioned Permeability Diagram: an innovative way to estimate Fault Zones hydraulics.

2021 ◽  
Author(s):  
Irène Aubert ◽  
Juliette Lamarche ◽  
Philippe Leonide
Keyword(s):  
Fault Zone ◽  
Damage Zone ◽  
Fault Zones ◽  
Matrix Permeability ◽  
Vertical Fault ◽  
Core Permeability ◽  
The Impact ◽  
Evolution With Time ◽  
Fluid Pathway

<p>Understanding the impact of fault zones on reservoir trap properties is a major challenge for a variety of geological ressources applications. Fault zones in cohesive rocks are complex structures, composed of 3 components: rock matrix, damage zone fractures and fault core rock. Despite the diversity of existing methods to estimate fault zone permeability/drain properties, up to date none of them integrate simultaneously the 3 components of fracture, fault core and matrix permeability, neither their evolution with time. We present a ternary plot that characterizes the fault zones permeability as well as their drainage properties. The ternary plot aims at (i) characterizing the fault zone permeability between the three vertices of matrix, fractures and fault core permeability ; and at (ii) defining the drain properties among 4 possible hydraulic system: (I) good horizontal and vertical, fault-perpendicular and -parallel; (II) moderate parallel fluid pathway; (III) good parallel fault-core and (IV) good parallel fractures. The ternary plot method is valid for 3 and 2 components fault zones. The application to the Castellas Fault case study show the simplicity and efficiency of the plot for studying underground and/or fossil, simple or polyphase faults in reservoirs with complete or limited permeability data.</p>

Faults and magma reservoirs along the Southern Andes Volcanic zone (SAVZ): linking oservations and numerical models of stress change controlling magmatic and hydrothermal fluid flow

2020 ◽  
Author(s):  
Javiera Ruz ◽  
Muriel Gerbault ◽  
José Cembrano ◽  
Pablo Iturrieta ◽  
Camila Novoa Lizama ◽  
...  
Keyword(s):  
Fluid Flow ◽  
Stress Field ◽  
Fault Zone ◽  
Numerical Models ◽  
Fluid Pressure ◽  
Seismic Inversion ◽  
Fault Zones ◽  
Fault System ◽  
Far Field ◽  
Trigger Failure

<p> The Chilean margin is amongst the most active seismic and volcanic areas on Earth. It hosts active and fossil geothermal and mineralized systems of economic interest documenting significant geofluid migration through the crust. By comparing numerical models with field and geophysical data, we aim at pinning when and where fluid migration occurs through porous domains, fault zone conduits, or remains stored at depth awaiting a more appropriate stress field. <span>Dyking and volcanic activity occur within fault zones</span> <span>along the S</span><span>A</span><span>VZ, linked with stress field variations</span> <span>in spatial and temporal association with</span> –<span>short therm-</span> <span>seismicity</span> <span>and -long term- oblique </span><span>plate </span><span>convergence.</span> <span>Volcanoes and geothermal domains are mostly located along or at the intersection of margin-oblique fault zones (Andean Transverse Faults), and along margin-parallel strike slip zones, some which may cut the entire lithosphere (Liquiñe-Ofqui fault system). Wh</span><span>ereas</span><span> the big picture displays</span> <span>fluid flow straight to the surface, at close look significant offsets between crustal structures occur. 3D numerical models using conventional elasto-plastic rheology provide insights on the interaction of (i) an inflating magmatic cavity, (ii) a slipping fault zone, and (iii) regional tectonic stresses. Applying either (i) a magmatic overpressure or (ii) a given fault slip can trigger failure of the intervening rock, and generate either i) fault motion or ii) magmatic reservoir failure, respectively, but only for distances less than the structures' breadth even at low rock</span> <span>strength. However, at greater inter-distances the bedrock domain in between the fault zone and the magmatic cavity undergoes dilatational strain of the order of 1-5x10-5. This dilation opens the bedrock’s pore space and forms «pocket domains» that may store up-flowing over-pressurized fluids, which may then further chemically</span> interact<span> with the bedrock, for the length of time</span> <span>that</span> <span>these pockets remain open. These porous pockets</span> <span>can reach kilometric size, questioning their parental link with outcropping plutons along the margin. Moreover, bedrock permeability may also increase as fluid flow diminishes effective bedrock friction and cohesion. Comparison with rock experiments indicates that such stress and fluid pressure changes may eventually trigger failure at the intermediate timescale (repeated slip or repeated inflation). Finally, incorporating far field compression (iii)</span> <span>loads the bedrock to</span> <span>a state of stress at the verge of failure. Then, failure around the magmatic </span><span>reservoir</span><span> or </span><span>at</span> <span>the fault zone occurs for lower load</span><span>ing</span><span>.</span> <span>Permanent tectonic loading on the one hand, far field episodic seismic inversion of the stress field on the other, and localized failure all together promote a transient stress field, thus explaining the occurrence of transient fluid pathways on seemingly independent timescales. These synthetic models are then discussed with regards to specific cases along the SVZ, particularly the Tatara-San Pedro area (~36°S), where magnetotelluric profiles </span><span>document</span><span> conductive volumes at different depths underneath active faults, volcanic edifices and geothermal vents. We discuss the mechanical link between these deep sources and surface structures</span>.</p>

Discrete and continuous conceptual models of fault zones.

2021 ◽  
Author(s):  
Tom Manzocchi
Keyword(s):  
Numerical Model ◽  
Conceptual Model ◽  
Fault Zone ◽  
Large Scale ◽  
Damage Zone ◽  
Fault Zones ◽  
Zone Model ◽  
Multiple Model ◽  
Modelling Studies ◽  
Scale Properties

<p>Faults can control the large-scale properties of rock volumes through their behaviour as flow conduits and/or barriers or by localising geomechanical effects. Hence, often the fidelity of a numerical model of faulted site relies on the accuracy with which the fault zone is represented.  There are two distinct factors that must be considered in a modelling study: first, does the model contain the most relevant characteristics of the fault that influence the behaviour of interest; and second, are these characteristics assigned realistic and representative values that capture both their natural variability and the uncertainty with which they can be determined for the specific case of interest. These two factors are contained in the conceptual fault model and choice of modelling proxy-properties, respectively.</p><p>In recent years, two classes of conceptual fault zone model have dominated the description of fault zones, broadly characterised by either a continuous or a discrete approach. Continuous fault zone properties (e.g. fault core and damage zone thickness, displacement partitioning statistics) often show high variability which many modelling studies attempt to capture by running multiple model containing property values sampled from the distribution. Discrete descriptions focus on the presence of individual fault zone elements (e.g. shale smears, relay zones), and models guided by a discrete conceptual model attempt to place representative frequencies of elements. A single discrete model might contain the same property distributions as an ensemble of continuous models yet, because it contains a representative frequency of different elements, its behaviour might lie beyond the extreme behaviour of the continuous ensemble. Hence, the manner in which a geologist’s conceptual model is represented in a modeller’s numerical model can be hugely important for the outcome of the study, and it is in the interest of both modellers and geologists to ensure that they have a correct understanding of the other’s part of the process.</p>

scholarly journals Heterogeneous stresses and deformation mechanisms at shallow crustal conditions, Hungaroa Fault Zone, Hikurangi Subduction Margin, New Zealand

2020 ◽  
Author(s):  
Carolyn Boulton ◽  
Marcel Mizera ◽  
Maartje Hamers ◽  
Inigo Müller ◽  
Martin Ziegler ◽  
...  
Keyword(s):  
Friction Coefficient ◽  
Fault Zone ◽  
Mixed Mode ◽  
Fluid Pressure ◽  
Fault Zones ◽  
Deformation Mechanisms ◽  
Multiple Observations ◽  
Calcite Veins ◽  
Dynamic Recrystallisation ◽  
Clumped Isotope

<p>The Hungaroa Fault Zone (HFZ), an inactive thrust fault along the Hikurangi Subduction Margin, accommodated large displacements (~4–10 km) at the onset of subduction in the early Miocene. Within a 40 m-wide high-strain fault core, calcareous mudstones and marls display evidence for mixed-mode viscous flow and brittle fracture, including: discrete faults; extensional veins containing stretched calcite fibers; shear veins with calcite slickenfibers; calcite foliation-boudinage structures; calcite pressure fringes; dark dissolution seams; stylolites; embayed calcite grains; and an anastomosing phyllosilicate foliation.</p><p>Multiple observations indicate a heterogeneous stress state within the fault core. Detailed optical and electron backscatter diffraction-based texture analysis of syntectonic calcite veins and isoclinally folded limestone layers within the fault core reveal that calcite grains have experienced intracrystalline plasticity and interface mobility, and local subgrain development and dynamic recrystallisation. The recrystallized grain size in two calcite veins of 6.0±3.9 µm (n=1339; 1SD; HFZ-H4-5.2m_A;) and 7.2±4.2µm (n=406; 1SD; HFZ-H4-19.9m) indicate high differential stresses (~76–134 MPa). Hydrothermal friction experiments on a foliated, calcareous mudstone yield a friction coefficient of μ≈0.35. Using this friction coefficient in the Mohr-Coulomb failure criterion yields a maximum differential stress of 55 MPa at 4 km depth, assuming a minimum principal stress equal to the vertical stress, an average sediment density of 2350 kg/m<sup>3</sup>, and hydrostatic pore fluid pressure. Interestingly, calcareous microfossils within the foliated mudstone matrix are undeformed. Moreover, calcite veins are oriented both parallel to and highly oblique to the foliation, indicating spatial and/or temporal variations in the maximum principle stress azimuth.</p><p>To further constrain HFZ deformation conditions, clumped isotope geothermometry was performed on six syntectonic calcite veins, yielding formation temperatures of 79.3±19.9°C (95% confidence interval). These temperatures are well below those at which dynamic recrystallisation of calcite is anticipated and exclude shear heating and the migration of hotter fluids as an explanation for dynamic recrystallisation of calcite at shallow crustal levels (<5 km depth).</p><p>Our results indicate that: (1) stresses are spatiotemporally heterogeneous in crustal fault zones containing mixtures of competent and incompetent minerals; (2) heterogeneous deformation mechanisms, including frictional sliding, pressure solution, dynamic recrystallization, and mixed-mode fracturing accommodate slip in shallow crustal fault zones; and (3) brittle fractures play a pivotal role in fault zone deformation by providing fluid pathways that promote fluid-enhanced recovery and dynamic recrystallisation in the deforming calcite at remarkably low temperatures. Together, field geology, microscopy, and clumped isotope geothermometry provide a powerful method for constraining the multiscale slip behavior of large-displacement fault zones.</p>

The importance of phase morphology for rheology of ferropericlase-bridgmanite mixtures

2021 ◽  
Author(s):  
Marcel Thielmann ◽  
Gregor Golabek ◽  
Hauke Marquardt
Keyword(s):  
Lower Mantle ◽  
Mantle Convection ◽  
Effective Viscosity ◽  
Large Scale ◽  
Numerical Models ◽  
Phase Morphology ◽  
Strong Impact ◽  
Analytical Approximations ◽  
The Impact ◽  
Strong Control

<p>The rheology of the Earth’s lower mantle is poorly constrained due to a lack of knowledge of the rheological behaviour of its constituent minerals. In addition, the lower mantle does not consist of only a single, but of multiple mineral phases with differing deformation behaviour. The rheology of Earth’s lower mantle is thus not only controlled by the rheology of its individual constituents (bridgmanite and ferropericlase), but also by their interplay during deformation. This is particularly important when the viscosity contrast between the different minerals is large. Experimental studies have shown that ferropericlase may be significantly weaker than bridgmanite and may thus exert a strong control on lower mantle rheology.</p><p>Here, we thus explore the impact of phase morphology on the rheology of a ferropericlase-bridgmanite mixture using numerical models. We find that elongated ferropericlase structures within the bridgmanite matrix significantly lower the effective viscosity, even in cases where no interconnected network of weak ferropericlase layers has been formed. In addition to the weakening, elongated ferropericlase layers result in a strong viscous anisotropy. Both of these effects may have a strong impact on lower mantle dynamics, which makes is necessary to develop upscaling methods to include them in large-scale mantle convection models. We develop a numerical-statistial approach to link the statistical properties of a ferropericlase-bridgmanite mixture to its effective viscosity tensor. With this approach, both effects are captured by analytical approximations that have been derived to describe the evolution of the effective viscosity (and its anisotropy) of a two-phase medium with aligned elliptical inclusions, thus allowing to include these microscale processes in large-scale mantle convection models.</p>

Vulnerability of coastal areas to increased aquifer saturation due to climate change

2021 ◽  
Author(s):  
Alexandre Gauvain ◽  
Ronan Abhervé ◽  
Jean-Raynald de Dreuzy ◽  
Luc Aquilina ◽  
Frédéric Gresselin
Keyword(s):  
Climate Change ◽  
Sea Level ◽  
Large Scale ◽  
Regional Scale ◽  
Drainage System ◽  
Hydrological Regime ◽  
Coastal Areas ◽  
Western Coast ◽  
Study Sites ◽  
The Impact

<p>Like in other relatively flat coastal areas, flooding by aquifer overflow is a recurring problem on the western coast of Normandy (France). Threats are expected to be enhanced by the rise of the sea level and to have critical consequences on the future development and management of the territory. The delineation of the increased saturation areas is a required step to assess the impact of climate change locally. Preliminary models showed that vulnerability does not result only from the sea side but also from the continental side through the modifications of the hydrological regime.</p><p>We investigate the processes controlling these coastal flooding phenomena by using hydrogeological models calibrated at large scale with an innovative method reproducing the hydrographic network. Reference study sites selected for their proven sensitivity to flooding have been used to validate the methodology and determine the influence of the different geomorphological configurations frequently encountered along the coastal line.</p><p>Hydrogeological models show that the rise of the sea level induces an irregular increase in coastal aquifer saturations extending up to several kilometers inland. Back-littoral channels traditionally used as a large-scale drainage system against high tides limits the propagation of aquifer saturation upstream, provided that channels are not dominantly under maritime influence. High seepage fed by increased recharge occurring in climatic extremes may extend the vulnerable areas and further limit the effectiveness of the drainage system. Local configurations are investigated to categorize the influence of the local geological and geomorphological structures and upscale it at the regional scale.</p>

Iron Oxide (U–Th)/He Thermochronology: New Perspectives on Faults, Fluids, and Heat

Elements ◽  
2020 ◽  
Vol 16 (5) ◽  
pp. 319-324
Author(s):  
Emily H. G. Cooperdock ◽  
Alexis K. Ault
Keyword(s):  
Fluid Flow ◽  
Iron Oxide ◽  
Heat Exchange ◽  
Iron Oxides ◽  
Fault Zone ◽  
Frictional Heating ◽  
Fault Zones ◽  
Chemical Behavior ◽  
Rock Interaction ◽  
Dynamic Motion

Fault zones record the dynamic motion of Earth’s crust and are sites of heat exchange, fluid–rock interaction, and mineralization. Episodic or long-lived fluid flow, frictional heating, and/or deformation can induce open-system chemical behavior and make dating fault zone processes challenging. Iron oxides are common in a variety of geologic settings, including faults and fractures, and can grow at surface-to magmatic temperatures. Recently, iron oxide (U–Th)/He thermochronology, coupled with microtextural and trace element analyses, has enabled new avenues of research into the timing and nature of fluid–rock interactions and deformation. These constraints are important for understanding fault zone evolution in space and time.

Structural evolution of a major fault zone: the Cape Ray Fault Zone, southwestern Newfoundland, Canada

1996 ◽  
Vol 33 (2) ◽  
pp. 199-215 ◽  
Author(s):  
Benoît Dubé ◽  
Kathleen Lauzière
Keyword(s):  
Fault Zone ◽  
Large Scale ◽  
Regional Scale ◽  
Point Group ◽  
Structural Evolution ◽  
Strain Partitioning ◽  
Two Phase ◽  
Detailed Structural Analysis ◽  
Grand Bay ◽  
Geodynamic Processes

The Cape Ray Fault Zone is a major Paleozoic structure in southwestern Newfoundland, and occurs at or close to the boundary between two major continental blocks, Laurentia and Avalonia. A detailed structural analysis demonstrates that the fault records early reverse-sinistral thrusting of the Grand Bay Complex at amphibolite grade (D2), followed by a protracted event (D3) characterized by reverse-dextral thrusting of the Grand Bay Complex rocks on top of the supracrustal rocks of the Windsor Point Group and retrogression to greenschist facies, as well as a pre-384 Ma orogen-parallel dextral transcurrent mylonite (D4) during the later stages of the collision. Regional-scale strain partitioning induced heterogeneity of strain both along and across the strike of the Cape Ray Fault Zone. The east–west-oriented segment of the Cape Ray Fault Zone is a tear fault that accommodated differential displacement along the length of the fault. Later stages of the deformation include post-384 Ma sinistral transcurrent reactivation of the dextral mylonite and extension. The reverse-sinistral thrusting and the reverse-dextral motion occurred between 415 and 386 Ma and correspond to the two-phase Acadian orogeny recognized at the scale of the orogen and believed to be related to collision between Laurentia and Avalonia. The Cape Ray Fault Zone preserves evidence of large-scale geodynamic processes affecting rocks where the kinematics and the timing are well constrained.

Sign in / Sign up

Export Citation Format

Share Document