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 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Imaging Complex Subsurface Structures for Geothermal Exploration at Pirouette Mountain and Eleven-Mile Canyon in Nevada

Accurate imaging of subsurface complex structures with faults is crucial for geothermal exploration because faults are generally the primary conduit of hydrothermal flow. It is very challenging to image geothermal exploration areas because of complex geologic structures with various faults and noisy surface seismic data with strong and coherent ground-roll noise. In addition, fracture zones and most geologic formations behave as anisotropic media for seismic-wave propagation. Properly suppressing ground-roll noise and accounting for subsurface anisotropic properties are essential for high-resolution imaging of subsurface structures and faults for geothermal exploration. We develop a novel wavenumber-adaptive bandpass filter to suppress the ground-roll noise without affecting useful seismic signals. This filter adaptively exploits both characteristics of the lower frequency and the smaller velocity of the ground-roll noise than those of the signals. Consequently, this filter can effectively differentiate the ground-roll noise from the signal. We use our novel filter to attenuate the ground-roll noise in seismic data along five survey lines acquired by the U.S. Navy Geothermal Program Office at Pirouette Mountain and Eleven-Mile Canyon in Nevada, United States. We then apply our novel anisotropic least-squares reverse-time migration algorithm to the resulting data for imaging subsurface structures at the Pirouette Mountain and Eleven-Mile Canyon geothermal exploration areas. The migration method employs an efficient implicit wavefield-separation scheme to reduce image artifacts and improve the image quality. Our results demonstrate that our wavenumber-adaptive bandpass filtering method successfully suppresses the strong and coherent ground-roll noise in the land seismic data, and our anisotropic least-squares reverse-time migration produces high-resolution subsurface images of Pirouette Mountain and Eleven-Mile Canyon, facilitating accurate fault interpretation for geothermal exploration.

15 GEOTHERMAL ENERGY↗

DOE OSTI2013

Microscale Metabolic, Redox and Abiotic Reactions in Hanford 300 Area Subsurface Sediments

The Hanford 300 Area is a unique site due to periodic hydrologic influence of river water resulting in changes in groundwater elevation and flow direction. This area is also highly subject to uranium remobilization, the source of which is currently believed to be the region at the base of the vadose zone that is subject to period saturation due to the changes in the water levels in the Columbia River. We found that microbial processes and redox and abiotic reactions which operate at the microscale were critical to understanding factors controlling the macroscopic fate and transport of contaminants in the subsurface. The combined laboratory and field research showed how microscale conditions control uranium mobility and how biotic, abiotic and redox reactions relate to each other. Our findings extended the current knowledge to examine U(VI) reduction and immobilization using natural 300 Area communities as well as selected model organisms on redox-sensitive and redox-insensitive minerals. Using innovative techniques developed specifically to probe biogeochemical processes at the microscale, our research expanded our current understanding of the roles played by mineral surfaces, bacterial competition, and local biotic, abiotic and redox reaction rates on the reduction and immobilization of uranium.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Organizational Principles of Hyporheic Exchange Flow and Biogeochemical Cycling in River Networks Across Scales

Hyporheic zones increase freshwater ecosystem resilience to hydrological extremes and global environmental change. However, current conceptualizations of hyporheic exchange, residence time distributions, and the associated biogeochemical cycling in streambed sediments do not always accurately explain the hydrological and biogeochemical complexity observed in streams and rivers. Specifically, existing conceptual models insufficiently represent the coupled transport and reactivity along groundwater and surface water flow paths, the role of autochthonous organic matter in streambed biogeochemical functioning, and the feedbacks between surface-subsurface ecological processes, both within and across spatial and temporal scales. While simplified approaches to these issues are justifiable and necessary for transferability, the exclusion of important hyporheic processes from our conceptualizations can lead to erroneous conclusions and inadequate understanding and management of interconnected surface water and groundwater environments. This is particularly true at the landscape scale, where the organizational principles of spatio-temporal dynamics of hyporheic exchange flow (HEF) and biogeochemical processes remain largely uncharacterized. This article seeks to identify the most important drivers and controls of HEF and biogeochemical cycling based on a comprehensive synthesis of findings from a wide range of river systems. We use these observations to test current paradigms and conceptual models, discussing the interactions of local-to-regional hydrological, geomorphological, and ecological controls of hyporheic zone functioning. This improved conceptualization of the landscape organizational principles of drivers of HEF and biogeochemical processes from reach to catchment scales will inform future river research directions and watershed management strategies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Importance of river water recharge to the San Joaquin Valley groundwater system

Groundwater is not a sustainable resource, unless abstraction is balanced by recharge. Identifying the sources of recharge in a groundwater basin is critical for sustainable groundwater management. In this paper, we studied the importance of river water recharge to groundwater in the south-eastern San Joaquin Valley (24,000 km 2 , population 4 million). We combined dissolved noble gas concentrations, stable isotopes, tritium, and carbon-14 analyses to analyse the sources, mechanisms, and timescales of groundwater recharge. Area-representative groundwater sampling and numerical model input data enabled a stable isotope mass balance and quantitative estimates of river and local recharge. River recharge, identified by a lighter stable isotope signature, represents 47 ± 4% of modern groundwater in the San Joaquin Valley (recharged after 1950) but only 26 ± 4% of premodern groundwater (recharged before 1950). This implies that the importance of river water recharge in the San Joaquin valley has nearly doubled and is likely the result of a 40% increase in total recharge, caused by river water irrigation return flows and increased stream depletion and river recharge due to groundwater pumping. Compared with the large and long-duration capacity for water storage in the subsurface, storage of water in rivers is limited in time and volume, as evidenced by cold river recharge temperatures resulting from fast infiltration and recharge. Finally, groundwater banking of seasonal surface water flows and expansion of managed aquifer recharge practices therefore appear to be a natural and promising method for increasing the resilience of the San Joaquin Valley water supply system.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗

DOE OSTI2018

Nonnegative tensor factorization for contaminant source identification

Unsupervised Machine Learning (ML) is becoming increasingly popular for solving various types of data analytics problems including feature extraction, blind source separation, exploratory analyses, model diagnostics, etc. In this work, we have developed a new unsupervised ML method based on Nonnegative Tensor Factorization (NTF) for identification of the original groundwater types (including contaminant sources) present in geochemical mixtures observed in an aquifer. Frequently, groundwater types with different geochemical signatures are related to different background and/or contamination sources. The characterization of groundwater mixing processes is a challenging but very important task critical for any environmental management project aiming to characterize the fate and transport of contaminants in the subsurface and perform contaminant remediation. This task typically requires solving complex inverse models representing groundwater flow and geochemical transport in the aquifer, where the inverse analysis accounts for available site data. Usually, the model is calibrated against the available data characterizing the spatial and temporal distribution of the observed geochemical types. Numerous different geochemical constituents and processes may need to be simulated in these models which further complicates the analyses. Additionally, the application of inverse methods may introduce biases in the analyses through the assumptions made in the model development process. Here, we substitute the model inversion with unsupervised ML analysis. The ML analysis does not make any assumptions about underlying physical and geochemical processes occurring in the aquifer. Our ML methodology, called NTFk, is capable of identifying (1) the unknown number of groundwater types (contaminant sources) present in the aquifer, (2) the original geochemical concentrations (signatures) of these groundwater types and (3) spatial and temporal dynamics in the mixing of these groundwater types. These results are obtained only from the measured geochemical data without any additional site information. In general, the NTFk methodology allows for interpretation of large high-dimensional datasets representing diverse spatial and temporal components such as state variables and velocities. NTFk has been tested on synthetic and real-world site three-dimensional datasets. Finally, the NTFk algorithm is designed to work with geochemical data represented in the form of concentrations, ratios (of two constituents; for example, isotope ratios), and delta notations (standard normalized stable isotope ratios).

58 GEOSCIENCES↗

DOE OSTI2014

Full waveform inversion for time-distance helioseismology

Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions and is likely for supergranular flows. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the subsurface structure of large-amplitude perturbations. Defining the penalty functional as the L {sub 2} norm of wave travel-time deviations, we compute the total misfit gradient of this functional with respect to the relevant model parameters at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or the quasi-Newton limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Performing nonlinear iterations requires privileging pixels (such as those in the near field of the scatterer), a practice that is not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a two-dimensional inverse problem with only a sound-speed perturbation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Coupled Geomechanical Modeling to Assess Cap Rock Integrity and Mechanical Fault Stability: Application to Farnsworth Field Unit Project

Abstract The SWP project is located in a mature waterflood undergoing conversion to CO2-WAG operations in the Farnsworth, Texas, USA. Anthropogenic CO2 is sourced from a fertilizer and an ethanol plant. This work utilizes Farnsworth’s full-field, history-matched, compositional hydrodynamic coupled geomechanical model for assessing the impact of stress changes observed through the history matched field life. Production and injection induced stress changes, fault stability and caprock integrity investigations are performed to project the potential for fault reactivation and the loss of caprock integrity under shear failure. A static mechanical earth model (MEM) was constructed for use in transient coupled geomechanical model based on the existing Southwest Regional Partnership (SWP) geological model. The static MEM inherits the stratigraphic and structural features of the geologic model and incorporates additional overburden, underburden, and sideburden formations required to impose mechanical boundary conditions. Mechanical properties were distributed in the 3D MEM through integration of geophysical logs and 3-dimensional seismic elastic inversion properties using a combination of Bayesian and stochastic interpolation methods. These data are further enhanced by lab derived strength and failure criteria for the caprock interval. Additionally, interpreted faults and other geological features were included as part of the static structural framework to fully represent subsurface hydraulic and mechanical systems and appropriately integrate heterogeneity. Two-way coupling of hydrodynamic flow and geomechanical simulations incorporates the Kozeny-Carman relationship for updating permeability and is history matched through primary, secondary (waterflooding) and tertiary (CO2 WAG) recovery periods before performing a twenty year forecast. Two-way simulations are performed to understand the effective stress perturbations imposed by field operations: water injection and the more recent implementation of WAG. Evaluation of Mohr circles with liquid production and field pressure charts as well as slip tendency and distance to failure metrics indicate that neither our faults or cap rock are critically stressed. This paper presents the results of the Farnsworth initial attempts to integrate the seismic driven 3D MEM into coupled hydrodynamic geomechanical history match simulation workflow.

02 PETROLEUM↗

DOE OSTI2025

Integration and Demonstration of Monitoring, Modeling, and Prediction of DV-1 Amendment Performance at the Bench Scale: DV-1 Amendment Demonstration

During fiscal years 2024 and 2025, the U.S. Department of Energy’s Hanford Field Office commissioned Pacific Northwest National Laboratory to conduct applied research aimed at reducing the cost, time, and uncertainty associated with in situ treatment of vadose zone contaminants at the Hanford Site. This report outlines the integration of three key research efforts into a meso-scale demonstration designed to advance field-scale solutions that aim to (1) optimize the delivery of chemical amendments to contaminated soils, (2) reduce uncertainty in amendment delivery performance assessment using advanced monitoring techniques, and (3) provide real-time insights into when and where amendment-induced precipitation reactions occur in the subsurface. To achieve these objectives, the tank-scale (~ 1 cubic meter) Geophysical Imaging of Flow and Transport (GIFT) system was developed. GIFT enables experimental testing of amendment delivery while incorporating automated multi-modal monitoring approaches, including pressure measurements, direct fluid sampling, and remote time-lapse geophysical imaging. The data generated from these monitoring techniques will serve as inputs for a generative artificial-intelligence-driven digital twin – a numerical simulation model designed to honor observed data while quantifying uncertainty in simulation accuracy. Using this simulator, researchers will refine an amendment injection strategy to maximize delivery efficiency within a low-permeability soil zone. Monitoring data will be interpreted through simulated outputs to enhance understanding of the injection process. The efficacy of this integrated approach will be evaluated through direct sampling at the conclusion of the experiment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

System and method for optimizing hydrocarbon production from subsurface reservoirs

A method is described for subsurface hydrocarbon reservoir characterization including receiving a time-lapse electromagnetic (EM) dataset and a flow dataset; inverting the time-lapse EM dataset using a parametric inversion that models steel well casings to determine a volume of fluid-changed reservoir; inverting the time-lapse EM dataset and the flow dataset using a joint inversion that honors the volume of the fluid-changed reservoir to determine relative permeability and capillary pressure; and characterizing flow characteristics in the volume of the fluid-changed reservoir. The method may be executed by a computer system.

Hoversten, Gary Michael↗

DOE OSTI2013

Development of a coupled thermo-hydro-mechanical model in discontinuous media for carbon sequestration

Geomechanical alteration of porous media is generally ignored for most shallow subsurface applications, whereas CO 2 injection, migration, and trapping in deep saline aquifers will be controlled by coupled multifluid flow, energy transfer, and geomechanical processes. The accurate assessment of the risks associated with potential leakage of injected CO 2 and the design of effective injection systems requires that we represent these coupled processes within numerical simulators. The objectives of this study were to develop a coupled thermal-hydro-mechanical model into a single software, and to examine the coupling of thermal, hydrological, and geomechanical processes for simulation of CO 2 injection into the subsurface for carbon sequestration. A numerical model is developed to couple nonisothermal multiphase hydrological and geomechanical processes for prediction of multiple interconnected processes for carbon sequestration in deep saline aquifers. The geomechanics model was based on Rigid Body-Spring Model (RBSM), one of the discrete methods to model discontinuous rock system. Poisson’s effect that was often ignored by RBSM was considered in the model. The simulation of large-scale and long-term coupled processes in carbon capture and storage projects requires large memory and computational performance. Global Array Toolkit was used to build the model to permit the high performance simulations of the coupled processes. The model was used to simulate a case study with several scenarios to demonstrate the impacts of considering coupled processes and Poisson’s effect for the prediction of CO 2 sequestration.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Two-Phase Flow Mechanisms Controlling CO2 Intrusion into Shaly Caprock

Abstract Geologic carbon storage in deep saline aquifers has emerged as a promising technique to mitigate climate change. CO 2 is buoyant at the storage conditions and tends to float over the resident brine jeopardizing long-term containment goals. Therefore, the caprock sealing capacity is of great importance and requires detailed assessment. We perform supercritical CO 2 injection experiments on shaly caprock samples (intact caprock and fault zone) under representative subsurface conditions. We numerically simulate the experiments, satisfactorily reproducing the observed evolution trends. Simulation results highlight the dynamics of CO 2 flow through the specimens with implications to CO 2 leakage risk assessment in field practices. The large injection-induced overpressure drives CO 2 in free phase into the caprock specimens. However, the relative permeability increase following the drainage path is insufficient to provoke an effective advancement of the free-phase CO 2 . As a result, the bulk CO 2 front becomes almost immobile. This implies that the caprock sealing capacity is unlikely to be compromised by a rapid capillary breakthrough and the injected CO 2 does not penetrate deep into the caprock. In the long term, the intrinsically slow molecular diffusion appears to dominate the migration of CO 2 dissolved into brine. Nonetheless, the inherently tortuous nature of shaly caprock further holds back the diffusive flow, favoring safe underground storage of CO 2 over geological time scales.

58 GEOSCIENCES↗

DOE OSTI2016

Impact of mineral precipitation on flow and mixing in porous media determined by microcomputed tomography and MRI

Here, precipitation reactions in porous media influence transport properties of the environment and can control advective and dispersive transport. In subsurface environments, mixing of saline groundwater or injected solutions for remediation with fresh groundwater can induce supersaturation of constituents and drive precipitation reactions. Magnetic resonance imaging (MRI) and micro-computed tomography (µ-CT) were employed as complimentary techniques to evaluate advection, dispersion and formation of precipitate in a 3D porous media flow cell. Two parallel fluids were flowed concentrically through the porous media under two flow rate conditions with Na 2 CO 3 and CaCl 2 in the inner and outer fluids, respectively. Upon mixing, calcium carbonate became supersaturated and formed a precipitate at the interface of the two fluids. Spatial maps of changing local velocity fields and dispersion in the flow cell were generated from MRI, while high resolution imaging of the precipitate formed in the porous media was achieved via µ-CT imaging. Formation of a precipitate layer minimized dispersive and advective transport between the two fluids and the shape of the precipitation was influenced by the flow rate condition.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Fetch-trap Pairs: Exploring definition of carbon storage prospects to increase capacity and flexibility in areas with competing uses

As CCS grows beyond isolated first-mover projects, the question of how to identify potential storage sites (prospects) is key to siting projects amidst the constraints of neighboring storage projects, hydrocarbon fields and other considerations. To date, prospects have been defined either by the extent of a buoyant trap or by the predicted maximum extent of the CO 2 or pressure plume. Both approaches have worked, but they may struggle in places like the Gulf of Mexico where complex geologic structure can invalidate simple assumptions of radial plume spread and buoyant traps are commonly unattractive due to numerous legacy well penetrations. Here, we propose identifying prospective storage sites by fetch area (drainage cell) rather than by buoyant closure alone. Doing so offers 1) greater freedom in siting injectors to avoid surface and subsurface constraints, including legacy wells that tend to cluster on structural high; and 2) a coherent flow regime in which buoyancy drives all injected CO 2 toward a common high. This strategy gives space to dissipate injection pressure and minimize the number of legacy wells needing review. It may also offer 1) better injectivity by bringing synclines and associated channel axes into play; and 2) potentially improved capacity by tapping a larger rock volume and taking advantage of migration losses.

58 GEOSCIENCES↗

DOE OSTI2015

3D Seismic Experimentation and Advanced Processing/Inversion Development for Investigations of the Shallow Subsurface

The work plan for this project was to develop and apply advanced seismic reflection and wide-angle processing and inversion techniques to high resolution seismic data for the shallow subsurface to seismically characterize the shallow subsurface at hazardous waste sites as an aid to containment and cleanup activities. We proposed to continue work on seismic data that we had already acquired under a previous DoE grant, as well as to acquire additional new datasets for analysis. The project successfully developed and/or implemented the use of 3D reflection seismology algorithms, waveform tomography and finite-frequency tomography using compressional and shear waves for high resolution characterization of the shallow subsurface at two waste sites. These two sites have markedly different near-surface structures, groundwater flow patterns, and hazardous waste problems. This is documented in the list of refereed documents, conference proceedings, and Rice graduate theses, listed below.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Experimental Investigation of Two-Phase Flow in Consolidating and Intact Salt: Analysis of Capillary Pressure Measurements from Core Labs

Two-phase flow of brine and air through very low permeability porous media is the dominant fluid flow state in and around geologic repositories for nuclear waste disposal. Constitutive models developed for other subsurface media (shale, sandstone, limestone, etc.) have been adopted but never experimentally parameterized for geologic salt. This report summarizes results of scoping capillary pressure and brine saturation measurements on reconsolidated run-of-mine salt and intact geologic salt samples. These measurements were obtained using mercury injection tests and yielded values of threshold pressure, pore-size distribution, and mercury saturation. The brine capillary pressure curve was then inferred from the raw mercury capillary pressure data. This report also summarizes estimates of the constitutive model parameters from the brine saturation and capillary pressure data (Brooks and Corey, 1964; van Genuchten, 1980). The results suggests that the model of van Genuchten (1980) is better suited for brine-air-salt systems. Additional testing is needed for the low brine saturation range of the capillary pressure curve where the present suite of tests did not yield sufficient data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗