土壤-作物系统模拟优化方法简介06.pptVIP

陈研 chenyan@cau.edu.cn Tel: 新学科综合楼 4-201 中国农业大学资源和环境学院 2006年1月 动态规划模型： 阶段变量： 状态变量： 将冬小麦的生育期划分为六个阶段（stage=6） 1. t 阶段初始时可用于分配的灌水量 （mm/ha）； 2. t 阶段初始时可用于分配的氮肥施用量 （kg/ha）； 3. t 阶段初始时冬小麦的生物量（包括叶、根、茎、穗分别 的生物量） （kg/ha）。 播种-出苗、出苗-返青，返青-拔节， 拔节-抽穗，抽穗-灌浆，灌浆-成熟。 决策变量： （mm/ha）：表示 t 阶段初的灌水决策。 （kg/ha）:表示 t 阶段初的施氮决策。 状态转移方程（系统方程） 1. 水量分配方程: 2. 肥量分配方程: 3.本模型属于土壤-作物系统中的一个子模块，因此，第三个 状态转移函数的角色由土壤模块和作物生长模块共同担任。 目标函数： 作物各器官干物质的增加量同水分胁迫系数和氮胁迫系数 的乘积成正比。而对于特定作物品种来说，粒穗比是一定的， 即贮藏器官的干物重与经济产量成正比。 具体为： 1. 初始阶段初始时可用于分配的灌水量 （mm/ha）； 2. 初始阶段初始时可用于分配的氮肥施用量 （kg/ha）； 3. 初始阶段初始时冬小麦的生物量（包括叶、根、茎、穗分别 的生物量） （kg/ha）。 初始条件： 约束条件： 每次进行田间管理时，一次性灌水量最大值为160 mm/ha， 即灌水允许决策集合为0~160 mm/ha；一次性施氮量最大值 为220 kg/ha，施肥允许决策集合为0~220 kg/ha。 存在的问题：状态的设计是否满足无后效性？ 进一步的研究：降雨对决策的影响 优化部分流程图： 5. 随机规划 带有随机参数的数学规划模型一般来说可以描述如下： 式中, X 为n 维决策向量; 为随机向量; 为目标函数; 是事先给定的约束条件的置信水平; P r{·}为表示{·}中事件成立的概率。 常见的约束方程的表达形式： i=1,2,…m 线性约束条件 基于作物生长模型的田间随机优化管理模型： 多目标非线性随机规划 二、反求参数 对土壤－作物系统的定量描述应包括以下三方面的内容 （1）过程的概化－模型的建立 （2）模型的求解－计算方法的选用 （3）各类参数的准确测定与计算 反求参数法 反求参数法就是从已知实测结果和计算结果来反求参 数，也叫做求逆问题或反问题。一般情况下，通过室内实 验求得的参数，和田间实验求得的参数相差很多，这样就 给实际计算带来很大的误差；无法通用。近年来引入了求 解微分方程逆问题的方法，并结合实验来反求参数。 逆问题的类型： 1.单纯求方程的参数； 2.求源汇项 3.求初始条件 4.求边界条件 5.以上几种类型的混合 使 初始参数有现有的实验或资料决定，曲线拟合时采用优化方法 用优化方法反求参数的基本思想是：求一组参数 单纯形法 单纯形法是对n维空间中的（n + 1）个点（它们构成 一个单纯形的顶点，如平面上一个三角形的三顶点）上的 函数值进行比较，丢掉其中最坏点代之以新点，从而构成 一个新的单纯形，以此逐步逼近极小点。 初始单纯形的选取： 取n + 1个顶点 ………………………… 的正规单纯形 其中 a 为单纯形边长 令 单纯形法的迭代过程，包括下述四种运算： * * 土壤－作物系统模拟的优化方法简介 现代农业生产的目标是多方面的，不仅要求具有很好 的经济效益，同时也需要一定的社会及生态效益。 在生产过程中仅追求高的产量常导致资源的低效利用， 并伴随环境污染问题； 由此，按农业可持续发展的要求，应在达到一定目标 产量的前提下，要求资源利用效率尽可能提高，从而实现农业生产的高产高效。 反之，若只重视资源的高效利用又可能存在低产的现象。 拟建整合模型的主要内容和结构 作物生长发育 （生理响应等） 根冠形态发育 3维模型 根系吸收水分 模型 根系吸收养分 模型 土壤水分、热、养分运移和转化模型 调控途径 逆境条件 土壤－作物系统模拟的优化问题 一、管理措施的优化 选择合适的灌水、施肥量，使得作物产量、水分利用 效率(WUE)和氮素利用效率(NUE)都达到最优。 常用的优化方法: 1.线性规划 例1：农业用水最优分配的线性规划模型 农业用水对象为小麦、玉米、经济作物、果树等四项， 其用水对策拟定为畦灌、喷灌、滴灌、雨养（即不灌） 四项。经过两层次二元组合，得出 14个决策变x1~