We compute the fluid flow time-correlation functions of incompressible, immiscible two-phase flow in porous media using a 2D network model. Given a properly chosen representative elementary volume, the flow rate distributions are Gaussian and the integrals of time correlation functions of the flows are found to converge to a finite value. The integrated cross-correlations become symmetric, obeying Onsager's reciprocal relations. These findings support the proposal of a non-equilibrium thermodynamic description for two-phase flow in porous media.

We perform more than 6000 steady-state simulations with a dynamic pore network model, corresponding to a large span in viscosity ratios and capillary numbers. From these simulations, dimensionless quantities such as relative permeabilities, residual saturations, mobility ratios and fractional flows are computed. Relative permeabilities and residual saturations show many of the same qualitative features observed in other experimental and modeling studies. However, while other studies find that relative permeabilities converge to straight lines at high capillary numbers we find that this is not the case when viscosity ratios are different from 1. Our conclusion is that departure from straight lines occurs when fluids mix rather than form decoupled flow channels. Another consequence of the mixing is that computed fractional flow curves, plotted against saturation, lie closer to the diagonal than they would otherwise do. At lower capillary numbers, fractional flow curves have a classical S-shape. Ratios of average mobility to their high-capillary number limit values are also considered. These vary, roughly, between 0 and 1, although values larger than 1 are also observed. For a given saturation and viscosity ratio, the mobilities are not always monotonically increasing with the pressure gradient. While increasing the pressure gradient mobilizes more fluid and activates more flow paths, when the mobilized fluid is more viscous, a reduction in average mobility may occur.

Immiscible fluids flowing at high capillary numbers in porous media may be characterized by an effective viscosity. We demonstrate that the effective viscosity is well described by the Lichtenecker-Rother equation. The exponent α in this equation takes either the value 1 or 0.6 in two-and 0.5 in three-dimensional systems depending on the pore geometry. Our arguments are based on analytical and numerical methods.PACS numbers:

In the original publication of the article, the Electronic Supplementary Material was missed.