Abstract To explore how the microscopic geometry of a pore space affects the macroscopic characteristics of fluid flow in porous media, we have used approximate solutions of the Navier-Stokes equations to calculate the flow of two fluids in random networks. The model pore space consists of an array of pores of variable radius connected to a random number of nearest neighbors by throats of variable length and radius. The various size and connectedness distributions may be arbitrarily assigned, as are the wetting characteristics of the two fluids in the pore space. The fluids are assumed to be incompressible, immiscible, Newtonian, and of equal viscosity. In the calculation, we use Stokes flow results for the motion of the individual fluids and incorporate microscopic capillary force by using the Washburn approximation. At any time, the problem is mathematically identical to a random electrical network of resistors, batteries, and diodes. From the numerical solution of the latter, we compute the fluid velocities and saturation rates of change and use a discrete timestepping procedure to follow the subsequent motion. The scale of the computation has restricted us so far to networks of hundreds of pores in two dimensions (2D). Within these limitations, we discuss the dependence of residual oil saturations and interface shapes on network geometry and flow conditions.

We discuss the low Reynolds number flow of one or two immiscible Newtonian fluids in network models of microscopically random porous media. For the case of a single fluid, we reduce the flow problem to an analog random electrical resistor problem and use an "effective medium theory" to express the permeability of such networks in terms of the pore space geometry. For the flow of two fluids we use the Washburn approximation to incorporate capillary pressure differences, and show that this problem may also be formulated as a random electrical network. In this case, the capillary menisci correspond to moving batteries, and we follow the motion of the fluid-fluid interface (the ensemble of analog batteries) by a time-step procedure. We study the time evolution of the interface and the dynamics of blobs of one fluid contained in the other. as a function of the network geometry.