The modeling on fracture networks and coupled thermal-hydraulic-mechanical processes in the near field rock surrounding a hypothetical nuclear waste repository are studied by means of tests, theoretical analysis and numerical simulation in this project. Field measurement of fracture system and long term heating-cooling THM test on surounding rock of YuccaMountain underground repository are utilized in our study. Multiple-point-geostastic method and Self avoiding walk method are introduced to estimate spatial distribution of fracture and its connectivity and roughness. An modeling method on rock mass continuum is established, considering long term evolution of permeability with mechanical behavior and phase transition and compressibility of fluid and non-equilibrium states of the processes in dual-continuum. Finally, the evolutionary processes of the stress, temperature, hydraulic behaviour in surrounding rock of the repository in 1000000 years are numercally simulated. And the uncertainty of both model and parameters are also studied.