Physicochemical Composition Variability and Hydraulic Conditions in a Geothermal Borehole—The Latest Study in Podhale Basin, Poland

In deep geothermal boreholes, an effect of temperature (so-called thermal lift) is observed, which results in the volumetric expansion of the fluid extracted. This process results in increased wellhead pressure values being recorded; in the absence of an appropriate correction, hydraulic properties of the reservoir layer cannot be properly determined. As an example of this effect, the Chochołów PIG-1 (CH PIG-1) geothermal borehole situated in Podhale Basin in Poland was used. Hydrodynamic tests including two pumping phases were carried out in the well in order to establish the basic hydraulic properties related to the determination of its operational resources (maximum water extraction value–capacity) and permissible groundwater level. Particular attention was paid to the thermal lift effect in the borehole. The conductivity, which depends on the pumping level, could be two to three times higher with temperature correction than results without any correction. The goal was to analyse the variability of the observed physiochemical properties of the exploited geothermal waters and to determine the correlation between the properties analysed and the temperature of the geothermal water. For the relationship between temperature and the observed pressure at the wellhead, the value of the correlation coefficient was negative (a negative linear relationship was determined), which means that as the temperature increases, the wellhead pressure decreases. The hydrodynamic tests carried out in the CH PIG-1 borehole and the analysis of variability of selected ions and parameters in exploited water were necessary to assess the possibility of increasing the efficiency (Q) of the CH PIG-1 borehole and to determine the water quality and its natural variability. Such information is crucial for the functioning of the recreational complex based on the use of geothermal water. A study of the phenomena affecting the exploitation of hot water from deep boreholes enables their effective exploitation and the use of resources in accordance with the expectations of investors.


Introduction
The physical process of volumetric expansion of fluid is observed in a borehole's column used to extract thermal waters from the depths, and this is known as thermal lift. Details of the process are described below. During the exploitation of geothermal boreholes, this effect may result in elevated wellhead pressure values being recorded. Only thorough knowledge of this phenomenon and appropriate corrections make it possible to properly determine the hydraulic properties of the reservoir layer. The impact of this phenomenon on the correct interpretation of hydrodynamic test results increases with borehole depth and formation temperature [1,2]. In shallow hydrogeological boreholes exploiting ordinary (not thermal) water, the pressure head in the borehole is practically independent of the temperature of the aquifer. The observed fluctuations of the water table in such boreholes are mainly a function of borehole efficiency and the filtration parameters of the aquifer. In deep geothermal boreholes, the effect of temperature on the observed water table level or measured wellhead pressure can be significant and increases with the deposit temperature depth of the borehole. The borehole is heated during exploitation. The temperature inside the casing pipes approaches the temperature at the bottom of the borehole with the duration of the exploitation, thus contributing to lower water heat loss during intake. During the heating of the borehole, a volumetric expansion of the water due to the reduction of its density can be observed. During the exploitation of the thermal borehole, this effect is manifested by an increased water table level or increased wellhead pressure compared to the cold borehole, although the bottom pressure is the same in both cases. This affects the design filtration parameters, mainly hydraulic conductivity. Since it is often impossible to measure the bottom pressure during a hydrodynamic test, the results obtained from measuring the water table or the wellhead pressure should be adjusted for the thermal effect by calculating the so-called reduced pressure [2,3], thus allowing the determination of absolute drawdown in the borehole, undisturbed by thermal conditions of the medium. By adjusting the results of hydrodynamic tests for the effect of borehole heating, more accurate values of filtration parameters of the aquifer, including hydraulic conductivity, are obtained [1,2]. The studies carried out so far in a borehole with a depth of approx. 2 km indicate that the conductivity calculated for reduced pressure (taking into account the thermal effect) is more than two times, and in the case of depression in the first phase of pumping even three times lower than the one calculated for the observed pressure [1]. Not taking thermal lift into account when interpreting hydrodynamic test results invariably leads to aquifer conductivity being overstated. The correct determination of aquifer filtration properties only becomes possible after the thermal lift effect has been adjusted for, i.e., the value of the recorded pressure in artesian/subartesian systems has been reduced to the value of pressure under thermal equilibrium with the surrounding formation in the borehole [4].
The key parameter for correcting the pressure readings at the wellhead is the average density of the water column in the borehole when it is in thermal equilibrium without exploitation. The quality of the data obtained is crucial here. The expected temperature of water in the reservoir and at the wellhead is usually the first criterion determining the direction of further conceptual work for selecting the perspective areas for thermal water production [5]. Aliyu and Chen [6] present an approach for predicting the long-term performance of a deep geothermal reservoir using multiple combinations of various reservoir parameters. The finite element method and factorial experimental design were applied to forecast which parameters had the greatest influence on long-term reservoir productivity. The solver employed was validated using a known analytical solution and experimental measurements with good agreement. After the validation, an investigation was then performed based on the Soultz lower geothermal reservoir. The results showed that the fluid injection temperature was the parameter that influenced the experiment the most during exploitation involving production temperature, whereas the injection pressure rate had a more significant impact on reservoir cooling [6]. The actual geothermal fluid temperature is a distinctive parameter determining the selection of the type of geothermal power plant, which is conditioned by thermodynamic and economic performance [7,8]. Temperature and pressure are also the most significant parameters to determine the feasibility of geothermal power generation using low-flow-rate downhole heat exchangers [7]. The utilisation of geothermal energy, which is a renewable and environmentally friendly source of energy, brings noticeable improvements in air quality as evidenced by the significant ecological effects achieved after commissioning geothermal district heating plants [9][10][11]. Utilization of geothermal energy depends to the largest degree on Energies 2020, 13, 3882 3 of 18 the main parameters (e.g., temperature, pressure, efficiency) so a correct assessment of them is necessary. Geothermal waters in southern Poland (especially in the Podhale Basin) at about a 3-3.5 km depth have temperatures reaching almost 90 • C and low mineralisation values below 3 g dm −3 [12]. These parameters make their technical exploitation economically viable, ensuring relatively good operating conditions [5,10,13].
This study presents the results and analysis of hydrogeological studies carried out in the CH PIG-1 borehole, which provide the most up-to-date information on the thermal water properties accessed through this borehole in the context of the thermal lift effect observed during its operation. The methodology presented may provide inputs for conducting standard hydrodynamic tests-so-called "good practice" for similar studies in geothermal boreholes. The variability of physico-chemical properties of water exploited by the CH PIG-1 in the 2016-2019 multiannual period in the context of system operation was also analysed. Water samples were taken every time from the test cock directly at the intake. Analyses were performed in the Hydrogeochemical Laboratory of the Hydrogeology and Engineering Geology Department at the AGH University of Science and Technology in Krakow. Laboratory was accredited by the Polish Centre for Accreditation on 14 June 2009 (Certificate No. AB 1050) and has implemented a quality management system in accordance with the requirements of the ISO 17025 standard [14]. The stationary observations available with respect to the water extracted made it possible to analyse correlations between the pH value and total dissolved solids (TDS), wellhead pressure, and extraction temperature, and to determine the possibility of relationships observed in the context of the measured wellhead pressure. This paper highlights a number of factors of key importance in the assessment and forecasting of operational resources of thermal water intake. Furthermore, the analysis of the correlation between hydrodynamic and physico-chemical parameters of exploited water is presented, which also plays an important role in assessing the possibility of safe water exploitation.

Description of the Studied Area
The CH PIG-1 well is located in Witów in southern Poland. It was constructed during the 1989-1990 period on the basis of a design which envisaged the drilling of nine deep boreholes in order to extract thermal water in the Podhale Basin in southern Poland [15] (Figure 1). Of these, six boreholes were ultimately completed, including the CH PIG-1 borehole [16]. In the following years, further boreholes were drilled to extract geothermal water. For many years, the Podhale Basin has been commonly known as one of the most promising structures in Poland from the point of view of recreation and balneotherapy [17][18][19][20][21][22][23][24]. The extraction of thermal waters in the area is also an important attraction boosting the development of tourism [25]. The Paleogene Podhale Basin, within which the CH PIG-1 borehole was drilled, is situated in the northern part of the Inner Carpathians, sandwiched between the Tatra Mountains in the south and the Pieniny Rock Belt in the north [27]. The area of the basin within Polish borders is approximately 490 km 2 and it extends latitudinally for about 40 km from the west to the east. There are two main geological and structural units in the Tatra region: the Tatra Massif (lower-order unit) and the Podhale Basin (upper-order unit) [28].
The waters of the Quaternary aquifer are associated with the presence of alluvial formations which fill beds of streams or fluvioglacial formations.
The waters of the Tertiary (Paleogene) aquifer are associated with Podhale Flysch formations The Paleogene Podhale Basin, within which the CH PIG-1 borehole was drilled, is situated in the northern part of the Inner Carpathians, sandwiched between the Tatra Mountains in the south and the Pieniny Rock Belt in the north [27]. The area of the basin within Polish borders is approximately 490 km 2 and it extends latitudinally for about 40 km from the west to the east. There are two main geological and structural units in the Tatra region: the Tatra Massif (lower-order unit) and the Podhale Basin (upper-order unit) [28].
The waters of the Quaternary aquifer are associated with the presence of alluvial formations which fill beds of streams or fluvioglacial formations.
The waters of the Tertiary (Paleogene) aquifer are associated with Podhale Flysch formations and carbonate sub-Flysch formations (Tatra Eocene). Podhale Flysch formations (shale and sandstone) are structures with a low water content. On a regional scale, as a whole, they can be considered isolating formations. They are fractured aquifers, and less often fractured-porous ones. The water table is most often confined [18].
Mesozoic aquifers are mostly of the karstic-fissured and karstic types and are associated with carbonate series. The precipitation recharge area of sub-Flysch aquifers is situated in the Tatra Mountains. Rainwater infiltrating in this area moves through a system of fissures and karstic voids towards the north, i.e., in the direction in which the Tatra series subside ( Figure 2). The circulation of some of this water is fairly short, since it appears on the surface in the form of springs or is drained through valleys. The remaining part penetrates deep into the ground, is heated up, and extracted through boreholes as thermal water. Early isotope studies indicated a relatively young age of the water within 100-2000 years, which indicates high intensity of water exchange in the Podhale Basin [29]. A detailed study [30] indicates that the waters in the eastern part of the Podhale Basin are in the zone of slower flows (e.g., due to the existence of sealing faults) and may be much older than the Holocene waters located closer to the outcrop.
The CH PIG-1 borehole reached a final depth of 3572 m (Figures 2 and 3). The composition of the extracted water is 0.11% SO 4 -Ca-Mg-Na. The extracted water mineral content at the outflow point ranges from 1050 to 1300 mg/dm 3 and this is natural long-term variability. The water temperature at the wellhead is 89.8 • C at Q = 160 m 3 /h. The CH PIG-1 borehole reached a final depth of 3572 m (Figures 2 and 3). The composition of the extracted water is 0.11% SO4-Ca-Mg-Na. The extracted water mineral content at the outflow point ranges from 1050 to 1300 mg/dm 3 and this is natural long-term variability. The water temperature at the wellhead is 89.8 °C at Q = 160 m 3 /h. The CH PIG-1 geothermal borehole extracts water from a confined aquifer. Due to favorable hydrodynamic conditions, the well is operated without any mechanical devices as an artesian well. Construction of the CH PIG-1 thermal water well is presented in Figure 3.  The CH PIG-1 geothermal borehole extracts water from a confined aquifer. Due to favorable hydrodynamic conditions, the well is operated without any mechanical devices as an artesian well. Construction of the CH PIG-1 thermal water well is presented in Figure 3.

Hydrodynamic Testing
The tests related to the assessment of operating parameters for the CH PIG-1 borehole in the context of the thermal lift effect were carried out in two phases on 27 and 28 November 2017. Details of the tests are presented in Table 1.

Hydrodynamic Testing
The tests related to the assessment of operating parameters for the CH PIG-1 borehole in the context of the thermal lift effect were carried out in two phases on 27 and 28 November 2017. Details of the tests are presented in Table 1.
Phase I was conducted with an average pumping capacity of 160.4 m 3 /h and included 2 h of pumping and 2 h and 5 min of stabilisation. Phase II was conducted with an average pumping capacity of 151.9 m 3 /h, with a pumping time of 2 h 15 min and a stabilisation time of 1 h 15 min. The well in which the tests were carried out is continuously operated in connection with the needs of the recreational complex. The purpose of the designed works was to determine the current operating parameters and the possibility of increasing the borehole capacity to 160 m 3 /h, which is dictated by further investment plans, i.e., extension of the aqua park (construction of, among others, a hotel facility). Prior to conducting proper hydrodynamic tests on 27 November 2017, the well was shut down for 2 h. In the course of the tests, during artesian flow of water from the well, continuous automatic measurements of the temperature of the flowing geothermal water, pressure at the wellhead, and artesian flow rate were conducted. Individual elements of the measuring system were connected to a computer system Energies 2020, 13, 3882 6 of 18 which recorded the parameters tested on an ongoing basis. In order to determine reduced wellhead pressure, it was necessary to adjust for the thermal lift effect. Temperature measurements in the CH PIG-1 borehole were carried out under thermally stabilised conditions. The value of the reduced wellhead pressure was determined using Equation (1) [1,2]: The reduced wellhead pressure value was used during the interpretation of hydrodynamic tests. To interpret test results, the Theis method as modified by Agarwal [33,34] was used with respect to the borehole pressure recovery phase. The choice of this method was determined by the fact that the flow rate was constant at that time (Q = 0), which is the basic assumption underlying tests carried out under transient inflow conditions [3,35,36]. The time here is defined as t/∆t , where t is the total duration of the test at a given stage (pumping + pressure recovery), and t is the pressure recovery time (the time that elapses from the start of the pressure recovery phase).
Non-linear flow resistance, which is irrelevant when considering the pressure recovery period (Q = 0), was not taken into account in the calculations. The hydraulic conductivity values obtained were converted into filtration coefficient values (k). For calculation purposes, the thickness of the reservoir layer accessed was assumed to be 329 m [37].

Hydrogeochemical Testing
Water samples from the borehole were taken at the outflow point in accordance with the Polish Standard PN-ISO 5667-11:2004 [38], the guidebook commonly used in Poland in similar analyses [39] and the practical guidelines described by Zdechlik et al. [40]. Once a month, an analysis of the main and other selected ions and chemical compounds was carried out at the accredited laboratory of the AGH University of Science and Technology in Krakow. The database includes results from the year 2016 (the beginning of the borehole operation) until the end of year 2019. In 2019, electrical conductivity (EC), pH, wellhead pressure, and temperature were measured continuously using gauges installed at the outflow. The reading of the pressure value was made on the basis of the indications of a pressure gauge installed on the wellhead. Individual elements of the measuring chain were connected to a computer system recording pumping parameters on an ongoing basis.
The basic analysis of the variation in concentrations of selected ions was carried out using the STATISTICA platform and was presented in the form of typical box charts. Additionally, correlations between the parameters observed were determined. The convention proposed by Guilford [41] was used to interpret the correlation coefficients of the sample.

Results and Discussion
The hydrodynamic tests carried out in the CH PIG 1 borehole and the analysis of variability of selected ions and parameters in exploited water were necessary to assess the possibility of increasing the efficiency of the CH PIG-1 borehole and to determine the water quality and its natural variability. Such information is crucial for the functioning of the recreational complex based on the use of geothermal water. A study of the phenomena affecting the exploitation of hot water from deep boreholes enables their effective exploitation and enables the use of resources in accordance with the expectations of investors.

Hydrodynamic Testing
Changes in pressure, flow rate, and temperature during hydrodynamic testing in the CH PIG-1 well are shown in Figures 4 and 5. In Figure 4, the brown line marks the wellhead pressure value which would have been recorded if there had been no volumetric expansion of water in the borehole due to the thermal lift effect. This assumption forms the basis for using existing analytical solutions to estimate the hydraulic parameters of the aquifer. The thermal lift phenomenon depends on the pressure condition and formation temperature-the higher the values of these properties, the more significant the thermal lift effect and the greater its impact on the correct interpretation of results. This distortion is clearly visible in Figure 4 (orange line-pressure measured), where despite pressure recovery being recorded in the borehole, the graph shape characteristic of this process could not be recorded. Figure 4 also demonstrates the fairly rapid stabilisation of pressure after each of the individual pumping stages in both hydrodynamic tests.
As a result of the hydrodynamic test carried out, the results shown in Table 2 were obtained. Semi-logarithmic graphs were used to interpret the results [3,42]. In the case of the calculation results of the hydrogeological properties, the maximum water extraction rate and the permissible lowering of the water table (head pressure) were determined on  table  in the CH PIG-1 well was assumed to correspond to a water table lowered by 145.5 m in relation to the  established water table, i.e., +155.46 m a.g.l. (above ground level). This means that depression had to be maintained above the 799.5 m a.s.l. (above sea level) (+10.0 m a.g.l.) as shown in Figure 6. Static water level fluctuations are caused by exploitation from many others boreholes in the Podhale Basin.   As a result of the hydrodynamic test carried out, the results shown in Table 2 were obtained. Semi-logarithmic graphs were used to interpret the results [3,42].   As a result of the hydrodynamic test carried out, the results shown in Table 2 were obtained. Semi-logarithmic graphs were used to interpret the results [3,42].  The hydraulic conductivity values obtained by using the Theis method [33] as modified by Agarwal [34] were converted into filtration coefficient values (k). For the CH PIG-1 well operated at 160 m 3 /h at a depression of s = 145.5 m, the range of the hydrodynamic changes (R-depression cone) theoretically determined using Formula (2) is 421 m. The average filtration coefficient is 9.32·10 −7 m/s.
Results of the tests conducted to determine the basic hydraulic parameters of the aquifer accessed via the CH PIG-1 well using the Theis method [33] as modified by Agarwal [34] are presented in Table 3. Interpretation results are graphically presented in Figures 7 and 8. Energies 2020, 13, 3882 9 of 18 extraction rate and the permissible lowering of the water table (head pressure) were determined on the basis of the results of hydrodynamic tests. The interpretations and calculations of exploitable water resources conducted with respect to the CH PIG-1 borehole made it possible to determine its operating flow rate: Qe = 160 m 3 /h = 3840 m 3 /d. The maximum decrease in the dynamic water table in the CH PIG-1 well was assumed to correspond to a water table lowered by 145.5 m in relation to the established water table, i.e., +155.46 m a.g.l. (above ground level). This means that depression had to be maintained above the 799.5 m a.s.l. (above sea level) (+10.0 m a.g.l.) as shown in Figure 6. Static water level fluctuations are caused by exploitation from many others boreholes in the Podhale Basin. The hydraulic conductivity values obtained by using the Theis method [33] as modified by Agarwal [34] were converted into filtration coefficient values (k). For the CH PIG-1 well operated at 160 m 3 /h at a depression of s = 145.5 m, the range of the hydrodynamic changes (R-depression cone) theoretically determined using Formula (2) is 421 m. The average filtration coefficient is 9.32·10 −7 m/s.
Results of the tests conducted to determine the basic hydraulic parameters of the aquifer accessed via the CH PIG-1 well using the Theis method [33] as modified by Agarwal [34] are presented in Table 3. Interpretation results are graphically presented in Figures 7 and 8.      Hydrogeological surveys carried out in the CH PIG-1 borehole made it possible to verify its exploitable resources. This parameter was first determined in 1997 as 190 m 3 /h [36] and was subsequently corrected in 2009 to 120 m 3 /h (estimated based only on the analysis of pressure recorded since 2006 in the wellhead) [30]. In 1997, hydrodynamic tests were carried out as part of wider research in all boreholes within the Podhale Basin. At that time, the production resources for the CH-PIG 1 borehole were established at 190 m 3 /h. The decision was issued for 10 years, so it expired in 2007. In the new documentation, the resources were corrected to 120 m 3 /h only on the basis of a few pressure measurements from 2007-2009, as there was no exploitation and no other results. CH PIG-1 was not in operation until 2016. Normal exploitation only started in spring 2016, and it was then found that after over one year of borehole operation, the static pressure in the heated borehole was 5 bar above the level recorded in 2009 in the cooled down borehole. Therefore, the production resource figures were adjusted again. Earlier interpretations did not take the thermal lift effect into account and therefore were burdened with a significant error. At present, exploitable resources have been established at 160 m 3 /h at a depression of 145.5 m, which corresponds to an elevation of H s = 799.5 m a.s.l. If the borehole is operated at the designated flow rate, the value of hydrodynamic changes determined is 421 m. The calculated value of the filtration coefficient for the aquifer accessed through the borehole is 9.32 × 10 −7 m/s, which is within the range typical of the parameters of the Podhale Basin geothermal reservoir in question [18,21].   Table 4 presents the physico-chemical characterisation of geothermal water samples from the CH PIG-1 well with the parameters analysed shown. Table 4 only shows those parameters for which the number of determinations enabled reliable statistical analysis in the next stage.

Hydrogeochemical Testing
The available database for 2016-2018 included monthly measurements of the levels of main ions, TDS, SO 4 , and H 2 S performed in the accredited Hydrogeochemical Laboratory of the Hydrogeology and Engineering Geology Department at the AGH University of Science and Technology in Krakow. H 2 SiO 3 and potassium contents were measured approximately once in two months; ferrumiron and manganese were measured quarterly in the same laboratory. Using the STATISTICA software, basic statistical analysis was performed with respect to the variability of the observed main ions and parameters of the geothermal waters analysed, and distribution histograms were plotted to present the frequency distribution in a graphical manner (Figure 9).  The analysis of Figure 9 indicates that the variability of chemical indicators is characterised by different empirical distributions, which are presented graphically in the form of histograms. For The analysis of Figure 9 indicates that the variability of chemical indicators is characterised by different empirical distributions, which are presented graphically in the form of histograms. For most indicators, the distributions are unimodal, except for H 2 S where the sample may be assumed to be heterogeneous. The greatest variability was observed for cations (Na + , K + , Ca 2+ ions) as well as for the hydrogen sulphide content. The results obtained may have been affected by the relatively small size of the data set shown in Table 4, since the indicators in question were determined once a month.
With respect to the data of observations of physical indicators, which were recorded in situ and for which a much larger sample of 267 results was available in 2019 (Table 4), a similar statistical variability analysis was performed using STATISTICA. The available database included the data from stationary studies from 2017-2019 carried out by employees of the "Chochołowskie Termy" Thermal Water Mining Plant: determination of the wellhead pressure-the reading of the pressure value was made on the basis of the indications of a pressure gauge installed on the wellhead, -measurement of the water temperature at the outlet point by means of an automatic thermometer with a measurement accuracy of ±0.1 • C, -measurement of the conductivity (EC) of water at the outlet point using an automatic probe with a measurement accuracy of ±1 µS/cm at 20 • C.
The data are presented in Figure 10.
Energies 2020, 13, x FOR PEER REVIEW 14 of 20 most indicators, the distributions are unimodal, except for H2S where the sample may be assumed to be heterogeneous. The greatest variability was observed for cations (Na + , K + , Ca 2+ ions) as well as for the hydrogen sulphide content. The results obtained may have been affected by the relatively small size of the data set shown in Table 4, since the indicators in question were determined once a month. With respect to the data of observations of physical indicators, which were recorded in situ and for which a much larger sample of 267 results was available in 2019 (Table 4), a similar statistical variability analysis was performed using STATISTICA. The available database included the data from stationary studies from 2017-2019 carried out by employees of the "Chochołowskie Termy" Thermal Water Mining Plant: determination of the wellhead pressure-the reading of the pressure value was made on the basis of the indications of a pressure gauge installed on the wellhead, -measurement of the water temperature at the outlet point by means of an automatic thermometer with a measurement accuracy of ±0.1 °C, -measurement of the conductivity (EC) of water at the outlet point using an automatic probe with a measurement accuracy of ±1 µS/cm at 20 °C.
The data are presented in Figure 10.  The analysis of Figure 10 demonstrates that the parameters of stationary observations in 2019 exhibit unimodal left-skewed distributions. The lowest variation is observed with respect to the range of pH values recorded.
Analysing the results of the water tests performed in the CH PIG-1 well, it should be stated that the water extracted is expected to be characterised by a constant chemical composition and good bacteriological status. In the previous investigation, Operacz [43] showed that the appropriate operation of the CH PIG-1 borehole allows a proper temperature to be maintained, which underpins the functioning of the Chochołowskie Termy recreation and balneological complex. Analysing the results of the water tests performed in the CH PIG-1 well, it should be stated that the water extracted is expected to be characterised by a constant chemical composition and good bacteriological status. In the previous investigation, Operacz [43] showed that the appropriate operation of the CH PIG-1 borehole allows a proper temperature to be maintained, which underpins the functioning of the Chochołowskie Termy recreation and balneological complex.
Since the temperature of exploited thermal waters is the essence of this article, it was decided to check whether the temperature in the studied borehole had an effect on other important parameters of the water from the intake. Relationships between the parameters of the geothermal waters extracted Since the temperature of exploited thermal waters is the essence of this article, it was decided to check whether the temperature in the studied borehole had an effect on other important parameters of the water from the intake. Relationships between the parameters of the geothermal waters extracted by the CH PIG-1 well determined during stationary observations in 2019 are presented under the assumption that temperature values are independent variables and the other parameters are dependent variables. The relationships were determined using Pearson's linear correlation coefficient. The values of the correlation coefficient obtained in scatter graphs are shown in Figure  11.   Pearson's correlation coefficient ranged from the lowest absolute value of −0.08 for total dissolved solids to the highest absolute value of −0.67 for the wellhead pressure. For the pH value, the correlation coefficient was 0.15.
As a result of the analysis, on the basis of the convention proposed by Guilford [41], a relationship was found only between the temperature measured at the outflow and the observed rare wellhead pressure. For this correlation, the coefficient was r = −0.67, which means moderate correlation. For TDS and pH, a practically insignificant relationship was observed.
For the relationship between temperature and the observed pressure at the wellhead, the value of the correlation coefficient was negative (a negative linear relationship was determined), which means that as the temperature increased, the wellhead pressure decreased This relationship can be presented using the following formula ( Figure 11): The trend line formula presented in Figure 10 indicates that with each degree Celsius increase in temperature in the CH PIG-1 well, there is a pressure drop at the wellhead equal to approx. 0.2 bar. The intercept of the trend line equation does not affect the pressure variation because it is not statistically significant. The significance of the intercept was analysed by Student's t-test.
The procedure presented in the article and the results obtained should constitute good practice in the evaluation of exploitable resources of geothermal water wells, of which evaluation should precede the subsequent stages of, inter alia, designing geothermal heat and power plants. In both cases, it is extremely important to analyse specific hydrogeothermal and geological conditions to assess the possibility of using geothermal resources and effective utilization [44,45]. Bujakowski et al. [18] argue that correct estimation of the operational parameters of the thermal water intake plays a key role in ensuring the possibility of long-term exploitation of thermal waters for heating purposes but also for recreational purposes. The importance of the stability of the chemical components of thermal water, used for recreational and balneological purposes, is emphasised by Kępińska i Ciągło [22]. The stability of chemical indicators in the water is important, not only because of the therapeutic qualities of thermal water (temperature, hydrogen sulphide, iodine, methacic acid contents), but also due to the issues related to the exploitation and disposal of water. According to Tomaszewska et al. [12], factors such as temperature, pH and the content of dissolved components in water, as well as operating parameters, including water pressure, have a decisive effect on the thermodynamic state of water and the tendency of precipitation of secondary (mineral) deposits from water (scaling). Scaling of the geothermal system, including the piping network, can cause immense problems with the operation of the system, as pointed out byŞimşek et al. [46]. Furthermore, the main parameters determining the power rating and the amount of energy available are the temperature of the geothermal water at the wellhead and the flow rate of geothermal water at the ground level [7,15,19,47,48]. As indicated by the authors of the aforementioned publications, the temperature of geothermal water is not the sole parameter determining the suitability of geothermal waters for energy purposes, although it is a crucial one. As already mentioned, the temperature of geothermal water feeding, e.g., a geothermal power plant, is important for the plant's power rating and efficiency, as both parameters increase with the temperature of the geothermal water [49]. Regarding changes in the flow rate of geothermal water at the ground level, the amount of heat supplied to the technological system increases with this value. This directly increases the amount of energy that can be obtained, e.g., using heat exchangers in a geothermal heat plant or an evaporator in a geothermal power plant. For this reason, the presented issues together with the proper estimation of the geothermal intake parameters and the exploitation capacity of the intake are of great importance and constitute the key engineering issues in the field.
Both the plant's power rating and the amount of heat or electricity that can be generated are parameters that have a direct bearing on the assessment of the economic viability of the potential investment project as well as on the environmental effect which can be achieved compared to methods based on conventional fuels [50,51]. This is particularly important if sustainable development is to be achieved in areas where no district heating networks are present and access to infrastructure such as gas mains is limited. Such absence of infrastructure results in increased environmental pollution (mainly air pollution) caused by the burning of conventional fuels. In such circumstances, geothermal energy may be the only prospect of improving this unfavourable state of affairs, since it offers a stable, environmentally friendly, and renewable energy source [52,53]. However, it is necessary to correctly estimate both the economic and environmental effects, and this is impossible without having precise information about the temperature and usable flow rate of geothermal water at the ground level. Furthermore, recreational water use and assessment of the stability of the physico-chemical and exploitation properties of the intake should always be correlated in order to ensure long-term and trouble-free operation.
One of the primary purposes of hydrodynamic tests is to determine the hydraulic parameters of the aquifer. The results of the tests are often flawed as a consequence of rejection of thermal heating of a well during pumping which leads to aberrant readings of the water level or wellhead pressure. A product of this is an erroneous higher hydraulic conductivity coefficient, which may lead to incorrect assessment of the admissible volume of extracted groundwater. The effect of thermal heating of a well, also called the thermal lift effect, is more important given a greater depth of a well and a greater temperature difference in the well's profile. Bielec and Miecznik [1] presented equations that allow us to calculate the above affect, including the sample analysis of the data from a hydrodynamic test carried out in a 2000 m deep geothermal well. In this case, the maximum difference between the observed and reduced wellhead pressure was 0.172 MPa, i.e., approximately 17.3 m of the water column. The difference in the hydraulic conductivity coefficient calculated for both cases ranged from 216.7% (steps II and III of pumping) to 319.4% (step I of pumping). To summarize the analyses shown above, the conductivity of exploited geothermal water depends on the pumping level and could be two to three times higher with temperature correction compared to results without any correction.

Conclusions
In deep boreholes used to access geothermal waters, it is extremely important to take the thermal lift effect into account. Disregarding this effect invalidates analysis results and results in incorrect determinations of basic hydrogeological parameters. The article points out the most important aspects of the effect in question and provides an example of correct analysis for the CH PIG-1 well in southern Poland.
The thermal effect is very important in the practice of documenting the exploitation resources of thermal waters in Poland and is widely used and recommended by legislators. Thus, in applicable procedures and methodological guides, an obligation to properly define the thermal effect is imposed. This scientific article can be an excellent example of conduct where scientific considerations can and should be used in practical activities when documenting thermal water intakes. Hydrogeological properties of the thermal water reservoir, including primarily hydraulic conductivity, calculated on the basis of hydrodynamic tests without taking into account the thermal correction are overstated. This is also of practical importance when determining the extent of the mining impact (demarcation of the mining area). The final thermal lift effect does not directly affect production from the CH PIG-1 borehole but is essential for the proper determination of its deposit parameters. Thus, under the conditions prior to starting the exploitation of geothermal waters, i.e., at the stage of preparing hydrogeological documentation and obtaining a concession for the extraction of a mineral, deposit parameters were correctly established, inter alia, the resource area of the borehole and area of influence constituting the scope of the R. On this basis, the operator has correctly designated the mining area covered by the concession, where operator is the sole entity with the authorization to exploit geothermal waters.
In order to make the study more thorough, the basic statistical variability of ions, chemical compounds, and the parameters observed for the extracted geothermal waters was analysed. The results obtained showed a natural hydrogeochemical variability of the water intake. The value ranges obtained are typical of geothermal waters, especially of hot waters extracted in southern Poland (in the Podhale Basin).
Another extremely important matter is confirming the relationship between temperature and the observed borehole pressure (measured at the wellhead) on the basis of correlation analysis. This relationship was described by a regression equation. The relationship confirmed and supplements Energies 2020, 13, 3882 16 of 18 the existing information-as the temperature of water in the borehole increases, the pressure observed in the borehole decreases. The application of Pearson's linear correlation method made it possible not just to confirm the phenomenon as such but also to determine a specific numerical relationship.
The procedure presented in this article, which was followed to correctly determine the basic parameters of the extracted geothermal waters, should be mandatory for all similar installations. Only the application of appropriate methodology and a thorough analysis of both hydrodynamic test results and water physico-chemical parameters allows for proper operation of such wells. Temperature, which is often regarded as the most important parameter of geothermal waters, determines their potential uses. Therefore, its precise determination and knowledge of the correlation between temperature and reservoir pressure are extremely important issues in geothermal studies.