Direct Quantification of Dynamic Effects in Capillary Pressure for Drainage–Wetting Cycles

The constitutive relationship between capillary pressure ( P c) and wetting fluid saturation ( S w), or retention curve, is needed to model multiphase flow in porous media. This relationship is usually measured under static conditions; however, transient flow is governed by a dynamic relationship between the P c and S w. Differences in P c measured under static and dynamic conditions are due to dynamic effects typically defined as a product of a dynamic coefficient (τ) and the rate of change in S w. To date, relatively few experimental studies have been conducted to directly quantify the magnitude of this effect. In this study, the magnitude of τ was quantified by measuring both static and dynamic retention curves in repeated drainage and wetting experiments using a field sand. The 95% confidence intervals for the static retention curves showed that the dynamic retention curves were statistically different. The measured τ for primary drainage generally increased with decreasing S w. The measured τ values were also compared with those estimated using a different approach based on redistribution time. The measured and estimated τ were in close agreement when the redistribution times were 146 s for the wetting cycle and 509 s for primary and main drainage cycles. The shape of the τ– S w relationship was largely controlled by the slope of the static retention curve. Numerical modeling demonstrated that a log-linear model relating τ and S w yielded the best match to experimental outflow results.

[1]  T. Phelan,et al.  Exploring dynamic effects in capillary pressure in multistep outflow experiments , 2005 .

[2]  Multi-step and two-step experiments in heterogeneous porous media to evaluate the relevance of dynamic effects , 2008 .

[3]  G. I. Barenblatt,et al.  The Mathematical Model of Nonequilibrium Effects in Water-Oil Displacement , 2003 .

[4]  Jack C. Parker,et al.  Determining Soil Hydraulic Properties from One-step Outflow Experiments by Parameter Estimation: I. Theory and Numerical Studies1 , 1985 .

[5]  Michael A. Celia,et al.  Dynamic Effect in the Capillary Pressure–Saturation Relationship and its Impacts on Unsaturated Flow , 2002 .

[6]  William G. Gray,et al.  Toward an improved description of the physics of two-phase flow , 1993 .

[7]  S. Hassanizadeh,et al.  Dynamic capillary pressure effects in two-phase flow through heterogeneous porous media , 2004 .

[8]  R. Hilfer,et al.  Dimensional analysis of two-phase flow including a rate-dependent capillary pressure-saturation relationship , 2008 .

[9]  William G. Gray,et al.  Mechanics and thermodynamics of multiphase flow in porous media including interphase boundaries , 1990 .

[10]  J. H. Dane,et al.  Capillary Pressure–Saturation Relations for Saprolite: Scaling With and Without Correction for Column Height , 2004 .

[11]  Experimental quantification of bulk sampling volume of ECH2O soil moisture sensors , 2009 .

[12]  A. Klute,et al.  TENSIOMETER RESPONSE TIME , 1962 .

[13]  Ole Wendroth,et al.  Unsaturated hydraulic conductivity from transient multistep outflow and soil water pressure data , 1994 .

[14]  N. T. Burdine Relative Permeability Calculations From Pore Size Distribution Data , 1953 .

[15]  Olaf Ippisch,et al.  Dynamic nonequilibrium in unsaturated water flow , 1999 .

[16]  Peter Droogers,et al.  Inverse method to determine soil hydraulic functions from multistep outflow experiments. , 1994 .

[17]  T. Patzek,et al.  On Barenblatt's Model of Spontaneous Countercurrent Imbibition , 2004 .

[18]  F.J-M. Kalaydjian,et al.  Dynamic Capillary Pressure Curve for Water/Oil Displacement in Porous Media: Theory vs. Experiment , 1992 .

[19]  P. A. Gorry General least-squares smoothing and differentiation by the convolution (Savitzky-Golay) method , 1990 .

[20]  J. Schembre,et al.  Estimation of Dynamic Relative Permeability and Capillary Pressure from Countercurrent Imbibition Experiments , 2006 .

[21]  Jack C. Parker,et al.  Determining Soil Hydraulic Properties from One‐step Outflow Experiments by Parameter Estimation: II. Experimental Studies , 1985 .

[22]  John S. Selker,et al.  Characterization of Miller‐Similar Silica Sands for Laboratory Hydrologic Studies , 1996 .

[23]  M. Celia,et al.  Bundle-of-Tubes Model for Calculating Dynamic Effects in the Capillary-Pressure- Saturation Relationship , 2005 .

[24]  D. Zhou,et al.  An Investigation of Counter-Current Imbibition Processes in Diatomite , 2001 .

[25]  Anthony R. Kovscek,et al.  Scaling of counter-current imbibition processes in low-permeability porous media , 2002 .

[26]  Jan W. Hopmans,et al.  Optimization of Hydraulic Functions from Transient Outflow and Soil Water Pressure Data , 1993 .

[27]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[28]  R. H. Brooks,et al.  Hydraulic properties of porous media , 1963 .

[29]  Susan E. Powers,et al.  Estimating unique soil hydraulic parameters for sandy media from multi-step outflow experiments , 2003 .

[30]  J. H. Dane,et al.  Capillary pressure-saturation relations for saprolite: scaling with and without correction for column height , 2004 .

[31]  Tissa H. Illangasekare,et al.  Comparison of height‐averaged and point‐measured capillary pressure–saturation relations for sands using a modified Tempe cell , 2007 .

[32]  Ruben Juanes,et al.  Nonequilibrium effects in models of three-phase flow in porous media , 2008 .

[33]  Malgorzata Peszynska,et al.  Numerical methods for unsaturated flow with dynamic capillary pressure in heterogeneous porous media , 2008 .

[34]  A. Kovscek,et al.  Nonequilibrium Effects During Spontaneous Imbibition , 2006 .

[35]  Mark E. Grismer,et al.  Parameter estimation of two-fluid capillary pressure–saturation and permeability functions , 1999 .