Efficient Numerical Framework for Geothermal Energy Production Optimization in Fracture‐Controlled Reservoirs

We describe a numerical framework designed to automate the search for the placements of injection and production wells in hot fracture‐controlled reservoirs that sustainably optimize geothermal energy production, where we consider deviated multiwell layouts (smart multiwell arrangement), and the framework is based exclusively on open‐source software. This search is carried out via 3D simulations of groundwater flow and heat transfer. We model the reservoirs as discrete fracture networks (DFN) in which the fractures are 2D manifolds with polygonal boundaries embedded in a 3D porous medium. The wells are modeled via the immersed boundary method. The flow and heat transport in the DFN‐matrix system are modeled by solving the balance equations for mass and energy, while expressing the momentum balance by the Darcy law. The fully developed numerical framework combines the finite element method with semi‐implicit time‐stepping, algebraic flux correction, and approximation of the wells via the non‐matching approach. To perform the optimization, we use various gradient‐free algorithms. We present the results of verification and validation tests with DFNs of simplified structure and realistic physical parameter values.