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

Generative modeling of three-dimensional multiphase flow in porous media 多孔介质中三维多相流的生成建模

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

Advances in Water Resources

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

Serveh Kamrava

Understanding multiphase flow in porous media requires models that capture complex pore geometries and reproduce realistic fluid distributions under diverse boundary conditions. In this paper, we propose a generative framework that combines variational machine learning with a denoising diffusion approach to build 3D multiphase pore structures directly from micro-CT images. The model is trained on small subvolumes to balance resolution with computational feasibility, enabling efficient training and rapid generation of new realizations. Quantitative comparisons show good agreement between generated and experimental samples across morphological metrics and statistical functions, while permeability distributions are reproduced within the range of variability. The variational component of the proposed method provides a compact latent representation that accelerates sampling and also increases the diversity of multiphase configurations that the model can generate, capturing a wider range of pore-scale fluid distributions. The framework allows tiling strategies to produce larger domains that remain consistent with pore-scale data, to scale from local heterogeneity to larger representative volumes. This method generates ensembles of realistic multiphase structures that can be used as input for sampling, conditioning, and optimization workflows in porous media applications.

理解多孔介质中的多相流需要能够捕捉复杂孔隙几何形状的模型，并在不同的边界条件下重现真实的流体分布。在本文中，我们提出了一个结合变分机器学习和去噪扩散方法的生成框架，直接从微ct图像中构建3D多相孔隙结构。该模型在小子体积上进行训练，以平衡分辨率和计算可行性，从而实现高效训练和快速生成新实现。定量比较表明，在形态度量和统计函数上，生成样品和实验样品之间具有良好的一致性，而渗透率分布在变异性范围内重现。该方法的变分成分提供了一个紧凑的潜在表示，加速了采样，也增加了模型可以生成的多相配置的多样性，捕获了更大范围的孔隙尺度流体分布。该框架允许平铺策略产生与孔隙尺度数据保持一致的更大域，从局部异质性扩展到更大的代表性体积。该方法可生成实际多相结构的集成，可作为多孔介质应用中采样、调节和优化工作流程的输入。

{"title":"Generative modeling of three-dimensional multiphase flow in porous media","authors":"Serveh Kamrava","doi":"10.1016/j.advwatres.2025.105196","DOIUrl":"10.1016/j.advwatres.2025.105196","url":null,"abstract":"<div><div>Understanding multiphase flow in porous media requires models that capture complex pore geometries and reproduce realistic fluid distributions under diverse boundary conditions. In this paper, we propose a generative framework that combines variational machine learning with a denoising diffusion approach to build 3D multiphase pore structures directly from micro-CT images. The model is trained on small subvolumes to balance resolution with computational feasibility, enabling efficient training and rapid generation of new realizations. Quantitative comparisons show good agreement between generated and experimental samples across morphological metrics and statistical functions, while permeability distributions are reproduced within the range of variability. The variational component of the proposed method provides a compact latent representation that accelerates sampling and also increases the diversity of multiphase configurations that the model can generate, capturing a wider range of pore-scale fluid distributions. The framework allows tiling strategies to produce larger domains that remain consistent with pore-scale data, to scale from local heterogeneity to larger representative volumes. This method generates ensembles of realistic multiphase structures that can be used as input for sampling, conditioning, and optimization workflows in porous media applications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105196"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731387","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

Impact of observation and surrogate-model noises on deep learning-based subsurface heterogeneous structure identification through monitoring network optimization 观测噪声和替代模型噪声对基于深度学习的地下非均质结构监测网络优化识别的影响

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

Advances in Water Resources

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

Yuzhou Xia , Chuanjun Zhan , Zhenxue Dai , Jichun Wu , Xiaoying Zhang , Huichao Yin , Jiahe Yan , Junjun Chen , Zihao Wang , Mohamad Reza Soltanian , Kenneth C. Carroll

Accurate identification of subsurface structures is essential for enhancing geoscientific modeling and monitoring, thereby deepening the understanding of earth system. Although recent advances in deep learning have greatly enhanced the computational efficiency of structure identification, the identification accuracy can still be degraded by the spatial stochasticity of monitoring locations, observation noise and surrogate model noise. To address these challenges, this study develops a transition probability and entropy-based monitoring network optimization method, which can determine monitoring locations with high data worth. Based on this method, a synthetic contaminant transport experiment was conducted to obtain dynamic observations of hydraulic head, solute concentration, and electrical resistivity tomography (ERT). The impacts of observation and model noise on inversion results for different data types were then quantitatively analyzed. The results indicate that ERT data exhibit strong noise resistance, with average identification accuracy reductions of only 1.59% and 2.24% under 15% observation noise and model noise, respectively. Moreover, incorporating ERT data into data fusion scenarios enhances the robustness of other observation types, particularly solute concentration data, for which the accuracy decrease is limited to 1.24% even under the highest observation noise level. Model noise not only reduces the accuracy of structure identification but also increases uncertainty, and this negative effect is further amplified in data fusion scenarios.

准确识别地下构造对于加强地球科学建模和监测，从而加深对地球系统的认识至关重要。尽管近年来深度学习的发展大大提高了结构识别的计算效率，但由于监测位置的空间随机性、观测噪声和替代模型噪声的影响，识别精度仍然会降低。针对这些挑战，本研究开发了一种基于转移概率和熵的监测网络优化方法，该方法可以确定具有高数据价值的监测位置。基于该方法，进行了合成污染物输运实验，获得了水头、溶质浓度和电阻率层析成像（ERT）的动态观测结果。然后定量分析了观测噪声和模型噪声对不同数据类型反演结果的影响。结果表明，ERT数据具有较强的抗噪能力，在15%的观测噪声和模型噪声下，ERT数据的平均识别准确率分别仅降低1.59%和2.24%。此外，将ERT数据纳入数据融合场景，增强了其他观测类型的鲁棒性，特别是溶质浓度数据，即使在最高观测噪声水平下，其精度下降也限制在1.24%以内。模型噪声不仅降低了结构识别的精度，而且增加了不确定性，这种负面影响在数据融合场景中进一步放大。

{"title":"Impact of observation and surrogate-model noises on deep learning-based subsurface heterogeneous structure identification through monitoring network optimization","authors":"Yuzhou Xia , Chuanjun Zhan , Zhenxue Dai , Jichun Wu , Xiaoying Zhang , Huichao Yin , Jiahe Yan , Junjun Chen , Zihao Wang , Mohamad Reza Soltanian , Kenneth C. Carroll","doi":"10.1016/j.advwatres.2025.105204","DOIUrl":"10.1016/j.advwatres.2025.105204","url":null,"abstract":"<div><div>Accurate identification of subsurface structures is essential for enhancing geoscientific modeling and monitoring, thereby deepening the understanding of earth system. Although recent advances in deep learning have greatly enhanced the computational efficiency of structure identification, the identification accuracy can still be degraded by the spatial stochasticity of monitoring locations, observation noise and surrogate model noise. To address these challenges, this study develops a transition probability and entropy-based monitoring network optimization method, which can determine monitoring locations with high data worth. Based on this method, a synthetic contaminant transport experiment was conducted to obtain dynamic observations of hydraulic head, solute concentration, and electrical resistivity tomography (ERT). The impacts of observation and model noise on inversion results for different data types were then quantitatively analyzed. The results indicate that ERT data exhibit strong noise resistance, with average identification accuracy reductions of only 1.59% and 2.24% under 15% observation noise and model noise, respectively. Moreover, incorporating ERT data into data fusion scenarios enhances the robustness of other observation types, particularly solute concentration data, for which the accuracy decrease is limited to 1.24% even under the highest observation noise level. Model noise not only reduces the accuracy of structure identification but also increases uncertainty, and this negative effect is further amplified in data fusion scenarios.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105204"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880963","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

Multiscale projection-based embedded discrete fracture modeling approach for CO2 storage in deep saline aquifers 基于多尺度投影的深部含盐含水层CO2封存嵌入离散裂缝建模方法

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

Advances in Water Resources

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

Mengjie Zhao , Ryan Haagenson , Marc Gerritsma , Hadi Hajibeygi

This study introduces a multiscale simulation framework, termed Projection-based Embedded Discrete Fracture Modeling with Algebraic Dynamic Multilevel method (pEDFM-ADM), which integrates an embedded discrete fracture network representation with a fully algebraic, front-tracking-based mesh adaptation strategy. Incorporating a fully implicit scheme, compositional thermodynamics, and algebraic multilevel operators, the framework captures essential subsurface processes such as buoyancy-driven migration, convective dissolution, phase partitioning, and fracture-matrix interactions under geologically realistic conditions. The method constructs a hierarchy of multilevel grids and localized multiscale basis functions that introduce fine-scale heterogeneities at each coarse level. Adaptive mesh refinement and coarsening are driven by local variations in CO₂ mass fraction and executed through algebraic prolongation and restriction operators, enabling efficient projection between grid levels. The framework is systematically evaluated across a sequence of test cases with increasing complexity, including systems with low-permeability flow barriers, highly conductive fractures, striking a trade-off between computational resource and detailed simulation accuracy. Overall, the pEDFM-ADM framework provides a scalable, fully algebraic, and physically adaptive modeling tool for large-scale CO₂ storage simulations in fractured porous media, supporting predictive simulation and risk assessment for long-term carbon sequestration.

本研究引入了一种多尺度模拟框架，称为基于投影的嵌入式离散裂缝建模与代数动态多级方法（pEDFM-ADM），该框架将嵌入式离散裂缝网络表示与完全代数的、基于前端跟踪的网格自适应策略集成在一起。该框架结合了完全隐式方案、成分热力学和代数多层运算符，捕捉了地质现实条件下的基本地下过程，如浮力驱动的运移、对流溶解、相分配和裂缝-基质相互作用。该方法构建了多层次网格和局部多尺度基函数的层次结构，在每个粗层次上引入了精细尺度的异质性。自适应网格细化和粗化是由二氧化碳质量分数的局部变化驱动的，并通过代数扩展和限制算子执行，从而实现网格层之间的有效投影。该框架通过一系列复杂的测试案例进行系统评估，包括具有低渗透流障、高导流裂缝的系统，在计算资源和详细的模拟精度之间进行权衡。总体而言，pEDFM-ADM框架为裂缝性多孔介质中大规模二氧化碳封存模拟提供了一个可扩展的、完全代数的、物理自适应的建模工具，支持长期碳封存的预测模拟和风险评估。

{"title":"Multiscale projection-based embedded discrete fracture modeling approach for CO2 storage in deep saline aquifers","authors":"Mengjie Zhao , Ryan Haagenson , Marc Gerritsma , Hadi Hajibeygi","doi":"10.1016/j.advwatres.2025.105200","DOIUrl":"10.1016/j.advwatres.2025.105200","url":null,"abstract":"<div><div>This study introduces a multiscale simulation framework, termed Projection-based Embedded Discrete Fracture Modeling with Algebraic Dynamic Multilevel method (pEDFM-ADM), which integrates an embedded discrete fracture network representation with a fully algebraic, front-tracking-based mesh adaptation strategy. Incorporating a fully implicit scheme, compositional thermodynamics, and algebraic multilevel operators, the framework captures essential subsurface processes such as buoyancy-driven migration, convective dissolution, phase partitioning, and fracture-matrix interactions under geologically realistic conditions. The method constructs a hierarchy of multilevel grids and localized multiscale basis functions that introduce fine-scale heterogeneities at each coarse level. Adaptive mesh refinement and coarsening are driven by local variations in CO<sub>2</sub> mass fraction and executed through algebraic prolongation and restriction operators, enabling efficient projection between grid levels. The framework is systematically evaluated across a sequence of test cases with increasing complexity, including systems with low-permeability flow barriers, highly conductive fractures, striking a trade-off between computational resource and detailed simulation accuracy. Overall, the pEDFM-ADM framework provides a scalable, fully algebraic, and physically adaptive modeling tool for large-scale CO<sub>2</sub> storage simulations in fractured porous media, supporting predictive simulation and risk assessment for long-term carbon sequestration.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105200"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823478","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 statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials 基于孔隙尺寸分布估算饱和岩土渗透率的统计和模拟模型

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

Advances in Water Resources

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

Mengwei Liu, Yongkoo Seol

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. This study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ_d², where ϕ is porosity, μ_d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ_d/μ_d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ_dand a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

准确的渗透率估计在地下工程应用中是必不可少的，但由于天然岩土材料的复杂孔隙结构，仍然具有挑战性。传统的经验方法和简化的理论模型往往不能充分反映孔径分布和连通性的作用。本研究开发了一个基于统计和模拟的渗透率模型，该模型将孔隙尺度复杂性分解为k = αϕμd2的紧凑尺度，其中φ为孔隙度，μd为平均孔径，α为弱变化系数。通过将孔隙网络模拟与单峰和双峰孔径分布的统计分析相结合，我们发现了三个关键发现：(1)渗透率对平均孔径比孔隙度更敏感；（ii）在广泛的数据集中，σd/μd（标准差与平均值）的比值在特征值~ 0.4附近聚集，使得整个孔径分布的影响可以用μd和狭窄变化的α≈0.05表示；（iii）对于双峰系统，存在一个小孔隙的临界分数~ 0.78，在此分数以上，流动以小孔隙为主，从而可以定义有效的控制流动的孔隙种群，并简化此类系统的渗透率估计。该模型只需要孔隙度和具有代表性的平均孔径作为输入，并通过综合实验数据集（>；1700个样本）对不同土壤和岩石进行了验证，取得了良好的预测精度。总的来说，这项工作为地下工程、环境保护和资源管理应用提供了一种物理接地但实际上简单的渗透率估算器。

{"title":"A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials","authors":"Mengwei Liu, Yongkoo Seol","doi":"10.1016/j.advwatres.2025.105197","DOIUrl":"10.1016/j.advwatres.2025.105197","url":null,"abstract":"<div><div>Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. This study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form <em>k</em> = <em>αϕμ<sub>d</sub></em><sup>2</sup>, where <em>ϕ</em> is porosity, <em>μ<sub>d</sub></em> is mean pore size, and <em>α</em> is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio <em>σ<sub>d</sub></em>/<em>μ<sub>d</sub></em> (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by <em>μ<sub>d</sub></em>and a narrowly varying <em>α</em> ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105197"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880964","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

Enhancing hydrological prediction in snow-dominant basins through multivariate calibration-assimilation framework 利用多元校准-数据同化框架加强积雪优势流域水文预测

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

Advances in Water Resources

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

Fatemeh Rezaei Aderyani , Keyhan Gavahi , Hamid Moradkhani

This study introduces a novel multivariate data assimilation framework for snow-dominant basins to enhance hydrologic prediction capabilities. Applied to 94 snow-dominant watersheds from the CAMELS dataset over a 30-year period, our approach first implements a multivariate calibration strategy that optimizes Snow-17 model parameters using observed snow water equivalent (SWE) and SACSMA model parameters using observed streamflow. This calibration strategy improved SWE prediction while maintaining streamflow performance. Following this calibration, we implemented a multivariate data assimilation using the Evolutionary Particle Filter with Markov Chain Monte Carlo (EPFM) algorithm. This system simultaneously integrates both SWE and streamflow observations with the models, enhancing prediction accuracy according to multiple performance measure. The framework demonstrates exceptional seasonal performance, with peak SWE prediction skill in winter and improved accuracy in predicting peak flow during spring snowmelt. Over 90 % of the basins showed improved peak flow prediction, critical for water resources planning and management. A comprehensive assessment across diverse geographic regions, hydroclimatic conditions, and multiple metrics confirms the framework's robustness and versatility.

为提高积雪优势流域的水文预报能力，提出了一种新的多变量数据同化框架。该方法应用于骆驼数据集30年期间的94个雪主导流域，首先实现了一种多变量校准策略，利用观测到的雪水当量（SWE）优化snow- 17模型参数，利用观测到的流量优化SACSMA模型参数。该校准策略在保持水流性能的同时提高了SWE预测。在此校准之后，我们使用带有马尔可夫链蒙特卡罗（EPFM）算法的进化粒子滤波实现了多变量数据同化。该系统同时将SWE和流量观测数据与模型相结合，根据多种性能指标提高了预测精度。该框架表现出优异的季节性能，在冬季具有峰值SWE预测能力，而在春季融雪期间预测峰值流量的准确性有所提高。超过90%的流域的峰值流量预测得到改善，这对水资源规划和管理至关重要。对不同地理区域、水文气候条件和多种指标的综合评估证实了该框架的稳健性和通用性。

{"title":"Enhancing hydrological prediction in snow-dominant basins through multivariate calibration-assimilation framework","authors":"Fatemeh Rezaei Aderyani , Keyhan Gavahi , Hamid Moradkhani","doi":"10.1016/j.advwatres.2025.105192","DOIUrl":"10.1016/j.advwatres.2025.105192","url":null,"abstract":"<div><div>This study introduces a novel multivariate data assimilation framework for snow-dominant basins to enhance hydrologic prediction capabilities. Applied to 94 snow-dominant watersheds from the CAMELS dataset over a 30-year period, our approach first implements a multivariate calibration strategy that optimizes Snow-17 model parameters using observed snow water equivalent (SWE) and SACSMA model parameters using observed streamflow. This calibration strategy improved SWE prediction while maintaining streamflow performance. Following this calibration, we implemented a multivariate data assimilation using the Evolutionary Particle Filter with Markov Chain Monte Carlo (EPFM) algorithm. This system simultaneously integrates both SWE and streamflow observations with the models, enhancing prediction accuracy according to multiple performance measure. The framework demonstrates exceptional seasonal performance, with peak SWE prediction skill in winter and improved accuracy in predicting peak flow during spring snowmelt. Over 90 % of the basins showed improved peak flow prediction, critical for water resources planning and management. A comprehensive assessment across diverse geographic regions, hydroclimatic conditions, and multiple metrics confirms the framework's robustness and versatility.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105192"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145689886","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

Convection behavior of interlayer disturbance during salt cavern leaching process 盐洞浸出过程中层间扰动的对流行为

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.105183

Huan Li , Jie Yang , Hao Wang , Haoran Zhang , Liangliang Ye , Hongling Ma , Xiangsheng Chen , Xilin Shi

Insoluble interlayers in bedded salt formations cause irregular cavern morphology during solution mining, yet the underlying hydrodynamics remain poorly quantified. This study employs a synergistic approach, combining Particle Image Velocimetry (PIV) with numerical modeling to quantitatively link interlayer attributes to disturbances in the brine flow field. The effects of interlayer collapse on brine flow field dynamics, including range, thickness, number, spacing, and tubing position, are systematically analyzed. The results indicate that the interlayer collapse range, rather than its thickness, is the dominant control on flow segmentation. For a single interlayer, a high overhanging degree (η > 0.4) combined with positioning the tubing string above the interlayer results in severe partitioning of the flow field, confining over 80% of the convective energy to the lower cavern. This partitioning effect is significantly mitigated when the tubing is placed below the interlayer or when the overhanging ratio is reduced below a critical threshold of 0.3. Furthermore, we demonstrate that multi-interlayer structures correct the primary flow deviation but induce localized, high-vorticity recirculation zones between layers, explaining the formation of wall irregularities. Mechanistically, we identified a dual-scale erosion process: a global turbulent flow drives bulk dissolution, while localized deflected jets sculpt the deep cavities. These findings provide critical theoretical guidance for optimizing the construction and ensuring the morphological stability of salt caverns in complex, bedded salt formations.

在溶采过程中，层状盐层中的不溶性夹层导致了不规则的洞穴形态，但底层的水动力学仍然很难量化。本研究采用协同方法，将颗粒图像测速（PIV）与数值模拟相结合，定量地将层间属性与盐水流场中的扰动联系起来。系统分析了层间塌缩对卤水流场动力学的影响，包括范围、厚度、数量、间距、油管位置等。结果表明，层间塌陷范围而非层间厚度是控制流动分割的主要因素。对于单个夹层，高悬垂度（η > 0.4）加上将管柱置于夹层上方，会导致流场的严重分配，将80%以上的对流能限制在下层洞穴中。当油管位于夹层下方或悬垂比低于0.3的临界阈值时，这种分割效果会显著减轻。此外，我们证明多层间结构纠正了初级流动偏差，但在层间诱导了局部的高涡度再循环区，解释了壁面不规则的形成。在力学上，我们确定了一个双尺度的侵蚀过程：全球湍流驱动整体溶解，而局部偏转射流雕刻深空腔。这些发现为复杂层状盐层中盐洞的优化构造和形态稳定性提供了重要的理论指导。

{"title":"Convection behavior of interlayer disturbance during salt cavern leaching process","authors":"Huan Li , Jie Yang , Hao Wang , Haoran Zhang , Liangliang Ye , Hongling Ma , Xiangsheng Chen , Xilin Shi","doi":"10.1016/j.advwatres.2025.105183","DOIUrl":"10.1016/j.advwatres.2025.105183","url":null,"abstract":"<div><div>Insoluble interlayers in bedded salt formations cause irregular cavern morphology during solution mining, yet the underlying hydrodynamics remain poorly quantified. This study employs a synergistic approach, combining Particle Image Velocimetry (PIV) with numerical modeling to quantitatively link interlayer attributes to disturbances in the brine flow field. The effects of interlayer collapse on brine flow field dynamics, including range, thickness, number, spacing, and tubing position, are systematically analyzed. The results indicate that the interlayer collapse range, rather than its thickness, is the dominant control on flow segmentation. For a single interlayer, a high overhanging degree (<em>η</em> > 0.4) combined with positioning the tubing string above the interlayer results in severe partitioning of the flow field, confining over 80% of the convective energy to the lower cavern. This partitioning effect is significantly mitigated when the tubing is placed below the interlayer or when the overhanging ratio is reduced below a critical threshold of 0.3. Furthermore, we demonstrate that multi-interlayer structures correct the primary flow deviation but induce localized, high-vorticity recirculation zones between layers, explaining the formation of wall irregularities. Mechanistically, we identified a dual-scale erosion process: a global turbulent flow drives bulk dissolution, while localized deflected jets sculpt the deep cavities. These findings provide critical theoretical guidance for optimizing the construction and ensuring the morphological stability of salt caverns in complex, bedded salt formations.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105183"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145657973","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

The impact of bare gaps on meadow-scale flows within fragmented emergent aquatic vegetation 裸地间隙对破碎化突生水生植被草地尺度流动的影响

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

Advances in Water Resources

Pub Date : 2026-01-01 Epub Date: 2025-11-27 DOI: 10.1016/j.advwatres.2025.105182

T.D. da Silva , J.C. Mullarney , C.A. Pilditch , G. Coco

Aquatic vegetation has been shown to shape landscape morphology by interacting with flows and influencing sedimentological processes. These interactions are modulated by the spatial distribution and density of the vegetation. We conducted idealised numerical simulations using a depth-averaged hydrodynamic model to explore how circular gaps (of diameters 0.4 to 9 m) within otherwise continuous beds of emergent aquatic vegetation alter in-canopy horizontal flow patterns. Our findings indicate that gaps reroute and concentrate horizontal flow, creating areas of accelerated flows within the canopy. As gap size increases, horizontal flow speeds also increase, though the rate of increase diminishes substantially for gap diameters larger than 5 m. In cases with paired gaps, two distinct interaction modes were identified, one occurring when the gaps were aligned with the flow, leading to mutually reinforced flow intensification within the gaps, and another when they were positioned side by side, resulting in diminished flow acceleration associated with each bare gap. The maximum flow change within the canopy occurred at streamwise gap separations of three gap diameters, with minimal changes for gap separations of greater than four diameters. The increase in horizontal flow speeds was directly proportional to the percentage of bare bed cover. For a given percentage of bare bed cover ranging from 8 to 40 %, flow changes were also modulated by the fragmentation level, which is determined by the number and size of gaps. The results reveal that fragmentation levels control in-canopy flow speeds, whereby fewer, larger gaps exerted a more pronounced impact on in-canopy horizontal flow speeds compared to numerous smaller gaps. The interconnected nature of vegetation coverage and gap distribution underscores the role of fragmentation levels in ecohydrological dynamics.

水生植被已被证明通过与水流相互作用和影响沉积过程来塑造景观形态。这些相互作用受到植被的空间分布和密度的调节。我们使用深度平均水动力模型进行了理想的数值模拟，以探索突生水生植被连续床内的圆形间隙（直径0.4至9米）如何改变冠层内的水平流动模式。我们的研究结果表明，间隙改变和集中了水平流动，在冠层内创造了加速流动的区域。随着间隙尺寸的增大，水平流速度也会增加，但当间隙直径大于5 m时，水平流速度的增加幅度会明显减小。在成对间隙的情况下，确定了两种不同的相互作用模式，一种发生在间隙与流动对齐时，导致间隙内相互增强的流动加剧，另一种发生在它们并排放置时，导致与每个裸间隙相关的流动加速减弱。冠层内流量变化最大的是3个直径的顺流间隙，大于4个直径的间隙变化最小。水平流动速度的增加与裸床覆盖的百分比成正比。对于一定比例的裸床覆盖物，流量变化也受到破碎程度的调节，而破碎程度是由间隙的数量和大小决定的。结果表明，破碎程度控制着冠层内的流动速度，相对于众多的小间隙，较少的大间隙对冠层内水平流动速度的影响更为显著。植被覆盖度和林隙分布的相互联系强调了破碎化水平在生态水文动力学中的作用。

{"title":"The impact of bare gaps on meadow-scale flows within fragmented emergent aquatic vegetation","authors":"T.D. da Silva , J.C. Mullarney , C.A. Pilditch , G. Coco","doi":"10.1016/j.advwatres.2025.105182","DOIUrl":"10.1016/j.advwatres.2025.105182","url":null,"abstract":"<div><div>Aquatic vegetation has been shown to shape landscape morphology by interacting with flows and influencing sedimentological processes. These interactions are modulated by the spatial distribution and density of the vegetation. We conducted idealised numerical simulations using a depth-averaged hydrodynamic model to explore how circular gaps (of diameters 0.4 to 9 m) within otherwise continuous beds of emergent aquatic vegetation alter in-canopy horizontal flow patterns. Our findings indicate that gaps reroute and concentrate horizontal flow, creating areas of accelerated flows within the canopy. As gap size increases, horizontal flow speeds also increase, though the rate of increase diminishes substantially for gap diameters larger than 5 m. In cases with paired gaps, two distinct interaction modes were identified, one occurring when the gaps were aligned with the flow, leading to mutually reinforced flow intensification within the gaps, and another when they were positioned side by side, resulting in diminished flow acceleration associated with each bare gap. The maximum flow change within the canopy occurred at streamwise gap separations of three gap diameters, with minimal changes for gap separations of greater than four diameters. The increase in horizontal flow speeds was directly proportional to the percentage of bare bed cover. For a given percentage of bare bed cover ranging from 8 to 40 %, flow changes were also modulated by the fragmentation level, which is determined by the number and size of gaps. The results reveal that fragmentation levels control in-canopy flow speeds, whereby fewer, larger gaps exerted a more pronounced impact on in-canopy horizontal flow speeds compared to numerous smaller gaps. The interconnected nature of vegetation coverage and gap distribution underscores the role of fragmentation levels in ecohydrological dynamics.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105182"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145611783","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

Ephemeral reservoirs of low-head wadi dams drained by seepage, evaporation and pumping: the Pavlovsky/Polubarinova-Kochina/Abel analytical legacy redux 渗流、蒸发和抽水排干的低水头瓦底大坝的短暂水库：Pavlovsky/Polubarinova-Kochina/Abel分析遗产Redux

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

Advances in Water Resources

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

Anvar Kacimov

Transient 2-D seepage through a porous bed of a wadi, where a small-head dam is constructed, is studied with the help of conformal mappings (the Schwarz-Christoffel formula), based on the analytical solution obtained by Pavlovsky and Polubarinova-Kochina and used by geotechnical engineers a century ago. The dam has a vertical cutoff wall, placed at an arbitrary distance between the head water and tail water. Evaporation from the reservoir horizontal water surface, as well as pumping (or injection) are contributing to the mass balance in the head water that gives an ODE, which is solved either exactly (if no pumping or a constant rate pumping that gives either a linear or the Abel ODE, solved two centuries ago), or by a computer algebra routine of numerical integration (as the Cauchy problem with a “time front” indicating a complete depletion of the reservoir). Optimal design problems are solved, aiming at smart management of the reservoir water storage that will help mini-farms or camel-goat ranches on the wadi banks, which use the reservoir water for irrigation and animal husbandry in harsh environment of Arabia (or other hyperarid regions). Numerical modeling in HYDRUS2D demonstrates kinematics of saturated-unsaturated seepage for the cases of unclogged and clogged reservoir bed. Subsurface hydrology of cascades of dams, including groundwater recharged ones, is also discussed from the perspective of smart water retention and flow in ephemeral and intermittent streams of drylands.

在Pavlovsky和Polubarinova-Kochina的解析解的基础上，利用共形映射（Schwarz-Christoffel公式）研究了小水头坝在河道多孔河床中的瞬态二维渗流，该解析解在一个世纪前被岩土工程师使用。大坝有一个垂直的截水墙，在头水和尾水之间任意距离放置。水库水平水面的蒸发，以及抽水（或注入）都有助于水头水的质量平衡，从而得到ODE，这可以精确地解决（如果没有抽水或恒定速率抽水，得到线性或阿伯ODE，两个世纪前解决），或者通过计算机代数程序的数值积分（如“时间前沿”的柯西问题，表明水库完全枯竭）。解决了优化设计问题，旨在对水库蓄水进行智能管理，帮助小型农场或河堤上的骆驼山羊牧场在阿拉伯（或其他极端干旱地区）恶劣的环境中使用水库水进行灌溉和畜牧业。在HYDRUS2D中进行数值模拟，模拟了未堵塞和堵塞储层情况下饱和-非饱和渗流的运动学。本文还从旱地短流和间歇流的智能蓄水和流动的角度，讨论了坝级（包括地下水补给）的地下水文。

{"title":"Ephemeral reservoirs of low-head wadi dams drained by seepage, evaporation and pumping: the Pavlovsky/Polubarinova-Kochina/Abel analytical legacy redux","authors":"Anvar Kacimov","doi":"10.1016/j.advwatres.2025.105201","DOIUrl":"10.1016/j.advwatres.2025.105201","url":null,"abstract":"<div><div>Transient 2-D seepage through a porous bed of a wadi, where a small-head dam is constructed, is studied with the help of conformal mappings (the Schwarz-Christoffel formula), based on the analytical solution obtained by Pavlovsky and Polubarinova-Kochina and used by geotechnical engineers a century ago. The dam has a vertical cutoff wall, placed at an arbitrary distance between the head water and tail water. Evaporation from the reservoir horizontal water surface, as well as pumping (or injection) are contributing to the mass balance in the head water that gives an ODE, which is solved either exactly (if no pumping or a constant rate pumping that gives either a linear or the Abel ODE, solved two centuries ago), or by a computer algebra routine of numerical integration (as the Cauchy problem with a “time front” indicating a complete depletion of the reservoir). Optimal design problems are solved, aiming at smart management of the reservoir water storage that will help mini-farms or camel-goat ranches on the wadi banks, which use the reservoir water for irrigation and animal husbandry in harsh environment of Arabia (or other hyperarid regions). Numerical modeling in HYDRUS2D demonstrates kinematics of saturated-unsaturated seepage for the cases of unclogged and clogged reservoir bed. Subsurface hydrology of cascades of dams, including groundwater recharged ones, is also discussed from the perspective of smart water retention and flow in ephemeral and intermittent streams of drylands.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105201"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145823466","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

Effects of wave forcing on groundwater and solute transport in lakeshore zone 波浪强迫对湖滨带地下水和溶质运移的影响

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

Advances in Water Resources

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

Xuesong Ding , Xiayang Yu , Xinghua Xu , Pei Xin

Waves are common in lakes and produce high-frequency water level fluctuations, which likely affect groundwater flow and solute transport in lakeshore zones. However, the underlying mechanisms remain unclear. This study used laboratory experiments and numerical simulations to investigate wave-driven surface water-groundwater exchange and solute transport. Results show that waves alter the hydraulic gradient to modify groundwater flow and drive water circulation in the lakeshore zone. Waves enhance water exchange between the lake and inland, while suppressing the net discharge of groundwater into the lake. Waves also alter subsurface solute pathways, shifting discharge zones toward the lake and increasing solute travel time. Sensitivity analysis reveals that the circulation zone extent and the lacustrine groundwater discharge rate are sensitive to wave height, wave period and inland groundwater level. Larger waves expand circulation areas and increase lacustrine groundwater discharge, whereas higher inland groundwater levels reduce circulation zones. These findings emphasize the role of waves in regulating hydrological and solute transport processes in lakeshore zones.

波浪在湖泊中很常见，并产生高频水位波动，这可能影响湖岸地区的地下水流动和溶质运输。然而，潜在的机制仍不清楚。本文采用室内实验和数值模拟相结合的方法对波浪驱动地表水-地下水交换和溶质运移进行了研究。结果表明，波浪通过改变水力梯度来改变湖滨带地下水流量，驱动水循环。波浪加强了湖泊与内陆之间的水交换，同时抑制了地下水进入湖泊的净排放量。波浪也改变了地下的溶质路径，使排放区向湖泊移动，增加了溶质的传播时间。敏感性分析表明，环流带范围和湖泊地下水流量对波高、波周期和内陆地下水位敏感。较大的波浪扩大了环流面积，增加了湖泊地下水排放量，而较高的内陆地下水水位则减少了环流面积。这些发现强调了波浪在调节湖滨带水文和溶质运输过程中的作用。

{"title":"Effects of wave forcing on groundwater and solute transport in lakeshore zone","authors":"Xuesong Ding , Xiayang Yu , Xinghua Xu , Pei Xin","doi":"10.1016/j.advwatres.2025.105199","DOIUrl":"10.1016/j.advwatres.2025.105199","url":null,"abstract":"<div><div>Waves are common in lakes and produce high-frequency water level fluctuations, which likely affect groundwater flow and solute transport in lakeshore zones. However, the underlying mechanisms remain unclear. This study used laboratory experiments and numerical simulations to investigate wave-driven surface water-groundwater exchange and solute transport. Results show that waves alter the hydraulic gradient to modify groundwater flow and drive water circulation in the lakeshore zone. Waves enhance water exchange between the lake and inland, while suppressing the net discharge of groundwater into the lake. Waves also alter subsurface solute pathways, shifting discharge zones toward the lake and increasing solute travel time. Sensitivity analysis reveals that the circulation zone extent and the lacustrine groundwater discharge rate are sensitive to wave height, wave period and inland groundwater level. Larger waves expand circulation areas and increase lacustrine groundwater discharge, whereas higher inland groundwater levels reduce circulation zones. These findings emphasize the role of waves in regulating hydrological and solute transport processes in lakeshore zones.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105199"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796276","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

Handling fracture intersections in the enriched-embedded discrete fracture model (nEDFM) 富嵌入离散裂缝模型（nEDFM）中裂缝相交的处理

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

Advances in Water Resources

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

Kaituo Jiao , Meng Yu , Yongwei Li , Dongxu Han , Yujie Chen , Bofeng Bai , Bo Yu

The embedded discrete fracture model (EDFM) with non-conforming grids has been widely used for simulating flow in fractured porous media; however, it suffers from limited accuracy in computing matrix-fracture mass exchange and fails to handle naturally cemented (blocking) fractures. To address these limitations, the enriched-embedded discrete fracture model (nEDFM) is proposed. The core idea of nEDFM lies in defining the local shape function that approximates the pressure distribution within the interaction region (IR). This function can capture discontinuities in pressure and its gradient across both conductive and blocking fractures and is applicable to media with isotropic or anisotropic matrix permeability. Nonetheless, the original local shape function is limited to cases with at most one fracture in the IR and cannot represent the complex discontinuous pressure fields arising from multiple intersecting fractures. This paper extends nEDFM to fracture intersection scenarios through two key innovations: 1) based on two modification criteria for local shape functions under multiple-fracture configurations, enriched degrees of freedom are reallocated to intersecting fracture cells, and corresponding discontinuous functions are designed to construct new local shape functions that can handle arbitrary intersections of conductive/blocking fractures; 2) to address the inaccuracy of the star-delta transformation in computing flow interaction between fractures with large permeability contrast, the permeability of the intersection area between conductive/blocking fractures is incorporated to establish an improved star-delta formulation. Four numerical cases involving intersecting or closely spaced fractures are conducted to evaluate the relative error and computational cost of nEDFM in comparison with EDFM, projected EDFM, and other non-conforming grid-based methods. The results show that nEDFM can accurately capture the pressure distribution within both conductive and blocking fractures, and effectively represent discontinuities when these fractures intersect or are in close proximity. Overall, nEDFM achieves an excellent balance among accuracy, computational efficiency, and general applicability, demonstrating strong potential for further development and practical implementation.

非均匀网格嵌入离散裂缝模型（EDFM）已广泛应用于裂缝性多孔介质的渗流模拟。然而，它在计算基质-裂缝质量交换方面的精度有限，并且无法处理自然胶结（阻塞）裂缝。针对这些局限性，提出了富嵌入离散裂缝模型（nEDFM）。nEDFM的核心思想在于定义近似于相互作用区域（IR）内压力分布的局部形状函数。该函数可以捕捉传导裂缝和阻断裂缝的压力不连续及其梯度，适用于各向同性或各向异性基质渗透率的介质。然而，原始的局部形状函数仅限于IR中最多有一条裂缝的情况，不能代表多个相交裂缝产生的复杂不连续压力场。本文通过两个关键创新将nEDFM扩展到裂缝相交场景：1)基于多裂缝构型下的两种局部形状函数修正准则，将丰富自由度重新分配到相交的裂缝单元中，并设计相应的不连续函数，构建新的局部形状函数，可以处理导电/阻塞裂缝的任意相交；2)为解决星δ变换在计算渗透率对比大的裂缝间流动相互作用时的不准确性问题，引入导/堵裂缝相交区域渗透率，建立改进的星δ公式。通过4个涉及相交裂缝或紧密裂缝的数值案例，对比EDFM、投影EDFM和其他非一致性网格方法，评估了nEDFM的相对误差和计算成本。结果表明，nEDFM可以准确地捕获导电裂缝和封堵裂缝内的压力分布，并有效地表示裂缝相交或靠近时的不连续面。总的来说，nEDFM在准确性、计算效率和一般适用性之间取得了很好的平衡，显示出进一步发展和实际实施的强大潜力。

{"title":"Handling fracture intersections in the enriched-embedded discrete fracture model (nEDFM)","authors":"Kaituo Jiao , Meng Yu , Yongwei Li , Dongxu Han , Yujie Chen , Bofeng Bai , Bo Yu","doi":"10.1016/j.advwatres.2025.105198","DOIUrl":"10.1016/j.advwatres.2025.105198","url":null,"abstract":"<div><div>The embedded discrete fracture model (EDFM) with non-conforming grids has been widely used for simulating flow in fractured porous media; however, it suffers from limited accuracy in computing matrix-fracture mass exchange and fails to handle naturally cemented (blocking) fractures. To address these limitations, the enriched-embedded discrete fracture model (nEDFM) is proposed. The core idea of nEDFM lies in defining the local shape function that approximates the pressure distribution within the interaction region (IR). This function can capture discontinuities in pressure and its gradient across both conductive and blocking fractures and is applicable to media with isotropic or anisotropic matrix permeability. Nonetheless, the original local shape function is limited to cases with at most one fracture in the IR and cannot represent the complex discontinuous pressure fields arising from multiple intersecting fractures. This paper extends nEDFM to fracture intersection scenarios through two key innovations: 1) based on two modification criteria for local shape functions under multiple-fracture configurations, enriched degrees of freedom are reallocated to intersecting fracture cells, and corresponding discontinuous functions are designed to construct new local shape functions that can handle arbitrary intersections of conductive/blocking fractures; 2) to address the inaccuracy of the star-delta transformation in computing flow interaction between fractures with large permeability contrast, the permeability of the intersection area between conductive/blocking fractures is incorporated to establish an improved star-delta formulation. Four numerical cases involving intersecting or closely spaced fractures are conducted to evaluate the relative error and computational cost of nEDFM in comparison with EDFM, projected EDFM, and other non-conforming grid-based methods. The results show that nEDFM can accurately capture the pressure distribution within both conductive and blocking fractures, and effectively represent discontinuities when these fractures intersect or are in close proximity. Overall, nEDFM achieves an excellent balance among accuracy, computational efficiency, and general applicability, demonstrating strong potential for further development and practical implementation.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"207 ","pages":"Article 105198"},"PeriodicalIF":4.2,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785703","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