Closed-Loop Systems for Geothermal Energy Exploitation from Hydrocarbon Wells: An Italian Case Study

: Geothermal energy resources associated with disused hydrocarbon wells in Italian oilﬁelds represent a considerable source of renewable energy. Using the information available on Italian hydrocarbon wells and on-ﬁeld temperatures, two simpliﬁed closed-loop-type systems models were implemented in the Python environment and applied to a selected hydrocarbon well (Trecate4) located inside the Italian Villafortuna–Trecate ﬁeld (Northwestern Italy). Considering the maximum extracted working ﬂuid temperatures, Coaxial WBHE turned out to be a better performing technology than the U-tube version. The obtained outﬂow temperatures of the working ﬂuid at the wellhead for Coaxial and U-tube WBHEs of 98.6 ◦ C and 84 ◦ C, respectively, are both potentially exploitable for ensuring a multi-variant and comprehensive use of the resource through its application in sectors such as the food industry, horticultural and ﬂower ﬁelds.


Introduction
Energy production using affordable, sustainable and reliable resources has become one of the central topics of European development policy vision while concerns over climate change and energy independence have given increasing importance to the implementation of non-fossil fuels in the national energy mix. Advances in technological capabilities and substantial cost reductions in renewables are needed to facilitate the integration of these energy sources in end-use applications for the scaled transformation of existing energy systems.
Of the available environmentally friendly energy resources, geothermal energy can be used directly for heating and cooling buildings as well as to produce electricity. Moreover, depending on the specific geological, hydrogeological and thermo-physical conditions such as temperature, flow rate and geothermal water mineralisation grades, prevailing in an area to be analysed, there are several possibilities for the industrial and economic utilisation of heat energy potentially accumulated.
Direct use technologies are currently growing in applications including district heating, space heating with the use of heat pumps, agricultural purposes (e.g., greenhouses), fish farming, milk pasteurization and other purposes. By cascading exploitations of available heat, realising a multi-variant and comprehensive use of the resource is possible according to the corresponding temperature demand [1].
In [2], the authors presented a review of worldwide applications of geothermal energy for direct utilization, updating the previous survey carried out in 2015 and comparing data from the previous versions. The total installed capacity, reported at the end of 2019, for geothermal direct utilization worldwide turns out to be equal to 107,727 MWt, with a percentage increase of 52% over 2015. The total annual energy use is 1,020,887 TJ (283,580 GWh), indicating a 72.3% increase and a compound annual growth rate of about 11.5%. Over time, energy companies and researchers have begun to put greater effort into developing various strategies for harnessing this type of deep geothermal energy resource. The majority of studies that have been carried out on existing abandoned petroleum wells have focused on open-loop systems designed to repurpose petroleum fields as geothermal reservoirs [12,13].
Although geothermal open-loop systems in sedimentary formations may provide a sustainable solution where the geothermal potential and heat demand coincide [14,15], these technologies are subject to some technical problems, including groundwater recession, corrosion and scaling. Further important issues are represented by the reinjection of fluids [16,17] and by the complex additional exploration activities [18].
An effective and optimally dimensionalized alternative for relatively low energy demand is represented by the use of closed-loop deep geothermal systems. Unlike conventional open-loop geothermal systems, heat carrier fluids in closed-loop systems circulate inside of a wellbore heat exchanger (WBHE), while no ground fluids are extracted from the surrounding rocks. Moreover, corrosion and scaling problems are limited. Due to their proven advantages, a large number of studies that deal with developing closed-loop system technologies have been carried out [19,20].
Despite the recent success of some theoretical oilfield geothermal closed-loop system experiments, several challenges remain related to harnessing geothermal resources at a larger scale in oilfields using closed-loop type technologies, including low levels of thermal energy recovery and low energy conversion efficiency [21]. Moreover, researchers have indeed primarily focused their attention on analysing the impacts on the energy performance of changes in working fluid-related parameters, such as initial temperature and injection flow rate values.
Appl. Sci. 2021, 11, 10551 3 of 18 In the elaborations reported in [22][23][24], the authors fixed the values of different thermophysical parameters (thermal conductivity, volumetric heat capacity and rock density) as weighted mean values. However, due to the continuous spatial and vertical variability of geological formations associated with deep wells in oilfields, the thermophysical parameters of geological strata, as well as the depth of the strata and their thickness, need to be properly considered in order to achieve increasingly accurate and realistic estimates of heat exchanger performance.
Given the above considerations, the proposed study aims primarily at contributing to the discussion by encouraging reflection on the potential benefits of using low-temperature to medium-temperature geothermal energy resources associated with dismissed Italian hydrocarbon wells. Using the information about Italian hydrocarbon wells and on-field temperatures available from both the National Mining Office of the Italian Ministry for Economic Development (MISE) and the Italian National Geothermal Database, the attention is focused on two different geothermal closed-loop-type technologies: U-tube and coaxial wellbore heat exchangers (WBHEs). Two different simplified heat exchange models (Coaxial and U-tube WBHE) and thermal resistance systems, described in references [24,25], respectively, were implemented in the Python environment and applied to a selected hydrocarbon well (Trecate4) located inside the Italian Villafortuna-Trecate field (Northwestern Italy).
The main purposes of the study were to analyse heat exchange mechanisms by examining the main factors that affect extraction efficiency from a geological perspective, proposing a simplified and automated tool for the identification of the best closed-loop configuration that allows maximal heat recovery from a selected hydrocarbon well.
The assumption that the thermophysical parameters (thermal conductivity, volumetric heat capacity and rock density) are constant values has been overcome by considering and implementing detailed stratigraphic data in the proposed models. The final use of the potentially accumulated heat energy was considered as possible in direct applications through a cascade plant system, which provides specific thermal energy amounts to different production cycles in manufacturing, agricultural and recreational districts near the oilfield.

Test Site: The Villafortuna-Trecate Oilfield
A variety of petroleum systems has been identified in Italy as a result of complex geological and sedimentary history. According to [26], Italian hydrocarbon occurrences can be classified as associated with three main tectono-stratigraphic systems: the carbonate Mesozoic substratum of the foredeep/foreland area and the external thrust belts, thrusted terrigenous Oligo-Miocene foredeep wedges (Southern Alps, Northern Apennines, Calabria and Sicily) and terrigenous Pliocene-Pleistocene successions of the late foredeep basins of the Apennines, in both the Central and Northern Adriatic Sea and the Po Plain.
In this study, attention is focused on an analysis of the geothermal resources potentially associated with a disused hydrocarbon well named Trecate4, which is located within the Villafortuna-Trecate oilfield ( Figure 1).
The Villafortuna-Trecate system represents one of the largest oil accumulations of the Italian Middle Triassic carbonate petroleum system. It involves dolomitized platform units of Late Triassic-Early Jurassic that were charged by Middle Triassic carbonate source rocks deposited in the confined basins created by rifting. Due to its depth, it can be pursued only in the outer sector of the foredeep and in the foreland regions (the Piedmont area) [27]. The main reservoir associated with the Villafortuna-Trecate field was identified at a depth between 5800 m and 6100 m [28].
The selection of the proposed case study was guided by the preliminary results presented in [29], which identified the wells located within the Villafortuna-Trecate oilfield as the most suitable to be converted into geothermal ones and also in function of their deep and constructional features. Considering the related thermophysical parameters, the amount of thermal energy that they would allow recovering turned out to be larger in comparison with the wells located in comparable Italian sedimentary contexts.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 4 of 18 area) [27]. The main reservoir associated with the Villafortuna-Trecate field was identified at a depth between 5800 m and 6100 m [28]. The selection of the proposed case study was guided by the preliminary results presented in [29], which identified the wells located within the Villafortuna-Trecate oilfield as the most suitable to be converted into geothermal ones and also in function of their deep and constructional features. Considering the related thermophysical parameters, the amount of thermal energy that they would allow recovering turned out to be larger in comparison with the wells located in comparable Italian sedimentary contexts. Detailed litho-stratigraphic information and temperature data visualisation related to the Trecate4 hydrocarbon well have been found on the Italian National Geothermal Database (BDNG), the largest collection of Italian geothermal data [8]. Moreover, information regarding productive and dismissed oil and gas wells in Italy was also provided by the National Mining Office of the Italian Ministry for Economic Development (MISE) and by the website VIGOR promoted by MISE-DGRME (Direzione Generale Risorse Minerarie ed Energetiche) and the Italian Geological Society and the Assomineraria Association.
For properly defining thermophysical parameters, values reported in [31,32] were used. Notably, in [31], by measuring its value in both dry and wet conditions, the authors investigated and defined the thermal conductivity of 200 rock samples collected from four different regions of Southern Italy (Calabria, Campania, Apulia and Sicily). Moreover, in [32], the framework of the MIUR 2008 project "Geothermal resources of the Mesozoic basement of the Po Basin: groundwater flow and heat transport" was utilized in order to accurately estimate the thermophysical properties of a wide variety of sedimentary and intrasedimentary volcanic rocks from the Po Basin.
As demonstrated by information related to the litho-stratigraphic units reported in Table 1, the stratigraphic succession associated with the Trecate4 well is mainly composed of carbonate platform rocks, with the exception of a surface unit consisting of terrigenous deposits. The analysed well has a maximum depth of 6282 m and bottom-hole temperatures reaching 168 °C ( Figure 2). Detailed litho-stratigraphic information and temperature data visualisation related to the Trecate4 hydrocarbon well have been found on the Italian National Geothermal Database (BDNG), the largest collection of Italian geothermal data [8]. Moreover, information regarding productive and dismissed oil and gas wells in Italy was also provided by the National Mining Office of the Italian Ministry for Economic Development (MISE) and by the website VIGOR promoted by MISE-DGRME (Direzione Generale Risorse Minerarie ed Energetiche) and the Italian Geological Society and the Assomineraria Association.
For properly defining thermophysical parameters, values reported in [31,32] were used. Notably, in [31], by measuring its value in both dry and wet conditions, the authors investigated and defined the thermal conductivity of 200 rock samples collected from four different regions of Southern Italy (Calabria, Campania, Apulia and Sicily). Moreover, in [32], the framework of the MIUR 2008 project "Geothermal resources of the Mesozoic basement of the Po Basin: groundwater flow and heat transport" was utilized in order to accurately estimate the thermophysical properties of a wide variety of sedimentary and intrasedimentary volcanic rocks from the Po Basin.
As demonstrated by information related to the litho-stratigraphic units reported in Table 1, the stratigraphic succession associated with the Trecate4 well is mainly composed of carbonate platform rocks, with the exception of a surface unit consisting of terrigenous deposits. The analysed well has a maximum depth of 6282 m and bottom-hole temperatures reaching 168 • C ( Figure 2).

Closed-Loop Geothermal Energy Systems: WBHEs
In current practice, two main types of closed-loop systems have been described fo harnessing geothermal energy resources by taking advantage of disused boreholes in oi fields: U-tube and coaxial double-pipe WBHE technologies [33,34].
In U-tube heat exchangers, fluid is pumped through one tube string and comes ou of the other. It is by this action of flowing through the well that the fluid in the U-tube ca gain heat energy from the surrounding geological formations. The gaps between intern pipes are filled with an insulating material, and the bottom hole is sealed. On the othe hand, the coaxial heat exchanger comprises two concentric pipes. Circulating workin fluid is injected into the outer pipe (the injection pipe), flows down to the lower part o the exchanger and is gradually warmed by heat from the rocks. After the fluid reaches th bottom hole of the well, it flows upwards through a thinner pipe, which acts as the inne pipe (the extraction pipe). Both pipes are thermally insulated, and the bottom hole sealed. Heat exchange occurs between the geological formation and the fluid in the inje tion pipe and between the fluid in the injection pipe and the fluid in the extraction pip [10].
For both closed-loop geothermal systems, the consequences of corrosion have to b considered in the selection phase of the working fluid. Due to its low cost, heat transfe and stored capacity, water represents one of the most used working fluids. Moreove

Closed-Loop Geothermal Energy Systems: WBHEs
In current practice, two main types of closed-loop systems have been described for harnessing geothermal energy resources by taking advantage of disused boreholes in oilfields: U-tube and coaxial double-pipe WBHE technologies [33,34].
In U-tube heat exchangers, fluid is pumped through one tube string and comes out of the other. It is by this action of flowing through the well that the fluid in the U-tube can gain heat energy from the surrounding geological formations. The gaps between internal pipes are filled with an insulating material, and the bottom hole is sealed. On the other hand, the coaxial heat exchanger comprises two concentric pipes. Circulating working fluid is injected into the outer pipe (the injection pipe), flows down to the lower part of the exchanger and is gradually warmed by heat from the rocks. After the fluid reaches the bottom hole of the well, it flows upwards through a thinner pipe, which acts as the inner pipe (the extraction pipe). Both pipes are thermally insulated, and the bottom hole is sealed. Heat exchange occurs between the geological formation and the fluid in the injection pipe and between the fluid in the injection pipe and the fluid in the extraction pipe [10].
For both closed-loop geothermal systems, the consequences of corrosion have to be considered in the selection phase of the working fluid. Due to its low cost, heat transfer and stored capacity, water represents one of the most used working fluids. Moreover, operating parameters such as the fluid flow rate and pipe diameter have to be selected in order to guarantee transient turbulent flow conditions in the subsurface, thus facilitating heat transfer from the ground to the fluid [33].
For properly evaluating the temperature profiles of the selected fluid (water) associated with both U-tube and coaxial double-pipe WBHE technologies, the models proposed below were implemented in a Python environment, identifying the best geothermal system configuration. Then, they were applied to the selected case study represented by the Trecate4 hydrocarbon well.

Heat Transfer in Coaxial WBHEs
In coaxial WBHEs, the steel downward is in contact with the hole in the well. The energy balance of the fluid in the injection pipe can be expressed with the following equation [14,35]: where A 0 and v f are the outer pipe area and fluid velocity, respectively, T fo is the fluid temperature in the outer pipe and dQ/dz is the heat extraction from the formation at unit well depth (W/m) ( Table 2). Although insulation is used to prevent heat loss from the inner-pipe fluid, heat is partly transferred between the two pipes; thus, dQ io /dz represents the heat flux from the inner pipe to the outer pipe. Therefore, the energy equation for the inner pipe can be given as follows.
By assuming steady heat transfer and constant heat flux in wellbore components (e.g., insulation, casing and cement), heat extraction from the formation dQ/dz can be assumed to be equal to the heat flux through the outside surface of the wellbore (interface of the wellbore and rock formation) to the injected fluid [14,35]: where T w is the temperature at the interface of the wellbore and the formation, k w is the heat transfer coefficient between the outer pipe fluid and wellbore exterior and R w is the resistance between the outer pipe and surrounding rocks. At the well bottom, the heated fluid is forced to enter and flow through the internal pipe of the coaxial WBHE. Proceeding upwards to the wellhead, heat transfer occurs only through the wall of the internal pipe. Thus, dQ i0 /dz is determined by considering the temperature difference between the outer-pipe and inner-pipe fluids, as well as the estimated thermal resistance of the insulation: where T fi is the fluid temperature in the inner pipe, k i0 is the heat transfer coefficient between the outer pipe and inner pipe and R i0 is the thermal resistance between the outer pipe and inner pipe.
2.3.1. Coaxial WBHE: Coefficient of Heat Exchange between Outer-Pipe Fluid and the Wellbore Exterior In an analysis of the energy balance equation for the fluid in the outer pipe (injection pipe) of a coaxial WBHE, a proper estimate of the parameter k w is fundamental for a correct evaluation of the heat exchange between the outer-pipe fluid and the drilled geological formations.
For a coaxial WBHE, the heat exchange coefficient for the injection pipe can be expressed as the sum of heat transfer components, expressed in terms of thermal resistance values (R w ) [24]: where R s is a function of time that represents thermal resistance due to conductive heat transfer in the rock, R a is the thermal resistance due to convective heat transfer in the pipe and R c is the thermal resistance due to conductive heat transfer through the casings of the well.
In the evaluation of total thermal resistance, the conductive term prevails; consequently, thermal exchange is directly proportional to the convective transfer coefficient. Conductive thermal resistance (R s ) can be expressed as follows: where λ s (W/mK) is the thermal conductivity of the rock, and α s (m/s) is the thermal diffusivity of the rock. In Equation (6), the numerator of the argument of the natural logarithm represents the time-dependent radius of the thermal influence of the well (r s ). This parameter considers the change, over time, of the heat flux into WBHE surroundings geological formations [24]. Convective thermal resistance (R a ) can be determined by the following equation: where r c is the external radius of the external casing, and h f is the convective heat transfer coefficient, which was calculated using the Nusselt number (Nu) and a form of the Dittus-Boelter equation that assumes turbulent flow inside the tubes (Reynolds number ≥ 104) [36]: with Pr = ρc f µ λ f and Re = ρv f 2r c µ . Finally, thermal resistance to heat conduction through the casings of the well can be determined as follows: where λ i is the thermal conductivity of the rock in correspondence with the different casings of the well. Generally, due to the high thermal conductivity of the steel piping, the total thermal resistance of the casing is negligible compared to the rock's thermal resistance. As a result, the heat exchange coefficient k w can be correctly determined as follows [22]: where r c = r w as the thickness of the external tube is negligible.

Coaxial WBHE: Coefficient of the Heat Exchange between the Outer-Pipe Fluid and the Inner Pipe
Unlike in the injection pipe, the total heat flux in the upward pipe (extraction pipe) is determined by a conductive component of the composite pipe and by two convective components: one on the internal wall and one on the external wall of the WBHE.
Consequently, the total heat exchange coefficient k i0 for the extraction pipe can be calculated as follows [22]: where r 0 is the radius of the inner pipe, d is the thicknesses of the pipe exchanger, h 0 and h i are the coefficients of convective heat transfer to the inner and outer wall, respectively, and λ i is the thermal conductivity of the pipe material (air and steel).

U-Tube WBHE: Thermal Resistances Model
In U-tube WBHE, the hole between the tubes and the well is filled with grout (bentonite) in order to prevent direct leakage between the ground and tubes, avoiding a connection between the ground and Earth's subsurface.
Several approaches have been proposed for reproducing the thermal behaviour of different U-tube WBHE configurations [37]. In this work, to compute the temperature profile in the U-tube configuration using a one-dimensional heat exchange, the model with a set of equivalent thermal resistances described in [25] was applied. Six thermal resistances were considered at each depth, including the thermal properties of the ground, the grout and the pipes (Figure 3).
Appl. Sci. 2021, 11, x FOR PEER REVIEW  The thermal resistances between the grout (bentonite) and pipe depend on the overall borehole thermal resistance R bhe . This parameter is the average thermal resistance between the fluid in the pipe and the borehole wall, and it is usually determined after experimental tests. The grout area is divided into two zones according to the pipe numbers; thus, R bhe was also considered as divided into two parallel resistances that connected each pipe with the corresponding grout zone. This parameter can be further divided into a convective (R h ) and a conductive (R c ) term (Equation (13)): where R c represents conductive thermal resistance, considering the total conductive resistance between the pipes and the borehole wall. As the grout node can be located at a certain distance D x , R c is divided into two different resistances (Equations (14) and (15)): where R b represents the conductive resistance between the pipe and the considered node. R h is calculated by following Equation (16), where r i is the internal pipe radius, h f is the convective heat coefficient and Nu is the dimensionless Nusselt number that is calculated using Equation (9).
The R c value is indeed estimated using a calculation method that requires the estimation of the pipe equivalent surface S eq and its diameter D eq (Equation (17)) ( Table 3).
The pipe equivalent surface (S eq ) can be calculated by following the approach proposed in [38], considering the sum of the areas S gg and S p as indicated in Figure 4. Consequently, D eq is obtained by following Equation (18).
Using the equivalent diameter, the conductive thermal resistances R b and R x are calculated by considering a semi-cylindrical conductive heat transfer: where k b is the borehole grout thermal conductivity. Usually, when the pipes are near the borehole wall, D x can be located at the same distance from the borehole wall; thus D x = D b , and, as a consequence, the R x part of the conductive resistance can be neglected.
Finally, the correlation as reported in Equation (6) is used for the estimation of ground thermal resistances. The total resistance value related to the analysed single tube includes convective (R a ), grout conduction (R b + R x ) and ground resistance (R s ) terms.    [22]. Copyright Elsevier, 2014).

WBHE Model Assumptions
Python is a high-level programming language with an object-oriented approach described in [39]. As the Python software allows algorithms to be solved easily, it was used to perform the analysis of the WBHEs by implementing the described models (coaxial and U-tube-WBHE models) in different codes that have the structure illustrated in Figure 5.
According to both selected models, the propagation of heat in the reservoir occurs through conduction; the propagation of heat inside the wellbore tubes takes place through conduction and convection phenomena. The reservoir model was built by assuming a single well located at the center of a circular reservoir.
The temperature profile in the radial direction was assumed to be constant. Therefore, there was no temperature gradient in the annulus or the inner tube. Due to the turbulent flow, enhanced mixing phenomena occurred, which decreased the radial gradient. The temperature changed only in the annulus and in the vertical direction of the inner tube. Consequently, the temperature profile was unidirectional (vertical temperature profile).
The properties of the heat carrier fluid were assumed to be constant. As the fluid used in this study was water (100 °C , 2 bar), no variations occurred due to pressure or/and temperature gradients.
The analytical models were built under steady-state conditions; there were no temperature variations over time, with each point in the tubes (annulus and inner tube) maintaining the same temperature for the lifecycle of the system. In addition, the models considered the resistance associated with tube thickness to be negligible. The tube material had very high conductivity (15 W/mK), so its resistance could be considered small compared to the other resistances in the system. For the estimation of resistance associated with the rock (see Equation (6)), the time value used was 3 years. The analysed system, in the period preceding the three years (1-3 years), turns out to ensure larger heat exchange phenomena with the possibility of causing overestimation in energy performances.

WBHE Model Assumptions
Python is a high-level programming language with an object-oriented approach described in [39]. As the Python software allows algorithms to be solved easily, it was used to perform the analysis of the WBHEs by implementing the described models (coaxial and U-tube-WBHE models) in different codes that have the structure illustrated in Figure 5. Considering the U-tube configuration, heat exchange was assumed to take place in an area that was half the area of the casing pipe, while the interaction between the downward and upward tubes was neglected.
Both proposed analytical models follow the path of the working fluid with an approach that could be called step-by-step. In detail, they considered intervals of length dz in which the inlet and outlet temperatures were calculated by solving the energy balance equation for each considered volume dv. For the estimation of the energy exchange in the radial direction, the mean value of the temperature in the volume dv was used and calculated by using the arithmetic mean. All energy exchanged in the radial direction in the volume dv was absorbed by water (selected working fluid).

Influence of Ground Thermal Properties on Heat Exchange Phenomena
During a preliminary phase of analysis to understand the extent of the variation in thermal outputs associated with the considered WBHE technologies, a sensitivity analysis was performed on the different underground thermal properties. Starting from the analysis of Equation 6, describing the thermal resistance due to conductive heat transfer in the rock, we proceeded to the graphic elaboration of curves of the variation of the resistance parameter over time with the values of density (ρ), volumetric heat capacity (ρcs) and the According to both selected models, the propagation of heat in the reservoir occurs through conduction; the propagation of heat inside the wellbore tubes takes place through conduction and convection phenomena. The reservoir model was built by assuming a single well located at the center of a circular reservoir.
The temperature profile in the radial direction was assumed to be constant. Therefore, there was no temperature gradient in the annulus or the inner tube. Due to the turbulent flow, enhanced mixing phenomena occurred, which decreased the radial gradient. The temperature changed only in the annulus and in the vertical direction of the inner tube. Consequently, the temperature profile was unidirectional (vertical temperature profile).
The properties of the heat carrier fluid were assumed to be constant. As the fluid used in this study was water (100 • C, 2 bar), no variations occurred due to pressure or/and temperature gradients.
The analytical models were built under steady-state conditions; there were no temperature variations over time, with each point in the tubes (annulus and inner tube) maintaining the same temperature for the lifecycle of the system. In addition, the models considered the resistance associated with tube thickness to be negligible. The tube material had very high conductivity (15 W/mK), so its resistance could be considered small compared to the other resistances in the system. For the estimation of resistance associated with the rock (see Equation (6)), the time value used was 3 years. The analysed system, in the period preceding the three years (1-3 years), turns out to ensure larger heat exchange phenomena with the possibility of causing overestimation in energy performances.
Considering the U-tube configuration, heat exchange was assumed to take place in an area that was half the area of the casing pipe, while the interaction between the downward and upward tubes was neglected.
Both proposed analytical models follow the path of the working fluid with an approach that could be called step-by-step. In detail, they considered intervals of length dz in which the inlet and outlet temperatures were calculated by solving the energy balance equation for each considered volume dv. For the estimation of the energy exchange in the radial direction, the mean value of the temperature in the volume dv was used and calculated by using the arithmetic mean. All energy exchanged in the radial direction in the volume dv was absorbed by water (selected working fluid).

Influence of Ground Thermal Properties on Heat Exchange Phenomena
During a preliminary phase of analysis to understand the extent of the variation in thermal outputs associated with the considered WBHE technologies, a sensitivity analysis was performed on the different underground thermal properties. Starting from the analysis of Equation 6, describing the thermal resistance due to conductive heat transfer in the rock, we proceeded to the graphic elaboration of curves of the variation of the resistance parameter over time with the values of density (ρ), volumetric heat capacity (ρc s ) and the thermal conductivity of rock (λ s ) change (Figure 6a-c). Moreover, an evaluation of the impact of thermal diffusivity value (α) variation in the R s parameter was reported in Figure 6d.
Different input values associated with thermal properties subject to the performed sensitivity analysis were chosen in order to properly cover a significant range for the types of geological formations involved (sedimentary rocks) [28,29]. The proposed analysis has been useful for understanding how thermal outputs change as a consequence of the variations in the input data of ground thermal properties. The ground thermal properties directly influence ground resistance values (R s ), as highlighted in Equation (6). The highest wellhead temperature variation is caused by changes in thermal conductivity values. The implementation in the proposed models of proper values of geological formation thermal properties (case studies-specific stratigraphy) turns out to be fundamental for improving the quality of the performed heat exchange analysis.

WBHE Configurations
The temperature profiles associated with the described Coaxial and U-tube WBHE system configurations were obtained by making use of the specific ground properties of the selected case study (Trecate4 hydrocarbon well) reported in Table 1. Values relating to the thermal properties of the different rock formations have been attributed by following those proposed in [31,32] (Table 1). Moreover, the thermal conductivity value of the insulating material was set to 0.025 W/mK. Working fluid inlet flow rate (water) and temperature values were considered as 3.0 kg/s and 50 °C , respectively. Selected inlet temperature represents a typical value for direct applications such as production cycles in manufacturing and agricultural districts [1]. Subsequently, an analysis was conducted on the term vector fluid temperature at the outlet as the inlet flow rate varies.
For the coaxial WBHE configuration, the sizings of the inner and outer tubes, as well as the final casing size, were fixed according to the values proposed in Table 4a. The characteristic configuration of the U-Tube WBHE is shown in Table 4b.

WBHE Configurations
The temperature profiles associated with the described Coaxial and U-tube WBHE system configurations were obtained by making use of the specific ground properties of the selected case study (Trecate4 hydrocarbon well) reported in Table 1. Values relating to the thermal properties of the different rock formations have been attributed by following those proposed in [31,32] (Table 1). Moreover, the thermal conductivity value of the insulating material was set to 0.025 W/mK. Working fluid inlet flow rate (water) and temperature values were considered as 3.0 kg/s and 50 • C, respectively. Selected inlet temperature represents a typical value for direct applications such as production cycles in manufacturing and agricultural districts [1]. Subsequently, an analysis was conducted on the term vector fluid temperature at the outlet as the inlet flow rate varies.
For the coaxial WBHE configuration, the sizings of the inner and outer tubes, as well as the final casing size, were fixed according to the values proposed in Table 4a. The characteristic configuration of the U-Tube WBHE is shown in Table 4b.
For the Coaxial WBHE configuration, over the first 1200 m downwards, the selected thermo-vector fluid decreases in temperature. As the downward profile line crosses the underground temperature line, the real heating process begins, and the ground contribution became positive. Due to the presence of the insulating material, the heat exchange coefficient between tubes turns out to be low, and the increase in working fluid temperature can be associated with the ground (Figure 7a). underground temperature line, the real heating process begins, and the ground contribution became positive. Due to the presence of the insulating material, the heat exchange coefficient between tubes turns out to be low, and the increase in working fluid temperature can be associated with the ground (Figure 7a). In coaxial WBHE, the maximum recorded temperature is equal to 98.6 °C . The inlet flow rate strongly influences the temperature of the wellhead thermal fluid and, consequently, the heat power amount. As observed in Figure 8a, by considering inlet flow rate values between 0.5 and 1 kg/s, the output fluid temperature increases up to about 155 °C. Consequently, for the specific case study and associated plant configuration, it is possible to identify an inlet flow rate value that potentially allows an optimization of the wellhead temperature.
Considering the U-tube WBHE configuration and its associated temperature profile, it was possible to identify how the ground's contribution was responsible for a sizable temperature variation, both in the downward and upward tube. Over the first borehole section (approximately 1200 m), the fluid was cooled by the ground. Subsequently, as the temperature profile line crossed the underground temperature line, the trend was inverted, and the heat carrier fluid temperature began to increase, reaching a maximum temperature value of 84 °C (Figure 7b).
As for the considered coaxial-WBHE configuration, the wellhead thermal fluid temperature in U-tube WBHE changes as the input flow rate parameter varies; these were estimated (Figure 8b).  In coaxial WBHE, the maximum recorded temperature is equal to 98.6 • C. The inlet flow rate strongly influences the temperature of the wellhead thermal fluid and, consequently, the heat power amount. As observed in Figure 8a, by considering inlet flow rate values between 0.5 and 1 kg/s, the output fluid temperature increases up to about 155 • C. Consequently, for the specific case study and associated plant configuration, it is possible to identify an inlet flow rate value that potentially allows an optimization of the wellhead temperature.

Conclusions
Geothermal energy resources associated with disused hydrocarbon wells in Italian oilfields represent a considerable source of renewable energy. As a type of energy stored in subsurface geological formations and associated with hydrocarbons in oilfields, geothermal energy needs to be extracted before final utilisations. Closed-loop technologies represent an effective alternative in contrast to conventional open-loop geothermal systems as heat carrier fluids circulate inside wellbore heat exchangers (WBHEs), while no ground fluids are extracted from surrounding rocks.
In this paper, with the main aim to propose a simplified tool for the definition of the best closed-loop system that allows maximizing heat recovered from a selected hydrocarbon well, two different simplified heat exchange models (U-tube and coaxial WBHE) have been described and implemented in a Python programming language. Both the site-specific geological and thermo-physical properties of a selected study case (Trecate4 hydrocarbon well) were properly considered for the computer analyses carried out.
The implementation of site-specific stratigraphy in the performed models was fundamental as it contributes to improving the amount of heat exchange phenomena. Moreover, by analysing the results obtained as output from the developed models (U-tube and Considering the U-tube WBHE configuration and its associated temperature profile, it was possible to identify how the ground's contribution was responsible for a sizable temperature variation, both in the downward and upward tube. Over the first borehole section (approximately 1200 m), the fluid was cooled by the ground. Subsequently, as the temperature profile line crossed the underground temperature line, the trend was inverted, and the heat carrier fluid temperature began to increase, reaching a maximum temperature value of 84 • C ( Figure 7b).
As for the considered coaxial-WBHE configuration, the wellhead thermal fluid temperature in U-tube WBHE changes as the input flow rate parameter varies; these were estimated (Figure 8b).
For a flow rate value of 1.5 kg/s, the fluid reaches the surface at a maximum temperature of 89.4 • C. For higher flow rate values, the wellhead temperatures recorded are progressively lower. Using the fixed inlet working fluid temperature (50 • C) and the estimated maximum fluid temperature at the outlet for the different configurations, thermal power values were evaluated for 616.7 KW (98.6 • C-Coaxial WBHE) and 427.9 KW (84 • C-U-tube WBHE) (Figure 9a,b).
Considering a cascading exploitation mode of the heat accumulated by the working fluid in Trecate4 WBHEs, it is possible to hypothesize a multi-variant and comprehensive use of the resource. The outflow temperatures of working fluids at the wellhead for Coaxial and U-tube WBHEs are 98.6 • C and 84 • C, respectively, which allows it to progressively be used for greenhouse heating (100-80 • C), food industry (80-70 • C), animal breeding (60 • C), biomass and agricultural cultures (<50 • C).
The utilization of the geothermal resources at cascade levels is indeed identified as a series of sequential operations by integrating different technologies for distribution and use of thermal energy, drying and dehydration processes, recreational uses and any other direct use [40].

Conclusions
Geothermal energy resources associated with disused hydrocarbon wells in Italian oilfields represent a considerable source of renewable energy. As a type of energy stored in subsurface geological formations and associated with hydrocarbons in oilfields, geothermal energy needs to be extracted before final utilisations. Closed-loop technologies represent an effective alternative in contrast to conventional open-loop geothermal systems as heat carrier fluids circulate inside wellbore heat exchangers (WBHEs), while no ground fluids are extracted from surrounding rocks.
In this paper, with the main aim to propose a simplified tool for the definition of the best closed-loop system that allows maximizing heat recovered from a selected hydrocarbon well, two different simplified heat exchange models (U-tube and coaxial WBHE) have been described and implemented in a Python programming language. Both the site-specific geological and thermo-physical properties of a selected study case (Trecate4 hydrocarbon well) were properly considered for the computer analyses carried out.
The implementation of site-specific stratigraphy in the performed models was fundamental as it contributes to improving the amount of heat exchange phenomena. Moreover, by analysing the results obtained as output from the developed models (U-tube and Coaxial), the coaxial WBHE technology turned out to be better performing. In detail, from the comparison between the temperature profiles obtained by fixing an inlet mass flow Geothermal energy derived from a converted geothermal well can contribute to reducing the consumption of fuel used for heating and is strongly dependent on the production techniques adopted, the species and varieties cultivated and the control strategies of the climatic conditions adopted.

Conclusions
Geothermal energy resources associated with disused hydrocarbon wells in Italian oilfields represent a considerable source of renewable energy. As a type of energy stored in subsurface geological formations and associated with hydrocarbons in oilfields, geothermal energy needs to be extracted before final utilisations. Closed-loop technologies represent an effective alternative in contrast to conventional open-loop geothermal systems as heat carrier fluids circulate inside wellbore heat exchangers (WBHEs), while no ground fluids are extracted from surrounding rocks.
In this paper, with the main aim to propose a simplified tool for the definition of the best closed-loop system that allows maximizing heat recovered from a selected hydrocarbon well, two different simplified heat exchange models (U-tube and coaxial WBHE) have been described and implemented in a Python programming language. Both the site-specific geological and thermo-physical properties of a selected study case (Trecate4 hydrocarbon well) were properly considered for the computer analyses carried out.
The implementation of site-specific stratigraphy in the performed models was fundamental as it contributes to improving the amount of heat exchange phenomena. Moreover, by analysing the results obtained as output from the developed models (U-tube and Coaxial), the coaxial WBHE technology turned out to be better performing. In detail, from the comparison between the temperature profiles obtained by fixing an inlet mass flow value of 3 kg/s, it can be found that the outlet fluid temperature for the coaxial configuration (98.6 • C) is higher than that estimated for the U-tube version (84 • C). Even for variable flow rate values, an ever-higher output temperature for the coaxial configuration is recorded. Thermal powers are consequently higher. The outflow temperatures of geothermal water at the wellhead for both Coaxial and U-tube WBHEs could progressively be used for greenhouse heating (100-80 • C), food industry (80-70 • C), animal breeding (60 • C), biomass and agricultural cultures (<50 • C).
An analysis of the possibility of having phase change (evaporation) of the working fluid and the role in heat transfer and performance of extracting heat from abandoned wells of intraformational flows needs to be performed in order to further improve the accuracy of models. The proposed approach covering the described simplified models represents a useful methodological tool that allows for preliminary studies of the possibility of a selected Italian hydrocarbon well to be converted into a geothermal one by using WBHE technologies.