Method for upscaling of CO2 migration in 3D heterogeneous geological models