土壤饱和导水率空间预测的不确定性分析

当土壤转换函数应用于土壤水力性质估计时,对于预测值的不确定性往往容易被忽视。为了有针对性地提出减少这种不确定性的方法和措施,提高土壤转换函数的实际应用能力,以两种现有的土壤转换函数(Vereecken和HYPRES模型)为例,将其应用于山东省平度市土壤饱和导水率的空间预测,并利用拉丁超立方抽样(LHS)方法对预测结果的不确定性进行了分析。结果表明,饱和导水率空间预测的不确定性主要来源于土壤基本性质的空间插值误差和土壤转换函数自身的预测误差。当Vereecken模型应用于饱和导水率空间预测时,预测结果的不确定性主要由土壤基本性质空间插值误差所决定,土壤转换函数预测误差的影响较小,而HYPRES模型则是受二者的双重影响。     

土壤饱和导水率的多尺度预测模型与转换关系

运用联合多重分形方法研究不同土层土壤饱和导水率与土壤基本物理特性的多尺度相关性,建立不同土层土壤饱和导水率的多尺度预测模型,构建不同土层土壤饱和导水率之间的转换关系。结果表明:不同土层土壤饱和导水率与土壤基本物理特性的相关程度排序不同;在单一尺度和多尺度上,0~20 cm土层土壤饱和导水率与土壤基本物理特性的相关程度排序相同,20~40 cm土层土壤饱和导水率与土壤基本物理特性的相关程度排序不同;土壤饱和导水率多尺度预测模型的预测精度较高,0~20 cm土层和20~40 cm土层拟合值的均方根误差分别为0.035 0和0.029 0;0~20 cm土层和20~40 cm土层土壤饱和导水率转换关系的计算精度较高,拟合值的均方根误差为0.037 5。     

土壤非饱和导水率模型中参数的敏感性分析

针对邵明安根据再分布过程得出的非饱和导水率模型,对其中参数的敏感性进行分析,通过实测资料计算非饱和导水率对不同参数的敏感度,对比参数对非饱和导水率的影响程度。研究表明采用线性关系表示土壤湿润剖面的平均湿度和湿润锋处湿度之间的关系时,非饱和导水率对参数的敏感性比其他参数高得多。     

基于PCA-GWR的包气带不同土层饱和导水率的传递函数及其回归克里金估计

饱和导水率是重要的土壤水力学参数,对渗流和溶质运移研究起着至关重要的作用.土壤传递函数可以代替大规模采样和室内外试验间接预测该参数,但由于土壤的空间变异性以及方法的局限性等原因,以往的传递函数往往精度有限.本文以黄河下游沿岸兰考县闫楼乡包气带不同土层为研究对象,基于64个钻孔数据,将土壤粒径分布（黏粒、粉粒和砂粒含量）、土壤结构分形维数、干容重、总孔隙度、pH值、有机质和电导率等9个基本理化特性参数作为影响因素,基于逐步回归、主成分回归和主成分-地理加权回归（PCA-GWR）三种方法,分别对研究区包气带不同土层的饱和导水率进行回归分析,比较精度后基于PCA-GWR对饱和导水率进行地理加权回归克里金插值.结果显示,除表层粉土以外,PCA-GWR法的预测精度有显著优势.不同土层传递函数和饱和导水率分布的差异性表明,由于黄河下游地区快速多变的沉积过程、黄河的决口改道以及人类活动等因素,包气带土层相变剧烈,土壤饱和导水率在平面上亦具有明显的非平稳特征.而局部地质作用和人类活动是第二层粉土层回归克里金估计结果欠佳的深层原因.     

大孔隙对农田耕作层饱和导水率的影响

介绍了土壤中大孔隙流的研究,并用一种简单的方法对土壤内的大孔隙和中小孔隙的饱和导水率进行分析,同时讨论了农田耕作层大孔隙流和Darcy流饱和导水率的变异情况。     

包气带增厚区土壤水力参数及其对入渗补给的影响

为探讨包气带深部增厚区土壤水力参数变化对入渗补给过程的影响,采用压力膜仪对河北正定深部包气带(8.0~21.0 m)10个原状土样进行水分特征曲线测试,利用RETC软件中Mualem-van Genuchten导水率模型对其拟合,获取含水率与非饱和导水率的关系曲线,并根据达西法对其进行分析讨论.结果表明:场地包气带深埋区的非饱和导水率为25~240 mm/a.当某一埋深历史水位下降速度越快,该埋深处相同含水率情况下土壤非饱和导水率越大,说明对应土层的入渗补给强度越大;因包气带厚度增大使原来位于饱水带的层状非均质土层转变为包气带,潜水位波动下降过程中深部包气带土层因排水压密作用,使得土壤水力特性发生变化,进而影响垂向入渗补给过程.      

多孔介质非饱和导水率预测的分形模型

多孔介质非饱和导水率是地下水污染预测与评价的重要参数。根据分形几何的基本原理和方法,推导出了与Campbell经验公式在形式上完全一致的多孔介质非饱和导水率的预测公式。公式中的幂指数为介质孔隙分维和随机行走分维的函数,分别体现了多孔介质的静态性质与动态性质对其中水分运动的影响,但静态性质的影响是主要的,即导水率主要受多孔介质的结构控制。根据文献中报道的大量数据,利用笔者推导的预测公式计算得到的幂指数的统计值与试验测定的幂指数的统计值基本一致,说明推导的理论公式预测多孔介质非饱和导水率是较为可靠的。     

邯郸地区土壤水入渗规律实验及数值模拟方法研究

在对土壤入渗理论进行总结和分析的基础上,依据邯郸市平原区的水文水情的基本特点,选定2~3处地点进行实验,完成了静态数据的收集和动态数据的观测,确定了实验区土壤质地类型、土壤容重等物理特性,结合实验测定的饱和导水率K、土壤水分特征曲线θ（ h）推算非饱和导水率K（θ）、扩散率 D（θ）等土壤水分运动参数。运用上述数据对Van Genuchten模型加以改进,并运用改进后的模型对实验中的水分变化过程进行数值模拟。结合实测成果和数值模拟成果,对邯郸市平原区已选定实验点土壤水入渗规律进行分析。     

应用连续混合器模型模拟非饱和土壤盐分淋洗过程

尽管连续混合器模型已经广泛应用于土壤中溶质迁移转化规律计算,但却仅限于饱和土壤条件。通过采用Kostiakov公式描述非稳定入渗过程,并结合活塞假设计算入渗过程中土壤水分的运动和分布规律,将连续混合器模型扩展应用于非饱和土壤的盐分淋洗过程计算,并求出了模型的解析解。分析发现,在一定的土壤质地、初始含水率剖面和初始含盐量剖面下,模型的计算结果会受到土层划分厚度的影响,而合理的土层厚度又与饱和导水率和溶质弥散系数之间存在量纲一的函数关系。当土层厚度选取合理时,连续混合器模型与HYDRUS-1D的计算结果十分接近,采用该模型模拟试验土壤盐分淋洗过程,显示了较好的计算精度。     

基于大地电磁法的水文地质结构及参数推断

大地电磁法已被广泛用于矿产资源勘探和工程勘查,但该方法在水文地质参数定量方面的成果文献尚不多见。以科尔沁沙丘—草甸相间地区巨厚第四系沉积层为研究对象,通过EH-4大地电磁法的测试结合地质钻探试验数据,划定出研究区埋深105 m范围内砂黏土层在空间上的分布;经检验大地电磁法探测总体误差在5.1%~7.2%,其精度足以满足水文地质勘察的需求。然后利用研究区已有的土壤传递函数(Pedo-transfer functions,PTFs)模型对测区埋深105 m范围内的饱和导水系数进行了估算,进而首次利用大地电磁法所得数据建立了视电阻率与饱和导水系数的拟合模型。利用该模型推断出研究区埋深200 m范围内含水层的富水性,刻画出饱和导水系数的空间分布,为研究区提供了更加丰富的水文地质信息,同时也极大地提升了水文地球物理方法在水力参数预测方面的尺度。     

随机单裂隙饱和/非饱和渗透系数研究

裂隙易在降雨作用下诱发滑坡等灾害,裂隙的饱和/非饱和渗透特性是研究此类问题的关键。利用精密数控机床制作随机粗糙裂隙面,并研制了一套仪器进行此随机粗糙裂隙的渗流试验,得到了裂隙的饱和渗透系数,然后通过间接方法预测此裂隙的非饱和渗透系数。研究发现,当裂隙平均开度为0.4 mm时,其饱和渗透系数为0.1 m/s。通过立方定律得到的水力等效隙宽为0.35 mm,小于其平均隙宽。同时裂隙的渗透系数与平均隙宽的平方成正比,这与立方定律的趋势相一致。研究得到了不同隙宽裂隙的非饱和渗透系数函数,当基质吸力小于进气值时,渗透系数为一常数,即为饱和渗透系数;当基质吸力大于进气值时,裂隙板的渗透系数急剧减小。当裂隙板的基质吸力达到其残余含水率对应的吸力值时,裂隙板的渗透系数基本稳定。在此情况下,基质吸力的继续增加对渗透系数的影响非常小,很难使渗透系数减小。     

Influence of Unsaturated Soil Properties Uncertainty on Moisture Flow Modeling

The quality of a numerical modeling solution of moisture flow through unsaturated soil, in part, depends on properly described unsaturated soil properties. The variability of the Soil Water Characteristic Curve, SWCC, is attributed to hysteresis and reproducibility of measurement. Because the unsaturated conductivity function is rarely directly measured, the variability of the unsaturated soil hydraulic conductivity function is attributed to the uncertainty associated with the estimation of this parameter with currently available fitting functions, and hence a range of reasonable variation was considered. One-dimensional modeling of expansive soil under dry initial conditions (suction of 1,500 kPa) was performed; both potential evaporation and infiltration boundary conditions were considered. It was found that small variations in the unsaturated soil hydraulic conductivity function result in significantly different modeling outputs, as expected, while substential variation in SWCC alone (assuming the same unsaturated soil hydraulic conductivity for all SWCCs) produced almost identical soil response in terms of soil suction when the slope of the SWCC is similiar. Thus, proper characterization of the slope of the SWCC is important to proper suction profile determination.     

Experimental and analytical model studies on leachate volume computation from solid waste

An extensive use of solid-waste landfills for disposal of municipal and industrial wastes have prompted increased attention to groundwater pollution caused by leachate generated in such landfills. The potential for groundwater contamination by leachate has necessitated engineering designs for landfills. The quantity of leachate generated from the solid waste and the movement of water through the solid waste depends on water input and the solid-waste characteristics. This paper dealt with the experimental investigations using the laboratory solid-waste leaching column to estimate the total leachate volume/leachate flow for unsaturated and saturated conditions. The hydraulic properties of the solid waste like initial moisture content, field capacity, permanent wilting point, saturation moisture content, effective void ratio, saturation hydraulic conductivity and saturation suction pressure were determined from the small-scale laboratory experiments, which are the input for analytical model study of leachate flow/total leachate volume for both unsaturated and saturated conditions. The result of analytical model study was compared with the results of experimental investigations. Comparisons of measured and computed total leachate volume/leachate flow using Darcy’s law showed reasonable agreement.     

晚更新世黄土渗透性的各向异性及其机制研究

非饱和黄土的渗透性是非饱和黄土性质的重要组成部分。研究黄土不同方向的渗透性对确定其湿陷范围和由于水的渗透引起的黄土滑坡具有很重要的理论意义。研究了黄土渗透性的各向异性特征及其机制。以具有明显各向异性的西安Q3原状黄土为研究对象,用TEN型张力计测量了黄土试样不同方向的、不同含水率下的吸力;用变水头渗透试验测量了黄土竖直和水平方向的饱和渗透系数。结果表明,当体积含水率在23%～41%时,张力计沿不同方向插入土样所测吸力相差不大;竖直方向的饱和渗透系数是水平方向的4.02倍。在吸力测量的基础上,根据土-水特征曲线,确定了竖直和水平方向的非饱和黄土的渗透系数。得出在黄土不同方向,随着吸力的增大或减小,渗透系数减小或增大;竖直方向的渗透系数普遍地大于水平方向的渗透系数;当吸力小于57 kPa时,随着吸力的增大,竖向渗透系数与水平向渗透系数的差值减小。通过观测黄土的结构,得出黄土结构对其渗透性有重要影响。     

垃圾土渗透性和持水性的试验研究

通过室内常水头试验测定了垃圾土的饱和渗透系数,通过室内压力板仪试验测得垃圾土的土-水特征曲线。基于土-水特征曲线预测了垃圾土的非饱和渗透系数,并通过室内入渗试验进行了初步验证。常水头试验得到深层、中层和浅层垃圾土的饱和渗透系数分别为3.56×10-4、3.50×10-3、4.81×10-2 cm/s。土-水特征曲线试验表明,垃圾土饱和含水率和残余含水率较高,进气值很小,土-水特征曲线在低基质吸力时存在陡降段,其中浅层垃圾土的陡于深层垃圾。验证试验表明,预测得到的非饱和渗透系数与实测结果接近,基于土-水特征曲线预测垃圾土非饱和渗透系数的方法基本可行。     

Predicting saturated hydraulic conductivity using particle swarm optimization and genetic algorithm

Soil saturated hydraulic conductivity (K_s) is considered as soil basic hydraulic property, and its precision estimation is a key element in modeling water flow and solute transport processes both in the saturated and vadose zones. Although some predictive methods (e.g., pedotransfer functions, PTFs) have been proposed to indirectly predict K_s, the accuracy of these methods still needs to be improved. In this study, some easily available soil properties (e.g., particle size distribution, organic carbon, calcium carbonate content, electrical conductivity, and soil bulk density) are employed as input variables to predict K_s using a fuzzy inference system (FIS) trained by two different optimization techniques: particle swarm optimization (PSO) and genetic algorithm (GA). To verify the derived FIS, 113 soil samples were taken, and their required physical properties were measured (113 sample points?×?7 factors?=?791 input data). The initial FIS is compared with two methods: FIS trained by PSO (PSO-FIS) and FIS trained by GA (GA-FIS). Based on experimental results, all three methods are compared according to some evaluation criteria including correlation coefficient (r), modeling efficiency (EF), coefficient of determination (CD), root mean square error (RMSE), and maximum error (ME) statistics. The results showed that the PSO-FIS model achieved a higher level of modeling efficiency and coefficient of determination (R²) in comparison with the initial FIS and the GA-FIS model. EF and R² values obtained by the developed PSO-FIS model were 0.69 and 0.72, whereas they were 0.63 and 0.54 for the GA-FIS model. Moreover, the results of ME and RMSE indices showed that the PSO-FIS model can estimate soil saturated hydraulic conductivity more accurate than the GA-FIS model with ME?=?10.4 versus 11.5 and RMSE?=?5.2 versus 5.5 for PSO-FIS and GA-FIS, respectively.     

考虑参数空间变异性的非饱和土坡可靠度分析

在考虑多个土体参数空间变异性的基础上,提出了基于拉丁超立方抽样的非饱和土坡稳定可靠度分析的非侵入式随机有限元法。利用Hermite随机多项式展开拟合边坡安全系数与输入参数间的隐式函数关系,采用拉丁超立方抽样技术产生输入参数样本点,通过Karhunen-Loève展开方法离散土体渗透系数、有效黏聚力和内摩擦角随机场,并编写了计算程序NISFEM-KL-LHS。研究了该方法在稳定渗流条件下非饱和土坡可靠度分析中的应用。结果表明：非侵入式随机有限元法为考虑多个土体参数空间变异性的非饱和土坡可靠度问题提供了一种有效的分析工具。土体渗透系数空间变异性和坡面降雨强度对边坡地下水位和最危险滑动面位置均有明显的影响。当降雨强度与饱和渗透系数的比值大于0.01时,边坡失效概率急剧增加。当土体参数变异性或者参数间负相关性较大时,忽略土体参数空间变异性会明显高估边坡失效概率。     

快速测量非饱和土渗透系数的湿润锋前进法适用性研究

非饱和土渗透系数函数是非饱和土体渗流分析中必不可少的参数。然而非饱和土渗透系数函数测量非常困难,其中稳态法测量范围过于狭小,而瞬时截面法耗时动辄数月,难以应用。为此,提出了湿润锋前进法,有望在一周内测得渗透系数函数。采用数值模拟手段检验湿润锋前进法的测量精度和适用范围。数值模拟结果显示：对于砂土、粉土、黏土,湿润锋前进法得到的渗透系数函数均与输入的理论渗透曲线一致,但测量范围限于高吸力范围（大于进气值）。此外,还采用瞬时截面法计算了渗透系数函数,结果表明：瞬时截面法计算的渗透系数在理论渗透曲线附近波动,其精度取决于监测截面的间距,间距越大,精度越低。在给出的20、60、100 mm这3种间距中,20 mm间距时计算得到的渗透系数精度最高,但与湿润锋前进法的精度仍有一定差距。因此,湿润锋前进法是一种可靠的能够快速得到渗透系数函数的测量手段。     

Unsaturated hydraulic properties of cemented tailings backfill that contains sodium silicate

Recycling the mine waste (tailings) into cemented tailings backfill has economical and environmental advantages for the mining industry. One of the most recent types of cemented tailings backfill is gelfill (GF), a backfill that contains sodium silicate as chemical additive. GF is typically made of tailings, water, binder and chemical additives (sodium silicate gel). It is a promising mine tailings backfill technology. From a design point of view, the environmental performance or durability of GF structures is considered as a key factor. Due to the fact that GF structures are cementitious tailings, their durability and environmental performance depend on their ability to resist the flow of aggressive elements (water and oxygen). Thus, understanding the unsaturated hydraulic properties of GF is essential for a cost-effective, environmentally friendly and durable design of GF structures. However, there is a lack of information with regards to unsaturated hydraulic properties of GF, the factors that affect them and their evolution with time. Hence, the unsaturated hydraulic properties (water retention curve (WRC) or water characteristic curve, air entry value (AEV), residual water content, unsaturated hydraulic conductivity) of GF are investigated in this paper. GF samples of various compositions and cured in room temperature for different times (3, 7, 28, and 90 days) are considered. Saturated hydraulic conductivity and microstructural tests have been conducted; WRCs are measured by using a WP4-T dewpoint potentiameter and the saline solution method. Unsaturated hydraulic conductivity is predicted using the van Genuchten (1980) equation. The water retention curve (WRC) is determined as the relationship between volumetric water content and suction for each GF mix and curing time. The van Genuchten (1980) equation is used to simulate the WRC to best-fit the experimental data. AEV and residual water content are also computed for each mix and curing time. Furthermore, functions are developed to predict the evolution of AEVs, residual water content and fitting parameters of the van Genuchten model with degree of hydration. Important outcomes have been achieved with regards to unsaturated hydraulic properties. The unsaturated hydraulic conductivity of GF was calculated to decrease when the suction, binder content, and degree of hydration increase. The effects of binder content and degree of hydration are more obvious at low suction ranges. The obtained results would contribute to a better design and assessment of the durability and environmental performance of GF structures.