SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Three-dimensional geologic mapping to assess geothermal potential: examples from Nevada and Oregon

Geologic structure plays an important role in controlling fluid flow in geothermal systems. In particular, very complex structural settings, consisting of many closely spaced and intersecting faults, host many geothermal systems. To elucidate the key geologic factors that affect fault-controlled geothermal circulation, it is critical to precisely characterize the structural and stratigraphic geometries in these complex settings. Here, we present a methodology and the results of 3D geologic analyses of two geothermal systems in the Basin and Range, USA. This methodology is a quantitative and geologically focused technique that can be used to precisely characterize geothermal areas, in a time when future geothermal growth demands increased exploration precision and efficiency. Surficial and subsurface geologic and geophysical data are synthesized in the construction of detailed 3D geologic maps of geothermal areas. Based on these 3D geologic maps, we examine several geologic attributes that control permeability development and geothermal fluid flow along faults. We use the stress state of faults and the distribution of structural discontinuities (i.e., fault intersections and fault terminations) to identify locations of upflow along faults in these geothermal systems. These results and the methodology presented herein are directly applicable to structurally controlled geothermal fields in the Basin and Range and worldwide. As development focus shifts toward blind geothermal resources, integration of precisely characterized subsurface structural information into exploration methods will be increasingly critical to continued growth in geothermal exploration and development.

15 GEOTHERMAL ENERGY↗

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 OSTI2019

Exploratory Project: A vegetative facies-based multiscale approach to modeling nutrient transport in the Columbia river Basin

The exploratory project titled “A vegetative facies-based multiscale approach to modeling nutrient transport in the Columbia river Basin” was a collaborative effort between Stanford University (PI: Ilenia Battiato) and University of Southern California (co-PI: Felipe de Barros). The primary objective of this project was to develop a predictive understanding of the impact that submerged vegetation has on flow and transport in riverine systems. This is based on the recognition that vegetation-mediated transport and its temporal and spatial coupling with hydrological and ecological responses is a critical component in achieving a predictive understanding of the future impact of, and the risks associated to, anthropogenic and climate change stressors on nutrient cycling in the Subsurface Interaction Zone (SIZ). Although the full project framework was focused on developing top-down and bottom-up coupling between processes at different scales as outlined in Figure 1, its exploratory counterpart was organized in two major thrusts: the first related to quantifying flow resistance in vegetated rivers (Level II) and the second connected to quantifying and predicting the coupling between vegetation growth and hydrodynamics in riverine systems (Level III). In this regard, we successfully used porous media theory to model flow and transport in vegetated channels at the canopy scale (Level II) and have demonstrated that the aforementioned models have predictive capabilities to model both mean flow and transport at the canopy scale. At Level III, we have developed a computational framework in OpenFOAM, CladoFOAM, that is able to model and predict submerged vegetation growth at the km-scale and to study the impact that climate stressors can have on algae blooms from basic conservation laws at the appropriate scale. In Section 2, we briefly outline the project objectives; in Section 3, we present our scientific findings of all manuscripts already published and in preparation. We conclude in Section 4 with comments about future lines of investigation that may be pursued in the context of a larger-scope project.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Geochemical Exports to River From the Intrameander Hyporheic Zone Under Transient Hydrologic Conditions: East River Mountainous Watershed, Colorado

To understand how redox processes influence carbon, nitrogen, and iron cycling within the intrameander hyporheic zone, we developed a biotic and abiotic reaction network and incorporated it into the reactive transport simulator PFLOTRAN. Additionally, two-dimensional reactive flow and transport simulations were performed (1) to evaluate how transient hydrological conditions control the lateral redox zonation within an intrameander region of the East River in Colorado and (2) to quantify the impact of a single meander on subsurface exports of carbon and other geochemical species to the river. The meander's overall contribution to the river was quantified by integrating geochemical outfluxes along the outside of the meander bend. The model was able to capture the field-observed trends of dissolved oxygen, nitrate, iron, pH, and total inorganic carbon along a 2-D transect. Consistent with field observations, simulated dissolved oxygen and nitrate decreased along the intrameander flow paths while iron (Fe 2+ ) concentration increased. The simulation results further demonstrated that the reductive potential of the lateral redox zonation was controlled by groundwater velocities resulting from river stage fluctuations, with low-water conditions promoting reducing conditions. The sensitivity analysis results showed that permeability had a more significant impact on biogeochemical zonation compared to the reaction pathways under transient hydrologic conditions. Lastly, the simulation results further indicated that the meander acted as a sink for organic and inorganic carbon as well as iron during the extended baseflow and high-water conditions; however, these geochemical species were released into the river during the falling limb of the hydrograph.

58 GEOSCIENCES↗

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 OSTI2018

Integrating geophysical monitoring data into multiphase fluid flow reservoir simulation

Simulation of multiphase flow systems are of critical importance in managing hydrological systems. Flow simulations are affected by a number of factors including structure and flow properties including porosity and permeability as well as the anisotropy and heterogeneity of these properties. In many cases traditional hydrological and reservoir data are highly affected by these parameters, but are not directly sensitive to them. As such modellers often adjust these parameters in an ad hoc manner until solutions numerically converge. Simulation models are generally based on structural data from reflection seismics whose physical flow properties are then populated using geostatistical extrapolation techniques utilizing a sparse number of borehole logs and core analysis. In multiphase systems including enhanced oil recovery and carbon capture and sequestration uncertainties regarding phase-dependent physical properties confounds this challenge further. Geophysical methods provide a means by which to gain an improved understanding of phase distributions in the subsurface. In this paper we will look at applications from active carbon capture and sequestration and enhanced oil recovery applications, as well as synthetic examples. Geophysical data including electromagnetic and gravity are inverted using structural constraints from the reservoir model. Inversions are then mapped into flow properties using calibrated relations such as Archie’s Equation. The coupled models can then be used to both verify and improve on the reservoir flow model which improves it’s predictive power and utility as a management tool.

02 PETROLEUM↗

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 OSTI2018

Casing-integrity mapping using top-casing electrodes and surface based EM fields

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage and hazardous waste disposal. We introduce a simple non-invasive technology for casing integrity screening, based on the continuity of electrical current flow. Applying low frequency current to a wellhead, with a distant return electrode, produces a casing current dependent on the properties and depth extent of the well casing as well as the background formation. These currents in-turn generate surface electrical fields that can be captured in a radial profile and be used to analyze properties of the well casing. Numerical modeling results reveal a strong relation of the electric field to the casing properties and depth extent of the well. A small breakage in the casing produces a profile coincident to a cased well with a completion depth above the break. A corroded patch, where the casing conductivity is reduced, also alters the field profiles and its depth may be estimated by comparing to the profile expected from the well completion diagrams. The electric field profiles are also strongly dependent on background resistivity distributions and on whether the well was drilled using water or oil-based drilling fluids. We validate the proof of concept in a field experiment, where we applied currents at the wellheads of two wells with different casing lengths. The two profiles were similar in appearance but offset in amplitude by more than a factor of 5, consistent with the theoretical analysis as well as the 3D modeling results. These results demonstrate that our proposed approach has promise for mapping the general casing condition without well intervention. This approach can be a practical and effective tool for rapidly screening a number of wells before expensive logging-based technologies are employed for casing inspection in detail.

58 GEOSCIENCES↗

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↗

NASA NTRS1984

Interpretation of the lunar microwave brightness temperature spectrum - Feasibility of orbital heat flow mapping

The feasibility of an orbital mapping of heat flow by means of microwave radiometers is examined on the basis of a detailed model of the lunar regolith. The results obtained are discussed in terms of contributions of physical properties and regolith structure to the spectral signature of lunation-averaged brightness temperatures. Even for regions for which remote heat flow detection is not feasible, detailed interpretations of remote measurements in terms of physical properties and regolith structure could be made by means of the model presented. Surface and subsurface scattering effects are considered in evaluating possible wavelength dependencies of the microwave emissivity.

Keihm, S. J.↗