In order to alleviate the environmental crisis and improve energy structure, countries from all over the world have focused on the hot dry rock geothermal resources with great potential and with little pollution. The geothermal heat production from enhanced geothermal system comes with complex multi-field coupling process, and it is of great significance to study the temporal and spatial evolution of geothermal reservoir. In this work, a practical numerical model is established to simulate the heat production process in EGS, and the comparison of thermal-hydraulic (TH), thermal-hydraulic-mechanical (THM) and thermal-hydraulic-mechanical-chemical (THMC) coupling in geothermal reservoir is analyzed. The results show that the stable production stage of the three cases is approximately 5 years; however, compared with TH and THMC coupling, the service-life for THM coupling decreased by 1140 days and 332 days, respectively. The mechanical enhanced effects are offset by the chemical precipitation, and the precipitation from SiO2 is much larger than the dissolution of calcite.