Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic velocity”

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 OSTI2022

Utah FORGE - Development of a Reservoir Seismic Velocity Model and Seismic Resolution Study

This is data from and a final report on the development of a 3D velocity model for the larger FORGE area and on the seismic resolution in the stimulated fracture volume at the bottom of well 16A-32. The velocity model was developed using RMS velocities of the seismic reflection survey and seismic velocity logs from borehole measurements as an input model. To improve the accuracy of the model in the shallow subsurface, travel times phase arrivals of the direct propagating P-waves were determined from the seismic reflection data, using PhaseNet, a deep-neural-network-based seismic arrival time picking method. The travel times were subsequently inverted using the input velocity model. The seismic resolution study used borehole and surface seismic sensors as well as the seismicity observed during the April 2022 stimulation experiment to estimate the seismic resolution in the activated fracture reservoir. The data contain a 3D P- and S-wave velocity model for the larger FORGE area.

15 GEOTHERMAL ENERGY↗
DOE OSTI2019

Temporal changes of seismic velocities in the San Jacinto Fault zone associated with the 2016 M w 5.2 Borrego Springs earthquake

We study temporal changes of seismic velocities associated with the 10 June 2016 Mw 5.2 Borrego Springs earthquake in the San Jacinto fault zone, using nine component Green's function estimates reconstructed from daily cross correlations of ambient noise. The analysed data are recorded by stations in two dense linear arrays, at Dry Wash (DW) and Jackass Flat (JF), crossing the fault surface trace ~3 km northwest and southeast of the event epicentre. The two arrays have 9 and 12 stations each with instrument spacing of 25–100 m. Relative velocity changes (δv/v) are estimated from arrival time changes in the daily correlation coda waveforms compared to a reference stack. The obtained array-average δv/v time-series exhibit changes associated with the Borrego Springs event, superposed with seasonal variations. The earthquake-related changes are characterized by a rapid coseismic velocity drop followed by a gradual recovery. This is consistently observed at both arrays using time- and frequency-domain δv/v analyses with data from different components in various frequency bands. Larger changes at lower frequencies imply that the variations are not limited to the near surface material. A decreasing coseismic velocity reduction with coda wave lapse time indicates larger coseismic structural perturbations in the fault zone and near-fault environment compared to the surrounding rock. Observed larger changes at the DW array compared to the JF array possibly reflect the northwestward rupture directivity of the Borrego Springs earthquake.

Geochemistry & Geophysics↗
DOE OSTI2016

Anomalous density and elastic properties of basalt at high pressure: Reevaluating of the effect of melt fraction on seismic velocity in the Earth's crust and upper mantle

Independent measurements of the volumetric and elastic properties of Columbia River basalt glass were made up to 5.5 GPa by high-pressure X-ray microtomography and GHz-ultrasonic interferometry, respectively. The Columbia River basalt displays P and S wave velocity minima at 4.5 and 5 GPa, respectively, violating Birch’s law. These data constrain the pressure dependence of the density and elastic moduli at high pressure, which cannot be modeled through usual equations of state nor determined by stepwise integrating the bulk sound velocity as is common practice. We propose a systematic variation in compression behavior of silicate glasses that is dependent on the degree of polymerization and arises from the flexibility of the aluminosilicate network. Likewise, this behavior likely persists into the liquid state for basaltic melts resulting in weak pressure dependence for P wave velocities perhaps to depths of the transition zone. By modeling the effect of partial melt on P wave velocity reductions it is suggested that melt fraction determined by seismic velocity variations may be significantly overestimated in the crust and upper mantle.

58 GEOSCIENCES↗
DOE OSTI2024

Utah FORGE: Composite 3D Seismic Velocity Model

This is a composite 3D seismic velocity that was constructed from compiled information from several local studies regarding seismic velocities and structural information. This seismic velocity model is provided in NonLinLoc format (slow_len), which is readily usable in NonLinLoc software. Other model formats and versions of the model can be produced using the Python script provided with this data set. Details on how the model was created and prior velocity and structural information was used is provided in the accompanying documentation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

EGS Collab Experiment 1: 3D Seismic Velocity Model and Updated Microseismic Catalog from Double-Difference Seismic Tomography

This package contains a 3D Seismic velocity model and an updated microseismic catalog obtained for a double-difference seismic tomography study. The 3D_seismic_velocity_model text file contains x (m), y(m), z(m), P-wave velocity (km/s), P-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained), S-wave velocity (km/s), and S-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained). The Updated_MEQ_catalog text file contains event origin time, x(m), y(m), z(m), error in x (m), error in y (m), error in z (m), and RMS misfit (millisecond). The 3D_seismic_P-wave_velocity_model animation file shows slices of the 3D P-wave velocity model. The 3D_seismic_S-wave_velocity_model animation file shows slices of the 3D S-wave velocity model. The Interactive_MEQ_location_comparison API file is an interactive visualization of the updated microseismic event locations and the original seismic catalog. The visualization allows users to view and compare the event locations by dragging, rotating, and zooming in. An updated version of 3D_seismic_velocity_model and associated animations were included, which were calculated with a more strict assumption for quality indicators.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Seeking Repeating Anthropogenic Seismic Sources: Implications for Seismic Velocity Monitoring at Fault Zones

Abstract Seismic velocities in rocks are highly sensitive to changes in permanent deformation and fluid content. The temporal variation of seismic velocity during the preparation phase of earthquakes has been well documented in laboratories but rarely observed in nature. It has been recently found that some anthropogenic, high‐frequency (>1 Hz) seismic sources are powerful enough to generate body waves that travel down to a few kilometers and can be used to monitor fault zones at seismogenic depth. Anthropogenic seismic sources typically have fixed spatial distribution and provide new perspectives for velocity monitoring. In this work, we propose a systematic workflow to seek such powerful seismic sources in a rapid and straightforward manner. We tackle the problem from a statistical point of view, considering that persistent, powerful seismic sources yield highly coherent correlation functions (CFs) between pairs of seismic sensors. The algorithm is tested in California and Japan. Multiple sites close to fault zones show high‐frequency CFs stable for an extended period of time. These findings have great potential for monitoring fault zones, including the San Jacinto Fault and the Ridgecrest area in Southern California, Napa in Northern California, and faults in central Japan. However, extra steps, such as beamforming or polarization analysis, are required to determine the dominant seismic sources and study the source characteristics, which are crucial to interpreting the velocity monitoring results. Train tremors identified by the present approach have been successfully used for seismic velocity monitoring of the San Jacinto Fault in previous studies.

58 GEOSCIENCES↗
DOE OSTI2015

Factors Affecting Seismic Velocity in Alluvium

Yucca Flat at the Nevada National Security Site has been selected as the Source Physics Experiment (SPE) Dry Alluvium Geology Phase II site. The alluvium in this part of Yucca Flat is typical of desert basin fill, with discontinuous beds that are highly variable in clast size and provenance. Detailed understanding of the subsurface geology will be needed for interpretation of the SPE seismic data. A 3D seismic velocity model, created for Yucca Flat using interval seismic velocity data, shows variations in velocity within alluvium near the SPE Phase II site beyond the usual gradual increase of density with depth due to compaction. In this study we examined borehole lithologic logs, geophysical logs, downhole videos, and laboratory analyses of sidewall core samples to understand which characteristics of the alluvium are related to these variations in seismic velocity. Seismic velocity of alluvium is generally related to its density, which can be affected by sediment provenance, clast size, gravel percentage, and matrix properties, in addition to compaction. This study presents a preliminary subdivision of the alluvial strata in the SPE Phase II area into mappable units expected to be significant to seismic modeling. Further refinements of the alluvial units may be possible when seismic data are obtained from SPE Phase II tests.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring Seismic Velocity Changes Across the San Jacinto Fault Using Train‐Generated Seismic Tremors

Abstract Microseismic noise has been used for seismic velocity monitoring. However, such signals are dominated by low‐frequency surface waves that are not ideal for detecting changes associated with small tectonic processes. Here we show that it is possible to extract stable, high‐frequency body waves using seismic tremors generated by freight trains. Such body waves allow us to focus on small velocity perturbations in the crust with high spatial resolution. We report on 10 years of seismic velocity temporal changes at the San Jacinto Fault. We observe and map a two‐month‐long episode of velocity changes with complex spatial distribution and interpret the velocity perturbation as produced by a previously undocumented slow‐slip event. We verify the hypothesis through numerical simulations and locate this event along a fault segment believed to be locked. Such a slow‐slip event stresses its surroundings and may trigger a major earthquake on a fault section approaching failure.

58 GEOSCIENCES↗
DOE OSTI2022

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗
DOE OSTI2021

Regional seismic velocity changes following the 2019 M w 7.1 Ridgecrest, California earthquake from autocorrelations and P / S converted waves

SUMMARY We examine regional transient changes of seismic velocities generated by the 2019 Mw 7.1 Ridgecrest earthquake in California, using autocorrelations of moving time windows in continuous waveforms recorded at regional stations. We focus on traveltime differences in a prominent phase generated by an interface around 2 km depth, associated with transmitted Pp waves and converted Ps waves from the ongoing microseismicity. Synthetic tests demonstrate the feasibility of the method for monitoring seismic velocity changes. Taking advantage of the numerous aftershocks in the early period following the main shock, we obtain a temporal resolution of velocity changes up to 20 min in the early post-main-shock period. The results reveal regional coseismic velocity drops in the top 1–3 km with an average value of ∼2 per cent over distances up to 100 km from the Ridgecrest event. These average velocity drops are likely dominated by larger changes in the shallow materials and are followed by rapid recoveries on timescales of days. Around the north end of the Ridgecrest rupture and the nearby Coso geothermal region, the observed coseismic velocity drops are up to ∼8 per cent. The method allows monitoring temporal changes of seismic velocities with high temporal resolution, fast computation and precise spatial mapping of changes. The results suggest that significant temporal changes of seismic velocities of shallow materials are commonly generated on a regional scale by large events.

58 GEOSCIENCES↗
DOE OSTI2020

EGS Collab Experiment 1: 3D Seismic Velocity Model and Updated Microseismic Catalog Using Transfer-Learning Aided Double-Difference Tomography

This package contains a 3D Seismic velocity model and an updated microseismic catalog associated with a proceedings paper (Chai et al., 2020) published in the 45th Workshop on Geothermal Reservoir Engineering. The 3D_seismic_velocity_model text file contains x (m), y(m), z(m), P-wave velocity (km/s), P-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained), S-wave velocity (km/s), and S-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained). The Updated_MEQ_catalog text file contains event origin time, x(m), y(m), z(m), error in x (m), error in y (m), error in z (m), and RMS misfit (millisecond). The 3D_seismic_P-wave_velocity_model animation file shows slices of the 3D P-wave velocity model. The 3D_seismic_S-wave_velocity_model animation file shows slices of the 3D S-wave velocity model. The Interactive_MEQ_locations API file is an interactive visualization of the updated microseismic event locations. The visualization allows users to view the event locations by dragging, rotating, and zooming in. References: Chai, C., Maceira, M., Santos-Villalobos, H. J., Venkatakrishnan, S. V., Schoenball, M., and EGS Collab Team, 2020, Automatic Seismic Phase Picking Using Deep Learning for the EGS Collab Project, in PROCEEDINGS, 45th Workshop on Geothermal Reservoir Engineering, edited, Stanford University, Stanford, California, 45, 1266-1276.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗
DOE OSTI2016

Seismic velocity changes associated with aseismic deformations of a fault stimulated by fluid injection

Fluid pressure plays an important role in the stability of tectonic faults. However, the in situ mechanical response of faults to fluid pressure variations is still poorly known. To address this question, we performed here a fluid injection experiment in a fault zone in shales while monitoring fault movements at the injection source and seismic velocity variations from a near-distance (<10 m) monitoring network. We measured and located the P and S wave velocity perturbations in and around the fault using repetitive active sources. We observed that seismic velocity perturbations dramatically increase above 1.5 MPa of injection pressure. This is consistent with an increase of fluid flow associated with an aseismic dilatant shearing of the fault as shown by numerical modeling. We find that seismic velocity changes are sensitive to both fault opening by fluid invasion and effective stress variations and can be an efficient measurement for monitoring fluid-driven aseismic deformations of faults.

58 GEOSCIENCES↗
DOE OSTI2019

Science Plan for Improving 3-D Seismic Velocity Models in the San Francisco Bay Region, 2019-24

This five-year science plan outlines short-term and long-term goals for improving 3-D seismic velocity models in the greater San Francisco Bay region as well as how to foster a community effort in reaching those goals. The short-term goals focus on improving the current USGS San Francisco Bay region geologic and seismic velocity model using existing data. The long-term goals focus on acquiring new data and leveraging better analysis tools to improve the model and characterize the uncertainty. The plan describes opportunities for contributions by members of the community to development of these seismic velocity models, provides current and potential users with general information on where efforts will likely be focused to improve these models and how new versions will be released, and outlines funding needs and obstacles for improving and maintaining such models. Several aspects of this plan, including how to foster a community effort, are independent of the geographic region and apply to other, similar efforts.

58 GEOSCIENCES↗
DOE OSTI2015

Imaging of 3-D seismic velocity structure of Southern Sumatra region using double difference tomographic method

Southern Sumatra region has a high level of seismicity due to the influence of the subduction system, Sumatra fault, Mentawai fault and stretching zone activities. The seismic activities of Southern Sumatra region are recorded by Meteorological Climatological and Geophysical Agency (MCGA’s) Seismograph network. In this study, we used earthquake data catalog compiled by MCGA for 3013 events from 10 seismic stations around Southern Sumatra region for time periods of April 2009 – April 2014 in order to invert for the 3-D seismic velocities structure (Vp, Vs, and Vp/Vs ratio). We applied double-difference seismic tomography method (tomoDD) to determine Vp, Vs and Vp/Vs ratio with hypocenter adjustment. For the inversion procedure, we started from the initial 1-D seismic velocity model of AK135 and constant Vp/Vs of 1.73. The synthetic travel time from source to receiver was calculated using ray pseudo-bending technique, while the main tomographic inversion was applied using LSQR method. The resolution model was evaluated using checkerboard test and Derivative Weigh Sum (DWS). Our preliminary results show low Vp and Vs anomalies region along Bukit Barisan which is may be associated with weak zone of Sumatran fault and migration of partial melted material. Low velocity anomalies at 30-50 km depth in the fore arc region may indicated the hydrous material circulation because the slab dehydration. We detected low seismic seismicity in the fore arc region that may be indicated as seismic gap. It is coincides contact zone of high and low velocity anomalies. And two large earthquakes (Jambi and Mentawai) also occurred at the contact of contrast velocity.

58 GEOSCIENCES↗
NASA NTRS1977

The distribution of seismic velocities and attenuation in the earth

Estimates of the radial distribution of seismic velocities and density and of seismic attenuation within the earth are obtained through inversion of body wave, surface wave, and normal mode data. The effect of attenuation related dispersion on gross earth structure, and on the reliability of eigenperiod identifications is discussed. The travel time baseline discrepancies between body waves and free oscillation models are examined and largely resolved.

Hart, R. S.↗
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↗
NASA NTRS1983

Seismic velocity structure of the lunar mantle

The recently completed set of seismic arrival times from the Apollo lunar seismic network are inverted to estimate the average seismic velocities in three sections of the lunar mantle: two for the upper mantle and one for the middle mantle. The method used is a variation of the linearized least squares inversion where the inversion is accomplished in steps. The estimated average velocities in the upper mantle decrease from Vp = 7.74 km/sec and Vs = 4.49 km/sec in the section above 270-km depth to Vp = 7.46 km/sec and Vs = 4.25 km/sec in the section between 270- and 500-km depth, confirming the earlier finding of negative gradients based on seismic amplitude variations. The average velocities in the middle mantle between the depths of 500 km and 1000 km of Vp = 8.26 km/sec and Vs = 4.65 km/sec are significantly higher than those in the upper mantle, contradicting earlier estimates based on more limited data. The higher velocities may suggest initial melting of the moon down to at least 1000-km depth.

Nakamura, Y.↗