Statistics of Lens galaxies in the Illustris simulation: power-law models and the bias of the Hubble constant from time delays