In tunnel construction in western China, a vast amount of carbonaceous slate is encountered. High in situ stress and foliation structures cause the rock mass to exhibit pronounced anisotropic creep, readily inducing a series of engineering disasters like collapses and lining cracks. Investigating the anisotropic time-dependent characteristics of carbonaceous slate is beneficial to the long-term stability of tunnel construction and operation. In view of this, carbonaceous slate specimens with different angles, beta, between the foliation plane and loading direction were studied using a graded loading method through uniaxial compression creep tests. The results show that the instantaneous axial strain, epsilon(i), the axial creep strain, epsilon(c), the duration time of decelerating creep stage, t(d), and the steady creep strain rate, epsilon(center dot)(s), increased with the rise in the loading ratio, k. Their variations followed a power law relationship, with the R-2 (Coefficient of Determination) values all exceeding 0.95. The value of epsilon(center dot)(s) was observed to be less than 1.5 x 10(-4)/h when beta < 45 degrees, while it was found to exceed 1.5 x 10(-4)/h in the cases of beta >= 45 degrees. The long-term strength, sigma(L), of carbonaceous slate showed a U-shaped pattern with the variation in beta. The maximum sigma(L) occurred at beta = 90 degrees and the minimum was observed at beta = 15 degrees. A fractional nonlinear creep model (FNC model) was developed. The sensitivity analysis reveals that the larger the fractional order n is, the t(d) and epsilon(center dot)(s) increase. eta(2) and E-2 primarily affect the decelerated creep stage, while the epsilon(center dot)(s) exhibits a rapid increase with the rise of eta(1). To further validate the FNC model, a comparison is made with the traditional Nishihara model. The R-2 of the FNC model is larger than 0.965, which is higher than that of the Nishihara model (R-2 <= 0.911). The FNC model can effectively cope with the impact of the sudden increase in strain and well describe the characteristics of the decelerating, steady-state, and accelerating creep stages at any stress level and any angle. The results provide a reference for the study of the creep mechanism of layered rocks.