Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic data”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

On the time-varying properties of the lunar seismic meteoroid population

Strong short-term fluctuations of meteroid impact rates are evident in a list of 1557 events derived from Apollo lunar seismic data. Times of fall and seismic signal amplitudes are considered in relation to the possible orbits and identification of the impacting objects. The entire lunar surface is the effective collector and the Apollo network data, gathered between 1970 and 1977, reflect the recognized abundance of large meteorite falls from early April through July. There is also some indication that briefly increased counts represent fragments of about 100 g or greater belonging to several meteor streams: Quadrantids, Aquarids, Perseids, Orionids, Leonids, Geminids, and possibly others as well.

Dorman, J.↗
DOE OSTI2020

The Sensitivity of Joint Inversions of Seismic and Geodynamic Data to Mantle Viscosity

Seismic tomography has revealed the existence of large-scale velocity heterogeneities in the mantle. The interpretation of seismic velocity anomalies in terms of temperature and chemical composition is nonunique. We use geodynamic observations including gravity, plate motions, dynamic topography, and excess ellipticity of the core-mantle boundary combined with seismic observations to investigate the thermo-chemical structure of the mantle through joint inversions. An outstanding issue, however, is the physical connection between mantle density anomalies and the surface geodynamic observations, which requires knowledge of the mantle viscosity structure. Here we perform joint inversions assuming different viscosity profiles and examine the dependence of the results on the viscosity. We first assume that mantle heterogeneity is due to thermal variations, which places a constraint on the relation between seismic velocity and density, and we subsequently relax the constraint to allow for potential nonthermal effects. In all of our joint inversions, a nonthermal origin of density anomalies is required to explain the geodynamic data, though the amount varies with the assumed viscosity structure. A common observation is a high-density chemical signal in the center of the large low-shear-velocity provinces at the base of the mantle resulting in a near neutral or slightly dense overall buoyancy there. Using the derived density models and their corresponding viscosity profiles, we also calculate instantaneous mantle flow fields. The predicted flow fields derived from joint inversions are generally similar but are quite different from flow fields using density models derived from a posteriori scaling of pure seismic tomography models.

58 GEOSCIENCES↗
NASA NTRS1973

Digital interactive image analysis by array processing

An attempt is made to draw a parallel between the existing geophysical data processing service industries and the emerging earth resources data support requirements. The relationship of seismic data analysis to ERTS data analysis is natural because in either case data is digitally recorded in the same format, resulting from remotely sensed energy which has been reflected, attenuated, shifted and degraded on its path from the source to the receiver. In the seismic case the energy is acoustic, ranging in frequencies from 10 to 75 cps, for which the lithosphere appears semi-transparent. In earth survey remote sensing through the atmosphere, visible and infrared frequency bands are being used. Yet the hardware and software required to process the magnetically recorded data from the two realms of inquiry are identical and similar, respectively. The resulting data products are similar.

Sabels, B. E.↗
DOE OSTI2020

Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. Here, we obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.

15 GEOTHERMAL ENERGY↗
NASA NTRS1974

Study of the characteristics of seismic signals generated by natural and cultural phenomena

Seismic data recorded at the Tonto Forest Seismological Observatory in Arizona and the Uinta Basin Seismological Observatory in Utah were used to compare the frequency of occurrence, severity, and spectral content of ground motions resulting from earthquakes, and other natural and man-made sources with the motions generated by sonic booms. A search of data recorded at the two observatories yielded a classification of over 180,000 earthquake phase arrivals on the basis of frequency of occurrence versus maximum ground velocity. The majority of the large ground velocities were produced by seismic surface waves from moderate to large earthquakes in the western United States, and particularly along the Pacific Coast of the United States and northern Mexico. A visual analysis of raw film seismogram data over a 3-year period indicates that local and regional seismic events, including quarry blasts, are frequent in occurrence, but do not produce ground motions at the observatories comparable to either the large western United States earthquakes or to sonic booms. Seismic data from the Nevada Test Site nuclear blasts were used to derive magnitude-distance-sonic boom overpressure relations.

Goforth, T. T.↗
DOE OSTI2023

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗
NASA NTRS1992

A test of present-day plate geometries for northeast Asia and Japan

Alternative geometries for the present-day configuration of plate boundaries in northeast Asia and Japan are tested using NUVEL-1 and 256 horizontal earthquake slip vectors from the Japan and northern Kuril trenches. Statistical analysis of the slip vectors is used to determine whether the North American, Eurasian, or Okhotsk plate overlies the trench. Along the northern Kuril trench, slip vectors are well-fit by the NUVEL-1 Pacific-North America Euler pole, but are poorly fit by the Pacific-Eurasia Euler pole. Results for the Japan trench are less conclusive, but suggest that much of Honshu and Hokkaido are also part of the North American plate. The simplest geometry consistent with the trench slip vectors is a geometry in which the North American plate extends south to 41 deg N, and possibly includes northern Honshu and southern Hokkaido. Although these results imply that the diffuse seismicity that connects the Lena River delta to Sakhalin Island and the eastern Sea of Japan records motion between Eurasia and North America, onshore geologic and seismic data define an additional belt of seismicity in Siberia that cannot be explained with this geometry. Assuming that these two seismic belts constitute evidence for an Okhotsk block, two published kinematic models for motion of the Okhotsk block are tested. The first model, which predicts motion of up to 15 mm/yr relative to North America, is rejected because Kuril and Japan trench slip vectors are fit more poorly than for the simpler geometry described above. The second model gives a good fit to the trench slip vectors, but only if Okhotsk-North America motion is slower than 5 mm/yr.

Demets, Charles↗
DOE OSTI2022

A Methodological Overview of Seismic Analysis for Nuclear Event Detection

Underground explosions generate potentially detectable signatures, including energy waves that travel through the Earth’s subsurface (i.e., seismic waves), low-frequency sound waves (i.e., infrasound and hydroacoustic waves), and radioactive gases and/or particles that might leak from the test cavity (if the event was nuclear). There can also be intelligence indicators of a test, such as observations of modified patterns of life and activity at a suspected test site. If all of these detectable signatures and intelligence indicators are present and self-consistent, then analysts have high confidence in classifying a signature generating event as an explosion. However, because only partial information about an event is likely to be available, determining whether an event was natural (e.g., an earthquake or landslide) or manmade (e.g., a chemical or nuclear explosion) is much more challenging. This primer describes how one category of event signatures—seismic signatures—can augment event analyses. While universities and government organizations have generated detailed technical descriptions of seismic analytic techniques, we seek to translate seismic event analysis for a broad, non-technical audience. When the geologic conditions near an event are well-characterized, seismic data can be used to calculate critical information, such as event location and depth, with relatively high accuracy. Moreover, specific features within seismic datasets can help determine whether an event was an explosion. However, a key challenge in seismic analysis is that geologic site conditions are often poorly characterized, complicating the ability to discern the true nature of the event. To overcome this challenge, geologists answer a series of questions (discussed in section 1) to guide seismic event analysis and determine the most probable nature of an event. As more information is gathered during each analytic step, confidence grows regarding the nature of the event. Section 2 addresses uncertainties in seismic analysis and the vital nature of high-fidelity geologic data for accurate seismic event analysis.

58 GEOSCIENCES↗
DOE OSTI2023

A Multimodal Event Catalog and Waveform Data Set That Supports Explosion Monitoring from Nevada, U.S.A.

Multimodal, curated data sets and nuisance event catalogs remain rare in the explosion monitoring community relative to curated seismic data sets. The source of this relative absence is the difficultly in deploying multimodal receivers that sense the seismic, acoustic, and other modalities from multiphysics sources. We provide such a data set in this study that delivers seismic, infrasound, and electromagnetic (magnetometer) sensor records collected over a two–week period, within 255 km of a 10 ton buried chemical explosion called DAG–4 that was located at 37.1146°, –116.0693° on 22 June 2019 21:06:19.88 UTC. This catalog includes 485 seismic, seismoacoustic, and infrasound–only events that an expert analyst manually built by reviewing waveforms from 29 seismic and infrasound sensors. Our data release includes waveforms from these 29 seismic, infrasound, and seismoacoustic stations and two magnetometer stations and their station metadata. We deliver these waveforms in NNSA KB Core CSS.w format (i4) with a corresponding wfdisc table that provides the header information. Here, we expect that this data set will provide a valuable, benchmark resource to develop signal processing algorithms and explosion monitoring methods against manual, human observations.

58 GEOSCIENCES↗
DOE OSTI2022

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗
DOE OSTI2025

TCCSP Site Characterization and Feasibility Report

The Tulare County Carbon Storage Project (TCCSP) focused on establishing the geological and commercial viability of a 50 million metric ton (MMT) carbon storage complex located in Tulare County, California. The geologic scope of work included collecting new and existing subsurface data, such as two-dimensional (2D) seismic data, advanced geophysical logs, core and fluid samples, and specialized core-analysis datasets, to support the integration and interpretation of findings from earlier tasks. These efforts contributed to the comprehensive feasibility study for commercial carbon capture and storage (CCS) development at the TCCSP site and in the surrounding region. The TCCSP project team evaluated the feasibility of the Lower Monterey Group (Santa Margarita Sand Member or its equivalent), Temblor Formation sand members (Olcese, Jewett, Vedder), and Lodo-Martinez Zone (Domengine Sandstone, Walker Formation, Lodo Formation, and the Martinez Sand) to serve as commercial CCS reservoirs along with potential for the Upper Monterey Group (Reef Ridge, Antelope Members) and the Kreyenhagen Formation to vertically seal underlying reservoirs from shallow above-zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Dense Seismic Array Study of a Legacy Underground Nuclear Test at the Nevada National Security Site

The complex postdetonation geologic structures that form after an underground nuclear explosion are difficult to constrain because increased heterogeneity around the damage zone affects seismic waves that propagate through the explosion site. Generally, a vertical rubble-filled structure known as a chimney is formed after an underground nuclear explosion that is composed of debris that falls into the subsurface cavity generated by the explosion. Compared with chimneys that collapse fully, leaving a surface crater, partially collapsed chimneys can have remnant subsurface cavities left in place above collapsed rubble. The 1964 nuclear test HADDOCK, conducted at the Nevada test site (now the Nevada National Security Site), formed a partially collapsed chimney with no surface crater. Understanding the subsurface structure of these features has significant national security applications, such as aiding the study of suspected underground nuclear explosions under a treaty verification. In this study, we investigated the subsurface architecture of the HADDOCK legacy nuclear test using hybrid 2D–3D active source seismic reflection and refraction data. The seismic data were acquired using 275 survey shots from the Seismic Hammer (a 13,000 kg weight drop) and 65 survey shots from a smaller accelerated weight drop, both recorded by ~1000 three-component 5 Hz geophones. First-arrival, P-wave tomographic modeling shows a low-velocity anomaly at ~200 m depth, likely an air-filled cavity caused by partial collapse of the rock column into the temporary postdetonation cavity. A high-velocity anomaly between 20 and 60 m depth represents spall-related compaction of the shallow alluvium. Hints of low velocities are also present near the burial depth (~364 m). The reflection seismic data show a prominent subhorizontal reflector at ~300 m depth, a short-curved reflector at ~200 m, and a high-amplitude reflector at ~50 m depth. Comparisons of the reflection sections to synthetic data and borehole stratigraphy suggest that these features correspond to the alluvium–tuff contact, the partial collapse cavity, and the spalled layer, respectively.

58 GEOSCIENCES↗
DOE OSTI2020

Methods for the robust computation of the long-period seismic spectrum of broad-band arrays

SUMMARY We describe array methods to search for low signal-to-noise ratio (SNR) signals in long-period seismic data using Fourier analysis. This is motivated by published results that find evidence of solar free oscillations in the Earth's seismic hum. Previous work used data from only one station. In this paper, we describe methods for computing spectra from array data. Arrays reduce noise level through averaging and provide redundancy that we use to distinguish coherent signal from a random background. We describe two algorithms for calculating a robust spectrum from seismic arrays, an algorithm that automatically removes impulsive transient signals from data, a jackknife method for estimating the variance of the spectrum, and a method for assessing the significance of an entire spectral band. We show examples of their application to data recorded by the Homestake Mine 3-D array in Lead, SD and the Piñon Flats PY array. These are two of the quietest small aperture arrays ever deployed in North America. The underground Homestake data has exceptionally low noise, and the borehole sensors of the PY array also have very low noise, making these arrays well suited to finding very weak signals. We find that our methods remove transient signals effectively from the data so that even low-SNR signals in the seismic background can be found and tested. Additionally, we find that the jackknife variance estimate is comparable to the noise floor, and we present initial evidence for solar g-modes in our data through the T2 test, a multivariate generalization of Student's t-test.

Caton, Ross C.↗
NASA NTRS1985

HF Doppler observations of acoustic waves excited by the earthquake

Ionospheric disturbances caused by the earthquake of a relatively small and large epicentral distance have been detected by a network of HF-Doppler sounders in central Japan and Kyoto station, respectively. The HF-Doppler data of a small epicentral distance, together with the seismic data, have been used to formulate a mechanism whereby ionospheric disturbances are produced by the Urakawa-Oki earthquake in Japan. Comparison of the dynamic spectra of these data has revealed experimentally that the atmosphere acts as a low-pass filter for upward-propagating acoustic waves. By surveying the earthquakes for which the magnitude M is larger than 6.0, researchers found the ionospheric effect in 16 cases of 82 seismic events. As almost all these effects have occurred in the daytime, it is considered that it may result from the filtering effect of the upward-propagating acoustic waves.

Ichinose, T.↗
NASA NTRS2013

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗
DOE OSTI2023

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗
NASA NTRS2011

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms.

Knapmeyer, Martin↗