Abstract: The elastoplastic model and the relative 2D FEM code were extended and improved to analyze coupled thermo-hydro-mechanical process in saturated-unsaturated porous media through introducing the percolation and migration equation,so that the code could be used for solving the temperature field,flow field,stress field and nuclide concentration field at the same time.Then through a computation example of an assumed nuclear waste repository the distribution and variation of temperature,pore pressure and nuclide concentration in the near field were investigated under the condition of thermo-hydro-mechanical-migratory coupling.