It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius fo...It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius for reducing the risk of sand accumulation.In this study,three-dimensional models of desert highway embankments with different vertical curve radii were constructed,and Fluent software was used to simulate the wind-sand flow field and sand accumulation distribution of vertical curve embankments.The results show that:(1)Along the direction of the road,the concave and the convex vertical curve embankments have the effect of collecting and diverging the wind-sand flow,respectively.When the radius of the concave vertical curve is 3000 m,5000 m,8000 m,10000 m and 20000 m,the wind velocity in the middle of the vertical curve is 31.76%,22.58%,10.78%,10.53%and 10.44%,higher than that at both ends.When the radius of the convex vertical curve is 6500 m,8000 m,10000 m,20000 m and 30000 m,the wind velocity at both ends of the vertical curve is 14.06%,9.99%,6.14%,3.22%and 2.41%,higher than that in the middle.The diversion effect also decreases with the increase of the radius.(2)The conductivity of the concave and convex vertical curve embankments with different radii is greater than 1,which is the sediment transport roadbed.The conductivity increases with the increase of radius and gradually tends to be stable.When the radius of the concave and convex vertical curves reaches 8000 m and 20000 m respectively,the phenomenon of sand accumulation is no longer serious.Under the same radius condition,the concave vertical curve embankment is more prone to sand accumulation than the convex one.(3)Considering the strength of the collection and diversion of the vertical curve embankment with different radii,and the sand accumulation of the vertical curve embankment in the desert section of Wuma Expressway,the radius of the concave vertical curve is not less than 8000 m,and the radius of the convex vertical curve is not less than 20000 m,which can effectively reduce the sand accumulation of the vertical curve embankment.In the desert highway area,the research results of this paper can provide reference for the design of vertical curve to ensure the safe operation of desert highway.展开更多
The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Prov...The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.展开更多
The distributed permutation flow shop scheduling problem(DPFSP)has received increasing attention in recent years.The iterated greedy algorithm(IGA)serves as a powerful optimizer for addressing such a problem because o...The distributed permutation flow shop scheduling problem(DPFSP)has received increasing attention in recent years.The iterated greedy algorithm(IGA)serves as a powerful optimizer for addressing such a problem because of its straightforward,single-solution evolution framework.However,a potential draw-back of IGA is the lack of utilization of historical information,which could lead to an imbalance between exploration and exploitation,especially in large-scale DPFSPs.As a consequence,this paper develops an IGA with memory and learning mechanisms(MLIGA)to efficiently solve the DPFSP targeted at the mini-malmakespan.InMLIGA,we incorporate a memory mechanism to make a more informed selection of the initial solution at each stage of the search,by extending,reconstructing,and reinforcing the information from previous solutions.In addition,we design a twolayer cooperative reinforcement learning approach to intelligently determine the key parameters of IGA and the operations of the memory mechanism.Meanwhile,to ensure that the experience generated by each perturbation operator is fully learned and to reduce the prior parameters of MLIGA,a probability curve-based acceptance criterion is proposed by combining a cube root function with custom rules.At last,a discrete adaptive learning rate is employed to enhance the stability of the memory and learningmechanisms.Complete ablation experiments are utilized to verify the effectiveness of the memory mechanism,and the results show that this mechanism is capable of improving the performance of IGA to a large extent.Furthermore,through comparative experiments involving MLIGA and five state-of-the-art algorithms on 720 benchmarks,we have discovered that MLI-GA demonstrates significant potential for solving large-scale DPFSPs.This indicates that MLIGA is well-suited for real-world distributed flow shop scheduling.展开更多
The multi-terminal direct current(DC)grid has extinctive superiorities over the traditional alternating current system in integrating large-scale renewable energy.Both the DC circuit breaker(DCCB)and the current flow ...The multi-terminal direct current(DC)grid has extinctive superiorities over the traditional alternating current system in integrating large-scale renewable energy.Both the DC circuit breaker(DCCB)and the current flow controller(CFC)are demanded to ensure the multiterminal DC grid to operates reliably and flexibly.However,since the CFC and the DCCB are all based on fully controlled semiconductor switches(e.g.,insulated gate bipolar transistor,integrated gate commutated thyristor,etc.),their separation configuration in the multiterminal DC grid will lead to unaffordable implementation costs and conduction power losses.To solve these problems,integrated equipment with both current flow control and fault isolation abilities is proposed,which shares the expensive and duplicated components of CFCs and DCCBs among adjacent lines.In addition,the complicated coordination control of CFCs and DCCBs can be avoided by adopting the integrated equipment in themultiterminal DC grid.In order to examine the current flow control and fault isolation abilities of the integrated equipment,the simulation model of a specific meshed four-terminal DC grid is constructed in the PSCAD/EMTDC software.Finally,the comparison between the integrated equipment and the separate solution is presented a specific result or conclusion needs to be added to the abstract.展开更多
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.展开更多
The flow field near a spur dike such as down flow and horseshoe vortex system(HVS)are susceptible to the topographic changes in the local scouring process,resulting in variation of the sediment transport with time.In ...The flow field near a spur dike such as down flow and horseshoe vortex system(HVS)are susceptible to the topographic changes in the local scouring process,resulting in variation of the sediment transport with time.In this study,large eddy simulations with fixed-bed at different scouring stages were conducted to investigate the changes in flow field.The results imply that the bed deformation leads to an increase in flow rate per unit area,which represent the capability of sediment transportation by water,in the scour hole.Moreover,the intensity of turbulent kinetic energy and bimodal motion near the sand bed induced by the HVS were also varied.However,the peak moments between the two sediment transport mechanisms were different.Hence,understanding the complex feedback mechanism between topography and flow field is essential for the local scour problem.展开更多
The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study are...The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study area which is extending along Karakorum Highway(KKH) from Besham to Chilas. Intense seismicity, deep gorges, steep terrain and extreme climatic events trigger multiple mountain hazards along the KKH, among which debris flow is recognized as the most destructive geohazard. This study aims to prepare a field-based debris flow inventory map at a regional scale along a 200 km stretch from Besham to Chilas. A total of 117 debris flows were identified in the field, and subsequently, a point-based debris-flow inventory and catchment delineation were performed through Arc GIS analysis. Regional scale debris flow susceptibility and propagation maps were prepared using Weighted Overlay Method(WOM) and Flow-R technique sequentially. Predisposing factors include slope, slope aspect, elevation, Topographic Roughness Index(TRI), Topographic Wetness Index(TWI), stream buffer, distance to faults, lithology rainfall, curvature, and collapsed material layer. The dataset was randomly divided into training data(75%) and validation data(25%). Results were validated through the Receiver Operator Characteristics(ROC) curve. Results show that Area Under the Curve(AUC) using WOM model is 79.2%. Flow-R propagation of debris flow shows that the 13.15%, 22.94%, and 63.91% areas are very high, high, and low susceptible to debris flow respectively. The propagation predicated by Flow-R validates the naturally occurring debris flow propagation as observed in the field surveys. The output of this research will provide valuable input to the decision makers for the site selection, designing of the prevention system, and for the protection of current infrastructure.展开更多
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.展开更多
基金The research described in this paper was financially supported by Youth Science Foundation Project’Research on Failure Mechanism and Evaluation Method of Sand Control Measures for Railway Machinery in Sandy Area’(12302511)Ningxia Transportation Department Science and Technology Project(20200173)Central guide local science and technology development funds(22ZY1QA005)。
文摘It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius for reducing the risk of sand accumulation.In this study,three-dimensional models of desert highway embankments with different vertical curve radii were constructed,and Fluent software was used to simulate the wind-sand flow field and sand accumulation distribution of vertical curve embankments.The results show that:(1)Along the direction of the road,the concave and the convex vertical curve embankments have the effect of collecting and diverging the wind-sand flow,respectively.When the radius of the concave vertical curve is 3000 m,5000 m,8000 m,10000 m and 20000 m,the wind velocity in the middle of the vertical curve is 31.76%,22.58%,10.78%,10.53%and 10.44%,higher than that at both ends.When the radius of the convex vertical curve is 6500 m,8000 m,10000 m,20000 m and 30000 m,the wind velocity at both ends of the vertical curve is 14.06%,9.99%,6.14%,3.22%and 2.41%,higher than that in the middle.The diversion effect also decreases with the increase of the radius.(2)The conductivity of the concave and convex vertical curve embankments with different radii is greater than 1,which is the sediment transport roadbed.The conductivity increases with the increase of radius and gradually tends to be stable.When the radius of the concave and convex vertical curves reaches 8000 m and 20000 m respectively,the phenomenon of sand accumulation is no longer serious.Under the same radius condition,the concave vertical curve embankment is more prone to sand accumulation than the convex one.(3)Considering the strength of the collection and diversion of the vertical curve embankment with different radii,and the sand accumulation of the vertical curve embankment in the desert section of Wuma Expressway,the radius of the concave vertical curve is not less than 8000 m,and the radius of the convex vertical curve is not less than 20000 m,which can effectively reduce the sand accumulation of the vertical curve embankment.In the desert highway area,the research results of this paper can provide reference for the design of vertical curve to ensure the safe operation of desert highway.
基金Under the auspices of National Natural Science Foundation of China (No 40571019)
文摘The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.
基金supported in part by the National Key Research and Development Program of China under Grant No.2021YFF0901300in part by the National Natural Science Foundation of China under Grant Nos.62173076 and 72271048.
文摘The distributed permutation flow shop scheduling problem(DPFSP)has received increasing attention in recent years.The iterated greedy algorithm(IGA)serves as a powerful optimizer for addressing such a problem because of its straightforward,single-solution evolution framework.However,a potential draw-back of IGA is the lack of utilization of historical information,which could lead to an imbalance between exploration and exploitation,especially in large-scale DPFSPs.As a consequence,this paper develops an IGA with memory and learning mechanisms(MLIGA)to efficiently solve the DPFSP targeted at the mini-malmakespan.InMLIGA,we incorporate a memory mechanism to make a more informed selection of the initial solution at each stage of the search,by extending,reconstructing,and reinforcing the information from previous solutions.In addition,we design a twolayer cooperative reinforcement learning approach to intelligently determine the key parameters of IGA and the operations of the memory mechanism.Meanwhile,to ensure that the experience generated by each perturbation operator is fully learned and to reduce the prior parameters of MLIGA,a probability curve-based acceptance criterion is proposed by combining a cube root function with custom rules.At last,a discrete adaptive learning rate is employed to enhance the stability of the memory and learningmechanisms.Complete ablation experiments are utilized to verify the effectiveness of the memory mechanism,and the results show that this mechanism is capable of improving the performance of IGA to a large extent.Furthermore,through comparative experiments involving MLIGA and five state-of-the-art algorithms on 720 benchmarks,we have discovered that MLI-GA demonstrates significant potential for solving large-scale DPFSPs.This indicates that MLIGA is well-suited for real-world distributed flow shop scheduling.
基金supported in part by Natural Science Foundation of Jiangsu Province under Grant BK20230255Natural Science Foundation of Shandong Province under Grant ZR2023QE281.
文摘The multi-terminal direct current(DC)grid has extinctive superiorities over the traditional alternating current system in integrating large-scale renewable energy.Both the DC circuit breaker(DCCB)and the current flow controller(CFC)are demanded to ensure the multiterminal DC grid to operates reliably and flexibly.However,since the CFC and the DCCB are all based on fully controlled semiconductor switches(e.g.,insulated gate bipolar transistor,integrated gate commutated thyristor,etc.),their separation configuration in the multiterminal DC grid will lead to unaffordable implementation costs and conduction power losses.To solve these problems,integrated equipment with both current flow control and fault isolation abilities is proposed,which shares the expensive and duplicated components of CFCs and DCCBs among adjacent lines.In addition,the complicated coordination control of CFCs and DCCBs can be avoided by adopting the integrated equipment in themultiterminal DC grid.In order to examine the current flow control and fault isolation abilities of the integrated equipment,the simulation model of a specific meshed four-terminal DC grid is constructed in the PSCAD/EMTDC software.Finally,the comparison between the integrated equipment and the separate solution is presented a specific result or conclusion needs to be added to the abstract.
基金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.
基金supported by Shenzhen Science and Technology Program(Grant No.JCYJ20220818102012024)Hong Kong Research Grants Council(Grant Nos.T21–602/16-R and RGC R5037–18)。
文摘The flow field near a spur dike such as down flow and horseshoe vortex system(HVS)are susceptible to the topographic changes in the local scouring process,resulting in variation of the sediment transport with time.In this study,large eddy simulations with fixed-bed at different scouring stages were conducted to investigate the changes in flow field.The results imply that the bed deformation leads to an increase in flow rate per unit area,which represent the capability of sediment transportation by water,in the scour hole.Moreover,the intensity of turbulent kinetic energy and bimodal motion near the sand bed induced by the HVS were also varied.However,the peak moments between the two sediment transport mechanisms were different.Hence,understanding the complex feedback mechanism between topography and flow field is essential for the local scour problem.
基金financially supported by the Higher Education Commission of Pakistan (HEC) grant under National Research Program for Universities (NRPU) with No: (20-14681/NRPU/R&D/HEC/20212021)。
文摘The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study area which is extending along Karakorum Highway(KKH) from Besham to Chilas. Intense seismicity, deep gorges, steep terrain and extreme climatic events trigger multiple mountain hazards along the KKH, among which debris flow is recognized as the most destructive geohazard. This study aims to prepare a field-based debris flow inventory map at a regional scale along a 200 km stretch from Besham to Chilas. A total of 117 debris flows were identified in the field, and subsequently, a point-based debris-flow inventory and catchment delineation were performed through Arc GIS analysis. Regional scale debris flow susceptibility and propagation maps were prepared using Weighted Overlay Method(WOM) and Flow-R technique sequentially. Predisposing factors include slope, slope aspect, elevation, Topographic Roughness Index(TRI), Topographic Wetness Index(TWI), stream buffer, distance to faults, lithology rainfall, curvature, and collapsed material layer. The dataset was randomly divided into training data(75%) and validation data(25%). Results were validated through the Receiver Operator Characteristics(ROC) curve. Results show that Area Under the Curve(AUC) using WOM model is 79.2%. Flow-R propagation of debris flow shows that the 13.15%, 22.94%, and 63.91% areas are very high, high, and low susceptible to debris flow respectively. The propagation predicated by Flow-R validates the naturally occurring debris flow propagation as observed in the field surveys. The output of this research will provide valuable input to the decision makers for the site selection, designing of the prevention system, and for the protection of current infrastructure.
基金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.
