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

工程建筑学相关资料ENGINEERING ARCHITECTURE RELATED INFORMATION 工程建筑学相关资料 ENGINEERING ARCHITECTURE RELATED INFORMATION PAGE PAGE 1 作物生长田间水分平衡的系统模拟  摘要：本文以作物生长模型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)，因此有：式中 h(cmH2O)为压力水头，K(h)为导水率(cm d-1-1)， C=dθ/dh为比水容量，Z为深度座标(cm)，t为时间(d)，Sw为根系吸水函数.　田间条件下，土壤水分运动有多种边界条件，而它们又反复交替变化.本模型有：　初始条件： h=hi， Z＞0， t=0.　上边界条件：(1)表层土壤接近饱和或有极薄积水时，h=0， Z=0， t＞0. (2)土表向土壤供水，供水强底低于入渗率，K(h)(h)/(Z)-K(h)=-i， Z=0， t＞0. (3)土表蒸发而土壤供水充足，K(h)(h)/(Z)-K(h)=E， Z=0， t＞0. (4)土壤蒸发而土体向土表供水速率低于蒸发速率，土表趋于风干状态，h=had， Z=0， t＞0; 其中脚标ad表示风干土.　　下边界条件：下边界条件主要考虑地下水位埋深ZW(m)的影响　　(1) ZW＞3.0; h=hi， Z=Zmax， t＞0. (2)3.0≥ZW≥1.0; h=ZW-Zmax， Z=Zmax， t＞0. (3)ZW＜1.0; h=0， Z=ZW， t＞0. 其中，当ZW＜1.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个线性方程求解.　　对土壤中的非饱和流，考虑到土壤水分运动模型与作物生长动态模型的综合模拟计算量较大，故采用基质流势(Matric flow potential)［8，9］的概念求解水流方程.定义U为基质流势U≡-∫hK(h)dh，则水流通量方程可改写为：　　　　　　　　　　　　　　q=-(dU)/(dZ)+K(h) ，(3)假设两相邻土层中心的流量随深度的分布不变，微薄土层中心之间的水势变化为线性的，则式(3)可以导出水流通量方程的积分形式为：