动态规划法对灌区水资源最优化分配.pdfVIP

维普资讯 安徽农业科学，JournalofAnhuiA ．Sci．2008，36 15 ：6147，6176 责任编辑 李玮 责任校对 张敏 动态规划法对灌区水资源最优化分配 王松林，周文 华北水利水电学院，河南郑州450o∞ 摘要 建立了各种作物全生育期的水分生产函数，以此为基础在农业种植结构不变的条件下，采用动态规划法进行有限水资源的合理 分配，以获得最大效益。 关键词 生产函数；动态规划法；合理分配；效益最大 中图分类号 $274．1 文献标识码 A 文章编号 0517—6611 2008 15—06147—01 我国北方地区由于降雨量少，水资源总量明显低于南方 产量即可得出灌水量与作物增产量之间的函数关系。 地区。随着工农业生产的发展，用水量增加，水资源短缺问 以相同农业生产投入为依据，根据农业生产过程中的投 题越来越严重。对于北方绝大部分灌区，一般只能对作物进 入产出关系，以不同频率年与灌水条件下的作物产出为对 行非充分灌溉。因此，研究在已定农作物种植 比例条件下如 照，利用如下公式即可取得不同频率年各种作物对应不同灌 何将有限的水资源分配给农作物，以获得最大效益，是灌 区 水量等级的增产效益及灌溉效益。 管理的主要研究内容之一。 B 一 变+ 不变 一 j0一 不变 一Yo 一变 1 水分生产函数的建立及参数推求 式中，B为由于灌水而获得的增产效益；k为作物价格 ；不变 为了进行水资源优化分配，达到最高效益的目标，必须 为单位面积上的固定投资，研究中一般将不同频率年和不 同 建立各种作物全生育期耗水量 ～产量函数。针对不 同频率 灌水量水平的 不蛮均取为定值；Y为单位面积不同灌水量的 年，在耗水量中扣除作物全生育阶段的有效降雨量和所分析 作物产量；Y。为单位面积不灌水条件下的作物产量；变为可 土层 内的土壤水变化量，即可获得作物对应频率年的灌水量 变投资。 ～ 作物产量函数，再利用作物耗水量 ～灌水量～产量 效益 2 灌区水资源优化分配 之间的间接关系，得到灌水量与净效益函数 ¨lj。 2．1 优化方法的选择 在灌区水资源管理中，采用的系统 1．1 水分生产函数的建立 在具体求解过程 中，首先建立 分析方法主要有线性规划、非线性规划、大系统优化、动态规 不同频率年灌水量～粮食产量函数关系式，然后利用投入产 戈等方法。在灌区有限水资源条件下的农作物种植 比例优 出关系建立灌水量与灌溉效益函数。对灌水量 ～灌溉效益 化多采用线性和非线性规划，多灌区间的水资源优化多采用 函数求导，可获得作物单位灌水量的效益变化率 边际效 大系统优化，而同一时期有限水资源优化多采用线性规划。 益 。根据边际效益，可推测单一作物有限水资源的优化灌 对于单一水库水资源优化调度、工程线路选择和有限水资源 溉方式。 在不进行农业种植结构调整的长时段，水资源优化分配多采 水分生产函数有线性关系，抛物线关系，幂指数关系等 用动态规划。 形式 2¨j。参照大量的研究成果，并分析以往灌溉试验的数 2．2 水资源优化分配计算 据，发现水分生产函数一般采用二次抛物线形式[引。 2．2．1 数学模型的建立。采用动态规划法进行水资源优化 1．2 水分生产 函数求解原理 设公式基本形式如下： 分配，首先应建立数学模型，其数学模型主要包括以下 5个 Y ax2+b +c 部分： 式中，y为单位面积产量 kg／hm2 ； 为单位面积耗水量 1 阶段变量。将每种作物作为一个用水单位，即看作 rn3／hm2 ；abc为回归分析系数。 — 个阶段 ，因此 ，具有 Ⅳ种作物的灌区，将共有 Ⅳ个阶段 若令Q y—Y1 乏 ax+bx+c—Y 要求a、b 变量。 2 状态变量和决策变量。状态变量为各阶段可用于分 和c的值 ，可以令 ：0、 ：0、 ：0，得到一个三元一次 配的有效水量，以 表示，决策变量为供给每种作物的水 方程组，该方程组的解即为所求，方程组形式如下： 量，即各阶段供水量以 表示。 ~Xi+b∑ 3+ C X 2 ~2 3 系统方程。根据状态变量和决策变量之间的关系， i Xi 1 i 1 i l 可推得系统方程为 ： +1 ％一 +b∑ 2l+c~Xi 锄 4 目标函数。设 为以总水量 分配给Ⅳ种作物而 i：1 l 获得的最大总净效益 ，则 ： + cn 式中，托、分别表示实测点的耗水量及其对应的作物产量； max ／n n 表示实测资料点数。在作物产量中扣除不灌水条件下的 5 约束条件。①供给各种作物水量之和不超过水资源 总量。②供给第 Ⅳ种作物 的水量 ％不能超过第 Ⅳ种作物 基金项 目 河南省科技厅攻关实验项 目 000291 。 在需供水阶段渠道最大过水能力、