Advances in Water Resources最新文献

{"title":"Global sensitivity analysis of mass transfer and reaction dynamics for electrokinetic transport in porous media","authors":"R. Sprocati ,&nbsp;A. Guadagnini ,&nbsp;L. Ceresa ,&nbsp;A. Gallo ,&nbsp;M. Rolle","doi":"10.1016/j.advwatres.2025.104887","DOIUrl":"10.1016/j.advwatres.2025.104887","url":null,"abstract":"<div><div>Electrokinetic (EK) technologies are promising solutions for the remediation of contaminated sites, particularly in low-permeability porous media. However, their widespread adoption is hindered by the challenge of predicting the complex, coupled physico-chemical processes triggered by the application of external electric fields in the subsurface. Numerical models therefore represent essential tools to interpret system behavior. Uncertainties in experimental data, as well as in the formulation of conceptual models, still pose a challenge to develop robust predictive tools. In this context, our work addresses the impact of various sources of uncertainty on model-based predictions of EK transport in porous media. We employ Monte Carlo-based global sensitivity analyses (GSA) within both single-model (SM-GSA) and multi-model (MM-GSA). The multi-model approach relies on a theoretical framework encompassing different models capable of interpreting a set of EK transport scenarios. This allows us to address the impact of model formulation besides parametric uncertainty on mass transfer and reaction dynamics of EK transport. All candidate models in our set are based on a 2D dipole electrode configuration and each model incorporates a different combination of physico-chemical processes to explore different EK remediation scenarios dominated by electromigration or electroosmosis, for both conservative and reactive transport settings. We also investigate the influence of background electrolytes, charge interactions, reactant mobility and degradation reaction kinetics on system dynamics. To overcome the computational burden of process-based modeling and GSA implementations, we develop machine learning-based surrogate models. The latter are employed within both SM- and MM-GSA frameworks, using Sobol’ and AMAE sensitivity indices, respectively. This work provides a comprehensive quantification of how multiple sources of uncertainty impact electrokinetic transport behavior in porous media.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104887"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094795","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A hybrid CNN-transformer surrogate model for the multi-objective robust optimization of geological carbon sequestration","authors":"Zhao Feng ,&nbsp;Bicheng Yan ,&nbsp;Xianda Shen ,&nbsp;Fengshou Zhang ,&nbsp;Zeeshan Tariq ,&nbsp;Weiquan Ouyang ,&nbsp;Zhilei Han","doi":"10.1016/j.advwatres.2025.104897","DOIUrl":"10.1016/j.advwatres.2025.104897","url":null,"abstract":"<div><div>The optimization of well controls over time constitutes an essential step in the design of cost-effective and safe geological carbon sequestration (GCS) projects. However, the computational expense of these optimization problems, due to the extensive number of simulation evaluations, presents significant challenges for real-time decision-making. In this paper, we propose a hybrid CNN-Transformer surrogate model to accelerate the well control optimization in GCS applications. The surrogate model encompasses a Convolution Neural Network (CNN) encoder to compress high-dimensional geological parameters, a Transformer processor to learn global patterns inherent in the well controls over time, and a CNN decoder to map the latent variables to the target solution variables. The surrogate model is trained to predict the spatiotemporal evolution of CO₂ saturation and pressure within 3D heterogeneous permeability fields under dynamic CO₂ injection rates. Results demonstrate that the surrogate model exhibits satisfactory performance in the context of prediction accuracy, computation efficiency, data scalability, and out-of-distribution generalizability. The surrogate model is further integrated with Multi-Objective Robust Optimization (MORO). Pareto optimal well controls are determined based on Non-dominated Sorting-based Genetic Algorithm II (NSGA-II), which maximize the storage efficiency and minimize the induced over-pressurization across an ensemble of uncertain geological realizations. The surrogate-based MORO reduces computational time by 99.99 % compared to simulation-based optimization. The proposed workflow not only highlights the feasibility of applying the CNN-Transformer model for complex subsurface flow systems but also provides a practical solution for real-time decision-making in GCS projects.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"196 ","pages":"Article 104897"},"PeriodicalIF":4.0,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143094797","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":"Modeling and design of air injection-based hydraulic barriers","authors":"Ilan Ben-Noah","doi":"10.1016/j.advwatres.2025.104902","DOIUrl":"10.1016/j.advwatres.2025.104902","url":null,"abstract":"<div><div>Hydraulic barriers are useful for manipulating groundwater flow to mitigate and contain harmful environmental effects. Injecting air into the aquifer has been suggested as a cost-efficient, sustainable, and reversible hydraulic barrier. In this, the injected air reduces the conductance of the aquifer to water flow. However, this practice is not commonly used despite its potential, probably due to a lack of a design tool and field-scale demonstrations. Evaluating the effect of air injection on the water flow is severely hindered by the ability to simulate the unstable coupled multiphase flow problem. Multiphase models are computationally expensive and unstable, require many (sometimes unattainable parameters), and generally provide poor predictive capabilities. Here, a simplified approach is suggested that decouples the phases flows, and uses analytical (fast, stable, parameters-parsimonious) solutions of the air injection problem to evaluate its effect on the water permeability field. In this manuscript, we describe the framework and the solutions of the air flow, discuss its limitations, and demonstrate its usability to evaluate the sensitivity to the media’s and system design parameters in three case studies: (i) point sources in a confined aquifer, (ii) line source in a confined aquifer, and (iii) line source in an unconfined aquifer.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104902"},"PeriodicalIF":4.0,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143158003","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Inferring experimental colloid removal with an inverse two-population model linking continuum scale data to nanoscale features","authors":"Sabrina N. Volponi ,&nbsp;Giovanni Porta ,&nbsp;Bashar M. Al-Zghoul ,&nbsp;Diogo Bolster ,&nbsp;William P. Johnson","doi":"10.1016/j.advwatres.2025.104905","DOIUrl":"10.1016/j.advwatres.2025.104905","url":null,"abstract":"<div><div>Models of colloid transport in porous media that assume constant fractional loss per grain passed fail in the presence of repulsive barriers to attachment, under which condition experiments produce profiles of colloid concentrations with distance from source that are nonexponential. Nonexponential removal is hypothesized to arise from variable likelihood of encountering nanoscale regions of attraction (heterodomains) on grain surfaces that allow attachment. Implementing heterodomains in mechanistic simulations of pore scale trajectories generates continuum-scale rate coefficients that produce experimentally-observed breakthrough-elution curves (BTEC) and retention profiles (RP). However, current one-directional simulation across scales is inefficient in finding a heterodomain surface coverage that yields observed RP and BTEC. In this work, we develop an inverse two-population model approach that not only estimates colloid rate coefficients from experimental BTEC and RP data but also quantifies the associated uncertainty, thereby allowing the problem to be worked from both ends via comparison of (1) rate coefficients upscaled from mechanistic pore scale simulations incorporating heterodomains with (2) rate coefficients inverted from continuum-scale BTEC and RP. We validate our inverse model using synthetic data with known removal rates and subsequently demonstrate its applicability to experimental data with multiexponential and nonmonotonic RPs. We moreover derive what is likely the first analytical expression for the RP under repulsive conditions, revealing that the hypoexponential distribution can be used to reproduce multiexponential and non-monotonic shapes. By addressing key limitations in present models, our inverse approach offers a valuable tool for advancing colloidal transport predictions in natural environments.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104905"},"PeriodicalIF":4.0,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143342377","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":"Simplifications of macroscopic models for heat and mass transfer in porous media","authors":"Didier Lasseux ,&nbsp;Francisco J. Valdés-Parada","doi":"10.1016/j.advwatres.2025.104899","DOIUrl":"10.1016/j.advwatres.2025.104899","url":null,"abstract":"<div><div>When performing upscaling of transport phenomena in multiscale systems it is not uncommon that terms of different physical nature than those present at the underlying scale arise in the resulting averaged differential equations. For diffusive species mass transfer with heterogeneous reaction and conductive heat transfer, additional terms result from upscaling using the volume averaging method, which are classically discarded by means of orders of magnitude estimates. In this work, these two cases are revisited and it is shown that, for single and two-species diffusive mass transfer with heterogeneous nonlinear reaction, the additional term is exactly zero using Green’s formula. This conclusion is shown to also be applicable when using the periodic homogenization method. Nevertheless, for heat conduction, with and without considering interfacial resistance, only the dominant conduction-corrective terms are shown to be zero also using Green’s formula. In contrast, the contribution of the co-conduction-corrective terms may be relevant depending on the systems characteristics, the properties of the phases and the macroscopic boundary conditions. This is exemplified by performing numerical simulations in a non-symmetric unit cell.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"197 ","pages":"Article 104899"},"PeriodicalIF":4.0,"publicationDate":"2025-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143342378","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":"On the flow behaviour of unconfined dual porosity aquifers with sloping base","authors":"Konstantinos N. Moutsopoulos ,&nbsp;John N.E. Papaspyros ,&nbsp;Antonis D. Koussis ,&nbsp;Frederick Delay ,&nbsp;Marwan Fahs","doi":"10.1016/j.advwatres.2024.104856","DOIUrl":"10.1016/j.advwatres.2024.104856","url":null,"abstract":"<div><div>The flow in an unconfined double-porosity aquifer with a sloping base is investigated and equations are developed for its description. In order to obtain a relatively simple description of the problem, similar assumptions as in Moutsopoulos (2021) have been adopted; the pressure in the unsaturated zone, especially in the fractures’ network, is considered to be atmospheric, and the Dupuit-Forchheimer approximation is invoked, reducing the dimensionality by eliminating the vertical direction. The derived equations have been linearized and solved analytically for a problem involving interaction between an inclined aquifer and an adjacent surface water body similar to the one examined by Akylas and Koussis (2007). The analytical solution has been checked against results obtained with state-of-the-art numerical codes. The agreement between the two approaches is excellent. The solution tools were used to gain insight in the influence of the aquifer's base inclination on the flow quantities.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"195 ","pages":"Article 104856"},"PeriodicalIF":4.0,"publicationDate":"2025-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143167151","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":"Numerical investigation of multiphase flow through self-affine rough fractures","authors":"Prakash Purswani,&nbsp;Javier E. Santos,&nbsp;Jeffrey D. Hyman,&nbsp;Eric J. Guiltinan","doi":"10.1016/j.advwatres.2024.104852","DOIUrl":"10.1016/j.advwatres.2024.104852","url":null,"abstract":"<div><div>Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"195 ","pages":"Article 104852"},"PeriodicalIF":4.0,"publicationDate":"2025-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143167149","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":"On transient qanat discharge in an unconfined aquifer underlain by a fractured aquifer","authors":"Mohammad M. Sedghi ,&nbsp;Hongbin Zhan","doi":"10.1016/j.advwatres.2024.104874","DOIUrl":"10.1016/j.advwatres.2024.104874","url":null,"abstract":"<div><div>Fractured aquifers sometimes underlie unconfined aquifers to form unconfined-fractured aquifer systems. The discharges of qanats constructed in such aquifer systems depend not only on the hydraulic parameters of the overlying aquifer but also on the hydraulic parameters of the underlying aquifer, where qanats can be approximated as constant-head slightly inclined horizontal wells. The available analytical solutions of the transient discharge variation of a qanat usually consider a single-layer unconfined aquifer and cannot be used to simulate the transient discharge variation of a qanat in an unconfined-fractured aquifer system. Thus, the goal of this work is to develop such an analytical model. The solution of the drawdown due to a point sink and the solution of the groundwater mound due to an areal recharge in an unconfined-fractured aquifer system are obtained using Laplace and Fourier transforms. These solutions are utilized to obtain the discharge variation of qanat in an unconfined-fractured aquifer system. By removing the inter-porosity flow term from the underlying fractured aquifer, the solution can be used to simulate the discharge variation of qanat in a two-layer aquifer or an aquifer-aquitard system. The effects of the hydraulic and geometric parameters of the underlying fractured aquifer on the discharge variation of qanat in an overlying, unconfined aquifer are explored. The solution can be utilized to predict the transient discharge of qanat due to rainfall or constant head horizontal well in an unconfined-fractured aquifer system; to estimate the hydraulic parameters of the underlying fractured aquifer using the discharge data of qanat or a horizontal constant-head well; to evaluate the response of qanat to areal recharge due to rainfall among other applications, thereby demonstrating the practical relevance and potential impact of our research on real-world hydrological problems.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"195 ","pages":"Article 104874"},"PeriodicalIF":4.0,"publicationDate":"2025-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142825420","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":"Fully coupled morphological modelling under the combined action of waves and currents","authors":"Wei Li ,&nbsp;Yiming Zhang ,&nbsp;Peng Hu","doi":"10.1016/j.advwatres.2024.104875","DOIUrl":"10.1016/j.advwatres.2024.104875","url":null,"abstract":"<div><div>The non-linear interactions of tides, waves and sediment transport in shallow coastal waters require coupled modelling of waves and currents, as well as the consideration of feedbacks between flow, sediment transport and morphology. However, most of the previous work about the wave-current morphological modelling lacks full consideration of the feedbacks of bed deformation on flow or/and assumes instant sediment adaption, which may be subjected to uncertainties. To bridge this gap, a two-dimensional depth-averaged (2DH) coastal model that couples an advanced fully coupled flow-morphological model and the phase-averaged SWAN model by a two-way coupling method has been developed. The feedback impacts of bed deformation on the flow, and the spatial and temporal scales required for sediment adaptation to the capacity regime are fully considered in the fully coupled flow-morphological model. For efficient and accurate modelling, a hybrid local time step/global maximum time step (LTS/GMaTS) method and the Open Multi-Processing (OpenMP) technique are used together with the shock-capturing finite volume method and the Harten-Latex-van Leer-Contact (HLLC) approximate Riemann solver. The applicability and accuracy of the present coastal model has been tested by four cases including the wave-induced currents at a tri-cuspate beach, longshore currents and sediment transport over transverse bars, channel infilling and migration under wave-current actions, and surge-wave driven morphology in a bay-inlet system. Reasonable agreement between the computation and measurement can be obtained if empirical and numerical parameters of the model are sensibly specified. The fully coupled morphological approach is suggested especially for extreme wave-current conditions due to its completeness in physics and negligible extra cost in computation. In addition, the hybrid LTS/GMaTS method can also effectively accelerate the computation in the wave-current environment.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"195 ","pages":"Article 104875"},"PeriodicalIF":4.0,"publicationDate":"2025-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143166700","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":"Visualization experiment for characterizing the influence of an open fracture on CO2 behavior in fractured porous media","authors":"Seung-Wook Ha ,&nbsp;Ji-Young Baek ,&nbsp;In-Woo Park ,&nbsp;Kang-Kun Lee","doi":"10.1016/j.advwatres.2024.104871","DOIUrl":"10.1016/j.advwatres.2024.104871","url":null,"abstract":"<div><div>While gas behavior influenced by the intricate interplay between fractures and the surrounding porous media can create complex flow paths and distribution, there is a scarcity of experimental research focusing on the influence of a fracture on gas behavior due to the limitation of implementing an open structure within porous media. In this study, we have implemented transparent porous structures embedding an open fracture in porous media by using a 3D printer. The study encompasses a series of experiments to visualize and quantify the spatiotemporal distribution of gaseous and dissolved CO₂ concentration by high-resolution pixel-wise image analysis. Results show that the presence of an open fracture significantly influences the behavior of CO₂, affecting both the migration and distribution of gas and dissolved CO₂ concentration. In addition, capillary fingering can be the trigger of early breaching across the fracture into the upper matrix which can create a highly discontinuous and heterogeneous spatiotemporal CO₂ distribution. The implications of this research are not limited to CO₂ behavior but extend to other gases and immiscible contaminants such as non-aqueous phase liquid whose behavior is different from that of groundwater. Our results suggest the applicability of the 3D printing technique in visualization experiments and the importance of incorporating detailed geological structures into environmental monitoring and predictive models for an effective environmental management.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"195 ","pages":"Article 104871"},"PeriodicalIF":4.0,"publicationDate":"2025-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143167155","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}