Flow Zone Research Articles

A Microscopic Experimental Study on the Dominant Flow Channels of Water Flooding in Ultra-High Water Cut Reservoirs

The water drive reservoir in Shengli Oilfield has entered a stage of ultra-high water cut development, forming an advantageous flow channel for the water drive, resulting in the inefficient and ineffective circulation of injected water. Therefore, the distribution characteristics of water drive flow channels and their controlled residual oil in ultra-high water cut reservoirs are of great significance for treating water drive dominant flow channels and utilizing discontinuous residual oil. Through microscopic physical simulation of water flooding, color mixing recognition and image analysis technology were used to visualize the evolution characteristics of water flooding seepage channels and their changes during the control process. Research has shown that during the ultra-high water content period, the shrinkage of the water drive seepage channel forms a dominant seepage channel, forming a “seepage barrier” at the boundary of the dominant seepage channel, and dividing the affected area into the water drive dominant seepage zone and the seepage stagnation zone. The advantage of water flooding is that the oil displacement efficiency in the permeable zone is as high as 80.5%, and the remaining oil is highly dispersed. The water phase is almost a single-phase flow, revealing the reason for high water consumption in this stage. The remaining oil outside the affected area and within the stagnant flow zone accounts for 89.8% of the remaining oil, which has the potential to further improve oil recovery in the later stage of ultra-high water cut. For the first time, the redundancy index was proposed to quantitatively evaluate the control effect of liquid extraction and liquid flow direction on the dominant flow channels in water flooding. Experimental data showed that both liquid extraction and liquid flow direction can regulate the dominant flow channels in water flooding and improve oil recovery under certain conditions. Microscopic physical simulation experiments were conducted through the transformation of well network form in the later stage of ultra-high water content, which showed that the synergistic effect of liquid extraction and liquid flow direction can significantly improve the oil recovery effect, with an oil recovery rate of 68.02%, deepening the understanding of improving oil recovery rate in the later stage of ultra-high water content.

Integrated Subsurface Hydrologic Modeling for Agricultural Management Using HYDRUS and UZF Package Coupled with MODFLOW

The present work aims to compare two different subsurface hydrological models, namely HYDRUS and MODFLOW UZF package, in terms of groundwater recharge; thus, both models were coupled with MODFLOW. The study area is an experimental kiwifruit orchard located in the Arta plain in the Epirus region of Greece. A novel conceptual framework is introduced in order to (i) use in situ and laboratory measurements to estimate parameter values for both sub-surface flow models; (ii) couple the developed models with MODFLOW to estimate groundwater recharge; and (iii) compare and evaluate the performance of both approaches, with differences stemming from the distinctive equations describing the flow in the unsaturated zone. Detailed soil investigation was conducted in two soil horizons in the research field to identify soil texture zones, along with infiltration experiments implementing both double-ring and single-ring infiltrometers. The results of the field measurements indicate that fine-textured soils are predominant within the field, affecting several hydrological processes, such as infiltration, drainage, and root water uptake. Field measurements were incorporated in unsaturated zone flow modeling and the infiltration fluxes were simulated with the application of both the UZF package of MODFLOW and the HYDRUS code. The two codes presented acceptable agreement between the simulated and observed hydraulic head values with a similar performance in terms of statistics; however, they produced different results regarding recharge rates in the aquifer as simulated by MODFLOW. HYDRUS produced higher hydraulic head values in the aquifer throughout the simulation, related to higher recharge rates arising from the root water uptake and the capillary effects that are computed by HYDRUS but neglected by the UZF package of MODFLOW.

IMPROVING THE RELIABILITY OF THE THREE-DIMENSIONAL GEOLOGICAL BASIS OF COMPLEX DEVELOPMENT TARGETS

This article proposes a new method for identifying different groups of reservoirs when constructing a three-dimensional geological and hydrodynamic model of the target. Productive formations of the Bajocian and Bathonian stages of the Jurassic system of the Kalamkas field were selected as the object of study. This field is tectonically associated with the northern part of the North Buzachinsky arch, located at the northwestern end of the Turan plate. It was noted that these productive deposits have a complex geological structure, due in lithological terms to uneven alternation of sand-siltstone and clay rocks. In addition, four different reservoir classes are identified within each tier, which are characterized by their own permeability-porosity correlation relationships. As a result, for reliable construction of the filtration field, it is necessary to determine the distribution zones of each of the identified reservoir classes within the reservoir. This is due to the fact that due to the different filtration field, the development of reserves from each group of reservoirs will occur unevenly both in area and in section. In this regard, the integrated Flow Zone Indicator (FZI) was evaluated using the results of laboratory core studies. In order to obtain a continuous FZI parameter for wells without core sampling and laboratory analysis, the results of laboratory studies were compared with the results of well geophysical studies (well logging). As a result, analytical dependence of FZI and log data was obtained. Based on this analytical dependence, new three-dimensional models were constructed and zones of distribution of different groups of reservoirs were identified. To verify the obtained results, the synthetic FZI curve calculated from the obtained analytical relationship was compared with the FZI estimate based on the results of laboratory core studies. It is concluded that the synthetic FZI curve is consistent with the FZI estimate based on the results of laboratory core studies and, therefore, the proposed approach can be used when performing work related to 3D modeling.

Study on the Control Effect of Borehole Gas Extraction in Coal Seams Based on the Stress–Seepage Coupling Field

In order to determine the reasonable parameters of high-gas and extra-thick coal seam drainage, considering the factors of the coal seam metamorphic degree, stress condition, gas occurrence state, and permeability dynamic change, the gas desorption, diffusion, and transport process of coal seam gas are analyzed. A secondary distribution model of coal around the borehole, a porosity variation model of coal around the borehole, a stress–seepage coupling model, a pore flow model of the pressure-driven transition flow zone, and a free molecular flow zone are established. Taking the gas drainage of Zhangcun Coal Mine of Lu’an Group as the research object, the influence of drilling hole diameter, coal seam permeability, gas original pressure, and other factors on the control range of coal seam drainage drilling is simulated by ANSYS Fluent 6.3.26. The results show that secondary stress distribution occurs in the coal seam drill hole under the action of lead stress, which leads to the change in porosity; the seepage zone, transition zone, molecular flow zone, and original rock stress zone are presented around the drill hole; and the range of influence of the drill hole is mainly based on the seepage zone and the transition zone, supplemented by the molecular flow zone. The control range of the drill hole is in a positive proportional relationship to the diameter of the drill hole, the porosity of the coal seam, and the original pressure of the gas.

Time-lapse 3D imaging and modelling of acid-rock interactions within a mudstone: Implications for cap-rock behaviour during carbon storage

Assessing caprock integrity is crucial for safe geological carbon storage. This study addresses this by utilising time-lapse micro-CT imaging to capture real-time, dynamic changes in mudstone during acid injection, providing a detailed analysis of the 3D geometry and evolution of minerals and fractures. Additionally, the integration of reactive transport modelling offers insights into the potential changes occurring immediately after acid injection, enhancing our understanding of fluid-rock interactions and their impact on caprock stability. Different acid concentrations were used to mimic a range of possible acid concentrations during CO2 injection and storage. Changes in mudstones due to acid interaction include initial pre-existing fracture closure, followed by fracture growth and sample swelling. Acid predominantly followed pathways like fractures or laminations. Reactive transport models showed most dissolution occurring near the inlet, driven by a decrease in H+ concentration along fluid pathways. Carbonate distribution controlled dissolved areas, with calcite dissolution and acid migration rates decreasing over time. Simulation indicated reduced shear stress after acid injection, especially with low-pH acids. This correlates to decreased fluid velocity, resulting in slower fluid flow. Slower acid movement and prolonged presence in immobile fluid flow zones led to more extensive calcite dissolution compared to mobile zones. This study provides valuable insights into microstructural changes in caprock following CO2 injection, emphasizing the potential risk of fracture development and leakage.

Flow behavior and heat transfer characteristics of liquid film on vertical hot surface by inclined jet impingement

In this study, the flow behavior and heat transfer characteristics of the liquid film on the hot wall by inclined jet impingement are experimentally studied in detail. The effects of jet parameters, such as jet inclination angle, impingement distance, jet Reynolds number (Rej), and nozzle diameter are explored. The liquid film flow is visualized using a high-speed camera, and the surface temperature and heat flux are obtained by solving the inverse heat conduction problem. The results indicate that the liquid film shape is strongly affected by jet inclination angle but is almost unaffected by other parameters. As the inclination angle increases, the liquid film shape changes from elliptical to fusiform. In addition, the onset, enhancement, and disappearance of boiling cause the expansion and contraction of liquid film. The splashing rate is barely affected by jet parameters and remains within the range of 80–95 % under all conditions. The propagation of wetting fronts exhibits anisotropy. Except for the impingement distance, the position of wetting fronts along y-axis direction displays a high dependence on all parameters. The Rej and nozzle diameter have a significant effect on the heat flux at the impingement point and parallel flow zone, while the jet inclination angle and impingement distance only effect the impingement point. The visualization image proves that the droplet impingement pattern is the main reason for the increase in heat flux at higher impingement distances. Optimizing jet parameters can promote the wall to enter the rapid cooling stage in advance and increase maximum heat flux, thereby improving the maximum cooling capacity of the jet. Finally, an empirical equation is proposed to predict the maximum Nusselt number.

The Zone of Proximal Flow in K-12 Creative Collaborative eLearning: A synthesis of the zone of proximal development and creative flow

Student groups coexist on individual and team levels toward learning symmetry within two zones: the zone of proximal development (ZPD) for new knowledge co-construction and the zone of proximal flow (ZPF), which occurs once a balance between the learner?s skill and challenge of the task emerges. Thus, each group functions within the ZPF, especially for K-12 students engaging in social behavior and collaborative learning activities. eLearning must support the dual persona of learners as users engage and remain within the ZPF via diverse pedagogical methodologies and tools, such as multimodal learning analytics, to support and enrich students? learning curves in self-organized and collaborative eLearning.

Gene Flow Between Populations With Highly Divergent Mitogenomes in the Australian Stingless Bee, Tetragonula hockingsi.

Coadaptation of mitochondrial and nuclear genes is essential for proper cellular function. When populations become isolated, theory predicts that they should maintain mito-nuclear coadaptation in each population, even as they diverge in genotype. Mito-nuclear incompatibilities may therefore arise when individuals from populations with divergent co-evolved mito-nuclear gene sets are re-united and hybridise, contributing to selection against inter-population hybrids and, potentially, to speciation. Here, we explored genetic divergence and gene flow between populations of a stingless bee (Tetragonula hockingsi) that have highly divergent mitogenomes. We identified three distinct populations across the species' 2500 km range on the east coast of Queensland (Australia): 'Cape York', 'Northern', and 'Southern'. The mitogenomes of each population showed > 12% pairwise nucleotide divergence from each other, and > 7% pairwise amino acid divergence. Based on nuclear SNPs from reduced representation sequencing, we identified at least two zones of gene flow between populations: a narrow natural zone between Northern and Southern populations (coinciding with a biogeographic barrier, the Burdekin Gap), and an artificial zone at the southern edge of the species' distribution, where Cape York, Northern, and Southern mito-lineages have been brought together in recent decades due to beekeeping. In the artificial hybrid zone, we also confirmed that males of all three mito-lineages were attracted to the mating aggregations of Southern queens, consistent with inter-population hybridisation. Populations of T. hockingsi thus appear to be in the 'grey zone' of the speciation continuum, having strong genetic differentiation but incomplete reproductive isolation. Among the nuclear SNPs most differentiated between Northern and Southern populations, several were associated with genes involved in mitochondrial function, consistent with populations having co-diverged mito-nuclear gene sets. Our observations suggest that coadapted sets of mitochondrial and nuclear genes unique to each population of T. hockingsi may play a role in maintaining population boundaries, though more study is needed to confirm the fitness costs of mito-nuclear incompatibilities in hybrid individuals.

Modelling Seasonal Variability in Parameters Defining Volumetric Water Content in a Low Permeability Soil in Central Illinois: An Application of MODFLOW‐6 and the Unsaturated Zone Flow Package

ABSTRACTIncreasing interest in solute transport phenomena in agricultural systems on a sub‐annual basis necessitates a better understanding of seasonal changes in natural systems and how these changes can be incorporated into modelling. A better understanding of the seasonal timing of nutrient loading in tile drained agricultural systems in particular is essential for efforts trying to replicate or predict the occurrence of harmful algal blooms. Literature exists showing there are seasonal dynamics (freeze–thaw, plant‐root processes, land management practices, etc.) that may cause changes in the hydraulic properties of the soil zone including hydraulic conductivity and porosity. To test whether these changes are important in an agricultural system, a MODFLOW‐6 model using the unsaturated zone flow package was constructed. The simulation was comprised of separate, seasonal models to be run sequentially with each year being broken into a winter and summer seasons. As part of this architecture, model parameters representing soil hydraulic properties were allowed to vary by season. The model was calibrated against soil moisture observations at multiple depths using a genetic algorithm machine learning technique. The parameters of the sub‐models were compared for the winter and summer seasons. Brook‐Corey epsilon, saturated vertical conductivity, saturated volumetric water content and residual volumetric water content were found to be consistently different between the modelled summer and winter periods. A more traditional model which did not allow hydraulic properties to vary seasonally was also run and compared to the seasonal architecture and the seasonal architecture was found to improve simulation results. The hydrologic dynamics of the unsaturated zone—particularly in tile drained agricultural systems—control the residence time for water and solutes, which is critical for in‐field chemical processes such as denitrification. This work has important implications for seasonal transport phenomena in agricultural systems and improving the simulation and prediction of harmful algal blooms.

Estimation of the Hydraulic Parameters of a Stratified Alluvial Soil in the Region of El Haouareb—Central Tunisia. Experiments, Empirical, Analytical and Inverse Models

Modeling water flow and contaminant transport in the unsaturated zone is a difficult task that relies heavily on good hydrodynamic soil characterization. This article presents a complementarity between experimentation, direct modeling and inverse modeling in order to provide a better estimate of the hydrodynamic parameters of stratified alluvial soil in the El Haouareb region of the Kairouane plain in Tunisia. A field sampling campaign was carried out. The samples collected underwent particle size analysis, bulk density measurements and infiltration tests using a mini-Muntz. In parallel, simple evaporation tests were applied to separate strata. In addition, a 2 m soil column was reconstituted and fitted with sensors to monitor water content, tension, temperature and electrical conductivity. An internal drainage test was performed on this monolith. Three methods were applied using experimental data to estimate soil hydrodynamic parameters. In the first method, pedotransfer functions were used (Rosetta platform) based on granulometric results and bulk density. In the second, water tension and water content monitored during the simple evaporation test were used to plot the soil–water retention curve (SWRC) using SWRC-Fit. In the third method, inverse modeling was applied to the internal drainage test. A comparison of the results showed that the inverse method had the lowest RMSE. Uncertainty analysis has been implemented for both the experimental and numerical set up.

Hybridization has localized effect on genetic variation in closely related pine species

BackgroundHybridization is a known phenomenon in nature but its genetic impact on populations of parental species remains less understood. We investigated the evolutionary consequences of the interspecific gene flow in several contact zones of closely related pine species. Using a set of genetic markers from both nuclear and organellar genomes, we analyzed four hybrid zones (384 individuals) and a large panel of reference allopatric populations of parental taxa (2104 individuals from 96 stands).ResultsWe observed reduced genetic diversity in maternally transmitted mitochondrial genomes of pure pine species and hybrids from contact zones compared to reference allopatric populations. The distribution of mtDNA haplotypes followed geographic rather than species boundaries. Additionally, no new haplotypes emerged in the contact zones, instead these zones contained the most common local variants. However, species diverged significantly at nuclear genomes and populations in contact zones exhibited similar or higher genetic diversity compared to the reference stands. There were no signs of admixture in any allopatric population, while clear admixture was evident in the contact zones, indicating that hybridization has a geographically localized effect on the genetic variation of the analyzed pine species.ConclusionsOur results suggest that hybrid zones act as sinks rather than melting pots of genetic diversity. Hybridization influences sympatric populations but is confined to contact zones. The spectrum of parental species ancestry in hybrids reflects the old evolutionary history of the sympatric populations. These findings also imply that introgression may play a crucial role in the adaptation of hybrids to specific environments.

Hemodynamic characteristics of pulsatile blood flow through bifurcated stenosed carotid artery

PurposeCarotid artery is often associated with plaque deposition because of its shape and associated flow features. The shape of stenosed bifurcation is characterised by bifurcation angle (ß), planarity angle (α) and severity of stenosis (b). In the present work, three-dimensional numerical computations have been performed to analyse the effect of these geometrical parameters of carotid bifurcation on the characteristics of flow.Design/methodology/approachGoverning equations of this study were solved using ANSYS Fluent 20.1 and the blood flow was considered as laminar, pulsatile and non-Newtonian. Instantaneous flow behaviour has been illustrated using vorticity, velocity and helicity contours, whereas the time-averaged wall shear stress (τw¯) and oscillatory shear index (OSI) quantify the time-averaged behaviour.FindingsThe recirculation zone and secondary flow are ascertained to be stronger for higher bifurcation angle as compared to the lower bifurcation angle. Strength of the secondary flow is found to reduce with increase in α from 0° to 10°, whereas it grows as α varies from 10° to 20°. For higher bifurcation angles, τw¯ is lower than 2 Pa and OSI is greater than 0.2 on the outer walls. Similar observations were made for τw¯ and OSI distribution on bottom wall in non-planar cases, which predicted atherogenic locations.Originality/valueThe values for ß were taken as 30°, 45°, 60° and 75°, whereas for α, range of 0°–20° was chosen. The stenosis was considered on the outer wall of internal carotid artery and its severity was considered within the range of 0%–60%.

Effects of mono- and multicomponent nonaqueous-phase liquid on the migration and retention of pollutant-degrading bacteria in porous media

The successful implementation of in-situ bioremediation of nonaqueous-phase liquid (NAPL) contamination in soil-groundwater systems is greatly influenced by the migration performance of NAPL-degrading bacteria. However, the impact and mechanisms of NAPL on the migration/retention of pollutant-degrading bacteria remain unclear. This study investigated the migration/retention performance of A. lwoffii U1091, a strain capable of degrading diesel while producing surfactants, in porous media without and with the presence of mono- and multicomponent NAPL (n-dodecane and diesel) under environmentally relevant conditions. The results showed that under all examined conditions (5 and 50 mM NaCl solution at flow rates of 4 and 8 m/d), the presence of n-dodecane/diesel in porous media could reduce the migration and enhance retention of A. lwoffii in quartz sand columns. Moreover, comparing with mutlicomponent NAPLs of n-dodecane, the monocomponent NAPLs (diesel) exhibited a greater reduction effect on the retention of A. lwoffii in porous media. Through systemically investigating the potential mechanisms via tracer experiment, visible chamber experiment, and theoretical calculation, we found that the reduction in porosity, repulsive forces and movement speeds, the presence of stagnant flow zones in porous media, particularly the biosurfactants generated by A. lwoffii contributed to the enhanced retention of bacteria in NAPL-contaminated porous media. Moreover, owing to presence of the greater amount of hydrophilic components in diesel than in n-dodecane, the available binding sites for the adsorption of bacteria were lower in diesel, resulting in the slightly decreased retention of A. lwoffii in porous media containing diesel than n-dodecane. This study demonstrated that comparing with porous media without NAPL contamination, the retention of strain capable of degrading NAPL in porous media with NAPL contamination was enhanced, beneficial for the subsequent biodegradation of NAPL.

Exploring comparative heterogeneity management for precise water saturation assessment in carbonate formations: Dean-Stark measurements and beyond

Accurate determination of water saturation is a critical aspect of reservoir characterization and development potential of a hydrocarbon field. The inherent heterogeneity in carbonate reservoirs significantly influences the Archie equation coefficients and, consequently, water saturation calculations. In this study, 160 direct core sample water saturation measurements using the Dean-Stark method have been used. The primary goal was to identify the most suitable reservoir heterogeneity management approach for water saturation calculation. The research focuses on the Dalan and Kangan formations (equivalent to Khuff Formation), known as some of the most hydrocarbon-rich (gas) reservoirs globally, located in the western part of the Persian Gulf. To manage reservoir heterogeneity, we employed various techniques, including Winland R35, flow zone indicator, Lucia, pore types, and electrofacies, utilizing porosity and permeability core data, thin section studies, and well-logging data to classify the studied intervals into more homogeneous categories. Subsequently, we measured Archie parameters and, for the first time, water saturations derived from Dean-Stark measurements were compared with Archie saturations from various heterogeneity management methods. Additionally, the influence of geological and petrophysical parameters on the accuracy of water saturation calculations in different rock types is discussed. Our results highlighted the significance of geometrical attributes of pore types, pore throat radius, and permeability as key factors affecting the precision of water saturation calculations. Moreover, our investigation revealed that Lucia's rock typing methodology exerted a substantial influence on the control of permeability and the distribution of water saturation within the reservoir. Our analysis revealed that the electrofacies method provided the most precise estimates for water saturation, with a mean difference of 0.07 (v/v) units compared to the Dean-Stark method. Subsequently, the accuracy of predicted water saturation in rock types determined by the Lucia method, with a mean error of 0.09 (v/v), ranked second among the rock typing methods. Winland R35 and flow zone indicator methods exhibited the highest uncertainty in water saturation estimation, with mean differences of 0.23 and 0.21 (v/v), respectively, compared to the Dean-Stark method.

Analysis of droplet motion in cavity zone of rotating packed bed

Droplet motion in outer cavity zone of rotating packed bed is a complex phenomenon that requires analysis of liquid and gas flow to clarify influences of factors such as rotation speed, liquid flow rate, or gas flow rate. In this study, hydrodynamics is described using laser Doppler velocimetry and high-resolution visualization for droplet motion characterization. The results reveal that droplet motion initially depends on the peripheral velocity of the packing, while further from the packing, it loses kinetic energy mainly due to the drag force caused by the surrounding gas. Two-way coupling is investigated as both phases interact with each other. Rankine vortex theory is applied to simulate gas flow in the outer cavity zone. A model, based on Newton’s second law, is developed to predict droplet dynamics, including impingement velocity, trajectory length, and total residence time. Validation of the model shows that it can predict the liquid velocity in the outer cavity zone accurately. Upon impingement of droplets on the casing, part of the liquid rebounds as secondary droplets into the cavity, where they circulate with the gas flow. Overall, this study provides insights into droplet dynamics in the rotating packed bed and offers a basis for optimizing process efficiency and design in various applications, including chemical processing and environmental engineering.

Induction of Controllable Vortical Flow in a Dual-Stenosis Aorta Model: A Replication of Disordered Eddies Flow in Aneurysms.

This paper presents a two-stenosis aorta model mimicking vortical flow in vascular aneurysms. More specifically, we propose to virtually induce two adjacent stenoses in the abdominal aorta to develop various vortical flow zones post stenoses. Computational fluid dynamics (CFD) simulations were conducted for the virtual two-stenosis model based on physiological and anatomical data (i.e., diameters, flow rate waveforms) from adult rabbits. The virtual model includes adult rabbits' infra-renal portion of the aorta and iliac arteries. 3D CFD simulations in five different dual-stenosis configurations were performed using a commercial CFD package (FLUENT). In-house software assessed the evolution of flow vortices. Notably, spatial-temporally averaged wall shear stress (STA-WSS) and oscillatory shear index (OSI), the total volume of vortex flow, the number of vortices, and the phase-to-phase overlap of vortex flow within each region were evaluated. In all models, we found consistent patterns of the vortex flow parameters, indicating that the adjacent stenoses induced three different hemodynamic zones, namely, stable vortical flow (after the first stenosis), transient vortical flow (after the second stenosis), and unstable vortical flow (further distal to the second stenosis). Also, different degrees of flow disturbance can be achieved in these three zones. It is significant to note that, although the 'dual-stenosis' geometry is completely hypothetical, it allows us to create various vortical flows in consecutive vessel segments for the first time. As a result, if implemented as a pre-clinical model, the proposed two-stenosis model offers an attractive, tunable environment to investigate the interplays between subject-specific hemodynamics and vascular remodeling. This aspect remains in our future directions.

Sensitivity of surface water and groundwater contributions to streamflow in a tropical glacierized basin under climate change scenarios

Abstract While mountain water faces threats posed by climate change, particularly in snow-dominated and glacierized systems, the role of groundwater in sustaining streamflow in these systems remains elusive. Changing mountain headwaters, marked by reduced snowpacks, retreating glaciers, shifting precipitation patterns, and rising temperatures, pose a crucial question: what is the resilience of streamflow in these mountains, and what role does groundwater play in this resilience? This is particularly uncertain in tropical high mountains where the seasonality of precipitation and glacier melt govern streamflow generation. A glacio-hydrological model was created using the Cold Regions Hydrological Modelling platform to investigate cryosphere-surface water-groundwater interactions in the Quilcayhuanca Basin, in Peru’s Cordillera Blanca. The model was forced by in-situ meteorological observations and parameterized using numerous data sources and process-based studies in the basin. Model results show that during the dry season, 37% of streamflow is generated from groundwater discharge, increasing to 56% during the lowest flows. Evapotranspiration is the largest mass flux from the basin at the peak of the dry season. Precipitation, temperature, and glacier change scenarios were used to assess the sensitivity of basin hydrology to climate change and glacier retreat. In a warmer, wetter, and nearly deglaciated future, Quilcayhuanca basin streamflow is expected to decrease by 4-19% annually, with a larger volumetric change in overland and vadose zone flow than in groundwater flow. The range in values is more closely linked to uncertainty in precipitation change than temperature change. Despite a strong reduction in snow and ice contribution to streamflow with warming and deglaciation, the concomitant increase in precipitation can limit the changes in streamflow and groundwater flow, showcasing the resilience of the system to shifts in climate and glacier cover.

Precision polishing force control technology of abrasive belt-wheel polishing servo control system based on HGO-SMC

Due to cylinder friction, valve dead zone effect, flow nonlinearity, and measurement noise exist in pneumatic system. In the actual machining, it is difficult to achieve precise control of the pneumatic servo polishing system using the traditional PID algorithm, so the machining accuracy of the blades cannot be guaranteed. A sliding mode control (SMC) method based on high gain observer (HGO) is proposed. The HGO observes the output signal of the system and feeds the estimated signal back to the SMC to ensure that the observation error is uniform ultimate boundedness. Simulation and experiments show that compared with the PID algorithm, the steady-state error of HGO-SMC can be reduced by 66.4%. Compared with SMC(sgn), the high-frequency chattering of HGO-SMC can be reduced by 72.5%. Moreover, the polishing processing experiment shows that HGO-SMC can improve the form and position accuracy by 52.6% and reduce the surface roughness by 55.6%, compared with the original PID control algorithm of the polishing machine tool.

Numerical simulation of two-phase oil–water flow in fractured-vuggy reservoirs based on the coefficient of porous medium proportion and coupled regions

Based on the flow characteristics of fluids in various reservoir media, fractured-vuggy oil reservoirs can be classified into seepage zones and conduit flow zones. An interface exists between these two regions, where the movement of formation fluid near this interface is characterized by a coupling or transitional phenomenon between seepage and conduit flow. However, the complexity of this coupling interface poses challenges for traditional numerical simulations in accurately representing the intricate fluid dynamics within fractured-vuggy oil reservoirs. This limitation impacts the development planning and production adjustment strategies for fractured-vuggy oil reservoirs. Consequently, achieving accurate characterization and numerical simulation of these systems remains a critical challenge that requires urgent attention. A new mathematical model for oil-water two-phase flow in fractured-vuggy oil reservoirs is presented, which developed based on a novel coupling method. The model introduces the concept of the proportion coefficient of porous media within unit grids and defines a coupling region. It employs an enhanced Stokes–Brinkman equation to address the coupling issue by incorporating the proportion coefficient of porous media, thereby facilitating a more accurate description of the coupling interface through the use of the coupling region. Additionally, this proportion coefficient characterizes the unfilled cave boundary, simplifying the representation of model boundary conditions. The secondary development on the open-source fluid dynamics software is conducted by using matrix & laboratory (MATLAB). The governing equations of the mathematical model are discretized utilizing finite volume methods and applying staggered grid techniques along with a semi-implicit calculation format for pressure coupling—the Semi-Implicit Method for Pressure Linked Equations algorithm—to solve for both pressure and velocity fields. Under identical mechanism models, comparisons between simulation results from this two-phase flow program and those obtained from Eclipse reveal that our program demonstrates superior performance in accurately depicting flow states within unfilled caves, thus validating its numerical simulation outcomes for two-phase flow in fractured cave reservoirs. Utilizing the S48 fault-dipole unit as a case study, this research conducted numerical simulations to investigate the water-in-place (WIP) behavior in fractured-vuggy oil reservoirs. The primary focus was on analyzing the upward trend of WIP and its influencing factors during production across various combinations of fractures and dipoles, thereby validating the feasibility of the numerical modeling approach in real-world reservoirs. The simulation results indicated that when multiple dissolution cavities at different locations communicated with the well bottom sequentially, the WIP in the production well exhibited a staircase-like increase. Furthermore, as the distance between bottom water and well bottom increased, its effect on water intrusion into the well diminished, leading to a slower variation in the WIP curve. These characteristics manifested as sudden influxes of water flooding, rapid increases in water levels, and gradual rises—all consistent with actual field production observations. The newly established numerical simulation method for fractured-vuggy oil reservoirs quantitatively describes two-phase flow dynamics within these systems, thus effectively predicting their production behaviors and providing guidance aimed at enhancing recovery rates typically observed in fractured-vuggy oil reservoirs.

An improved detached eddy simulation method for cavitation multiphase flow

The evolution of cavitation flow involves intense transient characteristics and complex multi-scale motions, significantly increasing the difficulty and computational cost of solving in separation zones. This study proposes an Improved Detached Eddy Simulation (IDES) method, based on a single-equation eddy-viscosity model, aimed at enhancing the resolution of small-scale cavitation flows. By incorporating the broad-spectrum von Karman scale concept, the method effectively improves the resolution capability for small-scale flows in separated zones. The performance of the IDES model was validated through three computational cases, and the main conclusions are as follows: The results of the triangular cylinder flows show that the IDES model demonstrates a resolution capability for large-scale turbulence comparable to the Large Eddy Simulation (LES) model. The orifice flow results indicate that activating the LES method with insufficient grid resolution is highly discouraged, whereas the IDES model performs more robustly under such conditions. For hydrofoil flows, the IDES model more accurately captures the temporal evolution of cavitation, avoiding the premature cavity shedding predicted by LES under coarse grid conditions. The decomposition results of the vorticity transport equation show that both the dilatation term and viscous term are significant contributors in the vorticity transport process, with their average contributions to vorticity reaching 49.25% and 39.99%, respectively. Vortices can be considered the primary controllers of cavitation dynamics, while the dilatation term and viscous term can be regarded as the main controllers of vorticity. Overall, the energy compensation mechanism of the IDES model, based on local flow information, gives it unique advantages in handling small-scale turbulence and cavitation phenomena, providing both high computational efficiency and accuracy.