非饱和膨胀土的土-水特征曲线研究

在干旱和半干旱地区 ,土体中含水量的变化常会引发各种工程问题。研究表明 ,非饱和土的工程性质不仅取决于土的组成、结构和应力状态 ,还与土中的吸力密切相关。非饱和土的土 -水特征曲线表达了土体中含水量与吸力的关系 ,是非饱和土研究的重要内容之一。     

Identifying and quantifying urban recharge: a review

The sources of and pathways for groundwater recharge in urban areas are more numerous and complex than in rural environments. Buildings, roads, and other surface infrastructure combine with man-made drainage networks to change the pathways for precipitation. Some direct recharge is lost, but additional recharge can occur from storm drainage systems. Large amounts of water are imported into most cities for supply, distributed through underground pipes, and collected again in sewers or septic tanks. The leaks from these pipe networks often provide substantial recharge. Sources of recharge in urban areas are identified through piezometry, chemical signatures, and water balances. All three approaches have problems. Recharge is quantified either by individual components (direct recharge, water-mains leakage, septic tanks, etc.) or holistically. Working with individual components requires large amounts of data, much of which is uncertain and is likely to lead to large uncertainties in the final result. Recommended holistic approaches include the use of groundwater modelling and solute balances, where various types of data are integrated. Urban recharge remains an under-researched topic, with few high-quality case studies reported in the literature. Electronic Publication     

Drainage in finite-sized unsaturated zones

After the initiation of gravity drainage, water is often assumed to be either (a) draining under unit gradient, or (b) at capillary/gravity equilibrium. Both of these simplifications can be useful, but the regimes of validity of each assumption must be delineated. Water pressures are measured versus time and distance as water drains out of a 1.6 m long sand column to determine the relative effects of capillary and gravitational forces during drainage. For medium sized sands (0.15–0.3 mm in diameter), the capillary pressure is constant in space in a large region of the column for over 12 days, and the water continues to flow under unit gradient for relatively long time scales. Similar results are seen for finer sands, but with a much faster approach to equilibrium. Numerical simulations and analytical estimates are presented and compare favorably to the measurements. Together, the experimental, theoretical and analytical results are used to calculate when capillary/gravity equilibrium is reached as a function of porous media properties and length of the unsaturated zone. The ratio of the length of the unsaturated zone to the bubbling pressure is a key parameter in determining the drainage regime, and that even for relatively short unsaturated zones the equilibrium time scale can be on the order of years.     

Parameters affecting maximum fluid transport in large aperture fractures

We present results of laboratory experiments to study the behavior of liquids moving in unsaturated wide-aperture fractures. A 5-mm-thick glass plate cut with a 1.7-mm aperture was used as a fractured rock analog to study behavior of film and capillary droplet flow modes. Flow rates ranged between 0.6 and 6.0 ml/min. Variability in the ambient barometric pressure, resulting from weather conditions, seemed to play a role in the natural selection of flow mode. For droplet mode, constant input conditions resulted in highly variable transport properties within the fracture. The advancing meniscus exhibited a dynamic contact angle that was a function of the droplet speed and much larger than the static contact angle. Flow rate was reduced due to the larger contact angle. Analytical expressions for droplet velocity and flow capacity are presented as a function of the dynamic rather than the static contact angle.     

An experimental investigation of LNAPL migration in an unsaturated/saturated sand

Accidental spills of hydrocarbons, such as Light Non-Aqueous Phase Liquids (LNAPLs), are one of the most common sources of subsurface contamination. Migration of LNAPL in a porous medium is influenced by various factors such as the number of fluids present in the unsaturated/saturated zones and the proportion of pores occupied by each fluid. The results for relationship between matric suction and degree of saturation are presented in this paper for water–air, water–LNAPL and LNAPL–air systems in a sand. A simple and reliable setup using Buchner funnel was designed to obtain these relations. It was found that the relationship between matric suction head and degree of saturation is hysteretic for all the fluid systems (water–air, water–LNAPL and LNAPL–air). Furthermore, the amount of hysteresis depended upon the fluid system, with the maximum hysteresis occurring for water–air system. The results suggest that the amount of trapped air depends upon the reversal degree of saturation from drying to wetting.     

Suction-controlled experiments on Boom clay

The understanding of the thermo-hydro-mechanical behaviour of a clay barrier is needed for the prediction of its final in situ properties after the hydration and thermal transient in a radioactive waste repository.

As part of the CEC 1990–1994 R&D programme on radioactive waste management and storage, the CEA (Fr), CIEMAT (Sp), ENRESA (Sp), SCK · CEN (B), UPC (Sp) and UWCC (UK) have carried out a joint project on unsaturated clay behaviour (Volckaert et al., 1996). The aim of the study is to analyse and model the behaviour of a clay-based engineered barrier during its hydration phase under real repository conditions. The hydro-mechanical and thermo-hydraulic models developed in this project have been coupled to describe stress/strain behaviour, moisture migration and heat transfer. A thermo-hydraulic model has also been coupled to a geochemical code to describe the migration and formation of chemical species.

In this project, suction-controlled experiments have been performed on Boom clay (B), FoCa clay (Fr) and Almeria bentonite (Sp). The aim of these experiments is to test the validity of the interpretive model developed by Alonso and Gens (Alonso et al., 1990), and to build a database of unsaturated clay thermo-hydro-mechanical parameters. Such a database can then be used for validation exercises in which in situ experiments are simulated.

The Boom clay is a moderately swelling clay of Rupellian age. It is studied at the SCK · CEN in Belgium as a potential host rock for a radioactive waste repository. In this paper, suction-controlled experiments carried out on Boom clay by SCK · CEN are described. SCK · CEN has performed experiments to measure the relation between suction, water content and temperature and the relation between suction, stress and deformation. The applied suction-control techniques and experimental setups are detailed. The results of these experiments are discussed in the perspective of the model of Alonso and Gens. The influence of temperature on water uptake was rather small. The measured swelling-collapse behaviour can be explained by the Alonso and Gens model.     

Soil moisture potential and water content in the unsaturated zone within the arid Ejina Oasis in Northwest China

Three soil profiles were selected in the Ejina Oasis, northwest China, to determine water content profiles and evolution of soil moisture potentials in the unsaturated zone within the arid area. The total soil moisture potentials have been monitored for about 3 months in 2001 at different depths in the soil profiles. The occurrence and movement of water in the unsaturated zone was analyzed using the zero flux plane (ZFP) method. It is shown that convergent ZFPs and divergent ZFPs may occur at depths between 0.5 and 3.0 m, and that the depth of the ZFPs was controlled by the root zone of plants growing on the land surface. Profiles of the total soil moisture potentials were observed to be coincident with those of the water contents at the three experimental sites. The total soil moisture potential showed a slight increasing trend and the ZFPs tend to vanish from summer to winter as the water extraction by roots decreased. Evapotranspiration through vegetation has an important bearing on the water content and the total potential in the unsaturated zone.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of a new method for efficiently calculating of evaporation from the phreatic aquifer in variably saturated flow modeling

In this paper, we proposed a new method that has been developed based on the surface soil moisture content (SSMC) to more efficiently calculate the groundwater evaporation in variably saturated flow modeling. In this method, the empirical formula to calculate evaporation was modified and the value of the formula varies from zero to one as a closed interval. In addition, the simulation code for calculating the groundwater evaporation based on the SSMC method was incorporated into the EOS9 module of Tough2, a variably saturated flow modeling code. Finally, two numerical tests and a case simulation were conducted to verify the feasibility and accuracy of the SSMC method. Simulation results indicate that the SSMC method is capable of appropriately simulating the characteristics of water flow in vadose zone and the amount of evaporation with the variable water table. And such results are in coincidence with the value calculated by the logistic function method, and fit well with the measured data globally rather than locally.