作物生长田间水分平衡的系统模拟

1、作物生长田间水分平衡的系统模拟 06-02-21 10:33:00 作者：吕 军编辑：studa9ngns摘要：本文以作物生长模型MACROS(L1D)为基础，建立了土壤水分运动和作物生长耦合关系综合模拟系统. 着重介绍了SPAC体系水分运动子模型，并对土壤水分运动模拟的结果和各类参数的灵敏度等进行了分析.最后，应用该模拟系统对南方多雨地区冬小麦田土壤水分平衡和土壤渍害的发生过程进行了模拟分析. 关键词 关键词：SPAC水分运动 水分与作物生长耦合模型 水分平衡 渍害 地表排水 近年来数值模拟方法在农田土壤水分运动研究中被广泛采用1,4，并逐步从实验室引向田间5,6，也考虑了农田作物状况对土壤

2、水分运动的影响2,3,7,9.但是，国内大多数水分运动模型没有涉及土壤水分条件对作物生长过程的动态影响，因而不能很好地反映土壤水分运动与作物生长动态之间所具有的强烈的反馈和耦合作用，实际上，不良水分条件限制作物生长，而作物生长状况的改变又对农田水分平衡过程产生重要影响.同时，作物的不同种植方式也将影响农田水分平衡状况.在作物生长动态模型MACROS(L1D)8和田间水分平衡模型的基础上，建立并引进了农田水分平衡与作物生长动态耦合子模型，本文着重讨论水分运动子模型的模拟方法，分析模型模拟结果的可靠性及其对主要参数的敏感性，并应用该模型研究了多雨地区冬小麦生长田间土壤水分动态，讨论多雨地区土壤渍、

3、涝的发生和消退过程.1土壤水分运动数学模型以近作物冠层的大气为上边界，土体一米深处为下边界，构成作物生长田间农田水分平衡系统.作物生长模型MACROS(L1D)以1d(24h)为时间步长，在作物出苗后开始模拟，并向土壤水分平衡模型提供所需作物参数.而土壤水分运动模型则以可变时间步长(0.1-0.001d)进行模拟.因此，土壤水分状态经24h循环计算后才取值与作物发生相互作用.当不良土壤水分条件发生时，作物生长受到限制，从而也改变了向土壤水分平衡模型提供的有关作物参数，这样就构成了土壤水分平衡与作物生长动耦合关系.模型中，仅考虑作物特性、土壤水分特性和气象要素，并认为其它条件均满足作物生长要求.

4、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)为导水率(cm d-1-1), C=d/dh为比水容量，Z为深度座标(cm)，t为时间(d)，Sw为根系吸水函数.田间条件下，土壤水分运动有多种边界条件，而它们又反复交替变化.本模型有：初始条件： h=hi, Z0, t=0.上边界条件：(1)表层土壤接近饱和或有极薄积水时，h=0, Z=0, t0. (2)土表向土壤供水，供

5、水强底低于入渗率，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时，地下水位埋深即为模拟系统下边界，流入地下水位中的水

6、亦即为土壤水渗漏损失.土壤中发生饱和流时， /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个线性方程求解.对土壤中的非饱和流，考虑到土壤水分运动模型与作物生长动态模型的综合模拟计算量较大，故采用基质流势(Matric flow potential)8,9的概念求解水流方程.定义U为基质流势U-hK(h)dh,则水流通量方程可改写为：q=-(dU)/(dZ)+K(h) ,(3)假设两相邻土层中心的流量随深度的分布不变，微薄土层中心之间的水势变化为线性的，则式(3)可

7、以导出水流通量方程的积分形式为：(4)当水流发生在非均质土层之间时，式(4)中的积分项由下式修正.(5)作物根系吸水函数Sw, 先由作物冠层蒸腾量E确定其总量，当土壤供水充足时，各土层中根系吸水量Swi，根据已确定的根系总量和由实验测定的根系在土壤剖面上的分布比例i计算，即Swi=iE .(6)不良土壤水分条件发生时，根据旱害、渍害程度对式(6)进行修正，修正参数由实验确定9.1.2水分汽态流Ep=ESR+E (7)其中(8)式(8)中，SR为辐射吸收量，为水的汽化潜热，S为饱和水汽压随温度变化曲线的斜率，P和*分别为空气干燥力和表观温度计常数项，它们又分别为：(9)式中es和ea分别为饱和水

8、汽压和水汽压实测值，Cp为空气的体积热容量，为湿度计常数，而h、b和s分别为边界层热交换阻力、边界层水汽交换阻力和气孔阻力.式(8)中，ESR和E受作物冠层叶面积状况的影响，因此作物潜在蒸腾速率Ecp,应由叶面积系数ALV进行校正8，则式(8)可改写为：(10)式中ES为土壤蒸发率，bL和bS分别为作物遮阴修正系数和土壤湿度修正系数，即(12)式中m和n分别经验系数，分别与作物种类和土壤性质有关，ALV为包括死叶和活叶的总叶面积系数，a和f则分别为风干土壤湿度和田间持水量.考虑到实际土壤蒸发并不都发生在一个峰面上，则应由各土层蒸发贡献率mi将ES值分配到各土层中去：(13)式中脚标i为土壤层次，d为土层厚度，L为土层深度，K为经验系数.此外，为避免地下侧渗流的引入而造成模型过于复杂，地下水埋深按实际观测值输入.1.3模型软件及其输入和输出整个模型模拟软件用CSMP语言编制(连续系统模拟程序语言).模型所需输入资料包括作物资料(品种特性和生长特性资料)、气象资料(日最高和最低