Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic”

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 253 records · Page 14

Seismic anisotropy, dominant slip systems and phase transitions in the lowermost mantle

SUMMARY The presence of seismic anisotropy at the base of the Earth's mantle is well established, but there is no consensus on the deformation mechanisms in lower mantle minerals that could explain it. Strong anisotropy in magnesium post-perovskite (pPv) has been invoked, but different studies disagree on the dominant slip systems at play. Here, we aim to further constrain this by implementing the most recent results from atomistic models and high-pressure deformation experiments, coupled with a realistic composition and a 3-D geodynamic model, to compare the resulting deformation-induced anisotropy with seismic observations of the lowermost mantle. We account for forward and reverse phase transitions from bridgmanite (Pv) to pPv. We find that pPv with either dominant (001) or (010) slip can both explain the seismically observed anisotropy in colder regions where downwellings turn to horizontal flow, but only a model with dominant (001) slip matches seismic observations at the root of hotter large-scale upwellings. Allowing for partial melt does not change these conclusions, while it significantly increases the strength of anisotropy and reduces shear and compressional velocities at the base of upwellings.

58 GEOSCIENCES↗

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↗

Preliminary Nuclear Containment Vessel Modeling for Multi-Hazard Probabilistic Risk Assessment under Seismic Hazards and Concrete Degradation

The current practice for natural phenomena hazards (NPH) risk assessment of nuclear facilities is to compute the risk for each hazard independently and then compound the total risk as a combination of single hazard risks. This state of practice does not consider correlations between hazards and the cascading impacts to structures, systems, and components (SSCs), and could thus underestimate the NPH risk or overestimate the nuclear facility safety. Events such as the Fukushima Daiichi accident have highlighted the importance of multi-hazard risk considerations to nuclear power plants (NPPs) that quantify the cascading damage effects to SSCs in the risk models. Moreover, the current fleet of NPPs in the United States is aging; these NPPs are now expected to operate well beyond their initially planned design life. Aging-related deterioration can potentially decrease the capacity of critical structures such as containment vessels to withstand NPH. Such aging considerations may not be adequately accounted for by the current NPH risk assessment guidelines. This paper presents a preliminary modeling and simulation of a representative reinforced concrete containment vessel subjected to seismic mainshock and aftershock considering concrete degradation due to alkali silica reaction. The broader aim is to develop multi-hazard time-dependent fragility functions that could be subsequently used in the probabilistic risk assessment (PRA) model. The multi-hazard component comes into play due to the consideration of damage to the containment vessel under seismic loads and concrete degradation. Consideration of concrete degradation also brings into play the time-dependent nature of the containment vessel response. The response of the containment vessel under varying degrees of concrete degradation to seismic loads is investigated. The results presented are simulated using the Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) software for seismic analysis and the Blackbear software for concrete degradation and damage modeling. Both software are open source and developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE).

42 ENGINEERING↗

Seismic Survey 2016 Data at Crescent Valley, Nevada

In September 2016, 989 vertical-component seismic instruments were deployed for 75 hours at the Crescent Valley greenfield geothermal play area in Nevada. Data were recorded 12/9/16 - 12/15/16. Data are stored in individual files in one-minute increments in SEGD formats. See the metadata in Crescent Valley Station Locations.sp1, Crescent Valley Observer Logs.xls, and the attached Metadata link for details about the seismic station locations, seismic data logger specifications, instrumentation specifications, descriptions of data, a fracture finding summary, and the final report for the 2016 seismic survey done in Crescent Valley, Nevada. Note that no discrete events were identified, and a 1750m x 2250m x 300m PSET volume of acoustic energy was delivered.

15 GEOTHERMAL ENERGY↗

Utah FORGE 5-2419: Seismicity-Permeability Relationships Probed via Nonlinear Acoustic Imaging - 2024 Annual Workshop Presentation

This is a presentation on Seismicity-Permeability Relationships Probed via Nonlinear Acoustic Imaging by Pennsylvania State University, presented by Derek Elsworth. This video slide presentation discussed controls and acoustic signatures of aseismic through seismic evolution of relocation-stability-permeability on fractures in shear-reactivation as a precursor to, and a key predictor of, seismic moment magnitude of prospective triggered-seismicity. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 13, 2024.

15 GEOTHERMAL ENERGY↗

Seismic Isolation of Major Advanced Reactor Systems for Economic Improvement and Safety Assurance

Advanced reactor concepts currently under development in the United States [e.g., the High-Temperature Gas Reactor (HTGR) of X-Energy and the Molten Chloride Fast Reactor (MCFR) of TerraPower] are striving to improve the economics of nuclear power primarily by using passively and ‘walk-away’ safe technologies in the reactor. The recently published report, ‘Future of Nuclear Energy in a Carbon-Constrained World’ by the Massachusetts Institute of Technology Energy Institute (MIT EI) has shown that a significant contributor to the capital cost of NPPs is the ‘civil works’ i.e., the construction of the balance of plant (BoP: all SSCs except those involved in power generation, e.g., reactor vessel) such as the buildings, containment dome, foundation, etc., which amount to almost half of the overnight capital cost. The report suggests that using advanced construction technologies such as seismic isolation and high-performance concrete to significantly reduce the capital costs of advanced reactors. This report demonstrates the application of two cost-cutting solutions on advanced reactor designs that are currently under development. These solutions include, (1) seismic base isolation, and (2) risk- and cost-based seismic design optimization including seismic isolation of individual components.

42 ENGINEERING↗

Seismic Observations of the OSIRIS-REx Sample Return Capsule Reentry: Deployment, Signal Characteristics, and Wavefield Phenomenology

Controlled spacecraft reentries from interplanetary trajectories provide rare, well-characterized hypersonic sources for advancing seismoacoustic observation techniques. Here we present seismic observations of the OSIRIS-REx sample return capsule (SRC) reentry on 24 September 2023, recorded by 16 three-component nodal seismometers deployed near Eureka, Nevada, at ground distances of 7–20 km from the capsule trajectory. Air-to-ground coupled signals are detected at all stations, exhibiting impulsive onsets consistent with ballistic shock arrivals from the descending Mach cone. We characterize the seismic wavefield through signal amplitude, period, waveform cross-correlation, and array processing. Signal periods decrease systematically with increasing distance from the trajectory within the airport array, indicating that higher-frequency content becomes more prominent at greater offsets, opposite to expectations from geometric spreading and atmospheric absorption. Seismic array processing identifies frequency-dependent back-azimuth variations whose origin remains unresolved; possible contributing factors include source geometry, scattering by fine-scale layered structure in the stratosphere, and near-surface effects. These observations document a spatially complex seismic wavefield from a well-characterized hypersonic line source and provide constraints for future modeling of atmospheric propagation and air-to-ground coupling.

OSIRIS-REx↗

Mass Wasting In Planetary Environments: Implications For Seismicity

On Earth, mass wasting events such as rock falls and landslides are well known conse-quences of seismic activity. Through a variety of re-mote sensing techniques, tectonic faults have been pos-itively identified on all four of the inner planets, Earth's Moon, several outer planet satellites, and aster-oids. High-resolution imaging has furthermore ena-bled positive identification of mass wasting events on many of these bodies. On Mars, it has been suggested that fallen boulders may be indicative of pale-omarsquakes. On the Moon, meteor impacts and moonquakes have likewise been suggested as potential triggering mechanisms for mass wasting. Indeed, we know from the Apollo era that the Moon experienc-es a wide variety of seismicity. Seismicity estimates play an important role in creat-ing regional geological characterizations, which are useful not only for understanding a planet's formation and evolution, but also of key importance to site selec-tion for landed missions. Here we investigate the re-gional effects of seismicity in planetary environments with the goal of determining whether surface features such as landslides and boulder trails on the Moon, Mars, and Mercury could be triggered by fault motion. We attempt to quantify the amount of near-source ground shaking necessary to mobilize the mate-rial observed in various instances of mass wasting.

seismicity↗

Cardinal: Seismic and Geoacoustic Array Processing

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Adaptive Array↗

In-situ observation of pre-, co- and post-seismic shear slip at 1.5 km depth

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and transient afterslip. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. < p >The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment.

Schoenball, Martin↗

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗

Seismic source characterization of the Arabian Peninsula and Zagros Mountains from regional moment tensor and coda envelopes

Reliable estimates of moment magnitude and source mechanism for seismic events in the Middle East can be challenging due to a small number of openly available stations, the complex tectonic setting, and regions of high attenuation. Access to high-quality waveform data from well-calibrated regional seismic stations is fundamental in producing robust and stable estimates of earthquake source parameters, particularly when measurements of absolute ground motion amplitudes are required. Earthquake source mechanism and moment magnitude are invaluable information in the assessment of seismic hazard, plate motions, and the characterization of faults and regional stress field. The expansion of the Kingdom of Saudi Arabia’s national seismic network in the past several decades provides an opportunity to develop the capabilities of routine focal mechanism and moment magnitude estimations. In this study, we use time-domain full waveform moment tensor inversion and coda envelope–derived amplitude measurements to solve for earthquake source mechanism, moment magnitudes, and their source type. Furthermore, we compared the moment magnitudes calculated from the two methods and publicly available earthquake catalogs and discuss the implications of the obtained source parameters.

58 GEOSCIENCES↗

A general approach to seismic inversion with automatic differentiation

Imaging Earth structure or seismic sources from seismic data involves minimizing a target misfit function, and is commonly solved through gradient-based optimization. The adjoint-state method has been developed to compute the gradient efficiently; however, its implementation can be time-consuming and difficult. We develop a general seismic inversion framework to calculate gradients using reverse-mode automatic differentiation. The central idea is that adjoint-state methods and reverse-mode automatic differentiation are mathematically equivalent. Here, the mapping between numerical PDE simulation and deep learning allows us to build a seismic inverse modeling library, ADSeismic, based on deep learning frameworks, which supports high performance reverse-mode automatic differentiation on CPUs and GPUs. We demonstrate the performance of ADSeismic on inverse problems related to velocity model estimation, rupture imaging, earthquake location, and source time function retrieval. ADSeismic has the potential to solve a wide variety of inverse modeling applications within a unified framework.

58 GEOSCIENCES↗

Seismic monitoring of a small CO 2 injection using a multi-well DAS array: Operations and initial results of Stage 3 of the CO 2 CRC Otway project

Active time-lapse seismic is widely employed for monitoring CO 2 geosequestration due to its ability to track the distribution of fluids in space and time. However, standard 4D seismic monitoring suffers from several challenges, including high cost, disruption to other land uses, and, consequently, relatively large intervals between monitor surveys. Some of these challenges can be mitigated using permanently installed sources and receivers. Such an approach was tested at the CO 2 CRC Otway site by continuous offset VSP monitoring of 15,000 t of supercritical CO 2 injected into an aquifer 1,500 m deep with nine permanent seismic sources (surface orbital vibrators or SOVs) and five downhole fibre-optic receivers. This continuous monitoring is complemented by multi-well 4D VSP using a mobile vibroseis source and the same DAS receivers, which included one baseline and two monitor surveys after injection of 4,000 and 12,000 t of CO 2 . The continuous DAS-SOV monitoring detected an abrupt increase of travel times below the injection interval on the second day of injection (after injection of 300 t of CO 2 ) and tracked the growth of the areal CO 2 plume by mapping changes of reflection amplitudes. The plume is also detected by time-lapse changes of reflection amplitudes in multi-well 4D VSPs. The plume images obtained from continuous offset VSP and 4D VSP are broadly consistent with each other but with some differences due to differences in illumination, lateral variations of velocities and seismic anisotropy. Furthermore, these differences also serve as a measure of uncertainty of 4D VSP images.

4D VSP↗

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗

Seismic Ambient Noise Analyses Reveal Changing Temperature and Water Signals to 10s of Meters Depth in the Critical Zone

The critical zone sustains terrestrial life, but we have few tools to explore it efficiently beyond the first few meters of the subsurface. Using analyses of high-frequency ambient seismic noise from densely spaced seismometers deployed in the forested Shale Hills subcatchment of the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), we show that temporal changes in seismic velocities at depths from ~1 m to tens of m can be detected. Furthermore, these changes are driven by variations at the land surface. The Moving-Window Cross-Spectral (MWCS) method was employed to measure seismic-velocity changes in coda waves at hourly resolution in 10 different frequency bands. We observed a diurnal signal, a seasonal signal, and a meteorological-event-based signal. These signals were compared to time-series measurements of precipitation, well water levels, soil moisture, soil temperature, air temperature, latent heat flux, and air pressure in the heavily instrumented catchment. Most of the velocity changes can be explained by variations in temperature that result in thermoelastic strains that propagate to depth. But some double minima in seismic velocity time-series observed after large rain events were attributed in part to the effects of water infiltration. These results show that high-frequency ambient noise data may in some locations be used to detect changes in the critical zone from ~1 to ~100 m or greater depth with hourly resolution. But interpretation of such data requires multiple environmental data sets to deconvolve the complex interrelationships among thermoelastic and hydrological effects in the subsurface critical zone.

54 ENVIRONMENTAL SCIENCES↗

Seismic Imaging of a Shale Landscape Under Compression Shows Limited Influence of Topography–Induced Fracturing

We used seismic refraction to image the P-wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill-valley-hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. Here, this contradiction could be explained by the greater importance of infiltration-driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.

58 GEOSCIENCES↗