(6.3 Distance between multiple downhole throttles)
 
(25 intermediate revisions by the same user not shown)
Line 50: Line 50:
 
{| style="text-align: center; margin:auto;"  
 
{| style="text-align: center; margin:auto;"  
 
|-
 
|-
| <math display="inline">{q}_{sc}=\displaystyle \frac{0.408{p}_{1}{d}^{2}}{\sqrt{{\gamma }_{g}{T}_{1}{Z}_{1}}}\sqrt{\displaystyle \frac{k}{k-1}}\left[ {\left( \displaystyle \frac{{p}_{2}}{{p}_{1}}\right) }^{\displaystyle \frac{2}{k}}-\right. </math><math>\left. {\left( \frac{{p}_{2}}{{p}_{1}}\right) }^{\frac{k+1}{k}}\right]</math>
+
| <math>{q}_{sc}=\displaystyle \frac{0.408{p}_{1}{d}^{2}}{\sqrt{{\gamma }_{g}{T}_{1}{Z}_{1}}}\sqrt{\displaystyle \frac{k}{k-1}}\left[ {\left( \displaystyle \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)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
 
|}
 
|}
  
 +
where, <math display="inline">{q}_{sc}\,</math> is natural gas production (under standard conditions), <math>10^4</math>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;  <math display="inline"> d </math> is the diameter of the orifice, mm; <math display="inline">{T}_{1}</math> is the upstream inlet temperature of the throttle, <math display="inline"> K </math>; <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; <math display="inline">k</math> is the natural gas adiabatic index.
  
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 Eq.(1), the maximum output is
 
+
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;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 65: Line 64:
 
{| style="text-align: center; margin:auto;"  
 
{| 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>
+
| <math>{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}}- {\left( \frac{2}{k+1}\right) }^{\frac{k+1}{k-1}}\right]</math>
 
|}
 
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
Line 71: Line 70:
  
  
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>.
+
When the output of natural gas, the temperature and pressure at the upstream inlet of the throttle are known, substituting them into Eq.(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===
 
===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.
+
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;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 82: Line 81:
 
{| style="text-align: center; margin:auto;"  
 
{| 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>
+
| <math>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)
 
| 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 Eq.(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;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 95: Line 93:
 
{| style="text-align: center; margin:auto;"  
 
{| 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>
+
| <math>{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}}- {\left( \frac{2}{k+1}\right) }^{\frac{k+1}{k-1}}\right]</math>
 
|}
 
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
 
|}
 
|}
  
 
+
where, <math display="inline">{F}_{w}</math> is water cut correction coefficient, and <math display="inline">{\gamma }_{w}</math> is the relative density of the composite gas and the other parameters are the same as before.
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===
 
===4.3 Downhole throttling temperature drop===
Line 112: Line 109:
 
{| style="text-align: center; margin:auto;"  
 
{| 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>
+
| <math>{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)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
 
|}
 
|}
  
 
+
where, <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, and <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.
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==
 
==5. Analysis of applicable conditions for multistage downhole throttling==
Line 124: Line 120:
 
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,
 
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)
+
{| class="formulaSCP" style="width: 100%; text-align: left;"
 +
|-
 +
|
 +
{| style="text-align: center; margin:auto;width: 100%;"
 +
|-
 +
| style="text-align: center;" | <math>\left\{\begin{matrix}\displaystyle p_{wf} -\sum_{i=1}^N \Delta p_{ci} - \Delta p_{w} \le \Delta p_{wsmax}\\\displaystyle T_{ws} \le T_{hp} + 3^\circ\mbox{C} \end{matrix}    \right.</math>
 +
|}
 +
| style="width: 5px;text-align: right;white-space: nowrap;" |(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, ''℃''.
+
where, <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">^\circ\mbox{C}</math>; <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, MPa; <math display="inline">N</math> is the number of downhole throttles. <math display="inline">{T}_{ws}</math> is the wellhead flow temperature after multistage throttling of the well bore, <math display="inline">^\circ\mbox{C}</math>.
  
 
==6. Optimization design of downhole multistage throttling parameters==
 
==6. Optimization design of downhole multistage throttling parameters==
Line 141: Line 145:
 
{| style="text-align: center; margin:auto;"  
 
{| style="text-align: center; margin:auto;"  
 
|-
 
|-
| <math display="inline">N=\frac{ln{p}_{wsmax}-ln{p}_{wf}}{ln{p}_{rc}}</math>
+
| <math>N=\frac{ln{p}_{wsmax}-ln{p}_{wf}}{ln{p}_{rc}}</math>
 
|}
 
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
 
|}
 
|}
  
 +
where <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.
  
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.
 
+
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===
 
===6.2 Minimum lowering depth of the first downhole throttle===
Line 160: Line 163:
 
{| style="text-align: center; margin:auto;"  
 
{| 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>
+
| <math>{L}_{min}\geq G{T}_{wf}\left[ \left( {T}_{h1}+273\right) {p}_{rc}^{-{Z}_{1}\left( k-1\right) /k}+ \left( {T}_{0}+273\right) \right]</math>
 
|}
 
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
 
|}
 
|}
  
 
+
where <math display="inline">{T}_{0}</math> is the average surface temperature, <math display="inline">^\circ\mbox{C}</math>.
Here, <math display="inline">\, {T}_{0}</math> is the average surface temperature, ''℃''.
+
  
 
===6.3 Distance between multiple downhole throttles===
 
===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
+
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 <math display="inline"> \Delta x </math> 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;"  
 
{| class="formulaSCP" style="width: 100%; text-align: center;"  
Line 177: Line 179:
 
{| style="text-align: center; margin:auto;"  
 
{| 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>
+
| <math>{\left. \frac{\partial T\left( {Q}_{g},p,GLR,K\right) }{\partial x}\right| }_{x=\Delta x}=G{T}_{wf}</math>
 
|}
 
|}
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
 
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
 
|}
 
|}
  
 +
where <math display="inline">G{T}_{wf}</math> is the average temperature gradient, m/<math display="inline">^\circ\mbox{C}</math>; <math display="inline">{Q}_{g}</math> is the reasonable distribution of production for gas wells, <math display="inline">10^4 \mbox{m}^3/d</math>; <math display="inline">GLR</math> is the gas and liquid ratio, m<math display="inline">^3</math>/m<math display="inline">^3</math>; <math display="inline">K</math> is the total heat transfer coefficient of well bore which can be obtained according to the well structure, W/(m<math display="inline">^2</math>.K).
  
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<math display="inline">^\circ\mbox{C}</math>/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.
 
+
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===
 
===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,
+
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.
+
# The production and bottom hole pressure of gas well are determined according to the reasonable production distribution method of gas well.
 +
# The number of downhole throttles is determined preliminarily according to the pressure of surface gathering system.
 +
# 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.
 +
# Calculate the temperature and pressure upstream of the first downhole throttle according to the well bore dynamic analysis method.
 +
# Call the downhole throttle pressure and temperature drop model to calculate the pressure and temperature downstream of the first-stage throttle.
 +
# Calculate hydrate formation temperature, <math display="inline">T_{h1}</math>, under downstream pressure of the throttle, if the temperature behind the nozzle (<math display="inline">T_2 > T_{h1} + 3^\circ\mbox{C}</math>), 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.
 +
# 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.
 +
# 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.
 +
# 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==
 
==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.
+
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 [[#tab-1|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;">
+
<div class="center" style="font-size: 75%;">'''Table 1'''.  Application of multi-stage downhole choke</div>
Table 1 Application of multi-stage downhole choke</div>
+
  
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;"  
+
<div id='tab-1'></div>
 +
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"  
 +
|-style="text-align:center"
 +
! No. !! Production <br> (10<sup>4</sup>m<sup>3</sup>/d) !! Mean annual temperature<br>  (℃) !! Choke series !! Choke depth<br> (m) !! Pressure before choke<br> (MPa) !! Temperature before choke<br> (℃) !! Pressure  after choke <br> (MPa) !!Temperature  after choke <br> (℃) !! Hydrate formation temperature <br> (℃)
 
|-
 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">No.</span>
+
rowspan='2' style="text-align: center;"|SuX1
|  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>
+
rowspan='2' style="text-align: center;"|7.00
|  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>
+
rowspan='2' style="text-align: center;"|8
|  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="text-align: center;"|1<sup>st</sup> level
|  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="text-align: center;"|2501
|  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="text-align: center;"|20.88
|  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="text-align: center;"|89.45
|  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="text-align: center;"|11.40
|  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="text-align: center;"|42.21
|  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>
+
|  style="text-align: center;"|22.23
 
|-
 
|-
rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">SuX1</span>
+
|  style="text-align: center;"|2<sup>nd</sup> level
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.00 </span>
+
|  style="text-align: center;"|1999
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
+
|  style="text-align: center;"|10.93
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1<sup>st</sup> level</span>
+
|  style="text-align: center;"|74.89  
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2501 </span>
+
|  style="text-align: center;"|5.97
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.88 </span>
+
|  style="text-align: center;"|29.55
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">89.45 </span>
+
|  style="text-align: center;"|18.34
|  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>
+
rowspan='2' style="text-align: center;"|SuX2
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1999 </span>
+
|  rowspan='2' style="text-align: center;"|6.00
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">10.93 </span>
+
rowspan='2' style="text-align: center;"|8
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">74.89 </span>
+
style="text-align: center;"|1<sup>st</sup> level
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.97 </span>
+
|  style="text-align: center;"|2500
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">29.55 </span>
+
style="text-align: center;"|17.99
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">18.34 </span>
+
|  style="text-align: center;"|87.72
 +
|  style="text-align: center;"|9.82
 +
|  style="text-align: center;"|40.70
 +
|  style="text-align: center;"|21.34  
 
|-
 
|-
rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">SuX2</span>
+
|  style="text-align: center;"|2<sup>nd</sup> level
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.00 </span>
+
|  style="text-align: center;"|2005
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
+
|  style="text-align: center;"|9.43
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1<sup>st</sup> level</span>
+
|  style="text-align: center;"|73.36
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2500 </span>
+
|  style="text-align: center;"|5.15
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.99 </span>
+
|  style="text-align: center;"|28.22
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">87.72 </span>
+
|  style="text-align: center;"|17.45
|  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>
+
rowspan='2' style="text-align: center;"|SuX3
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2005 </span>
+
|  rowspan='2' style="text-align: center;"|10.00
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9.43 </span>
+
rowspan='2' style="text-align: center;"|8
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">73.36 </span>
+
style="text-align: center;"|1<sup>st</sup> level
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.15 </span>
+
|  style="text-align: center;"|2408
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">28.22 </span>
+
style="text-align: center;"|20.77
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.45 </span>
+
|  style="text-align: center;"|90.24
 +
|  style="text-align: center;"|11.34
 +
|  style="text-align: center;"|42.90
 +
|  style="text-align: center;"|22.20
 
|-
 
|-
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>
+
|  style="text-align: center;"|2<sup>nd</sup> level
|  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>
+
|  style="text-align: center;"|2008  
|  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="text-align: center;"|10.97  
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1<sup>st</sup> level</span>
+
|  style="text-align: center;"|78.64  
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2408 </span>
+
|  style="text-align: center;"|5.99  
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.77 </span>
+
|  style="text-align: center;"|32.81  
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">90.24 </span>
+
|  style="text-align: center;"|18.37  
|  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==
 
==8. Conclusion==
Line 296: Line 279:
 
==References==
 
==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.
+
<div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;">
 +
 
 +
[1] Duncan K.W., Musaffar A.K.,  Anas M.Q. 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, 2011.
  
[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.
+
[2] Orta D., Ramanchandran S., Yang J. et al. 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, 2007.
  
[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.
+
[3] Siddiqui S.,  Yang J. 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, 1999.
  
[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)''':37-42.
+
[4] Wang F., Liang Z., Deng X. Automation of liquid removal from natural gas wells. Proceedings of Institution of Civil Engineers - Energy, 170(1):37-42, 2017.
  
[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.
+
[5] William H. 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, 2010.
  
[6] Xie C, Liu Y, Li X (2022) A novel comprehensive model for predicting production of downhole choke wells. Fuel '''313''':1-23.
+
[6] Xie C., Liu Y., Li X. A novel comprehensive model for predicting production of downhole choke wells. Fuel, 313:1-23, 2022.
  
[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.
+
[7] Zheng L., Fei W., Xiong D. 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, 2014.

Latest revision as of 14:25, 8 June 2022

Abstract

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.

Keywords: 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 , is the non-critical state. Among them, is the pressure in front of the throttle, is the pressure at the end of the throttle, is the critical pressure after the throttle. Downstream of the throttle, the gas velocity decreases with the increase of , and the downstream pressure is equal to the back pressure.
  2. When , Is the critical state, Downstream of the throttle, the gas velocity reaches its maximum and the downstream pressure is equal to the critical pressure . 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.

Draft Yijun 989359957-image1.png
Figure 1. Schematic diagram of well bore pressure and temperature variation


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,

(1)

where, is natural gas production (under standard conditions), Nm3/d; and are the the upstream and downstream pressures of the throttle, respectively, MPa; is the diameter of the orifice, mm; is the upstream inlet temperature of the throttle, ; is the gas compressibility coefficient at the inlet of the throttle; is the relative density of natural gas; is the natural gas adiabatic index.

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

(2)


When the output of natural gas, the temperature and pressure at the upstream inlet of the throttle are known, substituting them into Eq.(1) and calculating after deformation, the pressure after throttling can be determined, the pressure drop can be obtained as .

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:

(3)

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

(4)

where, is water cut correction coefficient, and 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,

(5)

where, and are the gas compressibility coefficients at the inlet and outlet conditions of the throttle respectively, and 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,

(6)

where, is the bottom hole pressure under reasonable production allocation of gas wells, MPa; is the highest safety gas gathering pressure at the wellhead of gas well, MPa; is the hydrate formation temperature corresponding to the highest safe gas gathering pressure, ; is well bore pressure loss, MPa; is the throttling pressure drop of the downhole throttling device, MPa; is the number of downhole throttles. 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,

(7)

where 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

(8)

where 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 when the temperature after the throttles returns to the geothermal gradient as the lowering position of the next-level throttle, namely

(9)

where is the average temperature gradient, m/; is the reasonable distribution of production for gas wells, ; is the gas and liquid ratio, m/m; is the total heat transfer coefficient of well bore which can be obtained according to the well structure, W/(m.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, , under downstream pressure of the throttle, if the temperature behind the nozzle (), 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.

Table 1. Application of multi-stage downhole choke
No. Production
(104m3/d)
Mean annual temperature
(℃)
Choke series Choke depth
(m)
Pressure before choke
(MPa)
Temperature before choke
(℃)
Pressure after choke
(MPa)
Temperature after choke
(℃)
Hydrate formation temperature
(℃)
SuX1 7.00 8 1st level 2501 20.88 89.45 11.40 42.21 22.23
2nd level 1999 10.93 74.89 5.97 29.55 18.34
SuX2 6.00 8 1st level 2500 17.99 87.72 9.82 40.70 21.34
2nd level 2005 9.43 73.36 5.15 28.22 17.45
SuX3 10.00 8 1st level 2408 20.77 90.24 11.34 42.90 22.20
2nd level 2008 10.97 78.64 5.99 32.81 18.37

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.

References

[1] Duncan K.W., Musaffar A.K., Anas M.Q. 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, 2011.

[2] Orta D., Ramanchandran S., Yang J. et al. 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, 2007.

[3] Siddiqui S., Yang J. 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, 1999.

[4] Wang F., Liang Z., Deng X. Automation of liquid removal from natural gas wells. Proceedings of Institution of Civil Engineers - Energy, 170(1):37-42, 2017.

[5] William H. 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, 2010.

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

[7] Zheng L., Fei W., Xiong D. 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, 2014.
Back to Top

Document information

Published on 09/06/22
Accepted on 07/06/22
Submitted on 06/04/22

Volume 38, Issue 2, 2022
DOI: 10.23967/j.rimni.2022.06.001
Licence: CC BY-NC-SA license

Document Score

0

Views 51
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?