Research on Temperature Distribution in Wellbore of Deep Well Drilling

The wellbore temperature has an important effect on design and drilling of deep well. Based on energy conservation equations and actual drilling data of one deep well, the wellbore temperature distribution was simulated and the influence of different parameters on the wellbore temperature was revealed using the software of Hydraulics Analysis System. The results show that, while drilling, the mud temperature in wellbore gradually decreases from the formation temperature to the stable temperature, and it is higher than the mud inlet temperature on ground, the annular temperature is higher than the temperature in drill string, and the bottom hole temperature is higher than the ground temperature. The effect of geothermal gradient on wellbore temperature is great, while the mud density is negligible. The bottom hole temperature increases with the increase of mud inlet temperature, geothermal gradient, mud thermal conductivity and decrease of mud flow rate, mud specific heat and mud density.

focused more on finding oil, and the geothermal gradient was low. In the future, more attention will be paid to deep gas exploration, and the problem of high temperature and high pressure will become prominent. The high temperature and high pressure will affect the mud property in the wellbore [3] [4], which makes the bottom hole pressure calculation have large errors and easily cause drilling accidents. Meanwhile, the excessive temperature also restricts the use of many downhole tools.
In 1962, Ramey proposed the classical method of predicting wellbore fluid temperature distribution [5]. He took the fluid temperature as a function of well depth and time, deduced the calculation formula of wellbore temperature and put forward the concept of comprehensive heat transfer coefficient. In 1967, Willhite gave the classical formula for calculating the total heat transfer coefficient [6]. Later, many scholars proposed new calculation models of wellbore temperature based on Ramey method [7] [8] [9]. In China, in 1996, aiming to solve the problem of heavy oil thermal recovery, Li Zifeng established a prediction model of well temperature distribution in the process of circulation [10], and solved the model by finite volume method. Other scholars subsequently developed some similar or improved models [11] [12]. These studies have shown that the temperature in the wellbore is not exactly the same as the geothermal field and is a dynamic process that requires a dedicated model for calculation.
For this reason, the HAS (Hydraulics Analysis System) software developed by Chongqing University of Science and Technology was imported and the temperature distribution in wellbore of deep well drilling was analyzed. The analysis result can be useful for the design of deep well drilling.

Model Equations
The well structure and heat transfer diagram in deep well is shown in Figure 1.
The energy conservation equations of various media in the wellbore system [11], such as drilling fluid, drill string, annulus and formation, can be expressed as follows: In the above equations, t is time, z is the axial position along wellbore, r is the radial position in borehole, T is the temperature in drill string or annular, ρ is the density of drilling fluid, c is the specific heat of drilling fluid, v is the drilling fluid flow velocity, T is the temperature of drilling fluid, k is the convective heat  string and the annular space, or between the annular space and the formation.

Model Solver
Apparently the temperature distribution is related not only with heat conduction but also with fluid pressure transfer. Due to the complexity the model cannot be Coupling the above model of pressure loss and calculation model [14]  Based on these models, HAS can be also used for hydraulics design before drilling, hydraulics monitoring during drilling and hydraulics analysis and optimization before or after drilling. It can be used for hydraulics analysis, monitoring, design and optimization of complex structural wells drilling such as conventional directional wells, horizontal wells, sidetrack wells, small wellbore wells, deep wells, ultra-deep wells, HTHP wells, deepwater wells, etc.
The temperature analysis module of HAS is shown in Figure 2.

Basic Data
A vertical deep well of 7000m was drilled with Φ215.9 mm bit, Φ139 mm drill pipe. The oil-based mud density was 1.75 g/cm 3 , mud inlet temperature was 27˚C, mud flow rate was 30 L/s, nozzle combination was Φ22 mm + Φ22 mm + Φ33 mm, the specific heat of mud, steel and formation were 1675, 400 and 837 J/(Kg·˚C) respectively, and the thermal conductivity was 1.73, 44 and 2.25 W/(m·˚C) respectively, the geothermal gradient was 3˚C/100m.  with time when drilling to 7000 m. It can be seen from the figures that with the progress of drilling, the mud temperature in wellbore gradually decreases from the geothermal temperature to a stable temperature. After 16 hours, the temperature change begins to slow down whether in annular or drill string. The bottom hole temperature is about 118˚C, which is 91˚C higher than the mud inlet temperature, and the annular temperature is greater than the temperature in drill string. If the geothermal gradient is higher, for example 5˚C/100m in Southwest oil and gas field in China, the bottom hole temperature will become a big problem for drilling.

Discussion
In general, the temperature in annular is always higher than drill string. Therefore, the sensitivity analysis is performed with the annular temperature distribution when drilling to 7000 m after 16 hours. The factors affecting the temperature distribution include mud inlet temperature (tin), geothermal gradient (tg), mud specific heat (Cf), mud thermal conductivity (Kf), mud density (ρf) and mud flow rate (Qf). The influence of each factor on the annulus temperature is shown in Figures 4-6. It can be seen from these figures that the bottom hole temperature increases with the increase of mud inlet temperature, geothermal gradient and mud thermal conductivity, and decreases with the increase of mud flow rate, mud specific heat and mud density. Geothermal gradient and mud specific heat have a greater impact on the annulus temperature field than mud thermal conductivity, mud inlet temperature and mud flow rate, while mud density has negligible impact.    2) With the progress of drilling, the mud temperature in wellbore gradually decreases from the geothermal temperature to a stable temperature. The mud B. Y. Xia Open Journal of Applied Sciences temperature in wellbore is higher than the mud inlet temperature on the ground, the annulus temperature is higher than the temperature in the drill string, and the bottom hole temperature is higher than the ground temperature.

Conclusions
3) The factors affecting the temperature field include mud inlet temperature, geothermal gradient, mud specific heat, mud thermal conductivity, mud density and mud flow rate, etc. The bottom hole temperature increases with the increase of mud inlet temperature, geothermal gradient and mud thermal conductivity, and decreases with the increase of mud flow rate, mud specific heat and mud density.
4) Geothermal gradient and mud specific heat have a greater impact on the annulus temperature field than mud thermal conductivity, mud inlet temperature and mud flow rate, while mud density has negligible impact.

Conflicts of Interest
The author declares no conflicts of interest regarding the publication of this paper.