The permeability of the Earth's crust commonly varies over many orders of magnitude. Flow velocity can range over several orders of magnitude in structures of interest that vary in scale from centimeters to kilometers. To accurately and efficiently model multiphase flow in geologic media, we introduce a fully conservative node-centered finite volume method coupled with a Galerkin finite element method on an unstructured triangular grid with a complementary finite volume subgrid. The effectiveness of this approach is demonstrated by comparison with traditional solution methods and by multiphase flow simulations for heterogeneous permeability fields including complex geometries that produce transport parameters and lengths scales varying over four orders of magnitude.