Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysics for exploration”

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

CarbonSAFE Phase II: Optimizing Alabama’s CO 2 Storage in Shelby County, Alabama (Project OASIS), Milestone M3 - Site Specific Drilling Report

The Phase II Storage Complex Feasibility project, entitled "Optimizing Alabama's CO 2 Storage in Shelby County, Alabama (Project OASIS)," is a part of the DOE/NETL's CarbonSAFE initiative. The Project is managed by the Southern States Energy Board (SSEB), and includes participation from Advanced Resources International, Inc. (ARI), Crescent Resource Innovation, Southern Company, Alabama A&M University, Auburn University, and Oklahoma State University. Project OASIS is working to establish the foundation for a commercial-scale geologic storage complex for CO 2 captured from Plant Gaston (home of the National Carbon Capture Center) and surrounding industrial sources of CO 2 located in Shelby County, Alabama. The Project objectives are: • Demonstrate that the subsurface saline formations at the storage complex can store commercial volumes of CO 2 safely and permanently. • Develop a comprehensive Community Benefits Plan. • Develop the infrastructure framework for a CO 2 storage hub. • Develop a rigorous risk registry and to conduct a comprehensive risk assessment. • Develop a monitoring plan. • Develop a comprehensive site characterization plan to support an Underground Injection Control Class VI Permit in a future Phase III program. • Evaluate the commercial viability of the project. Project OASIS is about 30 miles southeast of Birmingham, Alabama within a geologic province called the Valley and Ridge (Figure 1.1). The Valley and Ridge Province comprises a sequence of Paleozoic carbonate and clastic rocks that underwent structural deformation during the Alleghanian Orogeny. Storage prospects occur in relatively flat lying structural panels located between thrust faults. Available geologic studies related to hydrocarbon exploration suggest that Cambro-Ordovician carbonates and Cambrian clastic units offer multiple potential storage intervals, and that regional confining systems are present, such as the tectonically thickened Floyd-Parkwood Shale. The surface property is owned by a timber and land stewardship company, The Westervelt Company, Inc., who worked with the Project Team to select and prepare adequate sites for geologic assessment. The purpose of drilling the Westover Stratigraphic Test Well #2 was to collect geologic data to model the feasibility of commercial scale CO 2 injection and storage. This includes geological and geophysical evaluations, reservoir engineering analyses, and risk assessments.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Utah FORGE 5-2565: Hydrothermal Evolution of Fracture Properties - Workshop Presentation

This is a presentation on the Evolution of Permeability and Strength Recovery of Shear Fractures Under Hydrothermal Conditions project by the U.S. Geological Survey, presented by Dr. David Lockner. The project's objective was to determine how thermal, hydraulic, mechanical, and chemical processes affect the sustainability of fracture networks in geothermal reservoirs and provide strategies for improved EGS techniques that maximize thermal coupling and increase reservoir longevity. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering

Discovery of hidden geothermal resources is challenging. It requires the mining of large datasets with diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used play fairway analysis approach typically incorporates subject-matter expertise to analyze regional data to estimate geothermal characteristics and favorability. We demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset from southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. Several of these systems are not well characterized because of insufficient existing data and limited past explorative work. This study discovers hidden patterns and relations in the SWNM geothermal dataset to improve our understanding of the regional hydrothermal conditions and energy-production favorability. This understanding is obtained by applying an unsupervised ML algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden signatures characterizing analyzed datasets, (2) the optimal number of these signatures, (3) the dominant data attributes associated with each signature, and (4) the spatial distribution of the extracted signatures. Here, in this study, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, and geothermal attributes at 44 locations in SWNM. Using NMFk, we find data patterns and identify the spatial associations of hydrothermal signatures within two physiographic provinces (Colorado Plateau and Basin and Range) and two sub-regions of these provinces (the Mogollon-Datil volcanic field and the Rio Grande rift) in SWNM. The ML algorithm extracted five hydrothermal signatures in the SWNM datasets that differentiate between low (<90°C) and medium (90-150°C)-temperature hydrothermal systems. The algorithm also suggests that the Rio Grande rift and northern Mogollon-Datil volcanic field are the most favorable regions for future geothermal resource discovery. NMFk also identified critical attributes to identify medium-temperature hydrothermal systems in the study area. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored regions. The code to execute the performed analyses as well as the corresponding data can be found at https://github.com/SmartTensors/GeoThermalCloud.jl.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗
DOE OSTI2026

FY 2026 Midyear Report: Seismic Monitoring of Underground Vibration Sources Using Distributed Acoustic Sensing and Seismometers

Safeguards-relevant temporal changes in underground facilities can be observed using geophysical monitoring techniques. Seismic waves, in particular, provide valuable insights into subsurface activities and can serve as an important tool for detecting anomalous events that may indicate containment breaches at geological repositories. This midyear report summarizes ongoing efforts to automatically and rapidly detect and locate anomalous vibration signals that could be indicative of potential containment breaches. Previous work during FY25 focused on compiling continuous seismic datasets from two underground sites and developing a database of continuous waveforms and ground-truth event data derived from multiple sensing modalities. Building on this foundation, we are adapting anomaly detection and geolocation algorithms to explore methods for monitoring underground activities using two relatively low-maintenance sensing technologies: a dense surface geophone array deployed at the Pleasant Gap mine in Pennsylvania, and a three-dimensional fiber-optic cable array for distributed acoustic sensing (DAS) installed in the subsurface at the Sanford Underground Research Facility (SURF) in South Dakota. This report summarizes work conducted during the first two quarters of FY26, during which we refined a dynamic power spectral density (PSD)-based detector, applied it independently to each geophone station, and then combined the per‑station detections with density-based spatial clustering of applications with noise (DBSCAN) to cluster events and produce spatial maps over a nine‑day interval. In addition, we outline plans for a field trial at the Waste Isolation Pilot Plant (WIPP) in New Mexico to compare traditional seismic monitoring approaches with DAS techniques and to evaluate the benefits of combined data analysis. Activities during the past two quarters have included the preparation and submission of a Field Test Plan to WIPP for approval, as well as submission to headquarters for review and feedback.

58 GEOSCIENCES↗
DOE OSTI2026

Geophysical Impacts and Spectroscopic Identification of a Hydrous Iron Sulfate on Icy Worlds

Over geologic time-scales, large volumes of exogenic sulfur ions from Io's plasma torus have been supplied to the surface of Europa and Ganymede, which, combined with recent interpretations of orbiter images, dynamical modeling, and surface-subsurface exchange, suggests further sulfur transport into the interior of the icy worlds. These observations motivate mixed-phase spectral modeling for interpreting orbiter spectroscopy data and determination of hydration states of candidate surface materials including hydrous sulfates. In this work, we present a combined experimental and theoretical study of the low temperature and high pressure vibrational spectral signature of the iron-sulfate monohydrate endmember, szomolnokite (FeSO 4 ·H 2 O). By employing synchrotron Fourier-transform infrared spectroscopy (FTIR) in the diamond anvil cell up to 23 GPa and down to 20 K, we explore the extreme range of pressure-temperature domains relevant to icy environments throughout our solar system and beyond. Combined with our density-functional theory quantum-mechanics molecular dynamics results, we demonstrate that experimentally observed infrared features in the O-H stretching region commonly associated with nH 2 O (n > 1) hydration states can be attributed to a pure monohydrate without the need for pressure-induced exsolved ice, other coexisting hydrous iron sulfates, or strong overtone and combination modes. We further discuss the possibility of lateral variations in density and shear properties on icy worlds associated with temperature variations and the high-pressure phases of kieserite group monohydrated sulfates.

Geosciences↗
DOE OSTI2023

Exploratory analysis of machine learning techniques in the Nevada geothermal play fairway analysis

Play fairway analysis (PFA) is commonly used to generate geothermal potential maps and guide exploration studies, with a particular focus on locating and characterizing blind geothermal systems. This study evaluates the application of machine learning techniques to PFA in the Great Basin region of Nevada. Following the evaluation of various techniques, we identified two approaches to PFA that produced promising results, 1) supervised Bayesian probabilistic neural networks to generate geothermal potential maps with confidence intervals, and 2) unsupervised principal component analysis paired with k-means clustering to generate both cluster maps to help identify spatial patterns, as well as new combined feature inputs. We applied these techniques to perform a comparative analysis between two principal sets of geological and geophysical features related to permeability and heat and a set of positive (known geothermal resources) and negative training sites (known drill sites with unsuitable geothermal conditions). We found that these methods constrain previously unrecognized feature controls on geothermal favorability, many of which are spatially organized within the extent of cluster groups and the major structural-hydrologic domains of the study area. Furthermore, we utilized exploratory unsupervised modeling to highlight spatial relationships between input data and predictive output results of our supervised modeling. As a result, we demonstrate how our models compare to the previous Nevada PFA and how the rapid insights these machine learning techniques offer may support future assessments of both known and undiscovered blind geothermal systems in the Great Basin region of Nevada and beyond.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Utah FORGE 3-2417: Fiber-Optic Geophysical Monitoring of Reservoir Evolution - Workshop Presentation

This is a presentation on the Fiber-Optic Geophysical Monitoring of Reservoir Evolution at the Utah FORGE Milford Site project by Rice University, presented by Prof. Jonathan Ajo-Franklin. The project's objective was to develop an end-to-end fiber-optic sensing approach for EGS to track the initial zone of fracture creation, zones of connected mechanically compliant fractures, and zones of flowing fractures linking the injection/production well pair. This approach would also be used to integrate data into an improved thermo-hydro-mechanical model (THM). This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Enhanced Multi-Dimensional Inversion Through Target-Specific Inversion Parameter Bounds With an Application to Crosswell Electromagnetic for Sequestration Monitoring

In geophysical inversions, lower and upper model parameter bounds are a means of solution stabilization. Further, constraints that intend to let only geologically plausible inverse solutions pass are amenable to lower and upper bounds. Reliable prior information is paramount to construct such bound constraints. It is common practice to narrow and widen bound intervals for regions of, respectively, more and less certain prior information. Contrary to this practice, we experiment with widened bound intervals in zones that are poorly resolved by a given survey configuration but where prior information would suggest structural anomalies of interest. The purpose of enlarged parameter bounds that correlate spatially with predefined targets is to let the inversion explore a larger solution space, thus increasing the potential to resolve otherwise hidden anomalies. Application of the method is based on a carbon-sequestration baseline (pre-injection) crosswell electromagnetic (EM) field survey at the Containment and Monitoring Institute Field Research Station (Alberta, Canada), where impeded measurements led to generally reduced sensitivities for the interwell region. Synthetic-data proofs of concept use augmented bounds designed to boost the resolution of artificial plume targets, indicating an enhanced illumination compared to constant bounds. Comparative field data inversions with spatially variable bounds constructed from prior resistivity and velocity information highlight non-horizontal baseline structures.

3-D electrical resistivity imaging↗
DOE OSTI2021

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗
DOE OSTI2024

Seismic Monitoring near Ithaca, New York, Reveals Nonuniform Distribution of Microseismicity in an Intraplate Region

Abstract Cornell University intends to use a deep direct-use geothermal system to heat its Ithaca, New York, campus. In preparation for this project, the Cornell Seismic Network has been monitoring the background seismicity in this intraplate region since 2019. From January 2020 to June 2023, 95 events were detected within 20 km of the proposed geothermal well site, with local magnitudes ranging from −1.02 to 0.56. None of these events appear in regional or national catalogs. Events locate in a narrow geographic band, with one-fourth exhibiting multimodal hypocentral probability peaks both near the surface and at 1–4 km depth. We relocate events with a joint hypocenter and 1D velocity model inversion, in addition to a fully nonlinear method, and then compare observations with synthetic waveforms. Together, these approaches provide strong evidence for >95% of events locating at the surface or within the 3-km-thick sedimentary sequence. We explore how anthropogenic activity and regional topographic stress may contribute to frequent surficial events. This information is critical for characterizing the background microseismicity for comparison during future geothermal operations. Ithaca’s geology of Paleozoic sediments overlying Precambrian crystalline basement is typical of many continental interiors, so these results also provide insight into intraplate microseismicity patterns.

Geochemistry & Geophysics↗
DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Electromagnetic Induction (EMI) Data, 2024, Trail Creek, Colorado

This dataset contains Electromagnetic Induction (EMI) data collected at Trail Creek, Colorado, in 2024. EMI surveys were conducted to investigate the spatial distribution of electrical conductivity in the subsurface, providing insights into soil moisture and subsurface geological features. The surveys were performed along multiple transects to capture variations in conductivity influenced by changes in soil composition, moisture content, and underlying geological structures. This dataset complements other geophysical data collected in the region, including Electrical Resistivity Tomography (ERT) and Terrestrial LiDAR Scanning (TLS), providing a detailed understanding of the subsurface and its impact on surface vegetation and hydrological processes. The data are valuable for environmental geophysics, ecological research, and hydrological modeling in mountainous ecosystems. The files include: - data.zip: the raw EMI data (.csv) - inversion.zip: the inverted resistivity model (.csv and .kml) - kriging.zip: the kriging resistivity model (.csv, .tif, .kmz) - flmd.csv: file level metadata file describing all files within this dataset - dd.csv: data dictionary file describing the column headers within CSV files This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

CMD Mini-Explorer↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

High-pressure behavior of 3.65 Å phase: Insights from Raman spectroscopy

Abstract The 3.65 Å phase [MgSi(OH)6] is a hydrous phase that is predicted to be stable in a simplified MgO-SiO2-H2O (MSH) ternary system at pressures exceeding 9 GPa. Along cold subduction zones, it is likely to transport water, bound in its crystalline lattice, into the Earth’s interior. The 3.65 Å phase consists of Mg and Si octahedral sites attached to the hydroxyl group that forms a hydrogen bond and is predicted to undergo pressure-induced symmetrization of the hydrogen bond. Therefore, in this study, we investigate the high-pressure behavior of the 3.65 Å phase using Raman spectroscopy. We have conducted five distinct compressions up to ~60 GPa using two different pressure-transmitting media—alcohol mixture and neon. At ambient conditions, we identified vibrational modes using complementary first-principles simulations based on density functional perturbation theory. Upon compression, we note that the first derivative of the vibrational modes in the lattice region stiffens, i.e., b1lattice > 0. In contrast, the hydroxyl region softens, i.e., b1OH > 0. This is indicative of the strengthening of hydrogen bonding upon compression. We noticed a significant broadening of vibrational modes related to hydroxyl groups that are indicative of proton disorder. However, within the maximum pressures explored in this study, we did not find evidence for pressure-induced symmetrization of the hydrogen bonds. We used the pressure derivative of the vibrational modes to determine the ratio of the bulk moduli and their pressure derivative. We note that the smaller bulk moduli of hydrous phases compared to the major mantle phases are compensated by significantly larger pressure derivatives of the bulk moduli for the hydrous phases. This leads to a significant reduction in the elasticity contrast between hydrous and major mantle phases. Consequently, the detection of the degree of mantle hydration is likely to be challenging at greater depths.

Geochemistry & Geophysics↗
DOE OSTI2025

Nonlinear Ensemble Filtering with Diffusion Models: Application to the Surface Quasigeostrophic Dynamics

The intersection between classical data assimilation methods and novel machine learning techniques has attracted significant interest in recent years. Here, we explore another promising solution in which diffusion models are used to formulate a robust nonlinear ensemble filter for sequential data assimilation. Unlike standard machine learning methods, the proposed ensemble score filter (EnSF) is completely training free and can efficiently generate a set of analysis ensemble members. Here, in this study, we apply the EnSF to a surface quasigeostrophic model and compare its performance against the popular local ensemble transform Kalman filter (LETKF), which makes Gaussian assumptions in the analysis step. Numerical tests demonstrate that EnSF maintains stable performance in the absence of localization and for a variety of experimental settings. We find that while LETKF maintains optimal performance in the case of linear observations of the entire state and a perfect model, EnSF shows improvements over LETKF when nonlinear observations are assimilated and the system is subject to unexpected model errors. A spectral decomposition of the analysis results in this nonlinear observation regime shows that the largest improvements over LETKF occur at large scales (small wavenumbers), where LETKF lacks sufficient ensemble spread. Overall, this initial application of EnSF to a geophysical model of intermediate complexity motivates further development of the algorithm for more realistic problems.

Artificial intelligence↗
DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗
DOE OSTI2021

Data Arrays for Microearthquake (MEQ) Monitoring using Deep Learning for the Newberry EGS Sites

The 'Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties' project looks to apply machine learning (ML) methods to Microearthquake (MEQ) data for imaging geothermal reservoir properties and forecasting seismic events, in order to advance geothermal exploration and safe geothermal energy production. As part of the project, this submission provides data arrays for 149 microearthquakes between the year 2012 and 2013 at the Newberry EGS Site for use with the Deep Learning Algorithm that has been developed. The data provided includes raw waveform data, location data, normalized waveform data, and processed waveform data. Penn State Geothermal Team has shared the following files from the project: - 149 microearthquakes (MEQs) between 2012 and 2013 at Newberry EGS sites, 'Normalized Waveform Inputs.npz' are normalized waveforms. - labels of 149 MEQs: Processed Waveform Inputs.npz - location labels of 149 MEQs: Location Data.npz Note: .npz is the python file format by NumPy that provides storage of array data.

15 GEOTHERMAL ENERGY↗