Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to ...Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.展开更多
Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typica...Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.展开更多
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.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement m...Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.展开更多
Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fre...Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.展开更多
Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid pa...Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.展开更多
Oil-water two-phase flow is ubiquitous in shale strata due to the existence of connate water and the injection of fracturing fluid.In this work,we propose a relative permeability model based on a modified Hagen-Poiseu...Oil-water two-phase flow is ubiquitous in shale strata due to the existence of connate water and the injection of fracturing fluid.In this work,we propose a relative permeability model based on a modified Hagen-Poiseuille(HP)equation and shale reconstruction algorithm.The proposed model can consider the nanoconfined effects(slip length and spatially varying viscosity),oil-water distribution,pore size distribution(PSD),total organic matter content(TOC),and micro-fracture.The results show that the increasing contact angles of organic matters(OM)and inorganic minerals(iOM)increase the relative permeability of both oil and water.As the viscosity ratio increases,the relative permeability of oil phase increases while that of water phase decreases,due to the different water-oil distribution.The effective permeability of both oil and water decreases with the increasing TOC.However,the relative permeability of water phase increases while that of oil phase decreases.The increasing number and decreasing deviation angle of micro-fracture increase the effective permeability of oil and water.However,microfracture has a minor effect on relative permeability.Our model can help understand oil-water twophase flow in shale reservoirs and provide parameter characterization for reservoir numerical simulation.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
The existence of water phase occupies oil flow area and impacts the confined oil flow behavior at the solid substrate in inorganic nanopores of shale oil reservoirs,resulting in a completely different flow pattern whe...The existence of water phase occupies oil flow area and impacts the confined oil flow behavior at the solid substrate in inorganic nanopores of shale oil reservoirs,resulting in a completely different flow pattern when compared with the single oil phase flow.This study proposes an analytical model to describe the water-oil two-phase flow.In this model,water slippage at the solid substrate is considered while oil slip is introduced to calculate the oil movement at the solid-oil boundary in dry conditions.It is proven that the oil flow profiles of both the two-phase model and single-phase model show parabolic shapes,but the oil flow capacity drops when water takes up the flow space and the impact of water is more significant when the pore dimension is smaller than 30 nm.Also,the oil flow velocity at a pore center is found to drop linearly given a larger water saturation in wet conditions.The effects of surface wettability and oil properties on water-oil flow are also discussed.Compared with the existing singlephase models,this model describes oil flow pattern in the wet condition with the incorporation of the influence of nanopore properties,which better predicts the oil transport in actual reservoir conditions.Water-oil relative permeability curves are also obtained to improve oil yield.展开更多
The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study o...The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.展开更多
Electromagnetic Computer Tomography （ECT） is a method to probe the interior of an inhomogeneous medium via surface measurement in a non-linear way. Due to the great differences in conductivity and permittivity betwe...Electromagnetic Computer Tomography （ECT） is a method to probe the interior of an inhomogeneous medium via surface measurement in a non-linear way. Due to the great differences in conductivity and permittivity between oil and water in the well, Electromagnetic Tomography Well Logging （ETWL）, a new flow imaging measurement system, is proposed to describe the distribution and movement of oil/water two-phase flow in the well by scanning the detected region and applying a suitable data processing algorithm. The results of the numerical simulation and physical modeling show that the system could provide a clear image of the flow profile.展开更多
Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of t...Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of the most complex systems. In this paper, we use complex networks to study the inclined oil water two-phase flow. Two different complex network construction methods are proposed to build two types of networks, i.e. the flow pattern complex network （FPCN） and fluid dynamic complex network （FDCN）. Through detecting the community structure of FPCN by the community-detection algorithm based on K-means clustering, useful and interesting results are found which can be used for identifying three inclined oil-water flow patterns. To investigate the dynamic characteristics of the inclined oil-water two-phase flow, we construct 48 FDCNs under different flow conditions, and find that the power-law exponent and the network information entropy, which are sensitive to the flow pattern transition, can both characterize the nonlinear dynamics of the inclined oil-water two-phase flow. In this paper, from a new perspective, we not only introduce a complex network theory into the study of the oil-water two-phase flow but also indicate that the complex network may be a powerful tool for exploring nonlinear time series in practice.展开更多
The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under th...The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.展开更多
We generate a directed weighted complex network by a method based on Markov transition probability to represent an experimental two-phase flow. We first systematically carry out gas-liquid two-phase flow experiments f...We generate a directed weighted complex network by a method based on Markov transition probability to represent an experimental two-phase flow. We first systematically carry out gas-liquid two-phase flow experiments for measuring the time series of flow signals. Then we construct directed weighted complex networks from various time series in terms of a network generation method based on Markov transition probability. We find that the generated network inherits the main features of the time series in the network structure. In particular, the networks from time series with different dynamics exhibit distinct topological properties. Finally, we construct two-phase flow directed weighted networks from experimental signals and associate the dynamic behavior of gas-liquid two-phase flow with the topological statistics of the generated networks. The results suggest that the topological statistics of two-phase flow networks allow quantitative characterization of the dynamic flow behavior in the transitions among different gas-liquid flow patterns.展开更多
Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleto...Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives（normals and curvatures） and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.展开更多
In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The la...In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace＇s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.展开更多
Two-phase flow in a horizontal pipe was investigated by using numerical and experimental visualization methods.A horizontal pipe was built for qualitative and quantitative flow visualization.The length of horizontal p...Two-phase flow in a horizontal pipe was investigated by using numerical and experimental visualization methods.A horizontal pipe was built for qualitative and quantitative flow visualization.The length of horizontal pipe flow system was 9.5 m and the inner diameter was 51 mm.High-speed video method was used for the qualitative visualization and PIV method was applied for the quantitative visua-lization.The same geometry model was used for the numerical study.Three flow regimes including stratified flow,elongated bubble and slug flow field were generated and visualized by using numerical and experimental methods.The results show that the numerical simulation results are qualitatively si-milar to that of the experimental results.In addition,more quantitative results can be analyzed by numerical method.Development and decay process of slug flow was investigated,showing that the decay of slug heavily depends on the magnitude of nose velocity and its lasting time.It can also be found that the liquid superficial velocity plays a significant role in affecting the slug frequency.When keeping the gas superficial velocity constant,the frequency will increase with the liquid superficial velocity.展开更多
We extend the complexity entropy causality plane（CECP） to propose a multi-scale complexity entropy causality plane（MS-CECP） and further use the proposed method to discriminate the deterministic characteristics of ...We extend the complexity entropy causality plane（CECP） to propose a multi-scale complexity entropy causality plane（MS-CECP） and further use the proposed method to discriminate the deterministic characteristics of different oil-in-water flows. We first take several typical time series for example to investigate the characteristic of the MS-CECP and find that the MS-CECP not only describes the continuous loss of dynamical structure with the increase of scale, but also reflects the determinacy of the system. Then we calculate the MS-CECP for the conductance fluctuating signals measured from oil–water two-phase flow loop test facility. The results indicate that the MS-CECP could be an intrinsic measure for indicating oil-in-water two-phase flow structures.展开更多
We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eu...We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.展开更多
基金the support of the Opening Fund of State Key Laboratory of Multiphase Flow in Power Engineering(SKLMF-KF-2102)。
文摘Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.
基金supported by the National Natural Science Foundation of China Joint Fund Project (Grant/Award Number: U20B6003)National Natural Science Foundation of China (Grant/Award Number: 52304054)。
文摘Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.
基金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.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
基金supported by the National Natural Science Foundation of China(Grant Nos.22275092,52102107 and 52372084)the Fundamental Research Funds for the Central Universities(Grant No.30923010920)。
文摘Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.
文摘Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.
基金supported by the Fund of Innovation Research Group of National Natural Science Foundation of China (Grant NO.5052160450323001)Major Program of National Natural Science Foundation of China (Grant No.50536020)
文摘Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.
基金supported by the National Natural Science Foundation of China(51804328,51974348)
文摘Oil-water two-phase flow is ubiquitous in shale strata due to the existence of connate water and the injection of fracturing fluid.In this work,we propose a relative permeability model based on a modified Hagen-Poiseuille(HP)equation and shale reconstruction algorithm.The proposed model can consider the nanoconfined effects(slip length and spatially varying viscosity),oil-water distribution,pore size distribution(PSD),total organic matter content(TOC),and micro-fracture.The results show that the increasing contact angles of organic matters(OM)and inorganic minerals(iOM)increase the relative permeability of both oil and water.As the viscosity ratio increases,the relative permeability of oil phase increases while that of water phase decreases,due to the different water-oil distribution.The effective permeability of both oil and water decreases with the increasing TOC.However,the relative permeability of water phase increases while that of oil phase decreases.The increasing number and decreasing deviation angle of micro-fracture increase the effective permeability of oil and water.However,microfracture has a minor effect on relative permeability.Our model can help understand oil-water twophase flow in shale reservoirs and provide parameter characterization for reservoir numerical simulation.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
基金the NSERC/Energi Simulation and Alberta Innovates Chairs for providing research fundingthe Beijing Natural Science Foundation(2204093)Science Foundation of China University of Petroleum,Beijing(No.2462018YJRC033)for providing research funding。
文摘The existence of water phase occupies oil flow area and impacts the confined oil flow behavior at the solid substrate in inorganic nanopores of shale oil reservoirs,resulting in a completely different flow pattern when compared with the single oil phase flow.This study proposes an analytical model to describe the water-oil two-phase flow.In this model,water slippage at the solid substrate is considered while oil slip is introduced to calculate the oil movement at the solid-oil boundary in dry conditions.It is proven that the oil flow profiles of both the two-phase model and single-phase model show parabolic shapes,but the oil flow capacity drops when water takes up the flow space and the impact of water is more significant when the pore dimension is smaller than 30 nm.Also,the oil flow velocity at a pore center is found to drop linearly given a larger water saturation in wet conditions.The effects of surface wettability and oil properties on water-oil flow are also discussed.Compared with the existing singlephase models,this model describes oil flow pattern in the wet condition with the incorporation of the influence of nanopore properties,which better predicts the oil transport in actual reservoir conditions.Water-oil relative permeability curves are also obtained to improve oil yield.
基金supported by the National Natural Science Foundation of China(11722104,11671150)supported by the National Natural Science Foundation of China(11571280,11331005)+3 种基金supported by the National Natural Science Foundation of China(11331005,11771150)by GDUPS(2016)the Fundamental Research Funds for the Central Universities of China(D2172260)FANEDD No.201315
文摘The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.
基金This work was supported by the National Natural Science Foundation of China(60472019).
文摘Electromagnetic Computer Tomography （ECT） is a method to probe the interior of an inhomogeneous medium via surface measurement in a non-linear way. Due to the great differences in conductivity and permittivity between oil and water in the well, Electromagnetic Tomography Well Logging （ETWL）, a new flow imaging measurement system, is proposed to describe the distribution and movement of oil/water two-phase flow in the well by scanning the detected region and applying a suitable data processing algorithm. The results of the numerical simulation and physical modeling show that the system could provide a clear image of the flow profile.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 50674070 and 60374041)the National High Technology Research and Development Program of China (Grant No 2007AA06Z231)
文摘Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of the most complex systems. In this paper, we use complex networks to study the inclined oil water two-phase flow. Two different complex network construction methods are proposed to build two types of networks, i.e. the flow pattern complex network （FPCN） and fluid dynamic complex network （FDCN）. Through detecting the community structure of FPCN by the community-detection algorithm based on K-means clustering, useful and interesting results are found which can be used for identifying three inclined oil-water flow patterns. To investigate the dynamic characteristics of the inclined oil-water two-phase flow, we construct 48 FDCNs under different flow conditions, and find that the power-law exponent and the network information entropy, which are sensitive to the flow pattern transition, can both characterize the nonlinear dynamics of the inclined oil-water two-phase flow. In this paper, from a new perspective, we not only introduce a complex network theory into the study of the oil-water two-phase flow but also indicate that the complex network may be a powerful tool for exploring nonlinear time series in practice.
基金financially supported by the National Natural Science Foundation of China (No. 42102149)the Fundamental Research Funds for the Central Universities (No. 2462021YXZZ005)。
文摘The relationship between capillary pressure and saturation plays a critical role in the characterization of two-phase flow and transport in aquifers and oil reservoirs. This relationship is usually determined under the static condition, where capillary pressure is the only function of saturation. However,considerable experiments have suggested that the dependence of capillary pressure on desaturation rate is under the dynamic condition. Thus, a more general description of capillary pressure that includes dynamic capillary effect has been approved widely. A comprehensive understanding of the dynamic capillary effect is needed for the investigation of the two-phase flow in porous media by various methods. In general, dynamic capillary effect in porous media can be studied through the laboratory experiment, pore-to macro-scale modeling, and artificial neural network. Here, main principle and research procedures of each method are reviewed in detail. Then, research progress, disadvantages and advantages are discussed, respectively. In addition, upscaling study from pore-to macro-scale are introduced, which explains the difference between laboratory experiment and pore-scale modeling. At last, several future perspectives and recommendations for optimal solution of dynamic capillary effect are presented.
基金Project supported by the National Natural Science Foundation of China ( Grant Nos. 61104148, 41174109, and 50974095)the National Science and Technology Major Project of the Ministry of Science and Technology of China (Grant No. 2011ZX05020-006)the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20110032120088)
文摘We generate a directed weighted complex network by a method based on Markov transition probability to represent an experimental two-phase flow. We first systematically carry out gas-liquid two-phase flow experiments for measuring the time series of flow signals. Then we construct directed weighted complex networks from various time series in terms of a network generation method based on Markov transition probability. We find that the generated network inherits the main features of the time series in the network structure. In particular, the networks from time series with different dynamics exhibit distinct topological properties. Finally, we construct two-phase flow directed weighted networks from experimental signals and associate the dynamic behavior of gas-liquid two-phase flow with the topological statistics of the generated networks. The results suggest that the topological statistics of two-phase flow networks allow quantitative characterization of the dynamic flow behavior in the transitions among different gas-liquid flow patterns.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10972010,11028206,11371069,11372052,11402029,and 11472060)the Science and Technology Development Foundation of China Academy of Engineering Physics(CAEP),China(Grant No.2014B0201030)the Defense Industrial Technology Development Program of China(Grant No.B1520132012)
文摘Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method（GAHM） is presented for incompressible two-phase flow.A front tracking method（FTM） is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives（normals and curvatures） and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10872222 and 50921063)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20110191110037)the Fundamental Research Funds for the Central Universities,China(Grant Nos.CDJXS11240011 and CDJXS10241103)
文摘In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace＇s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.
基金supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) under the Ministry of Science,ICT&Future Planning(2018R1A2B6003623)。
文摘Two-phase flow in a horizontal pipe was investigated by using numerical and experimental visualization methods.A horizontal pipe was built for qualitative and quantitative flow visualization.The length of horizontal pipe flow system was 9.5 m and the inner diameter was 51 mm.High-speed video method was used for the qualitative visualization and PIV method was applied for the quantitative visua-lization.The same geometry model was used for the numerical study.Three flow regimes including stratified flow,elongated bubble and slug flow field were generated and visualized by using numerical and experimental methods.The results show that the numerical simulation results are qualitatively si-milar to that of the experimental results.In addition,more quantitative results can be analyzed by numerical method.Development and decay process of slug flow was investigated,showing that the decay of slug heavily depends on the magnitude of nose velocity and its lasting time.It can also be found that the liquid superficial velocity plays a significant role in affecting the slug frequency.When keeping the gas superficial velocity constant,the frequency will increase with the liquid superficial velocity.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.41174109 and 61104148)the National Science and Technology Major Project of China(Grant No.2011ZX05020-006)the Zhejiang Key Discipline of Instrument Science and Technology,China(Grant No.JL130106)
文摘We extend the complexity entropy causality plane（CECP） to propose a multi-scale complexity entropy causality plane（MS-CECP） and further use the proposed method to discriminate the deterministic characteristics of different oil-in-water flows. We first take several typical time series for example to investigate the characteristic of the MS-CECP and find that the MS-CECP not only describes the continuous loss of dynamical structure with the increase of scale, but also reflects the determinacy of the system. Then we calculate the MS-CECP for the conductance fluctuating signals measured from oil–water two-phase flow loop test facility. The results indicate that the MS-CECP could be an intrinsic measure for indicating oil-in-water two-phase flow structures.
基金supported by the National Natural Science Foundation of China(Grant No.11572062)the Fundamental Research Funds for the Central Universities,China(Grant No.CDJZR13248801)+2 种基金the Program for Changjiang Scholars and Innovative Research Team in University,China(Grant No.IRT13043)Key Laboratory of Functional Crystals and Laser Technology,TIPCChinese Academy of Sciences
文摘We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.
