应急平台中一维洪水演进模型研究.pdfVIP

维普资讯 第 33卷 第 5期 武 汉 大 学 学 报 ·信 息 科 学 版 V0l_33No．5 2008年 5月 GeomaticsandInformationScienceofWuhanUniversity M ay 2008 文章编号：1671—8860(2008)05—0542—04 文献标志码 ：A 应急平台中一维洪水演进模型研究 吴迪军 孙海燕 黄全义 。 熊金宝 (1 武汉大学测绘学院 ，武汉市珞喻路 129号 ，430079) (2 中铁 大桥勘测设计院有限公司，武汉 市汉 阳大道 34号，430050) (3 清华大学工程物理系 ，北京市海淀 区清华园，100084) 摘 要：根据应急平 台的特点和需求，建立 了一维洪水演进数学模型。利用简便易行的基 于特征线法的库 朗 格式求解方程 ，应用 Matlab编制 了洪水演进计算程序 ，并利用其强大的绘 图功能实现 了计算结果的可视化显 示和分析 。通过对长江某河段实测资料 的建模分析和应用研究，验证 了模型的实用性和可靠性 。 关键词 ：应急平台；一维洪水演进模 型 ；Matlab 中图法分类号：P208 一 维洪水演进 的计算方法主要包括两大类 ， 水流速；R为水力半径 ；C为谢才 系数 ；g为重力 一 类是水文学方法 ]，另一类是水力学方法 ]。 系数 ；q为单位长度旁侧人流。 随着现代计算机存储技术和运算速度的飞速发 上述方程组为一双 曲型偏微分方程组，可采 展 ，目前多采用水力学方法 ，而其中又以具有无条 用特征线法求出其数值解。由于无法得到两个常 件稳定特性的Preissmann隐式差分格式最为常 微分方程的解析解 ，因此 ，通常将两个常微分方程 用 ]。国内外学者对一维洪水演进模型的研究 各变为有限差分形式方程 ，再根据给定的初始条 主要集 中在水利领域 ，且更多地侧重于计算精度 件及边界条件求得任何时刻每一河道 断面上的 方面的研究，而基于应急平台的集成研究尚不多 、 7d或Q 的数值解 。 见。因此 ，本文结合应急平台中洪水演进计算 的 如图 1所示 ，在需求解 的 ￡平面上沿时间t 快速、高效、实时、可视化查询和显示等特殊要求， 坐标按时间步长 △￡划分若干段，并以 表示 时间 对一维洪水演进模型的建立方法及其关键技术进 分段序号；沿流程 5划分若干个断面，以 i表示各 行了研究。 个断面序号 ，△s一5… 一s为流程距离步长 。矩 形网格的交点称为结点，其坐标可写为 (s，t，)，简 1 一维洪水演进数学模型的建立 写为(i，)，相应的因变量为z(s，t)、Q( t)，简 写为 、 。 1．1 一维水流运动方程及其求解 j+2 沿河道 的一维洪水演进计算的基本方程一般