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

2021 ◽  
Vol 118 (51) ◽  
pp. e2023433118
Author(s):  
Marcello Gori ◽  
Vito Rubino ◽  
Ares J. Rosakis ◽  
Nadia Lapusta

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.

SPE Journal ◽  
2015 ◽  
Vol 20 (04) ◽  
pp. 689-700 ◽  
Author(s):  
S.. Ameen ◽  
A. Dahi Taleghani

Summary Injectivity loss is a common problem in unconsolidated-sand formations. Injection of water into a poorly cemented granular medium may lead to internal erosion, and consequently formation of preferential flow paths within the medium because of channelization. Channelization in the porous medium might occur when fluid-induced stresses become locally larger than a critical threshold and small grains are dislodged and carried away; hence, porosity and permeability of the medium will evolve along the induced flow paths. Vice versa, flowback during shut-in might carry particles back to the well and cause sand accumulation inside the well, and subsequently loss of injectivity. In most cases, to maintain the injection rate, operators will increase injection pressure and pumping power. The increased injection pressure results in stress changes and possibly further changes in channel patterns around the wellbore. Experimental laboratory studies have confirmed the presence of the transition from uniform Darcy flow to a fingered-pattern flow. To predict these phenomena, a model is needed to fill this gap by predicting the formation of preferential flow paths and their evolution. A model based on the multiphase-volume-fraction concept is used to decompose porosity into mobile and immobile porosities where phases may change spatially, evolve over time, and lead to development of erosional channels depending on injection rates, viscosity, and rock properties. This model will account for both particle release and suspension deposition. By use of this model, a methodology is proposed to derive model parameters from routine injection tests by inverse analysis. The proposed model presents the characteristic behavior of unconsolidated formation during fluid injection and the possible effect of injection parameters on downhole-permeability evolution.


2020 ◽  
Author(s):  
Valère Lambert ◽  
Nadia Lapusta

Abstract. Substantial insight into earthquake source processes has resulted from considering frictional ruptures analogous to cohesive-zone shear cracks from fracture mechanics. This analogy holds for slip-weakening representations of fault friction that encapsulate the resistance to rupture propagation in the form of breakdown energy, analogous to fracture energy, prescribed in advance as if it were a material property of the fault interface. Here, we use numerical models of earthquake sequences with enhanced weakening due to thermal pressurization of pore fluids to show how accounting for thermo-hydro-mechanical processes during dynamic shear ruptures makes breakdown energy rupture-dependent. We find that local breakdown energy is neither a constant material property nor uniquely defined by the amount of slip attained during rupture, but depends on how that slip is achieved through the history of slip rate and dynamic stress changes during the rupture process. As a consequence, the frictional breakdown energy of the same location along the fault can vary significantly in different earthquake ruptures that pass through. These results suggest the need for re-examining the assumption of pre-determined frictional breakdown energy common in dynamic rupture modeling and for better understanding of the factors that control rupture dynamics in the presence of thermo-hydro-mechanical processes.


Science ◽  
2019 ◽  
Vol 364 (6439) ◽  
pp. 464-468 ◽  
Author(s):  
Pathikrit Bhattacharya ◽  
Robert C. Viesca

Earthquake swarms attributed to subsurface fluid injection are usually assumed to occur on faults destabilized by increased pore-fluid pressures. However, fluid injection could also activate aseismic slip, which might outpace pore-fluid migration and transmit earthquake-triggering stress changes beyond the fluid-pressurized region. We tested this theoretical prediction against data derived from fluid-injection experiments that activated and measured slow, aseismic slip on preexisting, shallow faults. We found that the pore pressure and slip history imply a fault whose strength is the product of a slip-weakening friction coefficient and the local effective normal stress. Using a coupled shear-rupture model, we derived constraints on the hydromechanical parameters of the actively deforming fault. The inferred aseismic rupture front propagates faster and to larger distances than the diffusion of pressurized pore fluid.


SPE Journal ◽  
2019 ◽  
Vol 24 (04) ◽  
pp. 1839-1855 ◽  
Author(s):  
Bing Hou ◽  
Zhi Chang ◽  
Weineng Fu ◽  
Yeerfulati Muhadasi ◽  
Mian Chen

Summary Deep shale gas reservoirs are characterized by high in-situ stresses, a high horizontal-stress difference (12 MPa), development of bedding seams and natural fractures, and stronger plasticity than shallow shale. All of these factors hinder the extension of hydraulic fractures and the formation of complex fracture networks. Conventional hydraulic-fracturing techniques (that use a single fluid, such as guar fluid or slickwater) do not account for the initiation and propagation of primary fractures and the formation of secondary fractures induced by the primary fractures. For this reason, we proposed an alternating-fluid-injection hydraulic-fracturing treatment. True triaxial hydraulic-fracturing tests were conducted on shale outcrop specimens excavated from the Shallow Silurian Longmaxi Formation to study the initiation and propagation of hydraulic fractures while the specimens were subjected to an alternating fluid injection with guar fluid and slickwater. The initiation and propagation of fractures in the specimens were monitored using an acoustic-emission (AE) system connected to a visual display. The results revealed that the guar fluid and slickwater each played a different role in hydraulic fracturing. At a high in-situ stress difference, the guar fluid tended to open the transverse fractures, whereas the slickwater tended to activate the bedding planes as a result of the temporary blocking effect of the guar fluid. On the basis of the development of fractures around the initiation point, the initiation patterns were classified into three categories: (1) transverse-fracture initiation, (2) bedding-seam initiation, and (3) natural-fracture initiation. Each of these fracture-initiation patterns had a different propagation mode. The alternating-fluid-injection treatment exploited the advantages of the two fracturing fluids to form a large complex fracture network in deep shale gas reservoirs; therefore, we concluded that this method is an efficient way to enhance the stimulated reservoir volume compared with conventional hydraulic-fracturing technologies.


Author(s):  
H. Lim ◽  
K. Deng ◽  
Y.H. Kim ◽  
J.‐H. Ree ◽  
T.‐R. A. Song ◽  
...  

2020 ◽  
Vol 110 (5) ◽  
pp. 2350-2365 ◽  
Author(s):  
Dawid Szafranski ◽  
Benchun Duan

ABSTRACT In this work, we integrate a fluid-flow model of 3D deformable porous media with a dynamic rupture model of earthquakes in 3D heterogeneous geologic medium. The method allows us to go beyond fault failure potential analyses and to examine how big an earthquake can be if part of a fault reaches failure due to fluid injection. We apply the method to the 17 May 2012 Mw 4.8 Timpson, Texas, earthquake as a case study. The simulated perturbations of pore pressure and stress from wastewater injection at the time of the mainshock are high enough (several MPa) to trigger an earthquake. Dynamic rupture modeling could reproduce the major observations from the Mw 4.8 event, including its size, focal mechanism, and aftershock sequence, and thus building a more convincing physical link between fluid injection and the Mw 4.8 earthquake. Furthermore, parameter space studies of dynamic rupture modeling allow us to place some constraints on fault frictional properties and background stresses. For the Timpson case, we find that a dynamic friction coefficient of ∼0.3, a value of ∼0.1  m for the critical slip distance in the slip-weakening friction law, and uniform effective normal stress are associated with the Timpson earthquake fault. By reproducing main features of the aftershock sequence of the mainshock, we also demonstrate that the method has potential to become a predictive tool for fluid injection design in the future.


2020 ◽  
Vol 39 (12) ◽  
pp. 893-900
Author(s):  
Inga Berre ◽  
Ivar Stefansson ◽  
Eirik Keilegavlen

Hydraulic stimulation of geothermal reservoirs in low-permeability basement and crystalline igneous rock can enhance permeability by reactivation and shear dilation of existing fractures. The process is characterized by interaction between fluid flow, deformation, and the fractured structure of the formation. The flow is highly affected by the fracture network, which in turn is deformed because of hydromechanical stress changes caused by the fluid injection. This process-structure interaction is decisive for the outcome of hydraulic stimulation, and, in analysis of governing mechanisms, physics-based modeling has potential to complement field and experimental data. Here, we show how recently developed simulation technology is a valuable tool to understand governing mechanisms of hydromechanical coupled processes and the reactivation and deformation of faults. The methodology fully couples flow in faults and matrix with poroelastic matrix deformation and a contact mechanics model for the faults, including dilation because of slip. Key elements are high aspect ratios of faults and strong nonlinearities in highly coupled governing equations. Example simulations using our open-source software illustrate direct and indirect hydraulic fault reactivation and corresponding permeability enhancement. We investigate the effect of the fault and matrix permeability and the Biot coefficient. A higher matrix permeability leads to more leakage from a permeable fault and thus suppresses reactivation and slip of the fault compared to the case with a lower matrix permeability. If a fault is a barrier to flow, increase of pressure because of the fluid injection results in stabilization of the fault; the situation is opposite if the fault is highly permeable compared to the matrix. For the given setup, lowering the Biot coefficient results in more slip than the base case. While conceptually simple, the examples illustrate the strong hydromechanical couplings and the prospects of physics-based numerical models in investigating the dynamics.


2019 ◽  
Vol 109 (6) ◽  
pp. 2168-2186 ◽  
Author(s):  
Paul Peshette ◽  
Julian Lozos ◽  
Doug Yule ◽  
Eileen Evans

Abstract Investigations of historic surface‐rupturing thrust earthquakes suggest that rupture can jump from one fault to another up to 8 km away. Additionally, there are observations of jumping rupture between thrust faults ∼50  km apart. In contrast, previous modeling studies of thrust faults find a maximum jumping rupture distance of merely 0.2 km. Here, we present a dynamic rupture modeling parameter study that attempts to reconcile these differences and determines geometric and stress conditions that promote jumping rupture. We use the 3D finite‐element method to model rupture on pairs of thrust faults with parallel surface traces and opposite dip orientations. We vary stress drop and fault strength ratio to determine conditions that produce jumping rupture at different dip angles and different minimum distance between faults. We find that geometry plays an essential role in determining whether or not rupture will jump to a neighboring thrust fault. Rupture is more likely to jump between faults dipping toward one another at steeper angles, and the behavior tapers down to no rupture jump in shallow dip cases. Our variations of stress parameters emphasize these toward‐orientation results. Rupture jump in faults dipping away from one another is complicated by variations of stress conditions, but the most prominent consistency is that for mid‐dip angle faults rupture rarely jumps. If initial stress conditions are such that they are already close to failure, the possibility of a long‐distance jump increases. Our models call attention to specific geometric and stress conditions where the dynamic rupture front is the most important to potential for jumping rupture. However, our models also highlight the importance of near‐field stress changes due to slip. According to our modeling, the potential for rupture to jump is strongly dependent on both dip angle and orientation of faults.


Sign in / Sign up

Export Citation Format

Share Document