Effect of air mobility in a soil pipe on water flow through a model hillslope after pipe clogging

Soil pipes (continuous macropores expanding laterally in the soil subsurface) are a key factor controlling hillslope water cycles and sediment transport. Soil pipes usually enhance slope stability under rainfall events through their high water drainage ability, and pipe clogging by sediments is regarded as a risk for slope failure. In this study, we conducted a bench-scale pipe clogging experiment to clarify the effect of air mobility in soil pipes on water flow and water pressure build-up in the slope at the clogged point. Before pipe clogging, the soil pipe drained rainwater effectively and lowered the groundwater table. After the pipe clogging event, the mobility of air in the soil pipe before the clogging determined the water flow in the slope. When the air in the soil pipe connected to the atmosphere and moved freely, the water level in the soil pipe increased at the pipe clogging, and water pressure build-up was limited near the pipe outlet. On the other hand, when air in the soil pipe was entrapped by the clogging, water pressure suddenly increased, and the groundwater table of the whole slope rose correspondingly. This study clearly demonstrated the importance of pipe morphology with respect to air connectivity between the pipe and atmosphere to elucidate the water flow and slope stability during the pipe clogging event. © 2019 John Wiley & Sons, Ltd.     

Clogging of saturated porous media by silt‐sized suspended solids under varying physical conditions during managed aquifer recharge

Managed aquifer recharge is an effective method for utilizing excess flood flows, but clogging of porous media is a limiting factor in the implementation of this water storage technique. In recent years, much research on the physical clogging of porous media during artificial recharge has been conducted. However, the understanding of clogging due to silt‐sized suspended solids (SS) is still inadequate, especially under varying physical conditions. Here, we subjected sand columns to controlled rates of flow and SS suspensions to investigate the influence of media size, SS size, SS concentration, and flow velocity on the clogging of porous media by silt‐sized SS. The results show that the diameter ratio of SS particles to sand grains is the dominant factor influencing the position of physical clogging. As pore velocity increased, the mobility of silt‐sized SS was enhanced and retention in the porous media decreased noticeably. The spatial retention profiles in the porous media were found to vary greatly at different flow velocities. The SS concentration of the infiltrating suspension also dramatically influenced the mobility and deposition of silt‐sized SS particles, such that high concentrations accelerated the clogging process. As the different physical factors changed, the breakthrough curves and retention profiles of silt‐sized SS particles changed obviously and the mechanisms of retention differed. On the whole, clogging position is mainly determined by particle size ratio, but clogging rate is dominated by a variety of factors including particle size ratio, SS concentration, and flow velocity.     

垂直流人工湿地填料的淤堵机理初探

应用人工湿地技术处理污水时,填料的淤堵往往是制约其处理功效和使用寿命的主要因素.本研究采用垂直流人工湿地单元模型,通过测定填料有效孔隙率与截留SS总量、COD降解量、生物膜的增长量之间的关系,初步探讨填料淤堵机理.结果表明:进水中悬浮颗粒全部为无机颗粒的系统较进水中全部为溶解态污染物质的系统更易造成淤堵,进水中全部为溶解态污染物的系统至实验结束时未出现淤堵现象.淤堵发生后,填料上层的渗透系数明显低于中下层,表明淤堵主要发生在填料上层部分;截留的不可生物降解无机SS、被截留但未被降解的有机SS、生物膜的生长及老化脱落、以及被上述物质吸附的水是导致淤堵的主要原因.淤堵发生后,湿地日污染物去除总量降低,湿地功效受到影响.     

Clogging process by suspended solids during groundwater artificial recharge: Evidence from lab simulations and numerical modelling

Permeability reduction of infiltration media due to suspended solid (SS) clogging is the bane of groundwater artificial recharge. To overcome the clogging problem and advance the understanding of the process‐based spatial‐temporal evolution of SS clogging, a 1D laboratory column simulation was carried out, followed by numerical modelling of the experimental data in this study. It was found that clogging caused a reduction in the hydraulic conductivity (K) in the upper layer at the beginning and extended deeper to approximately 50 cm, and no reduction in K was detected below 52 cm throughout the experimental period of 129 hr. The most clogged layer spanned from the surface to a depth of 11 cm, and the middle 11–52 cm was characterized by a slight decrease in K. The clogging rates of the different layers decreased with the depth, which was based on data analysis, with the largest value of 0.038 hr^?1 in the upper 1 cm. The overall K began to decrease from the surface layer and was increasingly affected by clogging with time. A mathematical model was established to simulate the SS clogging process evolution based on considerations of the attachments and detachments of particles. Then the model was applied to perform several scenario analyses after calibration and validation using the data obtained in the experiment. The simulation results indicated that the SS concentration was much more sensitive than the groundwater depth below the land surface, and 10 days of constant recharge is recommended as the disposal cycle of the clogged layer under the given conditions.     

Biological and Physical Clogging in Infiltration Wells: Effects of Well Diameter and Gravel Pack

Gravity-driven infiltration into the shallow subsurface via small-diameter wells (SDWs), i.e., wells with an inner diameter smaller than 7.5 cm (3 inches) and no gravel pack) has proven to be a cost-efficient and flexible tool for managed aquifer recharge (MAR), as it provides relatively high recharge rates with minimal construction effort. SDWs have a significantly smaller open filter area than larger diameter wells with gravel pack, making the infiltration of low-quality waters through these wells more at risk clogging. To investigate their susceptibility for biological and physical clogging, 24 physical models with different well setups were evaluated by infiltrating either nutrient-poor but turbid water or nutrient-rich but clear water. The experiments showed that smaller diameters and the lack of a gravel pack increase the well's susceptibility to both kinds of clogging. However, this effect was observed to be much more pronounced for physical than for biological clogging. Our conclusion is that SDWs show severe disadvantages with respect to the infiltration of highly turbid waters in comparison to large diameter wells with a gravel pack. Nevertheless, this disadvantage is much less severe when it comes to the infiltration of clear but nutrient-rich waters (e.g., treated wastewater). Depending on the economic and geological circumstances of a MAR-project, this disadvantage could be outweighed by the significantly lower construction costs of SDWs.     

Mechanism of Suspended Kaolinite Particle Clogging in Porous Media During Managed Aquifer Recharge

Managed aquifer recharge is an effective strategy for urban stormwater management. Chemical ions are normally retained in stormwater and groundwater and may accelerate clogging during the recharge process. However, the effect of water chemistry on physical clogging has not previously been investigated. In this study, we investigated the hydrogeochemical mechanism of saturated porous media clogging in a series of column experiments. The column was packed with river sand and added suspensions of kaolinite particles. Calcium chloride and sodium chloride are used as representative ions to study chemical effects. We found that an increase in ionic strength resulted in retention of kaolinite solids in the column, with a breakthrough peak of C/C₀ value of 1 to 0.2. The corresponding hydraulic conductivity decreased with increased solids clogging. Divalent cations were also found to have a greater influence on kaolinite particle clogging than monovalent cations. The enhanced hydrochemical-related clogging was caused by kaolinite solids flocculating and increasing the deposition rate coefficient by 1 to 2 times in high ionic strength conditions. Three clogging mechanisms of kaolinite solids are proposed: surface filtration, inner blocking, and attachment. This study further deepens the understanding of the mechanisms of solids clogging during aquifer recharge and demonstrates the significance of ionic strength on recharge clogging risk assessments.     

Coupled Aquifer-Borehole Simulation

A model coupling fluid hydraulics in a borehole with fluid flow in an aquifer is developed in this paper. Conservation of momentum is used to create a one-dimensional steady-state model of vertical flow in an open borehole combined with radially symmetric flow in an aquifer and with inflow to the well through the wellbore screen. Both laminar and turbulent wellbore conditions are treated. The influence of inflow through the wellbore screen on vertical flow in the wellbore is included, using a relation developed by Siwoń (1987) . The influence of inflow reduces the predicted vertical variation in head up to 15% compared to a calculation of head losses due to fluid acceleration and the conventional Colebrook-White formulation of friction losses in a circular pipe. The wellbore flow model is embedded into the MODFLOW-2000 ground water flow code. The nonlinear conservation of momentum equations are iteratively linearized to calculate the conductance terms for vertical flow in the wellbore. The resulting simulations agree favorably with previously published results when the model is adjusted to meet the assumptions of the previous coupled models.     

River-induced flow dynamics in long-screen wells and impact on aqueous samples

Previously published field investigations and modeling studies have demonstrated the potential for sample bias associated with vertical wellbore flow in conventional monitoring wells constructed with long-screened intervals. This article builds on the existing body of literature by (1) demonstrating the utility of continuous (i.e., hourly measurements for ～1 month) ambient wellbore flow monitoring and (2) presenting results from a field experiment where relatively large wellbore flows (up to 4 L/min) were induced by aquifer hydrodynamics associated with a fluctuating river boundary located approximately 250 m from the test well. The observed vertical wellbore flows were strongly correlated with fluctuations in river stage, alternating between upward and downward flow throughout the monitoring period in response to changes in river stage. Continuous monitoring of ambient wellbore flows using an electromagnetic borehole flowmeter allowed these effects to be evaluated in concert with continuously monitored river-stage elevations (hourly) and aqueous uranium concentrations (daily) in a long-screen well and an adjacent multilevel well cluster. This study demonstrates that when contaminant concentrations within the aquifer vary significantly over the depth interval interrogated, river-induced vertical wellbore flow can result in variations in measured concentration that nearly encompass the full range of variation in aquifer contaminant concentration with depth.     

Analysis of Transient Flow to an Extended Fully Penetrating Well at Constant-Head Pumping

Vertical wells with radial extension at the well bottom can improve the rate of water production. No study has yet investigated the effects of the transient state and anisotropy in directional hydraulic conductivities on the wellbore flux rate for this type of well. This study derives a semianalytical transient drawdown solution for constant-head pumping at a fully penetrating well radially extended at the bottom of a confined, anisotropic aquifer by applying Laplace transform and separation of variables as well as conducting a Fourier analysis. The results of this new solution indicate that transient and steady-state wellbore flux rates can be increased by a factor of two for greater radial extension of the well. Compared with an isotropic aquifer (a ratio of vertical and horizontal hydraulic conductivities equal to one), an anisotropic aquifer with the ratio less than one may produce a higher transient wellbore flux rate and lower steady-state wellbore flux rate. Moreover, the time required to achieve the steady-state wellbore flux rate can be substantially affected by anisotropy of the aquifer.     

Distribution of clay-sized sediments in streambeds and influence of fine sediment clogging on hyporheic exchange

In this work, the deposition of clay-sized fine particles (d₅₀ = 0.006 mm) and its subsequent influence on the dune-induced hyporheic exchange are investigated. Fine sand (D₅₀ = 0.28 mm), coarse sand (D₅₀ = 1.7 mm), and gravel (D₅₀ = 5.5 mm) grains were used to form homogenous model streambeds; one control - no clay input, and two treatments - increasing clay inputs for each grain type. The results indicate that the clogging profiles of clay-sized sediments may not be predicted accurately using the previously proposed metric based on the relative sizes of infiltrating and substrate sediments. Further, the depositional patterns vary with the initial concentration of clay particles in the surface water. The assessment of clogging profiles in coarse-grained model streambeds also reveals a preferential infiltration of the clay particles in the hyporheic downwelling regions. The results from the dye tracer test suggest that the accumulation of clay particles altered the exchange characteristics in the treatment flumes. For each grain size, the treatment flumes exhibit lower hyporheic flux and higher median residence times compared to their respective control flumes. The dye penetration depths were lower in treatment flumes with fine and coarse sand compared to their respective control flumes. Interestingly, higher penetration depths were observed in treatment flumes with gravel compared to their respective control flume potentially due to the generation of preferential flow paths in the partially clogged gravel beds. The clogging altered the hyporheic fluxes and residence times in the coarse-grained model beds to a greater degree in comparison to the fine sand beds. Overall, our findings indicate that the properties of both fine and substrate sediments influence the clogging patterns in streambeds, and the subsequent influence of fine sediment clogging on hyporheic exchange and associated processes may vary across stream ecosystems.     

DNAPL accumulation in wells and DNAPL recovery from wells: Model development and application to a laboratory study

Dense nonaqueous phase liquid (DNAPL) accumulation and recovery from wells cannot be accurately modeled through typical pressure or flux boundary conditions due to gravity segregation of water and DNAPL in the wellbore, the effects of wellbore storage, and variations of wellbore inflow and outflow rates with depth, particularly in heterogeneous formations. A discrete wellbore formulation is presented for numerical modeling of DNAPL accumulation in observation wells and DNAPL removal from recovery wells. The formulation includes fluid segregation, changing water and DNAPL levels in the well and the corresponding changes in fluid storage in the wellbore. The method was added to a three-dimensional finite difference model (CompSim) for three phase (water, gas, DNAPL) flow. The model predictions are compared to three-dimensional pilot scale experiments of DNAPL (benzyl alcohol) infiltration, redistribution, recovery, and water flushing. Model predictions match experimental results well, indicating the appropriateness of the model formulation. Characterization of mixing in the extraction well is important for predicting removal of highly soluble organic compounds like benzyl alcohol. A sensitivity analysis shows that the incorporation of hysteresis is critical for accurate prediction. Among the multiphase flow and transport parameters required for modeling, results are most sensitive to soil intrinsic permeability.     

基于3D-FEM的洞穴型地层双侧向测井数值模拟及响应特征研究

洞穴型碳酸盐岩储层非均质性强、电测井响应复杂、测井识别和表征难度大,利用数值模拟方法明确该类储层的电测井响应特征,可为洞穴识别与评价提供理论依据.本文基于三维有限元素法(3D-FEM),引入边界局部加密技术,实现对复杂球形洞穴的精细刻画与响应精确计算;对比分析井眼钻穿型洞穴和井旁洞穴的双侧向测井响应特征和敏感性,进而考察双侧向测井对两者敏感范围的差异.结果表明:洞穴的存在导致双侧向测井响应明显降低,受洞穴边界及仪器探测深度等影响,井眼钻穿型洞穴双侧向测井曲线复杂,而井旁洞穴曲线呈"抛物线"型;双侧向测井对井眼钻穿型洞穴的敏感性远大于对井旁洞穴的敏感性,深侧向和浅侧向测井最大可对洞穴边界距井壁0.5 m和0.3 m的井旁洞穴敏感.     

Stochastic evaluation of the impact of sewer inlets’ hydraulic capacity on urban pluvial flooding

Sewer inlet structures are vital components of urban drainage systems and their operational conditions can largely affect the overall performance of the system. However, their hydraulic behaviour and the way in which it is affected by clogging is often overlooked in urban drainage models, thus leading to misrepresentation of system performance and, in particular, of flooding occurrence. In the present paper, a novel methodology is proposed to stochastically model stormwater urban drainage systems, taking the impact of sewer inlet operational conditions (e.g. clogging due to debris accumulation) on urban pluvial flooding into account. The proposed methodology comprises three main steps: (i) identification of sewer inlets most prone to clogging based upon a spatial analysis of their proximity to trees and evaluation of sewer inlet locations; (ii) Monte Carlo simulation of the capacity of inlets prone to clogging and subsequent simulation of flooding for each sewer inlet capacity scenario, and (iii) delineation of stochastic flood hazard maps. The proposed methodology was demonstrated using as case study design storms as well as two real storm events observed in the city of Coimbra (Portugal), which reportedly led to flooding in different areas of the catchment. The results show that sewer inlet capacity can indeed have a large impact on the occurrence of urban pluvial flooding and that it is essential to account for variations in sewer inlet capacity in urban drainage models. Overall, the stochastic methodology proposed in this study constitutes a useful tool for dealing with uncertainties in sewer inlet operational conditions and, as compared to more traditional deterministic approaches, it allows a more comprehensive assessment of urban pluvial flood hazard, which in turn enables better-informed flood risk assessment and management decisions.     

Near‐Well Nonlinear Flow Identified by Various Displacement Well Response Testing

The impact of nonlinear flow phenomena on well response tests is still not completely understood today. With the present paper, a set of 10 well response tests is investigated. The tests were conducted in a fractured Devonian limestone formation close to the western national border of Germany. The test design incorporates a packer as well as different solid cylinders to initiate a series of slug‐injection and slug‐withdrawal tests by various initial displacements. Nonlinear response characteristics were observed in the course of the tests, which cannot be explained by tubing‐controlled flow inside the cased well. The analysis shows that the nonlinear response characteristics are likely to be either formation controlled according to non‐Darcian flow developing in a high‐conductivity fracture compartment of the tested limestone formation or a consequence of a severe well inefficiency caused by some sort of screen clogging. This inference is obtained from analyzing the data by a nonlinear well response test model, which differentiates between wellbore internal hydraulic head losses and a generalized rate‐dependent skin effect accounting for nonlinear flow processes in the vicinity of the well. The potential of identifying near‐well nonlinear flow by various displacement well response testing may indicate this methodology to be a valuable complement to modern high‐resolution borehole imaging techniques used when characterizing fractured reservoirs and the tightness of fractured reservoir cap rocks.     

Numerical Modeling of Unsaturated Flow in Wastewater Soil Absorption Systems

It is common practice in the United States to use wastewater soil absorption systems (WSAS) to treat domestic wastewater. WSAS are expected to provide efficient, long-term removal of wastewater contaminants prior to ground water recharge. Soil clogging at the infiltrative surface of WSAS occurs due to the accumulation of suspended solids, organic matter, and chemical precipitates during continued wastewater infiltration. This clogging zone (CZ) creates an impedance to flow, restricting the hydraulic conductivity and rate of infiltration. A certain degree of clogging may improve the treatment of wastewater by enhancing purification processes, in part because unsaturated flow is induced and residence times are significantly increased. However, if clogging becomes excessive, the wastewater pond height at the infiltrative surface can rise to a level where system failure occurs. The numerical model HYDRUS-2D is used to simulate unsaturated flow within WSAS to better understand the effect of CZs on unsaturated flow behavior and hydraulic retention times in sandy and silty soil. The simulations indicate that sand-based WSAS with mature CZs are characterized by a more widely distributed flow regime and longer hydraulic retention times. The impact of clogging on water flow within the silt is not as substantial. For sand, increasing the hydraulic resistance of the CZ by a factor of three to four requires an increase in the pond height by as much as a factor of five to achieve the same wastewater loading. Because the degree of CZ resistance directly influences the pond height within a system, understanding the influence of the CZ on flow regimes in WSAS is critical in optimizing system design to achieve the desired pollutant-treatment efficiency and to prolong system life.     

Factors Governing Sustainable Groundwater Pumping near a River

The objective of this paper was to provide new insights into processes affecting riverbank filtration (RBF). We consider a system with an inflatable dam installed for enhancing water production from downstream collector wells. Using a numerical model, we investigate the impact of groundwater pumping and dam operation on the hydrodynamics in the aquifer and water production. We focus our study on two processes that potentially limit water production of an RBF system: the development of an unsaturated zone and riverbed clogging. We quantify river clogging by calibrating a time‐dependent riverbed permeability function based on knowledge of pumping rate, river stage, and temperature. The dynamics of the estimated riverbed permeability reflects clogging and scouring mechanisms. Our results indicate that (1) riverbed permeability is the dominant factor affecting infiltration needed for sustainable RBF production; (2) dam operation can influence pumping efficiency and prevent the development of an unsaturated zone beneath the riverbed only under conditions of sufficient riverbed permeability; (3) slow river velocity, caused by dam raising during summer months, may lead to sedimentation and deposition of fine‐grained material within the riverbed, which may clog the riverbed, limiting recharge to the collector wells and contributing to the development of an unsaturated zone beneath the riverbed; and (4) higher river flow velocities, caused by dam lowering during winter storms, scour the riverbed and thus increase its permeability. These insights can be used as the basis for developing sustainable water management of a RBF system.     

Two‐dimensional modelling of large wood transport during flash floods

Large woody material (LWM) transported by rivers may be entrapped at critical stream geometry configurations (e.g. bridges) and therefore dramatically increase the destructive power of floods. This was the case in a Spanish mountain river where a flood event with a high degree of LWM transport took place in 1997. The aim of this study was to simulate a bridge clogging process and reconstruct the wood deposit patterns, modelling individual pieces of wood moving with the water flow and interacting among them and with the bridge. A two‐dimensional numerical model was developed to simulate the transport of LWM and its effect on hydrodynamics. Different scenarios for the wood transport rate allowed us to study the influence of inlet boundary conditions on bridge clogging. For the studied event, the scenario which best reproduced the bridge clogging effect and flood characteristics was one in which 60% of the total wood entered before the peak discharge. This dropped to 30% at the peak itself, and finally fell to 10% during the recession curve. In addition, the accumulation patterns of LWM along the reach were computed and compared with post‐event field photographs, showing that the model succeeded in predicting the deposition patterns of wood and those areas prone to form wood jams. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparing MODFLOW simulation options for predicting intra‐meander flux

During the evolution of meander bends, the intra‐meander groundwater head gradients steepen and generate zones of accelerated water and nutrient intra‐meander fluxes important for ecosystem processes. This paper compares and contrasts three MODFLOW groundwater model packages based on their simulation of intra‐meander flux for two stages of meander evolution observed in a sandbox river table and one level of river bed clogging, where the hydraulic conductivity in the river bed is lower than in the adjacent aquifer. These packages are the Time‐Variant Specified Head package [constant head (CHD)], River package (RIV), and Streamflow‐Routing package (SFR2), each controlling the groundwater or river head bounding the intra‐meander region. The RIV and SFR2 packages fix river stage and allow for variation in groundwater head below the river, which is suggested for simulating intra‐meander flux for all sinuosities with and without river bed clogging whenever river bed parameters are available. The CHD package fixes below river groundwater head and fails to simulate intra‐meander head loss and flux in meanders with high sinuosity or river bed clogging. In low sinuosity meanders and in cases without river bed clogging, there were no significant differences between MODFLOW packages for simulating river intra‐meander head loss and flux. This research demonstrates why MODFLOW users need to consider the limitations of each package when simulating intra‐meander flux in reaches with river bed clogging, high sinuosity, or similarly steep hydraulic gradients. Copyright © 2014 John Wiley & Sons, Ltd.     

川06井水位固体潮效应变化初探

通过分析承压井孔储水效应对井水位固体潮效应(井潮)的影响, 重点研究了含水层储水系数S和导水系数T对潮汐排水响应的影响, 以探讨井潮变化形成的初步原因. 以川06井为例, 通过对实测井孔水位数据进行潮汐分析, 得到井孔储水效应引起的相对振幅A和相位移η, 进而分析A和η随时间变化的规律以及两者之间的关系. 结果表明, 井潮变化主要受井孔储水效应变化影响; 在潮汐排水响应阶段, 井潮主要受含水层导水系数T(或渗透系数K)影响, 而受储水系数S影响较小.      

Concentration Profiles in Screened Wells under Static and Pumped Conditions

Purge and pump samples from screened wells reflect concentration averaging and contaminant redistribution by wellbore flow. These issues were assessed in a screened well at the Hanford Site by investigating the vertical profile of a technetium-99 plume in a conventional well under static and pumped conditions. Specific conductance and technetium-99 concentrations were well correlated, and this enabled measurement of specific conductance to be used as a surrogate for technetium-99 concentration. Time-series measurements were collected during purging from three specific conductance probes installed in the well at 1.2, 3.1, and 4.9 m below the static water level in a 7.7-m-deep screened well. The vertical contaminant profile adjacent to the well in the aquifer was calculated using the concentration profile in the well during pumping, the pumping flow rate, and a wellbore flow and mixing model. The plume was found to be stratified in the aquifer—the highest concentrations occurred adjacent to the upper part of the screened interval. The purge and pump sample concentrations were 41% to 58% of the calculated peak concentration in the aquifer. Plume stratification in the aquifer adjacent to the well screen became more pronounced as pumping continued. Extended pumping may have partially reversed the effect of contaminant redistribution in the aquifer by wellbore flow and allowed the stratification of the plume to be more observable. It was also found that the vertical profile of contamination in the well under static (i.e., nonpumping conditions) was not representative of the profile in the aquifer. Thus, passive or micropurge sampling techniques, which sample the wellbore water at different depths, would not yield results representative of the aquifer in this well.