Sand transport and deposition behaviour in subsea pipelines for flow assurance

: Sand transport through tubing and pipeline could cause a series of problems to ﬂow assurance, if not properly managed or controlled. The most serious problem is the accumulation and erosion in multiphase ﬂow pipelines and the surface equipment. Therefore, the importance of understanding the transport and deposition behaviour of sands through multiphase ﬂow pipelines cannot be overemphasized. This study presents the sand transport and deposition characteristics in the complicated multiphase ﬂow pipeline. The numerical result shows that the slurry velocity presents a uniform distribution in the multiphase ﬂow pipeline at the sand concentration of 5% and the sand diameter of 50 µ m. However, the slurry velocity at the bottom of the pipeline is signiﬁcantly smaller than that at the top when the sand concentration and diameter reach 30% and 300 µ m, respectively. It indicates that the sand deposition at the bottom of the pipe declines the slurry velocity and transport capacity. The deposition thickness is approximately 10% of the pipe diameter even at the low concentration of 5% sand with a small sand diameter of 50 µ m and a high slurry velocity of 1.8 m / s. The sand deposition reaches about 30% of the pipe diameter at the same low concentration and high slurry velocity when the sand diameter increases to 300 µ m.


Introduction
The demand for energy supply is a key issue for the global economy. The petroleum supply plays an important role in primary energy. In the petroleum supply networks, sand transport is generally a severe problem for most oil and gas reservoirs all over the world [1][2][3]. It is more important especially for the subsea production systems considering the flow assurance. The deposition of solid sand in the pipe bottom and equipment can cause a series of risks. For instance, the formation of the sand bed in the pipeline decreases the flow area and correspondingly increases the frictional pressure losses. Moreover, sand transport can cause erosion and corrosion problems as a result of the impact on the pipe wall and the microbial attack under the sand bed. In addition, the slurry flow containing sands could lead to equipment failure. Therefore, the deposition and transport of the sand in slurry flows attract much attention, and it is also a very significant issue in the petroleum industry.
A series of computer-controlled experiments were conducted in a thin fluid (like water) and the key parameter distribution was obtained at various flow rates, injected sand concentrations, and fluid viscosities [4][5][6]. Doron et al. [7] experimentally studied the transport behaviour of the coarse solids in a horizontal pipeline. Gillies and Shook [8] employed the traversing gamma-ray gauges to measure the sand concentration in a horizontal pipeline with coarse solids nominally from 180 to 2400 µm.

of 12
Their results indicated that the maximal concentration gradient (0.47 in their studies) was irrelevant to the mean velocity. Gillies et al. [9] experimentally studied the sand transport characteristics in a horizontal pipeline with viscous Newtonian fluids. Their results showed that a pressure gradient of 2 kPa/m was required to transport sand in laminar flow. Kaushal et al. [10][11][12][13] conducted a series of experiments to study the sand transport characteristic in multiphase flow pipelines. The measurement of the pressure drop and concentration was experimentally observed in the liquid-solid slurry flow. Duan et al. [14] investigated experimentally the critical resuspension velocity and the critical deposition velocity for different size sands in slurry flows. It indicated that the critical resuspension velocity was two to three times smaller than the critical deposition velocity. Najmi et al. [15] experimentally studied the critical flow rate of gas and liquid to keep particles moving in a horizontal pipeline in stratified and intermittent flow with a small liquid hold-up of about 0.02. The volume fraction of the solid phase (0.01-0.1%) was experimentally measured in a horizontal pipeline at a very low concentration. Rice et al. [16] utilized two kinds of particles including the glass and plastic solids to understand the sand transport behaviour in the multiphase flow pipeline.
A computational model was developed by Roco and Shook [17] to predict the quasi-uniform slurry flow with the range of sand volume concentration below 40%. Their simulated results of concentration and velocity were compared with experimental data and showed a good agreement. The sand transport process was investigated in a horizontal well by means of numerical simulations [18]. Ebadian et al. [19,20] developed an algebraic slip mixture model to predict the sand transport behaviour within a horizontal pipeline. Their comparison of the pressure drop between the computational and experimental results showed a good agreement. Danielson [21] considered the slip velocity to develop a two-fluid model for the slurry flows, and a good agreement between the simulation and experiment was obtained in the sand bed height. An analytical near-wall model was developed to describe the liquid-solid two-phase flow, and the numerical result showed a good agreement with the experimental concentration [22]. A two-fluid model was proposed to predict the fully suspended slurry flow, and it was validated by experimental data [23][24][25]. In another study, they proposed an improved mathematical model for the slurry flow within a horizontal pipeline [26]. The comparison between the simulation and experiment showed that the new model was more accurate than the previous one. Hadinoto [27] adapted the two-fluid flow model to predict the turbulence modulation characteristic in a vertical pipeline at a low solid volume fraction from 0.5% to 4.0%. The numerical and experimental results showed that the drag correlation significantly affected the accuracy of the computational fluid dynamics (CFD simulation). The Euler-granular model was employed by Kaushal et al. [28] to predict the sand transport behaviour in a horizontal pipeline with fine particles at high concentrations. The concentration distributions of the numerical and experimental results were in good agreement. Soepyan et al. [29] developed a CFD model to predict the fluid velocity of the slurry flow in a horizontal pipeline at a low concentration by volume. Table 1 lists recent studies on sand transport and deposition in multiphase flow pipelines with a focus on low concentration conditions. These studies demonstrated the importance of the effect of sands on the flow process even in low concentrations for multiphase flow pipelines.  [31] 20-300 0-1% 50 Najmi et al. 2016 [32] 150-300 0.01-0.1% 50-100 Dabirian et al. 2018 [33] 150-600 0-1% 50 Dabirian et al. 2018 [34] 45-600 0.025-1% 97 Tebowei et al. 2018 [35] 255 0.04% 100 Fajemidupe et al. 2019 [36] 212-800 0-0.006% 50.4 Leporini et al. 2019 [37] 100-1100 0.0065-0.056% 63 Leporini et al. 2019 [38] 45-600 0.00161-0.0538 50-100 The purpose of this study is to evaluate the transportation and deposition process of the sand in multiphase flows, which will contribute to the flow assurance of the subsea pipeline system. The novelty of this work includes understanding the sand transportation and deposition behaviour in multiphase flows under extreme conditions in a complicated pipeline system. This was done in this study with the aim to introduce a numerical simulation to predict the sand transport process with the sand concentration up to 30%. The effect of the particle size and velocity on the sand deposition were carried out with the Eulerian-Eulerian two-fluid model. The deposition and transport of sands are described in detail in different cross-sections of the multiphase flow pipeline, which contributes to the understanding of the flow behaviour of this liquid-solid two-phase flow system for oil and gas industry. The key contribution of this work can be summarized as (a) Development and application of a computational fluid dynamics modelling that captures the sand transport process in complex multiphase flow pipelines, (b) the baseline simulation result of liquid-sand two-phase flows have been validated against experimental data, (c) case studies for the effect of the particle size and the inlet velocity has been carried out with the detailed analysis of flow behaviour.

Governing Equations
The Eulerian two-fluid model was employed to predict the sand transport behaviour in subsea multiphase flow pipelines. In this simulation, the continuous phase was assumed as a Newtonian fluid and the water was utilized as the working fluid. Correspondingly, the mass and the momentum conservation equations were solved individually as follows.
The conservation of mass equation for phase q is [39]: where α is the phase volume fraction, ρ q is the physical density of phase q, u q is the velocity of phase q, . m sq and . m qs represents the mass transfer between the phase s and phase q. t is time. The conservation of momentum equation for phase q is described as follows [39]: ∂ ∂t (α q ρ q u q ) + ∇ · (α q ρ q u q u q ) = −α q ∇p + ∇ · τ q + α q ρ q g+ + n s=1 F sq + . m sq u sq − . m qs u qs where τ is the phase stress-strain tensor for phase q: where g is the acceleration of gravity, µ q and λ q are the shear and bulk viscosity of phase q, I is unit tensor, F sq is an interaction force between phases, F q is an external body force, F lift,q is a lift force, F wl,q is a wall lubrication force, F vm,q is a virtual mass force, and F td,q is a turbulent dispersion force, p is the pressure shared by all phases. u qs and u sq are the interphase velocity. The liquid-sand two-phase flow presents high turbulence because of the large velocity in multiphase flow pipelines. In this paper, the standard k-ε turbulence model was employed to solve the liquid-solid flow as a result of the robustness, economy, and reasonable accuracy. The turbulence kinetic energy, k, and its rate of dissipation, ε, are obtained from the following transport equations [39]: where G k is the generation of turbulence kinetic energy due to the mean velocity gradients, G b represents the generation of turbulence kinetic energy due to buoyancy, Y M represents the contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate. where C u , C ε1 , C ε2 and C ε3 are constants. σ k and σ ε are the turbulent Prandtl numbers for k and ε, respectively. µ t is the turbulent viscosity.

Numerical Schemes
The ANSYS FLUENT [39] software was used as a computational platform for the numerical study. The Eulerian multiphase flow model was employed to calculate the liquid-solid flow. The water (liquid) is the primary phase while the sand (solid) is the secondary phase in this study. The pressure-velocity coupling we used was the phase-coupled SIMPLE algorithm. The continuity equation, momentum equations, and turbulence equations were discretized with the second-order upwind scheme [40,41].
The numerical simulation was conducted in a horizontal pipeline. The pipe diameter and length were 200 mm and 8000 mm, respectively. The structured mesh was employed for the numerical study. The pipe geometry and grid system were shown in Figure 1. The fine mesh near the wall region was generated for exactly capturing the sand deposition and transport information. The velocity and outflow conditions were assigned to the pipe inlet and outlet, respectively. During the simulation, the following error formula is used to estimate the rightness of calculations and has been added in the revised manuscript.
where φ is a general variable, P is the computational cell, R φ is the computational residual, a P is the centre coefficient, a nb is the influence coefficient for the neighbouring cell, and b is the contribution of the constant part of the source term.
The resolution of the grid is a significant factor to influence the numerical result during the computational study of the sand transport and deposition in multiphase flow pipelines. The coarse mesh of 121,080 cells, the medium mesh of 203,490 cells and the fine mesh of 466,820 cells were employed here to test the sensitivity of the grid resolution on the numerical simulation. The slurry velocity at the cross-section of 20 times the pipe diameter was selected as a rational parameter to evaluate the effect of the grid resolution on the flow behaviour. The comparison results were shown in Figure 2 and Table 2, respectively. From the contours of the slurry velocity in Figure 2, it can be observed that these three different mesh cells predict almost the same distribution of the slurry velocity at the cross-section of 20 times the pipe diameter. Table 2 lists the computed slurry velocity at the central point of the cross-section of 20 times the pipe diameter under different gird cells. We can see that the relative errors for coarse mesh and medium mesh are 0.40% and 2.24%, respectively. Therefore, the medium mesh of 203,490 cells was used for the numerical simulation of sand transport and deposition in multiphase flow pipelines considering the computational cost and accuracy.

Model Validation
The developed numerical method for sand transport in multiphase flow pipelines was validated with experimental data [42]. Figure 3 shows the variation of the hydraulic gradient (Im) along with the Reynolds number (Re) in the sand-slurry flow. The sand concentration is given by the volume ratio. The Reynolds number, Re, was calculated as Re = vDρ/μ, where D was the diameter of the pipe, v was the mean velocity of the slurry flow, ρ and μ were the density and the viscosity of the slurry flow, respectively. The hydraulic gradient, Im, was defined as Im = Δp/(ρgl), where p was the pressure, l was the pipe length, and g is the acceleration of gravity. It could be seen that our simulation agreed well with experimental data when the Reynolds number is from (0.5-1.3) × 10 6 . It demonstrated that the CFD method could accurately predict the sand transport in multiphase flow pipelines.

Model Validation
The developed numerical method for sand transport in multiphase flow pipelines was validated with experimental data [42]. Figure 3 shows the variation of the hydraulic gradient (Im) along with the Reynolds number (Re) in the sand-slurry flow. The sand concentration is given by the volume ratio. The Reynolds number, Re, was calculated as Re = vDρ/μ, where D was the diameter of the pipe, v was the mean velocity of the slurry flow, ρ and μ were the density and the viscosity of the slurry flow, respectively. The hydraulic gradient, Im, was defined as Im = Δp/(ρgl), where p was the pressure, l was the pipe length, and g is the acceleration of gravity. It could be seen that our simulation agreed well with experimental data when the Reynolds number is from (0.5-1.3) × 10 6 . It demonstrated that the CFD method could accurately predict the sand transport in multiphase flow pipelines.

Model Validation
The developed numerical method for sand transport in multiphase flow pipelines was validated with experimental data [42]. Figure 3 shows the variation of the hydraulic gradient (I m ) along with the Reynolds number (Re) in the sand-slurry flow. The sand concentration is given by the volume ratio. The Reynolds number, Re, was calculated as Re = vDρ/µ, where D was the diameter of the pipe, v was the mean velocity of the slurry flow, ρ and µ were the density and the viscosity of the slurry flow, respectively. The hydraulic gradient, I m , was defined as I m = ∆p/(ρgl), where p was the pressure, l was the pipe length, and g is the acceleration of gravity. It could be seen that our simulation agreed well with experimental data when the Reynolds number is from (0.5-1.3) × 10 6 . It demonstrated that the CFD method could accurately predict the sand transport in multiphase flow pipelines.

Slurry Velocity
The sand content in multiphase flow pipelines is usually very small. However, the sand content of up to 30% is employed in this simulation for accelerated transport and deposition considering some extreme circumstances.  Table 3.

Slurry Velocity
The sand content in multiphase flow pipelines is usually very small. However, the sand content of up to 30% is employed in this simulation for accelerated transport and deposition considering some extreme circumstances.  Table 3. Table 3. Detailed inputs for the numerical simulation of sand transport and deposition.

Slurry Velocity
The sand content in multiphase flow pipelines is usually very small. However, the sand content of up to 30% is employed in this simulation for accelerated transport and deposition considering some extreme circumstances.  Table 3.  Under the condition of the lower concentration of the fine sand (d = 50 µ m, c = 5%), we could see from Figure 4 that the slurry velocities generally presented a uniform distribution in the turbulence area. Correspondingly, the vertical slurry velocity was symmetrical about the central axis. However, it could be seen from Figure 5 that the slurry velocities showed an extraordinarily nonuniform distribution in the higher concentration of the coarse sand (d = 300 µ m, c = 30%). In this higher concentration condition, the slurry velocities near the pipe top area were obviously larger than those close to the bottom area of the pipe. One of the major reasons was that the sands deposited in the pipe bottom in the slurry flow because of the large density difference between the liquid and the sand. In addition, the slurry velocities near the pipe wall area declined sharply due to the strong viscous shear stress in the turbulent boundary layer and non-slip boundary condition. This indicated that the liquids consumed more energy to carry these solid sands. Figure 6 shows the sand concentration at different cross-sections away from the inlet of the pipeline. The initial computational condition included the sand concentration of 5%, the sand diameter of 50 µ m and the slurry velocity of 1.8 m/s, respectively. It can be seen from these concentration contours that the sand can hardly settle in the pipe bottom within the distance of 10 times the pipe diameter from the pipe inlet (Figure 6a). After that, the sand begins to deposit as a result of the effect of the gravitational force along with the flow direction. In these simulation cases, the sand deposition was observed at the cross-section of 20 times the pipe diameter (Figure 6b). The deposition of the solid sand can be obviously found when the distance from the pipe inlet reaches 30 times the pipe diameter (Figure 6c). In other words, the thickness of the sand deposition is very small and it accordingly could be ignored, when the sand travelled less than 10 times the pipe diameter. Nevertheless, the deposition thickness reached approximately 10% of the pipe diameter, when the sand travelled about 30 times the pipe diameter in this simulation case. Under the condition of the lower concentration of the fine sand (d = 50 µm, c = 5%), we could see from Figure 4 that the slurry velocities generally presented a uniform distribution in the turbulence area. Correspondingly, the vertical slurry velocity was symmetrical about the central axis. However, it could be seen from Figure 5 that the slurry velocities showed an extraordinarily nonuniform distribution in the higher concentration of the coarse sand (d = 300 µm, c = 30%). In this higher concentration condition, the slurry velocities near the pipe top area were obviously larger than those close to the bottom area of the pipe. One of the major reasons was that the sands deposited in the pipe bottom in the slurry flow because of the large density difference between the liquid and the sand. In addition, the slurry velocities near the pipe wall area declined sharply due to the strong viscous shear stress in the turbulent boundary layer and non-slip boundary condition. This indicated that the liquids consumed more energy to carry these solid sands. Figure 6 shows the sand concentration at different cross-sections away from the inlet of the pipeline. The initial computational condition included the sand concentration of 5%, the sand diameter of 50 µm and the slurry velocity of 1.8 m/s, respectively. It can be seen from these concentration contours that the sand can hardly settle in the pipe bottom within the distance of 10 times the pipe diameter from the pipe inlet (Figure 6a). After that, the sand begins to deposit as a result of the effect of the gravitational force along with the flow direction. In these simulation cases, the sand deposition was observed at the cross-section of 20 times the pipe diameter (Figure 6b). The deposition of the solid sand can be obviously found when the distance from the pipe inlet reaches 30 times the pipe diameter (Figure 6c). In other words, the thickness of the sand deposition is very small and it accordingly could be ignored, when the sand travelled less than 10 times the pipe diameter. Nevertheless, the deposition thickness reached approximately 10% of the pipe diameter, when the sand travelled about 30 times the pipe diameter in this simulation case.

Effect of the Particle Size
The effect of the particle size on the sand deposition was numerically simulated at the inlet velocity of 1.8 m/s with a sand concentration of 5%. The detailed inputs for the numerical stet-up are listed in Table 4. The sand diameter is assumed to be 50 μm, 150 μm, and 300 μm, respectively. The results of the sand concentration are shown in Figure 7 at the cross-section of 30 times the pipe diameter. We can see that the height of the sand bed in the pipe bottom is small when the particle diameter is less than 50 μm. However, the height of the sand bed increases rapidly and it is even approximately equal to one-third of the pipe diameter when the sand diameter reaches 300 μm. Under this condition, the slurry flow cannot carry the coarse sand and correspondingly the severe deposition occurs in this low velocity of 1.8 m/s. It indicates that the sands cannot be carried by the liquid and the deposition will inhibit the flow assurance, especially for the subsea multiphase flow pipeline system.

Effect of the Particle Size
The effect of the particle size on the sand deposition was numerically simulated at the inlet velocity of 1.8 m/s with a sand concentration of 5%. The detailed inputs for the numerical stet-up are listed in Table 4. The sand diameter is assumed to be 50 µm, 150 µm, and 300 µm, respectively. The results of the sand concentration are shown in Figure 7 at the cross-section of 30 times the pipe diameter. We can see that the height of the sand bed in the pipe bottom is small when the particle diameter is less than 50 µm. However, the height of the sand bed increases rapidly and it is even approximately equal to one-third of the pipe diameter when the sand diameter reaches 300 µm. Under this condition, the slurry flow cannot carry the coarse sand and correspondingly the severe deposition occurs in this low velocity of 1.8 m/s. It indicates that the sands cannot be carried by the liquid and the deposition will inhibit the flow assurance, especially for the subsea multiphase flow pipeline system. Table 4. Detailed inputs for the numerical simulation of sand transport and deposition under the condition of different particle sizes.

Effect of the Particle Size
The effect of the particle size on the sand deposition was numerically simulated at the inlet velocity of 1.8 m/s with a sand concentration of 5%. The detailed inputs for the numerical stet-up are listed in Table 4. The sand diameter is assumed to be 50 μm, 150 μm, and 300 μm, respectively. The results of the sand concentration are shown in Figure 7 at the cross-section of 30 times the pipe diameter. We can see that the height of the sand bed in the pipe bottom is small when the particle diameter is less than 50 μm. However, the height of the sand bed increases rapidly and it is even approximately equal to one-third of the pipe diameter when the sand diameter reaches 300 μm. Under this condition, the slurry flow cannot carry the coarse sand and correspondingly the severe deposition occurs in this low velocity of 1.8 m/s. It indicates that the sands cannot be carried by the liquid and the deposition will inhibit the flow assurance, especially for the subsea multiphase flow pipeline system.  Table 5. The sand concentration is 5% and the sand diameter is 50 µm for the numerical condition. The inlet velocity is about 0.8 m/s, 1.8 m/s and 3.0 m/s for the simulation, respectively. The sand deposition is obviously observed in the pipe bottom while the inlet velocity is approximate 0.8 m/s. In this condition, the slurry velocity is not large enough to carry the solid sand, and correspondingly the deposition occurs. If the inlet velocity increases to 3.0 m/s, the sand can hardly deposit in the pipe bottom. Therefore, the slurry flow velocity is a major influence factor for sand deposition and transport.   Figure 8 describes the effect of the inlet slurry velocity on the sand deposition and transport behaviour. The detailed inputs for the numerical stet-up are listed in Table 5. The sand concentration is 5% and the sand diameter is 50 μm for the numerical condition. The inlet velocity is about 0.8 m/s, 1.8 m/s and 3.0 m/s for the simulation, respectively. The sand deposition is obviously observed in the pipe bottom while the inlet velocity is approximate 0.8 m/s. In this condition, the slurry velocity is not large enough to carry the solid sand, and correspondingly the deposition occurs. If the inlet velocity increases to 3.0 m/s, the sand can hardly deposit in the pipe bottom. Therefore, the slurry flow velocity is a major influence factor for sand deposition and transport.

Discussion
Al-lababidi et al. [43] reported the experimental result of the sand transportation and deposition behaviour in a pipe of 0.05 m inner diameter under different sand concentrations and inlet velocities. If the inlet velocity was fixed at 1 m/s, most sands transported in the liquid core with the sand concentration of 0.00161%, while sand streaks were observed on the pipe bottom with the sand concentration of 0.0538%. This indicated that higher sand concentration resulted in the possibility of sand deposition in multiphase flows. Our numerical results in Figures 4 and 5 are in accord with these experimental observations. Furthermore, their experiments showed that increasing inlet velocity decreased the deposition of the sand in multiphase flows. For instance, the sand formed dunes on the bottom of the pipe and few sand particles moved on the top of the dunes with the sand concentration of 0.0108% when the inlet velocity was about 0.1 m/s. If the inlet velocity increased to

Discussion
Al-lababidi et al. [43] reported the experimental result of the sand transportation and deposition behaviour in a pipe of 0.05 m inner diameter under different sand concentrations and inlet velocities. If the inlet velocity was fixed at 1 m/s, most sands transported in the liquid core with the sand concentration of 0.00161%, while sand streaks were observed on the pipe bottom with the sand concentration of 0.0538%. This indicated that higher sand concentration resulted in the possibility of sand deposition in multiphase flows. Our numerical results in Figures 4 and 5 are in accord with these experimental observations. Furthermore, their experiments showed that increasing inlet velocity decreased the deposition of the sand in multiphase flows. For instance, the sand formed dunes on the bottom of the pipe and few sand particles moved on the top of the dunes with the sand concentration of 0.0108% when the inlet velocity was about 0.1 m/s. If the inlet velocity increased to 1 m/s, most sands transported in the liquid core at the same condition of the sand concentration. These experimental results are demonstrated in our simulations in Figure 8. Ofei and Ismail [44] numerically investigated the sand transport and deposition in multiphase flows in a pipe of 103 mm diameter. Their numerical results showed that the maximum sand concentrations deposited on the pipe bottom were 15.1%, 21.5%, 25.8%, and 30.7% for the particle size of dp = 90 µm, 150 µm, 210 µm, and 270 µm, respectively. This illustrated that larger sand particles could accumulate on the pipe bottom and lead to the safety issue of the pipe operation. Our results in Figure 7 agree well with this trend of the particle size effect in point view of the qualitative analysis.

Conclusions
The Eulerian two-fluid model was developed to predict the transport and deposition behaviour of sands in multiphase flow pipelines. The numerical result shows that the slurry velocities at the vertical cross-section presented an extraordinarily non-uniform distribution in a higher concentration of the coarse sand compared to the condition of lower concentration of the fine sand. The sand could hardly settle in the pipe bottom at the cross-section of 10 times the pipe diameter, while the deposition was obviously observed at the cross-section of 30 times the pipe diameter. If the sand diameter was 300 µm, the height of the sand bed would increase rapidly and it even approximately equal to one-third of the pipe diameter. For the sand diameter of 50 µm with a concentration of 5% slurry flow, the deposition was obviously observed in the pipe bottom with the inlet velocity of 0.8 m/s. The sand could hardly deposit in the pipe bottom if the inlet velocity increased to 3.0 m/s.
For the future work, we are carrying out the numerical and experimental studies on gas-liquid-solid multiphase flows in inclined pipelines, which will focus on the sand transport and deposition behaviour in this complicated pipeline system. The fluid flow analysis will be integrated to study the particle erosion issue, which is expected to contribute to the flow assurance of the complex oil and gas pipeline system.