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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Deep compressed seismic learning for fast location and moment tensor inferences with natural and induced seismicity

Fast detection and characterization of seismic sources is crucial for decision-making and warning systems that monitor natural and induced seismicity. However, besides the laying out of ever denser monitoring networks of seismic instruments, the incorporation of new sensor technologies such as Distributed Acoustic Sensing (DAS) further challenges our processing capabilities to deliver short turnaround answers from seismic monitoring. In response, this work describes a methodology for the learning of the seismological parameters: location and moment tensor from compressed seismic records. In this method, data dimensionality is reduced by applying a general encoding protocol derived from the principles of compressive sensing. The data in compressed form is then fed directly to a convolutional neural network that outputs fast predictions of the seismic source parameters. Thus, the proposed methodology can not only expedite data transmission from the field to the processing center, but also remove the decompression overhead that would be required for the application of traditional processing methods. An autoencoder is also explored as an equivalent alternative to perform the same job. We observe that the CS-based compression requires only a fraction of the computing power, time, data and expertise required to design and train an autoencoder to perform the same task. Implementation of the CS-method with a continuous flow of data together with generalization of the principles to other applications such as classification are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Seismic Spatial Gradients and Machine Learning-Based Classifiers for Explosion Monitoring (LDRD 218327)

This final report summarizes the work completed under the Laboratory Directed Research and Development (LDRD) project “Seismic Spatial Gradients as a Machine Learning-Based Classifier for Explosion Monitoring.” The overarching goal of the project was to explore the efficacy of using machine learning-based classification algorithms where the input data are the spatial gradient of the seismic wavefield collected at a single point on the Earth’s surface. The methods that I describe here are in direct contrast to conventional methods of seismic discrimination which typically rely on a spatially extended network of instruments and physics-based wavefield attributes such as, for example, the ratio between $\textit{P}$ and $\textit{S}$ waves. Rather, we use the spatial gradient of the seismic wavefield observed at a single point on the Earth’s surface and data processing approaches inspired by the machine learning community. We tested two algorithms, a neural network and a modified version of principal component analysis termed Spectrally Filtered Principal Component Analysis (SFPCA). To test these algorithms, we first conducted a series of numerical tests using synthetic data and then conducted a small-scale controlled field experiment. The tests using synthetic data showed that both algorithms had high success rates on gradiometric data, even when simulated noise was added to the signal. Furthermore, we found that using seismic spatial gradients increased the performance of our discrimination algorithms when compared to using just the traditional translational motion seismic data. The tests with field data also showed a high degree of discriminative success.

58 GEOSCIENCES↗
DOE OSTI2022

Review: Induced Seismicity During Geoenergy Development—A Hydromechanical Perspective

The basic triggering mechanism underlying induced seismicity traces back to the mid-1960s that relied on the process of pore-fluid pressure diffusion. The last decade has experienced a renaissance of induced seismicity research and data proliferation. An unprecedent opportunity is presented to us to synthesize the robust growth in knowledge. The objective of this article is to provide a concise review of the triggering mechanisms of induced earthquakes with a focus on hydro-mechanical processes. Four mechanisms are reviewed: pore-fluid pressure diffusion, poroelastic stress, Coulomb static stress transfer, and aseismic slip. For each, an introduction of the concept is presented, followed by case studies. Diving into these mechanisms sheds light on several outstanding questions. For example, why did some earthquakes occur far from fluid injection or after injection stopped? Our review converges on the following conclusions: (a) Pore-fluid pressure diffusion remains a basic mechanism for initiating inducing seismicity in the near-field. (b) Poroelastic stresses and aseismic slip play an important role in inducing seismicity in regions beyond the influence of pore-fluid pressure diffusion. (c) Coulomb static stress transfer from earlier seismicity is shown to be a viable mechanism for increasing stresses on mainshock faults. (d) Multiple mechanisms have operated concurrently or consecutively at most induced seismicity sites. (e) Carbon dioxide injection is succeeding without inducing earthquakes and much can be learned from its success. Future research opportunities exist in deepening the understanding of physical and chemical processes in the nexus of geoenergy development and fluid motion in the Earth’s crust.

58 GEOSCIENCES↗
DOE OSTI2021

Brady Geodatabase for Geothermal Exploration Artificial Intelligence

These files contain the geodatabases related to Brady's Geothermal Field. It includes all input and output files for the Geothermal Exploration Artificial Intelligence. Input and output files are sorted into three categories: raw data, pre-processed data, and analysis (post-processed data). In each of these categories there are six additional types of raster catalogs which are titled Radar, SWIR, Thermal, Geophysics, Geology, and Wells. These inputs and outputs were used with the Geothermal Exploration Artificial Intelligence to identify indicators of blind geothermal systems at the Brady Hot Springs Geothermal Site. The included zip file is a geodatabase to be used with ArcGIS and the tar file is an inclusive database that encompasses the inputs and outputs for the Brady Hot Springs Geothermal Site.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Upper crustal structure beneath East Java from ambient noise tomography: A preliminary result

East Java has a fairly complex geological structure. Physiographically East Java can be divided into three zones, i.e. the Southern Mountains zone in the southern part, the Kendeng zone in the middle part, and the Rembang zone in the northern part. Most of the seismic hazards in this region are due to processes in the upper crust. In this study, the Ambient Noise Tomography (ANT) method is used to image the upper crustal structure beneath East Java. We have used seismic waveform data recorded by 8Meteorological, Climatological and Geophysical Agency (BMKG) stationary seismographic stations and 16 portable seismographs installed for 2 to 8 weeks. The data were processed to obtain waveforms fromnoise cross-correlation between pairs of seismographic stations. Our preliminary results indicate that the Kendeng zone, an area of low gravity anomaly, is associated with a low velocity zone. On the other hand, the southern mountain range, which has a high gravity anomaly, is related to a high velocity anomaly as shown by our tomographic images.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

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 OSTI2023

Exploring for Superhot Geothermal Targets in Magmatic Settings: 2022 Field Campaign at Newberry Volcano

This paper presents preliminary results from a subset of work carried out as part of a multinational research project entitled DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN), supported by the U.S. Department of Energy (DOE) and Geothermica, a joint effort by EU member states and associated countries. The DEEPEN project will develop a customized approach to exploration for supercritical and superhot geothermal plays in magmatic systems, which will be applied to two demonstration sites. This paper summarizes field activities carried out at the U.S. demonstration site, Newberry Volcano in central Oregon. The objective of this work effort is to refine the subsurface model of Newberry Volcano, with special focus on deeper zones including the magmatic plumbing system and other key geologic elements. New data collection included gravity and wideband magnetotelluric (MT) surveys, as well as reinstallation of a seismic network. The National Renewable Energy Laboratory (NREL) and Enthalpion Energy LLC (Enthalpion) worked with the Deschutes National Forest Fort Rock District to use a low ground disturbance method of MT deployment to collect MT data inside the caldera and other restricted areas inside the National Volcanic Monument. This opened these areas to geophysical exploration for the first time in decades. Sites along and adjacent to the south rim of the volcano constituted the primary survey objectives. A team from Lawrence Berkeley National Laboratory (LBNL), the U.S. Geological Survey (USGS), and AltaRock also began the process of reinstalling the seismic network from the AltaRock enhanced geothermal system (EGS) demonstration in anticipation of further development activities at the site. The data ingestion, reduction, and analysis phase of the project is ongoing. We are currently processing the MT and gravity data and are developing a new, highly GPU-accelerated, 3D joint MT and gravity inversion to better localize the south rim/south flank conductive target and better understand its relationship to deep heat, fluid sources, and surface extrusive features. Joint inversions, which have not yet been undertaken at Newberry, will allow us to obtain constraints on the geologic model that cannot be determined from each method in isolation, improving our ability to image key geologic features at depth.

geophysics↗
DOE OSTI2017

Final Scientific/Technical Report of Gas Hydrate Dynamics on the Alaskan Beaufort Continental Slope: Modeling and Field Characterization

Methane Hydrates, a solid form of methane and water, exist at high pressures and low temperatures, occurs on every continental margin on Earth, represents one of the largest reservoirs of carbon on the planet, and, if destabilized, may play an important role in both slope stability and climate change. For decades, researchers have studied methane hydrates with the hope of determining if methane hydrates are destabilizing, and if so, how this destabilization might impact slope stability and ocean/atmosphere carbon budgets. In the past ~5 years, it has become well established that the upper “feather-edge” of methane hydrate stability (intermediate water depths of ~200-500 meters below sea level) represents an important frontier for methane hydrates stability research, as this zone is most susceptible to destabilization due to minor fluctuations in ocean temperature in space and time. The Arctic Ocean—one of the fastest warming regions on Earth—is perhaps the best place to study possible changes to methane hydrate stability due to ocean warming. To address the stability of methane hydrates at intermediate ocean depths, Southern Methodist University in partnership with Oregon State University and The United State Geological Survey at Woods Hole began investigating methane hydrate stability in intermediate water depths below both the US Beaufort Sea and the Atlantic Margin, from 2012-2017. The work was funded by the Department of Energy’s (DOE) National Energy Technology Laboratory (NETL). The key goal of the SMU component of this study was to collect the first ever heat flow data in the Beaufort Sea and compare measured shallow (probe-based1) heat flow values with deeper (BSR-derived2) heat flow values, and from this, determine whether hydrates were in thermal equilibrium. In September 2016, SMU/OSU collected the first ever heat flow measurements in the US Beaufort Sea. Despite poor weather and rough seas, the cruise was a success, with 116 heat flow measurements acquired across the margin, spanning 4 transects separated by more than 400 km. Useable heat flow data exists for 97% (113) of probe heat flow measurements, revealing a clear picture of regional heat flow across the basin. During the past 8 months since the cruise, SMU researchers have processed the heat flow and thermal conductivity measurements and compared results to deeper heat flow estimates obtained from seismic data. The analysis reveals clear, consistent trends: All probe heat flow measurements in depths greater than 800 mbsl are consistent with BSR-derived values; heat flow measurements obtained in water depths between ~250-750 mbsl are systematically lower than those estimated from BSRs; and heat flow estimates in water depths shallower than ~250 mbsl are systematically warmer than deeper estimates. The consistency between shallow (probe) and deep (BSR) heat flow measurements at depths greater than ~750 m where ocean temperature changes are minimal supports the premise that the hydrates consist primarily of methane and represent a valuable tool for estimating heat flow. The anomalous cooling trend observed in the upper 250 m is consistent with expected seasonal effects observed in shallow ocean buoy measurements in the arctic, when cold, less dense melting sea ice cools the upper 200 m of the ocean during the summer as ice melting occurs. The discrepancy in heat flow at intermediate water depths is best explained via recent intermediate ocean temperature warming, where long-term (annual or longer) warming intermediate ocean bottom waters result in an anomalously low heat flow in shallow heat flow measurements. Using the characteristic 1D time-length scale for diffusion, we estimate that ocean temperature warming began no later than ~1200 years ago but arguably much more recently as results are limited by seismic resolution. More importantly, our analysis indicates methane hydrate is destabilizing not only in the upper feather edge (200-500 mbsl) but at depths as great as 750 mbsl. The intermediate ocean warming rate supports previous studies suggesting geologically rapid warming (>0.1 deg C/decade) at intermediate ocean depths in the Beaufort Sea. Assuming no further changes or additional warming, our analysis indicates methane hydrates will destabilize at seafloor depths shallower than 750 mbsl in the Beaufort Sea within the next ~3000 years. 1 Probe outfitted with sensors inserted into the seafloor sediment 2 Bottom-simulating reflector (BSR) seismic data indicates presence of hydrate deposits

58 GEOSCIENCES↗
DOE OSTI2019

Integrating distributed acoustic sensing, borehole 3C geophone array, and surface seismic array data to identify long-period long-duration seismic events during stimulation of a Marcellus Shale gas reservoir

Microseismic monitoring by downhole geophones, surface seismic, fiber-optic distributed acoustic sensing (DAS), and distributed temperature sensing (DTS) observations were made during the hydraulic fracture stimulation of the MIP-3H well in the Marcellus Shale in northern West Virginia. DAS and DTS data measure the fiber strain and temperature, respectively, along a fiber-optic cable cemented behind the casing of the well. The presence of long-period long-duration (LPLD) events is evaluated in the borehole geophones, DAS data, and surface seismic data of one of the MIP-3H stimulated stages. LPLD events are generally overlooked during the conventional processing of microseismic data, but they represent significant nonbrittle deformation produced during hydraulic fracture stimulation. In a single stage that was examined, 160 preexisting fractures and two faults of suboptimal orientation are noted in the image logs. We identified two low-frequency ([Formula: see text]) events of large temporal duration (tens of seconds) by comparing the surface seismic data, borehole geophone data, and DAS amplitude spectra of one of the MIP-3H stages. Spectrograms of DAS traces in time and depth reveal that the first low-frequency event might be an injection noise that has footprints on all DAS channels above the stimulated stage. However, the surface seismic array indicates an LPLD event concurrent with the first low-frequency event on DAS. The second LPLD event on DAS data and surface seismic data is related to a local deformation and does not have footprints on all DAS channels. The interpreted events have duration less than 100 s with frequencies concentrated below 10 Hz, and are accompanied by microseismic events.

Geochemistry & Geophysics↗
DOE OSTI2020

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Hot Springs and Geysers: Exploring Historical and Modern Impacts of Geothermal Energy Production on Associated Natural Surface Systems and Standardizing Management Practices

Surface thermal features, most notably hot springs and geysers are increasingly being recognized for their importance to ecosystems, indigenous cultures, and in some cases agriculture, recreation, and tourism. Geothermal project development poses a potential risk to these natural features but current regulatory requirements for assessing and managing these risks during exploration, permitting and monitoring are somewhat inconsistent and unpredictable across different geothermal fields. This has resulted in uncertainty and increases in exploration risk for geothermal energy developers that have led to costly project delays, cancellations, or hesitation to commit. Varying regulatory requirements may also influence public perception, fostering confusion, distrust and ultimately opposition to geothermal projects, further contributing to project delays or cancellations. At a time when there is an increasing urgency for reliable baseload clean energy, geothermal is a net-zero, renewable solution that additionally provides access to more equitable and environmentally just clean power. Continued integration of geothermal energy into the national energy roadmap can be facilitated through consistent and predictable permitting, providing regulators the framework they need, developers a clear path forward, and transparency that the public deserves. This project, currently in its beginning phases, seeks to address this important issue by providing a technical basis from which to build a preliminary protocol for assessing and managing potential impacts from new or existing geothermal energy projects to surface thermal features and their associated ecosystems. Development of this preliminary protocol will be informed by (1) a literature review of well-documented case studies in the western U.S. and New Zealand to understand the range of conditions that exemplify geothermal-surface thermal systems; (2) development of generic illustrative conceptual-numerical models to quantify, understand, and predict the first-order controls (e.g., pressure and permeability) on surface flows; and (3) additional independent and scientifically rigorous evaluations of geothermal-surface thermal system case studies from the Basin and Range Province that incorporate publicly available data as well as data provided by industry through data-sharing agreements. Learning from the successes of the process used to develop the Induced Seismicity Management Protocol (ISMP), we ultimately aim to use these initial efforts as a springboard for establishing a surface thermal feature management working group that will work collaboratively to finalize the protocol as well as co-create recommended best practices for implementation. We envision that the working group will primarily be composed of representatives from regulatory entities, government agencies, Tribes, academia, national laboratories, and industry, and will include early and regular engagement with community organizations and environmental groups. This will help ensure broad acceptance and implementation of the protocol, which will facilitate a more consistent, predictable, and standardized regulatory process, and help to ensure that geothermal energy continues to provide a reliable source of clean energy, and a pathway to achieving greater energy equity in the U.S.

Best Practices↗
DOE OSTI2019

Work Plan for the Site-wide SSHAC Level 3 Probabilistic Seismic Hazard Analysis at INL

This work plan provides the overall framework, key activities, key tasks, and schedule for the execution of a probabilistic seismic hazard analysis at Idaho National Laboratory (INL). The probabilistic seismic hazard analysis will follow processes and procedures appropriate for the Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 study provided in U. S. Nuclear Regulatory Commission guidance, NUREG-2213. The SSHAC Level 3 probabilistic seismic hazard analysis (PSHA) is being conducted as an INL site-wide effort. The work plan was developed and revised by participants who will perform the SSHAC Level 3 PSHA. The work plan also includes the quality assurance requirements commensurate with American Society of Mechanical Engineers Nuclear Quality Assurance NQA 1 2008/2009 and 2011 Addendum for design nuclear facilities under the U.S. Department of Energy and U. S. Nuclear Regulatory Commission. The scope of the SSHAC Level 3 PSHA is designed to produce the mean-centered definition of the seismic hazard with the appropriate treatment of uncertainties supported by the associated technical justifications. The SSHAC Level 3 study will properly implement and completely document the activities of evaluation, integration, and documentation. The evaluation process will consider data, models, and methods proposed by the larger technical community that are relevant to the PSHA at INL. The integration process will build seismic-source and ground-motion characterization models that represent the center, body, and range of technically defensible interpretations, taking into consideration the outcomes of the evaluation process. Key tasks and activities—such as a project kickoff meeting, three workshops, four formal working meetings, and a briefing meeting with the Participatory Peer Review Panel—will occur for the seven-member technical integration teams to construct seismic-source and ground-motion-characterization models. Both preliminary and final hazard-input documents will be developed to perform sensitivity analyses and the final PSHA, respectively. The entire study will be thoroughly documented in draft and final reports. The Participatory Peer Review Panel will review the SSHAC Level 3 study from start to finish, including the draft and final INL Site-wide PSHA reports, so they can provide a closure letter documenting the study’s adherence to SSHAC processes and technical adequacy of the final products.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

A comparison of earthquake backprojection imaging methods for dense local arrays

Backprojection imaging has recently become a practical method for local earthquake detection and location due to the deployment of densely sampled, continuously recorded, local seismograph arrays. While backprojection sometimes utilizes the full seismic waveform, the waveforms are often pre-processed and simplified to overcome imaging challenges. Real data issues include aliased station spacing, inadequate array aperture, inaccurate velocity model, low signal-to-noise ratio, large noise bursts and varying waveform polarity. We compare the performance of backprojection with four previously used data pre-processing methods: raw waveform, envelope, short-term averaging/long-term averaging and kurtosis. Our primary goal is to detect and locate events smaller than noise by stacking prior to detection to improve the signal-to-noise ratio. The objective is to identify an optimized strategy for automated imaging that is robust in the presence of real-data issues, has the lowest signal-to-noise thresholds for detection and for location, has the best spatial resolution of the source images, preserves magnitude, and considers computational cost. Imaging method performance is assessed using a real aftershock data set recorded by the dense AIDA array following the 2011 Virginia earthquake. Our comparisons show that raw-waveform backprojection provides the best spatial resolution, preserves magnitude and boosts signal to detect events smaller than noise, but is most sensitive to velocity error, polarity error and noise bursts. On the other hand, the other methods avoid polarity error and reduce sensitivity to velocity error, but sacrifice spatial resolution and cannot effectively reduce noise by stacking. Of these, only kurtosis is insensitive to large noise bursts while being as efficient as the raw-waveform method to lower the detection threshold; however, it does not preserve the magnitude information. For automatic detection and location of events in a large data set, we therefore recommend backprojecting kurtosis waveforms, followed by a second pass on the detected events using noise-filtered raw waveforms to achieve the best of all criteria.

58 GEOSCIENCES↗
DOE OSTI2017

Imaging Fracture Networks Using Joint Seismic and Electrical Change Detection Techniques (Final Report)

During the initial phase of this SubTER project, we conducted a series of high resolution seismic imaging campaigns designed to characterize induced fractures. Fractures were emplaced using a novel explosive source, designed at Sandia National Laboratories, that limits damage to the borehole. This work provided evidence that fracture locations could be imaged at inch scales using high-frequency seismic tomography but left many fracture properties (i.e. permeability) unresolved. We present here the results of the second phase of the project, where we developed and demonstrated emerging seismic and electrical geophysical imaging technologies that characterize 1) the 3D extent and distribution of fractures stimulated from the explosive source, 2) 3D fluid transport within the stimulated fracture network through use of a contrasting tracer, and 3) fracture attributes through advanced data analysis. Focus was placed upon advancing these technologies toward near real-time acquisition and processing in order to help provide the feedback mechanism necessary to understand and control fracture stimulation and fluid flow. Results from this study include a comprehensive set of 4D cross-hole seismic and electrical data that take advantage of change detection methodologies allowing for perturbations associated with the fracture emplacement and particulate tracer to be isolated. During the testing the team also demonstrated near real-time 4D electrical resistivity tomography imaging and 4D seismic tomography using the CASSM approach with a temporal resolution approaching 1 minute. All of the data collected were used to develop methods of estimating fracture attributes from seismic data, develop methods of assimilating disparate and transient data sets to improve fracture network imaging resolution, and advance capabilities for near real-time inversion of cross-hole tomographic data. These results are illustrated here. Advancements in these areas are relevant to all situations where fracture emplacement is used for reservoir stimulation (e.g. Enhanced Geothermal Systems (EGS) and tight shale gases).

42 ENGINEERING↗
DOE OSTI2015

Citronelle 2013 DAS VSP

The data in this archive are the distributed acoustic sensor (DAS) and well log data used in Daley, et al, 2016. Both raw and processed data are included

borehole↗
DOE OSTI2022

Subsurface stress criticality associated with fluid injection and determined using machine learning

Machine-learning methods and apparatus are disclosed to determine critical state or other parameters related to fluid-driven failure of a terrestrial locale impacted by anthropogenic activities such as hydraulic fracturing, hydrocarbon extraction, wastewater disposal, or geothermal harvesting. Acoustic emission, seismic waves, or other detectable indicators of microscopic processes are sensed. A classifier is trained using time series of microscopic data along with corresponding data of critical state or failure events. In disclosed examples, random forests and artificial neural networks are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to live data from a fluid injection locale in order to assess a frictional state, assess seismic hazard, assess permeability, make predictions regarding a future fluid-driven failure event, or drive engineering solutions for mitigation or remediation. Variations are disclosed.

Johnson, Paul Allan↗
DOE OSTI2019

Relaxation algorithms for matrix completion, with applications to seismic travel-time data interpolation

Travel time tomography is used to infer the underlying three-dimensional wavespeed structure of the Earth by fitting seismic travel time data collected at surface stations. Data interpolation and denoising techniques are important pre-processing steps that use prior knowledge about the data, including parsimony in the frequency and wavelet domains, low-rank structure of matricizations, and local smoothness. In this work, we show how local smoothness structure can be combined with low rank constraints using level-set optimization formulations, and develop a new relaxation algorithm that can efficiently solve these joint problems. In the seismology setting, we use the approach to interpolate missing stations and de-noise observed stations. The new approach is competitive with alternative algorithms, and offers new functionality to interpolate observed data using both smoothness and low rank structure in the presence of data fitting constraints.

constrained optimization↗