Three-dimensional upscaling of fault damage zones for reservoir simulation