应用最优控制理论自动率定二维浅水方程的糙率参数-国家自然科学.PDF

维普资讯 第 19卷 第 3期 水 科 学 进 展 V01．19。No．3 2008年 5月 ADVANCESIN rAIER SCIENCE May，2008 应用最优控制理论 自动率定二维浅水方程的糙率参数 赖锡军，姜加虎，黄 群 (中国科学院南京地理与湖泊研究所湖泊与环境 国家重点实验室，江苏 南京 210008) 摘要：应用偏微分方程最优化控制理论，导出了二维浅水方程糙率 自动率定的最优化系统，实现了二维水力学模型 计算中分布式糙率参数的自动率定算法。结合实际应用的可操作性 ，按区块设定糙率，对恒定流和非恒定流两种流 态分别进行了常数糙率和随水深变化的糙率 自动率定实验。数值实验表明，在观测数据包含足够的信息时，该方法 不仅能够准确地 自动率定出相关的糙率参数，还能成功地将糙率随水深的变化关系式识别出来。但是，糙率水深变 化关系识别迭代收敛速度显著变慢，因此，建议采用相应的实现对策来考虑变糙率的洪水预测 问题。 关 键 词：糙率；自动率定；参数识别；最优化控制 ；浅水方程 中图分类号：TVI33．I 文献标识码 ：A 文章编号：1001．6791(2008)03．0383．06 计算水力学从 20世纪 70—80年代兴起 ，以有限差分格式为主 【；20世纪 90年代以来，基于有限体积方 法的新型高分辨率格式的建立和迅速发展／2J，使应用数值手段研究复杂特征水流 (如溃坝、水跃等)的能力也有 了质的飞跃。水力学模型应用，要给定各种参数，如糙率、边界条件等，当这些参数准确时可以获取理想的结 果 。其中，糙率值是一个间接量 ，不能直接测得 ，实际应用通常是通过不断地调整数值 ，将计算结果与水位、 流量等观测资料进行 比较 ，从而取定一个较好的值。这个过程带入了很多的主观性 ，得到的并不是最优解 ，而 且具有很强的经验性。 人工率定过程通常通过 比较一些观测点的计算和观测值差异 ，人为调整相关参数 ，力图设定最 “优”的糙 率值 ，使得差异最小 ，这实际上是人脑对计算系统的识别过程 。当问题变得复杂时，很难估出准确的参数。偏 微分方程的最优控制理论l3J为这类受控于动力方程的系统的最优化控制问题提供了数学理论基础。其基本思想 是引入代价函数，借助伴随方程计算代价函数关于参数的梯度 ，利用最小化算法 自动调整参数使代价函数最小 化 ，从而获取最优的参数解。利用该方法 ，本文实现了二维分布式糙率参数的自动率定(亦称为识别)。在二维 浅水方程基础上导出了糙率参数 自动率定的伴随方程 ，利用建立的模型开展相关 的率定实验。 1 水力学模型的基本方程 以二维浅水方程作为基本方程 的水力学模型在环境、水利等涉水领域的科学和工程问题的研究中得到了广 泛的应用。其常用的守恒型方程可写为 瓦 ++— 一 ++— 一 ：丑(u) (l1) 式 中 、Y为笛卡儿直角坐标系；t为时间；U为方程的状态变量 向量 ；F、G分别为 、Y方 向的通量向量 ；B 为源项 向量。表达式如下 ： U = (h，hu，hv) (2) F：(，h+专g，删) (3) G：( 聊， +告g) (4) 收稿 日期 ：2007．06．15 基金项 目：国家 自然科学基金资助项 目 40571028)；江苏省基础研究计划 (自然科学基金)(BK2007611) 作者简介：赖锡军(1977一)，男，浙江遂 昌人，助理研究员，博士，主要从事环境水力学研究。Em·sJlzxjlai@mg~ ．ac．cn