Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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 OSTI2025

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗
DOE OSTI2013

Recovery Act: High-Temperature Circuit Boards for use in Geothermal Well Monitoring Applications

The U.S. Department of Energy is leading the development of alternative energy sources that will ensure the long-term energy independence of our nation. One of the key renewable resources currently being advanced is geothermal energy. To tap into the large potential offered by generating power from the heat of the earth, and for geothermal energy to be more widely used, it will be necessary to drill deeper wells to reach the hot, dry rock located up to 10 km beneath the earth’s surface. In this instance, water will be introduced into the well to create a geothermal reservoir. A geothermal well produced in this manner is referred to as an enhanced geothermal system (EGS). EGS reservoirs are typically at depths of 3 to 10 km, and the temperatures at these depths have become a limiting factor in the application of existing downhole technologies. These high temperatures are especially problematic for electronic systems such as downhole data-logging tools, which are used to map and characterize the fractures and high-permeability regions in underground formations. Information provided by these tools is assessed so that underground formations capable of providing geothermal energy can be identified, and the subsequent drilling operations can be accurately directed to those locations. The mapping of geothermal resources involves the design and fabrication of sensor packages, including the electronic control modules, to quantify downhole conditions (300°C temperature, high pressure, seismic activity, etc.). Because of the extreme depths at which these measurements are performed, it is most desirable to perform the sensor signal processing downhole and then transmit the information to the surface. This approach necessitates the use of high-temperature electronics that can operate in the downhole environment. Downhole signal processing in EGS wells will require the development and demonstration of circuit boards that can withstand the elevated temperatures found at these depths. At present, the highest-temperature commercially available circuit boards are based on polyimide materials, and those have maximum use temperatures of 200 to 250°C. In addition to thermal stability, downhole electronics must also be fabricated into high-aspect-ratio packages. For example, the multilayer assemblies produced at SNL were approximately 2.5 cm wide and 50 cm long. Because of this very high form factor, glass-fiber-reinforced polymers are much more desirable than multilayer ceramic modules (MCM). MCMs have many advantages for some applications, but are susceptible to damage induced by the mechanical and vibrational loads commonly experienced by data-logging tools. Thus, as EGS technology continues to advance, there is a strong need for multilayer electronics that can provide the necessary thermal performance while also being compatible with high-form-factor circuit designs. This project involved the design and development of high-temperature circuit materials, as well as the fabrication and testing of circuit components. The material development included the evaluation of various polymer/fiberglass composites, whereas the circuit components were tested using conventional microelectronic evaluation techniques. This effort targeted development of a new class of high-temperature multilayer circuit boards for use in downhole data-logging applications where temperatures are on the order of 300°C. This is consistent with DOE’s multiyear plan for advancing technologies for use in enhanced geothermal systems. Organic and inorganic polymer systems, both with glass reinforcements, were considered to provide the following performance at elevated temperatures: • Mechanical strength and durability • High dielectric strength and electrical resistivity • Thermal stability • Strong adhesion to copper to ensure the reliability of the multilayer assemblies • Processing characteristics that are consistent with state-of-the-art multilayer circuit board manufacturing practices

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

Sound velocity of Fe-S liquids at high pressure: Implications for the Moon's molten outer core

Sound velocities of Fe and Fe–S liquids were determined by combining the ultrasonic measurements and synchrotron X-ray techniques under high pressure–temperature conditions from 1 to 8 GPa and 1573 K to 1973 K. Four different liquid compositions were studied including Fe, Fe–10 wt% S, Fe–20 wt% S, and Fe–27 wt% S. Our data show that the velocity of Fe-rich liquids increases upon compression and decreases with increasing sulfur content, whereas temperature has negligible effect on the velocity of Fe–S liquids. The sound velocity data were combined with ambient-pressure densities to fit the Murnaghan equation of state (EOS). Compared to the lunar seismic model, our velocity data constrain the sulfur content at 4±3 wt%, indicating a significantly denser (6.5±0.5 g/cm 3 ) and hotter (1870-70+100 K) outer core than previously estimated. A new lunar structure model incorporating available geophysical observations points to a smaller core radius. Our model suggests a top–down solidification scenario for the evolution of the lunar core. Such “iron snow” process may have been an important mechanism for the growth of the inner core.

Jing, Zhicheng↗
DOE OSTI2022

Wintertime Brine Discharge at the Surface of a Cold Polar Glacier and the Unexpected Absence of Associated Seismicity

A subglacial groundwater system beneath Taylor Glacier, Antarctica, discharges hypersaline, iron-rich brine episodically at the glacier surface to create Blood Falls. However, the triggering mechanism for these brine release events is not yet understood. Identifying which fracture processes are observed seismically can help us better characterize the hydrological system at Taylor Glacier, and more generally, provide us with a broader understanding of englacial hydrologic activity in cold glaciers. We document wintertime brine discharge using time-lapse photography. Subfreezing air temperatures during the brine discharge indicate that surface melt-induced hydrofracture is an unlikely trigger for brine release. Further, we analyze local seismic data to test a hypothesis that fracturing generates elevated surface wave energy preceding and/or coinciding with brine release events. Our results show no discernible elevated Rayleigh wave activity prior to or during Blood Falls brine release. Instead, we find a pattern of seismic events dominated by a seasonal signal, with more Rayleigh events occurring in the summer than the winter from the Blood Falls source area. We calculate that the volumetric opening of cracks that would generate Rayleigh waves at our detection limits are of similar size to myriad cracks in glacier ice, lake ice, and frozen sediment in the terminus area. We therefore propose that any fracturing coincident with brine release activity likely consists of a series of smaller opening events that are masked by other seismicity in the local environment.

58 GEOSCIENCES↗
DOE OSTI2018

Microseismicity Analysis for HFTS Pad and Correlation With Completion Parameters

Microseismic data can be used as a tool to help understand fracturing behavior during hydraulic stimulation. Recent studies have validated complex fracture growth and interaction during fracturing process and there is a need to better utilize microseismic data as a way of improving our understanding of this complexity. Moreover, significant stress variations along the laterals and various formations of interest can be better interpreted by making use of b-value distribution as a proxy measure for stress. In this study, we tie the Frequency Magnitude Distributions (FMD) from microseismic with the 3D seismic and petrophysical data to understand stress variations within the upper and middle Wolfcamp formations. We also look at how temporal microseismic FMD variance can be used as a tool to help understand fracturing behavior with time. Our results highlight the utility of microseismic FMD as a valuable hydraulic fracturing diagnostic tool post fracturing operations and as a validation for observations made during treatment including relevant treatment parameters.

Maity, Debotyam↗
DOE OSTI2014

Advanced geophysical underground coal gasification monitoring

Underground Coal Gasification (UCG) produces less surface impact, atmospheric pollutants and greenhouse gas than traditional surface mining and combustion. Therefore, it may be useful in mitigating global change caused by anthropogenic activities. Careful monitoring of the UCG process is essential in minimizing environmental impact. In this paper, we first summarize monitoring methods that have been used in previous UCG field trials. We then discuss in more detail a number of promising advanced geophysical techniques. These methods – seismic, electromagnetic, and remote sensing techniques – may provide improved and cost-effective ways to image both the subsurface cavity growth and surface subsidence effects. Active and passive seismic data have the promise to monitor the burn front, cavity growth, and observe cavity collapse events. Electrical resistance tomography (ERT) produces near real time tomographic images autonomously, monitors the burn front and images the cavity using low-cost sensors, typically running within boreholes. Interferometric synthetic aperture radar (InSAR) is a remote sensing technique that has the capability to monitor surface subsidence over the wide area of a commercial-scale UCG operation at a low cost. It may be possible to infer cavity geometry from InSAR (or other surface topography) data using geomechanical modeling. The expected signals from these monitoring methods are described along with interpretive modeling for typical UCG cavities. They are illustrated using field results from UCG trials and other relevant subsurface operations.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Rapid Refreezing of a Marginal Ice Zone Across a Seafloor Distributed Acoustic Sensor

Abstract Seismic instruments such as broadband seismometers and distributed acoustic sensors (DAS) have a demonstrated potential for wide‐scale and continuous in situ monitoring of near‐surface environmental and anthropogenic processes. DAS is attractive for development as a multi‐geophysical observatory due to the prevalence of existing fiber infrastructure in regions with environmental, cultural, or strategic significance. Here we present results from a multi‐seasonal acquisition of DAS data on a seafloor cable in the Beaufort Sea, Alaska. During a November 2021 data collect we captured the rapid transition of ambient noise characteristics from an “ice‐free” state to an “ice‐bound” state. A sea ice formation front was plainly visible on the DAS record and was observed to propagate 20 km seaward over a period of 8 hr. Satellite‐based instrumentation were unable to record this event due to cloud cover, low light conditions, and orbital frequency.

Baker, Michael G.↗
DOE OSTI2023

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗
DOE OSTI2016

Evaluation and correction of timing inconsistencies from the Thor 1 and Thor 2 experiments at the Nevada National Security Site.

Active source seismic data was collected at the Nevada National Security Site using the Seismic Hammer(TM) (SH), under contract from HK Exploration. The SH generates a seismic pulse by dropping a 13 metric ton mass from a height of 1.5 m. Post-survey evaluation of collected data revealed inconsistencies in shot trigger time that required additional analysis and correction using cross-correlation and/or time shifts derived from manual picks of trigger times. While the primary analysis for which this data set was collected is independent of the knowledge of shot trigger time, other processing methods require highly precise knowledge of the trigger time. In order to make the Thor data set more usable to the larger community, additional work was undertaken. Results using the preferred method of cross-correlation were found to be satisfactory. An improved timing fiducial approach is required to reduce timing errors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

Comparison of emission properties of two homologous flares in AR 11283

Large, complex, active regions may produce multiple flares within a certain period of one or two days. These flares could occur in the same location with similar morphologies, commonly referred to as 'homologous flares'. In 2011 September, active region NOAA 11283 produced a pair of homologous flares on the 6th and 7th, respectively. Both of them were white-light (WL) flares, as captured by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory in visible continuum at 6173 Å which is believed to originate from the deep solar atmosphere. We investigate the WL emission of these X-class flares with HMI's seeing-free imaging spectroscopy. The durations of impulsive peaks in the continuum are about 4 minutes. We compare the WL with hard X-ray (HXR) observations for the September 6 flare and find a good correlation between the continuum and HXR both spatially and temporally. In absence of RHESSI data during the second flare on September 7, the derivative of the GOES soft X-ray is used and also found to be well correlated temporally with the continuum. We measure the contrast enhancements, characteristic sizes, and HXR fluxes of the twin flares, which are similar for both flares, indicating analogous triggering and heating processes. However, the September 7 flare was associated with conspicuous sunquake signals whereas no seismic wave was detected during the flare on September 6. Therefore, this comparison suggests that the particle bombardment may not play a dominant role in producing the sunquake events studied in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An open-access simulated earthquake ground-motion database for an M7 Hayward Fault earthquake in the San Francisco Bay Region

Comprehensive understanding of earthquake ground motions, particularly in the near-fault region of large-magnitude events, is limited by gaps in strong-motion data. This challenge is prominent in areas with high seismic hazard but infrequent large earthquakes where data is sparse and difficult to interpret. These data limitations lead to uncertainties in the development of site-specific ground motions, which are crucial for engineering risk assessments. To address these challenges, physics-based regional-scale ground-motion simulations have been developed. With the emergence of exaflop-scale computing ecosystems, it is now possible to simulate regional earthquake processes at unprecedented fidelity and generate the large number of fault rupture realizations necessary to characterize both intra- and inter-event ground-motion variability. This article introduces a new database of simulated earthquake ground motions, created for applications in earthquake engineering, earthquake planning, and emergency response. The inaugural version of the database features simulated ground motions for a magnitude 7 Hayward Fault earthquake in the San Francisco Bay Region (SFBR), using the EarthQuake SIMulation (EQSIM) simulation framework and the Graves–Pitarka kinematic rupture model. The aim is to provide high-fidelity, spatially dense, three-component motions generated on the Department of Energy’s (DOE) newest generation of graphics processing unit (GPU)-accelerated supercomputers. These motions are being made openly available to the engineering, scientific, and disaster planning communities. In addition, this work develops protocols for the efficient dissemination of these large data sets and emphasizes community engagement to build confidence in their application. This article discusses the methodology behind the data, underlying software verification and validation, scalable data management, and a user interface for data access. The goal is to facilitate widespread use and elicit expert feedback to maximize the utility and exploitation of simulated motions. While the initial focus is on the San Francisco Region, simulations for additional regions will be added as the DOE program progresses.

Simulated ground-motion database↗
DOE OSTI2021

Guidelines for Broadband Seismometer Testing

The Comprehensive Nuclear-Test-Ban Treaty mandates the four verification technologies to be used by the International Monitoring System (IMS) to monitor compliance to the Treaty. These technologies are seismic, hydroacoustic, infrasound, and radionuclide. Operational manuals for each of these technologies specify the requirements that equipment installed at IMS stations must meet for data from said stations to be accepted by the International Data Center. Following a model in which a revised set of infrasound sensor specifications and accompanying test procedures were developed in an international collaboration overseen by the Comprehensive Nuclear-Test-Ban Treaty Organization, the same process was recently followed for seismometers potentially of use by the IMS. This document, a product of that collaboration, defines key concepts and recommends test procedures to be followed to characterize broadband seismometers in a consistent and standardized way. Results from these test procedures can be used to verify that a seismometer meets the manufacturers stated specifications and IMS requirements. ACKNOWLEDGEMENTS The writers of this document would like to acknowledge Charles R. Hutt, John R. Evans, Fred Followill, Robert L. Nigbor, and Erhard Wielandt, the authors of the 2009 Guidelines for Standardized Testing of Broadband Seismometers and Accelerometers funded by the U.S. Department of Interior, on which this document is based. Although this document is limited in scope to seismometers and much of the document has been updated, the structure remains essentially the same.

58 GEOSCIENCES↗
DOE OSTI2020

EGS Collab Experiment 1: SIMFIP Notch-164 GRL Paper

Characterizing the stimulation mode of a fracture is critical to assess the hydraulic efficiency and the seismic risk related to deep fluid manipulations. We have monitored the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ~1.5 km depth in the crystalline rock of the Sanford Underground Research Facility (USA). The fracture initiates at 61% of the minimum horizontal stress by micro-shearing of the borehole on a foliation plane. As the fluid pressure increases further, borehole axial and radial displacements increase with injection time highlighting the opening and sliding of a new hydrofracture growing ~10 m away from the borehole, in accordance with the ambient normal stress regime and in alignment with the microseismicity. Our study reveals how fluid-driven fracture stimulation can be facilitated by a mixed-mode process controlled by the complex hydromechanical evolution of the growing fracture. The data presented in this submission refer to the SIMFIP measurements and analyses of the stimulation tests conducted on the 164 ft (50 m) notch of the Sanford Underground Research Facility (SURF), during the EGS-Collab test 1. In addition to the datafiles, there is the draft of a manuscript submitted to Geophysical Research Letters (GRL).

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

Monitoring CO 2 Storage at Cranfield, Mississippi with Time-Lapse Offset VSP – Using Integration and Modeling to Reduce Uncertainty

A time-lapse Offset Vertical Seismic Profile (OVSP) data set was acquired as part of a subsurface monitoring program for geologic sequestration of CO 2 . The storage site at Cranfield, near Natchez, Mississippi, is part of a detailed area study (DAS) site for geologic carbon sequestration operated by the U.S. Dept. of Energy’s Southeast Regional Carbon Sequestration Partnership (SECARB). The DAS site includes three boreholes, an injection well and two monitoring wells. The project team selected the DAS site to examine CO 2 sequestration multiphase fluid flow and pressure at the interwell scale in a brine reservoir. The time-lapse (TL) OVSP was part of an integrated monitoring program that included well logs, crosswell seismic, electrical resistance tomography and 4D surface seismic. The goals of the OVSP were to detect the CO 2 induced change in seismic response, give information about the spatial distribution of CO 2 near the injection well and to help tie the high-resolution borehole monitoring to the 4D surface data. The VSP data were acquired in well CFU 31-F1, which is the ~3200 m deep CO 2 injection well at the DAS site. A preinjection survey was recorded in late 2009 with injection beginning in December 2009, and a post injection survey was conducted in Nov 2010 following injection of about 250 kT of CO 2 . The sensor array for both surveys was a 50-level, 3-component, Sercel MaxiWave system with 15 m (49 ft) spacing between levels. The source for both surveys was an accelerated weight drop, with different source trucks used for the two surveys. Consistent time-lapse processing was applied to both data sets. Time-lapse processing generated difference corridor stacks to investigate CO 2 induced reflection amplitude changes from each source point. Corridor stacks were used for amplitude analysis to maximize the signal-to-noise ratio (S/N) for each shot point. Spatial variation in reflectivity (used to ‘map’ the plume) was similar in magnitude to the corridor stacks but, due to relatively lower S/N, the results were less consistent and more sensitive to processing and therefore are not presented. We examined the overall time-lapse repeatability of the OVSP data using three methods, the NRMS and Predictability (Pred) measures of Kragh and Christie (2002) and the signal-to-distortion ratio (SDR) method of Cantillo (2011). Because time-lapse noise was comparable to the observed change, multiple methods were used to analyze data reliability. The reflections from the top and base reservoir were identified on the corridor stacks by correlation with a synthetic response generated from the well logs. A consistent change in the corridor stack amplitudes from pre- to post-CO 2 injection was found for both the top and base reservoir reflections on all ten shot locations analyzed. In addition to the well-log synthetic response, a finite-difference elastic wave propagation model was built based on rock/fluid properties obtained from well logs, with CO 2 induced changes guided by time-lapse crosswell seismic tomography (Ajo-Franklin, et al., 2013) acquired at the DAS site. Time-lapse seismic tomography indicated that two reservoir zones were affected by the flood. The modeling established that interpretation of the VSP trough and peak event amplitudes as reflectivity from the top and bottom of reservoir is appropriate even with possible tuning effects. Importantly, this top/base change gives confidence in an interpretation that these changes arise from within the reservoir, not from bounding lithology. The modeled time-lapse change and the observed field data change from 10 shotpoints are in agreement for both magnitude and polarity of amplitude change for top and base of reservoir. Therefore, we conclude the stored CO 2 has been successfully detected and, furthermore, the observed seismic reflection change can be applied to Cranfield’s 4D surface seismic for spatially delineating the CO 2 /brine interface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

LLNL-G3Dv3: Global P wave tomography model for improved regional and teleseismic travel time prediction

[1] We develop a global-scale P wave velocity model (LLNL-G3Dv3) designed to accurately predict seismic travel times at regional and teleseismic distances simultaneously. The model provides a new image of Earth's interior, but the underlying practical purpose of the model is to provide enhanced seismic event location capabilities. The LLNL-G3Dv3 model is based on ~2.8 millionP and Pnarrivals that are re-processed using our global multiple-event locator called Bayesloc. We construct LLNL-G3Dv3 within a spherical tessellation based framework, allowing for explicit representation of undulating and discontinuous layers including the crust and transition zone layers. Using a multiscale inversion technique, regional trends as well as fine details are captured where the data allow. LLNL-G3Dv3 exhibits large-scale structures including cratons and superplumes as well numerous complex details in the upper mantle including within the transition zone. Particularly, the model reveals new details of a vast network of subducted slabs trapped within the transition beneath much of Eurasia, including beneath the Tibetan Plateau. We demonstrate the impact of Bayesloc multiple-event location on the resulting tomographic images through comparison with images produced without the benefit of multiple-event constraints (single-event locations). We find that the multiple-event locations allow for better reconciliation of the large set of direct P phases recorded at 0–97° distance and yield a smoother and more continuous image relative to the single-event locations. Travel times predicted from a 3-D model are also found to be strongly influenced by the initial locations of the input data, even when an iterative inversion/relocation technique is employed.

58 GEOSCIENCES↗
DOE OSTI2023

A case study of denoising microseismic data in the Farnsworth, TX CO 2 EOR site

As geologic carbon sequestration projects begin to be funded with a higher degree of frequency, microseismic monitoring will become more necessary to establish caprock integrity and induced seismicity risk. The Farnsworth, TX site, which hosts enhanced oil recovery operations, provides an ideal laboratory to test essential components of microseismic monitoring such as surface station placement, denoising, and different autopicker methods. We find that despite careful optimization, the surface station placement at the Farnsworth site was insufficient given the high level of background noise from the industrial operations at the oil field. Thus, we use only borehole geophone data in our processing. Two denoising techniques were examined: the DeepDenoiser and the continuous wavelet transform. The continuous wavelet transform was shown to be a valuable tool in converting between raw waveforms to processed denoised waveforms for microseismic monitoring. In the case of the noisy waveforms, only ten detections are found, compared with 90 in the denoised data for a two-hour window. The DeepDenoiser suffered from the fact that the training data was regional earthquake data. In addition, the PhaseNet machine-learning autopicker was applied to both the noisy and the denoised data, and this algorithm detected thousands of more arrivals in the data denoised with the continuous wavelet transform technique.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

Abstract Subsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity ( V S ) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface V S variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred V S variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The V S variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

58 GEOSCIENCES↗
DOE OSTI2026

Machine learning for seismic low-frequency extrapolation

The cycle-skipping problem that plagues full waveform inversion (FWI) can be at least partially mitigated if low frequencies (which encode the kinematics of wave propagation in seismic data) are recorded. However, seismic sources and receivers are band-limited, so seismic data does not generally include signals down to 0 Hz. To improve our ability to solve the seismic inverse problem, one can synthesize this missing low-frequency (LF) content from the recorded high-frequency (HF) data using machine learning (ML) models. Deep learning models such as convolutional neural networks (CNNs) demonstrate impressive ability to perform low frequency extrapolation. However, such models require powerful hardware (GPU machines) and careful training. We assess the extrapolation capabilities of three different ML models that do not require GPU machines, namely, random forest, Gaussian process regression and gradient boosting, on both synthetic and real data. Experimental results on two synthetic data sets (generated from a low velocity lens embedded in a homogeneous medium, and the Marmousi model) demonstrate that FWI applied to the extrapolated data consistently improves inversion accuracy relative to FWI applied to the original data sets that do not contain low frequencies. Application of low-frequency extrapolation to real data from the Northwest Shelf of Australia demonstrates that tree-based ML models such as gradient boosting can outperform CNNs in terms of both accuracy and computational cost on non-GPU architectures.

58 GEOSCIENCES↗