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.展开更多
In the municipal solid waste incineration process,it is difficult to effectively control the gas oxygen content by setting the air flow according to artificial experience.To address this problem,this paper proposes an...In the municipal solid waste incineration process,it is difficult to effectively control the gas oxygen content by setting the air flow according to artificial experience.To address this problem,this paper proposes an optimization control method of gas oxygen content based on model predictive control.First,a stochastic configuration network is utilized to establish a prediction model of gas oxygen content.Second,an improved differential evolution algorithm that is based on parameter adaptive and t-distribution strategy is employed to address the set value of air flow.Finally,model predictive control is combined with the event triggering strategy to reduce the amount of computation and the controller's frequent actions.The experimental results show that the optimization control method proposed in this paper obtains a smaller degree of fluctuation in the air flow set value,which can ensure the tracking control performance of the gas oxygen content while reducing the amount of calculation.展开更多
The multiphase flow characteristic is one of the most concerning problems during solid fluidization exploitation of marine natural gas hydrate reservoirs.In this research,a new transient gas-liquid-solid multiphase fl...The multiphase flow characteristic is one of the most concerning problems during solid fluidization exploitation of marine natural gas hydrate reservoirs.In this research,a new transient gas-liquid-solid multiphase flow model with hydrate phase transition was developed.Meanwhile,this model considered the coupling relationship among convective heat transfer,hydrate dynamic decomposition,and multi-phase flow.The model can simulate the change of flow pattern from solid-liquid to gas-liquid-solid flow,and describe the distribution character of volume fraction of phase,wellbore temperature and pressure,and hydrate decomposition rate during transportation.The simulation results indicate that the hydrate decomposition region in the wellbore gradually expands,but the hydrate decomposition rate gradually decreases during the solid fluidization exploitation of hydrate.When mining time lasts for 4 h,and the bottom hole pressure decreases by about 0.4 MPa.Increasing NaCl concentration in seawater helps expand hydrate decomposition regions and improves the wellbore hydrate decomposition rate.When the Nacl mass fraction in seawater reaches 15%,it will raise the hydrate decomposition regions to the whole wellbore.In addition,the higher the wellhead backpressure,the lower the decomposition area and decomposition rate of hydrate in the wellbore.When wellhead backpressure reaches 2 MPa,the volume fraction of gas near the wellhead will reduce to about 12%.This work is expected to provide a theoretical basis for the development of marine hydrate reservoirs.展开更多
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%.展开更多
Horizontal gas-liquid two-phase flows widely exist in chemical engineering,oil/gas production and other important industrial processes.Slug flow pattern is the main form of horizontal gas-liquid flows and characterize...Horizontal gas-liquid two-phase flows widely exist in chemical engineering,oil/gas production and other important industrial processes.Slug flow pattern is the main form of horizontal gas-liquid flows and characterized by intermittent motion of film region and slug region.This work aims to develop the ultrasonic Doppler method to realize the simultaneous measurement of the velocity profile and liquid film thickness of slug flow.A single-frequency single-channel transducer is adopted in the design of the field-programmable gate array based ultrasonic Doppler system.A multiple echo repetition technology is used to improve the temporal-spatial resolution for the velocity profile.An experiment of horizontal gas-liquid two-phase flow is implemented in an acrylic pipe with an inner diameter of 20 mm.Considering the aerated characteristics of the liquid slug,slug flow is divided into low-aerated slug flow,high-aerated slug flow and pseudo slug flow.The temporal-spatial velocity distributions of the three kinds of slug flows are reconstructed by using the ultrasonic velocity profile measurement.The evolution characteristics of the average velocity profile in slug flows are investigated.A novel method is proposed to derive the liquid film thickness based on the instantaneous velocity profile.The liquid film thickness can be effectively measured by detecting the position and the size of the bubbles nearly below the elongated gas bubble.Compared with the time of flight method,the film thickness measured by the Doppler system shows a higher accuracy as a bubble layer occurs in the film region.The effect of the gas distribution on the film thickness is uncovered in three kinds of slug flows.展开更多
In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM)...In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..展开更多
Dynamic performance on solids flow with water in deviated tubing is essential for the reliability of pump and normal operation of horizontal and directional wells.Compared with coal-water flow in vertical tubing and s...Dynamic performance on solids flow with water in deviated tubing is essential for the reliability of pump and normal operation of horizontal and directional wells.Compared with coal-water flow in vertical tubing and sand-oil flow with high production in deviated tubing,solids deposition with water shows obvious non-symmetric distributions in deviated tubing from simulations and experiments.The mathematical model of two phase flow was developed under coupling conditions of deviated tubing,low flow rate and viscosity based on the kinetic theory of granular flow and first-order discrete scheme.The results show that solid-water stratified flow in deviated tubing can be divided into two zones of suspension bed and the moving bed throughout the flow field.The solid flowing velocity with water is negative and particles slide down at the bottom of moving bed zone.The process of solids flow with water in deviated tubing will produce pressure loss and consume the kinetic energy.The thickness of deposited layer and the flowing velocity of solids flow downward with water at the moving bed zone enhance with the decreased inlet flow rate and the increased particle size,tubing inside diameter(ID)and inclination angle.Solids are easier into suspension from the upper part of moving bed zone to suspension bed zone and more solid particles flow with water towards the tubing outlet with the increase of inlet flowing velocity.The decision is made to reduce the screen width,tubing ID and inclination angle to maintain to be greater than critical deposition velocity in order to prevent solids settling.And it provides the theoretical basis and technical reserves for solid control and offers an effective approach to enhance tubing cleaning in deviated strings.展开更多
Cyclone separators are widely used in industrial applications. The separation efficiency and pressure drop are the most important parameters to evaluate the performance of processing system. In the simulations,the flo...Cyclone separators are widely used in industrial applications. The separation efficiency and pressure drop are the most important parameters to evaluate the performance of processing system. In the simulations,the flow behavior of gas and particles within a square cyclone separator is simulated by means of computational fluid dynamics. The RNG k- ε model and the Reynolds stress model( RSM) are used to model gas turbulence. The flow behavior is examined in the term of tangential velocity components,static pressure and pressure drop contour plots for flow field and solid volume fraction. The effects of the turbulence model and solid volume fraction on the square cyclone are discussed. The results indicate that the pressure drop increases with the increase of solid volume fraction,and increase with the increase of inlet velocities for two turbulence models, moreover,the simulations results are compared with pressure field. For all runs,the RSM model gives a higher pressure drop compared to the RNG k- ε model. The RSM model provides well the forced vortex and free vortex,and captures better the phenomena occurring during intense vortex flow in the presence of walls within cyclone separators.展开更多
Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of d...Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of dry coal under gas adsorption equilibrium,gas flow and gas diffusion within wet coal under the generally non-equilibrium state are often ignored in the process of gas recovery.In this study,an improved apparent permeability model is proposed which accommodates the water and gas adsorption,stress dependence,water film thickness and gas flow regimes.In the process of modeling,the water adsorption is only affected by water content while the gas adsorption is time and water content dependent;based on poroelastic mechanics,the effective fracture aperture and effective pore radius are derived;and then the variation in water film thickness for different pore types under the effect of water content,stress and adsorption swelling are modeled;the flow regimes are considered based on Beskok’s model.Further,after validation with experimental data,the proposed model was applied to numerical simulations to investigate the evolution of permeability-related factors under the effect of different water contents.The gas flow in wet coal under the non-equilibrium state is explicitly revealed.展开更多
Sand production is one of the main obstacles restricting gas extraction efficiency and safety from marine natural gas hydrate(NGH)reservoirs.Particle migration within the NGH reservoir dominates sand production behavi...Sand production is one of the main obstacles restricting gas extraction efficiency and safety from marine natural gas hydrate(NGH)reservoirs.Particle migration within the NGH reservoir dominates sand production behaviors,while their relationships were rarely reported,severely constrains quantitative evaluation of sand production risks.This paper reports the optical observations of solid particle migration and production from micrometer to mesoscopic scales conditioned to gravel packing during depressurization-induced NGH dissociation for the first time.Theoretical evolutionary modes of sand migration are established based on experimental observations,and its implications on field NGH are comprehensively discussed.Five particle migration regimes of local borehole failure,continuous collapse,wormhole expansion,extensive slow deformation,and pore-wall fluidization are proved to occur during depressurization.The types of particle migration regimes and their transmission modes during depressurization are predominantly determined by initial hydrate saturation.In contrast,the depressurization mainly dominates the transmission rate of the particle migration regimes.Furthermore,both the cumulative mass and the medium grain size of the produced sand decrease linearly with increasing initial methane hydrate(MH)saturation.Discontinuous gas bubble emission,expansion,and explosion during MH dissociation delay sand migration into the wellbore.At the same time,continuous water flow is a requirement for sand production during hydrate dissociation by depressurization.The experiments enlighten us that a constitutive model that can illustrate visible particle migration regimes and their transmission modes is urgently needed to bridge numerical simulation and field applications.Optimizing wellbore layout positions or special reservoir treatment shall be important for mitigating sand production tendency during NGH exploitation.展开更多
A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are s...A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are solved in the framework of a finite element method.The results are validated against those obtained by using the commercial software CMG(Computer Modeling Group software for advanced recovery process simulation).It is shown that the proposed method is reliable.It can capture the fracture rejection characteristics of tight gas reservoirs better than the CMG.A sensitivity analysis of various control factors(initial water saturation,reservoir parameters,and fracturing parameters)affecting the production in tight gas wells is conducted accordingly.Finally,a series of theoretical arguments are provided for a rational and effective development/exploitation of tight sandstone gas reservoirs.展开更多
The flow regimes of GLCC with horizon inlet and a vertical pipe are investigated in experiments,and the velocities and pressure drops data labeled by the corresponding flow regimes are collected.Combined with the flow...The flow regimes of GLCC with horizon inlet and a vertical pipe are investigated in experiments,and the velocities and pressure drops data labeled by the corresponding flow regimes are collected.Combined with the flow regimes data of other GLCC positions from other literatures in existence,the gas and liquid superficial velocities and pressure drops are used as the input of the machine learning algorithms respectively which are applied to identify the flow regimes.The choosing of input data types takes the availability of data for practical industry fields into consideration,and the twelve machine learning algorithms are chosen from the classical and popular algorithms in the area of classification,including the typical ensemble models,SVM,KNN,Bayesian Model and MLP.The results of flow regimes identification show that gas and liquid superficial velocities are the ideal type of input data for the flow regimes identification by machine learning.Most of the ensemble models can identify the flow regimes of GLCC by gas and liquid velocities with the accuracy of 0.99 and more.For the pressure drops as the input of each algorithm,it is not the suitable as gas and liquid velocities,and only XGBoost and Bagging Tree can identify the GLCC flow regimes accurately.The success and confusion of each algorithm are analyzed and explained based on the experimental phenomena of flow regimes evolution processes,the flow regimes map,and the principles of algorithms.The applicability and feasibility of each algorithm according to different types of data for GLCC flow regimes identification are proposed.展开更多
This article presents simulation results on the effects of neutral gas flow for nanoparticle transport in atmospheric-pressure,radio-frequency,capacitively-coupled,and acetylene discharge.The acetylene gas is set to f...This article presents simulation results on the effects of neutral gas flow for nanoparticle transport in atmospheric-pressure,radio-frequency,capacitively-coupled,and acetylene discharge.The acetylene gas is set to flow into the chamber from the upper showerhead electrode.The internal energy of the gas medium therein is transferred into kinetic energy so the gas advection can be triggered.This is represented by the pressure volume work term of the gas energy converse equation.The gas advection leads to the gas temperature sink at the gas inlet,hence a large gas temperature gradient is formed.The thermophoresis relies on the gas temperature gradient,and causes the profile of nanoparticle density to vary from a double-peak structure to a single-peak one.The gas advection influences the properties of electron density and temperature as well and causes the drift-ambipolar mode profile of electron density asymmetric.In the bulk region,i.e.away from the inlet,the gas advection is more like one isovolumetric compression,which slightly increases the temperature of the gas medium at consuming its kinetic energy.展开更多
Methanation is an effective way to efficiently utilize product gas generated from the pyrolysis and gasification of organic solid wastes.To deeply study the heat transfer and mass transfer mechanisms in the reactor,a ...Methanation is an effective way to efficiently utilize product gas generated from the pyrolysis and gasification of organic solid wastes.To deeply study the heat transfer and mass transfer mechanisms in the reactor,a successful three-dimensional comprehensive model has been established.Multiphase flow behavior and heat transfer mechanisms were investigated under reference working conditions.Temperature is determined by the heat release of the reaction and the heat transfer of the gas-solid flow.The maximum temperature can reach 951 K where the catalyst gathers.In the simulation,changes in the gas inlet velocity and catalyst flow rate were made to explore their effects on CO conversion rate and temperature for optimization purposes.As the inlet gas velocity increases from 2.78 to 4.79 m/s,the CO conversion rate decreases from 81.6%to 72.4%.However,more heat is removed from the reactor,and the temperature rise increases from 78.03 to 113.49 K.When the catalyst flow rate is increased from 7.18 to 17.96 kg/(m^(2)·s),the mass of the catalyst in the reactor is increased from 0.0019 to 0.0042 kg,and the CO conversion rate is increased from 66.8%to 81.5%.However,this increases the maximum temperature in the reactor from 940.0 to 966.4 K.展开更多
Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow...Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.展开更多
The study of liquid film characteristics in multiphase flow is a very important research topic, however,the characteristics of the liquid film around Taylor bubble structure in gas, oil and water three-phase flow are ...The study of liquid film characteristics in multiphase flow is a very important research topic, however,the characteristics of the liquid film around Taylor bubble structure in gas, oil and water three-phase flow are not clear. In the present study, a novel liquid film sensor is applied to measure the distributed signals of the liquid film in three-phase flow. Based on the liquid film signals, the liquid film characteristics including the structural characteristics and the nonlinear dynamics characteristics in three-phase flows are investigated for the first time. The structural characteristics including the proportion, the appearance frequency and the thickness of the liquid film are obtained and the influences of the liquid and gas superficial velocities and the oil content on them are investigated. To investigate the nonlinear dynamics characteristics of the liquid film with the changing flow conditions, the entropy analysis is introduced to successfully uncover and quantify the dynamic complexity of the liquid film behavior.展开更多
A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined wit...A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.展开更多
Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,...Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.展开更多
A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient o...A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of 1m and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Be-sides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macro-scopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.展开更多
A three-dimensional model for gas-solid flow in a circulating fluidized bed(CFB) riser was developed based on computational particle fluid dynamics(CPFD).The model was used to simulate the gas-solid flow behavior ...A three-dimensional model for gas-solid flow in a circulating fluidized bed(CFB) riser was developed based on computational particle fluid dynamics(CPFD).The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport(DPT) regime and a fast fluidization(FF) regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(62373017,62073006)and the Beijing Natural Science Foundation of China(4212032)。
文摘In the municipal solid waste incineration process,it is difficult to effectively control the gas oxygen content by setting the air flow according to artificial experience.To address this problem,this paper proposes an optimization control method of gas oxygen content based on model predictive control.First,a stochastic configuration network is utilized to establish a prediction model of gas oxygen content.Second,an improved differential evolution algorithm that is based on parameter adaptive and t-distribution strategy is employed to address the set value of air flow.Finally,model predictive control is combined with the event triggering strategy to reduce the amount of computation and the controller's frequent actions.The experimental results show that the optimization control method proposed in this paper obtains a smaller degree of fluctuation in the air flow set value,which can ensure the tracking control performance of the gas oxygen content while reducing the amount of calculation.
基金supported by the Youth Program of National Natural Science Foundation of China(Grant No.52104012)the Key Program of the National Natural Science Foundation of China(Grant No.51734010)the Key Natural Science Projects of Scientific Research Plan in Colleges and Universities of Xinjiang Uygur Autonomous Region(Grant No.XJEDU20211028).
文摘The multiphase flow characteristic is one of the most concerning problems during solid fluidization exploitation of marine natural gas hydrate reservoirs.In this research,a new transient gas-liquid-solid multiphase flow model with hydrate phase transition was developed.Meanwhile,this model considered the coupling relationship among convective heat transfer,hydrate dynamic decomposition,and multi-phase flow.The model can simulate the change of flow pattern from solid-liquid to gas-liquid-solid flow,and describe the distribution character of volume fraction of phase,wellbore temperature and pressure,and hydrate decomposition rate during transportation.The simulation results indicate that the hydrate decomposition region in the wellbore gradually expands,but the hydrate decomposition rate gradually decreases during the solid fluidization exploitation of hydrate.When mining time lasts for 4 h,and the bottom hole pressure decreases by about 0.4 MPa.Increasing NaCl concentration in seawater helps expand hydrate decomposition regions and improves the wellbore hydrate decomposition rate.When the Nacl mass fraction in seawater reaches 15%,it will raise the hydrate decomposition regions to the whole wellbore.In addition,the higher the wellhead backpressure,the lower the decomposition area and decomposition rate of hydrate in the wellbore.When wellhead backpressure reaches 2 MPa,the volume fraction of gas near the wellhead will reduce to about 12%.This work is expected to provide a theoretical basis for the development of marine hydrate reservoirs.
基金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(41974139,42274148,42074142)。
文摘Horizontal gas-liquid two-phase flows widely exist in chemical engineering,oil/gas production and other important industrial processes.Slug flow pattern is the main form of horizontal gas-liquid flows and characterized by intermittent motion of film region and slug region.This work aims to develop the ultrasonic Doppler method to realize the simultaneous measurement of the velocity profile and liquid film thickness of slug flow.A single-frequency single-channel transducer is adopted in the design of the field-programmable gate array based ultrasonic Doppler system.A multiple echo repetition technology is used to improve the temporal-spatial resolution for the velocity profile.An experiment of horizontal gas-liquid two-phase flow is implemented in an acrylic pipe with an inner diameter of 20 mm.Considering the aerated characteristics of the liquid slug,slug flow is divided into low-aerated slug flow,high-aerated slug flow and pseudo slug flow.The temporal-spatial velocity distributions of the three kinds of slug flows are reconstructed by using the ultrasonic velocity profile measurement.The evolution characteristics of the average velocity profile in slug flows are investigated.A novel method is proposed to derive the liquid film thickness based on the instantaneous velocity profile.The liquid film thickness can be effectively measured by detecting the position and the size of the bubbles nearly below the elongated gas bubble.Compared with the time of flight method,the film thickness measured by the Doppler system shows a higher accuracy as a bubble layer occurs in the film region.The effect of the gas distribution on the film thickness is uncovered in three kinds of slug flows.
基金supported by the National Natural Science Foundation of China(22078009)National Key Research and Development Program of China(2021YFC3001102,2021YFC3001100)。
文摘In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..
基金funded by National Natural Science Foundation of China(Grant No.52074161)National Science and Technology Major Project of China(Grant No.2016ZX05065-001)+2 种基金Taishan Scholar Project of Shandong Province(Grant No.tsqn202211177)Shandong Provincial Plan for Introduction and Cultivation of Young Pioneers in Colleges and Universities(Grant No.2021-QingChuang-30613019)Natural Science Foundation of Shandong Province(Grant No.ZR2022ME173).
文摘Dynamic performance on solids flow with water in deviated tubing is essential for the reliability of pump and normal operation of horizontal and directional wells.Compared with coal-water flow in vertical tubing and sand-oil flow with high production in deviated tubing,solids deposition with water shows obvious non-symmetric distributions in deviated tubing from simulations and experiments.The mathematical model of two phase flow was developed under coupling conditions of deviated tubing,low flow rate and viscosity based on the kinetic theory of granular flow and first-order discrete scheme.The results show that solid-water stratified flow in deviated tubing can be divided into two zones of suspension bed and the moving bed throughout the flow field.The solid flowing velocity with water is negative and particles slide down at the bottom of moving bed zone.The process of solids flow with water in deviated tubing will produce pressure loss and consume the kinetic energy.The thickness of deposited layer and the flowing velocity of solids flow downward with water at the moving bed zone enhance with the decreased inlet flow rate and the increased particle size,tubing inside diameter(ID)and inclination angle.Solids are easier into suspension from the upper part of moving bed zone to suspension bed zone and more solid particles flow with water towards the tubing outlet with the increase of inlet flowing velocity.The decision is made to reduce the screen width,tubing ID and inclination angle to maintain to be greater than critical deposition velocity in order to prevent solids settling.And it provides the theoretical basis and technical reserves for solid control and offers an effective approach to enhance tubing cleaning in deviated strings.
基金Sponsored by the National Natural Science Foundation of China(Grant No.21276056,20490202)
文摘Cyclone separators are widely used in industrial applications. The separation efficiency and pressure drop are the most important parameters to evaluate the performance of processing system. In the simulations,the flow behavior of gas and particles within a square cyclone separator is simulated by means of computational fluid dynamics. The RNG k- ε model and the Reynolds stress model( RSM) are used to model gas turbulence. The flow behavior is examined in the term of tangential velocity components,static pressure and pressure drop contour plots for flow field and solid volume fraction. The effects of the turbulence model and solid volume fraction on the square cyclone are discussed. The results indicate that the pressure drop increases with the increase of solid volume fraction,and increase with the increase of inlet velocities for two turbulence models, moreover,the simulations results are compared with pressure field. For all runs,the RSM model gives a higher pressure drop compared to the RNG k- ε model. The RSM model provides well the forced vortex and free vortex,and captures better the phenomena occurring during intense vortex flow in the presence of walls within cyclone separators.
基金the National Natural Science Foundation of China(No.52079077)the Natural Science Foundation of Shandong Province(No.ZR2021QE069).
文摘Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of dry coal under gas adsorption equilibrium,gas flow and gas diffusion within wet coal under the generally non-equilibrium state are often ignored in the process of gas recovery.In this study,an improved apparent permeability model is proposed which accommodates the water and gas adsorption,stress dependence,water film thickness and gas flow regimes.In the process of modeling,the water adsorption is only affected by water content while the gas adsorption is time and water content dependent;based on poroelastic mechanics,the effective fracture aperture and effective pore radius are derived;and then the variation in water film thickness for different pore types under the effect of water content,stress and adsorption swelling are modeled;the flow regimes are considered based on Beskok’s model.Further,after validation with experimental data,the proposed model was applied to numerical simulations to investigate the evolution of permeability-related factors under the effect of different water contents.The gas flow in wet coal under the non-equilibrium state is explicitly revealed.
基金supported by the Laoshan Laboratory(No.LSKJ LSKJ202203506)the Taishan Scholars Program,and the National Natural Science Foundation of China(Grant No.41976074).
文摘Sand production is one of the main obstacles restricting gas extraction efficiency and safety from marine natural gas hydrate(NGH)reservoirs.Particle migration within the NGH reservoir dominates sand production behaviors,while their relationships were rarely reported,severely constrains quantitative evaluation of sand production risks.This paper reports the optical observations of solid particle migration and production from micrometer to mesoscopic scales conditioned to gravel packing during depressurization-induced NGH dissociation for the first time.Theoretical evolutionary modes of sand migration are established based on experimental observations,and its implications on field NGH are comprehensively discussed.Five particle migration regimes of local borehole failure,continuous collapse,wormhole expansion,extensive slow deformation,and pore-wall fluidization are proved to occur during depressurization.The types of particle migration regimes and their transmission modes during depressurization are predominantly determined by initial hydrate saturation.In contrast,the depressurization mainly dominates the transmission rate of the particle migration regimes.Furthermore,both the cumulative mass and the medium grain size of the produced sand decrease linearly with increasing initial methane hydrate(MH)saturation.Discontinuous gas bubble emission,expansion,and explosion during MH dissociation delay sand migration into the wellbore.At the same time,continuous water flow is a requirement for sand production during hydrate dissociation by depressurization.The experiments enlighten us that a constitutive model that can illustrate visible particle migration regimes and their transmission modes is urgently needed to bridge numerical simulation and field applications.Optimizing wellbore layout positions or special reservoir treatment shall be important for mitigating sand production tendency during NGH exploitation.
基金supported by the China Postdoctoral Science Foundation(2021M702304)and Natural Science Foundation of Shandong Province(ZR2021QE260).
文摘A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are solved in the framework of a finite element method.The results are validated against those obtained by using the commercial software CMG(Computer Modeling Group software for advanced recovery process simulation).It is shown that the proposed method is reliable.It can capture the fracture rejection characteristics of tight gas reservoirs better than the CMG.A sensitivity analysis of various control factors(initial water saturation,reservoir parameters,and fracturing parameters)affecting the production in tight gas wells is conducted accordingly.Finally,a series of theoretical arguments are provided for a rational and effective development/exploitation of tight sandstone gas reservoirs.
文摘The flow regimes of GLCC with horizon inlet and a vertical pipe are investigated in experiments,and the velocities and pressure drops data labeled by the corresponding flow regimes are collected.Combined with the flow regimes data of other GLCC positions from other literatures in existence,the gas and liquid superficial velocities and pressure drops are used as the input of the machine learning algorithms respectively which are applied to identify the flow regimes.The choosing of input data types takes the availability of data for practical industry fields into consideration,and the twelve machine learning algorithms are chosen from the classical and popular algorithms in the area of classification,including the typical ensemble models,SVM,KNN,Bayesian Model and MLP.The results of flow regimes identification show that gas and liquid superficial velocities are the ideal type of input data for the flow regimes identification by machine learning.Most of the ensemble models can identify the flow regimes of GLCC by gas and liquid velocities with the accuracy of 0.99 and more.For the pressure drops as the input of each algorithm,it is not the suitable as gas and liquid velocities,and only XGBoost and Bagging Tree can identify the GLCC flow regimes accurately.The success and confusion of each algorithm are analyzed and explained based on the experimental phenomena of flow regimes evolution processes,the flow regimes map,and the principles of algorithms.The applicability and feasibility of each algorithm according to different types of data for GLCC flow regimes identification are proposed.
基金supported by National Natural Science Foundation of China(Nos.11805107 and 12275039)the Fundamental Research Funds in Heilongjiang Provincial Universities of China(No.135509124)the Graduate Innovation Foundation of Qiqihar University(No.YJSCX2022014)。
文摘This article presents simulation results on the effects of neutral gas flow for nanoparticle transport in atmospheric-pressure,radio-frequency,capacitively-coupled,and acetylene discharge.The acetylene gas is set to flow into the chamber from the upper showerhead electrode.The internal energy of the gas medium therein is transferred into kinetic energy so the gas advection can be triggered.This is represented by the pressure volume work term of the gas energy converse equation.The gas advection leads to the gas temperature sink at the gas inlet,hence a large gas temperature gradient is formed.The thermophoresis relies on the gas temperature gradient,and causes the profile of nanoparticle density to vary from a double-peak structure to a single-peak one.The gas advection influences the properties of electron density and temperature as well and causes the drift-ambipolar mode profile of electron density asymmetric.In the bulk region,i.e.away from the inlet,the gas advection is more like one isovolumetric compression,which slightly increases the temperature of the gas medium at consuming its kinetic energy.
基金supported by the National Key Research and Development Program of China[Grant Number 2019YFC1906802].
文摘Methanation is an effective way to efficiently utilize product gas generated from the pyrolysis and gasification of organic solid wastes.To deeply study the heat transfer and mass transfer mechanisms in the reactor,a successful three-dimensional comprehensive model has been established.Multiphase flow behavior and heat transfer mechanisms were investigated under reference working conditions.Temperature is determined by the heat release of the reaction and the heat transfer of the gas-solid flow.The maximum temperature can reach 951 K where the catalyst gathers.In the simulation,changes in the gas inlet velocity and catalyst flow rate were made to explore their effects on CO conversion rate and temperature for optimization purposes.As the inlet gas velocity increases from 2.78 to 4.79 m/s,the CO conversion rate decreases from 81.6%to 72.4%.However,more heat is removed from the reactor,and the temperature rise increases from 78.03 to 113.49 K.When the catalyst flow rate is increased from 7.18 to 17.96 kg/(m^(2)·s),the mass of the catalyst in the reactor is increased from 0.0019 to 0.0042 kg,and the CO conversion rate is increased from 66.8%to 81.5%.However,this increases the maximum temperature in the reactor from 940.0 to 966.4 K.
基金This research was conducted within the DECOVALEX-2023 projectDECOVALEX is an international research project comprising participants from industry,government,and academia,focusing on development of understanding,models and codes in complex coupled problems in sub-surface geological and engineering applications.DECOVALEX-2023 is the current phase of the project.The authors appreciate the DECOVALEX-2023 Funding Organisations Andra,BASE,BGE,BGR,CAS,CNSC,COVRA,US DOE,ENRESA,ENSI,JAEA,KAERI,NWMO,RWM,SÚRAO,SSM and Taipower for their financial and technical support of the work described in this paper.The statements made in the paper are,however,solely those of the authors and do not necessarily reflect those of the Funding Organisations.This work was further supported by the German Federal Ministry for Economic Affairs and Climate Action(BMWK).
文摘Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.
基金supported by the National Natural Science Foundation of China (42074142, 51527805, 41974139)China Postdoctoral Science Foundation (2020M680969, 2021T140099)the Fundamental Research Funds for the Central Universities (N2104013)。
文摘The study of liquid film characteristics in multiphase flow is a very important research topic, however,the characteristics of the liquid film around Taylor bubble structure in gas, oil and water three-phase flow are not clear. In the present study, a novel liquid film sensor is applied to measure the distributed signals of the liquid film in three-phase flow. Based on the liquid film signals, the liquid film characteristics including the structural characteristics and the nonlinear dynamics characteristics in three-phase flows are investigated for the first time. The structural characteristics including the proportion, the appearance frequency and the thickness of the liquid film are obtained and the influences of the liquid and gas superficial velocities and the oil content on them are investigated. To investigate the nonlinear dynamics characteristics of the liquid film with the changing flow conditions, the entropy analysis is introduced to successfully uncover and quantify the dynamic complexity of the liquid film behavior.
基金Project supported by the National Natural Science Foundation of China(No.12072158)。
文摘A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.
文摘Gas-solid two-phase turbulent flows,mass transfer,heat transfer and catalytic cracking reactions areknown to exert interrelated influences in commercial fluid catalytic cracking(FCC)riser reactors.In the presentpaper,a three-dimensional turbulent gas-solid two-phase flow-reaction model for FCC riser reactors was devel-oped.The model took into account the gas-solid two-phase turbulent flows,inter-phase heat transfer,masstransfer,catalytic cracking reactions and their interrelated influence.The k-V-k_P two-phase turbulence modelwas employed and modified for the two-phase turbulent flow patterns with relatively high particle concentration.Boundary conditions for the flow-reaction model were given.Related numerical algorithm was formed and a nu-merical code was drawn up.Numerical modeling for commercial FCC riser reactors could be carried out with thepresented model.
基金Supported by the National Key Program of Basic Research in China (No.2004CB217702, No.2005CB221202, No.2006CB20030201) and the National Natural Science Foundation of China (No.20590367, No.50676021, No.50606006).
文摘A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of 1m and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Be-sides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macro-scopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.
基金support by the National Basic Research Program (Grant No. 2010CB226906,and 2012CB215000)
文摘A three-dimensional model for gas-solid flow in a circulating fluidized bed(CFB) riser was developed based on computational particle fluid dynamics(CPFD).The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport(DPT) regime and a fast fluidization(FF) regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.