华北地区夏玉米田间水分转化规律研究

1、华北地区夏玉米田间水分转化规律研究作者：刘群昌 谢森传时间：2007-11-25 12:15:00 摘要：本文利用一维垂直非饱和土壤水运动计算模型，通过模拟计算夏玉米生育期各生长阶段不同水文年、不同灌水定额下田间水分转化规律，对棵间蒸发，深层渗漏与灌水定额之间的数量关系进行了定量研究，并建议华北地区夏玉米节水灌溉定额为40m/亩. 关键词：夏玉米 田间水分转化 灌水定额 1 田间水分转化的动力学模式田间水分的转化过程主要包括降雨或灌水的入渗，径流，根区土壤水分运移，根系吸水，蒸发蒸腾，根区以下土壤水运动等.描述土壤水运动的动力学方程是根据达西定律和连续方程结合推导出来的，由于方程是非线性的，所

2、以只有在特定的初始和边界条件下才能用解析法求解.自60年代以来，随着电子计算机技术的发展，人们能够借助计算机运用数值模拟的方法，对一般条件下的土壤水分运动问题求解，并得到了满意的结果.本文采用HYDRUS模型软件进行数值模拟.1.1 HYDRUS模型 HYDRUS模型软件8是美国盐碱实验室在Worm模型的基础上的改进版，用于模拟计算一维垂直非饱和流和溶质运移，考虑了作物根系吸水和土壤持水能力的滞后影响，适应于恒定或非恒定的边界条件，具有灵活的输入输出功能，可用来模拟非匀质土壤，最多可模拟5种20层土壤.模型中方程解法采用Gakerin线性有限元法.1.2 方程和边界条件 在忽略土壤侧向水流运动

3、，仅考虑一维垂向运移时，有根系吸水项的土壤水分运动方程为1：c(h)(h)/(t)=()/(Z)k(h)(h)/(Z)-k(h)-s(z,t)(1)式中c(h)比水容重L-1,c(h)=d/dh,h土壤压力水头L，体积含水量L3L-3,k(h)水力传导度LT-1，s(z,t)单位体积根系吸水率T-1,z土壤深度L、向下为正，t时间T.初始条件：h(z,t)=h0(z), t=0.上边界条件：h(0,t)=h0(t), (积水条件)-k(h)(h)/(z)+k(h)/z=0=q0(t)(2)式中h0(t)已知压力水头L，q0(t)净通量，正通量表示下渗，负通量表示蒸发(棵间蒸发).下边界条件：h

4、(l,t)=hi(t),-k(h)(h)/(z)+k(h)/z=l=q1(t)或(h)/(z)/zl=0, (3)2 作物根系吸水模式本文采用VanGenuchten模型8，即s(z,t)=Ep(t),(z)(h,h0)(4)式中S(z,t)单位根据吸水率，Ep(t)作物最大蒸腾率，(h,h0)盐分应力函数，反映土壤盐分对田间根系吸水的影响.(h,h0)=(1)/(1+(h+h0)/(h50)p(5)式中h压力水头，h0渗透压，与溶液浓度C有关，h0=a1c,a1简单换算系数，其值依赖于压力水头和浓度表达的单位.h50作物潜在蒸腾率减少50%的土水势，对玉米而言，h50为-0.25-0.65M

5、Pa，一般取-0.43MPa(Ehler,1983),P经验常数，P3；(z)根系密度分布函数相对值. (6)其中 Lr根层深度，(z)根系密度分布函数，(z)是根深Lr的函数，实际应用中将Lr分为若干层且认为每层内根系分布是均匀的.根系密度采用层内干根重占根区总干根重的比值表示，由田间实测获得.3 非饱和导水率的确定非饱和导水率K是土壤水分运动的重要参数，在对非饱和土壤水运动基本方程进行数学分析时，无论用解析解或数值解的方法，都要用到它.获得K的方法有两种：一是实验方法，如瞬时剖面法，垂直下渗通量法，垂直土壤稳定蒸发法，结壳法等；另一种是间接法，即根据已知的水分特征曲线和其他条件推导出的函数

6、形式，如Van Genunchten(1980)将土壤水分特征曲线的函数形式与Mualem(1976)导出的用来预测非饱和导水率的函数形式相结合得到如下函数关系5：(h)=r+(s-r)/(1+(h)nm),(7)k(l)=kss1/211-(1-s1/ml)m2,(8)sl=(-r)/(s-r),(9)式中s饱和含水量；r残余含水量；Ks饱和导水率；Sl相对饱和度，,n,m拟合参数，通过非线性最小二乘法对室内试验获得的土壤持水数据(,h)进行拟合求模糊.将式(7)代入式(8)得：k(h)=ks1+(h)n-2/m1-1-(1+h)n-m2(10)表1 VG非饱和导水率函数有关参数土层深度/c

7、m饱和含水量/cm/cm残余含水量/cm/cm/a/cmnR20800.490.120.0181.330.97801000.510.0960.005531.1890.981002200.480.080.00871.4790.982204000.430.0640.00231.3870.994 作物腾发量的计算作物腾发量又称作物耗水量，是指在作物生长季节，从生长面积上失去的水量，它包括从作物体蒸腾的水量和组成作物体内的水量(所占比例很小，可忽略不计)以及从种植面积上棵间蒸发的水量.目前最常用的最大作物腾发率的计算方法是先计算参考作物腾发率ET0(t),然后将其乘以一个作物系数Kc(t),即：ETc

8、(t)=Kc(t)ET0(t);作物系数Kc由河北望都灌溉试验站的灌溉资料确定.参考作物腾发率由FAO新近推荐的Penman-monteith方法计算6，FAO-Penman-monteith方法把参考作物腾发量重新定义为“作物高度0.12m，固定叶而阻力70sm-1，反射率0.23的假想参考作物的腾发量”，由此结合monteith方法可得出FAO-Penman-monteith方程下：ET0=(0.408(Rn-G)+(900)/(T+273)U2(ea-ed)/(+(1+0.34U2)(24h) (11)式中 ET0假想草的参考腾发量(mmd-1),Rn净辐射(mJm-2-2d-1),G土

9、壤热通量，干湿球常数，U22m高处风速，ea饱和水汽压(kPa)，ed实际水汽压(kPa)，压力曲线斜率，P大气压力(KPa).5 棵间蒸发与植株蒸腾的划分作物腾发量中叶面蒸腾与棵间蒸发的分摊是农田水分循环以及土壤-植物-大气连续体水分传输动态模拟研究中必不可少的工作之一，同时也是一件困难的事情.Richie 和 Burnett(1971)研究了棉花和谷类作物的分摊系数和叶面积指数LAI的关系，提出了充分供水条件下叶面蒸腾Tp与ETc关系式：Tp=(-0.21+0.7LAI1/2)ETc,(0.1LAI2.7)，Tp0. (LAI0.1)，(12)国内对这方面研究工作开展不多，西北农业大学康绍

10、忠通过研究认为4：作物蒸腾在总腾发中所占的比例依赖于提供到作物冠层和棵间土壤表面的净辐射以及各部分的传输阻力，若叶面蒸发与腾发用波纹比能量平衡法计算：ETc=Rn/(1+1),Tp=Rnc/(1+2)，在充分质水条件下12，由此可导出：Tp/ETc=Rnc/Rn.据田间连续几年观测，得到关系：(13)式中t为一日中的时间，从零点开始；K、A为经验系数，对于玉米K=0.4016，A=0.9872；Tp叶面蒸腾，ETc作物腾发量.本文利用实测的LAI数据分别用康绍忠公式和Chids公式计算了雄县试验站1995年夏玉米在不同生育阶段的叶面蒸腾Tp并进行了线性回归分析，其相关系数R2=0.98，说明两

11、种公式计算Tp有很好的一致性，应用于实际比较可靠，本文模拟计算中潜在蒸腾率采用康绍忠公式计算结果.6 实际棵间蒸发的估算在模拟计算中，上边界通量包括实际棵间蒸发量，但比较准确地确定实际棵间蒸发量目前仍很困难，为此本文采用一种试算方法，即根据实测叶面积指数按康绍忠的棵间蒸发，叶面蒸腾分配公式求出各个阶段中棵间蒸发与叶面蒸腾的比例值，然后在计算中先输入潜在棵间蒸发，则有一个叶面蒸腾值(等于根系吸水)的输出，得出一个模拟的棵间蒸发与根系吸水的比值Es/S，判断两个比值的大小，如果不一致，则改变输入的棵间蒸发值，再进行模拟计算，至直两个比值比较接近.由式(13)可推导出：(14) 式中Es棵间蒸发，s

12、叶面蒸腾.表2 各个生育阶段的Es/S值生育阶段播种-出苗出苗-拔节拔节-抽穗抽穗-灌浆成熟平均LAI0.42255.386.48Es/ETc0.9140.2100.062Es/S10.6280.2660.066图1 玉米田间土壤剖面负压实测值与模拟值比较8 田间水分转分分析8.1 不同水文年田间水分转化模拟 表3和4分别为特旱年，干旱年、平水年和丰水年以及不同灌溉处理的模拟结果，结果表明：(1)夏玉米生长期降水量虽然占全年的7085%，但通过模拟计算表明其不同年份的实际腾发量相差仍然较大，特旱年只占作物最大腾发量的55%，灌13次水即可提高到80%左右；干旱年占79%，灌1次水可提高到90%

13、，所以对夏玉米进行灌溉是提高其产量所必需的.(2)供水量越大，作物吸水作用越强，而棵间无效消耗(棵间蒸发)增加较缓慢.(3)夏玉米的生育阶段需水量是由棵间土壤蒸发与叶面蒸腾量组成，在不同生育阶段，它们之间的比例变化很大，从夏玉米播种到拔节，恰处于6月中下旬至7月上旬，气温高，大气干燥.此时植株矮小，叶面积指数小，叶面蒸腾量很低，棵间土壤蒸发量却占较大的比例，从表4看出：该阶段棵间蒸发量比例达90%以上，从全生育期看，玉米棵间土壤蒸发量占总需水量的40%左右.棵间土壤蒸发量对产量形成基本上无意义，应当采取适当的栽培措施，尽量降低它所占的比例.80年代以来，采用地膜覆盖进行节水的田间管理，基本上消

14、除了棵间土壤蒸发量，减少了农田耗水量，有利于提高水的利用效率，应尽量推广应用.(4)从表3看出棵间蒸发与深层渗漏之和占总供水一般在3545%之间，所以从供水角度来讲，降低二者所占比例，节水用水意义很大.表3 不同水文年田间水分转化模拟结果水文年处 理EsSETaETa/ETcDD/(P+I)WD+Es(D+Es)/(P+I)特旱年不灌82.11126.6208.75524.410.4106.545.0灌1水92.78178.9271.77229.29.48.8122.039.3灌2水94.13198.7292.97742.0210.951.4136.235.3灌3水99.91219.5319.48474.