Darcy’s law is widely used to describe the flow in porous media in which there is a linear relationship between fluid velocity and pressure gradient. However, it has been found that for high numbers of Reynolds this ...Darcy’s law is widely used to describe the flow in porous media in which there is a linear relationship between fluid velocity and pressure gradient. However, it has been found that for high numbers of Reynolds this law ceases to be valid. In this work, the Ergun equation is employed to consider the non-linearity of air velocity with the pressure gradient in casting sands. The contribution of non-linearity to the total flow in terms of a variable defined as a non-Darcy flow fraction is numerically quantified. In addition, the influence of the shape factor of the sand grains on the non-linear flow fraction is analyzed. It is found that for values of the Reynolds number less or equal than 1, the contribution of non-linearity for spherical particles is around 1.15%.展开更多
Water seepage in soil is a fundamental problem involving various scientific and engineering fields.According to the literature,low-velocity water seepage in low-permeability porous media,such as clay,does not follow D...Water seepage in soil is a fundamental problem involving various scientific and engineering fields.According to the literature,low-velocity water seepage in low-permeability porous media,such as clay,does not follow Darcy's law,also known as pre-Darcy flow.The formation of immovable water due to water adsorption on the pore wall is believed to be responsible for the formation of pre-Darcy flow.However,this view lacks direct solid evidence.To investigate the pre-Darcy water flow in clay,head permeability experiments are conducted on six clay samples with different densities.The results indicate that water seepage in clay at low hydraulic gradients does not follow Darcy's law.A clear nonlinear relationship between flow velocity and hydraulic gradient is observed.Water flow in clay can be divided into the pre-Darcy flow and Darcy flow regions by the critical hydraulic gradient,which is 10-12 for the Albic soil with dry density between 1.3 g/cm^(3)and 1.8 g/cm^(3).According to the disjoining pressure theory,immovable water due to water adsorption on the pore wall is the primary reason for water flow deviating from Darcy's law in clay.The results indicate that the percentage of movable water ranges from 39.7%to 59.3%for the six samples at a hydraulic gradient of 1.As the hydraulic gradient increases,the percentage of moveable water also increases.Additionally,there is a strong correlation between the percentage of movable water and the variation in hydraulic conductivity with the hydraulic gradient.Furthermore,a quantitative relationship between the percentage of movable water and the hydraulic conductivity has been established.The results of this study suggest that water adsorption on the pore wall not only affects the water movability,but is also closely related to the pre-Darcy flow phenomenon in clay.展开更多
In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)...In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)model that can simultaneously consider the fracture evolution and non-Darcy flow dynamics in studying the thermo-hydro-mechanical(THM)coupling processes for heat extraction in geothermal reservoir.We further employed the model on the Habanero enhanced geothermal systems(EGS)project located in Australia.First,our findings illustrate a clear spatial-temporal variation in the thermal stress and pressure perturbations,as well as uneven spatial distribution of shear failure in 3D fracture networks.Activated shear failure is mainly concentrated in the first fracture cluster.Secondly,channeling flow have also been observed in DFNs during heat extraction and are further intensified by the expansion of fractures driven by thermal stresses.Moreover,the combined effect of non-Darcy flow and fracture evolution triggers a rapid decline in the resulting heat rate and temperature.The NR-DFN model framework and the Habanero EGS's results illustrate the importance of both fracture evolution and non-Darcy flow on the efficiency of EGS production and have the potential to promote the development of more sustainable and efficient EGS operations for stakeholders.展开更多
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.展开更多
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.展开更多
A three-dimensional Darcy Forchheimer mixed convective flow of a couple stress hybrid nanofluid flow through a vertical plate by means of the double diffusion Cattaneo-Christov model is presented in this study.The inf...A three-dimensional Darcy Forchheimer mixed convective flow of a couple stress hybrid nanofluid flow through a vertical plate by means of the double diffusion Cattaneo-Christov model is presented in this study.The influence of highorder velocity slip flow,as well as a passive and active control,is also considered.The motive of the research is to develop a computational model,using cobalt ferrite(Co Fe_(2)O_(4))and copper(Cu)nanoparticles(NPs)in the carrier fluid water,to magnify the energy and mass communication rate and boost the efficiency and performance of thermal energy conduction for a variety of commercial and biological purposes.The proposed model becomes more significant,with an additional effect of non-Fick's mass flux and Fourier's heat model to report the energy and mass passage rate.The results are obtained through the computational strategy parametric continuation method.The figures are plotted to reveal the physical sketch of the obtained solution,while the statistical assessment has been evaluated through tables.It has been observed that the dispersion of Cu and Co Fe_(2)O_(4)NPs to the base fluid significantly enhances the velocity and thermal conductivity of water,which is the most remarkable property of these NPs from the industrial point of view.展开更多
Measurement of bloodflow velocity is key to understanding physiology and pathology in vivo.While most measurements are performed at the middle of the blood vessel,little research has been done on characterizing the in...Measurement of bloodflow velocity is key to understanding physiology and pathology in vivo.While most measurements are performed at the middle of the blood vessel,little research has been done on characterizing the instantaneous bloodflow velocity distribution.This is mainly due to the lack of measurement technology with high spatial and temporal resolution.Here,we tackle this problem with our recently developed dual-wavelength line-scan third-harmonic generation(THG)imaging technology.Simultaneous acquisition of dual-wavelength THG line-scanning signals enables measurement of bloodflow velocities at two radially symmetric positions in both venules and arterioles in mouse brain in vivo.Our results clearly show that the instantaneous bloodflow velocity is not symmetric under general conditions.展开更多
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.展开更多
Brazilian pre-salt reservoirs are renowned for their intricate pore networks and vuggy nature,posing significant challenges in modeling and simulating fluid flow within these carbonate reservoirs.Despite possessing ex...Brazilian pre-salt reservoirs are renowned for their intricate pore networks and vuggy nature,posing significant challenges in modeling and simulating fluid flow within these carbonate reservoirs.Despite possessing excellent petrophysical properties,such as high porosity and permeability,these reservoirs typically exhibit a notably low recovery factor,sometimes falling below 10%.Previous research has indicated that various enhanced oil recovery(EOR)methods,such as water alternating gas(WAG),can substantially augment the recovery factor in pre-salt reservoirs,resulting in improvements of up to 20%.Nevertheless,the fluid flow mechanism within Brazilian carbonate reservoirs,characterized by complex pore geometry,remains unclear.Our study examines the behavior of fluid flow in a similar heterogeneous porous material,utilizing a plug sample obtained from a vugular segment of a Brazilian stromatolite outcrop,known to share analogies with certain pre-salt reservoirs.We conducted single-phase and multi-phase core flooding experiments,complemented by medical-CT scanning,to generate flow streamlines and evaluate the efficiency of water flooding.Subsequently,micro-CT scanning of the core sample was performed,and two cross-sections from horizontal and vertical plates were constructed.These cross-sections were then employed as geometries in a numerical simulator,enabling us to investigate the impact of pore geometry on fluid flow.Analysis of the pore-scale modeling and experimental data unveiled that the presence of dead-end pores and vugs results in a significant portion of the fluid remaining stagnant within these regions.Consequently,the injected fluid exhibits channeling-like behavior,leading to rapid breakthrough and low areal swept efficiency.Additionally,the numerical simulation results demonstrated that,irrespective of the size of the dead-end regions,the pressure variation within the dead-end vugs and pores is negligible.Despite the stromatolite's favorable petrophysical properties,including relatively high porosity and permeability,as well as the presence of interconnected large vugs,the recovery factor during water flooding remained low due to early breakthrough.These findings align with field data obtained from pre-salt reservoirs,providing an explanation for the observed low recovery factor during water flooding in such reservoirs.展开更多
文摘Darcy’s law is widely used to describe the flow in porous media in which there is a linear relationship between fluid velocity and pressure gradient. However, it has been found that for high numbers of Reynolds this law ceases to be valid. In this work, the Ergun equation is employed to consider the non-linearity of air velocity with the pressure gradient in casting sands. The contribution of non-linearity to the total flow in terms of a variable defined as a non-Darcy flow fraction is numerically quantified. In addition, the influence of the shape factor of the sand grains on the non-linear flow fraction is analyzed. It is found that for values of the Reynolds number less or equal than 1, the contribution of non-linearity for spherical particles is around 1.15%.
基金funding support from the China National Science Foundation(Grant Nos.42072280,U2244215,41172205).
文摘Water seepage in soil is a fundamental problem involving various scientific and engineering fields.According to the literature,low-velocity water seepage in low-permeability porous media,such as clay,does not follow Darcy's law,also known as pre-Darcy flow.The formation of immovable water due to water adsorption on the pore wall is believed to be responsible for the formation of pre-Darcy flow.However,this view lacks direct solid evidence.To investigate the pre-Darcy water flow in clay,head permeability experiments are conducted on six clay samples with different densities.The results indicate that water seepage in clay at low hydraulic gradients does not follow Darcy's law.A clear nonlinear relationship between flow velocity and hydraulic gradient is observed.Water flow in clay can be divided into the pre-Darcy flow and Darcy flow regions by the critical hydraulic gradient,which is 10-12 for the Albic soil with dry density between 1.3 g/cm^(3)and 1.8 g/cm^(3).According to the disjoining pressure theory,immovable water due to water adsorption on the pore wall is the primary reason for water flow deviating from Darcy's law in clay.The results indicate that the percentage of movable water ranges from 39.7%to 59.3%for the six samples at a hydraulic gradient of 1.As the hydraulic gradient increases,the percentage of moveable water also increases.Additionally,there is a strong correlation between the percentage of movable water and the variation in hydraulic conductivity with the hydraulic gradient.Furthermore,a quantitative relationship between the percentage of movable water and the hydraulic conductivity has been established.The results of this study suggest that water adsorption on the pore wall not only affects the water movability,but is also closely related to the pre-Darcy flow phenomenon in clay.
基金funded by the National Natural Science Foundation of China (No.U22A20166)Science and Technology Foundation of Guizhou Province (No.QKHJC-ZK[2023]YB074)+2 种基金Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical EngineeringInstitute of Rock and Soil MechanicsChinese Academy of Sciences (No.SKLGME022009)。
文摘In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)model that can simultaneously consider the fracture evolution and non-Darcy flow dynamics in studying the thermo-hydro-mechanical(THM)coupling processes for heat extraction in geothermal reservoir.We further employed the model on the Habanero enhanced geothermal systems(EGS)project located in Australia.First,our findings illustrate a clear spatial-temporal variation in the thermal stress and pressure perturbations,as well as uneven spatial distribution of shear failure in 3D fracture networks.Activated shear failure is mainly concentrated in the first fracture cluster.Secondly,channeling flow have also been observed in DFNs during heat extraction and are further intensified by the expansion of fractures driven by thermal stresses.Moreover,the combined effect of non-Darcy flow and fracture evolution triggers a rapid decline in the resulting heat rate and temperature.The NR-DFN model framework and the Habanero EGS's results illustrate the importance of both fracture evolution and non-Darcy flow on the efficiency of EGS production and have the potential to promote the development of more sustainable and efficient EGS operations for stakeholders.
基金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.
基金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.
基金Deanship of Scientific Research at King Khalid University for funding this work through Large Groups Project under grant number(RGP.2/155/43)。
文摘A three-dimensional Darcy Forchheimer mixed convective flow of a couple stress hybrid nanofluid flow through a vertical plate by means of the double diffusion Cattaneo-Christov model is presented in this study.The influence of highorder velocity slip flow,as well as a passive and active control,is also considered.The motive of the research is to develop a computational model,using cobalt ferrite(Co Fe_(2)O_(4))and copper(Cu)nanoparticles(NPs)in the carrier fluid water,to magnify the energy and mass communication rate and boost the efficiency and performance of thermal energy conduction for a variety of commercial and biological purposes.The proposed model becomes more significant,with an additional effect of non-Fick's mass flux and Fourier's heat model to report the energy and mass passage rate.The results are obtained through the computational strategy parametric continuation method.The figures are plotted to reveal the physical sketch of the obtained solution,while the statistical assessment has been evaluated through tables.It has been observed that the dispersion of Cu and Co Fe_(2)O_(4)NPs to the base fluid significantly enhances the velocity and thermal conductivity of water,which is the most remarkable property of these NPs from the industrial point of view.
基金funded by the National Natural Science Foundation of China(Grant/Award Numbers 62075135 and 61975126)the Science and Technology Innovation Commission of Shenzhen(Grant/Award Numbers JCYJ20190808174819083 and JCYJ20190808175201640)Shenzhen Science and Technology Planning Project(ZDSYS 20210623092006020).
文摘Measurement of bloodflow velocity is key to understanding physiology and pathology in vivo.While most measurements are performed at the middle of the blood vessel,little research has been done on characterizing the instantaneous bloodflow velocity distribution.This is mainly due to the lack of measurement technology with high spatial and temporal resolution.Here,we tackle this problem with our recently developed dual-wavelength line-scan third-harmonic generation(THG)imaging technology.Simultaneous acquisition of dual-wavelength THG line-scanning signals enables measurement of bloodflow velocities at two radially symmetric positions in both venules and arterioles in mouse brain in vivo.Our results clearly show that the instantaneous bloodflow velocity is not symmetric under general conditions.
基金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.
基金the support of EPIC-Energy Production Innovation Center,hosted by the University of Campinas(UNICAMP)sponsored by FAPESP-Sao Paulo Research Foundation(2017/15736e3 process).
文摘Brazilian pre-salt reservoirs are renowned for their intricate pore networks and vuggy nature,posing significant challenges in modeling and simulating fluid flow within these carbonate reservoirs.Despite possessing excellent petrophysical properties,such as high porosity and permeability,these reservoirs typically exhibit a notably low recovery factor,sometimes falling below 10%.Previous research has indicated that various enhanced oil recovery(EOR)methods,such as water alternating gas(WAG),can substantially augment the recovery factor in pre-salt reservoirs,resulting in improvements of up to 20%.Nevertheless,the fluid flow mechanism within Brazilian carbonate reservoirs,characterized by complex pore geometry,remains unclear.Our study examines the behavior of fluid flow in a similar heterogeneous porous material,utilizing a plug sample obtained from a vugular segment of a Brazilian stromatolite outcrop,known to share analogies with certain pre-salt reservoirs.We conducted single-phase and multi-phase core flooding experiments,complemented by medical-CT scanning,to generate flow streamlines and evaluate the efficiency of water flooding.Subsequently,micro-CT scanning of the core sample was performed,and two cross-sections from horizontal and vertical plates were constructed.These cross-sections were then employed as geometries in a numerical simulator,enabling us to investigate the impact of pore geometry on fluid flow.Analysis of the pore-scale modeling and experimental data unveiled that the presence of dead-end pores and vugs results in a significant portion of the fluid remaining stagnant within these regions.Consequently,the injected fluid exhibits channeling-like behavior,leading to rapid breakthrough and low areal swept efficiency.Additionally,the numerical simulation results demonstrated that,irrespective of the size of the dead-end regions,the pressure variation within the dead-end vugs and pores is negligible.Despite the stromatolite's favorable petrophysical properties,including relatively high porosity and permeability,as well as the presence of interconnected large vugs,the recovery factor during water flooding remained low due to early breakthrough.These findings align with field data obtained from pre-salt reservoirs,providing an explanation for the observed low recovery factor during water flooding in such reservoirs.
