Evolution of Faulting Induced by Deep Fluid Injection, Paradox Valley, Colorado

2020 ◽  
Vol 110 (5) ◽  
pp. 2308-2327 ◽  
Author(s):  
Roger P. Denlinger ◽  
Daniel R. H. O’Connell

ABSTRACT High-pressure fluid injection into a subhorizontal confined aquifer at 4.3–4.6 km depth induced >7000 earthquakes between 1991 and 2012 within once seismically quiescent Paradox Valley in Colorado, with magnitudes up to Mw 3.9. Earthquake hypocenters expanded laterally away from the well with time, defining the margins of the aquifer pressurized by injection at the well. Within 5 km of the well, alignment of earthquake hypocenters defines strikes of nine vertical fault zones. Previous studies show that these fault zones predate injection, producing left-stepping offsets in the normal faults of the Wray-Mesa fault system that cradles Paradox Valley. Hypocenters, rakes, and strikes of 2041 well-constrained focal mechanisms show that most injection-related earthquakes occur where these vertical faults intersect the pressurized aquifer. Well-defined focal mechanisms show that this induced seismicity consists of Riedel shear faults at acute angles to the strikes of these fault zones. These small faults develop an anastomosing fault structure of focal planes along each planar fault zone, as fluid injection continues, even as their hypocenters define a single planar fault zone. Failure conditions at each hypocenter are found using a fully coupled poroelastic analysis of stress induced by fluid injection, and this analysis indicates a minimum Coulomb failure condition of 0.1 MPa. This failure condition is primarily a result of aquifer pore-fluid pressurization, as almost all well-located seismicity is within the pressurized aquifer. Reducing the rate of injection and frequent well shutdowns in the second decade nearly eliminated induced seismicity, except very near the well where gradients in pressurization are the largest. Despite these decreases in failure conditions and seismicity, some fault zones continued to produce earthquakes larger than M 3 as injection continued.

2021 ◽  
Author(s):  
Fabien Caroir ◽  
Frank Chanier ◽  
Virginie Gaullier ◽  
Julien Bailleul ◽  
Agnès Maillard-Lenoir ◽  
...  

<p>The Anatolia-Aegean microplate is currently extruding toward the South and the South-West. This extrusion is classically attributed to the southward retreat of the Aegean subduction zone together with the northward displacement of the Arabian plate. The displacement of Aegean-Anatolian block relative to Eurasia is accommodated by dextral motion along the North Anatolian Fault (NAF), with current slip rates of about 20 mm/yr. The NAF is propagating westward within the North Aegean domain where it gets separated into two main branches, one of them bordering the North Aegean Trough (NAT). This particular context is responsible for dextral and normal stress regimes between the Aegean plate and the Eurasian plate. South-West of the NAT, there is no identified major faults in the continuity of the NAF major branch and the plate boundary deformation is apparently distributed within a wide domain. This area is characterised by slip rates of 20 to 25 mm/yr relative to Eurasian plate but also by clockwise rotation of about 10° since ca 4 Myr. It constitutes a major extensional area involving three large rift basins: the Corinth Gulf, the Almiros Basin and the Sperchios-North Evia Gulf. The latter develops in the axis of the western termination of the NAT, and is therefore a key area to understand the present-day dynamics and the evolution of deformation within this diffuse plate boundary area.</p><p>Our study is mainly based on new structural data from field analysis and from very high resolution seismic reflexion profiles (Sparker 50-300 Joules) acquired during the WATER survey in July-August 2017 onboard the R/V “Téthys II”, but also on existing data on recent to active tectonics (i.e. earthquakes distribution, focal mechanisms, GPS data, etc.). The results from our new marine data emphasize the structural organisation and the evolution of the deformation within the North Evia region, SW of the NAT.</p><p>The combination of our structural analysis (offshore and onshore data) with available data on active/recent deformation led us to define several structural domains within the North Evia region, at the western termination of the North Anatolian Fault. The North Evia Gulf shows four main fault zones, among them the Central Basin Fault Zone (CBFZ) which is obliquely cross-cutting the rift basin and represents the continuity of the onshore Kamena Vourla - Arkitsa Fault System (KVAFS). Other major fault zones, such as the Aedipsos Politika Fault System (APFS) and the Melouna Fault Zone (MFZ) played an important role in the rift initiation but evolved recently with a left-lateral strike-slip motion. Moreover, our seismic dataset allowed to identify several faults in the Skopelos Basin including a large NW-dipping fault which affects the bathymetry and shows an important total vertical offset (>300m). Finally, we propose an update of the deformation pattern in the North Evia region including two lineaments with dextral motion that extend southwestward the North Anatolian Fault system into the Oreoi Channel and the Skopelos Basin. Moreover, the North Evia Gulf domain is dominated by active N-S extension and sinistral reactivation of former large normal faults.</p>


2021 ◽  
Author(s):  
Mohsen Goudarzi ◽  
René de Borst ◽  
Taras Gerya ◽  
Meng Li ◽  
van Dinther Ylona

<p>Accurate representation of fault zones is important in many applications in Earth sciences, including natural and induced seismicity. The framework developed here can efficiently model fault zone localization, evolution, and spontaneous fully dynamic earthquake sequences in a continuum plasticity framework. The geometrical features of the faults are incorporated into a regularized continuum framework, while the response of the fault zone is governed by a rate and state-dependent friction. Although a continuum plasticity model is advantageous to discrete approaches in representing evolving, unknown, or arbitrarily positioned faults, it is known that either non-associated plasticity or strain-softening can lead to mesh sensitivity of the numerical results in absence of an internal length scale. A common way to regularize the numerical model and introduce an internal length scale is by the adoption of a Kelvin-type visco-plasticity element. The visco-plastic rheological behavior for the bulk material is implemented along with a return-mapping algorithm for accurate stress and strain evolution. High slip rates (in the order of 1 m/s) are captured through numerical examples of a predefined strike-slip fault zone, where a detailed comparison with a reference discrete fault model is presented. Additionally, the regularization effect of the Kelvin viscosity parameter is studied on the fault slip velocity for a growing fault zone due to an initial material imperfection.  The model is consistently linearized leading to quadratic convergence of the Newton solver. Although the proposed framework is a step towards the modeling of earthquake sequences for induced seismicity applications, the numerical model is general and can be applied to all tectonic settings including subduction zones.</p><div> <div> <div> </div> <div> <div> <div> </div> <div> <p> </p> <p> </p> </div> </div> </div> </div> </div>


2020 ◽  
Author(s):  
Takahiro Tagami

<p>Thermal signatures as well as timing of fault motions can be constrained by thermochronological analyses of fault-zone rocks (e.g., Tagami, 2012, 2019).  Fault-zone materials suitable for such analyses are produced by tectocic and geochemical processes, such as (1) mechanical fragmentation of host rocks, grain-size reduction of fragments and recrystallization of grains to form mica and clay minerals, (2) secondary heating/melting of host rocks by frictional fault motions, and (3) mineral vein formation as a consequence of fluid advection associated with fault motions.  The geothermal structure of fault zones are primarily controlled by the following three factors: (a) regional geothermal structure around the fault zone that reflect background thermo-tectonic history of studied province, (b) frictional heating of wall rocks by fault motions and resultant heat transfer into surrounding rocks, and (c) thermal influences by hot fluid advection in and around the fault zone.  Geochronological/thermochronological methods widely applied in fault zones are K-Ar (<sup>40</sup>Ar/<sup>39</sup>Ar), fission-track (FT), and U-Th methods.  In addition, (U-Th)/He, OSL, TL and ESR methods are applied in some fault zones, in order to extract temporal information related to low temperature and/or recent fault activities.  Here I briefly review the thermal sensitivity of individual thermochronological systems, which basically controls the response of each method against faulting processes.  Then, the thermal sensitivity of FTs is highlighted, with a particular focus on the thermal processes characteristic to fault zones, i.e., flash and hydrothermal heating.  On these basis, representative examples as well as key issues, including sampling strategy, are presented to make thermochronological analysis of fault-zone materials, such as fault gouges, pseudotachylytes and mylonites, along with geological, geomorphological and seismological implications.  Finally, the thermochronological analyses of the Nojima fault are overviewed, as an example of multidisciplinary investigations of an active seismogenic fault system.</p><p> </p><p>References:</p><ol><li>Tagami, 2012. Thermochronological investigation of fault zones. Tectonophys., 538-540, 67-85, doi:10.1016/j.tecto.2012.01.032.</li> <li>Tagami, 2019. Application of fission track thermochronology to analyze fault zone activity. Eds. M. G. Malusa, P. G. Fitzgerald, Fission track thermochronology and its application to geology, 393pp, 221-233, doi: 10.1007/978-3-319-89421-8_12.</li> </ol>


2020 ◽  
Author(s):  
Javiera Ruz ◽  
Muriel Gerbault ◽  
José Cembrano ◽  
Pablo Iturrieta ◽  
Camila Novoa Lizama ◽  
...  

<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>


2020 ◽  
Author(s):  
Inbar Vaknin ◽  
Andy Nicol ◽  
Conrad Childs

<p>Fault surfaces and fault zones have been shown to have complex geometries comprising a range of morphologies including, segmentation, tip-line splays and slip-surface corrugations (e.g., Childs et al., 2009*). The three-dimensional (3D) geometries of faults (and fault zones) is difficult to determine from outcrop data which are typically 2D and limited in size. In this poster we examine the small-scale geometries of faults from normal faults cropping out in well bedded parts of the Mount Messenger and Mahakatino formations in Taranaki, New Zealand. We present two main datasets; i) measurements and maps of 2D vertical and horizontal sections for in excess of 200 faults and, ii) 3D fault model of a small-fault (vertical displacement ~1 cm) produced by serial fault-perpendicular sections of a block 10x10x13 cm. The sectioned block contains a single fault that offsets sand and silt layers, and comprises two main dilational bends; in the 3D model we map displacement, bedding and fault geometry for the sectioned fault zone. Faults in the 2D dataset comprise a range of geometries including, vertical segmentation, bends, splays and fault-surface corrugations. Although we have little information on the local magnitudes and orientations of stresses during faulting, geometric analysis of the fault zones provides information on the relationships between bed characteristics (e.g., thickness, induration and composition) and fault-surface orientations. The available data supports the view that the strike and dip of fault surfaces vary by up to 25° producing undulations or corrugations on fault surfaces over a range of scales from millimetres to metres and in both horizontal and vertical directions. Preliminary analysis of the available data suggests that these corrugations appear to reflect fault refractions due to changing bed lithologies (unexpectedly the steepest sections of faults are in mudstone beds), breaching of relays and development of conjugate fault sets. The relative importance of these factors and their importance for fault geometry will be explored further in the poster.</p><p> </p><p>*Childs, C., Walsh, J.J., Manzocchi, T., Bonson, C., Nicol A., Schöpfer, M.P.J. 2009. A geometric model of fault zone and fault rock thickness variations. Journal of Structural Geology 31, 117-127.</p>


2011 ◽  
Vol 62 (4) ◽  
pp. 345-359 ◽  
Author(s):  
Erman Özsayin ◽  
Kadir Dirik

The role of oroclinal bending in the structural evolution of the Central Anatolian Plateau: evidence of a regional changeover from shortening to extensionThe NW-SE striking extensional Inönü-Eskişehir Fault System is one of the most important active shear zones in Central Anatolia. This shear zone is comprised of semi-independent fault segments that constitute an integral array of crustal-scale faults that transverse the interior of the Anatolian plateau region. The WNW striking Eskişehir Fault Zone constitutes the western to central part of the system. Toward the southeast, this system splays into three fault zones. The NW striking Ilıca Fault Zone defines the northern branch of this splay. The middle and southern branches are the Yeniceoba and Cihanbeyli Fault Zones, which also constitute the western boundary of the tectonically active extensional Tuzgölü Basin. The Sultanhanı Fault Zone is the southeastern part of the system and also controls the southewestern margin of the Tuzgölü Basin. Structural observations and kinematic analysis of mesoscale faults in the Yeniceoba and Cihanbeyli Fault Zones clearly indicate a two-stage deformation history and kinematic changeover from contraction to extension. N-S compression was responsible for the development of the dextral Yeniceoba Fault Zone. Activity along this structure was superseded by normal faulting driven by NNE-SSW oriented tension that was accompanied by the reactivation of the Yeniceoba Fault Zone and the formation of the Cihanbeyli Fault Zone. The branching of the Inönü-Eskişehir Fault System into three fault zones (aligned with the apex of the Isparta Angle) and the formation of graben and halfgraben in the southeastern part of this system suggest ongoing asymmetric extension in the Anatolian Plateau. This extension is compatible with a clockwise rotation of the area, which may be associated with the eastern sector of the Isparta Angle, an oroclinal structure in the western central part of the plateau. As the initiation of extension in the central to southeastern part of the Inönü-Eskişehir Fault System has similarities with structures associated with the Isparta Angle, there may be a possible relationship between the active deformation and bending of the orocline and adjacent areas.


2021 ◽  
Author(s):  
Miriana Chinello ◽  
Michele Fondriest ◽  
Giulio Di Toro

<p>The Italian Central Apennines are one of the most seismically active areas in the Mediterranean (e.g., L’Aquila 2009, Mw 6.3 earthquake). The mainshocks and the aftershocks of these earthquake sequences propagate and often nucleate in fault zones cutting km-thick limestones and dolostones formations. An impressive feature of these faults is the presence, at their footwall, of few meters to hundreds of meters thick damage zones. However, the mechanism of formation of these damage zones and their role during (1) individual seismic ruptures (e.g., rupture arrest), (2) seismic sequences (e.g., aftershock evolution) and (3) seismic cycle (e.g., long term fault zone healing) are unknown. This limitation is also due to the lack of knowledge regarding the distribution, along strike and with depth, of damage with wall rock lithology, geometrical characteristics (fault length, inherited structures, etc.) and kinematic properties (cumulative displacement, strain rate, etc.) of the associated main faults.</p><p>Previous high-resolution field structural surveys were performed on the Vado di Corno Fault Zone, a segment of the ca. 20 km long Campo Imperatore normal fault system, which accommodated ~ 1500 m of vertical displacement (Fondriest et al., 2020). The damage zone was up to 400 m thick and dominated by intensely fractured (1-2 cm spaced joints) dolomitized limestones with the thickest volumes at fault oversteps and where the fault cuts through an older thrust zone. Here we describe two minor faults located in the same area (Central Apennines), but with shorter length along strike. They both strike NNW-SSE and accommodated a vertical displacement of ~300 m.</p><p>The Subequana Valley Fault is about 9 km long and consists of multiple segments disposed in an en-echelon array. The fault juxtaposes pelagic limestones at the footwall and quaternary deposits at the hanging wall. The damage zone is < 25 m  thick  and comprises fractured (1-2 cm spaced joints) limestones beds with decreasing fracture intensity moving away from the master fault. However, the damage zone thickness increases up to ∼100 m in proximity of subsidiary faults striking NNE-SSW. The latter could be reactivated inherited structures.</p><p>The Monte Capo di Serre Fault is about 8 km long and characterized by a sharp ultra-polished master fault surface which cuts locally dolomitized Jurassic platform limestones. The damage zone is up to 120 m thick and cut by 10-20 cm spaced joints, but it reaches an higher fracture intensity where is cut by subsidiary, possibly inherited, faults striking NNE-SSW.</p><p>Based on these preliminary observations, faults with similar displacement show comparable damage zone thicknesses. The most relevant damage zone thickness variations are related to geometrical complexities rather than changes in lithology (platform vs pelagic carbonates).  In particular, the largest values of damage zone thickness and fracture intensity occur at fault overstep or are associated to inherited structures. The latter, by acting as strong or weak barriers (sensu Das and Aki, 1977) during the propagation of seismic ruptures, have a key role in the formation of damage zones and the growth of normal faults.</p>


2020 ◽  
Author(s):  
Hagen Deckert ◽  
Steffen Abe ◽  
Wolfgang Bauer

<p>In the course of hydrocarbon or geothermal exploration the characterisation of fault zone architectures is of interest for fluid flow modelling and geomechanical studies. Seismic data normally offer the best information for the identification of fault zone geometries in sedimentary basins. However, the internal structure or the damage zone of a fault can be hardly resolved with seismic data as displacements along single fault strands or fractures are by far too small. Thus, it is not possible to directly map small scale faults with seismic methods, though these structures might significantly influence fluid flow. We try to examine the architecture of extensional fault zones in carbonate rocks at subseismic scales by using discrete element method (DEM) techniques to numerically simulate the evolution of fault zones including their associated damage zones.</p><p>As a case study we have analysed the geometry, displacement and fault width of normal faults in fine grained jurassic limestones in a quarry in Franconia, Germany. The quarry shows a rather simple set of conjugated 60deg dipping normal faults. Displacement is rather small and varies between c. 5cm up to c. 2m, some faults show almost no offset. The fault thickness varies between 2cm and c. 1m. A closer investigation of the fault geometries reveals, next to planar parts, sometimes complex fault zone structures including restraining and releasing bends, multiple fault strands as well as lenses and associated riedel shears. Analysis of high resolution photogrammetric data revealed a high number of small scale fractures between neighbouring discrete fault surfaces which are interpreted as highly fractured damage zones. Some faults with rather small displacement suggest that the overall inclination of the fault is a result of small subvertical sections which are connected in a staircase like appearance. </p><p>The DEM models simulate normal faulting in a layered marl-limestone sequence driven by the displacement of an underlying basement fault. Different layer geometries and effective vertical stresses in the range of 15-45 MPa, equivalent to an overburden thickness of c. 1000-3000m, have been used in the models. The stress range covers the maximum burial depth of the carbonates, which is assumed to be c. 1500m. Material properties used in the DEM were calibrated based on laboratory data, i.e. results of triaxial deformation tests on the studied limestones.</p><p>Results of the models show fault geometries which resemble those observed in the studied outcrop. In particularly under low stress, small offsets and with strongly decoupled layers we observe steeply dipping faults (>70deg) which also show staircase structures composed of sub-vertical fractures within each of the layers and horizontal offsets along the layer interfaces. We also observe the development of multiple fault strands and associated damage zones. </p><p>Our study shows that the DEM models are capable to reproduce observed fault geometries and damage zones. The results help to understand fault zone architectures and depict highly fractured areas in a sub-seismic scale.</p>


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-15
Author(s):  
Ziye Lu ◽  
Yingtao Li ◽  
Ning Ye ◽  
Shaonan Zhang ◽  
Chaojin Lu ◽  
...  

The exploration of deeply buried hydrocarbon is still a challenge for the petroleum geology. The Shunbei area is a newly discovered oil fields, located in the center of the Tarim Basin. The oil is mainly yielded from the Middle–Lower Ordovician carbonate reservoirs with depth > 7000  m in the Shunbei No. 1 and No. 5 fault zones. Calcite cements filled in vugs (v-calcite) and fractures (f-calcite) are identified in limestones and dolostones of the carbonate reservoirs. F-calcites in the Shunbei No. 1 fault zone trap secondary inclusions in trails, which comprise liquid-dominated biphase aqueous inclusions, liquid-dominated biphase oil inclusions, and/or oil-bearing triphase inclusions. F-calcite and v-calcite in the No. 5 fault zone trap secondary inclusions in trails, which consist of liquid-only monophase aqueous inclusions, liquid-dominated biphase aqueous inclusions, liquid-dominated biphase oil inclusions, liquid-only monophase oil inclusions, and/or oil-bearing triphase inclusions. The ranges of the homogenization temperature ( T h ) and ice-melting temperature ( T m − ice ) in the Shunbei No. 1 fault zone are, respectively, 130–150°C and -2.1–-1.5°C. The coexistence of liquid-only and liquid-dominated aqueous inclusions in the Shunbei No. 5 fault zone indicates that the aqueous inclusions are trapped at low temperatures. The aqueous inclusions in the Shunbei No. 5 fault zone show a range from -0.4 to -0.2°C in T m − ice which is very close to the meteoric fluid. In the context of the burial-thermal history and the Cambrian source rock evolution, the charging process of hydrocarbon in the Shunbei No. 1 and No. 5 fault zones corresponds to the Silurian and Middle Ordovician, respectively. Results of fluid inclusions indicate a tightly coupling relationship between the hydrocarbon charging process and fault system evolution in the Shunbei area. This study reveals the application of fluid inclusion under the systemically petrographic constraints to decipher the charging history of hydrocarbon, especially for the deeply buried reservoirs.


Geophysics ◽  
1956 ◽  
Vol 21 (2) ◽  
pp. 368-387 ◽  
Author(s):  
Hans P. Laubscher

Seismic reflections in normal fault zones in the Eastern Venezuelan basin usually appear distorted. Studies of reflections over fault structures delineated by drilling indicate that this is due to the similar effects of two entirely different phenomena: 1. True structural deformation of beds on the downthrown side. 2. Purely seismic distortion of reflections from underneath the fault. Analysis indicates that the structural deformations form an integral part of the fault zone; the purely seismic distortion is caused by passage of the wave through this zone of deformation.


Sign in / Sign up

Export Citation Format

Share Document