利用畦灌试验资料推求土壤入渗参数的改进方法.docxVIP

利用畦灌试验资料推求土壤入渗参数的改进方法 土壤渗透参数是确定土壤灌溉制度和评价土壤侵蚀质量的重要指标。土壤入渗参数的现场测定多采用筒测法。该方法需要的仪器设备比较简单,但试验耗时长,在田间土壤存在空间变异时,需进行多个试验,求其平均值来代表田间土壤的总体入渗性能。采用灌水资料推求入渗参数,可以避免筒测法的上述缺点,只需对一个或多个灌水过程进行完整的观测,记录田面水流推进过程,观测一个或多个点(视方法而定)的水深变化,通过计算便可得到满意的结果,因此受很多学者的推崇。利用灌水资料推求入渗参数的方法有一点法、两点法、Maheshwari法、Esfandiari法等,其中后两者的计算精度较高,因可计算Kostiakov、Philip等入渗模型而被广泛应用。但这两种方法都存在计算工作量大的缺点。为此,缴锡云根据Kostiakov入渗透模型的特点,引入了两点法的计算思路,对这两种方法进行了改进(本文分别称改进Maheshwar1法1和改进Esfandiari法1),减少了计算工作量,但改进后的方法只适用于Kostiakov入渗模型。本文根据水量平衡原理,提出了将Maheshwari法和Esfandiari法利用模式搜索技术进行参数计算改为利用非线性回归法进行多次逼近计算(分别称Maheshwarl改进法2和Esfandiari改进法2),使计算工作量大大减少。利用实测资料进行验证,改进后的方法具有更高的计算精度。 1 计算原理 1.1 中心管网流量计算 在畦灌过程中,水流推进到某一距离时的地表水面线与湿润范围如图1所示。 图中,x为水流前锋推进距离,L;h为地表水深,L;Z为累积入渗水量,以水深表示,L;Vh为水流前锋推进至X处时,单位畦宽上的地表储水总量,L2;VZ为对应的单位畦宽上的下渗水总量,L2;s为沿畦长方向距畦首的距离,L。 根据水量平衡原理,在水流推进过程中,对于单位畦宽,有 qt=Vh+VΖ(1)qt=Vh+VZ(1) 式中:q为畦灌单宽流量,L2T-1;t为灌水时间,与畦首入渗时间相同,T。 1.2 地表内风形系数 由图1可知,当水流前锋到达X处时,单位畦宽上的地表总储水量为 Vh=∫?x?0hds(2)Vh=∫?x?0hds(2) 式中:Vh单位畦宽上的地表储水总量,L2;s为沿畦长方向距畦首的距离,L;x为水流前锋的推进距离,L;h为地表水深,L。 在沿畦长设多个水深观测点的情况下,式(2)的近似计算式为 Vh≈m-1∑i=1(hi+hi+12)(xi+1-xi)(3) 式中:xi,xi+1分别为第i,i+1个水深观测点距畦首的距离,L;hi,hi+1为对应于xi,xi+1的田面水深,L;m为田面水深观测点个数。 在已知畦首水深h1情况下,引入地表储水形状系数σh=Vhh1x,式(2)可表示为 Vh=σhh1x(4) 式中:h1为畦首处水深,L;σh为地表储水形状系数,其它符号意义同前。 地表储水形状系数σh随水流推进过程而变化,一般取值在0.7～0.8之间,Maheshwari、缴锡云建议取恒定值0.75。 1.3 地表流推进模型及入渗参数估算方法 Ζ=Κτα(5) 式中:Z为单位畦田面积上的累积入渗量,以水深表示,L;τ为入渗历时,T;α为入渗指数,无因次;K为入渗系数,LT-α,如果不考虑灌水引起土壤质地等的变化,则K、α在一次灌水过程中应保持不变。 根据图1,当水流前锋到达X处时,单位畦宽上的总下渗水量为 Vz=∫?x?0Ζds(6) 式中:VZ单位畦宽上的渗水总量,L2;s为沿畦长方向距畦首的距离,L;x为水流前锋的推进距离,L;Z为沿畦长各点的入渗水深,L。 VZ的近似计算式为 {VΖ≈l-2∑i=1Κ((t-τi)α+(t-τi+1)α2)(xi+1-xi)+Κ(t-τl-1)α1+α(xl-xl-1)水流推进阶段(2≤l≤m)VΖ≈m-1∑i=1Κ((t-τi)α+(t-τi+1)α2)(xi+1-xi)水流消退阶段(7) 式中:xi,xi+1分别为第i,i+1个水深观测点距畦首的距离,L;τi,τi+1分别为水流推进至第i,i+1个水深观测点的时间,从开始灌水计时,T;t为观测时间,从开始灌水计时,T;l为观测时间t时水流前锋所在观测点的序号,l≤m;其它符号意义同前。 上式中11+α是水流前锋段(xl～xl-1段)总入渗水量计算的权重系数(即形状系数),是在假定水流前锋从距离xl-1推进至xl期间,推进距离x与时间t呈直线关系(即r=1)的情况下,通过积分推导而来,也可从式(10)直接得到。 引入下渗水形状系数σz=VzΖ1x,则利用畦首入渗水深求单位畦宽上下渗水总量的计算式为 Vz=σΖΖ1x=σΖΚtαx(8) 式中:Z1为畦首入渗水深,L;t为畦首入渗历时,T;其它符号意义同前。 对于