Regular Article
A model of wax deposition under oilgas twophase stratified flow in horizontal pipe
^{1}
Department of Fuel, Army Logistics University of PLA, Chongqing 401331, China
^{2}
Training Base, Army Logistics University of PLA, Xiangyang 441118, China
^{3}
National Engineering Laboratory for Pipeline Safety, China University of Petroleum (Beijing), Beijing 102249, China
^{*} Corresponding authors: duanjimiao@126.com, hsliu0820@163.com
Received:
25
May
2018
Accepted:
10
October
2018
A model of wax deposition based on molecular diffusion mechanism, for oilgas twophase stratified pipe flow is developed. In the model, unidirectional fully developed flow analyses of momentum, heat and mass transfer are presented. And, a cube cage model is used to describe the wax deposit structure considering the effect of oil flow shear on the deposit. Calculation of wax deposit is compared well with a flow loop experiment. In particular, the model could give the wax deposit forming a crescent shape at the cross section of pipe, which is observed in different experiments. Furthermore, the cause of forming a crescent shape is revealed, which is indicated by the nonuniform circumferential distribution of mass flux for wax deposition along the pipe wall wetted by the oil. The mass flux from oil bulk flow to the oildeposit interface is closely related to three parameters, diffusivity at oildeposit interface, the temperature gradient at the oildeposit interface at different time, and the slope of the wax solubility curve at oildeposit interface temperature.
© J. Duan et al., published by IFP Energies nouvelles, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Nomenclature
r_{ i } : Radial distance from the center to the oildeposit interface, m
h_{e} : Overall heat transfer coefficient, W/(m^{2} K)
T_{e} : Environment temperature, (°C)
F_{w} : Wax fraction in the deposit, wt %
ΔL: Pipe length of wax deposit buildup, m
D_{WO} : Diffusivity wax in crude oil, m^{2}/s
C : Wax concentration in crude oil, kg/m^{3}
μ : Kinematic viscosity coefficient, m^{2}/s
V_{A} : Molar volume of wax molecules
D_{e} : Effective diffusivity in the deposit, m^{2}/s
ρ_{dep} : Density of deposit, kg/m^{3}
a : Side length of cubic form, m
V_{total} : Total volume of the cubic cage structure, m^{3}
Z : Ratio of thickness of the frame and length of the cube
V_{solid} (%): Volume fraction of solid wax in cubic cage structure
V_{liquid} (%): Volume fraction of the wrapped liquid oil
V_{β,total} : Volume of the cage structure after deformation
V_{β,solid} (%): Volume fraction of the solid wax after deformation
V_{β, liquid} (%): Volume fraction of wrapped liquid oil after deformation
V_{solid,O} (%): Solid volume fraction of the initial deposition layer
G : Shear modulus of cubic cage structure, Pa
G_{O}: Initial shear modulus, Pa
w : Velocity in the axial direction, m/s
Γ_{C}: Effective mass transfer
1 Introduction
Recently, the rapid economic development in China with the need for oil keeping up, the oil import dependency rapid rises, up to 64.4% in 2017, and exceeds the internationally recognized warning level 50% [1, 2] . It makes economic security under potential threat because of the energy factor [3, 4]. As oil and gas exploiting onshore and offshore, deep water reserving abundant oil and gas will be key production area in a short time. The low temperature (4 °C) and intense heat transfer between oil and environment in deep water make a few kilometers of submarine pipeline to be against severe flow assurance.
As the wax crude oil flows through subsea pipeline on the ocean floor at approximately 4 °C, the temperature of oil will decrease because of the heat transfer to the surroundings. Then, once the pipe wall temperature decreases below the oil Wax Appearance Temperature (WAT), the wax in the oil will precipitate and deposit on the pipe wall [5]. The deposit makes the flow area and efficiency decrease. The deposit thickness grows with time, and even makes pipe blockage in some cases [6], as shown in Figure 1.
Fig. 1 Wax deposit and blockage in submarine pipeline. 
The oil exploited from the Bohai Sea and South China Sea in China is mostly waxy [7]. The intense heat transfer between oilgas pipeline and the subsea environment in the deep sea (4 °C) aggravates the wax deposit. Compared with the wax deposition in the single phase flow, due to the complexity of flow, heat and mass transfer of the oilgas twophase pipe flow, and the wax deposition process becomes much more complex. There are some flow patterns in oilgas twophase flow, such as bubble flow, stratified flow, slug flow, annular flow, and dispersed flow for different flow condition. The stratified flow is the fundament and known as one of the most common and dominant flow pattern.
Matzain et al. [8], Apte et al. [9] and Rittirong [10] observed that the wax deposit only buildup on the wall wetted by the oil, and there was no deposit on the wall wetted by gas, and the deposit emerged a crescent shape at the pipe cross section, through the wax deposition in oilgas stratified flow experiment. A semiempirical model for wax deposit in oilgas stratified flow was developed, based on the molecule diffusion mechanism and the assumption that all precipitated wax crystals on pipeline wall can deposit. The model takes the mutual effect of molecule diffusion and shear dispersion into consideration [9–11]. The wax deposit experiment for oilgas stratified flow was also carried out by Gong et al. [12] with Daqing crude oil and air. They also observed that wax deposit presenting a crescent shape at the pipe cross section in the experiment. A wax deposition model based on classical molecule diffusion was established, considering the influence of superficial gas velocity, superficial liquid velocity, liquid hold up and flow pattern Reynolds number.
Onedimensional hydrodynamic, heat and mass transfer correlations are usually applied to predict the wax deposit thickness in multiphase flow previously. However, these models can only obtain the equivalent thickness of wax deposition assuming the deposit distributing in the whole pipe inner wall, and just give a general description for the deposit distribution around the circumference of the inner wall. Furthermore, these models can’t give the reason for deposit forming a crescent shape. In this work, a threedimensional wax deposition model of oilgas stratified pipe flow is developed, and the reason of forming crescent shape is revealed by analyzing the wax deposit model, and it is indicated by the nonuniform circumferential distribution of mass flux for wax deposition along the pipe wall wetted by the oil.
2 Wax deposit model
For oilgas twophase stratified flow, as the wax deposit increases on the pipe wall, the oildeposit interface and the gasoil interface are two moving boundaries, which must be updated at each computational time step, as shown in Figure 2. Wax deposition buildups only at bottom part of the pipe wall wetted by the oil. There is not any deposit on the top part of the pipe wall connected with the gas. An assumption that there is no mass transfer between the oil phase and the gas phase or there is no phase change due to the variation of pressure and temperature, is adopted to simplify the calculation.
Fig. 2 Illustration of wax deposit in oilgas stratified pipe flow. 
2.1 Model development
According to molecular diffusion mechanism, the wax deposit thickness () and the wax fraction of deposit (F_{w}), are generated by the wax molecular mass flux (M) from pipe center to the oildeposit interface. The mass flux M is a function of wax deposit thickness, δ and the wax fraction of deposit F_{w}.(1)
The rate of wax molecular mass flux into the deposit in the total derivative form can be written as:(2)where is the inner radius of pipe, is the radial distance from the center to the oildeposit interface, is time.
The first term on the right hand side of equation (2) is the rate of wax mass flux, leading to the deposit thickness growth (buildup term), . The second term on the right hand side of equation (2) is the rate of wax mass flux, leading to wax fraction increase (aging term), .
The rate of wax mass flux into the oildeposit interface can be written as follows:(3)where ΔL is the pipe length of wax deposit buildup, D_{wo} is the wax diffusivity in oil, can be calculated by the HaydukMinhas relationship [13], as given(4)(5)where T is the oil temperature, μ is the oil dynamic viscosity, V_{A} is the molar volume of wax molecules.
The rate of mass increase leading to aging only depends on the wax mass flux that goes into the deposit. Thus, the aging term can be written as(6)
The buildup term depends on the difference between the wax mass flux from pipe center to the oildeposit interface and the wax mass flux diffusing into the deposit, and it is given as(7)
The wax mass in the deposit can be written as(8)where is the density of deposit, assumed to be same as oil.
Equation (8) can be given as(9)(10)then(11)(12)
By rearranging equations (7) and (11), the wax deposit thickness changes with time can be given as(13)
So, the wax fraction of the deposit can be calculated from equation (6) and equation (12):(14)where D_{e} is the effective wax diffusivity in oil as given by the Cussler relationship.(15)where α is the ratio of length to diameter of wax crystal particles in the deposit.
Singh et al. [14] assume that the equivalent crystal aspect ratio, α, varies linearly with the wax fraction, as shown in equation (16):(16)
Singh et al. [14] used k_{α} as the fitting parameter to fit the simulation with the experimental results. Thus, this value is unknown prior to conducting the actual experiments.
2.2 The effect of shear stress
Bidmus and Mehrotra [15] and Mehrotra and Bhat [16] show wax fraction in deposit increases with the increase of oil flow rate through experiment. A simple cubic cage model is used to describe the structure of the wax deposition layer, by adding the effect of shear on wax fraction in deposit. In order to reflect the effect of shear on wax fraction simply and vividly, quantitative relationship between solid wax, liquid crude oil and shear in the deposit layer is given by the structure.
The structure unit of deposit layer presents a cubic cage structure under static condition, as shown in Figure 3. The precipitated solid wax is located at the edge of the cage structure, forming its frame, and the liquid oil is wrapped in the cage structure. In order to simplify the calculation, it is assumed that each wax crystal cage structure presents a cubic form, the size is a × a × a. Cross section of the frame of the solid wax crystal at the edge is a rectangle, the size is ξ × ξ. Total volume of the cubic cage structure is V_{total} = a^{3}. Ratio of thickness of the frame and length of the cube is Z = ξ/a. Volume fraction of solid wax in cubic cage structure is V_{solid} (%) = 12Z^{2} − 16Z^{3}, and volume fraction of the wrapped liquid oil is V_{liquid} (%) = 1 − 12Z^{2} + 16Z^{3}.
Fig. 3 The cubic cage structure unit of deposit layer under static condition. 
As shown in Figure 3, when the cube cage structure is applied a shear along the flow direction, the cage structure will cause a certain inclination along shear direction and the inclination angle is β. The top, bottom, left and right sides of the cage cube are square, and the front and back sides are rhombus. Frame of wax crystals remains the same, the total volume is unchanged, and volume of the cage structure after deformation is , volume fraction of the solid wax is , volume fraction of wrapped liquid oil is . When there is no liquid oil in the inclination cage structure, V_{β,liquid} (%) = 0, the maximum value of the inclination angle β is , so .
All of the above are functions of Z and inclination angle β. Value of Z is calculated by the precipitated solid portion of wax molecule of the oil spread into the deposition layer with the bound oil. In the initial stage of wax deposition, solid volume fraction of the initial deposition layer V_{solid,O} (%) is calculated by the solidliquid equilibrium of the oil system with the pipe wall temperature T_{wall}:(17)where is the solid fraction of the initial deposition layer, f(T_{wall}) is the wax precipitation fraction under the temperature T_{wall}.
It is assumed that the solid fraction increasing completely due to the wax mass transfer inside the deposit. So, without shear, the solid fraction inside the deposit, V_{solid} (%) given as(18)
As shown in Figure 3, when the shear imposes on the cubicalcage, the relationship between the shear strain , and shear stress can be given as(19)where G is the shear modulus of cubic cage structure. Assumed that the shear modulus increases linearly with the edgethickness .(20)where is the initial shear modulus.
The shear stress at the interface can been given as(21)
And then the solid fraction of the deposit under the shear can be obtained by(22)(23)
Singh et al. [14] presented the formation of incipient gel layer on the cold surface and assumed that the thickness of incipient gel layer δ_{o} = 10 μm. So the deposit thickness can be given as(24)
And then the wax fraction inside the deposit F_{w} can be calculated by(25)where is the saturation concentration of wax dissolved in oil.
It is important to note that in equation (16), k_{α} is the fitting parameter to fit the deposit thickness with the experimental results. In equation (23), the optimum initial shear modulus G_{O} is selected by matching the predicted average wax fraction of the deposit layer to the experiment measurement.
2.3 Hydrodynamic, heat and mass transfer calculation
The hydrodynamic, heat and mass transfer for oilgas twophase flow should be fundamentally studied to predict wax deposition. During the process of oilgas twophase flow in a pipe, there are the heat transfer and the momentum transfer between oil and gas, gas and the environment, oil and the environment. As the temperature of the oil decreases along the pipe, the oil phase flows with smaller velocity and a larger cross sectional area and the gas phase flows with larger velocity and a smaller cross sectional area. The axial momentum equations for oil and gas phases are used to calculate the pressure gradient, the liquid holdup, and the velocity profile of oil in the pipe cross section. The energy equations for two phases are utilized to predict the fluid temperature profile on cross section and fluid temperatures along the pipe segment. In this research, a threedimensional (3D) model is established using a unidirectional and fully developed flow assumption, in order to reduce the computational intensity, to predict wax deposition in oil/gas stratified flow in the circular pipe. The details of solving nonisothermal hydrodynamics, heat transfer phenomena can be found in our previous works [17], and are not repeated.
The temperature and shear stress on the oildeposit interface in pipe can be obtained by solving momentum and energy transfer equations. And the radial concentration distribution of the wax molecule in pipe can be obtained by solving the wax molecular mass transfer equation.
The wax concentration profile is calculated using the mass transfer equation, equation (26), in which the precipitation of wax molecule in the bulk oil is neglected and it is valid for laminar flow [13]. The mass transfer equation of wax molecules, in (x, y, z) rectangular coordinate system, is shown as equation (26):(26)where C is the wax concentration of crude oil, kg/m^{3}; Γ_{C} is the turbulent effective mass diffusion coefficient of wax molecules, Γ_{C} = ν_{m}/Sc + ν_{t}/σ_{C}, m^{2}/s; σ_{C} is the mass transfer diffusion coefficient resulted by turbulence; Sc is the Schmidt number, Sc = ν_{m}/D_{WO}; D_{WO} is the wax molecular diffusion coefficient in crude oil, m^{2}/s.
In (x, y, z) rectangular coordinate system, the oil phase and gas phase have irregular calculation area for gasliquid smooth stratified flow. But the calculation area becomes regular in bipolar coordinate system (ξ, η, z)(27)
The boundary condition of mass transfer equation is shown in equation (28). Concentration of the wax molecule in the cross section of the pipe inlet (z = 0) is assumed to be the saturated wax molecular concentration (C_{b}) at the pipe inlet temperature, C = C_{b}. The wax molecule is assumed to be solidliquid equilibrium on the oildeposit interface (r = r_{i}) and the gasliquid interface, and the wax molecule concentration on the oildeposit interface is wax saturated concentration in crude oil C = C_{ws} (T). For vertical center boundary of the pipe, due to symmetric relation of the pipe cross section at the vertical center position, ∂C/∂η = 0.(28)
2.4 Calculation procedure
The calculation process of wax deposition thickness and wax fraction in deposition layer of oilgas stratified flow is as shown in Figure 4. Firstly, velocity, temperature and wax molecular concentration profile are calculated along the pipeline. Then temperature, concentration gradient of the wax molecule and shear stress at the oildeposit interface are calculated. And growth rate of wax deposit thickness and wax fraction in the deposition layer can be obtained from the axial direction. According to the wax deposition model, the RungeKuttaFehlberg (RKF45) method is used to solve the ordinary differential equation. The wax deposition thickness and wax fraction can be obtained by integration in the set time step. Finally, the computing grid is redevised according to the requirements, and the above steps are repeated for the next time step calculation until the set simulation time reached.
Fig. 4 The calculation process of wax deposit for oilgas stratified pipe flow. 
As the deposit grows on the pipe wall, the oilgas interface and oildeposit interface are two moving boundaries that need to be updated at each time. This quasisteadystate simulator is modular in structure, and each module makes computations for a particular aspect. The simulator consists of four main modules that perform the following tasks: two phase flow hydrodynamic calculations, two phase flow heat transfer calculations, two phase flow mass transfer calculations, and wax deposition calculations.
The wax deposit rate is slow, so the quasisteadystate simulator is generally used in wax deposit calculation. The deposition rate is expected to be fast initially and then slow down, inasmuch as there is a larger driving force initially. Hence, initially very small time step is used (Δt = 30 s, for t < 10 min), and then progressively bigger time steps are used (Δt = 2 min, for 10 min < t < 1 h, and Δt = 10 min thereafter).
3 Results and discussion
In order to test the feasibility of the wax deposition model, it has been applied to the oilgas stratified flow experiment conducted by Matzain et al. [8] with a flow loop. Matzain et al. [8] find that there is no wax deposition on the wall wetted by the gas phase, and the deposition appearing the crescent on the lower part of the wall wetted by oil, as shown in Figure 5. Matzain et al. [8] adopted two nonconcentric circles with different radius to describe the crescent shape. The position of A and B is determined by the liquid holdup in the oilgas stratified smooth flow, and the thickness of wax deposition δ at the bottom of the pipe is acquired from the Vernier Caliper.
Fig. 5 Wax thickness distribution along the pipe wall for oilgas stratified flow. 
The wax concentration at the inlet of the test section has a reducing trend because of the reutilization of crude oil in the experiment flow loop with more wax molecules diffusing to the inner wall, precipitating from the oil and depositing on the wall. According to mass conservation of wax molecules, the wax concentration at the inlet of the test section varies with time and the relationship, as shown in equation (29)(29)where C_{inlet} represents the wax concentration at different time points at the inlet of the pipe, C_{ini,inlet} represents the initial wax concentration at the inlet of the pipe without wax deposition occurring, m_{dep} and V_{dep} are the mass and the volume of wax deposition, respectively.
3.1 Comparison between the modeling and the Matzain’s experiment
Matzain et al. [8] carry out the wax deposition experiment under the condition of gasliquid stratified smooth flow, using the South Pelto crude oil having 6.6% wax content and the Tulsa City natural gas as the experimental medium. The experimental facility, in general, consists of: an oil system; a gas system; a multiphase system with an oil/gas temperature trimming section, a test section, and an oil/gas separator; a chilling system; and a heating system. The detailed description of the facility is given by Matzain et al. [8].
The properties of the experimental oil and gas and the operation conditions of the flow loop are listed in Table 1. The detailed parameters of the viscosity of temperature characteristic curve and the solubilitytemperature characteristic curve of wax molecules of the South Pelto crude oil are shown in Figures 6 and 7.
Fig. 6 Viscositytemperature curve of South Pelto crude oil. 
Fig. 7 Solubilitytemperature curve of South Pelto crude oil. 
The operating conditions of the experiment and the oil properties.
The relationship of viscosity and temperature is given as(30)where T_{o} is the oil temperature, °C; μ_{o} is the crude oil viscosity, mPa s.
The relationship between wax solubility and temperature:(31)where C_{wax} (wt%) is the weight fraction of wax molecules dissolved in oil.
In Figure 8, a comparison between wax deposition thickness at cross section from experiment and model prediction is demonstrated, and the deposition thickness at the bottom of pipe obtained from the experiment and the model prediction is provided. Obviously, the model can give an accurate prediction on the thickness, the shape and the distribution of wax deposition.
Fig. 8 Comparison of distribution of wax deposition in the ring section between the experiment and prediction. 
It must be pointed out that, liquid holdup of oilgas twophase stratified flow given by the wax deposition model is 0.53, and it is consistent with the liquid holdup of 0.55 measured by the Gamma density meter in the experimental process. Accurate prediction of the hydraulic condition is the basis of accurate prediction of growth rate of the wax deposition layer. Figure 9 shows the predicted angular distribution of wax deposition layer thickness at the pipe wall. Wax deposition thickness is the largest at the lowest point of pipeline section (γ = 0). The closer to contact points A and B at the gasliquid interface with the pipe wall, the smaller the wax deposition layer thickness is, and the thickness of wax deposition is from 1.86 mm gradually to zero. As shown in Figure 10, predicted angular distribution of wax fraction in deposition layer is given, which is similar to that of the thickness of wax deposition of the pipe inner wall. At the thickest position of wax deposition, maximum of wax fraction in the wax deposition can also be obtained. The wax fraction in deposition layer decreases from 18.52% to 14.78%. The hardness of wax deposition layer is mainly determined by the wax fraction in the deposition. It is of great significance to design and implementation of the paraffin removal method of pipeline.
Fig. 9 Predicted angular distribution of deposit thickness, γ = 0° at the bottom of the pipe. 
Fig. 10 Predicted angular distribution of wax content in the deposit, γ = 0° at the bottom of the pipe. 
In Figure 11, the thickness variation with time predicted by the model is compared with the experimental results. In Figure 12, there is a contrast of the wax fraction in the deposit layer at the bottom of the pipe between the experiment and the prediction. Obviously, the deposition thickness and the wax fraction calculated by the model are very close to the experiment when time reaches 24 h. The thickness at the bottom of the pipe of experiment and prediction are 2.0 mm and 1.86 mm, respectively, with a 7.0% relative error; the wax fraction at the bottom from experiment and prediction are 17.9% and 18.5%, respectively, with a relative error about 3.4%.
Fig. 11 Comparison between experimental deposit thickness in the bottom of pipe and the prediction. 
Fig. 12 Comparison between experimental wax content in deposit in the bottom of pipe and the prediction. 
3.2 Analysis the wax deposit crescent shape
The radial mass transfer of wax molecules is a direct driving force for wax deposition. Huang et al. [18] confirmed that the growth rate of wax deposition is mainly depended on J_{wax}, the convection diffusion mass flux of wax molecules in oil, as shown in equation (32):(32)
The wax concentration gradient at the oildeposit interface is obtained by the calculation of liquidsolid equilibrium. In other words, is represented by the slope of the solubilitytemperature characteristic curve of wax molecules in oil at the temperature of oildeposit interface. Then, the growth rate of wax deposition is codetermined by the three factors: (1) the diffusion coefficient of wax molecules at the temperature of oildeposit interface, ; (2) the temperature gradient at oildeposit interface, ; 3) the slope of the solubilitytemperature characteristic curve of wax molecules in oil at the temperature of oildeposit interface, . Two of the three factors, and are the function of the oildeposit interface temperature. If the oildeposit interface temperature changes, and will change the growth rate of wax deposition. To interpret the crescent shape of wax deposition in gasliquid smooth stratified flow, , the distribution of the temperature of oildeposit interface around the wall, and these three factors: , and should be obtained clearly.
(1) The distribution of the oildeposit interface temperature
The accurate prediction of thermodynamics is the foundation for estimating the growth rate of wax deposition. The wax deposition interface temperature can be accurately determined by solving the heat transfer equation. The distribution of oildeposit interface temperature around the circumference of the inner wall wetted by the oil phase is described in Figure 13, under the condition of gasliquid smooth stratified flow. As shown in Figure 13, the temperature at the bottom of the pipe is 18.91 °C, the minimum value for the temperature distribution; for the position getting close to point B where the gasliquid interface intersects with the wall, the interface temperature increases gradually. And the temperature at B is the maximum value, 24.43 °C.
Fig. 13 Circular distribution of oildeposit temperature. 
(2) The diffusion coefficient of wax molecules at the oildeposit interface
Figure 14 gives the distribution, around the circumference, of diffusion coefficient of wax molecules at the oildeposit interface temperature, which is similar to the interface temperature distribution around the wall. Along the wall, when the position gets close to point B, the diffusion coefficient increases with the growth of interface temperature and the decrease in oil viscosity, which promotes the increase in the wax deposition thickness.
Fig. 14 Circular distribution of wax molecular diffusion coefficient. 
(3) The temperature gradient at the oildeposit interface
The growth rate of wax deposition is an increasing function of the temperature gradient at the oildeposit interface. The larger the temperature gradient, the easier the thickness of wax deposition grows. According to the discretization of the energy transfer equation, the temperature gradient at the oildeposit interface can be described as:(33)
The temperature gradient distribution around the circumference of oildeposit interface is shown in Figure 15. With being close to B where the gasliquid interface intersects with the wall, has a slow decrease at first and then an obvious reduction. The gradual decrease leads to a little decrease in the deposition thickness, while the obvious reduction results in a great decrease in the deposition thickness.
Fig. 15 Circular distribution of temperature gradient of oildeposit. 
(4) The slope of wax solubility curve at the oildeposit interface
The solubility of wax in oil has a close relation to the oil temperature. The slope of wax solubility curve at the oildeposit interface is essential to calculate the growth rate of wax deposition. And the greater the slope, the faster the wax deposition grows. The can be obtained from T_{interface}, knowing the solubility curve of the oil.
The slope of wax solubility curve at different temperature dC/dT is calculated, based on the solubility curve of wax molecules in oil in Figure 7. As shown in Figure 16, dC/dT has a linear increasing trend with the oil temperature decreasing. As shown in Figure 17, the slope of wax solubility curve dC/dT diminishes, with the position approaching to B along the wall, which leads to a decrease in the thickness of wax deposition.
Fig. 16 The slope of solubility curve of wax in oil at different temperature. 
Fig. 17 Circular distribution of the slope of solubility curve of wax in oil. 
Through the above analysis, it is known that the growth rate of wax deposition is controlled by the three factors. When the position approaches to point B along the wall where the gasliquid interface intersects with the wall, , the diffusion coefficient of wax molecules in oil, gradually increases; , the temperature gradient at oildeposit interface, and , the slope of wax solubility curve, decrease. And the increasing rate of is slower than the summing decreasing rate of and . As a result, J_{wax}, the convection diffusion mass flux of wax molecules in oil, decreases, resulting in the reduction of the growth rate of wax deposition. Therefore, for the researches on the growth rate of wax deposition, it is essential to have a comprehensive analysis of heat transfer, mass transfer and oil properties during wax deposition process.
4 Conclusion
A threedimensional wax deposition model of oilgas stratified flow in a circular pipe where deposition occurs at the lower part of the pipe has been developed based on the molecular diffusion mechanism. The nonisothermal hydrodynamics and the heat/mass transfer for stratified flow in a circular pipe are studied to predict the growth of the wax deposit along the circumferential distribution as a function of time. Predictions of wax deposition thickness along the circumferential pipe wall showed acceptable agreement with Matzain’s experimental results.
In addition, an analysis was carried out using the wax deposition model to explain why the deposit appears crescent shaped with two identifiable zones in the oilgas stratified flow pattern. The reason is indicated by the circumferential distribution of characteristic mass flux for wax deposition, J_{wax}, along the pipe wall wetted by the oil. The effect of the mass flux includes D_{wo} (T_{interface}), the diffusion coefficient of wax molecules in oil, , the temperature gradient at oildeposit interface, and , the slope of wax solubility curve. The molecular diffusivity rises gradually from the bottom of the pipe to the intersection of gasliquid interface and the pipe wall across the circumference. But the oildeposit interface temperature gradient and the slope of wax solubility curve in oil diminish. Besides, the increasing rate of is slower than the summing decreasing rate of and , which leads to diffusion mass flux of wax molecules J_{wax} decreases. It indicates that the growth rate of wax deposition should diminish.
Acknowledgments
The work received the support of by Natural Science Foundation of China (51604282), Chongqing Research Program of Basic Research and Frontier Technology (cstc2015jcyjB90001 and cstc2016jcyjA0095), China Postdoctoral Science Foundation (2017M623394 and 2018T111156), Chongqing Postdoctoral Research Project Funding (Xm2017050).
References
 Wang W., Huang Q., Wang C., Li S., Qu W., Zhao J., He M. (2015) Effect of operating conditions on wax deposition in a laboratory flow loop characterized with DSC technique, J. Therm. Anal. Calorim. 119, 1, 471–485. [Google Scholar]
 Zhou J., Liang G., Deng T., Zhou S. (2017) Optimization design of coalbed methane pipeline networkcoupled wellbore/reservoir simulation, Adv. Mech. Eng. 9, 6, 2071939978. [Google Scholar]
 Duan J., Liu H., Gong J., Jiao G. (2015) Heat transfer for fully developed stratified wavy gasliquid twophase flow in a circular crosssection receiver, Sol. Energy 118, 338–349. [Google Scholar]
 Zhou J., Liang G., Deng T. (2018) Optimal design of startree oilgas pipeline network in discrete space, J. Pipeline Syst. Eng. Pract. 9, 1, 1–11. [CrossRef] [Google Scholar]
 Lei Y., Han S., Zhang J., Bao Y., Yao Z., Xu Y.N. (2014) Study on the effect of dispersed and aggregated asphaltene on wax crystallization, gelation, and flow behavior of crude oil, Energy Fuels 28, 4, 2314–2321. [Google Scholar]
 Bai C., Zhang J. (2013) Thermal, macroscopic, and microscopic characteristics of wax deposits in field pipelines, Energy Fuels 27, 2, 752–759. [Google Scholar]
 Quan Q., Wang W., Wang P., Duan J., Yang J., Yao H., Gong J. (2015) Wax deposition from emulsionwater in stratified flow, Petrol. Sci. Technol. 33, 6, 749–755. [CrossRef] [Google Scholar]
 Matzain A., Apte M.S., Zhang H.Q., Volk M., Brill J.P., Creek J.L. (2002) Investigation of paraffin deposition during multiphase flow in pipelines and wellbores – Part 1: Experiments, J. Energy Resour. Technol. 124, 3, 180–187. [Google Scholar]
 Apte M.S., Matzain A., Zhang H.Q., Volk M., Brill J.P., Creek J.L. (2001) Investigation of paraffin deposition during multiphase flow in pipelines and wellbores – Part 2: Modeling, J. Energy Resour. Technol. 123, 2, 150–158. [Google Scholar]
 Rittirong A. (2014) Paraffin deposition under twophase gasoil slug flow in horizontal pipes, PhD Dissertation, The University of Tulsa. [Google Scholar]
 Cem S., Ekarit P. (2012) Review of paraffin deposition research under multiphase flow conditions, Energy Fuels 26, 7, 3968–3978. [Google Scholar]
 Gong J., Yu Z., Liao L., Duan J., Wang P., Zhou J. (2011) Wax deposition in the oil/gas twophase flow for a horizontal pipe, Energy Fuels 25, 4, 1624–1632. [Google Scholar]
 Huang Z., Lee H.S., Senra M., Scott F.H. (2011) A fundamental model of wax deposition in subsea oil pipelines, AIChE J. 57, 11, 2955–2964. [Google Scholar]
 Singh P., Venkatesan R., Scott F.H., Nagarajan N. (2000) Formation and aging of incipient thin film waxoil gels, AIChE J. 46, 5, 1059–1074. [Google Scholar]
 Bidmus H.O., Mehrotra A.K. (2012) Comments on “The Effect of Operating Temperatures on Wax Deposition” by Huang et al., Energy Fuels 26, 6, 3963–3966. [Google Scholar]
 Mehrotra A.K., Bhat N.V. (2010) Deposition from “Waxy” Mixtures under turbulent flow in pipelines: Inclusion of a viscoplastic deformation model for deposit aging, Energy Fuels 24, 4, 2240–2248. [Google Scholar]
 Duan J., Gong J., Yao H., Deng T., Zhou J. (2014) Numerical modeling for stratified gasliquid flow and heat transfer in pipeline, Appl. Energy 115, 83–94. [Google Scholar]
 Huang Z., Lee H.S., Senra M., Scott F.H. (2011) A fundamental model of wax deposition in Subsea oil pipelines, AIChE J. 57, 11, 2955–2964. [Google Scholar]
All Tables
All Figures
Fig. 1 Wax deposit and blockage in submarine pipeline. 

In the text 
Fig. 2 Illustration of wax deposit in oilgas stratified pipe flow. 

In the text 
Fig. 3 The cubic cage structure unit of deposit layer under static condition. 

In the text 
Fig. 4 The calculation process of wax deposit for oilgas stratified pipe flow. 

In the text 
Fig. 5 Wax thickness distribution along the pipe wall for oilgas stratified flow. 

In the text 
Fig. 6 Viscositytemperature curve of South Pelto crude oil. 

In the text 
Fig. 7 Solubilitytemperature curve of South Pelto crude oil. 

In the text 
Fig. 8 Comparison of distribution of wax deposition in the ring section between the experiment and prediction. 

In the text 
Fig. 9 Predicted angular distribution of deposit thickness, γ = 0° at the bottom of the pipe. 

In the text 
Fig. 10 Predicted angular distribution of wax content in the deposit, γ = 0° at the bottom of the pipe. 

In the text 
Fig. 11 Comparison between experimental deposit thickness in the bottom of pipe and the prediction. 

In the text 
Fig. 12 Comparison between experimental wax content in deposit in the bottom of pipe and the prediction. 

In the text 
Fig. 13 Circular distribution of oildeposit temperature. 

In the text 
Fig. 14 Circular distribution of wax molecular diffusion coefficient. 

In the text 
Fig. 15 Circular distribution of temperature gradient of oildeposit. 

In the text 
Fig. 16 The slope of solubility curve of wax in oil at different temperature. 

In the text 
Fig. 17 Circular distribution of the slope of solubility curve of wax in oil. 

In the text 