AbstractTwo-phase flow in fractured rock is an important phenomenon for a range of practical problems, not the least of which is non-aqueous phase liquid (NAPL) contamination of groundwater. Although multiphase systems have long been studied in the petroleum field, in the hydrogeological field progress has only just reaches the point where models are being developed. Scale effect is one of the main issues of concern. Although models presented in this paper have the potential to provide useful predictions, they can only be used to investigate a variety of possible scenarios with parameters being specified in the form of distribution of values. The calibration and validation of all but the simplest of these models poses a formidable task, with great demands on hydrogeologists and geophysicists to provide adequate data.