Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettabi...Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettability alteration by surfactants. Although the interacting capillary bundle(ICB) model shows potential in characterizing imbibition rates in different pores during wettability alteration, the existing ICB models neglect the influence of wettability and viscosity ratio on the imbibition behavior, making it difficult to accurately describe the oil-water imbibition behavior within the porous media. In this work,a new ICB mathematical model is established by introducing pressure balance without assuming the position of the leading front to comprehensively describe the imbibition behavior in a porous medium under different conditions, including gas-liquid spontaneous imbibition and oil-water imbibition.When the pore size distribution of a tight rock is known, this new model can predict the changes of water saturation during the displacement process in the tight rock, and also determine the imbibition rate in pores of different sizes. The water saturation profiles obtained from the new model are validated against the waterflooding simulation results from the CMG, while the imbibition rates calculated by the model are validated against the experimental observations of gas-liquid spontaneous imbibition. The good match above indicates the newly proposed model can show the water saturation profile at a macroscopic scale while capture the underlying physics of the multiphase flow in a porous medium at a microscopic scale. Simulation results obtained from this model indicate that both wettability and viscosity ratio can affect the sequence of fluid imbibition into pores of different sizes during the multiphase flow, where less-viscous wetting fluid is preferentially imbibed into larger pores while more-viscous wetting fluid tends to be imbibed into smaller pores. Furthermore, this model provides an avenue to calculate the imbibition rate in pores of different sizes during wettability alteration and capture the non-Darcy effect in micro-and nano-scale pores.展开更多
A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten st...A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.展开更多
Considering the phase behaviors in condensate gas reservoirs and the oil-gas two-phase linear flow and boundary-dominated flow in the reservoir,a method for predicting the relationship between oil saturation and press...Considering the phase behaviors in condensate gas reservoirs and the oil-gas two-phase linear flow and boundary-dominated flow in the reservoir,a method for predicting the relationship between oil saturation and pressure in the full-path of tight condensate gas well is proposed,and a model for predicting the transient production from tight condensate gas wells with multiphase flow is established.The research indicates that the relationship curve between condensate oil saturation and pressure is crucial for calculating the pseudo-pressure.In the early stage of production or in areas far from the wellbore with high reservoir pressure,the condensate oil saturation can be calculated using early-stage production dynamic data through material balance models.In the late stage of production or in areas close to the wellbore with low reservoir pressure,the condensate oil saturation can be calculated using the data of constant composition expansion test.In the middle stages of production or when reservoir pressure is at an intermediate level,the data obtained from the previous two stages can be interpolated to form a complete full-path relationship curve between oil saturation and pressure.Through simulation and field application,the new method is verified to be reliable and practical.It can be applied for prediction of middle-stage and late-stage production of tight condensate gas wells and assessment of single-well recoverable reserves.展开更多
Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information whe...Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information when the reservoir exhibits reservoir heterogeneity and multiphase flow simultaneously.This work proposes a practical-yet-rigorous method to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,and improve the evaluation of reservoir/fracture properties.A new,general,semi-analytical model is proposed that explicitly accounts for multiphase flow,fractalbased reservoir heterogeneity,anomalous diffusion,and pressure-dependent fluid properties.This is achieved by introducing a new Boltzmann-type transformation,the exponent of which includes reservoir heterogeneity and anomalous diffusion.In order to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,the modified Boltzmann variable allows the conversion of three partial differential equations(PDE's)(i.e.,oil,gas and water diffusion equations)into ordinary differential equations(ODE's)that are easily solved using the Runge-Kutta(RK)method.A modified time-power-law plot is also proposed to estimate the reservoir and fracture properties,recognizing that the classical square-root-of-time-plot is no longer valid when various reservoir complexities are exhibited simultaneously.Using the slope of the straight line on the modified time-power-law plot,the linear flow parameter can be estimated with more confidence.Moreover,because of the new Boltzmann-type transformation,reservoir and fracture properties can be derived more efficiently without the need for defining complex pseudo-variable transformations.Using the new semi-analytical model,the effects of multiphase flow,reservoir heterogeneity and anomalous diffusion on rate-decline behavior are evaluated.For the case of approximately constant flowing pressure,multiphase flow impacts initial oil rate,which is a function of oil relative permeability and well flowing pressure.However,multiphase flow has a minor effect on the oil production decline exponent.Reservoir heterogeneity/anomalous diffusion affect both the initial oil production rate and production decline exponent.The production decline exponent constant is a function of reservoir heterogeneity/anomalous diffusion only.The practical significance of this work is the advancement of RTA techniques to allow for more complex reservoir scenarios,leading to more accurate production forecasting and better-informed capital planning.展开更多
A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a...A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.展开更多
The research on the multiphase flow characteristics of hydrate slurry is the key to implementing the risk prevention and control technology of hydrate slurry in deep-water oil and gas mixed transportation system.This ...The research on the multiphase flow characteristics of hydrate slurry is the key to implementing the risk prevention and control technology of hydrate slurry in deep-water oil and gas mixed transportation system.This paper established a geometric model based on the high-pressure hydrate slurry experimental loop.The model was used to carry out simulation research on the flow characteristics of gas-liquid-solid three-phase flow.The specific research is as follows:Firstly,the effects of factors such as slurry flow velocity,hydrate particle density,hydrate particle size,and hydrate volume fraction on the stratified smooth flow were specifically studied.Orthogonal test obtained particle size has the most influence on the particle concentration distribution.The slurry flow velocity is gradually increased based on stratified smooth flow.Various flow patterns were observed and their characteristics were analyzed.Secondly,increasing the slurry velocity to 2 m/s could achieve the slurry flow pattern of partial hydrate in the pipeline transition from stratified smooth flow to wavy flow.When the flow rate increases to 3 m/s,a violent wave forms throughout the entire loop.Based on wave flow,as the velocity increased to 4 m/s,and the flow pattern changed to slug flow.When the particle concentration was below 10%,the increase of the concentration would aggravate the slug flow trend;if the particle concentration was above 10%,the increase of the concentration would weaken the slug flow trend,the increase of particle density and liquid viscosity would weaken the tendency of slug flow.The relationship between the pressure drop gradients of several different flow patterns is:slug flow>wave flow>stratified smooth flow.展开更多
In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer...In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer were added in the apparatus to follow the real flow line.Virtue tracers were considered in the mathematical model and the algorithm of tracers was built.The comparison of the results between the experiment and numerical calculation shows that the time of the tracer flows out of stirring tube are 40 s in the experiment and 42 s in numerical calculated result.The transient diffusion process and the solution residence time of the numerical calculation are in good agreement with the experimental results,which indicates that the mathematical model is reliable and can be used to predict the flow field of the air-agitated seed precipitation tank.展开更多
The venturi meter has an advantage in its use,because it can measure flow without being much affected by the type of the measured fluid or flow conditions.Hence,it has excellent versatility and is being widely applied...The venturi meter has an advantage in its use,because it can measure flow without being much affected by the type of the measured fluid or flow conditions.Hence,it has excellent versatility and is being widely applied in many industries.The flow of a liquid containing air is a representative example of a multiphase flow and exhibits complex flow characteristics.In particular,the greater the gas volume fraction（GVF）,the more inhomogeneous the flow becomes.As a result,using a venturi meter to measure the rate of a flow that has a high GVF generates an error.In this study,the cause of the error occurred in measuring the flow rate for the multiphase flow when using the venturi meter for analysis by CFD.To ensure the reliability of this study,the accuracy of the multiphase flow models for numerical analysis was verified through comparison between the calculated results of numerical analysis and the experimental data.As a result,the Grace model,which is a multiphase flow model established by an experiment with water and air,was confirmed to have the highest reliability.Finally,the characteristics of the internal flow Held about the multiphase flow analysis result generated by applying the Grace model were analyzed to find the cause of the uncertainty occurring when measuring the flow rate of the multiphase flow using the venturi meter.A phase separation phenomenon occurred due to a density difference of water and air inside the venturi,and flow inhomogeneity happened according to the flow velocity difference of each phase.It was confirmed that this flow inhomogeneity increased as the GVF increased due to the uncertainty of the flow measurement.展开更多
It is very important to understand the annular multiphase flow behavior and the effect of hydrate phase transition during deep water drilling. The basic hydrodynamic models, including mass, momentum, and energy conser...It is very important to understand the annular multiphase flow behavior and the effect of hydrate phase transition during deep water drilling. The basic hydrodynamic models, including mass, momentum, and energy conservation equations, were established for annular flow with gas hydrate phase transition during gas kick. The behavior of annular multiphase flow with hydrate phase transition was investigated by analyzing the hydrate-forming region, the gas fraction in the fluid flowing in the annulus, pit gain, bottom hole pressure, and shut-in casing pressure. The simulation shows that it is possible to move the hydrate-forming region away from sea floor by increasing the circulation rate. The decrease in gas volume fraction in the annulus due to hydrate formation reduces pit gain, which can delay the detection of well kick and increase the risk of hydrate plugging in lines. Caution is needed when a well is monitored for gas kick at a relatively low gas production rate, because the possibility of hydrate presence is much greater than that at a relatively high production rate. The shut-in casing pressure cannot reflect the gas kick due to hydrate formation, which increases with time.展开更多
Multiphase flows are ubiquitous in our daily life and engineering applications. It is important to investigate the flow structures to predict their dynamical behaviors ef- fectively. Lagrangian coherent structures （...Multiphase flows are ubiquitous in our daily life and engineering applications. It is important to investigate the flow structures to predict their dynamical behaviors ef- fectively. Lagrangian coherent structures （LCS） defined by the ridges of the finite-time Lyapunov exponent （FTLE） is utilized in this study to elucidate the multiphase interactions in gaseous jets injected into water and time-dependent turbu- lent cavitation under the framework of Navier-Stokes flow computations. For the gaseous jets injected into water, the highlighted phenomena of the jet transportation can be observed by the LCS method, including expansion, bulge, necking/breaking, and back-attack. Besides, the observation of the LCS reveals that the back-attack phenomenon arises from the fact that the injected gas has difficulties to move toward downstream re- gion after the necking/breaking. For the turbulent cavitating flow, the ridge of the FTLE field can form a LCS to capture the front and boundary of the re-entraint jet when the ad- verse pressure gradient is strong enough. It represents a bar- rier between particles trapped inside the circulation region and those moving downstream. The results indicate that the FFLE field has the potential to identify the structures of mul- tiphase flows, and the LCS can capture the interface/barrier or the vortex/circulation region.展开更多
Abrasive water jet cutting technology is widely applied in the materials processing today and attracts great attention from scholars, but many phenomena concerned are not well understood, especially in the internal je...Abrasive water jet cutting technology is widely applied in the materials processing today and attracts great attention from scholars, but many phenomena concerned are not well understood, especially in the internal jet flow of the cutting head at the condition of ultra-high pressure. The multiphase flow in the cutting head is numerically simulated to study the abrasive motion mechanism and wear inside the cutting head at the pressure beyond 300 MPa. Visible predictions of the particles trajectories and wear rate in the cutting head are presented. The influences of the abrasive physical properties, size of the jewel orifice and the operating pressure on the trajectories are discussed. Based on the simulation, a wear experiment is carried out under the corresponding pressures. The simulation and experimental results show that the flow in the mixing chamber is composed of the jet core zone and the disturbance zone, both affect the particles trajectories. The mixing efficiency drops with the increase of the abrasive granularity. The abrasive density determines the response of particles to the effects of different flow zones, the abrasive with medium density gives the best general performance. Increasing the operating pressure or using the jewel with a smaller orifice improves the coherency of p articles trajectories but increases the wear rate of the jewel holder at the same time. Walls of the jewel holder, the entrance of the mixing chamber and the convergence part of the mixing tube are subject to wear out. The computational and experimental results give a qualitative consistency which proves that this numerical method can provide a reliable and visible cognition of the flow characteristics of ultra-high pressure abrasive water jet. The investigation is benefit for improving the machining properties of water jet cutting systems and the optimization design of the cutting head.展开更多
Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronauti...Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronautical and astronautical,transportation,hydraulic and nuclear engineering. In this paper,advances and re-search needs in fundamental studies of dispersed multiphase flows,including the particle/droplet/bubble dynamics,particle-particle,droplet-droplet and bubble-bubble interactions,gas-particle and bubble-liquid turbulence interac-tions,particle-wall interaction,numerical simulation of dispersed multiphase flows,including Reynolds-averaged modeling(RANS modeling),large-eddy simulation(LES) and direct numerical simulation(DNS) are reviewed. The research results obtained by the present author are also included in this review.展开更多
In this paper, a pseudopotential-based multiplerelaxation-time lattice Boltzmann model is proposed for multicomponent/multiphase flow systems. Unlike previous models in the literature, the present model not only enabl...In this paper, a pseudopotential-based multiplerelaxation-time lattice Boltzmann model is proposed for multicomponent/multiphase flow systems. Unlike previous models in the literature, the present model not only enables the study of multicomponent flows with different molecular weights, different viscosities and different Schmidt numbers, but also ensures that the distribution function of each component evolves on the same square lattice without invoking ad- ditional interpolations. Furthermore, the Chapman-Enskog analysis shows that the present model results in the correct hydrodynamic equations, and satisfies the indifferentiability principle. The numerical validation exercises further demonstrate that the favorable performance of the present model.展开更多
Volumetric fraction distribution measurement is a constituent part of process tomography system in oil-water-gas multiphase flow. With the technological development of nuclear radial inspection, dual-energy γ-ray tec...Volumetric fraction distribution measurement is a constituent part of process tomography system in oil-water-gas multiphase flow. With the technological development of nuclear radial inspection, dual-energy γ-ray techniques make it possible to investigate the concentration of the different components on the cross-section of oil-water-gas multiphase pipe-flow. The dual-energy gamma-ray technique is based on materials attenuation coefficients measurement comprised of two radioactive isotopes of 241Am and 241Cs which have emission energies at 59.5 keV and 662 keV in this project. Nuclear instruments and data acquisition system were designed to measure the material’s attenuation dose rate and a number of static tests were conducted at the Multiphase Laboratory, Institute of Mechanics, Chinese Academy of Sciences. Three phases of oil-water-gas media were inves- tigated for their possible use to simulate different media volumetric fraction distributions in experimental vessels. Attenuation intensities were measured, and the arithmetic of linear attenuation coefficients and the equations of volumetric fractions were studied. Investigation of an unexpected measurement error from attenuation equations revealed that a modified arithmetic was involved and finally the system achieved acceptable accuracy in experimental research.展开更多
The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectan...The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.展开更多
The hydrate formation or dissociation in deep subsea flow lines is a challenging problem in oil and gas transport systems. The study of multiphase flows is complex while necessary due to the phase changes （i.e., liqu...The hydrate formation or dissociation in deep subsea flow lines is a challenging problem in oil and gas transport systems. The study of multiphase flows is complex while necessary due to the phase changes （i.e., liquid, solid, and gas） that occur with increasing the temperature and decreasing the pressure. A one-dimensional multiphase flow model coupled with a transient hydrate kinetic model is developed to study the characteristics of the multiphase flows for the hydrates formed by the phase changes in the pipes. The multiphase flow model is derived from a multi-fluid model, while has been widely used in modelling multiphase flows. The heat convection between the fluid and the ambient through the pipe wall is considered in the energy balance equation. The developed multiphase flow model is used to simulate the procedure of the hydrate transport. The results show that the formation of the hydrates can cause hold-up oscillations of water and gas.展开更多
Most multiphase flow separation detection methods used commonly in oilfields are low in efficiency and accuracy,and have data delay.An online multiphase flow detection method is proposed based on magnetic resonance te...Most multiphase flow separation detection methods used commonly in oilfields are low in efficiency and accuracy,and have data delay.An online multiphase flow detection method is proposed based on magnetic resonance technology,and its supporting device has been made and tested in lab and field.The detection technology works in two parts:measure phase holdup in static state and measure flow rate in flowing state.Oil-water ratio is first measured and then gas holdup.The device is composed of a segmented magnet structure and a dual antenna structure for measuring flowing fluid.A highly compact magnetic resonance spectrometer system and intelligent software are developed.Lab experiments and field application show that the online detection system has the following merits:it can measure flow rate and phase holdup only based on magnetic resonance technology;it can detect in-place transient fluid production at high frequency and thus monitor transient fluid production in real time;it can detect oil,gas and water in a full range at high precision,the detection isn’t affected by salinity and emulsification.It is a green,safe and energy-saving system.展开更多
The structure of multiphase flow helico-axial pump＇s rotor and how to model the rotor, especially the blades of the rotor, based on the Solidworks software. More important, the principle of the blade design is mainly...The structure of multiphase flow helico-axial pump＇s rotor and how to model the rotor, especially the blades of the rotor, based on the Solidworks software. More important, the principle of the blade design is mainly introduced. Under the guide of the principle, the 3D coordinates of the blade data points can be got by matlab programming. In the paper, the design step and the modeling step are particularly described through a concrete example.展开更多
The inner flow field of a biogas plant can be optimized by agitating the feedstock to be evenly distributed for a rising biogas production rate. A hydraulic agitator can be installed in the digester with outlets far a...The inner flow field of a biogas plant can be optimized by agitating the feedstock to be evenly distributed for a rising biogas production rate. A hydraulic agitator can be installed in the digester with outlets far above the bottom. Hydraulic mixing is essential in a solid-liquid two-phase flow process, in which large solid particles can be found at the initial stage and turn to being high-concentration viscous liquid （non-Newtonian fluid）. A 0.75 m3 digester was taken as a case study with CFD （computational fluid dynamics） software. The basic pattern was simulated by using water as the medium and the pattern of pseudo plastic fluid state was simulated by the Euler-Euler Model, then the effect of optimized design with bottom inflow and high dispersed outlets could be verified. Viewed from the mixing effects, the velocity of 0.6 m/s is better than l m/s for water medium, while 1 m/s better than 0.6 m/s for pseudo plastic fluid medium.展开更多
When new types of hydrate chemical inhibitor （such as hydrate anti-agglomerant） are used in offshore fields, some difficulties will be encountered in the pipe transportation processing design. In view of these diffi...When new types of hydrate chemical inhibitor （such as hydrate anti-agglomerant） are used in offshore fields, some difficulties will be encountered in the pipe transportation processing design. In view of these difficulties, a horizontal flow experiment pipe loop of transparent polypropylene pipe （25.4 mm inner diameter, 20 m long） was constructed, and some experiments were conducted. The working fluids were the mixed paraffin hydrocarbons, water and condensate oil mixtures obtained from offshore oil field. Analyzing the experimental results and adopting relevant theories of liquid-solid two-phase flow, the computational method of hydrate slurry transportation in pipe was developed. For the operating conditions of a certain multiphase pipeline, by using OLGA software, this method can be used to calculate pressure drop of the multiphase flow pipelines using anti-agglomerant, which can provide support for hydrate anti-agglomerant application.展开更多
基金financially supported by the General Program Grant from the National Natural Science Foundation of China(52274051 and 52174045)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-01)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(51521063)。
文摘Surfactants are widely used in the fracturing fluid to enhance the imbibition and thus the oil recovery rate. However, current numerical models cannot capture the physics behind capillary imbibition during the wettability alteration by surfactants. Although the interacting capillary bundle(ICB) model shows potential in characterizing imbibition rates in different pores during wettability alteration, the existing ICB models neglect the influence of wettability and viscosity ratio on the imbibition behavior, making it difficult to accurately describe the oil-water imbibition behavior within the porous media. In this work,a new ICB mathematical model is established by introducing pressure balance without assuming the position of the leading front to comprehensively describe the imbibition behavior in a porous medium under different conditions, including gas-liquid spontaneous imbibition and oil-water imbibition.When the pore size distribution of a tight rock is known, this new model can predict the changes of water saturation during the displacement process in the tight rock, and also determine the imbibition rate in pores of different sizes. The water saturation profiles obtained from the new model are validated against the waterflooding simulation results from the CMG, while the imbibition rates calculated by the model are validated against the experimental observations of gas-liquid spontaneous imbibition. The good match above indicates the newly proposed model can show the water saturation profile at a macroscopic scale while capture the underlying physics of the multiphase flow in a porous medium at a microscopic scale. Simulation results obtained from this model indicate that both wettability and viscosity ratio can affect the sequence of fluid imbibition into pores of different sizes during the multiphase flow, where less-viscous wetting fluid is preferentially imbibed into larger pores while more-viscous wetting fluid tends to be imbibed into smaller pores. Furthermore, this model provides an avenue to calculate the imbibition rate in pores of different sizes during wettability alteration and capture the non-Darcy effect in micro-and nano-scale pores.
基金the National Key R&D Program(No.2023YFB3709900)the National Natural Science Foundation of China(Nos.U22A20171 and 52104343)the High Steel Central(HSC)at North China University of Science and Technology and Yanshan Univ ersity,China。
文摘A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.
基金Supported by National Natural Science Foundation of China(52104049)Young Elite Scientist Sponsorship Program by BAST(BYESS2023262)Science Foundation of China University of Petroleum,Beijing(2462022BJRC004).
文摘Considering the phase behaviors in condensate gas reservoirs and the oil-gas two-phase linear flow and boundary-dominated flow in the reservoir,a method for predicting the relationship between oil saturation and pressure in the full-path of tight condensate gas well is proposed,and a model for predicting the transient production from tight condensate gas wells with multiphase flow is established.The research indicates that the relationship curve between condensate oil saturation and pressure is crucial for calculating the pseudo-pressure.In the early stage of production or in areas far from the wellbore with high reservoir pressure,the condensate oil saturation can be calculated using early-stage production dynamic data through material balance models.In the late stage of production or in areas close to the wellbore with low reservoir pressure,the condensate oil saturation can be calculated using the data of constant composition expansion test.In the middle stages of production or when reservoir pressure is at an intermediate level,the data obtained from the previous two stages can be interpolated to form a complete full-path relationship curve between oil saturation and pressure.Through simulation and field application,the new method is verified to be reliable and practical.It can be applied for prediction of middle-stage and late-stage production of tight condensate gas wells and assessment of single-well recoverable reserves.
基金The authors would like to acknowledge financial support provided by National Natural Science Foundation of China(No.52074338)We are also grateful to the support of the National Key R&D Program of China(No.2019YFA0708700)+1 种基金National Key Basic Research Program of China(20CX06071A)Bin Yuan would like to thank for the support of Shandong Mountain Tai Scholar Program.Chris Clarkson would like to acknowledge funding support from an NSERC Alliance grant(ALLRP 548652-19)for research related to the topic of this paper.
文摘Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information when the reservoir exhibits reservoir heterogeneity and multiphase flow simultaneously.This work proposes a practical-yet-rigorous method to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,and improve the evaluation of reservoir/fracture properties.A new,general,semi-analytical model is proposed that explicitly accounts for multiphase flow,fractalbased reservoir heterogeneity,anomalous diffusion,and pressure-dependent fluid properties.This is achieved by introducing a new Boltzmann-type transformation,the exponent of which includes reservoir heterogeneity and anomalous diffusion.In order to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,the modified Boltzmann variable allows the conversion of three partial differential equations(PDE's)(i.e.,oil,gas and water diffusion equations)into ordinary differential equations(ODE's)that are easily solved using the Runge-Kutta(RK)method.A modified time-power-law plot is also proposed to estimate the reservoir and fracture properties,recognizing that the classical square-root-of-time-plot is no longer valid when various reservoir complexities are exhibited simultaneously.Using the slope of the straight line on the modified time-power-law plot,the linear flow parameter can be estimated with more confidence.Moreover,because of the new Boltzmann-type transformation,reservoir and fracture properties can be derived more efficiently without the need for defining complex pseudo-variable transformations.Using the new semi-analytical model,the effects of multiphase flow,reservoir heterogeneity and anomalous diffusion on rate-decline behavior are evaluated.For the case of approximately constant flowing pressure,multiphase flow impacts initial oil rate,which is a function of oil relative permeability and well flowing pressure.However,multiphase flow has a minor effect on the oil production decline exponent.Reservoir heterogeneity/anomalous diffusion affect both the initial oil production rate and production decline exponent.The production decline exponent constant is a function of reservoir heterogeneity/anomalous diffusion only.The practical significance of this work is the advancement of RTA techniques to allow for more complex reservoir scenarios,leading to more accurate production forecasting and better-informed capital planning.
基金financially supported by the National Natural Science Foundation of China(No.51704062)the Fundamental Research Funds for the Central Universities,China(No.N2025019)。
文摘A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.
基金supported by the National Natural Science Foundation of China(Grant No.52274061&52004039&51974037)China Postdoctoral Science Foundation(Grant No.2023T160717&2021M693908)+1 种基金The major project of universities affiliated to Jiangsu Province basic science(natural science)research(Grant No.21KJA440001)Jiangsu Qinglan Project,Changzhou Longcheng Talent Plan-Youth Science and Technology Talent Recruitment Project.
文摘The research on the multiphase flow characteristics of hydrate slurry is the key to implementing the risk prevention and control technology of hydrate slurry in deep-water oil and gas mixed transportation system.This paper established a geometric model based on the high-pressure hydrate slurry experimental loop.The model was used to carry out simulation research on the flow characteristics of gas-liquid-solid three-phase flow.The specific research is as follows:Firstly,the effects of factors such as slurry flow velocity,hydrate particle density,hydrate particle size,and hydrate volume fraction on the stratified smooth flow were specifically studied.Orthogonal test obtained particle size has the most influence on the particle concentration distribution.The slurry flow velocity is gradually increased based on stratified smooth flow.Various flow patterns were observed and their characteristics were analyzed.Secondly,increasing the slurry velocity to 2 m/s could achieve the slurry flow pattern of partial hydrate in the pipeline transition from stratified smooth flow to wavy flow.When the flow rate increases to 3 m/s,a violent wave forms throughout the entire loop.Based on wave flow,as the velocity increased to 4 m/s,and the flow pattern changed to slug flow.When the particle concentration was below 10%,the increase of the concentration would aggravate the slug flow trend;if the particle concentration was above 10%,the increase of the concentration would weaken the slug flow trend,the increase of particle density and liquid viscosity would weaken the tendency of slug flow.The relationship between the pressure drop gradients of several different flow patterns is:slug flow>wave flow>stratified smooth flow.
基金Project(07JJ4016) supported by the Natural Science Foundation of Hunan Procvince,China
文摘In order to check the validity of the mathematical model for analyzing the flow field in the air-agitated seed precipitation tank,a scaled down experimental apparatus was designed and the colored tracer and KCl tracer were added in the apparatus to follow the real flow line.Virtue tracers were considered in the mathematical model and the algorithm of tracers was built.The comparison of the results between the experiment and numerical calculation shows that the time of the tracer flows out of stirring tube are 40 s in the experiment and 42 s in numerical calculated result.The transient diffusion process and the solution residence time of the numerical calculation are in good agreement with the experimental results,which indicates that the mathematical model is reliable and can be used to predict the flow field of the air-agitated seed precipitation tank.
基金supported by the Industrial Infrastructure Program through The Korea Institute for Advancement of Technology(KIAT) Grant funded by the Korea government Ministry of Trade,Industry and Energy(Grant N0000502)
文摘The venturi meter has an advantage in its use,because it can measure flow without being much affected by the type of the measured fluid or flow conditions.Hence,it has excellent versatility and is being widely applied in many industries.The flow of a liquid containing air is a representative example of a multiphase flow and exhibits complex flow characteristics.In particular,the greater the gas volume fraction（GVF）,the more inhomogeneous the flow becomes.As a result,using a venturi meter to measure the rate of a flow that has a high GVF generates an error.In this study,the cause of the error occurred in measuring the flow rate for the multiphase flow when using the venturi meter for analysis by CFD.To ensure the reliability of this study,the accuracy of the multiphase flow models for numerical analysis was verified through comparison between the calculated results of numerical analysis and the experimental data.As a result,the Grace model,which is a multiphase flow model established by an experiment with water and air,was confirmed to have the highest reliability.Finally,the characteristics of the internal flow Held about the multiphase flow analysis result generated by applying the Grace model were analyzed to find the cause of the uncertainty occurring when measuring the flow rate of the multiphase flow using the venturi meter.A phase separation phenomenon occurred due to a density difference of water and air inside the venturi,and flow inhomogeneity happened according to the flow velocity difference of each phase.It was confirmed that this flow inhomogeneity increased as the GVF increased due to the uncertainty of the flow measurement.
基金supported by the China National 863 Program (Grant No.2006AA09A106)the Doctoral Program of Higher Education of China (Grant No.20060425502)+1 种基金the National Natural Science Foundation of China (Grant No.50874116)Shandong Province Natural Science Foundation(Grant No.Z2007A01)
文摘It is very important to understand the annular multiphase flow behavior and the effect of hydrate phase transition during deep water drilling. The basic hydrodynamic models, including mass, momentum, and energy conservation equations, were established for annular flow with gas hydrate phase transition during gas kick. The behavior of annular multiphase flow with hydrate phase transition was investigated by analyzing the hydrate-forming region, the gas fraction in the fluid flowing in the annulus, pit gain, bottom hole pressure, and shut-in casing pressure. The simulation shows that it is possible to move the hydrate-forming region away from sea floor by increasing the circulation rate. The decrease in gas volume fraction in the annulus due to hydrate formation reduces pit gain, which can delay the detection of well kick and increase the risk of hydrate plugging in lines. Caution is needed when a well is monitored for gas kick at a relatively low gas production rate, because the possibility of hydrate presence is much greater than that at a relatively high production rate. The shut-in casing pressure cannot reflect the gas kick due to hydrate formation, which increases with time.
文摘Multiphase flows are ubiquitous in our daily life and engineering applications. It is important to investigate the flow structures to predict their dynamical behaviors ef- fectively. Lagrangian coherent structures （LCS） defined by the ridges of the finite-time Lyapunov exponent （FTLE） is utilized in this study to elucidate the multiphase interactions in gaseous jets injected into water and time-dependent turbu- lent cavitation under the framework of Navier-Stokes flow computations. For the gaseous jets injected into water, the highlighted phenomena of the jet transportation can be observed by the LCS method, including expansion, bulge, necking/breaking, and back-attack. Besides, the observation of the LCS reveals that the back-attack phenomenon arises from the fact that the injected gas has difficulties to move toward downstream re- gion after the necking/breaking. For the turbulent cavitating flow, the ridge of the FTLE field can form a LCS to capture the front and boundary of the re-entraint jet when the ad- verse pressure gradient is strong enough. It represents a bar- rier between particles trapped inside the circulation region and those moving downstream. The results indicate that the FFLE field has the potential to identify the structures of mul- tiphase flows, and the LCS can capture the interface/barrier or the vortex/circulation region.
基金supported by National Natural Science Foundation of China (Grant No. 50806031)
文摘Abrasive water jet cutting technology is widely applied in the materials processing today and attracts great attention from scholars, but many phenomena concerned are not well understood, especially in the internal jet flow of the cutting head at the condition of ultra-high pressure. The multiphase flow in the cutting head is numerically simulated to study the abrasive motion mechanism and wear inside the cutting head at the pressure beyond 300 MPa. Visible predictions of the particles trajectories and wear rate in the cutting head are presented. The influences of the abrasive physical properties, size of the jewel orifice and the operating pressure on the trajectories are discussed. Based on the simulation, a wear experiment is carried out under the corresponding pressures. The simulation and experimental results show that the flow in the mixing chamber is composed of the jet core zone and the disturbance zone, both affect the particles trajectories. The mixing efficiency drops with the increase of the abrasive granularity. The abrasive density determines the response of particles to the effects of different flow zones, the abrasive with medium density gives the best general performance. Increasing the operating pressure or using the jewel with a smaller orifice improves the coherency of p articles trajectories but increases the wear rate of the jewel holder at the same time. Walls of the jewel holder, the entrance of the mixing chamber and the convergence part of the mixing tube are subject to wear out. The computational and experimental results give a qualitative consistency which proves that this numerical method can provide a reliable and visible cognition of the flow characteristics of ultra-high pressure abrasive water jet. The investigation is benefit for improving the machining properties of water jet cutting systems and the optimization design of the cutting head.
基金Supported by the Key Projects of National Natural Science Foundation of China （50736006 9587003-13） the State Key Development Program for Basic Research of China （G1999-0222-08） the National Pandeng Project of China （85-06-1-2）
文摘Dispersed multiphase flows,including gas-particle(gas-solid),gas-spray,liquid-particle(liquid-solid) ,liquid-bubble,and bubble-liquid-particle flows,are widely encountered in power,chemical and metallurgical,aeronautical and astronautical,transportation,hydraulic and nuclear engineering. In this paper,advances and re-search needs in fundamental studies of dispersed multiphase flows,including the particle/droplet/bubble dynamics,particle-particle,droplet-droplet and bubble-bubble interactions,gas-particle and bubble-liquid turbulence interac-tions,particle-wall interaction,numerical simulation of dispersed multiphase flows,including Reynolds-averaged modeling(RANS modeling),large-eddy simulation(LES) and direct numerical simulation(DNS) are reviewed. The research results obtained by the present author are also included in this review.
基金supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (62311)supported by the National Natural Science Foundation of China(51006040)+1 种基金the Hong Kong Scholar Programthe National Science Fund for Distinguished Young Scholars of China (51125024)
文摘In this paper, a pseudopotential-based multiplerelaxation-time lattice Boltzmann model is proposed for multicomponent/multiphase flow systems. Unlike previous models in the literature, the present model not only enables the study of multicomponent flows with different molecular weights, different viscosities and different Schmidt numbers, but also ensures that the distribution function of each component evolves on the same square lattice without invoking ad- ditional interpolations. Furthermore, the Chapman-Enskog analysis shows that the present model results in the correct hydrodynamic equations, and satisfies the indifferentiability principle. The numerical validation exercises further demonstrate that the favorable performance of the present model.
文摘Volumetric fraction distribution measurement is a constituent part of process tomography system in oil-water-gas multiphase flow. With the technological development of nuclear radial inspection, dual-energy γ-ray techniques make it possible to investigate the concentration of the different components on the cross-section of oil-water-gas multiphase pipe-flow. The dual-energy gamma-ray technique is based on materials attenuation coefficients measurement comprised of two radioactive isotopes of 241Am and 241Cs which have emission energies at 59.5 keV and 662 keV in this project. Nuclear instruments and data acquisition system were designed to measure the material’s attenuation dose rate and a number of static tests were conducted at the Multiphase Laboratory, Institute of Mechanics, Chinese Academy of Sciences. Three phases of oil-water-gas media were inves- tigated for their possible use to simulate different media volumetric fraction distributions in experimental vessels. Attenuation intensities were measured, and the arithmetic of linear attenuation coefficients and the equations of volumetric fractions were studied. Investigation of an unexpected measurement error from attenuation equations revealed that a modified arithmetic was involved and finally the system achieved acceptable accuracy in experimental research.
基金Project supported by the National Natural Science Foundation of China(No.11425210)the Fundamental Research Funds for the Central Universities(No.WK2090050025)
文摘The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface （DI） model, particularly in the com- plex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line （MCL） when a geometric formulation of the MCL model is included. Therefore, combination of the etch- ing method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical exper- iments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is vali- dated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.
基金Project supported by the National Natural Science Foundation of China(Nos.11272198 and11572183)
文摘The hydrate formation or dissociation in deep subsea flow lines is a challenging problem in oil and gas transport systems. The study of multiphase flows is complex while necessary due to the phase changes （i.e., liquid, solid, and gas） that occur with increasing the temperature and decreasing the pressure. A one-dimensional multiphase flow model coupled with a transient hydrate kinetic model is developed to study the characteristics of the multiphase flows for the hydrates formed by the phase changes in the pipes. The multiphase flow model is derived from a multi-fluid model, while has been widely used in modelling multiphase flows. The heat convection between the fluid and the ambient through the pipe wall is considered in the energy balance equation. The developed multiphase flow model is used to simulate the procedure of the hydrate transport. The results show that the formation of the hydrates can cause hold-up oscillations of water and gas.
基金Supported by the National Natural Science Foundation of China(51704327)
文摘Most multiphase flow separation detection methods used commonly in oilfields are low in efficiency and accuracy,and have data delay.An online multiphase flow detection method is proposed based on magnetic resonance technology,and its supporting device has been made and tested in lab and field.The detection technology works in two parts:measure phase holdup in static state and measure flow rate in flowing state.Oil-water ratio is first measured and then gas holdup.The device is composed of a segmented magnet structure and a dual antenna structure for measuring flowing fluid.A highly compact magnetic resonance spectrometer system and intelligent software are developed.Lab experiments and field application show that the online detection system has the following merits:it can measure flow rate and phase holdup only based on magnetic resonance technology;it can detect in-place transient fluid production at high frequency and thus monitor transient fluid production in real time;it can detect oil,gas and water in a full range at high precision,the detection isn’t affected by salinity and emulsification.It is a green,safe and energy-saving system.
基金Supported by National "863" Project of China (2007AA09Z318)
文摘The structure of multiphase flow helico-axial pump＇s rotor and how to model the rotor, especially the blades of the rotor, based on the Solidworks software. More important, the principle of the blade design is mainly introduced. Under the guide of the principle, the 3D coordinates of the blade data points can be got by matlab programming. In the paper, the design step and the modeling step are particularly described through a concrete example.
文摘The inner flow field of a biogas plant can be optimized by agitating the feedstock to be evenly distributed for a rising biogas production rate. A hydraulic agitator can be installed in the digester with outlets far above the bottom. Hydraulic mixing is essential in a solid-liquid two-phase flow process, in which large solid particles can be found at the initial stage and turn to being high-concentration viscous liquid （non-Newtonian fluid）. A 0.75 m3 digester was taken as a case study with CFD （computational fluid dynamics） software. The basic pattern was simulated by using water as the medium and the pattern of pseudo plastic fluid state was simulated by the Euler-Euler Model, then the effect of optimized design with bottom inflow and high dispersed outlets could be verified. Viewed from the mixing effects, the velocity of 0.6 m/s is better than l m/s for water medium, while 1 m/s better than 0.6 m/s for pseudo plastic fluid medium.
文摘When new types of hydrate chemical inhibitor （such as hydrate anti-agglomerant） are used in offshore fields, some difficulties will be encountered in the pipe transportation processing design. In view of these difficulties, a horizontal flow experiment pipe loop of transparent polypropylene pipe （25.4 mm inner diameter, 20 m long） was constructed, and some experiments were conducted. The working fluids were the mixed paraffin hydrocarbons, water and condensate oil mixtures obtained from offshore oil field. Analyzing the experimental results and adopting relevant theories of liquid-solid two-phase flow, the computational method of hydrate slurry transportation in pipe was developed. For the operating conditions of a certain multiphase pipeline, by using OLGA software, this method can be used to calculate pressure drop of the multiphase flow pipelines using anti-agglomerant, which can provide support for hydrate anti-agglomerant application.
