The geological setting is typically complex, with horizons, faults, and intertwined patterns of natural fractures, and the nearwell region is often hydraulically fractured to enhance well communication. Therefore, in order to fully utilize the potential of shallow geothermal heat storage, numerical simulations are imperative. In this work, we show how to practically model such systems, including generation of computational grids with a large number of wells and fractures, numerical discretizations with discrete fractures (DFN), and complex storage strategies with multiple wells working together under common group targets. We also discuss how adjoint-based methods can be used to tune model parameters (e.g., well injectivities, rock properties, and hydraulic fracture conductivities) so that the model fits observed data, and to find well controls (e.g., rates and temperatures) that optimize storage operations.
The methodology is demonstrated using a set of real geothermal heat storage projects currently under development, and we highlight important challenges and our suggested solutions related to each of them.