应用最优控制理论自动率定二维浅水方程糙率参数.pdfVIP

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