Hydrochemical and operational parameters driving carbonate scale kinetics at geothermal facilities in the Bavarian Molasse Basin

Abstract

The majority of scales observed at geothermal facilities exploring the Malm Aquifer in the Bavarian Molasse Basin are carbonates. They form due to a disruption of the lime–carbonic acid equilibrium during production caused by a reduction of the partial pressure of carbon dioxide due to pressure change and degassing. These scales are found at the pumps, production pipes, filters, heat exchangers, and occasionally in the injection pipes. In this study, scales of all sections of geothermal facilities were taken. The database consists of scale samples from 13 geothermal pumps, 6,000 m production pipe (sample interval 10 - 12 m), 11 heat exchanger revisions, 2 injection pipes, and numerous filter elements. The samples were analyzed by SEM-EDX, XRD, Raman spectroscopy, and acid digestion to assess their chemical and mineralogical composition. From direct gauge measurements at six facilities during pump changes, scale rates were determined along the production pipes. From indirect measurements (multifinger caliper measurements) scale rates are derived for the region below the pump. Hydrochemical analyses from the wellhead were taken from 13 sites to feed the hydrogeochemical models. The calcite scale rates in the production pipes increase from the pump to the wellhead, where they reach 1.5 - 4.1 μmol/(m2· s). Scale rates below the pump reach up to 1.5 μmol/(m2· s). Given the slight change of hydrochemistry on the rise through the production pipe, where < 4 % of dissolved calcium ions precipitate as scale, scale rates cannot be derived from water samples at the wellhead, but require direct gauge measurements. The small amount of precipitation, together with fully turbulent conditions suggests that all measured rates are controlled by the surface-reaction of calcite crystallization following the nomenclature of Appelo and Postma (2004). Two approaches are used for the modeling of the scale rates. The first approach is based on hydrogeochemical modeling with PHREEQC. Scale rates calculated by this method are one order of magnitude higher than the measured ones. The second approach is based on correlations between the measured scale rates at the wellhead at six facilities and identified thermodynamic scale drivers (Δ log (pCO 2), Δ total pressure, Δ pH, and SI calcite). The correlations allow linear regressions which are used for the prediction of the scale rate at the wellhead, along the whole production pipe, and below. The modeling results show that scale prediction based on the new regressions that rely on thermodynamic scale drivers works better than existing hydrogeochemical models, already without implementation of kinetic parameters (CO 2-stripping and magnesium inhibition).