Coupled Geoflow Processes in Subsurface: CO_{2}Sequestration and Geoenergy Focus
View this Special IssueResearch Article  Open Access
Bisheng Wu, Tianshou Ma, Guanhong Feng, Zuorong Chen, Xi Zhang, "An Approximate Solution for Predicting the Heat Extraction and Preventing Heat Loss from a ClosedLoop Geothermal Reservoir", Geofluids, vol. 2017, Article ID 2041072, 17 pages, 2017. https://doi.org/10.1155/2017/2041072
An Approximate Solution for Predicting the Heat Extraction and Preventing Heat Loss from a ClosedLoop Geothermal Reservoir
Abstract
Approximate solutions are found for a mathematical model developed to predict the heat extraction from a closedloop geothermal system which consists of two vertical wells (one for injection and the other for production) and one horizontal well which connects the two vertical wells. Based on the feature of slow heat conduction in rock formation, the fluid flow in the well is divided into three stages, that is, in the injection, horizontal, and production wells. The output temperature of each stage is regarded as the input of the next stage. The results from the present model are compared with those obtained from numerical simulator TOUGH2 and show firstorder agreement with a temperature difference less than 4°C for the case where the fluid circulated for 2.74 years. In the end, a parametric study shows that () the injection rate plays dominant role in affecting the output performance, () higher injection temperature produces larger output temperature but decreases the total heat extracted given a specific time, () the output performance of geothermal reservoir is insensitive to fluid viscosity, and () there exists a critical point that indicates if the fluid releases heat into or absorbs heat from the surrounding formation.
1. Introduction
With the rapid increasing demand for clean and renewable energy, geothermal energy has become one of the most promising alternatives for energy supplies due to its many advantages. First, geothermal energy is inexhaustible and renewable due to radioactive decay from below the earth [1] and there is a wide range of geothermal resources in the world including in the USA, Mexico, Germany, Italy, Finland, Norway, Iceland, Sweden, New Zealand, Australia, Indonesia, and China [1, 2]. Second, geothermal energy is much cleaner and environmentally friendly compared to the conventional fossil fuels. It has been reported that the geothermal energy based plants produce 1/12 carbon dioxide of the coal based plant [3]; Third, up to date, as only small amount of geothermal resource has been utilized for heating or power generation compared to its large reserve, there is a large potential for extraction.
Closedloop geothermal system represents an important mean to extract heat from below the earth and study on maximizing its output performance is very important [4, 5]. The basic principle of a closedloop system is simple: a cold fluid (CO_{2} or water) is injected into the well, absorbs heat from hotter surrounding rocks when flowing in a channel such as wellbore or fracture, and then is pumped out from the production well. In the present paper, the connecting channel is chosen to be a wellbore. In addition to extracting heat from the reservoir by fluid circulation, the closedloop system can prevent contaminants in fluids from leaking into the reservoir. As the output performance, including temperature and thermal power, of a geothermal reservoir is affected by many factors such as number and location of wells and fractures, flow rates, and velocity and direction of areal flow [6, 7], an efficient mathematical model is required for studying the relationship between reservoir performance and these parameters depending on specific system configuration.
Fluid flow and heat transfer in a wellbore/pipe have received a great deal of attention during the last several decades due to their importance in oil and gas industry. For example, the thermoelastic stress change is significant near the wellbore; it may change the fracture initiation and path during hydraulic fracturing [10–12]. The mathematical models can be mainly classified into two types, that is, analytical and numerical. Many studies belong to the first type. For example, Bullard [13] and Moss and White [14] used a line/source model to study the time for a wellbore to attain temperature equilibrium and the temperature evolution in a waterinjection well, respectively. Recently, the infinite, cylindrical, and finite line source theories [15] have been widely used for geothermal simulation. For example, Eskilson [16], Zeng et al. [17], Sutton et al. [18], Bandos et al. [19], Michopoulos and Kyriakis [20], MolinaGiraldo et al. [21], HechtMéndez et al. [4], Rivera et al. [22], and Zhou et al. [23] used line source models to predict the temperature evolution in the ground source heat pumps. In addition, Ramey [24] presented approximate solutions for the wellbore/reservoir (W/R) system involved in injection of hot or cold fluid. Dowdle and Cobb [25] applied the Horner temperature plot method, which is similar to conventional pressure buildup, to predict the static wellbore temperature from well logs. Edwardson et al. [26] and Tragesser et al. [27] studied the wellbore temperature during mud circulation by using an exact method which is based on the solution of differential equation of heat conduction. Because the above models do not consider the heat exchange between the fluid inside and that outside the drill pipe, they cannot be used as a general tool for predicting the thermal behaviour in a wellbore/reservoir systems.
As for fluid circulation in a wellbore/reservoir system, the work by Raymond has to be mentioned [28]. In addition to providing systematic derivation of the governing equations characterizing the transient fluid flow and heat transfer, Raymond [28] obtained the analytical solutions for the case with steady state heat transfer in the wellbore by using Laplace transformation and also numerically solved the problem with transient heat transfer in both wellbore and formation by using the finite difference method (FDM). Most of the models studied later on fluid circulation in a wellbore/reservoir system more or less borrow some concepts from Raymond’s work [28]. For example, Holmes and Swift [29] solved approximately similar governing equations under assumption of steady state heat transfer between the annular fluid and the formation. Keller et al. [30] considered the effects of heat generated by fluid friction, sliding casing strings, and other energy sources in the system. García et al. [31], EspinosaParedes et al. [32], Fomin et al. [33], and Izgec et al. [34] presented a fully transient FDM model by considering the transient heat behaviour during fluid circulation in a W/R system. Wu et al. [35] solved analytically the transient heat transfer both in the wellbore and formation during drilling with fluid circulation in the wellbore.
In addition to the abovementioned approaches, recently a large number of numerical open sources or commercial software are available for geothermal simulation. For example, HYDROTHERM by US Geological Survey [36] and TOUGH2 by Lawrence Berkeley National Laboratory [37] based on FDM, OpenGeoSys by Helmholtz Centre of Environmental Research, Germany [38], ROCKFLOW by Kolditz et al. [39], FEFLOW by Diersch et al. [40, 41], FEHM by Los Alamos National Laboratory [42] and COMSOL [43] based on finite element method (FEM), ECLIPSE by Schlumberger [44], and Fluent by ANSYS [45] based on finite volume method (FVM) have been applied successfully to model the heat extraction from geothermal reservoirs with complicated geometries. In addition, the boundary element method (BEM) was also used for modelling similar problems. For example, Ding [46], Ghassemi et al. [47], Kumar and Gutierrez [48], and McClure and Horne [49] used the BEM or displacement discontinuity element method (DDM) to simulate the heat and mass transfer in geothermal reservoir with single fracture or fracture networks.
Although the above numerical methods can be used to solve complex problems, they sometimes are computationally costly. The analytical methods cannot be directly used to solve the present closedloop geothermal model because of the complex model geometry. This is the main motivation for carrying out the present work.
The objective of the present work is to provide an approximate approach to predict the heat extraction and thermal power from a closedloop geothermal reservoir. This paper is organized as follows: in Section 2, the problem description and formulation, including governing equations and boundary and initial conditions, are presented; Section 3 provides the dimensionless formulation, followed by solution method in Section 4. The model validation and numerical results are given in Section 5 and conclusions are presented in Section 6.
2. Problem Formulation
2.1. Problem Description and Assumptions
The model geometry studied is shown in Figure 1. It contains one simplified well containing two vertical parts and one horizontal part. The initial temperature of the system is a function of depth; that is, T_{0} = A_{0}z + B_{0}, where A_{0} is the geothermal gradient and B_{0} is the surface ground temperature. At time , a fluid (water) with a constant temperature is injected into the system at a constant volumetric rate Q_{1}. The geometrical sizes of the wellbore shown in Figure 1 are defined as follows: the lengths of the injection, horizontal, and production well sections are L_{1}, L_{2}, and L_{3} (L_{3} = L_{1}), respectively; the radius of the wellbore is (), possibly including all layers such as steel casing and cementation if they exist as shown in Figure 1. The origin of the general coordinate system (GCS), denoted by the coordinates (), is located at the injection point of the vertical well and coincides with the origin of the local cylindrical coordinate system (LCCS), denoted by the coordinate (r_{1}, z), for the injection well. The origin of the LCCS, denoted by the coordinate (r_{2}, x), for the horizontal well is at the bottom of the injection well, that is, intersection . In a similar way, the origin of the local cylindrical coordinate system (LCCS), denoted by the coordinate (r_{3}, z_{1}), for the production well is located at the bottom of the production well, that is, (L_{2}, 0, L_{1}) in the GCS.
(a)
(b)
In order to make the problem tractable, some assumptions are made:(a)The fluid is single phase, incompressible, and Newtonian, and the rock is impermeable.(b)The material properties of the fluid and the rock are constants, independent of temperature.(c)When there are multiple layers around the wellbore, the heat transfer through the layers is characterized by one overall heat transfer coefficient.(d)The heat diffusion around the wellbore can be regarded to be axissymmetric.
Assumption (d) is the most important one based on which the present model is simplified. According to JiJi [50, page 242–246], for a semiinfinite domain with a homogenous initial temperature and one timedependent surface heat flux imposed at the boundary at time , the penetration depth of the thermal layer is in the order of , where /() is the thermal diffusivity with , , and being the thermal conductivity, mass density, and specific heat capacity, respectively, of the media. Take shale, for example, when W/(m·K), Kg/m^{3}, and J/(Kg·K) (pp 106, [8]), the thermal diffusivity m^{2}/s. This means that, after 60 days and 30 years, the thermal layer is about 3.16 m and 42.68 m, respectively, illustrating the slow movement of the thermal front in rock. Because the reservoir size is large enough (say in the order of 1 km), the interaction of heat transfer in the horizontal and vertical directions around the bottom of the vertical wells can be neglected.
Based on the above assumptions, the current model can be simplified greatly without causing large errors in predicting the temperature profiles.
2.2. Governing Equations
The governing equations can be written in the corresponding local cylindrical coordinate system (LCCS) for the vertical and horizontal parts of the well. Besides, temperature continuity is enforced at the intersection point at the heel of the horizontal well.
2.2.1. Heat Exchange along the Wellbore
As the heat flow is uniform along the well, the transient energy equations for the fluid flow are written as follows, according to their LCCSs [28]:for the vertical wells, and for the horizontal well, where and are the mass density and specific heat capacity, respectively, of the fluid and , , , , , and ( for injection, for horizontal, and for production) denote the radius of fluid flow channel, the areas of the wellbore cross section, overall heat transfer coefficients (OHTCs), fluid velocity in the wellbore, fluid temperature, and the temperature at the wellbore wall, respectively. The value of the radius of fluid flow channel will be determined based on the condition as follows. If fluid flows in a pipe, it is equal to the pipe radius; if the fluid contacts directly with the rock formation, it is equal to the wellbore radius; that is, .
The OHTCs can be calculated based on Willhite’s equation [51] which considers a multilayered wellbore structurewhere or 2, for two well parts. The crosssection geometry can be found in Figure 1(b) and the subscripts provide the layers. k denote the thermal conductivities for different layers. is the natural convection and conduction HTC through the annulus, and is the radiation HTC through the annulus; h_{ℓ} is the HTC between the fluid and tubing. If the effects of the annulus, casing, and cements are not taken into account, the first two terms on the right side of (3) are used to calculate U_{ℓ}. If the fluid contacts directly with the rock formation, and U_{ℓ} = h_{ℓ}. From (3) we also know that the OHTC is mainly determined by the minimum value of the denominators of the terms on the right side. This indicates that if the thermal conductivity of some layer is very small, this layer will work as a thermal insulator.
When fluid flows in a tubing, the HTC between the fluid and tubing, that is, , is obtained by using the relationship , where denotes the Nusselt number, is the thermal conductivity of the fluid, and D is the hydraulic diameter of the tubing. For fully developed laminar flow in a pipe with circular cross section, the Nusselt number , while, for transitional and turbulent flows, the Nusselt number is obtained by using the wellknown Gnielinski correlation [52]when and . In the present model, the Reynolds and Prandtl numbers related to fluid flowing in the tubing or wellbore and the rock formation are expressed as where and denote the mass density and heat capacity, respectively, of the fluid and and denote the thermal conductivity and heat capacity, respectively, of the rock formation.
2.2.2. Heat Conduction in the Rock Formation Surrounding the Wellbore
As the radial heat conduction dominates the thermal diffusion process around the wellbore, the equations are written as follows in their corresponding LCCS:where T_{r1}, T_{r2}, and T_{r3} denote the temperature of the formation around the injection, horizontal, and production wells, respectively. The thermal diffusivities of the rock formation are denoted by , with being the mass density of the rock formation. It should be noted that the radial spatial variables for the injection, horizontal, and production wells are r_{1}, r_{2}, and r_{3}, respectively.
The heat transfer conditions between the rock and the whole well system are given as
Based on the above equation, if the OHTC , , indicating no heat flowing into or out of the rock formation. From the point of view of heat extraction, a thermal insulation layer along the production well can enhance the heat production by preventing the heat loss along the well during fluid flowing upwards.
2.3. Boundary and Initial Conditions
The injection rate is prescribed and the injection temperature is
The bottomhole temperature (BHT) of the vertical well is used as the input conditions for the horizontal well and the output temperature of the horizontal well is used as the input conditions for the production well, that is, at the intersections A and B
The ground temperature at the surface is a constant and the initial temperature of the whole system is a function of depth
The heat extraction rate or thermal power output by the fluid with output temperature is expressed aswhere and from which the total heat extracted is
3. Dimensionless Formulation
The governing equations, boundary and initial conditions, are simplified with the following transformation:where , , and ( for injection, for horizontal, and for production) denote dimensionless temperature in the fluid, at the wall and in the formation, respectively.
By using the above variables, the simplified formulation for the three subproblems of the original model is listed as follows.
(a) Injection Well. The governing and heat balance equations for the injection well arewith the initial and boundary conditions
(b) Horizontal Well. The governing and heat balance equations for the horizontal well arewith the initial and boundary conditions becoming
(c) Production Well. The governing and heat balance equations for the production well arewith the initial and boundary conditionswhere denotes the distance of the point from the bottom of production well, . The subproblem for production well is similar to that for injection well as it can be regarded as a model with the origin of the CS at the bottom of production well and with a negative geothermal gradient .
4. Solution Method
By using the Laplace transformation, the analytical solutions for the three stages are obtained. In the following equations, the symbol denotes the variables which are Laplace transformed and is a complex number as the Laplace symbol. For example, denotes the Laplace transform of the function .
4.1. Injection Well
From the Laplace transform of the governing equation, that is, the first of (15), for the heat diffusion in the rock around the injection well, we obtain the rock temperaturewhere is the Laplace transform of the dimensionless temperature ; and are the modified Bessel functions of the first and second kind of order n, respectively; and F_{1}(s) and F_{2}(s) are unknowns to be determined by the boundary conditions. As is finite when , it is easy to know that . Therefore, the temperature at the wellbore wall is obtained
By substituting (21) and (22) into the second of (15) we have which, after being used in the third of (15), produces an ordinary differential equation with respect to the unknown ; that is,
The solution of the above equation is found to bewhere the unknown function C_{1}(s) can be determined by using (23), (25), and the injection boundary condition (16)
Therefore, the solutions for the temperatures of the rock, , and fluid, , are obtained and rewritten aswhere the functions A_{1}, A_{2}, and A_{3} are defined aswhere is the Laplace transform of the injection temperature.
The Laplace transformation of the bottomhole temperature (BHT) of the injection well iswhich will be used as the input condition for the horizontal well.
4.2. Horizontal Well
Initially, the temperature along horizontal well is identical due to the same depth. Based on the above calculations, the thermal front only penetrates less than 50 meters into the rock within 30 years. In the region between the upper and lower thermal front, the largest temperature is less than 4°C. Therefore, it is reasonable to assume a constant temperature in the rock formation around the horizontal well.
In a similar way, the solutions for the horizontal well are obtainedwhere the functions B_{1} and B_{3} are defined as
The Laplace transformation of the temperature at the end of the horizontal well, also BHT of the production well, is obtained by using ; that is,which will be used as the input condition for the production well.
4.3. Production Well
The temperatures for the fluid in the production well and neighboring rock formation can be obtained in the same way as those for injection and horizontal well and the calculation details are omitted here where the functions C_{1}, C_{2}, and C_{3} are defined as
The Laplace transformation of the temperature at the end of the horizontal well, also BHT of the production well, is obtained by using ; that is,which will be used as the input condition for the production well.
5. Validation and Numerical Results
5.1. Method Validation
In the present model, in order to obtain analytical solutions for temperature prediction in a closedloop geothermal reservoir, the wellbore for fluid flow and heat transfer is divided into three parts with the output of one well part as the input conditions of the next well part. In order to show that these assumptions are reasonable and to show the error resulting from decoupling of the whole process, the results from the present model are compared with those from numerical simulator TOUGH 2. An automatic Laplace inversion technique developed by D’Amore et al. [53] based on Fourier series is used in the present model to obtain the values in the time space.
Figure 2 shows spatial discretization of the rock formation for the numerical calculation by using TOUGH2. As the temperature of the fluid in the radial direction in the pipe is assumed to be the same, it is denoted by blue node. The node for the rock is located in the middle of the each element. Figure 3 compares the temperature variations at points A, B, and C, which corresponds to the output points of the injection, horizontal, and production well parts, respectively, between TOUGH2 [37] and the present model. The solid curves denote results from TOUGH2 and the dashed ones are results from the present model. Figures 3(a) and 3(b) are for temperatures at short (20 days) and long (2.74 years) times of circulation, respectively, for the case m^{3}/s, while Figures 3(c) and 3(d) are for small m^{3}/s and large m^{3}/s, respectively. Other parameters used here are m, m, B_{0} = 20°C, and °C.
(a)
(b)
(c)
(d)
It can be seen that the temperature results at points A, B, and C predicted from both methods do not match at short time period, say 10 days, as shown in Figure 3(a) for the case when m^{3}/s. The temperature difference is mainly caused by the neglect of thermal interaction at the heel of the injection and production wells (i.e., around points A and B) in the present model. However, for large time of circulation, the difference in the temperature results from both methods is very small, generally less than 4°C, as shown in Figure 3(b). The temperature differences after 365 days and 1000 days between the present model and TOUGH2 are −3.87°C and −3.65°C, respectively. In order to further show the accuracy of the present model, the temperatures at these three points in another two cases, that is, m^{3}/s and m^{3}/s, are also compared with those obtained from TOUGH2 and show good agreement. This means that the present model can be used to predict approximately the output temperature of the closedloop system without causing large errors.
5.2. Numerical Results
In this section, the case with fluid fully contacted (without any tubing, casing, and cement) with the wellbore is first investigated. Under this condition, the radius of the tubing for fluid flow, , is equal to the wellbore radius, . By doing so, two objectives are achieved. First, the important factors affecting the fluid and heat flow behaviour are identified; second, the critical location determining whether the heat is flowing into or out of the rock formation is identified for thermal isolation design, especially for production well. Water is chosen for the fluid; sandstone and shale are chosen for the rock around the vertical and horizontal wells, respectively, for the following calculations. The parameters used for the examples are listed in Table 1 unless otherwise specified.

5.2.1. Temperature along the Vertical Wells
Figure 4 displays the dimensionless temperature changes, which are defined by physical temperature changes with respect to the initial state divided by (A_{0}H + B_{0}), along the injection and production wells, respectively. The red curves refer to the initial state. For the injection well, we can see from Figure 4(a) that (1) before the cold fluid reaches the bottom of the wellbore, the temperature change of the fluid on the upper part of the wellbore decreases linearly with depth, while that on the lower part of the wellbore is almost a constant; (2) the temperature along the injection well decreases very quickly, with a temperature of 46.3°C at the bottom after 3 hours of circulation; and (3) as the injection temperature is lower than the initial ground temperature for most of the well length, most of the fluid flowing downwards absorbs heat from the surrounding rock except a small portion of fluid on the upper part of the wellbore.
(a)
(b)
Before the cold fluid reaches the production well, the temperature along the whole production well increases because of hot fluid flowing upwards, as shown by solid and dashed curves in Figure 4(b). The production temperature increases from 57°C to 152°C when the time varies from 10 mins to 1.0 hours and then decreases to 126°C after 3 hours. The small difference in temperature change between times days and years means that the heat transfer in the fluid approaches pseudosteady state after 10 days. Moreover, the critical position which indicates if the fluid absorbs heat from or releases heat to the rock formation is around or m.
5.2.2. Effect of Flow Rates on Output Performance
Figure 5 presents the output temperature for different injection rates ranging from m^{3}/s to 0.14 m^{3}/s. It can be seen that the injection rate plays a major influence on the output temperature. First, for each injection rate, the output temperature increases quickly to a maximum value, as shown in Figure 5(b), and then decreases steadily to a constant value, as shown in Figure 5(a). Second, with increasing the injection rate, the production temperature is decreased. This is due to shortened time as a result of larger fluid speed for heat exchange between fluid and the surrounding rock. In addition, the temperature difference after 30 years between the cases m^{3}/s and m^{3}/s is around 40°C. However, when the injection rate is increased to some value (0.05 m^{3}/s for the present case), there is no large difference in the final production temperature, as shown in Figure 5(a).
(a)
(b)
The thermal power and total heat extracted for the above cases are plotted in Figure 6. According to the definition of thermal power in (12), thermal power is linearly proportional to the injection rate and temperature. However, as lower production temperature is obtained when the injection rate is increased, the curves for thermal power shown in Figure 6(b) do not exhibit similar trends to the production temperature in Figure 5(b). It has to be mentioned that the longtime thermal power ranges from 1.71 MWs to 2.21 MWs when is changed from 0.01 m^{3}/s to 0.14 m^{3}/s.
(a)
(b)
(c)
Figure 6(c) displays the total heat extracted from the reservoir up to 30 years under different injection rates. The total heat extracted is expected to increase when the injection rate is increased. However, the total heat extracted after 30 years when m^{3}/s is only 32.75% larger than that when m^{3}/s. This means that the cost related to higher pressure required for larger flow rates has to be balanced with the total energy extracted.
5.2.3. Effect of Injection Temperatures on Output Performance
The effect of injection temperatures on output performance is displayed in Figures 7 and 8. As for output temperature, it increases evenly with increasing the injection temperature, as shown in Figure 7(a), and the maximum value of output temperature is not sensitive to the injection temperature, as shown in Figure 7(b).
(a)
(b)
(a)
(b)
(c)
Figure 8 presents the thermal power and total heat extracted up to 30 years. As the thermal power is linearly proportional to the temperature and the injection rate is kept constant, the thermal power exhibits a similar trend to the output temperature, as shown in Figures 8 and 7. In addition, it is interesting to find that although higher output temperature is obtained under larger injection temperature, the total heat extracted from the reservoir is reduced, as shown in Figure 8(c).
5.2.4. Effect of Fluid Viscosity on Output Performance
Figures 9 and 10 display the final output performance for different fluid viscosities. The water viscosity ranges from about 2.8 × 10^{−4} Pa·s to 1.3 × 10^{−3} Pa·s when the temperature changes from 100°C to 10°C, which can be found in the website link (https://en.wikipedia.org/wiki/viscosity). Here = 10^{−4}, 5 × 10^{−4}, 2.5 × 10^{−3}, and 1.2 × 10^{−2} Pa·s are chosen to show the effect of the viscosity on the output temperature. From Figure 9 the output temperature follows almost the same curve under these four cases and the decreasing trend for rapid temperature drop in less than 5 hrs is insensitive to the fluid viscosity.
(a)
(b)
(a)
(b)
(c)
5.2.5. Effect of Geothermal Gradient on the Output Performance
The output temperatures under different gradient geothermal gradients, A_{0}, are displayed in Figure 11 when other parameters are kept the same. It can be found that (1) as the depth of the vertical well is the same, larger geothermal gradient produces larger maximum output temperature (MOT), which is increased by around 42.2°C if has an increase of 1.5°C. For example, the MOTs for the cases with A_{0} = 0.020, 0.035, 0.050, 0.065, 0.080, and 0.095 K/m are 76.2, 118.4, 160.5, 202.7, 244.9, and 287.0°C, respectively. However, the output temperature after a long time (say 30 years) is found to be increased by around 4°C when A_{0} has an increase of 1.5°C, as shown in Figure 11(b), where the output temperature is 34.2, 38.2, 42.1, 46.0, 49.9, and 53.8°C, respectively, for the above six cases. This means that, compared to the injection conditions, the geothermal gradient plays much less important role in the output performance.
(a)
(b)
5.2.6. Fully Insulated Case
Figure 12 compares the production temperatures under fully contacted and fully insulated conditions. The parameters used for the fully insulated case are the same as those in Figure 4 except that the overall heat transfer coefficient along the production well is set to zero.
It is found from Figure 12 that there is no large difference in production temperature between these two cases and the production temperature predicted from fully insulated case is a little smaller than that predicted from fully contacted case. From the curves for days and years in Figure 4(b) without any casing and cementing, it can be found that (1) the final output temperature decreases very quickly and approaches some pseudosteady value which is dominated by the inject conditions and (2) the fluid along a larger part (around 80 percent) of wellbore length absorbs heat from the surrounding rock formation and the fluid along the rest of the wellbore releases heat into the formation. Compared with this openhole case, if an insulated layer is put along the whole wellbore, the final output temperature will become smaller although the heat loss into the formation, which occurs along a small part of wellbore length in the openhole case, is prevented in the insulated case, as confirmed in Figure 12.
6. Conclusion
This paper deals with the heat extraction from a closedloop geothermal system. Through reasonable assumptions and the temperature continuity conditions at the intersection, the whole well is divided into three portions, two vertical wellbores and one horizontal wellbore, which can be solved independently based on the input conditions. Conclusions that can be made based on our numerical results are as follows:(1)The difference in the output temperature predicted from the proposed model and TOUGH2 is very small, less than 4°C in the present numerical study, and thus can be used to approximately predict the heat extraction from the closedloop geothermal system.(2)The injection rate plays the dominant role in affecting the output performance; although it can increase the thermal power, the reduced temperature as a result of high flow rate compromises its production performance.(3)Higher injection temperature produces larger output temperature but decreases the total heat extracted given a specific time.(4)The output performance of geothermal reservoir is insensitive to fluid viscosity.(5)There exists a critical point that indicates if the fluid above and below this point releases heat into or absorbs heat from surrounding formation.(6)This approximation model proposed in this work runs in terms of seconds on a personal notebook computer and thus provides an efficient tool for reservoir optimization.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors thank CSIRO for permission to publish the research outcomes. This work was supported by the Scientific Research Foundation of International Cooperation and Exchanges of Sichuan Province (Grant no. 2017HH0061), the National Natural Science Foundation of China (Grant no. 51604230), the China Postdoctoral Science Foundation (Grant nos. 2016M600626 and 2017T100592), and the Fund of the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation of Southwest Petroleum University (Grant nos. G201604 and PLN201611).
References
 J. W. Lund and T. L. Boyd, “Direct utilization of geothermal energy 2015 worldwide review,” Geothermics, vol. 60, pp. 66–93, 2016. View at: Publisher Site  Google Scholar
 A. Bahadori, S. Zendehboudi, and G. Zahedi, “RETRACTED: A review of geothermal energy resources in Australia: Current status and prospects,” Renewable and Sustainable Energy Reviews, vol. 21, pp. 29–34, 2013. View at: Publisher Site  Google Scholar
 J. M. Sharp Jr., “Energy and momentum transport model of the Ouachita Basin and its possible impact on formation of economic mineral deposits,” Economic Geology, vol. 73, no. 6, pp. 1057–1068, 1978. View at: Publisher Site  Google Scholar
 J. HechtMéndez, M. De Paly, M. Beck, and P. Bayer, “Optimization of energy extraction for vertical closedloop geothermal systems considering groundwater flow,” Energy Conversion and Management, vol. 66, pp. 1–10, 2013. View at: Publisher Site  Google Scholar
 A. Casasso and R. Sethi, “Efficiency of closed loop geothermal heat pumps: A sensitivity analysis,” Renewable Energy, vol. 62, pp. 737–746, 2014. View at: Publisher Site  Google Scholar
 B. Wu, X. Zhang, A. P. Bunger, and R. G. Jeffrey, “An efficient and accurate approach for studying the heat extraction from multiple recharge and discharge wells,” in Effective and Sustainable Hydraulic Fracturing, A. Bunger, J. McLennan, and R. Jeffrey, Eds., p. 16, InTech, 2013. View at: Google Scholar
 B. Wu, G. Zhang, X. Zhang, R. G. Jeffrey, J. Kear, and T. Zhao, “Semianalytical model for a geothermal system considering the effect of areal flow between dipole wells on heat extraction,” Energy, vol. 138, pp. 290–305, 2017. View at: Publisher Site  Google Scholar
 L. Eppelbaum, I. Kutasov, and A. Pilchin, Applied Geothermics, Lecture Notes in Earth System Sciences, SpringerVerlag, Berlin, Germany, 2014.
 H. Quick, J. Michael, U. Arslan, and H. Huber, “Geothermal application in lowenthalpy regions,” Renewable Energy, vol. 49, pp. 133–136, 2013. View at: Publisher Site  Google Scholar
 B. Wu, X. Zhang, R. G. Jeffrey, and B. Wu, “A semianalytic solution of a wellbore in a nonisothermal lowpermeability porous medium under nonhydrostatic stresses,” International Journal of Solids and Structures, vol. 49, no. 13, pp. 1472–1484, 2012. View at: Publisher Site  Google Scholar
 T. Ma, Q. Zhang, P. Chen, C. Yang, and J. Zhao, “Fracture pressure model for inclined wells in layered formations with anisotropic rock strengths,” Journal of Petroleum Science and Engineering, vol. 149, pp. 393–408, 2017. View at: Publisher Site  Google Scholar
 T. Ma, B. Wu, J. Fu, Q. Zhang, and P. Chen, “Fracture pressure prediction for layered formations with anisotropic rock strengths,” Journal of Natural Gas Science and Engineering, vol. 38, pp. 485–503, 2017. View at: Publisher Site  Google Scholar
 E. C. Bullard, “The time necessary for a bore hole to attain temperature equilibrium,” Geophysical Journal International, vol. 5, pp. 127–130, 1947. View at: Publisher Site  Google Scholar
 J. T. Moss and P. D. White, “How to calculate temperature profiles in a water injection well,” The Oil and Gas Journal, vol. 57, pp. 174–177, 1959. View at: Google Scholar
 H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids, The Clarendon Press, Oxford, UK, 1st edition, 1946. View at: MathSciNet
 P. Eskilson, Thermal analysis of heat extraction boreholes [Ph.D. thesis], University of Lund, Lund, Sweden, 1987.
 H. Y. Zeng, N. R. Diao, and Z. H. Fang, “A finite linesource model for boreholes in geothermal heat exchangers,” Heat Transfer—Asian Research, vol. 31, no. 7, pp. 558–567, 2002. View at: Publisher Site  Google Scholar
 M. G. Sutton, D. W. Nutter, and R. J. Couvillion, “A ground resistance for vertical bore heat exchangers with groundwater flow,” Journal of Energy Resources Technology, Transactions of the ASME, vol. 125, no. 3, pp. 183–189, 2003. View at: Publisher Site  Google Scholar
 T. V. Bandos, Á. Montero, E. Fernández et al., “Finite linesource model for borehole heat exchangers: effect of vertical temperature variations,” Geothermics, vol. 38, no. 2, pp. 263–270, 2009. View at: Publisher Site  Google Scholar
 A. Michopoulos and N. Kyriakis, “Predicting the fluid temperature at the exit of the vertical ground heat exchangers,” Applied Energy, vol. 86, no. 10, pp. 2065–2070, 2009. View at: Publisher Site  Google Scholar
 N. MolinaGiraldo, P. Blum, K. Zhu, P. Bayer, and Z. Fang, “A moving finite line source model to simulate borehole heat exchangers with groundwater advection,” International Journal of Thermal Sciences, vol. 50, no. 12, pp. 2506–2513, 2011. View at: Publisher Site  Google Scholar
 J. A. Rivera, P. Blum, and P. Bayer, “Influence of spatially variable ground heat flux on closedloop geothermal systems: Line source model with nonhomogeneous Cauchytype top boundary conditions,” Applied Energy, vol. 180, pp. 572–585, 2016. View at: Publisher Site  Google Scholar
 G. Zhou, Y. Zhou, and D. Zhang, “Analytical solutions for two pile foundation heat exchanger models in a doublelayered ground,” Energy, vol. 112, pp. 655–668, 2016. View at: Publisher Site  Google Scholar
 H. Ramey, “Wellbore heat transmission,” Journal of Petroleum Technology, vol. 14, no. 04, pp. 427–435, 2013. View at: Publisher Site  Google Scholar
 W. L. Dowdle and W. M. Cobb, “Static formation temperature from well logs  an empirical method,” JPT, Journal of Petroleum Technology, vol. 27, pp. 1326–1330, 1975. View at: Publisher Site  Google Scholar
 M. Edwardson, H. Girner, H. Parkison, C. Williams, and C. Matthews, “Calculation of formation temperature disturbances caused by mud circulation,” Journal of Petroleum Technology, vol. 14, no. 04, pp. 416–426, 2013. View at: Publisher Site  Google Scholar
 A. Tragesser, P. B. Crawford, and H. R. Crawford, “A method for calculating circulating temperatures,” Journal of Petroleum Technology, vol. 19, no. 11, pp. 1507–1512, 2013. View at: Publisher Site  Google Scholar
 L. Raymond, “Temperature distribution in a circulating drilling fluid,” Journal of Petroleum Technology, vol. 21, no. 03, pp. 333–341, 2013. View at: Publisher Site  Google Scholar
 C. S. Holmes and S. C. Swift, “Calculation of circulating mud temperatures,” Journal of Petroleum Technology, vol. 22, no. 6, pp. 670–674, 1970. View at: Publisher Site  Google Scholar
 H. Keller, E. Couch, and P. Berry, “Temperature distribution in circulating mud columns,” Society of Petroleum Engineers Journal, vol. 13, no. 1, pp. 23–30, 2013. View at: Publisher Site  Google Scholar
 A. García, E. Santoyo, G. Espinosa, I. Hernández, and H. Gutiérrez, “Estimation of temperatures in geothermal wells during circulation and shutin in the presence of lost circulation,” Transport in Porous Media, vol. 33, no. 12, pp. 103–127, 1998. View at: Publisher Site  Google Scholar
 G. EspinosaParedes, A. Garcia, E. Santoyo, and I. Hernandez, “TEMLOPI/V.2: A computer program for estimation of fully transient temperatures in geothermal wells during circulation and shutin,” Computers and Geosciences, vol. 27, no. 3, pp. 327–344, 2001. View at: Publisher Site  Google Scholar
 S. Fomin, T. Hashida, V. Chugunov, and A. V. Kuznetsov, “A borehole temperature during drilling in a fractured rock formation,” International Journal of Heat and Mass Transfer, vol. 48, no. 2, pp. 385–394, 2005. View at: Publisher Site  Google Scholar
 B. Izgec, C. S. Kabir, D. Zhu, and A. R. Hasan, “Transient fluid and heat flow modeling in coupled wellbore/reservoir systems,” SPE Reservoir Evaluation and Engineering, vol. 10, no. 3, pp. 294–301, 2007. View at: Publisher Site  Google Scholar
 B. Wu, X. Zhang, and R. G. Jeffrey, “A model for downhole fluid and rock temperature prediction during circulation,” Geothermics, vol. 50, pp. 202–212, 2014. View at: Publisher Site  Google Scholar
 K. L. Kipp, P. A. Hsieh, and S. R. Charlton, “Transport Simulator: HYDROTHERM — Version 3,” U.S. Geological Survey, 2008. View at: Google Scholar
 K. Pruess, C. Oldenburg, and G. Moridis, “TOUGH2 User's Guide Version 2,” Tech. Rep. LBNL43134, Lawrence Berkeley National Laboratory, Berkeley, Calif, USA, 1999. View at: Publisher Site  Google Scholar
 N. Böttcher, N. Watanabe, U. J. Görke, and O. Kolditz, Modeling I: Geothermal Processes in Fractured Porous Media, Springer, 2017.
 O. Kolditz, A. Habbar, R. Kaiser et al., ROCKFLOW users manual release 3.5. Institute of fluid mechanics and computer applications in civil engineering, University of Hannover, 2001.
 H.J. G. Diersch, D. Bauer, W. Heidemann, W. Rühaak, and P. Schätzl, “Finite element modeling of borehole heat exchanger systems. Part 2. Numerical simulation,” Computers and Geosciences, vol. 37, no. 8, pp. 1136–1147, 2011. View at: Publisher Site  Google Scholar
 H.J. G. Diersch, FEFLOW: Finite Element Modeling of Flow, Mass and Heat Transport in Porous and Fractured Media, Springer, Berlin, Germany, 2014. View at: Publisher Site
 N. Tenma, K. Yasukawa, and G. Zyvoloski, “Model study of the thermal storage system by FEHM code,” Geothermics, vol. 32, no. 4, pp. 603–607, 2003. View at: Publisher Site  Google Scholar
 K. J. Bakhsh, M. Nakagawa, M. Arshad, and L. Dunnington, “Modeling thermal breakthrough in sedimentary geothermal system,” in Proceedings of the 41st Workshop on Geothermal Reservoir Engineering Stanford University, pp. 22–24, Stanford, Calif, USA, 2016. View at: Google Scholar
 Schlumberger, “Eclipse, reference manual 2008,” Tech. Rep., Schlumberger Information Solutions, Business Development Central and Eastern Europe, Hannover, Germany, 2008. View at: Google Scholar
 V. Hamm and B. Bazargan Sabet, “Modelling of fluid flow and heat transfer to assess the geothermal potential of a flooded coal mine in Lorraine, France,” Geothermics, vol. 39, no. 2, pp. 177–186, 2010. View at: Publisher Site  Google Scholar
 Y. Ding, “Using boundary integral methods to couple a semianalytical reservoir flow model and a wellbore flow model,” in Proceedings of the 15th Symposium on Reservoir Simulation, pp. 195–205, February 1999. View at: Google Scholar
 A. Ghassemi, S. Tarasovs, and A. H.D. Cheng, “An integral equation solution for threedimensional heat extraction from planar fracture in hot dry rock,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 27, no. 12, pp. 989–1004, 2003. View at: Publisher Site  Google Scholar
 D. Kumar and M. Gutierrez, “Threedimensional heat flow model for enhanced geothermal systems using boundary element method,” in Proceedings of the 38th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 2013. View at: Google Scholar
 M. W. McClure and R. N. Horne, “An investigation of stimulation mechanisms in Enhanced Geothermal Systems,” International Journal of Rock Mechanics and Mining Sciences, vol. 72, pp. 242–260, 2014. View at: Publisher Site  Google Scholar
 L. M. Jiji, Heat Conduction, Springer, Berlin, Germany, 3rd edition, 2009.
 G. Willhite, “Overall heat transfer coefficients in steam and hot water injection wells,” Journal of Petroleum Technology, vol. 19, no. 05, pp. 607–615, 2013. View at: Publisher Site  Google Scholar
 V. Gnielinski, “New equations for heat and mass transfer in turbulent pipe and channel flow,” in Proceedings of the International Chemical Engineering 16, vol. 16, pp. 359–368, 1976. View at: Google Scholar
 L. D'Amore, G. Laccetti, and A. Murli, “An implementation of a Fourier series method for the numerical inversion of the Laplace transform,” ACM Transactions on Mathematical Software, vol. 25, no. 3, pp. 279–305, 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Bisheng Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.