非恒定水流计算的最优控制问题及其变分求解.pdfVIP

第19卷第4期 水科学进展 V01．19．No．4 2008年7月 ADVANCESINWATERSCIENCE Jul．，2008 非恒定水流计算的最优控制问题及其变分求解 赖锡军1，傅国圣2，孙 波3 (1．中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室．江苏南京210008；2．扬泰水文水资源勘测局， 江苏扬州2251)00；3．山东省安丘市水利局．山东潍坊262100) 摘要：为综合利用各类可用的信息源(如水流运动规律，观测资料和统计信息)，在反问题的框架下构建了非恒定水 流计算的变分模型。针对模型预测的可能误差来源，给出了初始条件、边界条件、物理参数三类独立参量及综合控 制问题的定义。在此基础上，阐述了变分法求解这类偏微分方程最优化控制问题的基本原理及模型求解步骤。同 时，为便于应用参考，选取初始条件、边界条件和物理参数作为控制向量，以实际常用的水流运动微分方程为基 础。导出了一、二维非恒定水流计算的通用伴随方程。并以珠江黄浦至大虎段二维潮流计算为例，展示了变分模型 的应用。 关键词：非恒定流；反问题；变分法；最优控制；数据同化；数学模型 中图分类号：IVl33．2文献标识码：A 文章编号：1001-6791(2008)04．0537-09 尽管数值方法有了空前发展，但是非恒定流的分析与预测仍受三方面因素的制约：未知的初始条件，不完 整的边界条件以及模型自身的误差。一直以来，通过引入新的数值方法。计算精度得到了有效地提升，数学模 型的性能和对复杂问题的适应性有了明显改善。但是，数值方法的改进主要在于减少模型自身不完备的误差， 不能解决实际问题中存在的众多内在和外部的不确定性影响。为此，对实际计算问题，必须经过模型选取，初 边值和物理参数设定等预处理和历史实测资料的验证等一系列流程。这一经典的计算流程是按正问题(direct／ forward 架下，模型和观测资料两者是分离的。观测信息的利用程度受主观因素影响大。 观测资料蕴含着丰富的水流运动状态信息。自然地，如果有模型可以在运用动力方程分析时，有机地提取 观测资料中的有用信息用于水流分析与预测，那么就可进一步提高模型的计算能力，求出与实际最为吻合的 “真实”水流状态。基于这一想法，本文在反问题(inverse problem)的框架内，概要地论述了融合动力模型和观 测资料的非恒定水流反问题的计算模型，给出了实现客观地利用观测资料来改进数值计算的实际应用能力的实 现方法。反问题与正问题相对应，是“由果及因”的思路，即根据已知的结果，来反求未知参数，从而得到最 优的计算结果。本文的模型方法是基于偏微分方程最优化控制理论的变分同化方法(又称伴随同化法)，其在大 气…、海洋【2J、水文D,4J等一些自然科学领域已得到了发展应用。应用该方法，可有效地实现各类水文观测资 料的客观利用，识别一些不确定参数，从而可以更好地实现非恒定水流的分析与预测预报。 首先以概化的非线性方程组为例，将非恒定水流反问题的计算归结为最优控制问题，然后阐述了变分法求 解该类偏微分方程最优化控制问题的基本原理以及求解步骤。接着，导出了实际中常用的一维、二维非恒定水 流基本方程的伴随方程形式作为建模参考。最后，以珠江黄浦至大虎段潮流的计算为例展示了模型的应用。 收稿日期：2007．07．30 基金项目：国家自然科学基金资助项目40571028) 538 水科学进展 第19卷 1 非恒定水流计算的控制问题 一般地，非恒定水流计算所采用的微分方程组可表示成如下向量形式： 万3 U=r(v，P) (1) 式中t为时间；F为非线性算子；U为描述水流运动状态的向量，如水位、流速等；P为方程的参数向量。为 求解该方程组，还必须同时给定相应的初始条件UI。：o=矿和边界条件UIr=Ur。 对非恒定流计算