The Enhanced geothermal system in Soultz-sous-Forêts, located in the geothermal favorable Upper Rhine Graben, is a fracture-controlled reservoir that was highly investigated in the last decades generating a huge geoscientific database. Numerical reservoir models use this database to simulate the operation of the subsurface heat exchanger, yet suffer from simplifications regarding the transfer of experimental into model data, dimensional extension, and computational power and efficiency. The new extensive transient 3-D simulations, based on geophysical, geological and hydraulic data, highlight the hydraulic and transport feedback of the Soultz EGS due to convective and advective fluid flow. Developed with the goal of simulating the vast tracer test data during the reservoir-testing phase in 2005, the Finite Element Model is focusing on the main fractured zones, which connect the wells in the deep reservoir. It comprises 13 major hydraulically active faults and fractures in a 13x11x5 km extending model domain, as well as open-hole sections of the wells GPK1 to GPK4 and their casing leakages. The simulation of the tracer experiment confirms the strong heterogeneity of the reservoir and highlights the importance of a potential fractured zone, hydraulically separating the reservoir in a northern (GPK1 to 3) and southern section (GPK4). This zone tends to connect the reservoir to the main fault system by hydraulically separating GPK4 from the other wells. The calibration and sensitivity analyses provide a unique, broad understanding of the reservoir flow zones providing information on the extension of the Soultz reservoir in the future and on the fluid pathways in the deep subsurface of the Upper Rhine Graben.

This study presents a probabilistic analysis of 3D Navier‐Stokes (NS) fluid flow through 30 randomly generated sheared fractures with equal roughness properties (Hurst exponent = 0.8). The results of numerous 3D NS realizations are compared with the highly simplified local cubic law (LCL) solutions regarding flow orientations and regimes. The transition between linear and nonlinear flow conditions cannot be described with a generally valid critical Reynolds number (italicRecrit), but rather depends on the individual fracture's void geometry. Over 10% reduction in flow is observed for increased global Re (>100) due to the increasing impact of nonlinear conditions. Furthermore, the fracture geometry promotes flow anisotropy and the formation of channels. Flow perpendicular to the shearing leads to increased channeling and fluid flow (∼40% higher) compared to flow parallel to the shearing. In the latter case, dispersed flow and irregular flow paths cause a reduction of LCL validity.

HT-ATES (high-temperature aquifer thermal energy storage) systems are a future option to shift large amounts of high-temperature excess heat from summer to winter using the deep underground. Among others, water-bearing reservoirs in former hydrocarbon formations show favorable storage conditions for HT-ATES locations. This study characterizes these reservoirs in the Upper Rhine Graben (URG) and quantifies their heat storage potential numerically. Assuming a doublet system with seasonal injection and production cycles, injection at 140 °C in a typical 70 °C reservoir leads to an annual storage capacity of up to 12 GWh and significant recovery efficiencies increasing up to 82% after ten years of operation. Our numerical modeling-based sensitivity analysis of operational conditions identifies the specific underground conditions as well as drilling configuration (horizontal/vertical) as the most influencing parameters. With about 90% of the investigated reservoirs in the URG transferable into HT-ATES, our analyses reveal a large storage potential of these well-explored oil fields. In summary, it points to a total storage capacity in depleted oil reservoirs of approximately 10 TWh a−1, which is a considerable portion of the thermal energy needs in this area.

No abstract

Temperature logs have important applications in the geothermal industry such as the estimation of the static formation temperature (SFT) and the characterization of fluid loss from a borehole. However, the temperature distribution of the wellbore relies on various factors such as wellbore flow conditions, fluid losses, well layout, heat transfer mechanics within the fluid as well as between the wellbore and the surrounding rock formation, etc. In this context, the numerical approach presented in this paper is applied to investigate the influencing parameters/uncertainties in the interpretation of borehole logging data. To this end, synthetic temperature logs representing different well operation conditions were numerically generated using our newly developed wellbore simulator. Our models account for several complex operation scenarios resulting from the requirements of high-enthalpy wells where different flow conditions, such as mud injection with-and without fluid loss and shut-in, occur in the drill string and the annulus. The simulation results reveal that free convective heat transfer plays an important role in the earlier evolution of the shut-in-time temperature; high accuracy SFT estimation is only possible when long-term shut-in measurements are used. Two other simulation scenarios for a well under injection conditions show that applying simple temperature correction methods on the non-shut-in temperature data could lead to large errors for SFT estimation even at very low injection flow rates. Furthermore, the magnitude of the temperature gradient increase depends on the flow rate, the percentage of fluid loss and the lateral heat transfer between the fluid and the rock formation. As indicated by this study, under low fluid losses (< 30%) or relatively higher flow rates (> 20 L/s), the impact of flow rate and the lateral heat transfer on the temperature gradient increase can be ignored. These results provide insights on the key factors influencing the well temperature distribution, which are important for the choice of the drilling data to estimate SFT and the design of the inverse modeling scheme in future studies to determine an accurate SFT profile for the high-enthalpy geothermal environment.