^{1}

^{*}

^{2}

^{1}

The development of shale gas reservoir is mainly based on horizontal well production. Slug flow of gas-liquid two-phase is invariably encountered in inclined wells and horizontal wells of a producing environment. Due to gravitational differentiation, oil-water two-phase flow pattern, the local velocity and local phase holdup along the radial direction of pipe in near horizontal wells will perform complicatedly. This paper presented the results of an experimental study and a theoretical analysis of two-phase gas/water flow in horizontal and highly inclined systems. Extensive experiments were conducted using a test loop made of 124 mm diameter acrylic pipe with inclination angles from the horizontal of 0°, 5°, 15°, 45°, -2°, -5° and -10°, and with the total flow rate ranging from 50 to 800 m
^{3}/day. Based on the research on the law of slug flow dynamics model for gas-water two-phase flow in near horizontal pipeline, the theoretical analysis and experimental researches were done to propose the expressions of stable and exact production logging interpretation model for two-phase flow in near horizontal pipeline. The performance of the proposed method for estimating water holdup and water superficial velocity is in good agreement with our measurements. As a result, the slug flow dynamics model of gas-water two-phase flow in near horizontal wellbore was developed. The application effect of production logging in near horizontal wells had been improved.

As an unconventional energy source, the shale gas plays an important role in the global energy structure. The horizontal well mining method has obvious advantages in the development of shale gas reservoirs. Horizontal well production in shale gas reservoir is a mining technology developed after N2 fracturing, foam fracturing, gel fracturing, clear water fracturing and so on [

Accurate horizontal well production logging and interpretation evaluation technology undoubtedly provide a reliable basis for efficient development. The demand is also growing for well logging at horizontal and large inclined wells. Production profile logging is the main technical means for monitoring the production dynamics of production wells and is one of the important supporting technologies for the development of horizontal wells. The characteristics of gas-liquid two-phase flow in near horizontal wells are important for the establishment of the interpretation models of production profile logging data. And production wells are in the best state of production and the oil recovery is improved.

Slug flow is a common flow pattern in the horizontal or slightly inclined production oil-gas wells and oil-gas pipeline, which is the alternation of the liquid plug body and the long air bag in space and time, and that shows intermittent and unstable characteristics during the flow process [

Slug flow and stratified flow often occur in the normal operation of oil and gas well production and oil and gas mixed transportation, and the slug flow is more frequent. Slug flow is almost always inevitable in long-distance pipeline. Actually, many researches on multi-phase flow of oil and gas mixed transportation are to study slug flow characteristics [

In this paper, a dynamic model of gas-liquid slug flow was established on the basis of the dynamic characteristics of gas-liquid two phase flow. By using a ground multiphase flow loop simulation experiment device, the internal relationship between the liquid holdup and the two-phase characteristic parameters of the slug flow in the horizontal and large inclined tubes was studied on the basis of the simulated well gas water two phase dynamic experiment to get a more perfect and more comprehensive law. The law can better guide the logging interpretation of horizontal and large inclined well production profile, and improve the production logging application technology.

For the convenience of analysis, it takes a section of the continuous flow tube with typical characteristics as a slug body unit. The typical model characteristics of a single slug unit are shown in

For a slug body unit in which the slug body is stable, its translational velocity is equal to the slug head velocity and airbag velocity v_{T}. The average liquid mass flow rate over the time of the passage of a slug unit is [

W L = 1 T U ( v LS A H LS ρ L T S + ∫ 0 T F v LF A H LF ρ L d t ) (1)

where W_{L} is the input (or average) mass flow rate, kg/s; T_{U}, T_{S}, and T_{F} are the times for the passage of the slug unit, the liquid slug, and the liquid film/gas pocket, s; A is the pipe cross-sectional area, m^{2}; ρ_{L} is liquid phase density, kg/m^{3}.

T U = L U v T , T S = L S v T , T F = L F v T (2)

If the liquid film height is stable, the superficial velocity of the liquid phase and gas phase of the slug body unit are respectively as follows.

v SL = v LS H LS L S L U + L F L U v LF H LF (3)

v SG = v GS ( 1 − H LS ) L S L U + L F L U v GF ( 1 − H LF ) (4)

where v_{SL} is the superficial velocity of the liquid phase; v_{SG} is the superficial velocity of the gas phase.

Because of that T_{U} = T_{S} + T_{F}, the length of the whole slug element can be expressed as follows.

Then the total length of the slug unit can be derived from Equation (1):

L U = L S ( v LS H LS − v LF H LF v SL − v LF H LF ) = L S [ v GS ( 1 − H LS ) − v GF ( 1 − H LF ) v SG − v GF ( 1 − H LF ) ] (5)

For the liquid slug body, the head of the slug moves forward at a velocity v_{T}. And it wraps the slow flow of fluid in the liquid film at a volume flow Q_{in} (v_{T} > v_{LF}). At the tail of the liquid slug, the liquid in the liquid film moves at a velocity of v_{LS}. Slug tail moves at airbag velocity v_{T}. As v_{T} is greater than v_{LS}, the amount of the remaining liquid slug tail is Q_{out}. With v_{T} as the reference coordinate, the followings are the liquid volume flow into the slug head and the amount of the remaining liquid of the slug tail respectively [

Q in = ( v T − v LF ) H LF A (6)

Q out = ( v T − v LS ) H LS A (7)

When the liquid slug is stable, according to the law of conservation of mass, the mass flow rate of liquid into and out of the liquid slug is equal for the whole slug body unit.

ρ L Q in = ρ L Q out (8)

( v T − v LF ) H LF = ( v T − v LS ) H LS (9)

That is, the liquid mass flow in liquid slug and liquid film is equal.

Similarly, the mass flow rate of gas phase in liquid slug and air bag is also equal.

( v T − v GF ) ( 1 − H LF ) = ( v T − v GS ) ( 1 − H LS ) (10)

For the whole slug unit, its average liquid holdup H_{LU} is defined by

H L U = H L S L S + H L F L F L U (11)

From Equations (3), (4), (9) and (11), the average liquid holdup of the whole slug unit can be expressed by

H LU = H LS v T + ( 1 − H LS ) v GS − v SG v TB = H LF v T + ( 1 − H LF ) v GF − v SG v TB = H LF v T + v SL − H LF v LF − v SG v T (12)

If the height of the liquid film is stable, the liquid film and the airbag are similar to the stratified flow in the slug body unit, and the momentum equation is

τ WF S L A L − τ WG S G A G − τ I S I ( 1 A L + 1 A G ) + ( ρ L − ρ G ) g sin θ = 0 (13)

where A_{L} is cross-sectional area available to the liquid phase, m^{2}; A_{G} is cross-sectional area available to the gas phase, m^{2}; S_{L} is wetted periphery for the liquid phase, m; S_{G} is no-wetted periphery for the liquid phase, m; S_{I} is wetted periphery for the interface, m; g is gravitational acceleration, m/s^{2}; ρ_{G} is gas phase density, kg/m^{3}.

The momentum equation is an implicit equation of liquid holdup rate H_{LF} and gas/liquid superficial velocity v_{LF}, v_{GF}.

Evaluation of the correlations and models was based on the statistical parameters in the above equations. In the horizontal and vertical pipe flow, the length of the slug in the stable slug body unit is respectively as follows [

L S = 30 D ( θ = 0 ) (14)

L S = 20 D ( θ = 90 ) (15)

where D is the inner diameter of the pipe.

In near-level large-diameter pipe flow (θ = ±1˚, D > 2in ), the formula given by Scott et al. [

ln ( L S ) = − 25.4 + 28.5 [ ln ( D ) ] 0.1 (16)

In the inclined pipe section, the formula given by Zhang Hongquan et al. [

L S = 30 D cos 2 θ + 20 D sin 2 θ (17)

The prediction formula of liquid holdup in liquid slug is given by Gomez et al. [

H LS = 1.0 e − ( 7.85 × 10 − 3 θ + 2.48 × 10 − 6 R eSL ) , ( 0 ≤ θ ≤ 90 ∘ ) (18)

where R_{eSL} is the apparent Reynolds number of the slug body.

The apparent Reynolds number of the slug body in the formula is as follows [

R eSL = D v M ρ L μ L (19)

v T = C 0 v M + ( 0.542 g D cos θ + 0.351 g D sin θ ) (20)

where v_{M} is the average flow rate of the fluid [

The gas phase velocity in the liquid slug body is

v G S = C 1 v M + v 0 ∞ sin θ H L S 0 . 5 (21)

where C_{1} is the velocity distribution coefficient and Chokshi et al. [_{1} is equal to 1.15. v 0 ∞ is the bubble rising velocity. Harmathy [

v 0 ∞ = 1.53 ( g σ ( ρ L − ρ G ) ρ L 2 ) 0.25 (22)

where σ is the viscosity of the liquid.

The experiment was completed at the horizontal well and high inclination angle flow simulation test facility of Yangtze University. The two simulation wellbore holes of the test flow loop are 16 m organic transparent glasses with inner diameter of 159 mm and 124 mm respectively. The mobile loop can be set at any angle between the horizontal and the vertical. The experimental conditions are air and tap water. The experimental conditions are 16˚C and 1 atm, and the density of gas and water is 0.0012 g/cm^{3} and 0.9992 g/cm^{3} respectively. 5 flows were used from low to high in the whole experiment. v_{M} was 50, 100, 200, 400 and 800 m^{3}/d respectively, the inclination angle of corresponding horizontal direction was 0˚, 5˚, 15˚, 45˚, −2˚, −5˚, −10˚ respectively, and the moisture content was 0%, 10%, 30%, 50%, 70% and 90% respectively. The instantaneous shut well holdup was measured under each experiment. The fast cutting valve was installed in the two sections of the simulated wellbore. The closing time interval of the two valve was less than 0.5 s. After that, the vertical wellbore was simulated and the liquid height was measured. Thus the water holdup in the shut well was calculated, and the measured value was taken as the standard holdup.

In order to simulate the mixed flow of gas and water in the near horizontal well, the smooth stratified flow (SF), the wave stratified flow (SWF), the bubble flow (BF), the slug flow (SLF) and the annular flow (AF) appeared in the 124 mm transparent wellbore under the above experimental conditions. The flow pattern of three different parameters of well deviation angle, different water holdup and total flow rate were obtained by analyzing the experimental data, as shown in

According to the flow pattern analysis observed by the experiment, the inclined angles of the slug flow were +5˚, +15˚ and +45˚ under the above experimental conditions. The slug flow was not observed under the condition of horizontal and inclined downward angle.

The experimental gas and water apparent velocity of the slug flow was taken as the true theoretical gas and water velocity. And the water holdup was substituted by the theoretical model of liquid holdup (Equation (12)) as a prediction variable. It is calculated that the relationship between the calculated value of water holdup and the actual value of the slug flow angle of +5˚, +15˚ and +45˚ under various experimental conditions. As shown in Figures 5-7, the theoretical water holdup is in good agreement with the experimental results.

By using Equation (12) and the apparent velocity of gas and water, the comparison between the predicted and actual values of the slug flow theoretical model and the error analysis can be gotten as shown in

According to the analysis of the experimental and theoretical calculation results above, it can be seen that the water holdup predicted by the theoretical model of gas and water two-phase flow model corresponding to the flow pattern structure is in good agreement with the experimental results. The higher the total flow rate is, the higher the coincidence rate is. The average absolute error of predicted water holding capacity and actual water holding capacity is 0.048128, and the average relative error is 6.8313%. The theoretical model corresponding to the slug flow can better reflect the flow characteristics of the multiphase flow, and the prediction error of water holdup and experimental water holdup is within the allowable error range of production log interpretation.

In the same way, the total gas-water flow rate of the slug flow was used as the true theoretical air-water total flow rate, and the actual water hold-up rate was measured as the true theoretical water holdup. Air-water superficial velocity was used as a predictor to substitute for the momentum conservation model of Equation (13). The theoretical calculations of gas and water apparent velocities under the experimental conditions of 5˚, 15˚, and 4˚ gas-water slug flow were calculated and compared with the experimental values, as shown in Figures 9-12.

According to the results of the comparison, it can be seen that the apparent velocity of the water phase predicted by the theoretical model (Equation (13)) of the gas-water two-phase model based on the slug flow pattern structure is close to the experimental value. The average absolute error is 0.01217 m/s, and the average relative error is 13.8551%. The theoretical and experimental results are in good agreement.

Through the production logging interpretation theoretical model based on slug flow structure, the following conclusions were obtained by comparing the predicted results of water holding capacity and water surface velocity with the experimental values in the inclined well.

1) In this paper, the theoretical distribution models for production logging interpretation of slug flow were substituted into the theoretical values of experimental gas, the water apparent velocity and the water holdup as known values for positive and negative inversion calculations. The calculated results are in good agreement with the experimental results, demonstrating that the model is stable.

2) For the gas-water two-phase flow, the theoretical model for production log interpretation based on the flow pattern structure can better reflect the flow characteristics of the gas-water two-phase flow in oil wells, and theoretically predict the water holding capacity and water phase. The apparent velocity is close to the experimental result, and the error value is within the allowable range of the production log interpretation output industry. The theoretical model can be used to explain the production gas and water two-phase output profile.

3) The total flow rate and the phase flow rate calculated by participating in the theoretical model are experimental rations, and the water holdup is the instantaneous water retention rate measured by the shut-in well. There is a certain error between the experimental ration value and the well logging instrument response value. When it is applied to multiphase flow interpretation of horizontal wells and inclined well production logging, it is also necessary to make certain corrections to the logging instrument response values, and the calculation accuracy needs to be further verified.

The research is supported by Educational Commission of Hubei Province of China (D20141302), Open Fund of Key Laboratory of Exploration Technologies for Oil and Gas Resources (Yangtze University), Ministry of Education (No. K2018-02) and National Natural Science Foundation of China (41474115).

The authors declare no conflicts of interest regarding the publication of this paper.

Song, H.W., Guo, H.M. and Xu, S.H. (2019) Analysis of Gas and Liquid Two-Phase Slug Flow Production Logging Interpretation Model in near Horizontal Shale Gas Wells. Open Journal of Yangtze Gas and Oil, 4, 100-112. https://doi.org/10.4236/ojogas.2019.42008