Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Geochemistry & Geophysics↗
DOE OSTI2023

2D Seismic Surveys as an Alternative for CO 2 Monitoring in the North Dakota CarbonSAFE Project

Conference presentation for International Meeting for Applied Geoscience & Energy (IMAGE), Houston, TX, August 28 – September 1, 2023. As part of the North Dakota CarbonSAFE project, seismic, controlled-source electromagnetic, and microgravity baseline monitoring surveys were collected to assess the safe, permanent, commercial-scale geologic storage of CO 2 generated by the Milton R. Young coal-fired power plant. Data from five 2D seismic lines were assessed as an alternative to 3D seismic for CO 2 monitoring. 2D and 3D surveys were compared using wave equation-based (WEB)-AVO inversion. The encouraging WEB-AVO results demonstrated that 2D surveys could be included in CO 2 monitoring programs as a cost-effective method.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Distributed Acoustic Sensing (DAS) on Opportunistic Networks_ A Feasibility Study Utilizing Fiber Optic Infrastructure at NEON Sites for DAS

This study demonstrates the feasibility of using Distributed Acoustic Sensing (DAS) on the opportunistic optical fiber infrastructure presented by NEON (National Ecological Observatory Network) sites. The test took place at the NEON Terrestrial Site designated D10 CPER. Unlit portions of the local fiber network were successfully spliced together and optimized for length and a Silixa iDAS sensor was installed. The resulting array (composed of 12 spliced fiber runs) totaled 875m in length. Over a period of 2 weeks, we gathered a high spatial and temporal resolution data set to characterize discernable signals and understand background noise. The data set was reviewed to extract wave signatures of identifiable origin, both anthropogenic and otherwise. Continent wide distribution of these sites makes them an interesting candidate for DAS and identifying measurable signal sources informs future DAS application on opportunistic fiber networks. These observations validate the utility of both the NEON and similarly small fiber networks (<1km) for geoscience, environmental, and security DAS applications and indicate possibilities for further inquiry in methods for DAS data collection and analysis.

42 ENGINEERING↗
DOE OSTI2021

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗
DOE OSTI2026

Economic assessment of seismic monitoring for underground hydrogen storage

Underground hydrogen storage (UHS) plays a key role in the energy landscape. However, like other subsurface engineering technologies, UHS may cause leakage into the groundwater or atmosphere and possibly induce local seismicity. To reduce these risks, seismic monitoring could be a viable technique to track the UHS plume, detect leakages, and locate induced seismicity events. Seismic monitoring has been proposed to safely monitor UHS, but research in this area is still new and requires field studies. Lab and theoretical studies have demonstrated the validity of seismic monitoring for UHS. Therefore, it is imperative to analyze the economic feasibility of seismic monitoring for UHS. Hence, we develop a cost model and open-source Python code for seismic monitoring that considers types of seismometers, comprehensive operational scenarios, detection thresholds, and long-term leakage monitoring. A case study is further provided to validate the cost model on reservoir simulations of UHS. We find that the levelized cost for a 10-year operating UHS site will range on the order of ∼0.003 $\$$/kg. The methods developed in this study could also be applied to the monitoring of groundwater, gas, and/or wastewater injection.

08 HYDROGEN↗
DOE OSTI2021

Aftershock Identification Using a Paired Neural Network Applied to Constructed Data

This report is intended to detail the findings of our investigation of the applicability of machine learning to the task of aftershock identification. The ability to automatically identify nuisance aftershock events to reduce analyst workload when searching for events of interest is an important step in improving nuclear monitoring capabilities and while waveform cross - correlation methods have proven successful, they have limitations (e.g., difficulties with spike artifacts, multiple aftershocks in the same window) that machine learning may be able to overcome. Here we apply a Paired Neural Network (PNN) to a dataset consisting of real, high quality signals added to real seismic noises in order to work with controlled, labeled data and establish a baseline of the PNN's capability to identify aftershocks. We compare to waveform cross - correlation and find that the PNN performs well, outperforming waveform cross - correlation when classifying similar waveform pairs, i.e., aftershocks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Digital Tools for the Preventive Conservation of Built Heritage: The Church of Santa Ana in Seville

Historic Building Information Modelling (HBIM) plays a pivotal role in heritage conservation endeavours, offering a robust framework for digitally documenting existing structures and supporting conservation practices. However, HBIM’s efficacy hinges upon the implementation of case-specific approaches to address the requirements and resources of each individual asset and context. This paper defines a flexible and generalisable workflow that encompasses various aspects (i.e., documentation, surveying, vulnerability assessment) to support risk-informed decision making in heritage management tailored to the peculiar conservation needs of the structure. This methodology includes an initial investigation covering historical data collection, metric and condition surveys and non-destructive testing. The second stage includes Finite Element Method (FEM) modelling and structural analysis. All data generated and processed are managed in a multi-purpose HBIM model. The methodology is tested on a relevant case study, namely, the church of Santa Ana in Seville, chosen for its historical significance, intricacy and susceptibility to seismic action. The defined level of detail of the HBIM model is sufficient to inform the structural analysis, being balanced by a more accurate representation of the alterations, through linked orthophotos and a comprehensive list of alphanumerical parameters. This ensures an adequate level of information, optimising the trade-off between model complexity, investigation time requirements, computational burden and reliability in the decision-making process. Field testing and FEM analysis provide valuable insight into the main sources of vulnerability in the building, including the connection between the tower and nave and the slenderness of the columns.

Chaves, Estefanía↗
DOE OSTI2021

The Swan: Data-driven Inference of Stellar Surface Gravities for Cool Stars from Photometric Light Curves

Stellar light curves are well known to encode physical stellar properties. Precise, automated, and computationally inexpensive methods to derive physical parameters from light curves are needed to cope with the large influx of these data from space-based missions such as Kepler and TESS. Here we present a new methodology that we call “The Swan,” a fast, generalizable, and effective approach for deriving stellar surface gravity (logg) for main-sequence, subgiant, and red giant stars from Kepler light curves using local linear regression on the full frequency content of Kepler long-cadence power spectra. With this inexpensive data-driven approach, we recover logg to a precision of ~0.02 dex for 13,822 stars with seismic logg values between 0.2 and 4.4 dex and ~0.11 dex for 4646 stars with Gaia-derived logg values between 2.3 and 4.6 dex. We further develop a signal-to-noise metric and find that granulation is difficult to detect in many cool main-sequence stars (T {sub eff} ≲ 5500 K), in particular K dwarfs. By combining our logg measurements with Gaia radii, we derive empirical masses for 4646 subgiant and main-sequence stars with a median precision of ~7%. Finally, we demonstrate that our method can be used to recover logg to a similar mean absolute deviation precision for a TESS baseline of 27 days. Our methodology can be readily applied to photometric time series observations to infer stellar surface gravities to high precision across evolutionary states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Deep Learning-Aided Workflow for Decoding the Stress Regime of Southern Nevada

The Rock Valley fault zone in southern Nevada has a notable history of seismic activity and is the site of a future direct comparison experiment of explosion and earthquake sources. This study aims to gain insight into regional tectonic processes by leveraging recent advances in seismic monitoring capabilities to elucidate the local stress regime. A crucial step in this investigation is the accurate determination of P-wave first-motion polarities, which play a vital role in resolving earthquake focal mechanisms of small earthquakes. Here, we deploy a deep learning-based method for automatic determination of first-motion polarities to vastly expand the polarity dataset beyond what has been reviewed by human analysts. By the integrating P-wave polarities with new measurements of S/P amplitude ratios, we obtain robust focal mechanism estimates for 1306 earthquakes with a local magnitude of 1 and above occurring between 2010 and 2023 in southern Nevada. We then use the focal mechanism catalog to examine the regional stress orientation, confirming an overall trans-tensional stress regime with smaller scale complexities illuminated by individual earthquake sequences. These findings demonstrate how detailed analyses of small earthquakes can provide fundamental information for understanding earthquake processes in the region and inform future experiments at the Nevada National Security Site.

58 GEOSCIENCES↗
DOE OSTI2025

Advancing Tunnel Detection Via Vertical Acoustic Profiling

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

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗
DOE OSTI2026

Survey of Deep Learning and Physics-Based Approaches in Computational Wave Imaging

Computational wave imaging (CWI) extracts hidden structure and physical properties of a volume of material by analyzing wave signals that traverse that volume. Applications include seismic exploration of the Earth’s subsurface, acoustic imaging and nondestructive testing (NDT) in material science, and ultrasound computed tomography (USCT) in medicine. Current approaches for solving CWI problems can be divided into two categories: those rooted in traditional physics and those based on deep learning. Physics-based methods stand out for their ability to provide high-resolution and quantitatively accurate estimates of acoustic properties within the medium. However, they can be computationally intensive and are susceptible to ill-posedness and nonconvexity typical of CWI problems. Machine learning (ML)-based computational methods have recently emerged, offering a different perspective to address these challenges. Diverse scientific communities have independently pursued the integration of deep learning in CWI. This review discusses how contemporary scientific ML techniques, and deep neural networks in particular, have been developed to enhance and integrate with traditional physics-based methods for solving CWI problems. We present a structured framework that consolidates existing research spanning multiple domains, including computational imaging, wave physics, and data science. This study concludes with important lessons learned from existing ML-based methods and identifies technical hurdles and emerging trends through a systematic analysis of the extensive literature on this topic.

42 ENGINEERING↗
DOE OSTI2022

On the measurement of S diff splitting caused by lowermost mantle anisotropy

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (S diff ) are potentially strongly influenced by lowermost mantle anisotropy. S diff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. S diff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SH diff –SV diff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of S diff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which S diff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on S diff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SH diff waves. In particular, we investigate how SH diff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to S diff splitting that can be clearly distinguished from other waveform effects. Here, these results enable us to lay out a strategy for the analysis of S diff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of S diff and the influence of source-side anisotropy. We demonstrate our S diff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained S diff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Imperial Valley Dark Fiber Project Continuous DAS Data

The Imperial Valley Dark Fiber Project acquired Distributed Acoustic Sensing (DAS) seismic data on a ~28 km segment of dark fiber between the cities of Calipatria and Imperial in the Imperial Valley, Southern California. Dark fiber refers to unused optical fiber cables in telecommunications networks and is repurposed in this project for DAS applications. The objective, which is further detailed in the attached journal article from Ajo-Franklin et al., is to demonstrate dark fiber DAS as a tool for basin-scale geothermal exploration and monitoring. The included DAS data were recorded during two days at the beginning the project. Data is stored in the .h5 (HDF5) file format, readable using various software tools, including the 'h5read' and 'h5info' functions in Matlab. Provided here are examples of methods to read and use the data with the 'h5py' package in Python.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

HOSS: an implementation of the combined finite-discrete element method

Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.

58 GEOSCIENCES↗
DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗