Over the last two decades there has been an increase in activity on the pore-scale modeling of multiphase flow in porous media. Excellent progress has been made in many areas of pore-scale modeling, particularly in (1) the representation of the rock itself and (2) our description of the pore-scale displacement physics (in model pore geome-tries). Three-dimensional voxelized images of actual rocks can be generated either numerically (e.g. from 2D thin sections) or from micro-CT imaging. A simplified network involving more idealized nodes and bonds can then be extracted from this numerical rock model and this can be used in modeling pore-scale displacement processes. Much progress has also been made in understanding these pore-scale processes (i.e. piston-like displacement, snap-off events, layer formation/collapse, pore-body filling draining). These processes can be mathematically modeled accurately for pores of non uniform wettability, if the geometry of the pore is sufficiently simple. In fact, in recent years these various pore-level processes in mixed and fractionally wet systems have been classified as "events" in an entire capillary-dominated "phase space" which can be defined in a thermodynamically consistent manner. Advances in our understanding and ability to compute several two- (and three-) phase properties a priori have been impressive and the entire flooding cycle of primary drainage (PD), aging/wetting change, and imbibition can be simulated.
In this paper, we review the successes of pore-scale network modeling and explain how it can be of great use in understanding and explaining many phenomena in flow through porous media. However, we also critically examine the issue of how predictive network modeling is in practice. Indeed, one of our conclusions on pore-scale modeling in mixed-wet systems is that we cannot predict two-phase functions reliably in "blind" tests. Interestingly, we make this statement not because we do not understand the pore-scale physics of the process, but because we do understand the physics. It is hoped that our comments will stimulate a more critical debate on the role of pore-scale modeling and its use in core analysis.