MONTE-CARLO SIMULATIONS OF MULTIPHASE FLOW INCORPORATING SPATIAL VARIABILITY OF HYDRAULIC-PROPERTIES

To study the effect of spatial variability of sediment hydraulic properties on multiphase flow, oil infiltration into a hypothetical glacial outwash aquifer, followed by oil extraction, was simulated using a cross-sectional multiphase flow model. The analysis was simplified by neglecting capillary hysteresis. The first simulation used a uniform mean permeability and mean retention curve. This was followed by 50 Monte Carlo simulations conducted using 50 spatially variable permeability realizations and corresponding spatially variable retention curves. For the type of correlation structure considered in this study, which is similar to that of glacial outwash deposits, use of mean hydraulic properties reproduces the ensemble average oil saturation distribution obtained from the Monte Carlo simulations. However, spatial variability the oil saturation distribution in an individual oil lens to differ significantly from that of the mean lens. Oil saturations at a given location may be considerably higher than would be predicted using uniform mean properties. During cleanup by oil extraction from a well, considerably more oil may remain behind in the heterogeneous case than in the spatially uniform case.