Soil foundations exhibit significant creeping deformation, which may result in excessive settlement and failure of superstructures. Based on the theory of viscoelasticity and fractional calculus, a fractional Kelvin-Voigt model is proposed to account for the time-dependent behavior of soil foundation under vertical line load. Analytical solution of settlements in the foundation was derived using Laplace transforms. The influence of the model parameters on the time-dependent settlement is studied through a parametric study. Results indicate that the settlement-time relationship can be accurately captured by varying values of the fractional order of differential operator and the coefficient of viscosity. In comparison with the classical Kelvin-Voigt model, the fractional model can provide a more accurate prediction of long-term settlements of soil foundation. The determination of influential distance also affects the calculation of settlements.