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 展开更多
Tapered ring with thin wall and three high ribs(TRTWTHR),showing complicated geometry(wall thickness is less than 4 mm and rib height exceeds 20 mm),is extensively utilized to fabricate the critical structural parts o...Tapered ring with thin wall and three high ribs(TRTWTHR),showing complicated geometry(wall thickness is less than 4 mm and rib height exceeds 20 mm),is extensively utilized to fabricate the critical structural parts of aerospace equipment such as spacecraft cabin,rocket body and fuel tank because of light weight and high carrying capacity.How to fabricate TRTWTHR with high performance is a critical problem that aerospace area needs to solve.In this work,constraining ring rolling(CRR)technique is first adopted to form TRTWTHR.However.unreasonable metal streamlines(UMS)and uncoordinated growth of three ribs easily occur in CRR of TRTWTHR,which makes the forming quality of TRTWTHR difficult to be controlled.Faced with this difficulty,an analytical model that can predict UMS and the height of three ribs in CRR of TRTWTHR is established so as to guide the process design of CRR.Subsequently,the reliability of the established analytical model and the feasibility of CRR of TRTWTHR are confirmed by FE simulation and experiment.Then,using the established analytical model,the window of UMS occurring relevant to the tapered angle of TRTWTHR and the location of the rib of middle end is developed.Finally,three uncoordinated growth modes among three ribs are found when the width of three ribs is identical and UMS do not occur,and the mechanisms of three uncoordinated growth modes are revealed.展开更多
The two-dimensional spreading under gravity of a thin fluid film with suction (fluid leak-off) or blowing (fluid injection) at the base is considered. The thin fluid film approximation is imposed. The height of the th...The two-dimensional spreading under gravity of a thin fluid film with suction (fluid leak-off) or blowing (fluid injection) at the base is considered. The thin fluid film approximation is imposed. The height of the thin film satisfies a nonlinear diffusion equation with a source/sink term. The Lie point symmetries of the nonlinear diffusion equation are derived and exist, which provided the fluid velocity at the base, <em>v<sub>n</sub></em> satisfies a first order linear partial differential equation. The general form has algebraic time dependence while a special case has exponential time dependence. The solution in which <em>v<sub>n</sub></em> is proportional to the height of the thin film is studied. The width of the base always increases with time even for suction while the height decreases with time for sufficiently weak blowing. The streamlines of the fluid flow inside the thin film are plotted by first solving a cubic equation. For sufficiently weak blowing there is a dividing streamline, emanating from the stagnation point on the centre line which separates the fluid flow into two regions, a lower region consisting of rising fluid and dominated by fluid injection at the base and an upper region consisting of descending fluid and dominated by spreading due to gravity. For sufficiently strong blowing the lower region expands to completely fill the whole thin film.展开更多
As a high specific speed pump, the contra-rotating axial flow pump distinguishes itself in a rear rotor rotating in the opposite direction of the front rotor, which remarkably contributes to the energy conversion, the...As a high specific speed pump, the contra-rotating axial flow pump distinguishes itself in a rear rotor rotating in the opposite direction of the front rotor, which remarkably contributes to the energy conversion, the reduction of the pump size, better hydraulic and cavitation performances. However, with two rotors rotating reversely, the significant interaction between blade rows was observed in our prototype contra-rotating rotors, which highly affected the pump performance compared with the conventional axial flow pumps. Consequently, a new type of rear rotor was designed by the rotational speed optimization methodology with some additional considerations, aiming at better cavitation performance, the reduction of blade rows interaction and the secondary flow suppression. The new rear rotor showed a satisfactory performance at the design flow rate but an unfavorable positive slope of the head - flow rate curve in the partial flow rate range less than 40% of the design flow rate, which should be avoided for the reliability of pump-pipe systems. In the present research, to understand the internal flow field of new rear rotor and its relation to the performances at the partial flow rates, the velocity distributions at the inlets and outlets of the rotors are firstly investigated. Then, the boundary layer flows on rotor surfaces, which clearly reflect the secondary flow inside the rotors, are analyzed through the limiting streamline observations using the multi-color oil-film method. Finally, the unsteady numerical simulations are carded out to understand the complicated internal flow structures in the rotors.展开更多
Based on the working of Lighthill and Hunt et al., in the present paper the author has established the topological rules adapting to analysing the skin-friction lines and the section streamlines in cascades. These rul...Based on the working of Lighthill and Hunt et al., in the present paper the author has established the topological rules adapting to analysing the skin-friction lines and the section streamlines in cascades. These rules are (1) for a rotor cascade without shroud band, the total number of nodal points equals that the saddle points on the skin-friction line vector fields in eachpitch range; (2) for an annular or straight cascade with no-clearances at blade ends, the total number of saddle points is two more than that of nodal points on the skin-friction line fields in a pitch; (3) the total number of saddles in the secondary flow fields on cross-sections in cascade is one less than that of nodes; (4) in the section streamline vector fields on a meridian surface penetrating a flow passage, and on leading and trailing edge sections, the total number of nodes is equal to that of saddles; (5) on the streamline vector fields of a blade-to-blade surface, the total number of nodes is one less than that of saddles.展开更多
Vortices that develop over intakes are a hazardous hydraulic phenomenon.In this study, a 3D model was developed to study the flow field in air-core vortices.This model is based on the spiral pattern of streamlines and...Vortices that develop over intakes are a hazardous hydraulic phenomenon.In this study, a 3D model was developed to study the flow field in air-core vortices.This model is based on the spiral pattern of streamlines and the analytical solution of the momentum and continuity equations for deriving the three components of velocity.The model provides equations for free surface profiles and 3D patterns of the streamlines.Moreover, a new relationship was suggested for calculating effective viscosity and its distribution across the vortex flow field.The performance of the proposed analytical model was compared with existing experimental data and the results of previous analytical models.The outcomes indicated that the proposed model could predict characteristics of the vortex flow with good accuracy.展开更多
The aim of this investigation is to determine the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous in...The aim of this investigation is to determine the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous incompressible fluid by studying the evolution of the streamlines in the thin film. It is assumed that the normal component of the fluid velocity at the base is proportional to the spatial gradient of the height of the film. Lie symmetry methods for partial differential equations are applied. The invariant solution for the surface profile is derived. It is found that the thin fluid film approximation is satisfied for weak to moderate leak-off and for the whole range of fluid injection. The streamlines are derived and plotted by solving a cubic equation numerically. For fluid injection, there is a dividing streamline originating at the stagnation point at the base which separates the flow into two regions, a lower region consisting mainly of rising fluid and an upper region consisting mainly of descending fluid. An approximate analytical solution for the dividing streamline is derived. It generates an approximate V-shaped surface along the length of the two-dimensional film with the vertex of each section the stagnation point. It is concluded that the fluid flow inside the thin film can be visualised by plotting the streamlines. Other models relating the fluid velocity at the base to the height of the thin film can be expected to contain a dividing streamline originating at a stagnation point and dividing the flow into a lower region of rising fluid and an upper region of descending fluid.展开更多
The mechanism of the fluid flow in low permeability reservoirs is different from that in middle-high permeability reservoirs because of the existence of the Threshold Pressure Gradient （TPG）. When the pressure gradi...The mechanism of the fluid flow in low permeability reservoirs is different from that in middle-high permeability reservoirs because of the existence of the Threshold Pressure Gradient （TPG）. When the pressure gradient at some location is greater than the TPG, the fluid in porous media begins to flow. By applying the mirror image method and the principle of potential superposition, the steady-state pressure distribution and the stream function for infinite five-spot well patterns can be obtained for a low permeability reservoir with the TPG effect. Based on the streamlines distribution, the flowing and stagnant zones in five-spot well patterns can be clearly seen. By the definition of the effective startup coefficient （SUC）, the ratio of the flowing and stagnant zones can be calculated accurately. It is shown that the SUC for five-spot well patterns is not constant, but decreases with the increase of the di- mensionless TPG. By increasing the effective permeability of the formation （such as by the acid treatment and the hydraulic fracture）, in increasing the injection-production differential pressure or shortening the well space （such as by infilling well）, the SUC can be improved. The results of the sensitivity analysis show that a better choice for the SUC enhancement is to shorten the well spacing for small permeability reservoirs and to increase the pressure difference for large permeability reservoirs. This streamline approach can be used to determine the distribution of remaining oil and provide guidance for infilling well.展开更多
基金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
基金the National Natural Science Foundation of China (No. U2037204)the 111 Project (No. B17034)+1 种基金Innovative Research Team Development Program of Ministry of Education of China (No. IRT17R83)the National Natural Science Foundation of China (No. 52005375)
文摘Tapered ring with thin wall and three high ribs(TRTWTHR),showing complicated geometry(wall thickness is less than 4 mm and rib height exceeds 20 mm),is extensively utilized to fabricate the critical structural parts of aerospace equipment such as spacecraft cabin,rocket body and fuel tank because of light weight and high carrying capacity.How to fabricate TRTWTHR with high performance is a critical problem that aerospace area needs to solve.In this work,constraining ring rolling(CRR)technique is first adopted to form TRTWTHR.However.unreasonable metal streamlines(UMS)and uncoordinated growth of three ribs easily occur in CRR of TRTWTHR,which makes the forming quality of TRTWTHR difficult to be controlled.Faced with this difficulty,an analytical model that can predict UMS and the height of three ribs in CRR of TRTWTHR is established so as to guide the process design of CRR.Subsequently,the reliability of the established analytical model and the feasibility of CRR of TRTWTHR are confirmed by FE simulation and experiment.Then,using the established analytical model,the window of UMS occurring relevant to the tapered angle of TRTWTHR and the location of the rib of middle end is developed.Finally,three uncoordinated growth modes among three ribs are found when the width of three ribs is identical and UMS do not occur,and the mechanisms of three uncoordinated growth modes are revealed.
文摘The two-dimensional spreading under gravity of a thin fluid film with suction (fluid leak-off) or blowing (fluid injection) at the base is considered. The thin fluid film approximation is imposed. The height of the thin film satisfies a nonlinear diffusion equation with a source/sink term. The Lie point symmetries of the nonlinear diffusion equation are derived and exist, which provided the fluid velocity at the base, <em>v<sub>n</sub></em> satisfies a first order linear partial differential equation. The general form has algebraic time dependence while a special case has exponential time dependence. The solution in which <em>v<sub>n</sub></em> is proportional to the height of the thin film is studied. The width of the base always increases with time even for suction while the height decreases with time for sufficiently weak blowing. The streamlines of the fluid flow inside the thin film are plotted by first solving a cubic equation. For sufficiently weak blowing there is a dividing streamline, emanating from the stagnation point on the centre line which separates the fluid flow into two regions, a lower region consisting of rising fluid and dominated by fluid injection at the base and an upper region consisting of descending fluid and dominated by spreading due to gravity. For sufficiently strong blowing the lower region expands to completely fill the whole thin film.
文摘As a high specific speed pump, the contra-rotating axial flow pump distinguishes itself in a rear rotor rotating in the opposite direction of the front rotor, which remarkably contributes to the energy conversion, the reduction of the pump size, better hydraulic and cavitation performances. However, with two rotors rotating reversely, the significant interaction between blade rows was observed in our prototype contra-rotating rotors, which highly affected the pump performance compared with the conventional axial flow pumps. Consequently, a new type of rear rotor was designed by the rotational speed optimization methodology with some additional considerations, aiming at better cavitation performance, the reduction of blade rows interaction and the secondary flow suppression. The new rear rotor showed a satisfactory performance at the design flow rate but an unfavorable positive slope of the head - flow rate curve in the partial flow rate range less than 40% of the design flow rate, which should be avoided for the reliability of pump-pipe systems. In the present research, to understand the internal flow field of new rear rotor and its relation to the performances at the partial flow rates, the velocity distributions at the inlets and outlets of the rotors are firstly investigated. Then, the boundary layer flows on rotor surfaces, which clearly reflect the secondary flow inside the rotors, are analyzed through the limiting streamline observations using the multi-color oil-film method. Finally, the unsteady numerical simulations are carded out to understand the complicated internal flow structures in the rotors.
文摘Based on the working of Lighthill and Hunt et al., in the present paper the author has established the topological rules adapting to analysing the skin-friction lines and the section streamlines in cascades. These rules are (1) for a rotor cascade without shroud band, the total number of nodal points equals that the saddle points on the skin-friction line vector fields in eachpitch range; (2) for an annular or straight cascade with no-clearances at blade ends, the total number of saddle points is two more than that of nodal points on the skin-friction line fields in a pitch; (3) the total number of saddles in the secondary flow fields on cross-sections in cascade is one less than that of nodes; (4) in the section streamline vector fields on a meridian surface penetrating a flow passage, and on leading and trailing edge sections, the total number of nodes is equal to that of saddles; (5) on the streamline vector fields of a blade-to-blade surface, the total number of nodes is one less than that of saddles.
基金supported by the Iran National Science Foundation(INSF,Grant No.97008045)
文摘Vortices that develop over intakes are a hazardous hydraulic phenomenon.In this study, a 3D model was developed to study the flow field in air-core vortices.This model is based on the spiral pattern of streamlines and the analytical solution of the momentum and continuity equations for deriving the three components of velocity.The model provides equations for free surface profiles and 3D patterns of the streamlines.Moreover, a new relationship was suggested for calculating effective viscosity and its distribution across the vortex flow field.The performance of the proposed analytical model was compared with existing experimental data and the results of previous analytical models.The outcomes indicated that the proposed model could predict characteristics of the vortex flow with good accuracy.
文摘The aim of this investigation is to determine the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous incompressible fluid by studying the evolution of the streamlines in the thin film. It is assumed that the normal component of the fluid velocity at the base is proportional to the spatial gradient of the height of the film. Lie symmetry methods for partial differential equations are applied. The invariant solution for the surface profile is derived. It is found that the thin fluid film approximation is satisfied for weak to moderate leak-off and for the whole range of fluid injection. The streamlines are derived and plotted by solving a cubic equation numerically. For fluid injection, there is a dividing streamline originating at the stagnation point at the base which separates the flow into two regions, a lower region consisting mainly of rising fluid and an upper region consisting mainly of descending fluid. An approximate analytical solution for the dividing streamline is derived. It generates an approximate V-shaped surface along the length of the two-dimensional film with the vertex of each section the stagnation point. It is concluded that the fluid flow inside the thin film can be visualised by plotting the streamlines. Other models relating the fluid velocity at the base to the height of the thin film can be expected to contain a dividing streamline originating at a stagnation point and dividing the flow into a lower region of rising fluid and an upper region of descending fluid.
基金Project Supported by the National Natural Science Foundation of China(Grant No.51204148)
文摘The mechanism of the fluid flow in low permeability reservoirs is different from that in middle-high permeability reservoirs because of the existence of the Threshold Pressure Gradient （TPG）. When the pressure gradient at some location is greater than the TPG, the fluid in porous media begins to flow. By applying the mirror image method and the principle of potential superposition, the steady-state pressure distribution and the stream function for infinite five-spot well patterns can be obtained for a low permeability reservoir with the TPG effect. Based on the streamlines distribution, the flowing and stagnant zones in five-spot well patterns can be clearly seen. By the definition of the effective startup coefficient （SUC）, the ratio of the flowing and stagnant zones can be calculated accurately. It is shown that the SUC for five-spot well patterns is not constant, but decreases with the increase of the di- mensionless TPG. By increasing the effective permeability of the formation （such as by the acid treatment and the hydraulic fracture）, in increasing the injection-production differential pressure or shortening the well space （such as by infilling well）, the SUC can be improved. The results of the sensitivity analysis show that a better choice for the SUC enhancement is to shorten the well spacing for small permeability reservoirs and to increase the pressure difference for large permeability reservoirs. This streamline approach can be used to determine the distribution of remaining oil and provide guidance for infilling well.
