Search NASA⌕ Search

SEARCH · Search NASA

Results for “Seismic Data Compression”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Applying Compression-Based Metrics to Seismic Data in Support of Global Nuclear Explosion Monitoring

The analysis of seismic data for evidence of possible nuclear explosion testing is a critical global security mission that relies heavily on human expertise to identify and mark seismic signals embedded in background noise. To assist analysts in making these determinations, we adapted two compression distance metrics for use with seismic data. First, we demonstrated that the Normalized Compression Distance (NCD) metric can be adapted for use with waveform data and can identify the arrival times of seismic signals. Then we tested an approximation for the NCD called Sliding Information Distance (SLID), which can be computed much faster than NCD. We assessed the accuracy of the SLID output by comparing it to both the Akaike Information Criterion (AIC) and the judgments of expert seismic analysts. Our results indicate that SLID effectively identifies arrival times and provides analysts with useful information that can aid their analysis process.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Seismology in the Cloud: A New Streaming Workflow

Data-intensive research in seismology is experiencing a recent boom, driven in part by large volumes of available data and advances in the growing field of data science. However, there are significant barriers to processing large data volumes, such as long retrieval times from data repositories, complex data management, and limited computational resources. New tools and platforms have reduced the barriers to entry for scientific cluster computing, including the maturation of the commercial cloud as an accessible instrument for research. Here, we build a customized research cluster in the cloud to test a new workflow for large-scale seismic analysis, in which data are processed as a stream (retrieved on-the-fly and acted upon without storing), with data from the Incorporated Research Institutions for Seismology Data Management Center. We use this workflow to deploy a spectral peak detection algorithm over 5.6 TB of compressed continuous seismic data from 2074 stations of the USArray Transportable Array EarthScope network. Using a 50-node cluster in the cloud, we completed the noise survey in 80 hr, with an average data throughput of 1.7 GB per minute. By varying cluster sizes, we find the scaling of our analysis to be sublinear, due to a combination of algorithmic limitations and data center response times. The cloud-based streaming workflow represents an order-of-magnitude increase in acquisition and processing speed compared to a traditional download-store-process workflow, and offers the additional benefits of employing a flexible, accessible, and widely used computing architecture. It is limited, however, due to its reliance on Internet transfer speeds and data center service capacity, and may not work well for repeated analyses or those for which even higher data throughputs are needed. These research applications will require a new class of cloud-native approaches in which both data and analysis are in the cloud.

58 GEOSCIENCES↗
DOE OSTI2021

A Consensus Equilibrium Approach for 3-D Land Seismic Shots Recovery

Physical and budget constraints often result in inadequate sampling for accurate subsurface imaging. Preprocessing approaches, such as missing trace interpolation, are typically employed to enhance seismic data in such cases. The compressed sensing (CS) framework has been applied for modeling missing seismic data, which is estimated by sparsity-based computational algorithms. While existing work mainly focuses on recovering missing traces resulting from receiver subsampling, source subsampling has greater economical advantages, as sources are more expensive than receivers. Moreover, stronger image models different from sparsity have not been explored for source recovery. This work presents a consensus equilibrium (CE) approach to recover missing seismic shots, which enables to incorporate various regularization operators modeling different data priors. Here, simulation results from a real 3-D land seismic dataset demonstrate that the CE approach provides more accurate estimations of the linear and hyperbolic events in the recovered shots, compared with pure sparsity-based reconstructions.

58 GEOSCIENCES↗
DOE OSTI2022

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↗
DOE OSTI2019

Effect of Carbon on the Volume of Solid Iron at High Pressure: Implications for Carbon Substitution in Iron Structures and Carbon Content in the Earth’s Inner Core

Understanding the effect of carbon on the density of hcp (hexagonal-close-packed) Fe-C alloys is essential for modeling the carbon content in the Earth’s inner core. Previous studies have focused on the equations of state of iron carbides that may not be applicable to the solid inner core that may incorporate carbon as dissolved carbon in metallic iron. Carbon substitution in hcp-Fe and its effect on the density have never been experimentally studied. We investigated the compression behavior of Fe-C alloys with 0.31 and 1.37 wt % carbon, along with pure iron as a reference, by in-situ X-ray diffraction measurements up to 135 GPa for pure Fe, and 87 GPa for Fe-0.31C and 109 GPa for Fe-1.37C. The results show that the incorporation of carbon in hcp-Fe leads to the expansion of the lattice, contrary to the known effect in body-centered cubic (bcc)-Fe, suggesting a change in the substitution mechanism or local environment. The data on axial compressibility suggest that increasing carbon content could enhance seismic anisotropy in the Earth’s inner core. The new thermoelastic parameters allow us to develop a thermoelastic model to estimate the carbon content in the inner core when carbon is incorporated as dissolved carbon hcp-Fe. The required carbon contents to explain the density deficit of Earth’s inner core are 1.30 and 0.43 wt % at inner core boundary temperatures of 5000 K and 7000 K, respectively.

58 GEOSCIENCES↗
DOE OSTI2021

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

CODAS Data from Oliktok Point, Beaufort Sea, Alaska

Cryosphere/Ocean Distributed Acoustic Sensing (CODAS) data collected from the Beaufort Sea, Alaska, using ~37.4 km of dark telecommunications fiber located at Oliktok Point, Alaska. Data were collected with a Silixa iDAS, using 10 m gauge length, 2 m spatial resolution, and 1000 Hz sample rate. Provided here are the DAS-recorded time series for the rapid refreeze event described in Baker & Abbott (2022) (see link below). This covers a date range of 2021-11-10 15:00 UTC to 2021-11-11 17:00 UTC. Data have been decimated to 100 Hz and 20 m (i.e., every 10th channel for 1831 channels, total), as used in Baker & Abbott (2022). Data have been extracted from raw format into 1-hour long .sac* files and organized into directories by channel number, spanning channels 100 to 18400. Time series units are nano-strainrate (nm/m/s). For distribution, data have been compressed into .zip files containing all time series files for 100 channels. *For information on the Seismic Analysis Code (SAC) file format: https://seiscode.iris.washington.edu/projects/sac

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Screening IMS Seismic Detections using Dynamic Correlation Processor: Progress report for 2020

LLNL completed the FY20 workplan entailing application of the Dynamic Correlation Processor (DCP) to International Monitoring System (IMS) primary seismic data in 2019 and 2020. The effective analysis period was curtailed to January 1 through November 12, 2019 due to a data compression change that rendered later data unusable. Fixing the data compression error is a priority for our data source, and we plan to extend the period of data analysis when the fix is implemented. We found that the percentage overlap between detections generated by the International Data Centre (IDC) DFX software and detections generated by DCP varied widely by station. Overlap is 30-40% at some stations, suggesting that a large percentage of DFX detections could be screened (removed) from the detection-event association algorithm. However, detection overlap was much lower, a few percent, at many stations and the percentage overlap was low at stations known to have many mines and other local sources of seismicity. Unexpectedly low detection overlap at many stations prompted us to examine results for the ARCES station in detail. There is considerable mine activity near ARCES, but our initial FY20 analysis found only ~2% DCP/DFX detection overlap. Differences in the pre-processing (e.g. beam recipes) used by DCP and DFX appear to be the cause of low detection overlap at ARCES. DCP uses wideband filters to improve detector robustness and the IDC uses narrow-band pre-filtering to improve sensitivity. After reprocessing ARCES using the IDC beam recipe, DCP/DFX overlap increased to ~74%. This result shows that we must reprocess the IMS network using the IDC station-specific beam recipes if we are to effectively screen DFX detections and improve automatic event building performance at the IDC.

58 GEOSCIENCES↗
DOE OSTI2025

Seismic Tremors From Sea‐Landfast Ice Interactions Near Utqiaġvik, Alaska

The mechanical state of Arctic landfast sea ice remains poorly constrained due to limited observations. This study investigates interactions between drifting sea ice and the coastal landfast ice near Utqiaġvik, Alaska by integrating data from broadband seismometer, Distributed Acoustic Sensing, and marine radar. We find that decreases in sea ice velocity, marking transitions from drift to compressive contact, coincide with increased seismic energy. Tremor characteristics vary seasonally with ice conditions. In January, dense ice packs produced sustained harmonic tremors with gliding and U-shaped spectral features, consistent with repetitive stick-slip motion at the ice–ice or ice–ground interface under velocity-weakening friction. In April, smaller fragmented floes generated short-lived, chaotic tremors linked to brittle failure and spatially dispersed impacts. These findings demonstrate that seismic tremors encode the mechanical properties of interacting ice, offering a new tool to distinguish ice regimes and monitor evolving Arctic coastal dynamics under climate change.

58 GEOSCIENCES↗
DOE OSTI2020

The Complicated Link Between Material Properties and Microfracture Density for an Underground Explosion in Granite

Geophysical techniques are often implemented as quick and inexpensive ways to locate and characterize fractures in the subsurface, which is important for a number of geoscience fields. Seismic velocities are the most widely used proxies for identification of fractures, but the correlation is not always well-defined. In this study we present material property data: unconfined compressive strength (UCS), bulk density (ρ), Young's modulus (E), Poisson's ratio (ν), P wave velocity (V p ), and S wave velocity (V s ), in conjunction with microfracture densities measured on samples of granite collected before and after underground chemical explosions. Results indicate the relationship between fractures and material properties is complex, even in this single-lithology environment. We interpret that this complexity arises from varying fracture mechanisms (e.g. dilation-inducing fractures vs compression-inducing fractures) in different parts of the core, due to differences in stress conditions. Additional complexity may result from chemical interactions between the fresh fractures and the fluids in the area. Water content appears to have a significant, if not dominant, role in the unconfined compressive strength (UCS) of the samples. We suggest caution when using elastic property measurements as a proxy for fracturing in areas of explosion-induced damage, or in other areas where a variety of mechanisms induce fracturing.

58 GEOSCIENCES↗
DOE OSTI2018

Unconfined Compression Results on Core from Borehole U-15n#10, NNSS, in support of SPE

The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This memorandum reports on an initial phase of unconfined compression testing from corehole U-15n#10. Specimens tested came from the U-15n#10 core hole, which was drilled at -60° to the horizontal in granitic rock (quartz monzonite) after the second SPE shot (SPE-2). Figure 1 illustrates at the surface, the core hole was approximately 90 feet from the central SPE borehole. Corehole U 15n#10 intersects the central SPE borehole (U-15n) at a core depth of 170 feet (approximately 150 feet vertical depth) which is within the highly damaged zone of SPE-2. The U-15n#10 location is the site of the first, second and third SPE's, in Area 15 of the NNSS.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Unconfined Compression Results on Core from Boreholes U-15n#12 and U-15n#13, NNSS in support of SPE

The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This report summarizes unconfined compression testing (UCS) from coreholes U-15n#12 and U-15n#13 and compares those datasets to UCS results from coreholes U-15n and U-15n#10. U-15n#12 corehole was drilled at -60° to the horizontal and U-15n#13 was drilled vertically in granitic rock (quartz monzonite) after the third SPE shot. Figure 1 illustrates at the surface, U 15n#12 and U-15n#13 coreholes were approximately 30 meters and 10 meters from the central SPE borehole (U-15n) respectively. Corehole U-15n#12 intersects the central SPE borehole (U 15n) at a core depth of 174 feet (approximately 150 feet vertical depth). The location of U 15n#12 and U-15n#13 is the site of the first, second and third SPE's, in Area 15 of the NNSS.

58 GEOSCIENCES↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An autoencoder compression approach for accelerating large-scale inverse problems

Abstract Partial differential equation (PDE)-constrained inverse problems are some of the most challenging and computationally demanding problems in computational science today. Fine meshes required to accurately compute the PDE solution introduce an enormous number of parameters and require large-scale computing resources such as more processors and more memory to solve such systems in a reasonable time. For inverse problems constrained by time-dependent PDEs, the adjoint method often employed to compute gradients and higher order derivatives efficiently requires solving a time-reversed, so-called adjoint PDE that depends on the forward PDE solution at each timestep. This necessitates the storage of a high-dimensional forward solution vector at every timestep. Such a procedure quickly exhausts the available memory resources. Several approaches that trade additional computation for reduced memory footprint have been proposed to mitigate the memory bottleneck, including checkpointing and compression strategies. In this work, we propose a close-to-ideal scalable compression approach using autoencoders to eliminate the need for checkpointing and substantial memory storage, thereby reducing the time-to-solution and memory requirements. We compare our approach with checkpointing and an off-the-shelf compression approach on an earth-scale ill-posed seismic inverse problem. The results verify the expected close-to-ideal speedup for the gradient and Hessian-vector product using the proposed autoencoder compression approach. To highlight the usefulness of the proposed approach, we combine the autoencoder compression with the data-informed active subspace (DIAS) prior showing how the DIAS method can be affordably extended to large-scale problems without the need for checkpointing and large memory.

Mathematics↗
DOE OSTI2020

Pressure-stabilized divalent ozonide CaO 3 and its impact on Earth's oxygen cycles

High pressure can drastically alter chemical bonding and produce exotic compounds that defy conventional wisdom. Especially significant are compounds pertaining to oxygen cycles inside Earth, which hold key to understanding major geological events that impact the environment essential to life on Earth. Here we report the discovery of pressure-stabilized divalent ozonide CaO3 crystal that exhibits intriguing bonding and oxidation states with profound geological implications. Our computational study identifies a crystalline phase of CaO 3 by reaction of CaO and O 2 at high pressure and high temperature conditions; ensuing experiments synthesize this rare compound under compression in a diamond anvil cell with laser heating. High-pressure x-ray diffraction data show that CaO 3 crystal forms at 35 GPa and persists down to 20 GPa on decompression. Analysis of charge states reveals a formal oxidation state of -2 for ozone anions in CaO 3 . These findings unravel the ozonide chemistry at high pressure and offer insights for elucidating prominent seismic anomalies and oxygen cycles in Earth’s interior. We further predict multiple reactions producing CaO 3 by geologically abundant mineral precursors at various depths in Earth’s mantle.

58 GEOSCIENCES↗
DOE OSTI2023

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗