A New Method for Accurate and Practical Upscaling In Highly Heterogeneous Reservoir Models