The long‐term temporal and spatial changes in statistical, source, and stress characteristics of one cluster of induced seismicity recorded at The Geysers geothermal field (U.S.) are analyzed in relation to the field operations, fluid migration, and constraints on the maximum likely magnitude. Two injection wells, Prati‐9 and Prati‐29, located in the northwestern part of the field and their associated seismicity composed of 1776 events recorded throughout a 7 year period were analyzed. The seismicity catalog was relocated, and the source characteristics including focal mechanisms and static source parameters were refined using first‐motion polarity, spectral fitting, and mesh spectral ratio analysis techniques. The source characteristics together with statistical parameters (b value) and cluster dynamics were used to investigate and understand the details of fluid migration scheme in the vicinity of injection wells. The observed temporal, spatial, and source characteristics were clearly attributed to fluid injection and fluid migration toward greater depths, involving increasing pore pressure in the reservoir. The seasonal changes of injection rates were found to directly impact the shape and spatial extent of the seismic cloud. A tendency of larger seismic events to occur closer to injection wells and a correlation between the spatial extent of the seismic cloud and source sizes of the largest events was observed suggesting geometrical constraints on the maximum likely magnitude and its correlation to the average injection rate and volume of fluids present in the reservoir.

The spatiotemporal, kinematic, and source characteristics of induced seismicity occurring at different fluid injection rates are investigated to determine the predominant physical mechanisms responsible for induced seismicity at the northwestern part of The Geysers geothermal field, California. We analyze a relocated hypocenter catalog from a seismicity cluster where significant variations of the stress tensor orientation were previously observed to correlate with injection rates. We find that these stress tensor orientation changes may be related to increased pore pressure and the corresponding changes in poroelastic stresses at reservoir depth. Seismic events during peak injections tend to occur at greater distances from the injection well, preferentially trending parallel to the maximum horizontal stress direction. In contrast, at lower injection rates the seismicity tends to align in a different direction which suggests the presence of a local fault. During peak injection intervals, the relative contribution of strike-slip faulting mechanisms increases. Furthermore, increases in fluid injection rates also coincide with a decrease in b values. Our observations suggest that regardless of the injection stage, most of the induced seismicity results from thermal fracturing of the reservoir rock. However, during peak injection intervals, the increase in pore pressure may likewise be responsible for the induced seismicity. By estimating the thermal and hydraulic diffusivities of the reservoir, we confirm that the characteristic diffusion length for pore pressure is much greater than the corresponding length scale for temperature and also more consistent with the spatial extent of seismicity observed during different injection rates.

The Northwest Geysers Enhanced Geothermal System (EGS) demonstration project aims to create an EGS by directly and systematically injecting cool water at relatively low pressure into a known High Temperature (280-400°C) Zone (HTZ) located under the conventional (240°C) geothermal steam reservoir at The Geysers geothermal field in California. In this paper, the results of coupled thermal, hydraulic, and mechanical (THM) analyses made using a model developed as part of the prestimulation phase of the EGS demonstration project is presented. The model simulations were conducted in order to investigate injection strategies and the resulting effects of cold-water injection upon the EGS system; in particular to predict the extent of the stimulation zone for a given injection schedule. The actual injection began on October 6, 2011, and in this paper a comparison of pre-stimulation model predictions with micro-earthquake (MEQ) monitoring data over the first few months of a one-year injection program is presented. The results show that, by using a calibrated THM model based on historic injection and MEQ data at a nearby well, the predicted extent of the stimulation zone (defined as a zone of high MEQ density around the injection well) compares well with observed seismicity. The modeling indicates that the MEQ events are related to shear reactivation of preexisting fractures, which is triggered by the combined effects of injection-induced cooling around the injection well and small changes in steam pressure as far as half a kilometer away from the injection well. Pressure-monitoring data at adjacent wells and satellite-based groundsurface deformation data were also used to validate and further calibrate reservoirscale hydraulic and mechanical model properties. The pressure signature monitored from the start of the injection was particularly useful for a precise back-calculation of

We performed a series of 3D thermo-hydro-mechanical (THM) simulations to study the influences of hydromechanical and thermal processes on the development of an enhanced geothermal system, strongly influenced by a network of short fault zones. The model we developed was calibrated by comparing the simulated THM responses to field observations, including ground-surface deformations, well pressure, and microseismic activity. Of particular importance was the comparison between the observed temporal and spatial distribution of microseismic activity, and the calculated shear reactivation of preexisting fractures inferred from simulated elasto-plastic mechanical responses in the short fault zone. Using this approach, we could identify when fault zones were reactivated (as manifested in the field by a surge of local microseismic activity within the fault zone), and we could back-calculate the in situ stress field as being close to the stress conditions required for shear reactivation. Our results show that the main mechanisms of inducing seismicity are related to injection-induced pressure increase and cooling. During injection, the reservoir expansion caused by the pressure increase led to mechanical stress transfer through the reservoir, which prevented or delayed the reactivation of preexisting fractures. After injection stopped, there was an inversion of the mechanical stress transfers that favored shear reactivation, which may explain why microseismic activity occurred after the cessation of the injection.

No abstract

a b s t r a c tAn Enhanced Geothermal System (EGS) Demonstration Project is currently underway in the Northwest Geysers. The project goal is to demonstrate the feasibility of stimulating a deep high-temperature reservoir (HTR) (up to 400 • C, 750 • F). Two previously abandoned wells, Prati State 31 (PS-31) and Prati 32 (P-32), were reopened and deepened to be used as an injection and production doublet to stimulate the HTR. The deepened portions of both wells have conductive temperature gradients of 10 • F/100 ft (182 • C/km), produce connate native fluids and magmatic gas, and the rocks were isotopically unexchanged by meteoric water. The ambient temperature meteoric water injected into these hot dry rocks has evidently created a permeability volume of several cubic kilometers as determined by seismic monitoring. Preliminary isotopic analyses of the injected and produced water indicate that 50-75% of the steam from the created EGS reservoir is injection-derived.

We present a three‐dimensional thermohydromechanical numerical study of the evolution and distribution of the stress tensor within the northwest part of The Geysers geothermal reservoir (in California), including a detailed study of the region around one injection well from 2003 to 2012. Initially, after imposing a normal faulting stress regime, we calculated local changes in the stress regime around injection wells. Our results were compared with previously published studies in which the stress state was inferred from inverting the focal plane mechanism of seismic events. Our main finding is that changes in stress tensor orientation are caused by injection‐induced progressive cooling of the reservoir, as well as by the seasonal variations in injection rate. Because of the gravity flow and cooling around a liquid zone formed by the injection, the vertical stress reduction is larger and propagates far below the injection well. At the same time, the horizontal stress increases, mostly because of stress redistribution below and above the cooling area. These two phenomena cause the rotation of the stress tensor and the appearance of a strike‐slip regime above, inside, and below the cooling area. The cooling and the associated rotation of the stress regime can play a significant role in the observed long‐term deepening of the microseismicity below active injection wells.

In this paper, we use the Seismicity‐Based Reservoir Characterization approach to study the spatiotemporal dynamics of an injection‐induced microseismic cloud, monitored during the stimulation of an enhanced geothermal system, and associated with the Northwest Geysers Enhanced Geothermal System (EGS) Demonstration project (California). We identified the development of a seismically quiet domain around the injection well surrounded by a seismically active domain. Then we compare these observations with the results of 3‐D Thermo‐Hydro‐Mechanical simulations of the EGS, which accounts for changes in permeability as a function of the effective normal stress and the plastic strain. The results of our modeling show that (1) the aseismic domain is caused by both the presence of the injected cold water and by thermal processes. These thermal processes cause a cooling‐stress reduction, which prevent shear reactivation and favors fracture opening by reducing effective normal stress and locally increasing the permeability. This process is accompanied by aseismic plastic shear strain. (2) In the seismic domain, microseismicity is caused by the reactivation of the preexisting fractures, resulting from an increase in injection‐induced pore pressure. Our modeling indicates that in this domain, permeability evolves according to the effective normal stress acting on the shear zones, whereas shearing of preexisting fractures may have a low impact on permeability. We attribute this lack of permeability gain to the fact that the initial permeabilities of these preexisting fractures are already high (up to 2 orders of magnitude higher than the host rock) and may already be fully dilated by past tectonic straining.