Operational decisions play a critical role in governing fluid and heat transport in Enhanced Geothermal Systems (EGS), particularly in crystalline rocks where flow pathways are highly heterogeneous. This study combines block-scale laboratory experiments with high-resolution thermal-hydraulic simulations in CMG-STARS and optimization routines in CMOST to investigate the influence of injection rate and production pressure constraints on energy recovery. Two well configurations were examined: a five-spot fracture-caging layout and a binary well system. Laboratory observations of fracture development and injectivity were used to construct permeability fields that capture the dominant fracture-controlled flow architecture, while numerical simulations were used to evaluate thermal circulation scenarios under prescribed operational conditions. This workflow links experimentally constrained fracture geometry with systematic thermal-flow analysis, allowing the effects of operational control to be evaluated within physically grounded fractured domains. Preferential flow patterns developed in both systems as a result of fracture connectivity and local transmissibility, while flow imbalance and early thermal breakthrough were mitigated through dynamic pressure adjustments at producer wells. Optimization results showed that pressure-controlled production improved flow distribution and heat extraction in complex fracture networks. In the fracture caging model, the most favorable configuration recovered approximately 107.6 kJ of cumulated enthalpy, with an estimated exergy efficiency of 29.4%. In the binary well model, cumulative enthalpy recovery increased to approximately 186.8 kJ, with an exergy efficiency of 30%. The comparable exergy efficiencies suggest that once hydraulic connectivity between injection and production wells is established, the thermodynamic effectiveness of heat recovery is governed primarily by the imposed flow rate and pressure conditions, which control produced fluid volume and the thermal state of the produced water. Fracture architecture still controls flow partitioning, preferential pathways, and the timing of localized thermal breakthrough, but the integrated exergy efficiency reflects the global balance between connectivity, operational constraints and available thermal potential. These findings highlight the importance of adaptive operational strategies and demonstrate the value of integrated experimental-modeling workflows for evaluating fracture-controlled thermal recovery in EGS.