Xiangyang Lu , Liang Xu , Wenbin Gao , Yiyan Zhong , Qi Li
{"title":"Sensitivity analysis of key influencing factors for CO2 residual trapping within porous media","authors":"Xiangyang Lu , Liang Xu , Wenbin Gao , Yiyan Zhong , Qi Li","doi":"10.1016/j.advwatres.2025.105047","DOIUrl":"10.1016/j.advwatres.2025.105047","url":null,"abstract":"<div><div>In this study, a pore-scale two-phase flow model is developed by coupling the Navier–Stokes equations with the phase-field method to investigate the effects of wettability, viscosity ratio (<em>M</em>), and drainage–imbibition cycles on CO<sub>2</sub> residual trapping. Numerical simulations reveal that under unfavorable viscosity ratios (<em>M</em> < 1), the CO<sub>2</sub> drainage front becomes unstable, while wettability significantly influences the two-phase drainage patterns. When the CO<sub>2</sub> front breakthrough, strongly water-wet media exhibit lower initial CO<sub>2</sub> saturation and higher inlet drainage pressure. After breakthrough, the inlet pressure in these media shows multiple peaks, with corresponding fluctuations in CO<sub>2</sub> saturation at the pressure maxima. During forced imbibition, CO<sub>2</sub> in the strongly water-wet medium undergoes initial residual trapping followed by transport and escape due to water phase expansion, resulting in a stepwise decline in residual CO<sub>2</sub> saturation and indicating potential secondary transport behavior. The two-phase <em>M</em> in the reservoir environment is not the main factor dominating the drainage pattern and CO<sub>2</sub> residual trapping efficiency. Strongly water-wet systems tend to form sparse, large-scale residual CO<sub>2</sub> clusters, while weakly water-wet systems favor dense, small-scale distributions. An increased <em>M</em> suppresses jamming by enhancing viscous resistance in the non-wetting phase, thereby increasing the average droplet size. Simulations of the drainage-imbibition cycle show that residual CO<sub>2</sub> droplets reduce the effective permeability of the porous medium, elevate the inlet pressure, and delay breakthrough phenomenon, which is beneficial for geological CO<sub>2</sub> sequestration. However, the cycling process increases the escape risk of residual CO<sub>2</sub> phase at the outlet, and the residual CO<sub>2</sub> saturation tends to decrease with the number of cycles. The research results provide theoretical support for the capacity evaluation and risk management of geological CO<sub>2</sub> storage.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"204 ","pages":"Article 105047"},"PeriodicalIF":4.0,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144596699","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}
Achyut Mishra , Hailun Ni , Seyed Ahmad Mortazavi , Ralf R. Haese
{"title":"Performance assessment of graph theory towards predicting fluid flow in rocks across multiple spatial scales","authors":"Achyut Mishra , Hailun Ni , Seyed Ahmad Mortazavi , 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<sub>2</sub> 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, 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":"A dynamic Budyko framework incorporating water storage changes and antecedent precipitation to assess runoff response in typical karst regions","authors":"Hang Chen , Yu Li , 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><sub>0</sub>/<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><sub>0</sub>-based dscomposition framework specially isolating meteorological drivers. This secondary analysis demonstrated that relative humidity dominated runoff variation among climatic variables controlling <em>ET</em><sub>0</sub>. 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 , Yong Zhang , 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 , Wenjuan Zheng , 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}
Kerem Bozkurt , Christoph Lohrmann , Felix Weinhardt , Daniel Hanke , Raphael Hopp , Robin Gerlach , Christian Holm , Holger Class
{"title":"Intermittent flow paths in biofilms grown in a microfluidic channel","authors":"Kerem Bozkurt , Christoph Lohrmann , Felix Weinhardt , Daniel Hanke , Raphael Hopp , Robin Gerlach , Christian Holm , 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, 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}
Jesús Fernández , Amanda C.S.N. Pessoa, Clarice de Amorim, Jorge Avendaño, Márcio S. Carvalho
{"title":"From pore-scale physics to macroscopic flow dynamics in porous media micromodels with different wettability","authors":"Jesús Fernández , Amanda C.S.N. Pessoa, Clarice de Amorim, Jorge Avendaño, 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}
Yuanzhe Wang , Yifei Zong , James L. McCreight , Joseph D. Hughes , Michael Fienen , Alexandre M. Tartakovsky
{"title":"Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation","authors":"Yuanzhe Wang , Yifei Zong , James L. McCreight , Joseph D. Hughes , Michael Fienen , Alexandre M. Tartakovsky","doi":"10.1016/j.advwatres.2025.105024","DOIUrl":"10.1016/j.advwatres.2025.105024","url":null,"abstract":"<div><div>We evaluate the performance of the Karhunen–Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen–Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters.</div><div>We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"203 ","pages":"Article 105024"},"PeriodicalIF":4.0,"publicationDate":"2025-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144335446","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}