Abstract: A numerical calculation model of fractured rock mass artificial seepage freezing with two holes was set up，and the problem of coupling and phase change during the freezing process were simulated by using the COMSOL Multiphysics finite-element analysis software．The influence of the fracture on the distribution of freezing temperature field and seepage field was analyzed．The distribution law of freezing temperature field corresponding to the different flow velocities was researched．The results show that the existence of fracture reduces the thickness of the frozen wall in the upstream，and the isotherms outside the frozen wall show“heart-shaped”distribution．The seepage velocity has a great influence on both the frozen wall development and the closure time．If the seepage velocity is increasing，the development of freezing wall will become slower and the closure time will become longer accordingly．In this paper，when the seepage velocity is less than 7 m/d，the closure time shows a linear relationship with the seepage velocity basically．But when the seepage velocity is more than 7 m/d，the closure time mount up sharply，the closure time tends to be infinite once the velocity exceeds 10 m/d，namely reach the unable closure stage．

Key words: freezing temperature field；seepage；the fractured rock mass；freezing wall；closure time