Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SEISMICITY”

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 685 records · Page 38

Identifying High Potential Well Targets with 3D Seismic and Mineralogy

Seismic reflection the primary tool used in petroleum exploration and production, but use in geothermal exploration is less standard, in part due to cost but also due to the challenges in identifying the highly-permeable zones essential for economic hydrothermal systems [e.g. Louie et al., 2011; Majer, 2003]. Newer technology, such as wireless sensors and low-cost high performance computing, has helped reduce the cost and effort needed to conduct 3D surveys. The second difficulty, identifying permeable zones, has been less tractable so far. Here we report on the use of seismic attributes from a 3D seismic survey to identify and map permeable zones in a hydrothermal area.

15 GEOTHERMAL ENERGY↗

The use of multiwavelets for uncertainty estimation in seismic surface wave dispersion.

This report describes a new single-station analysis method to estimate the dispersion and uncer- tainty of seismic surface waves using the multiwavelet transform. Typically, when estimating the dispersion of a surface wave using only a single seismic station, the seismogram is decomposed into a series of narrow-band realizations using a bank of narrow-band filters. By then enveloping and normalizing the filtered seismograms and identifying the maximum power as a function of frequency, the group velocity can be estimated if the source-receiver distance is known. However, using the filter bank method, there is no robust way to estimate uncertainty. In this report, I in- troduce a new method of estimating the group velocity that includes an estimate of uncertainty. The method is similar to the conventional filter bank method, but uses a class of functions, called Slepian wavelets, to compute a series of wavelet transforms of the data. Each wavelet transform is mathematically similar to a filter bank, however, the time-frequency tradeoff is optimized. By taking multiple wavelet transforms, I form a population of dispersion estimates from which stan- dard statistical methods can be used to estimate uncertainty. I demonstrate the utility of this new method by applying it to synthetic data as well as ambient-noise surface-wave cross-correlelograms recorded by the University of Nevada Seismic Network.

58 GEOSCIENCES↗

Seismic Analysis of Spatio-Temporal Fracture Generation During EGS Resource Development

This DOE GTO funded project was a collaborative effort between Array Information Technology and the University of California at Berkeley. The project goal is to improve technology to assess in-situ stress magnitude and orientation, kinematic fracture parameters, rupture size, rupture orientation as well as temporal changes and volumetric distribution of the injected fluid during Enhanced Geothermal System (EGS) resource development. The project leveraged high-frequency seismic data recorded by the LBNL 34-station permanent geophone network and seismic broadband data recorded by a temporary 33-station seismometer network that operated in The Geysers during the injection phase of the DOE GTO funded EGS demonstration project at the well Prati-32. The operation of the broadband network with high station density and good azimuthal coverage to monitor the injection phase produced an unprecedented dataset that is typically unavailable for EGS operations. The two datasets offer the advantages of studying the parameters under investigation over a broader frequency band and the higher sensitivity of the broadband seismometers resulted in a richer dataset with a lower magnitude of completeness. Moment tensor (MT) inversion have been employed to develop a catalog of earthquake source solutions including moment tensor and finite source rupture area. This resulted in a 168-event waveform-based seismic moment tensor catalog for events 0.7≤Mw≤3.9 based on a semi-automated moment tensor analysis technique. The MT solutions were used for in-situ stress estimation during the EGS injection phases. We found an approximate 15-degree counterclockwise rotation of the least compressive stress sigma-3, and a rotation of the maximum compressive stress sigma-1 toward the vertical as the injected volume of water increased. We developed a rupture-area magnitude scaling relationship of The Geysers earthquakes obtained from finite-source inversion for fault slip that reveals a very high correlation to published results. Based on the rupture area-magnitude relationship we mapped Mw from the LBNL earthquake catalog to estimate the activated fracture area, its location and orientation within the injection volume of Prati-32. The double-difference Wadati (DDW) technique was applied to high-accuracy P- and S-wave differential travel times derived from waveform cross correlation to image fluid saturation in the subsurface based on Vp/Vs-ratio. The recently developed technique has been applied to estimate temporal changes of fluid distribution in fault zones and in volcanic regions but not yet in engineered geo-reservoirs. The improvement of the technology includes high-resolution waveform cross correlation applied to the high-quality datasets with lower magnitude of completeness resulting in a higher number of recorded events. The improvement allowed to estimate both, the spatial and temporal distribution of injected fluids and the state of the injected fluid during reservoir stimulation. Results of the DDW analysis revealed the development of an injection-derived steam plume during the early injection phase, particularly in the northern part of the Prati-32 study area.

15 GEOTHERMAL ENERGY↗

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↗

Near-Field Imaging of Shallow Chemical Explosions in Granite Using Change Detection Methods with Surface and Borehole Seismic Data

Explosions detonated in geologic media damage it in various ways via processes that include vaporization, fracturing, crushing of interstitial pores, etc. Seismic waves interact with the altered media in ways that could be important to the discrimination, characterization, and location of the explosions. As part of the Source Physics Experiment, we acquired multiple pre- and post-explosion near-field seismic datasets and analyzed changes to seismic P-wave velocity. Our results indicate that the first explosion detonated in an intact media can cause fracturing and, consequently, a decrease in P-wave velocity. After the first explosion, subsequent detonations in the pre-damaged media have limited discernible effects. We hypothesize this is due to the stress-relief provided by a now pre-existing network of fractures into which gasses produced by the explosion migrate. We also see an overall increase in velocity of the damaged region over time, either due to a slow healing process or closing of the fractures by subsequent explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Regolith, rock and fluid distributions at the Upper Colorado River Basin via a multicomponent seismic imaging approach

In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.

58 GEOSCIENCES↗

Evaluation of the PhaseNet Model Applied to the IMS Seismic Network

Producing a complete and accurate set of signal detections is essential for automatically building and characterizing seismic events of interest for nuclear explosion monitoring. Signal detection algorithms have been an area of research for decades, but still produce large quantities of false detections and misidentify real signals that must be detected to produce a complete global catalog of events of interest. Deep learning methods have shown promising capabilities in effectively characterizing seismic signals for complex tasks such as identifying phase arrival times. We use the PhaseNet model, a UNet-based Neural Network, trained on local distance data from northern California to predict seismic arrivals on data from the International Monitoring System (IMS) global network. We use an analyst-curated bulletin generated from this data set to compare the performance of PhaseNet to that of the Short-Term Average/Long-Term Average (STA/LTA) algorithm. We find that PhaseNet has the potential of outperforming traditional processing methods and recommend the training of a new model with the IMS data to achieve optimal performance.

58 GEOSCIENCES↗

Assessing the Accuracy of a Borehole-Controlled Velocity Model of Yucca Flat, Nevada, Using Large-N Seismic Data

The Source Physics Experiment (SPE) aims to improve nuclear test monitoring capabilities through a series of chemical explosions in Yucca Flat within the Nevada National Security Site (NNSS) (Snelson et al., 2013). With geologic data from over 950 boreholes, Yucca Flat has excellent borehole control on stratigraphy. These data were used to create a Geologic Framework Model (GFM) of the basin to aid in the study of seismic wave propagation through the Earth. However, the seismic velocities of the geologic units in the GFM are less well-defined. With the acquisition of five new active source Large-N datasets, we can now quantitatively assess the accuracy of the seismic model previously controlled by borehole data alone.

58 GEOSCIENCES↗

Time-Lapse Vertical Seismic Profiling (VSP) for CO 2 Storage in a Depleted Oil Field in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. This report describes the Dover 33 VSP study to test the effectiveness of time-lapse Vertical Seismic Profile (VSP) for detecting and delineating a plume of more than 271,000 tonnes of CO2 injected into the Brown Niagaran and A-1 Carbonate formations within the Dover 33 reef between March 2013 and September 2016. Five 2D walkaway VSP (WVSP) source lines were acquired by SIGMA3 in September 2016 to investigate the possible time-lapse response in both P-wave and PS-wave seismic data. The data was compared to the same survey geometry acquired in March 2013 by SR2020. The seismic energy was recorded into an 80-level, three-component geophone array deployed into the Dover 1-33 well and placed just above the reef. P-wave and PS-wave reflection images were produced for each of the three source lines from both the 2013 baseline and 2016 monitor surveys. The images were then compared to look for changes in the reflectivity at and around the injection location that might indicate how the CO 2 has moved over this time period. This technique revealed several localized areas with sizable impedance differences inside the reef where CO 2 would be expected; however, a large number of similar impedance “hotspots” were also detected outside the reef in areas where injected CO 2 would not be expected. Therefore, a second analysis that involved calculating P-wave and S-wave travel time differences between the 2013 and 2016 VSPs was conducted to look for a change that could be caused by the CO 2 plume.

01 COAL, LIGNITE, AND PEAT↗

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗

Auto-Curation of Seismic Event Data for Signal Denoising

Denoising contaminated seismic signals for later processing is a fundamental problem in seismic signals analysis. Neural network approaches have shown success denoising local signals when trained on short-time Fourier transform spectrograms. One challenge of this approach is the onerous process of hand-labeling event signals for training. By leveraging the SCALODEEP seismic event detector, we develop an automated set of techniques for labeling event data. Despite region specific challenges, training the neural network denoiser on machine curated events shows comparable performance to the neural network trained on hand curated events. We showcase our technique with two experiments, one using Utah regional data and one using regional data from the Korean peninsula.

58 GEOSCIENCES↗

Bridge Seismic Screening Tool (BSST), Version 2.0

The Regional Resiliency Assessment Program (RRAP) is a cooperative assessment of specific critical infrastructure within a designated geographic area and a regional analysis of the surrounding infrastructure that addresses a range of infrastructure resilience issues that could have regionally and nationally significant consequences. In 2018, DHS’s Cybersecurity and Infrastructure Security Agency (CISA) sponsored the Oregon Transportation Systems RRAP project in coordination with the Office of the Governor (under the oversight of the state resilience officer), the Oregon Office of Emergency Management (OEM), the Oregon Department of Transportation (ODOT), and other regional stakeholders (CISA 2021). This project focuses on assessing the impacts of a Cascadia Subduction Zone (CSZ) earthquake on state transportation systems and, in particular, how those impacts may affect the ability of emergency response efforts to move supplies into the region. The intended outcome of this analysis is the prioritization of transportation routes and modes for additional planning, investment, hardening, or other activities to enhance their resilience—and therefore, to enhance their ability to support response and recovery efforts following a CSZ earthquake. An important part of this transportation system-level assessment has been to assess the seismic vulnerability of the state highway system. In doing so, the RRAP project team used the Bridge Seismic Screening Tool (BSST) to assess, at a system-level, the potential impacts that a CSZ earthquake could have on state highway bridges (Bergerson et al. 2019).1 Argonne National Laboratory (Argonne), in collaboration with the Washington State Department of Transportation (WSDOT), originally developed the BSST as part of the 2017 Washington State Transportation Systems RRAP project, a sister project to the 2018 Oregon Transportation Systems RRAP project. Argonne updated the BSST during this more recent project in Oregon based on feedback from stakeholders and subject matter experts (SMEs) on the original version of the tool. The first step in the BSST is to assess the seismic vulnerability of roadway bridges following a CSZ earthquake to determine a projected or potential damage state. Damage states then help determine approximate reopening times for bridge crossings.2 This document provides details on the BSST methodology, the implementation of that tool to analyze the projected damage incurred in a CSZ earthquake scenario, and the determination of corresponding reopening times of interstate, state highway, and local bridges following such an event.

58 GEOSCIENCES↗

LLNL 5-2428: Fracture Permeability and Seismic Slip Behavior

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that also includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modelling approaches to reduce parameteruncertainty and better predict and mitigate seismic hazard at Utah FORGE and future EGS sites.We propose a novel approach that incorporates 3D physics-based Earthquake simulations in THMC models, herein referred to as “THMC+E” models. This capability will enable improvedengineering decisions at Utah FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. Utah FORGE management and future EGS operators may employ results of THMC+E models for decision making purposes throughout the lifetime of the field operation. For example, before production well installation, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite ofpossible well locations and flow rates, thus enabling evaluation of optimal production well placement. Such efforts will be conducted throughout the lifetime of the project, whereby additional laboratory experiments will constrain key model parameters and machine learning (ML) will reduce the size of the parameter space and the associated uncertainty. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗

Identification of Faults Susceptible to Induced Seismicity Final Report

The primary objective of this project (referred to as the State of Stress project or SoS project) was to predict the presence of faults susceptible to movement in the presence of fluid injection as a consequence of vertical pressure migration from the storage unit to the crystalline basement. These faults included those difficult to resolve with traditional interpretations of surface seismic: faults with no or minimal offset, faults in the crystalline basement where no seismic reflectors are present, subseismic faults (due to size), and faults undetected (e.g., due to unfavorable orientation with surface seismic survey).

25 ENERGY STORAGE↗

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modeling approaches to reduce parameter uncertainty and better predict/mitigate seismic hazard at EGS sites. Our novel approach couples 3D physics-based earthquake simulations with THMC models (THMc+E). This capability will enable improve engineering decisions at Utah-FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. For example, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite of flow rates and under uncertainty, thus enabling evaluation of optimal circulation strategy. Laboratory experiments will be performed to constrain key model parameters and Bayesian techniques will provide a probabilistic evaluation of parameters used in models. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗

Toward Global Regional Seismic Moment Tensor Inversion with Three-Dimensional Earth Models for Nuclear Explosion Monitoring with Sparse Networks: Demonstration of Reciprocity for Strain Greens Tensor Database Simulation with Salvus

Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).

58 GEOSCIENCES↗

5-2428: Fracture Permeability Impact on Seismic Slip Behavior

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modeling approaches to reduce parameter uncertainty and better predict/mitigate seismic hazard at EGS sites. Our novel approach couples 3D physics-based earthquake simulations with THMC models (THMc+E). This capability will enable improve engineering decisions at Utah-FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. For example, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite of flow rates and under uncertainty, thus enabling evaluation of optimal circulation strategy. Laboratory experiments will be performed to constrain key model parameters and Bayesian techniques will provide a probabilistic evaluation of parameters used in models. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗

American-Made Geothermal Geophone Prize: Ultra-High Temperature Seismic Tool for Geothermal Wells Summary Report

This work was conducted in support of the American Made Geothermal Prize. The following data summary report presents the testing conducted at Sandia National Labs to validate the performance of the Ultra-High Temperature Seismic Tool for Geothermal Wells. The goal of the testing was to measure the sensitivity of the device to seismic vibrations and reliability of the instrument at elevated temperatures. To this end, two tests were conducted: 1) Ambient Temperature Seismic Testing, which measured the response of the tool to a sweep of frequencies from 1 to 1000 Hz, and 2) Elevated Temperature Survivability Testing which measured the voltage response of the device at 225°C over a month-long testing window. The details of the testing methodology and summary of the tests are presented herein.

15 GEOTHERMAL ENERGY↗