非饱和透明土优先流迁移规律分析

为实现土的非饱和优先流迁移可视化,设计立柱装置进行非饱和渗流试验,采用透明土与数字图像处理技术,建立归一化像素强度与透明土饱和度的关系,在此基础上,通过室内模型试验,研究优先流路径连通性与相邻优先流路径旋转角对非饱和土优先流迁移的影响。试验结果表明：基于图像灰度像素强度表征非饱和透明土饱和度方法是可行的;全连通优先流（O-O型）与上连通优先流（O-C型）剖面呈现T型,中心轴剖面饱和度与边缘差异明显,下连通优先流路径（C-O型）中土体油压与基质势不足以使流体进入优先流路径形成优先流,入渗趋势与均匀流一致,O-C型优先流的稳定入渗率和湿润锋推移速度分别为C-O型的1.5倍和1.4倍;相邻O-C型优先流之间区域形成新的优先流,土体达到较高饱和度,增长速率随旋转角的增加而减小,优先流转角为90?、60?、30?时稳定入渗速率分别为均匀流的1.5、1.3、1.2倍;流体受重力影响,转角小的优先流仅沿路径一侧水平渗透,但湿润锋推移速度分别为均匀流的1.3、1.4、1.5倍,相邻优先流相互作用减弱,难以形成新的优先流。     

Single phase water flow through rock fractures

Flow analysis plays a major role in various geotechnical applications, and the understanding of flow mechanisms is essential for the development of a hydro-mechanical flow model suitable for underground excavations in rock. Discrete flow analysis through discontinuities is reviewed including empirical and analytical flow models. The influence of external loading on joint deformation and single-phase flow show that the surface roughness and aperture size are the prime factors influencing flow rate. Nevertheless, the idealization of natural fractures as smooth parallel plate joints is still followed in many numerical models, because of the simplicity of the cubic law when applied to fracture networks. A numerical study of water flow through a network of joints employing Universal Distinct Element Code (UDEC) is used to quantify the effects of joint orientation and external stress acting on idealized joints.It is found that, for the same joint spacing, the flow rate into an excavation depends on the boundary block size (A_b) relative to the excavation size (A_e). The inflow becomes excessive if A_b/A_e is less than 4, but becomes very small if A_b/A_e exceeds 8.     

Modeling preferential water flow and solute transport in unsaturated soil using the active region model

Preferential flow and solute transport are common processes in the unsaturated soil, in which distributions of soil water content and solute concentrations are often characterized as fractal patterns. An active region model (ARM) was recently proposed to describe the preferential flow and transport patterns. In this study, ARM governing equations were derived to model the preferential soil water flow and solute transport processes. To evaluate the ARM equations, dye infiltration experiments were conducted, in which distributions of soil water content and Cl⁻ concentration were measured. Predicted results using the ARM and the mobile–immobile region model (MIM) were compared with the measured distributions of soil water content and Cl⁻ concentration. Although both the ARM and the MIM are two-region models, they are fundamentally different in terms of treatments of the flow region. The models were evaluated based on the modeling efficiency (ME). The MIM provided relatively poor prediction results of the preferential flow and transport with negative ME values or positive ME values less than 0.4. On the contrary, predicted distributions of soil water content and Cl⁻ concentration using the ARM agreed reasonably well with the experimental data, with ME values higher than 0.8. The results indicated that the ARM successfully captured the macroscopic behavior of preferential flow and solute transport in the unsaturated soil.     

Micromechanical approach to unsaturated water flow in structured geomaterials by two-scale computations

This article presents a micromechanical approach to the problem of unsaturated water flow in heterogeneous porous media in transient conditions. The numerical formulation is based on the two-scale model obtained previously by periodic homogenization. It allows for a coupled solution of the non-linear flow equations at macroscopic and microscopic scales and takes into account the macroscopic anisotropy of the medium and the local non-equilibrium of the capillary pressure. The model was applied to simulate two-dimensional water infiltration at constant flux into an initially dry medium containing inclusions of square and rectangular shapes. The numerical results showed the influence of the inclusion–matrix conductivity ratio and the local geometry on the macroscopic behavior. The influence of the conductivity ratio manifested itself by the acceleration or retardation of the onset of the macroscopic water flux at the outlet, while the local geometry (anisotropy) significantly affected the macroscopic spatial distribution of the water flux. Such type of approach can be extended to simulate coupled phenomena (for example hydro-mechanical problems) with evolving local geometry.     

Determining the active region model parameter from dye staining experiments for characterizing the preferential flow heterogeneity in unsaturated soils

The active region model (ARM) has been developed as a practical and effective approach for characterizing and representing preferential flow patterns in unsaturated soils. However, studies on methods to determine the ARM parameter (γ) are very limited. The major objective of this work was to refine the methods for determining the ARM parameter (γ) using the data from field-scale dye staining experiments. For this purpose, 13 field-scale dye staining experiments were conducted in silty clay, loam and sand with various initial and boundary conditions. The distributions of soil water content and fraction of stained region were measured to determine the ARM parameter (γ). Three determination methods of ARM parameter (γ) were presented according to the different distribution patterns of the soil water content of stained region. The efficiency of these approaches was demonstrated with the results from 13 field-scale dye staining tests that were directly related to preferential flow patterns. Impacts of soil texture on soil water redistribution, and effects of initial soil water content and infiltration amount on the preferential flow heterogeneity were also discussed in this research.     

Preferential flow,diffuse flow,and perching in an interbedded fractured-rock unsaturated zone

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.     

基于流动电位正演模型的非饱和带水流过程监测

为了获取非饱和带水流过程的信息,借助流动电位正演模型,通过数值实验探讨非降雨和降雨两种条件下非饱和带流动电位和水流过程的关系,然后用南京中山植物园试验场地野外观测的流动电位和张力数据加以对比和验证。野外试验表明:流动电位可以有效地反映非饱和带水流过程。在夏季无降雨入渗的条件下,日周期变化的地表地下温度差导致水分的运动,流动电位准确地指示了非饱和带含水量和毛细压力的变化情况,从而指示出了水分运移的方向;在夏季有降雨入渗的条件下,降雨锋面推进之处,含水量和流动电位同时有明显的响应,进而根据不同位置的流动电位对降雨入渗响应的时刻差,直接求出入渗锋面的推进速度。     

Geochemical exploration for mineralized breccia pipes in northern Arizona,U.S.A.

Thousands of solution-collapse breccia pipe crop out in the canyons and on the plateaus of northern Arizona. Over 80 of these are known to contain U or Cu mineralized rock. The high-grade U ore associated with potentially economic concentrations of Ag, Pb, Zn, Cu, Co and Ni in some of these pipes has continued to stimulate mining and exploration activity in northern Arizona, despite periods of depressed U prices. Large expanses of northern Arizona are comprised of undissected high plateaus; recognition of pipes in these areas is particularly important because mining access to the plateaus is far better than to the canyons. The small size of the pipes, generally less than 600 ft (200 m) in diameter, and limited rock outcrop on the plateaus, compounds the recognition problem. Although the breccia pipes, which bottom in the Mississippian Redwall Limestone, are occasionally exposed on the plateaus as circular features, so are unmineralized near-surface collapse features that bottom in the Permian Kaibab and Toroweap Formations. The distinction between these two classes of circular features is critical during exploration for this unique type of U deposit.Various geochemical and geophysical exploration methods have been tested over these classes of collapse features. Because of the small size of the deposits, and the low-level geochemical signatures in the overlying rock that are rarely dispersed for distances in excess of several hundred feet, most reconnaissance geochemical surveys, such as hydrogeochemistry or stream sediment, will not delineete mineralized pipes.Several types of detailed geochemical surveys made over collapse features, located through examination of aerial photographs and later field mapping, have been successful at delineating collapse features from the surrounding host rock: (1) Rock geochemistry commonly shows low level Ag, As, Ba, Co, Cu, Ni, Pb, Se and Zn anomalies over mineralized breccia pipes; (2) Soil surveys appear to have the greatest potential for distinguishing mineralized breccia pipes from the surrounding terrane. Although the soil anomalies are only twice the background concentrations for most anomalous elements, traverses made over collapse features show consistent enrichment inside of the feature as compared to outside; (3) B. Cereus surveys over a known mineralized pipe show significantly more anomalous samples collected from within the ring fracture than from outside of the breccia pipe; (4) Helium soil-gas surveys were made over 7 collapse features with discouraging results from 5 of the 7 features.Geophysical surveys indicate that scaler audio-magnetotelluric (AMT) and E-field telluric profile data show diagnostic conductivity differences over mineralized pipes as compared to the surrounding terrane. These surveys, coupled with the geochemical surveys conducted as detailed studies over features mapped by field and aerial photograph examination, can be a significant asset in the selection of potential breccia pipes for drilling.     

Partially to fully saturated flow through smooth,clean, open fractures: qualitative experimental studies

Fractures are both rough and irregular but can be expressed by a simple model concept of two smooth parallel plates and the associated cubic law governing discharge through saturated fractures. However, in natural conditions and in the intermediate vadose zone, these assumptions are likely violated. This paper presents a qualitative experimental study investigating the cubic law under variable saturation in initially dry free-draining discrete fractures. The study comprised flow visualisation experiments conducted on transparent replicas of smooth parallel plates with inlet conditions of constant pressure and differing flow rates over both vertical and horizontal inclination. Flow conditions were altered to investigate the influence of intermittent and continuous influx scenarios. Findings from this research proved, for instance, that saturated laminar flow is not likely achieved, especially in nonhorizontal fractures. In vertical fractures, preferential flow occupies the minority of cross-sectional area despite the water supply. Movement of water through the fractured vadose zone therefore becomes a matter of the continuity principle, whereby water should theoretically be transported downward at significantly higher flow rates given the very low degree of water saturation. Current techniques that aim to quantify discrete fracture flow, notably at partial saturation, are questionable. Inspired by the results of this study, it is therefore hypothetically improbable to achieve saturation in vertical fractures under free-draining wetting conditions. It does become possible under extremely excessive water inflows or when not free-draining; however, the converse is not true, as a wet vertical fracture can be drained.     

Modeling root water uptake when calculating unsaturated flow in the vadose zone and groundwater recharge

This paper considers models of the transpiration root water uptake and principles for substantiating their parameters when calculating unsaturated flow in the vadose zone and estimating infiltration groundwater recharge. The possible use of generalized parameters for different landscape and soil conditions is substantiated based on analyzing the natural patterns of water balance formation in the vadose zone and the sensitivity of transpiration models. Recommendations are given for selecting the models of transpiration water uptake and substantiating their parameters, which enables one to minimize the errors in estimating infiltration groundwater recharge.     

Recharge processes: piston flow vs preferential flow in semi-arid aquifers of India

Study of groundwater recharge processes is vital for quantification of total natural recharge to the aquifers. One of the recharge processes demonstrated earlier by tracer experiments in the unsaturated zone is that of piston flow movement of soil moisture. Based on this recharge process, environmental tritium, chloride and injected tritium studies have been carried out extensively in various geological environs of India. The purpose of this paper is to evaluate the validity of the piston flow concept in different geological environs viz. consolidated fractured and weathered granites, semi-consolidated sandstones and unconsolidated alluvial tracts, and quantify the contribution from this process as well as that from the preferential flow mechanism using different tracers. Analysis of tracer data demonstrates that the preferential flow recharge process contributes very significantly (an average of 75% of total recharge) in the case of fractured granites and is important (an average of 33% of total recharge) for semi-consolidated sandstones, whereas the preferential flow recharge component is minimal in unconsolidated alluvial tracts (piston flow model is applicable). These findings necessitate re-evaluation of the total natural recharge potential of the above mentioned geological environs in view of the significant preferential flow recharge that is evidenced and estimated. Electronic Publication     

Evaluation of probabilistic flow in two unsaturated soils

A variably saturated flow model is coupled to a first-order reliability algorithm to simulate unsaturated flow in two soils. The unsaturated soil properties are considered as uncertain variables with means, standard deviations, and marginal probability distributions. Thus, each simulation constitutes an unsaturated probability flow event. Sensitivities of the uncertain variables are estimated for each event. The unsaturated hydraulic properties of a fine-textured soil and a coarse-textured soil are used. The properties are based on the van Genuchten model. The flow domain has a recharge surface, a seepage boundary along the bottom, and a no-flow boundary along the sides. The uncertain variables are saturated water content, residual water content, van Genuchten model parameters alpha (α) and n, and saturated hydraulic conductivity. The objective is to evaluate the significance of each uncertain variable to the probabilistic flow. Under wet conditions, saturated water content and residual water content are the most significant uncertain variables in the sand. For dry conditions in the sand, however, the van Genuchten model parameters α and n are the most significant. Model parameter n and saturated hydraulic conductivity are the most significant for the wet clay loam. Saturated water content is most significant for the dry clay loam. Electronic Publication     

Effect of macropore tortuosity and morphology on preferential flow through saturated soil: A Lattice Boltzmann study

Unlike micropores where water moves upward or downward based on hydraulic gradient, in macropores, water flows predominantly downward due to the gravity. Therefore, models based on capillary flow are not capable of simulating macropore flow. There are attempts to model the macropore flow using two domains, one for capillary flow and another one for macropores. These models use Richard’s equation for capillary flow and Poiseuille’s law for macropores in which the macropore is approximated to be cylindrical or planar. This study quantifies the magnitudes of the errors induced by this assumption. Influence of macropore shapes and tortuosity was quantified by using a 3D Lattice Boltzmann model, which is capable of simulating fluid flow in micropores as well as macropores of cracked clays. Artificial macropores of constant sectional area and volume, but different shapes were generated in 3D and the influence of macropore shapes, shape related parameters, and tortuosity were systematically investigated. Macropore flow rate decreases with different shapes compared to cylindrical macropores and increase in aspect ratio of sectional shape leads to decrease in macropore flow rate. The maximum effect of bends/turnings along the tortuous macropore was about 25% on overall decrease of flow rate due to tortuosity. However, more detailed study is required on the influence of bends on macropore flow rate. The macropore flow rate reduces by about 70% for tortuosity of 1.41. A prediction equation is verified to predict the flow rate of different shapes and tortuous macropores based on straight cylindrical macropore using aspect ratio and tortuosity factor.     

Geochemical and isotopic evidence for upward flow of saline fluid to the Mississippi River Valley alluvial aquifer,southeastern Arkansas,USA

Hydrogeology Journal - Groundwater from the Quaternary Mississippi River Valley Alluvial (MRVA) aquifer in southeastern Arkansas (SE AR), USA, has higher salinity compared to other MRVA...     

Geochemical evidence of aeolian deposits in European soils

Hafnium (Hf) and zirconium (Zr) concentrations measured in over 4100 agricultural soil samples from Europe were assessed with the focus on their relationship to the distribution of aeolian deposits, such as loess and coversands. Comparison of extractable (aqua regia; ICP‐MS) and total (XRFS) concentrations shows that only 1.0 to 1.7% of the total Hf and Zr is chemically extractable because of the resistant nature of their host minerals. Resistate minerals, such as zircon, are commonly found in the predominantly silty fraction of loess deposits. In this study a statistical analysis of total Hf and Zr soil data from areas with and without loess was carried out to derive threshold values of 10 mg kg^?1 (Hf) and 318 mg kg^?1 (Zr). These values were subsequently applied across the project area in an attempt to indicate the presence of aeolian deposits. The spatial distribution of above‐threshold concentrations suggests a more extensive and coherent loess belt across central and eastern Europe, providing additional evidence of loess across Brittany, Aquitaine and near the Vosges mountains in France as well as in the basins of central and northern Spain. Above‐threshold concentrations were also detected in loess regardless of its thickness, emphasising the importance of the abundance of zircon in the upper part of the soil profile rather than the actual thickness of the deposit. Soil data however, failed to indicate various loess facies within most of the Pannonian Basin, suggesting that this approach only works where deposits contain a sufficient amount of zircon. This may also explain why the extensive coversands across northern Germany and Poland were largely undetected. This study demonstrates that continental‐scale soil geochemical data can help identify and map the distribution of zircon‐rich loess and coversand, and subsequently enhance and improve current knowledge of the extent of these deposits.     

Micro- and macro-behaviour of fluid flow through rock fractures: an experimental study

Microscopic and macroscopic behaviour of fluid flow through rough-walled rock fractures was experimentally investigated. Advanced microfluidic technology was introduced to examine the microscopic viscous and inertial effects of water flow through rock fractures in the vicinity of voids under different flow velocities, while the macroscopic behaviour of fracture flow was investigated by carrying out triaxial flow tests through fractured sandstone under confining stresses ranging from 0.5 to 3.0 MPa. The flow tests show that the microscopic inertial forces increase with the flow velocity with significant effects on the local flow pattern near the voids. With the increase in flow velocity, the deviation of the flow trajectories is reduced but small eddies appear inside the cavities. The results of the macroscopic flow tests show that the linear Darcy flow occurs for mated rock fractures due to small aperture, while a nonlinear deviation of the flow occurs at relatively high Reynolds numbers in non-mated rock fracture (Re?>?32). The microscopic experiments suggest that the pressure loss consumed by the eddies inside cavities could contribute to the nonlinear fluid flow behaviour through rock joints. It is found that such nonlinear flow behaviour is best matched with the quadratic-termed Forchheimer equation.     

Estimation of actual evapotranspiration by numerical modeling of water flow in the unsaturated zone: a case study in Buenos Aires, Argentina

A method is presented to estimate actual evapotranspiration (ET_A) from potential evapotranspiration (ET_P) by numerical modeling of water flow in the unsaturated zone. Water flow is described by the Richards equation with a sink term representing the root water uptake. Evaporation is included in the model as a Neumann boundary condition at the soil surface. The Richards equation is solved in a one-dimensional domain using a mixed finite element method. The values of ET_A are obtained by applying a water stress factor to ET_P to account for soil moisture changes during the simulation period. The proposed numerical model is used to estimate ET_A in an experimental plot located in a flatland area in Buenos Aires (Argentina). Numerical results show that the proposed model is a useful tool for evaluating evapotranspiration under different scenarios.     

Tracing of water movement through the unsaturated zone of a coarse gravel aquifer by means of dye and deuterated water

To better understand the movement and transport of water and pollution through the coarse gravel unsaturated zone, the presented research was conducted to estimate water flow and transport processes with a tracing experiment in a lysimeter in the Selniška Dobrava. A combined tracing experiment was performed with deuterated water and the fluorescent dye—uranine. The fastest and dominant flow velocities were calculated based on injection time, the first tracer appearance time and the time of highest concentration. Mean flow velocity and vertical dispersion were estimated by an analytical best-fit method using one-dimensional convection–dispersion model. Deuterium was confirmed as an ideal conservative tracer and a more suitable tracer than dye (uranine) for the study of water flow in the unsaturated zone of a coarse gravel aquifer. The retardation factor of the dye as compared with deuterium was 1.13–1.75, which is in agreement with previously published results. Artificial tracers, especially deuterated water, were also identified as a very useful tool to assess other properties and differences in water flow in the unsaturated zone of a coarse gravel aquifer such as velocity and dispersion.     

Geochemical and isotopic characteristics of lower crustal xenoliths, San Francisco Volcanic Field, Arizona, U.S.A

A wide variety of xenoliths has been entrained in Miocene-to-Recent alkali olivine and hypersthene-normative basalts in the San Francisco Volcanic Field (SFVF), northern Arizona, U.S.A. Based on petrography, mineralogy, bulk rock chemistry and Sr-Nd isotopic characteristics, SFVF xenoliths can be divided into two major groups: cumulates and granulites. The cumulates are genetically related to the Cenozoic volcanic rocks and represent under- and/or intraplated additions to the crust of the Colorado Plateau. Assemblages are mafic to ultramafic and are dominated by clinopyroxene-orthopyroxene-plagioclase-spinel-amphibole-olivine. The granulites are probably Proterozoic in age, mafic-to-intermediate/felsic in bulk composition, either two pyroxene-plagioclase-spinel or plagioclase-alkali feldspar-quartz-magnetite-amphibole-biotite assemblages. Many of the granulites show evidence of partial melting. Some high SiO₂, very high Rb/Sr glasses are close in composition to erupted rhyolites, and probably represent end-member melts that have interacted with basalt to produce a variety of hybrid intermediate lavas. The major element, trace element and Sr-Nd isotope geochemistry is highly variable in the SFVF xenoliths. Extremely high Ba contents and Ba/Nb of a number of the granulites are equivalent to values characteristics of modern supra-subduction zone magmas. The considerable variation of chemical and isotopic composition depends upon mineral proportions, assemblages and chemistry. Isotopically, three end-members can be identified within the granulites: (i) lowest ⁸⁷Sr/⁸⁶Sr (0.702870) with low ¹⁴³Nd/¹⁴⁴Nd (0.511541, εNd-21.4); (ii) high ⁸⁷Sr/⁸⁶Sr (0.711069) with the lowest ¹⁴³Nd/¹⁴⁴Nd (0.511434, εNd-23.5); (iii) highest ⁸⁷Sr/{86}Sr (0.715306) with low ¹⁴³Nd/¹⁴⁴Nd (0.511793, εNd-16.5). Two important age ranges deduced from the isotopic data probably relate to episodes of crustal-growth beneath the SFVF (1.88 ± 0.33 Ga and Cenozoic). Thermobarometric calculations assuming equilibrium show that the xenoliths are derived from the lower crust (0.6–1.3 GPa, 850–1050 °C). The average SFVF lower crust is mafic in composition. In the absence of partial lithospheric delamination, the lower crust may become mafic with time due to under- and intraplating of continental crust by mafic magmas derived from the mantle.     

Movement of water through the thick unsaturated zone underlying Oro Grande and Sheep Creek Washes in the western Mojave Desert,USA

Previous studies indicate that a small quantity of recharge occurs from infiltration of streamflow in intermittent streams in the upper Mojave River basin, in the western Mojave Desert, near Victorville, California. Chloride, tritium, and stable isotope data collected in the unsaturated zone between 1994 and 1998 from boreholes drilled in Oro Grande and Sheep Creek Washes indicate that infiltration of streamflow occurs to depths below the root zone, and presumably to the water table, along much of Oro Grande Wash and near the mountain front along Sheep Creek Wash. Differences in infiltration at sites along each wash are the result of hydrologic variables such as proximity to the mountain front, quantity of streamflow, and texture of the subsurface deposits. Differences in infiltration between the washes are the result of large-scale geomorphic processes. For example, Oro Grande wash is incised into the Victorville fan and infiltration has occurred at approximately the same location over recent geologic time. In contrast, Sheep Creek Wash overlies an active alluvial fan and the stream channel can move across the fan surface through time. Infiltration does not occur to depths below the root zone at control sites outside of the washes. Electronic Publication