Search NASA⌕ Search

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.

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
DOE OSTI2017

Two-Relaxation-Time Lattice Boltzmann Method and its Application to Advective-Diffusive-Reactive Transport

The lattice Boltzmann method (LBM) based on single-relaxation-time (SRT) or multiple-relaxation-time (MRT) collision operators is widely used in simulating flow and transport phenomena. The LBM based on two-relaxation-time (TRT) collision operators possesses strengths from the SRT and MRT LBMs, such as its simple implementation and good numerical stability, although tedious mathematical derivations and presentations of the TRT LBM hinder its application to a broad range of flow and transport phenomena. This paper describes the TRT LBM clearly and provides a pseudocode for easy implementation. Various transport phenomena were simulated using the TRT LBM to illustrate its applications in subsurface environments. These phenomena include advection-diffusion in uniform flow, Taylor dispersion in a pipe, solute transport in a packed column, reactive transport in uniform flow, and bacterial chemotaxis in porous media. Finally, the TRT LBM demonstrated good numerical performance in terms of accuracy and stability in predicting these transport phenomena. Therefore, the TRT LBM is a powerful tool to simulate various geophysical and biogeochemical processes in subsurface environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

California GAMA Special Study: Importance of River Water Recharge to Selected Groundwater Basins

River recharge represents 63%, 86% and 46% of modern groundwater in the Mojave Desert, Owens Valley, and San Joaquin Valley, respectively. In pre-modern groundwater, river recharge represents a lower fraction: 36%, 46%, and 24% respectively. The importance of river water recharge in the San Joaquin valley has nearly doubled and is likely the result of a total increase of recharge of 40%, caused by river water irrigation return flows. This emphasizes the importance of recharge of river water via irrigation for renewal of groundwater resources. Mountain front recharge and local precipitation contribute to recharge of desert groundwater basins in part as the result of geological features focusing scarce precipitation promoting infiltration. River water recharges groundwater systems under lower temperatures and with larger water table fluctuations than local precipitation recharge. Surface storage is limited in time and volume, as evidenced by cold river recharge temperatures resulting from fast recharge, compared to the large capacity for subsurface storage. Groundwater banking of seasonal surface water flows therefore appears to be a natural and promising method for increasing the resilience of water supply systems. The distinct isotopic and noble gas signatures of river water recharge, compared to local precipitation recharge, reflecting the source and mechanism of recharge, are valuable constraints for numerical flow models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework

Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Facilitated Strontium Transport by Remobilization of Strontium-Containing Secondary Precipitates in Hanford Site Subsurface

Significantly enhanced immobilization of radionuclides (such as 90Sr and 137Cs) due to adsorption and coprecipitation with neo-formed colloid-sized secondary precipitates has been reported at the U.S. Department of Energy’s Hanford Site. However, the stability of these secondary precipitates containing radionuclides in the subsurface under changeable field conditions is not clear. Here, the authors tested the remobilization possibility of Sr containing secondary precipitates (nitrate-cancrinite) in the subsurface using saturated column experiments under different geochemical and flow conditions. The columns were packed with quartz sand that contained secondary precipitates (nitrate-cancrinite containing Sr), and leached using colloid-free solutions under different flow rates, varying pH, and ionic strength conditions. The results indicate remobilization of the neo-formed secondary precipitates could be possible given a change of background conditions. The remobility of the neo formed precipitates increased with the rise in the leaching solution flow rate and pH (in a range of pH 4 to 11), as well as with decreasing solution ionic strength. The increased mobility of Sr-containing secondary precipitates with changing background conditions can be a potential source for additional radionuclide transport in Hanford Site subsurface environments.

colloid↗
DOE OSTI2018

FEHM V3.0 - V3.X

The numerical background of the FEHM computer code can be traced to the early 1970s when it was used to simulate geothermal and hot dry rock reservoirs. The primary use over a number of years was to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the potential Yucca Mountain repository. Today FEHM is used to simulate groundwater and contaminant flow and transport in deep and shallow, fractured and un-fractured porous media throughout the US DOE complex. FEHM has proved to be a valuable asset on a variety of projects of national interest including Environmental Remediation of the Nevada Test Site, the LANL Groundwater Protection Program, geologic CO2 sequestration, Enhanced Geothermal Energy (EGS) programs, Oil and Gas production, Nuclear Waste Isolation, and Arctic Permafrost. Subsurface physics has ranged from single fluid/single phase fluid flow when simulating basin scale groundwater aquifers to complex multifluid/ multi-phase fluid flow that includes phase change with boiling and condensing in applications such as unsaturated zone surrounding nuclear waste storage facility or leakage of CO2/brine through faults or wellbores. The numerical method used in FEHM is the control volume method (CV) for fluid flow and heat transfer equations which allows FEHM to exactly enforce energy/mass conservation; while an option is available to use the finite element (FE) method for displacement equations to obtain more accurate stress calculations. In addition to these standard methods, an option to use FE for flow is available, as well as a simple Finite Difference scheme.

Pawar, Rajesh↗
DOE OSTI2016

Seasonal hyporheic dynamics control coupled microbiology and geochemistry in Colorado River sediments

Riverbed microbial communities play an oversized role in many watershed ecosystem functions, including the processing of organic carbon, cycling of nitrogen, and alterations to metal mobility. The structure and activity of microbial assemblages depend in part on geochemical conditions set by river-groundwater exchange or hyporheic exchange. In order to assess how seasonal changes in river-groundwater mixing affect these populations in a snowmelt-dominated fluvial system, vertical sediment and pore water profiles were sampled at three time points at one location in the hyporheic zone of the Colorado River and analyzed by using geochemical measurements, 16S rRNA gene sequencing, and ecological modeling. Oxic river water penetrated deepest into the subsurface during peak river discharge, while under base flow conditions, anoxic groundwater dominated shallower depths. Over a 70 cm thick interval, riverbed sediments were therefore exposed to seasonally fluctuating redox conditions and hosted microbial populations statistically different from those at both shallower and deeper locations. Additionally, microbial populations within this zone were shown to be the most dynamic across sampling time points, underlining the critical role that hyporheic mixing plays in constraining microbial abundances. Given such mixing effects, we anticipate that future changes in river discharge in mountainous, semiarid western U.S. watersheds may affect microbial community structure and function in riverbed environments, with potential implications for biogeochemical processes in riparian regions.

16S rRNA↗
DOE OSTI2020

Strong temporal variability in methane fluxes from natural gas well pad soils

We measured methane and carbon dioxide fluxes at natural gas well pad soils and undisturbed soils in the Rocky Mountain and Gulf Coast regions of the United States, including producing and gas storage wells. We collected both short-term (15 min) and multi-day (between 3 and 8), continuous measurements at 47 well pads and two undisturbed locations. Methane fluxes varied by more than an order of magnitude over periods as short as 30 min (e.g., 19–593 mg m -2 h -1 in one instance), and diurnal and seasonal variability was also significant (e.g., spring-to-fall change from 509 to 14174 mg m -2 h -1 ). We hypothesize that short-term flux variability was caused by pulsed flow of methane during its migration through the subsurface. Barometric pressure and well conditions likely impacted fluxes, but we found only weak evidence for this. Bacterial methanotrophy appeared to impact methane flux magnitude and variability. We injected methane into the subsurface at one well, and we found that, while fluxes of methane and carbon dioxide, and combustible soil gas concentrations, increased in response to the injection, the response was not uniform, and fluxes exhibited high hourly-scale variability, in spite of a constant injection rate. Methane fluxes tended to be higher at well pad soils compared to background soils (often much higher), and fluxes tended to be higher at well pad locations closer to the well head.

03 NATURAL GAS↗
DOE OSTI2020

Operational and geological controls of coupled poroelastic stressing and pore-pressure accumulation along faults: Induced earthquakes in Pohang, South Korea

Coupled poroelastic stressing and pore-pressure accumulation along pre-existing faults in deep basement contribute to recent occurrence of seismic events at subsurface energy exploration sites. Our coupled fluid-flow and geomechanical model describes the physical processes inducing seismicity corresponding to the sequential stimulation operations in Pohang, South Korea. Simulation results show that prolonged accumulation of poroelastic energy and pore pressure along a fault can nucleate seismic events larger than M w 3 even after terminating well operations. In particular the possibility of large seismic events can be increased by multiple-well operations with alternate injection and extraction that can enhance the degree of pore-pressure diffusion and subsequent stress transfer through a rigid and low-permeability rock to the fault. This study demonstrates that the proper mechanistic model and optimal well operations need to be accounted for to mitigate unexpected seismic hazards in the presence of the site-specific uncertainty such as hidden/undetected faults and stress regime.

58 GEOSCIENCES↗
DOE OSTI2019

Four Draft Single-Cell Genome Sequences of Novel, Nearly Identical Kiritimatiellaeota Strains Isolated from the Continental Deep Subsurface

ABSTRACT The recently proposed bacterial phylum Kiritimatiellaeota represents a globally distributed monophyletic clade distinct from other members of the Planctomycetes , Verrucomicrobia , and Chlamydiae (PVC) superphylum. Here, we present four phylogenetically distinct single-cell genome sequences from within the Kiritimatiellaeota lineage sampled from deep continental subsurface aquifer fluids of the Death Valley Regional Flow System in the United States.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Local Sensitivities of the Gulf Stream Separationa

Abstract Robust and accurate Gulf Stream separation remains an unsolved problem in general circulation modeling whose resolution will positively impact the ocean and climate modeling communities. Oceanographic literature does not face a shortage of plausible hypotheses that attempt to explain the dynamics of the Gulf Stream separation, yet a single theory that the community agrees on is missing. In this paper, the authors investigate the impact of the deep western boundary current (DWBC), coastline curvature, and continental shelf steepening on the Gulf Stream separation within regional configurations of the Massachusetts Institute of Technology General Circulation Model. Artificial modifications to the regional bathymetry are introduced to investigate the sensitivity of the separation to each of these factors. Metrics for subsurface separation detection confirm the direct link between flow separation and the surface expression of the Gulf Stream in the Mid-Atlantic Bight. It is shown that the Gulf Stream separation and mean surface position are most sensitive to the continental slope steepening, consistent with a theory proposed by Melvin Stern in 1998. In contrast, the Gulf Stream separation exhibits minimal sensitivity to the presence of the DWBC and coastline curvature. The implications of these results to the development of a “separation recipe” for ocean modeling are discussed. This study concludes adequate topographic resolution is a necessary, but not sufficient, condition for proper Gulf Stream separation.

58 GEOSCIENCES↗
DOE OSTI2024

Model data for Flood Frequency Analysis using Stochastic Storm Transposition and an Integrated Surface-Subsurface Hydrological Model

This archived provides scripts and input files used for the implementation of a novel approach to conduct process-based Flood Frequency Analysis using a Stochastic Storm Transposition (SST) and an Integrated Surface-Subsurface Hydrological Model (ISSHM). As a proof-of-concept, this study uses the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct flood frequency analysis by simulating the flood response to 5,000 annual synthetic storm events in a ~2000 km2 Southeast Texas watershed.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include TXT, CSV, DAT, SBATCH, SHP, TIF, NetCDF, and HDF5 files, which can be read through Python scripts. The input files for the ATS model and RainyDay model have .XML and .SST extensions, respectively, and can be edited in any commonly used text editors.This archive contains:* Scripts and data files essential for generating the ATS model input. It uses the Watershed Workflow package to produce both mesh and ATS input files. * Jupyter notebooks designated for the ATS model evaluation, covering both long-term simulations and 40 rainfall-runoff events.* Input files required to simulate SST storm events using RainyDay.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Radionuclide Migration through Sediment and Concrete: 16 Years of Investigations

The Waste Management Project provides safe, compliant, and cost-effective waste management services for the Hanford Site and the U.S. Department of Energy (DOE) complex. Part of these services includes safe disposal of low-level waste and mixed low-level waste at the Hanford Low-Level Waste Burial Grounds in accordance with the requirements of DOE Order 435.1, Radioactive Waste Management. To partially satisfy these requirements, performance assessment analyses were completed and approved. DOE Order 435.1 also requires continuing data collection to increase confidence in the critical assumptions used in these analyses to characterize the operational features of the disposal facility that are relied on to satisfy the performance objectives identified in the order. Cement-based solidification and stabilization is considered for hazardous waste disposal because it is easily done and cost-efficient. One critical assumption is that concrete will be used as a waste form or container material at the Hanford Site to control and minimize the release of radionuclide constituents in waste into the surrounding environment. Concrete encasement would contain and isolate the waste packages from the hydrologic environment and act as an intrusion barrier. Any failure of concrete encasement may result in water intrusion and consequent mobilization of radionuclides from the waste packages. The radionuclides iodine-129, selenium-75, technetium-99, and uranium-238 have been identified as long-term dose contributors (Mann et al. 2001; Wood et al. 1995). Because of their anionic nature in aqueous solutions, these constituents of potential concern may be released from the encased concrete by mass flow and/or diffusion and migrate into the surrounding subsurface environment (Serne et al. 1989; 1992; 1993a, b; 1995). Therefore, it is necessary to assess the performance of the concrete encasement structure and the ability of the surrounding soil to retard radionuclide migration. Each of the test methods performed throughout the lifetime of the project has focused on different aspects of the concrete waste form weathering process. Diffusion of different analytes [technetium-99 (Tc-99), iodine-125 (I-125), stable iodine (I), uranium (U), and rhenium (Re)] has been quantified from experiments under both saturated and unsaturated conditions. The water-saturated conditions provide a conservative estimate of the concrete’s performance in situ, and the unsaturated conditions provide a more accurate estimate of the diffusion of contaminants from the concrete.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Electric Power Generation from Low to Intermediate Temperature Resources

This project was designed to test the concept on the Eland-Lodgepole Field near Dickinson, North Dakota in the Williston Basin. The field is in secondary-recovery water-flood and consists of 12 producing oil wells, 5 water injection wells and one disposal well. Water production at the site averages approximately 320 gallons per minute (20.2 l s-1) and the temperature is 100 ⁰C. Engineers at Ormat estimated power production potential with the existing resource to be approximately 350 kWh. Unfortunately, ownership of the field was transferred from Encore, Inc., to Denbury, Inc., within the first week of the project. After two years of discussion and planning, Denbury decided not to pursue this project due to complications with the site location and its proximity to Patterson Lake. Attempts to find other partners operating in the Williston Basin were unsuccessful. Consequently, we were unable to pursue the primary objective of the project. However, during negations with Denbury and subsequent time spent contacting other potential partners, we focused on objectives 2 and 3 and developed a clear understanding of the potential for co-produced production in the Williston Basin and the best practices for developing similar projects. At least nine water bearing formations with temperatures greater than 90 ⁰C extend over areas of several 10s of km2. The total energy contained in the rock volume of those geothermal aquifers is 283.6 EJ (1 EJ = 1018 J). The total energy contained in the water volume, determined from porosities which range from 2 percent to 8 percent, is 6.8 EJ. The aquifers grouped by 10 ⁰C temperature bins (Table 1) include one or more formations due to the bowl-shape structure of the basin. Table 1. Summary of energy available in geothermal aquifers in the Williston Basin Analysis of overall fluid production from active wells, units, fields and formations in North Dakota showed that few sites co-produce sufficient fluid for significant power production with ORC technology. Average co-produced water for 10,480 wells is 3.2 gallons per minute (gpm). Even excluding the tight formations, Bakken and Three Forks, average co-produced water for the remaining 3,337 is only 5 gpm. The output of the highest producing well is 184 gpm and the average of the top 100 wells is 52 gpm. Due to the depth of the oil producing formations in the Williston Basin, typically 3 km or greater, pumps are operated slowly to prevent watering out thus total fluid production is purposefully maintained at low volumes. There remain potential possibilities for development of geothermal fluids in the Williston Basin. Unitized fields in which water production from several tens of wells is collected at a single site are good possibilities for development. Water production in the unitized fields is greater than 1000 gpm is several areas. A similar possibility occurs where infill-drilling between Bakken and Three Forks horizontal wells has created areas where large volumes of geothermal fluids are available on multi-well pads and in unitized fields. Although the Bakken produces small amounts of water, the water/oil ration is typically less than 1, the oil and water mix produced at the well head can be sent through the heat exchanger on an ORC. It is estimated that several tens of MWh of power could be generated by a distributed system of ORC engines in the areas of high-density drilling in the Bakken Formation. Finally, horizontal drilling in water bearing formations is the other possibility. Several secondary recovery water-flood projects in the basin are producing water above 100 ⁰C at rates of 300 gpm to 850 gpm. Those systems also could produce several tens of MWh of power with ORC technology. Objective 3 of the project was highly successful. The program has produced 5 PhDs, 7 MS, and 3 BS students with theses in geothermal energy. The team has involved 7 faculty in 4 different engineering and science disciplines, ChE, EE, GE, and Geol. The team has produced 26 peer-reviewed papers and 62 presentations at professional meetings. Faculty involved in the program developed five graduate level courses covering different elements in heat flow and geothermal energy that are now offered in the Harold Hamm School of Geology and Geological Engineering. Lessons learned – Keys to developing a successful project;1. Determine target formations; a. Data from oil and gas operators, state oil and gas regulatory agencies, and state geological surveys help to identify producing formations and their properties; 2. Determine the quantity of energy available in the target formations; a. A complete thermal analysis of the basin or region yields the most useful information; b. Critical data include: BHT, heat flow, stratigraphy, lithology, lithological properties, and thermal conductivity, subsurface structure; 3. Determine fluid production potential; a. State oil and gas regulatory agencies, and state geological surveys have data on oil, gas and water production. State Water Commission/Agencies have data on water quality, aquifers, and regulations; b. Consider single horizontal wells, multiple conventional wells, and unitized fields; 4. Calculate energy production capacity of each formation based on different well combination and power plant scenarios. This is a broad overview rather than a site specific analysis; 5. Research and understand the local electrical power industry. Obtain the PPA before committing to the project; 6. Work with the high-level personnel in the oil company partner. Obtain an MOU that addresses all issues in the project including what to expect if the company goes out of business, is bought out, changes management, etc; and 7. Be prepared for project delays.

15 GEOTHERMAL ENERGY↗
DOE OSTI2015

Revealing the Role of Microbes in Controlling Contaminants

In Rifle, Colorado, Berkeley Lab earth scientist, Kenneth Hurst Williams, highlights the role subsurface microbial communities can play in controlling the flow of contaminants in groundwater. The DOE Joint Genome Institute is a key collaborator in the research. Williams is Component Lead of Watershed Structure and Controls within Berkeley Lab's Genomes-to-Watershed Scientific Focus Area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advancing process-based flood frequency analysis for assessing flood hazard and population flood exposure

Recent studies have showcased the use of process-based hydrological models with Stochastic Storm Transposition (SST) techniques to conduct Flood Frequency Analysis (FFA). This framework, referred hereby FFA-SST, has proved to be a robust strategy to estimate peak flows of specific annual exceedance probability (e.g., 100-year peak flow) that can reflect natural and anthropogenic disturbances, including changes in land use and meteorological patterns. With the objective of advancing the FFA-SST framework, this study presents for the first time the use of an Integrated Surface-Subsurface Hydrological Model (ISSHM) to conduct FFA-SST by extending the analysis from peak flow responses to flood extent, enabling a unique view and analysis of flood hazard and population flood exposure at the basin scale. As a proof-of-concept, we used the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct FFA-SST by simulating the flood response to 5,000 annual synthetic storm events in a 2,227 $km^2$ Southeast Texas watershed. We demonstrate that ATS, without site-specific calibration, provides a robust process-based representation of peak flows, flood extent, streamflow, evapotranspiration, soil moisture content, and water storage changes. Our results and analyses, covering frequency curves up to a 500-year return period for peak flows, basin inundation fractions, and the number of people exposed to flooding, offer a unique perspective to analyze flood impacts across spatial scales. Overall, this study provides critical insights for flood risk management by extending the FFA-SST framework to include both flood hazard and population flood exposure analyses at the basin scale. Such an approach will empower stakeholders and disaster emergency agencies with a more comprehensive understanding of flood impacts across the entire basin domain, facilitating informed decision-making for flood risk assessment and management.

58 GEOSCIENCES↗
DOE OSTI2019

EGS Collab Project: Accomplishments and Plan

To better understand stimulation and inter-well flow in crystalline rock with the goal of developing and improving subsurface heat exchangers for enhanced geothermal energy recovery, the EGS (Enhanced Geothermal Systems) Collab project is developing highly monitored, well- characterized intermediate scale (~10-20 m) field test beds, stimulating, and performing flow tests. These tests provide the ability to evaluate concepts, data to test and validate models, and tools to be used for EGS development, which may be evaluated at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) project. Numerical simulation is key to the project, for experiment design and interpretation of results. Testbed characterization included optical and acoustic televiewers, full waveform seismic, electromagnetic, gamma, and temperature. The test block was further characterized using seismic tomography (p- and s-), electrical resistance tomography (ERT), and hydrologic characterization including tracer tests. Coupled process models are constrained by the detailed site characterization data. Installed monitoring systems include: 1) passive seismic/acoustic emissions; 2) continuous active source seismic; 3) ERT, with dynamic imaging using high contrast fluids; and 4) distributed acoustic, temperature, and strain fiber optic sensors. During flow and stimulation tests, fracture aperture strain monitoring is performed using the Step-rate Injection Method for Fracture In-situ Properties (SIMFIP) tool. Pressure and flow in the injection and production boreholes are monitored, and tracer tests and wavefield imaging and inversion are performed to provide additional constraints. Stimulations have been performed at several locations, connecting our injection and production wells. An overview of the EGS Collab project is presented to provide a framework for other EGS Collab presentations, along with results from numerical simulations and experiments performed at the Sanford Underground Research Facility in South Dakota.

Enhanced Geothermal Systems, EGS Collab, stimulati↗