Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic methods”

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 289 records · Page 16

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗

Back-Azimuth Estimation of Air-to-Ground Coupled Infrasound from Transverse Coherence Minimization [Slides]

Infrasonic signals are generated by a variety of natural and anthropogenic sources and can propagate long distances through a windy, layered atmosphere. When these signals impinge on the ground, the resulting deformation can be recorded by buried seismometers, we call this signal a ground coupled airwave (GCA). GCAs are relatively common, and several methods now exist to leverage an infrasound microphone to separate GCAs from pure seismic motion.

47 OTHER INSTRUMENTATION↗

Passive seismic experiment

The design, deployment, and operation of the Apollo 16 passive seismic experiment (PSE) are discussed. Since activation, all elements of the PSE have operated as planned, with the exception of the sensor thermal control system. Significant progress in the measurement of meteoroid flux in near-earth space has been made, along with dilineation of active moonquake source regions. The data obtained indicate that moonquakes are concentrated at great depth (800 to 1000 km) and that the apparent disparity between meteoroid flux estimtes based on lunar crater counts and those from earth-based observations can be resolved by seismic measurements in favor of the lower flux indicated by the crater count method. The results obtained from the PSE are summarized and their significance is discussed in detail.

Latham, G. V.↗

Applications of aerospace technology to petroleum exploration. Volume 1: Efforts and results

The feasibility of applying aerospace techniques to help solve significant problems in petroleum exploration is studied. Through contacts with petroleum industry and petroleum service industry, important petroleum exploration problems were identified. For each problem, areas of aerospace technology that might aid in its solution were also identified where possible. Topics selected for investigation include: seismic reflection systems; down-hole acoustic techniques; identification of geological analogies; drilling methods; remote geological sensing; and sea floor imaging and mapping. Specific areas of aerospace technology are applied to 21 concepts formulated from the topics of concern.

Jaffe, L. D.↗

Interferometric detection of freeze-thaw displacements of Alaskan permafrost using ERS-1 data

The possibility of making large scale (50 km) measurements of motions of the earth's surface with high resolution (10 m) and very high accuracy (1 cm) from multipass SAR interferometry was established in 1989. Other experiments have confirmed the viability and usefulness of the method. Work is underway in various groups to measure displacements from volcanic activity, seismic events, glacier motion, and in the present study, freeze-thaw cycles in Alaskan permafrost. The ground is known to move significantly in these cycles, and provided that freezing does not cause image decorrelation, it should be possible to measure both ground swelling and subsidence. The authors have obtained data from multiple passes of ERS-1 over the Toolik Lake region of northern Alaska of suitable quality for interferometry. The data are processed into images, and single interferograms are formed in the usual manner. Phase unwrapping is performed, and the multipass baselines are estimated from the images using both orbit ephemerides and scene tie points. The phases are scaled by the baseline ratio, and a double-difference interferogram (DDI) is formed. It is found that there is a residual 'saddle-shape' phase error across the image, which is postulated to be caused by a small divergence (10(exp -2) deg.) in the orbits. A simulation of a DDI from divergent orbits confirms the shape and magnitude of the error. A two-dimensional least squares fit to the error is performed, which is used to correct the DDI. The final, corrected DDI shows significant phase (altitude) changes over the period of the observation.

Werner, Charles L.↗

A New Approach to Predict Hydrogeological Parameters Using Shear Waves from the Multichannel Analysis of Surface Waves Method

For near-surface contaminant characterization, the accurate prediction of hydrogeological parameters in anisotropic and heterogeneous environments has been a challenge since the last decades. However, recent advances in near-surface geophysics have facilitated the use of geophysical data for hydrogeological characterization in the last few years. A pseudo 3-D high resolution P-wave shallow seismic reflection survey was performed at the P Reactor Area, Savannah River Site, South Carolina in order to delineate and predict migration pathways of a large contaminant plume including trichloroethylene. This contaminant plume originates from the northwest section of the reactor facility that is located within the Upper Atlantic Coastal Plain. The data were collected with 40 Hz geophones, an accelerated weight-drop as seismic source and 1 m receiver spacing with near- and far-offsets of 0.5 and 119.5 m, respectively. In such areas with near-surface contaminants, a detailed subsurface characterization of the vadose zone hydraulic parameters is very important. Indeed, an inexpensive method of deriving such parameters by the use of seismic reflection surveys is beneficial, and our approach uses the relationship between seismic velocity and hydrogeological parameters together with empirical observations relating porosity to permeability and hydraulic conductivity. Shear wave velocity ( V s ) profiles were estimated from surface wave dispersion analysis of the seismic reflection data and were subsequently used to derive hydraulic parameters such as porosity, permeability, and hydraulic conductivity. Additional geophysical data including core samples, vertical seismic profiling, surface electrical resistivity tomography, natural gamma and electrical resistivity logs allowed for a robust assessment of the validity and geological significance of the estimated V s and hydrogeological models. The results demonstrate the usefulness of this approach for the upper 15 m of shallow unconsolidated sediments even though the survey design parameters were not optimal for surface wave analysis due to the higher than desired frequency geophones.

Engineering↗

Impact-Seismic Investigations of the InSight Mission

Impact investigations will be an important aspect of the InSight mission. One of the scientific goals of the mission is a measurement of the current impact rate at Mars. Impacts will additionally inform the major goal of investigating the interior structure of Mars. In this paper, we review the current state of knowledge about seismic signals from impacts on the Earth, Moon, and laboratory experiments. We describe the generalized physical models that can be used to explain these signals. A discussion of the appropriate source time function for impacts is presented, along with spectral characteristics including the cutoff frequency and its dependence on impact momentum. Estimates of the seismic efficiency (ratio between seismic and impact energies) vary widely. Our preferred value for the seismic efficiency at Mars is 5 × 10−4, which we recommend using until we can measure it during the InSight mission, when seismic moments are not used directly. Effects of the material properties at the impact point and at the seismometer location are considered. We also discuss the processes by which airbursts and acoustic waves emanate from bolides, and the feasibility of detecting such signals. We then consider the case of impacts on Mars. A review is given of the current knowledge of present-day cratering on Mars: the current impact rate, characteristics of those impactors such as velocity and directions, and the morphologies of the craters those impactors create. Several methods of scaling crater size to impact energy are presented. The Martian atmosphere, although thin, will cause fragmentation of impactors, with implications for the resulting seismic signals. We also benchmark several different seismic modeling codes to be used in analysis of impact detections, and those codes are used to explore the seismic amplitude of impactinduced signals as a function of distance from the impact site. We predict a measurement of the current impact flux will be possible within the timeframe of the prime mission (one Mars year) with the detection of ∼ a few to several tens of impacts. However, the error bars on these predictions are large. Specific to the InSight mission, we list discriminators of seismic signals from impacts that will be used to distinguish them from marsquakes. We describe the role of the InSight Impacts Science Theme Group during mission operations, including a plan for possible night-time meteor imaging. The impacts detected by these methods during the InSight mission will be used to improve interior structure models, measure the seismic efficiency, and calculate the size frequency distribution of current impacts.

Ingrid Daubar↗

Measurements of mantle wave velocities and inversion for lateral heterogeneity and anisotropy. II - Analysis by the single-station method

In the present investigation, the single-station method reported by Brune et al. (1960) is utilized for an analysis of long-period Love(G) and Rayleigh(R) waves recorded on digital seismic networks. The analysis was conducted to study the lateral heterogeneity of surface wave velocities. The data set is examined, and a description is presented of the single-station method. Attention is given to an error analysis for velocity measurements, the estimation of geographical distribution of surface wave velocities, the global distribution of surface wave velocities, and the correlation between the surface wave velocities and the heat flow on the geoid. The conducted measurements and inversions of surface wave velocities are used as a basis to derive certain conclusions. It is found that the application of the single-station method to long-period surface waves recorded on digital networks makes it possible to reach an accuracy level comparable to great circle velocity measurements.

Nakanishi, I.↗

LocOO3D User's Manual

LocOO3D is a software tool that computes geographical locations for seismic events at regional to global scales. This software has a rich set of features, including the ability to use custom 3D velocity models, correlated observations and master event locations. The LocOO3D software is especially useful for research related to seismic monitoring applications, since it allows users to easily explore a variety of location methods and scenarios and is compatible with the CSS3.0 data format used in monitoring applications. The LocOO3D software, User's Manual, and Examples are available on the web at: https://github.com/sandialabs/LocOO3D For additional information on GeoTess, SALSA3D, RSTT, and other related software, please see: https://github.com/sandialabs/GeoTessJava, www.sandia.gov/geotess, www.sandia.gov/salsa3d, and www.sandia.gov/rstt

58 GEOSCIENCES↗

Analysis of upper mantle structure using wave field continuation of P waves

Wave field continuation theory, which allows transformation of the seismic record section data directly into velocity-depth space, is tested for upper mantle analysis using a large array-recorded data set obtained at the 200-station Caltech-USGS Southern California Seismic Network that is representative of the structure beneath the gulf of California. The method's resolution capability is illustrated by the comparison of the slant stacks and downward continuation of both synthetic and data record sections. It is stressed that when high-quality, densely sampled digital data are available, the technique is easy to implement, provides an inversion which contains all the data in the global format, and produces an objective estimate of depth resolution as a function of ray parameter.

Walck, M. C.↗

Plate motions and deformations from geologic and geodetic data

An analysis of geodetic data in the vicinity of the Crustal Dynamics Program (CDP) site at Vandenberg Air Force Base (VNDN) is presented. The utility of space-geodetic data in the monitoring of transient strains associated with earthquakes in tectonically active areas like California is investigated. Particular interest is in the possibility that space-geodetic methods may be able to provide critical new data on deformations precursory to large seismic events. Although earthquake precursory phenomena are not well understood, the monitoring of small strains in the vicinity of active faults is a promising technique for studying the mechanisms that nucleate large earthquakes and, ultimately, for earthquake prediction. Space-geodetic techniques are now capable of measuring baselines of tens to hundreds of kilometers with a precision of a few parts in 108. Within the next few years, it will be possible to record and analyze large-scale strain variations with this precision continuously in real time. Thus, space-geodetic techniques may become tools for earthquake prediction. In anticipation of this capability, several questions related to the temporal and spatial scales associated with subseismic deformation transients are examined.

Jordan, Thomas H.↗

Convolutional Neural Networks for Signal Detection

Currently, traditional methods such as short-term average/long-term average (STA/LTA) are used to detect arrivals in three-component seismic waveform data. Accurately establishing the identity and arrival of these waves is helpful in detecting and locating seismic events. Convolutional Neural Networks (CNNs) have been shown to significantly improve performance at local distances. This work will expand the use of CNNs to more remote distances and lower magnitudes. Sandia National Labs (SNL) will explore the advantages and limits of a particular approach and investigate requirements for expanding this technique to different types, distances, and magnitudes of events in the future. The team will describe detailed performance results of this method tuned on a curated dataset from Utah with its expert-defined arrival picks.

58 GEOSCIENCES↗

MASTODON: An Open-Source Software for Seismic Analysis and Risk Assessment of Critical Infrastructure

Seismic analysis and risk assessment of safety-critical infrastructure like hospitals, nuclear power plants, dams, and facilities handling radioactive materials involve computationally intensive numerical models and coupled multiphysics scenarios. They are also performed in a strict regulatory environment that requires high software quality assurance standards, and in the case of safety-related nuclear facilities, a conformance to the American Society of Mechanical Engineers Nuclear Quality Assurance (NQA-1) standard. This paper introduces the open-source finite-element software, MASTODON (Multi-hazard Analysis of Stochastic Time-Domain Phenomena), which implements state-of-the-art seismic analysis and risk assessment tools in a quality-controlled environment. MASTODON is built on MOOSE (Multi-physics Object-Oriented Simulation Environment), which is a highly parallelizable, NQA-1 conforming, coupled multiphysics, finite-element framework developed at Idaho National Laboratory. MASTODON is capable of fault rupture and source-to-site wave propagation using the domain reduction method, nonlinear site response, and soil-structure interaction analysis, implicit and explicit time integration, automated stochastic simulations, and seismic probabilistic risk assessment. When coupled with other MOOSE applications, MASTODON can also solve strongly and weakly coupled multiphysics problems. This paper presents a summary of the capabilities of MASTODON and some demonstrative examples.

42 ENGINEERING↗

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗

cjohnson-LANL/GRL_Kilauea

The python routines are outlined in detail to perform the methods and results in the manuscript under review in the journal Geophysical Research Letters titled “Seismic features predict ground motions during repeating caldera collapse sequence” with LA-UR-23-33345. All routines are written in open source python and were applied to publicly available data sets. The codes formats the data into the appropriate structure required to train a boosted tree regression model. Other codes produce figure results.

Johnson, Christopher W↗

The Strange Case of Ground-Coupled Airwaves on Seismoacoustic Stations at Local to Near-Regional Scales

Here we investigate the application of ground-coupled airwaves observed by seismoacoustic stations at local to near-regional scales to detect signals of interest and determine back-azimuth information. Ground-coupled airwaves are created from incident pressure waves traveling through the atmosphere that couple to the earth and transmit as a seismic wave with retrograde elliptical motion. Previous studies at sub-local scales (<10 km from a source of interest) found the back-azimuth to the source could be accurately determined from seismoacoustic signals recorded by acoustic and 3-component seismic sensors spatially separated on the order of 10 to 150 m. The potential back-azimuth directions are estimated from the coherent signals between the acoustic and vertical seismic data, via a propagation-induced phase shift of the seismoacoustic signal. A unique solution is then informed by the particle motion of the 3-component seismic station, which was previously found to be less accurate than the seismoacoustic-sensor method. We investigate the applicability of this technique to greater source-receiver distances, from 50-100 km and up to 400 km, which contains pressure waves with tropospheric and stratospheric ray paths, respectively. Specifically, we analyze seismoacoustic sources with ground truth from rocket motor fuel elimination events at the Utah Test and Training Range (UTTR) as well as a 2020 rocket launch in Southern California. From these sources we observe evidence that while coherent signals can be seen from both sources on multiple seismoacoustic station pairs, the determined ground-coupled airwave back-azimuths are more complicated than results at more local scales. Our findings suggest more complex factors including incidence angle, coupling location, subsurface material, and atmospheric propagation effects need to be fully investigated before the ground-coupled airwave back-azimuth determination method can be applied or assessed at these further distances.

58 GEOSCIENCES↗

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↗