Search USGS⌕ Search

Geology topics

Nicholas C. Davatzes

Publications and source records attributed to Nicholas C. Davatzes.

11 recordsLinked to original sources
USGS2017

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Compartmentalization of the Coso East Flank geothermal field imaged by 3-D full-tensor MT inversion

Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2–5 Ohm m) and appreciable lateral extent (>1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Z xx and Z yy data components leads to a robust east-dip (60†) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.

Coso Range↗
USGS2016

Analysis and interpretation of stress indicators in deviated wells of the Coso Geothermal Field

Characterizing the tectonic stress field is an integral part of the development of hydrothermal systems and especially for enhanced geothermal systems (EGS). With a well characterized stress field the propensity of fault slip on faults with known location and orientation can be identified. Faults that are critically oriented for faulting with respect to the stress field are known to provide natural fluid pathways. A high slip tendency makes a fault a likely candidate for reactivation during the creation of an EGS. Similarly, the stress state provides insight for the potential of larger, damaging earthquakes should extensive portions of well-oriented, larger faults be reactivated. The analysis of stress indicators such as drilling-induced fractures and borehole breakouts is the main tool to infer information on the stress state of a geothermal reservoir. The standard procedure is applicable to sub-vertical wellbore sections and highly deviated sections have to be discarded. However, in order to save costs and reduce the environmental impact most recent wells are directionally drilled with deviations that require appropriate consideration of the deviated trajectory. Here we present an analysis scheme applicable to arbitrary well trajectories or a combination of wells to infer the stress state. Through the sampling of the stress tensor along several directions additional information on the stress regime and even relative stress magnitudes can be obtained. We apply this method on image logs from the pair of wells 58-10 and 58A-10 that were drilled from the same well pad. Both wells have image logs of about 2km of their trajectories that are separated by less than 300m. For both wells we obtain a mean orientation of S Hmax of N23° with large standard deviations of locations of stress indicators of 24° and 26°, respectively. While the local stress direction is highly variable along both wells with dominant wavelengths from around 50 to 500m, the mean directions are very consistent and also agree with previous stress estimates in the eastern part of the Coso Geothermal Field. In order to obtain a reliable estimation of the stress orientation in this setting, it is necessary to sample the stress field on an interval long to capture several of the dominant wavelengths.

California↗
USGS2015

Differentiating induced and natural seismicity using space-time-magnitude statistics applied to the Coso Geothermal field

A remarkable characteristic of earthquakes is their clustering in time and space, displaying their self-similarity. It remains to be tested if natural and induced earthquakes share the same behavior. We study natural and induced earthquakes comparatively in the same tectonic setting at the Coso Geothermal Field. Covering the preproduction and coproduction periods from 1981 to 2013, we analyze interevent times, spatial dimension, and frequency-size distributions for natural and induced earthquakes. Individually, these distributions are statistically indistinguishable. Determining the distribution of nearest neighbor distances in a combined space-time-magnitude metric, lets us identify clear differences between both kinds of seismicity. Compared to natural earthquakes, induced earthquakes feature a larger population of background seismicity and nearest neighbors at large magnitude rescaled times and small magnitude rescaled distances. Local stress perturbations induced by field operations appear to be strong enough to drive local faults through several seismic cycles and reactivate them after time periods on the order of a year.

California↗
USGS2015

Natural or Induced: Identifying Natural and Induced Swarms from Pre-production and Co-production Microseismic Catalogs at the Coso Geothermal Field

Increased levels of seismicity coinciding with injection of reservoir fluids have prompted interest in methods to distinguish induced from natural seismicity. Discrimination between induced and natural seismicity is especially difficult in areas that have high levels of natural seismicity, such as the geothermal fields at the Salton Sea and Coso, both in California. Both areas show swarm-like sequences that could be related to natural, deep fluid migration as part of the natural hydrothermal system. Therefore, swarms often have spatio-temporal patterns that resemble fluid-induced seismicity, and might possibly share other characteristics. The Coso Geothermal Field and its surroundings is one of the most seismically active areas in California with a large proportion of its activity occurring as seismic swarms. Here we analyze clustered seismicity in and surrounding the currently produced reservoir comparatively for pre-production and co-production periods. We perform a cluster analysis, based on the inter-event distance in a space-time-energy domain to identify notable earthquake sequences. For each event j, the closest previous event i is identified and their relationship categorized. If this nearest neighbor’s distance is below a threshold based on the local minimum of the bimodal distribution of nearest neighbor distances, then the event j is included in the cluster as a child to this parent event i. If it is above the threshold, event j begins a new cluster. This process identifies subsets of events whose nearest neighbor distances and relative timing qualify as a cluster as well as a characterizing the parent-child relationships among events in the cluster. We apply this method to three different catalogs: (1) a two-year microseismic survey of the Coso geothermal area that was acquired before exploration drilling in the area began; (2) the HYS_catalog_2013 that contains 52,000 double-difference relocated events and covers the years 1981 to 2013; and (3) a catalog of 57,000 events with absolute locations from the local network recorded between 2002 and 2007. Using this method we identify 10 clusters of more than 20 events each in the pre-production survey and more than 200 distinct seismicity clusters that each contain at least 20 and up to more than 1000 earthquakes in the more extensive catalogs. The cluster identification method used yields a hierarchy of links between multiple generations of parent and offspring events. We analyze different topological parameters of this hierarchy to better characterize and thus differentiate natural swarms from induced clustered seismicity and also to identify aftershock sequences of notable mainshocks. We find that the branching characteristic given by the average number of child events per parent event is significantly different for clusters below than for clusters around the produced field.

Conference Paper↗
USGS2014

Micro-seismicity and seismic moment release within the Coso Geothermal Field, California

We relocate 16 years of seismicity in the Coso Geothermal Field (CGF) using differential travel times and simultaneously invert for seismic velocities to improve our knowledge of the subsurface geologic and hydrologic structure. We expand on our previous results by doubling the number of relocated events from April 1996 through May 2012 using a new field-wide 3-D velocity model. Relocated micro-seismicity sharpens in many portions of the active geothermal reservoir, likely defining large-scale fault zones and fluid pressure compartment boundaries. However, a significant fraction of seismicity remains diffuse and does not cluster into sharply defined structures, suggesting that permeability is maintained within the reservoir through distributed brittle failure. The seismic velocity structure reveals heterogeneous distributions of compressional (Vp) and shear (Vs) wave speed, with Vs generally higher in the Main Field and East Flank and Vp remaining relatively uniform across the CGF, but with significant local variations. The Vp/Vs ratio appears to outline the two main producing compartments of the reservoir at depths below mean ground level of approximately 1 to 2.5 km, with a ridge of relatively high Vp/Vs separating the Main Field from the East Flank. Detailed analyses of spatial and temporal variations in earthquake relocations and cumulative seismic moment release in the East Flank reveal three regions with persistently high rates of seismic activity. Two of these regions exhibit sharp, stationary boundaries at the margins of the East Flank that likely represent barriers to fluid flow and advective heat transport. However, seismicity and moment release in a third region at the northern end of the East Flank spread over time to form an elongated NE to SW structure, roughly parallel both to an elongated cluster of seismicity at the southern end of the East Flank and to regional fault traces mapped at the surface. Our results indicate that high-precision relocations of micro-seismicity and simultaneous velocity inversions in conjunction with mapping of seismic moment release can provide useful insights into subsurface structural features and hydrologic compartmentalization within the Coso Geothermal Field.

California↗
USGS2012

Linear complementarity formulation for 3D frictional sliding problems

Frictional sliding on quasi-statically deforming faults and fractures can be modeled efficiently using a linear complementarity formulation. We review the formulation in two dimensions and expand the formulation to three-dimensional problems including problems of orthotropic friction. This formulation accurately reproduces analytical solutions to static Coulomb friction sliding problems. The formulation accounts for opening displacements that can occur near regions of non-planarity even under large confining pressures. Such problems are difficult to solve owing to the coupling of relative displacements and tractions; thus, many geomechanical problems tend to neglect these effects. Simple test cases highlight the importance of including friction and allowing for opening when solving quasi-static fault mechanics models. These results also underscore the importance of considering the effects of non-planarity in modeling processes associated with crustal faulting.

Computational Geosciences↗
USGS2012

Using micro-seismicity and seismic velocities to map subsurface geologic and hydrologic structure within the Coso geothermal field, California

Geothermal reservoirs derive their capacity for fluid and heat transport in large part from faults and fractures. Micro-seismicity generated on such faults and fractures can be used to map larger fault structures as well as secondary fractures that add access to hot rock, fluid storage and recharge capacity necessary to have a sustainable geothermal resource. Additionally, inversion of seismic velocities from micro-seismicity permits imaging of regions subject to the combined effects of fracture density, fluid pressure and steam content, among other factors. We relocate 14 years of seismicity (1996-2009) in the Coso geothermal field using differential travel times and simultaneously invert for seismic velocities to improve our knowledge of the subsurface geologic and hydrologic structure. We utilize over 60,000 micro-seismic events using waveform cross-correlation to augment to expansive catalog of P- and S-wave differential travel times recorded at Coso. We further carry out rigorous uncertainty estimation and find that our results are precise to within 10s of meters of relative location error. We find that relocated micro-seismicity outlines prominent, through-going faults in the reservoir in some cases. We also find that a significant portion of seismicity remains diffuse and does not cluster into more sharply defined major structures. The seismic velocity structure reveals heterogeneous distributions of compressional (Vp) and shear (Vs) wave speed, with Vp generally lower in the main field when compared to the east flank and Vs varying more significantly in the shallow portions of the reservoir. The Vp/Vs ratio appears to outline the two main compartments of the reservoir at depths of -0.5 to 1.5 km (relative to sea-level), with a ridge of relatively high Vp/Vs separating the main field from the east flank. In the deeper portion of the reservoir this ridge is less prominent. Our results indicate that high-precision relocations of micro-seismicity can provide useful insight into: 1) prominent structural features, faults and fractures that contribute to the flow of fluid and heat in the reservoir; 2) diffuse seismicity throughout the reservoir representing fractures that likely contribute to the overall permeability, storage and heat exchange capacity of the reservoir, but which are not confined to prominent faults; and 3) seismic velocities that outline the major hydrologic compartments within the Coso geothermal field.

Coso geothermal field↗
USGS2010

Stress, fracture, and fluid-flow analysis using acoustic and electrical image logs in hot fractured granites of the Coso geothermal field, California, U.S.A.

Acoustic and electrical image logs in fractured granitic rocks penetrated by U.S. Navy well 58A-10, Coso Wash, in the eastern margin of the Coso geothermal field, California, were compared to evaluate their relative ability to characterize fractures and fault rock textures and to measure stress orientations from borehole failure. Electrical image logs are sensitive to variations in mineralogy or porosity, which affect conductivity. Thus, they capture both open and healed natural fractures as well as rock foliation.

California↗
USGS2009

Fractures, stress and fluid flow prior to stimulation of well 27-15, Desert Peak, Nevada, EGS project

A suite of geophysical logs has been acquired for structural, fluid flow and stress analysis of well 27-15 in the Desert Peak Geothermal Field, Nevada, in preparation for stimulation and development of an Enhanced Geothermal System (EGS). Advanced Logic Technologies Borehole Televiewer (BHTV) and Schlumberger Formation MicroScanner (FMS) image logs reveal extensive drilling-induced tensile fractures, showing that the current minimum compressive horizontal stress, Shmin, in the vicinity of well 27-15 is oriented along an azimuth of 114±17°. This orientation is consistent with the dip direction of recently active normal faults mapped at the surface and with extensive sets of fractures and some formation boundaries seen in the BHTV and FMS logs. Temperature and spinner flowmeter surveys reveal several minor flowing fractures that are well oriented for normal slip, although over-all permeability in the well is quite low. These results indicate that well 27-15 is a viable candidate for EGS stimulation and complements research by other investigators including cuttings analysis, a reflection seismic survey, pressure transient and tracer testing, and micro-seismic monitoring.

Nevada↗