Abstract: In order to study the laws of surface subsidence after filling, we establish a mathematical model about fractional order creep and seepage. This model is solved by multi-physics numerical simulation software COMSOL and compared with classical creep-seepage model and field data. The calculation results show that the classical seepage coupling model has big difference with actual measured data with the increase of time, while the results calculated by fractional model are more close to the actual measured data.