# A Unified Model of Oil/Water Two-Phase Flow through the Complex Pipeline

Authors: Qingchun Gao; Zhiming Wang; Quanshu Zeng
Journal: Geofluids
 (2021)

---

## Abstract

Oil-water two-phase flow through the complex pipeline, consisting of varying pipes and fittings in series or parallel, is commonly encountered in the petroleum industry. However, the majority of the current study is mainly limited to single constant-radius pipe. In this paper, a unified model of oil-water two-phase flow in the complex pipeline is developed based on the combination of pipe serial-parallel theory, flow pattern transformation criterion, two-fluid model, and homogenous model. A case is present to verify the unified model and compare with CFD results. The results show that the proposed unified model can achieve excellent performance in predicting both the flow distributions and pressure drops of oil-water two-phase flow in the complex pipeline. Compared with CFD results for water volumetric fractions ranging from 0% to 100%, the highest absolute percentage error of the proposed model is 14.4% and the average is 9.8%.


---

## Body

## 1. Introduction

Oil-water two-phase flow in the complex pipeline, consisting of varying pipes and fittings in series or parallel, is commonly encountered in the petroleum industry [1, 2]. Study on oil-water two-phase flow law in complex pipeline, especially the diversion, and the pressure drop, can not only solve the important problems in the process of oil production but also has important significance to improve the multiphase flow theory. However, the study is mainly limited to single constant-radius pipe at present [3, 4], and the oil-water two-phase flow in complex pipelines still needs further study.Up to date, research focused on oil-water two-phase flow performance can be categorized by the following three types, encompassing theoretical mathematical model, numerical simulation, and laboratory experiments. Thanks to the rapid advances in computational capacity in recent years, the numerical simulation approach has developed at an impressive speed [5, 6]. Several physics fully-coupled commercial numerical simulator have been massively utilized for analysis, like CFD and FLUENT [7]. However, it is still challenging to completely mimic oil-water flow behavior inside realistic complex pipelines with the use of numerical simulation. Also, the simulation model requires total re-establishment when the realistic case evolves, which aggravates its complexity [8, 9]. As for the laboratory experiments, relevant results can be directly examined and collected, which is evidently its advantage. Notably, the accuracy or reliability of the experimental data is heavily dependent on the quality of apparatus and examination method, suffering strong uncertainty [10–12]. In contrast, the mathematical model derived from classic theory not only provides instanton reliable results but also facilitates to identify the effect of each influential factor [13]. Motivated by the above perspective, a theoretical model is proposed in this work and expects to provide a new available option to investigate oil-water flow behavior in complex pipelines. And the numerical simulation is used to clarify the reliability of the proposed model based on regular scenarios.In this paper, a unified model of oil-water two-phase flow in the complex pipeline is developed based on the combination of pipe serial-parallel theory, flow pattern transformation criterion, two-fluid model, and homogenous model. A case is present to verify the unified model, to compare with CFD results, to analyze the diversion of oil-water mixture under different oil viscosity and water content, and to predict the pressure drop.


## 2. Model Establishment

In this paper, four kinds of oil-water stratified flow patterns including o & w, w/o & w, w/o & o/w, and o & o/w were considered [14, 15], in which two kinds of oil-water dispersed flow patterns such as w/o and o/w were considered.The basic assumptions of the model are shown below: oil and water are Newton incompressible fluid; no heat conduction or working phenomenon between fluid and the outside world, the fluid system is under the isothermal flow state; steady flow; the oil-water interface is assumed to be planar; the interaction between the pipeline and pipe fittings is neglectful; the fluid properties in different channels are consistent.


### 2.1. Flow Regime Transition Criterion

Because the density difference between oil and water is small, and their interfacial tension is low, the oil-water interface is easy to fluctuate and translate into other flow patterns. The transition from stratified flow to dispersed flow depends on the equilibrium relationship between the total turbulent kinetic energy and the total free energy of the dispersed phase. If the turbulent kinetic energy is larger, the dispersed phase will be dispersed in the continuous phase in the form of spherical droplets and collide with each other in the turbulent motion [16, 17]. At the same time, if the droplet size is too large, under the influence of the turbulent force, the droplet will break. Therefore, the maximum dispersion of each flow unit depends on the turbulence intensity of the continuous phase. According to Atmaca et al. [18], when the oil and water mixture velocity is larger than a certain critical value, the fluid will translate the flow pattern from the oil-water two-phase flow to the oil-water two-phase dispersed flow.
(1)vm>6.325CeφIntσowρw−ρog1/2fmρm1/2.Among them,vm represents the speed of the mixture, m/s; Ce represents coefficient, dimensionless; φInt is the internal phase volume fraction, dimensionless; σow is the oil-water surface tension; N/m; ρw represents water density, kg/m3; ρo represents oil density, kg/m3; g is the acceleration of gravity, m/s2; fm represents the friction coefficient between the wall and the mixture, dimensionless; ρm represents the density of liquid mixture; kg/m3.
(2)Ce=2.5−sinθ2,(3)fm=64RemRem≤2320,0.0122Rem−23201680+0.02762320<Rem<4000,0.3164Rem0.254000<Re<105,(4)Rem=ρmvmdμm,(5)ρm=ρc1−φInt+ρdφInt,(6)vm=vc+vd,(7)μm=μc1−φInt−2.5ηc>ηcr,(8)ηcr=11+μo/μw0.4ρo/ρw0.6,where θ represents the angle deviating from the horizontal, degrees; Rem is the mixture Reynolds number, dimensionless; d is wellbore diameter; m; μm is the mixed liquid viscosity, Pa∗s; ρc represents the continuous phase density, kg/m3; ρd represents the dispersion phase density, kg/m3; vc represents the continuous phase velocity, m/s; vd is dispersion phase velocity, m/s; μc is the continuous phase viscosity, Pa s; ηcr represents the critical water content, dimensionless.


### 2.2. Fluid Flow in the Horizontal Pipe

The control body of the oil-water two-phase flow in the horizontal straight pipe is shown in Figure1.Figure 1
Control volume of oil-water two-phase flow in the horizontal pipe.A two-phase flow model in the horizontal straight pipe is established based on the homogeneous flow model, and the momentum balance equation of the control body showed in Figure1 is shown in equation (9).
(9)‐dpdx=4τmD,(10)τm=18fmρmvm2.Among them,τm is the friction stress between mixed fluid and wellbore wall, Pa; D represents the flow area diameter, m. By integrating the formula (9), the general equation of pressure drop of oil and water two-phase flow in the horizontal pipe can be obtained [19].
(11)Δpp=4τmDL.Among them,L is the horizontal straight pipe length, m.However, when the fluid flows from the dispersed flow to the stratified flow, the separation phenomenon will occur, as a result, the above equation is not applicable. The pressure drop equation needs to be derived based on the two-fluid model. The momentum equations of the w/o layer and the o/w layer column are shown in formulas (12) and (13), respectively.
(12)−A1dpdx−τ1S1−τiSi=0,(13)−A2dpdx−τ2S2+τiSi=0.Among them,A1 and A1 represent the pipe cross-sectional area of w/o layer and o/w layer, respectively, m2. P is the average pressure in X section, Pa; τ1 and τ2 are the wall friction stress between w/o layer or o/w layer and wellbore, respectively, Pa; S1 and S2 represent the wet week of w/o layer or o/w layer, m; τi is the interface friction stress between w/o layer and o/w layer under oil-water two-phase stratified flow condition, Pa; Si is the section chord length for w/o layer and o/w layer, m.Formulas (12) and (13) are added and integrated in order to obtain the general pressure drop equation of the oil-water two-phase stratified flow.
(14)Δpp=τ1S1+τ2S2AL.However, several key parameters must be determined for the above model, including the cross-section geometry parameters, the friction stress parameters, and the dispersed phase holdup parameters.The section geometric parameters include the wetted perimeter and the cross-sectional area of the upper and lower two layers. The wetted perimeter of the upper and lower two layers of the oil-water two-phase flow is shown in formulas (15) and (16), and the cross-sectional area is shown in formulas (17) and (18).
(15)S1=πD−S2,(16)S2=γD,(17)A1=D24π−γ+sinγcosγ,(18)A2=D24γ−sinγcosγ.Among them,γ represents the wetter perimeter of the lower o/w layer corresponding to the semicircular, dimensionless.The section chord of w/o layer and o/w layer can be obtained in the same way.(19)Si=Dsinγ.The friction stress is calculated by the concept of friction coefficient.(20)τ1=18f1ρ1v12,(21)τ2=18f2ρ2v22,(22)ρ1=ρo1−αw1+ρwαw1,(23)ρ2=ρoαo2+ρw1−αo2.Among them,f1 and f2 represent the Moody wall friction coefficient of w/o layer and o/w layer, respectively, dimensionless; ρ1 and ρ2 represent the effective density of w/o layer and o/w layer, respectively, kg/m3; v1 and v2 represent the actual flow rate of w/o layer and o/w layer, respectively, m/s; aw1 and ao2 represent the dispersed phase concentration of w/o layer and o/w layer, respectively, dimensionless.The effective Reynolds numbers of the upper and lower two layers of oil-water two-phase flow are shown in formulas (24) and (25), respectively [20, 21].
(24)Re1=ρ1v1d1μ1,(25)Re2=ρ2v2d2μ2.Among them,d1 and d2 represent the hydraulic diameter of w/o layer and o/w layer, respectively, m; μ1 and μ2 represent the effective viscosity of w/o layer and o/w layer, respectively, Pa∗s.According to the velocity of each layer, the hydraulic diameter of w/o layer and o/w layer is different.(26)d1=4A1S1+S2v1>v2,4A1S1v1≤v2,(27)d2=4A2S2+Siv1<v2,4A2S2v1≥v2.Using Roscoe-Brinkman formula to calculate the effective viscosity of w/o and o/w layers.(28)μ1=μo1−αw1−2.5,(29)μ2=μw1−αo2−2.5.


### 2.3. Fluid Flow in the Elbow Pipe

The control body of the oil-water two-phase flow is shown in Figure2.Figure 2
Control volume of oil/water two-phase flow in the elbow.It is similar to the horizontal oil-water two-phase flow that the general pressure drops equation of oil-water two-phase disperse flow and stratified flow in the elbow pipe can be obtained by integrating the pressure gradient equation along the radial direction.(30)Δpef=4τmDRθ,(31)Δpef=τ1S1+τ2S2RAθ,where R is the radius of the elbow pipe, m; θ is the circular measure of the elbow pipe, rad.


### 2.4. Fluid Flow through the Sudden Expansion Joint

The control body of the oil-water two-phase flow in the sudden expansion joint is shown in Figure3. The friction stress and gravity are neglected, and it is assumed that there is no mass exchange between the upper and lower two layers.Figure 3
Control volume of oil-water two-phase flow in the sudden expanded fitting.Based on the homogeneous flow model, a general pressure drop equation for oil-water two-phase flow in the sudden expansion joint is established, as shown in equation (32).
(32)Δpsef=1−AUAD2vU22ρm.Similarly, the separation phenomenon occurs when the dispersed flow is transferred to the stratified flow, and the above equation is no longer applicable. The momentum equation of oil-water two-phase stratified flow in the sudden expansion joint can be derived based on the two-fluid model.(33)pUAU−pDAD+p′AD−AU=ρ1q1vD1−vU1+ρ2q2vD2−vU2,where q1 and q2 represent the rate of flow of w/o layer and o/w layer, respectively.
(34)q1=AU1vU1=AD1vD1,(35)q2=AU2vU2=AD2vD2.According to Che and Li’s research [22], the pressure of the U-U cross-section is equal to the pressure of the smaller cross-section before the sudden expansion joint. Take this equal relation into consideration, the general pressure drop equation of oil-water two-phase flow in the sudden expansion joint is obtained.
(36)Δpsef=ρ1q1vD1−vU1+ρ2q2vD2−vU2AD.


### 2.5. Fluid Flow through the Sudden Shrunken Joint

The control body of the oil-water two-phase flow in the sudden shrunken joint is shown in Figure4. The friction stress and gravity are neglected, and it is assumed that there is no mass exchange between the upper and lower two layers.Figure 4
Control volume of oil-water two-phase flow in the sudden shrunken fitting.When the fluid flows into the U-U section and flows out from the C-C section, the fluid accelerates and the pressure energy is converted to kinetic energy, in this process, the friction dissipation is small and even little. After the fluid flows through the contraction C-C section, the flow state is similar to that of the sudden expansion joint and is accompanied by friction dissipation. Based on the homogeneous flow model, the general pressure drop equation of oil-water two-phase flow in the sudden shrunken joint can be established, as shown in equation (37).
(37)Δpssf=1+1Cv2Cc2−2CcvD22ρm,(38)Cc=ACAD,(39)Cv=v¯cv¯.Among them,Cc represent the contraction coefficient, dimensionless; Cv represent the velocity coefficient, dimensionless; v¯c is the actual average flow velocity in the contraction section, m/s; v¯ is the average velocity in the contraction section, m/s.Based on the two-fluid model, the momentum equation of the oil-water two-phase flow is(40)pCAC−pDAD+p′AD−AC=ρ1q1vD1−vC1+ρ2q2vD2−vC2,where
(41)q1=AU1vU1=AC1vC1=AD1vD1,(42)q2=AU2vU2=AC2vC2=AD2vD2.Similarly, the pressure of the C-C section is equal to the pressure at the outlet whose cross-section is larger than others of the former sudden shrunken joint, then substitute into the formula (40), recognize the formula then obtain the general pressure drop equation for oil-water two-phase stratified flow.
(43)Δpssf=ρ1q1vD1−vC1+ρ2q2vD2−vC2AD.


### 2.6. Serial-Parallel Pipe Theory

There are two important rules when the straight pipe is connected with various pipe fittings in series [23–25]:
(44)Q=Q=1Q2=⋯=Qn,(45)ΔP=ΔP1+ΔP2+⋯+ΔPn.These two rules change when they are connected in parallel:(46)Q=Q+1Q2+⋯+Qn,(47)ΔP=ΔP1=ΔP2=⋯=ΔPn.Therefore, by means of simultaneous (44)–(47), the split-flow and pressure drop of oil-water two-phase flow can be obtained.


## 2.1. Flow Regime Transition Criterion

Because the density difference between oil and water is small, and their interfacial tension is low, the oil-water interface is easy to fluctuate and translate into other flow patterns. The transition from stratified flow to dispersed flow depends on the equilibrium relationship between the total turbulent kinetic energy and the total free energy of the dispersed phase. If the turbulent kinetic energy is larger, the dispersed phase will be dispersed in the continuous phase in the form of spherical droplets and collide with each other in the turbulent motion [16, 17]. At the same time, if the droplet size is too large, under the influence of the turbulent force, the droplet will break. Therefore, the maximum dispersion of each flow unit depends on the turbulence intensity of the continuous phase. According to Atmaca et al. [18], when the oil and water mixture velocity is larger than a certain critical value, the fluid will translate the flow pattern from the oil-water two-phase flow to the oil-water two-phase dispersed flow.
(1)vm>6.325CeφIntσowρw−ρog1/2fmρm1/2.Among them,vm represents the speed of the mixture, m/s; Ce represents coefficient, dimensionless; φInt is the internal phase volume fraction, dimensionless; σow is the oil-water surface tension; N/m; ρw represents water density, kg/m3; ρo represents oil density, kg/m3; g is the acceleration of gravity, m/s2; fm represents the friction coefficient between the wall and the mixture, dimensionless; ρm represents the density of liquid mixture; kg/m3.
(2)Ce=2.5−sinθ2,(3)fm=64RemRem≤2320,0.0122Rem−23201680+0.02762320<Rem<4000,0.3164Rem0.254000<Re<105,(4)Rem=ρmvmdμm,(5)ρm=ρc1−φInt+ρdφInt,(6)vm=vc+vd,(7)μm=μc1−φInt−2.5ηc>ηcr,(8)ηcr=11+μo/μw0.4ρo/ρw0.6,where θ represents the angle deviating from the horizontal, degrees; Rem is the mixture Reynolds number, dimensionless; d is wellbore diameter; m; μm is the mixed liquid viscosity, Pa∗s; ρc represents the continuous phase density, kg/m3; ρd represents the dispersion phase density, kg/m3; vc represents the continuous phase velocity, m/s; vd is dispersion phase velocity, m/s; μc is the continuous phase viscosity, Pa s; ηcr represents the critical water content, dimensionless.


## 2.2. Fluid Flow in the Horizontal Pipe

The control body of the oil-water two-phase flow in the horizontal straight pipe is shown in Figure1.Figure 1
Control volume of oil-water two-phase flow in the horizontal pipe.A two-phase flow model in the horizontal straight pipe is established based on the homogeneous flow model, and the momentum balance equation of the control body showed in Figure1 is shown in equation (9).
(9)‐dpdx=4τmD,(10)τm=18fmρmvm2.Among them,τm is the friction stress between mixed fluid and wellbore wall, Pa; D represents the flow area diameter, m. By integrating the formula (9), the general equation of pressure drop of oil and water two-phase flow in the horizontal pipe can be obtained [19].
(11)Δpp=4τmDL.Among them,L is the horizontal straight pipe length, m.However, when the fluid flows from the dispersed flow to the stratified flow, the separation phenomenon will occur, as a result, the above equation is not applicable. The pressure drop equation needs to be derived based on the two-fluid model. The momentum equations of the w/o layer and the o/w layer column are shown in formulas (12) and (13), respectively.
(12)−A1dpdx−τ1S1−τiSi=0,(13)−A2dpdx−τ2S2+τiSi=0.Among them,A1 and A1 represent the pipe cross-sectional area of w/o layer and o/w layer, respectively, m2. P is the average pressure in X section, Pa; τ1 and τ2 are the wall friction stress between w/o layer or o/w layer and wellbore, respectively, Pa; S1 and S2 represent the wet week of w/o layer or o/w layer, m; τi is the interface friction stress between w/o layer and o/w layer under oil-water two-phase stratified flow condition, Pa; Si is the section chord length for w/o layer and o/w layer, m.Formulas (12) and (13) are added and integrated in order to obtain the general pressure drop equation of the oil-water two-phase stratified flow.
(14)Δpp=τ1S1+τ2S2AL.However, several key parameters must be determined for the above model, including the cross-section geometry parameters, the friction stress parameters, and the dispersed phase holdup parameters.The section geometric parameters include the wetted perimeter and the cross-sectional area of the upper and lower two layers. The wetted perimeter of the upper and lower two layers of the oil-water two-phase flow is shown in formulas (15) and (16), and the cross-sectional area is shown in formulas (17) and (18).
(15)S1=πD−S2,(16)S2=γD,(17)A1=D24π−γ+sinγcosγ,(18)A2=D24γ−sinγcosγ.Among them,γ represents the wetter perimeter of the lower o/w layer corresponding to the semicircular, dimensionless.The section chord of w/o layer and o/w layer can be obtained in the same way.(19)Si=Dsinγ.The friction stress is calculated by the concept of friction coefficient.(20)τ1=18f1ρ1v12,(21)τ2=18f2ρ2v22,(22)ρ1=ρo1−αw1+ρwαw1,(23)ρ2=ρoαo2+ρw1−αo2.Among them,f1 and f2 represent the Moody wall friction coefficient of w/o layer and o/w layer, respectively, dimensionless; ρ1 and ρ2 represent the effective density of w/o layer and o/w layer, respectively, kg/m3; v1 and v2 represent the actual flow rate of w/o layer and o/w layer, respectively, m/s; aw1 and ao2 represent the dispersed phase concentration of w/o layer and o/w layer, respectively, dimensionless.The effective Reynolds numbers of the upper and lower two layers of oil-water two-phase flow are shown in formulas (24) and (25), respectively [20, 21].
(24)Re1=ρ1v1d1μ1,(25)Re2=ρ2v2d2μ2.Among them,d1 and d2 represent the hydraulic diameter of w/o layer and o/w layer, respectively, m; μ1 and μ2 represent the effective viscosity of w/o layer and o/w layer, respectively, Pa∗s.According to the velocity of each layer, the hydraulic diameter of w/o layer and o/w layer is different.(26)d1=4A1S1+S2v1>v2,4A1S1v1≤v2,(27)d2=4A2S2+Siv1<v2,4A2S2v1≥v2.Using Roscoe-Brinkman formula to calculate the effective viscosity of w/o and o/w layers.(28)μ1=μo1−αw1−2.5,(29)μ2=μw1−αo2−2.5.


## 2.3. Fluid Flow in the Elbow Pipe

The control body of the oil-water two-phase flow is shown in Figure2.Figure 2
Control volume of oil/water two-phase flow in the elbow.It is similar to the horizontal oil-water two-phase flow that the general pressure drops equation of oil-water two-phase disperse flow and stratified flow in the elbow pipe can be obtained by integrating the pressure gradient equation along the radial direction.(30)Δpef=4τmDRθ,(31)Δpef=τ1S1+τ2S2RAθ,where R is the radius of the elbow pipe, m; θ is the circular measure of the elbow pipe, rad.


## 2.4. Fluid Flow through the Sudden Expansion Joint

The control body of the oil-water two-phase flow in the sudden expansion joint is shown in Figure3. The friction stress and gravity are neglected, and it is assumed that there is no mass exchange between the upper and lower two layers.Figure 3
Control volume of oil-water two-phase flow in the sudden expanded fitting.Based on the homogeneous flow model, a general pressure drop equation for oil-water two-phase flow in the sudden expansion joint is established, as shown in equation (32).
(32)Δpsef=1−AUAD2vU22ρm.Similarly, the separation phenomenon occurs when the dispersed flow is transferred to the stratified flow, and the above equation is no longer applicable. The momentum equation of oil-water two-phase stratified flow in the sudden expansion joint can be derived based on the two-fluid model.(33)pUAU−pDAD+p′AD−AU=ρ1q1vD1−vU1+ρ2q2vD2−vU2,where q1 and q2 represent the rate of flow of w/o layer and o/w layer, respectively.
(34)q1=AU1vU1=AD1vD1,(35)q2=AU2vU2=AD2vD2.According to Che and Li’s research [22], the pressure of the U-U cross-section is equal to the pressure of the smaller cross-section before the sudden expansion joint. Take this equal relation into consideration, the general pressure drop equation of oil-water two-phase flow in the sudden expansion joint is obtained.
(36)Δpsef=ρ1q1vD1−vU1+ρ2q2vD2−vU2AD.


## 2.5. Fluid Flow through the Sudden Shrunken Joint

The control body of the oil-water two-phase flow in the sudden shrunken joint is shown in Figure4. The friction stress and gravity are neglected, and it is assumed that there is no mass exchange between the upper and lower two layers.Figure 4
Control volume of oil-water two-phase flow in the sudden shrunken fitting.When the fluid flows into the U-U section and flows out from the C-C section, the fluid accelerates and the pressure energy is converted to kinetic energy, in this process, the friction dissipation is small and even little. After the fluid flows through the contraction C-C section, the flow state is similar to that of the sudden expansion joint and is accompanied by friction dissipation. Based on the homogeneous flow model, the general pressure drop equation of oil-water two-phase flow in the sudden shrunken joint can be established, as shown in equation (37).
(37)Δpssf=1+1Cv2Cc2−2CcvD22ρm,(38)Cc=ACAD,(39)Cv=v¯cv¯.Among them,Cc represent the contraction coefficient, dimensionless; Cv represent the velocity coefficient, dimensionless; v¯c is the actual average flow velocity in the contraction section, m/s; v¯ is the average velocity in the contraction section, m/s.Based on the two-fluid model, the momentum equation of the oil-water two-phase flow is(40)pCAC−pDAD+p′AD−AC=ρ1q1vD1−vC1+ρ2q2vD2−vC2,where
(41)q1=AU1vU1=AC1vC1=AD1vD1,(42)q2=AU2vU2=AC2vC2=AD2vD2.Similarly, the pressure of the C-C section is equal to the pressure at the outlet whose cross-section is larger than others of the former sudden shrunken joint, then substitute into the formula (40), recognize the formula then obtain the general pressure drop equation for oil-water two-phase stratified flow.
(43)Δpssf=ρ1q1vD1−vC1+ρ2q2vD2−vC2AD.


## 2.6. Serial-Parallel Pipe Theory

There are two important rules when the straight pipe is connected with various pipe fittings in series [23–25]:
(44)Q=Q=1Q2=⋯=Qn,(45)ΔP=ΔP1+ΔP2+⋯+ΔPn.These two rules change when they are connected in parallel:(46)Q=Q+1Q2+⋯+Qn,(47)ΔP=ΔP1=ΔP2=⋯=ΔPn.Therefore, by means of simultaneous (44)–(47), the split-flow and pressure drop of oil-water two-phase flow can be obtained.


## 3. Case Study

Figure5 is the schematic diagram of the complex pipeline, which contains 2 parallel channels, each of which is made up of a series of pipes and fittings, which includes the common fittings and connections in the petroleum industry. The model is verified by the example, and the results are compared with the numerical simulation results.Figure 5
Diagram of the complex pipeline.The structure parameters of the complex pipeline are shown in Table1, which is generated in the professional CAD software SOLIDWORKS, transferred to the GAMBIT software in the FLUENT through the CAD interface, then through Boolean operation get the internal flow model and mesh division. When using the FLUENT software to calculate, the inlet is set as the velocity inlet (velocity-inlet) condition, the exit is set as the outflow (outflow) condition, and the other is the default (wall) condition. The laminar flow and turbulent flow are simulated by the laminar flow model and the standard κ−ε model, respectively. The mixed model and VOF model are used to simulate the oil-water two-phase dispersion and stratified flow, respectively. And the inlet flow rate is set as 30 m3/D, and the velocity streamline diagram of oil phase and water phase through the complex pipeline can be obtained, as shown in Figure 6.Table 1
Structure parameters for the present complex pipeline.
Structure parametersValueUnitChannel 1Sudden shrunken joints number1DimensionlessSudden expansion joints number1DimensionlessExpansion channel length0.1633mExpansion channel radius0.005mShrunken channel length0.1633mShrunken channel radius0.0025mChannel 2Elbow pipe number12DimensionlessBend angle90°Bending radius0.03mFigure 6
Path line colored by velocity magnitude.According to Figure6, it can be seen that the oil phase and water phase through the complex pipe after shunting vary greatly, this is due to the fact that the pressure drop compositions of the two parallel lines are different, line 1 is mainly of local resistance and has large water resistance; the main pressure drop of line 2 is friction, so its resistance to the oil phase is greater. At the same time, as the great property difference of oil-water mixture in different oil phase viscosity and water content, in order to accurately describe the influence of water cut and oil phase viscosity on distributary situation and pressure drop, moisture content (%) values are as follows: 0, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100; the values of the oil-phase viscosity (cP) are as follows: 1, 10, 100, and other parameters are shown in Table 2.Table 2
Fluid parameters for oil-water two-phase flow in the complex pipeline.
Fluid parametersValueUnitRate of flow30m3/DWater content0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100%Oil phaseDensity850kg/m3Viscosity1, 10, 100mPa·sWater phaseDensity998.2kg/m3Viscosity1.003mPa·sSubstituted the above parameters into the model and the numerical simulation, then the distributary situation and the pressure drop of oil-water two-phase flow through the complex pipeline can be obtained, as shown in Figures7 and 8.Figure 7
Flow distribution through the complex pipeline with varying fluid properties.Figure 8
Pressure drop through the complex pipeline with varying fluid properties.According to Figures7 and 8, the results show that the model is in good agreement with the numerical simulation results. The highest average absolute error is about 14.4%, and the average error is 9.8% when the oil-water two-phase volumetric water content is 0% and 100%, respectively. And these curves have 3 major turning point: the first is reverse point, if the oil viscosity increases, the critical volume water content will decrease, then the position of reverse point will move forward; the second one and the 3rd are the turning point making the fluid flow from laminar transition to turbulent transition zone and from turbulent transition zone to the turbulent zone. In addition, if the flow rate is greater than a certain critical value, the flow pattern will also change, which will affect the flow distribution. In addition, if the flow rate is greater than a certain critical value, the flow pattern will also change, which will affect the flow distribution.Therefore, if the mixed liquid water content is larger than the critical moisture content, because the water phase is dominant, the distributary and the pressure drop are basically same under different oil-phase viscosity; in other cases, as the oil phase dominates, the distributary and the pressure drop vary greatly. With the increase of the oil-phase viscosity, the diversion effect will be better. With the increase of water content, the diversion effect will become good at first then go bad, and it will reach the best result at the critical moisture content. Correspondingly, the pressure drops of the oil-water two-phase fluid passing through the device increases with the oil phase-viscosity; with the increasing of water cut, the pressure drops become good at first then go bad and reach the best result at the critical moisture content.In addition, the close relationship between oil/water two phase flow inside pipelines and huge economic interests from the development of conventional or unconventional geological hydrocarbon should be highlighted, which is an indispensable branch of economic geology. Notably, the inherent scientific question required to be addressed regarding the production of geological hydrocarbon is the multiphase flow inside pipeline. The accurate prediction of flux or pressure field not only relies on the theoretical multiphase flow model but also heavily dependent on surrounding geological condition. In the other words, the research content lays the solid theoretical background for efficient production and tightly ties with economical geology. Also, from the academic perspective, the research content can be regarded as the cross-region between mechanical engineering as well as petroleum engineering. Meanwhile, based on the application motivation, the majority of potential readership will come from petroleum engineering. As a result, the paper content is demonstrated to be closely linked to economic geology, rather than mechanical engineering.


## 4. Conclusions

Based on the theory of series-parallel pipeline, the flow regime transition criterion, two-fluid model, and homogeneous flow model, the unified mechanism model of oil-water two-phase flow through the complex pipeline is established, the model was validated through an example, then compared with numerical simulations, and the main conclusions are as follows:(1)
The model has a performance in predicting the distributary situation and pressure drop of the oil-water two-phase mixture through complex pipeline, the highest average absolute error is about 14.4%, and the average error is 9.8% when the oil-water two-phase volumetric water content is 0% and 100%, respectively(2)
The shunt situation of the oil-water two-phase mixture through the complex pipeline will have 3 transitions, due to the transition of antiphase phenomena, the transition from laminar flow zone to the turbulent transition zone, and the transition from the transition zone to the turbulent flow zone. With the increase of water content, the diversion effect will become good at first then go bad, and it will reach the best result at the critical moisture content(3)
Corresponding to the diversion situation, the pressure drop of oil-water mixture fluid through the example device will also have 3 transitions, and with the increasing of water cut, the pressure drop become good at first then go bad and reach the best result at the critical moisture content



---

*Source: 3756577-2021-07-22.xml*