Advances in Water Resources最新文献

{"title":"Pore-scale insights on mixed-wettability and its impact on underground hydrogen storage in aquifers","authors":"Mansour Nazari,&nbsp;Hassan Mahani,&nbsp;Shahab Ayatollahi","doi":"10.1016/j.advwatres.2025.105044","DOIUrl":"10.1016/j.advwatres.2025.105044","url":null,"abstract":"<div><div>The efficient storage and extraction of hydrogen from underground porous formations are essential for advancing a sustainable energy transition. However, the pore-scale effects of wettability distribution on hydrogen displacement, trapping, and recovery remain insufficiently characterized. Spatial distribution of wettability (contact angle)—not just the average value—strongly influences fluid configurations and flow pathways, and neglecting these differences can lead to unrealistic storage predictions. This study employs direct numerical simulations to systematically examine how uniform, random, and correlated wettability distributions affect compressible flow of hydrogen in porous media. Results show that decreasing wettability toward a less water-wet state increases hydrogen saturation during drainage across nearly all capillary numbers, enhancing pore space utilization. Moreover, it is found that drainage is primarily governed by invasion percolation, while imbibition is dominated by I₁ and I₂ mechanisms. In strongly water-wet systems, optimal recovery occurs at low capillary numbers (on order of 10⁻⁷) in drainage, and high capillary numbers (on order of 10⁻⁵) in imbibition. Weakly water-wet systems require higher capillary numbers (on order of 10⁻⁶ in drainage, and 10⁻⁴ in imbibition) to suppress capillary fingering. Random wettability distributions reduce hydrogen recovery compared to uniform systems, highlighting the negative impact of wettability heterogeneity. In correlated wettability models, where wettability is linked to pore size, hydrogen trapping significantly increases during imbibition despite similar saturation levels in drainage. These findings enhance our understanding of how mixed wettability affects hydrogen storage and recovery in aquifers, revealing novel insights for identifying suitable reservoirs through wettability classification and predicting flow dynamics.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"204 ","pages":"Article 105044"},"PeriodicalIF":4.0,"publicationDate":"2025-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144633098","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Performance assessment of graph theory towards predicting fluid flow in rocks across multiple spatial scales","authors":"Achyut Mishra ,&nbsp;Hailun Ni ,&nbsp;Seyed Ahmad Mortazavi ,&nbsp;Ralf R. Haese","doi":"10.1016/j.advwatres.2025.105045","DOIUrl":"10.1016/j.advwatres.2025.105045","url":null,"abstract":"<div><div>Fluids tend to migrate along preferential flow paths in rocks which depend on several factors including pore size distribution at sub-Darcy scales and heterogeneity in flow and petrophysical properties in Darcy scale domains. Typically, fluid migration in rocks across scales is determined using numerical simulations which can be computationally expensive. This study presents a graph theory based reduced physics approach to determine potential fluid flow pathways in pore and continuum scale domains. This is an improvement over the existing methods based on graph theory which have largely been focussed on the analysis of the minimum resistance faced by fluid when flowing in rocks. The presented method utilises rock properties as edge weights which are then used to calculate the probability of fluid invading a particular node in the domain. Predictions from the algorithm were validated against results from full physics numerical simulations at pore scale as well as experimentally measured data at pore, core and sand-tank scales. Some of the datasets used for this exercise include Ketton and Estaillades Limestones, and Berea and Bentheimer Sandstones. Further, the algorithm was applied on a suite of reservoir models of the Parasequence-2 of the Paaratte Formation, Otway Basin, Australia. This application was aimed at assessing the influence of grid size resolution and rock property distribution on the uncertainty in CO₂ migration. The algorithm showed computational advantages as it was possible to achieve model runs on some scenarios which are typically not possible using conventional numerical simulations.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"204 ","pages":"Article 105045"},"PeriodicalIF":4.0,"publicationDate":"2025-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144563410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Semi-analytical solutions for flow with space or time-varying viscosity in fractured porous media","authors":"Benjamin Belfort,&nbsp;Anis Younes","doi":"10.1016/j.advwatres.2025.105041","DOIUrl":"10.1016/j.advwatres.2025.105041","url":null,"abstract":"<div><div>Some polymer solutions can form a gel when injected into a porous medium. These gel-forming solutions are of great interest for specific applications such as improving the geotechnical characteristics of soils and enhancing the recovery of contaminants or oil from fractured reservoirs. During its travel through porous media, the viscosity of the gelling solution increases, which can significantly affect its hydraulic properties. This work considers the injection of a gelling polymer, whose viscosity varies in space or time, into a system of two fractures with different apertures. Semi-analytical solutions are developed to describe the evolution of velocity as well as the position of the gelling solution front inside each fracture. The developed semi-analytical solutions are useful to investigate the effect of different parameters as demonstrated by the global sensitivity analysis carried out. GSA results show that fluid velocity inside fractures has different influencing parameters over time. Increasing the viscosity of the gelling fluid reduces the gap between the arrival times of the gelling fluid passing through the large and thin fractures.</div><div>The semi-analytical solutions can also serve for the validation of numerical models simulating variable viscosity flow in fractured media, as we demonstrate with an in-house code.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"204 ","pages":"Article 105041"},"PeriodicalIF":4.0,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144566169","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An integrated model of scale inhibitor (DETPMP) transport, coupled adsorption/precipitation (Γ/Π) and reaction in carbonate systems","authors":"M. Kalantari Meybodi,&nbsp;K.S. Sorbie,&nbsp;O. Vazquez,&nbsp;E.J. Mackay","doi":"10.1016/j.advwatres.2025.105042","DOIUrl":"10.1016/j.advwatres.2025.105042","url":null,"abstract":"&lt;div&gt;&lt;div&gt;Scale inhibitors (SI) are the most widely used materials to prevent inorganic scale formation in flow assurance. SIs operate at threshold concentrations to prevent the formation of inorganic scale; therefore, it is paramount to maintain these levels of SI as long as possible. SIs are usually directly injected (i.e. “bullheaded”) in a “squeeze” treatment into a producing well at very high concentration, where they are retained by adsorption (Γ) and precipitation (Π) mechanisms in the near-well rock formation. Retained SIs are gradually released to the aqueous production stream when the well is put back into production. The efficiency of the SI treatment will depend on the level and type of retention (Γ/Π) in the system. Therefore, a detailed model of the system chemistry and fluid dynamics is required to build a realistic design and for the implementation of more efficient SI field treatments.&lt;/div&gt;&lt;div&gt;The modelling of squeeze treatments in reactive formations requires a detailed model that can consider the SI chemical system, its binding to divalent cations (Ca&lt;sup&gt;2+&lt;/sup&gt; and Mg&lt;sup&gt;2+&lt;/sup&gt;), the precipitation (Π) of the formed SI_M&lt;sup&gt;2+&lt;/sup&gt; complexes, and the adsorption of the solution SI species, coupled to the reactive carbonate rock chemistry. No such model currently exists in the literature which models &lt;em&gt;all&lt;/em&gt; parts of this process, and this is addressed in this study. This paper presents a fully integrated geochemical transport model that can simulate the SI squeeze treatments in reactive formations, such as carbonates, or core flood tests used to support the design of such treatments. The general model is developed to simulate the reactive transport of multiple components, including the SI itself through carbonate formations while reactions are occurring between the bulk fluid and rock matrix, as well as the homogenous reactions in the bulk fluid. For modelling purposes, the SI is treated as a weak polyacid (H&lt;sub&gt;n&lt;/sub&gt;A), which on dissociation yields phosphonate ions along with the neutral species; &lt;em&gt;n&lt;/em&gt; = 10 for DETPMP modelled in this work, but it varies for other phosphonate. The reaction of the SI and carbonate substrate (here taken simply as calcite, CaCO&lt;sub&gt;3&lt;/sub&gt;), is fully described by a recently developed geochemical model for the SI-brine-carbonate system (Kalantari Meybodi et al., 2024a), which considers all the reactions occurring in such system to fully characterize the concentration of all engaged species under any conditions (Kalantari Meybodi et al., 2024a). Moreover, the reaction model includes the new concept of the coupled adsorption/precipitation (Γ/Π) isotherm, which can be used to determine the quantitative partitioning of SI into the adsorption and/or precipitation phases (Kalantari Meybodi et al., 2024a). The isotherm is constructed based on SI static bottle tests, measuring “apparent adsorption”, which are common in any SI application. In the geochemical model, separate parts of th","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"204 ","pages":"Article 105042"},"PeriodicalIF":4.0,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144632203","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A dynamic Budyko framework incorporating water storage changes and antecedent precipitation to assess runoff response in typical karst regions","authors":"Hang Chen ,&nbsp;Yu Li ,&nbsp;Lihua Chen","doi":"10.1016/j.advwatres.2025.105043","DOIUrl":"10.1016/j.advwatres.2025.105043","url":null,"abstract":"<div><div>Karst regions, characterized by dual-domain aquifers and rapid water infiltration, pose challenges when applying traditional hydrological models designed for stable catchments. Understanding these dynamics is critical for advancing water resource management in ecologically fragile areas. The traditional Budyko framework fails to capture the dynamic processes at finer temporal scales that are more acute in karst systems. This study addresses this limitation by proposing a dynamic Budyko framework that explicitly incorporates changes in soil water storage, groundwater storage and antecedent precipitation as key drivers of water availability. Through a multi-temporal analysis of eight karst basins in Southwest China (1980–2020), we demonstrated that the modification effectively constrained annual and monthly data within the Budyko domain, particularly resolving the phenomenon where the monthly distributions of <em>ET</em>/<em>P</em> and <em>ET</em>₀/<em>P</em> significantly exceeded the Budyko boundaries. Attribution analysis showed that runoff was most sensitive to rainfall, followed by landscape characteristics and potential evaporation. To further dissect climatic influences, we proposed a <em>ET</em>₀-based dscomposition framework specially isolating meteorological drivers. This secondary analysis demonstrated that relative humidity dominated runoff variation among climatic variables controlling <em>ET</em>₀. Furthermore, the monthly variation of Budyko parameter ω in karst regions exhibited distinct trends compared to non-karst regions, which was likely attributed to region’s unique hydrogeological structures. This study offers an efficient tool for karst systems, enabling the precise separation of climatic and anthropogenic impacts across temporal scales. These findings demonstrate the applicability of Budyko’s theory and provide critical insights into sustainable water management under intensifying climatic extremes.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105043"},"PeriodicalIF":4.0,"publicationDate":"2025-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144536072","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A conditioned backward fractional-derivative model for identifying super-diffusive pollutant source in aquatic systems with multiple observation data","authors":"Feng Zhang ,&nbsp;Yong Zhang ,&nbsp;HongGuang Sun","doi":"10.1016/j.advwatres.2025.105039","DOIUrl":"10.1016/j.advwatres.2025.105039","url":null,"abstract":"<div><div>Backward models for super-diffusion have been developed to identify pollutant source locations, but they are limited to a single observation and disregard field-measured concentrations. To overcome these limitations, this study derives the adjoint of the space-fractional advection–dispersion equation, incorporating measured concentrations from multiple observation data. Backward probabilities, such as the backward location probability density function (PDF), describe the likely source location(s) at a fixed time before sampling, offering a comprehensive modeling approach for source identification. By applying Bayes’ theorem, the individual PDFs from each observation and its corresponding concentration are combined into a joint PDF, enhancing both the information and reliability compared to the previous single PDF. Field applications show that the improved model enhances accuracy (with PDF peak locations closer to the actual source) and precision (with reduced variance) of backward PDFs for identifying point sources in a natural river and aquifer. The model’s performance is affected by observation count and measurement errors, with double peaks in the backward location PDF possible due to source mass uncertainty. Future refinements, such as incorporating backward travel time analysis and extending applications to reactive pollutants, could further enhance the utility of the conditioned backward fractional-derivative model developed in this study.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105039"},"PeriodicalIF":4.0,"publicationDate":"2025-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144503858","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Vertical textural contract (VTC) increases seawater evaporation from heterogeneous sand","authors":"Jing Yan ,&nbsp;Wenjuan Zheng ,&nbsp;Yan Jin","doi":"10.1016/j.advwatres.2025.105040","DOIUrl":"10.1016/j.advwatres.2025.105040","url":null,"abstract":"<div><div>The coupled soil texture and salinity effect on evaporation in porous systems is complex. This study investigates how vertical textural contrast (VTC) influences evaporation of artificial seawater (SW) compared to freshwater/deionized (DI) water. We conducted experiments using Hele-Shaw cells packed with a fine sand, a coarse sand, and their mixture, respectively, and monitored changes in evaporation rate and water/salt solution dynamics over 20 days. Using light transmittance to detect 2D distribution of water content and surface salt precipitation, we determined how fine/coarse textural contrast affects water transport and salt accumulation. Our findings demonstrated that VTC prolonged high-level water saturation through lateral water transfer from coarse to fine sand, which increased the evaporation rate in both SW and DI water systems. This lateral transfer was driven by differences in capillary pressure between coarse and fine sands, and SW was found to have limited effects on changing the capillary pressure difference compared to DI water systems. Notably, the increase in evaporation due to VTC was more pronounced in SW systems compared to DI systems, which was attributed to the delayed onset of salt precipitation in the fine fraction of heterogeneous chambers compared to homogeneous systems. Our results clarify the intricate mechanisms underlying saltwater evaporation in porous media and emphasize the critical control of soil textural heterogeneity and salinity in evaporation dynamics.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105040"},"PeriodicalIF":4.0,"publicationDate":"2025-06-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144503859","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Intermittent flow paths in biofilms grown in a microfluidic channel","authors":"Kerem Bozkurt ,&nbsp;Christoph Lohrmann ,&nbsp;Felix Weinhardt ,&nbsp;Daniel Hanke ,&nbsp;Raphael Hopp ,&nbsp;Robin Gerlach ,&nbsp;Christian Holm ,&nbsp;Holger Class","doi":"10.1016/j.advwatres.2025.105018","DOIUrl":"10.1016/j.advwatres.2025.105018","url":null,"abstract":"<div><div>Biofilms exposed to flow experience shear stress, which leads to a competitive interaction between the growth and development of a biofilm and shearing. In this study, <em>Pseudomonas fluorescens</em> biofilm was grown in a microfluidic channel and exposed to forced flow of an aqueous solution of variable velocity. It can be observed that under certain conditions preferential flow paths form with a dynamic, but quasi-steady state interaction of growth, detachment, and re-attachment. We find that the regimes for preferential flow path development are determined by nutrient availability and the ratio of shear stress versus the biofilm’s ability to resist shear forces. The intermittent regime of flow paths is mainly driven by the supply with nutrients, which we confirm by comparison with a numerical model based on coarse-grained molecular dynamics and Lattice Boltzmann hydrodynamics.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105018"},"PeriodicalIF":4.0,"publicationDate":"2025-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144335471","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Homogenization-informed convolutional neural network to predict permeability and dispersion in porous media","authors":"Ross M. Weber,&nbsp;Ilenia Battiato","doi":"10.1016/j.advwatres.2025.105022","DOIUrl":"10.1016/j.advwatres.2025.105022","url":null,"abstract":"<div><div>Understanding the transport properties of fluids through porous media is crucial in a wide range of scientific and engineering applications. Accurately predicting key parameters, such as permeability and effective dispersion, is essential for optimizing these processes. These parameters depend not only on the pore-scale geometry but also on flow conditions, and are traditionally expensive to compute since they are generally determined by solving direct numerical simulations on macroscopic pore-scale domains. Such computational costs limit the effectiveness of data-driven approaches in terms both of predictive accuracy and/or types of geometries that can be accurately handled. This is because the computational cost for training over a broad set of topologies and dynamic conditions is prohibitive. In this work, we propose an approach that combines deep learning with multiscale modeling techniques, and exploits the computational efficiency of homogenization theory for periodic domains to support a data-driven technique. By using only a unit cell for training purposes, we are able to generate a large dataset of porous media images and corresponding permeability and dispersion tensors at a significantly reduced computational cost, while spanning an unprecedented range of the geometric and dynamic parameter space. The dataset is composed of 10,000 images, is designed to include a wide variety of morphological properties and serves as the training set for a Convolutional Neural Network (CNN) that estimates permeability and dispersion tensors from both microstructural images and input flow conditions described by the Péclet number. The CNN can quickly and accurately characterize effective properties (permeability and dispersion tensors) spanning more than three orders of magnitude for a wide range of pore-scale topologies and flow regimes. These results highlight the potential to enhance porous media characterization and prediction in various fields.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105022"},"PeriodicalIF":4.0,"publicationDate":"2025-06-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144335445","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"From pore-scale physics to macroscopic flow dynamics in porous media micromodels with different wettability","authors":"Jesús Fernández ,&nbsp;Amanda C.S.N. Pessoa,&nbsp;Clarice de Amorim,&nbsp;Jorge Avendaño,&nbsp;Márcio S. Carvalho","doi":"10.1016/j.advwatres.2025.105038","DOIUrl":"10.1016/j.advwatres.2025.105038","url":null,"abstract":"<div><div>Multiphase flow through porous media involves complex pore-scale phenomena. In this work, we investigated water-oil flow in a PDMS-based porous media micromodel under different wettability conditions (water-wet, hybrid-wet, and oil-wet). Using fluorescence microscopy combined with in-situ pressure drop measurements, we analyzed fluid distributions in the pore space and flow intermittency dynamics as a function of capillary number and fractional flow rate. The effect of wettability conditions on the flow dynamics was correlated with the phase distribution in the pore space. Oil-wet media favors the continuity of the oleic phase. In contrast, water-wet media preferentially favors the flow of the aqueous phase through the porous medium, leaving behind a larger amount of oil. These findings provide a detailed perspective on two-phase flow behavior in micromodels of porous media, significantly improving our understanding of the underlying mechanisms that govern the impact of wettability on phase distribution, macroscopic flow pattern and intermittency dynamics.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105038"},"PeriodicalIF":4.0,"publicationDate":"2025-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144321373","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}