Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unsaturated 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Integrated systems-based approach to monitoring environmental remediation

The U.S. Department of Energy is responsible for risk reduction and cleanup of its nuclear weapons complex. Remediation strategies for some of the contamination may include techniques that mitigate risk, but leave contaminants in place. Monitoring to verify remedy performance and long-term mitigation of risk is key to implementing these strategies and can be a large portion of the total cost of remedy implementation. Especially in these situations, there is a need for innovative monitoring approaches that move away from the cost- and labor-intensive point-source monitoring. In this paper, alternative approaches for monitoring are presented for vadose zone, groundwater, groundwater/surface water interface, and surface water. To illustrate integrated, systems-based monitoring, this paper focuses on vadose zone contaminant remediation to mitigate impact to groundwater. In this context, vadose zone contamination is a source, or potential source, to groundwater plumes. The monitoring design uses a systems-based approach focused on developing a conceptual site model that highlights key features that control contaminant flux to groundwater. These features are derived considering the unsaturated flow and contaminant transport processes in the vadose zone and the nature of the waste discharge. Diagnostic properties and/or parameters related to both short- and long-term contaminant flux to groundwater can be identified and targeted for monitoring. The resolution of monitoring data needed to correspond to a functionally useful indicator of flux to groundwater can be estimated using quantitative analyses and the associated unsaturated flow properties relevant to the targeted site and vadose zone features. This monitoring design approach follows the process of developing a quantitative conceptual model suitable for supporting projections of future flux to groundwater. Support for such projections is important because it is likely that, in many cases, remediation decisions for the vadose zone will need to be made based all or in part on projected impacts to groundwater, and monitoring will then be applied to verify that remedy goals are being met. (authors)

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

Numerical Evaluation of Unsaturated-Zone Flow and Transport Pathways at Rainier Mesa, Nevada

Sixty-one underground nuclear tests were conducted at Rainier Mesa at the Nevada National Security Site between 1957 and 1992. The mesa includes the highest-elevation areas at the Nevada National Security Site and experiences some of the highest estimated infiltration rates. Perched water is ubiquitous in many of the tunnels used for testing. The presence of sloping layers with strongly varying degrees of fracturing and overall permeability suggests that some infiltrating water may be diverted laterally along higher-permeability layers, thereby reducing overall percolation flux and radionuclide transport through the low-permeability testing horizons but also potentially transporting radionuclides laterally away from the test locations. A high-resolution, two-dimensional, cross-sectional model was created to simulate water and particle movement along potential lateral flow paths and to investigate flow and transport paths within and surrounding the perched water bodies. The model results indicate that particles with starting locations within the low-permeability perching layers will follow dominantly vertical trajectories through the perching horizons and reach the water table in the underlying carbonate aquifer only after many thousands of years. In contrast, particles starting near ground surface above the testing horizon are mostly diverted laterally before reaching the testing horizon, except where nearby faults through the perching layers are present. Where the rock above the perching horizon is fractured, particles originating near ground surface can move laterally with relatively high transport velocities. Transport calculations for 3 H show that transport extent is limited by slow matrix flow, radioactive decay, and matrix diffusion in the fractured units.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Coupling of Nuclear Waste Form Corrosion and Radionuclide Transports in Presence of Relevant Repository Sediments

Assessments of waste form and disposal options start with the degradation of the waste forms and consequent mobilization of radionuclides. Long-term static tests, single-pass flow-through tests, and the pressurized unsaturated flow test are often employed to study the durability of potential waste forms and to help create models that predict their durability throughout the lifespan of the disposal site. These tests involve the corrosion of the material in the presence of various leachants, with different experimental designs yielding desired information about the behavior of the material. Though these tests have proved instrumental in elucidating various mechanisms responsible for material corrosion, the chemical environment to which the material is subject is often not representative of a potential radioactive waste repository where factors such as pH and leachant composition will be controlled by the near-field environment. Near-field materials include, but are not limited to, the original engineered barriers, their resulting corrosion products, backfill materials, and the natural host rock. For an accurate performance assessment of a nuclear waste repository, realistic waste corrosion experimental data ought to be modeled to allow for a better understanding of waste form corrosion mechanisms and the effect of immediate geochemical environment on these mechanisms. Additionally, the migration of radionuclides in the resulting chemical environment during and after waste form corrosion must be quantified and mechanisms responsible for migrations understood. The goal of this research was to understand the mechanisms responsible for waste form corrosion in the presence of relevant repository sediments to allow for accurate radionuclide migration quantifications. The rationale for this work is that a better understanding of waste form corrosion in relevant systems will enable increased reliance on waste form performance in repository environments and potentially decrease the need for expensive engineered barriers.Our current work aims are 1) quantifying and understanding the processes associated with glass alteration in contact with Fe-bearing materials; 2) quantifying and understanding the processes associated with glass alteration in presence of MgO (example of engineered barrier used in WIPP); 3) identifying glass alteration suppressants and the processes involved to reach glass alteration suppression; 4) quantifying and understanding the processes associated with Saltstone and Cast Stone (SRS and Hanford cementitious waste forms) in various representative groundwaters; 5) investigating positron annihilation as a new tool for the study of glass alteration; and 6) quantifying and understanding the processes associated with glass alteration under gamma irradiation.

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

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Soil Water Retention and Relative Permeability for Conditions from Oven-Dry to Full Saturation

Common conceptual models for unsaturated flow often rely on the oversimplified representation of medium pores as a bundle of cylindrical capillaries and assume that the matric potential is attributed to the capillary force only. The adsorptive surface forces are ignored. It is often assumed that aqueous flow is negligible when a soil is near or at the residual water content. These models are successful at high and medium water contents but often give poor results at low water contents. These models do not apply to conditions at which the water content is less than the residual water content. We extend the lower bound of existing water-retention functions and conductivity models from residual water content to the oven-dry condition (i.e., zero water content) by defining a state-dependent, residual-water content for a soil drier than a critical value. Furthermore, a hydraulic conductivity model for smooth uniform spheres was modified by introducing a correction factor to describe the film flow-induced hydraulic conductivity for natural porous media. The total unsaturated hydraulic conductivity is the sum of those due to capillary and film flow. The extended retention and conductivity models were verified measurements. Results show that, when the soil is at high and intermediate water content, there is no difference between the un-extended and the extended models; when the soil is at low water content, the un-extended models overestimate the water content but underestimate the conductivity. The extended models match the retention and conductivity measurements well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Performance of the Fluidized Bed Steam Reforming Product Under Hydraulically Unsaturated Conditions

Currently, several candidates for secondary waste immobilization at the Hanford site in the State of Washington, USA are being considered. To demonstrate the durability of the product in the unsaturated Integrated Disposal Facility (IDF) at the site, a series of tests have been performed one of the candidate materials using the Pressurized Unsaturated Flow (PUF) system. The material that was tested was the Fluidized Bed Steam Reformer (FBSR) granular product and the granular product encapsulated in a geopolymer matrix. The FBSR product is composed primarily of an insoluble sodium aluminosilicate matrix with the dominant phases being feldspathoid minerals mostly nepheline, sodalite, and nosean. The PUF test method allows for the accelerated weathering of materials, including radioactive waste forms, under hydraulically unsaturated conditions, thus mimicking the open-flow and transport properties that most likely will be present at the IDF. The experiments show a trend of decreasing tracer release as a function of time for several of the elements released from the material including Na, Si, Al, and Cs. However, some of the elements, notably I and Re, show a steady release throughout the yearlong test. This result suggests that the release of these minerals from the sodalite cage occurs at a different rate compared with the dissolution of the predominant nepheline phase.

Durability testing↗
DOE OSTI2022

Enhanced Transport of TiO 2 in Unsaturated Sand and Soil after Release from Biodegradable Plastic during Composting

Biodegradable plastics can reach full degradation when disposed of appropriately and thus alleviate plastic pollution caused by conventional plastics. But additives can be released into the environment during degradation, and the fate of these additives can be affected by the degradation process. In this study, we characterized TiO 2 particles released from a biodegradable plastic mulch during composting and studied the transport of the released TiO 2 particles in inert sand and agricultural soil columns under unsaturated flow conditions. TiO 2 microparticles were released from the biodegradable plastic mulch in both single particle and cluster forms. The released TiO 2 microparticles were fully retained in unsaturated soil columns due to its ne texture, slow flow rate, and favorable attachment conditions. But in unsaturated sand columns, the released TiO 2 microparticles were highly mobile. Additionally, comparison with the pristine TiO 2 microparticles revealed that the enhanced mobility of the released TiO 2 microparticles was caused by the presence of compost residues which blocked attachment sites on the sand surface. This study demonstrates that TiO 2 microparticles can be released during composting of biodegradable plastics, and the transport potential of the released TiO 2 microparticles in the terrestrial environment can be enhanced by the presence of compost material.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Performance of the Fluidized Bed Steam Reforming product under hydraulically unsaturated conditions

Several candidates for supplemental low-activity waste (LAW) immobilization at the Hanford site in Washington State, USA are being considered. One waste sequestering technology considered is Fluidized Bed Steam Reforming (FBSR). The granular product resulting from the FBSR process is composed primarily of an insoluble sodium aluminosilicate matrix with the dominant phases being feldspathoid minerals with a 1:1:1 molar ratio of Na, Al and Si. To demonstrate the durability of the product, which can be disposed of at the unsaturated Integrated Disposal Facility (IDF) at Hanford, a series of tests has been performed using the Pressurized Unsaturated Flow (PUF) system, which allows for the accelerated weathering of the solid materials. The system maintains hydraulically unsaturated conditions, thus mimicking the open-flow and transport properties that will be present at the IDF. Two materials were tested using the system: 1) the FBSR granular product and 2) the FBSR granular product encapsulated in a geopolymer to form a monolith. Results of the experiments show a trend of relatively constant effluent concentration of Na, Si, Al, and Cs as a function of time from both materials. The elements I and Re show a steady release throughout the yearlong test from the granular material but their concentrations seem to be increasing at one year from the monolith material. This result suggests that these two elements may be present in the sodalite cage structure rather than in the predominant nepheline phase because their release occurs at a different rate compared to nepheline phase. Also, these elements to not seem to reprecipitate when released from the starting material. Calculated one-year release rates for Si are on the order of 10 6 g/(m2 d) for the granular material and 10 5 g/(m2 d) for the monolith material while Re release is seen to be two orders of magnitude higher than Si release rates. SEM imaging and XRD analysis show how the alteration of the two materials is dependent on their depth in the column. This phenomenom is a result of depth-dependent solution concentrations giving rise chemical environments that may be supersaturated with respect to a number of mineral phases.

Fluidized Bed Steam Reforming↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
NASA NTRS1987

Programming a hillslope water movement model on the MPP

A physically based numerical model was developed of heat and moisture flow within a hillslope on a parallel architecture computer, as a precursor to a model of a complete catchment. Moisture flow within a catchment includes evaporation, overland flow, flow in unsaturated soil, and flow in saturated soil. Because of the empirical evidence that moisture flow in unsaturated soil is mainly in the vertical direction, flow in the unsaturated zone can be modeled as a series of one dimensional columns. This initial version of the hillslope model includes evaporation and a single column of one dimensional unsaturated zone flow. This case has already been solved on an IBM 3081 computer and is now being applied to the massively parallel processor architecture so as to make the extension to the one dimensional case easier and to check the problems and benefits of using a parallel architecture machine.

Devaney, J. E.↗
DOE OSTI2014

Assessment of Contaminated Brine Fate and Transport in MB139 at WIPP

Following the radionuclide release event of February 14, 2014 at the Waste Isolation Pilot Plant (WIPP), actinide contamination has been found on the walls and floor in Panel 7 as a result of a release in Room 7 of Panel 7. It has been proposed to decontaminate Panel 7 at the WIPP by washing contaminated surfaces in the underground with fresh water. A cost-effective cleanup of this contamination would allow for a timely return to waste disposal operations at WIPP. It is expected that the fresh water used to decontaminate Panel 7 will flow as contaminated brine down into the porosity of the materials under the floor – the run-of-mine (ROM) salt above Marker Bed 139 (MB139) and MB139 itself – where its fate will be controlled by the hydraulic and transport properties of MB139. Due to the structural dip of MB139, it is unlikely that this brine would migrate northward towards the Waste-Handling Shaft sump. A few strategically placed shallow small-diameter observation boreholes straddling MB139 would allow for monitoring the flow and fate of this brine after decontamination. Additionally, given that flow through the compacted ROM salt floor and in MB139 would occur under unsaturated (or two-phase) conditions, there is a need to measure the unsaturated flow properties of crushed WIPP salt and salt from the disturbed rock zone (DRZ).

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

Simulator for Hydrologic Unstructured Domains (SHUD v1.0): numerical modeling of watershed hydrology with the finite volume method

Abstract. Hydrologic modeling is an essential strategy for understanding and predicting natural flows, particularly where observations are lacking in either space or time or where complex terrain leads to a disconnect in the characteristic time and space scales of overland and groundwater flow. However, significant difficulties remain for the development of efficient and extensible modeling systems that operate robustly across complex regions. This paper introduces the Simulator for Hydrologic Unstructured Domains (SHUD), an integrated, multiprocess, multiscale, flexible-time-step model, in which hydrologic processes are fully coupled using the finite volume method. SHUD integrates overland flow, snow accumulation/melt, evapotranspiration, subsurface flow, groundwater flow, and river routing, thus allowing physical processes in general watersheds to be realistically captured. SHUD incorporates one-dimensional unsaturated flow, two-dimensional groundwater flow, and a fully connected river channel network with hillslopes supporting overland flow and baseflow. The paper introduces the design of SHUD, from the conceptual and mathematical description of hydrologic processes in a watershed to the model's computational structures. To demonstrate and validate the model performance, we employ three hydrologic experiments: the V-catchment experiment, Vauclin's experiment, and a model study of the Cache Creek Watershed in northern California. Ongoing applications of the SHUD model include hydrologic analyses of hillslope to regional scales (1 m2 to 106 km2), water resource and stormwater management, and interdisciplinary research for questions in limnology, agriculture, geochemistry, geomorphology, water quality, ecology, climate and land-use change. The strength of SHUD is its flexibility as a scientific and resource evaluation tool where modeling and simulation are required.

58 GEOSCIENCES↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Deep subsurface drip irrigation using coal-bed sodic water: Part I. Water and solute movement

Water co-produced with coal-bed methane (CBM) in the semi-arid Powder River Basin of Wyoming and Montana commonly has relatively low salinity and high sodium adsorption ratios that can degrade soil permeability where used for irrigation. Nevertheless, a desire to derive beneficial use from the water and a need to dispose of large volumes of it have motivated the design of a deep subsurface drip irrigation (SDI) system capable of utilizing that water. Drip tubing is buried 92 cm deep and irrigates at a relatively constant rate year-round, while evapotranspiration by the alfalfa and grass crops grown is seasonal. We use field data from two sites and computer simulations of unsaturated flow to understand water and solute movements in the SDI fields. Combined irrigation and precipitation exceed potential evapotranspiration by 300–480 mm annually. Initially, excess water contributes to increased storage in the unsaturated zone, and then drainage causes cyclical rises in the water table beneath the fields. Native chloride and nitrate below 200 cm depth are leached by the drainage. Some CBM water moves upward from the drip tubing, drawn by drier conditions above. Chloride from CBM water accumulates there as root uptake removes the water. Year over year accumulations indicated by computer simulations illustrate that infiltration of precipitation water from the surface only partially leaches such accumulations away. Field data show that 7% and 27% of added chloride has accumulated above the drip tubing in an alfalfa and grass field, respectively, following 6 years of irrigation. Maximum chloride concentrations in the alfalfa field are around 45 cm depth but reach the surface in parts of the grass field, illustrating differences driven by crop physiology. Deep SDI offers a means of utilizing marginal quality irrigation waters and managing the accumulation of their associated solutes in the crop rooting zone.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2012

TOUGHREACT-Pitzer V1.21

TOUGHREACT is a numerical simulation program for chemically reactive non-isothermal flows of multiphase fluids in porous and fractured media. The program is written in Fortran 77 and was developed by introducing reactive chemistry into the multiphase flow code TOUGH2 V2. Interactions between mineral assemblages and fluids can occur under local equilibrium or kinetic rates. The gas phase can be chemically active. Precipitation and dissolution reactions can change formation porosity and permeability, and can also modify the unsaturated flow properties of the rock. The code is distributed with a comprehensive user's guide that includes sample problems addressing geothermal reservoirs and hydrothermal systems, nuclear waste isolation, groundwater quality, sequestration of carbon dioxide in saline aquifers, and supergene copper enrichment. The TOUGHREACT-Pitzer V1.21 package supports the following fluid property modules packaged into TOUGH2 V2.0 (CR-1574): EOS1, EOS2, EOS3, EOS4, EOS7, EOS9, and ECO2N. This version incorporates a Pitzer ion interaction model for application to saline solutions. It is NOT compatible with the most recent TOUGHREACT V2 and relies on the previous TOUGHREACT V1.21.

Zhang, Guoxiang↗
DOE OSTI2011

TOUGHREACT-Pitzer V1.21

TOUGHREACT is a numerical simulation program for chemically reactive non-isothermal flows of multiphase fluids in porous and fractured media. The program is written in Fortran 77 and was developed by introducing reactive chemistry into the multiphase flow code TOUGH2 V2. Interactions between mineral assemblages and fluids can occur under local equilibrium or kinetic rates. The gas phase can be chemically active. Precipitation and dissolution reactions can change formation porosity and permeability, and can also modify the unsaturated flow properties of the rock. The code is distributed with a comprehensive user?s guide that includes sample problems addressing geothermal reservoirs and hydrothermal systems, nuclear waste isolation, groundwater quality, sequestration of carbon dioxide in saline aquifers, and supergene copper enrichment.

Source record↗
DOE OSTI2017

Current Efforts in Immobilized Low-Activity Waste Glass Testing to Support the Hanford Site Integrated Disposal Facility - 17119

Cleanup plans for Hanford site high-level radioactive mixed tank waste calls for separating the waste into high-level waste (HLW) and low-activity waste (LAW) fractions. The fractions will be vitrified separately at the Hanford Tank Waste Treatment and Immobilization Plant (WTP). The HLW glass will be sent to an offsite geological disposal facility while the immobilized low-activity waste (ILAW) glass will be disposed at the onsite Integrated Disposal Facility (IDF). In order to inform models used to predict the performance of ILAW glass at IDF, it is necessary to understand the corrosion behavior of the glass. In this paper, we discuss tests that account for the effects of 1) solution chemistry, 2) the ion exchange reaction, and 3) secondary phase formation on the glass dissolution rate. To quantify the effects of solution chemistry, specifically changes in pH, T, and Si concentration, the single-pass flow-through (SPFT) test method is applied. A designed set of these tests results in a set of kinetic rate law parameters that are used in a glass dissolution rate law based on transition-state theory. Furthermore, for SPFT tests with relatively high concentrations of Si, an ion exchange rate is provided. Additionally, through tests such as the pressurized unsaturated flow (PUF) test method, a suite of secondary phases can be compiled. The collection of kinetic rate law parameters, the ion exchange term, and the suite of secondary phases are provided to IDF performance assessment (PA) modelers. Lastly, we present a description of the efforts made to verify and validate the use of these parameters and secondary phases through the use of a field lysimeter present at the Hanford site. (authors)

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

Water Saturation Relations and Their Diffusion-Limited Equilibration in Gas Shale: Implications for Gas Flow in Unconventional Reservoirs

Large volumes of water are used for hydraulic fracturing of low permeability shale reservoirs to stimulate gas production, with most of the water remaining unrecovered and distributed in a poorly understood manner within stimulated regions. Because water partitioning into shale pores controls gas release, we measured the water saturation dependence on relative humidity (rh) and capillary pressure (P c ) for imbibition (adsorption) as well as drainage (desorption) on samples of Woodford Shale. Experiments and modeling of water vapor adsorption into shale laminae at rh = 0.31 demonstrated that long times are needed to characterize equilibrium in larger (5 mm thick) pieces of shales, and yielded effective diffusion coefficients from 9 × 10 -9 to 3 × 10 -8 m 2 s -1 , similar in magnitude to the literature values for typical low porosity and low permeability rocks. Most of the experiments, conducted at 50°C on crushed shale grains in order to facilitate rapid equilibration, showed significant saturation hysteresis, and that very large P c (~1 MPa) are required to drain the shales. These results quantify the severity of the water blocking problem, and suggest that gas production from unconventional reservoirs is largely associated with stimulated regions that have had little or no exposure to injected water. Finally, gravity drainage of water from fractures residing above horizontal wells reconciles gas production in the presence of largely unrecovered injected water, and is discussed in the broader context of unsaturated flow in fractures.

58 GEOSCIENCES↗