<p>The thermal consolidation theory of soils is a prominent research topic in modern geotechnical engineering. Thermal consolidation involves coupled interaction of thermal, hydraulic, and mechanical fields within the soil matrix. Among these interactions, the temperature-dependent permeability of saturated soils plays a crucial role in controlling water flow. However, all the existing models predominantly based on Darcy’s law often neglect the impact of non-Darcian flow and nonlinear compression on the thermal consolidation. In this context, a non-Darcian law of Hansbo’s flow is introduced into the one-dimensional thermal consolidation of soft soil for the first time in this study. In this way, the governing equations for this problem are firstly derived by further considering the nonlinearity of soil compression. Subsequently, numerical solutions are provided using finite difference methods to investigate the influence of non-Darcian flow in thermal consolidation of soft soils. The results show that non-Darcian flow in soft soils amplifies the impact of thermal gradients on the consolidation rate. The peak value of excess pore water pressure in soil with non-Darcian flow is significantly higher than that with Darcy flow, but it decreases as the soil heats up. Non-Darcian flow exerts less effects on the thermal diffusion process within the soil; varying trends of pore water flow velocity at different depths may be opposite due to the influence of temperature-dependent permeability within the process of consolidation.</p>