Wellbore Temperature and Pressure Calculation of Offshore Gas Well Based on Gas–Liquid Separated Flow Model

: Compared with land wells, the production environment and reservoir depth of offshore oil and gas wells are more complex and shallower. Further, HPHT production ﬂuid there will produce strong temperature and pressure disturbance that affects the wellbore, which easily generates wellbore safety problems, such as wellhead growth and leakage caused by the incompatible deformation of casing and cement sheath. Therefore, obtaining an accurate wellbore temperature and pressure ﬁeld is the key to implementing a wellbore safety assessment. Based on the gas–liquid two-phase separated method, this paper established an improved calculation model of wellbore temperature and pressure ﬁeld for offshore HPHT wells. This model also takes into account the heat transfer environment characteristics of “formation-seawater-air” and the inﬂuence of well structure. Compared with the measured data of the case well, the error of temperature and pressure calculation results of the improved model are only 0.87% and 2.46%. Further, its calculation accuracy is greatly improved compared to that of the traditional gas–liquid homogeneous ﬂow calculation model. Based on this model, the inﬂuencing factors of wellbore temperature and pressure in offshore gas wells are analyzed. The results show that forced convection heat exchange between seawater–air and wellbore is stronger than that between wellbore and formation. Reducing the gas–liquid ratio of the product can effectively reduce wellbore temperature and increase wellbore pressure. The gas production has a signiﬁcant impact on the wellbore temperature. When the gas production rises from 10 × 10 4 · m 3 /d to 60 × 10 4 · m 3 /d, the wellhead temperature rises from 63 ◦ C to 99 ◦ C. However, due to the mutual inﬂuence of friction pressure drop and hydrostatic pressure drop, wellbore pressure increases ﬁrst and then decreases with the increase in gas production. The improved model can provide a more accurate estimate of the time to reach the rated wellhead temperature. Meanwhile, this model displays accurate theoretical support for the rational formulation of the production plan after the well opening, so as to avoid excessive restrictions on the initial production rate.


Introduction
The offshore oil and gas field will become the main arena in which to improve domestic oil and gas production capacity and ensure national energy security. However, the harsh and complex offshore environment and operating conditions are the realistic challenges that offshore oil and gas development and production must face. The complex geological genesis of offshore oil and gas results in some HPHT reservoirs. For example, there are HPHT wells with an absolute reservoir temperature of more than 200 • C and a formation pressure coefficient up to 2.3 in Yingqiong Basin and the Ledong area [1,2]. Compared with Tarim Basin, where reservoir temperature and pressure are comparable despite buried depths of more than 7500 m [3,4], the buried depth of offshore oil and gas is usually only about 3000 m. The HPHT fluid in offshore oil and gas reservoirs will reach the wellhead in a shorter time, causing a stronger thermal-pressure coupling effect on the entire wellbore, which leads to difficulties in operation and production. Further there are even safety problems such as wellhead growth, well leakage, and well integrity damage [5].
In 1962, Ramey established a temperature distribution model for single-phase fluid flowing in the wellbore, which transformed heat transfer in the wellbore into steady state heat transfer and formation heat transfer into unsteady state heat transfer. This model is still widely used today [6]. Shiu and Beggs developed semi-empirical expressions for the relaxation distance parameters in the Ramey model based on the analysis of field data and gas and fluid flow characteristics, which improved the utility of the model [7]. Hasan and Kabir transformed single-phase flow into two-phase flow based on the Ramey model and established a semi-empirical analytical solution model based on gas-liquid two-phase flow by using the calculation method of homogeneous-phase flow [8]. Spindler obtained the exact analytical solution of the Hasan model by establishing dimensionless intermediate variables, which improved the accuracy of calculation [9]. Wu and Xu established the pressure and temperature coupling model of HPHT gas wells under single-phase gas flow and analyzed the influence of geothermal gradient on wellbore temperature field [10,11]. Martínez found that the change of fluid viscosity would have a significant impact on the fluid flow pattern and pressure drop, and the change of wellbore temperature and pressure would also have a great impact on the wellbore fluid flow pattern, and then change the heat-transfer ability of the wellbore fluid [12]. Based on this, Yin and Gao, combined with the correction of compressibility and gas-liquid mixing density, established a pseudo single-phase flow wellbore-temperature field model considering the change of wellbore fluid flow pattern [13,14]. Hou studied the influence of gas-liquid two-phase flow pattern change on fluid heat-transfer ability through experiments [15]. According to the change of fluid flow characteristics after well shut-in, Song established a transient model of the temperature field of ocean wellbore under the condition of well shut-in and studied the heattransfer mechanism during well shut-in [16]. Ma simulated the heat-transfer experiment of gas-liquid two-phase flow and gave the calculation formula of convective heat transfer between seawater and string [17]. Sivaramkrishnan developed a simulation algorithm for an unsteady state non-isothermal two-phase wellbore model to predict the downward flow of a wet steam based on a drift-flux model. However, the frictional multiplier of gas-liquid two-phase flow is not considered [18].
The above studies mainly focus on the influence of physical parameters of production fluid and the change of external heat-transfer environment on wellbore temperature and pressure. There are few studies on gas-liquid two-phase flow in wellbore, and these studies usually use the homogeneous flow method, treating gas-liquid two-phase flow as quasisingle flow, in order to calculate wellbore temperature and pressure. The homogenous flow model assumes that the two phases are well mixed and move with identical velocities. However, the separated flow model assumes the two phases to flow separately and share a definite and continuous interface. This method is more in line with the actual situation of oil and gas production. Therefore, this paper establishes a new gas-liquid two-phase flow wellbore temperature and pressure mathematical model by using the calculation method of gas-liquid-phase flow. Further, this model also takes into account the heat-transfer environment characteristics of "formation-seawater-air" and the influence of well structure. Figure 1 shows the wellbore structure. According to different wellbore structures, regions are divided, and a heat-transfer model between the wellbore and the external environment is established. The temperature field and pressure-field coupling should be calculated first, followed by multiple annular temperatures in semi-steady mode. Due to fluid migration from the bottom-up along the tubing, fluid temperature is significantly 1.

Calculation Model of Wellbore Temperature Field
Assume that gas-liquid two-phase fluid flows in a one-dimensional steady direction in the tubing. The fluid flow is turbulent flow, and there is no phase change.

2.
Tubing is concentric with the riser, casing, formation, and cement sheath. Further, the cement sheath is well-bonded to the formation and casing. 3.
The temperature profile of the formation and seawater is approximately linear. Further, the air temperature is considered the same as the temperature of the sea level.

4.
When the fluid flows in the tubing, only the heat-transfer process occurs, and no energy exchange with the outside world occurs.

5.
Radial heat transfer occurs between the fluid in the tubing and the environment, and axial heat transfer in the flow direction is ignored. Figure 1 shows the wellbore structure. According to different wellbore structures, regions are divided, and a heat-transfer model between the wellbore and the external environment is established. The temperature field and pressure-field coupling should be calculated first, followed by multiple annular temperatures in semi-steady mode. Due to fluid migration from the bottom-up along the tubing, fluid temperature is significantly higher than environment temperature in the same plane, resulting in heat transfer between them. Before establishing the calculation model, the following assumptions must be made:

Calculation Model of Wellbore Temperature Field
1. Assume that gas-liquid two-phase fluid flows in a one-dimensional steady direction in the tubing. The fluid flow is turbulent flow, and there is no phase change. 2. Tubing is concentric with the riser, casing, formation, and cement sheath. Further, the cement sheath is well-bonded to the formation and casing. 3. The temperature profile of the formation and seawater is approximately linear. Further, the air temperature is considered the same as the temperature of the sea level. 4. When the fluid flows in the tubing, only the heat-transfer process occurs, and no energy exchange with the outside world occurs. 5. Radial heat transfer occurs between the fluid in the tubing and the environment, and axial heat transfer in the flow direction is ignored.  Figure 2 depicts wellbore heat transfer. Taking a micro-segment of the tubing for example, the arrow indicates the flow direction. According to the law of energy conservation, the energy expression of this micro-segment is obtained:

Tubing Heat-Transfer Model
where H is the fluid of enthalpy, J; g is gravitational acceleration, m/s 2 ; v is fluid velocity, m/s; and Q is heat loss in a micro-segment, J.
where H is the fluid of enthalpy, J; g is gravitational acceleration, 2 m s ; v is fluid velocity, m s ; and Q is heat loss in a micro-segment, J. The differential formula can be obtained by deriving Equation (2): where q is the heat loss of the micro-segment in unit time, ( ) J m s ⋅ ; h is the specific enthalpy of gas, J kg ; and θ is the well inclination angle, °.
Fluid flows axially from bottom to wellhead, so the loss of the wellbore micro-segment in unit time-namely, heat transmitted from fluid to well wall-is [19,20]: where po r is the outer diameter of the tubing, m; T U is the total heat-transfer coefficient of the wellbore, ; f T is the fluid temperature of the tubing, o C ; to T is wellbore temperature, o C ; and ω is mass flow rate, kg/s. The total thermal resistance of radial heat transfer in the wellbore is composed of the thermal resistance of fluid in tubing (Rf), string thermal resistance (Rp and Rc), annulus fluid thermal resistance (Ran), and cement sheath thermal resistance (Rcem) in a series. The specific calculation equation of the total heat transfer coefficient of the wellbore, according to the heat-transfer principle of thick wall cylinder, is [21]: The differential formula can be obtained by deriving Equation (2): where q is the heat loss of the micro-segment in unit time, J/(m · s); h is the specific enthalpy of gas, J/kg; and θ is the well inclination angle, • . Fluid heat flowing into the tubing (Q in ) equals the sum of heat flowing out of the micro-segment per unit time (Q out ) and radial heat loss (Q loss ).
Fluid flows axially from bottom to wellhead, so the loss of the wellbore micro-segment in unit time-namely, heat transmitted from fluid to well wall-is [19,20]: where r po is the outer diameter of the tubing, m; U T is the total heat-transfer coefficient of the wellbore, J/ s · m 2 · • C ; T f is the fluid temperature of the tubing, • C; T to is wellbore temperature, • C; and ω is mass flow rate, kg/s. The total thermal resistance of radial heat transfer in the wellbore is composed of the thermal resistance of fluid in tubing (R f ), string thermal resistance (R p and R c ), annulus fluid thermal resistance (R an ), and cement sheath thermal resistance (R cem ) in a series. The specific calculation equation of the total heat transfer coefficient of the wellbore, according to the heat-transfer principle of thick wall cylinder, is [21]: where R T is the total heat-transfer coefficient of the wellbore, J/ s · m 2 · • C ; k is the number of annulus in the wellbore; h kanc is the convective heat-transfer coefficient of the kth annulus, J/ s · m 2 · • C ; h kanr is the radiative heat-transfer coefficient of the kth annulus, J/ s · m 2 · • C ; r kco is the outer radius of the outer casing of the kth annulus, m; r kci is the inner radius of the outer casing of the kth annulus, m; and k ck is the thermal conductivity of the outer casing of the kth annulus, J/ s · m 2 · • C .
The convection heat-transfer coefficient of the production fluid is [22]: I 0.25 (6) where h L is the liquid convection heat-transfer coefficient of the tubing, J/ s · m 2 · • C ; k L is the thermal conductivity of the tubing liquid, J/(s · m · • C); Re L is the in situ liquidphase Reynolds number of the tubing liquid, dimensionless; Pr is the Prandtl constant, dimensionless; µ w is the dynamic viscosity of the production casing wall, Pa · s; and F p is the fluidity factor, dimensionless. The calculation formula of the in situ Reynolds number of the tubing liquid is: The calculation formula of the fluidity factor is: where ϕ g is gas holdup, dimensionless; F s is the shape factor, dimensionless. The shape factor calculation formula is: The calculation formula of the inclination coefficient is: where σ p is surface tension, N/m. Heat transfer from well-wall to formation in a wellbore micro-segment: where f (t d ) is the dimensionless time function of formation, dimensionless; k d is the thermal conductivity of the formation, J/(s · m · • C); and T d is the formation temperature at any depth of well, • C. According to the heat-transfer model and the transient heat-transfer function improved by Ramey and Hasan in the steady heat-transfer process, the dimensionless time function of the formation can be obtained as follows [23]: where λ D is the thermal diffusion coefficient of the formation, m 2 /s; t D is dimensionless time, dimensionless; and t is the production time, s. When drilling and completing a well on an offshore platform, the seawater section of the mud line is run into a riser to isolate the seawater. Below the mudline, the wellbore temperature gradually decreases with the decrease in formation depth. Above the mud line, the wellbore temperature is affected by the seawater temperature, which drops from the sea level temperature to the temperature at the mud line (the mudline temperature is 0-4 • C).
The heat-transfer process of the riser section wellbore is similar to the heat-transfer model below the mudline, including the heat transfer of the oil casing wall, the heat transfer between the casing wall and the riser wall, and the convective heat transfer between the riser wall and the seawater.
Heat transfer from well wall to seawater in a wellbore micro-segment: where h sea is the convective heat-transfer coefficient between seawater and casing wall, J/ s · m 2 · • C ; T sea is the formation temperature at any depth of well, • C. The convective heat-transfer coefficient of seawater: where Nu sea is the Nusselt coefficient of seawater, dimensionless; k sea is the thermal conductivity of seawater, J/(s · m · • C); and r ro is the outer diameter of the riser, m. The Nusselt relation of the riser under different seawater flow states is as follows [17]: Between the sea level and the offshore platform, the convective heat transfer between the casing and the external air environment is not negligible. The heat transfer in the air section is similar to that in the seawater section, except that the external heat-transfer medium is different. Therefore, the seawater heat-transfer model is still adopted, and the convective heat-transfer coefficient of seawater can be changed into air.
Heat transfer from well wall to air in a wellbore micro-segment: where h air is the convective heat-transfer coefficient between air and casing wall, J/ s · m 2 · • C ; T air is the air temperature, • C. The air convective heat-transfer coefficient: where Nu air is the Nusselt coefficient of the air, dimensionless; k air is the thermal conductivity of air, J/(s · m · • C). The Nusselt relation of the riser under different air-flow states is as follows [24]: Laminar flow : Transitional flow : 1.08 × 10 4 ≤ Gr d ≤ 6.9 × 10 Nu air = 2.9Gr 1 where Ra h is the Rayleigh number of the air, dimensionless.
Radial heat in and out of the micro-segment per unit time can be considered equal for the second contact surface.
Combined with Equations (4), (12), (13), (16) and (20), the well wall temperature in different heat-transfer environments is: The relationship between specific enthalpy and temperature and pressure is obtained from a specific heat capacity at a constant pressure and the Joule-Thomson coefficient, according to the first law of thermodynamics: where C J is the Joule-Thomson coefficient, K/Pa; C p is the specific heat capacity of the fluid at constant pressure, J/(kg · K). The pressure gradient formula of gas-liquid two-phase flow is: where ρ f is the mixed fluid density, kg/m 3 ; f is the friction coefficient, dimensionless; and r pi is the tubing radius, m. The gas-phase flow rate in the tubing is: where B f is the volume coefficient of natural gas, K/MPa; v st is the gas-phase flow rate under standard state, m/s; Z is the deviation factor of natural gas, dimensionless; P f is the pressure of the fluid acting on the tubing, MPa; and Q st is the gas flow under standard state, m 3 /d. The liquid-phase flow rate in the tubing is: where v l is the liquid flow rate, m/s; Q l is the liquid flow of the tubing, m 3 /d.
The actual velocity of gas-liquid two-phase is [25]: Processes 2022, 10, 2043 8 of 18 The gas density in the tubing is: where s g is the gas relative density of the tubing. By the volume-weighted calculation of gas-liquid two-phase density, the mixed fluid density in the tubing is: where ρ l is the liquid density of the tubing, kg/m 3 . Liquid holdup reflects the ratio of liquid to production fluid on the section of tubing. The calculation formula is: where C o is a distribution parameter, dimensionless; v gm is the drift speed, m/s. The calculation formula of the distribution parameters is: The calculation formula of drift speed is [22]: The calculation formula of gas volume flow fraction λ is: Bring Equations (4), (21) and (23) into Equation (2) to obtain: Fluid temperature at the outlet of a micro-segment in formation is obtained by solving the first-order differential Equation (33): Fluid temperature at the outlet of a micro-segment in seawater is obtained using the same solution process: Similarly, fluid temperature at the outlet of a micro-segment in seawater is obtained by: The wellbore is divided into countless micro-segments, and the outlet temperature of each micro-segment is calculated as the inlet temperature of the next section, which is successively calculated until the wellhead fluid temperature is obtained, and finally, the tubing temperature profile is obtained. Figure 3a illustrates the single-layer annulus heat transfer. The heat flowing into the annulus Q ano is equal to the heat flowing out of the annulus per unit of time. Further, the single-layer annulus temperature can be calculated by [26]:

Wellbore Annulus Temperature Calculation
where: The wellbore is divided into countless micro-segments, and the outlet temperatur of each micro-segment is calculated as the inlet temperature of the next section, which i successively calculated until the wellhead fluid temperature is obtained, and finally, th tubing temperature profile is obtained. Figure 3a illustrates the single-layer annulus heat transfer. The heat flowing into th annulus a n o Q is equal to the heat flowing out of the annulus per unit of time. Further the single-layer annulus temperature can be calculated by [26]:    According to Fourier's law of thermal conductivity, calculate the annulus temperature. Then, the heat transfer of the kth annulus can be expressed as:

Wellbore Annulus Temperature Calculation
2πk an r ank h anck T n+1 ank−1 − T n+1 ank ∆r ank + 2πk an r ank+1 h anck+1 T n+1 ank − T n+1 ank+1 ∆r ank+1 = C pan ρ an π∆r 2 ank h anck T n+1 ank − T n ank ∆t Where k an is thermal conductivity of fluid in the annulus, J/(s · m · • C); ∆r ank is the difference between the outer surface radius of the kth annulus and the outer surface radius of the (k + 1)th annulus, m.

Wellbore Pressure Field Calculation Model
In the process of offshore HPHT gas-well exploitation, a small amount of liquid phase such as water and oil and natural gas are extracted to the wellhead at the same time, so the calculation model of separated flow pressure drop can be adopted. Different from the homogeneous flow model, it considers the slip at the gas-liquid interface.
The pressure drop of a flow pipeline mainly consists of three parts: hydrostatic pressure drop caused by the axial flow of pipeline fluid, pressure drop caused by friction of the pipe wall and gas-liquid interface, and pressure drop caused by acceleration. The contribution of each of these sections to the total pressure drop depends on the flow mode, fluid characteristics, pipe orientation, and pipe diameter.
The hydrostatic pressure drop of gas-liquid two-phase flow is [27]: The pressure drop caused by the acceleration of gas-liquid two-phase flow is: β is the correction coefficient of gas holdup, and the calculation formula is: E is the liquid entrainment fraction. Due to the relative motion between the gas and liquid phases, a special phenomenon called entrainment can be observed. The entrainment process is characterized by tiny droplets flowing into a rapidly moving gas phase at the center. Liquid entrainment is due to the obvious shear phenomenon at the gas-liquid interface, which depends on the fluid flow rate, pipe diameter and direction, and the surface tension at the gas-liquid interface. The liquid entrainment fraction E is defined as the ratio of the mass flow rate of droplets entering the gas phase to the total mass flow rate of the liquid phase.
The calculation formula of liquid entrainment fraction is: where We c is the Weber number, dimensionless; ρ c is the fluid core density, kg/m 3 . The calculation formula of the Weber number is: The fluid core density calculation formula is: In this paper, the friction pressure drop is calculated by selecting two-phase frictional multiplier Φ 2 lo , which is summarized by Ghajar and Bhagwat through a large number of experiments. The pressure drop caused by interface friction for gas-liquid two-phase flow is: The calculation formula of two-phase frictional multiplier Φ 2 lo is: where B 1 , B 2 , B 3 , Π 1 , Π 2 , Π 3 , ζ, Bo, and Y are the calculation parameters and correction coefficient of two-phase frictional multiplier Φ 2 lo (for the specific formula, please refer to the literature [21,25]).
For laminar flow and turbulent flow in a circular pipe, this calculation method can be used to calculate the friction coefficient of single-phase string [28]: where e is the absolute roughness of the tubing wall, m.

Model Verification
In order to ensure the reliability of the results of the calculation model established in this paper, a gas well in the China sea is used for an example. This case well is a directional well. The basic parameters are shown in Table 1 and Figure 1, and the numerical calculation process is shown in Figure 4.     Figure 5 shows the comparison of the wellbore temperature and pressure field results of the case wells, which were obtained by different calculation methods. It can be seen from the figure that the results obtained by the two calculation methods are relatively close, and the variation of the wellbore temperature and pressure curve is the same. However, the wellbore data obtained by the homogeneous flow calculation are larger than that of the separated flow calculation. The wellhead temperature and pressure are selected to compare with the measured data, as shown in Figure 6. The comparison between the calculated results of homogeneous flow and the measured data shows that there is an error of 2.33% in the wellhead temperature and an error of 17.45% in the wellhead pressure. Compared with the measured data, the calculated results of the separated flow in this  Figure 5 shows the comparison of the wellbore temperature and pressure field results of the case wells, which were obtained by different calculation methods. It can be seen from the figure that the results obtained by the two calculation methods are relatively close, and the variation of the wellbore temperature and pressure curve is the same. However, the wellbore data obtained by the homogeneous flow calculation are larger than that of the separated flow calculation. The wellhead temperature and pressure are selected to compare with the measured data, as shown in Figure 6. The comparison between the calculated results of homogeneous flow and the measured data shows that there is an error of 2.33% in the wellhead temperature and an error of 17.45% in the wellhead pressure. Compared with the measured data, the calculated results of the separated flow in this paper show an error of 0.87% in wellhead temperature and only 2.46% in wellhead pressure. This shows that the calculation model established in this paper has higher accuracy and is more consistent with the measured data, which greatly improves the accuracy of wellbore temperature and pressure.   Figure 7 shows the wellbore temperature curve obtained by using the calculation model in this paper. We can see from the drawings: (1) in the formation section, the fluid temperature in the tubing decreases slowly with the decrease in well depth; (2) in the seawater-air section, forced convection heat exchange between seawater-air and wellbore is stronger than that between wellbore and formation. Further, production fluid temperature decreases significantly faster. The annulus temperature of the formation section declines gradually with depth decreases. However, there is a sharp change near the mudline, and the annulus temperature drops rapidly. After entering the seawater section completely, the temperature of each annulus gradually decreases with the decrease in seawater depth, and the temperatures between annulus temperatures are close to each other. The maximum temperature difference is only about 7 • C. As a result, the impact of the seawater section and air section on wellbore temperature cannot be ignored for offshore wells.   Figure 8 depicts the influence of gas-liquid ratio on wellbore temperature and pressure under the conditions of gas production of 4 3 31.5 10 m d × ⋅

Gas-Liquid Ratio Production Fluid
, production time of 84 h, and the same oil-water ratio of production fluid (14:3). As can be seen from the figure, since the specific heat capacity of liquid phase is much larger than that of gas phase, the wellbore temperature becomes increasingly high with the increase in liquid oil-water con-  Figure 8 depicts the influence of gas-liquid ratio on wellbore temperature and pressure under the conditions of gas production of 31.5 × 10 4 · m 3 /d, production time of 84 h, and the same oil-water ratio of production fluid (14:3). As can be seen from the figure, since the specific heat capacity of liquid phase is much larger than that of gas phase, the wellbore temperature becomes increasingly high with the increase in liquid oil-water content in the mixed fluid. However, the increase in liquid content will increase the density of gas-liquid mixture, which makes hydrostatic pressure drop and friction pressure drop increase. As a result, the wellbore pressure drop decreases continuously, showing a trend of wellbore pressure decreasing with the decrease in the gas-liquid ratio (the increase in liquid content). It can also be seen from Equation (33) that wellbore temperature and pressure change in opposite trends. Therefore, by controlling the gas-liquid ratio reasonably, the wellbore pressure can be adjusted without changing the production process, so as to avoid the risk of increasing casing failure due to the high wellbore pressure.  Figure 8 depicts the influence of gas-liquid ratio on wellbore temperature and p sure under the conditions of gas production of 4 3 31.5 10 m d × ⋅

Gas-Liquid Ratio Production Fluid
, production time of 8 and the same oil-water ratio of production fluid (14:3). As can be seen from the fig since the specific heat capacity of liquid phase is much larger than that of gas phase, wellbore temperature becomes increasingly high with the increase in liquid oil-water c tent in the mixed fluid. However, the increase in liquid content will increase the den of gas-liquid mixture, which makes hydrostatic pressure drop and friction pressure d increase. As a result, the wellbore pressure drop decreases continuously, showing a tr of wellbore pressure decreasing with the decrease in the gas-liquid ratio (the increas liquid content). It can also be seen from Equation (33) that wellbore temperature and p sure change in opposite trends. Therefore, by controlling the gas-liquid ratio reasona the wellbore pressure can be adjusted without changing the production process, so a avoid the risk of increasing casing failure due to the high wellbore pressure.   Figure 9 shows the influence of gas production on wellbore temperature and pressure. Only gas production is changed, and other conditions are consistent with the case well. As can be seen from the figure, the gas production has a significant impact on the wellbore temperature. When the gas production rises from 10 × 10 4 · m 3 /d to 60 × 10 4 · m 3 /d, the wellhead temperature rises from 63 • C to 99 • C, which is increased by 1.6 times. However, the wellbore pressure increases first and then decreases with the increase in production because of the mutual influence of friction pressure drop and hydrostatic pressure drop. Figure 10 depicts the pressure drop in the micro-segment of the wellhead. Since the accelerated pressure drop is very small, only the friction pressure drop and hydrostatic pressure drop are shown. As can be seen from the figure, when the gas production is small, the liquid phase in the wellbore occupies a relatively high proportion, and the mixed fluid density is relatively large, resulting in a rapid hydrostatic pressure drop, and it accounts for a large proportion of the total pressure drop. However, with the increase in gas production, the fluid mixing density decreases, and the hydrostatic pressure drop slows down. At the same time, the friction at the gas-liquid interface gradually increases with the increase in the gas phase, leading to a significant decrease in the wellbore friction pressure drop, and the final wellbore pressure increases first and then decreases. tively large, resulting in a rapid hydrostatic pressure drop, and it accounts for a large p portion of the total pressure drop. However, with the increase in gas production, the f mixing density decreases, and the hydrostatic pressure drop slows down. At the same ti the friction at the gas-liquid interface gradually increases with the increase in the gas ph leading to a significant decrease in the wellbore friction pressure drop, and the final w bore pressure increases first and then decreases.   Figure 11 shows the influence of production time on wellbore temperature and pressure under the production conditions of the case well. It can be seen from the figure that the wellhead temperature increased from 91.6 o C to 96.4 o C during the period of 10 d to 100 d, while the wellhead temperature only increased by 1.8 o C during the period of 101 d to 300 d. This indicates that wellbore temperature tends to be stable with the increase in production time. The wellbore pressure does not change with the increase in production time in the 10~300 d production range.  Figure 11 shows the influence of production time on wellbore temperature and pressure under the production conditions of the case well. It can be seen from the figure that the wellhead temperature increased from 91.6 • C to 96.4 • C during the period of 10 d to 100 d, while the wellhead temperature only increased by 1.8 • C during the period of 101 d to 300 d. This indicates that wellbore temperature tends to be stable with the increase in production time. The wellbore pressure does not change with the increase in production time in the 10~300 d production range. Figure 11 shows the influence of production time on wellbore temperature and pressure under the production conditions of the case well. It can be seen from the figure that the wellhead temperature increased from 91.6 o C to 96.4 o C during the period of 10 d to 100 d, while the wellhead temperature only increased by 1.8 o C during the period of 101 d to 300 d. This indicates that wellbore temperature tends to be stable with the increase in production time. The wellbore pressure does not change with the increase in production time in the 10~300 d production range. Figure 11. The influence of production time on wellbore temperature and pressure. Figure 12 compares the variation trend of the wellhead temperature changing with the production time that was obtained by the homogeneous method model and the separated method. As can be seen from the figure, it takes 28 days for the wellhead temperature to reach 95 o C in the model based on homogeneous phase flow, which is much faster Figure 11. The influence of production time on wellbore temperature and pressure. Figure 12 compares the variation trend of the wellhead temperature changing with the production time that was obtained by the homogeneous method model and the separated method. As can be seen from the figure, it takes 28 days for the wellhead temperature to reach 95 • C in the model based on homogeneous phase flow, which is much faster than 48 days in the separated model. Therefore, the improved model can provide a more accurate estimate of the time taken to reach the rated wellhead temperature and accurate theoretical support for the rational formulation of the production plan after the well opening, so as to avoid excessive restrictions on the initial production rate.

Production Conditions
Processes 2022, 10, x FOR PEER REVIEW 18 of 20 than 48 days in the separated model. Therefore, the improved model can provide a more accurate estimate of the time taken to reach the rated wellhead temperature and accurate theoretical support for the rational formulation of the production plan after the well opening, so as to avoid excessive restrictions on the initial production rate.

Conclusions
Based on the gas-liquid two-phase separated method, an improved calculation model of wellbore temperature and pressure for gas-liquid two-phase flow in offshore gas wells is established in this paper. Further, the influencing factors of wellbore temperature and pressure are analyzed. The following conclusions are obtained: 1. Although the results of homogeneous flow and separated phase flow have the same trend, the errors of wellhead temperature and pressure obtained by homogeneous

Conclusions
Based on the gas-liquid two-phase separated method, an improved calculation model of wellbore temperature and pressure for gas-liquid two-phase flow in offshore gas wells