Trapping and hysteresis in two‐phase flow in porous media: A pore‐network study

[1] Several models for two-phase flow in porous media identify trapping and connectivity of fluids as an important contribution to macroscale hysteresis. This is especially true for hysteresis in relative permeabilities. The trapping models propose trajectories from the initial saturation to the end saturation in various ways and are often based on experiments or pore-network model results for the endpoints. However, experimental data or pore-scale model results are often not available for the trajectories, that is, the fate of the connectivity of the fluids while saturation changes. Here, using a quasi static pore-network model, supported by a set of pore-scale laboratory experiments, we study how the topology of the fluids changes during drainage and imbibition including first, main and scanning curves. We find a strong hysteretic behavior in the relationship between disconnected nonwetting fluid saturation and the wetting fluid saturation in a water-wet medium. The coalescence of the invading nonwetting phase with the existing disconnected nonwetting phase depends critically on the presence (or lack thereof) of connected nonwetting phase at the beginning of the drainage process as well as on the pore geometry. This dependence involves a mechanism we refer to as ‘‘reversible corner filling.’’ This mechanism can also be seen in laboratory experiments in volcanic tuff. The impact of these pore-network model results on existing macroscopic models is discussed.

[1]  Ruben Juanes,et al.  Relative-Permeability Hysteresis: Trapping Models and Application to Geological CO2 Sequestration , 2005 .

[2]  G. R. Jerauld,et al.  The effect of pore-structure on hysteresis in relative permeability and capillary pressure: Pore-level modeling , 1990 .

[3]  Jack C. Parker,et al.  A parametric model for constitutive properties governing multiphase flow in porous media , 1987 .

[4]  Vahid Joekar-Niasar,et al.  Specific interfacial area: The missing state variable in two‐phase flow equations? , 2011 .

[5]  Vahid Joekar-Niasar,et al.  Network model investigation of interfacial area, capillary pressure and saturation relationships in granular porous media , 2010 .

[6]  Carlon S. Land,et al.  Calculation of Imbibition Relative Permeability for Two- and Three-Phase Flow From Rock Properties , 1968 .

[7]  Vahid Joekar-Niasar,et al.  Non-equilibrium effects in capillarity and interfacial area in two-phase flow: dynamic pore-network modelling , 2010, Journal of Fluid Mechanics.

[8]  Rudolf Hilfer,et al.  Capillary pressure, hysteresis and residual saturation in porous media , 2006 .

[9]  M. Blunt Flow in porous media — pore-network models and multiphase flow , 2001 .

[10]  Y. Mualem,et al.  A conceptual model of hysteresis , 1974 .

[11]  Jack C. Parker,et al.  A model for hysteretic constitutive relations governing multiphase flow: 2. Permeability‐saturation relations , 1987 .

[12]  Jason I. Gerhard,et al.  Measurement and prediction of the relationship between capillary pressure, saturation, and interfacial area in a NAPL‐water‐glass bead system , 2010 .

[13]  R. Hilfer,et al.  Numerical solutions of a generalized theory for macroscopic capillarity. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  E. M. Braun,et al.  Relative permeability hysteresis: Laboratory measurements and a conceptual model , 1995 .

[15]  Jack C. Parker,et al.  A model for hysteretic constitutive relations governing multiphase flow: 1. Saturation-pressure relations , 1987 .

[16]  William G. Gray,et al.  Thermodynamic basis of capillary pressure in porous media , 1993 .

[17]  Roger L. Nielsen,et al.  TRACE. FOR: a program for the calculation of combined major and trace-element liquid lines of descent for natural magmatic systems , 1988 .

[18]  S. M. Hassanizadeh,et al.  A model for two‐phase flow in porous media including fluid‐fluid interfacial area , 2008 .

[19]  G. Jerauld General Three-Phase Relative Permeability Model for Prudhoe Bay , 1997 .

[20]  Florian Doster,et al.  Generalized Buckley–Leverett theory for two-phase flow in porous media , 2011 .

[21]  Florian Doster,et al.  Percolation as a Basic Concept for Macroscopic Capillarity , 2010 .

[22]  Tadeusz W Patzek,et al.  Verification of a Complete Pore Network Simulator of Drainage and Imbibition , 2001 .

[23]  Nicholas Kouwen,et al.  Hysteretic effects on net infiltration , 1983 .

[24]  J. Parker,et al.  Comparing Simulated and Experimental Hysteretic Two‐Phase Transient Fluid Flow Phenomena , 1991 .

[25]  V. Duijn,et al.  A mathematical model for hysteretic two-phase flow in porous media , 2001 .

[26]  S. M. Hassanizadeh,et al.  Insights into the Relationships Among Capillary Pressure, Saturation, Interfacial Area and Relative Permeability Using Pore-Network Modeling , 2008 .

[27]  G. Miller,et al.  Rapid Pseudo-Static Measurement of Hysteretic Capillary Pressure-Saturation Relationships in Unconsolidated Porous Media , 2007 .

[28]  D. H. Everett,et al.  A general approach to hysteresis , 1952 .

[29]  M. L. Porter,et al.  Image analysis algorithms for estimating porous media multiphase flow variables from computed microtomography data: a validation study , 2010 .

[30]  S. M. Hassanizadeh,et al.  Analysis of Fundamentals of Two-Phase Flow in Porous Media Using Dynamic Pore-Network Models: A Review , 2012 .