Aqueous Relative Permeability Options (GT)
Relative Permeability Options
Constant
The relative permeability is constant regardless of saturation.
Mualem
Modified Mualem
The pore scale parameter (aka tortuosity-connectivity coefficient) can be specified to be any value. The default value is 0.5.
Irreducible Mualem
There is a minimum saturation for the porous medium.
Mualem with Polmann
The Mualem (1976) relative permeability model is used to calculate the vertical permeability and the Polmann (1990) model is used to calculate horizontal permeability.
Polmann (1990) used a stochastic model to evaluate tension-dependent anisotropy. The Gardner (1958) relationship was used to describe unsaturated hydraulic conductivity as a function of saturated hydraulic conductivity and tension. Using a linear correlation model between the zero-tension intercept and β, Polmann (1990) presented a generalized model that accounts for the cross-correlation of the local soil property (i.e., lnKs and β) residual fluctuations. Compared against the uncorrelated lnKs and β model, the partial correlation of the properties was shown to have a significant impact on the magnitude of the effective parameters derived from the stochastic theory. The Polmann (1990) equations for deriving the effective parameters for strongly stratified porous media are as shown below.
A modified version of the original Polmann model is implemented in STOMP. This version uses a single segment of the log-linear K(h) function. Because the Polmann model is based on the log-linear hydraulic function, it may be applicable only to a relatively narrow pressure head range over which the parameters are derived. Extrapolating the derived hydraulic function to other pressure head conditions may cause significant error or even physically incorrect results (e.g., the predicted hydraulic conductivity in the horizontal direction may increase with decreasing saturation, which is physically incorrect). It calculates an anisotropy coefficient (Can) and the hydraulic conductivity in the horizontal direction (Kh). The modified model also requires user specified maximum (Canmax) and minimum (Canmin) limits of the anisotropy coefficient. If computed values of Can are outside these limits, the values are truncated as shown in the equations below:
Symbols
|
---|
|
equivalent unsaturated hydraulic conductivity in the horizontal direction, m/s |
|
equivalent unsaturated hydraulic conductivity in the vertical direction, m/s |
|
unsaturated hydraulic conductivity (pressure dependent), m/s |
|
saturated hydraulic conductivity, m/s |
|
mean of LnKs |
|
capillary head, m |
|
standard deviation of the residuals in the LnKs vs. β regression for all the samples considered |
|
standard deviation of LnKs |
|
variance of LnKs |
|
slope between LnKs and h |
|
slope of LnKs vs. β regression line for the samples considered |
|
mean slope of β for LnK vs. h regression |
|
vertical correlation length for LnKs |
|
anisotropy coefficient |
|
maximum value of anisotropy coefficient |
|
minimum value of anisotropy coefficient |
The Directional Aqueous Relative Permeability Card can be used to represent anisotropy for several aqueous relative permeability models. Zhang et al. 2003 and Raats et al. 2004 demonstrate the use of Directional Aqueous Relative Permeability with the Modified Mualem model.
Fatt and Klikoff
Corey Model
Free Corey
Haverkamp
Touma and Vauclin
Tabular
This option accepts tabulated relative permeability data. The default is the data of aqueous saturation and aqueous relative permeability data pairs and linear interpolation is used between data points; a cubic spline interpolation scheme can also be specified. Alternately, capillary head vs relative permeability can be provided using the keyword "head." Other interpolation schemes that can be specified for capillary head vs. aqueous relative permeability are log capillary head vs relative permeability, cubic spline, or cubic spline for log capillary head vs relative permeability.
Sub-Options
Dual Porosity/Permeability Model for Fractured Systems
IJK, JKI or KIJ Indexing
If IJK Indexing, JKI Indexing, or KIJ Indexing is specified as the Rock/Soil Name in the Rock/Soil Zonation Card, then this must also be specified as the Rock/Soil Name in the Aqueous Relative Permeability Card.
IJK Indexing
If the IJK Indexing option is specified in the Rock/Soil Zonation Card, then the relative permeability models and any or all associated parameters can be specified either as a single value that will be applied to each node in the domain, or in an external file with the values for every grid-cell ordered according to the IJK indexing scheme. Units shown in the input line will be applied to all parameters in the external file.
JKI Indexing
If the JKI Indexing option is specified in the Rock/Soil Zonation Card, then the relative permeability models and any or all associated parameters can be specified either as a single value that will be applied to each node in the domain, or in an external file with the values for every grid-cell ordered according to the JKI indexing scheme. Units shown in the input line will be applied to all parameters in the external file.
KIJ Indexing
If the KIJ Indexing option is specified in the Rock/Soil Zonation Card, then the relative permeability models and any or all associated parameters can be specified either as a single value that will be applied to each node in the domain, or in an external file with the values for every grid-cell ordered according to the KIJ indexing scheme. Units shown in the input line will be applied to all parameters in the external file.
Effective Aqueous Liquid Saturation
Mobile saturations are scaled by the residual aqueous liquid saturation to determine effective saturations for use with common relative permeability models.
Symbols
|
---|
|
effective aqueous liquid saturation |
|
actual aqueous liquid saturation |
|
actual aqueous liquid residual saturation |
References:
Burdine, NT. 1953. "Relative Permeability Calculation from Size Distribution Data," Trans. AIME, 198:71-78.
Corey, AT. 1977. Mechanics of Heterogeneous Fluids in Porous Media, Fort Collins, Colorado.
Fatt , I and WA Klikoff Jr. 1959. "Effect of Fractional Wettability on Multiphase Flow through Porous Media," AIME Trans, 215:426-429.
Gardner, WR. 1958. "Some Steady-State Solutions of the Unsaturated Moisture Flow Equation with Applications to Evaporation from a Water Table," Soil Science, 85:228-232.
Haverkamp, R, M Vauclin, J Touma, PJ Wierenga, and G Vachaud. 1977. "A Comparison of Numerical Simulation Models for One-Dimensional Infiltration," Soil Sci. Soc. Am. J., 41:285-294.
Hoffmann-Riem, H, MT Van Genuchten, and H Fluhler. 1999. A General Model of the Hydraulic Conductivity of Unsaturated Soils. In M.T. Van Genuchten, L.J. Leij and L. Wu Proceedings of Int. Workshop, Characterization And Measurements Of The Hydraulic Properties of UnsaturatedPorous Media, University of California Riverside, Riverside, CA.
Klavetter, EA and RR Peters. 1986. Estimation of Hydrologic Properties of Unsaturated Fractured Rock Mass, SAND84-2642, Sandia National Laboratories, Albuquerque, NM.
Mualem, Y. 1976. "A New Model for Predicting the Hydraulic Conductivity of Unsaturated Porous Media," Water Resources Research, 12:513-522.
Nitao, JJ. 1988. Numerical Modeling of the Thermal and Hydrological Environment around a Nuclear Waste Package Using the Equivalent Continuum Approximation: Horizontal Emplacement, UCID-2144, Lawrence Livermore National Laboratory, Livermore, CA.
Parker, JC and RJ Lenhard. 1987. "A Model for Hysteretic Constitutive Relations Governing Multiphase Flow 1. Saturation-Pressure Relations," Water Resources Research, 23(12):2187-2196.
Polmann, DJ. 1990. "Application of Stochastic Methods to Transient Flow and Transport in Heterogeneous Unsaturated Soils," PhD, Massachusetts Institute of Technology, Cambridge, MA.
Raats, PaC, ZF Zhang, AL Ward, and GW Gee. 2004. "The Relative Connectivity–Tortuosity Tensor for Conduction of Water in Anisotropic Unsaturated Soils," Vadose Zone Journal, 3:1471-1478.
Zhang, ZF, AL Ward, and GW Gee. 2003. "A Tensorial Connectivity–Tortuosity Concept to Describe the Unsaturated Hydraulic Properties of Anisotropic Soils," Vadose Zone Journal, 2:313-321.