Reuse of Decommissioned Hydrocarbon Wells in Italian Oilﬁelds by Means of a Closed-Loop Geothermal System

: Geological and geophysical exploration campaigns have ascertained the coexistence of low to medium-temperature geothermal energy resources in the deepest regions of Italian sedimentary basins. As such, energy production based on the exploitation of available geothermal resources associated with disused deep oil and gas wells in Italian oilﬁelds could represent a considerable source of renewable energy. This study used information available on Italian hydrocarbon wells and on-ﬁeld temperatures to apply a simpliﬁed closed-loop coaxial Wellbore Heat Exchanger (WBHE) model to three different hydrocarbon wells located in different Italian oilﬁelds (Villafortuna-Trecate, Val d’Agri ﬁeld, Gela ﬁelds). From this study, the authors have highlighted the differences in the quantity of potentially extracted thermal energy from different analysed wells. Considering the maximum extracted working ﬂuid temperature of 100 ◦ C and imagining a cascading exploitation mode of the heat accumulated, for Villafortuna 1 WBHE was it possible to hypothesise a multi-variant and comprehensive use of the resource. This could be done using existing infrastructure, available technologies, and current knowledge.


Introduction
The policy visions of the 2030 UN Agenda for Sustainable Development and the Paris Agreement on Climate Change were both approved by the member states of the UN in 2015 and represent two fundamental contributions that guide the transition towards an economic model that aims not only for profitability but also for social progress and environmental protection. To achieve increased energy efficiency, all nations must change the ways in which they produce and manage natural energy resources in order to create more sustainable and environmentally resilient communities.
The current Italian urban energy paradigm relies heavily on fossil fuels; given the air pollution and resource depletion caused by fossil-fuel use, this is unsustainable. Public and private investments in energy must be increased, especially for innovative technological models that transform the energy system, integrating renewable energy into end-use applications in buildings, transportation and industry.
Geoscience offers solutions to this issue through the development of various options that could encourage decarbonisation and the transition to renewable energy sources at local and regional scales: electricity production from renewable sources, domestic heating/cooling using low-enthalpy geothermal energy resources and larger-scale technologies that target harmful emissions, such as bioenergy production and carbon capture and storage [1].
Geothermal energy is a weather-independent, environmentally friendly and currently available renewable resource; it represents an effective solution for power generation, heating and cooling, and direct-use applications. Specifically, energy production from temperature geothermal energy resources associated with dismissed Italian hydrocarbon wells as a renewable source of energy.
Using the information on 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 simplified heat exchange model (coaxial WBHE) described in [13] was applied to three hydrocarbon wells. The main purpose was to properly analyse heat exchange mechanisms in three different Italian oil and gas fields, emphasising that the quantity of the potentially extracted thermal energy can change based on geological and depositional context. Wells from the Villafortuna-Trecate field, the Val d'Agri field and the Gela field were selected.
The assumption that the thermophysical parameters (thermal conductivity, volumetric heat capacity and rock density) are constant values, as described in [13], has been overcome by applying the elaborated model to the detailed stratigraphic data of each case study.
The final use of the potentially extracted heat would be in possible direct applications by means of a cascade plant system, which provides specific thermal energy amounts to production cycles in manufacturing, agricultural and recreational districts near the oilfields. Considering specifically the examples of disused hydrocarbon wells analysed, located in the Italian territory, it is possible to state that they are located far from inhabited areas or buildings: it is therefore not generally possible to hypothesize the use of the extracted thermal energy for heating purposes.

Oilfields in Italian Sedimentary Basins
To introduce the analysis of the selected case studies, this paper provides an overview of the main sedimentary systems and associated hydrocarbon fields.
The Italian landscape is dominated by the Southern Alps and the Apennines mountain chains. The two chains are bordered in their outer margins by well-developed sedimentary foreland basins, especially along the Adriatic sectors, and by relatively wide foreland areas (e.g., the Po Plain, Adriatic Sea and Hyblean Basin). Due to Italy's complex geological and sedimentary history, a variety of petroleum systems have been identified. Authors provided overviews of the Italian peninsula's geological evolution and described how Italian hydrocarbon occurrences can be classified by their association with three main tectono-stratigraphic systems ( Figure 1) [14][15][16]: (1) carbonate Mesozoic substratum of the foredeep/foreland area and of the external thrust belts; (2) thrusted terrigenous Oligo-Miocene foredeep wedges (Southern Alps, Northern Apennines, Calabria and Sicily); and (3) terrigenous Pliocene-Pleistocene successions of the late foredeep basins of the Apennines, in both the central and northern Adriatic Sea and the Po Plain.
According to [17], at least five important source rocks have been recognised, with ages ranging from the Mesozoic through the Pleistocene. Three of the source rocks were deposited during Mesozoic crustal extension and are mainly oil-prone. Hydrocarbon occurrences associated with these sources are usually found in complex carbonate structures along the Apennines thrust-and-fold belt and in the foreland; the Villafortuna-Trecate (Po Plain), Val d'Agri-Tempa Rossa (southern Apennines) and Gela (Sicily) fields represent the largest oil accumulations pertaining to these systems.
The two other source rocks were deposited in the foredeep terrigenous units of the foreland basins, which formed during the Cenozoic orogenesis. The older of the two sources is thermogenic gas-prone and is found in the highly tectonised Oligo-Miocene foredeep wedges; gas occurrences associated with this source are mainly concentrated along the northern Apennines margin (e.g., the Cortemaggiore field), in Calabria (e.g., the Luna field) and in Sicily (e.g., the Gagliano field). The younger source is biogenic gas-prone and is situated in the outer Plio-Pleistocene foredeeps.
Reference [17] clearly showed the geographic limitation of the oil-prone Villafortuna-Trecate Middle Triassic and Val d'Agri-Tempa Rossa Cretaceous systems, as opposed to the wide distribution of the Late Triassic-Early Jurassic system and of the biogenic gas-prone Plio-Pleistocene systems ( Figure 1). This study focuses on the analysis of geothermal energy potential associated with disused hydrocarbon wells of the different petroleum systems of the Mesozoic carbonate succession. Three deep hydrocarbon wells located inside of the Villafortuna-Trecate oilfield, the Tempa Rossa Field and the Gela Field are discussed in detail.
Detailed information related to the litho-stratigraphic units and temperature data visualisation can be found on The Italian National Geothermal Database (BDNG), the largest collection of Italian geothermal data. BDNG was established in the 1980s and implemented by the Institute of Geosciences and Earth Resources (IGG) of the National Research Council (CNR) of Italy [9].
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 of the VIGOR, promoted by the MISE-DGRME (Direzione Generale Risorse Minerarie ed Energetiche), the Italian Geological Society and the Assomineraria Association.

The Villafortuna-Trecate Field
First, data were analysed for the geothermal potential associated with the disused Villafortuna 1 hydrocarbon well, 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. The petroleum system is wholly developed inside the Triassic succession and consists of two main reservoirs, composed of dolomitised carbonate platform rocks and a set of source rocks deposited in the adjacent anoxic intra-platform basins.
Generally, the depth of the main reservoir associated with the Villafortuna-Trecate field is between 5800 m and 6100 m deep, with temperatures of approximately 160-170 • C [18]. Because of its depth, the main reservoir can be pursued only in the outer sector of the foredeeps and in foreland regions (the Piedmont area); areas along the thrust belt (the Po Plain) are generally too deep.
As demonstrated by information related to the litho-stratigraphic units and temperature data visualisation reported in Table 1 and Figure 2, the stratigraphic succession is mainly composed of clastic sedimentary and carbonate rocks. The analysed well has a maximum depth of 6202 m and temperatures reaching 165 • C.

The Val d'Agri-Tempa Rossa Field
The second case study, the Tempa Rossa 1D hydrocarbon well, is located within the Tempa Rossa oilfield ( Figure 1). The Tempa Rossa system lies in the Mesozoic carbonate substratum of the foredeep and foreland areas and in the external thrust belt of the southern Apennines. The system bears the largest oil and gas accumulations in Italy, including the Val d'Agri and Tempa Rossa oil fields.
The reservoirs are composed of fractured limestones from the buried Apulia Platform, which extends from the Cretaceous to the Miocene. The majority of the oil column exceeds 1000 m and sometimes exceeds 2000 m. The seal is composed of Lower Pliocene shales. The source rocks, identified in a few deep wells of the area, are mainly Albian-Cenomanian and are marine anoxic carbonates facies containing sulphur [19].
Unlike the Villafortuna 1 hydrocarbon well, the Tempa Rossa 1D well features lithostratigraphic units mainly composed of sandstone and associated shales ( Table 2). The maximum depth of the analysed well is 5042 m; temperatures reach 107 • C ( Figure 3).

The Gela Field
The third case study is the Gela 38 hydrocarbon well, located in the Gela oilfield in Sicily. The Gela field is part of the Late Triassic-Early Jurassic petroleum system and is linked to the main phase of the Tethyan rifting. It is the most explored of the three systems, both in the foreland and in the thrust belt, and from Lombardy to Sicily. The source rocks are terrigenous or mixed carbonate-terrigenous and were deposited during the anoxic stage that preceded the extension of the Jurassic basins.
The Ragusa-Gela fields were discovered in the 1950s and have been the largest source of Italian oil for some of the past decades. The reservoir is composed of fractured, massive dolomites of the Upper Triassic Gela Formation. The traps are large-scale, probably polyphase anticlines bounded by high-angle normal faults.
Based on the available lithological and temperature data reported in Table 3 and Figure 4, the area's stratigraphic succession is composed of marl, calcareous marl and clays. The well has a maximum depth of 3446 m and temperatures that reach 85 • C. Marls, clays and gypsum

Closed-Loop Geothermal Energy Systems: WBHEs
Exploiting the geothermal energy resources associated with decommissioned hydrocarbon wells requires that the borehole to be retrofitted with a heat exchanger. In current practice, two main types of closed-loop systems are used to harness geothermal energy resources by taking advantage of disused boreholes in oilfields: U-tube and coaxial doublepipe WBHE technologies [20,21]. Both kinds of systems allow for the extraction of heat from the ground without extracting or re-injecting any geothermal fluids.
In U-tube heat exchangers, fluid is pumped through one tube string and comes out the other, while coaxial heat exchangers are composed of two concentric pipes, as shown in Figure 5 and Table 4. In coaxial WBHEs 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). The gap between internal pipes is filled with an insulating material, 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.  Corrosion must be considered when selecting this type of system; it may be ideal for heat exchange, but it may also reduce the system's operational life. Due to its low cost and its heat transfer and storage capacity, water is still one of the most commonly used fluids. The operating parameters of subsurface closed-loop systems, such as fluid flow rate and pipe diameter, should be selected to guarantee transient turbulent flow conditions, since these conditions facilitate heat transfer, and low hydraulic head losses, since this indicates lower energy expenditure on circulation pumping.
Compared to U-tube heat exchangers, coaxial heat exchangers have a higher surface area and a higher volume of the working fluid through which heat exchange occurs. As a result, under the same injection rate conditions (q), the fluid flow velocity in the coaxial pipe system and the hydraulic pressure required for fluid circulation can be lower, resulting in decreased energy consumption from pumping [21]. The coaxial geometry of a doublepipe heat exchanger also has the advantage of reducing thermal resistance between the circulating fluid and the wellbore.
In this paper, in order to evaluate and analyse the temperature profiles of the selected fluid (water) associated with a coaxial double-pipe WBHE technology, the simplified models proposed below were implemented in MATLAB and applied to the selected case studies: the Villafortuna 1, Tempa Rossa 1D and Gela 38 hydrocarbon wells.
For the elaborated model, two main assumptions were considered: the propagation of heat in the reservoir occurs by means of conduction (convection phenomena are neglected), and the propagation of heat inside the wellbore tubes takes place through both conduction and convection.

Heat Transfer in Coaxial WBHEs
In coaxial WBHEs, the steel downward tube is cemented to the rock wall and so 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: 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). Although insulation is used to prevent heat loss from the inner-pipe fluid, heat is partly transferred between the two pipes, so 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: 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 [22]: 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. Going up 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.

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 careful estimate of the parameter k w is fundamental for proper 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 ) [22]: where R s is a function of time that represents the thermal resistance due to conductive heat transfer in the rock, R a is the thermal resistance due to convective heat transfer into 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, the 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 relationship in the numerator of the second term represents the time-dependent radius of the thermal influence of the well (r s ). 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 by using the Nusselt number (Nu) and a form of the Dittus-Boelter equation that assumes turbulent flow inside the tubes (Reynolds number ≥ 104) [23]: with Pr = ρc f µ λ f and Re = ρv f 2r c µ . Finally, the 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 thermal resistance [13]. As a result, the heat exchange coefficient k w can be correctly determined as follows: where rc = rw 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: 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.

WBHE Model Assumptions
MATLAB is a software for numerical and statistical calculations that is written in the C programming language. MATLAB R2018b version was used to perform the analysis of WBHEs by implementing the proposed model, with consideration of the following assumptions and approximations.
The reservoir model was built by assuming a single well-positioned at the centre 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 in the inner tube. Due to turbulent flow, enhanced mixing phenomena occurred, decreasing the radial gradient. The temperature changed only in the annulus and in the vertical direction of the inner tube, so the temperature profile was vertically unidirectional.
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 model was 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. For the elaborations, r c (external radius of the external casing) was considered equal to r w (radius of outside wellbore).
In addition, the model considered the resistance associated with tube thickness to be negligible. The tube material had very high conductivity, so its resistance could be considered small compared to the other resistances in the system. For estimation of the resistance associated with the rock (see Equation (6)), the time value used was 3 years after the starting of the system. This assumption made the method a conservative estimate: the system in the years before (1-2 years) turns out in fact to work better. As consequence, the heat exchange phenomena are bigger with the possibility of causing overestimations.
The proposed model followed the path of the fluid using a step-by-step approach. It considered intervals of length dz in which the inlet and outlet temperatures were calculated by solving the energy balance equation for each volume dv. For estimating the energy exchange in the radial direction, the mean value of the temperature in each volume dv was used.

Results and Discussion
The temperature profiles associated with the analysed coaxial WBHE system configurations were obtained by making use of the specific ground properties of the selected case studies (the Villafortuna 1, Gela 38 and Tempa Rossa 1D hydrocarbon wells). Values relating to the thermal properties of the different rock formations have been attributed in accordance with references [11,12] (Tables 1-3).
For all cases analysed, the flow rate of the fluid and the inlet temperature were first considered as 3 kg/s and 50 • C, respectively. This temperature is a typical value for direct applications like production cycles in manufacturing and agricultural districts [2]. Subsequently, an analysis was conducted on the thermovector fluid temperature at the outlet as the inlet flow rate varied.
The sizing of the inner and outer tubes, as well as the final casing size, were fixed according to the values proposed in Figure 6 and Table 5. The thermal conductivity value of the insulating material was set to 0.025 W/mK.  In the first section of an external pipe of a coaxial WBHE system, the downward fluid is in thermal contact with both the ground on one side and the upward tube on the other.
Because of its thermal properties, the ground in contact with the external piping provides a negative heat contribution, while the inner tube of a coaxial WBHE provides a positive one. As the negative contribution is usually larger, the water temperature (working fluid temperature) decreases slightly. This behaviour can be observed by analysing the obtained temperature profiles associated with the Gela 38 and Tempa Rossa 1D hydrocarbon wells, in which the thermovector fluid was cooled at depths of up to 1800 m and 1200 m, respectively. Unlike in the Gela 38 and Tempa Rossa 1D hydrocarbon wells, the thermovector fluid maintained a constant temperature (50 • C) for the first portion of the borehole (1400 m) in the Villafortuna 1 hydrocarbon well (Figures 7a, 8a and 9a). The presence of a thick stratigraphic horizon made up of terrigenous sedimentary deposits (Table 1), which are characterised by very low conductivity and specific heat values, negatively influenced the heat exchange and so limited heat dispersion.   However, when the water and ground temperatures approached similar values in all three hydrocarbon wells, the negative contribution decreased, and water temperature increased. As the downward water profile line crossed the ground temperature line, the ground's heating process began, and the ground contribution became positive. Due to the presence of the insulating material, the heat exchange coefficient between the annulus and the inner tube was quite low, and the increase in the working fluid temperature could be mainly associated with the contribution of the ground.
Considering a cascading exploitation mode of the heat accumulated by the working fluid water in Villafortuna 1 WBHE, it is correctly possible to hypothesize a multi-variant and comprehensive use of the resource. The outflow temperature of geothermal water at the wellhead is 100 • C, which allows it to progressively be used for greenhouse heating (100-80 • C), domestic hot water and food industry (80-70 • C), animal breeding (60 • C), biomass and agricultural culture (50 • C). On the contrary, the thermal load accumulated in correspondence with Tempa Rossa 1D and Gela 38 wells turns out to be sufficient neither to justify the costs of plant retrofitting nor to plan a cascading exploitation of the geothermal fluid produced.
As can be seen in Figures 7b, 8b and 9b, the inlet flow rate strongly influenced the temperature of the wellhead thermal fluid. In each case, it is possible to identify an inlet flow rate value (kg/s) in order to obtain a higher fluid temperature at the outlet and optimize the quantity of extracted thermal power. When considering inlet flow rate values between 0.5 and 0.8 kg/s, the output fluid temperature increased up to about 122 • C (0.8 kg/s, Villafortuna 1), 65.5 • C (0.7 kg/s, Gela 38) and 72 • C (0.5 kg/s Tempa Rossa 1D). For higher values, the trend was inverted in all three cases.
However, for coaxial WBHE systems like the one analysed, such low inlet flow rates may not be technically appropriate. The heat-exchange modalities are in fact not the only aspect that must be considered to carry out the correct analysis of the heat transfer mechanisms associated with coaxial WBHE. Pressure loss phenomena need to be analysed as they affect pumping costs and are not negligible in the management of a closed-loop geothermal system. In coaxial WBHE, the inlet flow rate strongly influences the temperature of the wellhead thermal fluid. As pressure losses are proportional to the velocity, an increase in the inlet flow rate values will cause an increase in the required pumping power.

Conclusions
Over time, the petroleum systems in Italian sedimentary basins have been explored for oil and gas extraction. Since 1985, 7246 wells have been drilled for hydrocarbon extraction. Of these, 898 wells are located onshore with various operational statuses. Geological and geophysical exploration campaigns have ascertained the coexistence of low-to medium-temperature geothermal energy resources in the deepest regions of such geological contexts. As such, energy production based on the exploitation of available lowto medium-temperature geothermal resources associated with disused deep oil and gas wells in Italian oilfields could represent a considerable source of renewable energy.
This study used information on Italian hydrocarbon wells and on-field temperatures to apply a simplified coaxial WBHE model to three hydrocarbon wells. The main purpose was to analyse heat exchange mechanisms in three different Italian oil and gas fields (the Villafortuna-Trecate, Val d'Agri and Gela fields), emphasizing differences in the quantity of extracted thermal energy and considering different geological and depositional contexts.
All calculations considered the detailed stratigraphic data and related thermophysical parameters (e.g., thermal conductivity, volumetric heat capacity and rock density) of each case study. The results indicate a substantial difference in the potential amount of extracted thermal energy between analysed sites, located in different Italian sedimentary contexts. With a fixed inlet working fluid flow rate of 3 kg/s and a fixed temperature of 50 • C, thermal power values were evaluated at 627.9 KW for Villafortuna 1, 75.3 KW for the Tempa Rossa 1D well and 100.5 KW for the Gela 38 hydrocarbon well. The Villafortuna 1 WBHE recorded a maximum extracted fluid temperature of 100 • C, which allows us to hypothesize a cascading exploitation mode of the heat accumulated. Unlike Villafortuna 1, Tempa Rossa and Gela fields had thermal load values that would be of no practical use. As such, the implementation of a coaxial borehole heat exchanger in a such hydrocarbon well may not be energetically or economically worthwhile.
Improving the accuracy of the proposed models by means of future analysis is required: the basic assumption related to the constancy of the properties of the water as working fluid must be overcome by properly analysing the possibility of having phase change (evaporation) in the well, which would change the proposed models. The role in heat transfer and performance of extracting heat from abandoned wells of intraformational flows also needs to be properly considered.
Despite the above considerations, the analysis approach with the associated simplified model proposed in this paper could represent a useful simplified methodological tool to allow the preliminary definition of the possibility of a selected Italian hydrocarbon well to be converted into a geothermal one by means of a coaxial-WBHE technology.
After a preliminary analysis about the presence of industries and agricultural districts, it will be useful to produce a detailed evaluation of the industrial plants available in the area near the Villafortuna 1 well: a technical feasibility and cost-benefit analysis of the selected configuration in this proposed case study could represent the subject of future research work.