With the energy demand arising globally, geothermal recovery by Enhanced Geothermal Systems (EGS) becomes a promising option to bring a sustainable energy supply and mitigate CO2 emission. However, reservoir management of EGS primarily relies on reservoir simulation, which is quite expensive due to the reservoir heterogeneity, the interaction of matrix and fractures, and the intrinsic multi-physics coupled nature. Therefore, an efficient optimization framework is critical for the management of EGS. We develop a general reservoir management framework with multiple optimization options. A robust forward surrogate model fl is developed based on a convolutional neural network, and it successfully learns the nonlinear relationship between input reservoir model parameters (e.g., fracture permeability field) and interested state variables (e.g., temperature field and produced fluid temperature). fl is trained using simulation data from EGS coupled thermal-hydro simulation model by sampling reservoir model parameters. As fl is accurate, efficient and fully differentiable, EGS thermal efficiency can be optimized following two schemes: (1) training a control network fc to map reservoir geological parameters to reservoir decision parameters by coupling it withfl ; (2) directly optimizing the reservoir decision parameters based on coupling the existing optimizers such as Adam withfl. The forward model fl performs accurate and stable predictions of evolving temperature fields (relative error1.27±0.89%) in EGS and the time series of produced fluid temperature (relative error0.26±0.46%), and its speedup to the counterpart high-fidelity simulator is 4564 times. When optimizing withfc, we achieve thermal recovery with a reasonable accuracy but significantly low CPU time during inference, 0.11 seconds/optimization. When optimizing with Adam optimizer, we achieve the objective perfectly with relatively high CPU time, 4.58 seconds/optimization. This is because the former optimization scheme requires a training stage of fc but its inference is non-iterative, while the latter scheme requires an iterative inference but no training stage. We also investigate the option to use fc inference as an initial guess for Adam optimization, which decreases Adam's CPU time, but with excellent achievement in the objective function. This is the highest recommended option among the three evaluated. Efficiency, scalability and accuracy observed in our reservoir management framework makes it highly applicable to near real-time reservoir management in EGS as well as other similar system management processes.