Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method     

Tao Chen  Christoph Clauser  Gabriele Marquart  Karen Willbrand  Thomas Hiller 《Hydrogeology Journal》2018,26(6):1903-1916

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.  相似文献   

Tracing submarine groundwater discharge flux in Tolo Harbor,Hong Kong (China)     

Yi Liu  Jiu Jimmy Jiao  Ho Kwan Cheng 《Hydrogeology Journal》2018,26(6):1857-1873

Submarine groundwater discharge (SGD) is an important pathway for groundwater and associated chemicals to discharge to the sea. Groundwater levels monitored along a transect perpendicular to the shoreline are used to calculate SGD flux from the nearshore aquifer to Tolo Harbor, Hong Kong (China). The calculated SGD flux—recharge/discharge measured with Darcy’s Law methods—agrees well with estimates based on geo-tracer techniques and seepage meter in Tolo Harbor during previous studies. The estimated freshwater SGD is 1.69–2.0 m²/d at the study site and 0.3?±?0.04 cm/d for the whole of Tolo Harbor, which is comparable to the river discharge (0.25?±?0.07 cm/d) and precipitation (0.45?±?0.15 cm/d). The tide-driven SGD in the intertidal zone is 13.98–17.59 m²/d at the study site and 2.42?±?0.56 cm/d for the whole of Tolo Harbor. The SGD occurring in the subtidal zone and the bottom of Tolo Harbor is 3.12?±?4.63 cm/d. Fresh SGD accounts for ~5% of the total SGD, while the rest (~95%) is contributed by saline SGD driven by various forces. About 96% of the tide-driven SGD in the intertidal zone occurs in the ebbing tide period because the head difference between the groundwater level and sea level is great during this period. Tide-driven SGD in the spring tide is ~1.2 times that during neap tide. The tidal fluctuation amplitude and tide-driven SGD in the intertidal zone are positively correlated to each other; thus, a spring neap variation of the tide-driven SGD is observed.  相似文献   

Characterising the vertical separation of shale-gas source rocks and aquifers across England and Wales (UK)     

Sian E. Loveless  John P. Bloomfield  Robert S. Ward  Alwyn J. Hart  Ian R. Davey  Melinda A. Lewis 《Hydrogeology Journal》2018,26(6):1975-1987

Shale gas is considered by many to have the potential to provide the UK with greater energy security, economic growth and jobs. However, development of a shale gas industry is highly contentious due to environmental concerns including the risk of groundwater pollution. Evidence suggests that the vertical separation between exploited shale units and aquifers is an important factor in the risk to groundwater from shale gas exploitation. A methodology is presented to assess the vertical separation between different pairs of aquifers and shales that are present across England and Wales. The application of the method is then demonstrated for two of these pairs—the Cretaceous Chalk Group aquifer and the Upper Jurassic Kimmeridge Clay Formation, and the Triassic sandstone aquifer and the Carboniferous Bowland Shale Formation. Challenges in defining what might be considered criteria for ‘safe separation’ between a shale gas formation and an overlying aquifer are discussed, in particular with respect to uncertainties in geological properties, aquifer extents and determination of socially acceptable risk levels. Modelled vertical separations suggest that the risk of aquifer contamination from shale exploration will vary greatly between shale–aquifer pairs and between regions and this will need to be considered carefully as part of the risk assessment and management for any shale gas development.  相似文献   

Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer     

Elena Tamayo-Mas  Marco Bianchi  Majdi Mansour 《Hydrogeology Journal》2018,26(6):1917-1933

This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings.  相似文献   

Spatial resolution of transport parameters in a subtropical karst conduit system during dry and wet seasons     

Anna Ender  Nadine Goeppert  Nico Goldscheider 《Hydrogeology Journal》2018,26(7):2241-2255

Karst aquifers are characterized by a high degree of hydrologic variability and spatial heterogeneity of transport parameters. Tracer tests allow the quantification of these parameters, but conventional point-to-point experiments fail to capture spatiotemporal variations of flow and transport. The goal of this study was to elucidate the spatial distribution of transport parameters in a karst conduit system at different flow conditions. Therefore, six tracer tests were conducted in an active and accessible cave system in Vietnam during dry and wet seasons. Injections and monitoring were done at five sites along the flow system: a swallow hole, two sites inside the cave, and two springs draining the system. Breakthrough curves (BTCs) were modeled with CXTFIT software using the one-dimensional advection-dispersion model and the two-region nonequilibrium model. In order to obtain transport parameters in the individual sections of the system, a multi-pulse injection approach was used, which was realized by using the BTCs from one section as input functions for the next section. Major findings include: (1) In the entire system, mean flow velocities increase from 183 to 1,043 m/h with increasing discharge, while (2) the proportion of immobile fluid regions decrease; (3) the lowest dispersivity was found at intermediate discharge; (4) in the individual cave sections, flow velocities decrease along the flow direction, related to decreasing gradients, while (5) dispersivity is highest in the middle section of the cave. The obtained results provide a valuable basis for the development of an adapted water management strategy for a projected water-supply system.  相似文献   

Hybrid Gaussian-cubic radial basis functions for scattered data interpolation     

Pankaj K. Mishra  Sankar K. Nath  Mrinal K. Sen  Gregory E. Fasshauer 《Computational Geosciences》2018,22(5):1203-1218

Scattered data interpolation schemes using kriging and radial basis functions (RBFs) have the advantage of being meshless and dimensional independent; however, for the datasets having insufficient observations, RBFs have the advantage over geostatistical methods as the latter requires variogram study and statistical expertise. Moreover, RBFs can be used for scattered data interpolation with very good convergence, which makes them desirable for shape function interpolation in meshless methods for numerical solution of partial differential equations. For interpolation of large datasets, however, RBFs in their usual form, lead to solving an ill-conditioned system of equations, for which, a small error in the data can cause a significantly large error in the interpolated solution. In order to reduce this limitation, we propose a hybrid kernel by using the conventional Gaussian and a shape parameter independent cubic kernel. Global particle swarm optimization method has been used to analyze the optimal values of the shape parameter as well as the weight coefficients controlling the Gaussian and the cubic part in the hybridization. Through a series of numerical tests, we demonstrate that such hybridization stabilizes the interpolation scheme by yielding a far superior implementation compared to those obtained by using only the Gaussian or cubic kernels. The proposed kernel maintains the accuracy and stability at small shape parameter as well as relatively large degrees of freedom, which exhibit its potential for scattered data interpolation and intrigues its application in global as well as local meshless methods for numerical solution of PDEs.  相似文献   

Convergence analysis of fixed stress split iterative scheme for anisotropic poroelasticity with tensor Biot parameter     

Saumik Dana  Mary F. Wheeler 《Computational Geosciences》2018,22(5):1219-1230

We perform a convergence analysis of the fixed stress split iterative scheme for the Biot system modeling coupled flow and deformation in anisotropic poroelastic media with tensor Biot parameter. The fixed stress split iterative scheme solves the flow subproblem with all components of the stress tensor frozen using a multipoint flux mixed finite element method, followed by the poromechanics subproblem using a conforming Galerkin method in every coupling iteration at each time step. The coupling iterations are repeated until convergence and Backward Euler is employed for time marching. The convergence analysis is based on studying the equations satisfied by the difference of iterates to show that the fixed stress split iterative scheme for anisotropic poroelasticity with Biot tensor is contractive. We also demonstrate that the scheme is numerically convergent using the classical Mandel’s problem solution for transverse isotropy.  相似文献   

Impact of recharge water temperature on bioclogging during managed aquifer recharge: a laboratory study     

Lu Xia  Zongjun Gao  Xilai Zheng  Jiuchuan Wei 《Hydrogeology Journal》2018,26(7):2173-2187

To investigate the effect of recharge water temperature on bioclogging processes and mechanisms during seasonal managed aquifer recharge (MAR), two groups of laboratory percolation experiments were conducted: a winter test and a summer test. The temperatures were controlled at ~5±2 and ~15±3 °C, and the tests involved bacterial inoculums acquired from well water during March 2014 and August 2015, for the winter and summer tests, respectively. The results indicated that the sand columns clogged ~10 times faster in the summer test due to a 10-fold larger bacterial growth rate. The maximum concentrations of total extracellular polymeric substances (EPS) in the winter test were approximately twice those in the summer test, primarily caused by a ~200 μg/g sand increase of both loosely bound EPS (LB-EPS) and tightly bound EPS (TB-EPS). In the first half of the experimental period, the accumulation of bacteria cells and EPS production induced rapid bioclogging in both the winter and summer tests. Afterward, increasing bacterial growth dominated the bioclogging in the summer test, while the accumulation of LB-EPS led to further bioclogging in the winter test. The biological analysis determined that the dominant bacteria in experiments for both seasons were different and the bacterial community diversity was ~50% higher in the winter test than that for summer. The seasonal inoculums could lead to differences in the bacterial community structure and diversity, while recharge water temperature was considered to be a major factor influencing the bacterial growth rate and metabolism behavior during the seasonal bioclogging process.  相似文献   

Hydrogeology of a montane headwater groundwater system downgradient of a coal-mine waste rock dump: Elk Valley,British Columbia,Canada     

Jakub T. Szmigielski  S. Lee Barbour  Sean K. Carey  John Kurylo  A. F. McClymont  M. Jim Hendry 《Hydrogeology Journal》2018,26(7):2341-2356

Steelmaking-coal waste rock placed in mountain catchments in the Elk Valley, British Columbia, Canada, drain constituents of interest (CIs) to surface water downgradient of the waste rock dumps. The role of groundwater in transporting CIs in the headwaters of mountain catchments is not well understood. This study characterizes the physical hydrogeology of a portion of a 10-km² headwater catchment (West Line Creek) downgradient of a 2.7-km² waste rock dump placed over a natural headwater valley-bottom groundwater system. The study site was instrumented with 13 monitoring wells. Drill core samples were collected to determine subsurface lithology and geotechnical properties. The groundwater system was characterized using field testing and water-level monitoring. The valley-bottom sediments were composed of unconsolidated glacial and meltwater successions (<64 m thick) deposited as a series of cut and fill structures overlying shale bedrock. An unconfined basal alluvial aquifer located above fractured bedrock was identified as the primary conduit for groundwater flow toward Line Creek (650 m from the toe of the dump). Discharge through the basal alluvial aquifer was estimated using the geometric mean hydraulic conductivity (±1 standard deviation). These calculations suggest groundwater discharge could account for approximately 15% (ranging from 2 to 60%) of the total water discharged from the watershed. The residence time from the base of the waste rock dump to Line Creek was estimated at <3 years. The groundwater system was defined as a snowmelt (i.e., nival) regime dominated by direct recharge (percolation of precipitation) across the catchment.  相似文献   

Groundwater resources evaluation in calcareous limestone using geoelectrical and VLF-EM surveys (El Salloum Basin,Egypt)     

Zarif  Fardous  Slater  Lee  Mabrouk  Mohamed  Youssef  Ahmed  Al-Temamy  Ayman  Mousa  Salah  Farag  Karam  Robinson  Judy 《Hydrogeology Journal》2018,26(4):1169-1185

Understanding and developing groundwater resources in arid regions such as El Salloum basin, along the northwestern coast of Egypt, remains a challenging issue. One-dimensional (1D) electrical sounding (ES), two-dimensional (2D) electrical resistivity imaging (ERI), and very low frequency electromagnetic (VLF-EM) measurements were used to investigate the hydrogeological framework of El Salloum basin with the aim of determining the potential for extraction of potable water. 1D resistivity sounding models were used to delineate geoelectric sections and water-bearing layers. 2D ERI highlighted decreases in resistivity with depth, attributed to clay-rich limestone combined with seawater intrusion towards the coast. A depth of investigation (DOI) index was used to constrain the information content of the images at depths up to 100 m. The VLF-EM survey identified likely faults/fractured zones across the study area. A combined analysis of the datasets of the 1D ES, 2D ERI, and VLF-EM methods identified potential zones of groundwater, the extent of seawater intrusion, and major hydrogeological structures (fracture zones) in El Salloum basin. The equivalent geologic layers suggest that the main aquifer in the basin is the fractured chalky limestone middle Miocene) south of the coastal plain of the study area. Sites likely to provide significant volumes of potable water were identified based on relatively high resistivity and thickness of laterally extensive layers. The most promising locations for drilling productive wells are in the south and southeastern parts of the region, where the potential for potable groundwater increases substantially.

  相似文献