Efficient approximation of the exponential operator for discrete 2D advection-diffusion problems