水利工程论文-作物生长田间水分平衡的系统模拟.doc

水利工程论文-作物生长田间水分平衡的系统模拟摘要：本文以作物生长模型MACROS(L1D)为基础，建立了土壤水分运动和作物生长耦合关系综合模拟系统.着重介绍了SPAC体系水分运动子模型，并对土壤水分运动模拟的结果和各类参数的灵敏度等进行了分析.最后，应用该模拟系统对南方多雨地区冬小麦田土壤水分平衡和土壤渍害的发生过程进行了模拟分析.关键词关键词：SPAC水分运动水分与作物生长耦合模型水分平衡渍害地表排水近年来数值模拟方法在农田土壤水分运动研究中被广泛采用1,4，并逐步从实验室引向田间5,6，也考虑了农田作物状况对土壤水分运动的影响2,3,7,9.但是，国内大多数水分运动模型没有涉及土壤水分条件对作物生长过程的动态影响，因而不能很好地反映土壤水分运动与作物生长动态之间所具有的强烈的反馈和耦合作用，实际上，不良水分条件限制作物生长，而作物生长状况的改变又对农田水分平衡过程产生重要影响.同时，作物的不同种植方式也将影响农田水分平衡状况.在作物生长动态模型MACROS(L1D)8和田间水分平衡模型的基础上，建立并引进了农田水分平衡与作物生长动态耦合子模型，本文着重讨论水分运动子模型的模拟方法，分析模型模拟结果的可靠性及其对主要参数的敏感性，并应用该模型研究了多雨地区冬小麦生长田间土壤水分动态，讨论多雨地区土壤渍、涝的发生和消退过程.1土壤水分运动数学模型以近作物冠层的大气为上边界，土体一米深处为下边界，构成作物生长田间农田水分平衡系统.作物生长模型MACROS(L1D)以1d(24h)为时间步长，在作物出苗后开始模拟，并向土壤水分平衡模型提供所需作物参数.而土壤水分运动模型则以可变时间步长(0.1-0.001d)进行模拟.因此，土壤水分状态经24h循环计算后才取值与作物发生相互作用.当不良土壤水分条件发生时，作物生长受到限制，从而也改变了向土壤水分平衡模型提供的有关作物参数，这样就构成了土壤水分平衡与作物生长动耦合关系.模型中，仅考虑作物特性、土壤水分特性和气象要素，并认为其它条件均满足作物生长要求.1.1水分液态流土壤水分液态流为有根系吸水条件下的垂直一维流，即/t=-q/Z-Sr,其中q为水流通量，q=-K(h)h/Z+K(h),因此有：C(h)(h)/(t)=()/(Z)K(h)(h)/(Z)-K(h)-Sw,(1)式中h(cmH2O)为压力水头，K(h)为导水率(cmd-1-1),C=d/dh为比水容量，Z为深度座标(cm)，t为时间(d)，Sw为根系吸水函数.田间条件下，土壤水分运动有多种边界条件，而它们又反复交替变化.本模型有：初始条件：h=hi,Z0,t=0.上边界条件：(1)表层土壤接近饱和或有极薄积水时，h=0,Z=0,t0.(2)土表向土壤供水，供水强底低于入渗率，K(h)(h)/(Z)-K(h)=-i,Z=0,t0.(3)土表蒸发而土壤供水充足，K(h)(h)/(Z)-K(h)=E,Z=0,t0.(4)土壤蒸发而土体向土表供水速率低于蒸发速率，土表趋于风干状态，h=had,Z=0,t0；其中脚标ad表示风干土.下边界条件：下边界条件主要考虑地下水位埋深ZW(m)的影响(1)ZW3.0；h=hi,Z=Zmax,t0.(2)3.0ZW1.0；h=ZW-Zmax,Z=Zmax,t0.(3)ZW1.0；h=0,Z=ZW,t0.其中，当ZW1.0m时，地下水位埋深即为模拟系统下边界，流入地下水位中的水亦即为土壤水渗漏损失.土壤中发生饱和流时，/t=0,有：q=-Ki(hi+1-hi)/(Zi+1-Zi)i,(h2-h1)+(h3-h2)+(hn+1-hn)=hn+1-h1,(2)i=n,为土壤层次数，得到n+1个线性方程求解.对土壤中的非饱和流，考虑到土壤水分运动模型与作物生长动态模型的综合模拟计算量较大，故采用基质流势(Matricflowpotential)8,9的概念求解水流方程.定义U为基质流势U-hK(h)dh,则水流通量方程可改写为：q=-(dU)/(dZ)+K(h),(3)假设两相邻土层中心的流量随深度的分布不变，微薄土层中心之间的水势变化为线性的，则式(3)可以导出水流通量方程的积分形式为：(4)当水流发生在非均质土层之间时，式(4)中的积分项由下式修正.(5)作物根系吸水函数Sw,先由作物冠层蒸腾量E确定其总量，当土壤供水充足时，各土层中根系吸水量Swi，根据已确定的根系总量和由实验测定的根系在土壤剖面上的分布比例i计算，即Swi=iE.(6)不良土壤水分条件发生时，根据旱害、渍害程度对式(6)进行修正，修正参数由实验确定9.1.2水分汽态流土壤蒸发和作物蒸腾由彭曼-蒙特斯公式计算1.2.1作物蒸腾设Ep为作物潜在蒸腾速率，ESR和E分别为太阳辐射力蒸发项和空气干燥力蒸发项，有：Ep=ESR+E(7)其中(8)式(8)中，SR为辐射吸收量，为水的汽化潜热，S为饱和水汽压随温度变化曲线的斜率，P和*分别为空气干燥力和表观温度计常数项，它们又分别为：(9)式中es和ea分别为饱和水汽压和水汽压实测值，Cp为空气的体积热容量