scholarly journals Thermo-Poroelastic Analysis of Induced Seismicity at the Basel Enhanced Geothermal System

2019 ◽  
Vol 11 (24) ◽  
pp. 6904 ◽  
Author(s):  
Sandro Andrés ◽  
David Santillán ◽  
Juan Carlos Mosquera ◽  
Luis Cueto-Felgueroso

Geothermal energy has emerged as an alternative to ensure a green energy supply while tackling climate change. Geothermal systems extract the heat stored in the Earth’s crust by warming up water, but the low rock permeability at exploitation depths may require the hydraulic stimulation of the rock fracture network. Enhanced Geothermal Systems (EGS) employ techniques such as hydro-shearing and hydro-fracturing for that purpose, but their use promotes anthropogenic earthquakes induced by the injection or extraction of fluids. This work addresses this problem through developing a computational 3D model to explore fault reactivation and evaluating the potential for earthquake triggering at preexisting geological faults. These are included in the model as frictional contacts that allow the relative displacement between both of its sides, governed by rate-and-state friction laws and fully coupled with thermo-hydro-mechanical equations. We apply our methodology to the Basel project, employing the on-site parameters and conditions. Our results demonstrate that earthquakes which occurred in December 2006 in Basel (Switzerland) are compatible with the geomechanical and frictional consequences of the hydraulic stimulation of the rock mass. The application of our model also shows that it can be useful for predicting fault reactivation and engineering injection protocols for managing the safe and sustainable operation of EGS.

Author(s):  
Dustin Crandall ◽  
Goodarz Ahmadi ◽  
Grant Bromhal

Fractures in rocks enable the motion of fluids through the large, hot geologic formations of geothermal reservoirs. The heat transfer from the surrounding rock mass to the fluid flowing through a fracture depends on the geometry of the fracture, the fluid/solid properties, and the flow rate through the fracture. A numerical study was conducted to evaluate the changes in heat transfer to the fluid flowing through a rock fracture with changes in the flow rate. The aperture distribution of the rock fracture, originally created within Berea sandstone and imaged using a CT-scanner, is well described by a Gaussian distribution and has a mean aperture of approximately 0.6 mm. Water was used as the working fluid, enabling an evaluation of the efficiency of heat flux to the fluid along the flow path of a hot dry geothermal system. As the flow through the fracture was increased to a Reynolds number greater than 2300 the effect of channeling through large aperture regions within the fracture were observed to become increasingly important. For the fastest flows modeled the heat flux to the working fluids was reduced due to a shorter residence time of the fluid in the fracture. Understanding what conditions can maximize the amount of energy obtained from fractures within a hot dry geologic field can improve the operation and long-term viability of enhanced geothermal systems.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. B121-B130 ◽  
Author(s):  
Jared R. Peacock ◽  
Stephan Thiel ◽  
Graham S. Heinson ◽  
Peter Reid

Realization of enhanced geothermal systems (EGSs) prescribes the need for novel methods to monitor subsurface fracture connectivity and fluid distribution. Magnetotellurics (MT) is a passive electromagnetic (EM) method sensitive to electrical conductivity contrasts as a function of depth, specifically hot saline fluids in a resistive porous media. In July 2011, an EGS fluid injection at 3.6-km depth near Paralana, South Australia, was monitored by comparing repeated MT surveys before and after hydraulic stimulation. An observable coherent change above measurement error in the MT response was present and causal, in that variations in phase predict variations in apparent resistivity. Phase tensor residuals proved the most useful representation for characterizing alterations in subsurface resistivity structure, whereas resistivity tensor residuals aided in determining the sign and amplitude of resistivity variations. These two tensor representations of the residual MT response suggested fluids migrated toward the northeast of the injection well along an existing fault system trending north-northeast. Forward modeling and concurrent microseismic data support these results, although microseismic data suggest fractures opened along two existing fracture networks trending north-northeast and northeast. This exemplifies the need to use EM methods for monitoring fluid injections due to their sensitivity to conductivity contrasts.


2021 ◽  
Author(s):  
Kwang-Il Kim ◽  
Hwajung Yoo ◽  
Seheok Park ◽  
Juhyi Yim ◽  
Linmao Xie ◽  
...  

<p>Hydraulic stimulation for the creation of an Enhanced Geothermal System (EGS) reservoir could potentially reactivate a nearby fault and result in man-made earthquakes. In November 15, 2017, an M<sub>w</sub> 5.5 earthquake, the second largest after the initiation of the South Korean national instrumental monitoring system, occurred near an EGS project in Pohang, South Korea. The earthquake occurred on a previously unmapped fault, that is here denoted the M<sub>w</sub> 5.5 Fault. A number of previous studies to model the hydraulic stimulation in the Pohang EGS project have been carried out to identify the mechanism of seismic events. Those previous studies focused on coupled hydro-mechanical processes without the consideration of pre-existing fractures and thermal effects. This study presents an investigation of the mechanisms of induced and triggered seismicity in the Pohang EGS project through three-dimensional coupled thermo-hydro-mechanical numerical simulations. Fractures intersecting the open-hole sections of two deep boreholes, PX-1 and PX-2, clearly indicated by field observations are modeled along with the M<sub>w</sub> 5.5 Fault. Models of stress-dependent permeability models are calibrated based on the numerical reproduction of the pressure-time evolution during the field hydraulic stimulations. The Coulomb failure stress change at the M<sub>w</sub> 5.5 Fault is calculated to quantify the impact of five hydraulic stimulations. In the case of PX-2 stimulations, the pore pressure buildup results in a volumetric expansion of the reservoir and thereby the perturbation of stresses is transferred to the M<sub>w</sub> 5.5 Fault. The volumetric contraction of the reservoir by the temperature reduction could slightly perturb the stress distribution at the M<sub>w</sub> 5.5 Fault. In the case of PX-1 stimulations, shear slip of the PX-1 fracture is explicitly modeled. The modeling shows that transfer of the shear stress drop by the shear slip stabilizes the M<sub>w</sub> 5.5 Fault, which is consistent with the field observation that the seismicity was not induced at the M<sub>w</sub> 5.5 Fault by the PX-1 stimulations. The cooling-induced thermal stress additionally reduces the effective normal stress of PX-1 fracture. Thus, some additional shear slip of the PX-1 fracture is induced by the thermal effect. However, the modeling shows that for both PX-1 and PX-2 stimulations, thermally-induced stress perturbations are very small compared to pressure-induced stress perturbations.</p>


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>


2020 ◽  
Vol 205 ◽  
pp. 03007
Author(s):  
Yejin Kim ◽  
Seong Jun Ha ◽  
Tae sup Yun

Hydraulic stimulation has been a key technique in enhanced geothermal systems (EGS) and the recovery of unconventional hydrocarbon resources to artificially generate fractures in a rock formation. Previous experimental studies present that the pattern and aperture of generated fractures vary as the fracking pressure propagation. The recent development of three-dimensional X-ray computed tomography allows visualizing the fractures for further analysing the morphological features of fractures. However, the generated fracture consists of a few pixels (e.g., 1-3 pixels) so that the accurate and quantitative extract of micro-fracture is highly challenging. Also, the high-frequency noise around the fracture and the weak contrast across the fracture makes the application of conventional segmentation methods limited. In this study, we adopted an encoder-decoder network with a convolutional neural network (CNN) based on deep learning method for the fast and precise detection of micro-fractures. The conventional image processing methods fail to extract the continuous fractures and overestimate the fracture thickness and aperture values while the CNN-based approach successfully detects the barely seen fractures. The reconstruction of the 3D fracture surface and quantitative roughness analysis of fracture surfaces extracted by different methods enables comparison of sensitivity (or robustness) to noise between each method.


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.


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.


Sign in / Sign up

Export Citation Format

Share Document