Intrinsic Harnack estimates for non-negative local solutions of degenerate parabolic equations