Rheology of Immiscible Two-phase Flow in Mixed Wet Porous Media: Dynamic Pore Network Model and Capillary Fiber Bundle Model Results

Immiscible two-phase flow in porous media with mixed wet conditions was examined using a capillary fiber bundle model, which is analytically solvable, and a dynamic pore network model. The mixed wettability was implemented in the models by allowing each tube or link to have a different wetting angle chosen randomly from a given distribution. Both models showed that mixed wettability can have significant influence on the rheology in terms of the dependence of the global volumetric flow rate on the global pressure drop. In the capillary fiber bundle model, for small pressure drops when only a small fraction of the tubes were open, it was found that the volumetric flow rate depended on the excess pressure drop as a power law with an exponent equal to 3/2 or 2 depending on the minimum pressure drop necessary for flow. When all the tubes were open due to a high pressure drop, the volumetric flow rate depended linearly on the pressure drop, independent of the wettability. In the transition region in between where most of the tubes opened, the volumetric flow depended more sensitively on the wetting angle distribution function and was in general not a simple power law. The dynamic pore network model results also showed a linear dependence of the flow rate on the pressure drop when the pressure drop is large. However, out of this limit the dynamic pore network model demonstrated a more complicated behavior that depended on the mixed wettability condition and the saturation. In particular, the exponent relating volumetric flow rate to the excess pressure drop could take on values anywhere between 1.0 and 1.8. The values of the exponent were highest for saturations approaching 0.5, also, the exponent generally increased when the difference in wettability of the two fluids were larger and when this difference was present for a larger fraction of the porous network.

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

We present in detail a set of algorithms for a dynamic pore-network model of immiscible two-phase flow in porous media to carry out fluid displacements in pores. The algorithms are universal for regular and irregular pore networks in two or three dimensions and can be applied to simulate both drainage displacements and steady-state flow. They execute the mixing of incoming fluids at the network nodes, then distribute them to the outgoing links and perform the coalescence of bubbles. Implementing these algorithms in a dynamic pore-network model, we reproduce some of the fundamental results of transient and steady-state two-phase flow in porous media. For drainage displacements, we show that the model can reproduce the flow patterns corresponding to viscous fingering, capillary fingering and stable displacement by varying the capillary number and viscosity ratio. For steady-state flow, we verify non-linear rheological properties and transition to linear Darcy behavior while increasing the flow rate. Finally we verify the relations between seepage velocities of two-phase flow in porous media considering both disordered regular networks and irregular networks reconstructed from real samples.