三维交叉裂隙渗流传质特性数值模拟

对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。     

Effects of long-term irrigation with treated wastewater. Part I: Evolution of soil physico-chemical properties

The long-term impact of irrigation on a Mediterranean sandy soil irrigated with Treated wastewater (TWW) since 1980 was evaluated. The main soil properties (CEC, pH, size distribution, exchangeable cations and chloride, hydraulic conductivity) as well as the organic matter and Cu, Cr and Pb speciation in an irrigated soil and a non-irrigated control soil at various soil depths were monitored and compared during a 2 years experiment. In this first part, the evolution of the physico-chemical soil properties was described. The irrigation with TWW was beneficial with regard to water and nutrient supplying. All the exchangeable cations other than K⁺ were higher in the irrigated soil than in the reference one. A part of the exchangeable cations was not fixed on the exchange complex but stored as labile salts or in concentrated soil solution. Despite the very sandy soil texture, both saturated and unsaturated hydraulic conductivity exhibited a significant diminution in the irrigated soil, but remained high enough to allow water percolation during rainy periods and subsequent leaching of accumulated salts, preventing soil salinization. In the irrigated soil, exchangeable sodium percentage (ESP) exhibited high values (20% on average) and the soil organic C was lower than in the reference. No significant effect was noticed on soil mineralogical composition due to irrigation.     

系统维度对变密度溶质运移的影响研究

关于地下水变密度流与溶质运移的研究通常局限于二维系统中,三维系统中的研究非常有限。然而,仍有零星研究表明,三维系统中的对流过程与二维系统不同。文章通过SEAWAT-2000进行数值模拟,系统研究了二维和三维系统中的变密度溶质运移过程,利用舍伍德数、空间矩和稀释指数量化了系统的不稳定性和溶质扩散、稀释程度。结果表明:二维系统产生分散指流,而三维系统中因扩散程度的增强使得分散指流的产生受到抑制,但是,三维系统的不稳定性比二维系统更强,对流入渗更快,与指流的产生与否并无直接关联。另外,三维系统中溶质的稀释程度大于二维系统,但溶质在二维系统中将更快接近稀释的最大值,传统二阶中心矩可能会造成自由对流和不稳定性溶质运移过程中溶质扩散和稀释程度的错误估计。研究结果将有助于正确预测三维自然含水层中的自由对流和溶质运移。     

Assesing the effect of an olive mill wastewater evaporation pond in Sousse, Tunisia

Olive oil is a typical and valuable agro-industrial product in Mediterranean countries. In Tunisia, olive mill wastewaters (OMW) reach an amount of about 1,000,000 t year⁻¹ and constitute a serious organic pollution risk because of the high chemical oxygen demand values and the presence of phytotoxic and antibacterial polyphenols. OMW have been generally stored in pond sites to be eliminated by natural evaporation or valorised by spreading on cultivated soils or by composting. Many researches on the interactions of OMW with soils at laboratory scale (columns) have been reported, but less attention have been paid to the effect of OMW on soils at field scale. The aim of this work is to investigate an area used for >15 years as an uncontrolled OMW pond site. The transformations of soil properties and groundwater occurring during OMW storage were characterised by the pH, phenolic contents, electrical conductivity (EC), moisture content and organic contents. The soil samples were taken from two borings and compared to those of a control one located near the pond site. Groundwater samples were taken on the accessible and nearest water wells to the evaporation ponds. The permeable silty and sandy layers in the site support the infiltration of OMW near the evaporation ponds. This infiltration has reached a depth of 6 m at a distance of almost 50 m laterally. The results show that the OMW infiltration in the subsoil has affected the pH, EC, organic content, phenolic compounds and the moisture.     

Preferential water and solute transport through sandy soil containing artificial macropores

Macropores resulting from soil pedogenesis and biological activity play important roles in soil water and chemical transport. Numerous studies have examined individual macropores and the effects of their size on solute transport, but few have assessed the effects of macropore continuity and of neighboring macropores. This paper describes a laboratory investigation of the effects of macropores, with varying degrees and types of continuity, on the transport and distribution of solutes in a sandy soil from the northern Loess Plateau, China. Breakthrough curves were obtained from 60 cm tall, 2-D columns containing standardized artificial macropores using an input solution of 1,190 mg/L KBr and 100 mg/L FD&C Blue #1 under a constant hydraulic head of 8 cm. The types of macropore were: open at both the surface and bottom of the soil column (O–O); open at the surface, closed at the bottom (O–C); and closed at the surface, open at the bottom (C–O). Columns with no macropores served as a control. In the O–O columns the solution reached the bottom 10–50 times faster than in any other treatment, bypassing most of the soil matrix. The presence of an O–C macropore resulted in weak retardation and much deeper penetration of the bromide and FD&C Blue #1 solution than in the control columns. However, the C–O macropore had little effect on either breakthrough curves or solute distributions. In further experiments that considered neighboring macropores effects, an inclined macropore strongly affected solute concentrations in the profile around a nearby vertical macropore. It was concluded that the length, type and position of single macropores, and the presence of neighboring macropores, all affect soil water flow and solute infiltration parameters in a sandy loam soil.     

土壤水分与溶质运移机制的分形理论研究进展

土壤中水分和溶质运移一直是土壤-水环境系统中的研究热点,也是目前仍未得到很好解决的问题.将分形理论应用于土壤水分溶质运移机制的研究,探讨其领域中的众多复杂问题,是一种全新的思路和方法.在对土壤结构定量化的分形表征进行简要阐述的基础上,重点介绍了分形理论在水分特征曲线和水力传导度等土壤水分运动基本参数的分形模型、土壤水分运移过程模拟、土壤溶质运移的非费克现象、弥散度的尺度效应以及溶质运移机制研究方面所取得的一系列成果,并就分形理论今后在土壤水分、溶质运移研究中的应用作了展望.     

膜下滴灌微区环境对土壤水盐运移的影响

基于膜下滴灌特有的"膜中"、"膜间"、"膜边"、"膜外"微区环境,利用2011—2013年田间对比试验方法获取的5 960个数据,运用柯布-道格拉斯模型,构建膜下滴灌环境土壤层次、灌水定额、土壤水分、气温、蒸发综合因素与土壤水盐关系及影响效应分析模型.结果表明,在气候干旱、蒸发强烈灌区,地膜覆盖与滴灌结合的地表介面灌溉形式下,土壤水盐具有水平方向由"膜中"向"膜边"地表裸露区定向迁移,垂直方向土壤水盐则由下向上层运移且趋于"膜外"边界积累的趋势,尤其是气温与蒸发因素交互作用,推进膜下滴灌土壤水盐在地膜覆盖与土壤裸露区域空间运移,研究结果进一步揭示了膜下滴灌"土壤水盐定向迁移"形成机理,为膜下滴灌土壤水盐地表排放模式应用提供了依据.     

Effect of soil hydraulic properties model on slope stability analysis based on strength reduction method

Soil hydraulic properties models which have been proposed were derived based on the empirical fitting curve such as Brooks-Corey model (BC) and Van Genuchten model (VG), or based on soil pore radius distribution such as Lognormal model (LN). Each model has different accuracy for predicting soil moisture distribution. In the analysis of rainfall-induced slope failure, the soil hydraulic properties model was needed to describe the physical phenomena of behavior characteristic of water in unsaturated soil. As moisture content has an effect on soil strength, it is vital to select the suitable soil hydraulic properties model for predicting Factor of Safety (FOS) especially in forecasting landslide hazard. In this study, a numerical model of seepage finite element analysis using BC, VG, and LN model were used and compared in order to analyze the soil moisture distribution, water movement phenomenon, and slope stability characteristic in unsaturated soil slope based on the strength reduction method (SRM). The results showed that the variations of the parameters predicting the moisture content of soil leads to differences of FOS in some cases. The parametric study showed that for the unsaturated soil condition, BC model has the greatest FOS value than the other model, while VG model has the lowest. On the other hand, the FOS of all models have the same result for the saturated condition. Other than that, it was found that the increasing of ESP value in the surface layer has significant effect in the sub-surface layer.     

采用间歇灌溉进行土壤盐分淋洗的适用性

为比较一次性灌溉和间歇灌溉在盐分淋洗效率方面的优劣,采用HYDRUS-1D两区模型模拟分析了不同频次灌溉条件下的土壤水盐运移过程。结果表明存在一个间歇灌溉具有更高淋洗效率的最大深度,并将其称为临界深度。对砂土的模拟显示在灌水总量为20 cm时,随着潜在蒸发强度从0 mm/d升高至6 mm/d,临界深度从60 cm以上降为0 cm;而在潜在蒸发强度为2 mm/d时,随着灌水总量从10 cm增加至40 cm,临界深度相应从20 cm增加至接近80 cm。比较而言,透水性好、不动水体占比高、两部分孔隙之间水分和溶质交换能力差的土壤,临界深度更大;对黏土、黏壤土、壤土、砂壤土和砂土分别进行模拟显示偏砂性土壤具有更大的临界深度,在灌水总量和潜在蒸发强度分别为30 cm和2 mm/d时,黏土和砂土的临界深度分别为56 cm和大于100 cm。总的来看,对于间歇灌溉是否能提高盐分淋洗效率不能一概而论,引入临界深度可以在一定程度上解释研究者们对于间歇灌溉是否能提高盐分淋洗效率的不同认识。     

Simulations of solute transport in fractured porous media using 2D percolation networks with uncorrelated hydraulic conductivity fields

Two-dimensional percolation networks have been used to model a disordered and fractured porous medium. The advantage of percolation networks is that they allow the flow and transport properties of the system to be systematically studied as a function of the connectivity of the fractures and/or the permeable regions. The aim of this research is to study hydrodynamic dispersion in such networks, and to investigate the behavior of the longitudinal dispersion coefficient D_L with binary and log-normally distributed hydraulic conductivity fields. In particular, the study focuses on the behavior of D_L at the percolation threshold p_c, where the insufficiency of flow field homogenization and the limited number of tortuous paths for flow and transport force D_L to behave anomalously, i.e., to be scale- and time-dependent. The simulations indicate that the D_L population taken over a large number of the network realizations resembles a log-normal distribution, hence indicating that, unlike the hydraulic conductivity, D_L is not a self-averaged property whose variance should tend to zero when the size of the system tends to infinity. In addition, it was found that the power law that characterizes the scale dependence of D_L is contingent upon its computation method. Moreover, D_L is found to have a completely different behavior in networks with low and high connectivities.

---

Resumen Simulaciones de transporte de solutos en medio heterogéneo utilizando redes de precolación 2D con campos de conductividad no correlacionados. Se han utilizado redes de percolación en dos dimensiones para modelizar un medio poroso fracturado y desordenado. La ventaja de redes de percolación es que permiten estudiar sistemáticamente las propiedades de flujo y transporte del sistema en función de la conectividad de las fracturas y/o regiones permeables. El objetivo de esta investigación es estudiar la dispersión hidrodinámica en las redes mencionadas, e investigar el comportamiento del coeficiente de dispersión longitudinal DL con campos de conductividad hidráulica que tienen distribución log-normal y binaria. El estudio se enfoca particularmente en el comportamiento de DL cuando se alcanza el threshold de percolación pc, donde la insuficiencia de homogenización del campo de flujo y el número limitado de trayectorias tortuosas de flujo y fuerza de transporte DL se comportan de manera anómala, i.e., son dependientes del tiempo y la escala. Las simulaciones indican que la población DL tomada en un número grande en las redes realizadas se asemeja a una distribución log-normal, indicando por lo tanto que, a diferencia de la conductividad hidráulica, DL no es una propiedad auto-promediable cuya variación tienda a cero cuando el tamaño del sistema tiende al infinito. Se encontró además que la ley potencia que caracteriza la dependencia en escala de DL es dependiente del método de cálculo. Más aún, se encontró que DL tiene un comportamiento completamente distinto en las redes de alta y baja conectividad.

---

Résumé Des réseaux de percolation bi-dimensionnels ont été utilisés pour modéliser des milieux poreux fracturés et désordonnés. Lavantage de ces réseaux est quils permettent détudier les propriétés découlement et de transport en fonction de la connexité des fractures et/ou des zones perméables. Lobjectif est dappréhender la dispersion hydrodynamique et dinvestiguer le comportement du coefficient de dispersion longitudinal D_L pour des réseaux ayant une distribution des conductivités hydrauliques binaire et log-normale. En particulier, le comportement de D_L a été étudié au seuil de percolation p_c, là où linhomogénéité des vitesses du fluide et le nombre limité de chemins tortueux disponibles pour lécoulement et le transport entraînent un comportement anormal de D_L, à savoir, une dépendance vis-à-vis des échelles de temps et despace. Les simulations montrent que les populations de D_L définies pour un grand nombre de réseaux ressemblent à des distributions log-normales, indiquant que, contrairement à la conductivité hydraulique, D_L nest pas une propriété dont la variance tend vers zéro lorsque la taille du système tend vers linfini. Il a également été trouvé que les lois de puissance qui caractérisent la dépendance déchelle de D_L découlent directement de la méthode de calcul. Enfin, les simulations engendrent un comportement très différent de D_L dans des réseaux faiblement ou fortement connectés.

     

Influence of streambed hydraulic conductivity on solute exchange with the hyporheic zone

A conservative solute tracer experiment was conducted in Indian Creek, a small urban stream in Philadelphia, Pennsylvania to investigate the role of subsurface properties on the exchange between streamwater and the hyporheic zone (subsurface surrounding the stream). Sodium Bromide (NaBr) was used as a conservative tracer, and it was monitored in the surface water at two stations and in the upper bed sediments (shallow hyporheic zone extending from 7.5 to 10 cm below the streambed). The hydraulic conductivity (K) of the upper bed sediments and the lower bed sediments (10–12.5 cm below the streambed) was measured in situ. High tracer concentrations were observed in the upper layer at locations where the hydraulic conductivity of the upper layer was larger than that of the lower layer. Low concentrations in the upper layer were observed in the converse case. A statistically significant relationship between the mass retained in the upper layer and the difference of K values between layers was observed.     

多场耦合作用下溶质运移模型及机制研究

最新研究表明,多孔介质中溶质运移现象受化学渗透和土壤表面电荷的多重影响。随着城市化进程的加快,卫生填埋场衬垫系统的长期安全性和淤泥脱水问题也愈加严重。为了综合研究水头,电场和化学渗透等物理作用对离子在土壤介质中迁移特性的影响,建立了一个新的多场耦合模型。该模型用有限元法取得数值解,并且利用移动边界成功实现了对淤泥脱水和填埋场衬垫系统中水溶液和化学溶质迁移问题的数值模拟。对于衬垫系统而言,化学渗透现象和渗透电势两者均可减缓溶质击穿土层的时间,有助于延长衬垫的工程寿命。对于淤泥脱水问题的分析结果表明,水力渗透系数和电渗系数分别是主导压滤式和电渗式脱水处理的主要变量。孔隙淤堵是压滤脱水中最主要的问题,土工织物表面所积聚的滤饼会导致水流量急剧下降     

Effects of saline water irrigation on soil properties in northwest China

Due to the lack of freshwater, highly saline groundwater was the main irrigation source in the last few decades in the Minqin Basin, which is in northwest China. The study evaluates the effects of salt accumulation on the soil physical–chemicals properties. Undisturbed and disturbed soil samples were taken from the experiment site, which was irrigated with saline water at a concentration of 0.8, 2 and 5 g L⁻¹ (coded later as C₀₈, C₂ and C₅). Undisturbed soil samples, at depths of 0–45 and 45–60 cm were taken to determine the water retention curve (WRC). Moreover, in the same place, another set of undisturbed soil samples were taken to determine the porosity and pore-size distribution (PoSD). From the WRC, the water-holding capacity of the soil was estimated. Disturbed soil samples at depths of 0–20, 0–45, 45–60 and 80–100 cm were taken to determine the index of aggregates stability in water (IC). The electrical conductivity of the saturated paste (EC_e) was determined at depths of 0–30, 30–60 and 60–90 cm, during the irrigation season on C₀₈, C₂ and C₅ treatments. The results show that the total porosity and the index of aggregates stability in water decrease with the increasing salinity of irrigation water, and the EC_e increases with the increasing salinity of irrigation water especially in the surface soil. The water-holding capacity (WHC) of soil also increases with the increasing salinity of irrigation water.     

土壤分形结构对其水力性质的指示作用

研究土壤的水力性质是进行生态水文模拟、农业水分管理和环境监测的关键,然而强烈的空间变异造成土壤的物性特征异常复杂,特别是水力性质测定困难,试验结果随机性强,耗费大量人力物力,却难以准确描述。本文在综合分析了国内外30余种量化土壤水力性质的研究方法的优劣后,系统总结了分形理论在土壤物性特征研究中的应用,剖析了分形理论与土壤水力性质之间的关系,旨在探讨土壤分形结构在水力性质领域的发展前景,以及在测定水力特性参数时所具备的优势。结果表明：(1)采用分形方法定量研究土壤结构具备可行性;(2)分形结构方法能够指示水力性质,并能为快速准确刻画不同尺度下的土壤水分分布特征提供科学依据;(3)利用已知土壤水力性质建立分形模型可以有效反演土壤结构。     

Effects of coarse-grained material on hydraulic properties and shear strength of top soil

Sidewalk failures associated with top soil of low shear strength are a common problem in urban areas. Mixing top soil with granite chips can be used to increase its permeability and shear strength. The effects of mixing granite chips with top soils on the hydraulic properties and shear strength under saturated and unsaturated conditions were investigated in this study. The results showed that the mixing top soils with granite chips caused changes in several key parameters of the soil–water characteristic curve (e.g., the air-entry value, the residual matric suction, and the residual volumetric water content) and the unsaturated permeability of the top soils. The saturated permeability and shear strength of the soil mixture increased with increasing content of granite chips.     

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.     

Effect of treated municipal wastewater on bean growth, soil chemical properties, and chemical fractions of zinc and copper

The current study was carried out in order to investigate the short-term effect of different dilutions of wastewater on soil chemical properties, chemical fractions of zinc (Zn) and copper (Cu), and to assess the chemical buildup of heavy metal on two bean species. The experiment was performed as factorial based on completely randomized design with different dilutions of wastewater on two bean cultivars in two soil textures. The treatments consisted of irrigation with treated wastewater over all growing season, irrigation with wastewater and freshwater in equal proportions, and irrigation with well water only as control. The result showed that soil parameters are significantly affected by application of wastewater irrigation. Irrigation with wastewater increased the concentrations of organic matter, electrical conductivity, N, K, P, Fe, Cu, Zn, Mn, and Ni in soils compared to the control treatment. However, their values were all below international standards. Application of wastewater decreased soil pH and calcium carbonate equivalent. Plant tissue analysis showed increases in N, P, K, Zn, Cu, Zn, and Mn concentrations in grain and frond of beans in wastewater treatment as compared to the control. The concentrations of all elements in plants were lower than the toxic threshold. Chemical fractionations of Zn and Cu indicated that chemical forms of these metals were affected by irrigation with wastewater. Irrigating with wastewater resulted in the movement of Zn from the labile fractions towards the nonlabile fractions. In turn, mobility factor of Cu increased with application of wastewater. Maximum fresh and dry yields of beans were obtained from wastewater treatment compared to the control treatment. This study indicated that wastewater irrigation improves soil properties, plant growth, and yield without any contamination in soil and toxicity in plants.     

Impact of long-term wastewater irrigation on variability of soil attributes along a landscape in semi-arid region of Iran

Twelve soil profiles and 62 soil cores were described from two neighboring landscapes in the northwest of Iran to understand whether long-term wastewater irrigation has caused variations in morphological, physicochemical, and DTPA-heavy metals contents in soils. Along a gradual slope, different positions reflected heterogeneous response to flooding wastewater irrigation. Upper slope and mainly midslope were subject to gradual soil loss as characterized by thin Ap-horizons and shallow soil depth, whereas periodic cycles of soil accumulation and sediment along with local reducing condition were highlighted in lower slope. Compared to different positions of control land, wastewater irrigation increased significantly (P????0.05) soil fertility attributes (N, P, K, and organic matter) along with electrical conductivity. Significant correlations were observed between soil fine particles fraction and available P and also with available K, indicating that spatial distribution pattern of these elements is mainly dependent on distribution of the <0.002?mm fractions along the landscape as it is affected by wastewater irrigation. Midslope position was found to be more degraded compared to the other position regarding soil loss and available nutrients (N, P, and K) as a result of the processes of soil redistribution and soil transport in wastewater irrigation system. Mean concentrations of DTPA-extractable Zn, Cu, Cd, and Pb in the different geomorphological positions were 3.8?C7.8, 2.5?C3, 3?C4, and 1.7?C3.5 times greater, respectively, when compared to the adjacent control soil. With the exception of lower slope position, which indicated a remarkable increasing pattern of smectite in response to wastewater irrigation and unfavorable drainage condition, the relative abundances of clay minerals in view of the peak position and intensity followed almost the same pattern in both wastewater-irrigated and control landscape.     

崩塌体水力输移与塌岸淤床交互影响试验

以二元结构河岸崩塌模式为研究背景,在弯道水槽中展开系列试验,研究凹岸坡脚处成型崩塌体在水力作用下输移过程及其对岸坡稳定性与河床冲淤的交互影响。试验成果表明,崩塌体经水力分解破碎呈块状、片状或颗粒状起动,部分随水流带至凸岸或下游,堆积在坡脚附近的崩塌体残留量随水力作用大小及土体特性不同而变化。崩塌体临水面周围尤其上下游端水流紊动强烈易形成较大剪切力区,临坡面上下游端附近剪切力较无崩塌体时减小;崩塌体体积越大,对剪切力区特征的改变也越明显。崩塌体的存在虽不能制止附近岸坡的再次崩塌,但可能抑制崩岸发展及附近河床淤积的程度,崩塌体的粘性或体积越大,这种抑制作用越显著;相同崩塌体抑制附近河床淤积的程度较抑制岸坡崩塌的程度大。     

裂隙岩体溶质运移模型研究

核废料的地下深部贮存、垃圾填埋造成的污水下渗、海水入侵、输油管道老化引起渗漏等诸多应用领域的发展,要求从机理、试验和模型等方面发展裂隙岩体溶质运移理论。基于裂隙岩体系统溶质运移的多尺度概念模型,分析了目前描述裂隙岩体系统中溶质运移各种数学模型的适用性和优缺点,为选取合理的数学模型求解具体的问题提供了参考依据。最后提出了需要进一步研究的问题。