Regular Article
Establishment and analysis of temperature field of riserless mud recovery system
State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University (SWPU), Chengdu, Sichuan, China
^{*} Corresponding author: swpivip@163.com
Received:
15
September
2018
Accepted:
11
December
2018
Due to the drill pipe in the Riserless Mud Recovery (RMR) system is exposed to seawater, so the characteristics of the temperature changing are very different from the conventional offshore drilling. Considering temperature is an important factor affecting the annulus pressure, it is necessary to study the variation law of the temperature field of the RMR. In this paper, according to the physical process of the heat transfer in RMR, the mathematical model of the temperature field is established. The Computational Fluid Dynamics (CFD) software is used to simulate the temperature distribution in the drill pipe and the annulus, so that the law of the temperature changing can be observed more intuitively. In order to be more aware of the influencing factors of the temperature field changing, this paper analyzes the influence mechanism of the different discharge capacity and the different injection temperature on temperature changing. Moreover, this paper also analyzes the influence of the annulus temperature on the annulus pressure, which provides a theoretical basis for the well control of RMR.
© J. Zhang et al., published by IFP Energies nouvelles, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Riserless Mud Recovery (RMR) system is an advanced offshore drilling technology developed by Norwegian AGR Company (Stave et al., 2008). Compared with the conventional offshore drilling, RMR adopts the dualgradient control principle to control the annulus pressure more precisely, so as to effectively solve problems such as narrow mud density window, shallow gas and shallow flow that the offshore drilling engineering has always been facing (Michael and Michael, 2001; Myers, 2008; Smith et al., 2010; Stave et al., 2008). Because RMR abandons the riser used in the traditional offshore drilling and uses a seabed pump as a device for lifting drilling fluid, it reduces the cost of drilling and the demand for drilling platform (Claudey et al., 2016; Gao et al., 2009; Hannegan and Stave, 2006; Stave et al., 2014). From the above two points, RMR is suitable for drilling in deep water and ultradeep water, which is in line with the development trend of the current offshore drilling engineering.
As shown in Figure 1, the RMR consists of three modules: the suction module, the subsea pump and the return line (Cohen et al., 2010). The suction module is used to collect the mud returned from the annulus and seal the wellhead (Scanlon and Medeiros, 2012). The outlet of the suction module is connected to the inlet of the seabed pump through a hose. The subsea pump belongs to the disc pump and has been developed to the third generation (Alford et al., 2005; Beek et al., 2016). Its main function is to provide power for the lifting of the mud and by adjusting its pump speed, the pressure acting on the wellhead is equal to the static pressure of seawater at this depth, thus achieving dualgradient. The return line is the only channel for mud to return to the platform from the seabed (Wang, 2013). The pipe size and arrangement mode of the return line will greatly affect the lifting efficiency of mud (Stave et al., 2005; Thorogood et al., 2007).
Fig. 1
Composition of the RMR system. 
Because the drill pipe in RMR is directly exposed to seawater, the characteristics of the temperature changing of the wellbore are greatly different from the conventional offshore drilling due to the variation of the seawater temperature gradient. Temperature is an important factor affecting the annulus pressure. In order to make RMR successfully applied in the deep water and the ultradeep water, and achieve higher precision of the annulus pressure control, it is necessary to study the temperature variation characteristics in RMR in combination with the characteristics of the temperature changing of the marine environment and explore the influence mechanism of the temperature changing.
Since RMR is an emerging drilling technology, its related theory is not very mature so far, which limits the development of RMR to deep water and ultra deep water. Only Norwegian AGR Company has been successful in implementing it commercially, but only in the shallow water (Eikemo, 2015; Ziegler et al., 2013). Currently, RMR is not widely used in deep water and ultradeep water worldwide (Chen et al., 2007; Peyton et al., 2013). In order to strengthen the development of oil and gas resources in deep water and ultra deep water, it is necessary to carry out a deeper theoretical research based on the technical characteristics and advantages of RMR, which is crucial for the development of the offshore oil and gas engineering.
2 Mathematical model
2.1 Basic assumed conditions
Since the drill pipe and the annulus are symmetric models, they can be simplified to a twodimensional planar model for calculation. A drilling fluid control body of length dx is taken in the drill pipe and the annulus respectively, and the flow direction of drilling fluid is set to positive direction, and the following assumed conditions are made:

The temperature in any section of drill pipe and annulus perpendicular to the flow direction is uniform.

Ignore the heat conduction along the flow direction.

All physical properties are constants.

No insulation layer is set on the outer wall of drill pipe.
2.2 Mathematical model in drill pipe
The physical model of heat transfer process of drill pipe exposed to seawater is shown in Figure 2.
Fig. 2
Heat transfer model of drill pipe. 
The heat injected from the upper surface of the control body:(1)
The heat that flows from the lower surface of the control body:(2)
Heat transferred from the drilling fluid in the drill pipe to seawater:(3)
Because the law of seawater temperature decline in different sea areas is different, this paper uses the law of seawater temperature change in the South China Sea as the calculation method of seawater temperature (Vaid and Polito, 2016; Wang et al., 2010):(4)
According to the law of conservation of energy (Φ_{in} = Φ_{out}), the temperature profile of the drilling fluid in the drill pipe can be obtained from the above four equations:(5)
Combined with boundary condition t_{ x = 0} = t_{i}, the calculation equation of the temperature of the drilling fluid in drill pipe can be obtained:(6) (7) (8) (9) (10)where t_{0} is the temperature of sea surface, °C; t_{i} is the injection temperature of drilling fluid, °C; q_{m1} is the mass flow of drilling fluid through control body, kg/s; c_{p} is the specific heat capacity of drilling fluid, J/(kg °C); k_{1} is the heat transfer coefficient between drilling fluid in drill pipe and seawater, W/(m^{2} °C); d_{1} is the inner diameter of drill pipe, m; t_{ x } is the temperature of drilling fluid at x, °C; t_{ x+dx } is the temperature of drilling fluid at x + dx, °C.
2.3 Mathematical model in annulus
The physical model of heat transfer process of annulus is shown in Figure 3.
Fig. 3
Heat transfer model of annulus. 
The heat injected from the lower surface of the control body:(11)
The heat transfer from the formation to the drilling fluid in the annulus:(12)
The temperature of the formation at depth x:(13)
The heat that flows from the upper surface of the control body:(14)
The heat transfer from the drilling fluid in annulus to the drilling fluid in drill pipe:(15)
According to the law of conservation of energy, and combined with the boundary condition t_{ x=0} = t_{u}, the calculation equation of the temperature of drilling fluid in annulus can be obtained:(16) (17) (18) (19) (20)where k_{2} is the heat transfer coefficient between the formation and the drilling fluid in annulus, W/(m^{2} °C); k_{3} is the heat transfer coefficient between the drilling fluid in annulus and the drilling fluid in drill pipe, W/(m^{2} °C); t_{u} is the temperature of the drilling fluid at the bottom of well, °C; t_{e} is the temperature of the drilling fluid in drill pipe, °C; t_{s} is the surface temperature of formation, °C; m is geothermal gradient, °C/100 m; q_{m2} is the mass flow of drilling fluid through control body, kg/s; d_{3} is the diameter of annulus, m; d_{2} is the outer diameter of drill pipe, m.
3 CFD analysis
In order to observe the changing of temperature field intuitively, and to provide a certain basis for the verification of the mathematical model, this paper uses Computational Fluid Dynamics (CFD) software to carry out numerical simulation of some well sections. The model is scaled to a certain scale in this paper, but the heat exchange time is still the same as the actual situation. The basic data used in simulation is from a real well in the South China Sea (Gao et al., 2013).
3.1 CFD analysis for drill pipe
The temperature changing of the drilling fluid in drill pipe at water depth 0–1000 m was simulated. The basic parameters are shown in Table 1.
Input parameters for CFD.
The finite element model of drill pipe is established, as shown in Figure 4. Considering that the flow of the drilling fluid in drill pipe is turbulent for most of the time, it is necessary to add boundary layer near the wall when establishing the finite element model.
Fig. 4
The finite element model of drill pipe. 
In order to make the observation more clearly, this paper simulates the temperature changing of the drilling fluid in drill pipe at water depth 0–500 m and 500–1000 m respectively. The results of the CFD analysis are shown in Figures 5 and 6.
Fig. 5
CFD analysis results of 0–500 m. 
Fig. 6
CFD analysis results of 500–1000 m. 
As can be seen from Figures 5 and 6, the injection temperature of the drilling fluid at 0 m is 20 °C, the temperature at 500 m is about 14 °C, and the temperature at 1000 m is about 10 °C. At the entire stage of 0–1000 m, the temperature of drilling fluid is reduced by about 50% due to the influence of the temperature gradient of seawater. From the results of CFD analysis, the temperature distribution of the drilling fluid in drill pipe is largely affected by seawater temperature when the insulation layer is not installed.
3.2 CFD analysis for annulus
The temperature changing of the drilling fluid in annulus at formation depth 5000–6000 m and 3500–4000 m was simulated. The basic parameters are shown in Table 2.
Input parameters for CFD.
The finite element model of annulus is established, as shown in Figure 7. The flow of the drilling fluid in annulus is also in a turbulent state for most of the time. Therefore, it is necessary to set a boundary layer near the wall of the finite element model of annulus. And because annulus is a symmetric figure, so it can be simulated in half.
Fig. 7
The finite element model of annulus. 
The results of the CFD analysis are shown in Figures 8 and 9.
Fig. 8
CFD analysis results of 5500–6000 m. 
Fig. 9
CFD analysis results of 3500–4000 m. 
As can be seen from Figures 8 and 9, the temperature of the drilling fluid at formation depth 6000 m is about 63 °C, and it is about 78 °C at 5500 m, which is about 15 °C higher. The temperature of the drilling fluid at formation depth 4000 m is about 92 °C, and it is about 83 °C at 3500 m, which is about 9 °C lower. Therefore, when the drilling fluid returns in annulus, the temperature changing trend of drilling fluid is first increased and then decreased.
4 Case study
In this section, based on the basic data of a well in the South China Sea (Gao et al., 2013), the mathematical model is used for calculation, and the results are compared with the results of CFD analysis, so as to verify the feasibility of the mathematical model. The temperature variation characteristics of the drilling fluid in drill pipe and annulus under the influence of the different discharge capacity and different injection temperature are analyzed. The effect of temperature on annulus pressure is also analyzed in this section. The basic data of the well is shown in Table 3.
Basic data of a well in the South China Sea.
4.1 Influence of different discharge capacity
When the injection temperature of drilling fluid is 20 °C, the discharge capacity is set to 35 L/s, 45 L/s, 55 L/s and 65 L/s. The results of calculation are shown in Figures 10 and 11.
Fig. 10
Calculation results of temperature value in drill pipe in different discharge capacity. 
Fig. 11
Calculation results of temperature value in annulus in different discharge capacity. 
It can be seen from Figure 10 that when the injection temperature of drilling fluid is 20 °C and the discharge capacity is 35 L/s, the temperature of the drilling fluid in drill pipe at seawater depth 500 m is about 13 °C, and about 11 °C at 1000 m. As can be seen from Figure 11, the temperature of the drilling fluid at formation depth 6000 m is about 62 °C, about 78 °C at 5500 m, about 91 °C at 4000 m, and about 85 °C at 3500 m. This result is basically consistent with the result of CFD analysis. Therefore, the mathematical model of the temperature field established in this paper has certain accuracy.
As can be seen from Figures 10 and 11, as the discharge capacity increases, the temperature of the drilling fluid in drill pipe and annulus will increase. This is because when the displacement is increased, the heat exchange time between the drilling fluid and the seawater is shortened, so that the temperature decrease of drilling fluid is reduced, and the temperature at the bottom of well is higher. Therefore, in order to maintain the temperature of the drilling fluid, the discharge capacity can be increased to some extent in combination with the actual situation.
4.2 Influence of different injection temperature
When the discharge capacity is 35 L/s, the injection temperature of drilling fluid is set to 20 °C, 25 °C, 30 °C and 35 °C. The results of calculation are shown in Figures 12 and 13.
Fig. 12
Calculation results of temperature value in drill pipe in different injection temperature. 
Fig. 13
Calculation results of temperature value in annulus in different injection temperature. 
As can be seen from Figures 12 and 13, the higher the injection temperature of drilling fluid is, the higher the temperature of drilling fluid at the same water depth and the same formation depth is. This is because when the drilling fluid has the same properties and the discharge capacity is the same, when the injection temperature of drilling fluid increases, the heat Φ_{in1} and Φ_{in2} entering the control body increase. Combined with equations (1) and (11), it can be seen that as Φ_{in1} and Φ_{in2} increase, the temperature of the drilling fluid at the same water depth and formation depth will increase.
4.3 Analysis of variation characteristics
It can be seen from Figures 10 to 12 that when the water depth is less than 1000 m, as the water depth increases, the temperature of the drilling fluid in drill pipe becomes lower and lower, and the extent of decrease is larger. This is because the temperature of seawater at this stage is getting lower and lower as the water depth increases, and the amount of heat exchange between the drilling fluid and the seawater is increasing. However, when the water depth is greater than 1000 m, the temperature of drilling fluid remains substantially constant. This is because the temperature of seawater at this stage reaches the constant temperature layer (Vaid and Polito, 2016), and the temperature of seawater no longer decreases with the increase of water depth. At this time, the heat exchange between the drilling fluid and the seawater is extremely small, so the temperature of drilling fluid is basically stable.
It can be seen from Figures 11 to 13 that when the drilling fluid returns in annulus, the temperature of drilling fluid will rise first, but the increase will gradually become slower. Subsequently, the temperature will continue to decrease. When the drilling fluid returns to a position about 1500 m away from the wellhead, the temperature reduction of drilling fluid is intensified. This is because when the drilling fluid just begins to return in the annulus, Φ_{in3} > Φ_{k2} at this time, so the temperature of drilling fluid will increase. However, when the drilling fluid continues to rise, Φ_{k2} gradually increases, and Φ_{in3} gradually decreases, so the temperature rise of drilling fluid gradually becomes slower. When Φ_{k2} is increased to be larger than Φ_{in3}, the temperature of the drilling fluid in annulus begins to decrease. When the drilling fluid continues to flow upward, Φ_{k2} remains substantially unchanged, and Φ_{in3} decreases further, so the decrease in the temperature of the drilling fluid is further increased.
4.4 Effect of temperature on annulus pressure
The annulus pressure of the RMR consists of seabed pump inlet pressure, circulating liquid column pressure and cyclic pressure loss. It can be obtained as:(21)where P_{ann} is the pressure of annulus, MPa; P_{inlet} is the inlet pressure of subsea pump, MPa; ECD is the Equivalent Circulating Density of drilling fluid, g/cm^{3}; h_{L} is the well depth, m; ∆P_{f} is the cyclic pressure loss, MPa.
The ECD can be calculated via (Dokhani et al., 2016; Lin et al., 2016):(22)where ρ_{w} is the density of seawater, g/cm^{3}; h_{w} is the seawater depth, m; ρ_{m} is the static density of drilling fluid, g/cm^{3}; ρ_{s} is the density of cuttings, g/cm^{3}; C_{a} is cuttings concentration in annulus, dimensionless.
The static density of drilling fluid can be obtained as (Vaid and Polito, 2016):(23)where ρ_{o} is the drilling fluid density on the platform, g/cm^{3}; C_{T} is the thermal expansion coefficient, dimensionless; C_{P} is the elastic compression coefficient, dimensionless; ∆T is the temperature difference between the drilling fluid at a certain depth and the platform, °C; ∆P is the pressure difference between the drilling fluid at a certain depth and the platform, MPa.
Since the main research content of this paper is the variation law of the temperature field of the RMR, the influence of pressure change on the drilling fluid density is not considered in this section. Combined with equations (22) and (23), and the calculation results of ECD results shown in Figure 14 can be obtained.
Fig. 14
Calculation results of ECD in annulus. 
As can be seen from Figure 14, the ECD in annulus decreases first and then increases. The main reason is that the temperature of drilling fluid in annulus rises first and then decreases, which causes the volume of drilling fluid to expand first and then shrink. According to ρ = m/V, the static density of drilling fluid first decreases and then increases. According to equation (22), the same trend of ECD in the annulus can be obtained.
Based on equation (21), the annulus pressure is calculated. The calculation results are shown in Figure 15.
Fig. 15
Calculation results of annulus pressure. 
It can be seen from Figure 15 that the annulus pressure is still within the safe pressure window at a well depth 2500 m. When the well depth exceeds 2500 m, the annulus pressure is lower than the formation pressure, and the formation fluid is easy to invade the wellbore to cause the kick and blowout.
Under such conditions, the discharge capacity should be appropriately increased in combination with the actual situation to reduce the heat exchange time between the drilling fluid in the annulus and the formation, effectively suppressing the volume expansion of the drilling fluid, and making the value of ECD increases, the annulus pressure rises, allowing more well sections to be within the safe pressure window, further simplifying the casing program and reducing drilling costs.
5 Conclusion

If no insulation layer is set on the outer wall of drill pipe, the temperature of the drilling fluid in drill pipe will continue to decrease and is basically stable when the water depth is 1000 m. This will have a great impact on the control of wellbore pressure. Therefore, in the actual drilling, it should be considered to set an insulation layer on the outer wall of drill pipe.

The variation of temperature of the drilling fluid in annulus is mainly affected by the two values of Φ_{in3} and Φ_{k2}. As the two values constantly change, the temperature of the drilling fluid in annulus first increases and then decreases and the magnitude of the changing is large. Therefore, the calculation of annulus pressure needs to take into account the changing characteristics of drilling fluid temperature.

Under the condition of same injection temperature, the temperature of the drilling fluid at same water depth and same formation depth will increase with the increase of discharge capacity. This is very helpful in maintaining the temperature of drilling fluid. Therefore, it is necessary to increase the discharge capacity to a certain extent according to the actual situation when drilling.

With the same discharge capacity, the temperature of the drilling fluid at same water depth and same formation depth will increase with the increase of injection temperature. This also helps to keep the temperature of drilling fluid. Therefore, it is necessary to increase the injection temperature to a certain extent according to the actual situation when drilling.

According to the actual situation, the heat exchange between the drilling fluid in the annulus and the formation should be controlled by adjusting the discharge capacity, so as to effectively control the ECD value, so that more well sections are in the safety pressure window.
Acknowledgments
The financial support from the Natural Science Foundation of China (NSFC) (No. 51274168) and the National Key R&D Program of China (No. 2018YFC0310202) is gratefully acknowledged.
References
 Alford S.E, Asko A., Campbell M., Aston M.S., Kvalvaag E. (2005) Silicatebased fluid, mud recovery system combine to stabilize surface formations of Azeri Wells, SPE/IADC Drilling Conference, 23–25 February, Amsterdam, The Netherlands. [Google Scholar]
 Beek P.J.G.V., Pereboom H.P., Slot H.J. (2016) Evaluating vibration performance of a subsea pump module by fullscale testing and numerical modeling, ASME International Conference on Ocean, pp. V005T04A059. [Google Scholar]
 Chen G.M., Yin Z.M., Xu L.B., Jiang S.Q. (2007) Review of deepwater dual gradient drilling technology, Pet. Explor. Dev. 34, 2, 246–251. [Google Scholar]
 Claudey E., Maubach C., Ferrari S. (2016) Deepest deployment of riserless dual gradient mud recovery system in drilling operation in the North Sea, SPE Bergen One Day Seminar, 20 April, Bergen, Norway. [Google Scholar]
 Cohen J., Kleppe J., Grns T., Martin T.B. (2010) Gulf of Mexico’s first application of riserless mud recovery for tophole drilling – a case study, Offshore Technology Conference, 3–6 May, Houston, Texas, USA. [Google Scholar]
 Dokhani V., Ma Y., Yu M. (2016) Determination of equivalent circulating density of drilling fluids in deepwater drilling, J. Nat. Gas Sci. Eng. 34, 1096–1105. [Google Scholar]
 Eikemo B. (2015) Method and device for riserless drilling fluid recovery, US14356701. [Google Scholar]
 Gao B., Chen G., Yin Z., Liu S.J. (2009) Deepwater riserless mud recovery drilling technology, Oil Drill. Prod. Technol. 28, 3, 381–390. [Google Scholar]
 Gao D., Sun T., Zhang H., Tang H. (2013) Displacement and hydraulic calculation of the SMD system in ultradeepwater condition, Liquid Fuels Technol. 31, 11, 1196–1205. [Google Scholar]
 Hannegan D., Stave R. (2006) The time has come to develop riserless mud recovery technology’s deepwater capabilities, Drill. Contractor 62, 5, 106–111. [Google Scholar]
 Lin T., Wei C., Zhang Q., Sun T. (2016) Calculation of equivalent circulating density and solids concentration in the annular space when reaming the hole in deepwater drilling, Chem. Technol. Fuels Oils 52, 1, 70–75. [CrossRef] [Google Scholar]
 Michael J., Michael R. (2001) Riserless drilling technique saves time and money by reducing logistics and maximizing borehole stability, SPE Annual Technical Conference and Exhibition, 30 September–3 October, New Orleans, Louisiana. [Google Scholar]
 Myers G. (2008) Ultradeepwater riserless mud circulation with dual gradient drilling, Sci. Drill. 6, 6, 76–89. [CrossRef] [Google Scholar]
 Peyton J., Mcphee A., Eikemo B., Evans H. (2013) World first: Drilling with casing and riserless mud recovery, 6th International Petroleum Technology Conference, Beijing, China. [Google Scholar]
 Scanlon T., Medeiros F. (2012) Enhanced drilling solution for tophole sections on jackup wells with environmentallyimproved method and dualgradient drilling techniques, Offshore Technology Conference, Houston, Texas, USA. [Google Scholar]
 Smith D., Tarr B., Winters W.J., Ziegler R. (2010) Deepwater riserless mud return system for dual gradient tophole drilling, SPE/IADC Managed Pressure Drilling and Underbalanced Operations Conference and Exhibition, 24–25 February, Malaysia. [Google Scholar]
 Stave R., Farestveit R., Hyland S., Rochmann P.O. (2005) Demonstration and qualification of a riserless dual gradient system, Offshore Technology Conference, Houston, Texas, USA. [Google Scholar]
 Stave R., Brainard R.R., Cohen J.C. (2008) Riserless mud recovery moves into deepwater, Offshore 68, 5, 100–104. [Google Scholar]
 Stave R., Nordas P., Fossli B., French C. (2014) Safe and efficient tophole drilling using riserless mud recovery and managed pressure cementing, Offshore Technology Conference Asia, 25–28 March, Kuala Lumpur, Malaysia. [Google Scholar]
 Thorogood J., Rolland N., Brown J., Rolland N.L. (2007) Deployment of a riserless mudrecovery system offshore Sakhalin Island, SPE/IADC Drilling Conference, 20–22 February, Amsterdam. [Google Scholar]
 Vaid B.H., Polito P.S. (2016) Influence of the South China Sea biweekly sea surface temperature on the South China Sea summer monsoon especially during the Indian Ocean dipole, Atmos. Ocean 54, 1, 48–59. [CrossRef] [Google Scholar]
 Wang G. (2013) Mrl lectotype and parametric optimization for deepwater riserless drilling, China Pet. Mach. 48, 6, 97–102. [Google Scholar]
 Wang Z.M., Hao X.N., Wang X.Q., Xue L., Guo X.L. (2010) Numerical simulation on deepwater drilling wellbore temperature and pressure distribution, Liquid Fuels Technol. 4828, 9, 911–919. [Google Scholar]
 Ziegler R., Ashley P., Malt R.F., Stave R. (2013) Successful application of deepwater dual gradient drilling, IADC/SPE Managed Pressure Drilling and Underbalanced Operations Conference and Exhibition, 15–17 October, Moscow. [Google Scholar]
All Tables
All Figures
Fig. 1
Composition of the RMR system. 

In the text 
Fig. 2
Heat transfer model of drill pipe. 

In the text 
Fig. 3
Heat transfer model of annulus. 

In the text 
Fig. 4
The finite element model of drill pipe. 

In the text 
Fig. 5
CFD analysis results of 0–500 m. 

In the text 
Fig. 6
CFD analysis results of 500–1000 m. 

In the text 
Fig. 7
The finite element model of annulus. 

In the text 
Fig. 8
CFD analysis results of 5500–6000 m. 

In the text 
Fig. 9
CFD analysis results of 3500–4000 m. 

In the text 
Fig. 10
Calculation results of temperature value in drill pipe in different discharge capacity. 

In the text 
Fig. 11
Calculation results of temperature value in annulus in different discharge capacity. 

In the text 
Fig. 12
Calculation results of temperature value in drill pipe in different injection temperature. 

In the text 
Fig. 13
Calculation results of temperature value in annulus in different injection temperature. 

In the text 
Fig. 14
Calculation results of ECD in annulus. 

In the text 
Fig. 15
Calculation results of annulus pressure. 

In the text 