This paper investigates coupled thermal, mechanical and hydraulic phenomena in deep rock mass especially for underground heat storage system. Firstly, concepts of underground heat storage were presented and coupling phenomena in this area were illustrated. In order to understand the basic mechanism of thermal, hydraulic and deformation behavior in rock cavern disturbed by thermal gradient about 10$0^{circ}C$, various numerical experiments were conducted using several codes. The study involves the behavior of fractured rock mass including rock joint. In spite of the limitation of codes modelling fully coupled effects, these codes could be applied in analysis of underground heat storage. The heat loss in rock mass, which is a major factor in heat storage, is insignificant in all results.