In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.Bu...In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.But such an approach has not been proved rigorously,and its accuracy must be carefully considered.In this paper,we discuss the problem of a sphere moving in a non-uniform flow field,on the basis of the fundamental theory of hydrodynamics.We introduce two assumptions:i)the diameter of the sphere is much smaller than the linear dimension of the flow field,and ii)the velocity of the sphere relative to the surrounding water is very small.Using these two assumptions,we solve the linearized Navier-Stokes equations and equations of continuity by the method of Laplace transform,and finally we obtain a formula for the resistance acting on a sphere moving in a non-uniform flow field.展开更多
In the metallurgical industries, it is very important to characterize the flow of molten metals in open channels given that they are transported through these devices to different plant sections. Howeve...In the metallurgical industries, it is very important to characterize the flow of molten metals in open channels given that they are transported through these devices to different plant sections. However, unlike the flow of water which has been studied since ancient times, the flow of molten metals in open channels has received little attention. The unsteady non-uniform flow of blast furnace molten pig iron in a rectangular open channel is analyzed in this work by numerical solution of the Saint-Venant equations. The influence of mesh size on the convergence of molten metal height is studied to determine the proper mesh and time step sizes. A sinusoidal inflow pulse is imposed at the entrance of the channel in order to analyze the propagation of the resulting wave. The influence of the angle of inclination of the channel and the roughness coefficient of the walls on the amplitude and the dynamic behavior of the height of the molten metal are analyzed. Phase portraits of the channel state variables are constructed and interpreted. Numerical simulations show that as the angle of inclination of the channel increases, the amplitude of the formed wave decreases. From 10 degrees onwards, the peak of the wave descends even below the initial height. On the other hand, the roughness coefficient affects the molten pig iron height profiles in an inverse way than the angle of inclination. The amplitude of the formed wave increases as the roughness coefficient increases.展开更多
Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are...Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are based on the traditional quasi-steady method. In this work, an improved quasi-steady method along with the transient method was presented to compute the rotordynamic coefficients of a long seal. By comparisons with experimental data, the shortcomings of quasi-steady methods have been identified. Then, the effects of non-uniform incoming flow on seal dynamic coefficients were studied by transient simulations. Results indicate that the long seal has large cross stiffness k and direct mass M which are not good for rotor stability, while the transient method is more suitable for the long seal for its excellent performance in predicting M. When the incoming flow is non-uniform, the stiffness coefficients vary with the eccentric directions. Based on the rotordynamic coefficients under uniform incoming flow, the linearized fluid force formulas, which can consider the effects of non-uniform incoming flow, have been presented and can well explain the varying-stiffness phenomenon.展开更多
Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinni...Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.展开更多
The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a...The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.展开更多
A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes...A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its in...In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.展开更多
Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks th...Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.展开更多
In this study,the installation of an airlift pump with inner diameter of 102 mm and length of 5.64 m was utilized to consider the conveying process of non-spherical coal particles with density of 1340 kg/m3 and graini...In this study,the installation of an airlift pump with inner diameter of 102 mm and length of 5.64 m was utilized to consider the conveying process of non-spherical coal particles with density of 1340 kg/m3 and graining 25-44.5 mm.The test results revealed that the magnitude of increase in the solid transport rate due to the changes in the three tested parameters between compressed air velocity,submergence ratio,and feeding coal possibility was not the same,which are stand in range of 20%,75%,and 40%,respectively.Hence,creating the optimal airlift pump performance is highly dependent on submergence ratio.More importantly,we measured the solid volume fraction using the method of one-way valves in order to minimize the disadvantages of conventional devices,such as fast speed camera and conductivity ring sensor.The results confirmed that the volume fraction of the solid phase in the transfer process was always less than 12%.To validate present experimental data,the existing empirical correlations together with the theoretical equations related to the multiphase flow was used.The overall agreement between the theory and experimental solid delivery results was particularly good instead of the first stage of conveying process.This drawback can be corrected by omitting the role of friction and shear stress at low air income velocity.It was also found that the model developed by Kalenik failed to predict the performance of our airlift operation in terms of the mass flow rate of the coal particles.展开更多
A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hyb...A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hybrid Eulerian-Lagrangian approach is used with the rag being modeled as a set of interconnected particles.The influence of various parameters is considered,namely the inlet velocity(1.5,2.0,and 2.5 m/s,respectively),the angle formed by the initially straight rag with the flow direction(45°,60°and 90°,respectively),and the inlet position(90,100,and 110 mm,respectively).The results show that the flow rate has a significant impact on the permeability of the rag.The higher the flow rate,the higher the permeability and the rag speed difference.The angle has a minor effect on rag permeability,with 45°being the most favorable angle for permeability.The inlet position has a small impact on rag permeability,while reducing the initial distance between the rag an the cylinder makes it easier for rags to pass through.展开更多
文摘In hydraulics,when we deal with the problem of sand particles moving relative to the surrounding water,Stokes'formula of resistance has usually been used to render the velocity of sedimentation of the particles.But such an approach has not been proved rigorously,and its accuracy must be carefully considered.In this paper,we discuss the problem of a sphere moving in a non-uniform flow field,on the basis of the fundamental theory of hydrodynamics.We introduce two assumptions:i)the diameter of the sphere is much smaller than the linear dimension of the flow field,and ii)the velocity of the sphere relative to the surrounding water is very small.Using these two assumptions,we solve the linearized Navier-Stokes equations and equations of continuity by the method of Laplace transform,and finally we obtain a formula for the resistance acting on a sphere moving in a non-uniform flow field.
文摘In the metallurgical industries, it is very important to characterize the flow of molten metals in open channels given that they are transported through these devices to different plant sections. However, unlike the flow of water which has been studied since ancient times, the flow of molten metals in open channels has received little attention. The unsteady non-uniform flow of blast furnace molten pig iron in a rectangular open channel is analyzed in this work by numerical solution of the Saint-Venant equations. The influence of mesh size on the convergence of molten metal height is studied to determine the proper mesh and time step sizes. A sinusoidal inflow pulse is imposed at the entrance of the channel in order to analyze the propagation of the resulting wave. The influence of the angle of inclination of the channel and the roughness coefficient of the walls on the amplitude and the dynamic behavior of the height of the molten metal are analyzed. Phase portraits of the channel state variables are constructed and interpreted. Numerical simulations show that as the angle of inclination of the channel increases, the amplitude of the formed wave decreases. From 10 degrees onwards, the peak of the wave descends even below the initial height. On the other hand, the roughness coefficient affects the molten pig iron height profiles in an inverse way than the angle of inclination. The amplitude of the formed wave increases as the roughness coefficient increases.
基金Project(51276213)supported by the National Natural Science Foundation of ChinaProject(2013BAF01B00)supported by the National Science and Technology Support Program of China
文摘Currently, the flow field of annular seals disturbed by the circular whirl motion of rotors is usually solved using computational fluid dynamics(CFD) to evaluate the five rotordynamic coefficients. The simulations are based on the traditional quasi-steady method. In this work, an improved quasi-steady method along with the transient method was presented to compute the rotordynamic coefficients of a long seal. By comparisons with experimental data, the shortcomings of quasi-steady methods have been identified. Then, the effects of non-uniform incoming flow on seal dynamic coefficients were studied by transient simulations. Results indicate that the long seal has large cross stiffness k and direct mass M which are not good for rotor stability, while the transient method is more suitable for the long seal for its excellent performance in predicting M. When the incoming flow is non-uniform, the stiffness coefficients vary with the eccentric directions. Based on the rotordynamic coefficients under uniform incoming flow, the linearized fluid force formulas, which can consider the effects of non-uniform incoming flow, have been presented and can well explain the varying-stiffness phenomenon.
基金This work is supported by the National Natural Science Foundation of China(No.52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(No.BYESS2023262)Science Foundation of China University of Petroleum,Beijing(No.2462022BJRC004).
文摘Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.
基金Funding support from Heilongjiang"Open Competition"project(Grant No.DQYT2022-JS-758)is greatly acknowledgedFinancial support from the National Natural Science Foundation of China(Grant Nos.52304025 and 52174025)is acknowledged+1 种基金supports from Northeast Petroleum University and Guangdong Basic and Applied Basic Research Foundationsupport from the Heilongjiang Touyan Innovation Team Program.
文摘The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.
基金This study was supported by the National Natural Science Foundation of China(U22B2075,52274056,51974356).
文摘A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金supported by the Australian Research Council(Grant No.DP200101293)supported by the UWA-China Joint Scholarships(201906430030).
文摘In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.
基金supported by National Natural Science Foundation(52204050)Sichuan Science and Technology Program(2021ZHCG0013,22ZDYF3009)。
文摘Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.
基金supported by the European Research Council(Research Fund for Coal and Steel)under Grant Agreement number 800757.
文摘In this study,the installation of an airlift pump with inner diameter of 102 mm and length of 5.64 m was utilized to consider the conveying process of non-spherical coal particles with density of 1340 kg/m3 and graining 25-44.5 mm.The test results revealed that the magnitude of increase in the solid transport rate due to the changes in the three tested parameters between compressed air velocity,submergence ratio,and feeding coal possibility was not the same,which are stand in range of 20%,75%,and 40%,respectively.Hence,creating the optimal airlift pump performance is highly dependent on submergence ratio.More importantly,we measured the solid volume fraction using the method of one-way valves in order to minimize the disadvantages of conventional devices,such as fast speed camera and conductivity ring sensor.The results confirmed that the volume fraction of the solid phase in the transfer process was always less than 12%.To validate present experimental data,the existing empirical correlations together with the theoretical equations related to the multiphase flow was used.The overall agreement between the theory and experimental solid delivery results was particularly good instead of the first stage of conveying process.This drawback can be corrected by omitting the role of friction and shear stress at low air income velocity.It was also found that the model developed by Kalenik failed to predict the performance of our airlift operation in terms of the mass flow rate of the coal particles.
基金funded by the Zhejiang Provincial Natural Science Foundation of China(Grant Nos.LY21E060004,LGG22E060011)National Natural Science Foundation of China(Grant No.51976193).
文摘A coupled Computational Fluid Dynamics-Discrete Element Method(CFD-DEM)approach is used to calculate the interaction of a flexible rag transported by a fluid current with a fixed solid cylinder.More specifically a hybrid Eulerian-Lagrangian approach is used with the rag being modeled as a set of interconnected particles.The influence of various parameters is considered,namely the inlet velocity(1.5,2.0,and 2.5 m/s,respectively),the angle formed by the initially straight rag with the flow direction(45°,60°and 90°,respectively),and the inlet position(90,100,and 110 mm,respectively).The results show that the flow rate has a significant impact on the permeability of the rag.The higher the flow rate,the higher the permeability and the rag speed difference.The angle has a minor effect on rag permeability,with 45°being the most favorable angle for permeability.The inlet position has a small impact on rag permeability,while reducing the initial distance between the rag an the cylinder makes it easier for rags to pass through.