A 3-D numerical simulation-based heat production performance research for enhanced geothermal system with two horizontal wells and rectangular multiparallel fractures