Search NASA⌕ Search

SEARCH · Search NASA

Results for “SEISMICITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Local-Distance Seismic Event Relocation and Relative Magnitude Estimation, Applications to Mining Related Seismicity in the Powder River Basin, Wyoming

Recent efforts to characterize small ( M w < 3 ) seismic events at local distances have become more important because of the increased observation of human-triggered and induced seismicity and the need to advance nuclear explosion monitoring capabilities. The signals generated by low-magnitude seismic sources necessitate the use of nearby short-period observations, which are sensitive to local geological heterogeneity. Local to near-regional distance ( < 300 km ) surface and shear waves can dominate short-period observations from small, shallow seismic sources. In this work, we utilize these observations to estimate precise, relative locations and magnitudes of ~ 700 industrial mining events in Wyoming, using nearly 360,000 observations. The precise, relative location estimates (with formal location uncertainty estimates of less than 1 km) collapse a diffuse collection of mining events into discrete clusters associated with individual blasting operations. We also invert the cross-correlation amplitudes to estimate precise, relative moment magnitude estimates, which help validate and identify disparities in the event sizes reported by regional network catalogs. Joint use of multiple phases allows for the inclusion of more seismic events due to the increase in the number of observations. In some cases, using a single phase allowed us to relocate only 50% of the original reported seismic events within a cluster. Combining shear- and surface-wave phases increased the number of events to above 90% of the original events, allowing us to characterize a broader range of event sizes, source to station distances, and event distributions. This analysis takes a step toward making a fuller characterization of small industrial seismic events observed at local distances.

58 GEOSCIENCES↗

Assessing the need for an update of a probabilistic seismic hazard analysis using a SSHAC Level 1 study and the Seismic Hazard Periodic Reevaluation Methodology

A key decision for nuclear facilities is evaluating the need for an update of an existing seismic hazard analysis in light of new data and information that has become available since the time that the analysis was completed. We introduce the newly developed risk-informed Seismic Hazard Periodic Review Methodology (referred to as the SHPRM) and present how a Senior Seismic Hazard Analysis Committee (SSHAC) Level 1 probabilistic seismic hazard analysis (PSHA) was performed in an implementation of this new methodology. The SHPRM offers a defensible and documented approach that considers both the changes in seismic hazard and engineering-based risk information of an existing nuclear facility to assess the need for an update of an existing PSHA. The SHPRM has seven evaluation criteria that are employed at specific analysis, decision, and comparison points which are applied to seismic design categories established for nuclear facilities in United States. The SHPRM is implemented using a SSHAC Level 1 study performed for the Idaho National Laboratory, USA. The implementation focuses on the first six of the seven evaluation criteria of the SHPRM which are all provided from the SSHAC Level 1 PSHA. Finally, to illustrate outcomes of the SHPRM that do not lead to the need for an update and those that do, the example implementations of the SHPRM are performed for nuclear facilities that have target performance goals expressed as the mean annual frequency of unacceptable performance at 1x10 -4 , 4x10 -5 and 1x10 -5 .

58 GEOSCIENCES↗

Real-Time Seismic System for Monitoring, Imaging, and Characterization (RT-SEISMIC)

The goal of this Laboratory Directed Research and Development (LDRD) project was to develop a borehole seismic source and sensor array to enable real-time seismic imaging at scales and conditions relevant to the energy industry including both fossil-energy and geothermal. In FY21 and FY22, we designed, built, and tested both a prototype impulse source module for generating seismic energy and a sensing module for recording ground motions generated by the source module array. A pneumatically driven vibratory source was also designed. The source modules were fabricated with all high temperature components and the team has worked to incorporate the current RT-SEISMIC electronics design into a commercially available, high temperature silicon-on-insulator chip integrated circuit. Several issues were identified during fabrication and lab testing that led to redesign of several system components and subsequent retesting. The final round of testing showed that while metal/graphite-based seals worked quite well for static seals, they were unable to provide an adequate gas seal for dynamic, reciprocating part movements which necessitated a final redesign using Kalrez. This change will result in a continuous temperature rating of approximately 275 degrees C for the system. While a field test of the RT-SEISMIC system was targeted in FY22, due to the extended lab testing and redesign efforts, field testing was not achieved. As a result of this LDRD investment, several sponsors have expressed interest in RT-SEISMIC and we expect to continue towards a field demonstration of the full system in the future.

58 GEOSCIENCES↗

Seismicity Monitoring in North-Central New Mexico by the Los Alamos Seismic Network

The initial seismic stations of the Los Alamos Seismic Network (LASN) were installed in 1973, as a part of research on monitoring of nuclear testing. The extent of the network rapidly expanded by the late 1970s. By the middle 1980s, spatial coverage of the network was drastically reduced, due to a loss of funding. Since then, however, it has been possible to expand the coverage of the network with additional stations and to slowly upgrade network equipment to match contemporary instrument standards. These improvements will make it possible to keep recording the network’s data and to initiate real-time exchange of data with other institutions to improve earthquake monitoring throughout New Mexico and neighboring states. During more than 40 yr of operation, the network has provided a slow but steady increase in the volume of earthquake seismograms available to study seismicity and tectonics in north-central New Mexico. The current network covers an area from the Valles Caldera on the west, to within the Rio Grande rift on the east. LASN has yielded locations for about 900 earthquakes in north-central New Mexico between 1973 and 2013 (the most recent locations available). Epicenters of these have a complex pattern, with some that can be attributed to the deformation of the rift, though most are spread in areas west of the rift. A lack of seismicity in and near the Valles Caldera reinforces an earlier observation, based on far fewer earthquakes, which first called attention to this paucity of seismicity and attributed it to ductile deformation resulting from elevated crustal temperatures.

58 GEOSCIENCES↗

Mechanism of co-seismic water level change following four great earthquakes – insights from co-seismic responses throughout the Chinese mainland

In this study, we analyze the co-seismic groundwater level responses to four great earthquakes recorded by China's network of groundwater monitoring wells. The large number of operational wells (164 wells for the 2007 Mw 8.5 Sumatra earthquake, 245 wells for the Mw 7.9 Wenchuan earthquake, 228 wells for the Mw 9.0 Tohoku earthquake and 223 wells for 2012 Mw 8.6 Sumatra earthquake) and co-seismic responses provide an opportunity to test hypotheses on mechanisms for co-seismic water level changes. Overall, the co-seismic water level responses are complex over large spatial scales, and there is great variability both in the sign and amplitude of water level responses in the data set. As shown in previous studies, permeability change, rather than static strain, is a more plausible mechanism to explain most of the co-seismic responses. However, we find through tidal analysis of water level responses to solid Earth tide that only one third of these wells that showed a sustained post-seismic response can be explained by earthquake-induced permeability change in aquifers, and these wells had sustained (>30 days) water level changes. Lastly, wells that did not show sustained changes are more likely affected by permeability changes only immediately adjacent to the wellbore.

58 GEOSCIENCES↗

Deep compressed seismic learning for fast location and moment tensor inferences with natural and induced seismicity

Fast detection and characterization of seismic sources is crucial for decision-making and warning systems that monitor natural and induced seismicity. However, besides the laying out of ever denser monitoring networks of seismic instruments, the incorporation of new sensor technologies such as Distributed Acoustic Sensing (DAS) further challenges our processing capabilities to deliver short turnaround answers from seismic monitoring. In response, this work describes a methodology for the learning of the seismological parameters: location and moment tensor from compressed seismic records. In this method, data dimensionality is reduced by applying a general encoding protocol derived from the principles of compressive sensing. The data in compressed form is then fed directly to a convolutional neural network that outputs fast predictions of the seismic source parameters. Thus, the proposed methodology can not only expedite data transmission from the field to the processing center, but also remove the decompression overhead that would be required for the application of traditional processing methods. An autoencoder is also explored as an equivalent alternative to perform the same job. We observe that the CS-based compression requires only a fraction of the computing power, time, data and expertise required to design and train an autoencoder to perform the same task. Implementation of the CS-method with a continuous flow of data together with generalization of the principles to other applications such as classification are also discussed.

54 ENVIRONMENTAL SCIENCES↗

Isotropic seismic radiation from rock damage and dilatancy

SUMMARY We present semi-analytical and numerical results on isotropic radiation expected to be generated by seismic sources with damaged dilatant materials that sustain during failure an overall shear deformation. The results are derived for inclusions with different stiffness components, including those responsible for shear-volumetric coupling in dilatant material, and different geometries, in a surrounding elastic solid. The obtained relations indicate that the isotropic seismic potency and seismic moment components are proportional to the shear potency components multiplied by the aspect ratio between the thickness and length of the source zone. For a typical rupture zone with a thin tabular geometry, the expected small aspect ratio explains the general seismological observations of small isotropic radiation. However, isotropic radiation induced by dynamic dilation for short and thick sources is more pronounced and may reach for spherical sources sustaining shear deformation 18 per cent of the shear components. Sources that include also tensile components and explosions can produce larger isotropic radiation. The results provide rigorous relations that can be used to estimate different source properties in numerical simulations, experiments and field studies.

Geochemistry & Geophysics↗

The effects of earth model uncertainty on the inversion of seismic data for seismic source functions

SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.

58 GEOSCIENCES↗

Detection of random noise and anatomy of continuous seismic waveforms in dense array data near Anza California

We develop a methodology to separate continuous seismic waveforms into random noise (RN), not random noise (NRN) produced by earthquakes, wind, traffic and other sources of ground motions, and an undetermined mixture of signals. The analysis is applied to continuous records from a dense seismic array on the San Jacinto fault zone. To detect RN signals, we cut hourly waveforms into non-overlapping 1 s time windows and apply cross-correlations to separate RN candidates from outliers. The cross-correlation coefficients between different RN candidates fall into a tight range (i.e. 0.09–0.35), while cross-correlation coefficients of RN candidates with NRN signals (e.g. seismic or air-traffic events) are lower. The amplitude spectra of RN candidates have a well-defined level, while the amplitude spectra of other signals deviate from that level. Using these properties, we examine the amplitude spectra of moving time windows and cross-correlation coefficients with RN templates in each hour. The hourly RN is quasi-stationary and the results cluster tightly in the parameter space of cross-correlation coefficients and L2 norm deviations from the mean spectra of RN candidates. Time windows with parameters in this tight cluster are identified as RN, windows that deviate significantly from the RN cluster are identified as NRN and windows with values in between are identified as mixed signals. Several iterations on each hourly data are used to update and stabilize the selection of RN templates and mean noise spectra. For the days examined, the relative fractions of RN, NRN and mixed signals in local day (night) times are about 26 (42 percent), 40 (33 percent) and 34 percent (25 percent), respectively.

Geochemistry & Geophysics↗

Enhancing trap and fault seal analyses by integrating observations from HR3D seismic data with well logs and conventional 3D seismic data, Texas inner shelf

An understanding of trap and fault seal quality is critical for assessing hydrocarbon prospectivity. To achieve this, modern analytical techniques leverage well data and conventional industry-standard 3D seismic data to evaluate the trap, and any faults displacing the reservoir and top seal intervals. Above all, geological interpretation provides the framework of trap and fault seal analyses, but can be hindered by the data resolution, quality and acquisition style of the conventional seismic data. Furthermore, limiting the analysis to only the petroleum system at depth may lead to erroneous perceptions because interpreting overburden features, such as shallow faults or gas chimneys, can provide valuable observations with respect to container performance, and can to help validate trap and fault seal predictions. Here, a supplement to conventional 3D data are high-resolution 3D seismic (HR3D) data, which provide detailed images of the overburden geology. This study utilizes an HR3D seismic volume in the San Luis Pass area of the Texas inner shelf, where shallow fault tips and a sizeable gas chimney are interpreted over an unsuccessful hydrocarbon prospect. Static post-drill fault seal and trap analyses suggest that the primary fault displacing the structural closure could have withheld columns of gas c. 100 m high, but disagree with our HR3D seismic interpretations and dry-well analyses. From our results, we hypothesize that tertiary gas migration through fault conduits reduced the hydrocarbon column in the prospective Early Miocene reservoir, and may have resulted from continued movement along the intersecting faults. Overall, this study reinforces the importance of understanding the overburden geology and geohistory of faulted prospects, and demonstrates the utility of pre-drill HR3D acquisition when conducting trap and fault seal analyses.

3D seismic↗

Multiple plates subducting beneath Colombia, as illuminated by seismicity and velocity from the joint inversion of seismic and gravity data

Subduction beneath the northernmost Andes in Colombia is complex. Based on seismicity distributions, multiple segments of slab appear to be subducting, and arc volcanism ceases north of 5° N. Here, we illuminate the subduction system through hypocentral relocations and Vp and Vs models resulting from the joint inversion of local body wave arrivals, surface wave dispersion measurements, and gravity data. The simultaneous use of multiple data types takes advantage of the differing sensitivities of each data type, resulting in velocity models that have improved resolution at both shallower and deeper depths than would result from traditional travel time tomography alone. The relocated earthquake dataset and velocity model clearly indicate a tear in the Nazca slab at 5° N, corresponding to a 250-km shift in slab seismicity and the termination of arc volcanism. North of this tear, the slab is flat, and it comprises slabs of two sources: the Nazca and Caribbean plates. The Bucaramanga nest, a small region of among the most intense intermediate-depth seismicity globally, is associated with the boundary between these two plates and possibly with a zone of melting or elevated water content, based on reduced Vp and increased Vp/Vs. As a result, we also use relocated seismicity to identify two new faults in the South American plate, one related to plate convergence and one highlighted by induced seismicity.

58 GEOSCIENCES↗

Integrating distributed acoustic sensing, borehole 3C geophone array, and surface seismic array data to identify long-period long-duration seismic events during stimulation of a Marcellus Shale gas reservoir

Microseismic monitoring by downhole geophones, surface seismic, fiber-optic distributed acoustic sensing (DAS), and distributed temperature sensing (DTS) observations were made during the hydraulic fracture stimulation of the MIP-3H well in the Marcellus Shale in northern West Virginia. DAS and DTS data measure the fiber strain and temperature, respectively, along a fiber-optic cable cemented behind the casing of the well. The presence of long-period long-duration (LPLD) events is evaluated in the borehole geophones, DAS data, and surface seismic data of one of the MIP-3H stimulated stages. LPLD events are generally overlooked during the conventional processing of microseismic data, but they represent significant nonbrittle deformation produced during hydraulic fracture stimulation. In a single stage that was examined, 160 preexisting fractures and two faults of suboptimal orientation are noted in the image logs. We identified two low-frequency ([Formula: see text]) events of large temporal duration (tens of seconds) by comparing the surface seismic data, borehole geophone data, and DAS amplitude spectra of one of the MIP-3H stages. Spectrograms of DAS traces in time and depth reveal that the first low-frequency event might be an injection noise that has footprints on all DAS channels above the stimulated stage. However, the surface seismic array indicates an LPLD event concurrent with the first low-frequency event on DAS. The second LPLD event on DAS data and surface seismic data is related to a local deformation and does not have footprints on all DAS channels. The interpreted events have duration less than 100 s with frequencies concentrated below 10 Hz, and are accompanied by microseismic events.

Geochemistry & Geophysics↗

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗

Apollo 14 and 16 Active Seismic Experiments, and Apollo 17 Lunar Seismic Profiling

Seismic refraction experiments were conducted on the moon by Apollo astronauts during missions 14, 16, and 17. Seismic velocities of 104, 108, 92, 114 and 100 m/sec were inferred for the lunar regolith at the Apollo 12, 14, 15, 16, and 17 landing sites, respectively. These data indicate that fragmentation and comminution caused by meteoroid impacts has produced a layer of remarkably uniform seismic properties moonwide. Brecciation and high porosity are the probable causes of the very low velocities observed in the lunar regolith. Apollo 17 seismic data revealed that the seismic velocity increases very rapidly with depth to 4.7 km/sec at a depth of 1.4 km. Such a large velocity change is suggestive of compositional and textural changes and is compatible with a model of fractured basaltic flows overlying anorthositic breccias. 'Thermal' moonquakes were also detected at the Apollo 17 site, becoming increasingly frequent after sunrise and reaching a maximum at sunset. The source of these quakes could possibly be landsliding.

Source record↗

Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

Here we perform numerical simulations of sequences of earthquake and aseismic slip on planar rate and state faults separating dissimilar material within the 2-D plane strain approximation. We resolve all stages of the earthquake cycle from aseismic slip to fast ruptures while incorporating full inertia effects during seismic event propagation. We show that bimaterial coupling results in favorable nucleation site and subsequent asymmetric rupture propagation. We demonstrate that increasing the material contrast enhances this asymmetry leading to higher slip rates and normal stress drops in the preferred rupture propagation direction. The normal stress drop, induced by the bimaterial effect, leads to strong dynamic weakening of the fault and may destabilize the creeping region on a heterogeneous rate and state fault, resulting in extended rupture propagation. Such rupture penetration into creeping patches may lead to more frequent opening of earthquake gates, causing increased seismic hazard. Furthermore, bimaterial coupling may lead to irregular seismicity pattern in terms of event length, peak slip rates, and hypocenter location, depending on the properties of the creeping patches bordering the seismogenically active part of the fault. Our results highlight robust characteristics of bimaterial interfaces that persist over long sequence of events and suggest the need for further exploration of the role of material contrast in earthquake physics and models of seismic hazard.

58 GEOSCIENCES↗