基于改进maheshwari法的地表储水形状系数计算.docxVIP

基于改进maheshwari法的地表储水形状系数计算 在设计、评价和管理灌溉地面灌溉系统时，土壤入渗特性是一个非常重要的因素。根据maheshwari等人的研究，影响耕地灌溉水流促进的不同因素中，入渗特性的重要性仅为耕地的入渗量。因此，为了准确估计模型的入渗参数，选择合适的入渗模型是土壤灌溉技术的一个重要问题。 测定田间土壤入渗参数的专用设备有圆筒测渗计等.然而,由于土壤空间变异性的存在,在利用这些专用设备测定田间土壤入渗参数时,试验点数必须很多,否则难以取得足够的精度.此外,这些方法也不适用于某些特殊情况,例如在膜孔灌溉中,由于渗水界面是一个由多孔塑料薄膜覆盖的非均质渗水界面,其入渗性能与垂直一维入渗有着较大的差异,所以利用圆筒测渗计等专用设备在露地上所测的入渗参数不能被直接应用.近年来,国内外学者提出了一些利用灌溉水流的推进及消退过程、地表水深等资料来估算入渗参数的方法. Ellion\,Walker对Kostiakov\_Lewis入渗模型提出了估算入渗参数的两点法.该方法仅需要分别观测水流前锋推进到沟长中点、沟末端的时间及相应的沟首地表水深,便可计算入渗参数.Shepard等针对沟灌提出利用水流前锋推进到沟末端的时间和沟中的平均过水面积来估算Philip模型入渗参数的一点法.一点法、两点法的观测、计算简单,但由于观测的数据量较少所以精度较低,并且一点法只能对Philip入渗模型进行参数估算,适用性受到了限制. 王文焰等在研究波涌灌溉理论时,对Kostiakov入渗模型提出了利用两个畦田的水流推进和消退过程资料来估算入渗参数的方法.费良军在研究浑水波涌灌溉理论时,对Kostiakov入渗模型提出了通过测量两个畦田上沿畦长若干点的地表水深来估算入渗参数的方法.王文焰法、费良军法,必须要用两个以上畦田的试验数据来计算,精度较高,但试验工作量较大. Esfandiari等提出了利用水流推进过程和沿沟长上若干点地表水深的变化过程,根据水量平衡原理,采用模式搜索技术来估算沟灌入渗参数的方法.由于该方法采用类似于试算的模式搜索法进行计算,所以计算量大,计算效率较低.缴锡云等对Esfandiari法在数据处理方法上进行了改进,变原来的模式搜索法为直接计算. Maheshwari等利用水流推进过程和畦首地表水深的变化过程,引用地表储水形状系数并取其值为0.75,采用类似于试算法的模式搜索技术来估算畦灌入渗参数.该方法可以仅通过一个畦田的水流推进过程、畦首地表水深等试验数据来估算入渗参数,试验工作量较小,但由于该方法采用类似于试算的模式搜索法进行计算,所以计算量大,计算效率较低. 本文根据水量平衡方程,针对Kostiakov入渗模型,通过数学推导,对Maheshwari法加以改进,得到了利用水流推进过程、畦首的地表水深等观测资料来估算畦灌入渗参数的简捷计算方法. 1 单次河流生长方向上距两自关系 对于畦灌,在水流的推进过程中,地表储水和土壤入渗水的轮廓如图1所示.图中,0为畦首;s为沿畦长方向上距畦首的距离,L;x为水流前锋的推进距离,L;h为地表水深,L;Z为累积入渗水量,以水深表示,L.按照一般规律,畦首的地表水深和入渗水量最大,水流前锋处的地表水深和入渗水量为0,沿畦长方向愈接近水流前锋地表水深和入渗水量愈小. 1.1 地表下储水形状系数的确定 Z=Ktα(1)Ζ=Κtα(1) 式中:Z为单位畦田面积上的累积入渗量,以水深表示,L;t为入渗历时,T;α为入渗指数,无因次;K为入渗系数,LT-α.α、K统称为入渗参数.一般认为,在某一次连续灌水的整个过程中,入渗参数是稳定不变的. 由图1可以看出,当水流前锋推进到x处时,单位畦宽上的渗水总量为 VZ=∫x0Zds(2)VΖ=∫0xΖds(2) 式中:VZ为单位畦宽上的渗水总量,L2,s为沿畦长方向距畦首的距离,L;x为水流前锋的推进距离,L;Z为沿畦长各点的入渗水深,L. 引入地表下储水形状系数σZ=∫x0ZdsZ0xσΖ=∫0xΖdsΖ0x,则利用畦首入渗水深求单位畦宽上渗水总量的计算式为 VZ=σZZ0x=σZKtαx(3)VΖ=σΖΖ0x=σΖΚtαx(3) 式中:Z0为畦首入渗水深,L;t为畦首的入渗历时,等于放水时间,T;α、K为式(1)中的入渗参数;σZ为地表下储水形状系数;其它符号意义同前. Fok和Bishop在假定地表水流推进过程符合幂函数规律、土壤入渗规律符合Kostiakov模型的基础上,通过数学推导,得出地表下储水形状系数σZ的计算式为 σZ=α+r(1?α)+1(1+α)(1+r)(4)σΖ=α+r(1-α)+1(1+α)(1+r)(4) 式中:r为地表水流推进过程幂函数表达式x=ptr中的经验指数.在该幂函数表达式中,x为水流前锋的推进距离,L;α为式(1