Abstract:We have developed a new approach for generating pore throat cross sections of various shapes based on distributions of shape factors and radii of inscribed circles. These distributions are obtained from analysis of grains packing. General formulas for calculating geometrical properties and entry capillary pressure for given shape factor and inscribed circle radius are developed. These relationships are employed in a pore network, which has a number of special features. In particular, it is highly flexible in terms of location of pore bodies, variable coordination number, as well as variable cross‐sectional shapes. The pore network model is employed for simulating the equilibrium distribution of two fluids in a granular porous medium, under both drainage and imbibition conditions. The pore network model is verified by comparing simulation results with experimental data of quasi‐static drainage and imbibition experiments in a glass bead medium. The pore‐level topology and geometrical description of pore bodies and pore throats, essential for building the network, are rigorously extracted from experimental data using image analysis (3DMA‐Rock software). Calculated capillary pressure‐saturation (Pc − Sw) and specific interfacial area‐saturation (anw − Sw) curves show very good agreement with measured ones, for both drainage and imbibition. We show that the shape factor can significantly influence the form of macroscopic Pc − Sw and anw − Sw curves, if the length and volumes associated to the pore throats are considerable. Furthermore, using continuous generation of shape factor distribution, the model can be validated against the grain size distribution. After validating the model against experiments, in addition to primary and main curves, we simulate many scanning curves to generate Pc − Sw − anw surfaces for drainage and imbibition, separately. Results show that these two surfaces lie very close to each other, and the average normalized difference is small, in the range of simulations uncertainty. Our results illustrate that Pc − Sw − anw surfaces show very little hysteresis and, therefore, specific interfacial area can be considered as an essential variable for reducing or eliminating the hysteresis observed in Pc − Sw curves.

Generalized network modeling of capillary-dominated two-phase flow

Pore-scale modeling of two-phase flow: A comparison of the generalized network model to direct numerical simulation

Validating the Generalized Pore Network Model Using Micro-CT Images of Two-Phase Flow

Three-dimensional mixed-wet random pore-scale network modeling of two- and three-phase flow in porous media. I. Model description

Validation of Model Predictions of Pore-Scale Fluid Distributions During Two-Phase Flow

Pore network modeling of the effects of viscosity ratio and pressure gradient on steady-state incompressible two-phase flow in porous media

Pore-scale Network Modeling of Three-Phase Flow Based on Thermodynamically Consistent Threshold Capillary Pressures. II. Results

Modeling support of functional relationships between capillary pressure, saturation, interfacial area and common lines

Pore network models to determine flow statistics and structural controls in variably saturated porous media

Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method

Capillary Network Model: Capillary Power and Effective Permeability

Network model investigation of interfacial area, capillary pressure and saturation relationships in granular porous media

The Brooks and Corey Capillary Pressure Model Revisited from Pore Network Simulations of Capillarity-Controlled Invasion Percolation Process

A Two-Dimensional Network Simulator for Two-Phase Flow in Porous Media

An Improved Pore-Network Model Including Viscous Coupling Effects Using Direct Simulation by the Lattice Boltzmann Method

Modeling flow in porous media with double porosity/permeability: Mathematical model, properties, and analytical solutions

Pore-to-mesoscale Network Modeling of Heat Transfer and Fluid Flow in Packed Beds with Application to Process Design

Fluid Meniscus Algorithms for Dynamic Pore-Network Modeling of Immiscible Two-Phase Flow in Porous Media

Pore-Scale Modeling of Non-Newtonian Flow in Porous Media

APPLICATION OF THE PRECONDITIONED CONJUGATE GRADIENTS METHOD IN THE SIMULATION OF RELATIVE PERMEABILITY PROPERTIES OF POROUS MEDIA

Pore-by-pore modeling, analysis, and prediction of two-phase flow in mixed-wet rocks