Upscaling polymer flooding to model sub-grid block geological heterogeneity and compensate for numerical dispersion