Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface flow”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Clear Lake Volcanic Field Disasters: Creating a Deformation Record Using InSAR to Assess Hazards and Detect Volcanic Unrest in Clear Lake Volcanic Field

Clear Lake Volcanic Field (CLVF) in northern California is at a high threat potential for volcanic hazards. Eruptions leading to increased seismic activity could result in silicic domes, cinder cones, and flows that would be dangerous to the residential areas near the volcanic field. Remotely sensed Earth observations can reveal volcanic processes in the subsurface, which are essential to the timely monitoring of potential volcanic activity. In particular, Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Digital Elevation Model (DEM) data capture relative surface deformation at unprecedented high spatial and temporal resolutions. Leveraging C-SAR and DEMs, we conducted interferometric analysis from January 2016 to December 2023. Our results demonstrate 1) the mean surface displacement velocity of the Clear Lake Volcanic Field is measured to undergo 5 to 10-centimeter scale deformation and shows a strong relationship with the surrounding faults, 2) apparent seasonal differences in rates of surface change, and 3) seismic activity associated with the geyser geothermal field has a strong association with cumulative surface displacement, with active fault zones having 2 to 5 cm of additional displacement. Results indicate that deformation is linked to deep pressure sources causing stresses on the surficial environment that should be considered in hazard mitigation. This study provides a baseline of historic deformation, aiding hazard analysts in communication efforts and streamlining decision-making for potential risks to region residents.

remote sensing↗
NASA NTRSDate not supplied

Creating a Deformation Time Series Utilizing InSAR to Assess Hazards and Detect Volcanic Unrest in Clear Lake Volcanic Field, California

According to the U.S Geological Survey, the Clear Lake Volcanic Field (CLVF) in northern California is at a high threat potential for volcanic hazards. Eruptions leading to increased seismic activity could result in silicic domes, cinder cones, and flows that would be dangerous to the residential areas near the volcanic field. Remotely sensed Earth observations can reveal volcanic processes in the subsurface, which are essential to the timely monitoring of potential volcanic activity. In particular, Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Digital Elevation Model (DEM) data capture relative surface deformation at unprecedented high spatial and temporal resolutions. Leveraging C-SAR and DEMs, we conducted interferometric analysis from January 2016 to December 2023. Our results demonstrate 1) the mean surface displacement velocity of the Clear Lake Volcanic field is measured to undergo 5 to 10-centimeter scale deformation and shows a strong correlation with the surrounding faults, 2) apparent seasonal differences in rates of surface change, and 3) seismic activity associated with the geyser geothermal field has a strong correlation with cumulative surface displacement, with active fault zones having 2 to 5 cm of additional displacement. Results indicate that deformation is linked to deep pressure sources causing stresses on the surficial environment that should be considered in hazard mitigation. This study provides a baseline of historic deformation, aiding hazard analysts in communication efforts and streamlining decision-making for potential risks to region residents.

Ivan Tochimani-Hernandez↗
DOE OSTI2024

Modeling the Cost of Onshore CO2 Pipeline Transport and Onshore CO2 Saline Storage

This paper describes the FECM/NETL CO2 Transport Cost Model (CO2_T_COM), a technoeconomic model of CO2 transport by pipeline, and the FECM/NETL CO2 Saline Storage Cost Model (CO2_S_COM), a technoeconomic model of storage of CO2 in a deep, subsurface saline formation. The results of applying CO2_T_COM to calculate break-even CO2 prices to transport CO2 at different mass flow rates and different distances are presented. Similarly, CO2_S_COM is used to calculate the break-even CO2 price for storing CO2 in 314 potential storage formations across the US. These break-even prices are used to construct cost-supply curves for CO2 storage which are presented on a national and regional basis. CO2_T_COM and CO2_S_COM are the most comprehensive open source technoeconomic models available for analyzing CO2 pipeline transport and CO2 saline storage. To be presented at the SPE?AAPG/SEG Carbon Capture Utilization and Storage Conference in Houston, TX, March 11-13, 2024.

Morgan, David↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗
DOE OSTI2021

2019 Groundwater Monitoring Report Project Shoal Area: Subsurface Corrective Action Unit 447

The Project Shoal Area is a site in Nevada where an underground nuclear test was conducted in 1963. It later came to be known as the Shoal, Nevada, Site. Surface contamination at the site has been remediated, but investigation of groundwater contamination resulting from the test is still in the corrective action process. Annual sampling and water-level monitoring are conducted as part of the subsurface corrective action strategy, which has focused on revising the site conceptual model and evaluating the adequacy of the monitoring well network. It has also included enhancements to the monitoring well network to address uncertainties in the groundwater flow direction and the cause of rising water levels in site wells west of the shear zone since the first hydrologic characterization (HC) wells were installed in 1996. Revisions to the site conceptual model and enhancements to the monitoring strategy were provided to Nevada Division of Environmental Protection (NDEP) in the Addendum to: Corrective Action Decision Document/Corrective Action Plan (CADD/CAP) for the Subsurface Corrective Action Unit 447 Shoal, Nevada, Site. NDEP approved the addendum to the CADD/CAP, which included an expanded contaminant boundary and compliance boundary for the site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2022

Computational Analysis of Coupled Geoscience Processes in Fractured and Deformable Media

Prediction of flow, transport, and deformation in fractured and porous media is critical to improving our scientific understanding of coupled thermal-hydrological-mechanical processes related to subsurface energy storage and recovery, nonproliferation, and nuclear waste storage. Especially, earth rock response to changes in pressure and stress has remained a critically challenging task. In this work, we advance computational capabilities for coupled processes in fractured and porous media using Sandia Sierra Multiphysics software through verification and validation problems such as poro-elasticity, elasto-plasticity and thermo-poroelasticity. We apply Sierra software for geologic carbon storage, fluid injection/extraction, and enhanced geothermal systems. We also significantly improve machine learning approaches through latent space and self-supervised learning. Additionally, we develop new experimental technique for evaluating dynamics of compacted soils at an intermediate scale. Overall, this project will enable us to systematically measure and control the earth system response to changes in stress and pressure due to subsurface energy activities.

58 GEOSCIENCES↗
NASA NTRS1992

Regional variations in the stability and diffusion of water-ice in the Martian regolith

Geologic evidence suggests subsurface water-ice has played an important role in the formation of Martian landforms. Forms of mass-wasting such as debris aprons and flow patterns on valley floors suggest creep deformation of ice-laden soil, while thermokarst and chaotic terrain suggest once extensive deposits of ground ice that were later removed. The global distribution of ice-related morphology was mapped. The mapping showed regional variation, in both latitude and longitude, in the distribution of debris aprons, concentric fill craters, and 'softened' crater profiles.

Mellon, Michael T.↗
NASA NTRS2018

Physics-Based Approach to Predict the Solar Activity Cycles

Observations of the complex highly non-linear dynamics of global turbulent flows and magnetic fields are currently available only from Earth-side observations. Recent progress in helioseismology has provided us some additional information about the subsurface dynamics, but its relation to the magnetic field evolution is not yet understood. These limitations cause uncertainties that are difficult take into account, and perform proper calibration of dynamo models. The current dynamo models have also uncertainties due to the complicated turbulent physics of magnetic field generation, transport and dissipation. Because of the uncertainties in both observations and theory, the data assimilation approach is natural way for the solar cycle prediction and estimating uncertainties of this prediction. I will discuss the prediction results for the upcoming Solar Cycle 25 and their uncertainties and affect of Ensemble Kalman Filter parameters to resulting predictions.

Physics-Based↗
DOE OSTI2025

Monitoring river flow status using low-cost wildlife camera and image segmentation artificial intelligence

Continuous measurement and monitoring of surface water coverage in non-perennial streams are essential for understanding the exchange fluxes between surface and subsurface waters under both inundated and non-inundated conditions. In this study, a wildlife camera photo-based framework was developed to monitor small stream water inundation, depth, discharge, and velocity. Two advanced machine learning models, YOLOv8 and Mask2Former, were utilized to efficiently analyze images captured by wildlife cameras. The accuracy of the framework was validated against on-site depth measurements at six sites in the Yakima River Basin, along with the gage height, discharge, and velocity data from four USGS sites. This approach facilitates long-term, continuous monitoring and quantification of river intermittency and water availability with high precision and low cost, thereby advancing river ecosystem research and management.

machine learning↗
DOE OSTI2022

Transport of polymer-coated metal–organic framework nanoparticles in porous media

Abstract Injecting fluids into deep underground geologic structures is a critical component to development of long-term strategies for managing greenhouse gas emissions and facilitating energy extraction operations. Recently, we reported that metal–organic frameworks are low-frequency, absorptive-acoustic metamaterial that may be injected into the subsurface to enhance geophysical monitoring tools used to track fluids and map complex structures. A key requirement for this nanotechnology deployment is transportability through porous geologic media without being retained by mineral-fluid interfaces. We used flow-through column studies to estimate transport and retention properties of five different polymer-coated MIL-101(Cr) nanoparticles (NP) in siliceous porous media. When negatively charged polystyrene sulfonate coated nanoparticles (NP-PSS-70K) were transported in 1 M NaCl, only about 8.4% of nanoparticles were retained in the column. Nanoparticles coated with polyethylenimine (NP-PD1) exhibited significant retention (> 50%), emphasizing the importance of complex nanoparticle-fluid-rock interactions for successful use of nanofluid technologies in the subsurface. Nanoparticle transport experiments revealed that nanoparticle surface characteristics play a critical role in nanoparticle colloidal stability and as well the transport.

58 GEOSCIENCES↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide

The project Characterization of the Triassic Newark Basin of New York and New Jersey for Geologic Storage of Carbon Dioxide is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. Sandia Technologies, LLC, and co-investigator Conrad Geoscience Corporation, examined the potential for large-scale, permanent CO2 storage in sedimentary strata within the Newark Rift Basin. The Newark Rift Basin underlies an industrialized, developed region comprising parts of New York, New Jersey, and Pennsylvania. The project characterized and investigated the suitability of Triassic age sedimentary formations for potential geologic CO2 storage. The project team drilled and cored two test wells to define the sedimentary geologic formations underlying the basin and to document or reach basement rock. With this geologic characterization phase, an integration of seismic, geologic, borehole, and formation core results provided a higher resolution assessment of CO2 storage potential. The Stockton Formation is known to be a potentially favorable geologic storage formation in the basin. In 2011, the 1-NYSTA Tandem Lot stratigraphic test well was drilled to a depth of 6,855 feet in the northern portion of the Newark Basin in southern New York State. Approximately 9 miles south-southeast on the Lamont Doherty Campus, TW-4 was drilled and cored in 2013 to a depth of 1,802 feet and contacted apparent igneous basement at a depth of 1,712 feet. Both wells penetrated the Palisades Sill ranging from 800 feet thick in the eastern well to approximately 1,800 feet in thickness at the 1-NYSTA Tandem Lot deep drill site. A diabase sill can provide an excellent seal and dense confining layer for potential CO2 storage reservoirs and flow layers that are situated beneath it within the Stockton Sandstone. The Stockton Sandstone was encountered beneath the sill in the TW-4 well on the Lamont campus, and data integration suggests that it was likely observed near total depth in the deep 1-NYSTA Tandem Lot well. The test wells confirm and define reservoirs are present beneath the sill and offer CO2 storage potential. The integration of geologic and reservoir characterization of well logs, formation cores, and formation fluids indicated Triassic age lacustrine playa lake and mudbank shales of the Upper Passaic Group can provide an effective seal for the porous and permeable underlying sandstone reservoir layers. This project acquired seismic data, drilled borehole well logs, acquired core samples, and integrated these findings to provide a better understanding of the subsurface geologic formations in the Newark Rift Basin. These findings have contributed to a higher degree of accuracy in predicting potential geologic storage opportunities, while refining geologic storage capacity estimates for the indicated reservoirs and flow units.

.las↗
NASA NTRS2015

Restoration of the Apollo Heat Flow Experiments Metadata

Geothermal heat flow probes were deployed on the Apollo 15 and 17 missions as part of the Apollo Lunar Surface Experiments Package (ALSEP). At each landing site, the astronauts drilled 2 holes, 10-m apart, and installed a probe in each. The holes were 1- and 1.5-m deep at the Apollo 15 site and 2.5-m deep at the Apollo 17 sites. The probes monitored surface temperature and subsurface temperatures at different depths. At the Apollo 15 site, the monitoring continued from July 1971 to January 1977. At the Apollo 17 site, it did from December 1972 to September 1977. Based on the observations made through December 1974, Marcus Langseth, the principal investigator of the heat flow experiments (HFE), determined the thermal conductivity of the lunar regolith by mathematically modeling how the seasonal temperature fluctuation propagated down through the regolith. He also determined the temperature unaffected by diurnal and seasonal thermal waves of the regolith at different depths, which yielded the geothermal gradient. By multiplying the thermal gradient and the thermal conductivity, Langseth obtained the endogenic heat flow of the Moon as 21 mW/m(exp 2) at Site 15 and 16 mW/m(exp 2) at Site 17.

HFE metadata↗
DOE OSTI2023

Laboratory investigation of hydraulic fracturing in granitic rocks using active and passive seismic monitoring

SUMMARY Knowledge of the fracturing processes can be important for the optimization of pressurized fluid injection operations in the deep underground rock mass. Active and passive seismic monitoring techniques have been used in the field for tracking or mapping the propagating hydraulic fracture. Although both these monitoring techniques provide valuable information about the generated fracture network, it is difficult for either technique to comprehensibly identify the different processes associated with hydraulic fracturing. The combined active and passive monitoring has the potential for better characterization of the complex hydraulic fracturing phenomena. In this study, laboratory hydraulic fracturing experiments with combined active and passive seismic monitoring were conducted on true triaxially loaded Barre granite cubes with different fluid injection rates. The seismic inelastic fracturing was detected by 16 passive acoustic emission sensors, where 3678 and 2370 seismic source events were detected for the high and low injection rate experiments, respectively. For active monitoring, strong variations in the attributes of signals were observed which were transmitted through four source–receiver pairs, placed both perpendicular and parallel to the generated hydraulic fracture. Positive velocity changes were observed for active sensor pairs with ray paths passing through the generated hydraulic fracture indicating fluid permeation, whereas isolated dry deformation was characterized by a slight but permanent velocity decrease. Compared to velocity, the energy of the active signals was 1–2 orders of magnitude more sensitive to different hydraulic fracturing processes. However, the sensitivity and signatures of the active signal attributes were found to be dependent on the frequency range and direction of ray path with respect to the location of the generated fracture network. Using the coupled evaluation of the active and passive signals we were able to systematically identify various hydraulic fracturing processes including: (1) aseismic deformation, (2) fracture initiation and fluid permeation, (3) pressure build-up, (4) fracture propagation and (5) pressure release and leak-off. The results of this study showed that combining the respective advantages of active and passive seismic techniques and using both of them to monitor the failure processes can facilitate a more comprehensive understanding and better control of the hydraulic stimulations in subsurface operations.

Geochemistry & Geophysics↗
DOE OSTI2024

De-risking fault leakage risk and containment integrity for subsurface storage applications

The subsurface is pivotal in the energy transition, for the sequestration of CO 2 and energy storage. It is crucial to understand to what extent geological faults may form leakage pathways that threaten the containment integrity of these projects. Fault flow behavior has been studied in the context of hydrocarbon development, supported by observations from wells drilled through faults, but such observations are rare in geoenergy projects. Focusing on mechanical behavior as early indicator of potential leakage risks, a probabilistic Coulomb Failure Stress workflow is developed and demonstrated using data from the Decatur CO 2 sequestration project to rank faults based on their containment risk. The analysis emphasizes the importance of fault throw relative to reservoir thickness and pore pressure change in assessing reactivation risks. Integrating this mechanical assessment with geological and dynamic fault analyses contributes to derisking fault containment for geoenergy applications, providing valuable insights for the successful development of subsurface storage projects.

58 GEOSCIENCES↗
DOE OSTI2025

Subsurface hydrological controls on the short-term effects of hurricanes on nitrate–nitrogen runoff loading: a case study of Hurricane Ida using the Energy Exascale Earth System Model (E3SM) Land Model (v2.1)

When the nutrient level in the soil surpasses vegetation demand, nutrient losses due to surface runoff and subsurface leaching are the major reasons for the deterioration of water quality. The lower Mississippi River basin (LMRB) is one of the sub-basins that deliver the highest nitrogen loads to the Gulf of Mexico. Potential changes in episodic events induced by hurricanes may exacerbate water quality issue in the future. However, uncertainties in modeling the hydrologic response to hurricanes may limit the modeling of nutrient losses during such events. Using a machine learning approach, we calibrated the land component of the Energy Exascale Earth System Model (E3SM), or ELM, version 2.1, based on the water table depth (WTD) of a calibrated 3D subsurface hydrology model. While the overall performance of the calibrated ELM is satisfactory, some discrepancies in WTD remain in slope areas with low precipitation due to the missing lateral flow process in ELM. Simulations including biogeochemistry performed using ELM with and without model calibration showed important influences of soil hydrology, precipitation intensity, and runoff parameterization on the magnitude of nitrogen runoff loss and the leaching pathway. Despite such sensitivities, both ELM simulations produced reduced WTD and increased runoff and accelerated nitrate–nitrogen runoff loading during Hurricane Ida in August 2021, consistent with the observations. With observations suggesting more pronounced effects of Hurricane Ida on nitrogen runoff than the simulations, we identified factors for model improvement to provide a useful tool for studying hurricane-induced nutrient losses in the LMRB region.

54 ENVIRONMENTAL SCIENCES↗