Two-phase subsurface flows with capillary heterogeneity are studied.
A new method is developed to upscale relative permeability and capillary pressure.
The procedure is tested on a numerical model of CO2 injection into an aquifer.
The method provides accurate results and clearly improves upon simpler treatments.
Upscaled properties are sufficiently robust for reuse in related simulations.