具有砂质夹层土壤入渗过程的稳定性分析.docxVIP

具有砂质夹层土壤入渗过程的稳定性分析 根据研究，在混合土壤中的水合花过程中，随着时间的推移，随时间的推移，积累的入渗量从非线性入渗量逐渐转移到线性稳定入渗量。在这项入渗过程中，作者通常使用到达稳定入渗阶段时间p1的权限，并对其进行分段计算。这样的处理尽管可以满足工程设计计算的要求,但仍未能完全反映出客观实际,且没有合理的物理解释,同时由于表达两个阶段的函数在分界点t1时刻的不连续性,也给后期的进一步分析带来诸多不便。为了改进这一不足,本文在一个新的思路基础上,对具有砂质夹层土壤的入渗过程,建立一个完整的连续函数模型,该模型经实测资料验证,不仅更能符合具有砂质夹层土壤的入渗客观规律,而且具有较好的计算精度。 1 模型的构建 1.1 稳定入渗阶段的影响 关于在具有砂质夹层土壤的积水入渗过程中,其累积入渗量随时间的变化存在着由非线性入渗关系逐渐转入线性稳定入渗的机理,作者在文献中对此做了详尽的论述。研究分析表明,在入渗初期,当湿润锋尚未到达砂质夹层之前,其入渗过程实际上是一个均质土积水自由入渗过程,而当入渗锋面到达砂层上界面之后,由于砂层的阻水作用,水流将不再继续下渗而在上层土体内聚集,使上层土体剖面含水量逐渐增加,并趋于饱和,其累积入渗量F随着土壤水各势能梯度的逐渐减小表现为一个非线性的变化关系。此后在出现砂层上界面的土壤含水量接近饱和而使其能量水平逐渐大于砂质夹层水势时,入渗水流方能继续向砂层中流动,整个积水入渗过程开始进入一个线性稳定入渗阶段,此时对于砂层以上的均质土来讲,可视为近似饱和的土体的渗透问题。入渗全过程的累积入渗量F随时间T的变化关系如图1中曲线3所示,其中t1为开始进入稳定入渗阶段的时间。 基于以上分析,如果从整个入渗过程进入稳定入渗阶段后的角度来看待砂层以上土体的非线性入渗问题,则可将此阶段的累积入渗量F随时间T的变化假定为两部分入渗过程的叠加,即一部分为稳定入渗过程,其入渗率与进入稳定阶段后的量值fp相等,对于表示这一部分的累积入渗量Fi与时间T的关系在图1中是一条斜率与稳定入渗阶段相同且过零点的直线,见图1中曲线1;而另一部分可视为使上层土体水分逐渐增加接近饱和状态的入渗过程,而这一部分的入渗率则随时间T的增大而减小,直至到达稳定入渗时间t1时趋于零,对于表示这一部分的累积入渗量Fs与时间T的变化关系,在图1中是以Fb为渐近线的一条曲线,如图1中曲线2所示。根据这一思路,即可将表达前一部分入渗过程的函数称为稳渗函数Fi,而将表达后一部分入渗过程的函数称为储水函数Fs,为此,对于具有砂质夹层的土壤水分入渗模型将具有以下基本形式: F=Fi(T)+Ts(T)(1)F=Fi(Τ)+Τs(Τ)(1) 1.2 连续函数入渗模型的构建 1.2.1 入渗水流入渗率与水层稳定性的关系 由图1可见,稳渗函数Fi为一线性函数,它直接与入渗水流进入砂层的稳定入渗率fp有关,即: Fi=fpT+C(2)Fi=fpΤ+C(2) 式中:C——常数,且因该直线过零点,故C=0。 1.2.2 储水函数的渐进值 同样由图1 可见,储水函数Fs为一非线性函数,由上述思路可知,反映储水要求的这一部分累积入渗量Fs将随入渗时间T的增加而增大,而其入渗率将逐渐减小趋于零,并以进入稳定入渗时间t1时的入渗量Fb为渐近线,为此可采用以下函数表示这一变化关系: Fs=Fb(1?e?T/β)(3)Fs=Fb(1-e-Τ/β)(3) 式中:Fb为该储水函数Fs的渐进值;β为时间常数,一般可取达到稳定入渗时间t1的1/3,即β=t1/3。 如果将式(2)与式(3)代入式(1),则连续函数入渗模型即为: F=fpT+Fb(1?e?T/β)(4)F=fpΤ+Fb(1-e-Τ/β)(4) 1.3 模型中连续函数入渗的参数的分析和确定 由上式(4)可见,连续函数入渗模型中的参数,主要有稳定入渗率fp与储水函数Fs的渐近值Fb,两个参数的确定有以下两种方法。 1.3.1 储水系数的变化 该方法的前提是需要有砂质夹层土壤积水入渗的试验资料,在此基础上可以绘制出累积入渗量F随入渗时间T的变化过程线(参见图1中曲线3),则由图中曲线变化规律明显可见,稳定入渗阶段的直线斜率即为稳渗率fp,而该直线在纵坐标(T=0)上的截距即为储水系数Fb。因此该方法除依据试验资料采用图解法可取得参数fp与Fb值外,也可以根据入渗试验资料将稳渗阶段的数据进行线性拟合,即可通过计算得出。 1.3.2 砂层中sm—利用土壤基本资料进行计算 该方法的前提是需要有土、砂的基本物理参数及上层均质土壤的入渗参数,在此基础上即可通过以下计算得出fp与Fb值。 (1) 稳渗率fp的计算:根据作者在文献中提出的方法,稳渗率fp可由下式求出: fp=K′(θs)(H+Sm+Z)/Z(5)fp=Κ′(θs)(Η+S