**1. Introduction**

20 Will-be-set-by-IN-TECH

[30] Mayne, W. H. [1962]. Common reflection point horizontal data stacking techniques,

[31] Mitrionics [n.d.]. Languaje Mitrion-C, http://www.mitrionics.com/. Reviewed:

[32] Moore, G. E. [1975]. Progress in digital integrated electronics, *Electron Devices Meeting,*

[33] Necsulescu, P. I. [2011]. *Automatic Generation of Hardware for Custom Instructions*, PhD

[34] Necsulescu, P. I. & Groza, V. [2011]. Automatic Generation of VHDL Hardware Code from Data Flow Graphs, *6th IEEE International Symposium on Applied Computational*

[35] Nios II C-to-Hardware Acceleration Compilern [n.d.]. http://www.altera.com/.

[37] Owens, J. D., Houston, M., Luebke, D., Green, S., Stone, J. E. & Phillips, J. C. [2008]. Gpu

[38] Rieber, F. [1936]. Visual presentation of elastic wave patterns under various structural

[39] Riverside Optimizing Compiler for Configurable Computing: Roccc 2.0 [n.d.]. http://

[40] Sánchez Fernández, R. [2010]. *Compilación C a VHDL de códigos de bucles con reuso de datos*,

[43] Xilinx [n.d.b]. High-Level Synthesis: AutoESL, http://www.xilinx.com/

[44] Yankova, Y., Bertels, K., Vassiliadis, S., Meeuws, R. & Virginia, A. [2007]. Automated HDL Generation: Comparative Evaluation, *2007 IEEE International Symposium on Circuits*

[36] Onifade, A. [2004]. History of the computer, *Conference of History of Electronics*.

www.jacquardcomputing.com/roccc/. Revisado: Junio de 2011.

[41] Stolt, R. H. [1978]. Migration by fourier transform, *Geophysics* (43): 23–48. [42] *Xilinx* [n.d.a]. http://www.xilinx.com/. Reviewed: April 2012.

university/tools/autoesl/. Reviewed: April 2012.

*Geophysics* 27(06): 927–938.

thesis, Ottawa, Canada.

Revisado: Junio de 2011.

*and Systems* pp. 2750–2753.

*1975 International*, Vol. 21, pp. 11–13.

*Intelligence and Informatics* pp. 523–528.

conditions, *Geophysics* 01(02): 196–218.

Tesis, Universidad Politécnica de Cataluña.

computing, *Proceedings of the IEEE* 96(5): 879–899.

April 2012.

The slug is a multiphase flow pattern that occurs in pipelines which connect the wells in seabed to production platforms in the surface in oil industry. It is characterized by irregular flows and surges from the accumulation of gas and liquid in any cross-section of a pipeline. In this work will be addressed the riser slugging, that combined or initiated by terrain slugging is the most serious case of instability in oil/water-dominated systems [5, 15, 21].

The cyclic behavior of the riser slugging, which is illustrated in Figure 1, can be divided into four phases: (i) Formation: gravity causes the liquid to accumulate in the low point in pipeline-riser system and the gas and liquid velocity is low enough to enable for this accumulation; (ii) Production: the liquid blocks the gas flow and a continuous liquid slug is produced in the riser, as long as the hydrostatic head of the liquid in the riser increases faster than the pressure drop over the pipeline-riser system, the slug will continue to grow; (iii) Blowout: when the pressure drop over the riser overcomes the hydrostatic head of the liquid in the riser the slug will be pushed out of the system; (iv) liquid fall back: after the majority of the liquid and the gas has left the riser the velocity of the gas is no longer high enough to drag the liquid upwards, the liquid will start flowing back down the riser and the accumulation of liquid starts over again [15].

The slug flow causes undesired consequences in the whole oil production such as: periods without liquid or gas production into the separator followed by very high liquid and gas rates when the liquid slug is being produced, emergency shutdown of the platform due to the high level of liquid in the separators, floods, corrosion and damages to the equipments of the process, high costs with maintenance. One or all these problems cause significant losses in oil industry. The main one has been of economic order, due to reduction in oil production capacity [6, 8, 16–20].

Currently, control strategies are considered as a promising solution to handle the slug flow [4, 5, 7, 10, 15, 20]. An alternative to the implementation of control strategies is to make use of a mathematical model that represents the dynamic of slug flow in pipeline-separator system.

©2012 Sausen et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0),which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. ©2012 Sausen et al., licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

**Figure 1.** Illustration of a slug cycle.

In this chapter has been used the dynamic model for a pipeline-separator system under the slug flow, with 5 (five) Ordinary Differential Equations (ODEs) coupled, nonlinear, 6 (six) tuning parameters and more than 40 (forty) internal, geometric and transport equations [10, 13], denominated Sausen's model.

To carry out the simulation and implementation of control strategies in the Sausen's model, first it is necessary to calculate its tuning parameters. For this procedure, are used data from a case study performed by [18] in the OLGA commercial multiphase simulator widely used in the oil industry. Next it is important to check how the main variables of the model change their behavior considering a change in the model's tuning parameters. This testing, called sensitivity analysis, is an important tool to the building of the mathematical models, moreover, it provides a better understanding of the dynamic behavior of the system, for later implementation of control strategies.

In this context, from the sensitivity analysis, the Sausen's model has been an appropriate environment for application of the different feedback control strategies in the problem of the slug in oil industries through simulations. The model enables such strategies can be applied in consequence the slug, that is in the oil or gas output valve separator, as well as in their causes, in the top riser valve, or yet in the integrated system, in other words, in more than one valve simultaneously.

Therefore, as part of control strategies that can be used to avoid or minimize the slug flow, this chapter presents the application of the error-squared level control strategy Proportional Integral (PI) in the methodology by bands [10], whose purpose damping of the load flow rate oscillatory that occur in production's separators. This strategy is compared with the level controls strategy PI conventional [1], widely used in industrial processes; and with the level control strategy PI also in the methodology by bands.

The remainder of this chapter is organized as following. Section 2 presents the equation of the Sausen's model for a pipeline-separator system. Section 3 shows the simulation results of the Sausen's model. Section 4 presents the control strategies used to avoid or minimize the slug flow. Section 5 shows the simulation results and analysis of the control strategies applied the Sausen's model. And finally, in Section 6, are discussed the conclusions and future research directions.

### **2. The dynamic model**

#### **2.1. Introduction**

2 Will-be-set-by-IN-TECH

In this chapter has been used the dynamic model for a pipeline-separator system under the slug flow, with 5 (five) Ordinary Differential Equations (ODEs) coupled, nonlinear, 6 (six) tuning parameters and more than 40 (forty) internal, geometric and transport equations [10,

To carry out the simulation and implementation of control strategies in the Sausen's model, first it is necessary to calculate its tuning parameters. For this procedure, are used data from a case study performed by [18] in the OLGA commercial multiphase simulator widely used in the oil industry. Next it is important to check how the main variables of the model change their behavior considering a change in the model's tuning parameters. This testing, called sensitivity analysis, is an important tool to the building of the mathematical models, moreover, it provides a better understanding of the dynamic behavior of the system, for later

In this context, from the sensitivity analysis, the Sausen's model has been an appropriate environment for application of the different feedback control strategies in the problem of the slug in oil industries through simulations. The model enables such strategies can be applied in consequence the slug, that is in the oil or gas output valve separator, as well as in their causes, in the top riser valve, or yet in the integrated system, in other words, in more than one

Therefore, as part of control strategies that can be used to avoid or minimize the slug flow, this chapter presents the application of the error-squared level control strategy Proportional Integral (PI) in the methodology by bands [10], whose purpose damping of the load flow rate oscillatory that occur in production's separators. This strategy is compared with the level controls strategy PI conventional [1], widely used in industrial processes; and with the level

The remainder of this chapter is organized as following. Section 2 presents the equation of the Sausen's model for a pipeline-separator system. Section 3 shows the simulation results of the Sausen's model. Section 4 presents the control strategies used to avoid or minimize the slug

**Figure 1.** Illustration of a slug cycle.

13], denominated Sausen's model.

implementation of control strategies.

control strategy PI also in the methodology by bands.

valve simultaneously.

This section presents a mathematical model for the pipeline-separator system, illustrated in Figure 2 with biphase flow (gas-liquid). The model is the result of coupling the simplified dynamic model of Storkaas [15, 17, 20] with the model for a biphase horizontal cylindrical separator [22]. The new model has been called of Sausen's model.

The following are shown the modelling assumptions, the model equations, how the distribution of liquid and gas inside the pipeline for the separator occurs. Finally, are presented the simulation results for this model considering two settings for simulations: (i) the opening valve *Z* in top of the riser *z* = 20% (flow steady); and (ii) the opening valve *Z* in top of the riser *z* = 50% (slug flow).

**Figure 2.** Illustration of the pipeline-separator system with the slug formation.

#### **2.2. Model assumption**

The Sausen's model assumptions are presented as follow.


#### **2.3. Model equations**

The Sausen's model is composed of 5 (five) ODEs that are based on the mass conservation equations. The equations (1)-(3) represent the dynamics of the pipeline system and the equations (4)-(5) represent the dynamics of the separator:

$$
\dot{M}\_L(t) = m\_{L, \text{in}} - m\_{L, \text{out}}(t) \tag{1}
$$

$$
\dot{M}\_{G1}(t) = m\_{G,in} - m\_{Gint}(t) \tag{2}
$$

$$
\dot{M}\_{G2}(t) = m\_{G1}(t) - m\_{G,out}(t) \tag{3}
$$

$$\dot{N}(t) = \frac{\sqrt{r\_s^2 - (r\_s - N(t))^2}}{2H\_4 \rho\_L N(t) \left[3r\_s - 2N(t)\right]} \left[m\_{L,out}(t) - m\_{LS,out}(t)\right] \tag{4}$$

$$\dot{P}\_{\rm G1}(t) = \frac{\{\rho\_L \Phi[m\_{\rm G,out}(t) - m\_{\rm GS,out}(t)] + P\_{\rm G1}(t) \left[m\_{\rm L,out}(t) - m\_{\rm LS,out}(t)\right]\}}{\rho\_L[V\_{\rm S} - V\_{\rm LS}(t)]} \tag{5}$$

where: *ML*(*t*) is the liquid mass at low point in the pipeline, (*kg*); *MG*1(*t*) is the gas mass in the upstream feed section of pipeline, (*kg*); *MG*2(*t*) is the gas mass at the top of the riser, (*kg*); *N*(*t*) is the liquid level in the separator, (*m*); *PG*1(*t*) is the gas pressure in the separator, (*N*/*m*2); and the *M*˙ *<sup>L</sup>*(*t*), *M*˙ *<sup>G</sup>*1(*t*), *M*˙ *<sup>G</sup>*2(*t*), *N*˙ (*t*), *P*˙ *<sup>G</sup>*1(*t*) are their respective derivatives in relation to time; *mL*,*in* is the liquid mass flowrate that enters the upstream feed section of the pipeline, (*kg*/*s*); *mG*,*in* is the gas mass flowrate that enters in the upstream feed section of the pipeline, (*kg*/*s*); *mL*,*out*(*t*) is the liquid mass flowrate leaving through the valve at the top of the riser enters the separator, (*kg*/*s*); *mG*,*out*(*t*) is the gas mass flowrate leaving through the valve at the top of the riser enters the separator, (*kg*/*s*); *mGint*(*t*) is the internal gas mass flowrate, (*kg*/*s*); *mLS*,*out*(*t*) is the liquid mass flowrate that leaves the separator through the valve *Va*1, (*kg*/*s*); *mGS*,*out*(*t*) is the gas mass flowrate that leaves the separator through the valve *Va*2, (*kg*/*s*); *rs* is the separator ray, (*m*); *H*<sup>4</sup> is the separator length, (*m*); *ρ<sup>L</sup>* is the liquid density, (*kg*/*m*3); *VS* is the separator volume, (*m*3); *VLS*(*t*) is the liquid volume in the separator, (*m*3); Φ = *RT MG* is a constant; *R* is the ideal gas constant (8314 *<sup>J</sup> <sup>K</sup>*.*kmol*); *T* is the temperature, (*K*); *MG* is the gas molecular weight, (*kg*/*kmol*).

The stationary pressure balance over the riser is assumed to be given by

$$P\_1(t) - P\_2(t) = \mathfrak{g}\mathfrak{f}(t)H\_2 - \rho\_L\mathfrak{g}h\_1(t)$$

where: *P*1(*t*) is the gas pressure in the upstream feed section of the pipeline, (*N*/*m*2); *P*2(*t*) is the gas pressure at the top of the riser, (*N*/*m*2); *g* is the gravity (9.81*m*/*s*2); *ρ*¯(*t*) is the average mixture density in the riser, (*kg*/*m*3); *H*<sup>2</sup> is the riser height, (*m*); *h*1(*t*) is the liquid level at the decline, (*m*).

A simplified valve equation is used to describe the flow through the *Z* valve at the top of the riser that is given by

$$m\_{\rm mix,out}(t) = zK\_1 \sqrt{\rho\_T(t)(P\_2(t) - P\_{G1}(t))}\tag{6}$$

where: *z* is the valve position (0 − 100%); *K*<sup>1</sup> is the valve constant and a tuning parameter; *ρT*(*t*) is the density upstream valve, (*kg*/*m*3); *PG*1(*t*) is the gas pressure into the separator, (*N*/*m*2). It is possible to observe that the coupling between the pipeline and the separator occurs through a pressure relationship, in other words, the gas pressure into the separator *PG*1(*t*) is the pressure before the *Z* valve at the top of the riser, according to equation (6).

Considering the result that has been shown in equation (6), it is also possible to obtain the liquid mass flowrate given by

$$m\_{L,out}(t) = \alpha\_L^m(t) m\_{\text{mix},out}(t)$$

and the gas mass flowrate given by

4 Will-be-set-by-IN-TECH

**A5:** Simplified valve equation for gas and liquid mixture leaving the system at the top of the

**A9:** The portion of liquid mixed with the gas in the entrance of the separator is neglected.

The Sausen's model is composed of 5 (five) ODEs that are based on the mass conservation equations. The equations (1)-(3) represent the dynamics of the pipeline system and the

*<sup>G</sup>*1(*t*) = {*ρL*Φ[*mG*,*out*(*t*) <sup>−</sup> *mGS*,*out*(*t*)] + *PG*1(*t*)[*mL*,*out*(*t*) <sup>−</sup> *mLS*,*out*(*t*)]}

where: *ML*(*t*) is the liquid mass at low point in the pipeline, (*kg*); *MG*1(*t*) is the gas mass in the upstream feed section of pipeline, (*kg*); *MG*2(*t*) is the gas mass at the top of the riser, (*kg*); *N*(*t*) is the liquid level in the separator, (*m*); *PG*1(*t*) is the gas pressure in the separator, (*N*/*m*2);

time; *mL*,*in* is the liquid mass flowrate that enters the upstream feed section of the pipeline, (*kg*/*s*); *mG*,*in* is the gas mass flowrate that enters in the upstream feed section of the pipeline, (*kg*/*s*); *mL*,*out*(*t*) is the liquid mass flowrate leaving through the valve at the top of the riser enters the separator, (*kg*/*s*); *mG*,*out*(*t*) is the gas mass flowrate leaving through the valve at the top of the riser enters the separator, (*kg*/*s*); *mGint*(*t*) is the internal gas mass flowrate, (*kg*/*s*); *mLS*,*out*(*t*) is the liquid mass flowrate that leaves the separator through the valve *Va*1, (*kg*/*s*); *mGS*,*out*(*t*) is the gas mass flowrate that leaves the separator through the valve *Va*2, (*kg*/*s*); *rs* is the separator ray, (*m*); *H*<sup>4</sup> is the separator length, (*m*); *ρ<sup>L</sup>* is the liquid density, (*kg*/*m*3); *VS* is the separator volume, (*m*3); *VLS*(*t*) is the liquid volume in the separator, (*m*3); Φ = *RT*

*P*1(*t*) − *P*2(*t*) = *gρ*¯(*t*)*H*<sup>2</sup> − *ρLgh*1(*t*)

*<sup>M</sup>*˙ *<sup>L</sup>*(*t*) = *mL*,*in* <sup>−</sup> *mL*,*out*(*t*) (1)

*<sup>M</sup>*˙ *<sup>G</sup>*1(*t*) = *mG*,*in* <sup>−</sup> *mGint*(*t*) (2) *<sup>M</sup>*˙ *<sup>G</sup>*2(*t*) = *mG*1(*t*) <sup>−</sup> *mG*,*out*(*t*) (3)

[*mL*,*out*(*t*) − *mLS*,*out*(*t*)] (4)

*MG* is

*<sup>ρ</sup>L*[*VS* <sup>−</sup> *VLS*(*t*)] (5)

*<sup>G</sup>*1(*t*) are their respective derivatives in relation to

*<sup>K</sup>*.*kmol*); *T* is the temperature, (*K*); *MG* is the gas

**A6:** Stationary pressure balance over the riser (between pressures *P*1(*t*) and *P*2(*t*)). **A7:** There is not chemical reaction between the fluids (gas-liquid) in pipeline. **A8:** Each one of the fluid consists of a single component in the separator.

**A10:** Simplified valve equation for the gas and the liquid leaving the separator.

*<sup>s</sup>* − (*rs* − *<sup>N</sup>*(*t*))<sup>2</sup> 2*H*4*ρLN*(*t*)[3*rs* − 2*N*(*t*)]

riser.

**A11:** The liquid is incompressible. **A12:** The temperature is constant. **A13:** The gas has ideal behavior.

equations (4)-(5) represent the dynamics of the separator:

 *r*2

*N*˙ (*t*) =

and the *M*˙ *<sup>L</sup>*(*t*), *M*˙ *<sup>G</sup>*1(*t*), *M*˙ *<sup>G</sup>*2(*t*), *N*˙ (*t*), *P*˙

a constant; *R* is the ideal gas constant (8314 *<sup>J</sup>*

The stationary pressure balance over the riser is assumed to be given by

molecular weight, (*kg*/*kmol*).

**2.3. Model equations**

*P*˙

$$m\_{G,out}(t) = [1 - \alpha\_L^m(t)] m\_{mix,out}(t)$$

that leave through the *Z* valve at the top of the riser, where *α<sup>m</sup> <sup>L</sup>* (*t*) is the liquid fraction upstream valve.

The liquid mass flowrate that leaves the separator is represented by the *Va*<sup>1</sup> valve equation given by

$$m\_{LS,out}(t) = z\_L K\_4 \sqrt{\rho\_L [P\_{G1}(t) + g\rho\_L N(t) - P\_{OL2}]} \tag{7}$$

where: *zL* is the liquid valve opening (0 − 100%); *K*<sup>4</sup> is the valve constant and a tuning parameter; *POL*<sup>2</sup> is the downstream pressure after the *Va*<sup>1</sup> valve, (*N*/*m*2).

The gas mass flowrate that leaves the separator is represented by the *Va*<sup>2</sup> valve equation given by

$$m\_{GS,out}(t) = z\_G K\_5 \sqrt{\rho\_G(t)[P\_{G1}(t) - P\_{G2}]} \tag{8}$$

where: *zG* is the gas valve position (0− 100%); *K*<sup>5</sup> is the valve constant and a tuning parameter; *ρG*(*t*) is the gas density, (*kg*/*m*3); *PG*<sup>2</sup> is the downstream pressure after the *Va*<sup>2</sup> valve, (*N*/*m*2).

The boundary condition at the inlet (inflow *mL*,*in* and *mG*,*in*) can either be constant or dependent on the pressure. In this work they are constant and have been considered disturbances of the process. The most critical section of the model is the phase distribution and phase velocities of the fluids in the pipeline-riser system. The gas velocity is based on an assumption of purely frictional pressure drop over the low point and the liquid distribution is based on an entrainment model. Finally, the internal, geometric and transport equations for the pipeline system are found in [15, 17, 20].

#### **2.4. Displacement of the gas flow**

The displacement of gas in the pipeline system occurs through a relationship between the gas mass flow and the variation of the pressure inside the pipeline. The acceleration has been neglected for the gas phase, so that it is the difference of the pressure that makes the fluids outflow pipeline above. Its equation is given by

$$
\Delta P(t) = P\_1(t) - \left[P\_2(t) + \mathbf{g}\rho\_L \mathbf{a}\_L(t)H\_2\right],
$$

where: *αL*(*t*) is the average liquid fraction in riser.

It is considered that there are two situations in the riser: (i) *h*1(*t*) *> H*1, in this case the liquid is blocking the low point and the internal gas mass flowrate *mGint*(*t*) is zero; (ii) *h*1(*t*) *< H*1, in this case the liquid is not blocking the low point, so the gas will flow from *VG*<sup>1</sup> to *VG*2(*t*) with a internal gas mass flowrate *mGint*(*t*) �= 0, where *VG*<sup>1</sup> is the gas volume in upstream feed section of the pipeline, (*m*3) and *VG*<sup>2</sup> is the gas volume at the top of the riser, (*m*3).

From physical insight, the two most important parameters determining the gas flowrate are the pressure drop over the low point and the free area given by the relative liquid level

$$\xi(t) = (H\_1 - h\_1(t)) / H\_1$$

at the low point. This suggests that the gas transport could be described by a valve equation, where the pressure drop is driving the gas through a valve with opening *ξ*(*t*). Based on trial and error, the following valve equation has been proposed

$$m\_{G1}(t) = K\_2 f(h\_1(t)) \sqrt{\rho\_{G1}(t)[P\_1(t) - P\_2(t) - g\rho\_L a\_L(t)H\_2]} \tag{9}$$

where: *K*<sup>2</sup> is the valve constant and a tuning parameter; *f*(*h*1(*t*)) = *A*ˆ(*t*)*ξ*(*t*) e *A*ˆ(*t*) is the cross-section area at the low point, (*m*2); *h*1(*t*) is the liquid level upstream in the decline, (*m*); *H*<sup>1</sup> is the critical liquid level, (*m*); *ρG*1(*t*) is the gas density in the volume 1, (*kg*/*m*3). The internal gas mass flowrate from the volume *VG*<sup>1</sup> to volume *VG*2(*t*) is given by

$$m\_{Gint}(t) = \upsilon\_{G1}(t)\rho\_{G1}(t)\hat{A}(t)\tag{10}$$

where: *υG*1(*t*) is the gas velocity at the low point, *m*/*s*. Therefore, substituting equation (10) into equation (9), it has been found that the gas velocity is

$$w\_{G1}(t) = \begin{cases} K\_2 \tilde{\varsigma}(t) \sqrt{\frac{\overline{P\_1(t) - P\_2(t) - \underline{g}\rho\_I a\_L(t)}{\rho\_{G1}(t)}}} \,\forall h\_1(t) < H\_{1, \prime} \\ 0 & \forall h\_1(t) \ge H\_{1}. \end{cases} \tag{11}$$

#### **2.5. Entrainment equation**

The distribution of liquid occurs through an entrainment equation. It is considered that the gas pushes the liquid riser upward, then the volume fraction of liquid (*αLT*(*t*)) that is leaving through the *Z* valve at the top of the riser is modelled.

The volume fraction of liquid will lie between two extremes: (i) when the liquid blocks the gas flow (*υG*<sup>1</sup> = 0), there is no gas flowing through the riser and *αLT*(*t*) = *α*<sup>∗</sup> *LT*(*t*), in most cases there will be only gas leaving the riser, so *α*∗ *LT*(*t*) = 0, however, eventually the entering liquid may cause the liquid to fill up the riser and *α*∗ *LT*(*t*) will exceed zero; (ii) when the gas velocity is very high there will be no slip between the phases, so *αLT*(*t*) = *αL*(*t*), where *αL*(*t*) is the average liquid fraction in the riser.

The transition between these two extremes should be smooth and occurs as follows: when the liquid blocks the low point of the riser, the liquid fraction on top is *α*∗ *LT*(*t*) = 0, so the amount of liquid in the riser goes on increases until *α*∗ *LT*(*t*) *>* 0. At this moment the gas pressure and the gas velocity in the feed upstream section of the pipeline is very high, then the entrainment occurs. This transition depends on a parameter *q*(*t*). The entrainment equation is given by

$$\mathfrak{a}\_{LT}(t) = \mathfrak{a}\_{LT}^\*(t) + \frac{q^n(t)}{1 + q^n(t)} (\mathfrak{a}\_L(t) - \mathfrak{a}\_{LT}^\*(t)) \tag{12}$$

where

6 Will-be-set-by-IN-TECH

The displacement of gas in the pipeline system occurs through a relationship between the gas mass flow and the variation of the pressure inside the pipeline. The acceleration has been neglected for the gas phase, so that it is the difference of the pressure that makes the fluids

Δ*P*(*t*) = *P*1(*t*) − [*P*2(*t*) + *gρLαL*(*t*)*H*2]

It is considered that there are two situations in the riser: (i) *h*1(*t*) *> H*1, in this case the liquid is blocking the low point and the internal gas mass flowrate *mGint*(*t*) is zero; (ii) *h*1(*t*) *< H*1, in this case the liquid is not blocking the low point, so the gas will flow from *VG*<sup>1</sup> to *VG*2(*t*) with a internal gas mass flowrate *mGint*(*t*) �= 0, where *VG*<sup>1</sup> is the gas volume in upstream feed

From physical insight, the two most important parameters determining the gas flowrate are the pressure drop over the low point and the free area given by the relative liquid level

*ξ*(*t*)=(*H*<sup>1</sup> − *h*1(*t*))/*H*<sup>1</sup>

at the low point. This suggests that the gas transport could be described by a valve equation, where the pressure drop is driving the gas through a valve with opening *ξ*(*t*). Based on trial

where: *K*<sup>2</sup> is the valve constant and a tuning parameter; *f*(*h*1(*t*)) = *A*ˆ(*t*)*ξ*(*t*) e *A*ˆ(*t*) is the cross-section area at the low point, (*m*2); *h*1(*t*) is the liquid level upstream in the decline, (*m*); *H*<sup>1</sup> is the critical liquid level, (*m*); *ρG*1(*t*) is the gas density in the volume 1, (*kg*/*m*3). The

where: *υG*1(*t*) is the gas velocity at the low point, *m*/*s*. Therefore, substituting equation (10)

The distribution of liquid occurs through an entrainment equation. It is considered that the gas pushes the liquid riser upward, then the volume fraction of liquid (*αLT*(*t*)) that is leaving

*<sup>P</sup>*1(*t*)−*P*2(*t*)−*gρLαL*(*t*)*H*<sup>2</sup>

*ρG*1(*t*)[*P*1(*t*) − *P*2(*t*) − *gρLαL*(*t*)*H*2] (9)

*mGint*(*t*) = *υG*1(*t*)*ρG*1(*t*)*A*ˆ(*t*) (10)

<sup>0</sup> <sup>∀</sup>*h*1(*t*) <sup>≥</sup> *<sup>H</sup>*1. (11)

*<sup>ρ</sup>G*1(*t*) <sup>∀</sup>*h*1(*t*) *<sup>&</sup>lt; <sup>H</sup>*1,

section of the pipeline, (*m*3) and *VG*<sup>2</sup> is the gas volume at the top of the riser, (*m*3).

internal gas mass flowrate from the volume *VG*<sup>1</sup> to volume *VG*2(*t*) is given by

**2.4. Displacement of the gas flow**

outflow pipeline above. Its equation is given by

where: *αL*(*t*) is the average liquid fraction in riser.

and error, the following valve equation has been proposed

*mG*1(*t*) = *K*<sup>2</sup> *f*(*h*1(*t*))

into equation (9), it has been found that the gas velocity is

 *K*2*ξ*(*t*)

*υG*1(*t*) =

through the *Z* valve at the top of the riser is modelled.

**2.5. Entrainment equation**

$$q(t) = \frac{K\_3 \rho\_{G1}(t) v\_{G1}^2(t)}{\rho\_L - \rho\_{G1}(t)}$$

and *K*<sup>3</sup> and *n* are tuning parameters of the model. The details of the modelling of the equation (12) are found in Storkaas [15].

#### **3. Simulation and analysis results of the Sausen's model**

In this section are presented the simulation results of the Sausen's model for a pipeline-separator system. Initially the tuning parameters are calculated: *K*<sup>1</sup> in *Z* valve equation (6), *K*<sup>2</sup> in gas velocity equation (11), *K*<sup>3</sup> and *n* in entrainment equation (12), *K*<sup>4</sup> in *Va*<sup>1</sup> liquid valve equation (7), and *K*<sup>5</sup> in *Va*<sup>2</sup> gas valve equation (8). The calculation of these tuning parameters depends on the available data from a real system or an experimental loop, but a complete set of data is not found in the literature and is not provided by oil industries.

Therefore, to calculate the tuning parameters of the dynamic model are used the case study data carried out by Storkaas [15] through the multiphase commercial simulator OLGA [2] that accurately represents the pipeline system under slug flow [15] and the data of separator dimensioned from a tank of literature [10]. In this case study the transition of the steady flow to a slug flow occurs in the valve opening *z* = 13% (i.e., *zcrit* = 13%). Table 1 presents the data for the simulation of the dynamic model and Table 2 presents the values of the tuning parameters of the dynamic model.

Now are presented the simulation results considering the Z valve opening *z* = 12%. Figure 3 shows the varying pressures *P*1(*t*) in the upstream feed section and *P*2(*t*) at the top of the riser. Figure 4 shows the dynamics of the liquid mass flowrate (up-left) and the dynamics of the gas mass flowrate (down-left) that are entering the separator, and the dynamics of the liquid mass flowrate (up-right) and the dynamics of the gas mass flowrate (down-right) that are leaving the separator. Figure 5 shows the dynamics of the liquid level (left) and of the gas pressure (right) in the separator. It is possible to observe in all these simulation results that the varying pressures induce oscillations, but because the valve position is less than *zcrit*, these oscillations eventually die out characterizing the steady flow in pipeline-separator system.

#### 8 Will-be-set-by-IN-TECH 110 New Technologies in the Oil and Gas Industry


**Table 1.** Data for simulation dynamic model.


**Table 2.** Model tuning parameters.

**Figure 3.** Varying pressures in pipeline system with *z* = 12% (steady flow).

In the following section we are presenting the simulation results considering the Z valve opening *z* = 50%. Figure 6 shows the varying pressures throughout the pipeline system. Figure 7 presents the dynamics of the liquid mass flowrate (up-left) and the dynamics of the gas mass flowrate (down-left) that are entering the separator with peak mass flowrate of the 14 *kg*/*s* for the liquid and 2 *kg*/*s* for the gas, and the dynamics of the liquid mass flowrate (up-right) and the dynamics of the gas mass flowrate (down-right) that are leaving the separator. Figure 8 shows the dynamics of the liquid level (left) and of the gas pressure (right) in the separator. Finally, it has been observed that in all these simulation results the varying pressures induce periodical oscillations, characterizing the slug flow that happens in the pipeline-separator system. It has also been shown that the slug flow happens in intervals of 12 minutes.

8 Will-be-set-by-IN-TECH

Symbol/Value Description SI *mL*,*in* = 8.64 Liquid mass flowrate into system *kg*/*s mG*,*in* = 0.362 Gas mass flowrate into system *kg*/*s <sup>P</sup>*1(*t*) = 71.7 <sup>×</sup> 105 Gas pressure in the upstream feed section of the pipeline *<sup>N</sup>*/*m*<sup>2</sup> *<sup>P</sup>*2(*t*) = 53.5 <sup>×</sup> 105 Gas pressure at the top of the riser *<sup>N</sup>*/*m*<sup>2</sup> *r* = 0, 06 Pipeline ray *m H*<sup>2</sup> = 300 Height of riser *m L*<sup>1</sup> = 4300 Length of horizontal pipeline *m L*<sup>3</sup> = 100 Length of horizontal top section *m H*<sup>4</sup> = 4.5 Length of separator *m Ds* = 1.5 Diameter of separator *m Nt* = 0.75 Liquid level *m PG*<sup>1</sup> <sup>=</sup> <sup>50</sup> <sup>×</sup> 105 Pressure after Z valve at the top of the riser *<sup>N</sup>*/*m*<sup>2</sup> *POL*<sup>2</sup> <sup>=</sup> <sup>49</sup> <sup>×</sup> 105 Pressure after *Va*<sup>1</sup> liquid valve of separator *<sup>N</sup>*/*m*<sup>2</sup> *PGL*<sup>2</sup> <sup>=</sup> <sup>49</sup> <sup>×</sup> 105 Pressure after *Va*<sup>2</sup> gas valve of separator *<sup>N</sup>*/*m*<sup>2</sup>

> *ϕ K*<sup>1</sup> *K*<sup>2</sup> *K*<sup>3</sup> *K*<sup>4</sup> *K*<sup>5</sup> 2.55 0.005 0.8619 1.2039 0.002 0.0003

> > P1 [Bar] P2 [Bar]

<sup>0</sup> <sup>50</sup> <sup>100</sup> <sup>150</sup> <sup>200</sup> <sup>250</sup> <sup>300</sup> <sup>350</sup> <sup>50</sup>

time [min]

In the following section we are presenting the simulation results considering the Z valve opening *z* = 50%. Figure 6 shows the varying pressures throughout the pipeline system. Figure 7 presents the dynamics of the liquid mass flowrate (up-left) and the dynamics of the gas mass flowrate (down-left) that are entering the separator with peak mass flowrate of the 14 *kg*/*s* for the liquid and 2 *kg*/*s* for the gas, and the dynamics of the liquid mass flowrate (up-right) and the dynamics of the gas mass flowrate (down-right) that are leaving the separator. Figure 8 shows the dynamics of the liquid level (left) and of the gas pressure (right) in the separator. Finally, it has been observed that in all these simulation results the varying pressures induce periodical oscillations, characterizing the slug flow that happens in

**Table 1.** Data for simulation dynamic model.

55

**Figure 3.** Varying pressures in pipeline system with *z* = 12% (steady flow).

60

pressure [Bar]

65

70

75

**Table 2.** Model tuning parameters.

**Figure 4.** Input liquid (up-left) and input gas (down-left) mass flowrate in the separator and output liquid (up-right) and output gas (down-right) mass flowrate in the separator with *z* = 12% (steady flow).

**Figure 5.** Liquid level (left) and gas pressure (right) in the separator with *z* = 12% (steady flow).

**Figure 6.** Varying pressures in the pipeline system with *z* = 50% (slug flow).

**Figure 7.** Input liquid (up-left) and gas (down-left) mass flowrate in the separator and output liquid (up-right) and gas (below-right) mass flowrate in the separator with *z* = 50% (slug flow).

**Figure 8.** Liquid level (left) and gas pressure (right) in the separator with *z* = 50% (slug flow).
