1 Prediction of Hydrate Plugs in Gas Wells in Permafrost Edward Bondarev Leading Researcher Kira Argunova Senior Researcher Igor Rozhin Senior Researcher Institute of Oil and Gas Problems, Yakutsk, Russia Abstract An approach to predictions of position and size of hydrate plugs inside gas wells has been proposed. It is based on the mathematical model of steady non-isothermal flow of real gas in tubes and an algorithm of calculation of equilibrium conditions of hydrate formation. The proposed approach includes the following steps. 1) Numerically solve the system of ordinary differential equations to find the distributions of pressure and temperature along a particular well. 2) Represent the results of calculations as connection between pressure and temperature. 3) Find the intersection of this function with the calculated or experimental equilibrium curve for a particular natural gas. 4) Find the depth of well from the results of numerical solution. Keywords: hydrate; gas well; optimal regime; gas extraction; computational experiment. Introduction Analysis of factors determining reliability of gas supply to consumers located in the permafrost zone show that the first weak link in the technological chain is a well itself and the adjoining bottom-hole zone of the gas reservoir. Just here the intensive cooling of gas occurs due to throttle effect and heat exchange with the permafrost rocks surrounding the well. Due to the fact that many fields have sufficiently high reservoir pressures the danger of gas hydrate formation directly in the wells arises, which can cause either lowering of flow rate or their complete plugging. Moreover, risk of well plugging with gas hydrates can arise at their stoppage because of low temperature of surrounding rocks. To prevent hydrate formation in wells it is necessary to establish such regime of gas extraction at which its temperature will be higher than the equilibrium temperature of hydrate formation. Such regime can be provided for at the projecting stage of field development, constructive parameters of wells being chosen accordingly. Evidently, it is not always possible because of different engineering constraints. In this connection study of the possibility of controlling the gas temperature without changing the constructive parameters of wells is of a significant interest. Such possibility is based on the following peculiarities of gas well temperature conditions, which are influenced by two factors: external heat transfer and internal dissipative thermodynamic processes being the consequence of gas imperfection. Let us estimate the effect of mass flow rate on each of these factors n considering temperature at a fixed cross section of a well 1xx at two mass flow rates 12 MM and two input temperatures eTT 10 and eTT 20 (where eT is the ambient temperature). In case when gas is perfect the flow rate increase at eTT 10 causes temperature growth in the cross section 1xx while at eTT 20 temperature drops. One can estimate the influence of gas flow rate on dissipative processes by modeling the real gas flow in the heat-insulated pipe. For all cases the increase in mass flow rate results in temperature drop at the fixed cross section. Therefore at eTT 10 the increase in mass flow rate leads to two contrary trends: 1) gas temperature grows due to the decrease in heat transfer while 2) it drops because of throttling. At low flow rates the first trend prevails while at higher ones – the second. Therefore if one plots the curve of temperature at the fixed cross section versus mass flow rate it will have maximum at some optimum value *M . Formulation of the Problem To make the quantitative estimation of this effect let us use the mathematical model of stationary gas flow in pipes derived in the monograph (Bondarev et al. 1988) where it was shown that at the normal operation of gas velocity of gas flow is much lower than that of the sound. For example, at mass flow rate of 10 kg/s and at well depth of 3000 m it is approximately 5 m/s. This allows reducing the equation of flow and energy to the following: 5.224s in Mgdxdp , (1) s inpep cgTTMc DdxdpdxdT , (2) where – gas density; pс – gas specific heat capacity; g – gravity acceleration; D, – cross section and diameter of a pipe; x – coordinate along the pipe axis; p – pressure; – angle of the pipe slope counted from the fixed horizontal plane; – coefficient of hydraulic resistance; 2 NINTH INTERNATIONAL CONFERENCE ON PERMAFROST T – gas temperature; – total coefficient of heat transfer; M – mass flow rate of a gas being a constant. Density is connected with pressure and temperature by the equation of state zRTp, (3) Coefficient of throttling is determined by the formulaTzpcRTp 2 ; coefficient of gas compressibility z is the empirical function of pressure and temperature; R – gas constant. As stated above, one of the possibilities of gas temperature control in a well is based on its non-monotone dependence on mass flow rate. Let us note that similar dependence on the production rate will not have such peculiarity. This phenomenon is explained by the structure of the equation (2). Indeed when one transfers the second member from the left to the right part of the equation and notices the dependence of the pressure gradient on flow rate it becomes clear that the intensity of gas heat exchange with the environment (the first member in the right part of the equation (2)) is inversely proportional to mass flow rate while that of throttling is directly proportional to the square of this magnitude. Therefore, to determine the possibility of controlling the gas temperature it is necessary to find its dependence on M . Let us formulate the initial conditions for the system of equations (1)-(2) in the form: .)0(,)0( 00 TTpp (4) In case when the Bertolou equation of state is used, coefficient 211 TbTpkz, where 2026TTb c, 0007.0 Tp pTkcc, cT , cp – critical temperature and pressure which depend on the natural gas composition. Results For the first variant calculations have been made at the following values of the parameters: 82.5 W/(m2K); 1.0D m; 50 10290 p N/m2; 3230T K; 520R J/(kgK); 2300pc J/(kgK); 191сT K; 5105.45 сp N/m2; 02.0 ; xTT ee 0 ; 3230 eT K; 02.0 K/m; 2900L m. In the second variant – 50 10290 p N/m2; 3000L m; 3 0 0 02 5 9 3 K, 15.271 2 5 9 30,0 xxxTT ee , at the same values of the rest parameters. Thus, the first variant corresponds to the standard well, while the second one – to the well drilled in permafrost. Dependence of gas temperature at the wellhead on mass flow rate is presented in Fig. 1. As it has been expected gas temperature in the first variant (curve 1) is higher than that in the second variant (curve 2). In the case considered this difference makes approximately 12%. One must note that values of mass flow rates corresponding to the maximum temperature of gas at the wellhead practically do not depend on the temperature of rocks. The earlier investigations (Argunova & Bondarev 2005) show that depth and initial values of temperature and pressure have the basic effect on this parameter. The analysis of curves shows that the optimum conditions of gas extraction correspond to mass flow rate 5*M kg/s and temperatures at the wellhead 275.2 K (the first variant) and 273 K (the second variant). Fig. 1. Dependence of dimensionless temperature at the wellhead on mass flow rate (bold numbers at the curves correspond to variants of calculation). Similar dependences for the pressure at the wellhead are presented in Fig. 2. The optimum conditions of gas extraction correspond to 5.2M kg/s and pressures 216.3105 N/m2(the first variant) and 212.0105 N/m2(the second variant). Here non-monotone character of curves is explained by the fact that at low values of M the second member in the equation of motion (2) is very small and does not influence pressure distribution along the depth of the well. This distribution is determined only by the first member in this equation, which decreases with the temperature growth according to the equation of state (3). Due to the fact that at low flow rates the gas temperature increases with the well production rate (see Fig. 1) pressure in the wellhead will increase also until hydraulic resistance, described by the second member in the equation (1), becomes predominant. Thus, the above peculiarity of the wellhead pressure dependence on the well production rate reveals itself only at the non-isothermal regime of gas flow. EDWARD BONDAREV ET AL. 3 One must note particularly that the optimal regime of gas extraction providing the minimum loss of pressure corresponds to much less value of mass flow rate than in the case of the regime with the minimum heat loss. Fig. 2. Dependence of dimensionless pressure in the wellhead on gas mass flow rate (bold numbers at the curves correspond to variants of calculation). The presence of two different extrema indicates the possibility to state the problem of searching such gas extraction regime, which corresponds to the minimum loss of the full energy of gas flow. This energy is expended for resisting the gas friction against pipe walls and heat loss due to throttling and heat transfer to the environment. Comparative analysis of the results of these two variants allows drawing the following conclusion. The presence of the permafrost strata and slight increase in the well depth and bottom-hole pressure do not influence the location of temperature and pressure extrema. Yet pressure and temperature values at the wellhead decrease. To determine the danger of hydrate formation in the borehole one can use the following procedure suggested in (Argunova & Bondarev 2005). On the equilibrium curve of hydrate formation plotted according to the experimental data one lays on the dependence between pressure and temperature obtained by solving the problem (1) - (4). The coordinate x, above which gas will be cooled below the equilibrium temperature, is determined by the temperature corresponding to the intersection of these two curves. The example of using the above procedure is presented in Figs. 3, 5 (mass flow rate corresponds to the optimum temperature of a gas) and Figs. 4, 6 (mass flow rate corresponds to the optimum pressure of a gas). All calculations have been made for the second variant of the initial data. Fig. 3. Intersection of the experimental equilibrium curve of hydrate formation (dots and dotted line) and calculated connection between gas pressure and temperature in the well for M = 5 kg/s. 4 NINTH INTERNATIONAL CONFERENCE ON PERMAFROST Fig. 4. The same as in Fig. 3 for M = 2.5 kg/s. Fig. 5. Dependence of temperature and pressure on the dimensionless depth of a well for M = 5 kg/s. Fig. 6. The same as in Fig. 5 for M = 2.5 kg/s. EDWARD BONDAREV ET AL. 5 The data presented show that in those sections of wells, where depth changes from 0 to 1299300057.01 m (Fig. 5) and to 1330300056.01 m (Fig. 6), gas temperature will be lower than the equilibrium temperature of hydrate formation. It is seen that these depths are essentially larger than the permafrost thickness (the dotted vertical line in Figs. 5, 6). It should be noted in particular that for the variant where mass flow rate corresponds to the optimal gas pressure, the depth of the dangerous zone is larger than that in the variant where mass flow rate corresponds to the optimal gas temperature. The examples presented visually demonstrate possibilities of mathematical simulation and, particularly, computational experiment for solving the actual problems of improving the reliability of the gas production systems in the Northern regions. References Argunova, K.K. & Bondarev, E.A. 2005. Control of gas well operation: possibilities of mathematical simulation. Science and Education, 1: 41-45. Bondarev, E.A., Vasiliev, V.I., Voevodin, A.F., et al. 1988. Thermohydrodynamics of the systems of gas production and transport. Novosibirsk: Nauka. Siberian Branch, 272 pp.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- International Conference on Gas Hydrates (ICGH) (6th : 2008) /
- Prediction of Hydrate Plugs in Gas Wells in Permafrost
Open Collections
International Conference on Gas Hydrates (ICGH) (6th : 2008)
Prediction of Hydrate Plugs in Gas Wells in Permafrost Bondarev, Edward; Argunova, Kira; Rozhin, Igor 2008-07-31
pdf
Page Metadata
Item Metadata
Title | Prediction of Hydrate Plugs in Gas Wells in Permafrost |
Creator |
Bondarev, Edward Argunova, Kira Rozhin, Igor |
Contributor | University of British Columbia. Department of Chemical and Biological Engineering International Conference on Gas Hydrates (6th : 2008 : Vancouver, B.C.) |
Date Issued | 2008-07 |
Description | An approach to predictions of position and size of hydrate plugs inside gas wells has been proposed. It is based on the mathematical model of steady non-isothermal flow of real gas in tubes and an algorithm of calculation of equilibrium conditions of hydrate formation. The proposed approach includes the following steps. 1) Numerically solve the system of ordinary differential equations to find the distributions of pressure and temperature along a particular well. 2) Represent the results of calculations as connection between pressure and temperature. 3) Find the intersection of this function with the calculated or experimental equilibrium curve for a particular natural gas. 4) Find the depth of well from the results of numerical solution. |
Extent | 135415 bytes |
Subject |
hydrate gas well optimal regime gas extraction computational experiment International Conference on Gas Hydrates 2008 ICGH 2008 |
Genre |
Conference Paper |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2008-08-15 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0041086 |
URI | http://hdl.handle.net/2429/1384 |
Affiliation |
Non UBC |
Citation | Bondarev, Edward; Argunova, Kira; Rozhin, Igor. 2008. Prediction of Hydrate Plugs in Gas Wells in Permafrost. Proceedings of the 6th International Conference on Gas Hydrates (ICGH 2008), Vancouver, British Columbia, CANADA, July 6-10, 2008. |
Peer Review Status | Unreviewed |
Aggregated Source Repository | DSpace |
Download
- Media
- 59278-5315.pdf [ 132.24kB ]
- Metadata
- JSON: 59278-1.0041086.json
- JSON-LD: 59278-1.0041086-ld.json
- RDF/XML (Pretty): 59278-1.0041086-rdf.xml
- RDF/JSON: 59278-1.0041086-rdf.json
- Turtle: 59278-1.0041086-turtle.txt
- N-Triples: 59278-1.0041086-rdf-ntriples.txt
- Original Record: 59278-1.0041086-source.json
- Full Text
- 59278-1.0041086-fulltext.txt
- Citation
- 59278-1.0041086.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.59278.1-0041086/manifest