Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface system”

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 757 records · Page 42

On the source of cross-grain lineations in the central Pacific gravity field

The source of cross-grain lineations in marine gravity field observed in central Pacific was investigated by comparing multiple collinear gravity profiles from Geosat data with coincident bathymetry profiles, in the Fourier transform domain. Bathymetric data were collected by multibeam sonar systems operating from two research vessels, one in June-August 1985, the other in February and March 1987. The results of this analysis indicate that the lineations are superficial features that appear to result from a combination of subsurface and surface loads supported by a thin (2 km to 5 km) lithosphere.

Mcadoo, David C.↗

Hawaii Play Fairway Analysis: Lanai Resistivity and Density 3D Inversion Models

To prepare for its third phase, the Hawaii Play Fairway project conducted groundwater sampling and analyses in ten locations in the Hawaiian islands, magnetotelluric (MT) and gravity surveys, as well as calculations of 3D subsurface stress due to the weight of the rock underlying the topography of the volcano. The subsurface stresses were used to evaluate the potential for fracture-induced permeability. Inversions of the MT and gravity data produce 3D models of resistivity and density, respectively, on Lanai, across Haleakala's SW rift (Maui), and surrounding Mauna Kea (Hawaii Island). The project developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The project then incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. These results were used to identify target sites for exploratory drilling. Spreadsheet information: Each sheet contains data for a particular depth in kilometers. Positive depths are above sea level, and negative below. For more information, go to the Hawaii Groundwater and Geothermal Resources Center website linked in the resources.

15 GEOTHERMAL ENERGY↗

Data for Machado-Silva et al. (2024), "Short-Term Groundwater Level Fluctuations Drive Subsurface Redox Variability"

This dataset contains the analytical data reported in Machado-Silva et al. (2024) as part of the COMPASS-FME project, which seeks to advance a scalable, predictive understanding of the fundamental biogeochemical processes, ecological structure, and ecosystem dynamics that distinguish coastal terrestrial-aquatic interfaces from the purely terrestrial or aquatic systems to which they are coupled. The dataset consists of water quality parameters as well as redox potential, water content, and electrical conductivity. These data were collected in 2022 in Crane Creek (CRC), Portage River (PTR), and Old Woman Creek (OWC). Each of these sites included uplands (UP), transitions (TR), wetland-transition edge (WTE), and wetland (W) zones. The sites represent replicates of the Lake Erie terrestrial-aquatic interface under fluctuating water levels and are located in well-preserved areas with natural or restored marsh and forest cover.This dataset consists of a single data file (Machado_Silva_et_al_2024_EST_data.csv) that is in comma-separated value (CSV) format. No special software is required to read it.This dataset uses the ESS-DIVE Hydrologic Monitoring Reporting Format 1.0.

54 ENVIRONMENTAL SCIENCES↗

Inland and Near Shore Water Profiles Derived from the High Altitude Multiple Altimeter Beam Experimental Lidar (MABEL)

The Advanced Topographic Laser Altimeter System (ATLAS) on the Ice, Cloud, and Land Elevation Satellite (ICESat-2) mission is a six beam, low energy, high repetition rate, 532 nm laser transmitter with photon counting detectors. Although designed primarily for detecting height changes in icecaps, sea ice and vegetation, the polar-orbital satellite will observe global surface water during its designed three year life span, including inland water bodies, coasts, and open oceans. In preparation for the mission, an ICESat-2 prototype or the Multiple Altimeter Beam Experimental Lidar (MABEL), was built and flown on high altitude aircraft experiments over a range of inland and near-shore targets. The purpose was to test the ATLAS concept and to provide a database for developing an algorithm that detects along track surface water height and light penetration under a range of atmospheric and water conditions. The current analysis examines the datasets of three MABEL transects observed from 20 km above ground of coastal and inland waters conducted in 2012 and 2013. Transects ranged from about 2 to 12 km in length and included the middle Chesapeake Bay, the near shore Atlantic coast at Virginia Beach, and Lake Mead. Results indicate MABEL's high capability for retrieving surface water height statistics with a mean height precision of approximately 5-7 cm per 100m segment length. Profiles of attenuated subsurface backscatter, characterized using a Signal to Background Ratio written in Log10 base, or LSBR0, were observed over a range of 1.3 to 9.3 meters depending on water clarity and atmospheric background. Results indicate that observable penetration depth, although primarily dependent on water properties, was greatest when solar background rate was low. Near shore bottom reflectance was detected only at the Lake Mead site down to maximum of 10 m under a clear night sky and low turbidity of approximately 1.6 Nephelometric Turbidity Units (NTU). The overall results suggest that the feasibility of retrieving operational surface water height statistics from space-based photon counting systems such as ATLAS is very high for resolutions down to about 100m, even in partly cloudy conditions. The capability to observe subsurface backscatter profiles is achievable but requires much longer transects of several hundreds of meters.

Topographic↗

Maps of growing season gross primary production and net ecosystem exchange for Council Road Mile Marker 71, Seward Peninsula, Alaska, [2017-2023]

This data archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025a). Murphy et al. (2025a) evaluated whether incorporating observed Arctic vegetation heterogeneity into ELM, the land model of the Department of Energy’s Energy Exascale Earth System Model (E3SM), improved simulations of tundra carbon cycling. The associated model archive can be found at Murphy et al. (2025b). The study focused on the spatial patterns and net landscape-level growing season productivity and carbon uptake. As part of this evaluation, observationally derived maps of average growing season (June–August) net ecosystem exchange (NEE) and gross primary production (GPP) were developed for the same domain. These maps, which form the dataset described here, integrate eddy covariance flux tower, remote sensing, and vegetation community data to provide spatially explicit benchmarks for model evaluation. The maps provide spatially explicit estimates of average growing season NEE and GPP across 13 tundra vegetation communities within the study domain. By combining flux tower observations with Airborne Visible-Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) hyperspectral imagery and drone-based normalized difference vegetation index (NDVI), these maps capture the heterogeneity of carbon fluxes associated with different Arctic vegetation types. While they represent average seasonal conditions rather than interannual variability, the maps provide a unique dataset for evaluating model performance, comparing vegetation community contributions to landscape-scale carbon cycling, and supporting regional analyses of Arctic carbon dynamics. This data archive contains 5 m resolution maps of vegetation communities, vegetation community average growing season GPP, and vegetation community average growing season NEE (three *.tif files), a User’s Guide (*pdf file), and Table 1 of the User’s Guide displaying vegetation community coverage and average growing season NEE and GPP values (*.csv file).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗

Diverse and unconventional methanogens, methanotrophs, and methylotrophs in metagenome-assembled genomes from subsurface sediments of the Slate River floodplain, Crested Butte, CO, USA

We use metagenome-assembled genomes (MAGs) to understand single-carbon (C1) compound-cycling—particularly methane-cycling—microorganisms in montane riparian floodplain sediments. We generated 1,233 MAGs (>50% completeness and <10% contamination) from 50- to 150-cm depth below the sediment surface capturing the transition between oxic, unsaturated sediments and anoxic, saturated sediments in the Slate River (SR) floodplain (Crested Butte, CO, USA). We recovered genomes of putative methanogens, methanotrophs, and methylotrophs (n = 57). Methanogens, found only in deep, anoxic depths at SR, originate from three different clades (Methanoregulaceae, Methanotrichaceae, and Methanomassiliicoccales), each with a different methanogenesis pathway; putative methanotrophic MAGs originate from within the Archaea (Candidatus Methanoperedens) in anoxic depths and uncultured bacteria (Ca. Binatia) in oxic depths. Genomes for canonical aerobic methanotrophs were not recovered. Ca. Methanoperedens were exceptionally abundant (~1,400× coverage, >50% abundance in the MAG library) in one sample that also contained aceticlastic methanogens, indicating a potential C1/methane-cycling hotspot. Ca. Methylomirabilis MAGs from SR encode pathways for methylotrophy but do not harbor methane monooxygenase or nitrogen reduction genes. Comparative genomic analysis supports that one clade within the Ca. Methylomirabilis genus is not methanotrophic. The genetic potential for methylotrophy was widespread, with over 10% and 19% of SR MAGs encoding a methanol dehydrogenase or substrate-specific methyltransferase, respectively. MAGs from uncultured Thermoplasmata archaea in the Ca. Gimiplasmatales (UBA10834) contain pathways that may allow for anaerobic methylotrophic acetogenesis. Overall, MAGs from SR floodplain sediments reveal a potential for methane production and consumption in the system and a robust potential for methylotrophy.

58 GEOSCIENCES↗

A flow-through gamma detection system for sub-surface radioxenon tracer measurements in LYNM PE1

A continuous flow-through, real-time radioxenon sensor array has been designed to monitor either 127 Xe or 133 Xe tracer gas as part of a series of underground chemical explosions in the Low-Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1). The sensors are pairs of shielded thallium-doped sodium iodide gamma-ray detectors that surround a gas measurement chamber. Each measurement chamber is connected to a sampling location in the geology surrounding the chemical explosion or in the access tunnel. The first experiment in the series, PE1 A, occurred in 2023 and the system quantified the 127 Xe tracer gas in the gas sampling lines for four weeks following the experiment execution. The array generated a set of radioxenon concentration time series with a 15-minute resolution that will be used to refine subsurface gas transport models.

Gamma-ray coincidence detection↗

Achieving Ultra-High Selectivity to Hydrogen Production from Formic Acid on Pd–Ag Alloys

Palladium-silver based alloy catalysts have a great potential for CO-free hydrogen production from formic acid for fuel cell applications. However, the structural factors affecting the selectivity of formic acid decomposition is still debated. Herein, the decomposition pathways of formic acid on Pd-Ag alloys with different atomic configurations have been investigated to identify the alloy structures yielding high H 2 selectively. Several Pd x Ag 1-x surface alloys with various compositions were generated on a Pd(111) single crystal; their atomic distribution and electronic structure were determined by a combination of infrared reflection absorption spectroscopy (IRAS), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT). It was established that the Ag atoms with Pd neighbors are electronically altered, and the degree of alteration correlates with the number of nearest Pd. Temperature programmed reaction spectroscopy (TPRS) and DFT demonstrated that the electronically altered Ag domains create a new reaction pathway that selectively dehydrogenates formic acid. In contrast, Pd monomers surrounded by Ag are demonstrated to have a similar reactivity compared to pristine Pd(111), yielding CO and H 2 O in addition to the dehydrogenation products. Furthermore, they bind to the produced CO weaker than pristine Pd, demonstrating an enhancement in resistance to CO poisoning. This work therefore shows that surface Ag domains modified by interaction with subsurface Pd are the key active sites for selective decomposition of formic acid, while surface Pd atoms are detrimental for selectivity. Hence the decomposition pathways can be tailored for CO-free H 2 production on Pd-Ag alloy systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Software Developed for Analyzing High- Speed Rolling-Element Bearings

COBRA-AHS (Computer Optimized Ball & Roller Bearing Analysis--Advanced High Speed, J.V. Poplawski & Associates, Bethlehem, PA) is used for the design and analysis of rolling element bearings operating at high speeds under complex mechanical and thermal loading. The code estimates bearing fatigue life by calculating three-dimensional subsurface stress fields developed within the bearing raceways. It provides a state-of-the-art interactive design environment for bearing engineers within a single easy-to-use design-analysis package. The code analyzes flexible or rigid shaft systems containing up to five bearings acted upon by radial, thrust, and moment loads in 5 degrees of freedom. Bearing types include high-speed ball, cylindrical roller, and tapered roller bearings. COBRA-AHS is the first major upgrade in 30 years of such commercially available bearing software. The upgrade was developed under a Small Business Innovation Research contract from the NASA Glenn Research Center, and incorporates the results of 30 years of NASA and industry bearing research and technology.

Fleming, David P.↗

Unveiling Clues from Spacecraft Missions to Comets and Asteroids through Impact Experiments

The Deep Impact Spacecraft mission was the first to boldly face the challenge of impacting the surface of a comet, 9P/Tempel 1, to investigate surface and subsurface 'pristine' materials. The Stardust mission to Comet 81P/Wild 2 brought back an exciting surprise: shocked minerals which were likely altered during the comet's lifetime. Signatures of shock in meteorites also suggest that the violent past of the solar system has left our small bodies with signatures of impacts and collisions. These results have led to the question: How have impacts affected the evolutionary path taken by comets and asteroids, and what signatures can be observed? A future planetary mission to a near-Earth asteroid is proposing to take the next steps toward understanding small bodies through impacts. The mission would combine an ESA led AIM (Asteroid Impact Mission) with a JHU/APL led DART (Double Asteroid Redirect Mission) spacecraft to rendezvous with binary near-Earth asteroid 65803 Didymus (1996 G2). DART would impact the smaller asteroid, 'Didymoon' while AIM would characterize the impact and the larger Didymus asteroid. With these missions in mind, a suite of experiments have been conducted at the Experimental Impact Laboratory (EIL) at NASA Johnson Space Center to investigate the effects that collisions may have on comets and asteroids. With the new capability of the vertical gun to cool targets in the chamber through the use of a cold jacket fed by liquid nitrogen, the effects of target temperature have been the focus of recent studies. Mg-rich forsterite and enstatite (orthopyroxene), diopside (monoclinic pyroxene) and magnesite (Mg-rich carbonate) were impacted. Target temperatures ranged from 25 deg to -100 deg, monitored by connecting thermocouples to the target container. Impacted targets were analyzed with a Fourier Transform Infrared Spectrometer (FTIR) and Transmission Electron Microscope (TEM). Here we present the evidence for impact-induced shock in the minerals through both spectra and TEM imaging and compare with unshocked samples.

Lederer, Susan M.↗

EERC Final Topical Report – Findings on Subtask 3.2 – Produced Water Management Through Geologic Homogenization, Conditions, and Reuse

The Energy & Environmental Research Center (EERC) was awarded a contract by the North Dakota Industrial Commission (NDIC) Oil and Gas Research Program (NDIC No. G-051-101) to conduct a study on the recycling of water used in oil and gas operations, also known as produced water, from oil- and gas-producing regions of North Dakota as directed by Section 19 of North Dakota House Bill 1014. This final report provides a compilation of results of the study, which include regulatory, scientific, technological, and feasibility methods and considerations associated with North Dakota produced water management. The report provides a synopsis of this project’s previously submitted produced water assessment report entitled “Produced Water Management and Recycling Options in North Dakota” (Energy & Environmental Research Center, 2020), with updated values provided as appropriate. The report provides the results from the investigation of a novel produced water management strategy, referred to as geologic homogenization, conditioning, and reuse (GHCR), which aims to use a subsurface geologic formation as a natural medium for managing produced water recycling and reuse. Water management is a significant technical and economic challenge for sustainable oil and gas production, and water volumes are intrinsically linked to oil production volumes. North Dakota oil production rose to over 1.5 million barrels (MMbbl)/day in 2019, and despite a downturn in oil price in early 2020, North Dakota oil production has recovered to 1.1 MMbbl/day as of August 2021. Bakken petroleum system development between 2008 and 2020 has resulted in a nearly fourfold increase in produced water volumes to 642 MMbbl/yr in 2020 after peaking at 740 MMbbl/yr in 2019 and a fivefold increase in saltwater disposal (SWD) volumes to 565 MMbbl/yr in 2020 after a peak of 682 MMbbl/yr in 2019. Produced water and SWD volumes are forecasted to double by 2030. SWD is the primary method of produced water management used in North Dakota, with approximately 95% of the SWD volume occurring through subsurface injection into sandstones of the Dakota Group (Dakota). Localized pressurization of the Dakota resulting from SWD and projected increases in produced water volumes could impact the economics of North Dakota oil production. As a result, there is an emerging need to pursue alternative produced water management approaches, including recycling and reuse. While produced water recycling is not yet widespread, commercial operators are making strides in overcoming the technical challenges of using high salinity produced water in completion operations (Marathon Oil, 2020). As water management continues to be a key focal point in companies’ environmental, social, and governance (ESG) initiatives, focus on water management, including recycling, will likely continue to increase. Laboratory column testing, field sample collection, geologic modeling and numerical simulation, and techno-economic analysis all indicate that GHCR could feasibly be implemented as a potential water management option. Laboratory column testing and field sample collection indicate that the Inyan Kara sandstone and native formation fluid are capable of homogenizing with the Bakken produced water to a point where the fluid composition appears to stabilize. Extracting that stabilized fluid could be considered homogeneous and capable of providing individual batches of hydraulic fracturing fluid. Numerical simulation results indicate that extraction of fluids from the Inyan Kara in a GHCR implementation scenario is capable of reducing formation pressure, which would help ease localized pressurization of the Inyan Kara and extend the available capacity for nearby existing SWD wells. Economic analysis indicates that there are scenarios where GHCR implementation can be a competitive or even lower-cost option than a conventional water management approach. Site-specific conditions will dictate the economic potential of GHCR, but potentially attractive sites for GHCR implementation will be those that are located above a pressurized zone of the Inyan Kara, need six or more Bakken infill wells, and face high costs for conventional SWD and/or freshwater. Based on the regulatory review, drilling into the Inyan Kara for SWD and to harness as a source water for industrial use have precedent in North Dakota, and a workable regulatory solution for GHCR seems likely. However, restrictions in the state regarding surface storage and transport of produced fluids may limit some activities, which will affect how GHCR could ultimately be implemented. In summary, this study reveals pursuing GHCR can be a viable approach to water management in North Dakota. The GHCR concept addresses some of the challenges that hinder the more traditional approaches to recycling in the industry. Furthermore, an assessment of the current landscape of water management within the state reveals the ongoing trend of increasing volumes of produced water and SWD. Projections reveal that the volumes of produced water that need to be managed are expected to double over the next decade (Energy & Environmental Research Center, 2020). With the continued development of the Bakken and continuing driving factors related to ESG initiatives, implementing a practice such as GHCR is a feasible approach to adding recycling of produced water to industry within the state. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission Oil and Gas Research Program. References Energy & Environmental Research Center, 2020, Produced water management and recycling options in North Dakota: Final Report for North Dakota Legislative Management Energy Development and Transmission Committee and North Dakota Industrial Commission. Marathon Oil, 2020, Sustainability report: https://cdn.sanity.io/files/ghcnw9z2/website/ 91744eb6ef8fbe59505a911c6b8d2e8dd9a537fa.pdf?dl (accessed November 2021).

02 PETROLEUM↗

Geologic Assessment and Recommended New Well Locations to Improve the Area 25 Potable Water System, Nevada National Security Site

The Mission Support and Test Services, LLC (MSTS) Geoscience Operations Group performed a hydrogeologic study of an area of interest in Area 25 at the request of the MSTS Enterprise Infrastructure Program to identify potential locations for a new water well to replace J-14 Water Well (J-14 WW). Based on limited subsurface hydrogeologic data for the area, geologists developed two options for construction of a new well. Option 1 is to construct the well within approximately 250 feet (ft) of J-14 WW, and Option 2 is to construct the well approximately 1,550 ft north of J-11 Water Well (J-11 WW). Option 1 may be preferred, as the hydrogeology is better understood adjacent to an existing well. A well drilled at the Option 2 location would not have to be as deep as that at the Option 1 location, but the hydrogeologic setting is more uncertain at this greater distance from a verified subsurface geologic data point.

58 GEOSCIENCES↗

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape: Modeling Archive

This modeling archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025). This archive contains model input files and outputs from landscape-scale simulations conducted using ELM, the land model component of the Department of Energy’s Energy Exascale Earth System Model (E3SM), at the Council NGEE Arctic field site (Council Road mile marker 71) on Alaska’s Seward Peninsula. Input data and model output from two sets of ELM simulations are provided. The first set of simulations were conducted with the two default ELM Arctic plant functional types (PFTs; broadleaf deciduous boreal shrub and a C3 grass) and the second set of simulations were conducted with a set of nine Arctic-specific PFTs including nonvascular mosses and lichens, graminoids, forbs, evergreen dwarf shrubs, three height classes of deciduous shrubs (dwarf, low, and low to tall), and deciduous alder shrubs (Sulman et al., 2021). Parameter names and major parameter changes in the Arctic-specific PFT configuration are described in Sulman et al. (2021) and archived in the Sulman et al. (2021) dataset (see below). Simulations were spatially explicit, covering an approximately 6.4X3.3 km domain at the Council site with a spatial resolution of 100 m for a total of 2,112 simulated grid cells under each ELM PFT configuration. The modeling archive contains meteorological forcing (seven *.nc files and one *.txt file), a domain definition file (one *.nc files), land surface configuration files (two *.nc files), parameter files (two *.nc files), annual ELM output files spanning 1980-2014 (68 *.nc files), and a User’s Guide (*pdf file). Additional information on the provided files is in the “Modeling Archive Contents” section of the User’s Guide. Model outputs are aggregated to the column scale (i.e. PFT-specific outputs are not provided here).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗

Inland and Near-Shore Water Profiles Derived from the High-Altitude Multiple Altimeter Beam Experimental Lidar (MABEL)

The Advanced Topographic Laser Altimeter System (ATLAS) on the Ice, Cloud, and Land Elevation Satellite (ICESat-2) mission is a six beam, low energy, high repetition rate, 532-nanometer laser transmitter with photon counting detectors. Although designed primarily for detecting height changes in ice caps, sea ice, and vegetation, the polar-orbiting satellite will observe global surface water during its designed three-year life span, including inland waterbodies, coasts, and open oceans. In preparation for the mission, an ICESat-2 prototype, the Multiple Altimeter Beam Experimental Lidar (MABEL), was built and flown on high-altitude aircraft experiments over a range of inland and near-shore targets. The purpose was to test the ATLAS concept and to provide a database for developing an algorithm that detects along track surface water height and light penetration under a range of atmospheric and water conditions. The current analysis examines the data sets of three MABEL transects observed from 20 kilometers above ground of coastal and inland waters conducted in 2012 and 2013. Transects ranged from about 2 to 12 kilometers in length and included the middle Chesapeake Bay, the near-shore Atlantic coast at Virginia Beach, and Lake Mead. Results indicate MABEL's high capability for retrieving surface water height statistics with a mean height precision ofapproximately 5-7 centimeters per 100-meter segment length. Profiles of attenuated subsurface backscatter, characterized using a Signal to Background Ratio written in Log10 base, or LSBR (sub 0), were observed over a range of 1.3 to 9.3 meters, depending on water clarity and atmospheric background. Results indicate that observable penetration depth, although primarily dependent on water properties, was greatest when the solar background rate was low. Near-shore bottom reflectance was detected only at the Lake Mead site down to a maximum of 10 meters under a clear night sky and low turbidity of approximately 1.6 Nephelometric Turbidity Units (NTU). The overall results suggest that the feasibility of retrieving operational surface water height statistics from space-based photon counting systems such as ATLAS is very high for resolutions down to about 100 meters, even in partly cloudy conditions. The capability to observe subsurface backscatterprofiles is achievable but requires much longer transects of several hundreds of meters.

inland↗

Fiber-Optic Sensing for Earthquake Hazards Research, Monitoring, and Early Warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

58 GEOSCIENCES↗

High-Resolution 3D Acoustic Borehole Integrity Monitoring (Research Performance Progress Report Q2 2021)

Real-time, in-situ, high spatial resolution (sub-cm) imaging of the near-borehole environment would revolutionize wellbore diagnostics and integrity assessment by direct observation of defects. It is becoming increasingly apparent, that better understanding of the near-wellbore environment is required to meet the safety and operational needs in challenging environments such as those present in subsurface energy extraction (geothermal) and storage (CO2 sequestration) applications. Therefore, it is important to have a more robust ability to image the near-borehole and reliably detect defects. It was proposed to further develop and improve our advanced 3D imaging system to evaluate casing defects (e.g. corrosion) and cement quality in either open- or cased-borehole with the ultimate goal to develop a commercially deployable technology. The system consists of a unique acoustic source (LANL) and advanced inversion techniques for image processing (LANL, ORNL). This system will provide comprehensive borehole integrity monitoring with improved resolution over existing techniques. As an application of this imaging system, we will characterize the effectiveness of next-generation wellbore completion technology (NETL, SNL), and will demonstrate that, unlike current technology, the proposed approach can successfully characterize foamed cements.

15 GEOTHERMAL ENERGY↗

High-Resolution 3D Acoustic Borehole Integrity Monitoring (Final Report)

Real-time, in-situ, high spatial resolution (sub-cm) imaging of the near-borehole environment would revolutionize wellbore diagnostics and integrity assessment by direct observation of defects. It is becoming increasingly apparent, that better understanding of the near-wellbore environment is required to meet the safety and operational needs in challenging environments such as those present in subsurface energy extraction (geothermal) and storage (CO2 sequestration) applications. Therefore, it is important to have a more robust ability to image the near-borehole and reliably detect defects. It was proposed to further develop and improve our advanced 3D imaging system to evaluate casing defects (e.g. corrosion) and cement quality in either open- or cased-borehole with the ultimate goal to develop a commercially deployable technology. The system consists of a unique acoustic source (LANL) and advanced inversion techniques for image processing (LANL, ORNL). This system will provide comprehensive borehole integrity monitoring with improved resolution over existing techniques. As an application of this imaging system, we will characterize the effectiveness of next-generation wellbore completion technology (NETL, SNL), and will demonstrate that, unlike current technology, the proposed approach can successfully characterize foamed cements.

42 ENGINEERING↗

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗