Saturation Function Syntax (EOR-BO)

STOMP-EOR considers three mobile phases: 1) aqueous, 2) nonaqueous-liquid, and 3 gas, each with a unique relative permeability. Relative permeability often varies by orders of magnitude with saturation and often is the most nonlinear component of the governing equations and constitutive relationships.  Gas is considered to be the non-wetting fluid and the liquid phases are considered to be wetting.  The wettability of the aqueous and non-aqueous liquids can vary between fully aqueous-wet to fully nonaqueous-liquid-wet and all degrees of mixed wet.  The wettability is controlled via the contact angle.  Zero contact angle yields a fully aqueous-wet system, a 90˚ contact angle yields mixed wettability system, and a 180˚ contact angle yields a nonaqueous-liquid wet system.  The saturation function and associated parameters are considered to be intrinsic properties of the rock/soil.  As such the saturation functions are specified in reference to rock/soil types.  The current version of the simulator allows the user to select one of two base saturation versus capillary head functions: 1) van Genuchten^[1] or 2) Brooks and Corey^[2], and whether the function is extended below residual saturation via the Webb^[3] function.

Card Title

~Saturation Function Card,

 Combination Model, Binary Models and Parameters

Repeat

Repeat for each rock/soil type 

 

Info

IJK, JKI, or KIJ Indexing

Expand for Syntax

If Rock/Soil Name is IJK Indexing, JKI Indexing, or KIJ Indexing then a parameter value can replaced with an external file using the following formatting for ASCII files:

...,file:file name,...

or binary files:

...,binary file:file name,...

The external file will contain parameter values for every grid-cell ordered according to the indexing scheme (i.e., IJK, JKI, or KIJ). Units shown in the input line will be applied to all parameters in the external file.

Choose One

Chose one Base Saturation Function Model

van Genuchten^[1] Model

Rock/Soil Name,van Genuchten,α Parameter,Units (1/m),n Parameter,m Parameter,Aqueous Residual Saturation,Nonaqueous-Liquid Residual Saturation,Contact Angle,Mixed-Wet b Parameter,

 van Genuchten^[1] Model with the Webb^[3] Extension

Rock/Soil Name,van Genuchten with Webb,α Parameter,Units (1/m),n Parameter,m Parameter,Aqueous Residual Saturation,Nonaqueous-Liquid Residual Saturation,Contact Angle,Mixed-Wet b Parameter,

Brooks and Corey^[2] Model

Rock/Soil Name,Brooks and Corey,ψ Parameter,Units (m),λ Parameter,Aqueous Residual Saturation,Nonaqueous-Liquid Residual Saturation,Contact Angle,Mixed-Wet b Parameter,

Brooks and Corey^[2] Model with the Webb^[3] Extension

Rock/Soil Name,Brooks and Corey with Webb,ψ Parameter,Units (m),λ Parameter,Aqueous Residual Saturation,Nonaqueous-Liquid Residual Saturation,Contact Angle,Mixed-Wet b Parameter,

 

References

[1] van Genuchten, MTh. 1980. “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Science Society of America Journal, 44:892-898.

[2] Brooks, RH and AT Corey. 1964. “Hydraulic properties of porous media,” Hydraulic Paper No. 3, Colorado State University, 1964.

[3] Webb, SW, 2000. “A simple extension of two-phase characteristic curves to include the dry region.” Water Resources Research, 36(6):1425-1430.