Adaptive analysis of multistage down-hole throttling and optimization design of process parameters in Sulige gas field

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<big>'''Adaptive Analysis of Multistage Down-hole Throttling and Optimization Design of Process Parameters in Sulige Gas Field'''</big></div>
-->

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Abstract'''</div>

''The Sulige Gas Field mainly adopts down-hole throttling production technology. For gas Wells with high bottom-hole pressure and productivity and low pressure drop in the testing process, it is difficult to reduce the wellhead pressure to the ideal value by installing a single-stage down-hole throttle, which cannot meet the needs of existing middle and low-pressure gas gathering technology in Sulige Gas field. Therefore, it is proposed to adopt down-hole multi-stage throttling technology, which not only makes full use of the residual pressure and heat of the well bore, but also ensures the safety of the surface gathering and transportation process system and the normal production of gas wells. In this paper, the optimization design method and theory for the technological parameters of the down-hole multi-stage throttling system of gas wells are established. According to the constraint of low-pressure gas collection mode in a temperature and pressure of wellhead safety to determine the number of down-hole choke and the process parameters. Discuss the applied conditions of down-hole multi-stage throttling, and provide theoretical support for the selection of the down-hole choke technology of the Sulige gas well and the safe and stable production of the gas well. It has important guiding significance for gas field quality and efficiency improvement and fine management.''

'''Key words: '''Sulige gas field; multistage throttling; optimization design.

==1. Introduction==

Sulige gas field mainly adopts the gas production technology of "Down-hole throttling, wellhead emergency cut-off, liquid metering, wellhead data remote transmission, no heating, inter-well serial connection, middle and low-pressure gas collection, normal temperature separation, and two-stage pressurization". The application of down-hole throttling can not only reduce the wellhead pressure, make it conform to the pressure level of wellhead device and surface gas gathering and transportation system, inhibit the formation of hydrates, but also greatly reduce the investment and construction cost of surface water jacket furnaces and related equipment, and save the cost of single well production and operation [1-2].

Due to the high unblocked flow rate and high well bore pressure in some new wells, if a single-stage down-hole choke is used, the residual heat and pressure of the well bore cannot be fully utilized, which will affect the productivity of the gas well and the safety of the surface medium and low-pressure gas gathering system [3-5]. However, there are few studies on the down-hole multistage throttling technology for reducing wellhead pressure in new high-pressure gas Wells. Therefore, for the newly put into production of high-pressure gas Wells, adopting down-hole multistage throttling can solve the problems of overpressure of surface gathering and transportation system and energy utilization of well bore, which has higher technical and economic benefits.

==2. Mechanism and characteristics of down-hole flow throttling==

Downhole throttling technology is to install the throttling device to a certain depth of the well bore to achieve throttling and depressurization in the well, Meanwhile, the formation heat transfer is used to heat the gas after throttling, so that the gas will not form hydrate in the well bore and gas production pipeline and block the pipeline due to the temperature reduction. Therefore, downhole throttling technology can not only reduce the pressure level of the surface natural gas gathering and transportation system, but also simplify the surface flow, save investment and realize the rapid production of gas wells. A variety of downhole throttles have been developed at home and abroad, but their throttling mechanism is the same. When natural gas passes through the throttle, there are two possible flow states,

(1) When ''p<sub>c</sub>/p<sub>1</sub>≤p<sub>2</sub>/p<sub>1</sub>≤1'', Is the non-critical state. Among them,'' p<sub>1 </sub>''is the pressure in front of the throttle,'' p<sub>2 </sub>''is the pressure at the end of the throttle, ''p<sub>c</sub>'' is the critical pressure after the throttle. Downstream of the throttle, the gas velocity decreases with the increase of ''p<sub>B</sub>/p<sub>1</sub>'', and the downstream pressure'' p<sub>2 </sub>''is equal to the back pressure.

(2) When 0''≤p<sub>B</sub>/p<sub>1</sub>≤p<sub>C</sub>/p<sub>1</sub>'', Is the critical state, Downstream of the throttle, the gas velocity reaches its maximum and the downstream pressure'' p<sub>2 </sub>''is equal to the critical pressure'' p<sub>C</sub>''. Under critical flow conditions, if there is a choke nozzle downhole, the pressure wave cannot pass through the choke nozzle and affect the upstream (well bottom) pressure when any pressure fluctuation occurs downstream, that is, between the pipe string above the choke nozzle and the gas gathering system, which can ensure the stable production of the gas well. Therefore, experts and scholars of oil and gas fields are inclined to study and design choke models for critical flow conditions. At the same time, downhole choke nozzles have auxiliary functions such as slowing downhole excitement and reducing sand production in gas wells.

==3. Dynamic analysis of well bore==

When the downhole throttle is not installed, the well bore fluid pressure and temperature are affected by the temperature of the borehole wall and show a nonlinear trend along the depth of the well, shown as in Figure 1. When there is a downhole throttle, the fluid pressure and temperature at the downstream outlet of the throttles will decrease dramatically, and then the temperature will gradually rise due to the influence of geothermal heat. The change of well bore pressure and temperature during downhole throttling is shown by the dotted line in figure 1.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Yijun_989359957-image1.png|384px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Figure 1 Schematic diagram of well bore pressure and temperature variation</div>

For downhole throttling single phase or multiphase pipe flow in gas well, the core is to explore the laws of pressure and temperature changes, combining with the hydrate formation temperature prediction of nozzle flow analysis, optimized design of throttle mouth down to the position, the throttle nozzle diameter, the multistage downhole throttling process, also include the reasonable distance between various choke, to ensure the full use of geothermal, Prevent hydrate formation, avoid well bore fluid accumulation and gas well safety production, etc. This paper refers to the pressure and temperature analysis model in reference. Using the node system analysis method, the gas reservoirs in central as the ordinary nodes, choke as function node, considering well bore heat transfer process is the steady state heat transfer, and the heat transfer process in the formation of unsteady heat transfer, the accurate mathematical solution of stable heat source structure related to time heat relational expression into the equation of conservation of energy, coupling analysis the distribution law of pressure and temperature in the entire well bore.

==4. Dynamic analysis of nozzle flow==

===4.1 Single-phase gas throttling pressure drop===

According to the downhole throttling mechanism, after the high temperature and high-pressure gas well adopts downhole throttling technology, when the natural gas flows through the downhole throttling device, the natural gas velocity will increase rapidly due to the sudden reduction of the cross-section, and part of the pressure energy in the well bore will be converted into natural gas kinetic energy. Therefore, the flow pressure at the downstream outlet of the throttle will decrease significantly. According to the principle of isentropic expansion during natural gas throttling [6], in the subcritical state, the relationship between the flow rate of natural gas and the upstream and downstream pressure ratio of throttling is,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{q}_{sc}=\frac{0.408{p}_{1}{d}^{2}}{\sqrt{{\gamma }_{g}{T}_{1}{Z}_{1}}}\sqrt{\frac{k}{k-1}}\left[ {\left( \frac{{p}_{2}}{{p}_{1}}\right) }^{\frac{2}{k}}-\right. </math><math>\left. {\left( \frac{{p}_{2}}{{p}_{1}}\right) }^{\frac{k+1}{k}}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}


Here, <math display="inline">{q}_{sc}\,</math> is natural gas production (under standard conditions), ''10<sup>4</sup>Nm<sup>3</sup>/d''. <math display="inline">{p}_{1}</math> and <math display="inline">{p}_{2}</math> are the the upstream and downstream pressures of the throttle respectively, ''MPa''.  [[Image:Draft_Yijun_989359957-image2.png|12px]] is the diameter of the orifice, mm. <math display="inline">{T}_{1}</math> is the upstream inlet temperature of the throttle, ''K''. <math display="inline">{Z}_{1}</math> is the gas compressibility coefficient at the inlet of the throttle. <math display="inline">{\gamma }_{g}</math> is the relative density of natural gas. ''k'' is the natural gas adiabatic index.

In the critical flow state, when the critical pressure ratio is calculated according to equation (1), the maximum output is

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{q}_{max}=\frac{0.408{p}_{1}{d}^{2}}{\sqrt{{\gamma }_{g}{T}_{1}{Z}_{1}}}\sqrt{\frac{k}{k-1}}\left[ {\left( \frac{2}{k+1}\right) }^{\frac{2}{k-1}}-\right. </math><math>\left. {\left( \frac{2}{k+1}\right) }^{\frac{k+1}{k-1}}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}


When the output of natural gas, the temperature and pressure at the upstream inlet of the throttle are known, substituting them into equation (1) and calculating after deformation, the pressure <math display="inline">{p}_{2}</math> after throttling can be determined, the pressure drop can be obtained as <math display="inline">{\Delta p=p}_{1}{-p}_{2}</math>.

===4.2 Gas-liquid two-phase throttle pressure drop===

Most of the gas Wells in Sulige gas field produce liquid such as condensate oil and water. In the process of exploitation, the ratio of liquid to gas will gradually increase, which will obviously affect the productivity of the gas well, and even make the gas well water flooded and stop production. Therefore, the correct prediction of gas-water two-phase throttling of gas nozzle is of great significance to the optimization design of downhole gas nozzle parameters of water-producing gas wells, gas well dynamic analysis and drainage gas production process design [7]. For gas wells with high gas-liquid ratio mist flow, the condensate oil and gas are taken as compound gas, and then the gas-liquid two-phase throttling of compound gas-water mixture is considered, and the throttling calculation model for high gas-liquid ratio gas well is obtained.

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">q=\frac{0.408{p}_{1}{d}^{2}}{\sqrt{{{F}_{w}\gamma }_{w}{T}_{1}{Z}_{1}}}\sqrt{\frac{k}{k-1}}\left[ {\left( \frac{{p}_{2}}{{p}_{1}}\right) }^{\frac{2}{k}}-\right. </math><math>\left. {\left( \frac{{p}_{2}}{{p}_{1}}\right) }^{\frac{k+1}{k}}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}


In the critical flow state, the maximum output of gas-liquid two-phase flow when the critical pressure ratio is calculated according to formula (3).

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{q}_{max}=\frac{0.408{p}_{1}{d}^{2}}{\sqrt{{{F}_{w}\gamma }_{w}{T}_{1}{Z}_{1}}}\sqrt{\frac{k}{k-1}}\left[ {\left( \frac{2}{k+1}\right) }^{\frac{2}{k-1}}-\right. </math><math>\left. {\left( \frac{2}{k+1}\right) }^{\frac{k+1}{k-1}}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}


Here, <math display="inline">{F}_{w}</math> is water cut correction coefficient, <math display="inline">{\gamma }_{w}</math> is the relative density of the composite gas and the other parameters are the same as before.

===4.3 Downhole throttling temperature drop===

The natural gas temperature at the downstream outlet of the throttling device can be obtained based on the classical thermodynamic as,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{T}_{2}={T}_{1}\frac{{Z}_{1}}{{Z}_{2}}{\left( \frac{{p}_{2}}{{p}_{1}}\right) }^{\frac{k-1}{k}}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}


Here, <math display="inline">{Z}_{1}</math> and <math display="inline">{Z}_{2}</math> are the gas compressibility coefficients at the inlet and outlet conditions of the throttle respectively. <math display="inline">{T}_{2}</math> is the downstream outlet temperature of the throttle. The phase balance equation before and after gas throttling and the throttling temperature drop model established in literature are used for calculation.

==5. Analysis of applicable conditions for multistage downhole throttling==

The application purpose of downhole multistage throttling process is basically the same as that of conventional downhole single-stage throttling process, which is to simplify surface process flow, prevent hydrate, improve liquid carrying capacity of well bore, etc., and ensure the normal production of gas well and the safety of middle and low-pressure gas gathering process system. In general, the downhole single-stage choke well selection conditions are gas wells with perfect well bore, wellhead oil pressure higher than the starting pressure of emergency cut-off valve (generally 5.0MPa), and infrequent opening and closing. The downhole multi-stage throttling process is mainly considered for gas wells with high bottom hole pressure to ensure stable production for a long time and the wellhead pressure before throttling is higher than the safe gas gathering pressure, avoiding the use of single-stage downhole throttling process to cause surface gathering system overpressure abnormal shut-in, and the well bore conditions are consistent with the single stage throttling. The conditions for the implementation of the multi-stage downhole process are obtained as follows,

[[File:Review_371620629495_8516_6.jpg]]   （6）

Here, <math display="inline">{p}_{wf}</math> is the bottom hole pressure under reasonable production allocation of gas wells, ''MPa''. <math display="inline">\, {p}_{wsmax}</math> is the highest safety gas gathering pressure at the wellhead of gas well, ''MPa. '' <math display="inline">{T}_{hp}</math>'' ''is the hydrate formation temperature corresponding to the highest safe gas gathering pressure, ''℃''. <math display="inline">{\Delta p}_{w}</math> is well bore pressure loss, ''MPa''. <math display="inline">{\Delta p}_{ci}</math> is the throttling pressure drop of the downhole throttling device, [[Image:Draft_Yijun_989359957-image3.png|36px]] . ''N'' is the number of downhole throttles. <math display="inline">{T}_{ws}</math> is the wellhead flow temperature after multistage throttling of the well bore, ''℃''.

==6. Optimization design of downhole multistage throttling parameters==

Using gas well node system analysis method and taking downhole throttle as function node, the calculation model of well bore pressure and temperature distribution is combined with the downhole throttle temperature step-down and pressure-drop model, and the established well bore nozzle flow dynamic analysis model is used to predict the distribution law of pressure and temperature in multi-stage downhole throttle process. Optimize the design of the number of throttles, throttle down depth at all levels, throttling nozzle diameter and other parameters.

===6.1 the number of downhole throttles entering===

Assuming the highest gas gathering pressure at the wellhead as the critical condition and take the reasonable production distribution of gas well and the corresponding bottom hole flow pressure into consideration, the downhole throttling series can be preliminarily obtained as,

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">N=\frac{ln{p}_{wsmax}-ln{p}_{wf}}{ln{p}_{rc}}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}


Here, <math display="inline">{p}_{rc}</math> is the critical pressure ratio of the throttle nozzle. The calculated number of downhole throttles is rounded according to the analysis results to ensure that all levels of throttles are in critical state and avoid the change of surface production parameters affecting the gas well productivity.

Assuming that the adiabatic index of Sulige gas field is 1.3, the critical pressure ratio is 0.546. According to the take-off pressure of the wellhead emergency shutoff valve as 5.0''MPa'', the one-stage throttling should be adopted when the bottom hole flow pressure is 9.2-16.78''MPa'', and when bottom hole pressure is 16.78-30.72 ''MPa'', two-stage throttling should be adopted; when bottom hole pressure is 30.72-56.26 MPa, three-stage throttling should be adopted.

===6.2 Minimum lowering depth of the first downhole throttle===

From the perspective of facilitating well bore fluid carrying, hydrate prevention and making full use of ground temperature and considering the hydrate formation temperature, the formula for calculating the minimum lowering depth of the first downhole throttle is theoretically determined as

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{L}_{min}\geq G{T}_{wf}\left[ \left( {T}_{h1}+273\right) {p}_{rc}^{-{Z}_{1}\left( k-1\right) /k}+\right. </math><math>\left. \left( {T}_{0}+273\right) \right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}


Here, <math display="inline">\, {T}_{0}</math> is the average surface temperature, ''℃''.

===6.3 Distance between multiple downhole throttles===

The distance between the throttles uses the well bore temperature field analysis theory, and the temperature after each level of the throttle is used as the starting point, and the node system analysis method is used to analyze the distance [[Image:Draft_Yijun_989359957-image4.png|18px]] when the temperature after the throttles returns to the geothermal gradient as the lowering position of the next-level throttle, namely

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{\left. \frac{\partial T\left( {Q}_{g},p,GLR,K\right) }{\partial x}\right| }_{x=\Delta x}=</math><math>G{T}_{wf}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}


Here, <math display="inline">G{T}_{wf}</math> is the average temperature gradient, ''m/℃''. <math display="inline">{Q}_{g}</math> is the reasonable distribution of production for gas wells, ''10<sup>4</sup>m<sup>3</sup>/d''. ''GLR'' is the gas and liquid ratio, ''m<sup>3</sup>/m<sup>3</sup>''. ''K'' is the total heat transfer coefficient of well bore which can be obtained according to the well structure, ''W/(m<sup>2</sup>.K)''.

The average cryogenic gradient of Sulige gas field is 2.5-2.9''℃/100m'', considering the current production parameters and well structure of the gas well, the initial calculation of the spacing between each throttle is 300-500m.

===6.4 Calculation method and precedence===

According to the above proposed theory, the specific calculation method and precedence of multi-stage downhole throttle parameter optimization design are as list follows,

(1) The production and bottom hole pressure of gas well are determined according to the reasonable production distribution method of gas well.

(2) The number of downhole throttles is determined preliminarily according to the pressure of surface gathering system.

(3) Predict the hydrate formation temperature based on the bottom hole flowing pressure, and preliminarily determine the minimum lowering depth of the downhole first-stage throttle.

(4) Calculate the temperature and pressure upstream of the first downhole throttle according to the well bore dynamic analysis method.

(5) Call the downhole throttle pressure and temperature drop model to calculate the pressure and temperature downstream of the first-stage throttle.

(6) Calculate hydrate formation temperature, ''T<sub>h1</sub>'', under downstream pressure of the throttle, if the temperature behind the nozzle (T2>Th1+3℃), reduce the lowering depth of the first throttle, otherwise increase the lowering depth of the throttle, repeat steps (4)-(6) to optimize and determine the lowering depth of the first throttle.

(7) Using the well bore fluid analysis theory to calculate the pressure and temperature at the depth under the optimization of the first-stage throttle, and calculate the diameter of the first-stage throttling nozzle.

(8) Using the well bore fluid analysis theory and node analysis method to calculate the pressure and temperature distribution behind the throttle nozzle. According to the distance between the downhole throttles, the depth of the remaining throttles and the temperature and pressure of the multi-stage throttles were determined successively, optimizing and determining the lowering depth and the diameter of the remaining throttles.

(9) If the wellhead pressure does not meet the presented safe gas gathering pressure value, readjust the number of throttles and repeat step (2)-(9) until it meets the requirements of the surface gathering and transportation system.

==7. Application==

Based on the dynamic analysis in this paper and the theory and method of parameter optimization design of downhole multi-stage throttles, carrying out the parameter optimization design of three newly put into production gas wells of ''Su X1'', ''Su X2 ''and ''Su X3''. The design results and application conditions are shown in Table 1. After the downhole first-level throttling is adopted, the pressure behind the throttling nozzle is greater than the safety gathering and transportation pressure at the wellhead. If the production is direct, there is the risk of overpressure, so it is necessary to use the gas production tree needle valve to throttle again, because the pressure before the wellhead throttling is high, and it is easy to form hydrate and block the pipeline. After the downhole two-stage throttle is adopted, the pressure behind the throttle nozzle is reduced to below the design pressure of the gas pipeline. At the same time, the temperature behind the throttle nozzle is higher than the hydrate formation temperature, which improves the safety of the wellhead and meets the needs of stable gas well production.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Table 1 Application of multi-stage downhole choke</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">No.</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Production(10<sup>4</sup>m<sup>3</sup>/d)</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Mean annual temperature (℃) </span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">choke series</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Choke depth (m)</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Pressure before choke (MPa)</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Temperature before choke (℃)</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Pressure after choke (MPa)</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Temperature after choke (℃)</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Hydrate formation temperature (℃)</span>
|-
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">SuX1</span>
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.00 </span>
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1<sup>st</sup> level</span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2501 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.88 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">89.45 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">11.40 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">42.21 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22.23 </span>
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2<sup>ed</sup> level</span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1999 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">10.93 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">74.89 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.97 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">29.55 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">18.34 </span>
|-
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">SuX2</span>
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.00 </span>
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1<sup>st</sup> level</span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2500 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.99 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">87.72 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9.82 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">40.70 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">21.34 </span>
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2<sup>ed</sup> level</span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2005 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9.43 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">73.36 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.15 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">28.22 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.45 </span>
|-
|  rowspan='2' style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">SuX3</span>
|  rowspan='2' style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">10.00 </span>
|  rowspan='2' style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1<sup>st</sup> level</span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2408 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.77 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">90.24 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">11.34 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">42.90 </span>
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">22.20 </span>
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2<sup>ed</sup> level</span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2008 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">10.97 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">78.64 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.99 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">32.81 </span>
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">18.37 </span>
|}


==8. Conclusion==

This paper optimizes the design of downhole multi-stage choke process parameters for Sulige gas wells that have been put into production, aiming at adapt to the Sulige medium and low-pressure gas gathering mode and improve the production efficiency of gas wells, from the aspects of gas well bore dynamic analysis, mouth flow dynamic analysis, applicable conditions of multi-stage throttling and optimization design of process parameters, etc., a systematic method and theory for the optimization design of process parameters of gas well downhole multi-stage throttling system have been established. Through the practical application of typical wells, it is verified that the research theory in this paper can be used to guide the selection of throttling technology and the optimization design of production parameters for new wells and high-pressure and high-yield gas wells in the Sulige Gas Field, ensuring the safety of the gathering and transportation system and the efficient and stable production of gas wells.

<span id='_GoBack'></span>

==References==

1. Duncan KW, Musaffar AK and Anas MQ (2011) Foam sticks application to increase production and encapsulated inhibitors to mitigate corrosion and scale formation. In SPE/PAPG Annual Technical Conference, Islamabad, Pakistan. Society of Petroleum Engineers, Richardson, TX, USA, paper SPE 156213.

2. Orta D, Ramanchandran S, Yang J et al. (2007) A novel foamer for deliquification of condensate-loaded wells. In SPE Rocky Mountain Oil & Gas Technology Symposium, Denver, Colorado, U.S.A. Society of Petroleum Engineers, Richardson, TX, USA, paper SPE 107980.

3. Siddiqui S and Yang J (1999) Successful application of foam for lifting liquids from low-pressure gas wells. In SPE Mid-Continent Operations Symposium, Oklahoma City, Oklahoma. Society of Petroleum Engineers, Richardson, TX, USA, paper SPE 52122.

<span style="text-align: center; font-size: 75%;">4. Wang F, Liang Z, Deng X. (2017) Automation of liquid removal from natural gas wells. Proceedings of Institution of Civil Engineers: Energy '''170(1)''':</span> <span style="text-align: center; font-size: 75%;">37-42.</span>

5. William H (2010) Gas well deliquification. In SPE Abu Dhabi International Petroleum Exhibition and Conference, Abu Dhabi, UAE. Society of Petroleum Engineers, Richardson, TX, USA, paper SPE 138672.

6. Xie C, Liu Y, Li X (2022) A novel comprehensive model for predicting production of downhole choke wells. Fuel '''313''':1-23.

7. Zheng L, Fei W, Xiong D (2014) A novel technology of combining foam injection and compression to lift liquid in water flooded gas wells. Journal of Natural Gas Science and Engineering '''19(7)''': 147–151.