Study on Calculation Method for Wellbore Pressure in Gas Wells with Large Liquid Production

: In order to solve the inaccuracy in the calculation of the wellbore pressure distribution caused by large liquid production in the Donghai gas ﬁeld, the gas–liquid production conditions of the Donghai gas ﬁeld were simulated with indoor experiments, and the ﬂow patterns for different pipe diameters, different inclinations, and different ﬂow patterns were systematically analyzed using a ﬂow pattern discrimination method, liquid holdup calculation method, and pressure drop calculation method. Using the experimental data, the division methods for different ﬂow patterns were screened. Finally, based on the fact that the change trend for the ﬂow patterns was consistent, the Kaya–Sarica–Brill method was selected to establish the ﬂow pattern discrimination formula. According to the calculation method for the Mukherjee and Brill (M–B) liquid holdup, the M–B model was re-established according to a 75 mm pipe diameter and 60 mm pipe diameter using the instantaneous liquid holdup measured in the laboratory. Through the comparison and analysis of the measured data and the calculated data for the Beggs and Brill (B–B_ pressure drop model under the same working conditions, it was found that when the B–B model was applied to different angles and different gas–liquid ratios, the error decreased with an increase in the angle and increased with a decrease in the gas–liquid ratio. After verifying the correlation of different dimensionless numbers that can characterize the gas–liquid ratio, it was considered that the introduction of the Reynolds number into the original model could greatly improve the accuracy of the calculation, so a new pressure drop calculation model was established. The new pressure drop calculation model takes into account the two parameters of the well deviation angle and gas–liquid ratio. The accuracy was greatly improved, as veriﬁed by ﬁeld measurements in four wells.


Introduction
With the substantial increase in the proportion of natural gas exploitation, the monitoring of gas wells is becoming ever more important. The correct calculation of the gas wellbore pressure distribution is related to the optimization of the completion string, the optimization of gas-well production parameters, the judgment of the amount of fluid accumulation in the gas well, and the prediction of the height of the fluid accumulation. After much research in China and other countries, a variety of gas-liquid two-phase flow and pressure drop calculation models have been established for the calculation of the wellbore pressure distribution.
There are six pressure drop calculation models commonly used for two-phase pipe flows, namely, the Aziz, Govier, and Fogarasi model [1], Beggs and Brill model [2,3], Hasan and correct the appropriate model on this basis. Chen Dechun evaluated eight common models by sorting pressure measurement data in China and other countries in 2017 [18]. He found that the same model error fluctuated under different gas-liquid ratios, which also confirmed that the applicable conditions for different models have certain limitations. According to the pressure drop optimization idea of Liao's team and Mukherjee and Brill, a matched indoor experiment was established, the flow pattern and pressure drop calculation methods of predecessors were selected and optimized, and a new high-precision pressure drop calculation method was established according to the production conditions of the Donghai gas field.

Experimental Phenomena and Flow Pattern Discrimination
The experiment was carried out in the Key Laboratory of Multiphase Flow of PetroChina at Yangtze University. The experimental devices had 75 mm and 60 mm pipe diameters and a total length of 11.5 m. A pressure sensor, a temperature sensor, a differential pressure sensor, a quick closing valve, and other devices were installed on the pipe section. The distance between the two quick-closing valves was 9.5 m, including a 7 m plexiglass pipe and 2.4 m stainless-steel pipe [19]. The drain valve was used to place the liquid remaining in the wellbore after the quick-closing valve was closed into the measuring cylinder for a liquid-phase volume measurement to measure the liquid holdup, as shown in Figure 1.
Processes 2022, 10, x FOR PEER REVIEW 3 of 21 mathematical model is not rigorous. It is necessary to verify whether the views of Zhang Tingting and others are correct according to the physical model experimental data and correct the appropriate model on this basis. Chen Dechun evaluated eight common models by sorting pressure measurement data in China and other countries in 2017 [18]. He found that the same model error fluctuated under different gas-liquid ratios, which also confirmed that the applicable conditions for different models have certain limitations. According to the pressure drop optimization idea of Liao's team and Mukherjee and Brill, a matched indoor experiment was established, the flow pattern and pressure drop calculation methods of predecessors were selected and optimized, and a new high-precision pressure drop calculation method was established according to the production conditions of the Donghai gas field.

Experimental Phenomena and Flow Pattern Discrimination
The experiment was carried out in the Key Laboratory of Multiphase Flow of Petro-China at Yangtze University. The experimental devices had 75 mm and 60 mm pipe diameters and a total length of 11.5 m. A pressure sensor, a temperature sensor, a differential pressure sensor, a quick closing valve, and other devices were installed on the pipe section. The distance between the two quick-closing valves was 9.5 m, including a 7 m plexiglass pipe and 2.4 m stainless-steel pipe [19]. The drain valve was used to place the liquid remaining in the wellbore after the quick-closing valve was closed into the measuring cylinder for a liquid-phase volume measurement to measure the liquid holdup, as shown in Figure 1. The gas used in the experiment was air, and the liquid was tap water. The surface tension of water is about 72.8 mN/m at 20 °C, the densities of water and air are 1000 and 1.05 kg/m 3 , respectively, and the viscosities of water and air are 1 and 0.0181 mPa·s, respectively.
For different gas-liquid volumes and gas-liquid ratios, experiments on the variation law for the pressure drop under different conditions were carried out.
A total of 760 groups of gas-liquid flow pattern experiments with pipe section angles of 10°, 30°, 40°, 45°, 50°, 55°, 60°, 70°, and 90° were completed. By adjusting the gas-liquid ratio with a fixed amount of liquid, the physical experimental simulation of the flow pattern at different angles was carried out, and the bubble flow, slug flow, stirred flow, and annular flow were observed, respectively. The experimental gas-liquid range is shown in Table 1. The gas used in the experiment was air, and the liquid was tap water. The surface tension of water is about 72.8 mN/m at 20 • C, the densities of water and air are 1000 and 1.05 kg/m 3 , respectively, and the viscosities of water and air are 1 and 0.0181 mPa·s, respectively.
For different gas-liquid volumes and gas-liquid ratios, experiments on the variation law for the pressure drop under different conditions were carried out.
A total of 760 groups of gas-liquid flow pattern experiments with pipe section angles of 10 • , 30 • , 40 • , 45 • , 50 • , 55 • , 60 • , 70 • , and 90 • were completed. By adjusting the gas-liquid ratio with a fixed amount of liquid, the physical experimental simulation of the flow pattern at different angles was carried out, and the bubble flow, slug flow, stirred flow, and annular flow were observed, respectively. The experimental gas-liquid range is shown in Table 1.
The experimental picture is shown in Figure 2 below. Through the observation of the experimental phenomenon of a large gas-liquid twophase pipe flow, it was found that the change in flow pattern was slightly different from that for small flows. According to the experimental data for a large water volume and large air volume under experimental conditions, various commonly used flow pattern identification methods were screened. The preliminary screening result was that the Kaya-Sarica-Brill flow pattern identification method was more accurate in wellbore flow pattern identification for the Donghai gas field, as shown in Figure 3. The reason is that the change trend for the Kaya-Sarica-Brill model is basically consistent with the flow pattern boundary trend observed in the experiment. The experimental picture is shown in Figure 2 below.  Through the observation of the experimental phenomenon of a large gas-liquid twophase pipe flow, it was found that the change in flow pattern was slightly different from that for small flows. According to the experimental data for a large water volume and large air volume under experimental conditions, various commonly used flow pattern identification methods were screened. The preliminary screening result was that the Kaya-Sarica-Brill flow pattern identification method was more accurate in wellbore flow pattern identification for the Donghai gas field, as shown in Figure 3. The reason is that the change trend for the Kaya-Sarica-Brill model is basically consistent with the flow pattern boundary trend observed in the experiment.  Considering that the basic data for using superficial gas velocity (vsg) and superficial liquid velocity (vsl) for flow pattern judgment are not comprehensive, the Aziz-Govier-Fogarasi model and Duns-Ros method [10] based on the speed criterion were also considered, as shown in Figures 4 and 5. Since most of the applications of the Orkiszewski  Considering that the basic data for using superficial gas velocity (vsg) and superficial liquid velocity (vsl) for flow pattern judgment are not comprehensive, the Aziz-Govier-Fogarasi model and Duns-Ros method [10] based on the speed criterion were also considered, as shown in Figures 4 and 5. Since most of the applications of the Orkiszewski method use the Duns-Ros method for flow pattern discrimination, the Duns-Ros flow pattern discrimination method was verified using data. Considering that the basic data for using superficial gas velocity (vsg) and superficial liquid velocity (vsl) for flow pattern judgment are not comprehensive, the Aziz-Govier-Fogarasi model and Duns-Ros method [10] based on the speed criterion were also considered, as shown in Figures 4 and 5. Since most of the applications of the Orkiszewski method use the Duns-Ros method for flow pattern discrimination, the Duns-Ros flow pattern discrimination method was verified using data.
Processes 2022, 10, 685 6 of 20 where v sl is the liquid superficial velocity, m/s; v sg is the gas superficial velocity, m/s; ρ l is the liquid density, kg/m 3 ; ρ g is the gas density, kg/m 3 ; σ is the liquid-phase surface tension, mN/m; ρ a is the air density under standard conditions, kg/m 3 ; ρ w is the water density under standard conditions, kg/m 3 ; σ w is the water-phase surface tension, mN/m. From Figures 3-5, it can be observed that the change trend for the Kaya-Sarica-Brill model is roughly the same as that for the experimental data. Finally, the Kaya-Sarica-Brill model was still selected as the basis for the flow pattern discrimination formula for correction.
The experimental data were introduced to the Kaya-Sarica-Brill model and shown below. The flow pattern was small under the condition of large flow, and the prediction was often inaccurate.
After fitting, the new boundaries of each flow-pattern discrimination were as shown in   The new flow pattern discrimination formula is as follows: The boundary of the bubble slug flow pattern transition: The boundary of the flow pattern transition between slug flow and stirred flow:    The new flow pattern discrimination formula is as follows: The boundary of the bubble slug flow pattern transition: The boundary of the flow pattern transition between slug flow and stirred flow: The boundary of the flow pattern transition between stirred flow and annular flow: where is the liquid superficial velocity, m/s; The boundary of the flow pattern transition between stirred flow and annular flow: where v sl is the liquid superficial velocity, m/s; v sg is the gas superficial velocity, m/s; σ is the liquid-phase surface tension, mN/m; θ is the pipe and horizontal inclination angle, • .
The new flow-pattern-discrimination formula is as follows:

Calculation Formula for Liquid Holdup
Through the analysis of 760 groups of gas-liquid flow pattern experiments under different angles, the accuracy of the liquid holdup calculation formula under different flow patterns and flow patterns was verified, respectively.
By sorting and analyzing the experimental data, the law of change for the liquid holdup under different conditions was drawn, shown in Figures 10 and 11. It can be seen from the figure that, under the same liquid flow condition, the liquid holdup rate decreases with an increase in the gas flow rate]. Under the same gas flow rate condition, the liquid holdup rate gradually increases with an increase in the liquid flow rate, but the change is not significant.

Calculation Formula for Liquid Holdup
Through the analysis of 760 groups of gas-liquid flow pattern experiments under different angles, the accuracy of the liquid holdup calculation formula under different flow patterns and flow patterns was verified, respectively.
By sorting and analyzing the experimental data, the law of change for the liquid holdup under different conditions was drawn, shown in Figures 10 and 11. It can be seen from the figure that, under the same liquid flow condition, the liquid holdup rate decreases with an increase in the gas flow rate. Under the same gas flow rate condition, the liquid holdup rate gradually increases with an increase in the liquid flow rate, but the change is not significant.

Calculation Formula for Liquid Holdup
Through the analysis of 760 groups of gas-liquid flow pattern experiments under different angles, the accuracy of the liquid holdup calculation formula under different flow patterns and flow patterns was verified, respectively.
By sorting and analyzing the experimental data, the law of change for the liquid holdup under different conditions was drawn, shown in Figures 10 and 11. It can be seen from the figure that, under the same liquid flow condition, the liquid holdup rate decreases with an increase in the gas flow rate]. Under the same gas flow rate condition, the liquid holdup rate gradually increases with an increase in the liquid flow rate, but the change is not significant.  There are many methods for calculating the liquid holdup. The Mukherjee and Brill method was established to calculate the liquid holdup, and the Beggs and Brill model also puts forward a calculation method for the liquid holdup. However, when the B-B method is used to calculate the liquid holdup of different flow patterns, there is discontinuity, so it is difficult to ensure the accuracy of the liquid holdup calculation [19][20][21]. The empirical formula provided by the Mukherjee and Brill method has good continuity, considers the influencing factors comprehensively, and is suitable for all flow patterns at all angles.
where is the apparent velocity of the liquid phase, m/s; is the apparent velocity of the gas phase, m/s; is the surface tension of the liquid phase, N/m; is the viscosity of the liquid phase, Pa·s.
Based on the good continuity of the Mukherjee correlation equation for the liquid holdup and comprehensive consideration of factors, the empirical constant in the relationship was chosen to be revised. According to the liquid holdup formula for gas-liquid twophase flows provided by Mukherjee and Brill, the relationship between the experimental liquid holdup Hl and the correlation criteria Nvl, Nvg, and Nl provided by Mukherjee and Brill was fitted. A comparison between the corrected calculated values and the actual values is shown in Figures 12 and 13. There are many methods for calculating the liquid holdup. The Mukherjee and Brill method was established to calculate the liquid holdup, and the Beggs and Brill model also puts forward a calculation method for the liquid holdup. However, when the B-B method is used to calculate the liquid holdup of different flow patterns, there is discontinuity, so it is difficult to ensure the accuracy of the liquid holdup calculation [19][20][21]. The empirical formula provided by the Mukherjee and Brill method has good continuity, considers the influencing factors comprehensively, and is suitable for all flow patterns at all angles.
where v sl is the apparent velocity of the liquid phase, m/s; v sg is the apparent velocity of the gas phase, m/s; σ is the surface tension of the liquid phase, N/m; µ l is the viscosity of the liquid phase, Pa·s.

Pressure Calculation Formula
Based on the accurate establishment of the liquid holdup calculation model and flow pattern discrimination method, the pressure drop calculation formula was optimized and established.
In the Beggs and Brill pressure drop calculation model, the calculation of the friction coefficient is related to the flow pattern [3]. Therefore, this coefficient is used as the main correction parameter. During this test, the flow ranges of the gas and liquid phases were relatively large, and the experimental phenomena covered all the flow pattern areas; The blue line represents the fitted value, and the red line represents the experimental value. The x-axis represents the number of iterations, and the y-axis represents the value of the result.

Pressure Calculation Formula
Based on the accurate establishment of the liquid holdup calculation model and flow pattern discrimination method, the pressure drop calculation formula was optimized and established.
In the Beggs and Brill pressure drop calculation model, the calculation of the friction coefficient is related to the flow pattern [3]. Therefore, this coefficient is used as the main correction parameter. During this test, the flow ranges of the gas and liquid phases were relatively large, and the experimental phenomena covered all the flow pattern areas; The blue line represents the fitted value, and the red line represents the experimental value. The x-axis represents the number of iterations, and the y-axis represents the value of the result.
The new empirical constant obtained by fitting and the corrected liquid holdup error are shown in Tables 2 and 3, respectively.

Pressure Calculation Formula
Based on the accurate establishment of the liquid holdup calculation model and flow pattern discrimination method, the pressure drop calculation formula was optimized and established.
In the Beggs and Brill pressure drop calculation model, the calculation of the friction coefficient is related to the flow pattern [3]. Therefore, this coefficient is used as the main correction parameter. During this test, the flow ranges of the gas and liquid phases were relatively large, and the experimental phenomena covered all the flow pattern areas; corresponding calculation methods for drag coefficients along the way were provided for different flow patterns.
The drag coefficient of the gas-liquid two-phase flow was obtained according to the Beggs and Brill model.
In the formula λ is the "no slip" drag coefficient along the way, dimensionless; s is the index; Re is the Reynolds number without slip; µ l and µ g are the viscosities of the liquid phase and gas phase, Pa·s.
The index s in the formula can be calculated according to the following formula: After verifying the B-B model with the data under the condition of the gas-liquid ratio in the laboratory, it was found that the error was large, and the average absolute deviation is shown in Table 4. As shown in Figures 14-17, it can be observed from the error analysis diagram that, although the distribution of the data points is irregular, the overall trend is that the error of the B-B model decreases with an increase in the gas-liquid ratio. This phenomenon is also mentioned in Chen Dechun's research [18]. The error of the B-B model is better when the gas-liquid ratio is low. However, he mainly used the field data to evaluate the existing model and did not propose a solution. Moreover, the analysis of the existing model used the field-measured data, there was a lack of wellbore flow data for large water volume gas wells, and the evaluation focused on the data with a gas-liquid ratio of more than 1750. With an increase in the gas-liquid ratio, the error of the B-B model will decrease significantly. Since the error decreases with an increase in the gas-liquid ratio, the dimensionless number that can characterize the gas-liquid flow is considered to re-establish the formula.     The Froude number (FR) is a dimensionless scalar in fluid mechanics; it is the ratio of the fluid kinetic energy to potential energy. When fr > 1, the water depth is small, and the velocity is turbulent. Similarly, when fr < 1, the velocity is slow, and the water depth is large. Similarly, the Reynolds number (Re) is also a classical dimensionless number used to characterize fluid flows in fluid mechanics. The gas-liquid velocity ratio is a more intuitive reaction gas-liquid ratio. In addition to the gas-liquid velocity ratio, the Froude number and Reynolds number can also characterize the effect of fluid flows from the side [22][23][24][25][26]. In order to improve the calculation accuracy for the pressure drop calculation formula, the three parameters of the Froude number, Reynolds number, and gas-liquid velocity ratio were brought into the model for correlation verification, and the most appropriate parameters representing the gas-liquid ratio were selected for the establishment of a new model.  The Froude number (FR) is a dimensionless scalar in fluid mechanics; it is the ratio of the fluid kinetic energy to potential energy. When fr > 1, the water depth is small, and the velocity is turbulent. Similarly, when fr < 1, the velocity is slow, and the water depth is large. Similarly, the Reynolds number (Re) is also a classical dimensionless number used  The Froude number (FR) is a dimensionless scalar in fluid mechanics; it is the ratio of the fluid kinetic energy to potential energy. When fr > 1, the water depth is small, and the velocity is turbulent. Similarly, when fr < 1, the velocity is slow, and the water depth is large. Similarly, the Reynolds number (Re) is also a classical dimensionless number used According to the pressure drop data for gas-liquid two-phase pipe flows with large amounts of water and gas obtained under the experimental conditions, the gas-liquid ratio was introduced, combined with the pressure drop calculation formula of the B-B model, and the correlation was verified according to the experimental data. The correlation coefficient obtained using the new formula was r = 0.81. The correlation between the calculated value and the actual value for the new formula is shown in Figure 18.
tuitive reaction gas-liquid ratio. In addition to the gas-liquid velocity ratio, the Froude number and Reynolds number can also characterize the effect of fluid flows from the side [22][23][24][25][26]. In order to improve the calculation accuracy for the pressure drop calculation formula, the three parameters of the Froude number, Reynolds number, and gas-liquid velocity ratio were brought into the model for correlation verification, and the most appropriate parameters representing the gas-liquid ratio were selected for the establishment of a new model.
According to the pressure drop data for gas-liquid two-phase pipe flows with large amounts of water and gas obtained under the experimental conditions, the gas-liquid ratio was introduced, combined with the pressure drop calculation formula of the B-B model, and the correlation was verified according to the experimental data. The correlation coefficient obtained using the new formula was r = 0.81. The correlation between the calculated value and the actual value for the new formula is shown in Figure 18. Figure 18. When the gas-liquid ratio is substituted into the calculation formula, the correlation between the calculated value and the experimental value is as shown (xyz axis represents , gasliquid ratio and pressure drop (kPa), respectively).
By introducing the Froude number and combining it with the pressure drop-calculation formula of the B-B model, the correlation between the Froude number and the calculation results for the B-B model was verified according to the experimental data. The correlation coefficient of the Froude number was r = 0.83, and the correlation between the calculated value and the actual value for the new formula is shown in Figure 19. Figure 18. When the gas-liquid ratio is substituted into the calculation formula, the correlation between the calculated value and the experimental value is as shown (xyz axis represents v sg , gasliquid ratio and pressure drop (kPa), respectively).
By introducing the Froude number and combining it with the pressure drop-calculation formula of the B-B model, the correlation between the Froude number and the calculation results for the B-B model was verified according to the experimental data. The correlation coefficient of the Froude number was r = 0.83, and the correlation between the calculated value and the actual value for the new formula is shown in Figure 19.
By correlating the Reynolds number with experimental data and the B-B model, it was found that the Reynolds number had a better correlation after fitting than the gasliquid ratio and Froude number, as shown in Figure 20. Therefore, the Reynolds number was selected as the dimensionless number for characterizing the gas-liquid ratio for the optimization of the new model: The above figure shows the correlation verification based on the optimal relationship automatically fitted by the TableCurve software. The correlation verification can be used for the three parameters of the Froude number, gas-liquid ratio, and Reynolds number as the parameters for characterizing the gas-liquid ratio, which can be applied to the adaptability screening in the optimization of the B-B model (selecting the optimal parameters through the order of correlation R under the recommended optimal relationship). Finally, the Reynolds number was selected and brought into the B-B model, and it was refitted according to the relationship obtained using TableCurve. For different angles and different flow patterns, the relationship is as follows: Figure 19. When the Froude number is substituted into the calculation formula, the correlation between the calculated value and the experimental value is as shown (xyz axis represents , Froude number and pressure drop (kPa), respectively).
By correlating the Reynolds number with experimental data and the B-B model, it was found that the Reynolds number had a better correlation after fitting than the gasliquid ratio and Froude number, as shown in Figure 20. Therefore, the Reynolds number was selected as the dimensionless number for characterizing the gas-liquid ratio for the optimization of the new model: Figure 20. When the Reynolds number is substituted into the calculation formula, the correlation between the calculated value and the experimental value is as shown(XYZ axis represents , Reynolds number and pressure drop (kPa), respectively).  By correlating the Reynolds number with experimental data and the B-B model, it was found that the Reynolds number had a better correlation after fitting than the gasliquid ratio and Froude number, as shown in Figure 20. Therefore, the Reynolds number was selected as the dimensionless number for characterizing the gas-liquid ratio for the optimization of the new model:  When the Reynolds number is substituted into the calculation formula, the correlation between the calculated value and the experimental value is as shown (XYZ axis represents v sg , Reynolds number and pressure drop (kPa), respectively).
The above formula was fitted with experimental data, and the pressure drop formula under each flow pattern was fitted with the Levenberg Marquardt method + general global optimization method. The new parameters are shown in Table 5.

Field Data Verification
In the evaluation and optimization of the pressure drop calculation model for gasliquid two-phase pipe flows in gas wells, Chen Dechun used the pressure measurement data for 102 wells in China and other countries. He conclude that "the gray model is applied to domestic gas wells and the calculated value is too small, while the gray model is applied to foreign gas wells, the calculation result is more accurate". Since Chen Dechun has discussed the Gray model, the Gray model will not be discussed in this article. Comparing the gas well data selected by Chen Dechun, there is no significant change in gas production, and the gas-liquid ratio ranges from 444.4~75,000 m 3 /m 3 to 2119.5~200,796.0 m 3 /m 3 . This shows that some pressure drop calculation models will not be applicable at high liquid volumes.
The results of the comparison between the calculated and measured values of the wellbore pressure distributions of six models are shown in The reason the error of these models is too large is that the water production of the Donghai gas field is huge, and the daily liquid production of a single well can reach 300 square meters. As the influence of the gas-liquid ratio on the pressure drop was not considered, the calculation of the above pressure drop model will be inaccurate.       The errors between different models are shown in Table 6. The errors between different models are shown in Table 6.

Conclusions
Indoor experiments were used to simulate the production status of gas wells under the conditions of large water content and large gas production. It was verified from the aspects of the flow pattern, liquid holdup, and pressure drop that the current model has the problem of inaccurate prediction when applied to large-volume production gas wells. The flow pattern, flow pattern discrimination method, liquid holdup calculation method, and pressure drop calculation method for a large flow range were optimized using indoor experimental data. After analyzing the calculation error of the B-B pressure calculation model when applied to atmospheric liquid flow, it was found that the error of the B-B model increased with an increase in the gas-liquid ratio and increased with a decrease in the angle. In terms of innovation, through many experimental data, the distribution law for the convection type and the calculation method for liquid holdup were optimized, and it was found that the B-B model is not applicable under the production condition of a large water well. Based on the experimental data, the distribution law for the pressure drop calculation error was analyzed. The calculated value, actual value, and possible influencing factors were analyzed using the TableCurve software. It is concluded that the introduction of the Reynolds number can greatly improve the calculation accuracy. On the basis of considering the influence of the angle and gas-liquid ratio on the calculation accuracy for the pressure drop, a new pressure drop-calculation formula was established. Therefore, the Reynolds number and angle-correction coefficients were introduced to re-establish the formula, and the on-site measured data for the East China Sea gas field were used to verify that the newly established pressure drop calculation model had a higher calculation accuracy in the high flow range. The field-measured data for the Donghai gas field verify that the newly established pressure drop calculation model has a high calculation accuracy in a high flow range. The research results lay a solid foundation for improving the accuracy of the prediction of the wellbore fluid pressure drop distributions in water-bearing gas wells. Additionally, a new method for calculating the pressure distribution along the gas well in the East China Sea gas field is provided.