Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Experimental investigation of oil–water flow in the horizontal and vertical sections of a continuous transportation pipe

## Abstract

A series of experiments were conducted to investigate flow pattern transitions and concentration distribution during simultaneous pipe flow of oil–water two-phase flow through the horizontal and vertical sections. The flowing media applied were white mineral oil and distilled water. Superficial oil and water velocities were between 0 and 0.57 m/s. Flow pattern maps revealed that the horizontal and vertical sections of the pipe lead to different flow pattern characteristics under the same flow conditions. The original contributions of this work are that a transition mechanism for predicting the boundary between oil-in-water (O/W) flow and water-in-oil (W/O) in oil–water two-phase flow was obtained. The effects of input water cut, oil and water superficial velocities on the concentration distribution of the dispersed phase were studied. The empirical formulas for the phase holdup based on the drift-flux model were obtained. The predicted results agreed well with those of the experimental data, especially for the O/W flow pattern.

## Introduction

With the increase of oil well production life, the volume fraction of water in oil pipes may gradually increase1. The oil–water two-phase flow in the pipe is a crucial topic in the petroleum industry. The oil–water flow is characterized by a variety of diverse flow patterns. Different concentration and pressure drop distribution characteristics occurred under different flow patterns. The understanding of flow pattern and phase holdup is of great importance for the evaluation of oilfield operating and oil production. The change of the flow pattern has a vital impact on the pressure drop, spatial phase distribution, and safety of pipeline transportation2. The calculation of water holdup is helpful to predict the oil transportation in oil pipelines, and water holdup is an important parameter in identifying flow patterns3.

Many researchers have been studied the simultaneous oil and water flow in horizontal or vertical pipes. For horizontal oil–water flows, the flow patterns are mainly divided into the segregated flow and the dispersed flow. The segregated flow did not occur in pipes with an inclination of more than 33·. The oil–water flow often occurs as the dispersed flow4. Piela et al.5 studied the phase inversion of oil–water flow experimentally through a horizontal pipe loop. The experiments were carried in a pipe loop with an inner pipe diameter of 16 mm, consisting of two straight parts of 6 m connected via two bends. Experimental data showed that the phase inversion occurred at the oil volume fraction of 0.9 for the water-to-oil inversion, and the pressure drop increases at the inversion compared to the initial pressure drop. The volume fraction of the dispersed phase at the water-to-oil experiment transition point was higher than that for oil-to-water experiments. Gong et al.6 developed a model for predicting phase inversion in oil–water flow in a stainless steel pipe loop using the relations between the surface energies before and after phase inversion. The results of the model were in reasonable agreement with the experimental results. The experiments were conducted in a stainless steel pipe loop (25.7 mm inner diameter, 52 m long). The experimental temperature was 40 °C, 50 °C, and 60 °C, respectively. Kumara et al. measured an oil–water two-phase flow in a horizontal pipe with a diameter of 56 mm and length of 15 m with particle image velocimetry and laser Doppler anemometry7. Zhou et al. designed parallel-wire conductivity probes to measure water hold-up of oil–water two-phase flow in a near-horizontal pipe. The experimental data showed that the method can achieve a good measurement result of water hold-up8. Wang et al. conducted oil–water experiments in a small horizontal channel with 2 mm. Flow patterns and the plug shape/length taken by highspeed camera were investigated9.

For vertical oil–water systems, Hasan and Kablr10 adopted the drift-flux model to predict in situ oil volume fraction in oil–water flow. The model parameters were determined by the experimental data. The experimental rig consisted of two separate 5.5 m long transparent columns with 64 and 127 mm inner diameters. Xu et al.11 conducted an experimental study of a co-current upward and downward in a transparent vertical pipe. The phase inversion and frictional pressure gradients were investigated using white oil and water in a 50 mm inner diameter pipe. The flow patterns are usually recognized by visual inspection and the fluctuations in the volume fraction. Keska et al.12 carried out the comparative experimental research on four methods (capacitive, resistive, optical, and pressure) commonly adopted to identify the flow pattern, and the results proved that the best way to judge the flow pattern is to use the capacitance or resistance fluctuation of the flow cross-section. The values of resistance and capacitance represent the change in the volume fraction. Jones and Zuber13 demonstrated that the probability density function (PDF) of the fluctuations in the volume fraction can be used as a statistical analysis tool for flow pattern identification. The PDF is a function of signal amplitude and a method of measuring the probability that a signal has a range of values. In this study, the PDF of the volume fraction of the water phase is used to identify changes in the flow pattern. Du et al. employed convolutional neural networks to identify oil–water two-phase flow patterns. The different flow pattern images were collected through oil-in-water flow experiments in a vertical 20 mm inner diameter Plexiglas pipe. The flow pattern images were collected by the highspeed camera14.

However, few comparative studies of horizontal and vertical flows in a combined pipe have been reported. In the on-site oil pipeline layout, the pipeline is horizontal, vertical, or even inclined. Under the same flow parameters, pipelines in different directions may show different flow pattern characteristics, which brings hidden dangers to the stability of pipeline transportation. Under the same flow parameters, it is necessary to study the different flow behaviors in horizontal and vertical pipelines3.

The oil–water flow in the pipe is complex and related to pipe geometries (e.g., inner pipe diameter, pipe angle), fluid properties (e.g., viscosity, density, and surface tension), and boundary conditions (e.g., superficial input velocities). Previous studies mainly focused on the effect of a single parameter, but the coupled effect of the controlling parameters on the flow is not well understood. Hence, to understand the fundamental mechanisms of oil–water flows, the controlling dimensionless parameters were derived by dimensional analysis first.

Thus, the objective of this work was to investigate the flow pattern transition and phase holdup for an oil–water two-phase flow through the horizontal and vertical sections, a comparative study of the flow behavior has been conducted in a pipe loop by using different superficial input velocities. A series of experiments were conducted to investigate the oil–water flow pattern transition and the water holdup, considering both the horizontal and vertical sections of a transportation pipe.

The current paper is structured as follows. The experimental setup is described in “Experimental setup and procedure” section. The dimensional analysis is presented in “Dimensional analysis and data structure” section. The overall results in terms of flow pattern and phase holdup are discussed in “Results and discussion” section. The relationship between the transient water holdup and the change of the flow pattern in a transportation pipe with horizontal and vertical sections is established, and the empirical formulas for the phase holdup based on the drift-flux model are obtained. Finally, “Conclusion” section presents the conclusions of this study.

## Experimental setup and procedure

### Experimental facility

The experimental study for oil–water flow was conducted in a closed circuit loop. The schematic drawing of the flow loop is depicted in Fig. 1, which can be used to carry out both oil–water or oil–water–gas flow experiments. The test equipment can identify the oil–water flow patterns and record the water holdup (in situ volume fraction of water) and the pressure drop.

The flowing pipe was composed of a stainless steel pipe section and a plexiglass pipe section, and its inner diameter was 50 mm. The flow pattern development section was a horizontal stainless steel pipe with a length of 10 m from the entrance (corresponding to 200D, where D was the inner diameter). Based on the experimental observation, this length ensured that a fully developed flow can be achieved before measurement6. Due to the mixture velocities were low (u < 0.57 m/s) and oil phase viscosity high (32 mPa s), this length provided sufficient length to stabilize the flow. Flores2 also used the development length of about 200 diameters to study the characterization of oil–water flow patterns. The reasonableness of the development length is further confirmed in “Flow pattern maps” section.

The flow observation sections consisted of a 1 m long horizontal transparent pipe and a 2 m long vertical transparent pipe. With this arrangement, experiments on horizontal and vertical flows can be carried out simultaneously. The water was pumped first by a centrifugal pump (QABP160M2A, ABB), with a capacity of 12.5 m3/h. Then the oil phase was pumped into the clapboard mixer in which the oil was mixed with the water. A 6.99 KW gear pump (SNH440) with a capacity of 17 m3/h and accuracy of ± 0.1% of the reading was used for the oil phase. The oil and water phases entered the double-Y junction fitting mixer from the upper and lower layers of the mixer, respectively. The angle of the meeting point of the oil and water flow lines before entering the mixer is 60 degrees. The oil–water mixture flowed simultaneously to the test section and then flowed down the return section to the separating tank. Leaving the separating tank, the oil and water were directed to their original tanks. The liquid volume flow rate was monitored by a mass flowmeter (CMF100, Micro Motion), whose accuracy was ± 0.1%.

### Measuring equipment

Measurements included water holdup, pressure drop, and flow patterns, in addition to recording parameters such as the volume flow rate of each phase and the experimental temperatures. A differential pressure transducer was installed in the return section to measure the pressure drop at a frequency of 1 Hz. When the pressure drop was independent of time, it was deemed that a steady state of the system was reached. Then a total of 300 s of data was recorded. A digital video was used to observe the flow structures and identify the flow patterns. The working fluids used in the experiments were white mineral oil and distilled water. All experiments were carried out under atmospheric pressure, and the room temperature was controlled at about 20 °C. The physical parameters of the oil phase and water phase are shown in Table 1.

Two water holdup instruments equipped with a conductance probe were assembled in the horizontal and vertical sections of the test pipe to measure the cross-sectional average water holdup $$\alpha_{{\text{w}}}$$ at a sampling frequency of 1 Hz, respectively. The photograph of the water holdup instrument is shown in Fig. 23. The conductance probe measured the voltage between the two ends of the conductor, which allowed the mean conductivity of the mixture in the pipe to be calculated. Due to the poor conductivity of the oil phase, the voltage value measured when the pipe was filled with pure oil was much higher than that when the pipe was filled with pure water. The voltage values of pure oil and water were denoted by Vo and Vw, respectively. Then, the mean conductivity Vexp of the in-suit liquid can be determined by the following equation:

$$V_{\exp } = f\left( {\alpha_{{\text{w}}} ,\;V_{{\text{w}}} ,\;V_{{\text{o}}} } \right)$$
(1)

Taking Vw as a unit system produced the following dimensionless relationship:

$$\frac{{V_{\exp } }}{{V_{{\text{w}}} }} = f\left( {\alpha_{{\text{w}}} ,\frac{{V_{{\text{o}}} }}{{V_{{\text{w}}} }}} \right)$$
(2)

According to the mixture theory15, it was assumed that the voltage value of the mixture satisfied the following relationship:

$$\frac{{V_{\exp } }}{{V_{{\text{w}}} }} = \alpha_{{\text{w}}} + \left( {1 - \alpha_{{\text{w}}} } \right)\frac{{V_{{\text{o}}} }}{{V_{{\text{w}}} }}$$
(3)

The relationship between water holdup $$\alpha_{{\text{w}}}$$ and the mean conductivity of the mixture Vexp can be expressed as:

$$\alpha_{{\text{w}}} = \frac{{V_{{\text{o}}} - V_{\exp } }}{{V_{{\text{o}}} - V_{{\text{w}}} }}$$
(4)

Noting that the voltage values of pure water and pure oil were slightly different under different experimental conditions. Therefore, both Vo and Vw adopt the values measured under the current experimental condition when the voltage value signal was converted into the water holdup using Eq. (4).

## Dimensional analysis and data structure

The flow patterns of oil–water two-phase flow are classically organized into the separated flow and dispersed flow4. The dispersed flow frequently occurs in oil–water flow, where one phase is dispersed phase, and the other is continuous phase. These two types of dispersion spontaneously invert at some operational conditions, which is called phase inversion. A single uppercase subscript C refers to the continuous phase and D refers to the dispersed phase. The parameters controlling the oil–water two-phase dispersed flow are listed as follows:

Dispersed phase: density $$\rho_{D}$$, viscosity $$\mu_{D}$$, drop diameter $$d_{D}$$, interfacial tension $$\sigma$$.

Continuous phase: density $$\rho_{C}$$, viscosity $$\mu_{C}$$.

Geometric parameter: inner pipe diameter D.

Boundary condition: superficial input velocity $$u_{s}$$ which is defined as $$u_{s} = Q{/}A$$ where Q is the volume flow rate of each phase and A is the cross-sectional area of the pipe, $$A = \pi D^{2} /4$$.

Gravitational acceleration: g.

The eventual steady state of the oil–water flow system, characterized by the holdup of the dispersed phase $$\alpha_{D}$$, pressure drop $$\Delta P$$, and flow pattern, is a function of the above controlling parameters:

$$\left( {\alpha_{D} ,\;\Delta P,\;flow\;pattern} \right) = f\left( {\rho_{C} , \, \mu_{C} , \, u_{sC} ; \, \rho_{D} , \, \mu_{D} , \, u_{sD} , \, d_{D} , \, \sigma ; \, D, \, g} \right)$$
(5)

Equation (5) can be nondimensionalized as:

$$\left( {\alpha_{D} ,\;\frac{\Delta P}{{\rho_{C} u_{sC}^{2} }}, \, flow\;pattern} \right) = f\left( {\frac{{u_{sD} }}{{u_{sC} }}, \, \frac{{\rho_{C} u_{sC} D}}{{\mu_{C} }}, \, \frac{{u_{sC}^{2} }}{gD}, \, \frac{{\left| {\rho_{C} - \rho_{D} } \right|gD^{2} }}{\sigma }, \, \frac{{\mu_{D} }}{{\mu_{C} }}, \, \frac{{\rho_{D} }}{{\rho_{C} }}, \, \frac{{d_{D} }}{D}} \right)$$
(6)

where $$u_{sD} /u_{sC}$$ is the superficial velocity ratio, $$\rho_{C} u_{sC} D/\mu_{C}$$ is Reynolds number $$Re_{C}$$, $$u_{sC}^{2} /gD$$ is Froude number FrC. $${{\left( {\left| {\rho_{C} - \rho_{D} } \right|gD^{2} } \right)} \mathord{\left/ {\vphantom {{\left( {\left| {\rho_{C} - \rho_{D} } \right|gD^{2} } \right)} \sigma }} \right. \kern-\nulldelimiterspace} \sigma }$$ is Eötvös number Eo that represents the ratio of buoyancy force to surface tension force. For a wholly developed dispersed flow, the velocity difference between oil and water phases is minimal. Then, in the present paper, the dimensionless numbers are calculated instead of the superficial velocity $$u_{sC}$$, with the total superficial velocity $$u_{{s{\text{m}}}}$$.

$$u_{{s{\text{m}}}} = u_{sC} + u_{sD}$$
(7)

Defining input water cut $$\varepsilon_{{\text{w}}}$$ as the ratio of the water flow rate to the mixture flow rate.

$$\varepsilon_{{\text{w}}} = \frac{{Q_{{\text{w}}} }}{{Q_{{\text{w}}} + Q_{{\text{o}}} }} = \frac{{u_{{s{\text{w}}}} }}{{u_{{s{\text{m}}}} }}$$
(8)

Assuring that the fluid properties and drop size remain constant, i.e., Eötvös number Eo, density ratio $$\rho_{D} /\rho_{C}$$, viscosity ratio $$\mu_{D} /\mu_{C}$$ , and the ratio of drop size to pipe diameter $$d_{D} /D$$ are fixed. Therefore, the Eq. (6) can be modified as:

$$\left( {\alpha_{D} ,\;\frac{\Delta P}{{\rho_{C} u_{{s{\text{m}}}}^{2} }},\; \, flow\;pattern} \right) = f\left( {\varepsilon_{{\text{w}}} , \, \frac{{\rho_{C} u_{{s{\text{m}}}} D}}{{\mu_{C} }}, \, \frac{{u_{{s{\text{m}}}}^{2} }}{gD}} \right)$$
(9)

The effects of the dimensionless parameters in Eq. (9) on the flow characteristics of the oil–water two-phase flow are investigated. All the experimental values of $$u_{{s{\text{w}}}}$$, $$u_{{s{\text{o}}}}$$, and the corresponding dimensionless numbers are given in Table 2. Noting that $$Re_{{\text{m}}}$$ is not listed in this table and calculated according to which phase is continuous in the experiments. The material-related dimensionless parameters remain constant.

## Results and discussion

### Flow pattern maps

The experiments of wash-out type oil–water two-phase pipe flow were carried out by starting with pumping water first and then adding oil applying different flow rates. The identification of flow regimes was based on both visual observation and interpretation of transient flow signals (e.g., fluctuations of water holdup and pressure drop). Flores4 mainly identified water-dominated pattern (dispersion of oil in water, O/W) and oil-dominated pattern (dispersion of water in oil, W/O) for oil–water flow in vertical and horizontal pipes. The present work focused on the phase inversion between oil-in-water and water-in-oil dispersions.

Figure 3 shows examples of the flow patterns that were observed in this study in both the horizontal and vertical sections in a continuous transportation pipe. Figure 3a presents the O/W pattern. Figure 3b shows the W/O pattern, where a liquid film with a certain thickness formed between the gas column and the pipe wall. Figure 3c shows the stratified flow. The upper layer was the continuous oil phase and the lower layer was the continuous water phase. Based on the experimental observations, the development length can achieve the fully developed flow patterns.

Figures 4 and 5 shows two examples of time series of the cross-sectional average water holdup $$\alpha_{{\text{w}}}$$ and the corresponding PDF as a function of the water holdup for different experimental conditions. The PDF shows a single narrow peak in both the horizontal and vertical sections, which indicates a small variance in water holdup. This is consistent with the characteristics of the dispersed flow and stratified flow2,13. Figure 4 shows the data for $$u_{{s{\text{w}}}} = 0.16$$ m/s and $$u_{{s{\text{o}}}} = 0.07$$ m/s. O/W pattern was observed in the horizontal and vertical sections. The O/W pattern was characterized by one narrow peak in the PDF at $$\alpha_{{\text{w}}} = 0.72$$ in this case. Figure 5 shows the result for $$u_{{s{\text{w}}}} = 0.02$$ m/s and $$u_{{s{\text{o}}}} = 0.20$$ m/s. W/O pattern was observed in the horizontal and vertical sections. The PDF showed a single narrow at a low input water holdup ($$\alpha_{{\text{w}}} < 0.10$$).

The phase inversion point between O/W and W/O is usually defined as the critical holdup of the dispersed phase above which the dispersed phase becomes the continuous phase. Around the point of phase inversion, the effective mixture viscosity and the pressure drop fluctuate violently6,15. Figure 6 shows the scaled pressure drop (ratio of the measured pressure drop and the pressure drop of pure water) as a function of the oil holdup $$\alpha_{{\text{o}}}$$ for three different total superficial velocities. The pressure drop increases slightly to a maximum value with increasing the oil holdup and then reduces immediately after passes through the maximum. In this study, phase inversion occurs at an oil holdup between 0.75 and 0.95. It can also be seen from Fig. 6 that the pressure drop decreases and the critical oil holdup at phase inversion increases with increasing total superficial velocities.

The flow pattern maps in terms of superficial velocities of oil and water studied in this paper are demonstrated in Fig. 7. Stratified flow (SF) is observed at high water velocities and low oil velocities in the horizontal section. For decreasing water velocities and increasing oil velocities, the transition to O/W takes place. In the vertical section, only O/W and W/O flow patterns are observed. A comparison between Fig. 7a,b indicates there is a transition from horizontal stratified flow to vertical dispersed flow at the same superficial velocities. The region of O/W flow is larger than that of W/O flow, and the transition boundary from O/W and W/O is almost identical in both the horizontal and vertical sections. The flow patterns in this experiment are obtained with $$\mu_{{\text{o}}} /\mu_{{\text{w}}} = 32$$ and $$\rho_{{\text{o}}} /\rho_{{\text{w}}} = 0.84$$.

### Flow pattern transition mechanism

The boundary of the dispersed flow depends on the balance of the turbulent dispersive forces that tend to deform the drop and the surface tension that tend to resist fission. Brauner16 proposed a unified approach for predicting the transitions to dispersed flow patterns in gas–liquid and liquid–liquid systems. The transition criteria are based on the part of the work of Hinze17. In an oil–water system, the transition to O/W takes place when the water turbulence is strong enough to decompose the oil into droplets smaller than the critical size $$d_{{{\text{crit}}}}$$:

$$d_{\max } \le d_{{{\text{crit}}}}$$
(10)

where dmax is the maximum diameter that the droplet can reach under the combined action of shear force and surface tension. In dense dispersions, the maximum diameter is given by:

$$\left( {\frac{{d_{\max } }}{D}} \right) = 7.61C_{{\text{H}}}^{0.6} \left( {\frac{{\rho_{{\text{w}}} u_{{s{\text{m}}}}^{2} D}}{\sigma }} \right)^{ - 0.6} \left( {\frac{{\rho_{{\text{w}}} u_{{s{\text{m}}}} D}}{{\mu_{{\text{w}}} }}} \right)^{0.08} \left( {\frac{{u_{{s{\text{o}}}} }}{{u_{{s{\text{w}}}} }}} \right)^{0.6} \left( {1 + \frac{{\rho_{{\text{o}}} u_{{s{\text{o}}}} }}{{\rho_{{\text{w}}} u_{{s{\text{w}}}} }}} \right)^{ - 0.4}$$
(11)

where $$C_{{\text{H}}}$$ is a constant to be determined experimentally, $$d_{{{\text{crit}}}}$$ refers to the maximum diameter of the droplet without deformation, which can be expressed as:

$$\frac{{d_{{{\text{crit}}}} }}{D} = \frac{0.224}{{\left( {\Delta \rho gD^{2} /8\sigma } \right)^{0.5} }}$$
(12)

The scope of application of the above transition criterion is ρwusmD/µw ≥ 2100 and $$1.82\left( {{{\rho_{{\text{w}}} u_{{s{\text{m}}}} D} \mathord{\left/ {\vphantom {{\rho_{{\text{w}}} u_{{s{\text{m}}}} D} {\mu_{{\text{w}}} }}} \right. \kern-\nulldelimiterspace} {\mu_{{\text{w}}} }}} \right)^{ - 0.7} < {{d_{{{\text{crit}}}} } \mathord{\left/ {\vphantom {{d_{{{\text{crit}}}} } D}} \right. \kern-\nulldelimiterspace} D} < 0.1$$.

The solid line in Fig. 7 represents the boundary that corresponds to the results of the criterion (10). In our experiments, it can better represent the trend from O/W to W/O when $$C_{{\text{H}}} = 0.012$$. The expression of the critical curve is:

$$\left( {u_{{s{\text{o}}}} + u_{{s{\text{w}}}} } \right)^{1.12} \left( {\frac{{u_{{s{\text{w}}}} }}{{u_{{s{\text{o}}}} }}} \right)^{0.6} \left( {1 + \frac{{\rho_{{\text{o}}} u_{{s{\text{o}}}} }}{{\rho_{{\text{w}}} u_{{s{\text{w}}}} }}} \right)^{0.4} = 0.845\frac{{(\Delta \rho g)^{0.5} \sigma^{0.1} D^{0.48} }}{{\rho_{{\text{w}}}^{0.52} \mu_{{\text{w}}}^{0.08} }}$$
(13)

As can be seen in Fig. 7, Eq. (13) can provide somewhat reasonable predictions for phase inversion, although more experimental data are needed to improve the accuracy of the model prediction.

### Phase holdup

The input water cut has an essential influence on the water holdup in the test section. Figures 8, 9 and 10 shows the average water holdup of the time series in the horizontal and vertical sections at different total superficial velocities (corresponding to different $$Re_{{\text{m}}}$$ and $$Fr_{{\text{m}}}$$). The dotted lines in Figs. 8, 9 and 10 represent the boundary between W/O and O/W flow patterns. In the range of W/O, oil forms the continuous phase $$Re_{{\text{m}}} = \rho_{{\text{o}}} u_{{s{\text{m}}}} D/\mu_{{\text{o}}}$$, whereas, in the range of O/W, water forms the continuous phase $$Re_{{\text{m}}} = \rho_{{\text{w}}} u_{{s{\text{m}}}} D/\mu_{{\text{w}}}$$. For $$Fr_{{\text{m}}} = 0.105$$ and 0.216, there are no significant differences in water holdup between horizontal and vertical sections (Figs. 8 and 9). At lower input water cut corresponding to experimental data, the water holdup is slightly lower than the input water cut, i.e., there is a significant velocity difference between the oil and water phases. The range with a smaller water holdup corresponds to the W/O flow pattern range. In the range of the O/W flow pattern, the water holdup is larger than the input water cut. The transition point appears at an input water cut of about 0.2. For $$Fr_{{\text{m}}} = 0.654$$, the transition point of water holdup from less than input water cut to higher than the input water cut appears at an input water cut of about 0.35 in the vertical section. However, the water holdup is considerably lower than the input water cut in the horizontal section due to water flows faster at this condition (Fig. 10). It also can be seen that for dispersed O/W flow with high input water cut ($$\varepsilon_{{\text{w}}} = 0.9$$) the dispersed oil droplets travel at approximately the same velocity as the water. The velocity difference between oil and water is close to zero. The result is also similar to that in gas–liquid flow2.

The concentration variation of the dispersed phase in the pipe can be predicted by the drift-flux model proposed by Zuber and Findlay18. Hasan and Kablr10 successfully applied the drift-flux model to the oil–water two-phase flow to predict the concentration of the dispersed phase. For the O/W flow pattern, the drift-flux model calculates the oil holdup:

$$\frac{{u_{{s{\text{o}}}} }}{{\alpha_{{\text{o}}} }} = Cu_{{s{\text{m}}}} + \frac{{J_{d} }}{{\alpha_{o} }}$$
(14)

where $$J_{d}$$ is the drift flux that can be calculated by the following semi-theoretical relationship:

$$J_{d} = u_{\infty } \alpha_{{\text{o}}} \left( {1 - \alpha_{{\text{o}}} } \right)^{n}$$
(15)

where C and $$n$$ are parameters determined by experiments, $$u_{\infty }$$ is the terminal rise velocity of an oil droplet in the static water column, which is calculated by the Harmathy correlation19.

$$u_{\infty } = 1.53\left[ {{{\sigma \left( {\rho_{{\text{w}}} - \rho_{{\text{o}}} } \right)g} \mathord{\left/ {\vphantom {{\sigma \left( {\rho_{{\text{w}}} - \rho_{{\text{o}}} } \right)g} {\rho_{{\text{w}}}^{2} }}} \right. \kern-\nulldelimiterspace} {\rho_{{\text{w}}}^{2} }}} \right]^{0.25}$$
(16)

Substituting Eq. (15) into Eq. (14) yields the following correlation:

$${{u_{{s{\text{o}}}} } \mathord{\left/ {\vphantom {{u_{{s{\text{o}}}} } {\alpha_{{\text{o}}} }}} \right. \kern-\nulldelimiterspace} {\alpha_{{\text{o}}} }} = Cu_{{s{\text{m}}}} + u_{\infty } \left( {1 - \alpha_{{\text{o}}} } \right)^{n}$$
(17)

Similarly, for the W/O flow pattern, the drift-flux model is transformed into the following form:

$${{u_{{s{\text{w}}}} } \mathord{\left/ {\vphantom {{u_{{s{\text{w}}}} } {\alpha_{{\text{w}}} }}} \right. \kern-\nulldelimiterspace} {\alpha_{{\text{w}}} }} = Cu_{{s{\text{m}}}} + u_{\infty l} \left( {1 - \alpha_{{\text{w}}} } \right)^{n}$$
(18)

where $$u_{\infty l}$$ is the terminal settlement velocity of a water droplet in static oil column which is calculated by:

$$u_{\infty l} = 1.53\left[ {{{\sigma \left( {\rho_{{\text{w}}} - \rho_{{\text{o}}} } \right)g} \mathord{\left/ {\vphantom {{\sigma \left( {\rho_{{\text{w}}} - \rho_{{\text{o}}} } \right)g} {\rho_{{\text{o}}}^{2} }}} \right. \kern-\nulldelimiterspace} {\rho_{{\text{o}}}^{2} }}} \right]^{0.25}$$
(19)

The values of C and n can be obtained by curve fitting analysis of experimental data. The fitting results for the horizontal and vertical sections are presented in Table 3. The standard deviation (SD) of the predicted value can be expressed as20:

$${\text{SD}} = \sqrt {\frac{1}{n - 1}\sum\limits_{k = 1}^{n} {\left( {\frac{{\left( {\alpha_{D} } \right)_{{{\text{pred}}}} - \left( {\alpha_{D} } \right)_{\exp } }}{{\left( {\alpha_{D} } \right)_{\exp } }}} \right)^{2} } }$$
(20)

Figure 11 depicts a comparison of the predicted dispersed phase holdup with the experimental data. The square and circular symbols represent the data of the horizontal section and the vertical section, respectively. Better fitting results are obtained for the O/W flow pattern, and the standard deviation of fitting data is within 4.5%. In the range of the W/O flow pattern, a good agreement is obtained for the experimental data of the vertical section at a standard deviation of 6.49%. However, the data in the horizontal section cannot be predicted by the drift-flux model. The reason may be due to the measured holdup data for the horizontal section in the low input water cut are considered to be less reliable. Overall, the existing model provides good estimates of the oil holdup in the dispersion of oil in water flow.

## Conclusion

The oil–water two-phase flow experiments have been conducted in the horizontal and vertical sections of the pipe simultaneously. The difference of oil–water two-phase distribution in horizontal and vertical sections is investigated and discussed. The newly designed water holdup instrument equipped with a conductance probe is used to measure the cross-sectional average water holdup.

The dispersed flow (W/O and O/W) and stratified flow have been identified over the range of oil and water superficial velocities. The pressure drop reaches the maximum value at the phase transition point of O/W and W/O in this study. The phase transition occurs at the oil holdup between 0.75 and 0.95. Based on the model of Brauner16, a unified correlation for predicting the transition between O/W and W/O is suggested. The model can provide a reasonable prediction for phase inversion. The drift-flux model proposed by Zuber and Findlay18 is modified for predicting the concentration of the dispersed phase in the horizontal and vertical sections. The model can well predict the concentration distribution of the O/W flow patterns. For the W/O flow pattern, the model is suitable for predicting the concentration change in the vertical section but fails to predict the concentration change in the horizontal section.

## References

1. 1.

Ghorai, S., Suri, V. & Nigam, K. D. P. Numerical modeling of three-phase stratified flow in pipes. Chem. Eng. Sci. 60, 6637–6648. https://doi.org/10.1016/j.ces.2005.05.050 (2005).

2. 2.

Ibarra, R., Nossen, J. & Tutkun, M. Two-phase gas–liquid flow in concentric and fully eccentric annuli. Part I: Flow patterns, holdup, slip ratio and pressure gradient. Chem. Eng. Sci. 203, 489–500. https://doi.org/10.1016/j.ces.2019.01.064 (2019).

3. 3.

Ren, G. et al. The flow pattern transition and water holdup of gas–liquid flow in the horizontal and vertical sections of a continuous transportation pipe. Water 13, 2077. https://doi.org/10.3390/w13152077 (2021).

4. 4.

Flores, J. G. Characterization of oil-water flow patterns in vertical and deviated wells. Soc. Petrol. Eng. 38810 (1997).

5. 5.

Piela, K. et al. Experimental investigation of phase inversion in an oil–water flow through a horizontal pipe loop. Int. J. Multiph. Flow 32, 1087–1099. https://doi.org/10.1016/j.ijmultiphaseflow.2006.05.001 (2006).

6. 6.

Gong, J., Li, Q. P., Yao, H. Y. & Yu, D. A model for predicting phase inversion in oil-water two-phase pipe flow. J. Hydrodyn. 18, 310–314 (2006).

7. 7.

Kumara, W. A. S., Elseth, G., Halvorsen, B. M. & Melaaen, M. C. Comparison of particle image velocimetry and laser doppler anemometry measurement methods applied to the oil–water flow in horizontal pipe. Flow Meas. Instrum. 21, 105–117. https://doi.org/10.1016/j.flowmeasinst.2010.01.005 (2010).

8. 8.

Zhou, Y. Y., Jin, N. D., Zhang, H. X. & Zhai, L. S. Method based on parallel-wire conductivity probe for measuring water hold-up in nearhorizontal oil–water two-phase flow pipes. IET Sci. Meas. Technol. 14, 676–683. https://doi.org/10.1049/iet-smt.2019.0308 (2020).

9. 9.

Wang, C. L., Tian, M. C., Zhang, J. Z. & Zhang, G. M. Experimental study on liquid–liquid two-phase flow patterns and plug hydrodynamics in a small channel. Exp. Thermal Fluid Sci. 129, 110455. https://doi.org/10.1016/j.expthermflusci.2021.110455 (2021).

10. 10.

Hasan, A. R. & Kablr, C. S. A new model for two-phase oil/water flow: production log interpretation and tubular calculations. Soc. Petrol. Eng. SPE 5, 193–199 (1990).

11. 11.

Xu, J. Y., Li, D. H., Guo, J. & Wu, Y. X. Investigations of phase inversion and frictional pressure gradients in upward and downward oil–water flow in vertical pipes. Int. J. Multiph. Flow 36, 930–939. https://doi.org/10.1016/j.ijmultiphaseflow.2010.08.007 (2010).

12. 12.

Keska, J. K., Smith, M. D. & Williams, B. E. Comparison study of a cluster of four dynamic flow pattern discrimination techniques for multiphase flow. Flow Meas. Instrum. 10, 65–77 (1999).

13. 13.

Jones, O. C. & Zuber, N. The interrelation between void fraction fluctuations and flow patterns in two-phase flow. Int. J. Multiph. Flow 2, 273–306 (1975).

14. 14.

Du, M., Yin, H. G., Chen, X. Y. & Wang, X. Q. Oil-in-water two-phase flow pattern identification from experimental snapshots using convolutional neural network. IEEE Access 7, 6219–6225. https://doi.org/10.1109/access.2018.2888733 (2019).

15. 15.

Descamps, M. N., Oliemans, R. V. A., Ooms, G. & Mudde, R. F. Experimental investigation of three-phase flow in a vertical pipe: Local characteristics of the gas phase for gas-lift conditions. Int. J. Multiph. Flow 33, 1205–1221. https://doi.org/10.1016/j.ijmultiphaseflow.2007.06.001 (2007).

16. 16.

Brauner, N. The prediction of dispersed flows boundaries in liquid–liquid and gas–liquid systems. Int. J. Multiphase Flow (2001).

17. 17.

Hinze, J. O. Fundamentals of the hydrodynamic mechanism of splitting in dispersion processes. AIChE J. 1, 289–295. https://doi.org/10.1002/aic.690010303 (1955).

18. 18.

Zuber, N. & Findlay, J. A. Average volumetric concentration in 2-phase flow systems. J. Heat Transf. 87, 453. https://doi.org/10.1115/1.3689137 (1965).

19. 19.

Harmathy, T. Z. Velocity of lager drops and bubbles in media of infinite or restricted extent. AIChE J. 6, 281–288. https://doi.org/10.1002/aic.690060222 (1960).

20. 20.

Al-Wahaibi, T. Pressure gradient correlation for oil–water separated flow in horizontal pipes. Exp. Thermal Fluid Sci. 42, 196–203. https://doi.org/10.1016/j.expthermflusci.2012.04.021 (2012).

## Author information

Authors

### Contributions

All the authors have contributed equally to this manuscript.

### Corresponding author

Correspondence to Peng Li.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Yang, J., Li, P., Zhang, X. et al. Experimental investigation of oil–water flow in the horizontal and vertical sections of a continuous transportation pipe. Sci Rep 11, 20092 (2021). https://doi.org/10.1038/s41598-021-99660-8