Geology Reports⌕ Search

USGS · 70189623

Numerical modeling of injection, stress and permeability enhancement during shear stimulation at the Desert Peak Enhanced Geothermal System

Abstract

Creation of an Enhanced Geothermal System relies on stimulation of fracture permeability through self-propping shear failure that creates a complex fracture network with high surface area for efficient heat transfer. In 2010, shear stimulation was carried out in well 27-15 at Desert Peak geothermal field, Nevada, by injecting cold water at pressure less than the minimum principal stress. An order-of-magnitude improvement in well injectivity was recorded. Here, we describe a numerical model that accounts for injection-induced stress changes and permeability enhancement during this stimulation. In a two-part study, we use the coupled thermo-hydrological-mechanical simulator FEHM to: (i) construct a wellbore model for non-steady bottom-hole temperature and pressure conditions during the injection, and (ii) apply these pressures and temperatures as a source term in a numerical model of the stimulation. In this model, a Mohr-Coulomb failure criterion and empirical fracture permeability is developed to describe permeability evolution of the fractured rock. The numerical model is calibrated using laboratory measurements of material properties on representative core samples and wellhead records of injection pressure and mass flow during the shear stimulation. The model captures both the absence of stimulation at low wellhead pressure (WHP ≤1.7 and ≤2.4 MPa) as well as the timing and magnitude of injectivity rise at medium WHP (3.1 MPa). Results indicate that thermoelastic effects near the wellbore and the associated non-local stresses further from the well combine to propagate a failure front away from the injection well. Elevated WHP promotes failure, increases the injection rate, and cools the wellbore; however, as the overpressure drops off with distance, thermal and non-local stresses play an ongoing role in promoting shear failure at increasing distance from the well.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 39.884450178234395° to 40.6056120582602° latitude; -118.68530273437499° to -117.56469726562499° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

David Dempsey, Sharad Kelkar, Nick Davatzes, Stephen H. Hickman, Daniel Moos. 2015. Numerical modeling of injection, stress and permeability enhancement during shear stimulation at the Desert Peak Enhanced Geothermal System. https://doi.org/10.1016/j.ijrmms.2015.06.003

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Shear failure of a granite pin traversing a sawcut fault

Fault heterogeneities such as bumps, bends, and stepovers are commonly observed on natural faults, but are challenging to recreate under controlled laboratory conditions. We study deformation and microseismicity of a 76 mm-diameter Westerly granite cylinder with a sawcut fault with known frictional properties. An idealized asperity is added by emplacing a precision-ground 21 mm-diameter solid granite dowel that crosses the center of the fault at right angles. This intact granite ‘pin’ provides a strength contrast that resists fault slip. Upon loading to 80 MPa in a triaxial machine, we first observed a M -4 slip event that ruptured the sawcut fault, slipped 40 µm, but was halted by the granite pin. With continued loading, the pin failed in a swarm of thousands of M -6 to M -8 events known as acoustic emissions (AEs). Once the pin was fractured to a critical point, it permitted complete rupture events (M -3) on the sawcut fault (stick-slip instabilities). Subsequent slip events were preceded by clusters of foreshock-like AEs, all located on the fault plane, and the spatial extent of the foreshock clusters is consistent with our estimate of a critical nucleation dimension h*. We also identified an aseismic zone on the fault plane surrounding the fractured rock pin. A post-mortem analysis of the sample showed a thick gouge layer where the pin intersected the fault, suggesting that dilatancy of this gouge propped open the fault and prevented microseismic events in its vicinity. Recorded microseismicity separates into three categories: slip on the sawcut fault, fracture of the intact rock pin, and off-fault seismicity associated with pin-related rock joints. We found that pin fracture events were exclusively implosive (anticrack) even though the shear process zone was overall dilatant. This shows how aseismic effects can lead to unexpected seismic manifestations of certain faulting processes.

International Journal of Rock Mechanics and Mining↗

Sequential Gaussian co-simulation of rate decline parameters of longwall gob gas ventholes

Gob gas ventholes (GGVs) are used to control methane inflows into a longwall mining operation by capturing the gas within the overlying fractured strata before it enters the work environment. Using geostatistical co-simulation techniques, this paper maps the parameters of their rate decline behaviors across the study area, a longwall mine in the Northern Appalachian basin. Geostatistical gas-in-place (GIP) simulations were performed, using data from 64 exploration boreholes, and GIP data were mapped within the fractured zone of the study area. In addition, methane flowrates monitored from 10 GGVs were analyzed using decline curve analyses (DCA) techniques to determine parameters of decline rates. Surface elevation showed the most influence on methane production from GGVs and thus was used to investigate its relation with DCA parameters using correlation techniques on normal-scored data. Geostatistical analysis was pursued using sequential Gaussian co-simulation with surface elevation as the secondary variable and with DCA parameters as the primary variables. The primary DCA variables were effective percentage decline rate, rate at production start, rate at the beginning of forecast period, and production end duration. Co-simulation results were presented to visualize decline parameters at an area-wide scale. Wells located at lower elevations, i.e., at the bottom of valleys, tend to perform better in terms of their rate declines compared to those at higher elevations. These results were used to calculate drainage radii of GGVs using GIP realizations. The calculated drainage radii are close to ones predicted by pressure transient tests.

Pennsylvania↗

An analytical solution for transient flow of Bingham viscoplastic materials in rock fractures

We present below an analytical solution to model the one-dimensional transient flow of a Bingham viscoplastic material in a fracture with parallel walls (smooth or rough) that is subjected to an applied pressure gradient. The solution models the acceleration and the deceleration of the material as the pressure gradient changes with time. Two cases are considered: A pressure gradient applied over a finite time interval and an applied pressure gradient that is constant over time. The solution is expressed in dimensionless form and can therefore be used for a wide range of Bingham viscoplastic materials. The solution is also capable of capturing the transition that takes place in a fracture between viscoplastic flow and rigid plug flow. Also, it shows the development of a rigid central layer in fractures, the extent of which depends on the fluid properties (viscosity and yield stress), the magnitude of the pressure gradient, and the fracture aperture and surface roughness. Finally, it is shown that when a pressure gradient is applied and kept constant, the solution for the fracture flow rate converges over time to a steady-state solution that can be defined as a modified cubic law. In this case, the fracture transmissivity is found to be a non-linear function of the head gradient. This solution provides a tool for a better understanding of the flow of Bingham materials in rock fractures, interfaces, and cracks. ?? 2001 Elsevier Science Ltd. All rights reserved.

International Journal of Rock Mechanics and Mining↗