In typical EGS fracturing projects, the fracture network is denser in the region close to the fracturing wells. The distribution of fractures has a significant impact on the heat extraction performance of enhanced geothermal system (EGS). In order to avoid premature thermal breakthrough and improve the thermal extraction efficiency, this paper studies the heat extraction performance of EGS in heterogeneous fractured reservoirs under different injection-production methods. First, a heterogeneous discrete fracture network generation method based on relative pressure distribution is developed. A thermo-hydraulic-mechanical (THM) coupling model is established to simulate the fluid flow and heat transfer in reservoirs. The heat extraction performance under different production and injection modes is compared with different fracture densities. The results show that injection in the dense side of the anti-parallel well configuration will be more efficient for heat extraction. Under the condition of same injection mass flow rate, the heat extraction performance is positively correlated with the fracture density. Although sparse-side injection is not conducive to fluid flow, under the condition of same injection mass flow rate with anti-parallel well configuration, it will improve the thermal extraction performance.