Analysis of forced-convection boiling flow instabilities in a single-channel upflow system

A numerical model is developed to predict the steady-state and transient behaviour of forced-convection boiling two-phase flow in a single channel. The model is based on the assumption of homogeneous two-phase flow and thermodynamic equilibrium of the phases. Compressibility effects in the two-phase region, motion of the bulk boiling interface and the thermal capacity of the heater wall have been included in the analysis. The model is used to study the effects of heat input, inlet subcooling and flow rate on the system behaviour. For comparison purposes, an experimental investigation was conducted using a single-channel, electrically heated, forced-convection upflow system. Steady-state operating characteristics, and stable and unstable regions, are determined as a function of heat flux, inlet subcooling and mass flow rate. Different modes of oscillation and their characteristics have been investigated. The model's predictions are in good agreement with the experimental results.     

Turbulence structure of air-water bubbly flow—II. local properties

Microstructure was studied experimentally in air-water two-phase bubbly flow flowing upward in a vertical pipe of 60 mm diameter under atmospheric pressure. The results indicate that over a large portion of fully-developed bubbly flow the phases, the velocities of bubbles and water, and the ratio between the velocities of the phases have fairly flat radial profiles. In the wall region a maximum void fraction was observed. Spectra of the velocities of bubbles and water showed a Poisson distribution and a normal distribution function, respectively. The experimental evidence indicated a trend for the turbulent intensity to decrease first with increasing gas flow rate for constant water velocity and to increase again with a further increase in the gas flow rate. This phenomenon was more significant for a higher water velocity.     

A Mathematical Model for Determining Oil Migration Characteristics in Low-Permeability Porous Media Based on Fractal Theory

Shale reservoirs are characterized by very low permeability in the scale of nano-Darcy. This is due to the nanometer scale of pores and throats in shale reservoirs, which causes a difference in flow behavior from conventional reservoirs. Slip flow is considered to be one of the main flow regimes affecting the flow behavior in shale gas reservoirs and has been widely studied in the literature. However, the important mechanism of gas desorption or adsorption that happens in shale reservoirs has not been investigated thoroughly in the literature. This paper aims to study slip flow together with gas desorption in shale gas reservoirs using pore network modeling. To do so, the compressible Stokes equation with proper boundary conditions was applied to model gas flow in a pore network that properly represents the pore size distribution of typical shale reservoirs. A pore network model was created using the digitized image of a thin section of a Berea sandstone and scaled down to represent the pore size range of shale reservoirs. Based on the size of pores in the network and the pore pressure applied, the Knudsen number which controls the flow regimes was within the slip flow regime range. Compressible Stokes equation with proper boundary conditions at the pore’s walls was applied to model the gas flow. The desorption mechanism was also included through a boundary condition by deriving a velocity term using Langmuir-type isotherm. It was observed that when the slip flow was activated together with desorption in the model, their contributions were not summative. That, is the slippage effect limited the desorption mechanism through a reduction of pressure drop. Eagle Ford and Barnett shale samples were investigated in this study when the measured adsorption isotherm data from the literature were used. Barnett sample showed larger contribution of gas desorption toward gas recovery as compared to Eagle Ford sample. This paper has produced a pore network model to further understand the gas desorption and the slip flow effects in recovery of shale gas reservoirs.

     

Pore-Network Modeling of Isothermal Drying in Porous Media

In this paper we present numerical results obtained with a pore-network model for the drying of porous media that accounts for various processes at the pore scale. These include mass transfer by advection and diffusion in the gas phase, viscous flow in the liquid and gas phases and capillary effects at the liquid--gas interface. We extend our work by studying the effect of capillarity-induced flow in macroscopic liquid films that form at the pore walls as the liquid--gas interface recedes. A mathematical model that accounts for the effect of films on the drying rates and phase distribution patterns is presented. It is shown that film flow is a major transport mechanism in the drying of porous materials, its effect being dominant when capillarity controls the process, which is the case in typical applications.     

Rheological properties of some biphenyl liquid crystals

Rheological measurements have been carried out using a rotational viscometer with a system of coaxial cylinders on four liquid crystalline substances from the group of cyanobiphenyls. On the basis of results of these investigations it was found that in the investigated range of shear rates the nematic phases exhibited Newtonian flow behaviour, while the smectic phases exhibited non-Newtonian behaviour. For shear rates up to ca. 1000 s^–1 the dependence of the shear stress on shear rate is well described by a power-law model.     

The stratified flow of gas and non-newtonian liquid in horizontal pipes

Predictions of pressure drop and holdup are presented for the stratified flow of gas and non-Newtonian liquid obeying the Ostwald-de Waele power law model. The model of Taitel & Dukler (1976) for gas/Newtonian liquid flow is extended to liquids possessing either shear-thinning or shear-thickening laminar flow behaviour and computed results are given for flow behaviour indices in the range $\text{0.1 ≤ n ≤ 2}$. In particular, conditions are defined for drag reduction of the liquid flow by the presence of the gas. It is concluded that drag reduction occurs over the largest ranges of liquid and gas flow rates at the lowest $\text{n}$ values, provided that liquid flow remains laminar, but that maximum drag reduction may be expected for shear-thickening liquids with $\text{n}$ values of 2 or greater. Ratios of the liquid flow rate in the presence of gas to that for liquid flow alone under a constant pressure gradient are also presented. These ratios frequently exceed unity and are greatest for highly shear-thinning liquids.Although the Taitel & Dukler approach is consistent with experiments on gas/Newtonian liquid flow, and, in addition, appears to be valid for immiscible Newtonian liquid-liquid systems, provided that the viscosity ratio of the two phases is at least five, experiments are required to confirm its applicability for gas/non-Newtonian systems.     

Hydrodynamics and energy consumption studies in a three-phase liquid circulating three-phase fluid bed contactor

The hydrodynamics and energy consumption have been studied in a cold flow, bubbling and turbulent, pressurized gas–liquid–solid three-phase fluidized bed (0.15 m ID × 1 m height) with concurrent gas–liquid up flow is proposed with the intention of increasing the gas hold up. The hydrodynamic behaviour is described and characterised by some specific gas and liquid velocities. Particles are easily fluidized and can be uniformly distributed over the whole height of the column. The effect of parameters like liquid flow rate, gas flow rate, particle loading, particle size, and solid density on gas hold up and effect of gas flow rate, solid density and particle size on solid hold up, energy consumption and minimum fluidization velocity has been studied. At the elevated pressures a superior method for better prediction of minimum fluidization velocity and terminal settling velocities has been adopted. The results have been interpreted with Bernoulli’s theorem and Richardson–Zaki equation. Based on the assumption of the gas and liquid as a pretend fluid, a simplification has been made to predict the particle terminal settling velocities. The Richardson–Zaki parameter n′ was compared with Renzo’s results. A correlation has been proposed with the experimental results for the three-phase fluidization.     

Dynamic Relative Permeability and Simulation of WAG Injection Processes

In this work, we investigate the impact of mobility changes due to flow reversals from co-current to counter-current flow on the displacement performance of water alternating gas (WAG) injection processes. In WAG processes, the injected gas will migrate toward the top of the formation while the injected water will migrate toward the bottom of the formation. The segregation of gas, oil and water phases will result in counter-current flow occurring in the vertical direction in some portions of the reservoir during the displacement process. Previous experimental and theoretical studies of counter-current flow have shown that the relative mobility of each of the phases in a porous medium is considerably less when counter-current flow prevails as compared to co-current flow settings. A reduction of the relative permeability in the vertical direction results in a dynamic anisotropy in phase mobilities. This effect has, to the best of our knowledge, not previously been considered in the modeling and simulation of WAG processes. A new flow model that accounts for flow reversals in the vertical direction has been implemented and tested in a three-phase compositional reservoir simulator. In order to investigate the impact of flow reversals, results from the new flow model are compared to cases where counter-current flow effects on the phase mobilities are ignored. A range of displacement settings, covering relevant slug sizes, have been investigated to gauge the impact of mobility reductions due to flow reversals. Significant differences, in terms of saturation distribution, producing GOR and oil recovery, are observed between the conventional flow model (ignoring mobility reductions due to counter-current flow) and the proposed new model that accounts for reductions in phase mobility during counter-current flow. Accordingly, we recommend that an explicit representation of flow transitions between co-current and counter-current flow (and the related impact on phase mobilities) should be considered to ensure accurate and optimal design of WAG injection processes.     

Automaton Simulations of Dispersion in Porous Media

A thermodynamic lattice gas (automaton) model is used to simulate dispersion in porous media. Simulations are constructed at two distinctly different scales, the pore scale at which capillary models are constructed and large scale or Darcy scale at which probabilistic collision rules are introduced. Both models allow for macroscopic (pore scale) phase separation. The pore scale models clearly show the effect of pore structure on dispersion. The large scale (mega scale) simulations indicate that when the pressure difference between the displacing phase and displaced phase is properly chosen (representing the average pressure gradient between the phases). The simulation results are consistent with both theoretical predictions and experimental observations.     

Micro Bubbles in Solution–Gas Drive in Heavy Oil: Their Existence and Importance

This work investigates the existence and importance of the micro bubbles in heavy oil subjected to solution?Cgas drive. The term ??micro bubble?? will be used to refer to the free-gas phase that flows with the oil, no faster and no slower. Two types of experiments are reported here; slow and fast experiments. These experiments were previously reported by Sheikha and Pooladi-Darvish (SPE Res Eval Eng 12(3):390?C398, 2009) and were used to investigate the effects of pressure gradient and depletion rate of oil recovery. In this work, we investigate the nature of two-phase flow, and find that in the slow runs, the flow was characterized by single-phase flow of oil until a gas saturation of 2±1% was reached. Above this gas saturation, bulk flow of gas was observed at mobilities much higher than that of micro bubbles. Recovery factor of the slow tests was below 4%. In the fast runs, flow of bubbles is observed shortly after they are formed in porous media. The gas mobility and fractional flow remain low until a gas saturation of 7±1% is reached. Flow of gas between approximately 2 and 7% gas saturation is consistent with that of micro bubbles. Gas fractional flow increases sharply at gas saturations above approximately 7%. The results indicate that the attainment of high recovery values (12?C14%) observed in the fast experiments is partially as a result of low mobility of micro-gas bubbles. The pressure decline rate of each flow experiment was varied independent of its respective withdrawal rate. This did not alter the difference in recovery and mobility behaviour of the fast and slow experiments; the fast experiments exhibited a significant period of low mobility gas flow consistent with flow of micro bubbles. Regardless of the pressure decline rate, the slow experiments did not exhibit this period of low mobility gas flow.     

Experimental investigation of horizontal air–water bubbly-to-plug and bubbly-to-slug transition flows in a 3.81 cm ID pipe

The present study seeks to investigate horizontal bubbly-to-plug and bubbly-to-slug transition flows. The two-phase flow structures and transition mechanisms in these transition flows are studied based on experimental database established using the local four-sensor conductivity probe in a 3.81 cm inner diameter pipe. While slug flow needs to be distinguished from plug flow due to the presence of large number of small bubbles (and thus, large interfacial area concentration), both differences and similarities are observed in the evolution of interfacial structures in bubbly-to-plug and bubbly-to-slug transitions. The bubbly-to-plug transition is studied by decreasing the liquid flow rate at a fixed gas flow rate. It is found that as the liquid flow rate is lowered, bubbles pack near the top wall of the pipe due to the diminished role of turbulent mixing. As the flow rate is lowered further, bubbles begin to coalesce and form the large bubbles characteristic of plug flow. Bubble size increases while bubble velocity decreases as liquid flow rate decreases, and the profile of the bubble velocity changes its shape due to the changing interfacial structure. The bubbly-to-slug transition is investigated by increasing the gas flow rate at a fixed liquid flow rate. In this transition, gas phase becomes more uniformly distributed throughout the cross-section due to the formation of large bubbles and the increasing bubble-induced turbulence. The size of small bubbles decreases while bubble velocity increases as gas flow rate increases. The distributions of bubble size and bubble velocity become more symmetric in this transition. While differences are observed in these two transitions, similarities are also noticed. As bubbly-to-plug or bubbly-to-slug transition occurs, the formation of large elongated bubbles is observed not in the uppermost region of bubble layer, but in a lower region. At the beginning of transitions, relative differences in phase velocities near the top of the pipe cross-section to those near the pipe center become larger for both gas and liquid phases, because more densely packed bubbles introduce more resistance to both phases.     

页岩有机质纳米孔隙气体吸附与流动规律研究

页岩储层孔隙结构复杂, 气体赋存方式多样. 有机质孔隙形状对受限空间气体吸附和流动规律的影响尚不明确, 导致难以准确认识页岩气藏气体渗流机理. 为解决该问题, 本文首先采用巨正则蒙特卡洛方法模拟气体在不同形状有机质孔隙(圆形孔隙、狭长孔隙、三角形孔隙、方形孔隙)内吸附过程, 发现不同形状孔隙内吸附规律符合朗格缪尔单层吸附规律, 分析了绝对吸附量、过剩吸附浓量、气体吸附参数随孔隙尺寸、压力的变化, 研究了孔隙形状对气体吸附的影响. 在明确不同形状有机质孔隙内气体热力学吸附规律基础上, 建立不同形状有机质孔隙内吸附气表面扩散数学模型和考虑滑脱效应的自由气流动数学模型, 结合分子吸附模拟结果研究了不同孔隙形状、孔隙尺寸有机质孔隙内吸附气流动与自由气流动对气体渗透率的贡献. 结果表明, 狭长孔隙内最大吸附浓度和朗格缪尔压力最高, 吸附气表面扩散能力最弱. 孔隙半径5 nm以上时, 吸附气表面扩散对气体渗透率影响可忽略. 本文研究揭示了页岩气藏实际生产过程中有机质孔隙形状对页岩气吸附和流动能力的影响机制.      

Two-phase flow instabilities in a horizontal single boiling channel

A study of the stability of an electrically heated single channel, forced convection horizontal system was conducted by using Freon-11 as the test fluid. Two major modes of oscillations, namely, density-wave type (high frequency) and pressure-drop type (low frequency) oscillations have been observed. The steady-state operating characteristics and stable and unstable regions are determined as a function of heat flux, exit orifice diameter and mass flow rate. Different modes of oscillations and their characteristics have been investigated. The effect of the exit restriction on the system stability has also been studied.A mathematical model has been developed to predict the transient behavior of boiling two-phase systems. The model is based on homogenous flow assumption and thermodynamic equilibrium between the liquid and vapor phases. The transient characteristics of boiling two-phase flow horizontal system are obtained for various heat inputs, flow rates and exit orifice diameters by perturbing the governing equations around a steady state. Theoretical and experimental results have been compared.     

Modelling Foamy Oil Flow in Porous Media

This paper describes a dynamic model for the simulation of foamy oil flow in porous media. The model includes expressions for the rate processes of nucleation, bubble growth and disengagement of dispersed gas bubbles from the oil. The model is used to simulate experimental results pertaining to primary depletion tests conducted in a sand pack. Using the model to interpret experimental results indicated that, although the lifetimes of supersaturation and dispersed gas bubbles may be short, supersaturated conditions are likely to exist, and dispersed gas bubbles are likely to be present during the entire production period, as long as the pressure continues to decline at a high rate. The model developed in this paper gave better agreement with experimental data than other proposed models. The effect of foamy oil flow increases as the rate of pressure decline increases.     

Spontaneous Inertial Imbibition in Porous Media Using a Fractal Representation of Pore Wall Rugosity

Considering the separable phenomena of imbibition in complex fine porous media as a function of timescale, it is noted that there are two discrete imbibition rate regimes when expressed in the Lucas–Washburn (L–W) equation. Commonly, to account for this deviation from the single equivalent hydraulic capillary, experimentalists propose an effective contact angle change. In this work, we consider rather the general term of the Wilhelmy wetting force regarding the wetting line length, and apply a proposed increase in the liquid–solid contact line and wetting force provided by the introduction of surface meso/nanoscale structure to the pore wall roughness. An experimental surface pore wall feature size regarding the rugosity area is determined by means of capillary condensation during nitrogen gas sorption in a ground calcium carbonate tablet compact. On this nano size scale, a fractal structure of pore wall is proposed to characterize for the internal rugosity of the porous medium. Comparative models based on the Lucas–Washburn and Bosanquet inertial absorption equations, respectively, for the short timescale imbibition are constructed by applying the extended wetting line length and wetting force to the equivalent hydraulic capillary observed at the long timescale imbibition. The results comparing the models adopting the fractal structure with experimental imbibition rate suggest that the L–W equation at the short timescale cannot match experiment, but that the inertial plug flow in the Bosanquet equation matches the experimental results very well. If the fractal structure can be supported in nature, then this stresses the role of the inertial term in the initial stage of imbibition. Relaxation to a smooth-walled capillary then takes place over the longer timescale as the surface rugosity wetting is overwhelmed by the pore condensation and film flow of the liquid ahead of the bulk wetting front, and thus to a smooth walled capillary undergoing permeation viscosity-controlled flow.     

Understanding Oil and Gas Flow Mechanisms in Shale Reservoirs Using SLD–PR Transport Model

More and more attention has been paid to the oil and gas flow mechanisms in shale reservoirs. The solid–fluid interaction becomes significant when the pores are in the nanoscale. The interaction changes the fluid’s physical properties and leads to different flow mechanisms in shale reservoirs from those in conventional reservoirs. By using a Simplified Local Density–Peng Robinson transport model, we consider the density and viscosity profiles, which result from solid–fluid interaction. Gas rarefaction effect is negligible at high pressure, so we assume it is viscous flow. Considering the density- and viscosity-changing effects, we proposed a slit permeability model. The velocity profiles are obtained by this newly established model. This proposed model is validated by matching the density profile and velocity profile from molecular dynamic simulation. Then, the effects of pressure and pore size on gas and oil flow mechanisms are also studied in this work. The results show that both gas and oil exhibit enhanced flow rates in nanopores. Gas-phase flow in nanopores is dominated by the density-changing effect (adsorption), while the oil-phase flow is mainly controlled by the viscosity-changing effect. Both gas and oil permeability quickly decrease to the Darcy permeability when the slit aperture becomes large. The results reported in this work are representative and should significantly help us understand the mechanisms of oil and gas flow in shale reservoirs.     

Filtration Behaviour of Cement-Based Grout in Porous Media

Filtration behaviour of cement particles, especially under the high grouting pressure with a rapid grout flow velocity, has a significant effect on the grout injection. However, there have been few studies on this field where the governing equation of this behaviour remains unclear. In the present study, a novel experimental procedure for grout injection was adopted to acquire the spatial and temporal variations in porosity and viscosity of high-speed grout flow in coarse sand. Experimental observations showed that there were dramatic variations in viscosity and porosity during the grout penetration within the first 50 s, suggesting that the high velocity had a significant influence on the distribution of the filtration coefficient. A model based on the Stokes–Brinkman (S–B) equation and advection–filtration equations was established to describe the filtration of grout flow in porous media. Meanwhile, numerical solutions from both the proposed model and traditional Darcy’s law were compared with experimental results. The comparative results showed that the proposed approach can match the laboratory tests well; the analysis indicated that Darcy’s law was unable to properly describe high-speed grout flow in porous media due to the lack of a shear force and the inertial term. Nonuniform filtration behaviour of cement grout flowing in porous media was revealed. Due to the nonuniform distribution of the pore velocity isoline caused by Poiseuille flow, it led to a heterogenous distribution of porosity as well. Parametric studies on the applicability of Darcy’s law and S–B equation for grout flow were discussed, in which an error of less than 10% was calculated when the Reynolds number was less than 2.5.     

Formation of bubbles at submerged orifices - Experimental investigation and theoretical prediction

The present work reports an experimental investigation on bubble release through submerged orifices. Bubble frequency has been measured as a function of gas flow rate for three different orifice sizes at various pool heights. Needle type conductivity probe has been used for bubble count. Analysis of probe signal not only gives the bubble frequency but also indicates a transition from bubbling to jetting regime. Further, to validate the experimental observations a simple mechanistic model has been developed considering the evolution of non-spherical bubbles at the orifice mouth. Reasonable agreement between the model prediction and the experimental result has been observed.     

Pore Network Modelling of Electrical Resistivity and Capillary Pressure Characteristics

In order to model petrophysical properties of hydrocarbon reservoir rocks, the underlying physics occurring in realistic rock pore structures must be captured. Experimental evidence showing variations of wetting occurring within a pore, and existence of the so-called 'non-Archie' behaviour, has led to numerical models using pore shapes with crevices (for example, square, elliptic, star-like shapes, etc.). This paper presents theoretical derivations and simulation results of a new pore space network model for the prediction of petrophysical properties of reservoir rocks. The effects of key pore geometrical factors such as pore shape, pore size distribution and pore co-ordination number (pore connectivity) have been incorporated into the theoretical model. In particular, the model is used to investigate the effects of wettability and saturation history on electrical resistivity and capillary pressure characteristics. The petrophysical characteristics were simulated for reservoir rock samples. The use of the more realistic grain boundary pore (GBP) shape allows simulation of the generic behaviour of sandstone rocks, with various wetting scenarios. The predictions are in close agreement with electrical resistivity and capillary pressure characteristics observed in experiments.