考虑淋洗系数和优先流系数的土壤储水流通量计算.docxVIP

考虑淋洗系数和优先流系数的土壤储水流通量计算 灌溉。特别是由于优质淡水资源日趋紧张,在利用微咸水、咸水和污水灌溉的条件下,这些劣质水对土壤,特别是上部主要耕作层水盐动态的影响显得更加重要。为此,本文本着简单适用宜于推广的原则,从农田实际出发,提出了灌溉水入渗条件下的农田土壤水盐动态简化模型。 1 在灌溉条件下，耕地中的土壤水盐动态简化模型 1.1 计算模型假设 由于虫洞、根孔和干缩裂缝等大孔隙结构的发育,田间土壤一般不会发生均匀渗透,部分入渗水会较迅速地直接从这些较大孔隙中漏走,产生所谓的优先流。因此,在建立入渗水流模型时,首先引入优先流系数f1,表示优先流所占的比例系数。对于某一土层,设入渗率为i,则(1-f1)i发生均匀渗透,f1i为不均匀的大孔隙流即优先流。对于f1i部分快速水流,本文忽略其淋盐作用。 有关土壤表面入渗规律的经验、半经验半理论公式较多也比较成熟,但相应情况下该土层下边界水分渗漏规律的同类公式却很少。因此,对于(1-f1)i均匀入渗渗透水流,我们假设土层下边界水流通量q与土层储水量W呈线性关系,即: q(1-f1)i=w-w0wS-w0(1)q(1?f1)i=w?w0wS?w0(1) 又 dwdt=(1-f1)i-q(2)dwdt=(1?f1)i?q(2) 式中,WS—计算土层饱和储水量,mm;W0—计算土层初始储水量即灌前一定土层的储水量,mm;q—土层下边界水流通量,mm?s-1;i—上边界入渗率,mm?s-1;t—时间,s;f1—优先流系数,无量纲。当W=W0时,q=0;当W=WS时,q=(1-f1)i。 公式(1)表示土壤未饱和时地面入渗率i大于土层下部的渗漏率q,随着土壤含水量W的增加,q增大,当土壤饱和时q=(1-f1)i。为了计算土层下边界渗漏率q,本文提出了该线性假设。通过室内土柱和田间灌溉试验验证,该公式基本符合实际。 因为人们并不关心入渗过程中水盐的变化,主要是关心入渗后的结果。所以,设i为常量,则解上述微分方程,有: q={1-exp[-t(1-f1)iwS-w0]}(1-f1)i(3)w=wS-(wS-w0)?exp{-[t(1-f1)iws-w0]}(4) 设灌水量I全部入渗,则 w=wS-(wS-w0)?exp[-(1-f1)ΙwS-w0](5)Δw=(wS-w0)[1-exp(-(1-f1)ΙwS-w0)](6)dip=Ι-Δw(7)wS=10?γ?D?θS(8) 式中,W—入渗后土壤储水量,mm;ΔW—入渗后土壤储水量变化,mm;I—灌水量或降雨量,mm;dip—下边界渗漏量,mm;D—土层厚度,m;θS—计算土层饱和含水量,占干土重%;γ—土壤容量,T?m-3;其他意义同前。 1.2 土壤溶液初始盐量k 在入渗水流中,由于其淋洗作用,下边界渗漏水溶液浓度是灌溉水矿化度和原土壤溶液浓度的组合。为此,引入淋洗系数f2表示渗漏水溶液浓度中原土壤溶液浓度的贡献比例。设土层下边界渗漏量的溶液浓度: CΡ=(1-f2)CΙ+f2Ct(9)Ct=S0W0(10) 式中:f2—淋洗系数,无量纲;S0—计算土层初始盐量,g;CI—灌溉水矿化度,g?L-1;Ct—计算土层土壤溶液浓度,g?L-1。 其盐量变化公式如下: ΔS=(1-f1)ΙCΙ-[(1-f1)Ι-ΔW]?[(1-f2)CΙ+f2Ct](11)S=aS0+bΙCΙ(12)a=1+Δww0f2-Ι(1-f1)w0f2(13)b=f2+Δw(1-f2)Ι(1-f1)(14)Sd=ΙCΙ+S0-S(15) 式中:S—入渗后计算土层盐量,g;ΔW由式(6)计算;Sd—计算土层下边界漏盐量,g;其他意义同前。 1.3 模型因素的物理意义的讨论 1.3.1 tfdf11 优先流系数f1表示入渗时优先流所占的比例系数,数值在0～1之间变化,是无因次的比例系数。由于土壤上层疏松,大孔隙结构发育,随着土层深度的加厚,土壤结构趋于均匀。因此,f1一般随着土层的加厚而减少。公式(6)中: WS=10·D·γ·θS;W0=10·D·γ·θ0;W=10·D·γ·θ 可反求得: f1=1+10?γ?(θS-θ0)?DΙLn(1-θ-θ0θS-θ0)(16) 令β=-10?γ?(θS-θ0)Ln(1-θ-θ0θS-θ0),则f1=1-βDΙ;1-f1=βDΙ; 上式1-f1表示了土层的储水特性,并随土层厚度及灌水量而变化。而β值相对稳定,可代表土壤储水特性,β值越大,储水能力越大。 1.3.2 土壤溶液浓度ct和盐盐水构成ss 淋洗系数受许多因素影响,由于较多的灌溉水将通过较大的孔隙迅速地漏走,来不及溶解更多的盐分,故具有许多大孔隙的土壤将比无结构土壤的f2值低。由于大孔隙的数量随深度而减小,故浅根区将比深根区的f2值低。此外,灌水方法也影响淋洗系数,实际