Flow Regime Identification of Gas-liquid Two-phase Flow Based on HHT

A new method to identify flow regime in two-phase flow was presented, based on signal processing of differential pressure using Hilbert Huang transform （HHT）. Signals obtained from a Venturi meter were decomposed into different intrinsic mode functions （IMFs） with HHT, then the energy fraction of each intrinsic mode and the mean value of residual function were calculated, from which the rules of flow regime identification were summarized. Experiments were carried out on two-phase flow in the horizontal tubes with 50mm and 40mm inner diameter, while water flowrate was in the range of 1.3m^3.h^-1 to 10.5m^3.h^-1, oil flowrate was from 4.2m^3.h^-1 to 7.0m^3.h^-1 and gas flowrate from 0 to 15m^3.h^-1. The results show that the proposed rules have high precision for single phase, bubbly, and slug, plug flow regirne identification, which are independent of not only properties of two-phase fluid. In addition, the method can meet the need of industrial application because of its simple calculation.

Identification Method of Gas-Liquid Two-phase Flow Regime Based on Image Multi-feature Fusion and Support Vector Machine

The knowledge of flow regime is very important for quantifying the pressure drop, the stability and safety of two-phase flow systems. Based on image multi-feature fusion and support vector machine, a new method to identify flow regime in two-phase flow was presented. Firstly, gas-liquid two-phase flow images including bub- bly flow, plug flow, slug flow, stratified flow, wavy flow, annular flow and mist flow were captured by digital high speed video systems in the horizontal tube. The image moment invariants and gray level co-occurrence matrix texture features were extracted using image processing techniques. To improve the performance of a multiple classifier system, the rough sets theory was used for reducing the inessential factors. Furthermore, the support vector machine was trained by using these eigenvectors to reduce the dimension as flow regime samples, and the flow regime intelligent identification was realized. The test results showed that image features which were reduced with the rough sets theory could excellently reflect the difference between seven typical flow regimes, and successful training the support vector machine could quickly and accurately identify seven typical flow regimes of gas-liquid two-phase flow in the horizontal tube. Image multi-feature fusion method provided a new way to identify the gas-liquid two-phase flow, and achieved higher identification ability than that of single characteristic. The overall identification accuracy was 100%, and an estimate of the image processing time was 8 ms for online flow regime identification.

Multi-scale Chaotic Analysis of the Characteristics of Gas-Liquid Two-phase Flow Patterns

Using the high-speed camera the time sequences of the classical flow patterns of horizontal gas-liquid pipe flow are recorded, from which the average gray-scale values of single-frame images are extracted. Thus obtained gray-scale time series is decomposed by the Empirical Mode Decomposition (EMD) method, the various scales of the signals are processed by Hurst exponent method, and then the dual-fractal characteristics are obtained. The scattered bubble and the bubble cluster theories are applied to the evolution analysis of two-phase flow patterns. At the same time the various signals are checked in the chaotic recursion chart by which the two typical characteristics (diagonal average length and Shannon entropy) are obtained. Resulting term of these properties, the dynamic characteristics of gas-liquid two-phase flow patterns are quantitatively analyzed. The results show that the evolution paths of gas-liquid two-phase flow patterns can be well characterized by the integrated analysis on the basis of the gray-scale time series of flowing images from EMD, Hurst exponents and Recurrence Plot (RP). In the middle frequency section (2nd, 3rd, 4th scales), three flow patterns decomposed by the EMD exhibit dual fractal characteristics which represent the dynamic features of bubble cluster, single bubble, slug bubble and scattered bubble. According to the change of diagonal average lengths and recursive Shannon entropy characteristic value, the structure deterministic of the slug flow is better than the other two patterns. After the decomposition by EMD the slug flow and the mist flow in the high frequency section have obvious peaks. Anyway, it is an effective way to understand and characterize the dynamic characteristics of two-phase flow patterns using the multi-scale non-linear analysis method based on image gray-scale fluctuation signals.

Reconstructing bubble profiles from gas-liquid two-phase flow data using agglomerative hierarchical clustering method

The knowledge of bubble profiles in gas-liquid two-phase flows is crucial for analyzing the kinetic processes such as heat and mass transfer, and this knowledge is contained in field data obtained by surface-resolved computational fluid dynamics (CFD) simulations. To obtain this information, an efficient bubble profile reconstruction method based on an improved agglomerative hierarchical clustering (AHC) algorithm is proposed in this paper. The reconstruction method is featured by the implementations of a binary space division preprocessing, which aims to reduce the computational complexity, an adaptive linkage criterion, which guarantees the applicability of the AHC algorithm when dealing with datasets involving either non-uniform or distorted grids, and a stepwise execution strategy, which enables the separation of attached bubbles. To illustrate and verify this method, it was applied to dealing with 3 datasets, 2 of them with pre-specified spherical bubbles and the other obtained by a surface-resolved CFD simulation. Application results indicate that the proposed method is effective even when the data include some non-uniform and distortion.

An Experimental Study on the Void Fraction for Gas-Liquid Two-Phase Flows in a Horizontal Pipe

The flow patterns and the void fraction related to a gas-liquid two-phase flow in a small channel are experimentally studied.The test channel is a transparent quartz glass circular channel with an inner diameter of 6.68 mm.The working fluids are air and water and their superficial velocities range from 0.014 to 8.127 m/s and from 0.0238 to 0.556 m/s,respectively.The void fraction is determined using the flow pattern images captured by a high-speed camera,while quick closing valves are used for verification.Four flow patterns are analyzed in experiments:slug flow,bubbly flow,annular flow and stratified flow.For intermittent flows(bubbly flow and slug flow),the cross-sectional void fraction is in a borderline condition while its probability distribution function(PDF)image displays a bimodal structure.For continuous flows(annular flow and stratified flow)the cross-sectional void fraction behaves as a fluctuating continuous curve while the(PDF)image displays a single peak structure.The volumetric void fraction data are also compared with available predictive formulas,and the results show that the agreement is very good.An effort is also provided to improve the so-called Gregory and Scott model using the available data.

New numerical algorithm of gas-liquid two-phase flow considering characteristics of liquid metal during mold filling

A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Stokes equations and solved through Projection method. The Level set method is used to track the gas-liquid interface boundary. In order to demonstrate the correctness of this new program for simulation of gas-liquid two-phase mold filling in casting, a benchmark filling experiment is simulated （this benchmark test is designed by XU and the filling process is recorded by a 16-mm film camera）. The simulated results agree very well with the experimental results, showing that this new program can be used to properly predicate the gas-liquid two-phase mold filling simulation in casting.

Branch Quality Control of Gas-Liquid Two-Phase Flow Using a Novel T-Junction Type Distributor

In order to eliminate mal-distribution and ensure the side arm to produce desirable gas quality a special distributor is proposed. The experimental distributor mainly consists of a straight through section,a gas extraction line,a liquid extraction line and a side arm branch. A gas orifice and a liquid orifice are mounted at the gas and liquid extraction line respectively to control the outlet gas quality. The diameter of the liquid orifice was set to 2. 50 mm and three gas orifices with different size( dG= 2. 65,5. 00,10. 00 mm) were tested. The experiments were carried out at an air-water two-phase flow loop. The gas superficial velocity ranged from 6. 0 to 20. 0 m /s and the liquid superficial velocity was in the range of 0. 02- 0. 18 m /s. Flow patterns such as wave flow,slug flow and annular flow were observed. The gas quality of the side arm branch was found mainly determined by the flow area ratio of the gas orifice to the liquid orifice and independent of gas and liquid superficial velocity,flow patterns and extraction flux.

Local Flow Regime Transition Criteria of Gas-Liquid Two-phase Flow in Vertical Upward Tube with a Horizontal Rod

The upward multiphase cross flow and heat transfer in the vertical tube may occur in oil production and chemical facilities. In this study, the local flow patterns of an upward gas-water two phase cross flow in a vertical tube with a horizontal rod have been investigated with an optical probe and the digital high speed video system. The local flow patterns are defined as the bubble, slug, churn and annular flow patterns. Optical probe signals are ana- lyzed in terms of probability density function, and it is proved that the local flow patterns can be recognized by this method. The transition mechanisms between the different flow patterns have been analyzed and the corresponding transitional models are proposed. Finally, local flow pattern maps of the upward gas-water two-phase flow in the vertical tube with a horizontal rod are constructed.

Numerical Simulation and Experimental Analysis of the Influence of Asymmetric Pressure Conditions on the Splitting of a Gas-Liquid Two-Phase Flow at a T-Junction

Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid velocities.The flow rates in the two branches were measured accurately to determine how the two considered phases distribute in the two outlets.The experimental results have shown that when the two outlet pressures are asymmetric,the two-phase flow always tends to flow into the outlet which has a lower pressure.As the inlet liquid velocity increases,however,the two-phase flow gradually tends to split evenly.Compared with the experiment results,the pressure difference between the two outlets can be determined more accurately by means of numerical simulation.The trends of experimental results and simulations are in very good agreement.

A Hydrodynamic Modei for Slug Frequency in Horizontal Gas-Liquid Two-Phase Flow

The prediction of slug frequency has important significance on gas-liquid two-phase flow. A hydrody-namic modei was put forward to evaluate slug frequency for horizontal two-phase flow, based on the dependence of slug frequency on the frequency of unstable interfacial wave. Using air and water, experimental verification of the modei was carried out in a large range of flow parameters. Six electrical probes were installed at different positions of a horizontal plexiglass pipe to detect slug frequency development. The pipe is 30m long and its inner diameter is 24 mm. It is observed experimentally that the interfacial wave frequency at the inlet is about l to 3 times the frequency of stable slug. The slug frequencies predicted by the modei fit well with Tronconi (1990) modei and the experimental data. The combination of the hydrodynamic modei and the experimental data results in a conclusion that the frequency of equilibrium liquid slug is approximately half the miniraum frequency of interfacial wave.

Local Flow Regime Transition Criteria of Gas-Liquid Two-phase Flow in Vertical Upward Tube with a Horizontal Rod

The upward multiphase cross flow and heat transfer in the vertical tube may occur in oil production and chemical facilities. In this study, the local flow patterns of an upward gas-water two phase cross flow in a vertical tube with a horizontal rod have been investigated with an optical probe and the digital high speed video system. The local flow patterns are defined as the bubble, slug, churn and annular flow patterns. Optical probe signals are ana- lyzed in terms of probability density function, and it is proved that the local flow patterns can be recognized by this method. The transition mechanisms between the different flow patterns have been analyzed and the corresponding transitional models are proposed. Finally, local flow pattern maps of the upward gas-water two-phase flow in the vertical tube with a horizontal rod are constructed.

Study on Venturi nozzle for gas-liquid two-phase flow metering

Experiments were carried out to study gas-liquid two-phase flow through a horizontally mounted Venturi nozzle. The inner diameter of pipe is 25 mm and the throat diameter of Venturi nozzle is 5. 1 mm. The pressure difference and pressure loss across the nozzle were measured. It was found that the degree of pressure fluctuation strongly depends on the gas quality. However,the relationship between the standard deviation of pressure difference and the gas quality is not monotonous. Multiple solutions may occur when the relationship was used to determine gas quality. On the other hand,the standard deviation of pressure loss was found to be monotonously correlated to X. This phenomenon was applied to measured gas quality. Also a modified Lin's equation is proposed to calculate the two-phase flow rate. The experimental measurements agree well with the predicted values.

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines.Friction factors play an important role in the accurate calculation of pressure drop.Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall,gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions.In this paper by modification of a friction model available in the literature,an improved semiempirical model is proposed.The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved.Pressure gradient and velocity profiles are validated against experimental data.Using the improved model,the pressure gradient deviation from experiments diminishes by about 3%;the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.

Numerical Investigation on Dynamic Response Characteristics of Fluid-Structure Interaction of Gas-Liquid Two-Phase Flow in Horizontal Pipe

Fluid-structure interaction(FSI)of gas-liquid two-phase fow in the horizontal pipe is investigated numerically in the present study.The volume of fluid model and standard k-e turbulence model are integrated to simulate the typical gas-liquid two-phase fow patterns.First,validation of the numerical model is conducted and the typical fow patterns are consistent with the Baker chart.Then,the FSI framework is established to investigate the dynamic responses of the interaction between the horizontal pipe and gas-liquid two-phase fow.The results show that the dynamic response under stratified fow condition is relatively flat and the maximum pipe deformation and equivalent stress are 1.8 mm and 7.5 MPa respectively.Meanwhile,the dynamic responses induced by slug fow,wave fow and annular fow show obvious periodic fuctuations.Furthermore,the dynamic response characteristics under slug flow condition are maximum;the maximum pipe deformation and equivalent stress can reach 4mm and 17.5 MPa,respectively.The principal direction of total deformation is different under various flow patterns.Therefore,the periodic equivalent stress will form the cyclic impact on the pipe wall and affect the fatigue life of the horizontal pipe.The present study may serve as a reference for FSI simulation under gas-liquid two-phase transport conditions.

Numerical investigation on characteristics of interfacial wave of liquid film in gas-liquid two-phase flow using OpenFOAM

Liquid film cooling as an advanced cooling technology is widely used in space vehicles.Stable operation of liquid film along the rocket combustion inner wall is crucial for thermal protection of rocket engines.The stability of liquid film is mainly determined by the characteristics of interfacial wave,which is rarely investigated right now.How to improve the stability of thin film has become a hot spot.In view of this,an advanced model based on the conventional Volume of Fluid(VOF)model is adopted to investigate the characteristics of interfacial wave in gas-liquid flow by using OpenFOAM,and the mechanism of formation and development of wave is revealed intuitively through numerical study.The effects from gas velocity,surface tension and dynamic viscosity of liquid(three factors)on the wave are studied respectively.It can be found that the gas velocity is critical to the formation and development of wave,and four modes of droplets generation are illustrated in this paper.Besides,a gas vortex near the gas-liquid interface can induce formation of wave easily,so changing the gas vortex state can regulate formation and development of wave.What’s more,the change rules of three factors influencing on the interfacial wave are obtained,and the surface tension has a negative effect on the formation and development of wave only when the surface tension coefficient is above the critical value,whereas the dynamic viscosity has a positive effect in this process.Lastly,the maximum height and maximum slope angle of wave will level off as the gas velocity increases.Meanwhile,the maximum slope angle of wave is usually no more than 38°,no matter what happens to the three factors.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.