[1] The Dead Sea Transform (DST) comprises a boundary between the African and Arabian plates. During the last 15-20 m.y. more than 100 km of left lateral transform displacement has been accumulated on the DST and about 10 km thick Dead Sea Basin (DSB) was formed in the central part of the DST. Widespread igneous activity since some 20 Ma ago and especially in the last 5 m.y., thin (60-80 km) lithosphere constrained by seismic data and absence of seismicity below the Moho, seem to be quite natural for this tectonically active plate boundary. However, surface heat flow values of less than 50-60 mW/m 2 and deep seismicity in the lower crust (deeper than 20 km) reported for this region are apparently inconsistent with the tectonic settings specific for an active continental plate boundary and with the crustal structure of the DSB. To address these inconsistencies which comprise what we call the "DST heat-flow paradox," we have developed a numerical model that assumes an erosion of initially thick and cold lithosphere just before or during the active faulting at the DST. The optimal initial conditions for the model are defined using transient thermal analysis. From the results of our numerical experiments we conclude that the entire set of observations for the DSB can be explained within the classical pull-apart model assuming that the lithosphere has been thermally eroded at about 20 Ma and the uppermost mantle in the region have relatively weak rheology consistent with experimental data for wet olivine or pyroxenite.

Geothermal district heating development has been gaining momentum in Europe with numerous deep geothermal installations and projects currently under development. With the increasing density of geothermal wells, questions related to the optimal and sustainable reservoir exploitation become more and more important. A quantitative understanding of the complex thermo-hydraulic interaction between tightly deployed geothermal wells in heterogeneous temperature and permeability fields is key for a maximum sustainable use of geothermal resources. Motivated by the geological settings of the Upper Jurassic aquifer in the Greater Munich region, we develop a computational model based on finite element analysis and gradient-free optimization to simulate groundwater flow and heat transport in hot sedimentary aquifers, and numerically investigate the optimal positioning and spacing of multi-well systems. Based on our numerical simulations, net energy production from deep geothermal reservoirs in sedimentary basins by smart geothermal multi-well arrangements provides significant amounts of energy to meet heat demand in highly urbanized regions. Our results show that taking into account heterogeneous permeability structures and a variable reservoir temperature may drastically affect the results in the optimal configuration. We demonstrate that the proposed numerical framework is able to efficiently handle generic geometrical and geological configurations, and can be thus flexibly used in the context of multi-variable optimization problems. Hence, this numerical framework can be used to assess the extractable geothermal energy from heterogeneous deep geothermal reservoirs by the optimized deployment of smart multi-well systems. Keywords Porous and fractured geothermal reservoir modeling • Geothermal multi-well configurations • Finite element method • Thermo-hydraulic coupling • Optimization • Open-source software PACS 65M60 • 76S05 • 86-08 • 86A20 Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

No abstract

Abstract. This work studies intensively the flow in fractures with finite hydraulic conductivity intersected by a well injecting or producing at constant pressure, either during an injection or production well test or the operation of a production well. Previous investigations showed that for a certain time the reciprocal of flow rate is proportional to the fourth root of time, which is characteristic of the flow regime known as bilinear flow. Using a 2D numerical model, we demonstrated that during the bilinear flow regime the transient propagation of isobars along the fracture is proportional to the fourth root of time. Moreover, we present relations to calculate the termination time of bilinear flow under constant injection or production well pressure as well as an expression for the bilinear hydraulic diffusivity of fractures with finite hydraulic conductivity. To determine the termination of bilinear flow regime, two different methods were used: (a) numerically measuring the transient flow rate in the well and (b) analyzing the propagation of isobars along the fracture. Numerical results show that for low dimensionless fracture conductivities the transition from bilinear flow to another flow regime (e.g., pseudo-radial flow) occurs before the pressure front reaches the fracture tip, and for high dimensionless fracture conductivities it occurs when the pressure front arrives at the fracture tip. Hence, this work complements and advances previous research on the interpretation and evaluation of well test analysis under different reservoir conditions. Our results aim to improve the understanding of the hydraulic diffusion in fractured geologic media, and as a result they can be utilized for the interpretation of hydraulic tests, for example to estimate the fracture length. Highlights. The reciprocal of flow rate is proportional to the fourth root of time. The migration of isobars in the fracture is proportional to the fourth root of time. For low dimensionless fracture conductivities, bilinear flow ends before the pressure front reaches the fracture tip. For high dimensionless fracture conductivities, bilinear flow ends when the pressure front reaches the fracture tip. Isobars accelerate when they approach the fracture tip.