Applying Petroleum the Pressure Buildup Well Test Procedure on Thermal Response Test — A Novel Method for Analyzing Temperature Recovery Period

The theory of Thermal Response Testing (TRT) is a well-known part of the sizing process of the geothermal exchange system. Multiple parameters influence the accuracy of effective ground thermal conductivity measurement; like testing time, variable power, climate interferences, groundwater effect, etc. To improve the accuracy of the TRT, we introduced a procedure to additionally analyze falloff temperature decline after the power test. The method is based on a premise of analogy between TRT and petroleum well testing, since the origin of both procedures lies in the diffusivity equation with solutions for heat conduction or pressure analysis during radial flow. Applying pressure build-up test interpretation techniques to borehole heat exchanger testing, greater accuracy could be achieved since ground conductivity could be obtained from this period. Analysis was conducted on a coaxial exchanger with five different power steps, and with both direct and reverse flow regimes. Each test was set with 96 h of classical TRT, followed by 96 h of temperature decline, making for almost 2000 h of cumulative borehole testing. Results showed that the ground conductivity value could vary by as much as 25%, depending on test time, seasonal period and power fluctuations, while the thermal conductivity obtained from the falloff period provided more stable values, with only a 10% value variation.


Introduction
The number of installed systems utilizing shallow geothermal energy has constantly increased over the last decade.Although it is common practice to model borehole heat exchangers grid of power less than 30 kW [1] with various software that estimate thermogeological property of the soil, Thermal Response Test (TRT) is the only justifiable method for obtaining the correct properties for certain locations.
The catalogue properties of soil, which are the source of data for the software modeling phase, often lead to undersized systems and functional problems.Soil sample analysis in the laboratory is too expensive for everyday use.Considering this, TRT is also the most profitable method for modeling long-term operating systems.Inaccuracies in the estimation of thermal properties from TRT are mostly caused by outdoor effects.It is necessary to recognize and reduce the atmospheric influences, as well as changes in ambient temperature.Energy dissipation from the pipes above the ground results in a higher value of estimated thermal conductivity, and may lead to the undersizing of the future system installation.
Steep changes caused by the demand of consumers in the local electrical grid will cause fluctuation of the power provided for electric heaters in TRT equipment, leading to variation in the rate of heat transfer to the borehole fluid.This directly affects the entering source temperature (EST), as well as the orderliness of the results of the measurement, and the change of temperature over time, which is a source for calculating the thermal conductivity of the found and borehole thermal resistance.Therefore, TRT should be performed using 24 h cycles, because it is easier to link discrepancies in the measurement with the daily oscillation patterns of electricity use.Prerequisites for credible calculations based on TRT measurements are described in numerous pieces of research and handbooks.
The novel method for TRT presented in this paper could achieve more precise determination of the thermal properties of the ground.It is based on theoretical and experiential methods from petroleum engineering, which has been utilized in practice for almost century, and where it is possible to identify similarities resulting from its use of the same theoretical framework for describing pressure behavior and heat conduction based on the solution of the diffusivity equation.Therefore, our recommendation is to always conduct analysis of thermal recovery period after performing classical TRT.Such prolonged TRT procedures could decrease chance of errors in interpretation due to the fluctuation of voltage in the public electrical grid and the effect of test duration on the final calculations of ground thermal properties.Furthermore, the presented method is suitable for diminishing the effects of ambient temperature interference, which can affect the values of ground thermal conductivity.This is especially true during the winter season, when there is a large difference between the air temperature and fluid inside the pipes during TRT.However, during the TRT recovery period there is a significantly lower temperature difference, due to the smaller heat flux between the surface equipment and surrounding.

Literature Overview
The Duration of Thermal Response Test is one of the most discussed parameters in interpretation procedures.The ASHRAE Handbook [1] recommends that TRT be performed for at least 36 to 48 h.Bujok et al. [2] presented measurements conducted on experimental underground heat storage.TRT was performed on eight boreholes with various test-time durations on identical ground environments.Software simulation of thermal conductivity showed that data recorded for the first 24 h of TRT deviated by up to 7.8% compared to the results of TRT provided for 70 h.Furthermore, borehole thermal resistance value obtained from 24-h long TRT differs by up to 17.9% from values calculated based on 70-h long TRT.
The least discussed influence on the credibility of TRT is the variation of ambient temperature, which can strongly affect the measured temperature results of circulating fluid.Bandos et al. [3] first introduced a method whereby the effect of atmospheric conditions is subtracted by using air temperature data for the time at which the test was conducted.Application of this method reduces the oscillation of the thermal conductivity value from 30% to 10% of the mean value.It has been shown in research that the delay of the ambient to mean fluid temperature is about 3 h.
Borinaga-Trevino et al. [4] developed a method to reduce the influence of atmospheric conditions on the TRT.The main advantage is that it is not necessary to know its physical origin; rather, it is based on analyzing the influence of the chosen time interval in order to fit the data of the infinite line source theory (ILS) when predicting the ground thermal conductivity.Two TRTs were analyzed, each with different equipment and levels of isolation from the environment.They concluded that poorly isolated TRT equipment that is exposed to variable ambient conditions can lead to an error of ±33% in determining the thermal conductivity of the ground, depending on the climate and season.
Signorelli et al. [5] tested the duration of the TRT for the evaluation of ground thermal properties, comparing the results from 3D finite element numerical models with the results from synthetic line source solutions.Research has focused on the estimation of the time t0, which denotes the point in time after the start of the test after which the data will not be affected by the lower thermal conductivity of the tubes and grouting material and the effect of unsteady heat transfer.In this case, the simulated 200 h test in the numerical model was evaluated for t0 values of 10 h, 20 h, 40 h, and 60 h, and variable values for the end of the test, tE, to the extent t0 < tE < 200 h.According to the results, accurate definition of t0 has a more significant impact on values of calculated thermal conductivity than the total duration of test.
A novel method for the determination of the duration of the unsteady state was presented by Kurevija et al. [6], in which a derivation curve-a method from well testing in petroleum engineering-was applied to recorded temperature data.It is based on monitoring the change of temperature for certain small periods of time versus cumulative time (dT/dt vs. time).The dampening temperature change curve precisely shows the time after which the change of temperature drops below 0.25 • C per 5 min interval, as an arbitrarily established value.Such determination of the beginning time of semi-steady-state heat transfer is much more accurate than the standard method from Gehlin [7], where the value of thermal diffusivity needs to be assumed.
Badenes et al. [8] investigated the influence of operational parameters during TRT, distinguishing the importance of heat rejection control.Research shows that, when a proportional-integral-derivative (PID) controller is used, the impact of environmental temperature fluctuations is decreased, fluctuations of inlet and outlet temperature are negligible, and values of thermal conductivity and borehole thermal resistance are closer to the values corrected using the least-squares fitting algorithm method.
The ASHRAE Handbook [1] recommends that the standard deviation of input power, which directly affects the rejected heat to the ground from electric heaters, should be less than ±1.5% of the average value, with peaks less than ±10% of the average.Heaters can be powered either from the local electricity grid or from a generator.Demand for electricity changes during the day, according to households' habits and industrial consumption; but mostly, regular patterns can be found.Peak demand happens during the morning, where huge amounts of energy must be delivered to users, generally around 7-9 a.m.Another peak usually happens between 4 and 9 p.m. [9].
Witte [10] conducted research into the causes of errors in TRT interpretation.It was found that the average temperature of the fluid, as the arithmetic mean of the entering source temperature and leaving source temperature, was correct only when the heat flux was constant along the entire borehole.This fact confirms the importance of achieving semi-steady heat transfer during TRT.Most of the analysis was conducted using infinite line source theory, solely based on the conduction of the heat around the borehole.The ASHRAE Handbook [1] recommends heat rates from 40-80 W/m, as those are expected to be equal to peak loads of the actual heat pump system.Furthermore, below this threshold, thermally induced convection can occur in pipes, which does not fit with ILS assumptions.
Unexpected events, such as electrical power outage, can interrupt TRT before a time duration sufficient for properly estimating the thermal properties of the ground has been reached.ASHRAE [1] recommends that borehole temperature be allowed to return to within 0.3 • C of the initial ground temperature before conducting any new tests.It is assumed that it will require 10 to 12 days in mid-to-high-conductivity formations and 14 days in low-conductivity formations, in the case of a 48 h initial TR test.
Raymond [11] et al. analyzed TRT and improved it from a hydrogeological perspective, based on the concepts applied to pumping test analysis.The radius of influence was used to evaluate the duration of a TRT, prior to conducting the test.The sensitivity analysis showed that the uncertainty related to thermal properties is reduced by first using temperature recovery data to evaluate the subsurface thermal conductivity, and then using heat injection data, which can be extended with converted recovery measurements, to independently determine the borehole thermal resistance.

Solutions of the Diffusivity Equations for the Case of Infinite Medium and Line Source Well
TRT analysis, in most cases, is based on the assumption that a borehole heat exchanger of sufficient length and negligible radius can be regarded as a line source.In this case, Kelvin's line source theory analytical solution can be used to describe radial heat transfer and the behavior of temperature as a function of time and radius during the test [12][13][14].This theory originates from Fourier's law, which describes the heat conduction in solids, and which was later applied in solutions of pressure distribution during radial fluid flow in porous media (1), since both are derived from the diffusivity partial equation.For a case of radial flow into a well in a homogenous, isotropic and porous medium, and for a fluid of small and constant compressibility, constant viscosity and small pressure gradients, the diffusivity equation can be written as: where k φµc t = η is the hydraulic diffusivity factor in petroleum engineering.For a case of heat conduction: where λ cρ = α describes the thermal diffusivity factor in thermogeology.Due to their having the same theoretical origin, it is possible to correlate some practical methods from petroleum engineering with thermogeological engineering.Based on this assumption, the aims of the present research arise: to implement a novel improved TRT method based on petroleum engineering well testing methods.Since these methods are based on the same model of the diffusivity equation for describing the behavior of pressure as a function of time and radius during the pressure build-up test, the behavior of temperature in borehole heat exchanger during TRT temperature recovery phase should, analogously, also be able to be described.
The solution to the diffusivity equation for the line source theory is well known in petroleum engineering, and has been utilized for almost a century for pressure transient analysis in drawdown well testing: In geothermal applications, infinite line source theory is used to describe temperature change as a function of time and borehole radius during heat exchange between the borehole and the surrounding ground.In cases of heat extraction from the ground during the heat pump heating phase: In cases of heat rejection to the ground during the heat pump cooling phase or TRT: The term Ei represents the exponential integral function with the solution: In cases where x < 0.02, or term 4αt r w 2 > 50, the exponential integral can be approximated with a natural logarithmic function, with an error of less than 0.6% [14], as follows: where γ represents Euler's constant, which is equal to 0.5772; and e Euler's number, which is equal to 2.7183.Considering this, Equation ( 4) can be derived for heat extraction from the ground as follows: T(r w , t) ext = T i + q 4πλ ln e γ rw 2 4αt = T i + q 4πλ ln rw 2 4αt + 0.5772 = T i + q 4πλ ln rw 2 αt − 0.8097 (8) Analogously, for heat rejection to the ground, Equation ( 5) becomes: T(r w , t) rej = T i − q 4πλ ln e γ rw 2 4αt = T i − q 4πλ ln rw 2 4αt + 0.5772 = T i − q 4πλ ln rw 2 αt − 0.8097 (9) where κ = q 4πλ is the slope of the line when plotting T vs. ln(t), and is a standardized principle for obtaining effective ground thermal conductivity.

Application of Horner's Pressure Build-Up Method in Applied Thermogeology and TRT
The classic TRT procedure is analogous to the well testing in petroleum engineering practice.However, flow tests are rarely used for the determination of formation properties because of the long periods required to achieve a fully developed transient flow region.The pressure buildup test is the most frequently used test for determining formation permeability (equivalent to thermal conductivity) and near-well damage formation or skin factor (equivalent to borehole thermal resistance).Considering the similarities in the origin of the equations describing the behavior of pressure and temperature as functions of time and radius suggests the following hypothesis for the proposed novel method of TRT: "the thermal properties of the ground and borehole thermal resistance (skin) could also be estimated by observing the temperature recovery period after classical TRT, which is analogous to observing the pressure build-up (recovery) after achieving constant flow rate in petroleum well testing".
Petroleum pressure build-up test method imply constant rate of production and a stabilized bottom-hole pressure for a certain time, after which the well is closed from the surface, allowing pressure to build up/recover.During the recovery phase, bottom-hole pressure is monitored as a function of time.The acquired data is used to estimate the physical properties of the reservoir and skin, indicating any damage and disturbance of the initial near-wellbore formation permeability [14].
The interpretation of collected pressure data is most often performed with Horner's method [15], in which the measured pressure and (tp + ∆t)/ ∆t are plotted on semi-log paper (Figure 1).The term t p is the duration of the constant production rate before the well shut-in, while ∆t is the time following the well shut-in.A typical Horner's method procedure is to analyze the data and identify where the curve becomes linear on a semi-log graph, as this is the period of specific interest, from which the slope of the curve, m, can then be determined.This factor describes the pressure change (or temperature change during TRT) over one logarithmic cycle of time.The curve from semi-log paper can be arbitrary and logically divided into three regions: Early Time Region (ETR), Middle Time Region (MRT), and Late Time Region (LTR).In petroleum well testing, ETR is affected by altered (damaged) permeability of the near-well zone and the afterflow effect from the reservoir into the wellbore after shut-in.In thermogeology, and particularly TRT analysis, ETR describes the time necessary to achieve a semi-steady heat flow state from the borehole heat exchanger to the surrounding ground.Like in petroleum well testing, initial data distortion in the form of unsteady-state heat flow occurs under influence of borehole thermal resistances; this is primarily because of the low conductivity of pipe material and grout properties, as described in Section 3.4.
The MTR period starts when the investigation radius exceeds the altered zone influence near the test wells, or with the achievement of a semi-steady-state heat flow in TRT.Accurate determination of this period is crucial for the valid interpretation of the measured data.In MTR, the slope of the straight line, m, is directly dependent on the effective permeability of the formation in petroleum engineering; or, in the case of thermogeology, the slope of the line is dependent on the effective thermal conductivity of the ground.
Absolute change in pressure or temperature for one log cycle, in accordance with Figure 1, can be presented as: QBµ 4πkh in petroleum well testing (10) and, m = 2.303 q 4πλ ⇒ λ = 2.303 q 4πm in thermogeology and TRT interpretation (11) The constant factor 2.303 is the result of transferring the linear source solution with a natural logarithm into the decimal logarithm required for the use of semi-log paper in Horner's method: log e x = log 10 x log 10 e = log 10 x 0.4343 therefore ln(x) = 2.303 log(x) (12) The straight line from MRT will continue until the investigation radius reaches the natural reservoir boundaries, or pressure interference effects of other wells from the oil or geothermal reservoir appear.After this moment, deviations in pressure from a straight line can again be observed.In applied thermogeology and thermal response test, the MRT period will also cease following the appearance of interference due to the intersection of heat flow lines from other surrounding borehole heat exchangers.
The constant factor 2.303 is the result of transferring the linear source solution with a natural logarithm into the decimal logarithm required for the use of semi-log paper in Horner's method: The straight line from MRT will continue until the investigation radius reaches the natural reservoir boundaries, or pressure interference effects of other wells from the oil or geothermal reservoir appear.After this moment, deviations in pressure from a straight line can again be observed.In applied thermogeology and thermal response test, the MRT period will also cease following the appearance of interference due to the intersection of heat flow lines from other surrounding borehole heat exchangers.As seen from Figure 1, by extrapolating the linear part of the curve to the value of (tp + ∆t)/∆t = 1, the related pressure value will be equal to initial, undisturbed static pressure in the reservoir.By following the same investigative principles, the initial effective temperature alongside the borehole heat exchanger can be determined.As seen from Figure 1, by extrapolating the linear part of the curve to the value of (tp + ∆t)/ ∆t = 1, the related pressure value will be equal to initial, undisturbed static pressure in the reservoir.By following the same investigative principles, the initial effective temperature alongside the borehole heat exchanger can be determined.
Describing the pressure behavior curve during multi-flow rates of production in petroleum engineering is based on a mathematical technique called the principle of superposition.Considering the solution of line source theory for pressure behavior during pressure drawdown or flow test (3), and the constant flow rate, q, reduced to time of interest with an approximation of the value of exponential integral (7), the bottom-hole pressure drop is: after the end of production period t p .For shut-in of the well at a time ∆t after the production period t p , as required by Horner's method (Figure 2), the pressure drop at time ∆t can be obtained by the principle of superposition as: p i − p w f = (pressure drop caused by rate Q for the time t p + ∆t) + (pressure drop caused by rate change (−Q) for the time ∆t).
Pressure p wf can be replaced with static pressure p ws at a certain time ∆t after closing the well; therefore, according to Horner's method, the corresponding equation for pressure build-up becomes: which is the basic pressure build-up equation.
Considering the required decimal logarithm for plotting in Horner's semi-log paper, as shown in Section 3.2, final form of Equation ( 15) becomes: This equation confirms that if the pressure p ws observed during a shut-in period is plotted versus the logarithm of (t p + ∆t)/ ∆t, a straight line should be obtained (Figure 1).
Correspondingly, the TRT period is equal to the pressure drawdown test and the injector test (Figure 2, quadrant a1, a2 and c1/c2).After TRT is performed with a constant heat rate q for time t p , when turning off the heaters and changing the heat rate to (−q) for the time t p + ∆t, temperature recovery (equal to pressure buildup test and fall-off test) (Figure 2 quadrant b1, b2 and d1, d2) can be derived as: Energies 2018, 11, x FOR PEER REVIEW 8 of 23

Application of Petroleum Engineering Derivation Curves in TRT Interpretation
As mentioned, precise determination of the MTR period threshold is crucial for the valid interpretation of the measured data, either to determine formation permeability from the pressure data, or ground conductivity from the temperature data.In thermogeology, the currently most-used method for determining the beginning of semi-steady heat flow is by using the following formula

Application of Petroleum Engineering Derivation Curves in TRT Interpretation
As mentioned, precise determination of the MTR period threshold is crucial for the valid interpretation of the measured data, either to determine formation permeability from the pressure data, or ground conductivity from the temperature data.In thermogeology, the currently most-used method for determining the beginning of semi-steady heat flow is by using the following formula [13,[16][17][18]: The applicability of the ILS solution depends on the very large ratio between borehole depth and radius, as well as a sufficiently large Fourier number (F o = αt/r 2 b > 20).In the numerator, this equation contains a thermal diffusivity value, α, that is assumed from the drilling data or catalogue values for different types of soil.Since this is not the exact value, this method of determining the duration of the transition period can cause further error in interpretation, especially for highly heterogeneous ground.We suggest the much more accurate grapho-analytical method-the so-called derivation curve principle-which is widely used in petroleum engineering.In TRT interpretation, it could be used to precisely estimate the transition from the ETR into the MTR period, or the transition from unsteady-state heat flow to the relevant semi-steady-state heat-flow regime.The typical shape of the derivation curve for a TRT is shown in Figure 3, where each data point is derived as follows: As can be seen from Figure 3, it can be argued that the MTR period clearly appears for a value of ~0.25 • C/∆t.This threshold between ETR and MTR could be used for any type of TRT procedure, independently of power rate or test duration.The small distortions seen in Figure 3 are usual in TRT when the electric grid is used as a source of power, due to the daily fluctuations of voltage.The derivation curve could also be replaced with the simple ∆T/∆t principle.This method monitors the segment of change in the entering borehole fluid temperature against some small segment of time.The duration of this segment during TRT should not be longer than 5-10 min, as this is roughly the duration of one flow cycle trough for a coaxial or 2U heat exchanger of 100 m in depth.This method could be interpreted like a real derivation curve, because the derivative of a dot on a certain curve is a tangent to that curve.If we are looking at the change of ΔT vs. Δt on a curve, the line that connects these two dots is secant.If ΔT vs. Δt is small enough, then secant and tangent fall almost on the same line, therefore making the change of temperature against time so small as to be negligible.

Interpretation of Borehole Thermal Resistance Analogous to Well Testing Skin Effect
An inclined coaxial heat exchanger, as a pipe in a pipe-type system, was used as the geothermal The derivation curve could also be replaced with the simple ∆T/∆t principle.This method monitors the segment of change in the entering borehole fluid temperature against some small segment of time.The duration of this segment during TRT should not be longer than 5-10 min, as this is roughly the duration of one flow cycle trough for a coaxial or 2U heat exchanger of 100 m in depth.This method could be interpreted like a real derivation curve, because the derivative of a dot on a certain curve is a tangent to that curve.If we are looking at the change of ∆T vs. ∆t on a curve, the line that connects these two dots is secant.If ∆T vs. ∆t is small enough, then secant and tangent fall almost on the same line, therefore making the change of temperature against time so small as to be negligible.

Interpretation of Borehole Thermal Resistance Analogous to Well Testing Skin Effect
An inclined coaxial heat exchanger, as a pipe in a pipe-type system, was used as the geothermal source for this research.There are two principal cases of fluid circulation: Coaxial-annular or CXA, where inlet of fluid is through the annular space, and outlet is through the central pipe; and Coaxial-column or CXC, with inlet of fluid through the central pipe and outlet through the annular space.A coaxial borehole heat exchanger consists of three components (Figure 4): one pipe-in (marked as i1), one pipe-out (marked as o1), and grout material, which here is considered only as a single zone (marked as g1) [19].Heat exchange to the grout material g1 is only performed through the pipe-out o1, while the pipe-out o1 only exchanges heat with the pipe-in i1.
According to Diersch et al. [20], it is possible to divide thermal resistance, in the case of coaxial borehole heat exchangers, into three components: 1.
Thermal resistance due to the advective heat transport in the pipes and between the pipes (R ff ) 2.
Thermal resistance due to pipe wall material and grout transition (R fig , R fog ) 3.
Thermal resistance due to grout-soil exchange (R gs ) The distribution of thermal resistances is presented in Figure 4 for both coaxial borehole heat exchanger flow settings.In petroleum engineering, the permeability of the formation near the wellbore is reduced as a result of the invasion of drilling fluids into the permeable zone, dispersion of clays, and the presence of mudcake on the wellbore wall [21].The effect of the reduction in near-well permeability can be taken into account as an additional pressure drop ∆p proportional to the rate of production Q.The zone of reduced permeability has been called the "skin", and the resulting effect a "skin effect" [22].
If the value of the original formation permeability is known (from pressure build-up tests), it is possible to calculate the skin factor, which is a dimensionless factor for evaluating near-borehole damage and the causes of decreased permeability.Skin factor can be presented as an additional pressure drop [21]: If the same principle is linked to applied thermogeology, skin could represent cumulative thermal resistances to heat flow inside the borehole, as described above for the case of coaxial heat exchangers, and can be expressed as an initial temperature rise during TRT: After introducing Equations ( 19) and ( 20) into Equations ( 8) and ( 13), well pressure or borehole heat exchanger temperature could be expressed, following a production time tp, for a pressure drawdown test or extraction of heat from the ground: In petroleum engineering, the permeability of the formation near the wellbore is reduced as a result of the invasion of drilling fluids into the permeable zone, dispersion of clays, and the presence of mudcake on the wellbore wall [21].The effect of the reduction in near-well permeability can be taken into account as an additional pressure drop ∆p proportional to the rate of production Q.The zone of reduced permeability has been called the "skin", and the resulting effect a "skin effect" [22].
If the value of the original formation permeability is known (from pressure build-up tests), it is possible to calculate the skin factor, which is a dimensionless factor for evaluating near-borehole damage and the causes of decreased permeability.Skin factor can be presented as an additional pressure drop [21]: If the same principle is linked to applied thermogeology, skin could represent cumulative thermal resistances to heat flow inside the borehole, as described above for the case of coaxial heat exchangers, and can be expressed as an initial temperature rise during TRT: After introducing Equations ( 19) and (20) into Equations ( 8) and ( 13), well pressure or borehole heat exchanger temperature could be expressed, following a production time t p , for a pressure drawdown test or extraction of heat from the ground: For the case of an injector well test in petroleum engineering or heat rejection to the ground during TRT: T rej = T i − q 4πλ ln e γ rw 2 4αt In the case of heat rejection to the ground (classic TRT), the skin effect could be calculated from Equation (24), for the case of the initial conditions and the conditions just before heater shut-in.
The order of magnitude of the skin effect (or effective borehole thermal resistance) for a pressure build-up test or TRT recovery period can be estimated from the difference between the pressure or fluid temperature before shut-in and the one shortly after.By combining Equations ( 21) and (22), which give the pressure/temperature before shut-in, with Equations ( 16) and (17), which give the pressure/temperature after shut-in, it can be found that: For values of ∆t that are small compared to t p , it could be approximated that (t p + ∆t)/t p = 1.Rearranging Equation ( 25) and arbitrarily choosing a value of 1 h for ∆t, and replacing T i = T 1h , it could be rewritten [21]: This equation is one of the most important expressions in well testing.In this form, it could also be used to define the skin effect during the TRT recovery period, and to compare it with the skin derived from the classic TRT (Equation ( 25)).The value of T 1h is not strictly set, and changing the value of ∆t in the equation to values other than 1 h would merely change the constant 0.351 in Equation (27).
In an ideal case, these two values should be the same, but due to fluctuations in voltage during classic TRT and different ambient thermal interferences for the two periods, they are slightly different (different heat flux depending on the temperature difference between the air temperature and the borehole temperature during measurement).
To describe the entire curve for the two steps (TRT period + TRT recovery), including borehole thermal resistances, the principle of superposition applies.For the first TRT period, the general Equation ( 5) could be used with addition of the temperature rise term, due to the borehole skin: To fit the second part of the curve (TRT recovery period), the superposition technique is applied: where ∆t is the difference between cumulative time at a certain point and t p as the duration of the first period.

Experimental Site Setup
The ground TR tests were conducted on the coaxial borehole heat exchanger system in the city of Zagreb/Croatia, at the Faculty of Geology, Mining and Petroleum Engineering.The installation serves as the testing heat exchanger for students.
The system comprises two boreholes, each with a length of 50 m, hydraulically connected in series to provide effectively one borehole of 100 m in length.The boreholes were drilled with a standard diameter of 110 mm, and drilling was performed with specialized equipment that allows the drilling angle to be set from 35 • to 65 • , and in all directions (Figure 5b).Each borehole has an angle of 45 • , and these are placed opposite to one another inside a polyethylene shaft with a diameter of 1 m and depth of 1 m (Figure 5c).The coaxial heat exchangers comprise an outer polyethylene pipe of 63 mm with Standard Dimension Ratio of SDR11, while the inner polyethylene pipe is 32 with SDR11.Thermal response testing was conducted on both possible flow arrangements, CXC and CXA setup (see Section 3.4).Cementing with thermally enhanced grout was not possible due to the total losses of material into the high-permeability gravel layer.Therefore, the borehole was cemented with a mixture of water, bentonite and cement, with a somewhat lower thermal conductivity of 1.2 W/m K, measured with the needle probe method.
The detailed geological setting of the location, and the city of Zagreb in general, consists of Middle and Upper Pleistocene sediments, where lateral changes of gravel, silt, sand and clay are frequent, and Holocene sediments that consist of yellow-brown gravel, sand and limestone pebbles.The faculty is located in the northern part of the Zagreb aquifer, just near the outer boundary, with a thin aquifer layer present at a depth of between 6.5 and 12.0 m beneath the surface.The lithological profile, obtained from drilling data, is shown in Figure 5f.The undisturbed ground temperature and geothermal gradient at the Zagreb location were investigated in our previous research [7,[23][24][25], which demonstrated a geothermal gradient corresponding to 5.5 • C per 100 m of depth, and an undisturbed ground temperature of approximately 14.5 • C at a depth of 10 m.
As can be seen from Figure 5e, the geometry of the inclined coaxial system implies that the final depth of the boreholes is 35 m beneath the surface.Since aquifer is present in the shallow thin layer, which has a thickness of 5.5 m, the cumulative length of the pipes affected by additional convective heat transfer from groundwater is 16 m out of 100 m in total.Based on regional hydrogeological research data, the hydraulic conductivity of the aquifer near the outer boundary is only 0.3 cm/s; therefore, the convective component of heat transfer was ignored in the thermal response test.In cases of higher hydraulic conductivity, or a thicker groundwater layer in the lithological column, basic hydrogeological interpretation needs to be conducted to determine the Peclet number.In saturated coarse gravel, the Peclet number usually has a value greater than 500, which means that there is a dominant convective heat transfer; while for fine saturated sand, the value is around 1. Silts and clays give values lower than 0.01, which suggests a completely conductive heat transfer [26][27][28].
The measurements on the coaxial heat exchangers were conducted between September 2016 and June 2017, with a Geocube GC500 TRT apparatus (Precision Geothermal LLC, Maple Plain, MN, USA).The equipment has a maximum available power for the electric heaters of 9.0 kW @ 240 V.An internal logger collects 5-min interval data about inlet and outlet fluid temperature, air temperature, flow, voltage and electric current.Sealed temperature sensors (Onset Computer Corporation, Bourne, MA, USA) (resistance temperature detectors-RTD) on inlet and outlet connection have an accuracy of ±0.2 • C from 0 • C to 50 • C.
The testing procedure was organized as a classic TRT heat rejection step with a duration of 96 h, followed each time by a recovery period of an additional 96 h (only circulation).Five different heat steps were used-35, 42, 54, 61 and 71 W/m-for each of the two possible flow setups, CXC and CXA.Ultimately, the total testing time on the coaxial heat exchanger was 192 h for each of the ten different conditions, making for an accumulated 1920 h of data with 5 min logging intervals.After each testing condition, a pause was taken for a duration of seven days, allowing recovery of the ground and the borehole fluid temperature to static initial conditions.Unlike classic 2U-loop vertical heat exchangers with a depth of 100 m, inclined coaxial heat exchangers exploit shallow geothermal sources up to vertical depths of only 35 m.Hence, the initial ground temperature conditions due to different climate seasons are much more affected by the surface air temperature for this kind of installation.Additionally, in practical installations, coaxial heat exchangers are densely radially drilled from a single shaft, causing much more thermal interference between adjacent boreholes in the first few meters of depth.Table 1 presents the results of MS Excel (Supplementary Materials) descriptive statistics related to TRT heating power for each of the ten different heat steps; five for the CXC flow arrangement and five for the CXA arrangement.The acceptable maximum standard deviation for input power is +/−1.5% from the average power level, according to ASHRAE standard 1118-TRP.Furthermore, peak variations must be kept at less than 10% of the average power level.It can be seen from test data that Table 1 presents the results of MS Excel (Supplementary Materials) descriptive statistics related to TRT heating power for each of the ten different heat steps; five for the CXC flow arrangement and five for the CXA arrangement.The acceptable maximum standard deviation for input power is +/−1.5% from the average power level, according to ASHRAE standard 1118-TRP.Furthermore, peak variations must be kept at less than 10% of the average power level.It can be seen from test data that the entirety of the power data fits inside the industry guidelines provided, with the largest deviations seen for the highest power steps, as expected.Table 2 provides descriptive statistics results related to air temperature measurements for a 1-h step.The entire test procedure extended over a six-month period due to there being ~2000 h of testing, in addition to one week of waiting time for the ground to return to initial conditions.Therefore, it could be seen that the mean air temperature data are scattered for each of the heat steps, especially the minimum and maximum values.This also, to some degree, affects the initial temperature conditions for such a shallow geothermal installation, since solar perturbation affects ground up to a depth of 10 m at the test location, and the final depth of the boreholes is 35 m (Figure 5e).

Results and Discussion
As explained in Section 4, TRT + recovery period was performed for ten heat steps and two different flow arrangements.For each of the steps, a classic TRT analysis was performed under heat power conditions.Before turning on the electric heaters, the initial borehole temperature was recorded for 30 min of solely fluid circulation.The flow for each of the ten steps was set to 0.42 L/s, and the circulating fluid was pure water.Considering the coaxial pipe arrangement and dimensions (D63 SDR11 /D32 SDR11 ), the volume flow and velocity, the viscosity and density of water, and the pipe roughness, there was a fully developed turbulent regime in both the annular space (Re = 7100) and the column pipe (Re = 23,000).
Figure 6a presents the standard analysis for each step, where the entering source temperature (EST), leaving source temperature (LST), and unit heat power is plotted as a function of test time (96 h + 96 h).Since TRT is usually performed based on the principle of heat rejection into the ground, inverse mirror flow/return curves were charted to represent the cycle in which the heat pump was in heating mode (subcooling the ground).This approach better represents the working conditions during the extraction of heat from the ground (the heat pump heating cycle, for example).To determine the effective ground thermal conductivity for each of the ten analyzed cases, the emersion time of the semi-steady-state had be identified from the initial unsteady state.Based on the theory presented in Section 3.3, derivation curves were created for each of the ten TRT periods.The results are shown in Figure 7.It is evident that a semi-steady state appears after the ~10th h of investigation, where the change of the temperature per unit time reaches 0.25 °C per 5-min time step.Since most of the lithology column is made out of damp clay (except thin saturated gravel near surface), the effective thermal diffusivity, according to catalogue soil data, could be anywhere from 0.030 to 0.060 m 2 /d, depending on the moisture content.As a reasonable estimation, we assumed a value of 0.050 m 2 /d for further analysis.If the appearance of SS-state is deduced from the standardized equation by Mogensen presented in Section 3.3, then the corresponding time would be 7.5 h for a borehole diameter of 110 mm.To achieve maximum accuracy and nullify transient effects, the ground thermal conductivities were derived from intervals of 15-96 h for all test conditions in the TRT period and recovery period.To determine the effective ground thermal conductivity for each of the ten analyzed cases, the emersion time of the semi-steady-state heat flow had to be identified from the initial unsteady state.Based on the theory presented in Section 3.3, derivation curves were created for each of the ten TRT periods.The results are shown in Figure 7.It is evident that a semi-steady state appears after the ~10th h of investigation, where the change of the temperature per unit time reaches 0.25 • C per 5-min time step.Since most of the lithology column is made out of damp clay (except thin saturated gravel near surface), the effective thermal diffusivity, according to catalogue soil data, could be anywhere from 0.030 to 0.060 m 2 /d, depending on the moisture content.As a reasonable estimation, we assumed a value of 0.050 m 2 /d for further analysis.If the appearance of SS-state is deduced from the standardized equation by Mogensen presented in Section 3.3, then the corresponding time would be 7.5 h for a borehole diameter of 110 mm.To achieve maximum accuracy and nullify transient effects, the ground thermal conductivities were derived from intervals of 15-96 h for all test conditions in the TRT period and recovery period.In Figure 6a,b, the procedure of deriving ground thermal conductivities can be seen.For the case of the TRT period, semi-log axes were used, and borehole Tavg vs. tp was plotted for the Semi-Steady State or SSS interval.As stated by Equation ( 10), since a semi-log plot is used, there is the need to determine the change of temperature for a one-log cycle of time (m), and then calculate the effective ground conductivity.The same result could be obtained by using the standardized principle of plotting Tavg vs. ln(tp) on a normal graph, and then using Equation ( 9) and the slope of the line κ.
When analyzing the recovery period, a slightly different approach has to be used, as explained in Section 3.2.A reversed semi-log graph has to be used, and Tavg is plotted vs. (tp + Δt)/Δt, as shown by Figure 6c.Effective ground thermal conductivity is then derived from Equation ( 17) by knowing the log slope m, just like in the case for a standard TRT period.When the temperature recovery line is extended until (tp + Δt)/Δt = 1, initial conditions-i.e., undisturbed ground temperature-are reached.Figure 8 shows the entire Horner's procedure for all ten heat step conditions.It is important to note that testing times of 15-96 h were used in this analysis, in order to nullify transient effects, just like in the case of TRT.When extending the data trendline until (tp + Δt)/Δt = 1, it can be seen that the initial temperature value ranges between 14.0 °C and 15.2 °C.As mentioned before, the reason for this effect is the shallow installation of the coaxial system with a final depth of 35 m, while solar energy perturbation penetrates to a depth of 10 m.Therefore, the initial condition values are somewhat dependent on the climate and season in which the measurements are taking place.In Figure 6a,b, the procedure of deriving ground thermal conductivities can be seen.For the case of the TRT period, semi-log axes were used, and borehole T avg vs. t p was plotted for the Semi-Steady State or SSS interval.As stated by Equation ( 10), since a semi-log plot is used, there is the need to determine the change of temperature for a one-log cycle of time (m), and then calculate the effective ground conductivity.The same result could be obtained by using the standardized principle of plotting T avg vs. ln(t p ) on a normal graph, and then using Equation ( 9) and the slope of the line κ.
When analyzing the recovery period, a slightly different approach has to be used, as explained in Section 3.2.A reversed semi-log graph has to be used, and T avg is plotted vs. (t p + ∆t)/ ∆t, as shown by Figure 6c.Effective ground thermal conductivity is then derived from Equation ( 17) by knowing the log slope m, just like in the case for a standard TRT period.When the temperature recovery line is extended until (t p + ∆t)/ ∆t = 1, initial conditions-i.e., undisturbed ground temperature-are reached.Figure 8 shows the entire Horner's procedure for all ten heat step conditions.It is important to note that testing times of 15-96 h were used in this analysis, in order to nullify transient effects, just like in the case of TRT.When extending the data trendline until (t p + ∆t)/ ∆t = 1, it can be seen that the initial temperature value ranges between 14.0 • C and 15.2 • C. As mentioned before, the reason for this effect is the shallow installation of the coaxial system with a final depth of 35 m, while solar energy perturbation penetrates to a depth of 10 m.Therefore, the initial condition values are somewhat dependent on the climate and season in which the measurements are taking place.In Figure 6a,b, the procedure of deriving ground thermal conductivities can be seen.For the case of the TRT period, semi-log axes were used, and borehole Tavg vs. tp was plotted for the Semi-Steady State or SSS interval.As stated by Equation ( 10), since a semi-log plot is used, there is the need to determine the change of temperature for a one-log cycle of time (m), and then calculate the effective ground conductivity.The same result could be obtained by using the standardized principle of plotting Tavg vs. ln(tp) on a normal graph, and then using Equation ( 9) and the slope of the line κ.
When analyzing the recovery period, a slightly different approach has to be used, as explained in Section 3.2.A reversed semi-log graph has to be used, and Tavg is plotted vs. (tp + Δt)/Δt, as shown by Figure 6c.Effective ground thermal conductivity is then derived from Equation ( 17) by knowing the log slope m, just like in the case for a standard TRT period.When the temperature recovery line is extended until (tp + Δt)/Δt = 1, initial conditions-i.e., undisturbed ground temperature-are reached.Figure 8 shows the entire Horner's procedure for all ten heat step conditions.It is important to note that testing times of 15-96 h were used in this analysis, in order to nullify transient effects, just like in the case of TRT.When extending the data trendline until (tp + Δt)/Δt = 1, it can be seen that the initial temperature value ranges between 14.0 °C and 15.2 °C.As mentioned before, the reason for this effect is the shallow installation of the coaxial system with a final depth of 35 m, while solar energy perturbation penetrates to a depth of 10 m.Therefore, the initial condition values are somewhat dependent on the climate and season in which the measurements are taking place.Since the total duration of TRT is not strictly defined, only the minimum time is required, as explained in Section 2; in practical field testing, this always leads to a certain degree of analysis error.Since the final duration is arbitrarily chosen by the TRT operator, different test times could lead to quite different thermal conductivities.Figure 9 shows an analysis of thermal conductivity values with variable test duration between 36 and 96 h.When looking at the thick lines on both graphs, for CXC and CXA, it can be seen that choosing different TRT periods could lead to as much as a 20% difference in the final result.From a standpoint of modelling geothermal heat pumps with multiple boreholes, such discrepancies could have a significant impact on oversizing or undersizing the geoexchange system, and on technoeconomical benefit.
Since the total duration of TRT is not strictly defined, only the minimum time is required, as explained in Section 2; in practical field testing, this always leads to a certain degree of analysis error.Since the final duration is arbitrarily chosen by the TRT operator, different test times could lead to quite different thermal conductivities.Figure 9 shows an analysis of thermal conductivity values with variable test duration between 36 and 96 h.When looking at the thick lines on both graphs, for CXC and CXA, it can be seen that choosing different TRT periods could lead to as much as a 20% difference in the final result.From a standpoint of modelling geothermal heat pumps with multiple boreholes, such discrepancies could have a significant impact on oversizing or undersizing the geoexchange system, and on technoeconomical benefit.Such a difference in results is influenced primarily by two factors: fluctuations of voltage in the public electrical grid, which is the usual method of powering TRT; and interferences between the surface equipment and air temperature, especially during the winter months.As explained in Section 2, fluctuations in voltage from the public grid are regularly seen to be following recognizable patterns in a 24-h cycle, depending on the specific demand during the day and night.Therefore, to at least try to minimize the effects of this on the development of line slope, as shown in Figure 6b, TRT duration should follow a 24 h multiplication factor (i.e., 48, 72, 96 h or 36, 60, 84 h).The other major concern is that of ambient temperature interfering with the TRT equipment and header pipes.As seen from Figure 5e, we used a header pipe length of 3 m from the borehole shaft to the TRT equipment in order to purposely magnify this effect, although the entire setup was properly insulated with 12 mm of caoutchouc rubber insulation.The effect of this could be seen in terms of the higher conductivity deviation for the CXC, as opposed to the CXA flow arrangement.Theoretically, flow direction arrangement does not have an impact on the conductivity measurement procedure, but rather on borehole resistance; nevertheless, the CXC tests were carried out during autumn/winter months, and the CXA during the spring months.The data shown in Table 2 conclusively suggest this claim, explaining the reason for the higher measured thermal conductivity for the CXC compared to the CXA tests.
Prolonging the TRT procedure by conducting additional recovery period, and using Horner's method to interpret the data, could lead to higher accuracy in approaching the actual ground thermal conductivity.This is clearly seen from Figure 9, where thermal conductivity obtained by Horner's method from the recovery period gives more symmetrical results, for both the CXC and CXA setups.This is explained by the fact that the recorded data is smooth, since there is no heat power applied (coefficient of determination from 0.996 to 1.0, as seen from Figure 8), and the fact that during the winter months, the interference heat flux between ambient and TRT equipment is lower, due to the Such a difference in results is influenced primarily by two factors: fluctuations of voltage in the public electrical grid, which is the usual method of powering TRT; and interferences between the surface equipment and air temperature, especially during the winter months.As explained in Section 2, fluctuations in voltage from the public grid are regularly seen to be following recognizable patterns in a 24-h cycle, depending on the specific demand during the day and night.Therefore, to at least try to minimize the effects of this on the development of line slope, as shown in Figure 6b, TRT duration should follow a 24 h multiplication factor (i.e., 48, 72, 96 h or 36, 60, 84 h).The other major concern is that of ambient temperature interfering with the TRT equipment and header pipes.As seen from Figure 5e, we used a header pipe length of 3 m from the borehole shaft to the TRT equipment in order to purposely magnify this effect, although the entire setup was properly insulated with 12 mm of caoutchouc rubber insulation.The effect of this could be seen in terms of the higher conductivity deviation for the CXC, as opposed to the CXA flow arrangement.Theoretically, flow direction arrangement does not have an impact on the conductivity measurement procedure, but rather on borehole resistance; nevertheless, the CXC tests were carried out during autumn/winter months, and the CXA during the spring months.The data shown in Table 2 conclusively suggest this claim, explaining the reason for the higher measured thermal conductivity for the CXC compared to the CXA tests.
Prolonging the TRT procedure by conducting additional recovery period, and using Horner's method to interpret the data, could lead to higher accuracy in approaching the actual ground thermal conductivity.This is clearly seen from Figure 9, where thermal conductivity obtained by Horner's method from the recovery period gives more symmetrical results, for both the CXC and CXA setups.This is explained by the fact that the recorded data is smooth, since there is no heat power applied (coefficient of determination from 0.996 to 1.0, as seen from Figure 8), and the fact that during the winter months, the interference heat flux between ambient and TRT equipment is lower, due to the lower temperature difference between the air temperature and the fluid temperature.The entire analysis presented with the two-step TRT principle could provide geothermal engineers a more precise method for qualitatively interpreting field data.
Furthermore, as stated in Section 3.4, both the TRT period and the recovery period could be used to determine the equivalent borehole resistance or skin effect, an equally important parameter for the efficient design of geoexchange systems.By applying Equations ( 24) and (25) to TRT data for each heat step, and Equation ( 27) for Horner's method on the recovery period, the following values were obtained, as shown in Table 3.It can be seen that the CXA flow arrangement generally shows a higher skin effect or equivalent borehole heat resistance (approximately +5%), which is in line with our previous research in this field [23], in which the borehole resistances of inclined coaxial and vertical 2U-loop heat exchangers were compared.Figure 10 shows the analysis of the ground temperature change at a certain radius away from the borehole after a Thermal Response Test time of 72 h.Temperature change was calculated with the infinite line source solution presented in a Section 3.1 for the case of five heat rejection steps.Due to its having more equable values, the ground effective thermal conductivity was set as the average of all ten measurements conducted by the Horner method in the recovery test period (Table 3).As mentioned in Section 4, densely drilled and installed inclined coaxial heat exchangers from a single shaft are prone to thermal interference between adjacent pipes for the first few meters of depth.As seen in Figure 10, for a TR test time of 72 h, these interferences could be ignored, as temperature change is negligible for a radius larger than 1.0 m.As shown in Figure 5, the two inclined boreholes were closer than this value only very near the surface, since drilling was conducted at a 45 • angle and in opposite directions.This result also suggests that the infinite line source solution is reasonable to use for two boreholes connected in series.
Discrepancies between the thermal conductivity obtained for the TRT period and the recovery period were rather high (TRT analysis showed 10-20% higher values than the Horner method), as can be seen from Table 3 and Figure 9.We have already explained the causes for this phenomenon, but using simple statistical analysis, such as sum of squares of difference, could provide knowledge of exactly which thermal conductivity coefficients are of statistical significance.The Sum of Squares or simply variation (SUMXMY2 function in MS Excel = Σ(x − y) 2 ) is a statistical technique used in regression analysis to determine the dispersion of data points.In a regression analysis, the goal is to determine how well a data series (in this case, measured EST) can be fitted to a function that might help to explain how the data series was generated (in this case, ILS with Ei function with two different thermal conductivities).The sum of squares is used as a mathematical way of finding the function that best fits (varies least) from the measured data.
The procedure was carried out for every heat step separately, with two obtained values for the thermal conductivity factor: one from the TRT period and one from the recovery period.Then, the entire curve was fitted with Equation ( 28) for the first step period, and Equation (29) for the second step period.An example is shown in Figure 11 for the case of CXC and 71 W/m.The dark blue dotted line is an ILS exponential integral function curve that takes into account a thermal conductivity value of 2.62 W/m • C (obtained during TRT), while the light blue dotted line is for a thermal conductivity value of 2.25 W/m • C (obtained during the recovery period).The purpose of this procedure was to find which thermal conductivity value better described the recorded temperature during the test for both periods (TRT 96 h + Recovery 96 h).In this way, a total of 20 points were obtained on the basis of the same principle as that presented in Figure 11; each data set was fitted with two Ei function curves with two different thermal conductivities.The results are presented in Figure 12a.
The same kind of variation analysis was performed with the EST hourly measured data and hourly air temperatures for each period (TRT or recovery) and heat step.The assumption was that the measurements with higher differences between the air temperature and the circulating fluid temperature would have higher thermal interferences and, consequently, higher degrees of error.This is especially true for this kind of TRT setup, where equipment is placed 3 m from the borehole (Figure 5e).Some thermal interference phenomena are to be expected in such cases, although the surface pipes are insulated.The results are shown in Figure 12b, in which the associated thermal conductivities are plotted against the sum of squares of differences between EST and Tair.Since the air temperature was recorded for the entire ~2000 h of measurement, the hypothesis is that climate interference would be zero if the air temperature were the same value as the current EST of the borehole fluid.The higher the temperature difference between the air temperature and the EST of the borehole fluid at any moment, the higher the interference would be on the TRT ground conductivity analysis, since a certain heat transfer would take place between the environment and the surface collector pipes and equipment.In this way, a total of 20 points were obtained on the basis of the same principle as that presented in Figure 11; each data set was fitted with two Ei function curves with two different thermal conductivities.The results are presented in Figure 12a.Both statistical analyses in Figure 12a,b shows similar trendlines, where for a value of SUMXMY2 = 0, the ground thermal conductivity coefficient is in the range 2.0-2.2W/m K.This is, in fact, a value very near the average of the ground thermal conductivities obtained by ten tests conducted in the recovery period (lhorner avg = 2.13 W/m °C), as seen in Figure 9.
As previously stated, the soil thermal diffusivity was reasonably estimated based on the known geological column and the obtained drilling samples.The soil thermal diffusivity termaffects the value of the calculated borehole resistance or skin term, due to the interconnectivity of the two variables in solutions presented for the diffusivity equations for cases of infinite medium and line source well.The conducted sum of squares of differences analysis, as presented in Figures 11 and 12, would not be affected by altering the initially assumed soil thermal diffusivity term.Changing this value, in a range of 0.03-0.06m 2 /d as expected for clays, would only change the resistivity/skin term.Inverse analysis could also be conducted, by calculating the thermal resistance term based on a known well geometry and material thermodynamic characteristics, and determining the soil thermal diffusivity term using ILS equations.However, practical geometry and the completion of the  The same kind of variation analysis was performed with the EST hourly measured data and hourly air temperatures for each period (TRT or recovery) and heat step.The assumption was that the measurements with higher differences between the air temperature and the circulating fluid temperature would have higher thermal interferences and, consequently, higher degrees of error.This is especially true for this kind of TRT setup, where equipment is placed 3 m from the borehole (Figure 5e).Some thermal interference phenomena are to be expected in such cases, although the surface pipes are insulated.The results are shown in Figure 12b, in which the associated thermal conductivities are plotted against the sum of squares of differences between EST and T air.Since the air temperature was recorded for the entire ~2000 h of measurement, the hypothesis is that climate interference would be zero if the air temperature were the same value as the current EST of the borehole fluid.The higher the temperature difference between the air temperature and the EST of the borehole fluid at any moment, the higher the interference would be on the TRT ground conductivity analysis, since a certain heat transfer would take place between the environment and the surface collector pipes and equipment.
Both statistical analyses in Figure 12a,b shows similar trendlines, where for a value of SUMXMY2 = 0, the ground thermal conductivity coefficient is in the range 2.0-2.2W/m K.This is, in fact, a value very near the average of the ground thermal conductivities obtained by ten tests conducted in the recovery period (l horner avg = 2.13 W/m • C), as seen in Figure 9.
As previously stated, the soil thermal diffusivity was reasonably estimated based on the known geological column and the obtained drilling samples.The soil thermal diffusivity termaffects the value of the calculated borehole resistance or skin term, due to the interconnectivity of the two variables in solutions presented for the diffusivity equations for cases of infinite medium and line source well.The conducted sum of squares of differences analysis, as presented in Figures 11 and 12, would not be affected by altering the initially assumed soil thermal diffusivity term.Changing this value, in a range of 0.03-0.06m 2 /d as expected for clays, would only change the resistivity/skin term.Inverse analysis could also be conducted, by calculating the thermal resistance term based on a known well geometry and material thermodynamic characteristics, and determining the soil thermal diffusivity term using ILS equations.However, practical geometry and the completion of the borehole heat exchanger is often different from the one assumed in this project.This problem is often related to misaligned pipes in the wellbore (especially in inclined coaxial systems) and the quality of the grouting (especially perfect adherence of grout and pipe).Soil diffusivity could be more precisely determined by laboratory measurements of undisturbed soil samples, by determining density, specific heat capacity and conductivity.
Therefore, the presented results confirm the hypothesis of this paper that a prolonged TRT should be conducted wherever possible, since it provides a higher certainty of the real ground thermal conductivity value.

Conclusions
Ground thermal conductivity and borehole skin factor, obtained from TRT measurements, are a prerequisite for optimized design of the borehole heat exchanger field.Such systems are, in practice, often unnecessarily oversized, as a result of the insecurity factor, or undersized, as a result of poor engineering design and lack of understanding of ground properties.In both cases, this negatively reflects on the economic viability of the geothermal project itself.From this perspective, longer and more advanced ground thermal response measurements are more cost-effective, as they ensure the longevity of the system and the knowledge of the borehole temperature evolution throughout the entire year.
Therefore, practical borehole heat exchanger field tests should always consist of a classic TRT test and a corresponding recovery test of the same duration.Interpreting the temperature recovery period with the Horner technique, derived from petroleum well testing experiences, elevates certainty that the obtained ground thermal conductivity coefficient is as close as possible to the real value.Furthermore, such a procedure minimizes the negative effects of electric grid voltage fluctuations and interferences from the ambient climate when performing classic TR tests.

Figure 1 .
Figure 1.Typical plot of pressure response during well testing and application of Horner's method.

Figure 1 .
Figure 1.Typical plot of pressure response during well testing and application of Horner's method.

Figure 2 .
Figure 2. Display of corresponding behavior of pressure/temperature during different types of well tests/thermal response tests.

Figure 2 .
Figure 2. Display of corresponding behavior of pressure/temperature during different types of well tests/thermal response tests.

Figure 3 .
Figure 3. Determination of semi-steady-state heat flow period via the derivation curve technique.

Figure 3 .
Figure 3. Determination of semi-steady-state heat flow period via the derivation curve technique.

Figure 4 .
Figure 4. Components of borehole thermal resistance for (a) the CXA arrangement; and (b) the CXC arrangement.

Figure 4 .
Figure 4. Components of borehole thermal resistance for (a) the CXA arrangement; and (b) the CXC arrangement.

Energies 2018 , 23 Figure 6 .
Figure 6.Example of two-step TRT analysis carried out for each heat power condition (CXC 71 W/m presented here); (a) Recorded borehole and air temperatures and unit heat power; (b) Determination of conductivity from TRT period; (c) Determination of conductivity by Horner's method from recovery period.

Figure 6 .
Figure 6.Example of two-step TRT analysis carried out for each heat power condition (CXC 71 W/m presented here); (a) Recorded borehole and air temperatures and unit heat power; (b) Determination of conductivity from TRT period; (c) Determination of conductivity by Horner's method from recovery period.

Figure 7 .
Figure 7. EST derivation curves for CXA and CXC flow arrangement and each TRT heat step.

Figure 8 .
Figure 8.Horner's semi-log method for determining ground thermal conductivity applied to a TRT recovery period for each heat power step and flow regime.

Figure 7 .
Figure 7. EST derivation curves for CXA and CXC flow arrangement and each TRT heat step.

Energies 2018 , 23 Figure 7 .
Figure 7. EST derivation curves for CXA and CXC flow arrangement and each TRT heat step.

Figure 8 .
Figure 8.Horner's semi-log method for determining ground thermal conductivity applied to a TRT recovery period for each heat power step and flow regime.

Figure 8 .
Figure 8.Horner's semi-log method for determining ground thermal conductivity applied to a TRT recovery period for each heat power step and flow regime.

Figure 9 .
Figure 9. Results of the ground thermal conductivity values obtained for classic TRT and recovery period as a function of variable test time 36-96 h.

Energies 2018 ,
11,  x FOR PEER REVIEW 19 of 23 step period.An example is shown in Figure11for the case of CXC and 71 W/m.The dark blue dotted line is an ILS exponential integral function curve that takes into account a thermal conductivity value of 2.62 W/m °C (obtained during TRT), while the light blue dotted line is for a thermal conductivity value of 2.25 W/m °C (obtained during the recovery period).The purpose of this procedure was to find which thermal conductivity value better described the recorded temperature during the test for both periods (TRT 96 h + Recovery 96 h).

Figure 10 .
Figure 10.Ground temperature rise at a certain radius from the borehole for five different TRT heat steps.

Figure 10 . 23 Figure 11 .
Figure 10.Ground temperature rise at a certain radius from the borehole for five different TRT heat steps.Energies 2018, 11, x FOR PEER REVIEW 20 of 23 of: CXC -71 W/m Measured EST ILS Ei fitted curve with conductivity derived from TRT period (

Figure 11 .
Figure 11.Example of fitting measured EST data with the infinite line source equation and two different values of thermal conductivity (the first obtained from the TRT period, and the second from the recovery period) for a case of 71 W/m and CXC setup.

Figure 12 .
Figure 12. Results of sum of squares of differences analysis: (a) SUMXMY2 between ILS fitted curve and EST measured data; (b) SUMXMY2 between air temperature and EST measured data.

Figure 12 .
Figure 12. Results of sum of squares of differences analysis: (a) SUMXMY2 between ILS fitted curve and EST measured data; (b) SUMXMY2 between air temperature and EST measured data.

Table 1 .
TRT heat step descriptive statistics.

Table 2 .
Air temperature descriptive statistics during TRT.

Table 3 .
Obtained data from analysis of TRT and recovery periods.

Table 3
Figure 11.Example of fitting measured EST data with the infinite line source equation and two different values of thermal conductivity (the first obtained from the TRT period, and the second from the recovery period) for a case of 71 W/m and CXC setup.