Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic wave”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗
DOE OSTI2021

Toward Robust and Routine Determination of Mw for Small Earthquakes: Application to the 2020 Mw 5.7 Magna, Utah, Seismic Sequence

To better characterize seismic hazard, particularly, for induced seismicity, there is an increasing interest in methods to estimate moment magnitude (M w ) for small earthquakes. M w is generally preferred over other magnitude types, but, it is difficult to estimate M w for earthquakes with local magnitude (M L ) <3-3.5, using conventional moment tensor (MT) inversion. The 2020 M ww 5.7 Magna, Utah, seismic sequence provides an opportunity to illustrate and evaluate the value of spectral methods for this purpose. Starting with a high-quality seismic catalog of 2103 earthquakes (M L <5.6), we estimate M w using two independent spectral methods—one based on direct waves, yielding M w,direct , and the other based on coda waves, yielding M w,coda . For the direct-wave method, we present a non-parametric (NP) inversion scheme that solves for apparent geometrical spreading, G(R), and site effects (S), similar to other NP procedures that have been used to calibrate regional M L scales. The NP inversion is constrained using M ws derived from MTs for nine events in the Magna sequence. We recover statistically robust and physically reasonable G(R) and S and compute M w,direct for 635 Magna earthquakes down to M L 0.7. For the coda-wave method, we consider two separate calibration schemes involving previous MT solutions and compute M w,coda for 311 earthquakes down to M L 1.0. For 280 of the events that were processed with both methods—M w,direct and M w,coda —are strongly correlated (r = 0.98), with a mean difference of only 0.05. We compare M w,direct and M w,coda with M L and find reasonably good agreement for M L <3.6 with the theoretically predicted relationship of M w =(2/3)M L +C, in which C is a regional constant. Our results imply that seismic network operators can use spectral-based M w estimates to replace M L estimates for events with M L ≥1.0, and possibly smaller. The main requirement is the existence of a small number of MT solutions for calibration purposes.

58 GEOSCIENCES↗
DOE OSTI2023

Active‐Source Seismic Imaging of Fault Re‐Activation and Leakage: An Injection Experiment at the Mt Terri Rock Laboratory, Switzerland

Abstract We conducted a time‐lapse seismic experiment utilizing automated active seismic source and sensor arrays to monitor a reactivated fault within the Opalinus clay formation at the Mont Terri Rock Laboratory (Switzerland), an analog caprock for geologic carbon storage. A series of six brine injections were conducted into the so‐called Main Fault to reactivate it. Seismic instrumentation in five monitoring boreholes on either side of the fault was used to continuously probe changes in P‐wave travel‐times associated with fault displacement and leakage. We performed time‐lapse travel‐time tomography on five hundred sequential data sets; this revealed a zone of decreased P‐wave velocity, up to 16 m/s, during each injection cycle, followed by a velocity increase during shut‐in. These results demonstrate varying elastic property perturbations, both spatially and temporally, along the fault plane during reactivation. We then interpreted these velocity changes in terms of fault dilation induced by pressurized fluids along the fault.

58 GEOSCIENCES↗
DOE OSTI2022

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

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

Advancing Tunnel Detection Via Vertical Acoustic Profiling

This study demonstrates that placing an acoustic source at depth significantly enhances the detection of tunnels. By adapting the methodology of vertical seismic profiling to acoustic-range frequencies, we introduce a novel approach for tunnel detection. To our knowledge, this report documents the first application of this technique. The field study was conducted on a cleared site at the Oak Ridge National Laboratory; the subsurface was clay. A series of vertical shafts were bored immediately before and after the placement of a 40 ft long, 3 ft diameter steel pipe located 10 ft below the surface. Acoustic signatures were produced at various depths up to 30 ft and an array of surface geophones was used to collect the resulting signals. The innovation lies in using subsurface acoustic signatures to better enable the collection of wide-angle diffracted waves by surface geophones. These diffracted waves appear as subharmonic frequencies of the source signal and are recorded by a surface geophone array. While the concept draws from vertical seismic profiling used in oil and gas exploration, a key distinguishing feature is the use of acoustic frequencies, which are necessary for detecting small subsurface features. The findings have important implications. Critical information pertaining to diffracted source signals are not measured by current state-of-the-art methods. Although this study confirms enhanced detectability, accurate localization and imaging will require additional research. Future work should focus on leveraging diffracted-wave strength and time delays to enable full imaging and characterization.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

STILGAR: Subsurface Models for Graymont Pleasant Gap Mine

The detection, location, and monitoring of underground structures are of great importance to national and global security. Tunnels and voids generate seismic signatures detectable at the surface, but using non-invasive seismic data to image near-surface presents several challenges in real-world applications. In this report, we describe the use of a dense surface seismic deployment to generate subsurface models of the Graymont Pleasant Gap mine - a single-layer mine with a complex structure embedded in a high-velocity P-wave limestone bedrock. Our approach consists of three key methods. We use P-wave arrival times from local blast events to perform a tomography inversion with the tomoTD method, constructing a P-wave velocity model of the subsurface. We model the layer above the mine using Rayleigh wave ellipticity and inversion techniques. We leverage ongoing anthropogenic activities to identify and locate noise sources both on the surface and within the subsurface. With this integrated approach we aim to overcome the challenges and enhance our ability to non-invasively characterize underground structures, contributing to improved seismic monitoring techniques.

58 GEOSCIENCES↗
DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗
DOE OSTI2021

The Water-Fe-Pressure dependent single-crystal elastic properties of wadsleyite: Implications for the seismic anisotropy in the upper Mantle Transition Zone

Recent seismic studies suggest an anisotropic Mantle Transition Zone (MTZ) in areas adjacent to subducted slabs. Wadsleyite is the main anisotropy contributor in the upper MTZ, therefore the interpretation of these seismic observations requires the knowledge of single-crystal elastic moduli (C ij s) and the deformation-induced lattice preferred orientation (LPO) of wadsleyite. Wadsleyite can host up to 3 wt% water in its crystal structure as point defects in the form of hydroxyl groups, however, the combined effect of water content, Fe content, and pressure on the C ij s of wadsleyite remains unclear. In this study, we measured the high-pressure single-crystal elasticity of a synthetic hydrous Fe-bearing wadsleyite (0.14(5) wt% water, Fe#=9.4, Fe 3+ /ΣFe=0.3) up to 18.2 (2) GPa. In combination with previous experimental data, we separated the effects of pressure, water, and Fe contents on the Cijs and intrinsic elastic anisotropy of wadsleyite. Overall, our results suggest that the intrinsic elastic anisotropy of wadsleyite decreases with pressure, water, and Fe contents. At 15 GPa, increasing the water content by 0.1 wt% or Fe# by 1 decreases the V P and V S anisotropy of wadsleyite by ~1.1-1.3%, and ~0.8-1.3% in average, respectively. Combined with the LPO determined in previous deformation experiments, we modeled the seismic anisotropy in the upper MTZ generated by a sub-vertical mantle flow near cold subducted slabs and a sub-horizontal mantle flow in the ambient mantle. In both scenarios, the LPO of wadsleyite leads to V SV (vertically polarized shear wave velocity) > V SH (horizontally polarized shear wave velocity). Our results suggest that wadsleyite may account for a weak anisotropic MTZ (<1%) on the global scale. Considering the fact that water decreases the elastic anisotropy but promotes LPO of wadsleyite, seismic anisotropy may not be a good water sensor.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Internal Structure of the Central Garlock Fault Zone From Ridgecrest Aftershocks Recorded by Dense Linear Seismic Arrays

Abstract We provide high‐resolution seismic imaging of the central Garlock fault using data recorded by two dense seismic arrays that cross the Ridgecrest rupture zone (B4) and the Garlock fault (A5). Analyses of fault zone head waves and P ‐wave delay times at array A5 show that the Garlock fault is a sharp bimaterial interface with P waves traveling ∼5% faster in the northern crustal block. The across‐fault velocity contrast agrees with regional tomography models and generates clear P ‐wave reflections in waveforms recorded by array B4. Kirchhoff migration of the reflected waves indicates a near‐vertical fault between 2 and 6 km depth. The P ‐wave delay times imply a ∼300‐m‐wide transition zone near the Garlock fault surface trace beneath array A5, offset to the side with faster velocities. The results provide important constraints for derivations of earthquake properties, simulations of ruptures and ground motion, and future imaging studies associated with the Garlock fault.

58 GEOSCIENCES↗
DOE OSTI2021

Surface-wave tomography of the Emeishan large igneous province (China): Magma storage system, hidden hotspot track, and its impact on the Capitanian mass extinction

Large igneous provinces (LIPs) are commonly associated with mass extinctions. However, the precise relations between LIPs and their impacts on biodiversity is enigmatic, given that they can be asynchronous. It has been proposed that the environmental impacts are primarily related to sill emplacement. Therefore, the structure of LIPs' magma storage system is critical because it dictates the occurrence and timing of mass extinction. We use surface-wave tomography to image the lithosphere under the Permian Emeishan large igneous province (ELIP) in southwestern China. We find a northeast-trending zone of high shear-wave velocity (Vs) and negative radial anisotropy (Vsv > Vsh; v and h are vertically and horizontally polarized S waves, respectively) in the crust and lithosphere. We rule out the possibilities of rifting or orogenesis to explain these seismic characteristics and interpret the seismic anomaly as a mafic-ultramafic, dike-dominated magma storage system of the ELIP. We further propose that the anomaly represents a hidden hotspot track that was emplaced before the ELIP eruption. A zone of higher velocity but less-negative radial anisotropy, on the hotspot track but to the northeast of the eruption center in the Panxi region, reflects an elevated proportion of sills emplaced at the incipient stage of the ELIP. Liberation of poisonous gases by the early sill intrusions explains why the mid-Capitanian global biota crisis preceded the peak ELIP eruption by 2–3 m.y.

58 GEOSCIENCES↗
DOE OSTI2023

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗
DOE OSTI2024

A Pn magnitude scale m b ( Pn ) for earthquakes along the equatorial Mid-Atlantic Ridge

We developed a short-period Pn magnitude scale m b (Pn) for earthquakes along the equatorial Mid-Atlantic Ridge. Due to low signal-to-noise ratios, teleseismic body wave magnitude and long-period surface wave magnitude cannot be confidently determined for small earthquakes of m b < 4. Local magnitude scales are also not useful for these events because the oceanic environment does not allow the propagation of crustal phases. However, regional high-frequency Pn waves from these small- to moderate-size (m b 3–6) earthquakes are well recorded in the equatorial Atlantic region and can be used to assign magnitudes. We measured over 2041 Pn peak amplitudes on vertical records from about 21 stations in northeastern Brazil and 11 stations in western Africa in the distance range of 700–3700 km. We analysed data from 189 events from the global centroid moment tensor catalogue to tie our m b (Pn) scale to M w so that seismic moments can be readily estimated. Pn arrivals show apparent group velocity between 7.9 km s –1 at short ranges (~1000 km) and up to 9.1 km s –1 at 3500 km. The measured peak amplitudes have a frequency between 0.8 and 3 Hz at 1000–1800 km, but at greater distances, 1800–3700 km, they show a remarkably consistent frequency of about 0.8 Hz. The peak amplitude attenuates at a higher rate at short distances (~0.65 magnitude units between 700 and 2000 km) but attenuates at a lower rate at long distances (~0.35 magnitude units between 2000 and 3700 km). The low rate of amplitude decay with distance and nearly constant frequency content of the peak amplitudes suggest that Pn waves propagate efficiently in the lower part of the upper mantle in the equatorial Atlantic Ocean basins. Furthermore, these are important attributes of oceanic Pn waves that can be used to assign magnitude for small- to moderate-size earthquakes in the equatorial mid-Atlantic region. The estimated station corrections correlate well with upper mantle low-velocity anomalies, especially in Brazil.

58 GEOSCIENCES↗
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2021

Hydraulic fracture mapping using wavelet-based fusion of wave transmission and emission measurements

Fracture characterization is essential to hydrocarbon and geothermal exploration and production as well as carbon geo-sequestration. Active and passive measurements have been extensively applied to map fractures across various scales. Active and passive methods have inherent advantages and limitations that complement each other. Here, in this study, to improve the mapping of embedded fractures and the surrounding geomechanically altered regions, we integrate active shear-wave transmission measurement with passive acoustic-emission measurement collected during a lab-scale hydraulic-fracturing experiment. The proposed data fusion leverages the advantages of the two modalities of measurements, while minimizing their individual limitations. Polarized shear-wave transmission (active) measurement was collected before and after the hydraulic fracturing, and the acoustic emission (passive) measurement was collected during hydraulic fracturing. Two sets of two-dimensional maps of fracture and fracture-induced damage in axial, median and frontal planes were obtained by separately processing active and passive measurements. These 2D maps are then fed to wavelet-based image fusion technique to integrate the two sources of information for the reliable mapping/imaging of the embedded fractures and the surrounding geomechanically altered regions.

58 GEOSCIENCES↗
DOE OSTI2025

Microseismicity Modulation Due To Changes in Geothermal Production at San Emidio, Nevada, USA

Brief cessations of geothermal production can induce seismicity, a phenomenon that has drawn increasing attention in recent years. Such observations are rare, and the underlying mechanism requires careful analysis. In April 2022, a dense seismic and hydrologic monitoring system was deployed at the San Emidio geothermal field, Nevada, to accompany a planned power plant shutdown. Using the dense seismic array data, we detected and located ∼1,800 microseismic events (MSEs) and developed a high-resolution tomographic P-wave velocity model. We observed substantially increased microseismicity during shutdown. Most MSEs occurred on pre-existing normal faults, which are contained within extremely low-velocity zones that are likely damaged, fluid-filled, and hydraulically connected to nearby production wells. Hydrologic data show rapid fluid pressure increases of <60 kPa following the shutdown. We suggest that the cessation of production rapidly increased fluid pressures along pre-existing fault zones, activating critically stressed fault patches and fractures and producing microseismicity.

Guo, Hao [University of Wisconsin-Madison, WI (Uni↗
DOE OSTI2023

Infrasound Detections of Low-Magnitude Earthquakes: Preliminary Results of the West Texas Acoustic Experiment

Infrasound observations have grown increasingly important for the monitoring of earthquakes. While large earthquakes generate infrasound that can be detected thousands of kilometers away, there are few near-field observations of infrasound generated by low-magnitude events. We describe preliminary results of the West Texas Acoustic Experiment, during which infrasound sensors collected continuous data in the Permian Basin for a six-month period spanning January—June 2023. During this time, more than 1000 earthquakes with magnitudes between 1.2 and 4.2 occurred within 50 km of the network. We used spectral analysis, array processing, and manual inspection of waveforms to evaluate arrivals of infrasound signals following 84 events with magnitudes between 2.5 and 4.2. Here, we describe eight such events and the infrasound signals associated with each. We find detections of seismic-to-acoustic infrasound signals associated with seven events. We also find strong evidence of a laterally-propagating, purely acoustic wave generated by an M2.9 earthquake.

58 GEOSCIENCES↗