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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
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 OSTI2012

An Experimental Study of Diffusivity of Technetium-99 in Hanford Vadose Zone Sediments

One of the methods being considered at the Hanford site in Washington for safely disposing of low-level radioactive wastes (LLW) is to encase the waste in concrete and entomb the packages in the Hanford vadose zone sediments. The current plan for waste isolation consists of stacking low-level waste packages on a trench floor, surrounding the stacks with reinforced steel, and encasing these packages with concrete. Any failure of the concrete encasement may result in water intrusion and consequent mobilization of radionuclides from the waste packages. The mobilized radionuclides may escape from the encased concrete by mass flow and/or diffusion and move into the surrounding subsurface sediments. It is therefore necessary to conduct an assessment of the performance of the concrete encasement structure and the surrounding soil’s ability to retard radionuclide migration. The retardation factors for radionuclides contained in the waste packages can be determined from measurements of diffusion coefficients for these contaminants through concrete and fill material. Because of their anionic nature in aqueous solutions, the radionuclides, 99Tc and 129I were identified as long-term dose contributors in LLW. The leachability and/or diffusion of these radionuclide species must be measured in order to assess the long-term performance of waste grouts when contacted with vadose-zone porewater or groundwater. To measure the diffusivity, a set of experiments were conducted using 99Tc-spiked concrete (with 0 and 4% metallic iron additions) in contact with unsaturated soil half-cells that reflected the typical moisture contents of Hanford vadose zone sediments. The 99Tc diffusion profiles in the soil half cells were measured after a time lapse of ~1.9 yr. Using the concentration profiles, the 99Tc diffusivity coefficients were calculated based on Fick’s Second Law.

techetium, diffusion, vadose zone sediments↗
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↗
DOE OSTI2016

Integrated surface/subsurface permafrost thermal hydrology: Model formulation and proof-of-concept simulations

The need to understand potential climate impacts and feedbacks in Arctic regions has prompted recent interest in modeling of permafrost dynamics in a warming climate. A new fine-scale integrated surface/subsurface thermal hydrology modeling capability is described and demonstrated in proof-of-concept simulations. The new modeling capability combines a surface energy balance model with recently developed three-dimensional subsurface thermal hydrology models and new models for nonisothermal surface water flows and snow distribution in the microtopography. Surface water flows are modeled using the diffusion wave equation extended to include energy transport and phase change of ponded water. Variation of snow depth in the microtopography, physically the result of wind scour, is also modeled heuristically with a diffusion wave equation. The multiple surface and subsurface processes are implemented by leveraging highly parallel community software. Fully integrated thermal hydrology simulations on the tilted open book catchment, an important test case for integrated surface/subsurface flow modeling, are presented. Fine-scale 100-year projections of the integrated permafrost thermal hydrological system on an ice wedge polygon at Barrow Alaska in a warming climate are also presented. Finally, these simulations demonstrate the feasibility of microtopography-resolving, process-rich simulations as a tool to help understand possible future evolution of the carbon-rich Arctic tundra in a warming climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Center for Frontiers of Subsurface Energy Security (Final Report)

For the past two decades a slowly growing consensus has been that rising temperature is being caused by increases in carbon dioxide (CO 2 ) in the Earth’s atmosphere. What is less consensual is what to do about it. Solutions range from reducing emissions to various ways to capture CO 2 (from effluents of power plants, as plants, etc.) and store it away from the atmosphere. The problem is of such magnitude that it is likely that we will need many methods. Also some methods would cause draconian harm to global economies. This work is about the storage part of the reducing of CO 2 in the atmosphere. There have been many methods proposed for storing CO 2 : in vegetation, at the bottom of oceans, in land farms, and in subsurface formations. The latter is particularly prominent given the large estimated volume in porous formation in the Earth, and the 50-year old experience in injecting CO 2 for enhanced oil recovery. This report is about the science behind subsurface storage, called here geological carbon storage or GCS. Effective GCS requires that billions of metric tons of CO 2 be injected into storage reservoirs annually. GCS feasibility requires balancing the opposing impacts of large emplacement (injection) rates without compromising caprock seal integrity during injection. Once injected, the challenge is to assure the injected CO 2 remains permanently stored. Under most conditions, CO 2 is in a supercritical state, or scCO 2 . scCO 2 is immiscible with the subsurface brine and has a factor of ten smaller viscosity, leading to the development of two-phase fluid flow and scCO 2 bypassing. After the injection of large volumes of scCO 2 , the pore pressure, state of stress, chemical, thermal, and biological steady-state or equilibrium conditions in the subsurface are disrupted, which may lead to unpredictable behavior because the responses to these changes are often non-linear. Current challenges for efficient and reliable scCO 2 storage are: (1) sustaining large injection rates; (2) using available pore space efficiently, and (3) controlling undesired or unexpected behavior, such as scCO 2 leakage. Various aspects of these challenges are addressed by collaborative projects within Center for Frontiers in Subsurface Energy Security (CFSES). Challenge 1 – Sustaining large storage rates. To offset the annual CO 2 emissions of in the United States, ~7 billion metric tons of CO 2 must be injected into the storage reservoirs annually. Injected CO 2 forms a plume that starts spreading, displacing originally resident brine, though inefficiently so. Because of its low viscosity, and buoyant forces, a CO 2 plume rises to the reservoir rock-caprock interface. Injecting large volumes of scCO 2 at sufficiently high rates into subsurface storage formations leads to increases in pore pressures, potential expansion of the reservoir rock, and may result in fracturing of the reservoir rock and/or caprock. CFSES research uses existing and new experimental and modeling approaches to identify chemical-physical controls on the permeability and pore pressure dynamics, anticipating geomechanical (fracturing) events, predicting multi-phase flow patterns and trapping of CO 2 . Challenge 2 – Using pore space with unprecedented efficiency. Since large volumes of scCO 2 are to be stored, the volume available in storage reservoirs must be used efficiently. Storage efficiency is defined as a fraction of pore space occupied by carbon dioxide. Current estimates are that less than 5 percent of available pore volume is available for storage because scCO 2 , having a smaller density and viscosity than reservoir brine, tends to form viscous fingers, and bypass available pore space. This flow regime results in a CO 2 plume spreading over large areas in subsurface. The research goal in Challenge 2 is to advance the fundamental understanding of two-phase flow through porous media, including molecular-to meter- scale experimental and modeling approaches. This work leads to novel strategies to ensure that injected CO 2 occupies more than 5% of the pore space, ideally up to 50%. Challenge 3 – Controlling undesired or unexpected behavior. Following the injection of scCO 2 , the subsurface reservoir re-equilibrates to different physical and chemical states. This re-equilibration in space and time could lead to unexpected and/or undesired behavior: fracturing of reservoir rock and/or caprock, development of preferential flow paths, CO 2 leakage (to the surface or overlaying fresh water reservoirs), induced seismicity, and activation of faults. One of the main research goals of CFSES was to advance the fundamental understanding of how individual chemical and physical processes are coupled. We particularly focused on the chemical-mechanical coupling, and flow-mechanical coupling. We collected experimental and modeling data to determine the temporal and spatial scales of the coupled processes, which may lead to pore collapse caused by creep, activation of fractures and faults, chemically-induced fracture propagation, and self-focusing of scCO 2 flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
DOE OSTI2019

Nitrogen-Driven Chromatographic Separation During Gas Injection Into Hydrate-Bearing Sediments

Hydrates are solid phases composed of water cages enclosing gas molecules that may host large quantities of recoverable natural gas and may serve to sequester CO 2 on geological time scales. Most hydrate studies focus on hydrates containing a single gas component, such as CH 4 or CO 2 . Yet, there are several settings in which multiple components form hydrate mixtures, or mixed hydrates, including a subsurface injection technique that claims to simultaneously recover CH 4 and sequester CO 2 called “guest molecule exchange.” Here, we combine multicomponent phase behavior for hydrate-forming systems with a multiphase fluid flow simulator to understand the evolution of hydrate and nonhydrate phases during subsurface injection. We simulate various scenarios for systems composed of H 2 O, CH 4 , CO 2 , and N 2 . Our study probes the impact of injection composition, initial reservoir composition, and transport of each component through the model domain. We observe chromatographic separation from the combined effect of compositional partitioning in each phase, variable flow speed of each phase, and compositional dependence of phase stabilities. We show that N 2 drives chromatographic separation to create a CH 4 -free zone and a CO 2 -free zone that are connected by a continuous N 2 -dominated vapor phase. While our results are theoretical and should be validated experimentally, they imply that guest molecule exchange acts as two sequential processes rather than as a simultaneous process. Furthermore, they show that injections into reservoirs with and without free water have vastly different behaviors, which has implications for the interpretation of the guest molecule exchange field test and various laboratory studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗