In this paper, through establish the momentum equilibrium equation, the continuity equation for water, the energy balance equation and the elastoplastic matrix, one coupled thermo-hydro-mechanical model for saturated-unsaturated porous media is suggested. Taking the Prototype Repository Project (PRP) in-situ test for geological disposal of high level radioactive nuclear waste as the modeling object, the numerical analysis for the coupled thermo-hydro-mechanical processes in near field of a conceptual nuclear waste repository was carried out by using a 2D FEM code developed by the author, and the changes and distributions of temperature, saturation and stresses in the engineered barrier and surrounding rockmass were investigated. The part results obtained were compared with the computation data for the same problem obtained by THAMES code of Japan, and it was shown that the regulations of both results have certain identity qualitatively and quantitatively, but there exists some differences between them. The author also tried to analysis the reasons and got some realization.
