Advances in Water Resources最新文献

{"title":"A generalized analytical model for tidal wave propagation in multilayered, horizontally heterogeneous leaky aquifer systems","authors":"Chao Zhuang ,&nbsp;Desheng Bao ,&nbsp;Long Yan ,&nbsp;Hongbin Zhan ,&nbsp;Yun Yang ,&nbsp;Qiaona Guo ,&nbsp;Walter A. Illman ,&nbsp;Zhifang Zhou ,&nbsp;Jinguo Wang","doi":"10.1016/j.advwatres.2025.105156","DOIUrl":"10.1016/j.advwatres.2025.105156","url":null,"abstract":"<div><div>Existing analytical models for tidal wave propagation in coastal leaky aquifers are limited to simplified homogeneous or layered systems, failing to address complex heterogeneity such as aquitard preferential conduits. This study develops a generalized analytical model accommodating an arbitrary number of vertically stratified aquifer and aquitard layers and horizontally segmented zones with distinct hydraulic properties. The model’s generality is enhanced by incorporating four key features: non-negligible aquitard storativity, finite lateral extent, zone-specific tidal forcing at the top boundary, and zone-specific tidal loading effects on both aquifer and aquitard layers. The solution is derived by employing matrix eigenvalue analysis to obtain the zone-specific general solution, followed by a recursive formulation to enforce hydraulic head and flux continuity across aquifer zone interfaces. The newly derived analytical solution subsumes multiple existing analytical solutions as special cases. Theoretical investigations reveal that aquitard preferential conduits can reverse the typical landward amplitude decay and phase shift increase in the less permeable adjacent aquifer. This reversed tidal wave propagation pattern is pronounced under strong aquifer permeability contrasts, high conduit-to-aquitard permeability ratios, and large conduit widths, and it exhibits remarkable period dependence. For conduit identification through multi-frequency tidal analysis, greater weight should be given to higher-frequency components.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105156"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145412653","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":"Managing groundwater extraction to desalinate residual saltwater in coastal subsurface reservoirs: Transient behaviors and environmental implications","authors":"Yunhai Fang ,&nbsp;Lei Ma ,&nbsp;Tianyuan Zheng ,&nbsp;Jian Luo ,&nbsp;Huan Wang ,&nbsp;Jiazhong Qian ,&nbsp;Xilai Zheng ,&nbsp;Marc Walther","doi":"10.1016/j.advwatres.2025.105175","DOIUrl":"10.1016/j.advwatres.2025.105175","url":null,"abstract":"<div><div>The construction of subsurface reservoirs using physical barriers for controlling seawater intrusion triggers the retention of residual saltwater, raising concerns about contamination in regions intended for freshwater supply. This study investigated the temporal pollution behaviors of residual saltwater in coastal subsurface reservoirs. Groundwater extraction was analyzed to compare the effects of constant and intermittent pumping on the spread of residual saltwater. The results showed that residual saltwater pollution was a great threat to freshwater resources, contaminating up to 20 % of the reservoir area after 20,000 d under natural (no-pumping) conditions. Constant pumping exhibited positive environmental effects, inhibiting the inland intrusion of residual saltwater, reducing contamination, and accelerating desalination with high efficiency. Its contribution was more evident over time; for instance, the contaminated area was reduced to 4 % at 20,000 d with a desalination efficiency improvement of 28 %. The desalination mechanism under constant pumping was different from natural conditions, with minimal impact on high-salinity zones and its primary effect observed in reducing contamination in low-salinity zones. In contrast, intermittent pumping had adverse effects, increasing the contaminated area and delaying desalination, requiring caution in its application. The desalination efficiency and failure risk of pumping were related to well location and pumping rate. By optimizing well setup, an effective remediation strategy of residual saltwater was proposed, achieving at least 20 % improvement in desalination efficiency. These findings provide a general framework for understanding and addressing residual saltwater pollution and offer practical guidance for the management of coastal subsurface reservoirs.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105175"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145531170","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 tracer dispersion in a coupled system composed of a proppant-packed hydraulic fracture and a tight porous medium","authors":"Morteza Dejam","doi":"10.1016/j.advwatres.2025.105122","DOIUrl":"10.1016/j.advwatres.2025.105122","url":null,"abstract":"<div><div>Tracer tests are widely performed for the characterization of reservoir properties during the hydraulic fracturing operation. The dispersion of the tracer depends on the interaction of the proppant-packed hydraulic fracture and the tight porous medium through the naturally porous walls. However, the effects of the interaction of the porous walls and dynamics of flow in the proppant-packed hydraulic fracture on the tracer dispersion and reservoir dynamic mass/heat storage capacity have not yet been reported in the literature. In this work, the tracer dispersion in a proppant-packed hydraulic fracture surrounded by a tight porous medium is theoretically modeled and the dynamic storage capacity is evaluated. The Darcy-Brinkman equation is used to describe the fully developed laminar Stokes fluid flow in the proppant-packed hydraulic fracture. We used the Taylor dispersion theory and Reynolds decomposition approach to derive the exact equivalent transport parameters, including dispersion and advection coefficients, as well as the storage capacity of the tight porous medium. It is found that the tracer dispersion is controlled by the Darcy and the Peclet numbers in the proppant-packed hydraulic fracture. The results indicate that the ratio of tracer dispersion in the proppant-packed hydraulic fracture with porous walls to that with nonporous walls ranges from zero for very small Darcy numbers to 0.3 for large Darcy numbers. The ratio of the advection velocity in the proppant-packed hydraulic fracture with porous walls to that with nonporous walls ranges from unity for very small Darcy numbers to 7/5 for large Darcy numbers. The results also indicate that tracer mass storage capacity in the tight porous medium increases as the Peclet number for fluid flow in the proppant-packed hydraulic fracture increases. Conversely, storage decreases as the Darcy number in the proppant-packed hydraulic fracture rises. A comparison reveals that a flow transport model based on proppant-free hydraulic fracture may lead to the overestimation of the tracer mass/heat storage capacity. The findings of this study pave the way to advance our understanding of tracer tests for evaluating reservoir characteristics during fracturing operations in enhanced geothermal systems.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105122"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145094011","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":"Estimating the full anisotropy of the covariance function in geostatistical inversion using the pilot-point ensemble Kalman filter","authors":"Janek Geiger,&nbsp;Michael Finkel,&nbsp;Olaf A. Cirpka","doi":"10.1016/j.advwatres.2025.105103","DOIUrl":"10.1016/j.advwatres.2025.105103","url":null,"abstract":"<div><div>In geostatistical inversion, good prior knowledge about the covariance function is important in estimating hydraulic conductivity from hydraulic-head observations, but may be hampered by poor knowledge about anisotropy. In this study we propose an extension of the pilot-point ensemble Kalman filter (PP-EnKF) that can infer the full anisotropy of the covariance function based on attainable, initially random knowledge. We address the periodicity of rotation by incorporating the unique elements of the covariance transformation matrix into the set of parameters to be estimated. The filter is further modified by generating conditional realizations in each assimilation step, increasing the inherent variance of the ensemble and counteracting filter inbreeding. We demonstrate the methodology in a synthetic study of a 2-D groundwater-flow model where we estimate the full anisotropy of the covariance function and the hydraulic conductivity at pilot points via the assimilation of hydraulic-head data. The success of this method depends more on the configuration of pilot points than on the quality of prior knowledge, as ensembles initialized with faulty random priors successfully estimated the correct parameters of the covariance function, as well as the log-hydraulic conductivity values at the pilot points. The resulting parameter fields enabled accurate predictions of hydraulic heads during a verification period, with normalized root mean square errors reduced by up to 66% compared to ensembles with isotropic covariance functions. The methodology presented in this study mitigates the importance of informative prior knowledge of the covariance function in geostatistical parameter-inference methods, especially in highly anisotropic settings.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105103"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145263652","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":"Measurement as functional minimization: An image-based method for river surface velocity and discharge estimation","authors":"Kailin Huang ,&nbsp;Lionel Pénard ,&nbsp;Hua Chen ,&nbsp;Chong-Yu Xu","doi":"10.1016/j.advwatres.2025.105160","DOIUrl":"10.1016/j.advwatres.2025.105160","url":null,"abstract":"<div><div>Hydrological monitoring establishes the foundation for scientific studies of water resources. The introduction of image-based hydrometry methods offers a non-intrusive and highly efficient approach for collecting discharge data from rivers. However, challenges persist for image-based hydrometry methods, such as partial observation, light reflection, unknown surface-to-depth coefficients, and unphysical velocity estimates requiring post-processing. In this study we propose a new image-based functional optimization method to address these issues. This Lagrangian multipliers formulation transforms river surface velocity and discharge estimation into the minimization of a distance in the functional space. The algorithm effectively handles partial, large scale, and high reflection observations, and the velocity field obtained through this method adheres to the physical governing equations. Three experiments were conducted across rivers of different scales, lighting, and weather conditions in France and China to demonstrate the applicability and robustness of the proposed algorithm. Through optimization of an objective functional, the results reveal that the proposed method can directly estimate river surface velocity and discharge from an image sequence using numerical model without the need for post-processing outliers or missing values. The study demonstrates that the proposed method has the potential to be widely deployed in different hydrological monitoring scenarios to improve hydrologic monitoring efficiency.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105160"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145383870","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":"Competitive dissolution of binary minerals in porous media: A lattice Boltzmann study","authors":"Wenxin Yang,&nbsp;Hai Sun,&nbsp;Lei Zhang,&nbsp;Gloire Imani,&nbsp;Dongyan Fan,&nbsp;Junjie Zhong,&nbsp;Yongfei Yang,&nbsp;Jun Yao","doi":"10.1016/j.advwatres.2025.105132","DOIUrl":"10.1016/j.advwatres.2025.105132","url":null,"abstract":"<div><div>In the geological storage of hydrogen and carbon dioxide in underground salt caverns, the water-injected salt dissolution technology is widely applied in the expansion of salt caverns. During the salt dissolution, a variety of minerals with different properties are often involved, and there are obvious competitive interactions among them. This competition is primarily manifested in differences in diffusion rates, dissolution rates, and ionic concentration equilibrium constraints. In this study, we developed a lattice Boltzmann model, taking into account the competitive dissolution mechanisms of different minerals under the constraint of ion equilibrium. The model was used to investigate permeability changes of porous structure under various injection velocities and different combinations of minerals. The results reveal that whether the physical properties of binary minerals vary greatly or not, such as reaction rate and diffusion rate, the impact of competitive dissolution cannot be ignored. And as the injection rate increases, the influence of the competitive dissolution effect on the pore structure evolution becomes greater. This research provides theoretical insights into binary minerals' competitive dissolution mechanisms and references for its applications in fields such as environmental science, resource development, and chemical engineering.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105132"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145217975","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":"Statistical integro-differential fracture model (Sid-FM) for isolated fractures with variable apertures and lengths","authors":"Daniel Stalder,&nbsp;Shangyi Cao,&nbsp;Daniel W. Meyer,&nbsp;Patrick Jenny","doi":"10.1016/j.advwatres.2025.105125","DOIUrl":"10.1016/j.advwatres.2025.105125","url":null,"abstract":"<div><div>Flow in fractured porous media is associated with high uncertainty, particularly regarding fracture properties and their overall configuration within the domain. This is especially pronounced for disconnected fractures of smaller yet comparable size to the domain. Consequently, ensemble averages are often used to capture this statistical variability and predict the expected behavior. This leads to enormous computational costs, as flow simulations of single realizations with millions of fractures are extremely expensive; and much more so full Monte Carlo studies involving hundreds of realizations. Alternatively, a recently introduced model aims to directly estimate expected flow rates and pressure fields. The model involves few degrees of freedom, leading to low-cost computations. This is achieved by using integro-differential equations involving non-local kernel functions that encompass the statistical information of fractures. So far this statistical integro-differential fracture model (Sid-FM) considers only ensembles with identical fractures having constant aperture and lengths. In this paper Sid-FM is extended to account for arbitrary fracture aperture profiles and reservoirs with fractures following specified length distributions, which is a crucial step towards applications with realistic fractured reservoirs. In a series of numerical experiments, it is demonstrated that the Sid-FM’s predictions are in excellent agreement with Monte Carlo reference data, which are based on many fracture-resolving simulations. The applicability is demonstrated through statistically one-dimensional cases, laying crucial groundwork for 2D and 3D extensions. Future work will focus on further generalizations and extensions such as transport processes and 2D/3D applications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105125"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145217979","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 closed-form solution for pre-asymptotic dispersion with transient velocity","authors":"Sabrina N. Volponi,&nbsp;Rajitha Senavirathna,&nbsp;Timothy R. Ginn","doi":"10.1016/j.advwatres.2025.105172","DOIUrl":"10.1016/j.advwatres.2025.105172","url":null,"abstract":"<div><div>Advective–dispersive transport in natural systems is often modeled assuming steady one-dimensional flow and complete cross-sectional mixing, the latter allowing the use of a constant dispersion coefficient. In reality, however, many environments exhibit pre-asymptotic dispersion (PAD), where the effective dispersion coefficient increases as solutes progressively sample velocity heterogeneity. Time- and/or distance-dependent dispersion coefficients have historically been used to address PAD but can produce inconsistencies when solute pulses interact, particularly in nonlinear systems where superposition fails. Here we generalize our recent closed-form solution for one-dimensional advective–dispersive transport with PAD to the case of transient velocity. This new solution involves expressing dispersion as a function of both solute age in the flow (residence time), and physical time. These solutions apply to both initial and boundary value problems over an infinite domain, including cases with first-order reactions. We validate our boundary value solution by reproducing breakthrough curves from turbulent, accelerating pipe flow experiments. Our results show that a solute’s “memory” of reduced early-time dispersion continues to influence its transport even after complete mixing is achieved. More broadly, this framework fills a critical modeling gap by enabling physically consistent analysis of solute transport under PAD and transient flow—conditions common in natural systems but rarely studied together.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105172"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145485561","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":"Migration mechanisms of oil-water-sand three-phase flow in unconsolidated sandstones during water flooding: An image analysis approach resilient to threshold interference","authors":"Maozhu Li,&nbsp;Yiqiang Li,&nbsp;Zhipeng Wang,&nbsp;Xinji Du,&nbsp;Lei Wang,&nbsp;Jinxin Cao,&nbsp;Zheyu Liu","doi":"10.1016/j.advwatres.2025.105154","DOIUrl":"10.1016/j.advwatres.2025.105154","url":null,"abstract":"<div><div>Unconsolidated sandstone reservoirs exhibit high sand production susceptibility during water flooding, impairing fluid migration and well productivity. Elucidating three-phase (oil-water-sand) migration mechanisms and residual fluid distribution characteristics is critical for sanding rate prediction and fluid saturation estimation. This study examines sand migration and residual oil distribution during waterflooding in unconsolidated sandstone oil reservoirs. Using methylene blue-dyed water to flood artificial models, we replicated sand production patterns and quantified outcomes via novel image analysis. Results show that: In water flooding processes through heterogeneous unconsolidated sandstones, sand production rate exhibits a quadratic polynomial relationship with time. Greater fluctuation ranges in areal sweep variation rates indicate frequent flooding front diversions. Post-flood residual oil saturation measures 66.07–71.87 % in heterogeneous models versus 70.99–76.65 % in homogeneous counterparts. Residual oil saturation demonstrates broader distribution in swept zones of heterogeneous formations, whereas homogeneous systems show localized concentration. Strategic injection placement in high-permeability zones induces enhanced sweep efficiency when fluid propagates into low-permeability regions, establishing optimized drainage networks. These insights aid sand control and residual fluid management during fluid injection/production. The image analysis method enables pore-scale quantitative characterization for aquifer management and contaminant remediation.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105154"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145358153","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":"Experimental and numerical studies on the transport mechanism of clay particles in porous media considering double layer forces","authors":"Xiaoyi Zhao ,&nbsp;Xiangjun Pei ,&nbsp;Xiaochao Zhang ,&nbsp;Mao Qiu ,&nbsp;Renji Zhang","doi":"10.1016/j.advwatres.2025.105145","DOIUrl":"10.1016/j.advwatres.2025.105145","url":null,"abstract":"<div><div>Micron- and nanoscale particles can cause varying degrees of pore clogging in porous media due to inter-particle forces, yet their role in seepage erosion remains unclear. This study couples computational fluid dynamics (CFD) and the discrete element method (DEM), incorporating Coulomb’s law and the Johnson–Kendall–Roberts model to simulate electrostatic repulsion and van der Waals attraction, respectively, while accounting for surface roughness and force distance. Through simulations of micron-scale clay particle seepage at varying flow rates, we systematically analyzed particle erosion, hydraulic conductivity, clogging types, and micromechanical fields. A custom-designed microfluidic chip, matching the numerical model, validated the simulations. Results show that inter-particle force and fluid resistance nearly doubled, and the clogging rate increased by ∼10 % when non-contact forces and surface roughness were considered. However, their influence diminished significantly at flow velocities above 1.6 mm/s. Dependent clogging dominated, particularly when strong force chains formed mid-pore, resisting breakup. The Kozeny–Carman equation explained 91.336 % of hydraulic conductivity variation. The findings contribute to the understanding and prediction of erosion processes in porous media, providing a theoretical basis for related engineering.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"206 ","pages":"Article 105145"},"PeriodicalIF":4.2,"publicationDate":"2025-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145324351","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}