Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Rippled scales in a geothermal facility in the Bavarian Molasse Basin: a key to understand the calcite scaling process

Rippled scales in a geothermal facility in the Bavarian Molasse Basin: a key to understand the... bernhard.koehl@tum.de Institute of Hydrochemistry, Mineral precipitates, scales, are wide-spread in geothermal facilities and hamper Technische Universität technical and economic efficiency. Scales reported from earlier studies in the Bavar - München, Marchioninistrasse 17, 81377 Munich, Germany ian Molasse Basin consisted mostly of calcite, appeared rather homogeneous, grew Full list of author information flat and perpendicular to the surface, and showed little or no internal macrostructure is available at the end of the apart from small gradients in crystal size and chemical composition. In this study, we article investigated extraordinary calcite scales collected from the inlet chamber of a heat exchanger. These scales showed current ripples and a sequence of coarse and dense layers. By an extensive investigation of the ripples, the hydrochemistry, and by the application of existing empirical equations, the formation of the ripple scales could be reconstructed: Due to degassing of CO in the geothermal pump calcite scales precipitated on the pipes of the ground-level facilities. During maintenance works, these scales were partially dissolved, and remobilized as individual mobile grains. These particles subsequently settled again, and solidified as ripple scales in the horizontal pipes at surface level. For this ripple forming process, flow had to be small, because otherwise the ripples would have been altered or washed away. After deposition of the ripple structure, a dense new calcite layer precipitated on top of the ripples, and stabi- lized it due to transitory supersaturation. This sequence of ripples and a dense calcite cover layer is repeated two times. Finally, during start-up of regular plant operation, the whole composite scale was remobilized and got washed into the inlet chamber of a heat exchanger. The investigation presented here indicates that this type of calcite scale only forms in geothermal facilities in the Bavarian Molasse Basin if certain flow and hydrochemical conditions are present during the maintenance procedure. Keywords: Scales, Scaling, Precipitations, Calcite, CaCO , Ripples, Malm, Bavarian Molasse Basin, Carbonaceous aquifer, PhreeqC Introduction Mineral precipitation and build-up in pipes, also known as scales or scalings, are a com- mon problem for geothermal plants in the Bavarian Molasse Basin. Growth of scale in pipes reduces the technical and economic efficiency of geothermal plants and requires frequent, costly maintenance or even process modification of installed geothermal plants. Mineral build-up can occur on any surface exposed to geothermal fluid, including © The Author(s) 2020. 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/. Köhl et al. Geotherm Energy (2020) 8:23 Page 2 of 27 the submersible pump, its motor, the production pipes, and ground-level components, such as heat exchangers (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). In previous studies in the Bavarian Molasse Basin scales had a more or less even sur- face and the scales were interpreted to form during normal operation of the facility (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). The scale samples were analyzed by a variety of chemical and mineralogical methods (SEM/EDX, XRD, acid digestion, polarizing microscopy) (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). The mor - phology of the scales, however, received little attention, mainly because the observed scales grew evenly on the pipes. The scales presented in this study show a distinct ripple pattern, which was not reported anywhere else (Fig.  1). The macroscopic structure was used to gain insight into the formation process and to develop countermeasures against this type of min- eral precipitations. The geothermal facility investigated in this study is located in the Bavarian Molasse Basin, and exploits Upper Jurassic sedimentary rocks, also known as Malm aquifer. This formation, especially in its dolomitic karstified facies and its massive reef facies hosts high permeabilities and offers favorable conditions for geothermal use (Bauer et al. 2014). The formation is situated at a depth of 2000 m north of Munich and up to 5000 m in the South. Flow rates reach up to 150 L/s and temperatures up to 150 C (Agemar et al. 2014; Bauer et al. 2014). The formation of scalings is controlled by the hydrochemical conditions. Therefore, we give a brief overview of the conditions in the target reservoirs. Within the Bavarian Molasse Basin, six hydrogeochemical provinces can be distinguished (Birner et al. 2011). Following this classification, the Munich area belongs to the central part of the Molasse Basin and shows a Na–Ca–HCO –Cl-type or Na–HCO –Cl-type water chemistry after the classifica - 3 3 tion of Furtak and Langguth (1965). South of Munich, the influx of waters with higher min - eral content from the overburden increases (Mayrhofer 2013). In the central Molasse Basin, Fig. 1 Ripple scale sample collected from the inlet of the heat exchanger. The sample has a calcite composition. Some oil contained in the thermal water is incorporated in the scale causing the black color. The shape of the ripples with a luv and lee side points to current flow ripples Köhl  et al. Geotherm Energy (2020) 8:23 Page 3 of 27 the Malm aquifer has high gas contents up to 2 L gas per L water. The gas is mainly com - posed of CH , CO , and N . H S concentrations in the gas phase can reach up to 5 vol %. 4 2 2 2 Under reservoir conditions, due to the long residence time the water is in chemical and thermal equilibrium with the surrounding rock (Wanner et  al. 2017). Here, the key equi- librium is the carbonate equilibrium, and the driving factors are temperature and pres- sure. During the production of the thermal water, the equilibrium is easily disrupted. Local increase of the temperature and degassing of CO cause oversaturation with respect to car- bonates (Wanner et al. 2017). On the other side, the cooled water injected into the reservoir shows a benign behavior with respect to scales due to the reduced temperature. In contrast to any other scale samples reported from geothermal sites in the Bavarian Molasse (Wanner et  al. 2017, Herbrich 2017, Boch et  al. 2017), the scales in the ground- level facilities of the investigated geothermal power plant showed a distinct morphology (Fig.  1) with ripple structures usually known from river beds or sea shores. In general, ripples are formed when either wave action (wave or orbital ripples) or current flow (cur - rent ripples) over morphological perturbations of the bedform lead to abnormalities in local flow conditions, such as vortices (Ha and Chough 2003). This results in local areas of increased and decreased shear velocity where sediments are eroded from or re-deposited onto the bedform, respectively. An example of flow regimes along with erosion and deposi - tion mechanisms near ripple structures is shown in Fig. 2. Geometric indices can be used to determine whether ripples were formed by waves or currents (Tanner 1967). The ripple symmetry index (RSI) is the ratio of stoss length to lee length. A RSI < 2.5 is indicative of wave ripples, whereas RSI > 3.0 is pointing to current ripples (Allaby 2008). The boundary conditions at which the scales under investigation have been formed are ruling out oscillation ripples (wave ripples). The formation of ripples is governed by particle diameter, velocity, fluid viscosity, and relative density of the sediment to the fluid. Although ripples are not static but in con - stant motion, the shape of the ripples, that is the periodicity of the ripples (wavelength), the height of the ripple crests, and the RSI, to name a few reaches an equilibrium state. The ripple-forming parameters have different effects: While experiments have shown that the equilibrium morphology of current ripples depends on the particle grain size and their immersed density (Kraemer and Winter 2016), the time to reach a stable ripple pattern depends on fluid density and fluid velocity (Baas et al. 2019). Many studies have attempted to derive empirical relationships for transient and steady- state ripple morphology (ripple height, η (m), ripple length,  (m)) based on median max max particle diameter ( d ) ( Y alin 1985; Flemming 1988; Baas 1994; Coleman et al. 2003; Soulsby et al. 2012) (Eqs. 1–5). Median particle sizes used in these studies range from 0.045 mm to 0.87 mm, the particles were either glass spheres or sand, and experiments were often imple- mented in an open-channel rectangular flume. =1000 · d max (Yalin) 50 (1) η =0.0677 · 1000 · d max (Flemming) 50 (2) 0.75 =175 · (d ) (3) max (Coleman) 50 Köhl et al. Geotherm Energy (2020) 8:23 Page 4 of 27 Fig. 2 Schematic drawing of ripples with indication of ripple length, height and key formation processes modified after Ha and Chough (2003) The range of flow velocities that result in ripple formation for a given particle size has also been studied (Soulsby et al. 2012). The minimum flow velocity is given by the shear forces which have to exceed the force required to erode the sediment bed. At flow veloci - ties lower than the minimum the sediment grains simply stay put. There is also a maxi - mum flow velocity for ripples, again depending on the grain size. Above this velocity the particles are washed out and ripples are destroyed, and the bedform exhibits sheet flow behavior. If ripples have been formed at suitable flow velocities, and the velocity later drops below the minimum threshold, the ripple structure in the bedform is conserved (Soulsby et al. 2012). Flow velocities resulting in ripples can also be read-off from several bedform phase diagrams which confine the boundary conditions for ripples and other bedforms (Southard and Boguchwal 1990; Berg and Gelder 2009; Baas et al. 2016). For instance, the diagram of Southard and Boguchwal (1990) (Fig.  3) shows the bed sedi- ment size against flow velocity for flow depths between 0.25 and 0.4 m. At these condi - tions, as they occur in natural systems, and with pure non-cohesive sediment, ripples can form at velocities of roughly 0.1 - 0.9 m/s, depending on sediment particle size. For mixtures of cohesive and non-cohesive sediments, the equilibrium morphol- ogy of current ripples appears to follow the same empirical relationships that were found for non-cohesive sediments, but the time required to reach the equilibrium morphology is longer (Baas et  al. 2019; Schieber and Southard 2009; Schieber et  al. Köhl  et al. Geotherm Energy (2020) 8:23 Page 5 of 27 ◦ ◦ Fig. 3 Bedform diagram showing the 10 C-equivalent mean bed sediment size against 10 C-equivalent mean flow velocity for 10 C-equivalent flow depths between 0.25 and 0.4 m after Southard and Boguchwal (1990). Indication of the equivalent temperature is needed due to the effect of temperature on fluid viscosity and density. The diagram shows that current ripples are likely to occur at mean flow velocities between roughly 0.1 and 0.9 m/s and at a sediment size of roughly 0.1–0.8 mm 2013). When cohesive sediments form ripples, the particles flocculate, which enables them to migrate as bedload creating ripples. Floccules composed of silt and clay are typically in the range of 100 - 200µm (Yawar and Schieber 2017), sometimes several 100µm (Schieber et  al. 2007). To our knowledge the shortest ripples reported in lit- erature can be found in Baas et al. (2019) who observed wavelengths down to 77µm . These ripples consisted of a mixture of cohesive sand and 13  % mud fraction. Flow velocities leading to the formation of mixed cohesive and non-cohesive ripples are in the range of 0.1 - 0.7 m/s, but typically lower than velocities needed for pure non- cohesive sediments (Baas et al. 2019; Schieber et al. 2007; Yawar and Schieber 2017). If the cohesive portion exceeds about 30  %, no ripples form and a flat bed remains (Baas et al. 2016). Baas et al. (2016) published a review on the available equations for the prediction of rip- ples. Baas found that Soulsby et  al. (2012) created the most comprehensive predictor for current generated and combined-flow bedforms capable of predicting ripple geometry and ripple-forming velocity. Therefore, we are applying these equations (Soulsby et al. 2012) for the calculations in this paper. Equations 4 and 5 predict the maximum ripple height, η max (m), and maximum ripple length,  (m), for a given mean diameter d (m) and dimen- max 50 sionless relative particle density s, which is described by the dimensionless parameter D defined in Eq. 6 : −0.554 η =202d D | , (4) max 50 1.2<D <16 ∗ ∗ −1.5 =d (500 + 1881D )| , (5) max 50 1.2<D <16 ∗ ∗ Köhl et al. Geotherm Energy (2020) 8:23 Page 6 of 27 1/3 g(s − 1) D =d , (6) ∗ 50 where g is the gravitational constant ( m/s ), s = ρ /ρ is the dimensionless ratio of p f particle density (ρ ) to fluid density ( ρ ), and ν (m /s) is the kinematic viscosity of the p f fluid. Flow regimes favoring ripple formation can be determined by comparing the non- dimensional Shields number θ (Eq. 7) to the critical Shields number θ (Eq. 10) and the c cr washout Shields number θ (Eq. 11) (Soulsby et al. 2012): wo θ <θ no bedform movement, c cr θ < θ <θ ripple formation, cr c wo θ <θ ripples are washed out, resulting in either washed-out ripple or sheet bed wo morphology: C U θ = , (7) g(s − 1)d 0.40 C = , (8) ln(h/z ) − 1 z =d /12, (9) 0 50 0.3 θ = + 0.055[1 − exp(−0.02D )], cr ∗ (10) 1 + 1.2D 0.916, if D ≤ 1.58 ′ ∗ θ = (11) wo −1.3 1.66D , if D > 1.58, where U (m/s) is the average cross-sectional flow velocity in the pipe, C is the non- dimensional drag coefficient due to the particles (Eq.  8), h (m) is the water depth (for a full pipe this is the pipe diameter), and z (m) is the grain-related bed-roughness length (Eq. 9). Combining Eqs. 7, 8, and 9, and solving for U , one arrives at an expression for aver- age fluid velocity as a function of the Shields number (Eq.  12). Plugging values for θ cr and θ into Eq.  12 gives the lower and upper average fluid velocities, respectively, at wo which ripples would form for a given d and s. 1/2 U = 2.5(ln(12h/d ) − 1) θ g − 1 d (12) 50 c 50 For the characterization of the flow in the geothermal pipes investigated in this paper, Reynolds numbers can be used. Transition from laminar flow to turbulent flow occurs at Reynolds numbers of 2300 (below laminar) to 2900 (above fully turbulent) (Rotta 1956; Schade 1989; Avila et  al. 2011; Freimann 2014). The Reynolds numbers in the investi - gated facility (see Table  1) are on the order of Re = 10 and well in the fully turbulent flow regime. Köhl  et al. Geotherm Energy (2020) 8:23 Page 7 of 27 Table 1 Pipe dimensions and flow conditions during operation, and *) during maintenance ◦ ◦ Tubing ID, mm ˙ v, m/s Re 139 C Re 20 C V , L/s DN 245 250 255 110 2.15 - 2.51 · 10 DN 300 300 316 110 1.40 - 2.02 · 10 DN 350 350 348 110 1.16 - 1.84 · 10 DN 350 *) 350 348 4 - 8 0.04 - 0.08 - 1.5 −2.9 · 10 In this study, we investigate whether the framework of ripple formation provides access to the yet unknown flow conditions which led to scale formation in the surface- level pipes of the geothermal facility. A camera rover inspection of the pipes conducted shortly after the sampling date showed only scales at the bottom part of the horizontal or slightly slanted pipes indicating deposition. We assume that the shape of the scales will allow to differentiate between regular operation (110 L/s, respectively, 1.2 –2.2 m/s, depending on the pipe diameter) and maintenance operation (0–8  L/s, respectively, 0–0.08  m/s). We also have to assume at least one step to cement the particles which make up the scales and form a solid matrix. When this cementation took place is of great interest to the operators to prevent the formation and further displacement of this type of calcite scales in geothermal facilities in the Bavarian Molasse Basin. Materials and methods Site description and operation The site under investigation produces geothermal water from the carbonate reservoir in the sediments of the Upper Jurassic in the Bavarian Molasse Basin south of Munich. The temperature at the well head is 139 C and the typical flow rate is 110  L/s. During production the pressure decreases from about 300 bar to 12 bar ground-level pressure. Table 1 gives the dimensions, flow velocities, and Reynolds numbers for the major pipes in the facility. The pipes are entirely filled with water so that water depth of the horizon - tal pipes (DN 300 and DN 350) equals pipe diameter. Reynolds numbers were calculated ◦ ◦ for a pressure of 12 bar and temperatures of 139 C and 20 C. The calculated Reynolds numbers are on the order of 1 · 10 for full production in front of the heat exchanger. During maintenance, the Reynolds numbers are still in the 1 · 10 range. The numbers indicate turbulent flow in all pipes, and at all flow conditions. In order to prevent detached scales from the production pipes, rock fragments from the reservoir, and left-over cuttings from the drilling of the well to enter the produc- tion plant, particle filters with a cutoff of 50µm are installed close to the well head. Two automatic filters are operating in parallel. Any larger particles detected behind the filters have likely precipitated in the surface-level installations or grew from smaller particles. Samples were collected after a series of interruptions of the geothermal energy pro- duction (see Table  7) in Oct. 2015. One of the automatic particle filters was in service during the majority of this period. There were two maintenance cycles to mechanically remove scales from the inside of the heat exchanger (Aug 20, Sep 3). Here, the electric submersible pump (ESP) was kept running, bypassing the heat exchanger. All other inci- dents came with a shutdown of the ESP, hence the flow velocity in the facility was zero. Köhl et al. Geotherm Energy (2020) 8:23 Page 8 of 27 Unexpectedly, right after the yearly maintenance, huge amounts of scales blocked the heat exchanger, and caused another interruption (Oct 1). Cleaning of ground‑level facilities The procedure for the chemical cleaning of the ground-level pipes comprised the follow - ing steps: First, the ESP was shut-off, and the ground-level facilities were disconnected from the production and injection pipes. Second, the geothermal water in the surface- level pipes was displaced by tap water. Then, 70 % methanesulfonic acid and a surfactant solution with a pH around  7 were gradually added to the tap water in the pipes until the concentration of the acid in the fluid reached the desired 5  % at a dilution ratio of about 1:17, and a pH of 0 - 1. The surfactant was needed to emulsify oil, which is con - tained in the geothermal water, and to improve interaction of scales and acid. The clean - ing solution was continuously pumped through the facility at a flow rate of 4 - 8 L/s for approximately 8 h. The solution was then allowed to sit quietly overnight. The following day, the same solution was again continuously pumped through the facilities for 8 h of flow and allowed to sit overnight. Starting the next day the acid was drained from the facilities into a neutralization tank. The ground-level pipes were filled with tap water before normal production was restarted. The pressure during the entire cleaning oper - ation was held at 1.2 bar to 3 bar by a nitrogen gas pressure system. In total, 1,800  L of concentrated acid along with 2,000 L surfactant solution were consumed during the acidification. Methods Collection and analyses of water and gas samples Water and gas samples were taken at regular intervals before and after the cleaning procedure, and showed little variance (i.e. ± 4% TDS Herbrich 2017). The samples pre - sented here were taken after the incident and after full operation resumed. Samples were taken at a sampling port between well head and the automatic filter. An on-site cooling device was used to bring the temperature to 54 C to be able to handle the samples safely. Temperature, pH-value, electric conductivity, and redox potential were measured on-site with portable electrodes and sensors (InLab Expert Go-ISM, InLab 738-ISM, InLab Redox ORP, Mettler Toledo, Germany). Alkalinity and acidity were analyzed by mass analysis on-site with 0.1 M HCl (Carl Roth), 0.1 M NaOH (Carl Roth), and using colorimetric indicators phenolphthalein (Merck) and cooper (Sigma − − − 2− Aldrich). Anion ( F , Cl , Br , SO ) concentrations were quantified by IC (Dionex, IC25). Total sulfide in water was precipitated on-site with 10 mL zinc acetate (2 %). The concentration of the precipitated zinc sulfide was measured with the methylene blue colorimetric method (Deutsches Institut für Normung, 1995) using a spectrophotom- + + eter (Double Beam UV-190, Shimadzu). Na and K were quantified using a flame pho - tometer (Eppendorf ELEX 6361), calcium was determined by titration with EDTA. The other concentrations of the other cations were determined by AAS (PerkinElmer, 3300). The concentration of dissolved gases in the thermal water loading was quantified using an on-site degassing unit as described by Baumann (2016). The unit was composed of a 9−L plastic tank with two connectors in the lower third of the tank serving as inlet and outlet of the thermal water. A third connector at the top of the tank was used to transfer Köhl  et al. Geotherm Energy (2020) 8:23 Page 9 of 27 the accumulated gas into a gas cylinder (1  L gas). At the beginning the sampling tank was filled completely with thermal water with the outlets for water and gas open. After a constant flow of water was established and measured, the gas outlet was closed. The gas which then accumulated inside the container was measured as a function of time using an imprinted scale. The gas composition was later analyzed by GC–MS (Shimadzu GCMS-QP2010 Ultra, Kyoto, Japan) with helium as carrier gas. Sampling, preparation, and analysis of scales The scale samples (20 samples in total) were taken from the inlet chamber of the heat exchanger on October 1, 2015. The facility was in full operation for 6 days before shut down and opening of the heat exchanger. A total of 46, 700 m of thermal water had passed through the heat exchanger after the last chemical maintenance. In order to assess the mineralogical composition and the structure of the scales by optical transmission microscopy, the samples were rinsed with n-hexane (Chemsolute, min. 99 %). Then, thin sections of an approximate thickness of 30μm were cut parallel to the flow direction and perpendicular to the ripple crests. The sections were embedded in a polyurethane resin (Araldite ). For better visualization of the pores, a blue dye (Stru- ers, Bluedye Farbstoff für Imprägniermittel DE-L831203) was added to the resin. The samples were then visually examined with an optical transmission microscope. To prepare the samples for scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM/EDX), the scale fragments were also rinsed with n-hexane (Chem- solute, min. 99 %) about five times to remove any oil. The samples were then rinsed with water once or twice and dried in a desiccator for a few days to remove residual moisture. After sputtering a gold layer onto the samples the samples were visualized with a SEM (Tescan Vega, 3 LM, Software: Vega TC, Kohoutovice, Czech Republic) and measured with EDX (EDX Oxford Instruments, x-act, Software: Aztex, Austin TX 78744, USA). The SEM/EDX images were processed with ImageJ (Schindelin et al. 2012; Rueden et al. 2017) in order to measure the pore- and crystal-size distribution. The porosity and pore size distribution of scale samples were measured using micro computer tomogra- phy (micro-CT) (v|tome|x s 240, General Electrics) by C. Schulbert for a typical sam- ple cube with 0.5 × 0.5 × 0.5 cm at the University of Erlangen. The scan was run at 80 kV, 100µA and a distance between object and focal point of 34 mm. Data reconstruc- tion and processing was done with the GE-software datos, Version 2.3. and the software VGSTUDIO MAX 3.0 (VolumeGraphics) (GmbH 2019). Hydrogeochemical modeling In this study we also wanted to test our hypothesis that the observed ripple scale did not form by simple precipitation during normal operation, but by deposition of preexisting particles, which we assume due to the ripple morphology. For this test, a hydrogeochem- ical model for normal operation was set up with PhreeqC, v.3.4 (Parkhurst and Appelo 2013) following the procedure of Baumann (2016). The model accounts for the differ - ence in ambient conditions between the sampling location and the geothermal reservoir, and uses the database phreeqc.dat. The first step to simulate the water and gas composition throughout the facility was a back-calculation of the samples taken at surface level to reservoir conditions. Next, the Köhl et al. Geotherm Energy (2020) 8:23 Page 10 of 27 Table 2 Analysis of the geothermal water from 24 June 2016 Field parameters pH (water at surface) 6.31 Electric conductivity, µS/cm T = 25 C 639 ref Redox potential, mV -247 One liter of water contains: Mass Equivalents mg mmol % Cations 129 5.62 76.5 Sodium (Na ) Potassium (K ) 20.5 0.524 7.13 3.10 0.172 2.34 Ammonium (NH ) 2+ 16.7 0.833 11.3 Calcium ( Ca ) 2+ 2.25 0.185 2.52 Magnesium (M g ) 2+ 0.54 0.012 0.17 Strontium (S r ) 2+ 0.17 0.003 0.03 Barium (Ba ) 2+ 0.00 0.000 0.00 Manganese (M n ) 3+,2+ 0.01 0.001 0.01 Iron(II)+(III) as Iron(II) (Fe ) Sum of cations 172 7.35 100.0 Anions Fluoride (F ) 2.81 0.148 2.03 Chloride ( Cl ) 88.2 2.49 34.1 Bromide (Br ) 0.45 0.006 0.08 2− 12.4 0.258 3.54 Sulfate (S O ) Hydrogen sulfide ( HS ) 1.91 0.058 0.79 Hydrogencarbonate (HCO ) 265 4.33 59.4 Sum of anions 371 7.29 100.0 Undissociated matter Silicic acid (H SiO ) 105 1.34 2 3 Boric acid (H BO ) 9.20 0.149 3 3 Total dissolved solids ( TDS) 656 pressure and temperature were adjusted according to the defined simulation points. For each simulation point the model calculated the bubble pressure which is the pressure at which gas bubbles begin to form in the geothermal water (gas phase > 0.1 mL/L ) and the theoretical mass of precipitates. Here, the mass was given as the amount of calcite, that can precipitate to reach thermodynamic equilibrium (SI  =  0). This simplification was justified by the mineralogical assessment of the scales which consist primarily of calcite. Quantification of ripple features The ripple geometries were measured by sliding caliper measurements. The measured parameters were ripple length, ripple height, and sample thickness. Results and discussion Hydrochemical conditions A representative analysis of the geothermal water is given in Table 2. Following the clas- sification of Furtak and Langguth (1965), the water is a Na– HCO –Cl-type with a TDS 3 Köhl  et al. Geotherm Energy (2020) 8:23 Page 11 of 27 Table 3 Analysis of  the  gas contained in  the  geothermal facility normalized to  standard ◦ ◦ conditions (1013 hPa, 0 C ). Degassing took place at 54 C and p = 1030 hPa Gas CH N CO C H H S 4 2 2 2 6 2 Gas, Vol % 46.3 35.9 16.0 1.0 0.3 Dissolved gas @ reservoir conditions (298 bar, 2.1 1.4 3.1 0.04 0.17 140.2 C ), mmol/L Gas load (sum of all gases), NmL gas/ L water 86 Table 4 Pressure, temperature, and saturation index throughout the geothermal facility Pressure, bar Temperature , C SI Calcite Reservoir during production 298.0 140.2 0.0 Pump intake 30.7 140.2 0.24 Pump outlet 84.4 140.2 0.19 In front of heat exchanger 12.0 139.0 0.24 Behind heat exchanger 10.5 60.0 1.00 of 656 mg/L. With a (Na+K-Cl)/HCO ratio of 0.84 the water is strongly influenced by cation exchange. The Ca/Mg ratio of 4.5 of the water indicates a limestone reservoir matrix (Langmuir 1971). The water contains 86  mL/L gas (at standard conditions, Table  3), which is com- posed of methane (46.3  %), nitrogen (35.9  %), carbon dioxide (16.0  %), ethane (1.0  %) and hydrogen sulfide (0.3  %). At reservoir conditions (298  bar, 140.2 C ) this equals 2.1  mmol/L methane, 1.4  mmol/L nitrogen, 3.1  mmol/L carbon dioxide, 0.04  mmol/L ethane and 0.17  mmol/L hydrogen sulfide. Carbon dioxide is in carbonate equilibrium and is stripped from the solution if other gases form gas bubbles. The calculated saturation index for calcite along the geothermal cycle is given in Table  4. With an age of the thermal water exceeding 10,000 years we can safely assume equilibrium conditions in the reservoir. At the pump intake, the SI increases to 0.24 due to decompression, slightly decreasing to 0.19 at the outlet of the pump due to pressure build-up in the pump, and rising back again to 0.24 at the well head, again due to decom- pression (Table  4). The pressure dependence of calcite equilibrium constants was vali - dated by Hörbrand et al. (2018). As soon as the temperature decreases inside the heat exchanger, the SI turns to under- saturation (SI =−1.00 at 60 C ). This leads to a benign behavior of the geothermal water with respect to the precipitation of calcite, and no calcite scales are expected between the heat exchanger and the well head of the injection pipe. Description of scales Figure  1 shows a representative scale from 20 evaluated samples. The measured curva - ture of the samples resulting from the cylindric shape of the pipe is d ≥ 290mm for all samples. This rules out any smaller pipes, and indicates that the scales must have formed in the ground-level facilities, with pipe diameters between 316 and 348  mm (Table  1). At this point, we can safely rule out the production pipes ( 255mm ) as a source for the remobilized scales. Köhl et al. Geotherm Energy (2020) 8:23 Page 12 of 27 Fig. 4 Thin section of scale embedded with blue dye and observed with analysator. Whole thin section (left) and detail (right) with layers A and A , B and B 1 2 1 2 The ripples crests are oriented along the inner perimeter of the pipes, that is perpen - dicular to the flow in the pipes. This rules out accumulation of material which slides down the inner pipe walls at no flow conditions. The median of the measured ripple length is 9.7  mm (25  %-quartile (Q1): 9.0  mm, 75  %-quartile (Q3): 10.9  mm; N  =  31). The median of the measured ripple height is 1.2 mm (Q1: 1.1 mm, Q3: 2.0 mm; N = 31). The ripple morphology shows gently sloping stoss sides and steep lee sides (Figs.  1 and 4). This indicates that the ripples were formed in unidirectional current (current rip - ples). Also, the Ripple Symmetry Index (RSI) of the ripples is 2.9 (Q1: 2.5, Q3: 3.5) which indicates current ripples. This makes sense, because the facility is operating either in a unidirectional flow regime or no flow at all. Oscillation, like in wave ripples, does not occur even though some of the ripples show a RSI < 2.5 which is at the transition to orbital waves. The median bulk density of the ripple scales calculated from the geometry and weight 3 3 3 of the ripples is 1.3 g/cm (Q1: 1.2 g/cm , Q3: 1.4 g/cm ). Since the ripples consist of cal- cite ( 2.7 g/cm ), this indicates a large porosity in the scale fragments. More information is revealed by the thin sections of the ripple scales (Fig. 4). The scale samples have a median total thickness of 5.0 mm (Q1: 4.5 mm, Q3: 5.0 mm), and are all composed of four layers. The initial layer ( A , 2 mm thickness) consists of a fine grained calcite matrix with many pores which appear round or even circular. During prepara- tion of the samples oil came out of these pores. A dense layer made of calcite with large crystals ( B , 400µm thickness) covers layer A . It shows only very few pores which are 1 1 again filled with oil. Layers A (2 mm thickness) and B ( 400µm thickness) repeat layers 2 2 A B and but are phase-shifted against the underlying layer. 1 1 The ripple structure is created in layers A and A , respectively. The calcite layers  B 1 2 cover the ripple structures of the layers A with a constant thickness and like that stabi- lize the ripple structures. Scanning electron micrographs (Figs. 5, 6) show large pores in the A layers surrounded by small calcite crystals, and small pores but large calcite crystals in the B layers. Köhl  et al. Geotherm Energy (2020) 8:23 Page 13 of 27 Fig. 5 SEM image of the pipe side (layer A ) showing large pores surrounded by small calcite Fig. 6 SEM image of the water side (B ) showing small pores but large calcite crystals The SEM/EDX analysis showed that the porous matrix of the A  layers and the crys - talline layers B have the same elemental composition that points to a magnesium-rich calcite (Table 5). A few small (< 20µm ) pyrrhotite (FeS) and pyrite crystals ( FeS ) were found in the samples, particularly in the A layers. XRD analysis indicated magnesium calcite as dominant crystalline phase (48 at. %) fol- lowed by calcite (35  at.  %). Interestingly, kutnohorite CaMn(CO ) was detected with 3 2 17 at. %. The absence of kutnohorite in SEM/EDX might be due to a limited field-of-view or due to magnesium substitution of the manganese ions. For the crystal and pore size distribution, 100 - 180 measurements for each layer of the scales were conducted on the thin sections using ImageJ (Schindelin et al. 2012; Rueden Köhl et al. Geotherm Energy (2020) 8:23 Page 14 of 27 Table 5 Average EDX analysis of the ripple scales in atomic % A A B B 1 2 1 2 C 21.2 21.5 21.2 21.7 O 63.7 64.5 63.7 65.2 Na 0.1 0.0 0.3 0.0 Mg 0.5 0.2 0.9 0.5 Al 0.0 0.0 0.1 0.0 Si 0.0 0.0 0.1 0.1 S 0.0 1.5 0.0 0.5 Ca 14.4 11.8 13.7 12.1 Fe 0.0 0.3 0.0 0.0 Cu 0.0 0.2 0.0 0.0 Sum 100.0 100.0 100.0 100.0 Ca/Mg 24 32 15 23 Measurements 6 5 8 7 et al. 2017). The crystal size was determined by averaging the longest two visible crystal edges, and the pore size was determined by averaging the longest axis and its perpen- dicular axis. The results revealed that the pore sizes and crystal sizes differ between the A and B layers. The crystals of the A layers are significantly smaller (median: 22µm , Q1: 18µm , Q3: 28µm ) compared to the crystals of the B layers (median: 80µm , Q1: 59µm , Q3: 144 µm ). Here, some crystals grew distinctly larger (up to 500µm ) resulting a strong skewness of the distribution. The anisotropy of the crystals was calculated based on the evaluation of the transmis - sion micrographs, and is expressed as ratio of the shortest crystal edge to the longest crystal edge. With a median ratio of 0.8 (Q1: 0.7, Q3: 0.9) in the A layers and a median ratio of 0.7 (Q1: 0.6, Q3: 0.9) in the B layers the crystals show a similar, more or less iso- metric, only slightly prismatic crystal habit indicating spatially unconfined growth. The pores in the A layers from 180  measurements (median: 92µm , Q1: 71µm , Q3: 117µm ) are significantly larger compared to the pores in the B  layers (median: 57µm , Q1: 46µm , Q3: 68µm ). The largest pores in the A layers ( 250µm ) are twice as large as the largest pores in the B  layers. The micro-CT data indicate a median pore size of 50 - 60µm in all layers. Evaluation of the transmission micrographs (e.g. Fig.  7), allows a quantitative assess- ment of the porosity. The average area covered by matrix elements is 55 % leaving a pore space of 45 %. The 3-D micro-CT measurements yield a porosity for the A layers of 42 - 53 %, which is in good agreement with the 2-D measurements with TM. The porosity given for the B  layers is 0–20  %. Here, no measurements were possible with the trans- mission micrographs. The anisotropy of the pores expressed as the ratio of short and long axis of the ellipsoid shows a slight median anisotropy of 0.8 (Q1: 0.7, Q3: 0.9) for all layers. All pores look more or less round. Pores like this are known from fluid inclusions. Some degree of ani - sotropy can be attributed to static pressure, droplets coalescing and later growth of crys- tals into the pore space. The anisotropy is confirmed by the micro-CT measurements. Köhl  et al. Geotherm Energy (2020) 8:23 Page 15 of 27 Fig. 7 TM image showing spherical oil droplets surrounded by a thin calcite layer Mass balance for accumulated scales In order to understand the scaling process, the growth rates of the scales are important. The mass of the ripple scales will provide some evidence about the time frame in which the scales built up. The mass of scale present in the ground-level pipe can be calculated for each scale layer from the observed scale thickness and the known geometries of the pipes according to Eqs. 13 and 14: m = ρ · (1 − �) · V , calcite (13) −3 where m is the mass of scale, ρ is the density of calcite ( 2.71 gcm ),  is the poros- calcite ity with an average value of 45  % taken from the TM and micro-CT analysis and V is volume of scale given by V = x · d · π · 1/3 · l, (14) where x is the observed median scale thickness (5.0 mm, Q1 = 4.5 mm, Q3 = 5.0 mm), d is the ground-level pipe diameter (350 mm), and l is the pipe length between automatic 50 ± 5m filter and heat exchanger ( ). This is the section where all the scale fragments, which were later washed into the heat exchanger, originally formed, as proven by the scaling curvatures (see Section: Description of the scales). The factor of 1/3 accounts for the observation that only the lower third of the pipe showed scales as reported from a camera rover inspection of the pipes done shortly after the sampling date. Applying Eqs.  13 and  14 with a grain density of 2.71 ± 0.05 g/cm and a porosity of 0.45 ± 5 in the A  layers, and neglecting the minor porosity in the compact B  lay- ers, the mass of mineral scale attached to the pipe between the filter and the heat exchanger at the sampling date was roughly 154 ± 22 kg with the highest uncertain- ties resulting from porosity and scaled area. Furthermore, in case of a stillstand of the facility, the geothermal water volume of 4.8 m sitting in the ground-level facility could only precipitate 0.3 % of layer B , and 1 Köhl et al. Geotherm Energy (2020) 8:23 Page 16 of 27 is so excluding the possibility, that layer B formed during a short stillstand of the geo- thermal production. If we assume full thermodynamic equilibrium, the minimum volume of water required to precipitate the observed scales can be calculated. The analyses of the scales shows that the main component is calcite. With a molarity of calcium β of 0.416 mmol/L, and a SI of 0.24 the calcite scale potential, derived from the PhreeqC model during normal operation in front of the heat exchanger at 139.0 C and 12 bar is 0.130  mmol/L. If that total scale potential precipitated, a minimum volume of 4 3 1.2 · 10 m of geothermal water would be required to create the observed mass of scale of 154 kg on the substrate surface of 18.3 m . With the flow rate of 110 L/s dur - ing full operation, the time until this water volume is produced is about 30 hours. The time until layer B is formed would be about 4 hours. In reality, however, scale formation is limited by the kinetics of calcite precipita- tion and also controlled by the effective mineral surface. The average calcite scale rate at the well head for this specific site is 2.7 · µmol/(m · s ) as reported by Her- brich (2017). This rate is also in line with other sites (Wanner et al. 2017). From this real scale rate and from the production rate, it can be calculated that it would take 361 days during normal operation to create the observed 154 kg of calcite scale in the pipes, and the cementation of layer B would take 46 days. The conclusion of the mass calculation is that the observed ripple scales cannot have formed by new precipitation during regular operation within the 6 days since the last cleaning. Solid particle deposition versus in situ precipitation In this study the particles which built the ripples originally precipitated by heterogene- ous crystallization of calcite on the walls of the horizontal ground-level pipes. Earlier rover inspections have proven this scale growth on the pipe walls. This scale formed over the time of one year, between the last ground-level cleaning, and the scale sample collec- tion in October 2015. Precipitation was caused by degassing of CO in the pump which increases the pH and leads to supersaturation with respect to calcite as described in ear- lier studies (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). Evidence for degassing during operation was given by a glass window at the well head, in which tiny gas bubbles could be observed in the thermal water. It is known from various natural environments that ripples are the result of solid par- ticle deposition (Ha and Chough 2003; Kraemer and Winter 2016; Schieber et al. 2007; Soulsby et al. 2012; Tanner 1967; Yalin 1985). Therefore, the ripple shape itself which is visible in the A  layers of the scale directly implies that solid particle deposition of pre- existing particles was the driving process for the formation of these scale layers. In spite of that, it cannot be excluded that some minor additional in situ precipitation occurred during the settling of the particles. Furthermore, individual grains can easily be man- ually detached from the A  layers of the scale samples. This strongly suggests that the A  layers consist of weakly connected particles which settled down in a slow current, and did not precipitate in  situ. Additionally, the oil inclusions trapped in the A  layers indicate a fast deposition process. As shown in the mass balance section it is not pos- sible that scale of that thickness (2 mm A layers) precipitated fast enough in situ during the 6  days of regular operation since the last cleaning. Secondly, it is also unlikely that Köhl  et al. Geotherm Energy (2020) 8:23 Page 17 of 27 scale of that thickness precipitated during the maintenance works (2 - 3 days) since this would require very high supersaturation with respect to calcite. Maintenance includes cleaning with organic acid which makes long lasting supersaturation unlikely. Consider- ing all arguments together, the hypothesis that the rippled A  layers formed from pre- existing particles is much more evident than in  situ crystallization. This process is not yet described in literature for geothermal facilities, and is not accessible with current models. In contrast to the A layers, the B layers of the ripple scales are only 400µm thick, have a dense structure, and cover the A  layers evenly with constant thickness which means without superimposing another ripple structure on top of the A layers. These relatively thin B layers can already be created within a relatively short period of transient supersat- uration. Such short supersaturation can be reached at the end of the maintenance works, during flushing of the pipes with communal water of higher pH. This suggests that the B layers were actually formed by in situ precipitation during the maintenance works, in opposition to the A layers. Ripple formation—calculations The key to the understanding of the formation of the observed scalings is the presence of ripple structures. The establishment of current ripples requires sedimented but mobile particles and fluid flow. Sedimentation requires that the density of the mobile particles is significantly higher than the density of water or the particles will float or stay in suspen - sion. By the use of empirical equations, we reconstructed the conditions at which the ripples formed. With Eqs.  4 and  5, we calculated the theoretical ripple geometries with the observed features and fluid properties. For the latter, temperature-dependent changes of the viscosity have most impact. From the comparison to the measured ripple geometry, we could then derive the flow velocities in the system at which the ripples could have formed. In our system, regression equations  4 and  5 are valid for calcite particles ( ̺ = 2710 kg/m ) with a diameter between approx. 20 and 220µm at 139 C (operation temperature) and between 50 and 630µm at 20 C (maintenance temperature). Particles with lower particle density cause a shift to higher particle sizes. For an initial assessment 4  different scenarios for the ripple formation were defined (Table 8): Scenario 1 is standard operation at full production rate based on the observed calcite particle diameter ( d = 22µm ), scenario 2 is like scenario 1 but assumes floccu - lation/agglomeration to larger calcite particles ( d = 150µm ). Scenario 3 is represent- ative for the cleaning step with water and acid at a lower pumping rate assuming solid particle deposition of the observed calcite particle size ( d = 22µm ). Values for D 50 ∗ are 0.6 ± 0.2 for scenario 3, therefore, ripple height and length could not be calculated for this scenario. Scenario 4 assumes solid particle deposition of calcite-coated oil drop- lets ( d = 92µm ), also during maintenance. This scenario was introduced because the SEM images showed a thin layer of calcite surrounding each pore (Fig.  7). The process which could have created these coated oil droplets is known from literature as Picker- ing emulsion (Huang et al. 2019), where oil droplets serve as crystallization nuclei, onto which calcite precipitates. Then, the coated oil droplets would have become more dense Köhl et al. Geotherm Energy (2020) 8:23 Page 18 of 27 and experienced deposition. Also for this scenario D is < 1.2 and thus ripple height and length could not be calculated. Before applying Eqs.  4 and  5, we made a quick test by the use of Stokes’ law if par- ticles could have settled down at a reasonable time in the defined scenarios 1  -  4, or if the particles would have stayed in suspension. This test confirmed that in scenario  1 the sedimentation of calcite crystals from the top to the bottom of the pipe, would have taken only 3 min, so that solid particle deposition is reasonable. For larger particles (sce- nario 2), the sedimentation would have happened even faster. During maintenance (sce- nario 3) due to higher viscosity, the same sedimentation would have taken 13 min. If the ripples had been formed by calcite-coated oil droplets (scenario  4) ( ρ = 900 kg/m , oil 1.5µm calcite coating) sedimentation of these particles would have taken 19  min. The calculations for scenarios 1  -  4 by Stokes’ equation showed that sedimentation of the particles was fast enough to occur in all considered scenarios. Then, we applied Eqs.  4 and  5 to test which scenario is most likely. The results of the calculated ripple height and length are presented in Table 8 (upper part). With the meas- ured crystal sizes and composition the maximum ripple length for current ripples is around 31  mm for the operational scenario  1. Particle aggregation increases the ripple length to 83  mm (scenario  2). The ripple height in both cases is roughly 10  % of the ripple length. During maintenance (scenarios  3 and 4), the maximum ripple lengths are expected to be even longer, due to increased fluid viscosity, but equations for rip - ple length and height are not valid. To summarize the ripple geometries, all calculated values based on Eqs.  4 and  5 were either exceeding the measured values (  ≈ 10 mm , η ≈ 1.2 mm ) by at least a factor of 3, or out of calibration range. This indicates that the ripples were likely not formed from the individual crystals observed in layers A, or that the equations cannot be used in this setting. Why did the calculations according to Eqs.  4 and  5 predict higher and longer ripples than we observed, or in other words which parameters would need to be changed to get the observed small ripples? • Smaller ripples could be reached by lower viscosities (see Eqs. 4 -  6). However, lower viscosities can be excluded, and so viscosity does not explain the discrepancy. Vis- cosities were likely even higher than assumed, due to oil contained in the water. • Furthermore, a higher density contrast leads to smaller ripple dimensions, too. How- ever, for production scenarios 1 and 2, the density contrast was already at the upper limit, because density of non-porous calcite was assumed for the initial particles, and so this does not explain the small ripples. For maintenance scenario 4, an increase of the particle density to 1200 kg/m would allow using Eqs.  4 and  5, but the resulting ripple height and length would be 19  mm, respectively, 188  mm and therefore also way too long. The small ripples are therefore not explained by an erroneous density contrast, either. • However, the fine A  layers of the scale consist of mixture of cohesive and non- cohesive sediments. Such mixtures are reported to have significantly smaller ripple dimensions (height and length) due to cohesive forces (Baas et al. 2019, 2016). This could explain why our calculations based on experiments with pure non-cohesive material produced too large dimensions. For mixtures of cohesive and non-cohesive Köhl  et al. Geotherm Energy (2020) 8:23 Page 19 of 27 sediments somewhat smaller velocities for ripples formation would be expected (Baas et al. 2019; Schieber et al. 2007; Yawar and Schieber 2017). In contrast to the sophisticated Eqs. 4 – 6, the simple Eqs. 1 – 3 (see Table 8, lower part) which only depend on d showed relatively good agreement with the observed geom- etries even though they neglect other parameters. This suggests that these equations are still a good first guess if ripples are found in new settings. As a next step for the reconstruction of the ripple-forming conditions, the velocity was calculated according to Eqs.  6  -  12. This gave another hint that the ripples had not formed from individual crystals observed in layer  A: during full operation (scenario  1 and  2), the flow velocity in the pipes was well above the washout velocity of 0.52  m/s, and during maintenance (scenario 3) the critical velocity (0.21 m/s) for ripple formation from calcite particles was not reached (Table  8, upper part). For scenario  4, the calcu- lated flow velocities (0.04 – 0.23 m/s) matched the real flow conditions during mainte - nance, but since the equations for ripple length and height, which belong to the same regression of the data set, were not valid, also the velocity equations could strictly spo- ken not be applied for this scenario. Figures  8 and  9 show the calculated ripple forming velocities as a function of d for scenarios  1  –  3. The figure demonstrates the strong effect of temperature and viscos - ity on the shape of the curve U . While for 20 C , the washout velocity U washout washout increases with increasing particle size (Fig. 9), the washout velocity decreases with rising particle size, for a temperature of 139 C (Fig.  8). However, for the considered particle size ( 22µm ) of the ripple scale, the effect is small. Moreover, Figs.  8 and 9 show, which formation velocities of the ripples would be calculated for larger particles, as for instance created by flocculation. Since there is already a huge uncertainty for the viscosity in our system, the same applies for the calculated velocities. While oil on the surface of particles increases the viscosity, high temperatures lead to a drastic decrease in the viscosity. If the ripple scales had been formed during maintenance, low temperatures of about 20 C would have been present. Then, viscosity would have been 4 times higher than during normal operation ◦ ◦ at 139 C . Calculations for 20 C seemed in general better applicable than calculations at 139 C because the experiments with which the equations were set up were done at ambient temperatures. Discussion of the ripple‑forming scenarios Considering all calculations and aspects, when did the ripple scales, consisting of the rippled A  layers and the cemented B  layers, form? Table  6 summarizes the different options that are discussed in the following: • Ripple formation (A layers) could have happened during full production (scenario 1 and 2), even if it is not predicted by the calculations, but with the viscosity being out- side the calibration of Soulsby et al. (2012) ripples could possibly have formed even at the high flow velocity of 1–2 m/s. The observed small ripple dimensions can well be explained by cohesion of the silt size calcite crystals (Baas et al. 2019, 2016). This scenario seems unlikely, since no one reported ripples at these flow velocities before. Köhl et al. Geotherm Energy (2020) 8:23 Page 20 of 27 Table 6 Schematic overview of  different scale and  ripple formation mechanisms in  the  ground-level pipes between  filter and  heat exchanger with  indication of  the  calculated scenarios 1–4. No calculation scenarios were defined for  transition phases (start of  pump, shut down  of  pump) since  these could be ruled out  for  several reasons Flow rate, L/s Flow Pressure, bar Process velocity, m/s Start of production 110 1.40 12 Remobilization of existing scale particles due to turbulence. Pump rate is increased too fast to allow for ripple formation and stabilization Full production (scenarios 1 110 1.40 12 Slow scale precipitation due to and 2) degassing and oversaturation. Velocity > calculated ripple washout velocity, therefore no ripples predicted Shut-down of production 0.1 - 0.5 0.00 1 - 4 Solid particle deposition (rip- ples) of suspended particles followed by stabilization of scales due to sudden pressure drop and degassing. Too little particle supply Acid treatment (scenarios 3 5 0.06 1.8 - 2.6 Acid causes remobilization and 4) of partially dissolved scale particles and subsequent solid particle deposition (ripples). Oil prevents total dissolution of scales. Stabilization of ripples takes place when pH rises too high Flushing of facilities with tap 6.5 0.08 1.6 - 3.5 Stabilization of ripples by calcite water after acid treatment precipitation due to high cal- (scenarios 3 and 4) cium ion concentration in the cleaning solution and a neutral pH of the tap water Table 7 Course of  events during  the  3 months before  scale sample collection. Evaporator cleaning is  either  done by  flushing the  evaporator, or  by  opening it with  subsequent mechanical removal, or  by  chemical cleaning. T = temperature in  the  ground-level pre pipes before  interruption, T = temperature in  the  ground-level pipes at  the  end post of the interruption ◦ ◦ Date Cause Duration Flow rate, L/s T , C T , C Remark pre post Aug 3 Electric 11 d 0 140 25 Aug 17 Scales 3 h 110 139 126 Flushing of evaporator Aug 20 Scales 3 h 0 140 116 Mechanical removal Aug 24 Synchronization 2 h 0 141 120 Aug 28 Scales 2 h 110 139 132 Flushing of evaporator Sep 3 Scales 4 h + 2 h 0 140 118 Mechanical removal, 2 pump starts Sep 12 Revision 13 d 5 140 20 Chemical cleaning Oct 1 Scales 3 h 0 140 20 Sample collection If, despite of that, the A layers of the ripples formed during production, the B layers would have formed during different conditions (e.g. flushing with tap water during maintenance) due to their different compact appearance. Köhl  et al. Geotherm Energy (2020) 8:23 Page 21 of 27 Table 8 Comparison of  real observations and  conditions with  calculation results of  scenarios  1–4. The critical flow velocity U and  the  washout flow critical=min. velocity U were calculated by  Eq.  12 considering also  the  parameter ranges. The ripple height and  ripple length were calculated by  Eqs.  4 and  5 washout=max . considering also the parameter ranges. In the lower part of the table, ripple length and height predictions by the simple Eqs. 1–3 are indicated for comparison. Empty fields: unknown parameters, n.v. = equations not valid for the given input parameters (out of calibration) Observed ripple Observed conditions Observed conditions Scenario 1 (production Scenario 2 (production Scenario 3 (maintenance Scenario 4 geometries during production during maintenance 22μm calcite) 150μm flocc. calcite) 22μm calcite) (maintenance 92μm coated oil) T, C 139 20 139 139 20 20 ˙ 110 5 ± 2 110 110 5 ± 2 5 ± 2 −7 2 2.19 9.81 2.19 2.19 9.81 9.81 ν , 10 · m /s 2710 ± 0.14 2710 ± 0.14 2710 ± 0.14 2710 ± 0.14 2710 ± 0.14 1073 ± 28 ρ , kg/m 927 927 927 1020 ± 10 1020 ± 10 ρ , kg/m −6 22 (Q1: 18, Q3: 28) 22 (Q1: 18, Q3: 28) 150 22 (Q1: 18, Q3: 28) 92 (Q1: 71, Q3: 117) d , m · 10 D* 1.35 – 2.00 11.00 0.44 – 0.73 0.38 - 1.4 0.09 - 12 0.03 0.16 - 0.20 0.12 - 0.21 U , m/s 1.2 0 0.18 0.22 0.21 0.04 critical=min. U , m/s 2.2 0.08 0.52 0.33 0.58 0.23 washout=max . −3 1.2 3.1 - 3.9 8 n.v. n.v. η , m · 10 max −3 9.7 31 - 33 83 n.v. n.v. , m · 10 max −3 1.2 1.2 - 1.9 10 1.2 - 1.9 4.8 - 7.9 η , m · 10 max−flem. −3 9.7 8.6 - 12.0 42 8.6 - 12.0 24 - 35 , m · 10 max−colem. −3 9.7 18 - 28 150 18 - 28 71 - 117 , m · 10 max−yalin Köhl et al. Geotherm Energy (2020) 8:23 Page 22 of 27 Fig. 8 Boundaries for ripple-forming flow velocities U based on calcite particle size d (scenario 1 between ◦ −7 2 dashed lines). T = 139 C , ν = 2.19 · 10 m /s Fig. 9 Boundaries for ripple-forming flow velocities U based on calcite particle size d (scenario 3 between −7 2 dashed lines). T = 20 C , ν = 9.81 · 10 m /s Köhl  et al. Geotherm Energy (2020) 8:23 Page 23 of 27 • The second option is that ripples (A layers) formed during maintenance (scenario 3 and  4). The calculated velocities are higher than the maintenance velocities. At maintenance velocity, calcite particles would be below the critical threshold. Cal- cite-coated oil droplets fit to the flow velocities, but are unlikely regarding the rip - ple dimensions. However, ripple formation during maintenance is supported by the evidence of the 4  layer structure of the scale which fits well to the cleaning cycle consisting of 2 pumping phases, where the A layers could have formed, each one fol- lowed by a stillstand overnight, in which the B layers could form on top and stabilize the ripples at increasing pH due to degassing of CO . Alternatively, the B layers could have formed during flushing of the pipes at the end of the maintenance. Here, strong supersaturation was reached by high calcium ion concentrations in the fluid leading to calcite precipitation when the fluid was mixed to tap water with neutral pH. While this scenario assumed ripple formation (A  layers) during maintenance, the formation of the B layers could also have happened after the maintenance works, i.e. during normal production. However, this seems unlikely, since B  layer cementation would take 46  days during normal production. This fits neither to the 6  days since the last chemical cleaning, nor to the time span of one year lying between the yearly chemical cleanings. Furthermore, a cementation of the B layer during full production is unlikely because the B layers contain almost no oil inclusions which suggests grav- ity driven oil–water phase separation at quiescent conditions, not at high turbulence. • Another option for the formation of the ripple scales is the start phase of the pump. During the start phase, according to our calculations of the flow velocity, ripple formation was possible. Particles which had not fully been removed by earlier acid treatments of the facility, could have been remobilized and have formed ripples. This requires that the pump rate was only very slowly increased, so that the particles could stabilize by new precipitation of calcite before higher flow velocities could have lead to washout. Since the pump, at the investigated time period, was usually imme- diately turned on to almost final pump rate, this scenario seems unlikely. • Also the shut-down phase is a potential candidate for ripple scale formation. During the shut-down phase of the pump, suspended particles could have settle down in a ceasing current creating a rippled A layer. However, this scenario seems unlikely, as well, because there are not enough particles in suspension during normal produc- tion to create a continuous layer of the observed thickness (2 mm scaling layer A lay- ers in Fig.  4). The particle load in geothermal facilities in the Molasse Basin is only about 0.01  –  1  mg/L (Wolfgramm et  al. 2011), so that additional preexisting parti- cles (scales) would have needed to be remobilized to reach that thickness. Secondly, suspended particles during full operation consist of other particles, too (Wolfgramm et  al. 2011), and do not consist solely of calcite as observed in our study. Third, rip - ples need some time to be created (Soulsby et  al. 2012). During laboratory experi- ments, ripple formation only started after 40 min, while before that, rolling particles did not pile up significantly (Doucette and O’Donoghue 2006). Therefore ripple for - mation during shut-down phase is unlikely. To summarize, the calculated flow velocities contradict that the ripples were formed during full production and during maintenance. Additionally, the start- and shut-down Köhl et al. Geotherm Energy (2020) 8:23 Page 24 of 27 phase have the calculated flow velocities, but seem also unlikely for ripple formation for several above explained reasons. This suggests that the investigated ripple scalings in geothermal facilities did actually form at conditions which are not predicted by the empirical equations. This is possible since the ripples of this study were formed in a complex multiphase system (water, gas, oil) with simultaneous precipitation, dissolu- tion and degassing while the empirical equations are derived from controlled laboratory experiments. Major differences between the systems are that the empirical equations were derived from experiments with particles in the sand grain size while in this study particles are in the silt size. However, studies have reported ripples formed from silt and clay particles, as well. These fine grained ripples are formed by spontaneous floccula - tion of the particles which by that get large enough to settle down according to Stokes’ law (Schieber and Southard 2009). Furthermore, the particles in this study have a rhombohedral shape, while empirical experiments mostly studied ripples composed of rounded particles. The process of ripple formation comprising entrainment of particles, transient suspension, and reattachment as described by Ha and Chough (2003) is also possible with poorly rounded particles. Well rounded rolling particles are not required for ripple formation. Furthermore, rip- ples can also build from clay- and silt-size particles (Baas et al. 2019, 2016; Schieber and Southard 2009) which are typically also not rounded but rather platy. This shows that the principles of ripple formation can be applied to the rhombohedral particles of this study, as well. However, the different particle shape may explain some derivation from the quantita - tive theoretical prediction. The exposed surface area and the drag resistance of rhom - bohedral particles in a fluid current are different to rounded particles. This interplay of drag force and drag resistance due to the particle shape may influence the produced geometries and necessary ripple-forming velocities. However, this is yet little described in literature, and should therefore be addressed in future research. In literature (Soulsby et al. 2012; Baas et al. 2016) mostly only the particle diameter D is considered as criti- cal geometric parameter, but the shape of the particle is neglected in the equations. Conclusion Given all differences and uncertainties, we come to the conclusion that the ripples inves - tigated in this study formed during maintenance works at a very low flow rate (scenario 3 or 4). The silt size of the particles supports the small flow rates. Based on the evaluation, we conclude that the rippled scalings were created the following way. The facility was shut-down for the yearly maintenance. After displacement of the ther - mal water by tap water, acid was added to the solution which partially dissolved and destabilized the scaling, but did not reach all parts of the scale. During the acid- circula- tion, small scale fragments or individual crystals or rounded grains were transported, and underwent solid particle deposition creating ripples (layer A ). Later, layer B pre- 1 1 cipitated on top of layer A . There are two options, when and how layer  B formed: 1 1 • The first option is that layer  B formed overnight, at quiescent conditions, when the acid worked leading to degassing of CO and a rising pH exceeding the thermody- namical equilibrium, so that dissolution was followed by precipitation/cementation Köhl  et al. Geotherm Energy (2020) 8:23 Page 25 of 27 from the fluid with a very high calcium ion concentration. Like that, layer  B formed on top of layer A and stabilized it. The procedure was repeated the next pumping day and the following night resulting in layer A and B . It is striking that there are 2 2 2 pump phases and 2 stillstand phases during maintenance while there are also 2 rip- ple layers and 2 calcite cement layers. • The second option, how layer  B formed, is that it formed during flushing of the pipes at the end of the cleaning procedure. The fluid in the pipes has a high calcium ion concentration due to the dissolution of calcite. When this fluid gets mixed with fresh tap water of neutral pH, strong temporary supersaturation with respect to calcite is likely. If that is, what happened, then the observed ripple scales were created during 2 maintenance events, each event forming one A layer and one B layer, since flushing is only done once at the end of each cleaning cycle. Then, layer  A and B would have 1 1 formed before and during the cleaning one year ago, and subsequently would have withstood the last cleaning cycle. Rippled scales give insight into the hydrochemical conditions and fluid dynamics in the investigated geothermal facility, and might be a key to avoid cost-intensive shut-downs due to scales in the facility. According to our calculations based on Soulsby et al. (2012), ripple formation is nei- ther predicted for the production phase nor for the maintenance phase because the cal- culated velocities lie in between. However, major differences between this study and the experimental setup, from which the equations used for the calculations, were derived, suggest that ripples were created during maintenance, in spite of that. Particle diameter and fluid viscosity are major parameters that differ between this study and the empirical predictions. Therefore, we conclude that solid particle deposition of mobile calcite crystals at the bottom of the pipe created the rippled scalings during maintenance works. We are con- fident that due to a protective oil layer on an already existing scale, the acid used during the maintenance was not able to entirely dissolve the scale. However, the acid weakened and mobilized small particles from the scale which were redeposited as ripples in a weak current. Here, either pH increase overnight or a small oversaturation of the tap water at the end of the cleaning was sufficient to stabilize the ripples by precipitating a tiny layer of calcite on top. In general, oil plays an important role in the formation of calcite scale. It can mitigate the initial phase of scale formation by making the attachment to the pipe harder, but it can also somewhat protect an existing scale from acid. Furthermore, operational shut- downs tend to create oil layers by promoting water–oil phase separation. Finally, the importance of removing the oil from the scale before and during acid treatment must be emphasized. Only when you know that all scale can be accessed by the acid and a steady pH is reached with no further acid addition the facility can be considered clean. Abbreviations SEM/EDX: Scanning electron microscope with energy-dispersive X-ray spectroscopy; SI: Saturation Index. Acknowledgements The financial support by the Bayerisches Staatsministerium für Wissenschaft und Kunst is gratefully acknowledged. We thank the Stadtwerke München for providing us the access to their geothermal facility, to their scanning electron Köhl et al. Geotherm Energy (2020) 8:23 Page 26 of 27 microscope and laboratory and for supplying us with the operational data. We would very much like to thank Prof. Dr. Martin Elsner and Dr. Martina Ueckert for multiple readings and discussions of the manuscript. Furthermore we thank Christian Schulbert from the University of Erlangen for the micro-CT measurements and visualization. Vladimir Ruttner from the Chair of Engineering Geology at TUM is gratefully acknowledged for the preparation of the thin sections. We also want to thank Amanada Günther for the XRD-measurements and Christine Benning for additional SEM/EDX measurements. Authors’ contributions BK collected the water, gas and ripple scale samples, analyzed the scales, performed the PhreeqC calculations, checked the ripple reconstruction calculations of JG and discussed the formation scenarios with JG and TB. JG performed the rip- ple reconstruction calculations, discussed the scenarios and assisted during the water and gas sampling. TB had the idea, developed the hydrogeochemical modeling framework and discussed the calculations and scenarios with BK and JG. All authors read and approved the final manuscript. Funding This work was partially funded by the Bavarian Ministry of Science and Art in the project Geothermie-Allianz Bayern. Open access funding provided by Projekt DEAL. Availability of data and material The datasets used and analyzed during the current study are available from the corresponding author on reasonable request. Compliance with ethical standards Competing interests The authors declare that they have no competing interests. Author details 1 2 Institute of Hydrochemistry, Technische Universität München, Marchioninistrasse 17, 81377 Munich, Germany. Depart- ment of Civil, Architectural, and Environmental Engineering, The University of Texas at Austin, 301 E Dean Keeton St., 78712 Austin, TX, USA. Chair of Hydrogeology, Technische Universität München, Arcisstr. 21, Munich 80333, Germany. Received: 8 March 2020 Accepted: 4 August 2020 References Agemar T, Alten J-A, Ganz B, Kuder J, Kühne K, Schumacher S, Schulz R. The Geothermal Information System for Germany - GeotIS. Zeitschrift der Deutschen Gesellschaft für Geowissenschaften. 2014;165(2):129–44. https ://doi. org/10.1127/1860-1804/2014/0060. Allaby M, editor. A Dictionary of Earth Sciences. 3rd ed. : Oxford University Press, online version; 2008. www.oxfor drefe rence .com/view/10.1093/oi/autho rity.20110 80310 04223 07. Avila K, Moxey D, de Lozar A, Avila M, Barkley D, Hof B. The onset of turbulence in pipe flow. Science. 2011;333(6039):192– 6. https ://doi.org/10.1126/scien ce.12032 23. Baas JH. A flume study on the development and equilibrium morphology of current ripples in very fine sand. Sedimen- tology. 1994;41:185–209. https ://doi.org/10.1111/j.1365-3091.1994.tb014 00.x. Baas JH, Best JL, Peakall J. Predicting bedforms and primary current stratification in cohesive mixtures of mud and sand. J Geol Soc. 2016;173(1):12–45. https ://doi.org/10.1144/jgs20 15-024. Baas JH, Baker M, Malarkey J, Bass S, Manning A, Hope J, Peakall J, Lichtman ID, Ye L, Davies A, Parsons D, Paterson D, Thorne P. Integrating field and laboratory approaches for ripple development in mixed sand-clay-eps. Sedimentol- ogy. 2019;66:2749–68. https ://doi.org/10.1111/sed.12611 . Bauer M, Freeden W, Jacobi H, Neu T. Handuch Tiefe Geothermie. Berlin: Springer; 2014. Baumann T. Validation of hydrochemical analyses and gas concentrations of deep geothermal aquifers. In: 41st Workshop on Geothermal Reservoir Engineering. Stanford: Stanford University; 2016. p. 1269–77. https ://pange a.stanf ord.edu/ ERE/pdf/IGAst andar d/SGW/2016/Bauma nn.pdf. Berg JH, Gelder A. A new bedform stability diagram, with emphasis on the transition of ripples to plane bed in flows over fine sand and silt. God. 2009;. https ://doi.org/10.1002/97814 44303 995.ch2. Birner J, Mayr C, Thomas L, Schneider M, Baumann T, Winkler A. Hydrochemistry and evolution of deep groundwaters in the Malm aquifer in the bavarian part of the South German Molasse Basin. Z Geol Wiss. 2011;39(3/4):291–308. Boch R, Leis A, Haslinger E, Goldbrunner J, Mittermayr F, Froeschl H, Hippler D, Dietzel M. Scale-fragment formation impairing geothermal energy production: interacting H2S corrosion and CaCO3 crystal growth. Geotherm Energy. 2017;5:1–19. https ://doi.org/10.1186/s4051 7-017-0062-3. Coleman SE, Fedele JJ, Garcia MH. Closed-conduit bed-form initiation and development. J Hydraul Eng. 2003;129(12):956–65. https ://doi.org/10.1061/(ASCE)0733-9429(2003)129:12(956). Doucette JS, O’Donoghue T. Response of sand ripples to change in oscillatory flow. Sedimentology. 2006;53(3):581–96. https ://doi.org/10.1111/j.1365-3091.2006.00774 .x. Flemming B. On the classification of subaquatic flow-transverse bedforms. 1988;29. Freimann R, editor. Hydraulik Für Bauingenieure. Grundlagen und Anwendungen. München: Carl-Hanser; 2014. Furtak H, Langguth H-R. Zur hydrochemischen Kennzeichnung von Grundwässern und Grundwassertypen mittels Ken- nzahlen. IAH-Congress. 1965;7:89–96. GmbH VG. VGSTUDIO MAX 3.0. 2019. https ://www.volum egrap hics.com/en/produ cts/broch ure-downl oad.html. Köhl  et al. Geotherm Energy (2020) 8:23 Page 27 of 27 Ha HK, Chough SK. Intermittent turbulent events over sandy current ripples: a motion-picture analysis of flume experi- ments. Sediment Geol. 2003;161:295–308. https ://doi.org/10.1016/S0037 -0738(03)00132 -5. Herbrich M. Anlagen- und sicherheitsrelevante Prozesse im Thermalwasserkreislauf eines geothermischen Heizkraft- werks. PhD thesis, Technische Universität München. http://nbn-resol ving.de/urn/resol ver.pl?urn:nbn:de:bvb:91- diss-20161 209-13102 43-0-3 Hörbrand T, Baumann T, Moog HC. Validation of hydrogeochemical databases for problems in deep geothermal energy. Geotherm Energy. 2018;. https ://doi.org/10.1186/s4051 7-018-0106-3. Huang F, Liang Y, He Y. On the pickering emulsions stabilized by calcium carbonate particles with various morphologies. Coll Surf. 2019;580:123722. https ://doi.org/10.1016/j.colsu rfa.2019.12372 2. Kraemer K, Winter C. Predicted ripple dimensions in relation to the precision of in situ measurements. Ocean Sci Discuss. 2016;. https ://doi.org/10.5194/os-2016-20. Langmuir D. The geochemistry of some carbonate ground waters in central Pennsylvania. Geochimica et Cosmochimica Acta. 1971;35:1023–45. https ://doi.org/10.1016/0016-7037(71)90019 -6. Mayrhofer C. Hydrochemische Untersuchungen im Malmaquifer im bayerischen Molassebecken. Dissertation, Technis- che Universität München, München. 2013. https://mediatum.ub.tum.de?id=1128419 Parkhurst DL, Appelo CAJ. Description of input and examples for PHREEQC version 3 – a computer program for specia- tion, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Technical Report Techniques and Methods 6–A43, U.S. Geological Survey. 2013. http://pubs.usgs.gov/tm/06/a43/. Rotta J. Experimenteller Beitrag zur Entstehung turbulenter Strömung im Rohr. Ingenieur-Archiv. 1956;24:258–81. Rueden CT, Schindelin J. Imagej for the next generation of scientific image data. BMC Bioinform. 2017;. https ://doi. org/10.1186/s1285 9-017-1934-z. Schade H KE. Strömungslehre. Berlin New York: de Gruyter; 1989. Schieber J, Southard JB. Bedload transport of mud by floccule ripples-direct observation of ripple migration processes and their implications. Geology. 2009;37:483–6. https ://doi.org/10.1130/G2531 9A.1. Schieber J, Southard J, Thaisen K. Accretion of mudstone beds from migrating floccule ripples. Science. 2007;318(5857):1760–3. https ://doi.org/10.1126/scien ce.11470 01. Schieber J, Southard JB, Kissling P, Rossman B, Ginsburg R. Experimental deposition of carbonate mud from moving suspensions: Importance of flocculation and implications for modern and ancient carbonate mud deposition. J Sediment Res. 2013;83(11):1025–31. https ://doi.org/10.2110/jsr.2013.77. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Tinevez J-Y, White DJ, Hartenstein V, Eliceiri K, Tomancak P, Cardona A. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;. https ://doi.org/10.1038/nmeth .2019. Soulsby RL, Whitehouse RJS, Marten KV. Prediction of time-evolving sand ripples in shelf seas. Contin Shelf Res. 2012;38:47–62. https ://doi.org/10.1016/j.csr.2012.02.016. Southard JB, Boguchwal LA. Bed configuration in steady unidirectional water flows; Part 2, Synthesis of flume data. J Sedi- ment Res. 1990;60(5):658–79. https ://doi.org/10.1306/212F9 241-2B24-11D7-86480 00102 C1865 D. Tanner WF. Ripple mark indices and their uses. Sedimentology. 1967;9:89–104. https ://doi.org/10.1111/j.1365-3091.1967. tb013 32.x. Wanner C, Eichinger F, Jahrfeld T, Diamond LW. Causes of abundant calcite scaling in geothermal wells in the Bavarian Molasse Basin, Southern Germany. Geothermics. 2017;70:324–38. https ://doi.org/10.1016/j.geoth ermic s.2017.05.001. Wolfgramm M, Rauppach K, Thorwart K. Mineralneubildungen und Partikeltransport im Thermalwasserkreislauf geother- mischer Anlagen Deutschlands. Z Geol Wiss. 2011;39:213–39. Yalin MS. On the determination of ripple geometry. J Hydraul Eng. 1985;111(8):1148–55. https ://doi.org/10.1061/ (ASCE)0733-9429(1985)111:8(1148). Yawar Z, Schieber J. On the origin of silt laminae in laminated shales. Sediment Geol. 2017;360:22–34. https ://doi. org/10.1016/j.sedge o.2017.09.001. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Geothermal Energy Springer Journals

Rippled scales in a geothermal facility in the Bavarian Molasse Basin: a key to understand the calcite scaling process

Loading next page...
 
/lp/springer-journals/rippled-scales-in-a-geothermal-facility-in-the-bavarian-molasse-basin-Rx2dHBnsi1
Publisher
Springer Journals
Copyright
Copyright © The Author(s) 2020
eISSN
2195-9706
DOI
10.1186/s40517-020-00177-6
Publisher site
See Article on Publisher Site

Abstract

bernhard.koehl@tum.de Institute of Hydrochemistry, Mineral precipitates, scales, are wide-spread in geothermal facilities and hamper Technische Universität technical and economic efficiency. Scales reported from earlier studies in the Bavar - München, Marchioninistrasse 17, 81377 Munich, Germany ian Molasse Basin consisted mostly of calcite, appeared rather homogeneous, grew Full list of author information flat and perpendicular to the surface, and showed little or no internal macrostructure is available at the end of the apart from small gradients in crystal size and chemical composition. In this study, we article investigated extraordinary calcite scales collected from the inlet chamber of a heat exchanger. These scales showed current ripples and a sequence of coarse and dense layers. By an extensive investigation of the ripples, the hydrochemistry, and by the application of existing empirical equations, the formation of the ripple scales could be reconstructed: Due to degassing of CO in the geothermal pump calcite scales precipitated on the pipes of the ground-level facilities. During maintenance works, these scales were partially dissolved, and remobilized as individual mobile grains. These particles subsequently settled again, and solidified as ripple scales in the horizontal pipes at surface level. For this ripple forming process, flow had to be small, because otherwise the ripples would have been altered or washed away. After deposition of the ripple structure, a dense new calcite layer precipitated on top of the ripples, and stabi- lized it due to transitory supersaturation. This sequence of ripples and a dense calcite cover layer is repeated two times. Finally, during start-up of regular plant operation, the whole composite scale was remobilized and got washed into the inlet chamber of a heat exchanger. The investigation presented here indicates that this type of calcite scale only forms in geothermal facilities in the Bavarian Molasse Basin if certain flow and hydrochemical conditions are present during the maintenance procedure. Keywords: Scales, Scaling, Precipitations, Calcite, CaCO , Ripples, Malm, Bavarian Molasse Basin, Carbonaceous aquifer, PhreeqC Introduction Mineral precipitation and build-up in pipes, also known as scales or scalings, are a com- mon problem for geothermal plants in the Bavarian Molasse Basin. Growth of scale in pipes reduces the technical and economic efficiency of geothermal plants and requires frequent, costly maintenance or even process modification of installed geothermal plants. Mineral build-up can occur on any surface exposed to geothermal fluid, including © The Author(s) 2020. 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/. Köhl et al. Geotherm Energy (2020) 8:23 Page 2 of 27 the submersible pump, its motor, the production pipes, and ground-level components, such as heat exchangers (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). In previous studies in the Bavarian Molasse Basin scales had a more or less even sur- face and the scales were interpreted to form during normal operation of the facility (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). The scale samples were analyzed by a variety of chemical and mineralogical methods (SEM/EDX, XRD, acid digestion, polarizing microscopy) (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). The mor - phology of the scales, however, received little attention, mainly because the observed scales grew evenly on the pipes. The scales presented in this study show a distinct ripple pattern, which was not reported anywhere else (Fig.  1). The macroscopic structure was used to gain insight into the formation process and to develop countermeasures against this type of min- eral precipitations. The geothermal facility investigated in this study is located in the Bavarian Molasse Basin, and exploits Upper Jurassic sedimentary rocks, also known as Malm aquifer. This formation, especially in its dolomitic karstified facies and its massive reef facies hosts high permeabilities and offers favorable conditions for geothermal use (Bauer et al. 2014). The formation is situated at a depth of 2000 m north of Munich and up to 5000 m in the South. Flow rates reach up to 150 L/s and temperatures up to 150 C (Agemar et al. 2014; Bauer et al. 2014). The formation of scalings is controlled by the hydrochemical conditions. Therefore, we give a brief overview of the conditions in the target reservoirs. Within the Bavarian Molasse Basin, six hydrogeochemical provinces can be distinguished (Birner et al. 2011). Following this classification, the Munich area belongs to the central part of the Molasse Basin and shows a Na–Ca–HCO –Cl-type or Na–HCO –Cl-type water chemistry after the classifica - 3 3 tion of Furtak and Langguth (1965). South of Munich, the influx of waters with higher min - eral content from the overburden increases (Mayrhofer 2013). In the central Molasse Basin, Fig. 1 Ripple scale sample collected from the inlet of the heat exchanger. The sample has a calcite composition. Some oil contained in the thermal water is incorporated in the scale causing the black color. The shape of the ripples with a luv and lee side points to current flow ripples Köhl  et al. Geotherm Energy (2020) 8:23 Page 3 of 27 the Malm aquifer has high gas contents up to 2 L gas per L water. The gas is mainly com - posed of CH , CO , and N . H S concentrations in the gas phase can reach up to 5 vol %. 4 2 2 2 Under reservoir conditions, due to the long residence time the water is in chemical and thermal equilibrium with the surrounding rock (Wanner et  al. 2017). Here, the key equi- librium is the carbonate equilibrium, and the driving factors are temperature and pres- sure. During the production of the thermal water, the equilibrium is easily disrupted. Local increase of the temperature and degassing of CO cause oversaturation with respect to car- bonates (Wanner et al. 2017). On the other side, the cooled water injected into the reservoir shows a benign behavior with respect to scales due to the reduced temperature. In contrast to any other scale samples reported from geothermal sites in the Bavarian Molasse (Wanner et  al. 2017, Herbrich 2017, Boch et  al. 2017), the scales in the ground- level facilities of the investigated geothermal power plant showed a distinct morphology (Fig.  1) with ripple structures usually known from river beds or sea shores. In general, ripples are formed when either wave action (wave or orbital ripples) or current flow (cur - rent ripples) over morphological perturbations of the bedform lead to abnormalities in local flow conditions, such as vortices (Ha and Chough 2003). This results in local areas of increased and decreased shear velocity where sediments are eroded from or re-deposited onto the bedform, respectively. An example of flow regimes along with erosion and deposi - tion mechanisms near ripple structures is shown in Fig. 2. Geometric indices can be used to determine whether ripples were formed by waves or currents (Tanner 1967). The ripple symmetry index (RSI) is the ratio of stoss length to lee length. A RSI < 2.5 is indicative of wave ripples, whereas RSI > 3.0 is pointing to current ripples (Allaby 2008). The boundary conditions at which the scales under investigation have been formed are ruling out oscillation ripples (wave ripples). The formation of ripples is governed by particle diameter, velocity, fluid viscosity, and relative density of the sediment to the fluid. Although ripples are not static but in con - stant motion, the shape of the ripples, that is the periodicity of the ripples (wavelength), the height of the ripple crests, and the RSI, to name a few reaches an equilibrium state. The ripple-forming parameters have different effects: While experiments have shown that the equilibrium morphology of current ripples depends on the particle grain size and their immersed density (Kraemer and Winter 2016), the time to reach a stable ripple pattern depends on fluid density and fluid velocity (Baas et al. 2019). Many studies have attempted to derive empirical relationships for transient and steady- state ripple morphology (ripple height, η (m), ripple length,  (m)) based on median max max particle diameter ( d ) ( Y alin 1985; Flemming 1988; Baas 1994; Coleman et al. 2003; Soulsby et al. 2012) (Eqs. 1–5). Median particle sizes used in these studies range from 0.045 mm to 0.87 mm, the particles were either glass spheres or sand, and experiments were often imple- mented in an open-channel rectangular flume. =1000 · d max (Yalin) 50 (1) η =0.0677 · 1000 · d max (Flemming) 50 (2) 0.75 =175 · (d ) (3) max (Coleman) 50 Köhl et al. Geotherm Energy (2020) 8:23 Page 4 of 27 Fig. 2 Schematic drawing of ripples with indication of ripple length, height and key formation processes modified after Ha and Chough (2003) The range of flow velocities that result in ripple formation for a given particle size has also been studied (Soulsby et al. 2012). The minimum flow velocity is given by the shear forces which have to exceed the force required to erode the sediment bed. At flow veloci - ties lower than the minimum the sediment grains simply stay put. There is also a maxi - mum flow velocity for ripples, again depending on the grain size. Above this velocity the particles are washed out and ripples are destroyed, and the bedform exhibits sheet flow behavior. If ripples have been formed at suitable flow velocities, and the velocity later drops below the minimum threshold, the ripple structure in the bedform is conserved (Soulsby et al. 2012). Flow velocities resulting in ripples can also be read-off from several bedform phase diagrams which confine the boundary conditions for ripples and other bedforms (Southard and Boguchwal 1990; Berg and Gelder 2009; Baas et al. 2016). For instance, the diagram of Southard and Boguchwal (1990) (Fig.  3) shows the bed sedi- ment size against flow velocity for flow depths between 0.25 and 0.4 m. At these condi - tions, as they occur in natural systems, and with pure non-cohesive sediment, ripples can form at velocities of roughly 0.1 - 0.9 m/s, depending on sediment particle size. For mixtures of cohesive and non-cohesive sediments, the equilibrium morphol- ogy of current ripples appears to follow the same empirical relationships that were found for non-cohesive sediments, but the time required to reach the equilibrium morphology is longer (Baas et  al. 2019; Schieber and Southard 2009; Schieber et  al. Köhl  et al. Geotherm Energy (2020) 8:23 Page 5 of 27 ◦ ◦ Fig. 3 Bedform diagram showing the 10 C-equivalent mean bed sediment size against 10 C-equivalent mean flow velocity for 10 C-equivalent flow depths between 0.25 and 0.4 m after Southard and Boguchwal (1990). Indication of the equivalent temperature is needed due to the effect of temperature on fluid viscosity and density. The diagram shows that current ripples are likely to occur at mean flow velocities between roughly 0.1 and 0.9 m/s and at a sediment size of roughly 0.1–0.8 mm 2013). When cohesive sediments form ripples, the particles flocculate, which enables them to migrate as bedload creating ripples. Floccules composed of silt and clay are typically in the range of 100 - 200µm (Yawar and Schieber 2017), sometimes several 100µm (Schieber et  al. 2007). To our knowledge the shortest ripples reported in lit- erature can be found in Baas et al. (2019) who observed wavelengths down to 77µm . These ripples consisted of a mixture of cohesive sand and 13  % mud fraction. Flow velocities leading to the formation of mixed cohesive and non-cohesive ripples are in the range of 0.1 - 0.7 m/s, but typically lower than velocities needed for pure non- cohesive sediments (Baas et al. 2019; Schieber et al. 2007; Yawar and Schieber 2017). If the cohesive portion exceeds about 30  %, no ripples form and a flat bed remains (Baas et al. 2016). Baas et al. (2016) published a review on the available equations for the prediction of rip- ples. Baas found that Soulsby et  al. (2012) created the most comprehensive predictor for current generated and combined-flow bedforms capable of predicting ripple geometry and ripple-forming velocity. Therefore, we are applying these equations (Soulsby et al. 2012) for the calculations in this paper. Equations 4 and 5 predict the maximum ripple height, η max (m), and maximum ripple length,  (m), for a given mean diameter d (m) and dimen- max 50 sionless relative particle density s, which is described by the dimensionless parameter D defined in Eq. 6 : −0.554 η =202d D | , (4) max 50 1.2<D <16 ∗ ∗ −1.5 =d (500 + 1881D )| , (5) max 50 1.2<D <16 ∗ ∗ Köhl et al. Geotherm Energy (2020) 8:23 Page 6 of 27 1/3 g(s − 1) D =d , (6) ∗ 50 where g is the gravitational constant ( m/s ), s = ρ /ρ is the dimensionless ratio of p f particle density (ρ ) to fluid density ( ρ ), and ν (m /s) is the kinematic viscosity of the p f fluid. Flow regimes favoring ripple formation can be determined by comparing the non- dimensional Shields number θ (Eq. 7) to the critical Shields number θ (Eq. 10) and the c cr washout Shields number θ (Eq. 11) (Soulsby et al. 2012): wo θ <θ no bedform movement, c cr θ < θ <θ ripple formation, cr c wo θ <θ ripples are washed out, resulting in either washed-out ripple or sheet bed wo morphology: C U θ = , (7) g(s − 1)d 0.40 C = , (8) ln(h/z ) − 1 z =d /12, (9) 0 50 0.3 θ = + 0.055[1 − exp(−0.02D )], cr ∗ (10) 1 + 1.2D 0.916, if D ≤ 1.58 ′ ∗ θ = (11) wo −1.3 1.66D , if D > 1.58, where U (m/s) is the average cross-sectional flow velocity in the pipe, C is the non- dimensional drag coefficient due to the particles (Eq.  8), h (m) is the water depth (for a full pipe this is the pipe diameter), and z (m) is the grain-related bed-roughness length (Eq. 9). Combining Eqs. 7, 8, and 9, and solving for U , one arrives at an expression for aver- age fluid velocity as a function of the Shields number (Eq.  12). Plugging values for θ cr and θ into Eq.  12 gives the lower and upper average fluid velocities, respectively, at wo which ripples would form for a given d and s. 1/2 U = 2.5(ln(12h/d ) − 1) θ g − 1 d (12) 50 c 50 For the characterization of the flow in the geothermal pipes investigated in this paper, Reynolds numbers can be used. Transition from laminar flow to turbulent flow occurs at Reynolds numbers of 2300 (below laminar) to 2900 (above fully turbulent) (Rotta 1956; Schade 1989; Avila et  al. 2011; Freimann 2014). The Reynolds numbers in the investi - gated facility (see Table  1) are on the order of Re = 10 and well in the fully turbulent flow regime. Köhl  et al. Geotherm Energy (2020) 8:23 Page 7 of 27 Table 1 Pipe dimensions and flow conditions during operation, and *) during maintenance ◦ ◦ Tubing ID, mm ˙ v, m/s Re 139 C Re 20 C V , L/s DN 245 250 255 110 2.15 - 2.51 · 10 DN 300 300 316 110 1.40 - 2.02 · 10 DN 350 350 348 110 1.16 - 1.84 · 10 DN 350 *) 350 348 4 - 8 0.04 - 0.08 - 1.5 −2.9 · 10 In this study, we investigate whether the framework of ripple formation provides access to the yet unknown flow conditions which led to scale formation in the surface- level pipes of the geothermal facility. A camera rover inspection of the pipes conducted shortly after the sampling date showed only scales at the bottom part of the horizontal or slightly slanted pipes indicating deposition. We assume that the shape of the scales will allow to differentiate between regular operation (110 L/s, respectively, 1.2 –2.2 m/s, depending on the pipe diameter) and maintenance operation (0–8  L/s, respectively, 0–0.08  m/s). We also have to assume at least one step to cement the particles which make up the scales and form a solid matrix. When this cementation took place is of great interest to the operators to prevent the formation and further displacement of this type of calcite scales in geothermal facilities in the Bavarian Molasse Basin. Materials and methods Site description and operation The site under investigation produces geothermal water from the carbonate reservoir in the sediments of the Upper Jurassic in the Bavarian Molasse Basin south of Munich. The temperature at the well head is 139 C and the typical flow rate is 110  L/s. During production the pressure decreases from about 300 bar to 12 bar ground-level pressure. Table 1 gives the dimensions, flow velocities, and Reynolds numbers for the major pipes in the facility. The pipes are entirely filled with water so that water depth of the horizon - tal pipes (DN 300 and DN 350) equals pipe diameter. Reynolds numbers were calculated ◦ ◦ for a pressure of 12 bar and temperatures of 139 C and 20 C. The calculated Reynolds numbers are on the order of 1 · 10 for full production in front of the heat exchanger. During maintenance, the Reynolds numbers are still in the 1 · 10 range. The numbers indicate turbulent flow in all pipes, and at all flow conditions. In order to prevent detached scales from the production pipes, rock fragments from the reservoir, and left-over cuttings from the drilling of the well to enter the produc- tion plant, particle filters with a cutoff of 50µm are installed close to the well head. Two automatic filters are operating in parallel. Any larger particles detected behind the filters have likely precipitated in the surface-level installations or grew from smaller particles. Samples were collected after a series of interruptions of the geothermal energy pro- duction (see Table  7) in Oct. 2015. One of the automatic particle filters was in service during the majority of this period. There were two maintenance cycles to mechanically remove scales from the inside of the heat exchanger (Aug 20, Sep 3). Here, the electric submersible pump (ESP) was kept running, bypassing the heat exchanger. All other inci- dents came with a shutdown of the ESP, hence the flow velocity in the facility was zero. Köhl et al. Geotherm Energy (2020) 8:23 Page 8 of 27 Unexpectedly, right after the yearly maintenance, huge amounts of scales blocked the heat exchanger, and caused another interruption (Oct 1). Cleaning of ground‑level facilities The procedure for the chemical cleaning of the ground-level pipes comprised the follow - ing steps: First, the ESP was shut-off, and the ground-level facilities were disconnected from the production and injection pipes. Second, the geothermal water in the surface- level pipes was displaced by tap water. Then, 70 % methanesulfonic acid and a surfactant solution with a pH around  7 were gradually added to the tap water in the pipes until the concentration of the acid in the fluid reached the desired 5  % at a dilution ratio of about 1:17, and a pH of 0 - 1. The surfactant was needed to emulsify oil, which is con - tained in the geothermal water, and to improve interaction of scales and acid. The clean - ing solution was continuously pumped through the facility at a flow rate of 4 - 8 L/s for approximately 8 h. The solution was then allowed to sit quietly overnight. The following day, the same solution was again continuously pumped through the facilities for 8 h of flow and allowed to sit overnight. Starting the next day the acid was drained from the facilities into a neutralization tank. The ground-level pipes were filled with tap water before normal production was restarted. The pressure during the entire cleaning oper - ation was held at 1.2 bar to 3 bar by a nitrogen gas pressure system. In total, 1,800  L of concentrated acid along with 2,000 L surfactant solution were consumed during the acidification. Methods Collection and analyses of water and gas samples Water and gas samples were taken at regular intervals before and after the cleaning procedure, and showed little variance (i.e. ± 4% TDS Herbrich 2017). The samples pre - sented here were taken after the incident and after full operation resumed. Samples were taken at a sampling port between well head and the automatic filter. An on-site cooling device was used to bring the temperature to 54 C to be able to handle the samples safely. Temperature, pH-value, electric conductivity, and redox potential were measured on-site with portable electrodes and sensors (InLab Expert Go-ISM, InLab 738-ISM, InLab Redox ORP, Mettler Toledo, Germany). Alkalinity and acidity were analyzed by mass analysis on-site with 0.1 M HCl (Carl Roth), 0.1 M NaOH (Carl Roth), and using colorimetric indicators phenolphthalein (Merck) and cooper (Sigma − − − 2− Aldrich). Anion ( F , Cl , Br , SO ) concentrations were quantified by IC (Dionex, IC25). Total sulfide in water was precipitated on-site with 10 mL zinc acetate (2 %). The concentration of the precipitated zinc sulfide was measured with the methylene blue colorimetric method (Deutsches Institut für Normung, 1995) using a spectrophotom- + + eter (Double Beam UV-190, Shimadzu). Na and K were quantified using a flame pho - tometer (Eppendorf ELEX 6361), calcium was determined by titration with EDTA. The other concentrations of the other cations were determined by AAS (PerkinElmer, 3300). The concentration of dissolved gases in the thermal water loading was quantified using an on-site degassing unit as described by Baumann (2016). The unit was composed of a 9−L plastic tank with two connectors in the lower third of the tank serving as inlet and outlet of the thermal water. A third connector at the top of the tank was used to transfer Köhl  et al. Geotherm Energy (2020) 8:23 Page 9 of 27 the accumulated gas into a gas cylinder (1  L gas). At the beginning the sampling tank was filled completely with thermal water with the outlets for water and gas open. After a constant flow of water was established and measured, the gas outlet was closed. The gas which then accumulated inside the container was measured as a function of time using an imprinted scale. The gas composition was later analyzed by GC–MS (Shimadzu GCMS-QP2010 Ultra, Kyoto, Japan) with helium as carrier gas. Sampling, preparation, and analysis of scales The scale samples (20 samples in total) were taken from the inlet chamber of the heat exchanger on October 1, 2015. The facility was in full operation for 6 days before shut down and opening of the heat exchanger. A total of 46, 700 m of thermal water had passed through the heat exchanger after the last chemical maintenance. In order to assess the mineralogical composition and the structure of the scales by optical transmission microscopy, the samples were rinsed with n-hexane (Chemsolute, min. 99 %). Then, thin sections of an approximate thickness of 30μm were cut parallel to the flow direction and perpendicular to the ripple crests. The sections were embedded in a polyurethane resin (Araldite ). For better visualization of the pores, a blue dye (Stru- ers, Bluedye Farbstoff für Imprägniermittel DE-L831203) was added to the resin. The samples were then visually examined with an optical transmission microscope. To prepare the samples for scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM/EDX), the scale fragments were also rinsed with n-hexane (Chem- solute, min. 99 %) about five times to remove any oil. The samples were then rinsed with water once or twice and dried in a desiccator for a few days to remove residual moisture. After sputtering a gold layer onto the samples the samples were visualized with a SEM (Tescan Vega, 3 LM, Software: Vega TC, Kohoutovice, Czech Republic) and measured with EDX (EDX Oxford Instruments, x-act, Software: Aztex, Austin TX 78744, USA). The SEM/EDX images were processed with ImageJ (Schindelin et al. 2012; Rueden et al. 2017) in order to measure the pore- and crystal-size distribution. The porosity and pore size distribution of scale samples were measured using micro computer tomogra- phy (micro-CT) (v|tome|x s 240, General Electrics) by C. Schulbert for a typical sam- ple cube with 0.5 × 0.5 × 0.5 cm at the University of Erlangen. The scan was run at 80 kV, 100µA and a distance between object and focal point of 34 mm. Data reconstruc- tion and processing was done with the GE-software datos, Version 2.3. and the software VGSTUDIO MAX 3.0 (VolumeGraphics) (GmbH 2019). Hydrogeochemical modeling In this study we also wanted to test our hypothesis that the observed ripple scale did not form by simple precipitation during normal operation, but by deposition of preexisting particles, which we assume due to the ripple morphology. For this test, a hydrogeochem- ical model for normal operation was set up with PhreeqC, v.3.4 (Parkhurst and Appelo 2013) following the procedure of Baumann (2016). The model accounts for the differ - ence in ambient conditions between the sampling location and the geothermal reservoir, and uses the database phreeqc.dat. The first step to simulate the water and gas composition throughout the facility was a back-calculation of the samples taken at surface level to reservoir conditions. Next, the Köhl et al. Geotherm Energy (2020) 8:23 Page 10 of 27 Table 2 Analysis of the geothermal water from 24 June 2016 Field parameters pH (water at surface) 6.31 Electric conductivity, µS/cm T = 25 C 639 ref Redox potential, mV -247 One liter of water contains: Mass Equivalents mg mmol % Cations 129 5.62 76.5 Sodium (Na ) Potassium (K ) 20.5 0.524 7.13 3.10 0.172 2.34 Ammonium (NH ) 2+ 16.7 0.833 11.3 Calcium ( Ca ) 2+ 2.25 0.185 2.52 Magnesium (M g ) 2+ 0.54 0.012 0.17 Strontium (S r ) 2+ 0.17 0.003 0.03 Barium (Ba ) 2+ 0.00 0.000 0.00 Manganese (M n ) 3+,2+ 0.01 0.001 0.01 Iron(II)+(III) as Iron(II) (Fe ) Sum of cations 172 7.35 100.0 Anions Fluoride (F ) 2.81 0.148 2.03 Chloride ( Cl ) 88.2 2.49 34.1 Bromide (Br ) 0.45 0.006 0.08 2− 12.4 0.258 3.54 Sulfate (S O ) Hydrogen sulfide ( HS ) 1.91 0.058 0.79 Hydrogencarbonate (HCO ) 265 4.33 59.4 Sum of anions 371 7.29 100.0 Undissociated matter Silicic acid (H SiO ) 105 1.34 2 3 Boric acid (H BO ) 9.20 0.149 3 3 Total dissolved solids ( TDS) 656 pressure and temperature were adjusted according to the defined simulation points. For each simulation point the model calculated the bubble pressure which is the pressure at which gas bubbles begin to form in the geothermal water (gas phase > 0.1 mL/L ) and the theoretical mass of precipitates. Here, the mass was given as the amount of calcite, that can precipitate to reach thermodynamic equilibrium (SI  =  0). This simplification was justified by the mineralogical assessment of the scales which consist primarily of calcite. Quantification of ripple features The ripple geometries were measured by sliding caliper measurements. The measured parameters were ripple length, ripple height, and sample thickness. Results and discussion Hydrochemical conditions A representative analysis of the geothermal water is given in Table 2. Following the clas- sification of Furtak and Langguth (1965), the water is a Na– HCO –Cl-type with a TDS 3 Köhl  et al. Geotherm Energy (2020) 8:23 Page 11 of 27 Table 3 Analysis of  the  gas contained in  the  geothermal facility normalized to  standard ◦ ◦ conditions (1013 hPa, 0 C ). Degassing took place at 54 C and p = 1030 hPa Gas CH N CO C H H S 4 2 2 2 6 2 Gas, Vol % 46.3 35.9 16.0 1.0 0.3 Dissolved gas @ reservoir conditions (298 bar, 2.1 1.4 3.1 0.04 0.17 140.2 C ), mmol/L Gas load (sum of all gases), NmL gas/ L water 86 Table 4 Pressure, temperature, and saturation index throughout the geothermal facility Pressure, bar Temperature , C SI Calcite Reservoir during production 298.0 140.2 0.0 Pump intake 30.7 140.2 0.24 Pump outlet 84.4 140.2 0.19 In front of heat exchanger 12.0 139.0 0.24 Behind heat exchanger 10.5 60.0 1.00 of 656 mg/L. With a (Na+K-Cl)/HCO ratio of 0.84 the water is strongly influenced by cation exchange. The Ca/Mg ratio of 4.5 of the water indicates a limestone reservoir matrix (Langmuir 1971). The water contains 86  mL/L gas (at standard conditions, Table  3), which is com- posed of methane (46.3  %), nitrogen (35.9  %), carbon dioxide (16.0  %), ethane (1.0  %) and hydrogen sulfide (0.3  %). At reservoir conditions (298  bar, 140.2 C ) this equals 2.1  mmol/L methane, 1.4  mmol/L nitrogen, 3.1  mmol/L carbon dioxide, 0.04  mmol/L ethane and 0.17  mmol/L hydrogen sulfide. Carbon dioxide is in carbonate equilibrium and is stripped from the solution if other gases form gas bubbles. The calculated saturation index for calcite along the geothermal cycle is given in Table  4. With an age of the thermal water exceeding 10,000 years we can safely assume equilibrium conditions in the reservoir. At the pump intake, the SI increases to 0.24 due to decompression, slightly decreasing to 0.19 at the outlet of the pump due to pressure build-up in the pump, and rising back again to 0.24 at the well head, again due to decom- pression (Table  4). The pressure dependence of calcite equilibrium constants was vali - dated by Hörbrand et al. (2018). As soon as the temperature decreases inside the heat exchanger, the SI turns to under- saturation (SI =−1.00 at 60 C ). This leads to a benign behavior of the geothermal water with respect to the precipitation of calcite, and no calcite scales are expected between the heat exchanger and the well head of the injection pipe. Description of scales Figure  1 shows a representative scale from 20 evaluated samples. The measured curva - ture of the samples resulting from the cylindric shape of the pipe is d ≥ 290mm for all samples. This rules out any smaller pipes, and indicates that the scales must have formed in the ground-level facilities, with pipe diameters between 316 and 348  mm (Table  1). At this point, we can safely rule out the production pipes ( 255mm ) as a source for the remobilized scales. Köhl et al. Geotherm Energy (2020) 8:23 Page 12 of 27 Fig. 4 Thin section of scale embedded with blue dye and observed with analysator. Whole thin section (left) and detail (right) with layers A and A , B and B 1 2 1 2 The ripples crests are oriented along the inner perimeter of the pipes, that is perpen - dicular to the flow in the pipes. This rules out accumulation of material which slides down the inner pipe walls at no flow conditions. The median of the measured ripple length is 9.7  mm (25  %-quartile (Q1): 9.0  mm, 75  %-quartile (Q3): 10.9  mm; N  =  31). The median of the measured ripple height is 1.2 mm (Q1: 1.1 mm, Q3: 2.0 mm; N = 31). The ripple morphology shows gently sloping stoss sides and steep lee sides (Figs.  1 and 4). This indicates that the ripples were formed in unidirectional current (current rip - ples). Also, the Ripple Symmetry Index (RSI) of the ripples is 2.9 (Q1: 2.5, Q3: 3.5) which indicates current ripples. This makes sense, because the facility is operating either in a unidirectional flow regime or no flow at all. Oscillation, like in wave ripples, does not occur even though some of the ripples show a RSI < 2.5 which is at the transition to orbital waves. The median bulk density of the ripple scales calculated from the geometry and weight 3 3 3 of the ripples is 1.3 g/cm (Q1: 1.2 g/cm , Q3: 1.4 g/cm ). Since the ripples consist of cal- cite ( 2.7 g/cm ), this indicates a large porosity in the scale fragments. More information is revealed by the thin sections of the ripple scales (Fig. 4). The scale samples have a median total thickness of 5.0 mm (Q1: 4.5 mm, Q3: 5.0 mm), and are all composed of four layers. The initial layer ( A , 2 mm thickness) consists of a fine grained calcite matrix with many pores which appear round or even circular. During prepara- tion of the samples oil came out of these pores. A dense layer made of calcite with large crystals ( B , 400µm thickness) covers layer A . It shows only very few pores which are 1 1 again filled with oil. Layers A (2 mm thickness) and B ( 400µm thickness) repeat layers 2 2 A B and but are phase-shifted against the underlying layer. 1 1 The ripple structure is created in layers A and A , respectively. The calcite layers  B 1 2 cover the ripple structures of the layers A with a constant thickness and like that stabi- lize the ripple structures. Scanning electron micrographs (Figs. 5, 6) show large pores in the A layers surrounded by small calcite crystals, and small pores but large calcite crystals in the B layers. Köhl  et al. Geotherm Energy (2020) 8:23 Page 13 of 27 Fig. 5 SEM image of the pipe side (layer A ) showing large pores surrounded by small calcite Fig. 6 SEM image of the water side (B ) showing small pores but large calcite crystals The SEM/EDX analysis showed that the porous matrix of the A  layers and the crys - talline layers B have the same elemental composition that points to a magnesium-rich calcite (Table 5). A few small (< 20µm ) pyrrhotite (FeS) and pyrite crystals ( FeS ) were found in the samples, particularly in the A layers. XRD analysis indicated magnesium calcite as dominant crystalline phase (48 at. %) fol- lowed by calcite (35  at.  %). Interestingly, kutnohorite CaMn(CO ) was detected with 3 2 17 at. %. The absence of kutnohorite in SEM/EDX might be due to a limited field-of-view or due to magnesium substitution of the manganese ions. For the crystal and pore size distribution, 100 - 180 measurements for each layer of the scales were conducted on the thin sections using ImageJ (Schindelin et al. 2012; Rueden Köhl et al. Geotherm Energy (2020) 8:23 Page 14 of 27 Table 5 Average EDX analysis of the ripple scales in atomic % A A B B 1 2 1 2 C 21.2 21.5 21.2 21.7 O 63.7 64.5 63.7 65.2 Na 0.1 0.0 0.3 0.0 Mg 0.5 0.2 0.9 0.5 Al 0.0 0.0 0.1 0.0 Si 0.0 0.0 0.1 0.1 S 0.0 1.5 0.0 0.5 Ca 14.4 11.8 13.7 12.1 Fe 0.0 0.3 0.0 0.0 Cu 0.0 0.2 0.0 0.0 Sum 100.0 100.0 100.0 100.0 Ca/Mg 24 32 15 23 Measurements 6 5 8 7 et al. 2017). The crystal size was determined by averaging the longest two visible crystal edges, and the pore size was determined by averaging the longest axis and its perpen- dicular axis. The results revealed that the pore sizes and crystal sizes differ between the A and B layers. The crystals of the A layers are significantly smaller (median: 22µm , Q1: 18µm , Q3: 28µm ) compared to the crystals of the B layers (median: 80µm , Q1: 59µm , Q3: 144 µm ). Here, some crystals grew distinctly larger (up to 500µm ) resulting a strong skewness of the distribution. The anisotropy of the crystals was calculated based on the evaluation of the transmis - sion micrographs, and is expressed as ratio of the shortest crystal edge to the longest crystal edge. With a median ratio of 0.8 (Q1: 0.7, Q3: 0.9) in the A layers and a median ratio of 0.7 (Q1: 0.6, Q3: 0.9) in the B layers the crystals show a similar, more or less iso- metric, only slightly prismatic crystal habit indicating spatially unconfined growth. The pores in the A layers from 180  measurements (median: 92µm , Q1: 71µm , Q3: 117µm ) are significantly larger compared to the pores in the B  layers (median: 57µm , Q1: 46µm , Q3: 68µm ). The largest pores in the A layers ( 250µm ) are twice as large as the largest pores in the B  layers. The micro-CT data indicate a median pore size of 50 - 60µm in all layers. Evaluation of the transmission micrographs (e.g. Fig.  7), allows a quantitative assess- ment of the porosity. The average area covered by matrix elements is 55 % leaving a pore space of 45 %. The 3-D micro-CT measurements yield a porosity for the A layers of 42 - 53 %, which is in good agreement with the 2-D measurements with TM. The porosity given for the B  layers is 0–20  %. Here, no measurements were possible with the trans- mission micrographs. The anisotropy of the pores expressed as the ratio of short and long axis of the ellipsoid shows a slight median anisotropy of 0.8 (Q1: 0.7, Q3: 0.9) for all layers. All pores look more or less round. Pores like this are known from fluid inclusions. Some degree of ani - sotropy can be attributed to static pressure, droplets coalescing and later growth of crys- tals into the pore space. The anisotropy is confirmed by the micro-CT measurements. Köhl  et al. Geotherm Energy (2020) 8:23 Page 15 of 27 Fig. 7 TM image showing spherical oil droplets surrounded by a thin calcite layer Mass balance for accumulated scales In order to understand the scaling process, the growth rates of the scales are important. The mass of the ripple scales will provide some evidence about the time frame in which the scales built up. The mass of scale present in the ground-level pipe can be calculated for each scale layer from the observed scale thickness and the known geometries of the pipes according to Eqs. 13 and 14: m = ρ · (1 − �) · V , calcite (13) −3 where m is the mass of scale, ρ is the density of calcite ( 2.71 gcm ),  is the poros- calcite ity with an average value of 45  % taken from the TM and micro-CT analysis and V is volume of scale given by V = x · d · π · 1/3 · l, (14) where x is the observed median scale thickness (5.0 mm, Q1 = 4.5 mm, Q3 = 5.0 mm), d is the ground-level pipe diameter (350 mm), and l is the pipe length between automatic 50 ± 5m filter and heat exchanger ( ). This is the section where all the scale fragments, which were later washed into the heat exchanger, originally formed, as proven by the scaling curvatures (see Section: Description of the scales). The factor of 1/3 accounts for the observation that only the lower third of the pipe showed scales as reported from a camera rover inspection of the pipes done shortly after the sampling date. Applying Eqs.  13 and  14 with a grain density of 2.71 ± 0.05 g/cm and a porosity of 0.45 ± 5 in the A  layers, and neglecting the minor porosity in the compact B  lay- ers, the mass of mineral scale attached to the pipe between the filter and the heat exchanger at the sampling date was roughly 154 ± 22 kg with the highest uncertain- ties resulting from porosity and scaled area. Furthermore, in case of a stillstand of the facility, the geothermal water volume of 4.8 m sitting in the ground-level facility could only precipitate 0.3 % of layer B , and 1 Köhl et al. Geotherm Energy (2020) 8:23 Page 16 of 27 is so excluding the possibility, that layer B formed during a short stillstand of the geo- thermal production. If we assume full thermodynamic equilibrium, the minimum volume of water required to precipitate the observed scales can be calculated. The analyses of the scales shows that the main component is calcite. With a molarity of calcium β of 0.416 mmol/L, and a SI of 0.24 the calcite scale potential, derived from the PhreeqC model during normal operation in front of the heat exchanger at 139.0 C and 12 bar is 0.130  mmol/L. If that total scale potential precipitated, a minimum volume of 4 3 1.2 · 10 m of geothermal water would be required to create the observed mass of scale of 154 kg on the substrate surface of 18.3 m . With the flow rate of 110 L/s dur - ing full operation, the time until this water volume is produced is about 30 hours. The time until layer B is formed would be about 4 hours. In reality, however, scale formation is limited by the kinetics of calcite precipita- tion and also controlled by the effective mineral surface. The average calcite scale rate at the well head for this specific site is 2.7 · µmol/(m · s ) as reported by Her- brich (2017). This rate is also in line with other sites (Wanner et al. 2017). From this real scale rate and from the production rate, it can be calculated that it would take 361 days during normal operation to create the observed 154 kg of calcite scale in the pipes, and the cementation of layer B would take 46 days. The conclusion of the mass calculation is that the observed ripple scales cannot have formed by new precipitation during regular operation within the 6 days since the last cleaning. Solid particle deposition versus in situ precipitation In this study the particles which built the ripples originally precipitated by heterogene- ous crystallization of calcite on the walls of the horizontal ground-level pipes. Earlier rover inspections have proven this scale growth on the pipe walls. This scale formed over the time of one year, between the last ground-level cleaning, and the scale sample collec- tion in October 2015. Precipitation was caused by degassing of CO in the pump which increases the pH and leads to supersaturation with respect to calcite as described in ear- lier studies (Wanner et al. 2017, Herbrich 2017, Boch et al. 2017). Evidence for degassing during operation was given by a glass window at the well head, in which tiny gas bubbles could be observed in the thermal water. It is known from various natural environments that ripples are the result of solid par- ticle deposition (Ha and Chough 2003; Kraemer and Winter 2016; Schieber et al. 2007; Soulsby et al. 2012; Tanner 1967; Yalin 1985). Therefore, the ripple shape itself which is visible in the A  layers of the scale directly implies that solid particle deposition of pre- existing particles was the driving process for the formation of these scale layers. In spite of that, it cannot be excluded that some minor additional in situ precipitation occurred during the settling of the particles. Furthermore, individual grains can easily be man- ually detached from the A  layers of the scale samples. This strongly suggests that the A  layers consist of weakly connected particles which settled down in a slow current, and did not precipitate in  situ. Additionally, the oil inclusions trapped in the A  layers indicate a fast deposition process. As shown in the mass balance section it is not pos- sible that scale of that thickness (2 mm A layers) precipitated fast enough in situ during the 6  days of regular operation since the last cleaning. Secondly, it is also unlikely that Köhl  et al. Geotherm Energy (2020) 8:23 Page 17 of 27 scale of that thickness precipitated during the maintenance works (2 - 3 days) since this would require very high supersaturation with respect to calcite. Maintenance includes cleaning with organic acid which makes long lasting supersaturation unlikely. Consider- ing all arguments together, the hypothesis that the rippled A  layers formed from pre- existing particles is much more evident than in  situ crystallization. This process is not yet described in literature for geothermal facilities, and is not accessible with current models. In contrast to the A layers, the B layers of the ripple scales are only 400µm thick, have a dense structure, and cover the A  layers evenly with constant thickness which means without superimposing another ripple structure on top of the A layers. These relatively thin B layers can already be created within a relatively short period of transient supersat- uration. Such short supersaturation can be reached at the end of the maintenance works, during flushing of the pipes with communal water of higher pH. This suggests that the B layers were actually formed by in situ precipitation during the maintenance works, in opposition to the A layers. Ripple formation—calculations The key to the understanding of the formation of the observed scalings is the presence of ripple structures. The establishment of current ripples requires sedimented but mobile particles and fluid flow. Sedimentation requires that the density of the mobile particles is significantly higher than the density of water or the particles will float or stay in suspen - sion. By the use of empirical equations, we reconstructed the conditions at which the ripples formed. With Eqs.  4 and  5, we calculated the theoretical ripple geometries with the observed features and fluid properties. For the latter, temperature-dependent changes of the viscosity have most impact. From the comparison to the measured ripple geometry, we could then derive the flow velocities in the system at which the ripples could have formed. In our system, regression equations  4 and  5 are valid for calcite particles ( ̺ = 2710 kg/m ) with a diameter between approx. 20 and 220µm at 139 C (operation temperature) and between 50 and 630µm at 20 C (maintenance temperature). Particles with lower particle density cause a shift to higher particle sizes. For an initial assessment 4  different scenarios for the ripple formation were defined (Table 8): Scenario 1 is standard operation at full production rate based on the observed calcite particle diameter ( d = 22µm ), scenario 2 is like scenario 1 but assumes floccu - lation/agglomeration to larger calcite particles ( d = 150µm ). Scenario 3 is represent- ative for the cleaning step with water and acid at a lower pumping rate assuming solid particle deposition of the observed calcite particle size ( d = 22µm ). Values for D 50 ∗ are 0.6 ± 0.2 for scenario 3, therefore, ripple height and length could not be calculated for this scenario. Scenario 4 assumes solid particle deposition of calcite-coated oil drop- lets ( d = 92µm ), also during maintenance. This scenario was introduced because the SEM images showed a thin layer of calcite surrounding each pore (Fig.  7). The process which could have created these coated oil droplets is known from literature as Picker- ing emulsion (Huang et al. 2019), where oil droplets serve as crystallization nuclei, onto which calcite precipitates. Then, the coated oil droplets would have become more dense Köhl et al. Geotherm Energy (2020) 8:23 Page 18 of 27 and experienced deposition. Also for this scenario D is < 1.2 and thus ripple height and length could not be calculated. Before applying Eqs.  4 and  5, we made a quick test by the use of Stokes’ law if par- ticles could have settled down at a reasonable time in the defined scenarios 1  -  4, or if the particles would have stayed in suspension. This test confirmed that in scenario  1 the sedimentation of calcite crystals from the top to the bottom of the pipe, would have taken only 3 min, so that solid particle deposition is reasonable. For larger particles (sce- nario 2), the sedimentation would have happened even faster. During maintenance (sce- nario 3) due to higher viscosity, the same sedimentation would have taken 13 min. If the ripples had been formed by calcite-coated oil droplets (scenario  4) ( ρ = 900 kg/m , oil 1.5µm calcite coating) sedimentation of these particles would have taken 19  min. The calculations for scenarios 1  -  4 by Stokes’ equation showed that sedimentation of the particles was fast enough to occur in all considered scenarios. Then, we applied Eqs.  4 and  5 to test which scenario is most likely. The results of the calculated ripple height and length are presented in Table 8 (upper part). With the meas- ured crystal sizes and composition the maximum ripple length for current ripples is around 31  mm for the operational scenario  1. Particle aggregation increases the ripple length to 83  mm (scenario  2). The ripple height in both cases is roughly 10  % of the ripple length. During maintenance (scenarios  3 and 4), the maximum ripple lengths are expected to be even longer, due to increased fluid viscosity, but equations for rip - ple length and height are not valid. To summarize the ripple geometries, all calculated values based on Eqs.  4 and  5 were either exceeding the measured values (  ≈ 10 mm , η ≈ 1.2 mm ) by at least a factor of 3, or out of calibration range. This indicates that the ripples were likely not formed from the individual crystals observed in layers A, or that the equations cannot be used in this setting. Why did the calculations according to Eqs.  4 and  5 predict higher and longer ripples than we observed, or in other words which parameters would need to be changed to get the observed small ripples? • Smaller ripples could be reached by lower viscosities (see Eqs. 4 -  6). However, lower viscosities can be excluded, and so viscosity does not explain the discrepancy. Vis- cosities were likely even higher than assumed, due to oil contained in the water. • Furthermore, a higher density contrast leads to smaller ripple dimensions, too. How- ever, for production scenarios 1 and 2, the density contrast was already at the upper limit, because density of non-porous calcite was assumed for the initial particles, and so this does not explain the small ripples. For maintenance scenario 4, an increase of the particle density to 1200 kg/m would allow using Eqs.  4 and  5, but the resulting ripple height and length would be 19  mm, respectively, 188  mm and therefore also way too long. The small ripples are therefore not explained by an erroneous density contrast, either. • However, the fine A  layers of the scale consist of mixture of cohesive and non- cohesive sediments. Such mixtures are reported to have significantly smaller ripple dimensions (height and length) due to cohesive forces (Baas et al. 2019, 2016). This could explain why our calculations based on experiments with pure non-cohesive material produced too large dimensions. For mixtures of cohesive and non-cohesive Köhl  et al. Geotherm Energy (2020) 8:23 Page 19 of 27 sediments somewhat smaller velocities for ripples formation would be expected (Baas et al. 2019; Schieber et al. 2007; Yawar and Schieber 2017). In contrast to the sophisticated Eqs. 4 – 6, the simple Eqs. 1 – 3 (see Table 8, lower part) which only depend on d showed relatively good agreement with the observed geom- etries even though they neglect other parameters. This suggests that these equations are still a good first guess if ripples are found in new settings. As a next step for the reconstruction of the ripple-forming conditions, the velocity was calculated according to Eqs.  6  -  12. This gave another hint that the ripples had not formed from individual crystals observed in layer  A: during full operation (scenario  1 and  2), the flow velocity in the pipes was well above the washout velocity of 0.52  m/s, and during maintenance (scenario 3) the critical velocity (0.21 m/s) for ripple formation from calcite particles was not reached (Table  8, upper part). For scenario  4, the calcu- lated flow velocities (0.04 – 0.23 m/s) matched the real flow conditions during mainte - nance, but since the equations for ripple length and height, which belong to the same regression of the data set, were not valid, also the velocity equations could strictly spo- ken not be applied for this scenario. Figures  8 and  9 show the calculated ripple forming velocities as a function of d for scenarios  1  –  3. The figure demonstrates the strong effect of temperature and viscos - ity on the shape of the curve U . While for 20 C , the washout velocity U washout washout increases with increasing particle size (Fig. 9), the washout velocity decreases with rising particle size, for a temperature of 139 C (Fig.  8). However, for the considered particle size ( 22µm ) of the ripple scale, the effect is small. Moreover, Figs.  8 and 9 show, which formation velocities of the ripples would be calculated for larger particles, as for instance created by flocculation. Since there is already a huge uncertainty for the viscosity in our system, the same applies for the calculated velocities. While oil on the surface of particles increases the viscosity, high temperatures lead to a drastic decrease in the viscosity. If the ripple scales had been formed during maintenance, low temperatures of about 20 C would have been present. Then, viscosity would have been 4 times higher than during normal operation ◦ ◦ at 139 C . Calculations for 20 C seemed in general better applicable than calculations at 139 C because the experiments with which the equations were set up were done at ambient temperatures. Discussion of the ripple‑forming scenarios Considering all calculations and aspects, when did the ripple scales, consisting of the rippled A  layers and the cemented B  layers, form? Table  6 summarizes the different options that are discussed in the following: • Ripple formation (A layers) could have happened during full production (scenario 1 and 2), even if it is not predicted by the calculations, but with the viscosity being out- side the calibration of Soulsby et al. (2012) ripples could possibly have formed even at the high flow velocity of 1–2 m/s. The observed small ripple dimensions can well be explained by cohesion of the silt size calcite crystals (Baas et al. 2019, 2016). This scenario seems unlikely, since no one reported ripples at these flow velocities before. Köhl et al. Geotherm Energy (2020) 8:23 Page 20 of 27 Table 6 Schematic overview of  different scale and  ripple formation mechanisms in  the  ground-level pipes between  filter and  heat exchanger with  indication of  the  calculated scenarios 1–4. No calculation scenarios were defined for  transition phases (start of  pump, shut down  of  pump) since  these could be ruled out  for  several reasons Flow rate, L/s Flow Pressure, bar Process velocity, m/s Start of production 110 1.40 12 Remobilization of existing scale particles due to turbulence. Pump rate is increased too fast to allow for ripple formation and stabilization Full production (scenarios 1 110 1.40 12 Slow scale precipitation due to and 2) degassing and oversaturation. Velocity > calculated ripple washout velocity, therefore no ripples predicted Shut-down of production 0.1 - 0.5 0.00 1 - 4 Solid particle deposition (rip- ples) of suspended particles followed by stabilization of scales due to sudden pressure drop and degassing. Too little particle supply Acid treatment (scenarios 3 5 0.06 1.8 - 2.6 Acid causes remobilization and 4) of partially dissolved scale particles and subsequent solid particle deposition (ripples). Oil prevents total dissolution of scales. Stabilization of ripples takes place when pH rises too high Flushing of facilities with tap 6.5 0.08 1.6 - 3.5 Stabilization of ripples by calcite water after acid treatment precipitation due to high cal- (scenarios 3 and 4) cium ion concentration in the cleaning solution and a neutral pH of the tap water Table 7 Course of  events during  the  3 months before  scale sample collection. Evaporator cleaning is  either  done by  flushing the  evaporator, or  by  opening it with  subsequent mechanical removal, or  by  chemical cleaning. T = temperature in  the  ground-level pre pipes before  interruption, T = temperature in  the  ground-level pipes at  the  end post of the interruption ◦ ◦ Date Cause Duration Flow rate, L/s T , C T , C Remark pre post Aug 3 Electric 11 d 0 140 25 Aug 17 Scales 3 h 110 139 126 Flushing of evaporator Aug 20 Scales 3 h 0 140 116 Mechanical removal Aug 24 Synchronization 2 h 0 141 120 Aug 28 Scales 2 h 110 139 132 Flushing of evaporator Sep 3 Scales 4 h + 2 h 0 140 118 Mechanical removal, 2 pump starts Sep 12 Revision 13 d 5 140 20 Chemical cleaning Oct 1 Scales 3 h 0 140 20 Sample collection If, despite of that, the A layers of the ripples formed during production, the B layers would have formed during different conditions (e.g. flushing with tap water during maintenance) due to their different compact appearance. Köhl  et al. Geotherm Energy (2020) 8:23 Page 21 of 27 Table 8 Comparison of  real observations and  conditions with  calculation results of  scenarios  1–4. The critical flow velocity U and  the  washout flow critical=min. velocity U were calculated by  Eq.  12 considering also  the  parameter ranges. The ripple height and  ripple length were calculated by  Eqs.  4 and  5 washout=max . considering also the parameter ranges. In the lower part of the table, ripple length and height predictions by the simple Eqs. 1–3 are indicated for comparison. Empty fields: unknown parameters, n.v. = equations not valid for the given input parameters (out of calibration) Observed ripple Observed conditions Observed conditions Scenario 1 (production Scenario 2 (production Scenario 3 (maintenance Scenario 4 geometries during production during maintenance 22μm calcite) 150μm flocc. calcite) 22μm calcite) (maintenance 92μm coated oil) T, C 139 20 139 139 20 20 ˙ 110 5 ± 2 110 110 5 ± 2 5 ± 2 −7 2 2.19 9.81 2.19 2.19 9.81 9.81 ν , 10 · m /s 2710 ± 0.14 2710 ± 0.14 2710 ± 0.14 2710 ± 0.14 2710 ± 0.14 1073 ± 28 ρ , kg/m 927 927 927 1020 ± 10 1020 ± 10 ρ , kg/m −6 22 (Q1: 18, Q3: 28) 22 (Q1: 18, Q3: 28) 150 22 (Q1: 18, Q3: 28) 92 (Q1: 71, Q3: 117) d , m · 10 D* 1.35 – 2.00 11.00 0.44 – 0.73 0.38 - 1.4 0.09 - 12 0.03 0.16 - 0.20 0.12 - 0.21 U , m/s 1.2 0 0.18 0.22 0.21 0.04 critical=min. U , m/s 2.2 0.08 0.52 0.33 0.58 0.23 washout=max . −3 1.2 3.1 - 3.9 8 n.v. n.v. η , m · 10 max −3 9.7 31 - 33 83 n.v. n.v. , m · 10 max −3 1.2 1.2 - 1.9 10 1.2 - 1.9 4.8 - 7.9 η , m · 10 max−flem. −3 9.7 8.6 - 12.0 42 8.6 - 12.0 24 - 35 , m · 10 max−colem. −3 9.7 18 - 28 150 18 - 28 71 - 117 , m · 10 max−yalin Köhl et al. Geotherm Energy (2020) 8:23 Page 22 of 27 Fig. 8 Boundaries for ripple-forming flow velocities U based on calcite particle size d (scenario 1 between ◦ −7 2 dashed lines). T = 139 C , ν = 2.19 · 10 m /s Fig. 9 Boundaries for ripple-forming flow velocities U based on calcite particle size d (scenario 3 between −7 2 dashed lines). T = 20 C , ν = 9.81 · 10 m /s Köhl  et al. Geotherm Energy (2020) 8:23 Page 23 of 27 • The second option is that ripples (A layers) formed during maintenance (scenario 3 and  4). The calculated velocities are higher than the maintenance velocities. At maintenance velocity, calcite particles would be below the critical threshold. Cal- cite-coated oil droplets fit to the flow velocities, but are unlikely regarding the rip - ple dimensions. However, ripple formation during maintenance is supported by the evidence of the 4  layer structure of the scale which fits well to the cleaning cycle consisting of 2 pumping phases, where the A layers could have formed, each one fol- lowed by a stillstand overnight, in which the B layers could form on top and stabilize the ripples at increasing pH due to degassing of CO . Alternatively, the B layers could have formed during flushing of the pipes at the end of the maintenance. Here, strong supersaturation was reached by high calcium ion concentrations in the fluid leading to calcite precipitation when the fluid was mixed to tap water with neutral pH. While this scenario assumed ripple formation (A  layers) during maintenance, the formation of the B layers could also have happened after the maintenance works, i.e. during normal production. However, this seems unlikely, since B  layer cementation would take 46  days during normal production. This fits neither to the 6  days since the last chemical cleaning, nor to the time span of one year lying between the yearly chemical cleanings. Furthermore, a cementation of the B layer during full production is unlikely because the B layers contain almost no oil inclusions which suggests grav- ity driven oil–water phase separation at quiescent conditions, not at high turbulence. • Another option for the formation of the ripple scales is the start phase of the pump. During the start phase, according to our calculations of the flow velocity, ripple formation was possible. Particles which had not fully been removed by earlier acid treatments of the facility, could have been remobilized and have formed ripples. This requires that the pump rate was only very slowly increased, so that the particles could stabilize by new precipitation of calcite before higher flow velocities could have lead to washout. Since the pump, at the investigated time period, was usually imme- diately turned on to almost final pump rate, this scenario seems unlikely. • Also the shut-down phase is a potential candidate for ripple scale formation. During the shut-down phase of the pump, suspended particles could have settle down in a ceasing current creating a rippled A layer. However, this scenario seems unlikely, as well, because there are not enough particles in suspension during normal produc- tion to create a continuous layer of the observed thickness (2 mm scaling layer A lay- ers in Fig.  4). The particle load in geothermal facilities in the Molasse Basin is only about 0.01  –  1  mg/L (Wolfgramm et  al. 2011), so that additional preexisting parti- cles (scales) would have needed to be remobilized to reach that thickness. Secondly, suspended particles during full operation consist of other particles, too (Wolfgramm et  al. 2011), and do not consist solely of calcite as observed in our study. Third, rip - ples need some time to be created (Soulsby et  al. 2012). During laboratory experi- ments, ripple formation only started after 40 min, while before that, rolling particles did not pile up significantly (Doucette and O’Donoghue 2006). Therefore ripple for - mation during shut-down phase is unlikely. To summarize, the calculated flow velocities contradict that the ripples were formed during full production and during maintenance. Additionally, the start- and shut-down Köhl et al. Geotherm Energy (2020) 8:23 Page 24 of 27 phase have the calculated flow velocities, but seem also unlikely for ripple formation for several above explained reasons. This suggests that the investigated ripple scalings in geothermal facilities did actually form at conditions which are not predicted by the empirical equations. This is possible since the ripples of this study were formed in a complex multiphase system (water, gas, oil) with simultaneous precipitation, dissolu- tion and degassing while the empirical equations are derived from controlled laboratory experiments. Major differences between the systems are that the empirical equations were derived from experiments with particles in the sand grain size while in this study particles are in the silt size. However, studies have reported ripples formed from silt and clay particles, as well. These fine grained ripples are formed by spontaneous floccula - tion of the particles which by that get large enough to settle down according to Stokes’ law (Schieber and Southard 2009). Furthermore, the particles in this study have a rhombohedral shape, while empirical experiments mostly studied ripples composed of rounded particles. The process of ripple formation comprising entrainment of particles, transient suspension, and reattachment as described by Ha and Chough (2003) is also possible with poorly rounded particles. Well rounded rolling particles are not required for ripple formation. Furthermore, rip- ples can also build from clay- and silt-size particles (Baas et al. 2019, 2016; Schieber and Southard 2009) which are typically also not rounded but rather platy. This shows that the principles of ripple formation can be applied to the rhombohedral particles of this study, as well. However, the different particle shape may explain some derivation from the quantita - tive theoretical prediction. The exposed surface area and the drag resistance of rhom - bohedral particles in a fluid current are different to rounded particles. This interplay of drag force and drag resistance due to the particle shape may influence the produced geometries and necessary ripple-forming velocities. However, this is yet little described in literature, and should therefore be addressed in future research. In literature (Soulsby et al. 2012; Baas et al. 2016) mostly only the particle diameter D is considered as criti- cal geometric parameter, but the shape of the particle is neglected in the equations. Conclusion Given all differences and uncertainties, we come to the conclusion that the ripples inves - tigated in this study formed during maintenance works at a very low flow rate (scenario 3 or 4). The silt size of the particles supports the small flow rates. Based on the evaluation, we conclude that the rippled scalings were created the following way. The facility was shut-down for the yearly maintenance. After displacement of the ther - mal water by tap water, acid was added to the solution which partially dissolved and destabilized the scaling, but did not reach all parts of the scale. During the acid- circula- tion, small scale fragments or individual crystals or rounded grains were transported, and underwent solid particle deposition creating ripples (layer A ). Later, layer B pre- 1 1 cipitated on top of layer A . There are two options, when and how layer  B formed: 1 1 • The first option is that layer  B formed overnight, at quiescent conditions, when the acid worked leading to degassing of CO and a rising pH exceeding the thermody- namical equilibrium, so that dissolution was followed by precipitation/cementation Köhl  et al. Geotherm Energy (2020) 8:23 Page 25 of 27 from the fluid with a very high calcium ion concentration. Like that, layer  B formed on top of layer A and stabilized it. The procedure was repeated the next pumping day and the following night resulting in layer A and B . It is striking that there are 2 2 2 pump phases and 2 stillstand phases during maintenance while there are also 2 rip- ple layers and 2 calcite cement layers. • The second option, how layer  B formed, is that it formed during flushing of the pipes at the end of the cleaning procedure. The fluid in the pipes has a high calcium ion concentration due to the dissolution of calcite. When this fluid gets mixed with fresh tap water of neutral pH, strong temporary supersaturation with respect to calcite is likely. If that is, what happened, then the observed ripple scales were created during 2 maintenance events, each event forming one A layer and one B layer, since flushing is only done once at the end of each cleaning cycle. Then, layer  A and B would have 1 1 formed before and during the cleaning one year ago, and subsequently would have withstood the last cleaning cycle. Rippled scales give insight into the hydrochemical conditions and fluid dynamics in the investigated geothermal facility, and might be a key to avoid cost-intensive shut-downs due to scales in the facility. According to our calculations based on Soulsby et al. (2012), ripple formation is nei- ther predicted for the production phase nor for the maintenance phase because the cal- culated velocities lie in between. However, major differences between this study and the experimental setup, from which the equations used for the calculations, were derived, suggest that ripples were created during maintenance, in spite of that. Particle diameter and fluid viscosity are major parameters that differ between this study and the empirical predictions. Therefore, we conclude that solid particle deposition of mobile calcite crystals at the bottom of the pipe created the rippled scalings during maintenance works. We are con- fident that due to a protective oil layer on an already existing scale, the acid used during the maintenance was not able to entirely dissolve the scale. However, the acid weakened and mobilized small particles from the scale which were redeposited as ripples in a weak current. Here, either pH increase overnight or a small oversaturation of the tap water at the end of the cleaning was sufficient to stabilize the ripples by precipitating a tiny layer of calcite on top. In general, oil plays an important role in the formation of calcite scale. It can mitigate the initial phase of scale formation by making the attachment to the pipe harder, but it can also somewhat protect an existing scale from acid. Furthermore, operational shut- downs tend to create oil layers by promoting water–oil phase separation. Finally, the importance of removing the oil from the scale before and during acid treatment must be emphasized. Only when you know that all scale can be accessed by the acid and a steady pH is reached with no further acid addition the facility can be considered clean. Abbreviations SEM/EDX: Scanning electron microscope with energy-dispersive X-ray spectroscopy; SI: Saturation Index. Acknowledgements The financial support by the Bayerisches Staatsministerium für Wissenschaft und Kunst is gratefully acknowledged. We thank the Stadtwerke München for providing us the access to their geothermal facility, to their scanning electron Köhl et al. Geotherm Energy (2020) 8:23 Page 26 of 27 microscope and laboratory and for supplying us with the operational data. We would very much like to thank Prof. Dr. Martin Elsner and Dr. Martina Ueckert for multiple readings and discussions of the manuscript. Furthermore we thank Christian Schulbert from the University of Erlangen for the micro-CT measurements and visualization. Vladimir Ruttner from the Chair of Engineering Geology at TUM is gratefully acknowledged for the preparation of the thin sections. We also want to thank Amanada Günther for the XRD-measurements and Christine Benning for additional SEM/EDX measurements. Authors’ contributions BK collected the water, gas and ripple scale samples, analyzed the scales, performed the PhreeqC calculations, checked the ripple reconstruction calculations of JG and discussed the formation scenarios with JG and TB. JG performed the rip- ple reconstruction calculations, discussed the scenarios and assisted during the water and gas sampling. TB had the idea, developed the hydrogeochemical modeling framework and discussed the calculations and scenarios with BK and JG. All authors read and approved the final manuscript. Funding This work was partially funded by the Bavarian Ministry of Science and Art in the project Geothermie-Allianz Bayern. Open access funding provided by Projekt DEAL. Availability of data and material The datasets used and analyzed during the current study are available from the corresponding author on reasonable request. Compliance with ethical standards Competing interests The authors declare that they have no competing interests. Author details 1 2 Institute of Hydrochemistry, Technische Universität München, Marchioninistrasse 17, 81377 Munich, Germany. Depart- ment of Civil, Architectural, and Environmental Engineering, The University of Texas at Austin, 301 E Dean Keeton St., 78712 Austin, TX, USA. Chair of Hydrogeology, Technische Universität München, Arcisstr. 21, Munich 80333, Germany. Received: 8 March 2020 Accepted: 4 August 2020 References Agemar T, Alten J-A, Ganz B, Kuder J, Kühne K, Schumacher S, Schulz R. The Geothermal Information System for Germany - GeotIS. Zeitschrift der Deutschen Gesellschaft für Geowissenschaften. 2014;165(2):129–44. https ://doi. org/10.1127/1860-1804/2014/0060. Allaby M, editor. A Dictionary of Earth Sciences. 3rd ed. : Oxford University Press, online version; 2008. www.oxfor drefe rence .com/view/10.1093/oi/autho rity.20110 80310 04223 07. Avila K, Moxey D, de Lozar A, Avila M, Barkley D, Hof B. The onset of turbulence in pipe flow. Science. 2011;333(6039):192– 6. https ://doi.org/10.1126/scien ce.12032 23. Baas JH. A flume study on the development and equilibrium morphology of current ripples in very fine sand. Sedimen- tology. 1994;41:185–209. https ://doi.org/10.1111/j.1365-3091.1994.tb014 00.x. Baas JH, Best JL, Peakall J. Predicting bedforms and primary current stratification in cohesive mixtures of mud and sand. J Geol Soc. 2016;173(1):12–45. https ://doi.org/10.1144/jgs20 15-024. Baas JH, Baker M, Malarkey J, Bass S, Manning A, Hope J, Peakall J, Lichtman ID, Ye L, Davies A, Parsons D, Paterson D, Thorne P. Integrating field and laboratory approaches for ripple development in mixed sand-clay-eps. Sedimentol- ogy. 2019;66:2749–68. https ://doi.org/10.1111/sed.12611 . Bauer M, Freeden W, Jacobi H, Neu T. Handuch Tiefe Geothermie. Berlin: Springer; 2014. Baumann T. Validation of hydrochemical analyses and gas concentrations of deep geothermal aquifers. In: 41st Workshop on Geothermal Reservoir Engineering. Stanford: Stanford University; 2016. p. 1269–77. https ://pange a.stanf ord.edu/ ERE/pdf/IGAst andar d/SGW/2016/Bauma nn.pdf. Berg JH, Gelder A. A new bedform stability diagram, with emphasis on the transition of ripples to plane bed in flows over fine sand and silt. God. 2009;. https ://doi.org/10.1002/97814 44303 995.ch2. Birner J, Mayr C, Thomas L, Schneider M, Baumann T, Winkler A. Hydrochemistry and evolution of deep groundwaters in the Malm aquifer in the bavarian part of the South German Molasse Basin. Z Geol Wiss. 2011;39(3/4):291–308. Boch R, Leis A, Haslinger E, Goldbrunner J, Mittermayr F, Froeschl H, Hippler D, Dietzel M. Scale-fragment formation impairing geothermal energy production: interacting H2S corrosion and CaCO3 crystal growth. Geotherm Energy. 2017;5:1–19. https ://doi.org/10.1186/s4051 7-017-0062-3. Coleman SE, Fedele JJ, Garcia MH. Closed-conduit bed-form initiation and development. J Hydraul Eng. 2003;129(12):956–65. https ://doi.org/10.1061/(ASCE)0733-9429(2003)129:12(956). Doucette JS, O’Donoghue T. Response of sand ripples to change in oscillatory flow. Sedimentology. 2006;53(3):581–96. https ://doi.org/10.1111/j.1365-3091.2006.00774 .x. Flemming B. On the classification of subaquatic flow-transverse bedforms. 1988;29. Freimann R, editor. Hydraulik Für Bauingenieure. Grundlagen und Anwendungen. München: Carl-Hanser; 2014. Furtak H, Langguth H-R. Zur hydrochemischen Kennzeichnung von Grundwässern und Grundwassertypen mittels Ken- nzahlen. IAH-Congress. 1965;7:89–96. GmbH VG. VGSTUDIO MAX 3.0. 2019. https ://www.volum egrap hics.com/en/produ cts/broch ure-downl oad.html. Köhl  et al. Geotherm Energy (2020) 8:23 Page 27 of 27 Ha HK, Chough SK. Intermittent turbulent events over sandy current ripples: a motion-picture analysis of flume experi- ments. Sediment Geol. 2003;161:295–308. https ://doi.org/10.1016/S0037 -0738(03)00132 -5. Herbrich M. Anlagen- und sicherheitsrelevante Prozesse im Thermalwasserkreislauf eines geothermischen Heizkraft- werks. PhD thesis, Technische Universität München. http://nbn-resol ving.de/urn/resol ver.pl?urn:nbn:de:bvb:91- diss-20161 209-13102 43-0-3 Hörbrand T, Baumann T, Moog HC. Validation of hydrogeochemical databases for problems in deep geothermal energy. Geotherm Energy. 2018;. https ://doi.org/10.1186/s4051 7-018-0106-3. Huang F, Liang Y, He Y. On the pickering emulsions stabilized by calcium carbonate particles with various morphologies. Coll Surf. 2019;580:123722. https ://doi.org/10.1016/j.colsu rfa.2019.12372 2. Kraemer K, Winter C. Predicted ripple dimensions in relation to the precision of in situ measurements. Ocean Sci Discuss. 2016;. https ://doi.org/10.5194/os-2016-20. Langmuir D. The geochemistry of some carbonate ground waters in central Pennsylvania. Geochimica et Cosmochimica Acta. 1971;35:1023–45. https ://doi.org/10.1016/0016-7037(71)90019 -6. Mayrhofer C. Hydrochemische Untersuchungen im Malmaquifer im bayerischen Molassebecken. Dissertation, Technis- che Universität München, München. 2013. https://mediatum.ub.tum.de?id=1128419 Parkhurst DL, Appelo CAJ. Description of input and examples for PHREEQC version 3 – a computer program for specia- tion, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Technical Report Techniques and Methods 6–A43, U.S. Geological Survey. 2013. http://pubs.usgs.gov/tm/06/a43/. Rotta J. Experimenteller Beitrag zur Entstehung turbulenter Strömung im Rohr. Ingenieur-Archiv. 1956;24:258–81. Rueden CT, Schindelin J. Imagej for the next generation of scientific image data. BMC Bioinform. 2017;. https ://doi. org/10.1186/s1285 9-017-1934-z. Schade H KE. Strömungslehre. Berlin New York: de Gruyter; 1989. Schieber J, Southard JB. Bedload transport of mud by floccule ripples-direct observation of ripple migration processes and their implications. Geology. 2009;37:483–6. https ://doi.org/10.1130/G2531 9A.1. Schieber J, Southard J, Thaisen K. Accretion of mudstone beds from migrating floccule ripples. Science. 2007;318(5857):1760–3. https ://doi.org/10.1126/scien ce.11470 01. Schieber J, Southard JB, Kissling P, Rossman B, Ginsburg R. Experimental deposition of carbonate mud from moving suspensions: Importance of flocculation and implications for modern and ancient carbonate mud deposition. J Sediment Res. 2013;83(11):1025–31. https ://doi.org/10.2110/jsr.2013.77. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Tinevez J-Y, White DJ, Hartenstein V, Eliceiri K, Tomancak P, Cardona A. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;. https ://doi.org/10.1038/nmeth .2019. Soulsby RL, Whitehouse RJS, Marten KV. Prediction of time-evolving sand ripples in shelf seas. Contin Shelf Res. 2012;38:47–62. https ://doi.org/10.1016/j.csr.2012.02.016. Southard JB, Boguchwal LA. Bed configuration in steady unidirectional water flows; Part 2, Synthesis of flume data. J Sedi- ment Res. 1990;60(5):658–79. https ://doi.org/10.1306/212F9 241-2B24-11D7-86480 00102 C1865 D. Tanner WF. Ripple mark indices and their uses. Sedimentology. 1967;9:89–104. https ://doi.org/10.1111/j.1365-3091.1967. tb013 32.x. Wanner C, Eichinger F, Jahrfeld T, Diamond LW. Causes of abundant calcite scaling in geothermal wells in the Bavarian Molasse Basin, Southern Germany. Geothermics. 2017;70:324–38. https ://doi.org/10.1016/j.geoth ermic s.2017.05.001. Wolfgramm M, Rauppach K, Thorwart K. Mineralneubildungen und Partikeltransport im Thermalwasserkreislauf geother- mischer Anlagen Deutschlands. Z Geol Wiss. 2011;39:213–39. Yalin MS. On the determination of ripple geometry. J Hydraul Eng. 1985;111(8):1148–55. https ://doi.org/10.1061/ (ASCE)0733-9429(1985)111:8(1148). Yawar Z, Schieber J. On the origin of silt laminae in laminated shales. Sediment Geol. 2017;360:22–34. https ://doi. org/10.1016/j.sedge o.2017.09.001. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Journal

Geothermal EnergySpringer Journals

Published: Aug 20, 2020

There are no references for this article.