Abstract: The distribution of the deep geothermal field and the tectonic thermal control mode in the Yingshan area were not understood systematically, which would bring in constraints to the sustainable exploration of local geothermal resources. Based on the field survey on continuous temperature, fracture measurements in the outcrops of granite rock, slug-test, and pumping test in the boreholes, a three-dimensional model coupling fluid flow and heat transfer process in fractured rock reservoir was established. Then the distribution of fluid pressure temperature and the Darcy velocity in the whole geothermal reservoir were calculated considering different treatment combinations of main faults. The results show that the deep temperature field and fluid pressure distribution are greatly affected by the faults system. At the location of bedrock and impermeable fault, the fluid movement is very weak, and the heat transfer is mainly controlled by heat conduction, while in the high permeable fault, it is mainly controlled by heat convection. Low temperature and fluid pressure occurred in the deep, high-permeable fault zone whereas high fluid pressure occurred in the impermeable fault zone. At a depth larger than 2000 m, the zones with high temperatures and intense hydraulic exchange coincide with the conductive fault zones, which indicates that the high-permeable faults are the most important channel for seepage and heat transfer in the reservoir. Therefore, the northeast faults are the dominant channel for fluid migration and heat conduction in the Yinshan area. Hot water migrates upward along the northeast faults and is blocked by the north-south fault. The hot spring is formed at the intersecting position of two groups of faults. This study is helpful for scientific, sustainable exploration of geothermal resources in the Yingshan area.