Journal of Applied Mathematics

Volume 2014, Article ID 542343, 11 pages

http://dx.doi.org/10.1155/2014/542343

## Numerical Investigation of Gas Mixture Length of Nitrogen Replacement in Large-Diameter Natural Gas Pipeline without Isolator

^{1}State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu, Sichuan 610500, China^{2}Petroleum Engineering Institute, Southwest Petroleum University, Chengdu, Sichuan 610500, China

Received 6 December 2013; Accepted 15 February 2014; Published 25 March 2014

Academic Editor: M. Montaz Ali

Copyright © 2014 Hongjun Zhu and Qinghua Han. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Nitrogen replacement is a key process for natural gas pipeline before it is put into operation. A computational fluid dynamic model coupled to a species-transportation model has been used to investigate the gas mixture length of nitrogen replacement in large-diameter pipeline without isolator. A series of numerical simulations are performed over a range of conditions, including pipe length and diameter, inlet rate, and inclination angle of pipe. These affecting factors are analyzed in detail in terms of volume fraction of nitrogen, the maximum gas mixture length, and gas mixture length varied with time. Gas mixture length increases over time, and the maximum gas mixture length is present at outlet of pipe. Long and large-diameter pipe and fast speed of nitrogen lead to long length of mixed gas, while large inclination angle of pipe brings about short length. Several fitting formulas have been obtained, which can predict the maximum gas mixture length in gas pipelines. The used method of fitting formula is shown in the paper by examples. The results provide effective guidance for practical operation of nitrogen replacement.

#### 1. Introduction

Facing the increasing gas consumption, natural gas pipeline is progressing towards large diameter and long distance, resulting in high investment. It is very important to ensure the safety of gas pipeline operation. However, explosion is easily caused when the content of natural gas in air reaches 5~15% [1, 2]. Therefore, before the natural gas pipeline is put into operation, the air in pipe must be replaced by an inert gas. In practice, nitrogen, as the most readily available and least expensive inert gas, is usually used to replace air in pipe. And qualified replacement is performed when the content of air in pipe is less than 2% [3].

To ensure safety, isolator is commonly employed to separate the nitrogen and air. As shown in Figure 1(a), before pumping the nitrogen into pipe, isolator is firstly placed at the entrance of pipe. Then nitrogen flow pushes the isolator moving forward together [4]. However, a number of natural gas pipelines are laid in mountains and plateaus. The undulating terrain makes the passing through of isolator more difficult. And in larger diameter, the cost of isolator is doubled up. In addition, the isolator seal failure usually occurs due to wear. Therefore, nitrogen replacement without isolator is popular in today’s natural gas pipeline engineering as shown in Figure 1(b). After all, the receiving and transmitting cost of isolator can be saved. In this process, the gas mixture length becomes a key parameter to determine the displacement effect. According to the gas mixture length, appropriate operating parameters such as pumping volume and forward speed of nitrogen can be defined. However, in China, the amount of nitrogen is mainly determined by experience with great blindness [5–7]. And the mixture rule of nitrogen and air and affecting factors on mixture length are still unclear.

Several works in the literature focus on nitrogen injection to develop gas and gas condensate fields [8–10]. However, the mixture rule in formation pore is different with it in large-diameter pipeline, due to the great difference in scales. Therefore, an in-depth research for gas mixture length of nitrogen replacement in large-diameter pipeline without isolator is urgently needed in natural gas pipeline engineering.

Since numerical simulation can provide detailed information of flow field which is not easily obtained by physical experiments and has advantages of low cost and short research time, in present work, CFD model coupling with a species-transportation model has been employed to investigate the gas mixture length of nitrogen replacement in large-diameter pipeline without isolator. The gas mixture length has been examined in terms of volume fraction of nitrogen. By conducting a series of numerical simulations, effects of pipe length, pipe diameter, inlet rate, and inclination angle of pipe are examined. Then, the maximum gas mixture length is analyzed and the fitting formulas are obtained. Using these formulas, we can predict the maximum gas mixture length in gas pipelines. The used method of fitting formula is shown in the paper by examples. The results provide useful guidance for practical operation of nitrogen replacement.

The remaining part of this paper is organized as follows. In Section 2, the description of simulation problem is provided; Section 3 presents the governing equations and numerical method; Section 4 presents the numerical results and discussion; Section 5 is the concluding remarks.

#### 2. Problem Description

Straight pipe and undulating pipe with an inclined upward section are adopted in this study. Figure 2 shows a sketch of the geometry and numerical grid for computational domain. The straight pipe is placed horizontally. And the length of it () ranges from 50 m to 1000 m in comparing cases, in order to examine the effect of pipe length on gas mixture length, while the diameter of straight pipe is set to 346.0 mm, 647.2 mm, 851.2 mm, and 1000.2 mm to explore the effect of pipe diameter.

The undulating pipe consists of three sections: straight pipe section, bend section, and inclined upward section. The bend curvature (*R*/*D*) is defined as 3, the common one used in practice. The length of inclined upward section is fixed at 600 m, while the inclination angle changes from 10° to 40° in comparing cases conducted to observe the effect of pipe inclination.

Gas flow in pipeline is a symmetric problem, so two-dimensional flow simulation is accurate enough to capture the gas mixture length. In addition, three-dimensional simulation needs a higher CPU cost. Due to time limitations, 2D simulation is applied in this work. All geometry generation and meshing are performed using GAMBIT 2.3 mesh-generator. As shown in Figure 2, quadrilateral grids are used in the whole computational domain of both straight pipe and undulating pipe. And progressive mesh is used near the pipe wall. The computational domain of undulating pipe is divided into three blocks. Denser mesh is employed in the bend section. A suitable grid density is reached by repeating computations until a satisfactory independent grid is found. For example, the number of quadrilateral cells for straight pipe with diameter of 647.2 mm and length of 600 m is 480000 at last.

The velocity of gas flow is slow in pipe (2~5 m/s), so both nitrogen and air can be seen as incompressible fluids. In simulations, density and viscosity of nitrogen are defined as 1.138 kg/m^{3} and Pa*·*s, respectively. And the density and viscosity of air are set to 1.225 kg/m^{3} and Pa*·*s, respectively.

#### 3. Governing Equations and Numerical Method

##### 3.1. Governing Equations

The gas flows of nitrogen and air are governed by the Reynolds-Averaged-Navier-Stokes (RANS) equations, including continuity and momentum equations written as follows [11]: where where represents instantaneous velocity component in direction, for example, and are velocity in and direction, respectively, while is fluctuation velocity component in direction, is space coordinate in direction, is gravitational acceleration in direction, is time, is pressure, is density of mixed gas, and are density of nitrogen and air, respectively, is kinematic viscosity of mixed gas, and and are viscosity of nitrogen and air, respectively.

Reynolds number ranges from to in calculations. Thus, a realizable turbulence model [12–14] is employed to close the flow governing equations and describe the turbulent properties: where where and represent turbulent kinetic energy and turbulent kinetic energy dissipation rate per unit mass, respectively, is the relative strain parameter, is the strain rate, and represent production term of turbulent kinetic energy due to the average velocity gradient and production term of turbulent kinetic energy due to lift, respectively, is turbulent viscosity, is Prandtl number taken as 0.85, , , , and are empirical model constants taken as 0.09, 1.44, 1.9, and 0.9, respectively, and and are turbulent Prandtl numbers taken as 1.0 and 1.2, respectively.

There is no reaction between nitrogen and air in the replacement process, so a simplified species-transportation model is employed to capture the mixture of nitrogen and air written as [15] follows: where is nitrogen volume fraction and and are diffusion coefficient of nitrogen and air, respectively.

##### 3.2. Numerical Method

Finite volume method (FVM) is employed to discretize above equations. All the calculations are performed using a commercial software package FLUENT 14.5. Patankar’s well-known SIMPLE algorithm [16] is adopted to solve the pressure-velocity coupling. In order to ensure the accuracy of calculation, second-order upwind scheme and second-order central-differencing scheme are used for convective terms and diffusion terms, respectively. The convergent criteria for all calculations are set as that the residual in the control volume for each equation is smaller than 10^{−5}.

##### 3.3. Boundary and Initial Conditions

Velocity inlet boundary condition is used for the inlet, and the inlet rates of nitrogen are taken as 2 m/s, 3 m/s, 4 m/s, and 5 m/s in comparing cases. In the outlet of computational domain, pressure outlet boundary condition is employed, and the value is defined as 0 Pa in order to facilitate comparative analysis. No slip boundary condition is imposed on the pipe wall.

Since it is an unsteady problem, the whole computational domain of pipe is defined as being filled with air at the initial time. And the time step in simulations is set to 0.001 s. The information of simulation cases is listed in Table 1, in which pipe length, pipe diameter, and inlet rate and inclination angle are the four variables.

#### 4. Numerical Results and Discussion

##### 4.1. Standard Case Analysis

Case 6 is adopted as the standard case. Figure 3 shows the gas mixture length at different times, which is reflected by the volume fraction of nitrogen. The volume fraction of nitrogen more than 98% can be considered to meet the safe replacement requirement. And the content of nitrogen in air is 79%. Therefore, volume fraction of nitrogen ranging from 0.79 to 0.98 corresponds to the mixed gas. And the length of this mixed gas is called gas mixture length.

It is seen that the head of the mixed gas shows a bullet-shaped distribution. Frictional resistance at pipe wall is the major cause for this phenomenon. It is well known that there is a viscous sublayer near pipe wall. The closer the gas is to the wall, the greater the viscosity resistance gas gets. So the gas in the axis of pipe moves ahead. It is worth noting that the gas mixture length increases over time. At the time of 50 s, the gas mixture length is 9.5 m, while it increases to 18.5 m at the time of 200 s. In the mixed gas, the volume fraction of nitrogen gradually decreases along the pipe axis direction. And the decreasing rate is gradually reduced with time. This is because nitrogen is mixed with air more fully as time passes. However, the growth rate of mixture length is reducing over time. In the period from 50 s to 100 s, the mixture length is increased by 4 m, while it is increased by 3.5 m in the same interval of time (from 100 s to 150 s). And the increment is just 1.8 m for time increased from 150 s to 200 s. It can be explained that kinetic energy of nitrogen flow is consumed continuously as time goes. Therefore, the driving force and forward speed decrease gradually.

##### 4.2. Effect of Pipe Length

Since gas mixture length increases over time, the maximum mixture length will appear at the outlet of pipe. Figure 4 illustrates the maximum gas mixture length for different pipe lengths. It is observed that the maximum gas mixture length increases with the increase in pipe length. The maximum gas mixture length reaches 24 m in 1000 m long pipe, which is about six times longer than that in 50 m long pipe. The main reason is that the longer the pipe is, the more the time is required for gas flow passing through, resulting in more gas mixed. However, the increment of the maximum gas mixture length reduces for pipe with length ranging from 400 m to 1000 m. And the increment is not decreased linearly.

The gas mixture length curves, as shown in Figure 5, also show nonlinear relationship. Gas mixture length increases rapidly at the beginning, while the growth rate decreases after a period of time. This main reason is that, at the beginning of replacement, the volume fraction gradient and velocity gradient in pipe are both large, but the two gradients decrease gradually over time.

The longest pipeline used in simulation is just 1000 m. However, actual piping often has hundreds or thousands of kilometers. In order to predict the maximum gas mixture length in longer pipelines, the maximum gas mixture length () versus the pipe length () is analyzed, and we have obtained the fitting formula, as shown in Figure 6. The adjusted -squares reach 0.98722. Namely, the error in mixing length prediction is of the order of 1~2%. The fitting formula is written as

If the actual parameters are the same as that used in this paper such as mm and m/s, we can calculate the maximum gas mixture length in an arbitrary long pipe using this formula. For example, if the length of pipe is m, would be calculated as 1998.13 m. In order to detect the accuracy of this prediction equation, 1D simulation for pipe length m has been conducted. One order central difference method is employed to discrete one-dimensional governing equations. For m long pipe with diameter 647.2 mm and inlet rate 3 m/s, the calculated maximum mixing length by 1D simulation is 1959.62 m. So the error in mixing length prediction is less than 2%.

Therefore, the actual replacement time is not simply equal to pipe length divided by the flow speed. Qualified replacement is achieved until the mixed gas flows out thoroughly. So the replacement time should be more than the sum of pipe length divided by flow velocity and the maximum gas mixture length divided by flow rate.

##### 4.3. Effect of Pipe Diameter

Figure 7 presents the maximum gas mixture length for different pipe diameters. It can be seen that the maximum gas mixture length increases with the increase in pipe diameter. The main cause of this result is that the total kinetic energy of nitrogen flow is more in large-diameter pipe as the same inlet rate. Affected by the impact of friction, the smaller the pipe diameter is, the sharper the bullet-shaped head is.

The gas mixture length varied with time in pipes of different diameters is shown in Figure 8. The growth trend of mixture length is similar for different diameters. However, in small-diameter pipe, the growth rate of mixture length is relatively small. It is attributed to the joint action of kinetic energy and frictional resistance.

The actual pipe diameter is not only the four we have studied. So we also get the fitting formula of versus , as shown in Figure 9. The fitting line meets the linear distribution well, and the adjusted -squares are 0.99295. The fitting formula is written as If the pipe diameter is 800.2 mm, would be calculated as 20.8 m.

##### 4.4. Effect of Inlet Rate

Figure 10 depicts the maximum gas mixture length at different inlet rates. The faster the nitrogen forward, the longer the maximum gas mixture length. It is attributed to that high speed of nitrogen flow has more kinetic energy. Gas mixed more fuller at large inlet rate, resulting in the longer mixture length. As inlet rate increases 1 m/s, the increments of the maximum gas mixture length are basically the same, about 1.9 m.

The gas mixture length varied with time at different inlet rates is shown in Figure 11. The growth rate of mixture length presents the maximum value (0.2 m/s) at m/s, which is about 3.33 times than that at m/s. Less time required for high-speed flow to pass through the pipe is the main cause. Though faster replacement is reached by high inlet rate, the rate cannot be set to too large. Residue particles will be carried by high-speed flow and strike pipe wall, leading to electric spark. Therefore, inlet rate should be controlled in a certain range. In practice, 4 m/s is an appropriate choice.

As shown in Figure 12, the fitting formula of versus is obtained. This fitting line also meets the linear distribution, and the adjusted -squares are 0.99483. The fitting formula is as follow: Using this formula, would be calculated as 21 m if the inlet rate is 3.5 m/s.

##### 4.5. Effect of Inclination Angle of Pipe

Taking the terrain into account, four undulating pipes with different inclination angles are analyzed. Figure 13 provides the maximum gas mixture length at different inclination angles. We find an obvious reduction in mixture length as inclination angle increases. For , the maximum gas mixture length is 26 m, while it decreases to 24.7 m for . The main cause is that partial energy of gas flow is used to overcome gravity acting, and more energy is consumed for steep pipe.

Figure 14 shows the gas mixture length varied with time at different inclination angles. There is no difference in straight pipe section. The significant difference appears in inclined upward section. The larger the inclination angle, the smaller the growth rate of mixture length.

The fitting formula of versus is shown in Figure 15. This fitting line is a negative slope straight line, and the adjusted -squares reach 0.99691. The fitting formula is If the inclination angle is 50°, would be calculated by the formula as 24.25 m.

#### 5. Conclusions

A computational fluid dynamic model coupled to a species-transportation model has been used to investigate the gas mixture length of nitrogen replacement in large-diameter pipeline without isolator. Effects of pipe length and diameter, inlet rate, and inclination angle of pipe are examined by conducting a series of simulations. Based on the numerical results, the following conclusions can be drawn.(1)Gas mixture length increases over time, and the maximum gas mixture length is present at outlet of pipe. Long and large-diameter pipe and fast speed of nitrogen lead to long length of mixed gas, while large inclination angle of pipe brings about short length.(2)Four fitting formulas have been obtained, which can predict the maximum gas mixture length in gas pipelines. Besides that the formula for pipe length is a quadratic polynomial, the other formulas meet linear relationship. The calculation results provide effective guidance for practical operation of nitrogen replacement.

#### Conflict of Interests

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

#### Acknowledgments

Research work was supported by Open Fund (no. PLN1210) of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Southwest Petroleum University) and Key Project of Sichuan Provincial Education Department (No. 12ZA189). Without the support, this work would not have been possible.

#### References

- G. Liu and G. W. Wang, “Method to control gas displacement velocity in gas pipeline commissioning,”
*Oil & Gas Storage and Transportation*, vol. 26, no. 22, pp. 53–56, 2007. View at Google Scholar - W. Duan, J. Zhang, S. Zhang, and H. P. Liu, “Study on length of gas mixture segment for commissioning replacement of long-distance gas pipeline,”
*Oil & Gas Storage and Transportation*, vol. 30, no. 3, pp. 5–7, 2012. View at Google Scholar - L. W. Tan, J. Q. Jing, Z. X. Dai, Z. Lv, and H. M. Shao, “Study of natural gas pipeline replacement technology by nitrogen gas,”
*Pipeline Technique and Equipment*, vol. 15, no. 3, pp. 25–28, 2007. View at Google Scholar - K. F. Fan, W. Q. Wang, Y. Ma, H. Jiang, and Z. X. Pang, “The numerical simulation on mixture gas of nitrogen replacement in natural gas pipeline,”
*Journal of Petrochemical Universities*, vol. 26, no. 1, pp. 63–67, 2013. View at Google Scholar - L. W. Tan, J. Q. Jing, Z. X. Dai, Z. Lv, and H. Zhang, “Determination of process parameters of nitrogen replacement in gas pipeline,”
*Petroleum Planning & Engineering*, vol. 18, no. 5, pp. 43–45, 2007. View at Google Scholar - Y. M. Guo and J. J. Xue, “Numerical simulation on gas mixing rule in course of nitrogen isolation gas displacement technology for commissioning of gas transmission pipeline,”
*Oil Field Equipment*, vol. 35, no. 5, pp. 31–34, 2006. View at Google Scholar - J. Z. Huang and L. Fang, “Analysis of nitrogen replacement in gas pipeline,”
*Petroleum and Chemical Construction*, vol. 32, no. 3, pp. 73–75, 2010. View at Google Scholar - J. E. Hanssen, “Nitrogen as a low-cost replacement of natural gas reinjection offshore,” in
*Proceedings of the SPE Gas Technology Symposium*, SPE 17709, Dallas, Tex, USA, June 1988. View at Scopus - L. C. Berger and J. P. Arnoult, “Production of inert gas for partial replacement of natural gas trapped in an underground aquifer storage reservoir,” in
*Proceedings of the SPE Gas Technology Symposium*, SPE 19089, Dallas, Tex, USA, June 1989. View at Scopus - X. R. Wu, K. G. Liang, and D. X. Liu, “Nitrogen injection experience to development gas and gas condensate fields in rocky mountains,” in
*Proceedings of the International Petroleum Technology Conference*, IPTC 16830, Beijing, China, March 2013. - G. Alfonsi, “Reynolds-averaged Navier-Stokes equations for turbulence modeling,”
*Applied Mechanics Reviews*, vol. 62, no. 4, Article ID 040802, pp. 1–20, 2009. View at Publisher · View at Google Scholar · View at Scopus - K. van Maele and B. Merci, “Application of two buoyancy-modified k-
*ε*turbulence models to different types of buoyant plumes,”*Fire Safety Journal*, vol. 41, no. 2, pp. 122–138, 2006. View at Publisher · View at Google Scholar · View at Scopus - P. Rohdin and B. Moshfegh, “Numerical predictions of indoor climate in large industrial premises. A comparison between different k-
*ε*models supported by field measurements,”*Building and Environment*, vol. 42, no. 11, pp. 3872–3882, 2007. View at Publisher · View at Google Scholar · View at Scopus - M. Ji, T. Utomo, J. Woo, Y. Lee, H. Jeong, and H. Chung, “CFD investigation on the flow structure inside thermo vapor compressor,”
*Energy*, vol. 35, no. 6, pp. 2694–2702, 2010. View at Publisher · View at Google Scholar · View at Scopus - Y. Sun, J. N. Petersen, and T. P. Clement, “Analytical solutions for multiple species reactive transport in multiple dimensions,”
*Journal of Contaminant Hydrology*, vol. 35, no. 4, pp. 429–440, 1999. View at Publisher · View at Google Scholar · View at Scopus - S. V. Patankar,
*Numerical Heat Transfer and Fluid Flow*, McGraw-Hill Press, New York, NY, USA, 1980.