基于Morgan模型灌区灌溉水系统优化.docVIP

基于Morgan模型灌区灌溉水系统优化 　　摘要:从整个灌溉水系统着手,研究灌溉水系统管理技术,提出了Morgan模型,利用线性规划方法求出某种状态下的最优目标函数值及每次灌水决策问题;田间灌溉系统优化模型求解是由上面作出的灌水决策,以效益函数为目标函数,根据约束条件,按解线性规划原理,求解出灌水时间、入畦流量、畦宽、畦长等问题;渠系优化调度求解模型是采用决策数学中解整数规划的背包问题方法,针对灌区配水渠道“定流量,变时间”的运行方式,求解灌区渠系的优化配水问题。 　　关键词:灌区;Morgan模型;优化;模型求解 　　 　　1模型的建立 　　 　　针对不同的需要,建立了各种形式的模型,它们分为两类:静态和动态模型。静态模型是以作物耗水量为变量,寻找作物不同生育阶段不同程度水分亏缺与产量的关系;动态模型是跟踪作物生长过程,建立不同农田水分状况与干物质累积量之间的关系。可见,动态模型具有连续动态的特征,这里采用Morgan等人(1980)提出的动态产量模型。 　　1.1Morgan 模型的形式 　　在水分胁迫下,最终干物质产量可表达为: 　　 YD=YDO Γ(t)#8226;p(Amt) (1) 　　式中:YD――最终干物质产量; 　　YDO――计算起始时的干物质产量; 　　n――计算时段个数; 　　Γ(t)――干物质产量增加率; 　　p(Amt)――修正系数,随着土壤相对有效含水率而变。 　　1.2干物质产量与籽粒产量 　　一般而言,干物质累积总量愈大,籽粒产量愈高。经多年资料分析,夏玉米、冬小麦干物质产量与籽粒产量关系式如下: 　　冬小麦:Yg=8.53876+0.32258YD 　　夏玉米:Yg=-66.70+0.5862YD 　　1.3输配水损失 　　设变量Xij(0,1),其中,i表示斗渠序数,j表示出水口序数,当出水口j由斗渠i取水时,有Xij=1,反之Xij=0。配水渠配水损失总量为: 　　Z= Cifi( Xij),(M≤N) (2) 　　式中:Ci――序号为i的斗渠沿程输水损失。 　　1.4模型描述 　　灌溉水系统优化管理决策模型的设计过程,是一个多阶段决策过程,用动态规划的数学模型描述如下:①阶段变量。将作物的生长过程每日划为一个时段(n),n 为生育期天数序号,n=1,2,3……N,N为生育期天数。②状态变量。a各时段可利用的灌溉水量S(n)(m3/hm2)。b各时段根层内平均土壤含水量。③决策变量。a灌溉定额d(n),灌水时间由土壤含水率控制。b各时段作物腾发量ETn,它是土壤含水率、作物生长和气象因素的函数。c畦田宽度及畦田行水时间,与畦田土壤质地及平整状况有关。d灌区开启的支(斗)渠数目Xi(i),与轮灌制度有关。 　　④目标函数。 　　第一阶段:MaxYDΓ(t)#8226;P(mt)(3) 　　第二阶段:MaxPRO=YD#8226;PC -(WC+LC)-OC (4) 　　第三阶段: MinZ= Cifi( Xij) (5) 　　⑤边界条件。S(n)0;W(0)=W0 　　式中:W0――播前土壤含水率。 　　 　　2模型的求解 　　 　　2.1作物优化用水求解 　　①对状态变量进行离散。将占干土重百分数表示的调萎含水率(%)和田间持水率(%)用公式: 　　W0(i)=10γHW(i) (6) 　　转化为以水深表示土壤储水量WO(w) 和WO(C),然后对储水量和可供水量分别在区间[WO(w) ,WO(C)] 和[0,q(i)] 间以规定的间隔进行离散。 　　②利用直接搜索法寻求递推方程的解。对于n时段,求解结果得到不同可灌溉水量和含水率状态下的二维优化灌水量表,即有一对[WO(j),q(l)],就对应一个优化决策d(i,j,l)和相应的目标函数YD(i,j,l)。 　　③第i=n阶段的求解完成后,令i=n-1,以递推方程为目标函数进行求解,求解后得i=n-1时的二维优化灌水量表。 　　④令i=n-2,n-3…1,重复第(3)步骤,直至全生育期各时段的优化过程完成。 　　⑤顺序决策过程。利用播前含水率W0和可供灌溉水量S(n),依次查i=1,2,3…n 对应的结果,则可得到顺序决策结果。 　　根据土壤初始含水率及气象资料(或气象预报)。求出阶段初ETO及阶段末土壤含水率。判断诸阶段是否存在临界含水率过渡点,并据此进行灌水决策。用状态转移方程计算阶段末土壤含水率、下阶段初可供灌溉水量S(n)及相应的目标函数。根据可供灌溉水量,利用线性规划方法求出某种状态下的最优目标函数值及每次灌水决策。 　　2.2田间灌溉系统优化求解 　　根据上面作出的灌水决策,即灌溉定额及累积生物产量,把生物产量转化为籽粒产量,以效益函数