Advances in Water Resources最新文献

{"title":"Boosting the reconstruction performance of 3D Multi-porous media using double generative adversarial networks","authors":"Xiaoxiang Yin,&nbsp;Mingliang Gao,&nbsp;Ai Luo,&nbsp;Geling Xu","doi":"10.1016/j.advwatres.2024.104843","DOIUrl":"10.1016/j.advwatres.2024.104843","url":null,"abstract":"<div><div>With the continuous improvement of mathematical modeling technology, reconstructing the three-dimensional structure of media from two-dimensional reference images has become an important research method for the three-dimensional modeling of multi-porous media. Deep-learning-based methods are currently popular and form the focus of this research field. However, the performance of deep learning in reconstructing the three-dimensional structure of media from two-dimensional reference images still requires improvement. To enhance the diversity and generalization of network-generated three-dimensional models, this study proposed a preprocessing method that correlated two-dimensional reference images with Gaussian noise, a three-orthogonal random section constraint method, and a dual generative adversarial network (DGAN)-based model. Multiple sets of core samples, a set of building materials, and a set of battery-material samples were used to verify the performance of the proposed network. Both intuitive morphological and statistical feature comparisons showed that the DGAN model solved the problem of insufficient diversity and generalization when reconstructing three-dimensional porous media from a single image using deep-learning-based methods. The morphological and statistical features of the reconstructed three-dimensional structure also exhibited good consistency with the reference two-dimensional image, and the training efficiency of the network was greatly improved.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"194 ","pages":"Article 104843"},"PeriodicalIF":4.0,"publicationDate":"2024-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142573113","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An extended two-parameter mixed-dimensional model of fractured porous media incorporating entrance flow and boundary-layer transition effects","authors":"Marcio A. Murad ,&nbsp;Luciane A. Schuh ,&nbsp;Igor Mozolevski ,&nbsp;Josue Barroso","doi":"10.1016/j.advwatres.2024.104838","DOIUrl":"10.1016/j.advwatres.2024.104838","url":null,"abstract":"<div><div>We develop an enhanced reduced model for single-phase flow in fractured porous media capable of incorporating more realistic interface conditions at the fracture terminations. In addition to the traditional dimensional model reduction, where the elements of the discrete fracture network are treated as lower dimensional manifolds embedded in the porous matrix, we explore the microscale behavior of the boundary layer flow at the entrances of a fracture bounded by two parallel plates to construct a new set of interface conditions of Robin-type, giving rise to localized pressure jumps at the fracture edges. Within this enriched description, sharper reduced flow and tracer transport mixed-dimensional models are constructed in the asymptotic limit ruled by two small parameters related to the ratio between fracture aperture and entrance developing length and a macroscopic length scale. The discrete flow/transport mixed-dimensional model is discretized by a new discontinuous Galerkin(dG)-based formulation. An adequate version of the Galerkin-Newton method is developed for the numerical treatment of the non-linear Robin interface condition. Considering several fracture arrangements, numerical results illustrate the sharper description of the model proposed herein in predicting flow and tracer transport patterns in fractured media.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104838"},"PeriodicalIF":4.0,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535560","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media","authors":"Aleksei G. Sorokin ,&nbsp;Aleksandra Pachalieva ,&nbsp;Daniel O’Malley ,&nbsp;James M. Hyman ,&nbsp;Fred J. Hickernell ,&nbsp;Nicolas W. Hengartner","doi":"10.1016/j.advwatres.2024.104836","DOIUrl":"10.1016/j.advwatres.2024.104836","url":null,"abstract":"<div><div>Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution.</div><div>Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as <span><math><mrow><mi>O</mi><mrow><mo>(</mo><mi>n</mi><mo>log</mo><mi>n</mi><mo>)</mo></mrow></mrow></math></span> instead of the typical <span><math><mrow><mi>O</mi><mrow><mo>(</mo><msup><mrow><mi>n</mi></mrow><mrow><mn>3</mn></mrow></msup><mo>)</mo></mrow></mrow></math></span> rate in the number of samples <span><math><mi>n</mi></math></span> used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104836"},"PeriodicalIF":4.0,"publicationDate":"2024-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535559","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A physics-constraint neural network for CO2 storage in deep saline aquifers during injection and post-injection periods","authors":"Mengjie Zhao ,&nbsp;Yuhang Wang ,&nbsp;Marc Gerritsma ,&nbsp;Hadi Hajibeygi","doi":"10.1016/j.advwatres.2024.104837","DOIUrl":"10.1016/j.advwatres.2024.104837","url":null,"abstract":"<div><div>CO₂ capture and storage is a viable solution in the effort to mitigate global climate change. Deep saline aquifers, in particular, have emerged as promising storage options, owing to their vast capacity and widespread distribution. However, the task of proficiently monitoring and simulating CO₂ behavior within these formations poses significant challenges. To address this, we introduce the physics-constraint neural network for CO₂ storage (CO₂PCNet), a model specifically designed for simulating and monitoring CO₂ storage in deep saline aquifers during injection and post-injection periods. Recognizing the significant challenges in accurately modeling the distribution and movement of CO₂ under varying permeability conditions, the CO₂PCNet integrates the principles of physics with the robustness of deep learning, serving as a powerful surrogate model. The architecture of CO₂PCNet starts with an encoder that adeptly processes spatial features from overall mole fraction (<span><math><msub><mrow><mi>z</mi></mrow><mrow><msub><mrow><mtext>CO</mtext></mrow><mrow><mn>2</mn></mrow></msub></mrow></msub></math></span>) and pressure fields (<span><math><msub><mrow><mi>P</mi></mrow><mrow><mi>l</mi></mrow></msub></math></span>), capturing the complex dynamics of a CO₂ trajectory. By incorporating permeability information through a conditioning step, the network ensures a faithful representation of the influences on CO₂ behavior in subsurface conditions. A ConvLSTM module subsequently discerns temporal evolutions, reflecting the real-world progression of CO₂ plumes within the reservoir. Lastly, the decoder precisely reconstructs the predictive spatial profile of CO₂ distribution. CO₂PCNet, with its integration of convolutional layers, recurrent mechanisms, and physics-informed constraints, offers a refined approach to CO₂ storage simulation. This model offers the potential of utilizing advanced computational methods in advancing CCS practices.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104837"},"PeriodicalIF":4.0,"publicationDate":"2024-10-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535558","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multi-temporal drought rarity curves—A yearly classification of meteorological drought severity in France","authors":"Juliette Blanchet ,&nbsp;Baptiste Ainési ,&nbsp;Sandra Rome ,&nbsp;Jean-Dominique Creutin","doi":"10.1016/j.advwatres.2024.104829","DOIUrl":"10.1016/j.advwatres.2024.104829","url":null,"abstract":"<div><div>Droughts are recurrent phenomena that present a large variety of space and time patterns making rather difficult the assessment of their rarity and the comparison between events.</div><div>Our study focuses on the “memory effect” of meteorological drought over France using the daily raingauge network maintained by Météo-France over 1950–2022. The proposed easy tool of rarity curves analyses how drought events build and persist across time. The approach is purely statistic, assuming that drought consequences depend on the probability of non exceedance (“rarity”) of antecedent rainfall accumulations. In order to cover a large spectrum of “memory effects”, we consider a continuum of accumulation periods ranging from a few weeks to several years. The rarity curve of a given year displays the most severe rarity values encountered during the year as a function of the various accumulation periods (256 values ranging from 4 to 260 weeks in our case).</div><div>Over the study period of 1950–2022 we show how the shape of rarity curves discriminate short- and long-term historical droughts. A k-means algorithm proves to be useful to classify the different years into typical drought situations from continuous drying to continuous wetting over the antecedent five years. Looking at the chronology of the found clusters also allows studying how droughts build and persist across time.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104829"},"PeriodicalIF":4.0,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535561","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Laboratory and numerical investigations on land-sourced solute transport in coastal fractured aquifers","authors":"Yipeng Zeng ,&nbsp;Yifan Xie ,&nbsp;Yu Ye ,&nbsp;Chengji Shen ,&nbsp;Tongchao Nan ,&nbsp;Chunhui Lu","doi":"10.1016/j.advwatres.2024.104839","DOIUrl":"10.1016/j.advwatres.2024.104839","url":null,"abstract":"<div><div>Land-based pollutants threaten coastal aquifers, highlighting the need to protect groundwater and nearshore marine ecosystems. While aquifer heterogeneity has been recognized as a significant factor affecting solute behavior, the impact of fractures on land-sourced solute transport in coastal aquifers remains unclear. This study attempted to address this issue through laboratory experiments and discrete fracture matrix (DFM) models. The impact of horizontal fractures on the temporal and spatial characteristics of solute transport, spreading, and discharge under seawater intrusion was analyzed based on variations in fracture position and length. The results show that fractures/low-velocity zones (LVZ) can accelerate/delay solutes, dividing them into different transport modes and enhancing/prolonging their spreading/discharge duration. Changes in fracture position and length also affect its transport acceleration and path deviation abilities, which ultimately determine when solute discharge occurs. The mixing zone and unsaturated zone, in addition to the LVZ, hinder solute transport, reducing the rate and delaying the end of solute discharge. Meanwhile, fractures facilitate solute transport into the saltwater wedge, expanding the solute discharge zone. However, if solutes are initially within the LVZ, their entry into the fracture rely on their distance from the fracture's near-land edge and the size of the fracture's convergence zone.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104839"},"PeriodicalIF":4.0,"publicationDate":"2024-10-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535556","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Dynamic mode decomposition of GRACE satellite data","authors":"G. Libero ,&nbsp;V. Ciriello ,&nbsp;D.M. Tartakovsky","doi":"10.1016/j.advwatres.2024.104834","DOIUrl":"10.1016/j.advwatres.2024.104834","url":null,"abstract":"<div><div>Advancements in satellite technology yield environmental data with ever improving spatial coverage and temporal resolution. This necessitates the development of techniques to discern actionable information from large amounts of such data. We explore the potential of dynamic mode decomposition (DMD) to discover the dynamics of spatially correlated structures present in global-scale data, specifically in observations of total water storage anomalies provided by GRACE satellite missions. Our results demonstrate that DMD enables data compression and extrapolation from a reduced set of dominant spatiotemporal structures. The accuracy of its predictions of global system dynamics is preserved in its reconstruction of local time series. These findings suggest potential uses of DMD in analysis of remote-sensing data for hydrologic applications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104834"},"PeriodicalIF":4.0,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535555","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Semi-implicit schemes for modeling water flow and solute transport in unsaturated soils","authors":"Hamza Kamil ,&nbsp;Abdelaziz Beljadid ,&nbsp;Azzeddine Soulaïmani ,&nbsp;Yves Bourgault","doi":"10.1016/j.advwatres.2024.104835","DOIUrl":"10.1016/j.advwatres.2024.104835","url":null,"abstract":"<div><div>The coupled model of water flow and solute transport in unsaturated soils is addressed in this study. Building upon previous research findings by Keita, Beljadid, and Bourgault, we investigate a class of second-order time-stepping techniques where two free parameters are introduced, to identify the most stable and accurate scheme. The spatial discretization of the Richards equation is accomplished using the mixed finite element method. The proposed approach involves formulating noniterative schemes using an extrapolation formula and Taylor approximation in time to linearize nonlinear terms. Additionally, a specialized regularization technique is applied to ensure the convergence of the proposed numerical methods. Numerical simulations are conducted to determine the optimal scheme for solving the Richards equation, which is subsequently extended to the transport equation.</div><div>Numerical simulations of water flow reveal the good accuracy of three schemes—SBDF, MSBDF, and Richards-M2 for homogeneous and heterogeneous soils. Notably, the SBDF scheme stands out for its computational efficiency and stability, especially when gravity forces dominate over capillary forces. Through numerical analysis of the coupled semi-implicit schemes, our results affirm the SBDF scheme’s superior robustness, establishing it as the optimal choice among the proposed numerical methods. Therefore, the SBDF scheme is employed to solve the coupled model of water flow and solute transport. We conducted various numerical experiments to solve the coupled model, addressing scenarios including single and multispecies nitrogen transport, pore water electrical conductivity, and nitrate transport. The SBDF scheme’s accuracy was rigorously verified through comparisons with reference solutions and experimental data. This establishes the SBDF scheme as an efficient alternative to traditional implicit methods for modeling water flow and solute transport in unsaturated soils.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104835"},"PeriodicalIF":4.0,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535557","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Deep learning-based geological parameterization for history matching CO2 plume migration in complex aquifers","authors":"Li Feng ,&nbsp;Shaoxing Mo ,&nbsp;Alexander Y. Sun ,&nbsp;Dexi Wang ,&nbsp;Zhengmao Yang ,&nbsp;Yuhan Chen ,&nbsp;Haiou Wang ,&nbsp;Jichun Wu ,&nbsp;Xiaoqing Shi","doi":"10.1016/j.advwatres.2024.104833","DOIUrl":"10.1016/j.advwatres.2024.104833","url":null,"abstract":"<div><div>History matching is crucial for reliable numerical simulation of geological carbon storage (GCS) in deep subsurface aquifers. This study focuses on inferring highly complex aquifer permeability fields with multi- and intra-facies heterogeneity to improve the characterization of CO₂ plume migration. We propose a deep learning (DL)-based parameterization strategy combined with the ensemble smoother with multiple data assimilation (ESMDA) algorithm to formulate an integrated inverse framework. The DL model is employed to parameterize non-Gaussian permeability fields using low-dimensional latent variables in a Gaussian distribution, thereby mitigating the non-Gaussianity issue faced by the ensemble-based ESMDA inverse method and simultaneously alleviating the computational burden of high-dimensional inversion. The efficacy of the integrated DL-ESMDA inverse framework is demonstrated using a 3-D GCS model, where it estimates the non-Gaussian permeability field characterized by multi- and intra-facies heterogeneity. Results show that the DL model is able to represent the highly complex and high-dimensional permeability fields using low-dimensional latent vectors. The DL-ESMDA framework sequentially updates these low-dimensional latent vectors instead of the original high-dimensional permeability field to obtain posterior estimations of the permeability field. The resulting CO₂ plume migration closely matches historical measurements, suggesting a significantly improved model reliability after history matching. Additionally, a substantial reduction in uncertainty for future plume migration predictions beyond the history matching period is observed. The proposed framework provides an effective approach for reliable characterization of CO₂ plume migration in highly heterogeneous aquifers, enhancing GCS project operation and risk analysis.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104833"},"PeriodicalIF":4.0,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142535554","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Why do published models for fluvial and estuarine morphodynamics use unrealistic representations of the effects of transverse bed slopes?","authors":"Michael Tritthart ,&nbsp;Davide Vanzo ,&nbsp;Victor Chavarrías ,&nbsp;Annunziato Siviglia ,&nbsp;Kees Sloff ,&nbsp;Erik Mosselman","doi":"10.1016/j.advwatres.2024.104831","DOIUrl":"10.1016/j.advwatres.2024.104831","url":null,"abstract":"<div><div>The sediment transport direction is affected by the bed slope. This effect is of crucial importance for two- and three-dimensional modelling of the interaction between the flow of water and the alluvial bed. It is not uncommon to find applications of numerical morphodynamic models in the literature that exaggerate the effects of transverse bed slopes on sediment transport compared to results from laboratory experiments. We investigate mathematically the consequences of such an approach, and we analyse through numerical simulations different explanations for the need to apply deviating values. The study reveals that the reason often lies in the setup of the numerical models, such as the choice of mesh resolution or the necessity to comply with specific aspects of the numerical scheme. The missing or inadequate implementation of physical processes in the model is another cause. All of these effects can be compensated by artificial diffusion added through the bed slope effect coefficients. Since increased diffusion strongly alters the physical processes of self-formed bed morphology, we recommend that modellers address the root causes of inflated erosion and deposition. Bed slope effect coefficients should be applied within the range found in the original publications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"193 ","pages":"Article 104831"},"PeriodicalIF":4.0,"publicationDate":"2024-10-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142445457","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"环境科学与生态学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}