Pore‐scale analysis of trapped immiscible fluid structures and fluid interfacial areas in oil‐wet and water‐wet bead packs

The objective of this study is to obtain quantitative evidence of pore-scale immiscible fluid distribution in oil-wet and water-wet porous media using X-ray computed microtomography. Temporal and spatial saturation profiles, as well as surface and interfacial areas, are thoroughly analyzed through cycles of drainage and imbibition using samples with different wetting characteristics but similar pore structures. The population of individual immiscible fluid structures (‘blobs’) was also evaluated. The specific nonwetting phase surface areas of both porous media are found to be in close correlation with the specific solid surface area. On the other hand, the differing wetting strengths of the two porous media affect the curvature of the fluid–fluid interface and thus the specific meniscus interfacial area of the two porous media. Although the magnitude of the specific meniscus interfacial areas is different, they both trend toward a maximum at wetting phase saturations of 0.35–0.55. The differences in wetting characteristics are also apparent in the blob populations. The number of blobs in the oil-wet porous media is three times greater than that of the water-wet porous media at similar saturations; the increase in population is a result of the increase in the amount of smaller blobs inhabiting the smaller pore spaces. The surface areas of individual blobs as a function of the individual blob volumes are found to closely agree with the specific surface area of a sphere at blob volumes below the minimum individual grain volume and with the specific pore space surface area above this volume. These results show how wettability and saturation history influence the distribution of immiscible fluids within the pore space.

[1]  F. Leij,et al.  Estimating interfacial areas for multi-fluid soil systems , 1997 .

[2]  Markus Hilpert,et al.  Computation of the interfacial area for two-fluid porous medium systems. , 2002, Journal of contaminant hydrology.

[3]  M. Brusseau,et al.  Synchrotron X‐ray microtomography and interfacial partitioning tracer test measurements of NAPL‐water interfacial areas , 2008, Water resources research.

[4]  Nicholas J. Giordano,et al.  Measurement of interfacial area per volume for drainage and imbibition , 2007 .

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

[6]  Michael A. Celia,et al.  Modeling support of functional relationships between capillary pressure, saturation, interfacial area and common lines , 2001 .

[7]  Gurpreet Singh,et al.  Experimental investigation of trapped oil clusters in a water‐wet bead pack using X‐ray microtomography , 2010 .

[8]  Mohammad Piri,et al.  Pore-scale modeling of dispersion in disordered porous media. , 2011, Journal of contaminant hydrology.

[9]  Michael A. Celia,et al.  A Functional Relationship Between Capillary Pressure, Saturation, and Interfacial Area as Revealed by a Pore‐Scale Network Model , 1996 .

[10]  Kenneth Stuart Sorbie,et al.  Two-Phase Pore-Network Modelling: Existence of Oil Layers During Water Invasion , 2009 .

[11]  M. Piri,et al.  The effects of wettability and trapping on relationships between interfacial area, capillary pressure and saturation in porous media: A pore-scale network modeling approach , 2009 .

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

[13]  I. Chatzis,et al.  Dissolution of residual non-aqueous phase liquids in porous media : Pore-scale mechanisms and mass transfer rates , 2002 .

[14]  R. Al-Raoush,et al.  A pore-scale investigation of a multiphase porous media system. , 2005, Journal of contaminant hydrology.

[15]  William G. Gray,et al.  Pore-scale characteristics of multiphase flow in porous media: A comparison of air–water and oil–water experiments , 2006 .

[16]  Kenneth Stuart Sorbie,et al.  Prediction of Three-Phase Relative Permeabilities Using a Pore-Scale Network Model Anchored to Two-Phase Data , 2007 .

[17]  M. Blunt,et al.  Pore space reconstruction of vuggy carbonates using microtomography and multiple‐point statistics , 2007 .

[18]  R. Al-Raoush.  Impact of wettability on pore-scale characteristics of residual nonaqueous phase liquids. , 2009, Environmental science & technology.

[19]  Mohammad Piri,et al.  Prediction of fluid occupancy in fractures using network modeling and x-ray microtomography. II: results. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  M. Honarpour,et al.  Capillary end effects in coreflood calculations , 1998 .

[21]  Kenneth Stuart Sorbie,et al.  3D Stochastic Modelling of Heterogeneous Porous Media – Applications to Reservoir Rocks , 2006 .

[22]  M. Brusseau,et al.  Relationships among air‐water interfacial area, capillary pressure, and water saturation for a sandy porous medium , 2005 .

[23]  W. E. Soll,et al.  Pore level imaging of fluid transport using synchrotron X-ray microtomography , 1996 .

[24]  Mohammad Piri,et al.  Direct pore-level modeling of incompressible fluid flow in porous media , 2010, J. Comput. Phys..

[25]  Mohammad Piri,et al.  Prediction of fluid occupancy in fractures using network modeling and x-ray microtomography. I: data conditioning and model description. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Stanley R. Rotman,et al.  Comments on Picture thresholding using an iterative selection method , 1990, IEEE Trans. Syst. Man Cybern..

[27]  Steven Robert McDougall,et al.  Anchoring methodologies for pore-scale network models: Application to relative permeability and capillary pressure prediction , 2002 .

[28]  S. Bryant,et al.  A level set method for determining critical curvatures for drainage and imbibition. , 2006, Journal of colloid and interface science.

[29]  Z. Karpyn,et al.  A Study of Hydraulic Fracture Conductivity and Its Dependence on Proppant Wettability , 2010 .

[30]  T. Austad,et al.  “Smart water” as a wettability modifier in chalk: The effect of salinity and ionic composition , 2010 .

[31]  L. Girifalco,et al.  A Theory for the Estimation of Surface and Interfacial Energies. I. Derivation and Application to Interfacial Tension , 1957 .

[32]  K. Sorbie,et al.  Efficient flow and transport simulations in reconstructed 3D pore geometries , 2010 .

[33]  Edo S. Boek,et al.  Lattice-Boltzmann studies of fluid flow in porous media with realistic rock geometries , 2010, Comput. Math. Appl..

[34]  Clinton S. Willson,et al.  Comparison of Network Generation Techniques for Unconsolidated Porous Media , 2003 .

[35]  Ioannis Chatzis,et al.  Entrapment and Mobilization of Residual Oil in Bead Packs , 1988 .

[36]  S. L. Wellington,et al.  Tomographic imaging of three‐phase flow experiments , 1987 .

[37]  P. Meakin,et al.  Pore scale modeling of immiscible and miscible fluid flows using smoothed particle hydrodynamics , 2006 .

[38]  Mohammad Piri,et al.  Three-dimensional mixed-wet random pore-scale network modeling of two- and three-phase flow in porous media. II. Results. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[39]  Nicos Martys,et al.  Transport in sandstone: A study based on three dimensional microtomography , 1996 .

[40]  S. Bryant,et al.  Influence of wettability and saturation on liquid-liquid interfacial area in porous media. , 2003, Environmental science & technology.

[41]  Marcel G. Schaap,et al.  Comparison of pressure‐saturation characteristics derived from computed tomography and lattice Boltzmann simulations , 2007 .

[42]  S. Bryant,et al.  BULK AND FILM CONTRIBUTIONS TO FLUID/FLUID INTERFACIAL AREA IN GRANULAR MEDIA , 2004 .

[43]  M. L. Porter,et al.  Lattice-Boltzmann simulations of the capillary pressure–saturation–interfacial area relationship for porous media , 2009 .

[44]  M. Costanza-Robinson,et al.  X-ray microtomography determination of air-water interfacial area-water saturation relationships in sandy porous media. , 2008, Environmental science & technology.

[45]  William G. Gray,et al.  Interfacial area measurements for unsaturated flow through a porous medium , 2004 .

[46]  M. Blunt,et al.  A functional relation for field-scale nonaqueous phase liquid dissolution developed using a pore network model. , 2001, Journal of contaminant hydrology.

[47]  Mohammad Piri,et al.  Three-dimensional mixed-wet random pore-scale network modeling of two- and three-phase flow in porous media. I. Model description. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  J. Gerhard,et al.  Simulating the dissolution of a complex dense nonaqueous phase liquid source zone: 1. Model to predict interfacial area , 2007 .

[49]  Kenneth Stuart Sorbie,et al.  Efficient extraction of networks from three‐dimensional porous media , 2007 .

[50]  M. Brusseau,et al.  Measurement and estimation of organic-liquid/water interfacial areas for several natural porous media. , 2009, Environmental science & technology.

[51]  N. Wardlaw,et al.  Oil blob populations and mobilization of trapped oil in unconsolidated packs , 1985 .

[52]  J. Kearney,et al.  Functional Relationship Between , 1991 .

[53]  C. Lévi-Strauss,et al.  Experimental investigation , 2013 .