A new matrix propagation method for theoretical analysis of wave propagation in low-temperature gradient layered rock mass with parallel ice-filled joints