Research on Temperature Distribution in Wellbore of Deep Well Drilling ()
1. Introduction
With the increase of world energy demand and the development of petroleum industry, deep and ultra-deep wells have become an important way of oil-gas development [1] . In our country, deep and ultra-deep wells with depths of more than 4500 m to 7000 m are common in Xinjiang and Sichuan, especially in Sichuan and Chongqing, where there are a large number of deep and ultra-deep gas wells in the HTHP (high temperature and high pressure) gas fields [2] . Generally speaking, the deeper the well, the higher the downhole pressure and temperature, and the greater the technical difficulties. In the past, deep exploration 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.
2. Model Development
2.1. 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:
(1)
(2)
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 transfer coefficient,
is the heat generated by mechanical energy loss of rotary drilling tool, pressure drop of bit nozzle or fluid circulating pressure loss in per volume,
is the heat transfer quantity in per volume between the tube
Figure 1. Temperature transfer in deep well.
string and the annular space, or between the annular space and the formation.
2.2. 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 calculated directly. HAS solves the equations to get the temperature distribution in deep well using full implicit difference and Gauss-Seidel iteration methods.
To calculate borehole pressure loss accurately, HAS uses different mud rheological modes (powerlaw, binham, carson and H-B) and a variety of rheological parameter calculation and optimization methods to descript the mud rheological properties. The effects of drill pipe rotation, eccentricity, joints and the downhole motor pressure drop are also taken into account as well as the different hole cleaning conditions, which are the flow with single-phase fluid, the flow with full suspension of cuttings, the flow with general cleaning condition and the flow with poor cleaning condition in the calculation of borehole pressure loss, based on the advanced three-layer cuttings dynamic transport theoretical model of ERW, empirical calculation model of cuttings transport in whole section of ERW under multiple working conditions [13] .
Coupling the above model of pressure loss and calculation model [14] of heat transfer in deep well and ultra-deep well [15] , HAS can calculate and analyze the wellbore temperature distribution during drilling whether to consider the effects of temperature and pressure on mud performance, including mud density and mud rheological properties. Thus, the calculation results of wellbore pressure can be corrected to improve the accuracy of calculation and design of deep wells, HTHP wells and deep water wells. The calculation results are very accurate. A large number of field applications show that the calculation error can be controlled within 5%.
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.
3. Case Study
3.1. 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/cm3, 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.
3.2. Calculation Result
Figure 3 is the variation of temperature distribution in drill string and in annular
Figure 2. The temperature analysis module of HAS.
Figure 3. Temperature variation in drill string (left) and annular (right).
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.
4. 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.
5. Conclusions
1) Based on the theory of fluid mechanics and heat transfer, considering the characteristics of deep well drilling, a wellbore temperature model was proposed in drilling of deep wells. Using full implicit difference and Gauss-Seidel iteration
Figure 4. Annular temperature vs. mud inlet temperature (left) and geothermal gradient (right).
Figure 5. Annular temperature vs. mud flow rate (left) and mud specific heat (right).
Figure 6. Annular temperature vs. mud thermal conductivity (left) and mud density (right).
methods, the model was solved and the distribution of wellbore temperature can be calculated and analyzed.
2) With the progress of drilling, the mud temperature in wellbore gradually decreases from the geothermal temperature to a stable temperature. The mud 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.
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.