Refine
Has Fulltext
- no (24)
Document Type
- Article (24) (remove)
Language
- English (24)
Is part of the Bibliography
- yes (24)
Keywords
- induced seismicity (5)
- fault creep (2)
- fluid injection (2)
- fluid-induced seismicity (2)
- geothermal (2)
- maximum magnitude (2)
- pore pressure (2)
- reservoir characterization (2)
- source parameters (2)
- strain localization (2)
- torsion (2)
- (compliant) cracks (1)
- 2-D numerical model (1)
- 3D geomechanical numerical model (1)
- Alpine Fault (1)
- Anisotropy (1)
- Authigenic mineral formation (1)
- Brittle deformation (1)
- Brittle precursors (1)
- Broad-band seismometers (1)
- Bulk compressibility (1)
- Clay rock (1)
- Contact model (1)
- Coulomb failure stress (1)
- Dynamic bulk modulus (1)
- Earthquakes (1)
- Equant (stiff) pores (1)
- Fault healing (1)
- Fault-rock microstructures (1)
- Fluid Flow (1)
- Fluid-rock interaction (1)
- Fracture and flow (1)
- Geomechanics (1)
- Granular materials (1)
- High-temperature shear zones (1)
- Induced seismicity (1)
- Ion-Beam (FIB) (1)
- Micromechanical model (1)
- Microstructural deformation mechanisms (1)
- Narrow (1)
- Porous sandstone (1)
- Pressure-temperature and strain rate-dependent mechanical behaviour (1)
- Roughness (1)
- Sandy facies of Opalinus Clay (1)
- Sea of Marmara (1)
- Shear zones (1)
- Static bulk modulus (1)
- Strain localization Carrara marble (1)
- The Geysers (1)
- Triaxial deformation (1)
- Triaxial deformation experiments (1)
- acoustic emission (1)
- aseismic (1)
- central south island (1)
- deep-level mining (1)
- dislocation creep (1)
- earthquake nucleation (1)
- earthquake sequences (1)
- fault slip (1)
- fluid pressurization (1)
- focal mechanisms (1)
- fracking (1)
- frictional properties (1)
- hanging wall (1)
- internal structure (1)
- loading conditions (1)
- localization (1)
- marble (1)
- non-double-couple components (1)
- passive seismic monitoring (1)
- picoseismicity (1)
- rate (1)
- rheological weakening (1)
- rupture (1)
- san andreas fault (1)
- seismic (1)
- seismic hazard (1)
- seismic moment release (1)
- seismomechanics (1)
- slow slip (1)
- static stress change (1)
- stick slip (1)
- stick-slip (1)
- stick–slip tests (1)
- strain transient (1)
- strainmeter data (1)
- stress state (1)
- tabular mining (1)
- tensile opening (1)
- thermal effect (1)
- thickness-displacement relationships (1)
- transform faults (1)
- two phase aggregates (1)
- weakening mechanisms (1)
We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.
Preparatory mechanisms accompanying or leading to nucleation of larger earthquakes have been observed at both laboratory and field scales, but conditions favoring the occurrence of observable preparatory processes are still largely unknown. In particular, it remains a matter of debate why some earthquakes occur spontaneously without noticeable precursors as opposed to events that are preceded by an extended failure process. In this study, we have generated new high-resolution seismicity catalogs framing the occurrence of 20 M-L > 2.5 earthquakes at The Geysers geothermal field in California. To this end, a seismicity catalog of the 11 days framing each large event was created. We selected 20 sequences sampling different hypocentral depths and hydraulic conditions within the field. Seismic activity and magnitude frequency distributions displayed by the different earthquake sequences are correlated with their location within the reservoir. Sequences located in the northwestern part of the reservoir show overall increased seismic activity and low b values, while the southeastern part is dominated by decreased seismic activity and higher b values. Periods of high injection coincide with high b values and vice versa. These observations potentially reflect varying differential and mean stresses and damage of the reservoir rocks across the field. About 50% of analyzed sequences exhibit no change in seismicity rate in response to the large main event. However, we find complex waveforms at the onset of the main earthquake, suggesting that small ruptures spontaneously grow into or trigger larger events.
Triaxial high temperature (900 °C) deformation experiments were conducted at constant strain rate in a Paterson-type deformation apparatus on cylinders of Carrara marble with two right or left stepping, non-overlapping weak inclusions of Solnhofen limestone, oriented at 45° to the cylinders’ longitudinal axes. Applying different values of confinement (30, 50, 100 and 300 MPa) we induced various amounts of brittle deformation in the marble matrix and investigated the effect of brittle precursors on the initiation and development of heterogeneity-induced high temperature shear zones.
Viscosity contrast between the matrix and the inclusions induces local stress concentration at the tips of these latter. The initial arrangement of the inclusions results in either an overpressured (contractional) or underpressured (extensional) domain in the step-over region of the sample. At low confinement (30 and 50 MPa) abundant brittle deformation is observed, but the spatial distribution of microfractures is dependent on the kinematics of the step-over region: microcracks occur either along the shearing plane between inclusions (in extensional bridge samples), or broadly distributed outside the step-over region (contractional bridge samples). Accordingly, ductile deformation localizes along the inclusions plane in the extensional bridge samples as opposed to distributing over large areas of the matrix in the contractional bridge samples. If microcracking is suppressed (high confinement), strain is accommodated by viscous creep and strain progressively de-localizes in extensional bridge samples. Our experiments demonstrate that brittle precursors enhance the degree of localization in the ductile deformation regime, but only if the interaction of pre-existing heterogeneities induces an extensional mean stress regime in between.
We show that near-real-time seismic monitoring of fluid injection allowed control of induced earthquakes during the stimulation of a 6.1-km-deep geothermal well near Helsinki, Finland. A total of 18,160 m(3) of fresh water was pumped into crystalline rocks over 49 days in June to July 2018. Seismic monitoring was performed with a 24-station borehole seismometer network. Using near-real-time information on induced-earthquake rates, locations, magnitudes, and evolution of seismic and hydraulic energy, pumping was either stopped or varied-in the latter case, between well-head pressures of 60 and 90 MPa and flow rates of 400 and 800 liters/min. This procedure avoided the nucleation of a project-stopping magnitude M-W 2.0 induced earthquake, a limit set by local authorities. Our results suggest a possible physics-based approach to controlling stimulation-induced seismicity in geothermal projects.
The long-term temporal and spatial changes in statistical, source, and stress characteristics of one cluster of induced seismicity recorded at The Geysers geothermal field (U.S.) are analyzed in relation to the field operations, fluid migration, and constraints on the maximum likely magnitude. Two injection wells, Prati-9 and Prati-29, located in the northwestern part of the field and their associated seismicity composed of 1776 events recorded throughout a 7year period were analyzed. The seismicity catalog was relocated, and the source characteristics including focal mechanisms and static source parameters were refined using first-motion polarity, spectral fitting, and mesh spectral ratio analysis techniques. The source characteristics together with statistical parameters (b value) and cluster dynamics were used to investigate and understand the details of fluid migration scheme in the vicinity of injection wells. The observed temporal, spatial, and source characteristics were clearly attributed to fluid injection and fluid migration toward greater depths, involving increasing pore pressure in the reservoir. The seasonal changes of injection rates were found to directly impact the shape and spatial extent of the seismic cloud. A tendency of larger seismic events to occur closer to injection wells and a correlation between the spatial extent of the seismic cloud and source sizes of the largest events was observed suggesting geometrical constraints on the maximum likely magnitude and its correlation to the average injection rate and volume of fluids present in the reservoir.
The mechanical behavior of the sandy facies of Opalinus Clay (OPA) was investigated in 42 triaxial tests performed on dry samples at unconsolidated, undrained conditions at confining pressures (p(c)) of 50-100 MPa, temperatures (T) between 25 and 200 degrees C and strain rates (epsilon) (over dot ) of 1 x-10(-3)-5 x-10(-6) -s(-1). Using a Paterson-type deformation apparatus, samples oriented at 0 degrees, 45 degrees and 90 degrees to bedding were deformed up to about 15% axial strain. Additionally, the influence of water content, drainage condition and pre-consolidation was investigated at fixed p(c)-T conditions, using dry and re-saturated samples. Deformed samples display brittle to semi-brittle deformation behavior, characterized by cataclastic flow in quartz-rich sandy layers and granular flow in phyllosilicate-rich layers. Samples loaded parallel to bedding are less compliant compared to the other loading directions. With the exception of samples deformed 45 degrees and 90 degrees to bedding at p(c) = 100 MPa, strain is localized in discrete shear zones. Compressive strength (sigma(max)) increases with increasing pc, resulting in an internal friction coefficient of approximate to 0.31 for samples deformed at 45 degrees and 90 degrees to bedding, and approximate to 0.44 for samples deformed parallel to bedding. In contrast, pre-consolidation, drainage condition, T and epsilon(over dot )do not significantly affect deformation behavior of dried samples. However, sigma(max) and Young's modulus (E) decrease substantially with increasing water saturation. Compared to the clay-rich shaly facies of OPA, sandy facies specimens display higher strength sigma(max) and Young's modulus E at similar deformation conditions. Strength and Young's modulus of samples deformed 90 degrees and 45 degrees to bedding are close to the iso-stress Reuss bound, suggesting a strong influence of weak clay-rich layers on the deformation behavior.
Knowledge of pressure-dependent static and dynamic moduli of porous reservoir rocks is of key importance for evaluating geological setting of a reservoir in geo-energy applications. We examined experimentally the evolution of static and dynamic bulk moduli for porous Bentheim sandstone with increasing confining pressure up to about 190 MPa under dry and water-saturated conditions. The static bulk moduli (K-s) were estimated from stress-volumetric strain curves while dynamic bulk moduli (K-d) were derived from the changes in ultrasonic P- and S- wave velocities (similar to 1 MHz) along different traces, which were monitored simultaneously during the entire deformation. In conjunction with published data of other porous sandstones (Berea, Navajo and Weber sandstones), our results reveal that the ratio between dynamic and static bulk moduli (K-d/K-s) reduces rapidly from about 1.5 - 2.0 at ambient pressure to about 1.1 at high pressure under dry conditions and from about 2.0 - 4.0 to about 1.5 under water-saturated conditions, respectively. We interpret such a pressure-dependent reduction by closure of narrow (compliant) cracks, highlighting thatK(d)/K(s)is positively correlated with the amount of narrow cracks. Above the crack closure pressure, where equant (stiff) pores dominate the void space,K-d/K(s)is almost constant. The enhanced difference between dynamic and static bulk moduli under water saturation compared to dry conditions is possibly caused by high pore pressure that is locally maintained if measured using high-frequency ultrasonic wave velocities. In our experiments, the pressure dependence of dynamic bulk modulus of water-saturated Bentheim sandstone at effective pressures above 5 MPa can be roughly predicted by both the effective medium theory (Mori-Tanaka scheme) and the squirt-flow model. Static bulk moduli are found to be more sensitive to narrow cracks than dynamic bulk moduli for porous sandstones under dry and water-saturated conditions.
New Zealand's Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (M-w > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.
Rock deformation at depths in the Earth’s crust is often localized in high temperature shear zones occurring at different scales in a variety of lithologies. The presence of material heterogeneities is known to trigger shear zone development, but the mechanisms controlling initiation and evolution of localization are not fully understood. To investigate the effect of loading conditions on shear zone nucleation along heterogeneities, we performed torsion experiments under constant twist rate (CTR) and constant torque (CT) conditions in a Paterson-type deformation apparatus. The sample assemblage consisted of cylindrical Carrara marble specimens containing a thin plate of Solnhofen limestone perpendicular to the cylinder’s longitudinal axis. Under experimental conditions (900 °C, 400 MPa confining pressure), samples were plastically deformed and limestone is about 9 times weaker than marble, acting as a weak inclusion in a strong matrix. CTR experiments were performed at maximum bulk shear strain rates of ~ 2*10-4s-1, yielding peak shear stresses of ~ 20 MPa. CT tests were conducted at shear stresses of ~ 20 MPa resulting in bulk shear strain rates of 1-4*10-4s-1. Experiments were terminated at maximum bulk shear strains of ~ 0.3 and 1.0.Strain was localized within the Carrara marble in front of the inclusion in an area of strongly deformed grains and intense grain size reduction. Locally, evidences for coexisting brittle deformation are also observed regardless of the imposed loading conditions. The local shear strain at the inclusion tipis up to 30 times higher than the strain in the adjacent host rock, rapidly dropping to 5times higher at larger distance from the inclusion. At both investigated bulk strains, the evolution of microstructural and textural parameters is independent of loading conditions. Ourresults suggest that loading conditions do not significantly affect material heterogeneity-induced strain localization during its nucleation and transient stages.
In this paper, an underground experiment at the Aspo Hard Rock Laboratory (HRL) is described. Main goal is optimizing geothermal heat exchange in crystalline rock mass at depth by multistage hydraulic fracturing with minimal impact on the environment, that is, seismic events. For this, three arrays with acoustic emission, microseismicity and electromagnetic sensors are installed mapping hydraulic fracture initiation and growth. Fractures are driven by three different water injection schemes (continuous, progressive and pulse pressurization). After a brief review of hydraulic fracture operations in crystalline rock mass at mine scale, the site geology and the stress conditions at Aspo HRL are described. Then, the continuous, single-flow rate and alternative, multiple-flow rate fracture breakdown tests in a horizontal borehole at depth level 410 m are described together with the monitoring networks and sensitivity. Monitoring results include the primary catalogue of acoustic emission hypocentres obtained from four hydraulic fractures with the in situ trigger and localizing network. The continuous versus alternative water injection schemes are discussed in terms of the fracture breakdown pressure, the fracture pattern from impression packer result and the monitoring at the arrays. An example of multistage hydraulic fracturing with several phases of opening and closing of fracture walls is evaluated using data from acoustic emissions, seismic broad-band recordings and electromagnetic signal response. Based on our limited amount of in situ tests (six) and evaluation of three tests in Avro granodiorite, in the multiple-flow rate test with progressively increasing target pressure, the acoustic emission activity starts at a later stage in the fracturing process compared to the conventional fracturing case with continuous water injection. In tendency, also the total number and magnitude of acoustic events are found to be smaller in the progressive treatment with frequent phases of depressurization.