Harnack estimates for non-negative weak solutions of a class of singular parabolic equations