Dynamic motion monitoring of a 3.6 km long steel rod in a borehole during cold-water injection with 1 distributed fiber-optic sensing

16 Fiber-optic distributed acoustic sensing (DAS) data finds many applications in wellbore monitoring such as e.g. 17 flow monitoring, formation evaluation, and well integrity studies. For horizontal or highly deviated wells, wellbore 18 fiber-optic installations can be conducted by mounting the sensing cable to a rigid structure (casing/tubing) which 19 allows for a controlled landing of the cable. We analyze a cold-water injection phase in a geothermal well with a 20 3.6 km long fiber-optic installation mounted to a ¾” sucker-rod by using both DAS and distributed temperature 21 sensing (DTS) data. During cold-water injection, we observe distinct vibrational events (shock waves) which 22 originate in the reservoir interval and migrate up-and downwards. We use temperature differences from the DTS 23 data to determine the theoretical thermal contraction and integrated DAS data to estimate the actual deformation 24 of the rod construction. The results suggest that the rod experiences thermal stresses along the installation length 25 – partly in the compressional and partly in the extensional regime. We find strong evidence that the observed 26 vibrational events originate from the release of the thermal stresses when the friction of the rod against the borehole


Introduction
Fiber-optic distributed sensing in bore-hole applications has gained a lot of attention in the recent years.Distributed temperature sensing (DTS) has been used to assess rock thermal properties and locations of water-bearing fractures (e.g.Hurtig, 1994, Förster, 1997).DTS was used to perform cement job evaluations and wellbore integrity analysis during and after production tests (e.g.Pearce et al., 2009, Bücker andGroßwig, 2017).The performance of a borehole heat exchanger was monitored with DTS to evaluate the heat input along the wellbore and to measure the regeneration time after a heat extraction period (Storch et al., 2010).While DTS has found its way as a standard tool for wellbore monitoring over the last two decades, the utilization of distributed acoustic sensing (DAS) is still subject to many research questions.Johannessen et al., 2012 introduced the potential and capabilities for acoustic in-well monitoring applications based on DAS systems which range from e.g.flow measurements, sand detection, gas breakthrough, leak detection to vertical seismic profiling (VSP).Daley et al., 2013, Mateeva et al., 2014, Harris et al., 2016, Daley et al., 2016and Henninges et al., 2021 compare traditional geophone with DAS recordings acquired during a vertical seismic profiling campaign (VSP).Götz et al., 2018 report on a multi-well VSP campaign at a carbon dioxide storage site by using only one single DAS interrogator.Finfer et al., 2014 performed an experiment to study DAS applications for turbulent single-phase water flow monitoring in a steel pipe.Bruno et al., 2018 investigate the potential to use downhole DAS data for cross-hole monitoring between two adjacent wells by inducing low frequency pressure pulses to detect high conductivity zones by measuring characteristic vertical strain patterns.Naldrett et al., 2018 compare fiber-optic technology to traditional production logging tools and provides field data examples of flow monitoring based on both DTS and DAS with wirelinetype installations.Ghahfarokhi et al., 2019 analyze an extensive data set including borehole geophone and DAS during hydraulic fracturing (cable behind casing) to study micro-seismicity and low frequency events in the reservoir sandstone and silt caprock by injected CO2 can be quantitively evaluated using static distributed strain sensing over periods of 42 hours (cable behind casing).Zhang et al., 2020 provide an attempt to use distributed strain sensing to monitor elastic rock deformation during borehole aquifer testing to derive hydraulic parameter information.Miller et at., 2018 compare DTS and time-integrated DAS recordings from a borehole and finds a correlation between DTS recordings and very low frequent DAS strain recordings.In their work, they report on repeating "slip events" seen in the DAS data as short and confined vibrational events upon temperature changes in the well.
The study at hand observes similar "slip events" and shows their causal connection to the thermo-mechanical response of the borehole construction to water flow therein.
Installing a fiber-optic cable in a borehole requires specialized equipment.Depending on the aim of the fiber-optic monitoring campaign, different cable installation types are possible.One way is to permanently install the cable by mounting it to the outside of a casing and run it together with the casing into the well and cement it in place (e.g.Henninges et al. 2005, Reinsch et al., 2013, Lipus et al., 2021).A cemented fiber-optic cable generally provides a thorough mechanical coupling to the surrounding structure which is favorable for DAS data quality.
Due to its placement behind the casing, the fibers do not interfere with well operations and monitoring of the well can be performed at any time.However, the cemented annulus of a well is a crucial secondary barrier element for well integrity which is compromised by the installation of a fiber-optic cable.A fluid pathway could potentially be created along the cable.cases where the well completion design includes liner elements, a permanent cable installation behind casing to the end of the well is technically not possible, or at least, very challenging.In such cases, other installation types are available.A semi-permanent installation along e.g. a production tubing or a temporary installation via a wireline cable or coiled tubing allow cable placements inside the borehole after drilling is finished.Munn et al., 2017 present a field test of a novel "flexible borehole coupling technique" that allows deploying fiber-optic cables in boreholes after completion has finished with an improved mechanical coupling compared to lose installed fiber-optic cables.Due to physical constrains, this technology is best suited for shallow boreholes (< 425 m).Becker et al. 2017 provide an analysis of borehole fracture displacements using such kind of cable coupling technique.Another method to land a fiber-optic cable into a well is by mounting it to a rigid rod (e.g. a pump sucker-rod).The stiff sucker-rod acts as a centralizer and guides the flexible fiber preventing it from coiling up.Such type of installation is especially advantageous when the cable should be placed in a deep and deviated well.
To utilize acquired fiber-optic data from a free-hanging/free-lying rod with the highest possible confidence, it is important to understand the behavior of such a long and stiff structure inside a well.Heating and/or cooling of the well will lead to thermal stresses in the material which potentially result in contraction or expansion of the suckerrod and fiber-optic cable construction.As the fiber-optic cable is firmly attached to the rods, these dynamics influence the distributed strain and temperature sensing.From DTS monitoring, Schölderle et al., 2021 found that measurement equipment in the previously described setting does indeed contract upon the injection of cold water and that the points spatially sampled by the distributed sensing change their position.Besides a detailed analysis based on DAS and DTS data of the rod's dynamics in response to temperature changes during a cold-water injection, we show that the resulting thermal stresses are released by the observed vibrational events thus indicating stick-slip like behavior of the rod-borehole wall compound.

Well description and cable installation
The fiber-optic cable is installed within a production well at the geothermal site Schäftlarnstraße in Munich, Germany.A detailed description of the geothermal site and the cable installation procedure is presented in Schölderle et al., 2021.The well was completed with a 20" anchor casing, a 13 3/8", a 9 5/8" liner and a perforated 7" production liner.An overview of the landing depths is presented in Table 1.The design of the borehole completion is schematically shown in Figure 3 (right subplot).The well is vertical to a depth of 250 m.Below 250 m, the well is slightly inclined to 4° down to a depth of 879 m TVD (880 m MD).A number of kick-off-points (KOP) are located along the well path.These are also listed in Table 2.In the result section, a survey shows the well path.From a flow-meter log it is known, that the most prominent feed zone in the well is just below the transition from 9 5/8" liner to 7" perforated liner in the depth interval between 2825 -2835 m MD.The downhole fiber-optic cable is a tubing-encapsulated-fiber (TEF) that contains two multi-mode and two singlemode fibers.In this fiber-in-metal-tube (FIMT) construction, the sensing fibers are embedded in gel and placed in a metal tube.At elevated strain levels, the gel deforms plastically and allows for a relative motion between fiber and cable.Also, creep between cable construction and optical fibers can occur.Strain measurements with such a type of cable are typically applicable for dynamic strain changes (high frequencies) and low deformations (Reinsch et.al, 2017).For longer periods and higher deformations, fiber-optic strain sensing with FIMT cables is still possible but it becomes less localized due to deformation of the material.A laboratory experiment on the relative motion between cable structure and optical fiber in a FIMT cable at higher mechanical stress over time is presented in literature (Lipus et al., 2018).The cable has a total nominal diameter of 0.43 inch (1.1 cm) and the cable mantle is made of polypropylene.The cable was landed in the well after drilling was finished.To safely and effectively navigate the placement of the fiber-optic cable down to the end of the almost 3.6 km long well, the cable was strapped to steel rods (sucker rods) which were installed in the well together with the cable.The steel sucker rod also helps to retrieve the cable from the bore-hole when needed.Due to the high deviations in the well at depth, the cable needs to be gently pushed into the well.Therefore, the rigid sucker rod is used for the installation instead of a wireline-type installation.The final landing depth of the sucker rod construction is 3691 m (MD).Figure 1 depicts the configuration of the sucker rod/fiber-optic compound.Together with a number of cross-over elements and the final landing joint, more than 400 of individual sucker rod elements were installed in the well.In the following, we refer to the sucker rod / fiber-optic cable construction as "the rod".The depth reference for the DTS (spot warming) and DAS (tap test) are set to surface.
A fiber-optic pressure/temperature (p/T) gauge was installed with the rod and positioned at the top of the reservoir section at 2755 m (MD).

Monitoring campaign
The data shown in this study was measured before and during a cold-water injection test in a geothermal well.
Before the start of fluid injection, the well was shut-in for 29 days, so that the initial temperature profile is close to the natural geothermal gradient of the Bavarian Molasse basin (see Schölderle et al., 2021).The temperature at the well head was 17 °C and increasing up to 110 °C at the bottom of the well just before the injection start (see profile "00:48" in left panel in Figure 4).Cold-water fluid injection started on January 23, 2020 at 00:56 by pumping water through the wellhead which leads to a cooling of the well.With an initial water table at a depth of 170m below surface, water was injected from the surface without pressure built-up at the wellhead.The cold-water injection was maintained for 24h at a flow rate of 83 m3/h.In this study, we analyze the transient phase of well temperature change for the first 72 minutes of cold-water fluid injection.
Figure 1: Down-hole cable configuration of the sucker rod with a centralizer (black) and the fiber-optic cable (yellow)

Data Analysis
The analysis in this study is based on the comparison of strain derived from fiber-optic distributed temperature sensing (DTS) on the one hand and distributed acoustic sensing (DAS) on the other.

Derivation of strain from DTS
DTS uses each location of a glass fiber as a sensor for temperature (Hartog, 1983, Hartog andGamble, 1991).
This is achieved by coupling laser-light pulses into a glass fiber and analyzing the Raman spectrum of the backscattered light whose origin along the fiber is determined by the two-way travel time of the light.In this study, we use a system based on Raman backscatter.Temperature profiles were acquired every 10 minutes with a spatial sampling of 0.25 m.Detailed information about the performance of the fiber-optic system and the calibration procedure are presented in Schölderle et al., 2021.
We calculate the change in temperature from DTS at the start of fluid injection and the profile later during fluid injection.From DTS measurements we may predict thermo-mechanical deformation according to where αrod is the thermal expansion coefficient and ΔT(x) is the temperature difference at two subsequent points in time at some location x of the fiber.The rod construction as a whole consists of many different materials with different thermal expansion coefficients, such as the sensing fibers, gel filling, metal tubes, polypropylene mantle, steel rod and nylon centralizers.However, the steel of the sucker rod and the steel of the fiber-optic mantle are the dominant material by weight and the most relevant for any thermal stresses.The sucker rod consists of 4332 SRX Nickel Chromium Molybdenum steel with a thermal expansion coefficient of 10 -13 μϵ/K (Hidnert, 1931) and a modulus of elasticity of 200 GPa (T.E.Toolbox, 2012).The second most dominant material is the polypropylene cable mantle with a modulus of elasticity of 1.5-2 GPa (T.E.Toolbox, 2012).The proportion of steel on the thermal stresses in the rod construction are 99.8%.For simplicity, we assume that thermal expansion coefficient αrod = 13 μϵ/K for the sucker rod / fiber-optic cable construction and neglect the other materials.

Direct measurement of strain via DAS
Similar to DTS, DAS also analyzes the back scatter of light coupled into a fiber from one end.Upon contraction or dilatation, the strain-rate of the fiber, i.e. the temporal derivative of relative change of length, can be derived from the temporal change of the interference pattern of coherent light elastically scattered (Rayleigh scattering) from adjacent points within a certain interval of fiber called the gauge length (Masoudi et al., 2013).The centroid of the gauge length is defined as a sensor node.The location (x) of a sensor node along the fiber is again determined by the two-way travel time of light from its source to the node and back.In our study, DAS data is acquired at 10000 Hz and down-sampled to 1000 Hz.The gauge length and spatial samping are 10 m and 1 m, respectively.
No additional filtering was applied in post-processing (no high pass and no low pass filtering).
In contrast to DTS, DAS directly yields the temporal derivative of strain.In order to convert the measured strain rate ˙(x,t) data to strain εDAS(x) at each location, we integrate in time: (2) where t1 and t2 delineate the time window and ˙(x,t) the recorded strain rate at position x.In the following we speak of "measured strain" εDAS in contrast to "predicted or expected" strain εDTS.
We compare εDTS with εDAS measurements.We then use the εDTS data to compute the contractional forces along the rod which occur due to cooling.We compare the result with a static friction curve that was estimated from the sucker rod tally and borehole inclination.

Stick-slip approach
As the thermal contraction of the cooled sucker rod inflicts a sliding movement of the rods along the borehole wall, we must consider the friction of their relative motion.This friction would yield a stick-slip motion which is observed almost everywhere when two solid objects are moving relative to one another.A detailed review of the origins of stick-slip behavior in mechanical parts as well as an experimental and theoretical analysis on stick-slip characteristics is presented in literature (e.g.Berman et al., 1995).In the simplest case, a stick-slip motion appears when the static friction force Ff between two stationary solid bodies is overcome.A schematic drawing of the forces on an interval of the sucker-rod construction at a depth with borehole inclination is presented in Figure 2. The static friction force Ff can be calculated according to where Fn is the normal force and μ the static friction coefficient.The value for μ = 0.36 was obtained from a plateto-plate experimental analysis on the stick-slip behavior between steel and glass fiber-reinforced nylon specimen (Muraki et al., 2003).The force FN is calculated according to where Fw is the gravitational weight force and α the borehole inclination.Each sucker rod element is 9.1 m long, weights 15.7 kg and is equipped with four nylon centralizers and the fiber-optic cable (20 g/m).Therefore, the weight force for each contact point of the rod construction yields Fw = 9.81 m/s 2 * 15.9 kg / 4 = 39.0 N. Regarding the lowermost part of the rod construction as an example, this means that for the last nylon centralizer (borehole inclination of 54°), a static friction force of Ff = 0.36 • 39.0 N • sin (54°) = 11.3N is calculated.With respect to contraction of an initially unstressed rod construction, for each subsequent nylon centralizer towards the surface, the friction force of the rod at the given depth is calculated by the cumulative sum of all friction forces from the nylon centralizers below.The friction force increases with decreasing well depth.Two further weights are added to the friction force profile: the bottom end of the sucker rod is a 1.4 m long steel piece with a weight of 64 kg and the carrier of the pT gauge at 2755 m MD is a 2.2 m long steel piece with a weight of 105 kg.Here, we applied a static friction coefficient for steel on steel of μ = 0.8 (Lee and Polycarpou, 2007).
The expected thermal contraction εDTS can also be translated to a force.Assuming a Young's modulus for stainless steel of E = 200 GPa (Cardarelli, 2018) and given the cross-sectional area of the rod (Arod =2.9 cm 2 ), we can calculate the applied force Fapp at each location along the rod which was thermally induced within the investigated one-hour cold-water injection period: For simplicity, we assume that the elasticity from the fiber-optic cable and the nylon centralizers are neglectable and that the steel dominates the mechanical behavior of the structure.Furthermore, we make the assumption that no mechanical stresses are exerted on the rod prior to the cold-water injection.This allows us to set a zero-force baseline before injection start for the stick-slip analysis.

Stick-slip event detection and picking
In the DAS data we monitored repeating vibrational events with ongoing cold-water injection in the deeper part of the well.These events are characterized by a sudden DAS amplitude peak at some depth and an up-and downward directed move-out.With time, the spatio-temporal distribution of these vibrational events changes.To automate the detection of depth location and moveout of an event, we employ a short-term/long-term average trigger (Allen, 1978, Vaezi andv.d. Baan (2015)).The parameters used for the STA/LTA analysis can be found in Table 2:

Results
Figure 3 shows examples of raw and unprocessed strain rate data measured with the DAS unit in the well at the start of cold-water injection (1 st subplot), one hour after start of fluid injection (2 nd subplot) and shortly after the end of the 24 hours water injection period (3 rd subplot).Each subplot depicts 10 seconds of data with the same data color scaling.A number of features can be recognized in each of the data examples.At the depth marked with the arrow "A", there is a transition from a noisy depth interval above to a rather quiet one below.The transition marks the location of the water table in the well.From the wellhead, the water free-falls down to the water table at about 170 m below surface.In the cased hole section down to the depth of the transition to the perforated liner, high velocity tube waves (around 1500 m/s) are present which are reflected at the liner shoe of the 9 5/8" casing at ca. 2810 m MD (arrow "B" in first subplot).Below "B", the cable is located inside the perforated liner.The tube waves are not further guided in this interval and the noise level is rather low.In the uppermost 100 m of the perforated liner section (2810 -2900 m MD), a strong signal is present in the 2 nd and 3 rd subplot (arrow "C").The arrow "D" marks another common characteristic feature in the DAS data which was observed over the analyzed cold-water injection period.This abrupt and localized signal is interpreted as a sudden contraction of the sucker rod.Between 2900-3100 m MD, the temperature difference between the two DTS profiles rapidly decreases (see Figure 4, 1 st and 2 nd subplot).At 02:08, the DTS profile shows slightly increased temperatures (+1 °C) with a constant offset from 3100 m to the end of the cable compared to the DTS profile at 01:18.This leads to a constant offset of a positive expected strain εDTS.The measured strain εDAS shows no offset in this depth interval.

Event description
A close-up of raw DAS data is shown for the depth interval 2500-3300 m MD around the transition from cased hole to perforated liner 52 minutes after the start of the cold-water fluid injection (see Figure 5).At this time, the DAS records a transient strain-rate anomaly.Similar events are repeatedly observed in the course of the measurement during the cold-water injection periods.Using the event shown in Figure 5 as a representative example, we describe common features of these events in the following.Its origin lies at 600 ms and 3000 m MD and is characterized by an abrupt increase of the measured strain rate.The sudden increase of strain rate amplitude propagates both up-and downwards along the well with compressional and tensional sign of amplitude, respectively, where the propagation velocity upwards is approximately 3900 m/s (green Line A in Figure 5).In contrast, the downward propagation velocity is slower and shows irregularities from 650-1260 m/s.Most striking is the decay of the velocity from 3200 m MD onwards and the eventual stop of propagation slightly above 3300 m MD.In upward direction, this event is halted somewhere in the noisy interval where the reservoir section of the borehole begins.The event is followed by elastic reverberations that decay after approximately half a second.
Further examples of such kind of events are plotted in Figure 6 A, B, C and D. They all have in common, that they originate below 2900 m MD and trigger a contraction above and an extension below.The previously discussed event is characterized by a smaller precursor 100 ms before the origin of the large event at the same depth.
Precursors and successors can also be observed in the examples in Figure 6 (in particular in Figure 6 B), yet the events shown here are distinguished by the fact that their upwards propagation extends beyond the noisy reservoir section.All exemplary events except 6A whose downward propagation arrests rather sudden, have in common that the up-and downwards propagation slow down before coming to a halt.Another striking observation in all of the events is that the initial onset propagates slower than the reverberations in the coda.
While the exact shape of the spatial propagation and length varies (length between 20-1600 m), the duration of these events is mostly in the range of 0.5 s with some fading noise/reverberation afterwards.These events typically show a tensional signal at the energy front in the downward direction while the initial energy front upwards is mostly compressional.As the vibrational signal propagates along the rod, a succession of compressional and tensional waves is created which moves with a velocity of about 3900 m/s along the rod (as shown by the green line A in Figure 5).The downward propagation of the first arrival changes its velocity from the onset of the event towards the end of the vibrational event.In the first 50 ms, it increases in velocity, then it stays constant before it gradually decreases in velocity at around 700 ms below 3200 m MD.
The four black arrows on the left y-axis in Figure 5 indicate the timeseries for which the four spectrograms shown in Figure 7 were calculated with a moving window of 250 ms.The DAS strain-rate timeseries at 3000 and 3200 m MD show the onset of the slip event at 0.5 s with dominant frequencies of the first break between 30 and 75 Hz.
The slip only lasts approximately half a second but reverberations of different duration and different frequencies can be observed in band below 30 Hz depending on the rod segment.For instance, at 3000 m MD long lasting reverberations occur at ~10Hz whereas at 3200 m MD they occur at 20 Hz.As can be seen from the spectrogram from the DAS strain-rate recordings at 2700 and 2835 m MD, the slip event does not penetrate into and beyond the feed zone whose characteristic noise at 24 Hz remains undisturbed just as the low frequency pattern of the tube waves above.

Event detection over time
We applied a STA/LTA algorithm to automate the detection of the sudden vibrational events within the first 72 minutes of cold-water fluid injection.Three attributes are obtained for each event: a) the depth location where the event starts b) the lower end and c) the upper end of the event according to the STA/LTA algorithm.Figure 8 shows  relatively small number of bigger vibrational events occur, i.e. events which extend over more than 20m.Early events (within the first 5 minutes relative to injection start) appear in the depth region between 1250-1900 m MD.
Except for two large events (4 minutes: 2260-2730 m MD and 6 minutes: 2040-2700 m MD), the spatial extent of the vibrational events is rather small.One single event was recorded at a depth of 3540-3580 m MD close to the shoe of the installation.With time, the depth of vibrational events increases to 2900 m MD.From 17 minutes onwards, the occurrence of vibrational events is mostly constrained in the depth region from 2900-3100 m MD.
The maximum spatial extent of large vibrational events increases with time.From 01:18 (+22 minutes after injection start) onwards, most of the events extend into the depth region of 2835-3080 m MD.At 02:08 (+72 minutes to injection start), the spatial extent of the events is 2500-3470 m MD.
With time, the frequency of the occurrence of the events decreases.4-5 hours after injection start, large events (such as in Figure 6 C and D) appear every 10 -15 minutes.8 hours after injection start, large events appear approximately every 25 -40 minutes.

Friction force model
The static friction force Ff after formula (3) is compared to the applied force from thermal contraction of the rod Fapp after formula (5) which was evaluated for the period from injection start to 01:18 (+22 minutes after start of injection) and to 02:08 (+72 minutes after start of injection) (Figure 10).Ff and Fapp intersect at 17.0 kN and 9.3 kN, respectively.At the depth interval from 2732-2820 m MD, the applied force peaks at 22.0 kN (shown in Figure 10).For all estimates given above, it is assumed that the sucker rod did not move relative to the casing, i.e., thermal stresses can build up but will not be released by relative motion.

Discussion
With the help of distributed fiber-optic temperature and acoustic data, we monitored a cold-water injection period in a geothermal well at the site Schäftlarnstraße, Munich.The downhole monitoring data allows for an analysis of the deformation of the 3.6 km long sucker rod/fiber-optic cable construction due to cooling.We observe numerous localized episodes of large strain-rates that nucleate along the inclined stretch of the borehole and propagate both towards greater depth and the surface.Such events induce quickly declining elastic vibrations along the entire extent of the affected interval.The emergence of these vibrational events strongly correlates with the beginning of the fluid injection.In the following, we thus argue that the vibrational events are a result of the substantial temperature changes which the sucker rods with the optic fiber are exposed to.The contraction of the sucker rods upon cooling induces stress where the sucker rod is held to the borehole wall by frictional forces.On the basis of a simple mechanic model we show that accumulated stresses may eventually exceed the friction giving rise to sudden stress release and the observed strain changes.

Assessment of measuring errors
Our monitoring data analysis is based on a debatable approach of integrating DAS data over longer time periods.
To obtain the εDAS profile over the period of 1 hour, a total number of 3.6 million strain rate profiles are integrated (sample rate: 1000 Hz).Such kind of numerical operation has a high risk of creating numerical errors due to e.g.
rounding off or value truncation.In addition, the smallest systematic error in the DAS measurement system results in a significant drift over time which would misrepresent the strain profile measured by the sensing fiber.Also, it is well known that for FIMT type of installations, the gel filling allows for a creep and differential movement of fibers with respect to its surroundings which makes strain sensing unreliable for greater deformations and longer periods (e.g.Lipus et al., 2018, Becker et al., 2020).However, a creep over many meters or even kilometers is most likely improbable.To strengthen the meaningfulness of our integrated strain profile, we analyzed the εDAS for a deeper section of the well, where little temperature change (ca. 1 °C) was measured by the DTS.In 3500 m depth, we do not observe any strain accumulation after temporal integration of strain rate data over a period of 60 minutes.This indicates that the measured strain rate has no significant drift during the time of interest.For measurements with higher amplitudes such as within the depth interval 2800-2900 m, non-linear effects influencing the temporal integration of the data cannot be excluded.

DAS data integration
We integrated DAS data in time over 72 minutes to assess the absolute contraction of the rod construction prior to the cold-water injection start (see Figure 4).For the well interval from water table to the transition to the perforated liner, the results show a good match to the contraction that was theoretically assumed from the cooling of the well.
However, from 2800-2900 m MD, we obtain much higher deformation from the DAS data than what we expected.
We cannot give an unambiguous explanation for that but see two likely reasons for that observation.Firstly, the DAS integration process might result in a drift when integrating high amplitude DAS data.Especially from 2800-2900 m MD, constantly high energy is recorded by the system.The second explanation could be that the integrated DAS data measured a true deformation of the construction.In the depth region around 2800 m, the annular space of the borehole is rather irregular (transition to 7 " liner interval, localized increase in the borehole inclination (see Figure 4, 3 rd panel)).The repeating sudden sucker rod events might lead to an uneven distribution of the thermal stresses along the rod.Interestingly, the most prominent feed zone of the well coincides with the one single DAS interval which shows an extensional signal.
The sudden slip events presented in this study show some similarity to the "slip events" which were previously observed in FIMT-type fiber-optic installation in a geothermal well (Miller et al., 2018).In the reported DAS campaign, a fiber-optic cable was installed in a geothermal well and it is argued that repeated thermal cycles led to a loss of frictional coupling between fiber-optic cable and the borehole wall.Miller et al., 2018 reported that a sudden loss triggered a movement of the cable with a first arrival speed of 4600 m/s (we measured a first arrival speed of 4000 m/s).The integrated strain of the reported event shows a balance towards absolute contraction which we also observe in our events.Another similarity is given by the frequency content of these events.They recorded a dominant frequency of 45 Hz with some harmonics in both directions which we also observed in our data (see Figure 7 at 3000 m MD).

Stick-slip rod behaviour
We calculated the static friction force Ff along the rod construction by a cumulative sum of the friction of each nylon centralizer with the borehole inner wall.Independently of that, we computed the applied force Fapp on the rod construction by thermal contraction using the DTS monitoring data.By comparing both curves, we can distinguish depth regions where the rod remains immobile (Ff > Fapp) and depth regions where the applied forces overcome the static friction force (Ff < Fapp).The temperature difference in the course of the investigated time period is particularly high over the 9 5/8" liner interval (depth region from 2485-2890 mMD) which in consequence also means that Fapp is high.According to our model calculation, the contraction forces surpass the frictional forces at 2800 m MD around 01:18 (22 minutes after injection start).This result implies that after this time, the construction can contract in this depth interval.In other words, the thermal stresses on the rod construction in this depth region are high enough that the rod starts to move and to contract.Hence, the literature values assumed for the static friction between sucker rod and steel liner are assumed to approximate the real values.
With ongoing cold-water injection and further cooling of the well, the applied forces Fapp increase.This leads to a continuous growing of the depth interval where Fapp surpasses the static friction Ff of the rod.The STA/LTA detections match the predictions of the friction fore model.After a rather quiet initial phase of low energetic events (before 17 minutes in Figure 9) which could be caused by the relaxation of previously accumulated stress anomalies along the sucker rods, repeated vibrational events start to concentrate in the region 2800-3100 m MD.
As the region with Ff < Fapp increases, the length of the vibrational events increases.From our friction force model, we would expect vibrational events (more specifically: the contraction part of the movement) at 02:08 in the depth region 2150-2910 m MD.However, the observed events extend from 2500-3500 m MD.Regarding the upper limit, we can see in Figure 10 that there is a significant change in slope for Fapp at 02:08 at 2500 m MD.The friction force model is based on numerous assumptions (i.e.static friction coefficient nylon-steel, Young's modulus for stainless steel, neglecting fiber-optic cable, stress-free initial conditions) which might not accurately depict the downhole conditions.This could mean, that either the calculated applied force Fapp is too high and/or the static friction force Ff is too low.
With respect to the lower limit of the vibrational events, we predict the contraction part (Fapp(02:08) Figure 10) of the vibrational events down to a depth of 2912 m MD from our friction force model.However, we record vibrational events down to a depth of 3480 m MD.This discrepancy can partly be explained by the fact that the model prediction only shows the contraction part of the vibrational event.As seen in the cumulative strain εDAS (see Figure 6 event A and B), the lowest part of a vibrational event yields extension.The most likely reason is that the contraction above results in a pulling of the rod from a lower lying region to compensate for the missing rod length.Therefore, the events can be traced down to a greater depth than predicted.
The constant temperature offset by +1 °C in the DTS profiles from 02:08 (relative to 01:18) in the depth interval from 3100 m MD to the end of the cable is unlikely to be caused by any fluid movement.While DTS temperature measurements did show a variation, no additional offset was recorded from the measured strain εDAS.This could mean that the rod builds up thermal extensional stresses without actual movement taking place (εDTS > 0 εDAS = 0).
However, we speculate that the temperature anomaly is related to the processing of the DTS data.DTS temperature was measured in a double-ended configuration.A temperature profile is created by overlaying the DTS signal from both directions which are measured consecutively for both fiber branches.Close to the folding location (at the bottom of the well), an asymmetry in the temperature reading was observed between both fiber branches, which does not seem to be caused by any fluid motion.Averaging this difference between both branches led to a temperature offset.This offset was only visible if strong temperature changes were observed.

Conclusion
The field test at the geothermal site Schäftlarnstraße demonstrates that simultaneous recording of DTS and DAS data can be used for a detailed analysis of the deformation of a sucker rod type of fiber-optic cable installation in a 3.6 km deep well.By comparing the theoretical contraction of the rod structure from DTS with an estimated contraction from DAS, we can distinguish depth intervals with higher and lower thermal stresses in the material.
We introduce a friction force model which accurately predicts the onset and extent of sucker rod events releasing accumulated thermal stress.This is an important finding for DAS monitoring in geothermal settings because it shows that localized high-energetic vibrational events must not necessarily be related to microseismic events occurring in the rock formation but can originate in the subsurface construction and the way how the fiber-optic monitoring equipment is installed in the well.Moreover, the friction force model is useful to predict the data quality for DAS measurement campaigns for deep sucker-rod types of fiber-optic installations.Especially for the recordings of weak acoustic signals that are e.g.induced by fluid movement in the annulus, it is essential to know the potential sources of errors and artifacts in the data.During operations which introduce a temporal temperature gradient, thermo-mechanical response of freely hanging steel parts in the borehole may introduce stick-slip events that must be distinguished from any other relevant seismogenic source.Potentially, the vibrational energy from the sucker rod events can also be used to study the formation velocity in the near-field around the borehole.
Furthermore, the large-scale contraction along certain sucker rod and fiber intervals must be considered with respect to the location of the distributed sensor nodes.Our description also serves as a starting point for a more detailed dynamic description of the observed processes.This can be of use to predict onset and depth interval of such sucker rod events and to contain their destructive potential in case of too quick cooling of the construction.

Figure 2 :
Figure 2: Static friction force Ff and normal force FN applying on a sucker-rod contact point (nylon centralizer) as a function of the weight force Fw and the borehole inclination 90°-α

Figure 3 :
Figure 3: DAS raw data examples over the scope of the cold-water injection phase for (1) the onset of fluid injection (2) ongoing injection and (3) termination of fluid injection.Blue colours show relative compression and red colours relative expansion.The color ranges are the same for all subplots

Figure 4 :
Figure 4: Downhole monitoring data during the cold-water injection test. 1 st subplot: DTS temperature profiles.2 nd subplot: Comparison of strain profiles εDTS and εDAS.3 rd subplot: borehole inclination.4 th subplot: wellbore schematicIn general, a clear match is visible between εDTS and εDAS for the entire well which means that the strain the steel rod experiences (εDAS) follows the predicted thermal contraction (εDTS).However, there are depth intervals where the experienced strain (εDAS) exceeds and others where it falls short on the predicted strain (εDTS).Until 2825-2835 m MD where the most prominent injection interval is located, ΔT increases with increasing depth.At the injection interval ΔT rapidly increases.Below this zone, no thermal contraction is expected.

Figure 5 :
Figure 5: Sucker rod contraction event displayed by strain rate DAS data (left).The black arrows on the left y-axis mark the depth location of timeseries used for the spectrograms in Figure 7. Line "A" marks the moveout of the signal at a speed of 4000 m/s.The schematic drawing shows the inclination of the borehole with the fiber-optic cable (red) lying inside of the casing (right).The inflow profile from a wireline flowmeter measurement is shown by the blue graph one example of the automated detection with the STA/LTA trigger.The upper subplot shows an example trace of raw DAS data at a depth of 3120 m MD (marked by the black arrow in the lower subplot) and the corresponding STA/LTA characteristic function.Beginning and end of the detection are marked by the green and orange crosses, respectively.The lower subplot shows spatio-temporal DAS data and the detection of two vibrational events.

Figure 8 :
Figure 8: STA/LTA trigger algorithm applied as an automated detection method for vibrational events.Trigger start and end is marked with green and orange crosses All vibrational events which occurred within the first 72 minutes of cold-water fluid injection are plotted in Figure 9. Gray circles mark the spatio-temporal origin of vibrational events.The corresponding vertical black line indicates the spatial extent of the respective event.In this representation, events with a spatial extent of less than 20 m are neglected.Such small events occur between 4-10 times per minute in the depth region from 1250-2750 m MD over the entire investigated 72 minutes after fluid injection start.Within the first 15 minutes, only a

Figure 9 :
Figure 9: Gray circles and black vertical lines indicate the spatio-temporal origin and spatial extent of vibrational events in the well, respectively.The shown period comprises the first 72 minutes of cold-water fluid injection.
The gravitational weight force Fw per nylon centralizer is constant for every contact point of the rod.The force needed to overcome the cumulative static friction Ff is a function of the borehole inclination.Ff increases from the bottom of the rod installation at 3691 m MD towards 1000 m MD.The bottom end of the sucker rod and the carrier of the pT gauge at 2755 m MD create an additional static friction force of 0.4 kN and 0.5 kN, respectively.Above 1000 m MD, the well is nearly vertical and only little static friction is expected.The static friction Ff at 1000 m MD yields 26.1 kN.Fapp at 01:18 is lower than Ff for the entire installation length.Only in the depth interval 2731-2820 m MD, Fapp approaches a force of 10.5 kN which is close to Ff.This indicates that forces are sufficient to initiate relative motion between sucker rod and casing at that depth.With ongoing cold-water fluid injection, the applied forces Fapp increase with further decreasing temperatures.At 02:08, Fapp surpasses the frictional forces in the depth range from 2150-2912 m MD.

Figure 10 :
Figure 10: Comparison of static friction Ff with applied forces Fapp from thermal contraction of the rod within the first 72 minutes of cold-water fluid injection.The pale colors in Fapp originate from measured DTS data and the solid lines are constructed by a moving average over 15 m

Table 1 :
Well design at geothermal site Schäftlarnstraße, Munich (see also Figure3)

Table 2 :
Parameters used for the STA/LTA detection method