The paper explores the gravity-driven flow of the thin film of a viscoelastic-fluid-based nanofluids(VFBN)along an inclined plane under non-isothermal conditions and subjected to convective cooling at the free-surface...The paper explores the gravity-driven flow of the thin film of a viscoelastic-fluid-based nanofluids(VFBN)along an inclined plane under non-isothermal conditions and subjected to convective cooling at the free-surface.The Newton’s law of cooling is used to model the convective heat-exchange with the ambient at the free-surface.The Giesekus viscoelastic constitutive model,with appropriate modifications to account for non-isothermal effects,is employed to describe the polymeric effects.The unsteady and coupled non-linear partial differential equations(PDEs)describing the model problem are obtained and solved via efficient semi-implicit numerical schemes based on finite difference methods(FDM)implemented in Matlab.The response of the VFBN velocity,temperature,thermal-conductivity and polymeric-stresses to variations in the volume-fraction of embedded nanoparticles is investigated.It is shown that these quantities all increase as the nanoparticle volume-fraction becomes higher.展开更多
The focus of the study is to examine thermal radiation and viscous dissipative heat transfers of magnetohydrodynamics (MHD) stagnation point flow past a permeable confined stretching cylinder with non-uniform heat sou...The focus of the study is to examine thermal radiation and viscous dissipative heat transfers of magnetohydrodynamics (MHD) stagnation point flow past a permeable confined stretching cylinder with non-uniform heat source or sink. The formulated equation governing the flow is non-dimensional. The dimensionless momentum and energy equation are solved using shooting technique coupled with fourth-order Runge-kutta integrated scheme which satisfied smoothness conditions at the edge of the boundary layer. The result for the velocity and temperature distributions are presented graphically and discussed to portray the effects of some important embodiment parameters on the flow. The Nusselt number and skin friction were obtained and compared with the previous scholars’ results in others to validate the present research work.展开更多
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.展开更多
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.展开更多
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.展开更多
Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat ...Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.展开更多
Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 t...Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km<sup>3</sup> to 0.08 km<sup>3</sup>. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.展开更多
BACKGROUND The recently introduced ultrasonic flow ratio(UFR),is a novel fast computational method to derive fractional flow reserve(FFR)from intravascular ultrasound(IVUS)images.In the present study,we evaluate the d...BACKGROUND The recently introduced ultrasonic flow ratio(UFR),is a novel fast computational method to derive fractional flow reserve(FFR)from intravascular ultrasound(IVUS)images.In the present study,we evaluate the diagnostic performance of UFR in patients with intermediate left main(LM)stenosis.METHODS This is a prospective,single center study enrolling consecutive patients with presence of intermediated LM lesions(diameter stenosis of 30%-80%by visual estimation)underwent IVUS and FFR measurement.An independent core laboratory assessed offline UFR and IVUS-derived minimal lumen area(MLA)in a blinded fashion.RESULTS Both UFR and FFR were successfully achieved in 41 LM patients(mean age,62.0±9.9 years,46.3%diabetes).An acceptable correlation between UFR and FFR was identified(r=0.688,P<0.0001),with an absolute numerical difference of 0.03(standard difference:0.01).The area under the curve(AUC)in diagnosis of physiologically significant coronary stenosis for UFR was 0.94(95%CI:0.87-1.01),which was significantly higher than angiographic identified stenosis>50%(AUC=0.66,P<0.001)and numerically higher than IVUS-derived MLA(AUC=0.82;P=0.09).Patient level diagnostic accuracy,sensitivity and specificity for UFR to identify FFR≤0.80 was 82.9%(95%CI:70.2-95.7),93.1%(95%CI:82.2-100.0),58.3%(95%CI:26.3-90.4),respectively.CONCLUSION In patients with intermediate LM diseases,UFR was proved to be associated with acceptable correlation and high accuracy with pressure wire-based FFR as standard reference.The present study supports the use of UFR for functional evaluation of intermediate LM stenosis.展开更多
文摘The paper explores the gravity-driven flow of the thin film of a viscoelastic-fluid-based nanofluids(VFBN)along an inclined plane under non-isothermal conditions and subjected to convective cooling at the free-surface.The Newton’s law of cooling is used to model the convective heat-exchange with the ambient at the free-surface.The Giesekus viscoelastic constitutive model,with appropriate modifications to account for non-isothermal effects,is employed to describe the polymeric effects.The unsteady and coupled non-linear partial differential equations(PDEs)describing the model problem are obtained and solved via efficient semi-implicit numerical schemes based on finite difference methods(FDM)implemented in Matlab.The response of the VFBN velocity,temperature,thermal-conductivity and polymeric-stresses to variations in the volume-fraction of embedded nanoparticles is investigated.It is shown that these quantities all increase as the nanoparticle volume-fraction becomes higher.
文摘The focus of the study is to examine thermal radiation and viscous dissipative heat transfers of magnetohydrodynamics (MHD) stagnation point flow past a permeable confined stretching cylinder with non-uniform heat source or sink. The formulated equation governing the flow is non-dimensional. The dimensionless momentum and energy equation are solved using shooting technique coupled with fourth-order Runge-kutta integrated scheme which satisfied smoothness conditions at the edge of the boundary layer. The result for the velocity and temperature distributions are presented graphically and discussed to portray the effects of some important embodiment parameters on the flow. The Nusselt number and skin friction were obtained and compared with the previous scholars’ results in others to validate the present research work.
基金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 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 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(52304067,62273213)the Natural Science Foundation of Shandong Province of China(ZR2021QE073)+1 种基金the Natural Science Foundation of Shandong Province for Innovation and Development Joint Funds(ZR2022LZH001)the China Postdoctoral Science Foundation(2023M732111)。
文摘Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.
文摘Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km<sup>3</sup> to 0.08 km<sup>3</sup>. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.
基金supported by CAMS Innovation Fund for Medical Sciences(CIFMS)(2022–12M-C&TB-043).
文摘BACKGROUND The recently introduced ultrasonic flow ratio(UFR),is a novel fast computational method to derive fractional flow reserve(FFR)from intravascular ultrasound(IVUS)images.In the present study,we evaluate the diagnostic performance of UFR in patients with intermediate left main(LM)stenosis.METHODS This is a prospective,single center study enrolling consecutive patients with presence of intermediated LM lesions(diameter stenosis of 30%-80%by visual estimation)underwent IVUS and FFR measurement.An independent core laboratory assessed offline UFR and IVUS-derived minimal lumen area(MLA)in a blinded fashion.RESULTS Both UFR and FFR were successfully achieved in 41 LM patients(mean age,62.0±9.9 years,46.3%diabetes).An acceptable correlation between UFR and FFR was identified(r=0.688,P<0.0001),with an absolute numerical difference of 0.03(standard difference:0.01).The area under the curve(AUC)in diagnosis of physiologically significant coronary stenosis for UFR was 0.94(95%CI:0.87-1.01),which was significantly higher than angiographic identified stenosis>50%(AUC=0.66,P<0.001)and numerically higher than IVUS-derived MLA(AUC=0.82;P=0.09).Patient level diagnostic accuracy,sensitivity and specificity for UFR to identify FFR≤0.80 was 82.9%(95%CI:70.2-95.7),93.1%(95%CI:82.2-100.0),58.3%(95%CI:26.3-90.4),respectively.CONCLUSION In patients with intermediate LM diseases,UFR was proved to be associated with acceptable correlation and high accuracy with pressure wire-based FFR as standard reference.The present study supports the use of UFR for functional evaluation of intermediate LM stenosis.