firstname.lastname@example.org GFZ German Research Cementing operations in wellbores, especially for long casings, are often challenging Centre for Geosciences, and prone to deficiencies when not properly planned and executed. While exploring Telegrafenberg, 14473 Potsdam, Germany for and exploiting of geothermal resources at temperatures above the critical point of Full list of author information water was attempted in different drilling projects in recent years, the well design, and is available at the end of the especially the procedure to run and cement long production casings became a key article challenge for drilling engineers. For the first time, a reverse cementing job for a 2.97 km long production casing in a high-temperature geothermal well could be monitored and analyzed using a combination of permanently installed distributed fiber optic and electronic sensors as well as conventional well logging equipment. Data from the permanently installed sensors were used to monitor and evaluate the cementation process as well as the onset of the cement hydration. Based on the data, the under- standing of downhole fluid dynamics during cementation could be improved. Our analysis suggests that the cement was diluted during cement placement and partly lost into the formation. These findings can help to better prepare for future drilling ventures under similar downhole conditions. Keywords: Well cementing, Well integrity, Distributed temperature sensing (DTS), Iceland deep drilling project (IDDP-2) Introduction Cementing a casing string is considered one of the most important operations during the construction of a wellbore (Allouche et al. 2006). Not only does the cement pre- vent short mixing of fluids from different geological formations, it also gives the well construction its required strength to withstand hostile downhole environments. The completion of deep geothermal wells is challenging because of high pressure and tem- perature at depth and high radial thermal gradient during load changes (Teodoriu and Falcone 2009). High-temperature geothermal wells in Iceland are typically equipped with a conductor casing, a surface casing, an anchor casing, a production casing and a perforated liner. Such casings are typically cemented to the surface to withstand the extreme thermal and mechanical loads (except the liner which is usually just hang- ing). The practice to cement the entire casing gives rise to technical challenges during cementation of exceptionally long casings due to the differential pressure conditions © The Author(s) 2021. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/. Lipus et al. Geotherm Energy (2021) 9:5 Page 2 of 20 both at the casing and the formation (collapse and fracturing). The completion engineer must therefore be very selective in determining the correct cement weight (Salim and Amani 2012). Long cemented production casings (close to 3 km) were introduced in the high-temperature geothermal industry to explore for superhot or supercritical geother- mal fluids downhole, e.g., in the Descramble project (Bertani et al. 2018) and IDDP-1 (Pálsson et al. 2014). Long cemented sections require dedicated cementing strategies, e.g., reverse cementing in Iceland (RN-15/DEEPEGS/IDDP-2) (Fridleifsson et al. 2016) or two-stage cementation using a tieback in Italy (Venelle-2). The setting behavior of cement slurry is affected by the temperature occurring in the borehole and the cement additives. Depending on the cementing length and environment, a specific cement mixture is designed which allows for optimal pumping and a quick cement hardening onset once the cement is in place. A higher temperature reduces the setting time so the cement sets quicker compared to a cooler well temperature. A higher temperature also affects the cement rheology. The plastic viscosity and yield point decreases with increas - ing temperature (Ravi and Sutton 1990). Earlier, such resource conditions were encoun- tered (partly accidentally) and posed an insurmountable barrier to successfully complete and test such wells (Bertini et al. 1980; Fridleifsson et al. 2014; Reinsch et al. 2017). Cementing operations are typically poorly monitored and largely rely on measure- ments of density, volumetric flow rate and pressure, temperature at the surface. Down - hole displacement processes are assessed during post job cement evaluation, either by cement bond-logs (CBL) or an ultra-sonic imager tool (USIT). Fiber optic sensing technology allows to measure a variety of downhole parameters passively in the borehole in real-time. In previous studies, distributed temperature sens- ing (DTS) data from hot and cold water injection have been used to determine the loca- tion of permeable zones (fractures) and showed the propagation of the injected fluids into the formation (Hurtig et al. 1994). DTS was used to study thermal properties of the lithologies in a borehole (Förster et al. 1997; Henninges et al. 2005) and to monitor res- ervoir temperatures (Carnahan et al. 1999; Patterson et al. 2017). Moreover, fiber optic distributed data were used to perform cement job evaluations and wellbore integrity analysis during and after production tests (Pearce et al. 2009; Earles et al. 2011; Koelman et al. 2012; Reinsch et al. 2013; Bücker and Grosswig 2017; Lipus and Reinsch 2018; Raab et al. 2019). In this work, we report on DTS data acquired during cement placement along a permanently installed fiber optic cable. With depth reference to the ground surface, the fiber optic cable was installed down to 832 m behind the 9 5/8” (445–2932 m) and 9 7/8” (0–445 m) production casing in well RN-15/DEEPEGS/IDDP-2 in the Reykjanes geothermal field, south-west Iceland (Fig. 1). In addition to the fiber optic cable, eight thermocouples and one combined pressure gauge and thermocouple were installed. Table 1 lists the different installed sensors. For the first time, downhole data from a reverse cementing job in a geothermal well could be analyzed and hydraulic processes during the operation could be described. Well drilling and completion The depth information throughout the manuscript are referenced to Thor rig floor (9 m above surface), unless it is explicitly mentioned that the depth reference is to surface. The official well naming RN-15/DEEPEGS/IDDP-2 is replaced by “the well”. Lipus et al. Geotherm Energy (2021) 9:5 Page 3 of 20 Depth reference: Ground surface RN-15 –Surface casing Ø 22 ½'', Ø 571.5 x 12.0 mm, 117 lb/ft, Drill bit Ø 26'' St 37.2, welded RN-15 –Intermediate casing 1 Ø 18 5/8'', Ø 473.1 x 12.0 mm, Drill bit Ø 21'' (ID 444.0 mm) X56, welded RN-15 –Intermediate casing 2 Ø 13 3/8'', Ø 339.7 x 12.2 mm, 68 lb/ft, K55, BTC Drill bit Ø 17 ½'' IDDP-2 –Productioncasing Top 445 m: Ø 9 7/8'', Ø 250.8 x 15.9 mm, 62.8 lb/ft, T95, GEOCONN with Special Drift 8.5'' 445 –MD 2932.4 m: Ø 9 5/8‚‘, Drill bit Ø 12¼ Ø 244.5 x 12 mm, 47 lb/ft, L80, GEOCONN Fig. 1 Design of the well and landing depth of fiber optic cable (red line). The cementing section is displayed in green color Table 1 List of downhole sensors installed in RN-15/DEEPEGS/IDDP-2 with depth reference to surface Sensor Measured depth (m) Thermocouple No. 1 2632 Thermocouple No. 2 2332 Thermocouple No. 3 2132 Thermocouple No. 4 1832 Thermocouple No. 5 1532 Thermocouple No. 6 932 Thermocouple No. 7 632 Thermocouple No. 8 332 Pressure gauge and TC 1232 Fiber optic cable 832 The well is located in the Reykjanes geothermal field in south-west Iceland (Fridleif - sson et al. 2017). An already existing well (RN-15) in the geothermal field, which was drilled in 2004, was selected as the target well for the DEEPEGS/IDDP-2 project. The well was completed with a 13 3/8” production casing cemented from 0-794 m and a 12 1/4” open hole drilled to 2500 m (from surface) (Jónsson et al. 2010). The well was left as an open hole without being supported by, i.e., a perforated liner. The activity for the deepening of the well started with drilling the well vertically to 2750 m. Below that, the well was drilled directionally to the south-west (16 planned inclination) to a measured depth of 3000 m. Within the interval of the production casing, several lost circulation zones were observed. The well was drilled blind (from RN-15 to 3000 m) with a total loss of circulation without cuttings being retrieved at surface (Stefánsson et al. 2017). A cement job was performed to seal a zone near the KOP and multiple cycles of cleaning and reaming were performed. Polymer pills and bentonite were circulated in the well 3000 m 2932 m 793.8 m 292.8 m 84.4 m Lipus et al. Geotherm Energy (2021) 9:5 Page 4 of 20 before installing the 9 7/8” and 9 5/8” casing. During running the casing, the integrity of the TCs as well as the fiber optic cable was regularly assessed after each joint. As the casing string reached a depth of 2932 m (from surface), a damage on TC3 was detected. Instead of running the casing to TD, it was decided to land the casing at this location. In the next phase, the well was drilled with an 8 1/2” bit from 3000–4626 m. Additional 33 m were drilled with a 6” bit and coring bit to a measured depth of 4659 m. After drill- ing was finished, a 3 1/2” drill string was inserted into the well to thermally enhance the permeability in the lower part of the reservoir using cold fluid injection. During retrieval of the drill string, it was observed that the production casing was severely damaged between 2307 and 2380 m close to one of the main former feed zone for production well RN-15 (IDDP 2018). Cementing operation Prior to the cementing operation, it was known that there were loss zones near 2360 m depth and that the water table was at 600–800 m during pumping. Other older loss zones were also known. The following steps were planned for the reverse cementing operation of the 9 5/8” and the 9 7/8” casing: 3 3 • Start to batch mix 40 m at 1910 kg/m on 4 cementing trucks on location. • Fill lines and pressure test to 350 bar. 3 3 • Start cementing with two units for the first 80 m at 1.6 m /min. • Shut in valve on casing swage after 75 m are pumped. 3 3 • Reduce pumping rate to 0.5 m /min and pump 25 m . • Reduce pumping rate to 0.25 m /min. 3 3 3 • Stop pumping after 132 m and a slurry density between 1920 kg/m to 1970 kg/m . The cementing slurry is API G cement with 40 % silica flour and retarder in the first 127 ton. The laboratory analysis of the cement slurry can be found in the Appendix in Tables 4 and 5. The remaining cement slurry was without retarder. Cementing of the casing in well started around 22:45 on September 5th by pumping cement slurry through the kill-line from the top of the well. At 02:45 on the morning of September 6th, 132 m of cement slurry were already pumped into the annulus and the pressure on the casing was noted 15.6 psi (1.07 bar). After pumping the 132 m , the annulus needed to be filled up. So, after 6 h, the annulus was filled up with 17.3 m cement slurry with a density of 3 3 1970 kg/m at a rate of 0.34 m /min. The cementing operation was finished on Septem - ber 6th about 10:30 (Stefánsson et al. 2017, 2020). Cement job evaluation Two consecutive cement-bond logs (CBL) were carried out in order to monitor the evolution of the cement bonding of the casing in the well. Previous to each CBL a temperature profile was logged, both to get a temperature profile, and to see if con - ditions in the well were within the CBL-tool’s limitations (< 150 C). Additionally, a pressure tool was run simultaneously together with the temperature tool before the first CBL log was carried out. Figure 2 shows the pressure (black) and temperature profile (blue) logged on September 6th at 15:00 and on September 7th at 14:04 (red). Lipus et al. Geotherm Energy (2021) 9:5 Page 5 of 20 Fig. 2 Temperature (blue line) and pressure (black line) logged on September 6th at 15:00, Q = 0 l/s. Temperature (red line) logged on September 7th at 14:04 after completion of cementing job. Cement bond logs were carried out subsequently. Depth reference is Thor rig floor (9.0 m above surface). Figure modified from (Stefánsson et al. 2020) With respect to the pressure log, the water table inside the casing can be determined at 75 m. Below, the pressure gradually increases following the hydrostatic profile down to 2780 m. From there a slight decrease in the slope is observed where the tool probably entered the zone of cement slurry which fills the lowermost section of the well. The measured slope indicates a density of 1730 kg/m . This is lower than the potential slope for a cement slurry static profile (1940 kg/m , i.e., the used density of the pumped down cement slurry). The temperature curves exhibit a gradual increase in the uppermost 800 m with a con - stant offset of about 20–25 C. From 800–1600 m, the temperatures slightly decrease again. Below 1600 m temperatures strongly fluctuate for both logs in a similar way. Prominent cooling zones can be seen at 1600–1700 m and 2300–2400 m where most of the drilling fluid flowed out of the well during drilling. Some minor cooling zones, as part of the newly drilled section of the well (> 2500 m), are found at 2505, 2620, 2778 and 2920 m, though they rather indicate zones of higher permeability than feed zones. Lipus et al. Geotherm Energy (2021) 9:5 Page 6 of 20 The first CBL was run on September 6th at 19:40, approximately 10 h after cement - ing was finished and 22 hours after cementing was started. The second CBL was run on September 7th at 17:35. The first CBL was run to 2778 m depth. The second CBL was terminated at a depth of 2500 m because temperatures were reaching 150 C b elow that depth. The logs can be found in Fig. 3. Essentially for both cement bond logs, the recorded signal does not indicate a good cement bond. This may have different causes: e.g., the cement did not yet harden enough, or there is no cement at all. Starting with the first CBL, in the uppermost 800 m (where the outer casing is located) a good cement bond can only be seen at 780–850 m. Nevertheless, an overall low first amplitude may suggest a proper bonding except for some weak zones at 760–770, 815 and 850–860 m. It can be concluded that, throughout the most part of the uppermost casing, the cement bond is reasonably good. From 860-1400 m the cement bond is fairly good but a poor bond is observed from 895–925 m and from 970–980 m. A dramatic change in ampli- tudes can be observed below 1400 m. From there on, the cement bond (or hardening) seems generally weak. Thus, from 1400 to 2370 m, a non-uniform and poor cement bond is detected and the casing collars are readily seen in that interval. Previous tem- perature logs also infer a feed zone of good permeability around 2370 m depth. From 2370 to 2774 m the cement bonding conditions are poor and the CBL is indicative for Fig. 3 Two consecutive cement bond logs after the end of the cementing operation. CBL1: September 6th after 14–21 h curing. CBL2: September 7th 33-40 h curing. A well completion schematic with the caliper log is shown next to the CBL data. Depth reference is surface. Figure modified from (Stefánsson et al. 2020) Lipus et al. Geotherm Energy (2021) 9:5 Page 7 of 20 a nearly free pipe amplitude (unbonded casing). An exception is shown between 2710 and 2725 m indicating a good bond condition. The first CBL was limited to a depth of 2774 m. Below, the top of cementing slurry in the well was reached which is also indi- cated by the pressure log. The second CBL indicates some minor improvements in the cement bond especially in the uppermost part (< 800 m), but is still rather weak below 1400 m depth. In general, the cement bond is slightly better than in the first CBL. Com - paring the first and second CBL, some obvious changes can be observed showing more zones with improved bonding for the second CBL and casing collars not as clearly visible as in the first CBL. Cement droplet velocities in the annulus Because the water table is at around 700 m before cementing starts, we assume the fol- lowing model for the falling annular cement droplets during the reverse cementing oper- ation. We assume that the wall friction of the fluid is balancing the hydrostatic pressure gradient. This requires the fluid to maintain contact with the walls as the blobs are falling down the annulus, and for simplicity, only single-phase flow is considered while counter- current flow of air escaping upwards is neglected. If wall contact is not maintained, the blobs would accelerate indefinitely unless two-phase effects are accounted for. In the following, the annulus is assumed to be perfectly concentric. While water is a Newtonian fluid, the rheology of the cement slurries can typically be fitted to a Her - schel–Bulkley model, τ = τ + Kγ˙ ,for τ >τ , y y (1) where τ is the shear stress, γ˙ the shear rate and τ , K and n are the yield stress, consist- ency index and flow index, respectively. By modeling the annulus as a slot (Founargiota - kis et al. 2008), formulated a method to estimate the velocity of Herschel–Bulkley fluids in laminar, transitional and turbulent flow regimes. If the blobs are spanning the annular region and the axial lengths of each blob is sufficiently long, this method can be used to estimate the velocity of the falling fluid. Following their approach, the flow velocity V in laminar flow is given by m+1 m+1 dp/dL (h/2) (1 − ξ) (ξ + m + 1) V = , (2) K (m + 1)(m + 2) where m = 1/n, h is the annulus gap. The pressure gradient is the hydrostatic pressure gradient, dp/dL = ρg , (3) where ρ is the density and g the gravitational acceleration. The dimensionless shear stress ξ can then be calculated as 2τ ξ = . (4) hdp/dL To justify an assumption that the flow is laminar, the generalized Reynolds number needs to be below the transitional range, as outlined in Founargiotakis et al. 2008. They Lipus et al. Geotherm Energy (2021) 9:5 Page 8 of 20 provide a numerical procedure for predicting the velocity for transitional and turbulent flow via estimates for the friction factor, which is also used in the calculation below. Methods The fiber optic cable was installed behind the production casing (9 7/8” and 9 5/8”) of the well. A schematic of the cable design is presented in Fig. 4. The cable contains one multi-mode (MM) fiber and one single-mode (SM) fiber. The two fibers are located inside a stain-less steel tube. To reduce degradation of the optical fibers at elevated temperature, the tube is filled with dry nitrogen. Another resilient steel sheath covers the outside of the cable, giving the cable a nominal overall diameter of 6.35 mm. The maximum operating temperature is 300 C. All cables (thermocouples and fiber optic cable) were protected using three clamps per joint. Two different types of clamps were installed, two mid-joint clamps at 1/3 and 2/3 of each pipe and an over-collar clamp that also centralized the casing. Within the 13 3/8” casing, the former production cas- ing of RN-15, only a single mid-joint clamp was used. Using OTDR measurements, the fiber health was confirmed after each joint during run-in-hole. To ensure accurate temperature readings, 100 m of the fiber optic cable was fully submerged within an ice bath to correct for an eventual offset in the temperature data. A reference tem - perature gauge (VEMCO Minilog-II-T) was used to verify the constant temperature of the ice bath. Ice spray was used to locate individual traces along the optical cable. A second temperature logger was placed in the measurement container at the well site to monitor the temperature in vicinity of the DTS system. A third logger was placed on the outside of the reference bath to record the ambient temperatures. Distributed temperature sensing (DTS) Fiber optic distributed temperature sensing, or DTS, describes the technology of using each location of a glass fiber as a sensor for temperature (Hartog 1983). This is done by coupling optical pulses with a wavelength of 1064 nm into a glass fiber and measuring the backscattered light as a function of time. Given the two-way travel time and the speed of light in the glass fiber, a physical origin of the backscattering location can be determined. Different methods are available for DTS measurements. In this study, we use a system that is based on Raman backscatter. single-mode (SM) fiber multi-mode (MM) fiber strainless steel 304 Incoloy 825 (not to scale) Fig. 4 Fiber optic cable schematic Lipus et al. Geotherm Energy (2021) 9:5 Page 9 of 20 Table 2 DTS measurement overview Number of measurements 8470 Temporal resolution (s) 38 Spatial resolution (m) 0.5 Start measurement 03.09.2016 21:00 End measurement 07.09.2016 10:30 Fig. 5 Overview of downhole data which was acquired during the reverse cement job of the 9 5/8” x 9 7/8” production casing in RN-15/DEEPEGS/IDDP-2 Measurement overview The interrogation unit for distributed temperature sensing (DTS) was installed on September 3rd and continuously measured until September 7th. A total number of 8470 temperature profiles were recorded. The fiber optic cable covers the interval from surface to 832 m depth. A temperature profile was acquired every 38 s. The spa - tial resolution was 0.5 m (Table 2). Results Downhole data from pre-cementing observations and the cement pumping are shown first. After that, the analysis of the annular velocity of falling droplets is shown. Fig - ure 5 shows an overview of the DTS temperature readings gathered during the cemen- tation of the 9 7/8” and 9 5/8” production casing together with the flow rates pumped into the casing and/or the annulus as well as the readings of the downhole thermo- couples and the downhole pressure gauge. The depth interval in the DTS colormap shows the depth interval from 400–832 m. On the right, a schematic of the well with Lipus et al. Geotherm Energy (2021) 9:5 Page 10 of 20 the position of the fiber optic cable is shown together with the borehole geometry from a caliper log. Figure 6 shows a close-up of the depth interval from 550–832 m. A number of tem- perature isolines are plotted on top of the colormap to delineate subtle temperature changes along the well column. Pre‑cementing observations A pressure of 51.5 bar was measured at the downhole pressure gauge before the start of the cementation. At this time, water was injected into the annulus with a flow 3 ◦ rate of 0.672 m /min. Given a fresh water column at 20 C , the water table is esti- mated at 710 m. The pressure increase from 20:00 to 21:30 is 5.6 bar at an increas - ing injection rate of 0.567 m /min through the casing. After 21:30, both mud pumps (casing and kill-line) were stopped. This shows as a transient pressure decrease from 57 bar to 46 bar. Given the water injection rates and the downhole pressure informa- tion, an injectivity index can be calculated for the entire well column. As the pres- sure approaches 46 bar, the temperature at TC6 (932 m) and TC3 (2132 m) increases (21:50). In the same time, the DTS data show locally increased temperatures in the open-hole section at 802 m and 809 m. The DTS temperature increases to above ◦ ◦ 70 C within 5 min and then abruptly decreases back to around 20 C upon pres- sure increase at the pressure gauge. This temperature effect occurs two times, at 22:00 and 22:15. Prior to the abrupt temperature changes below 800 m, we observe a gen- tle negative temperature anomaly migrating downwards in the annulus. The pressure reading slightly increases until 22:50 by 3 bar. With the increase in pressure, no local DTS temperature peaks are present and the TC6 and TC3 readings remain relatively stable. Fig. 6 Close-up of downhole data which was acquired during the reverse cement job of the 9 5/8” and 9 7/8” casing in RN-15/DEEPEGS/IDDP-2. The pressure gauge reading is again displayed in white Lipus et al. Geotherm Energy (2021) 9:5 Page 11 of 20 Cement pumping At 22:50, the cement pumping starts. The DTS shows a cooling over the entire measure - ment length. The pressure increases to around 60 bar during the initial phase of cement pumping until 23:10. At the same time the pressure starts to increase, the temperature at TC3 (2132 m) decreases by 10 C until 23:00. The temperature then increases until the end of the initial cement pumping phase at 23:10. From 23:12 onwards, the pressure reading at 1232 m decreases to 8 bar at 23:40. In the meantime, TC3 (2132 m) shows a rapid increase in temperature from 85 to 163 C within 30 min after the main cementa- tion phase starts at 23:18. With ongoing cement pumping, the pressure increases grad- ually from 23:40 to 26:20. In the time window from hour 25:00 to 27:00, the pressure gauge reading is superimposed by a wavy pattern (indicated by the black arrow in the bottom subplot in Fig. 5). In the same time, the DTS temperature fluctuates strongly in the open-hole section at 794–832 m (see Figs. 5 and 6). From 25:50 onwards, the vertical columns marked by the DTS isolines are gradually replaced by horizontal lines which rise from the bottom upwards. As this pattern transition (white dotted line “B” in Fig. 6) reaches the shoe of the 13 3/8” casing (26:00), the depth interval from 797–805 m cools down for 6 min, before this pattern continues to rise in the cased hole interval. Below 805 m, no thermal perturbation was observed in this time period. A theoretical model of the location of the TOC was assessed by calculating annular volumes from the caliper log, casing diameters and the pump rate (white solid line “A” in Fig. 6). In the hypothetical case that no loss zones were present in the well, a total volume of approximately 141 m of cement slurry would be sufficient to fill the annu - lus to surface (neglecting any shrinkage of the cement). For comparison, line “A” was delayed along the time axis in Fig. 6 to match the location when line “B” continues to rise in the cased hole interval (26:06). The highest point of line “A” is placed at 573 m. This depth marks the location to which the well can be maximally filled in order to allow the top-up cement volume of 17.3 m (which was pumped 6 hours later) to fit into the annulus. Both in the open hole section as well as the cased hole section, the calculated TOC is significantly steeper compared to the transition line picked from DTS data. At the moment where line “A” and line “B” intersect, the pressure gauge at 1232 m reads 73 bar. After primary cement pumping stops, the pressure at the pressure gauge stabi- lizes at 57 bar. For the first 80 m , the standpipe valve was open but after closing the backpressure started to rise and ended at about 2 bar when cementing was halted. The top-up cement job after 33:00 shows a relative cooling of the upper 660 m of the well and no temperature change below, indicating an immobile cement column. The shallowest thermocouple (TC8) shows a gentle cooling during the top-up cement pumping. Annular velocities The annular fluid velocity from the DTS was determined manually with the help of contour plot (depth over time) displaying isotherms in the well (see Fig. 7). As cold water falls in the relatively warm annulus during water injection, the DTS cable cools down with time and increasing depth. During event A and B (left subplot in Fig. 7), cold water falls along the casing side with the DTS cable attached to it. This shows as a local elevation of the isotherms which propagates downwards in the well with time. Lipus et al. Geotherm Energy (2021) 9:5 Page 12 of 20 Fig. 7 Picking of droplet velocity from contour plot for event A and event B (water) and for event E and event F (cement) Table 3 Fluid speed in annulus (see events marked in Fig. 5) Event Depth start (m) Depth Time start (s) Time end (s) Velocity (m/s) Comment bottom (m) A 0 829 320 690 2.25 Water B 0 829 1850 2305 1.83 C 300 829 10405 10450 11.84 Cement D 0 829 12720 12800 10.41 E 0 700 14700 14770 10 F 0 700 15760 15820 11.66 Displayed times relative to 22:00 of September 5th The same approach was chosen to determine the fluid velocity of the cement slurry (right subplot in Fig. 7). The contour plots allow only for a rough velocity estimate because the temperature perturbation of falling cement is rather low and the fluid velocity exceeds the DTS sampling rate of 38 s. The velocity of the interpreted falling water in the DTS data in Table 3 is an order- of-magnitude greater than the expected bulk flow velocity as calculated from the pump rate, and for the falling cement it is greater by almost two orders-of-magnitude. By adopting the model outlined before, the resulting velocity for the falling cement blobs when assuming laminar flow is 12 m/s. However, this suggests that the flow is well into the turbulent regime, which adjusts the velocity to 6.4 m/s by following the procedure for transitional and turbulent flow. This places the predicted velocity for the cement blobs within the observed order-of-magnitude but lower by about 40% . It should be noted that the models used, assumes the fluid to be continuous, which in our case is a simplification that may result in inaccuracies. For the falling water, the approach predicts 10 m/s, which is higher than for cement and an order-of-mag- nitude higher than observed. On the other hand, since the falling water is not asso- ciated with DTS temperature fluctuations or the wavy pattern superimposed on the Lipus et al. Geotherm Energy (2021) 9:5 Page 13 of 20 pressure gauge as for the falling cement, it is not unlikely that these are two different phenomena and that the modeling approach is less suitable for water. Discussion In the RN-15/DEEPEGS/IDDP-2 well, a reverse cementing operation was performed for the 9 5/8” and 9 7/8” casing interval. The downhole monitoring data from distributed temperature sensing (DTS), thermocouples and pressure data allow an analysis of the downhole dynamics during the cementing operation. The data provides insight on the (1) assessment of the downhole formation injectivity, (2) falling cement droplets and (3) the impact on the outcome of the cementing operation. Annular injectivity and static formation pressure During the pre-cementing phase, the annular temperature was monitored in combina- tion with the flow rate of the mud pumps into the casing as well as the annulus of the casing. Given the water injection rates and the downhole pressure information, an injec- tivity index can be calculated for the entire annular well column if hydrostatic equilib- rium is assumed between casing and annulus. The assumption of hydrostatic equilibrium is supported by the fact that an open shoe is installed to allow for reverse cementing. The pressure increase p from 20:00 to 21:30 is 5.6 bar at an additional flow rate of 0.567 m / q˙ 3 m min (34.02 m /h). This gives an estimated injectivity index of II = = 6.07 for h·bar the open-hole section below the 13 3/8” casing shoe. Given this injectivity index and the 3 3 initial injection flow rate of 0.672 m /min (40.32 m /h), we can use the same equation to solve for pressure (5.6 bar). With this, we estimate the static formation pressure to be at p = 51.5 bar − 5.6 bar = 44.9 bar and the natural water table at 773 m in the well, assum- ing a fresh water column of 20 C. When the mud pumps were stopped at about 21:45, the annular pressure decreased to about 46 bar at the pressure gauge, close to the estimated formation pressure. The annu - lar temperature at a depth of 800 m increased significantly. As the temperature of the annular fluid is lower than the formation temperature, the density of the annular fluid is higher compared to the density of the formation fluid. Schematic profiles of pressure and fluid density for the depth interval above the pressure gauge before and after fluid injection are depicted in Fig. 8. The cooler and denser fluid column might cause under - balanced conditions in the annulus and possibly a local inflow of warmer formation fluid. Contemptuously, at the depth of thermocouples TC6 (932 m) and TC3 (2132 m) a similar temperature increase occurred, indicating a thermal anomaly likely caused by an interaction with formation fluid. From the observed temperature anomalies, we inter - pret underbalanced conditions down to TC3. In the time window between the stopping of both mud pumps and start of cementation, the pressure slightly increases by 3 bar indicating a heating up of annular fluid and an increase of the water table in the well. However, we observe a negative temperature anomaly, indicating a fluid flow down the annulus that cools the well along the entire monitored well column and which is espe- cially prominent below 800 m. The volume is unknown as there is no flow documented in the rig log. This might be the more likely reason for the pressure increase prior to the onset of cement pumping at 22:50. The complete picture of the underlying hydraulic process cannot be resolved by the limited amount of information available. Lipus et al. Geotherm Energy (2021) 9:5 Page 14 of 20 Fig. 8 Simplified representation of the pressure gradient in the well at different phases of the cementing operation Analysis of cement pumping and velocity of free‑falling cement After the initial phase on cement pumping into the annulus (from 22:50 to 23:10), the annular pressure drops to 8 bar at the pressure gauge, indicating that (a) the flow rate of the cement is not sufficient to build up a homogeneous cement plug injected into the annulus and (b) cement is flowing to below the pressure gauge. During pumping into the annulus, the cement column drives down the fluid column and counteracts the reservoir pressure. Due to its higher density, the TOC falls below the initial water table creating underbalanced conditions in the upper part of the annulus indicated by the pressure drop seen on the pressure gauge. A pressure reading of 8 bar at 1232 m can mean that the TOC is located at 1190 m at 23:50, given a cement slurry density of 1940 kg/m and assuming all annular water is displaced below the depth of the pressure gauge. The subsequent pressure increase (until 26:40) can be explained by a rising TOC above the pressure gauge. At the estimated static formation pressure of 44.9 bar, the TOC would be located at a depth of 996 m given the initial slurry density (see Fig. 8). In the DTS data, we observed a transition from a vertical pattern to a horizontal pat- tern which rises from bottom to top from 25:45 onwards (line “B” in Fig. 6). We interpret this interface as the location of the TOC. According to our understanding, the vertical pattern was created by cold falling individual cement blobs into the annulus while in the mean time, relatively hot formation fluid from the surrounding formation intersects the well. Once a continuous cement column formed in the well, a horizontal pattern in the DTS formed because no falling cement blobs occurred any more and the feed zone below the 13 3/8” casing shoe were partly filled with cement slurry. If we follow our interpre - tation that line “B” marks the location of the TOC, we can say that the TOC reaches the shoe of the 13 3/8” at 26:00 and that it continues to rise in the cased hole interval from 26:06 onwards. In the 6 min in between, we observe a cooling in the depth interval from 795 to 806 m. This might indicate that a part of the cold cement accumulated in the for - mation just below the shoe of the 13 3/8” casing. Giving the time (6 min) and the cement 3 3 pumping rate (0.3 m /min), an estimated volume of 1.8 m might be accumulated and possibly lost into the formation. At 26:06, the average cement slurry density ρ in the Lipus et al. Geotherm Energy (2021) 9:5 Page 15 of 20 interval between the casing shoe and the pressure gauge can be calculated according to 5 2 ρ = p/(g · h) = (73 · 10 Pa)/(9.81 m/s · (1232 m − 794 m)) = 1699 kg/m . The pump density from the cement report states an initial cement density of 1940 kg/m . The density differences could mean that the cement slurry was diluted by 34 % of formation water. The calculated density of 1699 kg/m is close to the density observed inside the casing during pressure and temperature logging after the end of cementation (1730 kg/ m ). Following line “B” to the end of the primary cementing phase (26:45), we interpret that a depth of about 700 m is reached in the well. In comparison, the thermal perturbation observed in the DTS data during the top-up cement job was observed down to 660 m and the volumetric calculation shows that a top-up cement slurry volume of 17.3 m should have allowed to fill an annular volume from rig floor down to 573 m. The dis - crepancy between 700 and 660 m can be partly explained by remaining cement slurry which is still sinking in the annulus just after stopping of the pumps and the inaccu- rate estimate of the TOC as the thermal anomaly caused by the cement is small. This is also shown by the DTS measurements in the time from 27:00 until the top-up cement- ing operation. We measured latent heat from cement hydration up to a depth of around 300 m in the well, indicating that some cement has set at the casing wall and at the fiber optic cable. This is also the reason, why the 17.3 m of cement slurry are sufficient to top- up the annulus while measuring the thermal perturbation down to 660 m. Contradictory to the visual TOC assessment of line “B” is the fact, that the calculated TOC (line “A”) is much steeper for both the open-hole and cased hole section. In the open- hole section, line “A” shows 4.5 m/min (line “B” is 2.0 m/min) and in the cased hole section, line “A” shows 9.5 m/min (line “B” is 2.7 m/min). From our analysis, we feel confident that the TOC increases in the cased hole from 26:06 (intersection of both line “A” and “B”). Despite the artificial temporal lag of line “A” to match line “B” at the 13 3/8” casing shoe, the projected final TOC after primary cementing is still too high for line “A” to allow for the top-up cement to fit the annulus. This indicates that cement slurry continues to sag into the open-hole section in the period where the TOC increases in the cased hole section. After the pressure reaches a maximum of 75 bar when the TOC reaches the shoe of the 13 3/8” casing, it decreases back to 57 bar. The gradual decrease coincides with a temperature increase within the annulus which is likely caused by the latent heat of hydration of the cement. The pressure decrease can either be attributed to the setting of the cement or a gradual sagging of part of the cement to below the pressure gauge, opening a hydraulic connection to the formation. While the CBL shows homogene- ous values between the shoe of the 13 3/8” casing and the pressure gauge, it is assumed that the cement does not further sag down in the annulus. This is also supported by the fact that there is no vertical temperature perturbation observed below the shoe after 26:00. The theory of free-falling cement blobs is further confirmed by the comparison of calculated velocity with the velocities that were interpreted from the DTS data. The fiber optic cable is in direct contact with the pumped fluids in the well annulus. The cable mantle material is stain-less steel with a high thermal conductivity. Therefore, we assume an almost instantaneous thermal effect on the temperature sensing glass fiber upon a temperature change caused by the moving fluid front. The simple modeling approach predicted a velocity within the observed order-of-magnitude, which is almost Lipus et al. Geotherm Energy (2021) 9:5 Page 16 of 20 two order-of-magnitudes greater than what would be expected from the injection of a homogeneous cement plug in the annulus. Impact on the outcome of the cementing operation From the CBL log evaluation we conclude a low cement bond. We interpret the low cement bond as a result of a non-optimal placement of cement slurry (e.g., dilution, free-falling of cement blobs). Analyzing the DTS temperature data acquired during the cementing operation in combination with the thermo-electrical sensors allowed to reveal some of the governing hydraulic processes. Not only did the analysis show that the limited flow rate led to a free fall of cement blobs within the annulus, it could also be shown that dilution of cement occurred that led to a dilution of the cement slurry which corresponds not only to a low CBL reading but also to the analysis of the TOC together with the pressure readings within the annulus. It is important to mention that there is a risk that the pressure gauge was malfunction- ing and that the recorded values do not display the actual downhole pressure conditions. The calibration or restriction to protect the gauge could also be affecting the read-out. The damage on TC3 which was observed during run-in-hole of the casing also imposes the risk that the temperature reading is incorrect. Conclusions For the first time, a reverse cementing job in a high-temperature geothermal well could be monitored and analyzed using a combination of permanently installed distributed fiber optic and electronic sensors as well as conventional well logging equipment. Data from the permanently installed sensors were used to interpret the hydraulic situation during cement pumping and the initial phase of cement setting. Based on this analysis, CBL data from the uppermost 832 m could be interpreted with more confidence. The low CBL reading in the depth interval below the 13 3/8” casing shoe is most likely related to the excessive amount of retarder. The amount of retarder was calculated based on the BHCT that did not consider the well cooling due to the loss of circulation. Never- theless, the low CBL reading can additionally be attributed to the observed injectivity of water and cement into the formation which is pointed out by our analysis. We observe an inflow of fluid (or steam) from the formation at times where the pressure gauge read - ing drops below the static formation pressure both in the DTS temperatures and ther- mocouple data. The inflow of formation fluid dilutes the cement slurry and reduces its density. From our calculation based on pressure gauge and DTS data, each cubic meter of cement slurry was diluted with an additional 34% of water in the depth interval from 13 3/8” casing shoe (794 m below surface) and the pressure gauge (1232 m below sur- face). Moreover, we find evidence from the DTS temperatures, that at least 1.8 m of cement slurry was pumped into the formation just below the shoe of the 13 3/8” casing. This finding is partly based on a cooling effect below the casing shoe and partly based on the change of the general DTS pattern for depth intervals where cement is falling (verti- cal pattern) and depth intervals where the cement begins to settle (horizontal pattern). We interpret this pattern transition in the DTS data as the TOC. Complementary to conventional monitoring data from cementing operations, fiber optic distributed temperature data can give valuable information about downhole Lipus et al. Geotherm Energy (2021) 9:5 Page 17 of 20 processes in real-time. From the analysis done in this study, we find evidence that cur - ing loss zones concurrently with cementing a long casing string might cause undesired results. The current practice of drilling blind, i.e., drilling with total circulation loss, as a standard procedure in high-temperature geothermal drilling should be reconsidered for such long cementing intervals. While drilling an interval with total circulation loss might be acceptable if no casing is planned to be cemented, drilling blind for an interval that should be cased and fully cemented needs to be carefully evaluated. The experience gathered in Iceland from this work and from the longer-term emphasis on research targeting deeper understanding of the deep roots of geothermal energy har- nessing, has been funded through both national efforts (GEORG 2020) and with support from European framework research programmes (Fridleifsson et al. 2017). This type of joint collaborative efforts involving research organizations, energy industry and energy agencies that collate and openly share results are needed to increase predictability and successful harnessing of geothermal energy. Abbreviations BHCT:: Bottom-hole-circulation-temperature; BHST:: Bottom-hole-static-temperature; bwoc:: By-weight-of-cement; Cal.x/ Cal.y:: Caliper log; CBL:: Cement-bond-log; DTS:: Distributed temperature sensing; IDDP:: Iceland Deep Drilling Project; II:: Injectivity index; KOP:: Kick-off-point; MD:: Measured depth; MM:: Multi-mode (fiber); OTDR:: Optical time domain reflectometry; PG:: Pressure gauge; pT logging:: Pressure and temperature logging; SM:: Single-mode (fiber); sk:: Sack, unit of measure for Portland cement; TC:: Thermocouple; TOC:: Top-of-cement; TD:: Total depth, maximum depth reached in the well. Acknowledgements The authors are thankful to Gudmundur Ómar Fridleifsson, Ómar Sigurdsson, as well as the staff from HS Orka for provid- ing access to the well and well data and support during and after the field campaign. The authors are also thankful to the staff of Iceland Drilling for the support during the field campaign. This work would not have been possible without the continuous support from our colleagues at ÍSOR and GEORG Geothermal Research cluster. The authors would also like to thank Adrian Desmarais from Petrospec Engineering for the fruitful discussions on the data from the electric downhole gauges. The authors are also grateful to the reviewers for their thorough comments and suggestions to this manuscript. Authors’ contributions MPL and TR planned the measurement campaign, collected the DTS data, interpreted the DTS data and analyzed the annular injectivity and static formation pressure. TBW performed the CBL log interpretation. SK provided the analysis on the cement droplet velocities in the annulus. AS provided data for the well drilling and completion and the cementing procedure. SGB managed the DEEPEGS project and commented on the final draft. All authors read and approved the final manuscript. Funding Open Access funding enabled and organized by Projekt DEAL. The research leading to these results has received fund- ing from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 654662 (Project GeoWell) and grant agreement No. 690771 (DEEPEGS). Also, this project has received funding from the European Union’s Seventh Programme for research, technological development and demonstration under grant agreement No: Data availability statement The acquired fiber optic DTS data are freely available under the Creative Commons Attribution International 4.0 Licence (CC BY 4.0) by GFZ Data Services (Lipus and Reinsch 2020). Data from the temperature and pressure gauges are also freely available (Eggertsson and Stefánsson 2020). Competing interests The authors declare that they have no competing interests. Author details 1 2 GFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany. Fraunhofer IEG, Fraunhofer Research Institution for Energy Infrastructures and Geothermal Systems IEG, Lennershofstrasse 140, 44801 Bochum, 3 4 Germany. Iceland GeoSurvey (ÍSOR), Grensásvegur 9, 108 Reykjavík, Iceland. NORCE Norwegian Research Centre AS, 5 6 Prof. Olav Hanssensvei 15, 4021 Stavanger, Norway. HS Orka hf, Svartsengi, 240 Grindavík, Iceland. GEORG Geothermal Research Cluster, Grensásvegur 9, 108 Reykjavík, Iceland. Appendix See Tables 4, 5 and 6 Lipus et al. Geotherm Energy (2021) 9:5 Page 18 of 20 Table 4 Laboratory analysis of cement slurry (Part 1/2) Chemical Concentration Unit Chemical composition Cement 94.00 lb/sk Silica flour 40.0 % bwoc Retarder 0.50 % bwoc Calculations Surface density 16.2 lb/gal Downhole density 16.2 lb/gal Sack weight 94.00 lb/sk Dry weight 132.07 lb/sk Dry yield 1.55 ft3/sk Water 5.87 gal/sk Total Fluids 5.87 gal/sk Surface Yield 1.49 ft3/sk Downhole Yield 1.49 ft3/sk Conditions Time to BHCT 33 min BHCT 126.6 BHST 260 C Initial pressure 459 psi Final pressur 3632 psi Simulated depth 6189 ft See Table 6 for unit conversion Table 5 Laboratory analysis of cement slurry (Part 2/2) Thickening time Free fluids Bearden units of consistency Time Surface (21 C) 0 (bc) 19 bc 0:00 Downhole (82 C) 0 40 bc 4:32:52 50 bc 4:34:30 60 bc 4:35:01 70 bc 4:35:22 99 bc 4:38:19 ◦ ◦ Surface rheology (21 C) Downhole rheology (82 C) RPM Dial reading RPM Dial reading 600 302.2 600 183.3 300 189.3 300 113.8 200 135.6 200 87.4 100 76.6 100 54.8 60 50.9 60 42.5 30 33.0 30 36.3 6 18.2 6 14.3 3 32.3 3 15.1 PV 144.1 cP PV 83.7 cP 2 2 YP YP 28.5 lbf/100 ft 24.6 lbf/100 ft The cementing contractor used a Fann Model 3530A to measure the rheological parameter Lipus et al. Geotherm Energy (2021) 9:5 Page 19 of 20 Table 6 Conversion table from field unit to SI unit Length 1 ft 0.3048 m Mass 1 lb 0.453 kg Volume 100 gal 0.3785 m 3 3 100 ft 2.83 m Density 1 lb/gal 119.8 kg/m Pressure 1 psi 6894 Pa 100 Pa 2.09 lbf/ft Viscosity 1 cP 1 mPa s Received: 20 August 2020 Accepted: 22 January 2021 References Allouche M, Guillot D, Hayman AJ, Butsch RJ, Morris CW. Well Cementing 2nd Edition. Schlumberger. 2006:562–563. Teodoriu C, Falcone G. Comparing completion design in hydrocarbon and geothermal wells: The need to evaluate the integrity of casing connections subject to thermal stresses. Geothermics. 2009;38:238–46. https ://doi.org/10.1016/j. geoth ermic s.2008.11.006 Salim P, Amani M. Special considerations in cementing high pressure high temperature wells. International Journal of Engineering and Applied Sciences. 2012;1.https ://doi.org/10.1007/978-94-009-9059-3 Bertani R, Büsing H, Buske S, Dini A, Hjelstuen M, Luchini M, et al. The first results of the Descramble project. In: Proceed- ings, 43rd workshop on geothermal reservoir engineering, Stanford University, Stanford, California. 2018. Pálsson B, Hólmgeirsson S, Gudmundsson A, Bóasson HA, Ingason K, Sverrisson H, et al. Drilling of the well IDDP-1. Geo- thermics. 2014. https ://doi.org/10.1016/j.geoth ermic s.2013.08.010 Fridleifsson GÓ, Bogason S, Ingolfsson H, Vergnes P, Thorbjörnsson I, Peter-Borie M, et al. Deployment of deep enhanced geothermal systems for sustainable energy business. Strasbourg, France: European Geothermal Congress; 2016. Ravi KM, Sutton DL. New rheological correlation for cement slurries as a function of temperature. SPE Annual Technical Conference and Exhibition. 1990.https ://doi.org/10.2118/20449 -MS Bertini G, Giovannoni A, Stefani GC, Gianelli G, Puxeddu M, Squarci P. Deep Exploration in Larderello Field: Sasso 22 Drill- ing Venture. Advances in European Geothermal Research. 1980.https ://doi.org/10.1007/978-94-009-9059-3 Fridleifsson GÓ, Elders WA, Albertsson A. The concept of the Iceland deep drilling project. Geothermics. 2014. https ://doi. org/10.1016/j.geoth ermic s.2013.03.004 Reinsch T, Dobson P, Huenges E, Poletto H, Sanjuan B. Utilizing supercritical geothermal systems: a review of past ven- tures and ongoing research activities. Geothermal Energy. 2017. https ://doi.org/10.1186/s4051 7-017-0075-y Hurtig E, Grosswig S, Jobmann M, Kühn K, Marschall P. Fibre-optic temperature measurements in shallow boreholes: experimental application for fluid logging. Geothermics. 1994;23. https ://doi.org/10.1016/0375-6505(94)90030 -2 Förster A, Schrötter J, Merriam DF, Blackwell DD. Application of optical-fiber temperature logging - an example in a sedimentary environment. Geophysics. 1997;62(4).https ://doi.org/10.1190/1.14442 11 Henninges J, Huenges E, Burkhardt H. In situ thermal conductivity of gas-hydrate-bearing sediments of the Mallik 5L–38 well. Journal of Geophysical Research. 2005;110. https ://doi.org/10.1029/2005J B0037 34 Carnahan BD, Clanton RW, Koehler KD, Harkins GO, Williams GR. Sensor Highway, Fiber Optic Temperature Monitoring Technology. SPE Western Regional Meeting held in Anchorage. 1999. Patterson JR, Cardiff M, Coleman T, Wang H, Feigl KL, Akerley J, et al. Geothermal reservoir characterization using distributed temperature sensing at Brady Geothermal Field. The Leading Edge: Nevada; 2017. https ://doi. org/10.1190/tle36 12102 4a1.1 Pearce JG, Rambow FHK, Shroyer W, Huckabee P, deJongh H, Dria DE, et al. High Resolution, Real-Time Casing Strain Imaging for Reservoir and Well Integrity Monitoring: Demonstration of Monitoring Capability in a Field Installa- tion. SPE Annual Technical Conference and Exhibition. 2009. https ://doi.org/10.2118/12493 2-MS Earles DM, Stoesz CW, Surveyor N, Pearce JG, DeJongh HA. Fiber Optic Strain Sensing at the Sand Face Enables Real- Time Flow Monitoring and Compaction Mitigation in Openhole Applications. SPE Annual Technical Conference and Exhibition held in Denver, Colorado, USA. 2011. Koelman JMVA, Lopez JL, Potters JHHM. Fiber Optic Technology for Reservoir Surveillance. Conference Proceedings, IPTC 2012:International Petroleum Technology Conference. 2012. Reinsch T, Henninges J, Ásmundsson R. Thermal, mechanical and chemical influences on the performance of optical fibres for distributed temperature sensing in a hot geothermal well. Environmental Earth Science. 2013;70. https :// doi.org/10.1007/s1266 5-013-2248-8 Bücker C, Grosswig S. Distributed temperature sensing in the oil and gas industry - insights and perspectives. Oil Gas European Magazine. 2017;43: Lipus M, Reinsch T, Schmidt-Hattenberger C, Henninges J, Reich M. Gravel Pack Monitoring With a Strain Sensing Fiber Optic Cable. Oil Gas European Magazine 4. 2018. https ://doi.org/10.19225 /18120 2 Lipus et al. Geotherm Energy (2021) 9:5 Page 20 of 20 Raab T, Reinsch T, Cifuentes SRA, Henninges J. Real-Time Well-Integrity Monitoring Using Fiber-Optic Distributed Acoustic Sensing. SPE Journal. 2019;24(05). https ://doi.org/10.2118/19567 8-PA Fridleifsson GÓ, Elders WA, Zierenberg RA, Stefánsson A, Fowler APG, Weisenberger TB, et al. The Iceland Deep Drilling Project 4.5km deep well, IDDP-2, in the seawater-recharged Reykjanes geothermal field in SW Iceland has success- fully reached its supercritical target. Scientific Drilling. 2017;23:1–12. https ://doi.org/10.5194/sd-23-1-2017 Jónsson SS, Sigurgeirsson MÁ, Sigurdsson Ó, Ingólfsson H. Reykjanes - Hola RN-15. 3. áfangi: Borun vinnsluhluta frá 804 m í 2507 m dypi. Iceland GeoSurvey, ÍSOR-2010/050 (in Icelandic). Report. 2010. Stefánsson A, Gíslason T, Sigurdsson Ó, Fridleifsson GO. The Drilling of RN-15/IDDP-2 Research Well at Reykjanes in SW Iceland. GRC Transactions. 2017;41. IDDP. IDDP-2 Way Forward Workshop (WFW ), Report of Workshop No. 11 of the Iceland Deep Drilling Project, 23 March 2018. https://iddpis/wp-content/uploads/2018/04/Way-Forward-Workshop-20-21-March-2018-SAGA-REPORT-No-11pdf. Stefánsson A, Fridleifsson GÓ, Sigurdsson Ó, Gíslason T. The IDDP-2 DEEPEGS Drilling Experience and Lesson Learned. Proceedings World Geothermal Congress 2020. 2020. Founargiotakis K, Kellessidis VC, Maglione R. Laminar, Transitional and Turbulent Flow of Herschel-Bulkley Fluids in Con- centric Annulus. The Canadian Journal of Chemical Engineering. 2008;86. https ://doi.org/10.1002/cjce.20074 Hartog AH. A Distributed Temperature Sensor Based on Liquid-Core Optical Fibers. Journal of Lightwave Technology. 1983. https ://doi.org/10.1109/JLT.1983.10721 46 GEORG. Legacy of project results, publications;. visited by Internet (https://georgclusteris/publications/), October 2020. 2020. Lipus M, Reinsch T. Data from distributed temperature sensing (DTS) measured along a fiber optic cable permanently installed behind casing in well RN-15/DEEPEGS/IDDP-2. GFZ Data Services: Iceland; 2020. https ://doi.org/10.5880/ GFZ.6.2.2018.0 Eggertsson GH, Stefánsson A. Temperature and pressure data from permanently installed sensors behind casing in well RN-15/DEEPEGS/IDDP-2. GFZ Data Services: Iceland; 2020. https ://doi.org/10.5880/GFZ.4.8.2020.004 Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Geothermal Energy – Springer Journals
Published: Feb 18, 2021