Simulating drainage and imbibition experiments in a high‐porosity micromodel using an unstructured pore network model

Development of pore network models based on detailed topological data of the pore space is essential for predicting multiphase flow in porous media. In this work, an unstructured pore network model has been developed to simulate a set of drainage and imbibition laboratory experiments performed on a two-dimensional micromodel. We used a pixel-based distance transform to determine medial pixels of the void domain of micromodel. This process provides an assembly of medial pixels with assigned local widths that simulates the topology of the porous medium. Using this pore network model, the capillary pressure-saturation and capillary pressure-interfacial area curves measured in the laboratory under static conditions were simulated. On the basis of several imbibition cycles, a surface of capillary pressure, saturation and interfacial area was produced. The pore network model was able to reproduce the distribution of the fluids as observed in the micromodel experiments. We have shown the utility of this simple pore network approach for capturing the topology and geometry of the micromodel pore structure.

[1]  Christoph H. Arns,et al.  Digital Core Laboratory: Properties of reservoir core derived from 3D images , 2004 .

[2]  Fluid topology, pore size and aspect ratio during imbibition , 1988 .

[3]  S. K. Jenson,et al.  Extracting topographic structure from digital elevation data for geographic information-system analysis , 1988 .

[4]  D. Denney Digital Core Laboratory: Reservoir-Core Properties Derived From 3D Images , 2004 .

[5]  M. Brady,et al.  Smoothed Local Symmetries and Their Implementation , 1984 .

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

[7]  Hans-Jörg Vogel,et al.  A new approach for determining effective soil hydraulic functions , 1998 .

[8]  Philippe Saint-Marc,et al.  B-spline Contour Representation and Symmetry Detection , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[9]  C. Montemagno,et al.  Porosity of natural fracture networks , 1995 .

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

[11]  Hans-Jörg Vogel,et al.  Morphological determination of pore connectivity as a function of pore size using serial sections , 1997 .

[12]  R. Glantz,et al.  Dual models of pore spaces , 2007 .

[13]  Shouxiang Ma,et al.  Effect of contact angle on drainage and imbibition in regular polygonal tubes , 1996 .

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

[15]  Ching Y. Suen,et al.  Thinning Methodologies - A Comprehensive Survey , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[16]  Pore-level modelling of wetting on correlated porous media , 1998 .

[17]  W. E. Soll,et al.  Developments in Synchrotron X-Ray Microtomography with Applications to Flow in Porous Media , 1996 .

[18]  I. Chatzis,et al.  A mercury porosimeter for investigating capillary phenomena and microdisplacement mechanisms in capillary networks , 1991 .

[19]  R. P. Mayer,et al.  Mercury porosimetry—breakthrough pressure for penetration between packed spheres , 1965 .

[20]  D. T. Lee,et al.  Medial Axis Transformation of a Planar Shape , 1982, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[21]  Q. Ye The signed Euclidean distance transform and its applications , 1988, [1988 Proceedings] 9th International Conference on Pattern Recognition.

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

[23]  Martin J. Blunt,et al.  Predictive pore‐scale modeling of two‐phase flow in mixed wet media , 2004 .

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

[25]  D. Wildenschilda,et al.  Using X-ray computed tomography in hydrology : systems , resolutions , and limitations , 2002 .

[26]  Azriel Rosenfeld,et al.  Computation of geometric properties from the medial axis transform in (O(nlogn) time , 1986, Comput. Vis. Graph. Image Process..

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

[28]  Hans-Jörg Vogel,et al.  A numerical experiment on pore size, pore connectivity, water retention, permeability, and solute transport using network models , 2000 .

[29]  I. Fatt The Network Model of Porous Media , 1956 .

[30]  David M. Mount,et al.  Parallel computational geometry of rectangles , 1992, Algorithmica.

[31]  Markus Ilg,et al.  Voronoi skeletons: theory and applications , 1992, Proceedings 1992 IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[32]  Matthew D. Jackson,et al.  Detailed physics, predictive capabilities and macroscopic consequences for pore-network models of multiphase flow. , 2002 .

[33]  Azriel Rosenfeld,et al.  Parallel computation of geometric properties from the medial axis transform , 1988, Comput. Vis. Graph. Image Process..

[34]  D. Nolte,et al.  Relating capillary pressure to interfacial areas , 2008 .

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

[36]  David D. Nolte,et al.  Linking pressure and saturation through interfacial areas in porous media , 2004 .

[37]  Á. Rodríguez-Rey,et al.  Microfractography of granitic rocks under confocal scanning laser microscopy , 1995 .

[38]  E. Bazzo,et al.  Manufacturing of metallic porous structures to be used in capillary pumping systems , 2003 .

[39]  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.

[40]  W. B. Lindquist,et al.  Pore and throat size distributions measured from synchrotron X-ray tomographic images of Fontaineble , 2000 .

[41]  Raymond W. Smith,et al.  Computer processing of line images: A survey , 1987, Pattern Recognit..

[42]  F. Dullien,et al.  Mercury porosimetry curves of sandstones. Mechanisms of mercury penetration and withdrawal , 1981 .

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

[44]  Roland Lenormand,et al.  Role Of Roughness And Edges During Imbibition In Square Capillaries , 1984 .

[45]  Theodosios Pavlidis,et al.  Feature Analysis Using Line Sweep Thinning Algorithm , 1999, IEEE Trans. Pattern Anal. Mach. Intell..

[46]  H. M. Princen Capillary phenomena in assemblies of parallel cylinders: II. Capillary rise in systems with more than two cylinders , 1969 .

[47]  Roland Glantz,et al.  Tight dual models of pore spaces , 2008 .

[48]  Madalena M. Dias,et al.  Network models for two-phase flow in porous media Part 1. Immiscible microdisplacement of non-wetting fluids , 1986, Journal of Fluid Mechanics.

[49]  Ugo Montanari,et al.  Continuous Skeletons from Digitized Images , 1969, JACM.

[50]  W. Mahmud,et al.  Effects of Snap-Off in Imbibition in Porous Media with Different Spatial Correlations , 2006 .

[51]  J. T. Fredrich,et al.  Imaging the Pore Structure of Geomaterials , 1995, Science.

[52]  J. W. Martin,et al.  Pore geometry and transport properties of Fontainebleau sandstone , 1993 .

[53]  A. ROSENFELD,et al.  Distance functions on digital pictures , 1968, Pattern Recognit..

[54]  R. Lenormand,et al.  Mechanisms of the displacement of one fluid by another in a network of capillary ducts , 1983, Journal of Fluid Mechanics.

[55]  R. Al-Raoush,et al.  Extraction of physically realistic pore network properties from three-dimensional synchrotron X-ray microtomography images of unconsolidated porous media systems , 2005 .