The impact of rheological deformation in soft soils on engineering construction is a critical consideration. To gain a comprehensive understanding of the variations in soft soils during the consolidation rheological process, the non-Darcy seepage model proposed by Hansbo was introduced, in conjunction with the fractional-order Merchant model, to rederive the one-dimensional consolidation equation for saturated clay. Subsequently, the finite difference method was employed to solve this equation, and a computer program was developed for implementation. Finally, the influence of different mechanical parameters on ground settlement and pore pressure dissipation was examined. The results indicate that the consolidation fractional derivative order, denoted as α, leads to intersections in both the settlement and pore pressure curves. A larger viscosity coefficient, represented as F, results in an extended duration for settlement.