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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Large-scale seismic waveform quality metric calculation using Hadoop

Here in this work we investigated the suitability of Hadoop MapReduce and Apache Spark for large-scale computation of seismic waveform quality metrics by comparing their performance with that of a traditional distributed implementation. The Incorporated Research Institutions for Seismology (IRIS) Data Management Center (DMC) provided 43 terabytes of broadband waveform data of which 5.1 TB of data were processed with the traditional architecture, and the full 43 TB were processed using MapReduce and Spark. Maximum performance of ~0.56 terabytes per hour was achieved using all 5 nodes of the traditional implementation. We noted that I/O dominated processing, and that I/O performance was deteriorating with the addition of the 5th node. Data collected from this experiment provided the baseline against which the Hadoop results were compared. Next, we processed the full 43 TB dataset using both MapReduce and Apache Spark on our 18-node Hadoop cluster. We conducted these experiments multiple times with various subsets of the data so that we could build models to predict performance as a function of dataset size. We found that both MapReduce and Spark significantly outperformed the traditional reference implementation. At a dataset size of 5.1 terabytes, both Spark and MapReduce were about 15 times faster than the reference implementation. Furthermore, our performance models predict that for a dataset of 350 terabytes, Spark running on a 100-node cluster would be about 265 times faster than the reference implementation. We do not expect that the reference implementation deployed on a 100-node cluster would perform significantly better than on the 5-node cluster because the I/O performance cannot be made to scale. Finally, we note that although Big Data technologies clearly provide a way to process seismic waveform datasets in a high-performance and scalable manner, the technology is still rapidly changing, requires a high degree of investment in personnel, and will likely require significant changes in other parts of our infrastructure. Nevertheless, we anticipate that as the technology matures and third-party tool vendors make it easier to manage and operate clusters, Hadoop (or a successor) will play a large role in our seismic data processing.

58 GEOSCIENCES↗
DOE OSTI2017

Faulting of natural serpentinite: Implications for intermediate-depth seismicity

The seismic potential of serpentinites at high pressure was investigated via deformation experiments on cored natural serpentinite samples, during which micro-seismicity was monitored by recording Acoustic Emissions (AEs). Deformation was performed at pressures of 3–5 GPa, using a Deformation-DIA device, and over a wide range of temperatures, both within and outside antigorite's stability field. Below 400 °C, serpentinite deformation involves “silent” semi-brittle mechanisms, even in cases where strain localization is observed. At high temperature (i.e., above 600 °C), despite conditions propitious to dehydration embrittlement (i.e., fast strain rates and reaction kinetics), joint deformation and dehydration lead to ductile shear, without generation of AEs. Brittle behavior was observed in a narrow temperature window ca. 500 °C. In this latter case, AEs are consistently observed upon faulting and extremely sharp strain localization is observed in recovered samples. The resulting microstructures are consistent with the inverse ductile-to-brittle transition proposed by Proctor and Hirth (2016) in antigorite. This may therefore be a source of seismicity in subducting slabs at mantle pressures and temperatures from 500 to 600 °C. However, the acoustic signal observed here is orders of magnitude weaker than what is obtained at low PT conditions with brittle failure, consistently with low radiation efficiency of serpentinite faulting (Prieto et al., 2013) and suggests that other mechanisms are responsible for large intermediate-depth earthquakes. In fact, the present results are in line with a recent study (Ferrand et al., 2017), that suggests that intermediate earthquakes are likely induced by mechanical instabilities due to dehydration in partly hydrated peridotites.

58 GEOSCIENCES↗
DOE OSTI2017

Bayesian seismic inversion based on rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies

We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to post-stack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rock-physics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Modeling and simulation of earthquake soil structure interaction excited by inclined seismic waves

Presented is an application of wave potential formulation (WPF) together with domain reduction method (DRM) to modeling earthquake soil structure interaction (ESSI) behavior in horizontally layered ground under inclined incident seismic waves. Wave potential formulation is used to develop a spatially varying, inclined seismic wave field from incident Primary (P) and Secondary (S) waves that propagate through layered ground. Developed seismic wave field is then used to develop effective forces for Domain Reduction Method that are then used for analyzing ESSI response of a soil structure system. Developed methodology, called WPF-DRM, is verified using analytic solution for a free field response of layered ground subjected to inclined incident waves. Developed WPF-DRM methodology is illustrated through analysis of an ESSI response of a deeply embedded structure, a small modular reactor (SMR) subjected to incident S wave polarized in vertical plane (SV) with variation in inclinations and frequencies. Presented example highlights the influences of incident wave inclination and frequency on ESSI response of analyzed SMR.

58 GEOSCIENCES↗
DOE OSTI2021

Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy)

In seismic modelling, fault planes are normally assumed to be flat due to the lack of data which can constrain fault morphology. However, incorporating 3D fault morphology is important for modelling several phenomena, for example calculating mainshock induced stress changes. Here we utilize a data-analytical method to unveil the 3D rupture morphology of faults using unsupervised clustering techniques applied to earthquake hypocenters in seismic sequences. We apply this method to the 2009 L'Aquila seismic sequence which involved a M W 6.1 mainshock on April 6th. We use a dataset of about 50,000 relocated events, mostly microearthquakes, reaching magnitude of completeness equal to 0.7. Clustering distinguishes the earthquakes as occurring in three main clusters along with other minor fault segments. We then represent the morphology of the main Paganica fault system (responsible for the largest mainshock) using splines. This method shows promise as a step toward robustly and quickly obtaining 3D rupture morphologies where earthquake sequences have been monitored. The 3D model is presented interactively online, and the processing is presented in an interactive Jupyter Notebook (https://bit.ly/2MnCFdj).

58 GEOSCIENCES↗
DOE OSTI2019

Global Seismic Nowcasting With Shannon Information Entropy

Abstract Seismic nowcasting uses counts of small earthquakes as proxy data to estimate the current dynamical state of an earthquake fault system . The result is an earthquake potential score that characterizes the current state of progress of a defined geographic region through its nominal earthquake “cycle.” The count of small earthquakes since the last large earthquake is the natural time that has elapsed since the last large earthquake (Varotsos et al., 2006, https://doi.org/10.1103/PhysRevE.74.021123 ). In addition to natural time, earthquake sequences can also be analyzed using Shannon information entropy (“information”), an idea that was pioneered by Shannon (1948, https://doi.org/10.1002/j.1538‐7305.1948.tb01338.x ). As a first step to add seismic information entropy into the nowcasting method, we incorporate magnitude information into the natural time counts by using event self‐information. We find in this first application of seismic information entropy that the earthquake potential score values are similar to the values using only natural time. However, other characteristics of earthquake sequences, including the interevent time intervals, or the departure of higher magnitude events from the magnitude‐frequency scaling line, may contain additional information.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Relocation of Light and Moderate-Magnitude ( M 4-6) Seismicity Along the Central Mid-Atlantic

In this study, we estimate precise relative epicentroid locations of earthquakes along the Mid-Atlantic Ridge (MAR) between 12°N and 19°N using intermediate-period Rayleigh and Love wave cross-correlation measurements. We relocate 152 earthquakes with magnitudes ranging between 4.0 and 6.0 and when possible use P waveform modeling to estimate centroid depths. We compare surface-wave locations to locations in the U.S. Geological Survey, Global Centroid Moment Tensor, and Pacific Marine Environmental Laboratory hydroacoustic catalogs (Smith et al., 2002, https://doi.org/10.1029/2001GL013912) and demonstrate that surface-wave based relative location techniques provide significant improvement on precision. Our analyses shift locations to more tectonically consistent alignments that correlate strongly with bathymetric features. Relocations along the Fifteen-Twenty and Marathon Fracture Zones (FZs) move to within a few kilometers of the transform bathymetry. Events along ridge segments scatter over a distance of about 10–20 km perpendicular to the local plate boundary and most ridge-segment seismicity clusters spatially. North of the Fifteen-Twenty FZ, ridge events along several segments locate preferentially in the North American Plate, suggesting recent asymmetric seismic strain release and perhaps asymmetric spreading. Along the ridge between 13°N and 14°N, moderate-magnitude seismic strain release concentrates in an area with two active oceanic core complexes (OCC) along the western rift boundary, but the largest events along this segment appear to occur along the eastern boundary rift valley faults and represent deformation opposite the OCC detachment surface processes. Finally, precise moderate-magnitude earthquake epicentroid locations for events that occurred over the last few decades complement the results of temporary seismometer deployments that have imaged microseismic activity along the MAR.

58 GEOSCIENCES↗
DOE OSTI2019

The Seismic Response to Injected Carbon Dioxide: Comparing Observations to Estimates Based Upon Fluid Flow Modeling

Time-lapse seismic amplitude differences and travel time shifts, obtained while monitoring enhanced oil recovery at Cranfield, Mississippi, reveal coherent changes that are associated with the injection of carbon dioxide. Rock physics modeling highlights the importance of the oil, brine, and gas content of pore fluids prior to the injection of carbon dioxide. For example, compressional velocity changes due to the injection of carbon dioxide can drop from 300 m/s to less than 100 m/s as the percentage of oil increases from 1% to 50%. Predictions based upon a new technique for modeling wave propagation in a poroelastic medium containing an arbitrary number of fluids, coupled with multicomponent numerical reservoir modeling at Cranfield, reproduce the general pattern of observed seismic amplitude changes and travel time shifts. In particular, time-lapse amplitude changes suggest a significant and widespread lowering of compressional velocities due to the injection of CO 2 into an aquifer bounding the oil rim of the reservoir. It appears that the large-scale variations in preexisting pore fluid content have a major influence on seismic velocity changes, even in the highly heterogeneous reservoir at Cranfield.

58 GEOSCIENCES↗
DOE OSTI2020

Urban Seismic Site Characterization by Fiber-Optic Seismology

ccurate ground motion prediction requires detailed site effect assessment, but in urban areas where such assessments are most important, geotechnical surveys are difficult to perform, limiting their availability. Distributed acoustic sensing (DAS) offers an appealing alternative by repurposing existing fiber-optic cables, normally employed for telecommunication, as an array of seismic sensors. We present a proof-of-concept demonstration by using DAS to produce high-resolution maps of the shallow subsurface with the Stanford DAS array, California. We describe new methods and their assumptions to assess H/V spectral ratio—a technique widely used to estimate the natural frequency of the soil—and to extract Rayleigh wave dispersion curves from ambient seismic field. These measurements are jointly inverted to provide models of shallow seismic velocities and sediment thicknesses above bedrock in central campus. The good agreement with an independent survey validates the methodology and demonstrates the power of DAS for microzonation.

Geochemistry & Geophysics↗
DOE OSTI2021

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗
DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗
DOE OSTI2021

Spectral Element Modeling of Acoustic to Seismic Coupling Over Topography

Acoustic waves in the atmosphere are commonly recorded on seismometers as they couple into the ground. These signals, here called ground coupled airwaves, are not commonly considered in numerical modeling of infrasound propagation, which often assumes a rigid unmeshed boundary. Starting from an analytically-tractable spherical wave model, we analyze how the coupling of an acoustic wave into a planar elastic halfspace changes over a wide range of scenarios. We use energy admittance to quantify acoustic to seismic coupling over both a planar elastic halfspace and meshed topography. Our spectral element and analytic calculations have different maxima as a function of incidence angle, with very high admittance values for near-vertical incidence (maximum ≈78%). Energy admittance calculations at shallow incidence angles are much smaller (less than 1%). In simulations over the complex topography of Sakurajima Volcano, we attribute the variable spatial pattern of energy admittance to changes in earth parameters between each model. The observed pressure difference over the simulated 15 km region appears to be < 2%. While this value is relatively small, the cumulative addition over 100s of km and multiple acoustic bounce points may be significant. Acoustic to seismic coupling along the propagation path may bias long distance yield estimates, particularly when infrasound propagates over regions with steep topography or particularly slow seismic velocities, such as alluvial planes.

58 GEOSCIENCES↗
DOE OSTI2021

Seasonal Seismicity in the Lake Biwa Region of Central Japan Moderately Modulated by Lake Water Storage Changes

Seasonal seismicity in Japan has been observed and attributed to variations induced by hydrological loading, as well as atmospheric and pore-fluid pressure. We examine the seasonal variations of crustal (<15 km depth) and low magnitude (M < 4.5) earthquakes in central Japan near the Biwako-Seigan Fault Zone (BSFZ), and analyze their correlation with the annual hydrological cycle of nearby Lake Biwa, the largest freshwater lake in Japan. The 2002–2018 lake water storage is estimated using water-level gauge and satellite altimetry data. Regional surface mass loading from soil moisture and snow predicted by the Global Land Data Assimilation System is used to augment the lake mass variation. We compute the stress and stressing-rate changes induced by the surface mass loading from lake storage, snow, and soil moisture at seismogenic depth (10 km) along the BSFZ using a finite element model incorporating a 3D Earth structure. The results show that seasonal surface mass loading generates 1–5 kPa peak-to-peak Coulomb stress change and 2.5–15 kPa/year stressing-rate variations on the fault segments along the BSFZ. Two of the four fault segments present significant correlation between local seismicity rate and loading-derived stress and/or stressing rate changes. Among the hydrological load components, lake water, snow, and soil moisture exhibit different annual phases, and the lake water storage dominates the total annual stress changes. Our analysis indicates that seasonal surface hydrological loads, especially lake water variations, modulate the seismicity rate along the BSFZ in central Japan.

58 GEOSCIENCES↗
DOE OSTI2023

The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves

Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.

58 GEOSCIENCES↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2025

Analysis of Seismic Wave Amplification in Sedimentary Basins Using 3D Wavefield Simulations: Nenana Basin, Central Alaska

We investigate seismic wave amplification in Nenana basin, central Alaska, using 3D seismic wavefield simulations. We quantify frequency-dependent amplification by comparing synthetic seismograms among four different seismic velocity models: (a) a tomographic model, (b) a tomographic model with the uppermost 6.5 km replaced by a 1D basin profile, (c) a tomographic model with an embedded ellipsoidal basin, and (d) a tomographic model with an embedded realistic basin. For each model we perform wavefield simulations accurate up to 2 Hz for a set of 10 different earthquakes, which provide a range of frequencies and incident angles for waves interacting with the basin. Average amplification ratios are four on the horizontal components and seven on the vertical component. The dominant influence of the amplification is the 3D geometry of the basin, not the slow-velocity profile near the surface. Our synthetic amplification ratios are in general agreement with amplification estimates obtained from 14 stations that recorded the same set of earthquakes. Our approach offers a general strategy for documenting frequency-dependent basin amplification for a region with realistic basin structures and local earthquakes.

Geosciences↗
DOE OSTI2020

An exponential build-up in seismic energy suggests a months-long nucleation of slow slip in Cascadia

Slow slip events result from the spontaneous weakening of the subduction megathrust and bear strong resemblance to earthquakes, only slower. This resemblance allows us to study fundamental aspects of nucleation that remain elusive for classic, fast earthquakes. We rely on machine learning algorithms to infer slow slip timing from statistics of seismic waveforms. We find that patterns in seismic power follow the 14-month slow slip cycle in Cascadia, arguing in favor of the predictability of slow slip rupture. Here, we show that seismic power exponentially increases as the slowly slipping portion of the subduction zone approaches failure, a behavior that shares a striking similarity with the increase in acoustic power observed prior to laboratory slow slip events. Our results suggest that the nucleation phase of Cascadia slow slip events may last from several weeks up to several months.

58 GEOSCIENCES↗
DOE OSTI2024

Crustal permeability generated through microearthquakes is constrained by seismic moment

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability (Δk) are linearly related to the seismic moment (M) of individual MEQs as Δk ∝ M. This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs.

15 GEOTHERMAL ENERGY↗