Abstract

Identifying fluid circulation in fracture zones (FZs) is a key challenge in the extraction of deep geothermal heat from natural reservoirs in the Upper Rhine Graben. This study focuses on permeable FZs present within the granitic basement penetrated by deep geothermal well GPK-1 at Soultz and GRT-1 and GRT-2 at Rittershoffen (France). The various temperature (T) log datasets acquired from these wells during production and at equilibrium, with the associated flow logs, allow for the unique opportunity to interpret fluid circulation at the borehole scale. All permeable FZs identified by permeability indicators measured during drilling operations and from image logs spatially coincide with positive or negative T anomalies observed in the T logs during production and/or at equilibrium. However, within the FZs, partially open fractures act as narrower paths for circulation at different temperatures. These temperatures can even be estimated with confidence if the associated flow log is available. The polarity of the T anomalies correlates with the state of equilibrium of the well and thus can change over the well history. During production, the temperature of the water inflow through the fractures can be estimated relative to the mixture of water circulating below the fractures. At thermal equilibrium, the water temperature is estimated with respect to the temperature of the surrounding rock formation. Because temperature fluxes and geothermal fluids are intimately linked, T logs are a useful, reliable, and very sensitive tool to localize the inflow of geothermal water through FZs.

1. Introduction

In the Upper Rhine Graben (URG), the underground temperature distribution is spatially heterogeneous, and a series of local anomalies with temperatures above 140°C at a 2 km depth are observed at the regional scale [14]. These temperature anomalies are interpreted as convective cells circulating along vertical local faults (Figure 1) [2, 5, 6]. Hot upwellings are mainly concentrated around N-S-striking local faults on the western side of the URG. These N-S-striking local faults, which originate in the Paleozoic granitic basement and extend through the sedimentary cover, channelize these deep fluids. Over more than 30 years, 18 deep geothermal wells have targeted these faults in Germany (Landau, Insheim, Bruchsal, and Brühl) and France (Soultz-sous-Forêts, Rittershoffen and, most recently, Illkirch and Vendenheim) (Figure 1) [712]. Temperature (T) logs obtained from these wells confirm the presence of hydrothermal convection and assert the roles of faults and fracture zones (FZs) as preferential pathways for geothermal brine [1215]. The top of the granitic basement, which is highly fractured and affected by hydrothermal alteration, presents a low geothermal gradient associated with the vertical flow of brine and is locally affected by temperature anomalies interpreted as the thermal signatures of the FZs (Figures 24) [1416].

Thus, the characterization of permeable fracture zones that act as primary potential paths for fluid flow is a key challenge for the development of geothermal reservoirs in granitic rocks. Experience shows that the characterization of permeable FZs by flow logs (obtained with a spinner tool) is difficult because flow logs are hard to obtain and interpret in these reservoirs. In contrast, T logs are cheap and easily acquired. T anomalies are interpreted as the thermal expressions of permeable FZs in several deep wells in geothermal systems worldwide [13, 14, 1719]. However, in this paper, the datasets of the T logs acquired from wells at equilibrium and during production, combined with flow logs when they are available, offer a unique chance to further interpret T logs at the borehole scale. The reinterpretation of T logs at equilibrium and during production and their comparison with the permeability indicators and structural data from the FZs in well GPK-1 at Soultz and wells GRT-1 and GRT-2 at Rittershoffen provide a new perspective to better understand brine circulation through FZs in the granitic basement.

After a structural analysis of the main permeable fracture network that has channelized geothermal brines at the first order, the attention is focused on a comparison of T logs obtained at equilibrium and during production with flow logs if they are available from the wells. Finally, a conceptual model of fluid circulations through the FZs and the associated T anomalies is proposed at the borehole scale.

2. Soultz and Rittershoffen Geothermal Sites

Both projects exploit the hot upwelling of geothermal brine concentrated around the Soultz and the Rittershoffen faults (Figure 1). Extensive studies have identified convection associated with these hot upwellings; these studies have used numerical models on both the regional scale [4, 20] and local scale [21]. The role of faults as preferential fluid pathways is evident in other geothermal systems [2225].

The pilot project of Soultz-sous-Forêts was initiated in the late 1980s. The target of the project was the development, hydraulic testing, and modeling of two Enhanced Geothermal System (EGS) heat exchangers within the granitic basement at depths of 3.5 and 5 km [15, 26, 27]. In 1987, exploration well GPK-1 was drilled to a depth of 2 km and extended to a depth of 3.6 km in 1992 after hydraulic stimulation (Figure 1). At Soultz, geothermal water was pumped from the production well (GPK-2) and reinjected at a lower temperature into the injection wells (GPK-3 and GPK-4) after delivering geothermal energy through a heat exchanger to a binary power plant [10].

The geothermal project in Rittershoffen was initiated in 2008 and advised by lessons learned in Soultz. It is located less than 10 km east of Soultz (Figure 1). Injection well GRT-1 was drilled to a depth of 2.5 km in 2012 and thermally, chemically, and hydraulically stimulated in 2013. Production well GRT-2 was drilled to a depth of 2.5 km and is highly deviated accordingly to the geometry of the Rittershoffen fault (Figure 1). Well GRT-2 was not stimulated and thus, is defined as an hydrothermal well. The geothermal reservoir is mainly located in the granitic basement.

Geochemical studies indicate that the same brine is circulating within the large-scale fractured reservoir between Soultz and Rittershoffen [28, 29]. The reservoir was thoroughly investigated with a large amount of data (obtained from core samples, image logs, geophysical logs, microseismic data, vertical seismic profile, seismic reflection, etc.); these data were used to develop structural models of the multiscale fracture network associated with the Soultz and Rittershoffen faults [13, 3033] as well as models of hydraulic circulation [34].

At the borehole scale, structural information from the granitic basement indicates that the highest density of fractures occurs in the first hundred meters of the hydrothermally altered granite. Data from acoustic image logs indicate values of 1,15 frac/m in GPK-1 in the first 500 m (data from a BoreHole TeleViewer, [35]) (Figure 2), 1,15 frac/m in GRT-1 in the first 200 m (Figure 3), and 0,75 frac/m in the deviated well GRT-2 (data from ultrasonic borehole images, [36]) (Figure 4). The fractures are mainly oriented NNW-SSE, with a steep dip mainly westward in GPK-1 (data from BHTV, [35]). In well GRT-1, the fractures strike mainly N10°E and dip approximately 60°W, whereas they are more scattered in GRT-2, with the main set oriented N160° to N-S and dipping 60°W (data from UBI, [36]). Hydraulic tests in GPK-1 indicate that 95% of the fluid flow enters the rock mass at only 10 discrete flow points [14]. These tests are consistent with cores from the permeable FZ, where a high density of partially open fractures promotes fluid circulation. These partially open fractures are centimeters thick and are mainly filled by geodic quartz [37, 38]. An investigation of a quartz vein in the exploration borehole EPS-1 revealed at least seven generations of quartz, indicating a pulsed system of circulation [39]. Surrounding the core is a damage zone, with a width that ranges from centimeters to tens of meters. In this altered and porous zone, primary minerals are partly dissolved, and illite precipitation dominates the granite matrix [38, 40]. The degrees of alteration in the granite were determined based on the observations of the illitization amount in the core samples or cuttings [37, 4143].

These permeable FZs mainly intersected in the granitic section, which is dominated by low temperature gradients associated with convection. The T gradients are approximately 20°C/km in GPK-1, nearly in GRT-1, and 18°C/km in GRT-2 (Figures 24). All permeable FZs have T anomalies that spatially correlate with null observations in the T logs (Figures 24) [1316]. Negative T anomalies are interpreted as the remnant cooling of porous damage zones after mud invasions during drilling operations and water injection during stimulation operations [15]. However, is this interpretation still true several months after drilling operations and hydraulic tests, when the well is at thermal equilibrium? Could we propose another interpretation based on the internal organization of the FZs? Is the intensity of the anomaly linked with the degree of permeability of the FZs? How could we explain the polarity and variations in the intensity of T anomalies over time?

3. Materials and Methods

3.1. Permeable FZs from Mud Logs and Image Logs

FZs that correspond to clusters of fractures partially open at the borehole scale and surrounded by a halo of hydrothermal alteration were identified in previous studies in Soultz and Rittershoffen [16, 44]. The core of the FZ is composed of opened fractures and quartz veins, and the damage zone is composed of illitic minerals [38]. Other secondary minerals could be present (carbonates, sulfates, sulfides, and Fe-oxides).

Permeable FZs were identified using permeability indicators from several mud logs as follows: (i)Drilling mud losses and brine outflow are linked to permeable FZs [1315, 33, 38, 45, 46](ii)The occurrence of alkanes, as well as other gases such as helium, CO2, and radon, indicates permeable fracture zones [47, 48](iii)The mineralogical investigation of rock cuttings or core samples, when available, is an efficient method for identifying hydrothermally altered zones in the granite [13, 49]. Brine circulation through a permeable fracture zone leads to partial sealing by secondary geodic quartz and clay mineral deposits, which are easily detectable in cuttings by visual inspection and in the laboratory by XRD [39, 40, 43, 50]. When cuttings are not available or reliable, gamma ray (GR) logging, which measures natural radioactivity, is a good indicator of hydrothermal alteration [42]. In GR, the occurrence of geodic quartz is associated with sharp localized negative anomalies, whereas clay minerals, which are K-bearing minerals, are associated with positive anomalies that can extend several meters

Finally, once the permeability indicators are identified, acoustic image logs allow for the identification of the geometry of the FZs and the acoustic aperture of the fractures at the borehole scale.

3.2. T Logs and Flow Logs

The T anomalies observed in a temperature profile are considered the most reliable of the permeability indicators [14, 16, 18, 24]. T logs can be acquired at thermal equilibrium, i.e., after several weeks without well operations, or during production. When T logs are acquired during production, a flow log can also be acquired, and the flow is measured based on the speed of the rotation of the helix (spinner tool). T logs acquired at equilibrium are acquired downward; T logs and flow logs acquired during production are acquired upward. T logs acquired during production allow for temperature estimation of the water circulating into the fracture zone with the following equation (assuming a first-order approximation) (Figure 5): where is the temperature of the water coming from below the fracture, is the flowrate of the water coming from below the fracture, is the temperature of the water coming from the fracture, is the flowrate of the water coming from the fracture, is the temperature of the mix of water above the fracture, and is the flowrate of the mix of water above the fracture.

At equilibrium, T anomalies were identified primarily with T logs acquired before stimulation (Figures 24). Thermal equilibrium must not be influenced by cold mud or cold water stored in fractures after drilling and hydraulic operations. Poststimulation T logs were compared to prestimulation T logs. It is possible to observe some vertical depth shifts between T anomalies and fracture zones from core samples and image logs. These shifts could have several explanations. All logs were not acquired at the same time, and the process of depth matching among all logs is complicated. Moreover, temperature variations led to cable elongations that were not corrected for the well-logging data. The logs were shifted downward or upward based on the depths of the main open fractures observed in the acoustic logs.

4. Results

4.1. Permeable FZs from T Logs at Thermal Equilibrium

Four permeable FZs are identified in GPK-1, and they are all associated with negative thermal anomalies at thermal equilibrium (Table 1 and Figure 2). The borehole images of GPK1-FZ2 and GPK1-FZ3 are available in Vidal et al. [16, 43, 44]. Another negative anomaly is visible at a depth of 2100 m, but based on hydraulic tests, the FZ was identified as permeable after stimulation [51].

In GRT-1, one FZ was identified as permeable after stimulation and is associated with a negative anomaly at thermal equilibrium (Table 1 and Figure 3). GRT1-FZ1 includes the two permeable fractures observed in the granitic basement, as described in Vidal et al. [16]. A poststimulation anomaly is observed in the T log at 2236 m MD, but it is not spatially correlated with permeability indicators obtained before stimulation. Four permeable FZs are identified in GRT-2; two are associated with positive anomalies (GRT2-FZ1 and GRT2–FZ2), and two are associated with negative anomalies (GRT2-FZ3 and GRT2-FZ4) (Figure 4). GRT2-FZ1 includes the two uppermost permeable fractures observed in the table 1 of Vidal et al. [16], and GRT2-FZ2 includes the five deepest fractures. The borehole images of GRT-1 and GRT-2 are available in Vidal et al. [16].

Permeable fractures are mainly oriented N170°E and dip westward (Figure 6). Fractures that dip westward are associated with positive or negative T anomalies. Fractures that dip eastward are also associated with T anomalies that are positive or negative. The polarity of T anomalies is not directly linked to the fracture orientation or the dip direction. Thus, it is not possible to anticipate the polarity of the anomalies with respect to the fracture geometry.

4.2. Fluid Circulation from T Logs during Production
4.2.1. GPK-1

The permeable FZ in GPK1-FZ4 is presented in Figure 7. In September 1996, during production at 10.8 L/s, a flow log and a T log were acquired. The first inflow of water into the well was observed at 3496 m MD, with a measured flow of 1.5 L/s. This inflow depth corresponds with a negative T anomaly (148°C) in the T logs acquired during production. These anomalies are spatially correlated with an open fracture, which is visible in the acoustic image log, striking N170°E and dipping westward 70° at 3496 m MD. Then, several small positive anomalies are observed and associated with a flow increase at the same depth as the open fracture, striking N160°E and dipping 60°W at 3489 m MD. All of these small anomalies are assumed to originate from the same fracture. Above this, the flow is 2.7 L/s, and the measured temperature is 149°C. Following the equation described in Section 3.2, the temperature of the water that arrives in this major fracture is approximately 150°C:

At equilibrium (five months after the hydraulic tests), the water inflow through GPK1-FZ4 is associated with a negative anomaly in the T log acquired in March 1993.

The subsequent T logs acquired from well GPK-1 during production are more affected by hydraulic circulation than the other wells and thus are less reliable.

4.2.2. GRT-1

In GRT-1, during hydraulic tests in January 2013, a T log and a flow log were acquired during airlift production at 12,3 L/s. However, the flow log is not entirely useful; anomalies are observed but are barely quantifiable. T logs at equilibrium were acquired pre- and poststimulation in April and December of 2013, respectively (Figure 8).

In the production T log, a small negative anomaly is observed at 2233 m MD, which is not correlated with the flow anomaly (Figure 8(a)). This T anomaly is not observed in the prestimulation T log at thermal equilibrium. However, a T anomaly is observed in the poststimulation T log at 2235.5 m MD. This anomaly corresponds with an open fracture striking N180°E and dipping 60° westward at 2236.5 m MD.

The main permeable FZ of GRT1-FZ1 extends from 2325 to 2368 m MD. Two major open fractures associated with mud losses were also identified. The first one, which is found at the top of the FZ, strikes N20°E and dips 50° westward and is located at 2326 m MD [16]. Although mud losses were observed during drilling operations, no T anomaly was observed during production at this depth. A small negative anomaly is observed in the poststimulation T log (Figure 3).

In the production T log, a small negative anomaly is observed at 2341 m MD that is not correlated with a flow anomaly (Figure 8(b)). This T anomaly is incorporated into a larger anomaly evident in the prestimulation T log at thermal equilibrium and in the poststimulation T log (Figure 3). These anomalies correlate with an open fracture striking N25°E and dipping 60° westward at 2345 m MD (Figure 8(b)).

The second open fracture of GRT1-FZ1 is associated with mud losses and is located at 2368 m MD. It is striking N175°E and dipping 65° westward. This fracture controls 2/3 of the total flowrate and is at the interface between highly altered and fractured granite above and low altered granite below [16]. Thus, this fracture permits the first inflow of water into the well. Surprisingly, the flow log is not consistent with this observation and is thus considered unreliable for the evaluation of the absolute flow associated with the fractures in the well. The flow anomaly is associated with a T anomaly (+3°C) at 2364 m MD (Figure 8(c)). A second T anomaly (+5°C) at 2350 m MD is observed. This positive T anomaly is not clearly associated with an open fracture but is associated with a cluster of thin fractures. However, the calculation of the water temperature at this depth is not possible because the flow log is not reliable. At equilibrium, the prestimulation log indicates a large negative T anomaly, whereas the poststimulation T log indicates two negative anomalies that correspond with those observed during production (Figure 3).

4.2.3. GRT-2

Only one temperature log at thermal equilibrium was acquired in this well in September 2014. Two positive T anomalies are observed for GRT2-FZ1 and GRT2-FZ2 (Figure 4). These anomalies are the only positive ones observed in the well, but as they are observed at thermal equilibrium, without previous injection into the well, they are considered reliable. They indicate water inflows hotter than the surrounding rock formation. Unfortunately, the temperature of this water cannot be precisely calculated because flow logs and T logs acquired during production are not available. Both anomalies are associated with open fractures at the borehole scale, which mainly strike N170°E and dip westward [16].

5. Discussion

5.1. Thermal Signature of Permeable FZs

GPK1-FZ4 corresponds to the first inflow of water into the well and thus marks the beginning of thermal convection. The first inflow of water into the well is at 148°C through a fracture at 3496 m MD (Figure 7). The negative anomaly in the production T log is explained by an inflow of water colder than the deeper reservoir during production. The positive anomaly described above is interpreted as the inflow of water from a fracture at 3489 m MD, which is hotter (150°C) than water from the fracture at 3496 m MD (148°C). At thermal equilibrium, GPK1-FZ4 is associated with a negative anomaly because the water inflows have temperatures of 148°C and 150°C and are probably colder than the temperature of the deep granite at equilibrium (Figure 2).

GRT1-FZ1 also corresponds to the first inflow of water into the well and thus to the beginning of thermal convection. Positive T anomalies are observed in the bottom part of the FZ at 2365 and 2350 m MD because geothermal water that comes from the fractures is hotter than the water coming from the underlying reservoir (Figure 8(c)). At thermal equilibrium, these hot water inflows are no longer observed, but the negative T anomaly does not suggest that hot water does not circulate. This means that the circulating water is probably colder than the reservoir at equilibrium. These fractures correlate with a flow anomaly. However, T anomalies during production at 2341 m MD and 2235 m MD indicate inflows of water colder than the mix of water below, a trend not seen in the flow log (Figures 8(a) and 8(b)). The T anomaly at 2231 m MD is associated with an FZ enhanced after stimulation at 2236 m MD because it is associated with a negative T anomaly in the poststimulation T log (Figure 8(a)).

GRT2-FZ1 and GRT2-FZ2 exhibit the only two positive T anomalies at equilibrium observed in the wells of Rittershoffen. Positive T anomalies are probably less common than negative T anomalies at equilibrium because the reservoir at equilibrium is generally hotter than the circulating water. Fractures are direct paths for cold mud or water during drilling and hydraulic operations and thus cool faster than the rest of the rock formation. The greater the quantity of cold water that is injected, the fewer positive anomalies visible. Because well GRT-2 was not stimulated, positive anomalies are probably more common in that well than in the stimulated well GRT-1 and GPK-1.

As observed by Evans et al. [14], temperature flux and geothermal fluid flow are intimately linked. T logs are better tools than flow logs for examining water circulation because T logs are more sensitive to small water inflows than flow logs, which require massive water inflow to produce an observable modification in the rotation speed of the helix. However, the estimation of the absolute water temperature circulating in a fracture requires a reliable associated flow log, which is often not available.

5.2. Water Circulation through Permeable FZs

Based on T and flow logs, most of the water circulation in the wells occurs in the first hundreds of meters of the granitic basement, where a high density of permeable FZs is observed: in the first 500 meters of the Soultz basement (1400 to 2100 m MD in GPK-1) and in the first 200 meters of the Rittershoffen basement (2200 to 2400 m MD in GRT-1 and 2500 to 2800 m MD (i.e., 2200 to 2400 m TVD) in the deviated well GRT-2) (Figures 24). The partially open fractures are several centimeters wide inside the FZs and are probably interconnected to each other over short distances, leading to a rather complex 3D organization that has channelized the water circulation [45, 46]. Inside a permeable FZ, geothermal water can circulate at different temperatures through different partially open fractures, representing multichannel pathways for geothermal fluids, as observed in GPK-1, where the water coming from the upper fracture is hotter (Figure 7). In this case, water coming from fractures less than 10 m apart circulates at different temperatures, which is surprising because circulations coming from narrow fractures should reach a state of equilibrium and a homogeneous temperature. Water may circulate through the FZs too fast to reach equilibrium. From the dataset presented in this article, the velocity and direction of these circulations cannot be determined.

The orientations and dips of the fractures cannot be linked to the polarity of the T anomalies at equilibrium (Figure 6). Negative T anomalies at equilibrium are observed more often than positive ones in the Soultz and Rittershoffen wells, but the negative anomalies are linked to the contrast between the temperature of the fluid circulating through the FZs and the temperature of the reservoir at equilibrium, which is generally hotter than the circulating water. During production, a positive T anomaly indicates an inflow of water hotter than the mix of water circulating below (Figure 9). At equilibrium, this anomaly could be positive if the geothermal water coming from the fracture is hotter than the surrounding rock formation or negative if the geothermal water is colder than the surrounding rock formation. T logs acquired during production are representative of the mix of water present during circulation, whereas T logs at thermal equilibrium reflect the temperature of the water leaching through the FZs relative to the temperature of the surrounding rock formations.

Contrary to the interpretation proposed by Genter et al. [15], the negative T anomalies in T logs at equilibrium are not a result of cooling of the FZs but only indicate that the circulation of water is colder than the surrounding formation. However, equilibrium could be difficult to identify depending on the volume of cold water previously injected and the time since the last injection.

The recurrence of negative anomalies could also be influenced by hydrothermal alteration (Figure 9). As modeled by Gentier et al. [52], during cold injection, thermal microcracking of quartz within the FZ is observed, which creates preferential flow paths and thus leads to preferential cooling in these fractures. All permeable fractures are associated with quartz veins, which could enhance the thermal effect of cold reinjection through the FZs (Figures 24). Moreover, Tung et al. [53] and Alevizos et al. [54] suggested that occurrences of illite could enhance fault creeping, which can release fluid and trigger endothermic chemical reactions. In the case of GPK-1, Sausse et al. [33] identified the correlation between GPK1-FZ4 and the fault zone imaged in the reservoir. The influence of clay minerals, particularly illite, on the fault movements at Soultz was already proposed by Meller et al. [49]. In the case of GRT-1 and GRT-2, preliminary geometrical modeling suggests a possible correlation between GRT1-FZ1 and GRT2-FZ2 and a fault zone imaged in the reservoir [7, 16]. The influence of these endothermic reactions on permeable FZs and the associated T anomalies should be investigated in further studies.

6. Concluding Remarks

In this study, a network of permeable FZs was identified in well GPK-1 at Soultz and wells GRT-1 and GRT-2 at Rittershoffen based on the reinterpretation of T logs obtained during production and at thermal equilibrium; these T logs were correlated with permeability indicators measured during drilling operations, acoustic image logs, and flow logs obtained during hydraulic tests. All permeable FZs present similar characteristics as follows: (1)They are composed of cm thick fractures partly filled by quartz(2)They are associated with T anomalies during production and/or at thermal equilibrium

This study presents a unique dataset of T logs and associated flow logs (if available) acquired during production and at thermal equilibrium at different times in well histories. T logs are useful and inexpensive tools for the localization of permeable FZs because they are more sensitive to small water inflows into wells than flow logs. Moreover, the comparison of T logs acquired in production and at thermal equilibrium provides a new interpretation of temperature anomalies. Based on these data, the temperature of water coming from partially open fractures was estimated. The study shows that geothermal water could circulate at different temperatures and in different fractures a few meters apart within the same FZ, which suggests a new interpretation of T anomalies. At temperature equilibrium, T anomalies are positive or negative with respect to the temperature of the surrounding rock formations. Because the surrounding formations are often hotter than the fractures, which are direct paths for fluid and thus cool faster than the surrounding rocks, negative anomalies are more often observed than positive ones. The nature of secondary deposits could also enhance the cooling effect observed in the T logs. In contrast, during production, anomalies are positive or negative with respect to the mix of water circulating below the fracture and not with respect to the surrounding formation. Therefore, anomaly polarity could change over time depending on the equilibrium state.

Data Availability

Please contact ES-Géothermie for data requests at [email protected].

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 792037 (MEET project). The study was partly funded by ADEME (French Agency of Energy and Environment) in the framework of the DEEP-EM project. The authors acknowledge the GEIE EMC and ECOGI for providing the data. The authors appreciate the helpful and constructive remarks of the two anonymous reviewers and the lead guest editor of this special issue, Dr Fabrizio Agosta.