#### Abstract

In order not to hinder gas production, we usually hope that the bottom hole effusion can be discharged to the surface with high-pressure natural gas. For the production data of high water content gas wells, the problems of insufficient water content and liquid-carrying capacity affecting gas well production should be considered. Based on the wellbore gas-liquid two-phase pipe flow theory and heat transfer theory, the temperature and pressure coupling prediction model of a high water-bearing gas well is established. Combined with the downhole throttling mechanism and gas-liquid two-phase homogeneous flow theory, the temperature and pressure field distribution model is established. The results show that compared with the Ramey model and Hassan and Kabir model, the temperature and pressure coupling prediction model of high water-bearing gas wells established in this study has the smallest coefficient of variation in the four groups of data tests. Based on this, the effects of different working conditions and choke diameter on downhole throttling characteristics of high water-bearing gas wells are analyzed. The findings of this study are helpful to better predict the wellbore temperature and pressure coupling of high water-bearing gas wells and provide more effective help for the smooth production of gas wells.

#### 1. Research Status of Gas-Liquid Two-Phase Well

The prediction of wellbore temperature and pressure field has always been a major problem of common concern to scientists [1] because it is related to the accurate construction of oil production technology [2, 3], the safety of pipe string [4–6], the normal use of various downhole tools [7], the reliability of production process [8], and the wax and scale prevention of oil and gas wells affected by temperature and pressure [9–12]. Nowadays, in the daily production process of gas wells, the phase change of natural gas in the wellbore is becoming more and more common [13, 14]. For example, affected by the formation and other environmental factors, formation water or other free water often appears in the wellbore of producing gas wells, and a large amount of accumulation will lead to fluid accumulation in the wellbore. The site often hopes that the water infiltrating from the formation to the bottom of the well can be carried to the ground by high-pressure natural gas fluid in time. In order to solve this engineering problem, many scholars have proposed new calculation methods to calculate and obtain the relevant data of gas well temperature, pressure field, and liquid carrying.

Abdelhafiz et al. [15] proposed a model for predicting the temperature distribution of vertical wellbore systems under circulation and shut-in conditions. The model can simulate the transient temperature disturbance of drilling fluid, drill string, casing string, cement behind the casing, and surrounding rock formation. However, the coupling effect of temperature and pressure is not considered in the model. Turner first considered the wellbore flow process and established the droplet model. Among them, the calculation formula of critical gas velocity and the droplet parameters carried by gas flow are obtained based on the droplet geometric size and force model [16]. However, the model has requirements for the gas-liquid ratio. When the gas-liquid ratio reaches a specific value (), the model can be applied to the gas well with annular fog flow. Based on Turner’s model, Gray considered the influence of temperature gradient, gas composition, fluid acceleration, and other factors during liquid carrying in gas wells and obtained the wellbore pressure drop of gas-liquid two-phase pipe flow by using an empirical model [17]. However, the Gray model has no different flow patterns, which is only applicable to stable wellbore flow and not applicable to production wells with effusion. In 2006, Ghalambor and Xu proposed a new liquid-carrying model for gas wells [18]. In view of the difficulty in calculating the pressure of multiphase flow gas wells, a simplified process is adopted. It is considered that the kinetic energy of the gas is only slightly higher than that required for liquid carrying in order to bring droplets from the bottom of the well to the surface. A method to predict critical flow is used to solve the defects of the Turner model, which cannot effectively calculate the pressure of multiphase flow gas wells, but the model only explains the droplet being brought out of the wellhead from the macro perspective, without considering the effects of droplet morphology and wellbore environment. Based on the formation process of droplets in gas wells, Li et al. proposed an ellipsoidal droplet model considering factors such as pressure and environment [19]. The calculation method of liquid-carrying critical flow is obtained. Compared with the Turner model, the calculated value is about 38% of its result, which is more accurate than the Turner model, but the model does not take into account the influence of wellbore pressure and temperature environment on liquid-carrying capacity. In 2010, Chaoyang adopted a new ellipsoid model [20], considered the influence of droplet physical model and size on stress mode, analyzed the liquid-carrying process of gas wells with liquid in the wellbore, and obtained a new ellipsoid calculation model. However, the application is limited by the field well conditions.

In 2012, Zhibin and Yingchuan considered the influence of droplet parameters, stress model, and other factors on the liquid-carrying capacity and gave a reasonable explanation for the difference in the calculation results of liquid-carrying capacity of each gas field [21], but the model ignored the change of fluid flow pattern in the wellbore. In 2012, Zhi et al. established the critical liquid-carrying model in the throttling process. It is qualitatively analyzed that downhole throttling can improve the liquid-carrying capacity of gas wells [22], but the model fails to solve the quantitative relationship between liquid-carrying flow and gas production. In 2014, Gang simplified the calculation of droplet surface tension, ignored the influence of surface tension on the critical flow of liquid carrying, and replaced it with quantitative constant. It is considered that the surface tension will change with different temperatures and pressure [23], but the quantitative relationship between surface tension and temperature and pressure is not described. In 2018, aiming at the problem of insufficient liquid-carrying capacity of low-pressure and low-yield gas wells, Yang calculated the changes of critical liquid-carrying flow at different depths [24] and obtained that high-pressure gas wells mostly occurred at the wellhead and low-pressure and low-yield gas wells mostly occurred at the bottom of the well. However, he did not establish an appropriate critical liquid-carrying flow prediction model but made a qualitative judgment. In 2020, Jinbao established the prediction model of critical liquid-carrying velocity of inclined gas wells based on the characteristics of gas-liquid two-phase flow [25] and carried out multiparameter sensitivity analysis. The prediction effect of the model based on an annular fog flow pattern is better.

Through the above data investigation, it is concluded that considering the throttling process, there is little research on the change of critical liquid carrying the flow of gas well, and the research on the change of flow pattern in wellbore on the change of critical liquid carrying the flow of gas well is rarely involved. Therefore, it is necessary to study the simulation of two-phase throttling temperature and pressure field in high water-bearing gas wells. In addition, the above methods do not mention the prediction of wellbore temperature and pressure of high water cut natural gas under high temperature and high pressure. Therefore, taking the gas wells in the Daning Jixian block as the research object, this paper not only considers the wellbore heat transfer but also considers the influence of formation temperature and wellbore pressure on the phase state of natural gas, establishing the wellbore pressure coupling model of water-bearing gas wells, and discusses the tubing diameter, gas production, and the relationship between gas water ratio and wellbore temperature field and pressure field. The technical process of this study is shown in Figure 1.

Firstly, the gas-liquid two-phase pipe flow model of a high water-bearing gas well is given. Secondly, the gas well temperature and pressure are calculated according to the model, and the relative errors are compared with the Ramey model, Hassan and Kabir model, and field measured data. Then, the downhole throttling model of high water content gas well is introduced to calculate the throttling pressure and temperature drop of high water content gas well, and an example of gas well throttling pressure and temperature drop is analyzed. Finally, the numerical simulation of downhole throttling is carried out.

#### 2. Temperature Pressure Coupling Calculation of High Water-Bearing Gas Well

Based on the wellbore heat transfer mechanism, considering the gas thermophysical parameters and Joule Thomson coefficient, a wellbore pressure-temperature coupling prediction model for high water-bearing gas wells is established. The relationship between the prediction model and Hassan and Kabir model, Ramey model, and field measured data is explored by using the coefficient of variation method, the accuracy of the model is verified, the sensitivity of wellbore temperature and pressure field is analyzed, and the relationship between tubing inner diameter, gas production, and gas water ratio and wellbore temperature and pressure field of high water cut gas wells is explored. The wellbore temperature and pressure coupling model of single-phase gas well and the temperature and pressure calculation model of high water-bearing gas well are established, respectively. The models are verified by Daji 14-1 and Daji 4-5 gas wells.

##### 2.1. Prediction Model of Wellbore Pressure Temperature Coupling Distribution in High Water-Bearing Gas Wells

Gas wells with low liquid content are generally treated as single-phase pure gas wells, but for gas wells with high water content, the pressure and temperature prediction model of single-phase gas wells cannot meet the accuracy requirements, so it is necessary to reestablish the pressure and temperature prediction model on the basis of considering water content.

In the daily production and management of natural gas wells, the gas phase plays a dominant role, and there is only a small amount of liquid. However, with the continuous extension of production time, the liquid content in the wellbore will also increase so that the flow characteristics in the whole wellbore are completely different from the flow of pure gas, reflecting the characteristics of two-phase flow. In addition, due to the continuous attenuation of temperature and pressure in the wellbore from the bottom to the ground, mass exchange between gas and liquid occurs, and the physical parameters of gas and liquid change at the same time. Therefore, it is necessary to establish a gas-liquid two-phase pipe flow model in high water-bearing gas wells to describe the dynamic model of gas and liquid wellbore flow.

The fluid flow in the wellbore is regarded as one-dimensional flow, that is, the flow parameters and physical parameters of gas and liquid phases on any section of the pipeline are uniform, which is the average value of the section. The homogeneous flow model in the two-phase flow research method is used for analysis.

###### 2.1.1. Continuity Equation

Take the one-dimensional micro element section to study; its diameter is and the pipe cross-sectional area is , as shown in Figure 2. Establish the flow model along the flow direction, and the continuity equation can be expressed as

Known

In the homogeneous flow model, the sliding velocity ratio is 1, that is, it is considered that there is no sliding between gas and liquid phases, and the volumetric void fraction is equal to the mass void fraction, so the following can be obtained:

###### 2.1.2. Momentum Equation

Similar to single-phase flow, the momentum equation of homogeneous flow can be expressed in the form of three pressure drop gradients, which can be expressed as

The pressure gradient generated by the gravity of homogeneous flow is

The friction gradient can be expressed as

Sorting:

The acceleration pressure gradient can be expressed as

Sorting:

###### 2.1.3. Energy Equation

In the homogeneous flow model, according to the principle of energy conservation, the energy conservation equation of micro element is

It can be expressed as

Substitute the basic parameters of two-phase flow:

Because there is no external work during gas-liquid two-phase flow in the wellbore, therefore

By introducing the specific enthalpy and taking the micro element section with the length of on the oil pipe, the energy conservation equation can be obtained as follows:

The specific enthalpy is a function of temperature and pressure, namely,

The heat transferred radially to the contact surface between the cement sheath and the formation is expressed as

The radial heat transfer from the cement layer to the surrounding stratum is

The heat transferred to the second contact surface is equal to the heat given to the surrounding formation by the second contact surface. The outlet temperature of each section can be obtained simultaneously:

The Joule Thomson coefficient can be expressed as [26]

where .

The wellbore pressure-temperature coupling model of the high water-bearing gas well needs to be solved by an iterative method. The wellbore is divided into sections. Assuming that the thermophysical parameters in each section are equal, the bottom hole temperature and pressure parameters are set as the boundary conditions, and the relevant physical parameters of the next section are calculated according to the temperature and pressure until the complete wellbore is calculated.

##### 2.2. Example Calculation of Gas Well Pressure and Temperature

Based on the established temperature and pressure coupling model of the high water-bearing gas well, the above model is verified in combination with the field actual wellbore data, and the relative errors between the Ramey model [27], Hassan and Kabir model, temperature-pressure coupling model of high water-bearing gas well, and field measured data are compared. The temperature and pressure distribution of wellbore temperature-pressure field with the change of tubing size, gas production, and gas water ratio is analyzed.

###### 2.2.1. Comparison of Wellbore Temperature and Pressure Field Models

Table 1 shows the wellbore structure parameters of Daji 14-1 and Daji 4-5 in a block.

For the above Daji 14-1 well, the surface casing running depth is 554.91 m, the gas reservoir casing running depth is 2471.75 m, the surface casing running depth of Daji 4-5 well is 554.64 m, and the gas reservoir casing running depth is 2486.47 m. As the overall GWR of the two wells is less than 2000, they cannot be treated as single-phase gas wells. Using the temperature pressure coupling theoretical model of high water-bearing gas wells, write the measured basic parameters of gas wells into the program (such as bottom-hole temperature, bottom hole pressure, well depth, and fluid density in the well) as the initial conditions and replace them into the temperature pressure coupling theoretical model of high water-bearing gas wells set above. The well is divided into micro element segments for an iterative solution. Because the fourth-order Runge Kutta method is more accurate to solve the differential equation, the fourth-order Runge Kutta method is used to solve the differential equation. To sum up, the calculation results through MATLAB programming are shown in Figure 3.

**(a) Comparison of wellbore temperature distribution of Daji 14-1 well**

**(b) Comparison of wellbore pressure distribution of Daji 14-1 well**

**(c) Comparison of wellbore temperature distribution of Daji 4-5 well**

**(d) Comparison of wellbore pressure distribution of Daji 4-5 well**

As shown in Figure 3, the wellbore temperature and pressure coupling prediction model of high water-bearing gas well established in this paper has higher accuracy and is closer to the field measured data than the Hassan and Kabir model and Ramey model. It can be seen from the figure that the Ramey model cannot correctly judge the wellbore temperature and pressure field of the high water-bearing gas well; the Hassan and Kabir model is more accurate than the Ramey model because it has the theoretical basis of a two-phase pipe flow model. However, when calculating the wellbore temperature and pressure field, it adopts the method of first calculating the pressure drop and independently considers the relationship between pressure change and temperature change, which has a large deviation from the field reality. For Daji 14-1 well model, the temperature prediction is more accurate and the pressure prediction error is large. There are two reasons for this phenomenon: (1) the model is caused by the discontinuity of formation pressure gradient in pressure calculation; (2) the actual measurement error of the measured pressure data will also lead to the deviation from the theoretical calculation value. For Daji 4-5 well, wellbore pressure prediction is more accurate and temperature error is relatively large. This may be due to the approximate continuous change of geothermal gradient in the calculation of formation temperature. The theoretical model can be applied to engineering practice.

###### 2.2.2. Model Error Analysis

Since the dimensions and well conditions of the four groups of data in Figure 3 are different from each other, a unified error cannot be used for comparative analysis. The coefficient of variation is introduced to describe the deviation between the models relative to the field measured data. The coefficient of variation can be expressed as

Through the coefficient of variation analysis of the four groups of data, the results are shown in Table 2.

The data discrete points in Table 2 are represented in two-dimensional coordinates, as shown in Figure 4.

As shown in Figure 4, the temperature-pressure coupling prediction model of high water-bearing gas well established in this paper has the smallest coefficient of variation in the four groups of data tests, indicating that the dispersion degree with the field measured data is the smallest and has the highest accuracy.

#### 3. Establishment of Downhole Throttling Model for High Water-Bearing Gas Wells

This paper introduces the pressure and temperature changes of natural gas fluid in the gas wellbore through the choke and puts forward the solution model of downhole throttling pressure and temperature. Mainly for single-phase gas wells and high water content gas wells in downhole throttling operation, combined with downhole throttling mechanism, the downhole throttling pressure and temperature drop model is established to predict the fluid flow parameters before and after throttling. Combined with the wellbore pressure and temperature prediction model, the whole wellbore pressure and temperature distribution under throttling conditions is obtained. Finally, the accuracy of the model is verified by comparing field data.

##### 3.1. Downhole Throttling Mechanism

For conventional gas well testing, a fixed nozzle or variable nozzle is usually used for throttling and depressurization. It is of great theoretical value and practical significance to study the variation of pressure and temperature with flow and orifice diameter before and after throttling [28].

###### 3.1.1. Thermodynamic Analysis of Downhole Throttling

The process of downhole high-pressure fluid passing through the choke belongs to isentropic expansion, accompanied by the conversion of pressure energy and kinetic energy. Natural gas does not exchange heat with the outside world during throttling, which is called adiabatic throttling [29]. Enthalpy is a function of pressure and temperature. Due to the isoenthalpy process, the pressure drop after throttling is caused by the change of temperature. Adiabatic throttling effect is usually expressed by the Joule Thomson coefficient.

The throttling coefficient is a physical parameter related to the specific state of gas, which can be calculated by looking up the diagram method and calculation formula. The average throttling coefficient can be calculated by means of the enthalpy pressure diagram and isoenthalpy calculation formula of natural gas, and the gas temperature after throttling can be obtained. If the temperature drop after throttling is too large and the gas temperature is too low, resulting in condensation or frost, it is advisable to reduce the gas adiabatic throttling coefficient, and measures can be taken to reduce the throttling pressure drop and increase the gas temperature before throttling.

###### 3.1.2. Downhole Throttling Critical Flow Conditions

When natural gas flows through the choke in the well, the flow process obeys the nozzle flow equation, that is, there are critical and subcritical areas, as shown in Figure 5.

When , the natural gas flow reaches the critical flow state, and the pressure wave generated before throttling will not be transmitted to after throttling.

When , the natural gas flow is subcritical, and the pressure wave generated before throttling will be transmitted to after throttling, resulting in wellhead pressure fluctuation [30].

Therefore, the fluid is generally required to be in a critical flow state for downhole throttling design.

##### 3.2. Calculation of Throttling Pressure and Cooling Drop in High Water-Bearing Gas Wells

###### 3.2.1. Calculation of Throttling Pressure Drop in High Water-Bearing Gas Wells

According to the theorem of momentum, there are

It is assumed that there is no mass exchange between gas and liquid phases and their respective densities do not change significantly. Continuity equations are available:

where is the change rate of temperature with pressure during natural gas throttling, if . After throttling, the temperature will be consistent with the change of pressure. is not only a function of pressure but also a function of temperature. Therefore, can only be solved approximately.

Approximate calculation

where is the average value of in the range of . The calculation method is

Substituting equations (22), (23), and (24) into equation (27) can obtain

For homogeneous flow, the section void fraction is equal to the mass void fraction, and the above formula can be simplified as

Obviously, through the above formula, on the premise of knowing the mass flow and fluid parameters, the gas-liquid two-phase throttling pressure drop can be calculated by formula (29).

###### 3.2.2. Calculation of Throttling Temperature Drop in High Water-Bearing Gas Wells

In the process of gas-liquid communication through the choke, because the specific heat capacity of the liquid is much larger than that of the gas, the temperature reduction rate of the liquid is much slower than that of the gas, so the energy transfer direction is from the liquid with higher temperature to the free gas with lower temperature. In addition, because the gas-liquid two-phase flow velocity in the throat of the choke is very fast, which is generally close to or exceeds the Mach number, the heat and mass balance of the gas-liquid two-phase cannot be established again in a very short time and can only be reestablished from the throttle outlet to the pressure recovery section. The basic equation satisfied by the heat exchange capacity is

Obviously, equation (30) expresses that the heat absorbed by the gas phase is equal to the heat released by the liquid phase when the gas and liquid are connected through the choke. Ignoring the changes of mass fraction and specific heat capacity with pressure and temperature conditions, equation (30) can be expressed as

The left end of equation (31) fully exchanges heat to obtain heat balance due to the high temperature of the liquid phase and the heat released by the liquid; the right end represents the heat absorbed by the gas due to the low temperature of the gas phase.

Generally, the enthalpy difference is used to calculate the energy change of material in the process. The enthalpy of material is the sum of ideal gas enthalpy and isothermal enthalpy difference at this temperature. The thermodynamic equation of isothermal enthalpy difference is

The equation of state of natural gas can be transformed into pressure in the form of

Substituting (32) into (31), the isothermal enthalpy difference formula can be obtained:

##### 3.3. Example Analysis of Gas Well Throttling Pressure and Temperature Drop

After obtaining the throttling pressure drop and temperature drop models of single-phase gas wells and high water content gas wells, combined with the wellbore temperature and pressure model, the wellbore temperature and pressure field under throttling condition can be analyzed. Taking Daji 4-5 gas well as an example, the wellbore temperature and pressure field under throttling conditions are analyzed and its sensitivity is also analyzed. The effects of water content, choke diameter, and gas production on the temperature-pressure field under throttling conditions are mainly considered. Combined with field data, the depth of the downhole choke is verified and analyzed.

Combined with the hydrate distribution curve of well Daji 4-5, °C and surface temperature °C can be obtained. Combined with the field data, the ground temperature gradient of well Daji 4-5 is 2.4°C/100 m, and m/°C can be obtained after conversion. Take and , and iteratively obtain at the corresponding temperature through the deviation coefficient of wellbore natural gas. Select 4 mm choke diameter and substitute formula (35) as follows:

By solving equation (37), the minimum running depth of well Daji 4-5 choke is 1800 m. The node analysis method is used to obtain the wellbore temperature and pressure distribution under throttling conditions. The steps are as follows: (1) the temperature and pressure distribution from the bottom of the well to the inlet of the choke is solved by the wellbore temperature and pressure program. (2) The temperature and pressure distribution in the throttling section is solved by using the throttling temperature and pressure drop model. (3) Taking the throttled temperature and pressure as the new initial value and substituting it into the temperature and pressure calculation program, the temperature and pressure distribution from the choke outlet to the wellhead is solved. The calculated results are shown in Figures 2 and 3.

It can be seen from Figure 6 that the wellbore temperature and pressure distribution change strongly after throttling. When the downhole choke is lowered at the position of 1800 m in the wellbore, the temperature and pressure have a sudden change at this position, and the temperature and pressure have decreased significantly. Due to the heating of the fluid in the wellbore by the ground temperature gradient, the temperature of the fluid in the wellbore will recover after throttling and then gradually decrease until the wellhead. The wellhead temperature after throttling is slightly lower than that before throttling, but the overall difference is small. After throttling, the pressure decreases greatly, making the wellhead in a low-pressure state compared with that before throttling, to achieve the purpose of preventing hydrate formation.

**(a) Temperature comparison of Daji 4-5 well before and after throttling**

**(b) Pressure comparison of Daji 4-5 well before and after throttling**

In the actual downhole throttling process, there are many factors affecting the wellbore temperature and pressure distribution after throttling, such as gas water ratio, choke diameter, and gas production. Through the field data, combined with the established throttling model, the sensitivity of the throttled wellbore is analyzed. The results are shown in Figures 7–9.

**(a) Effect of gas water ratio on throttle temperature distribution**

**(b) Effect of gas water ratio on throttle pressure distribution**

Combined with the throttling temperature and pressure distribution model of single-phase gas well and the throttling temperature and pressure distribution model of high water-bearing gas well, the throttling sensitivity of the Daji 4-5 gas well is analyzed. Figure 7 shows that (1) with the decrease of gas water ratio, the fluid water content in the wellbore increases and the throttling temperature decreases more; (2) on the premise of keeping the wellhead flow pressure unchanged, change the gas-water ratio to obtain the pressure change after throttling. With the decrease of gas-water ratio, the pressure before throttling becomes higher and higher, and the throttling pressure drop becomes larger and larger. Figure 8 shows that (1) with the increase of nozzle diameter, the throttling temperature drop becomes smaller and smaller. This is because the increase of nozzle diameter increases the amount of fluid passing through the nozzle at the unit interface, resulting in the weakening of the throttling effect, resulting in a smaller throttling temperature drop; (2) with the increase of orifice diameter, the throttling pressure drop decreases gradually. Figure 9 shows that (1) with the increase of gas production, the temperature before throttling gradually increases, and the throttling temperature drops gradually increases. The higher the production, the higher the wellhead temperature after throttling. (2) With the increase of gas production, the pressure before throttling increases gradually, and the pressure drop increases gradually.

**(a) Influence of choke nozzle diameter on throttle temperature distribution**

**(b) Influence of choke nozzle diameter on throttle pressure distribution**

**(a) Effect of gas production on throttling temperature distribution**

**(b) Effect of gas production on throttling pressure distribution**

##### 3.4. Numerical Simulation of Downhole Throttling

In order to obtain the variation of internal parameters of the throttling flow field, CFD simulation can be used. The gas-water mixing choke can be regarded as a steady and compressible flow problem of fluid through sudden contraction and sudden expansion pipeline. The purpose of the solution is to simulate the velocity field and temperature field of the downhole choke through Fluent. The working medium is gas-water two-phase fluid. Usually, the fluid flow state will change after passing through this pipeline. With reference to the actual parameters of the well example, the outer diameter is selected as *Φ*73 mm tubing movable downhole choke with an inner diameter of *Φ*62 mm, choke diameter , length , inlet section , and outlet section . The structural model is shown in Figure 10.

The structural grid in this paper is mainly generated by ICEM CFD 2020 R2, and the boundary layer at the throttle valve is encrypted. At the same time, in order to reduce the number of grids, a coarser grid is used in the part where the inlet and outlet are far away from the throttle valve. The grid diagram is shown in Figure 11. Through the grid independence analysis, the number of grids greater than 200,000 has no effect on the calculation results. Finally, about 200,000 grids are selected for calculation, as shown in Figure 12.

We use Fluent to run the model. First, open multiphase model, energy model, and viscous model, respectively. Then, in the multiphase flow model, phase 1 is defined as the main phase as methane, and phase 2 is defined as the dispersed phase as water-liquid. Throttle inlet pressure , . Throttle outlet pressure , . The flow of natural gas in the restrictor is considered as a steady state.

Using energy equation and standard model, the natural gas flow in the restrictor is regarded as a compressible fluid. The mixture homogenization model based on pressure is used to simulate the changes of pressure, velocity, and temperature in the internal flow field of the restrictor.

The general change trend of pressure is that the gas pressure before throttling does not change significantly. The cross-sectional area at the throttling inlet decreases sharply and the flow rate increases abruptly (as shown in Figure 13), resulting in the instantaneous reduction of the internal flow field from high pressure to low pressure (as shown in Figure 14). Throttling converts the pressure energy into kinetic energy. Near the throttling outlet, the flow rate of high-speed gas begins to decrease with the mixing with the surrounding gas. The pressure rises to a certain extent and then tends to be stable. It can be seen from the temperature distribution curve in Figure 15 that there is a large temperature drop process at the throttling position. The temperature of the fluid reaches the lowest value after passing through the throttling nozzle, and then, the low-temperature gas and the surrounding gas fuse and rise gradually under the action of the ground temperature gradient.

#### 4. Summary and Conclusions

Compared with the Ramey model and Hassan and Kabir model, the wellbore temperature and pressure coupling prediction model of the high water content gas well in this paper is closer to the field measured data. (1)With the decrease of gas water ratio, the fluid water content in the wellbore increases and the throttling temperature decreases more(2)On the premise of keeping the wellhead flow pressure unchanged, change the gas-water ratio to obtain the pressure change after throttling. With the decrease of gas-water ratio, the pressure before throttling becomes higher and higher, and the throttling pressure drop becomes larger and larger. With the increase of the orifice diameter, the throttling temperature drop becomes smaller and smaller. This is because the increase of the orifice diameter increases the amount of fluid passing through the orifice at the unit interface, resulting in the weakening of the throttling effect, resulting in the smaller throttling temperature drop(3)With the increase of orifice diameter, the throttling pressure drop decreases gradually. With the increase of gas production, the temperature before throttling increases gradually, and the throttling temperature drop increases gradually. The higher the production, the higher the wellhead temperature after throttling(4)With the increase of gas production, the pressure before throttling increases gradually, and the pressure drop increases gradually. There is a large temperature drop process at the throttling position. When the fluid passes through the throttling nozzle, the temperature reaches the lowest value, and then, the low-temperature gas merges with the surrounding gas and rises gradually under the action of the ground temperature gradient

In conclusion, through comparison, the model established in this paper has more accurate prediction data, but this research model is only applicable to the prediction and calculation of temperature and pressure field of high water cut two-phase flow gas wells.

#### Nomenclature

and : | Gas density and liquid density (kg/m^{3}) |

and : | Gas flow rate and liquid flow rate (m/s) |

: | Section air content |

: | Mass flow (kg/s) |

: | Mass flow rate (kg/(m^{2}·s)) |

: | Volumetric gas content |

: | Shear stress between fluid and pipe wall (N/m^{2}), |

: | The frictional resistance coefficient of two-phase flow |

: | Mass flow rate (kg/(m^{2}·s)) |

: | The specific volume of two-phase flow (kg/m^{3}) |

: | Specific enthalpy (J/kg) |

: | Flow rate of mixture (m/s) |

: | Heat of unit control body (J/m·s) |

: | Mass flow of wellbore fluid (kg/s) |

: | Specific heat capacity of fluid at constant pressure (J/(kg·K)) |

: | Joule Thomson coefficient of gas-liquid two-phase fluid |

: | Outer diameter of oil pipe (m) |

: | Total heat transfer coefficient () |

: | Wellbore fluid temperature (K) |

: | Temperature of the second contact surface (K) |

: | Formation thermal conductivity () |

: | Formation temperature at any depth (K) |

: | Dimensionless time function |

: | Fluid temperature at the outlet of each section (K) |

: | Outlet of each section (m) |

: | Entrance of each section (m) |

: | Fluid temperature at the inlet of each section (K) |

: | Formation temperature at the outlet of each section (K) |

: | Formation temperature at the inlet of each section (K) |

: | Mass flow of wellbore fluid (kg/s) |

: | Gas mass flow (kg/s) |

: | Deviation coefficient of natural gas |

: | Gas density (kg/m^{3}) |

: | Liquid density (kg/m^{3}) |

: | Coefficient of variation |

: | Standard deviation |

: | Average value |

: | Mass flow (kg/s) |

: | Flow rate (m/s) |

: | Throat sectional area (m^{2}); and , respectively, represent gas phase and liquid phase, and 1 and 2 represent throttle inlet and throttle throat, respectively |

: | Mass air content |

: | Section void content |

: | Mole fraction of component in the gas phase |

: | Temperature change before and after throttling (K) |

: | Pressure change before and after throttling (MPa) |

and : | Pressure value before and after throttling (MPa) |

: | Gas phase mass fraction |

: | Specific heat capacity of liquid phase at constant volume (kJ/(kg·K)) |

: | Specific heat capacity of gas phase at constant volume (kJ/(kg·K)) |

: | Temperature before throttling (K) |

: | Gas-liquid two-phase heat balance temperature after throttling (K) |

: | Temperature under downstream pressure of isoenthalpy process (K) |

: | Adiabatic index of natural gas with different components |

: | Comparison temperature of component (K) |

: | Gas constant (8314 kJ/(kmol·K)). |

#### Data Availability

The manuscript is a data self-contained article, whose results were obtained from the laboratory analysis, and the entire data are presented within the article.

#### Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

#### Acknowledgments

This project was supported by the National Natural Science Foundation of China (52004215, 12101482, and 51674199) and Graduate Innovation and Practical Ability Training Program of Xi’an Shiyou University (YCS211053).