Investigation of injection-induced seismicity using a coupled fluid flow and rate/state friction model

Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. WC181-WC198 ◽  
Author(s):  
Mark W. McClure ◽  
Roland N. Horne

We describe a numerical investigation of seismicity induced by injection into a single isolated fracture. Injection into a single isolated fracture is a simple analog for shear stimulation in enhanced geothermal systems (EGS) during which water is injected into fractured, low permeability rock, triggering slip on preexisting large scale fracture zones. A model was developed and used that couples (1) fluid flow, (2) rate and state friction, and (3) mechanical stress interaction between fracture elements. Based on the results of this model, we propose a mechanism to describe the process by which the stimulated region grows during shear stimulation, which we refer to as the sequential stimulation (SS) mechanism. If the SS mechanism is realistic, it would undermine assumptions that are made for the estimation of the minimum principal stress and unstimulated hydraulic diffusivity. We investigated the effect of injection pressure on induced seismicity. For injection at constant pressure, there was not a significant dependence of maximum event magnitude on injection pressure, but there were more relatively large events for higher injection pressure. Decreasing injection pressure over time significantly reduced the maximum event magnitude. Significant seismicity occurred after shut-in, which was consistent with observations from EGS stimulations. Production of fluid from the well immediately after injection inhibited shut-in seismic events. The results of the model in this study were found to be broadly consistent with results from prior work using a simpler treatment of friction that we refer to as static/dynamic. We investigated the effect of shear-induced pore volume dilation and the rate and state characteristic length scale, [Formula: see text]. Shear-induced pore dilation resulted in a larger number of lower magnitude events. A larger value of [Formula: see text] caused slip to occur aseismically.

Water ◽  
2020 ◽  
Vol 12 (11) ◽  
pp. 3171
Author(s):  
Liangliang Guo ◽  
Zihong Wang ◽  
Yanjun Zhang ◽  
Zhichao Wang ◽  
Haiyang Jiang

In order to study the mechanism of hydraulic fracturing in enhanced geothermal systems, we analyzed the influence of high temperatures and embedded fractures on the initiation and propagation of hydraulic fractures using a laboratory test and numerical simulation. The analysis was conducted via large-scale true triaxial hydraulic fracturing tests with acoustic emission monitoring. Moreover, we discussed and established the elastic-plastic criterion of hydraulic fracturing initiation. The corresponding fracturing procedure was designed and embedded into the FLAC3D software. Then, a numerical simulation was conducted and compared with the laboratory test to verify the accuracy of the fracturing procedure. The influence of high temperatures on hydraulic fracturing presented the following features. First, multi-fractures were created, especially in the near-well region. Second, fracturing pressure, extension pressure, and fracture flow resistance became larger than those at room temperature. 3D acoustic fracturing emission results indicated that the influence of the spatial distribution pattern of embedded fractures on hydraulic fracturing direction was larger than that of triaxial stress. Furthermore, the fracturing and extension pressures decreased with the increase of embedded fracture density. For hydraulic fracturing in a high temperature reservoir, a plastic zone was generated near the borehole, and this zone increased as the injection pressure increased until the well wall failed.


Geothermics ◽  
2007 ◽  
Vol 36 (3) ◽  
pp. 185-222 ◽  
Author(s):  
Ernest L. Majer ◽  
Roy Baria ◽  
Mitch Stark ◽  
Stephen Oates ◽  
Julian Bommer ◽  
...  

2020 ◽  
Author(s):  
Kyung Won Chang ◽  
Gungor Beskardes ◽  
Chester Weiss

<p>Hydraulic stimulation is the process of initiating fractures in a target reservoir for subsurface energy resource management with applications in unconventional oil/gas and enhanced geothermal systems. The fracture characteristics (i.e., number, size and orientation with respect to the wellbore) determines the modified permeability field of the host rock and thus, numerical simulations of flow in fractured media are essential for estimating the anticipated change in reservoir productivity. However, numerical modeling of fluid flow in highly fractured media is challenging due to the explosive computational cost imposed by the explicit discretization of fractures at multiple length scales. A common strategy for mitigating this extreme cost is to crudely simplify the geometry of fracture network, thereby neglecting the important contributions made by all elements of the complex fracture system.</p><p>The proposed “Hierarchical Finite Element Method” (Hi-FEM; Weiss, Geophysics, 2017) reduces the comparatively insignificant dimensions of planar- and curvilinear-like features by translating them into integrated hydraulic conductivities, thus enabling cost-effective simulations with requisite solutions at material discontinuities without defining ad-hoc, heuristic, or empirically-estimated boundary conditions between fractures and the surrounding formation. By representing geometrical and geostatistical features of a given fracture network through the Hi-FEM computational framework, geometrically- and geomechanically-dependent fluid flow properly can now be modeled economically both within fractures as well as the surrounding medium, with a natural “physics-informed” coupling between the two.</p><p>SNL is managed and operated by NTESS under DOE NNSA contract DE-NA0003525.</p>


2021 ◽  
Author(s):  
Saumik Dana

The effect of fluid pulse driven fractures (FPDF) propagating in poroelastic media on fault slip in the presence of natural fractures is a complicated interplay between fracture propagation, fracture-fracture interaction, fracture-fault interaction, friction model governing fault slip and wave propagation associated with pulsing injection. Furthermore, the problem is stochastic due to the uncertainty associated with the existing fracture-fault topology.


2014 ◽  
Author(s):  
Roland Gritto ◽  
◽  
Douglas Dreger ◽  
Oliver Heidbach ◽  
Lawrence Hutchings ◽  
...  

2012 ◽  
Author(s):  
Ernie Majer ◽  
James Nelson ◽  
Ann Robertson-Tait ◽  
Jean Savy ◽  
Ivan Wong

2020 ◽  
Author(s):  
Nathaniel Forbes Inskip ◽  
Tomos Phillips ◽  
Kevin Bisdom ◽  
Georgy Borisochev ◽  
Andreas Busch ◽  
...  

<p>Fractures are ubiquitous in geological sequences, and play an important role in the movement of fluids in the earth’s crust, particularly in fields such as hydrogeology, petroleum geology and volcanology. When predicting or analysing fluid flow, fractures are often simplified as a set of smooth parallel plates. In reality, they exhibit tortuosity on a number of scales: Fine-scale tortuosity, or roughness, is the product of the small-scale (µm – mm) irregularities in the fracture surface, whereas large-scale (> mm) tortuosity occurs as a result of anisotropy and heterogeneity within the host formation that leads to the formation of irregularities in the fracture surfaces. It is important to consider such tortuosity when analysing processes that rely on the movement (or hindrance) of fluids flowing through fractures in the subsurface. Such processes include fluid injection into granitic plutons for the extraction of heat in Engineered Geothermal Systems, or the injection of CO<sub>2</sub> into reservoirs overlain by fine-grained mudrocks acting as seals in Carbon Capture and Storage projects.</p><p>Although it is generally assumed that tortuosity is controlled by factors such as grain size, mineralogy and fracture mode, a systematic study of how these factors quantitatively affect tortuosity is currently lacking. Furthermore, in anisotropic rocks the fracture orientation with respect to any inherent anisotropy is also likely to affect tortuosity.</p><p>In order to address this gap, we have induced fractures in a selection of different rock types (mudrocks, sandstones and carbonates) using the Brazil disk method, and imaged the fracture surfaces using both a digital optical microscope and X-ray Computed Tomography. Using these methods we are able to characterise both the fine-scale (roughness) and large-scale tortuosity. In order to understand the effect of fracture orientation on tortuosity we have also analysed fractures induced at different angles to bedding in samples of a highly anisotropic mudrock taken from South Wales, UK. Results indicate that fine-scale tortuosity is highly dependent on the fracture orientation with regards to the bedding plane, with fractures normal to bedding being rougher than those induced parallel to bedding. Finally, in order to measure the effect of tortuosity on fluid flow, we have carried out a series of core flooding experiments on a subset of fractured samples showing that fracture transmissivity decreases with increasing tortuosity.</p>


Sign in / Sign up

Export Citation Format

Share Document