Anisotropic rate-dependent saturation functions for compositional simulation of sandstone composites