Due to the multi-well enhanced geothermal systems can significantly improve the heat transfer efficiency and the multi-branch radial well geothermal system can reduce the well construction cost, it is necessary to investigate the fracture propagation and the heat mining performance of the multi-branch radial well geothermal system. Based on the constructed THM-D coupling model, the fracture propagation rules of multi-branch radial wells under the co-function of low-temperature induced thermal stress and radial well induced stress is explored. Then the heat mining prediction study during the heat extraction stage is conducted based on the generated fracture. Furthermore, the influence of radial well characteristic parameters on fracture propagation and the corresponding heat mining performance of multi-branch radial well geothermal system is discussed. The results show that the smaller radial well angle in multi-branch radial well system can generate a larger scope of fracture network between the included angle under the co-function of radial well induced stress and low-temperature induced thermal stress, which will correspond to a larger heat extraction scope in heat mining stage. When the radial well included angle along the direction of maximum horizontal in-situ stress is 120°, the fracture network obtained by fracturing is the most complex. Excessive branching of radial wells can lead to significant induced stress in the plane of radial wells, which will affect the propagation of the fracture network in the direction of longitudinal connection between the doublet wells, thereby affecting the corresponding heat mining performance. With the increase of radial well length, the ability to induce the formation of branching fractures raises, resulting in more complex fracture network. In the heat mining stage of 10-year development, the outlet temperature and heat mining rate respectively increased by 5.83% and 6.9% when the radial well length is 75 m compared to 50 m.