Advances in Water Resources最新文献_第3页

Pore-scale imaging of H₂, N₂, and CO₂ behavior in sandstone: Wettability, trapping, pore occupancy, and displacement patterns 砂岩中H₂、N₂和CO₂行为的孔隙尺度成像：润湿性、圈闭、孔隙占用和位移模式

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-02-01 Epub Date: 2026-01-09 DOI: 10.1016/j.advwatres.2026.105213

Guanglei Zhang , Jinliang Chen , Yaxin Shao , Martin J. Blunt , Guowei Ma

Underground hydrogen storage (UHS) in porous aquifers is a promising solution for large-scale renewable energy storage. Despite growing interest in underground hydrogen storage, a fundamental gap remains in understanding whether hydrogen behaves distinctly from other subsurface gases at the pore scale. This study employs in situ micro-CT imaging to directly compare gas injection and withdrawal of H₂, N₂, and CO₂ in Bentheimer sandstone. The results confirm a water-wet system for all gases, with contact angles between 50° and 70°. However, significant differences in gas trapping and distribution emerged. After gas injection, the initial gas saturations were 33 % for H₂, 37 % for N₂, and 36 % for CO₂. After imbibition (withdrawal), H₂ exhibited the lowest residual saturation (14 %), compared to N₂ (25 %), while CO₂ was completely dissolved. Low residual saturation of H₂ indicate high recovery efficiency, which is a favorable for hydrogen storage. Analysis of gas ganglia revealed that H₂'s low viscosity and density promote capillary fingering and heterogeneous redistribution. In contrast, N₂'s higher viscosity facilitates stable displacement and higher residual saturation, making it a suitable candidate as a cushion gas. The capillary pressure for H₂ was found to be approximately 1.2 times higher than for N₂, a difference not fully explained by interfacial tension alone. These findings demonstrate that gas properties, particularly solubility and viscosity, govern pore-scale distribution and trapping mechanisms. Consequently, N₂ and CO₂ are poor proxies for predicting H₂ behavior, and accurate forecasting of UHS performance requires models that account for H₂'s mobility and its propensity for unstable flow.

多孔含水层地下储氢（UHS）是一种很有前途的大规模可再生能源储存方案。尽管人们对地下储氢的兴趣日益浓厚，但在了解氢气在孔隙尺度上是否与其他地下气体表现不同方面，仍然存在一个根本性的差距。本研究采用原位微ct成像技术，直接对比Bentheimer砂岩中H₂、N₂、CO₂的注、抽气情况。结果证实了水湿系统适用于所有气体，接触角在50°到70°之间。然而，在天然气捕获和分布方面出现了显著差异。注气后，初始气体饱和度H₂为33%，N₂为37%，CO₂为36%。在渗吸（抽吸）后，H₂的残余饱和度最低（14%），而N₂的残余饱和度最低（25%），而CO₂完全溶解。H 2残余饱和度低，表明回收效率高，有利于储氢。气体神经节的分析表明，H₂的低粘度和密度促进了毛细管指动和非均质再分布。相比之下，n2的高粘度有利于稳定的驱替和较高的残余饱和度，使其成为缓冲气体的合适人选。研究发现，H₂的毛细管压力大约是N₂的1.2倍，这种差异不能完全用界面张力来解释。这些发现表明，气体性质，特别是溶解度和粘度，控制着孔隙尺度分布和捕集机制。因此，N₂和CO₂是预测H₂行为的不良代理，而对UHS性能的准确预测需要考虑H₂的流动性及其不稳定流动倾向的模型。

{"title":"Pore-scale imaging of H₂, N₂, and CO₂ behavior in sandstone: Wettability, trapping, pore occupancy, and displacement patterns","authors":"Guanglei Zhang , Jinliang Chen , Yaxin Shao , Martin J. Blunt , Guowei Ma","doi":"10.1016/j.advwatres.2026.105213","DOIUrl":"10.1016/j.advwatres.2026.105213","url":null,"abstract":"<div><div>Underground hydrogen storage (UHS) in porous aquifers is a promising solution for large-scale renewable energy storage. Despite growing interest in underground hydrogen storage, a fundamental gap remains in understanding whether hydrogen behaves distinctly from other subsurface gases at the pore scale. This study employs <em>in situ</em> micro-CT imaging to directly compare gas injection and withdrawal of H₂, N₂, and CO₂ in Bentheimer sandstone. The results confirm a water-wet system for all gases, with contact angles between 50° and 70°. However, significant differences in gas trapping and distribution emerged. After gas injection, the initial gas saturations were 33 % for H₂, 37 % for N₂, and 36 % for CO₂. After imbibition (withdrawal), H₂ exhibited the lowest residual saturation (14 %), compared to N₂ (25 %), while CO₂ was completely dissolved. Low residual saturation of H₂ indicate high recovery efficiency, which is a favorable for hydrogen storage. Analysis of gas ganglia revealed that H₂'s low viscosity and density promote capillary fingering and heterogeneous redistribution. In contrast, N₂'s higher viscosity facilitates stable displacement and higher residual saturation, making it a suitable candidate as a cushion gas. The capillary pressure for H₂ was found to be approximately 1.2 times higher than for N₂, a difference not fully explained by interfacial tension alone. These findings demonstrate that gas properties, particularly solubility and viscosity, govern pore-scale distribution and trapping mechanisms. Consequently, N₂ and CO₂ are poor proxies for predicting H₂ behavior, and accurate forecasting of UHS performance requires models that account for H₂'s mobility and its propensity for unstable flow.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"208 ","pages":"Article 105213"},"PeriodicalIF":4.2,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956968","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}

引用次数: 0

Pore-scale dynamics of multiphase reactive transport in water-wet carbonates under co2-acidified brine injection: Dissolution patterns and reaction rates co2酸化盐水注入下水湿碳酸盐多相反应输运的孔隙尺度动力学：溶解模式和反应速率

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-02-01 Epub Date: 2026-01-05 DOI: 10.1016/j.advwatres.2025.105202

Qianqian Ma , Rukuan Chai , Sajjad Foroughi , Yanghua Wang , Martin J. Blunt , Branko Bijeljic

Depleted carbonate reservoirs are promising sites for geological CO

_{2}

storage, yet the presence of residual hydrocarbon introduces complex pore-scale interactions that influence the dynamics of solid dissolution. We combined time-resolved X-ray microtomography (micro-CT), core-flooding experiments, and pore-scale modeling to investigate how residual hydrocarbon affects dissolution patterns and effective reaction rates during CO

_{2}

-acidified brine injection into Ketton limestone under reservoir conditions. We find that the pore structure and fluid distribution control flow heterogeneity, reactive surface accessibility, dissolution patterns and the reaction rates. At low injection rate, two distinct dissolution patterns were observed: 1) a positive feedback loop of channel widening that efficiently enhanced transport properties; and 2) a suppressed regime in which heterogeneity and hydrocarbon blockage resulted in only a modest increase in permeability. At high injection rates, a more uniform dissolution occurred caused by re-mobilization of hydrocarbon that initially blocked the flow of brine. Effective reaction rates in two-phase flow were lower than in the equivalent single-phase case and up to two orders of magnitude lower than the batch rates due to persistent transport limitations. These findings provide mechanistic insights into multiphase reactive transport in carbonates and highlight the importance of accurately understanding the impact of the residual phase on reactions to improve predictions of CO

_{2}

storage efficiency.

枯竭的碳酸盐岩储层是地质CO2储存的理想场所，但残余烃的存在引入了复杂的孔隙尺度相互作用，影响了固体溶解的动力学。我们结合了时间分辨x射线微断层扫描（micro-CT）、岩心驱油实验和孔隙尺度模型，研究了油藏条件下二氧化碳酸化盐水注入Ketton石灰岩过程中残余烃对溶解模式和有效反应速率的影响。我们发现孔隙结构和流体分布控制着流动非均质性、反应表面可达性、溶解模式和反应速率。在低注入速率下，观察到两种不同的溶解模式：1)通道加宽的正反馈回路有效地增强了输运性质；2)非均质性和油气堵塞只导致渗透率适度增加的抑制状态。在高注入速率下，由于碳氢化合物的再动员，最初阻碍了盐水的流动，导致了更均匀的溶解。由于持续的输运限制，两相流的有效反应速率低于等效的单相流，并且比批处理速率低两个数量级。这些发现为碳酸盐中多相反应输运提供了机理见解，并强调了准确理解剩余相对反应的影响对提高二氧化碳储存效率预测的重要性。

{"title":"Pore-scale dynamics of multiphase reactive transport in water-wet carbonates under co2-acidified brine injection: Dissolution patterns and reaction rates","authors":"Qianqian Ma , Rukuan Chai , Sajjad Foroughi , Yanghua Wang , Martin J. Blunt , Branko Bijeljic","doi":"10.1016/j.advwatres.2025.105202","DOIUrl":"10.1016/j.advwatres.2025.105202","url":null,"abstract":"<div><div>Depleted carbonate reservoirs are promising sites for geological CO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span> storage, yet the presence of residual hydrocarbon introduces complex pore-scale interactions that influence the dynamics of solid dissolution. We combined time-resolved X-ray microtomography (micro-CT), core-flooding experiments, and pore-scale modeling to investigate how residual hydrocarbon affects dissolution patterns and effective reaction rates during CO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span>-acidified brine injection into Ketton limestone under reservoir conditions. We find that the pore structure and fluid distribution control flow heterogeneity, reactive surface accessibility, dissolution patterns and the reaction rates. At low injection rate, two distinct dissolution patterns were observed: 1) a positive feedback loop of channel widening that efficiently enhanced transport properties; and 2) a suppressed regime in which heterogeneity and hydrocarbon blockage resulted in only a modest increase in permeability. At high injection rates, a more uniform dissolution occurred caused by re-mobilization of hydrocarbon that initially blocked the flow of brine. Effective reaction rates in two-phase flow were lower than in the equivalent single-phase case and up to two orders of magnitude lower than the batch rates due to persistent transport limitations. These findings provide mechanistic insights into multiphase reactive transport in carbonates and highlight the importance of accurately understanding the impact of the residual phase on reactions to improve predictions of CO<span><math><msub><mrow></mrow><mrow><mn>2</mn></mrow></msub></math></span> storage efficiency.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"208 ","pages":"Article 105202"},"PeriodicalIF":4.2,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145902871","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}

引用次数: 0

A modified capillary pressure–saturation model for unsaturated flow in asphalt mixtures and hydraulic behavior analysis 沥青混合料非饱和流动的改进毛细管压力-饱和模型及水力特性分析

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-02-01 Epub Date: 2025-12-29 DOI: 10.1016/j.advwatres.2025.105203

Xinxing Bian , Huining Xu , Shiyuan Li

The capillary pressure–saturation (P_c–S) relationship is a critical constitutive model for modeling multiphase flow in unsaturated porous media. In the context of asphalt pavements, however, existing P_c–S models predominantly rely on direct adoption of empirical formulations in unsaturated soil, lack rigorous experimental validation, and often fail to capture the unsaturated hydraulic characteristics of pavement materials. This study aims to (1) develop a continuous mathematical framework for the P_c–S relationship in asphalt mixtures by integrating the physical constraints derived from experimental data with the functional form of classical models, and (2) elucidate the impact of this relationship on unsaturated seepage dynamics. A modified van Genuchten model was formulated to provide a concise, physically consistent P_c–S function. The proposed model demonstrates improved fidelity to experimental data compared to conventional soil-based hydraulic models, while preserving mathematical tractability. By implementing this P_c–S constitutive relationship within a numerical solver (SEEP/W) for the Richards equation, the hysteretic behavior of unsaturated flow was simulated in asphalt mixtures. The simulations reveal that the unsaturated flow analysis, incorporating the proposed P_c–S relationship, predicts a more complex and realistic moisture transport regime than the saturated flow approximation, and can reflect the observed in-service behavior of pavements. This work advances understanding of moisture transport in unsaturated asphalt mixtures and provides a useful tool for evaluating the long-term hydraulic performance in pavement–subsurface water systems.

毛细管压力-饱和度（Pc-S）关系是模拟非饱和多孔介质中多相流动的重要本构模型。然而，对于沥青路面，现有的Pc-S模型主要依赖于直接采用非饱和土的经验公式，缺乏严格的实验验证，往往无法捕捉路面材料的非饱和水力特性。本研究旨在(1)通过整合实验数据得出的物理约束与经典模型的函数形式，建立沥青混合料Pc-S关系的连续数学框架；(2)阐明该关系对非饱和渗流动力学的影响。一个改进的van Genuchten模型被制定，以提供一个简洁的，物理一致的Pc-S功能。与传统的基于土壤的水力模型相比，该模型在保持数学可追溯性的同时，提高了对实验数据的保真度。通过在Richards方程的数值求解器（SEEP/W）中实现Pc-S本构关系，模拟了沥青混合料中非饱和流动的滞后行为。模拟结果表明，结合Pc-S关系的非饱和流动分析比饱和流动近似预测更复杂、更真实的水分输运情况，并能反映观察到的路面在使用中的行为。这项工作促进了对非饱和沥青混合料中水分运移的理解，并为评估路面-地下水系统的长期水力性能提供了有用的工具。

{"title":"A modified capillary pressure–saturation model for unsaturated flow in asphalt mixtures and hydraulic behavior analysis","authors":"Xinxing Bian , Huining Xu , Shiyuan Li","doi":"10.1016/j.advwatres.2025.105203","DOIUrl":"10.1016/j.advwatres.2025.105203","url":null,"abstract":"<div><div>The capillary pressure–saturation (<em>P</em><sub>c</sub><em>–S</em>) relationship is a critical constitutive model for modeling multiphase flow in unsaturated porous media. In the context of asphalt pavements, however, existing <em>P</em><sub>c</sub><em>–S</em> models predominantly rely on direct adoption of empirical formulations in unsaturated soil, lack rigorous experimental validation, and often fail to capture the unsaturated hydraulic characteristics of pavement materials. This study aims to (1) develop a continuous mathematical framework for the <em>P</em><sub>c</sub><em>–S</em> relationship in asphalt mixtures by integrating the physical constraints derived from experimental data with the functional form of classical models, and (2) elucidate the impact of this relationship on unsaturated seepage dynamics. A modified van Genuchten model was formulated to provide a concise, physically consistent <em>P</em><sub>c</sub><em>–S</em> function. The proposed model demonstrates improved fidelity to experimental data compared to conventional soil-based hydraulic models, while preserving mathematical tractability. By implementing this <em>P</em><sub>c</sub><em>–S</em> constitutive relationship within a numerical solver (SEEP/W) for the Richards equation, the hysteretic behavior of unsaturated flow was simulated in asphalt mixtures. The simulations reveal that the unsaturated flow analysis, incorporating the proposed <em>P</em><sub>c</sub><em>–S</em> relationship, predicts a more complex and realistic moisture transport regime than the saturated flow approximation, and can reflect the observed in-service behavior of pavements. This work advances understanding of moisture transport in unsaturated asphalt mixtures and provides a useful tool for evaluating the long-term hydraulic performance in pavement–subsurface water systems.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"208 ","pages":"Article 105203"},"PeriodicalIF":4.2,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146035612","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}

引用次数: 0

Multi-mechanism modeling of gas-liquid hydraulic conductance 气液水力传导的多机理建模

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-02-01 Epub Date: 2026-01-09 DOI: 10.1016/j.advwatres.2026.105214

Qingyuan Zhu , Keliu Wu , Qiqi Liu , Fei Peng , Shengting Zhang , Jianlin Zhao , Jing Li , Japan Trivedi , Zhangxin Chen

Hydraulic conductance is crucial for determining fluid transport capacity in porous media during energy recovery and storage. Inspired by the quasi-steady-state assumption of seepage flows, we developed a numerical model coupling creeping flow, momentum transfer, and interfacial slip, enabling precise modeling of two-phase hydraulic conductance under various mechanisms. We generated tens of thousands of data samples and constructed three machine learning models for hydraulic conductance prediction. Validation results demonstrate that these models outperform traditional empirical models in accuracy. Viscous coupling enhances the wetting phase conductance more significantly than the non-wetting phase, with its effect controlled by fluid and geometric configuration. A dimensionless slip length exceeding 0.01 elevates wetting phase conductance but has limited impact on the non-wetting phase. Non-wetting phase slip boosts the hydraulic conductance of both phases at a Knudsen number exceeding 0.01. The wetting film profoundly reshapes the velocity profiles and amplifies viscous coupling effect, with wetting film flow contributing as much as 60∼100% to the total wetting phase flow. Wetting film also strengthens the impact of slip on both phases' conductance. These findings highlight the intricate mechanisms controlling gas-liquid hydraulic conductance, especially the powerful control of wetting film, which is often underestimated in current pore-scale to core-scale analyses.

在能量回收和储存过程中，水力导度是决定流体在多孔介质中传输能力的关键因素。受渗流准稳态假设的启发，我们建立了一个耦合蠕变流动、动量传递和界面滑移的数值模型，从而能够精确模拟各种机制下的两相水力导度。我们生成了数以万计的数据样本，并构建了三种用于水力导度预测的机器学习模型。验证结果表明，这些模型在精度上优于传统的经验模型。黏性耦合对润湿相电导率的提高比非润湿相更显著，其影响受流体和几何构型的控制。超过0.01的无因次滑移长度提高了润湿相的电导率，但对非润湿相的影响有限。非润湿相滑移提高了两相的水力导度，其无因次克努森数超过0.01。润湿膜深刻地重塑了速度分布，放大了粘性耦合效应，润湿膜流量对总润湿相流的贡献高达60 ~ 100%。湿膜也增强了滑移对两相电导的影响。这些发现强调了控制气液水力电导的复杂机制，特别是润湿膜的强大控制，这在目前的孔隙尺度到岩心尺度的分析中经常被低估。

{"title":"Multi-mechanism modeling of gas-liquid hydraulic conductance","authors":"Qingyuan Zhu , Keliu Wu , Qiqi Liu , Fei Peng , Shengting Zhang , Jianlin Zhao , Jing Li , Japan Trivedi , Zhangxin Chen","doi":"10.1016/j.advwatres.2026.105214","DOIUrl":"10.1016/j.advwatres.2026.105214","url":null,"abstract":"<div><div>Hydraulic conductance is crucial for determining fluid transport capacity in porous media during energy recovery and storage. Inspired by the quasi-steady-state assumption of seepage flows, we developed a numerical model coupling creeping flow, momentum transfer, and interfacial slip, enabling precise modeling of two-phase hydraulic conductance under various mechanisms. We generated tens of thousands of data samples and constructed three machine learning models for hydraulic conductance prediction. Validation results demonstrate that these models outperform traditional empirical models in accuracy. Viscous coupling enhances the wetting phase conductance more significantly than the non-wetting phase, with its effect controlled by fluid and geometric configuration. A dimensionless slip length exceeding 0.01 elevates wetting phase conductance but has limited impact on the non-wetting phase. Non-wetting phase slip boosts the hydraulic conductance of both phases at a Knudsen number exceeding 0.01. The wetting film profoundly reshapes the velocity profiles and amplifies viscous coupling effect, with wetting film flow contributing as much as 60∼100% to the total wetting phase flow. Wetting film also strengthens the impact of slip on both phases' conductance. These findings highlight the intricate mechanisms controlling gas-liquid hydraulic conductance, especially the powerful control of wetting film, which is often underestimated in current pore-scale to core-scale analyses.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"208 ","pages":"Article 105214"},"PeriodicalIF":4.2,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956967","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}

引用次数: 0

Chaotic advection in regional groundwater flow under periodic water table fluctuations: An analytical model with depth-dependent aquifer properties 地下水位周期性波动下区域地下水流动的混沌平流：一个含深度依赖含水层性质的解析模型

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-02-01 Epub Date: 2026-01-10 DOI: 10.1016/j.advwatres.2026.105215

Saurabh Maurya, Ratan Sarmah, Ickkshaanshu Sonkar

Reliable prediction of groundwater flow dynamics is often hindered by the assumption of uniform aquifer parameters, leading to biased estimates of hydraulic head. Despite their well-established depth-dependence relationship, most analytical solutions remain restricted to homogeneous or single-layer aquifer systems. This study develops a two-dimensional transient analytical model for a multilayer regional aquifer driven by a fluctuating water table, explicitly incorporating depth-dependent hydraulic conductivity and specific storage while accounting for pumping. The analytical solution is derived using the Generalized Integral Transform Technique (GITT), ensuring implicit continuity of head and flux across layer interfaces without requiring iterative eigenvalue estimation, providing a rigorous framework to capture the complexities of stratified aquifers. Model verification is performed against a benchmark single-layer solution, and independent validation is carried out with COMSOL Multiphysics simulations, showing excellent agreement. Parameter influence and reliability are further assessed through global sensitivity analysis and uncertainty evaluation. A novel contribution of this work is the identification of chaotic flow behaviour within a three-layer Tóthian basin, analyzed using the Finite-Time Lyapunov Exponent (FTLE). Results show that chaotic dynamics are most pronounced near the upper boundary, where FTLE values are significantly higher than in deeper zones and further enhancement is observed under periodic injection. Moreover, specific storage is found to enhance the swirly motion of particle trajectories, increasing residence times near pseudo-stagnation zones. The proposed analytical framework bridges a critical gap in multilayer transient flow modeling, advances the theoretical understanding of stratified aquifer systems, and provides a robust benchmark for numerical and field-scale groundwater investigations.

地下水流动动力学的可靠预测往往受到均匀含水层参数假设的阻碍，导致水头估计有偏差。尽管它们建立了深度依赖关系，但大多数解析解仍然局限于均质或单层含水层系统。本研究建立了一个由波动地下水位驱动的多层区域含水层的二维瞬态分析模型，明确地结合了与深度相关的水力传导性和特定储水量，同时考虑了抽水。利用广义积分变换技术（GITT）推导了解析解，保证了水头和通量跨层界面的隐式连续性，而不需要迭代特征值估计，为捕获分层含水层的复杂性提供了严格的框架。针对基准单层解决方案进行了模型验证，并通过COMSOL多物理场仿真进行了独立验证，结果显示出良好的一致性。通过全局敏感性分析和不确定度评估，进一步评估了参数的影响和可靠性。这项工作的一个新贡献是识别三层Tóthian盆地内的混沌流动行为，使用有限时间李雅普诺夫指数（FTLE）进行分析。结果表明，在上边界附近，混沌动力学最为明显，其中FTLE值明显高于深部区域，并且在周期性注入下进一步增强。此外，发现特定的存储增强了粒子轨迹的漩涡运动，增加了在伪停滞区附近的停留时间。所提出的分析框架弥补了多层瞬态流动模型的关键缺陷，推进了对分层含水层系统的理论认识，并为数值和现场尺度地下水调查提供了一个可靠的基准。

{"title":"Chaotic advection in regional groundwater flow under periodic water table fluctuations: An analytical model with depth-dependent aquifer properties","authors":"Saurabh Maurya, Ratan Sarmah, Ickkshaanshu Sonkar","doi":"10.1016/j.advwatres.2026.105215","DOIUrl":"10.1016/j.advwatres.2026.105215","url":null,"abstract":"<div><div>Reliable prediction of groundwater flow dynamics is often hindered by the assumption of uniform aquifer parameters, leading to biased estimates of hydraulic head. Despite their well-established depth-dependence relationship, most analytical solutions remain restricted to homogeneous or single-layer aquifer systems. This study develops a two-dimensional transient analytical model for a multilayer regional aquifer driven by a fluctuating water table, explicitly incorporating depth-dependent hydraulic conductivity and specific storage while accounting for pumping. The analytical solution is derived using the Generalized Integral Transform Technique (GITT), ensuring implicit continuity of head and flux across layer interfaces without requiring iterative eigenvalue estimation, providing a rigorous framework to capture the complexities of stratified aquifers. Model verification is performed against a benchmark single-layer solution, and independent validation is carried out with COMSOL Multiphysics simulations, showing excellent agreement. Parameter influence and reliability are further assessed through global sensitivity analysis and uncertainty evaluation. A novel contribution of this work is the identification of chaotic flow behaviour within a three-layer Tóthian basin, analyzed using the Finite-Time Lyapunov Exponent (FTLE). Results show that chaotic dynamics are most pronounced near the upper boundary, where FTLE values are significantly higher than in deeper zones and further enhancement is observed under periodic injection. Moreover, specific storage is found to enhance the swirly motion of particle trajectories, increasing residence times near pseudo-stagnation zones. The proposed analytical framework bridges a critical gap in multilayer transient flow modeling, advances the theoretical understanding of stratified aquifer systems, and provides a robust benchmark for numerical and field-scale groundwater investigations.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"208 ","pages":"Article 105215"},"PeriodicalIF":4.2,"publicationDate":"2026-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145956936","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}

引用次数: 0

Generative adversarial network-based super-resolution of subsurface rock images: Visual, petrophysical, and flow simulation assessment 基于生成对抗网络的地下岩石图像超分辨率：视觉、岩石物理和流体模拟评估

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-01-01 Epub Date: 2025-12-01 DOI: 10.1016/j.advwatres.2025.105184

Mohammad Hamidian, Rohaldin Miri, Hossein Fazeli

High-resolution (HR) micro-CT imaging of porous reservoir rocks plays a critical role in digital rock physics and flow simulations, yet it is limited by the resolution–field-of-view trade-off. To address this challenge, we propose a Pore-Preserving Denoised Generative Adversarial Network (PPD-GAN), trained on denoised, HR two-dimensional computed tomography (CT) slices to recover fine-scale pore structures from low-resolution (LR) images. The PPD-GAN model is systematically compared against seven convolutional neural network (CNN)-based SR methods and bicubic interpolation using both image-based metrics—including Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index Measure (SSIM) and Learned Perceptual Image Patch Similarity (LPIPS)—and a comprehensive set of petro-physical metrics, including porosity, connectivity, pore size distribution, skeleton topology, and grain boundary preservation. Results show that although CNN models yield high pixel-wise fidelity, the PPD-GAN model achieves superior perceptual quality and reconstructs structural features that are consistently closer to ground truth across all physical metrics. Furthermore, pore-scale transport simulations on three-dimensional (3D) core images confirm that PPD-GAN outputs accurately replicate ground-truth (GT) properties such as permeability, diffusivity, and relative permeability—substantially outperforming bicubic interpolation. These findings demonstrate that the proposed PPD-GAN model not only enhances visual resolution but also preserves physically meaningful characteristics, enabling reliable downstream simulation and analysis.

多孔储层岩石的高分辨率（HR）微ct成像在数字岩石物理和流动模拟中起着至关重要的作用，但它受到分辨率和视场权衡的限制。为了解决这一挑战，我们提出了一种保留孔隙的去噪生成对抗网络（PPD-GAN），该网络在去噪的HR二维计算机断层扫描（CT）切片上进行训练，以从低分辨率（LR）图像中恢复精细尺度的孔隙结构。利用基于图像的指标（包括峰值信噪比（PSNR）、结构相似指数测量（SSIM）和学习感知图像斑块相似度（LPIPS））和一套全面的岩石物理指标（包括孔隙度、连通性、孔径分布、骨架拓扑结构和晶粒边界保存），将PPD-GAN模型与7种基于卷积神经网络（CNN）的SR方法和双三次插值进行了系统比较。结果表明，尽管CNN模型产生了高像素保真度，但PPD-GAN模型实现了卓越的感知质量，并重建了在所有物理指标中始终更接近地面真相的结构特征。此外，三维（3D）岩心图像上的孔隙尺度输运模拟证实，PPD-GAN输出准确地复制了渗透率、扩散率和相对渗透率等地面真值（GT）属性，大大优于双三次插值。这些发现表明，所提出的PPD-GAN模型不仅提高了视觉分辨率，而且保留了物理上有意义的特征，从而实现了可靠的下游模拟和分析。

{"title":"Generative adversarial network-based super-resolution of subsurface rock images: Visual, petrophysical, and flow simulation assessment","authors":"Mohammad Hamidian, Rohaldin Miri, Hossein Fazeli","doi":"10.1016/j.advwatres.2025.105184","DOIUrl":"10.1016/j.advwatres.2025.105184","url":null,"abstract":"<div><div>High-resolution (HR) micro-CT imaging of porous reservoir rocks plays a critical role in digital rock physics and flow simulations, yet it is limited by the resolution–field-of-view trade-off. To address this challenge, we propose a Pore-Preserving Denoised Generative Adversarial Network (PPD-GAN), trained on denoised, HR two-dimensional computed tomography (CT) slices to recover fine-scale pore structures from low-resolution (LR) images. The PPD-GAN model is systematically compared against seven convolutional neural network (CNN)-based SR methods and bicubic interpolation using both image-based metrics—including Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index Measure (SSIM) and Learned Perceptual Image Patch Similarity (LPIPS)—and a comprehensive set of petro-physical metrics, including porosity, connectivity, pore size distribution, skeleton topology, and grain boundary preservation. Results show that although CNN models yield high pixel-wise fidelity, the PPD-GAN model achieves superior perceptual quality and reconstructs structural features that are consistently closer to ground truth across all physical metrics. Furthermore, pore-scale transport simulations on three-dimensional (3D) core images confirm that PPD-GAN outputs accurately replicate ground-truth (GT) properties such as permeability, diffusivity, and relative permeability—substantially outperforming bicubic interpolation. These findings demonstrate that the proposed PPD-GAN model not only enhances visual resolution but also preserves physically meaningful characteristics, enabling reliable downstream simulation and analysis.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105184"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145657972","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}

引用次数: 0

Decision-making under uncertainty for LT-ATES systems in complex subsurface settings: application to a low-transmissivity aquifer 复杂地下环境中LT-ATES系统的不确定性决策：在低透射率含水层中的应用

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-01-01 Epub Date: 2025-12-06 DOI: 10.1016/j.advwatres.2025.105193

Luka Tas , Jef Caers , Thomas Hermans

Low-temperature aquifer thermal energy storage systems (LT-ATES) store excess thermal energy in the subsurface and recover it when needed to heat and cool buildings sustainably. Complex subsurface settings are increasingly targeted for storage. Existing design guidelines cannot be readily applied because they were developed for homogeneous sandy aquifers. In complex settings, subsurface uncertainty translates into uncertainty about system feasibility or efficiency. This paper introduces a method to establish a decision tree for ATES based on uncertainty analysis. It combines Monte Carlo simulations of groundwater flow and heat transport with sensitivity analysis, joint probability distribution estimation and kernel density estimation. For illustration the method is applied to a low-transmissivity aquifer where hydraulic and thermal feasibility are the main prediction targets. Sensitivity analyses based on clustering reveal the impact of model parameters and engineering actions on both prediction targets. In the specific case study, we find for example that storage conditions with transmissivity below 20 m²/d lead to inefficient systems. Desirable storage conditions have transmissivity above 40 m²/d, while intermediate values require thorough risk analysis. Thermal breakthrough risk is higher when longitudinal dispersivity is above 3 m. Further, our approach results in some minimum requirements in terms of subsurface properties that have to be reached for which an investment is justified. On top, the decision tree proposes target engineering actions to decrease investment risk while optimizing return. As such, the proposed method effectively guides early investment decisions for LT-ATES in complex subsurface settings without requiring field data.

低温含水层热能储存系统（LT-ATES）将多余的热能储存在地下，并在需要时回收，以可持续地加热和冷却建筑物。复杂的地下环境越来越多地成为存储的目标。现有的设计准则不能轻易适用，因为它们是针对均质砂质含水层制定的。在复杂的环境中，地下的不确定性转化为系统可行性或效率的不确定性。介绍了一种基于不确定性分析的ATES决策树的建立方法。它将蒙特卡罗模拟地下水流动和热输运与灵敏度分析、联合概率分布估计和核密度估计相结合。为了说明该方法应用于低透射率含水层，其中水力和热可行性是主要预测目标。基于聚类的敏感性分析揭示了模型参数和工程行为对两个预测目标的影响。在具体的案例研究中，我们发现，例如，传输率低于20 m²/d的存储条件会导致系统效率低下。理想的储存条件下，透过率应高于40 m²/d，而中间值则需要进行彻底的风险分析。当纵向分散度大于3 m时，热突破风险较高。此外，我们的方法在地下属性方面产生了一些最低要求，必须达到这些要求才能证明投资是合理的。决策树提出目标工程行动，以降低投资风险，同时优化收益。因此，所提出的方法可以有效地指导LT-ATES在复杂地下环境中的早期投资决策，而无需现场数据。

{"title":"Decision-making under uncertainty for LT-ATES systems in complex subsurface settings: application to a low-transmissivity aquifer","authors":"Luka Tas , Jef Caers , Thomas Hermans","doi":"10.1016/j.advwatres.2025.105193","DOIUrl":"10.1016/j.advwatres.2025.105193","url":null,"abstract":"<div><div>Low-temperature aquifer thermal energy storage systems (LT-ATES) store excess thermal energy in the subsurface and recover it when needed to heat and cool buildings sustainably. Complex subsurface settings are increasingly targeted for storage. Existing design guidelines cannot be readily applied because they were developed for homogeneous sandy aquifers. In complex settings, subsurface uncertainty translates into uncertainty about system feasibility or efficiency. This paper introduces a method to establish a decision tree for ATES based on uncertainty analysis. It combines Monte Carlo simulations of groundwater flow and heat transport with sensitivity analysis, joint probability distribution estimation and kernel density estimation. For illustration the method is applied to a low-transmissivity aquifer where hydraulic and thermal feasibility are the main prediction targets. Sensitivity analyses based on clustering reveal the impact of model parameters and engineering actions on both prediction targets. In the specific case study, we find for example that storage conditions with transmissivity below 20 m²/d lead to inefficient systems. Desirable storage conditions have transmissivity above 40 m²/d, while intermediate values require thorough risk analysis. Thermal breakthrough risk is higher when longitudinal dispersivity is above 3 m. Further, our approach results in some minimum requirements in terms of subsurface properties that have to be reached for which an investment is justified. On top, the decision tree proposes target engineering actions to decrease investment risk while optimizing return. As such, the proposed method effectively guides early investment decisions for LT-ATES in complex subsurface settings without requiring field data.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105193"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145689898","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}

引用次数: 0

A Physics-Informed Adaptive-Wavelet Neural Network (PIAWNN) for reservoir simulation with two-phase flow 基于物理信息的自适应小波神经网络（PIAWNN）在两相流油藏模拟中的应用

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-01-01 Epub Date: 2025-11-14 DOI: 10.1016/j.advwatres.2025.105174

Haifeng Zhu , Zhiming Chen , Xin Gao , Xing Chen , Wei Yu , Kamy Sepehrnoori

The flow of oil and water represents the most common two-phase regime in reservoir development. The frontal position and diffusion law of water saturation and pressure are typically delineated using nonlinear partial differential equations. However, these equations exhibit strong nonlinearity. Conventional numerical approaches demand substantial computational resources, posing challenges in handling such equations. Recently, physics - informed neural networks (PINNs) have opened new avenues for solving partial differential equations, while existing PINNs methods still suffer from high computational expenses. To address this, we propose a physics - informed adaptive wavelet neural network (PIAWNN) as a simulator for oil - water two - phase flow. The solution of pressure at the subsequent time step is primarily achieved by iterating the loss function structured in the form of the finite volume method. Additionally, the implicit pressure explicit saturation method is employed to update saturation. For two-dimensional problems, results demonstrate that the mean absolute error and mean relative error reach a low level. When solve three-dimensional complex scenarios, we find the proposed PIAWNN remains applicable. Comparing training performances across different networks, a notable finding emerges: as grid size increases, PIAWNN's training time and accuracy do not deteriorate significantlyand it outperforms physics-informed convolutional neural network (PICNN) and PINN. These results suggest PIAWNN holds promise for yielding novel insights into solving nonlinear systems.

油水同时流动是油藏开发中最常见的两相流形式。含水饱和度和压力的锋面位置和扩散规律对地下资源的有效开发至关重要，这些现象通常用非线性偏微分方程来描述。但油水两相流的偏微分方程表现出较强的非线性。传统的数值方法，包括牛顿迭代法和有限元法，需要大量的计算资源来求解，这给处理这类方程带来了挑战。最近，基于物理的神经网络（PINN）为求解偏微分方程开辟了新的途径，而基于该框架的现有方法仍然存在相对较高的计算开销。为了解决这个问题，我们提出了一个物理信息自适应小波神经网络（PIAWNN）作为油水两相流的模拟器。后续时间步长的压力解主要通过迭代有限体积法形式的损失函数来实现。此外，采用隐式压力显式饱和方法更新饱和度，降低了全隐式算法的计算成本。对于二维问题，我们进行了非均质条件下一次注一次采和多次注一次采两种场景。结果表明，平均绝对误差和平均相对误差达到了较低的水平。进一步将研究扩展到三维单注多采，我们发现所提出的PIAWNN仍然是适用的。通过比较不同网络的训练性能并进行误差分析，一个值得注意的发现出现了：随着网格大小的增加，所提出的网络的训练时间和准确性并没有显着下降——这一特征优于当前流行的方法，如物理卷积神经网络（PICNN）和PINN。这些结果表明，所提出的PIAWNN有望为解决非线性系统引入新的见解。

{"title":"A Physics-Informed Adaptive-Wavelet Neural Network (PIAWNN) for reservoir simulation with two-phase flow","authors":"Haifeng Zhu , Zhiming Chen , Xin Gao , Xing Chen , Wei Yu , Kamy Sepehrnoori","doi":"10.1016/j.advwatres.2025.105174","DOIUrl":"10.1016/j.advwatres.2025.105174","url":null,"abstract":"<div><div>The flow of oil and water represents the most common two-phase regime in reservoir development. The frontal position and diffusion law of water saturation and pressure are typically delineated using nonlinear partial differential equations. However, these equations exhibit strong nonlinearity. Conventional numerical approaches demand substantial computational resources, posing challenges in handling such equations. Recently, physics - informed neural networks (PINNs) have opened new avenues for solving partial differential equations, while existing PINNs methods still suffer from high computational expenses. To address this, we propose a physics - informed adaptive wavelet neural network (PIAWNN) as a simulator for oil - water two - phase flow. The solution of pressure at the subsequent time step is primarily achieved by iterating the loss function structured in the form of the finite volume method. Additionally, the implicit pressure explicit saturation method is employed to update saturation. For two-dimensional problems, results demonstrate that the mean absolute error and mean relative error reach a low level. When solve three-dimensional complex scenarios, we find the proposed PIAWNN remains applicable. Comparing training performances across different networks, a notable finding emerges: as grid size increases, PIAWNN's training time and accuracy do not deteriorate significantlyand it outperforms physics-informed convolutional neural network (PICNN) and PINN. These results suggest PIAWNN holds promise for yielding novel insights into solving nonlinear systems.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105174"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145531299","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}

引用次数: 0

A simple derivation of the geometric theorems for averaging and Errata for the article Adv. Water Res., 32 (2013) 340-352 平均几何定理的简单推导及文章的勘误表。水利学报，32 (2013):340-352

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-01-01 Epub Date: 2025-12-04 DOI: 10.1016/j.advwatres.2025.105181

Francisco J. Valdés-Parada , Didier Lasseux , Brian D. Wood

The derivation of upscaled models using the volume averaging method involves the use of the spatial averaging theorem. From its use, further geometrical identities arise when it is applied to a position vector or its dyadic products at the

n

th-order, i.e., the

n

-adic product. The derivation of these geometrical identities is revised in this work in a concise manner that corrects some inconsistencies in the technical note by Wood (2013). These new derivations strengthen the use of the volume averaging method in many transfer problems, in particular in porous media.

用体积平均法推导上尺度模型涉及到空间平均定理的使用。从它的使用，进一步的几何恒等式产生，当它被应用到一个位置向量或其在n阶的并矢积，即n进积。这些几何恒等式的推导在这项工作中以简洁的方式进行了修订，纠正了Wood（2013）的技术说明中的一些不一致之处。这些新的推导加强了体积平均法在许多传递问题中的应用，特别是在多孔介质中。

引用次数: 0

Analytical investigation of pumping in finite coastal aquifers under heterogeneous recharge and various boundary conditions 非均质补给和不同边界条件下有限海岸含水层抽水分析研究

IF 4.2 2区环境科学与生态学 Q1 WATER RESOURCES

Advances in Water Resources

Pub Date : 2026-01-01 Epub Date: 2025-11-08 DOI: 10.1016/j.advwatres.2025.105171

Yuening Tang , Mengfan Wang , Chunhui Lu , Jian Luo

Effective pumping management is critical for coastal aquifers to address high freshwater demand in densely populated areas while mitigating seawater intrusion risks. Previous analytical studies have not considered the effects of heterogeneous recharge for pumping in coastal aquifers. We present the first analytical solutions for groundwater withdrawal under heterogeneous recharge conditions in finite rectangular aquifers with diverse boundary conditions. The solutions allow determination of toe positions for single- and multi-well systems, as well as the critical pumping rate for single wells. Results show that the maximum allowable pumping rate increases with the domain aspect ratio and eventually stabilizes, with stabilized values differing by <5 % across domains with varying numbers of sea boundaries for a centrally located well at an aspect ratio of 5. Boundary conditions exert a distinct influence on system behavior. The critical pumping rates of wells adjacent to sea boundaries show minimal sensitivity to heterogeneous recharge patterns, whereas those near no-flow boundaries exhibit high sensitivity. Uncertainty analysis confirms the contrasting influences of these two boundary types and further demonstrates that inter-percentile differences remain largely insensitive to recharge heterogeneity under the studied configurations. These analytical solutions extend classical potential theory to more realistic recharge and boundary conditions, offering a computationally efficient and physically transparent tool for preliminary design, sensitivity analysis, and sustainable groundwater management in coastal aquifers.

有效的抽水管理对于解决人口密集地区的高淡水需求，同时减轻海水入侵风险的沿海含水层至关重要。以往的分析研究没有考虑非均质补给对沿海含水层抽水的影响。本文首次提出了具有不同边界条件的有限矩形含水层在非均质补给条件下地下水抽取的解析解。该解决方案可以确定单井和多井系统的趾部位置，以及单井的临界泵速。结果表明，当长径比为5时，最大允许抽油速率随区域长径比的增加而增加，并最终趋于稳定，对于位于中心位置的井来说，在不同海域边界数量的区域，稳定值相差约5%。边界条件对系统行为有明显的影响。靠近海界的井的临界抽速对非均质补给模式的敏感性最低，而靠近无流界的井的临界抽速对非均质补给模式的敏感性很高。不确定性分析证实了这两种边界类型的不同影响，并进一步表明在研究配置下，百分位间差异对补给非均质性基本不敏感。这些解析解将经典电位理论扩展到更现实的补给和边界条件，为沿海含水层的初步设计、敏感性分析和可持续地下水管理提供了一种计算效率高、物理透明的工具。

{"title":"Analytical investigation of pumping in finite coastal aquifers under heterogeneous recharge and various boundary conditions","authors":"Yuening Tang , Mengfan Wang , Chunhui Lu , Jian Luo","doi":"10.1016/j.advwatres.2025.105171","DOIUrl":"10.1016/j.advwatres.2025.105171","url":null,"abstract":"<div><div>Effective pumping management is critical for coastal aquifers to address high freshwater demand in densely populated areas while mitigating seawater intrusion risks. Previous analytical studies have not considered the effects of heterogeneous recharge for pumping in coastal aquifers. We present the first analytical solutions for groundwater withdrawal under heterogeneous recharge conditions in finite rectangular aquifers with diverse boundary conditions. The solutions allow determination of toe positions for single- and multi-well systems, as well as the critical pumping rate for single wells. Results show that the maximum allowable pumping rate increases with the domain aspect ratio and eventually stabilizes, with stabilized values differing by <5 % across domains with varying numbers of sea boundaries for a centrally located well at an aspect ratio of 5. Boundary conditions exert a distinct influence on system behavior. The critical pumping rates of wells adjacent to sea boundaries show minimal sensitivity to heterogeneous recharge patterns, whereas those near no-flow boundaries exhibit high sensitivity. Uncertainty analysis confirms the contrasting influences of these two boundary types and further demonstrates that inter-percentile differences remain largely insensitive to recharge heterogeneity under the studied configurations. These analytical solutions extend classical potential theory to more realistic recharge and boundary conditions, offering a computationally efficient and physically transparent tool for preliminary design, sensitivity analysis, and sustainable groundwater management in coastal aquifers.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105171"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145473287","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}

引用次数: 0