《土壤-作物系统模拟的优化方法简介06》由会员分享，可在线阅读，更多相关《土壤-作物系统模拟的优化方法简介06（62页珍藏版）》请在金锄头文库上搜索。
1、土壤作物系统模拟的 优化方法简介 陈研 Tel: 62732959新学科综合楼 4-201 中国农业大学资源和环境学院中国农业大学资源和环境学院20062006年年1 1月月现代农业生产的目标是多方面的，不仅要求具有很好的经济效益，同时也需要一定的社会及生态效益。在生产过程中仅追求高的产量常导致资源的低效利用，并伴随环境污染问题；由此，按农业可持续发展的要求，应在达到一定目标产量的前提下，要求资源利用效率尽可能提高，从而实现农业生产的高产高效。反之，若只重视资源的高效利用又可能存在低产的现象。拟建整合模型的主要内容和结构拟建整合模型的主要内容和结构作物生长发育作物生长发育 （生理响应等）（生理
2、响应等）根冠形态发育根冠形态发育 3 3维模型维模型根系吸收水分根系吸收水分 模型模型根系吸收养分根系吸收养分 模型模型土壤水分、热、养分运移和转化模型土壤水分、热、养分运移和转化模型调控途径调控途径逆境条件逆境条件土壤作物系统模拟的优化问题一、管理措施的优化选择合适的灌水、施肥量，使得作物产量、水分利用效率(WUE)和氮素利用效率(NUE)都达到最优。常用的优化方法:1.线性规划例1：农业用水最优分配的线性规划模型农业用水对象为小麦、玉米、经济作物、果树等四项，其用水对策拟定为畦灌、喷灌、滴灌、雨养（即不灌）四项。经过两层次二元组合，得出 14个决策变x1x14，由于滴灌应用于大田作物是不适
3、宜的，所以不考虑小麦、玉米的滴灌方案。 目的：每万亩土地净收益最大决策变量 小 麦玉 米经济作物果 树畦 灌x1x4x7x11喷 灌x2x5x8x12滴 灌 x9x13雨 养x3x6x10x14用水对象 用水对策限制因素：1.（总面积限制） 2. x1 + x2 + x3 77 （小麦面积限制） 3.（小麦畦灌面积限制） 4.（小麦喷灌面积限制） 5. x4+x5 + x6 20 （玉米面积限制） 6.（玉米畦灌面积限制） 7.（玉米喷灌面积限制） 8. x7 + x8 + x9 + x10 26 （经济作物面积限制） 9. （经济作物畦灌面积限制） 10.（经济作物喷灌面积限制） 11. x
4、9 4.1 （经济作物滴灌面积限制） 12. x11 + x12 + x13 + x14 4.1 （果树面积限制） 13. x11 2 （果树畦灌面积限制）14. （果树喷灌面积限制） 15. x13 0.9 (果树滴灌面积限制）16. 210x1 + 110x2 + 190x4 + 100x5 + 200x7 + 100x8 + 70x9 + 180x11 + 95x12 + 60x13 11 707 （水资源量限制） xj 0 (j = 1, 2, , 6) （非负限制） (式中变量的系数为水资源消耗系数,如x1的系数210表示畦灌 每万亩小麦耗水210万立方米。)目标函数：Max f =
5、 60x1 + 55x2 + 35x3 + 52x4 + 47x5 + 30x6 + 100x7 + 90x8 + 85x9 +60x10 + 600x11 + 550x12 + 520x13 + 350x14 目标函数中各系数均为每万亩土地净收益（万元） 线性规划模型的一般形式求一组决策变量x1，x2，xn，满足下列约束条件：使目标函数取最大值（或最小值）。 线性规划问题的解法： 单纯形法 ！ 例 1：某工厂生产甲、乙两种产品，已知有关数据见下表：III拥有量原材料2111设 备1210利润 元/件810试求获利最大的生产方案。2. 目标规划 ：解：这是一个单目标问题，用线性规划模型描述为：
6、最优解为：元但实际上工厂在作决策时，要考虑市场等一系列其他条件：（1） 根据市场信息，产品I的销售量有下降的趋势，故考虑产品甲的产量不大于产品II。（2） 超过计划供应的原材料时，需要高价采购，这就使成本增加。（3）应尽可能利用设备I ，但不希望加班。（4） 应尽可能达到并超过计划利润指标56 元。1. 偏差变量设 为决策变量，此外，引进正、负偏差变量正偏差变量 表示决策值超过目标值的部分；负偏差变量表示决策值未达到目标值的部分。 2、绝对约束和目标约束绝对约束是指必须严格满足的等式约束和不等式约束,是硬约束； 目标约束是目标规划特有的，可把约束右端项看作要追求的目标值。在达到此目标值时允许发
7、生正或负偏差，因此在这些约束中加入正、负偏差变量，是软约束。 线性规划问题的目标函数，在给定目标值和加入正、负偏差变量后可变换为目标约束。也可根据问题的需要将绝对约束变换为目标约束。 3、优先因子（优先等级）与权系数一个规划问题常常有若干目标。但决策者在要求达到这些目标时，是有主次或轻重缓急的不同。凡要求第一位达到的目标赋于优先因子P1，次位的目标赋于优先因子P2，并规定PkPk+1,k = 1, 2, , K。表示Pk比Pk+1有更大的优先权。即首先保证P1级目标的实现，这时可不考虑次级目标；而P2级目标是在实现P1级目标的基础上考虑的；以此类推，若要区别具有相同优先因子的两个目标的差别，这
8、时可分别赋于它们不同的权系数j，这些都由决策者按具体情况而定。 4、目标规划的目标函数目标规划的目标函数(准则函数)是按各目标约束的正、负偏差变量和赋于相应的优先因子而构造的。当每一目标值确定后，决策者的要求是尽可能缩小偏离目标值。因此目标规划的目标函数只能是其基本形式有三种： （2）要求不超过目标值，既允许达不到目标值，就是正偏差变量尽可能地小，这时（3）要求超过目标值，既超过量不限，但必须是负偏差变量尽可能地小，这时（1）要求恰好达到目标值，既正、负偏差变量都要尽可能地小，这时目标规划的数学模型目标规划模型的一般形式为：使目标函数： 求x = (x1, , xn)T，满足条件：3. 非线性
9、目标规划非线性目标规划模型：在目标规划的基础上，其目标函数以及约束条件之一是非线性的。例：冬小麦夏玉米水氮管理措施的优化王凤仙，陈研，李韵珠； 2000这是一个多效应指标的综合优化分析，水氮效应函数多为非线性函数，从而选择非线性目标函数。优化目标（优先级）：1. 产量 （Y）； 2. 土壤水资源利用效率 （WUE）； 3. 土壤氮资源利用效率 （NUE）。约束条件：1. 冬小麦：灌水总量为40440mm; 施氮总量为0320kg/ hm22.夏玉米:灌水总量为40mm;施氮总量为0320kg/ hm23.冬小麦-夏玉米周年:灌水总量为80480mm; 施氮总量为0640kg/ hm2决策变量：
10、灌水总量 I，施氮总量N；优化设计方案：1. 冬小麦:干旱、中等、湿润3 种降雨年型;灌水处理为2. 5 个,即40 、140 、240 、340 和440mm;施氮量也是5个处理,0 、80 、160 、240 、320kg/ hm2 。2. 夏玉米: 以前茬冬小麦为干旱年型、灌水量240mm、施氮量160kg/ hm2 为基础条件,夏玉米降雨年型分为干旱、中等及湿润3 种;播前灌40mm ,无灌水处理;氮处理5 个,数量与小麦同。3.冬小麦-夏玉米周年:冬小麦、夏玉米均为中等降雨年型,两种作物的水氮处理同上,即周年灌水总量处理在80 至480mm 之间, 施氮总量在0 至640kg/ hm
11、2 之间。目标函数见原文。4. 动态规划动态规划是运筹学的一个分支，是求解多阶段决策过程（multistep decision process）的最优化方法。多阶段决策：有一类活动的过程，由于其特殊性，可以分为若干个互相联系的阶段，在它的每一个阶段都需要作出决策，从而使整个过程达到最好的活动效果。各个阶段的决策的选取不是任意确定的，它依赖于当前面临的状态，又影响以后的发展。当各个阶段决策确定后，就组成了一个决策序列，因而也就决定了整个过程的一条活动路线。这种把一个问题可看作是一个前后关联具有链状结构的多阶段过程就称为多阶段决策过程，这种问题就称为多阶段决策问题。动态规划的基本概念：(1) 阶段
12、把所给定问题的过程，恰当地分为若干个相互联系的阶段，以便能按一定的次序去求解，描述阶段的变量成为阶段变量。(2) 状态在动态规划模型中，状态表示每个阶段开始所处的自然状况或客观条件，它包含了描述系统情况所必需的信息，所以状态又叫做不可控因素。状态既是该阶段的起始状况，又是前一阶段末时的状况，通常一个阶段有若干个状态，描述过程状态的变量称为状态变量。 状态应具有下面的性质：如果某阶段状态给定后，则在这阶段以后过程的发展不受这阶段以前各段状态的影响。换句话说，过程的过去历史只能通过当前的状态去影响它未来的发展。这个性质称为无后效性。（马尔科夫性）(3) 决策当过程处于某一阶段的某个状态时，可以做出
13、不同的决定（或选择），从而确定下一阶段的状态，这种决定（或选择）称为决策，决策用决策变量来描述。 (4) 策略按顺序排列的决策组成的集合(5) 状态转移方程状态转移方程，又叫系统方程，是确定过程由一个状态到另一个状态的演变过程。 (6) 指标函数和最优值函数任何决策过程都必然有一个衡量其策略优劣的尺度，一个数量指标，称为指标函数（目标函数）。指标函数的最优值，称为最优值函数k=1k=2kk=n x1x2x3xk+1xn+1xku1u2ukunf1(x2)f2(x3)fk(xk+1)fn(xn+1)动态规划中各阶段、各变量之间关系xi : 状态ui : 决策变量是指第 k 阶段状态为 xk+1时
14、，相对于始点的最优指标函数值。而表示第阶段状态为取决策为时对本阶段的阶段效益值。在一些动态规划模型中，每个阶段内都有多于一个的目标函数，动态规划和目标规划的思想常常结合起来加以应用！ 例：基于作物生长模型的冬小麦田间优化管理冯凌，陈研。2005作物生长模型PS123模型PS123模型是荷兰Wageningen大学土壤地质系Driessen P.M.教授等人开发的用于定量化土地生产力评价的普适模型，可以用来进行小麦、玉米、粟类及棉花等多种作物生产的生长模型。PS123模型描述了生产水平的层次性。 PS1（最简单的生产水平）：是在作物生理可能达到的范围内把作物行为定量为只是有效的太阳辐射和温度的影
15、响；其它的土地质量因子则假设为充分满足相应的土地利用需求。 PS2 ：改变水分供给最适宜的假设，即认为作物生产是由所截获的太阳辐射量、温度和有效水分共同决定的，其它在实际农田生产中影响生产的土地产量因子都假设为最适宜。 PS3：有效养分成为附加考虑的因子 在原PS2模型中主要采用的是水量平衡模型计算实际的作物蒸腾（TR）和作物充分供水下的实际蒸腾速率TRM ，从而求得作物水分满足系数 。平衡模型在应用于区域均衡研究时，避开了一些土壤水分物理参数，计算相对简单；然而某些分量难以准确确定，因而精度不够高。1. 通过基于水动力学的土壤模块计算提供TR和TMR，在动态定量化地考虑土壤水转化和作物吸水等因素时更为直观、精确。 对原PS123 模型的改进：2.参照了Daisy模型对PS3模型进行了修正，并计算氮胁迫系数 动态规划模型：阶段变量：状态变量：将冬小麦的生育期划分为六个阶段（stage=6）1. t 阶段初始时可用于分配的灌
