Abstract: Based on the stress equilibrium equation, the seepage continuity equation and the energy conservation equation, with consideration of the effect of the thermal expansion coefficients between the soil particles and the pore water, a three-dimensional thermo-mechanical governing equation for saturated soil with consideration of the thermohydraulic coupling is established.Using the Fourier transform and the Laplace transform to derive the governing equations and obtain the solution in the transform domain, then to obtain the analytic solution of the point heat source in the real domain, as well as the analytic solution of the temperature, the pore pressure and the displacement of the thermally consolidated soil with two parallel cylindrical heat sources The results show that the increase of the span ratio will lead to a weakening of the interaction between the two heat sources. The decrease of the thermal consolidation coefficient will lead to the increase of the pore pressure and the displacement peak.