深水气井清井放喷过程井筒气液流动型态模拟研究

针对深水气井液气比变化剧烈的清井放喷过程,以井筒内部气液流动型态分布规律为目标,结合深水气井清井测试工况,开展了清井放喷瞬态数值仿真模拟与分析。结果表明:放喷后井口处和泥线处的温度呈现不同的温度变化特征,井口温度出现先降低后升高,而泥线处流温先升后降;低井口压力下井筒内更易出现段塞流、环状流,高井口压力下返排初期往往呈现泡状流,中后期逐渐转变为段塞和环状流型;同一放喷时刻,井身沿程也将呈现不同流型分布,井筒下部主要呈现泡状流、段塞流,中上部范围内主要为环状流型。

深水气井测试返排阶段水合物生成及防治研究

清井返排是深水气井测试流程的重要阶段,该阶段管柱内水合物的生成及堵塞问题一直是安全测试研究的重点。从水合物生成、生长、沉积、堵塞的微观机理分析出发,结合清井返排的具体阶段,首次将水合物形成与管柱气液流型结合研究不同流型下水合物形成流动障碍的风险大小。研究结果表明,泡状流、段塞流、搅拌流下水合物形成流动障碍的风险较小,环状流下仍然具有形成测试管柱堵塞的风险,须进行合理的测试参数设计,并在不同阶段实施防控措施,从而有效避免水合物堵塞事故。本文研究成果可为现场测试施工提供理论支撑。

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.

Effects of Orifice Orientation and Gas-Liquid Flow Pattern on Initial Bubble Size

In many gasliquid processes, the initial bubble size is determined by a series of operation parameters along with the sparger design and gasliquid flow pattern. Bubble formation models for variant gasliquid flow pat terns have been developed based on force balance. The effects of the orientation of gasliquid flow, gas velocity, liquid velocity and orifice diameter on the initial bubble size have been clarified. In ambient airwater system, thesultable gasllquid flow pattern is important to obtain smaller bubbles under the low velocity liquid crossflow con ditions with stainless steel spargers. Among the four types of gasliquid flow patterns discussed, the horizontal orifice in a vertically upward liquid flow produces the smallest initial bubbles. However the orientation effects of gas and liquid flow are found tobe insgnifican whenliq.uid velocity is.higher than. 3.2 m;sa or theorifice diameter is small enough.

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.

Computational fluid dynamics simulation of gas-liquid two phases flow in 320 m^3 air-blowing mechanical flotation cell using different turbulence models

According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.

Application of time–frequency entropy from wake oscillation to gas–liquid flow pattern identification

Gas–liquid two-phase flow abounds in industrial processes and facilities. Identification of its flow pattern plays an essential role in the field of multiphase flow measurement. A bluff body was introduced in this study to recognize gas–liquid flow patterns by inducing fluid oscillation that enlarged differences between each flow pattern. Experiments with air–water mixtures were carried out in horizontal pipelines at ambient temperature and atmospheric pressure. Differential pressure signals from the bluff-body wake were obtained in bubble, bubble/plug transitional, plug, slug, and annular flows. Utilizing the adaptive ensemble empirical mode decomposition method and the Hilbert transform, the time–frequency entropy S of the differential pressure signals was obtained. By combining S and other flow parameters, such as the volumetric void fraction β, the dryness x, the ratio of density φ and the modified fluid coefficient ψ, a new flow pattern map was constructed which adopted S（1–x）φ and （1–β）ψ as the vertical and horizontal coordinates, respectively. The overall rate of classification of the map was verified to be 92.9% by the experimental data. It provides an effective and simple solution to the gas–liquid flow pattern identification problems.

A Model of Turbulent-Laminar Gas-Liquid Stratified Flow

The time-dependent liquid film thickness and pressure drop are measured by using parallel-wire conduc tance probes and capacitance differential-pressure transducer. A mathematical model with iterative procedure to calculate holdup and pressure drop in horizontal and inclined gas-liquid stratified flow is developed. The predictions agree well with over a hundred experimental data in 0.024 and 0.04 m diameter pipelines.

A dual-scale turbulence model for gas–liquid bubbly flows

A dual-scale turbulence model is applied to simulate cocurrent upward gas-liquid bubbly flows and validated with available experimental data. In the model, liquid phase turbulence is split into shear-induced and bubble- induced turbulence. Single-phase standard k-e model is used to compute shear-induced turbulence and another transport equation is added to model bubble-induced turbulence. In the latter transport equation, energy loss due to interface drag is the production term, and the characteristic length of bubble-induced turbulence, simply the bubble diameter in this work, is introduced to model the dissipation term. The simulated results agree well with experimental data of the test cases and it is demonstrated that the proposed dual-scale turbulence model outperforms other models. Analysis of the predicted turbulence shows that the main part of turbulent kinetic en- ergy is the bubble-induced one while the shear-induced turbulent viscosity predominates within turbulent vis- cosity, especially at the pipe center. The underlying reason is the apparently different scales for the two kinds of turbulence production mechanisms： the shear-induced turbulence is on the scale of the whole pipe while the bubble-induced turbulence is on the scale of bubble diameter. Therefore, the model reflects the multi-scale phe- nomenon involved in gas-liquid bubbly flows.

A Model for Predicting Holdup and Pressure Drop in Gas- Liquid Stratified Flow

The time-dependent liquid film thickness and pressure drop were measured by using parallel-wire conductance probes and capacitance differential-pressure transducers. Applying the eddy viscosity theory and an appropriate correlation of interfacial sear stress,a new two-dimensional separated model of holdup and pressure drop of turbulent/turbulent gas-liquid stratified flow was presented. Prediction results agreed well with experimental data.

Parametric study of gas−liquid two-phase flow field in horizontal stirred tank

Based on Fluent software,the gas−liquid two-phase flow in the horizontal stirred tank was simulated with SST k−ωturbulence model,Eulerian−Eulerian two-fluid model,and multi-reference flame method.The mixing process in the tank was calculated by tracer method.The results show that increasing the rotating speed or gas flow is conducive to a more uniform distribution of the gas phase and accelerates the mixing of the liquid phase.When the rotating speed exceeds 93 r/min,the relative power demand remains basically constant.The change in the inclination angle of the upper impeller has minimal effect on the gas phase distribution.When the inclination angle is 50°,the relative power demand reaches the maximum.An appropriate increase in the impeller distance from the bottom improves the gas holdup and gas phase distribution but increases the liquid phase mixing time.

Numerical Simulations on Countercurrent Gas-Liquid Flow in a PWR Hot Leg with Air-Water Flow in a 1/15th Scale-Model

To clarify the countercurrent flow in a PWR hot leg under reflux condensation, numerical simulations of countercurrent air-water flow for a 1/15th scale model of the PWR hot leg were conducted using the two-fluid model implemented in CFD software. In this paper, the effect of expansion of the inclined pipe, which is the actual plant geometry, was evaluated. When increasing the air velocity, CCFL characteristics and the mechanism of flow pattern transition had significant differences between the case with and without expansion of the inclined pipe. CCFL characteristics were mitigated in the case with expansion. The effect of computational grid size was also discussed. When the supplied water velocity was small, the predicted flow pattern transition point agreed well with the measured data by increasing the number of cells. On the other hand, when the air velocity was decreasing, there were no significant differences in each case.

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.

Three-phase flow of submarine gas hydrate pipe transport

In the hydraulic transporting process of cutter-suction mining natural gas hydrate, when the temperature-pressure equilibrium of gas hydrate is broken, gas hydrates dissociate into gas. As a result, solid-liquid two-phase flow(hydrate and water) transforms into gas-solid-liquid three-phase flow(methane, hydrate and water) inside the pipeline. The Euler model and CFD-PBM model were used to simulate gas-solid-liquid three-phase flow. Numerical simulation results show that the gas and solid phase gradually accumulate to the center of the pipe. Flow velocity decreases from center to boundary of the pipe along the radial direction. Comparison of numerical simulation results of two models reveals that the flow state simulated by CFD-PBM model is more uniform than that simulated by Euler model, and the main behavior of the bubble is small bubbles coalescence to large one. Comparison of numerical simulation and experimental investigation shows that the values of flow velocity and gas fraction in CFD-PBM model agree with experimental data better than those in Euler model. The proposed PBM model provides a more accurate and effective way to estimate three-phase flow state of transporting gas hydrate within the submarine pipeline.