Zi Li , Zhenlei Yang , Sergio Galindo-Torres , Ling Li
{"title":"Lattice Boltzmann modelling of capillarity, adsorption and fluid retention in simple geometries: Do capillary and film water have equal matric suction or not?","authors":"Zi Li , Zhenlei Yang , Sergio Galindo-Torres , Ling Li","doi":"10.1016/j.advwatres.2025.104950","DOIUrl":"10.1016/j.advwatres.2025.104950","url":null,"abstract":"<div><div>The pore water retained in unsaturated soil includes film water attached on the solid surface and capillary water in corners or pores, which are mainly controlled by adsorptive force from the solid surface and capillary force from the water-gas interface, respectively. The soil water retention (SWR) curve represents the fundamental characteristic of unsaturated soil, in which the connected capillary and adsorptive water are conventionally presumed to have equivalent suction values as matric suction. Here, a long-range adsorptive fluid-solid interaction force is developed in the mesoscopic multiphase lattice Boltzmann model (LBM) framework to model the macroscopic processes of capillarity and adsorption. The pressure tensors of capillary and film water derived based on mechanical equilibrium and the results of numerical simulations combine to show that the adsorptive suction is much higher than the capillary suction, not following the classical relationship. We attribute this inequality to the different adsorptive interaction potentials incorporated in the capillary and film water pressures, due to the fluid density profiles varying differently with the separation distance from solid surface, and, from the perspective of thermodynamic equilibrium, the deviation of film water and capillary water densities from free water density. The film thickness almost does not change for the given radii of meniscus curvature in simple geometries (i.e., slits and corners). The adsorption effects on the matric suction upscaled from the intrinsic phase average method and on the equivalent pore size distribution are investigated for both single-sized slits and complex pore networks. The findings reveal the influences of capillarity and adsorption on the shape of SWR curve, and help establish the SWR function with accurate physical meanings in the field of soil physics and hydrology and measure the disjoining pressure isotherm properly in colloid and interface chemistry.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104950"},"PeriodicalIF":4.0,"publicationDate":"2025-03-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143620567","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":"Impact of saturation on continuum-scale conductivity and tracer dispersion in heterogeneous porous media","authors":"Doron Kalisman, Brian Berkowitz","doi":"10.1016/j.advwatres.2025.104953","DOIUrl":"10.1016/j.advwatres.2025.104953","url":null,"abstract":"<div><div>This study investigates the interplay among water saturation, hydraulic conductivity, and mechanical dispersion in heterogeneous porous media at the continuum scale. Mechanical dispersion of dissolved chemical tracers is influenced directly by water velocity variability, which is governed by the porous structure, the distribution of the water phase within it, and its corresponding conductivity field. Previous studies have either examined the relationships between these factors in fully saturated conditions, or in partial saturation but without considering continuum scale heterogeneity of the media. Through numerical simulations, the analysis here demonstrates how variations in saturation affect the hydraulic conductivity field and, consequently, mechanical dispersion. The study reveals that longitudinal spreading of the tracer plume, when scaled for varying transport times and velocities, shows a non-monotonic relationship with saturation, being least pronounced at an intermediate degree of saturation. These insights contribute to a more nuanced understanding of tracer transport in partially saturated, heterogeneous media, with implications for environmental and engineering applications.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104953"},"PeriodicalIF":4.0,"publicationDate":"2025-03-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143679608","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":"Derivation of weakly hydrodynamic models in the Dupuit–Forchheimer regime","authors":"Martin Parisot","doi":"10.1016/j.advwatres.2025.104933","DOIUrl":"10.1016/j.advwatres.2025.104933","url":null,"abstract":"<div><div>The current study is dedicated to the formal derivation of a hierarchic of asymptotic models that approximate the groundwater waves problem within the Dupuit–Forchheimer regime, over a regular, non-planar substratum. The derivation methodology employed bears resemblance to the techniques utilized in hierarchic of asymptotic models for approximating the water waves problem in the shallow water regime. Mathematically speaking, the asymptotic models manifest as nonlinear, non-local diffusion equations. We identify an energy dissipation law inherent to these models, thereby bolstering the physical validity and confidence in the proposed framework. A numerical strategy is proposed that preserved at the discrete level the energy dissipation.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104933"},"PeriodicalIF":4.0,"publicationDate":"2025-03-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143611035","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}
Alban Gilletta , Cyrille Bonamy , Marie Robert , Julien Chauchat
{"title":"Assessment of hybrid RANS/LES models for the prediction of the flow and scour around a wall-mounted cylinder","authors":"Alban Gilletta , Cyrille Bonamy , Marie Robert , Julien Chauchat","doi":"10.1016/j.advwatres.2025.104934","DOIUrl":"10.1016/j.advwatres.2025.104934","url":null,"abstract":"<div><div>In many environmental flow situations, a solid body emerges from a sediment bed. This may occur in natural systems, for example when a tree-root emerges from the bed, or around a man-made structure such as a bridge pier or a wind turbine foundation. When this situation occurs, various flow-structure interactions such as the horseshoe vortex or vortex shedding lead to the scour process, namely the erosion pattern around the solid obstacle. In this context, accurate numerical prediction of the flow-structure interactions is of primary importance. In this contribution, three-dimensional hydrodynamic simulations are performed using the single-phase incompressible flow solver pimpleFOAM (OpenFOAMv2206plus). Two hybrid RANS-LES turbulence models available in OpenFOAM are tested: the non-zonal <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-SAS model and the zonal <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-IDDES model. As a first benchmark, a plane channel flow configuration at large bulk Reynolds number (<span><math><mrow><mi>R</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>b</mi></mrow></msub><mo>=</mo></mrow></math></span>1.9 <span><math><mo>×</mo></math></span> 10<sup>4</sup>) is investigated to evaluate the mesh requirements and to determine a set of acceptable numerical parameters. These simulations show that while the <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-IDDES model allows to simulate some features of the fluctuating motions with fewer grid requirements than LES, the <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-SAS model behaves as a RANS model. Concerning the wall-mounted cylinder flow at high Reynolds number (<span><math><mrow><mi>R</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>b</mi></mrow></msub><mo>=</mo></mrow></math></span>1.8 <span><math><mo>×</mo></math></span> 10<sup>5</sup>), the conclusions are quite different, the <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-SAS model switches to LES in the region of interest namely the horseshoe-vortex region and the wake thanks to its non-zonal property. The <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-SAS model provides better predictions than the <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>SST-IDDES model in terms of amplitude of the bimodal oscillation and vortex merging processes in the horseshoe vortex at lower CPU cost than LES. After this benchmarking and assessment of the hybrid RANS-LES model, a new <span><math><mrow><mi>k</mi><mo>−</mo><mi>ω</mi></mrow></math></span>2006-SAS turbulence model has been developed for the two-phase flow solver <em>sedFOAM</em>. The model is then applied to a scour configuration around a wall-mounted cylinder and successfully represents the upstream scour depth. The simulations improve the resolution of the complex horseshoe-vortex dynamics which is directly in","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104934"},"PeriodicalIF":4.0,"publicationDate":"2025-03-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143643459","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":"Transfer learning for geological carbon storage forecasting using neural operator","authors":"Andres Nunez , Siddharth Misra , Yusuf Falola","doi":"10.1016/j.advwatres.2025.104948","DOIUrl":"10.1016/j.advwatres.2025.104948","url":null,"abstract":"<div><div>Geological carbon storage (GCS) is critical for sequestering CO<sub>2</sub> deep underground. GCS projects may face environmental challenges, such as leakage risks, adverse pressure buildup, and groundwater contamination. Numerical simulators play a vital role in accurate forecasting but can be computationally expensive. In this work, we leveraged an updated Fourier Neural Operator (FNO) which includes data sparsity management, to learn to rapidly forecast pressure and CO<sub>2</sub> phase saturation distributions in a geological carbon storage (GCS) reservoir. Compared to commercial reservoir simulators, FNO-based forecasting offers accurate prediction while reducing the computational time by a factor of 40, enabling high volume of forecasting in less time. Additionally, we applied transfer learning (TL) to further reduce the data and computational requirements of the FNO-based forecasting across a wide array of scenarios. Specifically, we demonstrated the usefulness of TL in accurately predicting the pressure and CO<sub>2</sub> saturation distributions for uncertain and variable geological and operational conditions. The results of this study indicate that the improved FNO workflow reduces the computational time by approximately 97 %, and the relative mean error for predicting both CO<sub>2</sub> saturation and pressure distributions is <1 %. Generally, the use of TL effectively transfers knowledge from a pre-existing model to other related tasks. TL significantly reduces the required training data by 78 % while maintaining a relative mean error below 5 %. Although, the results in this work can be further improved, this study demonstrates the potential of integrating FNO and TL to reduce computational time and data requirements for CO<sub>2</sub> forecasts during GCS projects, providing a more efficient and faster approach.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104948"},"PeriodicalIF":4.0,"publicationDate":"2025-03-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143620566","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}
Brandon Hilliard , William J. Reeder , Ralph Budwig , Vibhav Durgesh , Bishal Bhattarai , Benjamin T. Martin , Tao Xing , Daniele Tonina
{"title":"Unveiling surface-subsurface flow interactions of a salmon redd","authors":"Brandon Hilliard , William J. Reeder , Ralph Budwig , Vibhav Durgesh , Bishal Bhattarai , Benjamin T. Martin , Tao Xing , Daniele Tonina","doi":"10.1016/j.advwatres.2025.104947","DOIUrl":"10.1016/j.advwatres.2025.104947","url":null,"abstract":"<div><div>Female salmonids bury their eggs in streambed gravel by digging a pit where they lay their eggs, which they then cover with gravel from a second pit, forming a rough-surfaced dune-like structure called a redd. The interaction between a redd and the stream flow induces surface water to flow into the sediment, through egg pockets, and reemerge downstream of the crest. These downwelling and upwelling flows form the hyporheic exchange, which is vital for the embryos’ development because it regulates the egg pocket temperature regime and delivers oxygen-rich surface water to the embryos. Here, we experimentally investigate the effects of (1) redd surface roughness and (2) egg pocket presence on redd-induced hyporheic flows under different surface hydraulics by constructing a synthetic redd made of transparent grains that allow us to measure freestream, near-bed, and interstitial flow velocities with non-invasive techniques in a recirculating flume. Results indicate that flow through redds is proportional to the squared Froude number, but surface roughness causes this relation to be more complex. The egg pocket, having larger grains than the surrounding matrix, enhances mechanical dispersion, which increases water mixing within the egg pocket and causes convergence of downwelling flows towards the egg pocket. Overall, these results show that river modeling studies focused on salmonid preservation should consider surface roughness as well as varying hydraulic conductivities within the nesting environments.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104947"},"PeriodicalIF":4.0,"publicationDate":"2025-03-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143620568","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}
Muhammad Nasir, Shintaro Matsushita, Kailin Wang, Masayuki Osada, Shu Yamashita, Wilson Susanto, Sotheavuth Sin, Tetsuya Suekane
{"title":"Pore-scale relative permeability and saturation analysis under wide-ranging injection velocity and wettability during primary CO2 injection for geological carbon sequestration","authors":"Muhammad Nasir, Shintaro Matsushita, Kailin Wang, Masayuki Osada, Shu Yamashita, Wilson Susanto, Sotheavuth Sin, Tetsuya Suekane","doi":"10.1016/j.advwatres.2025.104938","DOIUrl":"10.1016/j.advwatres.2025.104938","url":null,"abstract":"<div><div>We performed two-dimensional (2D) pore-scale simulations of primary CO<sub>2</sub> injection using a weakly compressible scheme for geological carbon sequestration (GCS) applications. The aim was to analyze pore-scale relative permeability and saturation of CO<sub>2</sub> under wide-ranging injection velocities and wettabilities. The results show that saturation is highest for viscous fingering, lowest for crossover (− 5.82 < <em>logCa</em> < − 4.86; θ < 60°), and remains high in the capillary fingering regime even though the relative permeability of CO<sub>2</sub> is minimum. This trend occurs because saturation is influenced not only by the value of relative permeability but also by the frequency of relative permeability fluctuations. At a low injection velocity and contact angle, frequent permeability fluctuations due to Haines jumps result in high saturation despite the low relative permeability. At intermediate injection velocity and low contact angle, both the relative permeability and its fluctuations are moderate, leading to lower CO<sub>2</sub> saturation. The present work bridges the understanding of displacement-front advancement at the pore-network scale with relative permeability, which links the pore-scale meniscus dynamics with the large-scale Darcy-flow parameters. As the CO<sub>2</sub> flows away from the injection site in large-scale GCS applications, the displacement pattern exhibits crossover regime, resulting in minimal displacement efficiency. In a strongly wetting porous medium, this condition is severe because crossover regime spans a wide range of capillary numbers.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104938"},"PeriodicalIF":4.0,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143549864","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}
Anthony D. Miller , Amir Jazayeri , S. Cristina Solórzano-Rivas , Adrian D. Werner
{"title":"A Dupuit-Forchheimer solution to the extraction of seawater from coastal aquifers","authors":"Anthony D. Miller , Amir Jazayeri , S. Cristina Solórzano-Rivas , Adrian D. Werner","doi":"10.1016/j.advwatres.2025.104937","DOIUrl":"10.1016/j.advwatres.2025.104937","url":null,"abstract":"<div><div>Seawater intrusion can be mitigated by extracting saltwater, creating a negative hydraulic barrier that reduces the extent of saltwater in coastal aquifers. The effects of seawater extraction are analyzed in the current study through a semi-analytical methodology based on sharp-interface, steady-state conditions. The methodology is based on the Dupuit-Forchheimer approximation and applies a power series approach to obtain an exact solution to the seawater extent in the aquifer. Alternatively, the solution can be obtained by a Runge-Kutta method, thereby allowing for rapid assessment of the efficacy of seawater extraction for simple situations (e.g., uniform, homogeneous aquifer, and a continuous line sink well). Comparisons with numerical simulations using SEAWAT (including cases with nominally zero dispersion and with dispersion) display good agreement between the sharp-interface solution using a modified density and the 50 % seawater contour from SEAWAT. The results show that for a given well location there is an optimal (maximum) extraction rate of seawater that minimises the landward extent of seawater. Equivalently, for a given extraction rate, there is an optimal well location. These optima place a well at the seawater-freshwater interface, so in practice will likely lead to the partial extraction of freshwater. Nonetheless, they represent the limits defining the operating region for a negative hydraulic barrier in non-dimensional parameter space and thereby can inform initial decisions regarding applicability. The method presented provides a rapid assessment tool for examining the interplay between the extraction well location and the extraction rate.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104937"},"PeriodicalIF":4.0,"publicationDate":"2025-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143593686","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}
C. Nguyen, Y. Yang, G. Loi, T. Russell, P. Bedrikovetsky
{"title":"Salinity effects on fines migration in aquifers: Stochastic model and its upscaling","authors":"C. Nguyen, Y. Yang, G. Loi, T. Russell, P. Bedrikovetsky","doi":"10.1016/j.advwatres.2025.104932","DOIUrl":"10.1016/j.advwatres.2025.104932","url":null,"abstract":"<div><div>Colloidal-suspension-nano flows with varying ionic strength are widely present in nature and industry. The variation of brine salinity, which highly affects electrostatic particle-rock interaction, triggers fines detachment and consequent rock alteration. The microscale models for fines detachment at the pore-particle and at rock-reservoir scales are widely used to predict core and field behaviour under fines migration, while the relationship between those models hasn't been established. This includes the lack of upscaling and downscaling procedures, which prevents the determination of micro-scale parameters from lab corefloods or well production histories and calculation of large-scale model functions from SEM and microfluid tests. In this work for the first time, we derive the rock-scale detachment model as expressed as a function of maximum attached concentration versus brine salinity (maximum retention function MRF) from particle-scale torque balance of attaching and detaching forces. Reflecting micro heterogeneity of the pore space and attached particles, we consider the mutual probabilistic distributions of geometric and electrostatic coefficients from the torque balance of detaching drag and attaching electrostatic DLVO forces exerting the particle. This determines the cumulative distribution of critical salinity, which is calculated from the torque balance and the mutual distribution of microscale parameters and defines the salinity-dependent MRF. This upscaling procedure is performed by the Monte Carlo algorithm for MRF calculation. The corresponding downscaling comprises tuning the mean and variance values for some micro-scale parameters from the MRF. These algorithms are used to treat three coreflood data sets with varying salinity, determine the MRF, and calculate mean values of lever arm and aspect ratios; the match is high, and the obtained microscale parameters are within the common intervals. The upscaling technique developed allows for sensitivity analysis of detachment with respect to microscale parameters and velocity. We also developed the upscaling procedure for MRF recalculation from dependency of one flow parameter to another. The velocity-dependent MRF was recalculated from three coreflood-based salinity-dependent MRFs, yielding lab-based prediction of well behaviour for water injectors and producers.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104932"},"PeriodicalIF":4.0,"publicationDate":"2025-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143620565","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}
Shakhawat Hossain , Gary J. Hampson , Carl Jacquemyn , Matthew D. Jackson , Dmytro Petrovskyy , Sebastian Geiger , Julio D. Machado Silva , Sicilia Judice , Fazilatur Rahman , M. Costa Sousa
{"title":"Effective permeability of fluvial lithofacies in the Bunter Sandstone Formation, UK","authors":"Shakhawat Hossain , Gary J. Hampson , Carl Jacquemyn , Matthew D. Jackson , Dmytro Petrovskyy , Sebastian Geiger , Julio D. Machado Silva , Sicilia Judice , Fazilatur Rahman , M. Costa Sousa","doi":"10.1016/j.advwatres.2025.104936","DOIUrl":"10.1016/j.advwatres.2025.104936","url":null,"abstract":"<div><div>Understanding effective permeability is crucial for predicting fluid migration and trapping in subsurface reservoirs. The Bunter Sandstone of northwestern Europe hosts major groundwater and geothermal resources and is targeted for CO<sub>2</sub> storage projects. Here the effective permeability of fluvial facies within the Bunter Sandstone Formation was assessed using facies-scale models. Twelve lithofacies were modeled based on core and outcrop observations of their geometries and dimensions. Permeability values from minipermeameter measurements were assigned to low- and high-permeability lithologies in each facies. The dimensions of a Representative Elementary Volume (REV) in depositional dip, depositional strike and vertical directions were determined by extracting sub-volumes from the models at different scales, calculating values of effective permeability for each sub-volume, and identifying the sub-volume at which the values of effective permeability stabilise as the REV. The REV dimensions vary with facies type and flow direction, but are typically of order tens of centimetres to metres in size, significantly larger than a typical core plug. Having identified the REV, we analyze the effective permeabilities of the different facies types. Normalized values of effective permeabilities in depositional dip, strike and vertical directions (<em>k<sub>d</sub>, k<sub>s</sub>, k<sub>v</sub></em>), relative to the permeability of low- and high-permeability lithologies in each facies, display a positive linear correlation with the proportion of high-permeability lithology (clay-poor sandstone) for all facies. Therefore, the proportion of clay-poor sandstone, as measured in core data, can be used to predict facies-scale effective permeability in the Bunter Sandstone Formation, as well as in analogous fluvial deposits globally.</div></div>","PeriodicalId":7614,"journal":{"name":"Advances in Water Resources","volume":"199 ","pages":"Article 104936"},"PeriodicalIF":4.0,"publicationDate":"2025-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143549861","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}