Geological Society of America

T193 Approaches to Applied Modeling for Industry and the Public

Tuesday, 24 October 2017
8:00 AM - 12:00 PM
Room 613/614

IMPLEMENTATION OF A DYNAMIC CAPILLARY SUCTION MODEL FOR NUMERICAL MODELING OF UNSATURATED SYSTEMS WITH CHANGING POROSITY

Paper number: T193-11
Program: Technical Sessions
Time: 10:40 AM – 10:55 AM

Authors

Peter Johnson Presenting Author
Los Alamos National Laboratory, Los Alamos, NM
Peter Johnson Presenting Author
The University at Buffalo, Amherst, NY
Philip H. Stauffer
Los Alamos National Laboratory, Los Alamos, NM
Suzanne Bourret
Los Alamos National Laboratory, Los Alamos, NM
George A. Zyvoloski
Los Alamos National Laboratory, Los Alamos, NM

Abstract

Many geological situations involve fluid flow through the unsaturated zone during which porosity of the medium changes due to chemical, thermal, or mechanical processes. This effect can be seen in a variety of scenarios, such as hydrothermal systems, soil development, diagenesis, and others. In some cases, opening or closing of pore space can be extensive, resulting in local porosity that changes greatly over time. This can present a problem for numerical simulations of flow in these systems. One concern in such simulations pertains to functions calculating capillary pressure, which are commonly specified for initial conditions or some average value during the course of a model run. Neglecting the effects of porosity on capillary suction can lead to unreasonable model results if porosity varies greatly over time from initial conditions.

In support of experiments and studies relating to disposal of high level radioactive waste in salt, we have conducted numerous model runs using the coupled heat/stress/mass flow model Finite Element Heat and Mass (FEHM, https://fehm.lanl.gov). Capillary forces tend to pull brine towards heat sources where evaporating brine causes element undersaturation. Precipitation of salt as brine evaporates decreases porosity, which should result in greater capillary forces for comparable saturation values than the initial specified conditions. Vapor condensation elsewhere causes local dissolution of salt, increasing porosity, and should decrease the importance of capillary effects in those areas. This caused problems during model runs where increasingly porous areas had unphysically high suction effects which kept them more saturated than expected and led to a further dissolution of the salt. As a result, portions of the model domain erroneously show porosity approaching 1 while retaining high saturation. To address this, we have implemented capillary pressure functions that calculate capillary pressure for each node within each timestep as a function of both porosity and saturation. Residual saturation and the maximum saturation above which capillary pressures go to zero are calculated as end members and a fit applied between them. This approach allows for dynamic capillary suction calculations that ameliorate problematic pore-capillary behaviors in the model.

Geological Society of America Abstracts with Programs. Vol. 49, No. 6

doi: 10.1130/abs/2017AM-305255

© Copyright 2017 The Geological Society of America (GSA), all rights reserved. Permission is hereby granted to the author(s) of this abstract to reproduce and distribute it freely, for noncommercial purposes. Permission is hereby granted to any individual scientist to download a single copy of this electronic file and reproduce up to 20 paper copies for noncommercial purposes advancing science and education, including classroom use, providing all reproductions include the complete content shown here, including the author information. All other forms of reproduction and/or transmittal are prohibited without written permission from GSA Copyright Permissions.