Impact of injection rate ramp-up on nucleation and arrest of dynamic fault slip

Fluid injection into underground formations reactivates preexisting geological discontinuities such as faults or fractures. In this work, we investigate the impact of injection rate ramp-up present in many standard injection protocols on the nucleation and potential arrest of dynamic slip along a planar pressurized fault. We assume a linear increasing function of injection rate with time, up to a given time $$t_c$$ t c after which a maximum value $$Q_m$$ Q m is achieved. Under the assumption of negligible shear-induced dilatancy and impermeable host medium, we solve numerically the coupled hydro-mechanical model and explore the different slip regimes identified via scaling analysis. We show that in the limit when fluid diffusion time scale $$t_w$$ t w is much larger than the ramp-up time scale $$t_c$$ t c , slip on an ultimately stable fault is essentially driven by pressurization at constant rate. Vice versa, in the limit when $$t_c/t_w \gg 1$$ t c / t w ≫ 1 , the pressurization rate, quantified by the dimensionless ratio $$\dfrac{Q_m t_w}{t_c Q^*}$$ Q m t w t c Q ∗ with $$Q^*$$ Q ∗ being a characteristic injection rate scale, does impact both nucleation time and arrest distance of dynamic slip. Indeed, for a given initial fault loading condition and frictional weakening property, lower pressurization rates delay the nucleation of a finite-sized dynamic event and increase the corresponding run-out distance approximately proportional to $$\propto \left( \dfrac{Q_m t_w}{t_c Q^*}\right) ^{-0.472}$$ ∝ Q m t w t c Q ∗ - 0.472 . On critically stressed faults, instead, the ramp-up of injection rate activates quasi-static slip which quickly turn into a run-away dynamic rupture. Its nucleation time decreases non-linearly with increasing value of $$\dfrac{Q_m t_w}{t_c Q^*}$$ Q m t w t c Q ∗ and it may precede (or not) the one associated with fault pressurization at constant rate only.

Dynamic rupture initiation and propagation in a fluid-injection laboratory setup with diagnostics across multiple temporal scales

Fluids are known to trigger a broad range of slip events, from slow, creeping transients to dynamic earthquake ruptures. Yet, the detailed mechanics underlying these processes and the conditions leading to different rupture behaviors are not well understood. Here, we use a laboratory earthquake setup, capable of injecting pressurized fluids, to compare the rupture behavior for different rates of fluid injection, slow (megapascals per hour) versus fast (megapascals per second). We find that for the fast injection rates, dynamic ruptures are triggered at lower pressure levels and over spatial scales much smaller than the quasistatic theoretical estimates of nucleation sizes, suggesting that such fast injection rates constitute dynamic loading. In contrast, the relatively slow injection rates result in gradual nucleation processes, with the fluid spreading along the interface and causing stress changes consistent with gradually accelerating slow slip. The resulting dynamic ruptures propagating over wetted interfaces exhibit dynamic stress drops almost twice as large as those over the dry interfaces. These results suggest the need to take into account the rate of the pore-pressure increase when considering nucleation processes and motivate further investigation on how friction properties depend on the presence of fluids.

Comparative analysis of acoustic and electromagnetic emissions of rocks

Comparative analysis of acoustic and electromagnetic emissions recorded during the intact rock samples deformation and dynamic rupture of simulated crustal fault is presented. Specialized machines for uniaxial compression and shear tests of rock samples with identical data acquisition systems for both test cases were employed. Increase of acoustic emission was observed accompanied by significant rise of intensity and amplitude of electromagnetic signals at high stress of the rock samples under the uniaxial compression or dynamic failure in the spring-block model. Such correlation is consistent with the previous conclusions that an increase of electromagnetic emission may be considered as a rock failure precursor. Any specific characteristics of the detected electromagnetic signals to be used for prediction of impending rock failure or the earthquake fault rupture were not found. The similarity of electromagnetic signals and their spectra obtained at the press equipment and the spring-block model suggests that in both cases, the signals observed are generated by the crack formations and shear. The electromagnetic emission appeared only in dry samples. The samples saturated by water with the salinity of over 0.1% demonstrated no electromagnetic emission.

Unified Reservoir And Seismic Simulation With Explicit Representation Of Fractures And Faults

In the context of remote sensing, the vast disparity in characteristic scales between seismic deformation (e.g. milliseconds) and transient flow (e.g. hours) allows a "two-model paradigm" for geophysics and reservoir simulation. In the context of flow-induced geohazard risk mitigation and micro-seismic data integration, this paradigm breaks down. Under micro-seismic deformation, events occur with high-frequency, and over sustained duration during which the rock-fluid coupling is significant. In risk mitigation scenarios, the onset of seismic deformation is directly tied to quasi-static coupling periods. This work develops an approach to reservoir simulation modeling that allows simultaneous resolution of transient (inertial) poromechanics and multiphase fluid flow in the presence of fracture. A mixed discretization scheme combining the extended finite element method (XFEM) and the embedded discrete fracture model (EDFM) is extended using a second-order implicit Newmark time integration scheme for the inertial mechanics. A Lagrange multiplier method is developed to model pressure-dependent contact traction in fractures. The contact constraints are adapted to accommodate fracture opening. Slip-weakening fracture friction models are incorporated. Finally, a time-step controller is proposed to combine local discretization error with contact traction and slip-rate control along the fractures. This strategy allows automatic adaptation to resolve quasi-static, inter-seismic triggering, and co-seismic spontaneous rupture periods within one model. The model is verified to simulate complete induced earthquake sequences, including inter-seismic and dynamic rupture phases. The performance of the adaptive model is illustrated for cases with various set-ups of production and injection periods in a fractured reservoir with explicit fracture representation.

Effects of the Rock Bridge Ligament on Fracture and Energy Evolution of Preflawed Granite Exposed to Dynamic Loads

This paper aims to reveal the mechanical properties, energy evolution characteristics, and dynamic rupture process of preflawed granite under impact loading with different rock bridge angles and strain rates. A series of dynamic impact experiments were conducted along with the separate Hopkinson press bar (SHPB) testing system to analyze and study the overall rock fracture process. Under the impact load, the peak stress of granite increases with the increase of rock bridge angle and strain rate, but the increase gradually decreases. The peak strain also increases gradually with the increase of rock bridge angle, but there is an upper limit value; the total input strain energy increases with the increase of strain rate and rock bridge angle. It is shown that the higher the strain rate, the higher the unit dissipation energy, and the greater the degree of rock fragmentation. For rock under impact loads, the crack first initiates from the wing end of the prefabricated flaw, the preflaw closes gradually, and finally the crack propagates at the locking section leading to the coalescence of rock bridge. With the increase of strain rate, the fragmentation degree of the specimen increases asymptotically, and the average fragmentation size of the specimen decreases with the increase of strain rate. It is suggested that the stability of large rocked slopes is controlled by the locked section, and understanding the fracture evolution of the rock bridge is the key to slope instability prediction.

Fully Dynamic Earthquake Cycle Simulations on a Nonplanar Fault Using the Spectral Boundary Integral Element Method (sBIEM)

ABSTRACT One of the most suitable methods for modeling fully dynamic earthquake cycle simulations is the spectral boundary integral element method (sBIEM), which takes advantage of the fast Fourier transform (FFT) to make a complex numerical dynamic rupture tractable. However, this method has the serious drawback of requiring a flat fault geometry due to the FFT approach. Here, we present an analytical formulation that extends the sBIEM to a mildly nonplanar fault. We start from a regularized boundary element method and apply a small-slope approximation of the fault geometry. Making this assumption, it is possible to show that the main effect of nonplanar fault geometry is to change the normal traction along the fault, which is controlled by the local curvature along the fault. We then convert this space–time boundary integral equation of the normal traction into a spectral-time formulation and incorporate this change in normal traction into the existing sBIEM methodology. This approach allows us to model fully dynamic seismic cycle simulations on nonplanar faults in a particularly efficient way. We then test this method against a regular BIEM for both rough-fault and seamount-fault geometries and demonstrate that this sBIEM maintains the scaling between the fault geometry and slip distribution.

Assessing the role of selected constraints in Bayesian dynamic source inversion: Application to the 2017 Mw 6.3 Lesvos earthquake

Summary A dynamic finite-fault source inversion for stress and frictional parameters of the Mw 6.3 2017 Lesvos earthquake is carried out. The mainshock occurred on June 12, offshore the southeastern coast of the Greek island of Lesvos in the north Aegean Sea. It caused 1 fatality, 15 injuries, and extensive damage to the southern part of the island. Dynamic rupture evolution is modeled on an elliptic patch, using the linear slip-weakening friction law. The inversion is posed as a Bayesian problem and the Parallel Tempering Markov Chain Monte Carlo algorithm is used to obtain posterior probability distributions by updating the prior distribution with progressively more constraints. To calculate the first posterior distribution, only the constraint that the model should expand beyond the nucleation patch is used. Then, we add the constraint that the model should reach a moment magnitude similar to that obtained from our centroid moment tensor inversion. For the final posterior distribution, 15 acceleration records from Greek and Turkish strong motion networks at near regional distances ($\approx 30 - 150$ km) in the frequency range of 0.05–0.15 Hz are used. The three posterior distributions are compared to understand how much each constraint contributes to resolving different quantities. The most probable values and uncertainties of individual parameters are also calculated, along with their mutual trade-offs. The features best determined by seismograms in the final posterior distribution include the position of the nucleation region, the mean direction of rupture (towards WNW), the mean rupture speed (with 68 per cent of the distribution lying between 1.4–2.6 km/s), radiated energy (12–65 TJ), radiation efficiency (0.09–0.38), and the mean stress drop (2.2–6.5 MPa).

Multicycle Simulation of Strike-Slip Earthquake Rupture for Use in Near-Source Ground-Motion Simulations

ABSTRACT Realistic dynamic rupture modeling validated by observed earthquakes is necessary for estimating parameters that are poorly resolved by seismic source inversion, such as stress drop, rupture velocity, and slip rate function. Source inversions using forward dynamic modeling are increasingly used to obtain earthquake rupture models. In this study, to generate a large number of physically self-consistent rupture models, rupture process of which is consistent with the spatiotemporal heterogeneity of stress produced by previous earthquakes on the same fault, we use multicycle simulations under the rate and state (RS) friction law. We adopt a one-way coupling from multicycle simulations to dynamic rupture simulations; the quasidynamic solver QDYN is used to nucleate the seismic events and the spectral element dynamic solver SPECFEM3D to resolve their rupture process. To simulate realistic seismicity, with a wide range of magnitudes and irregular recurrence, several realizations of 2D-correlated heterogeneous random distributions of characteristic weakening distance (Dc) in RS friction are tested. Other important parameters are the normal stress, which controls the stress drop and rupture velocity during an earthquake, and the maximum value of Dc, which controls rupture velocity but not stress drop. We perform a parametric study on a vertical planar fault and generate a set of a hundred spontaneous rupture models in a wide magnitude range (Mw 5.5–7.4). We validate the rupture models by comparison of source scaling, ground motion (GM), and surface slip properties to observations. We compare the source-scaling relations between rupture area, average slip, and seismic moment of the modeled events with empirical ones derived from source inversions. Near-fault GMs are computed from the source models. Their peak ground velocities and peak ground accelerations agree well with the ground-motion prediction equation values. We also obtain good agreement of the surface fault displacements with observed values.