强冲积河流过程二维水沙耦合数值模拟-中科院寒旱所遥感与地理信息.PDF

中国科学 G 辑:物理学 力学 天文学 2008 年 第 38 卷 第 8 期 : 962 ~ 972 SCIENCE IN CHINA PRESS 强冲积河流过程二维水沙耦合数值模拟 ① ①* ② ② 岳志远 , 曹志先 , 李新 , 车涛 ① 武汉大学水资源与水电工程科学国家重点实验室, 武汉 430072; ② 中国科学院寒区旱区环境与工程研究所, 兰州 730000 * 联系人, E-mail: zxcao@ 收稿日期: 2007-06-11; 接受日期: 2008-04-18 国家重点基础研究发展计划(编号: 2007CB714106)、国家自然科学基金(批准号: 和中国科学院知识创新工 程重要方向项目(编号: KZCX3-SW-357-02)资助 摘要 强冲积河流过程泥沙运动非常活跃、河床变形快, 与水流之间存在 关键词 强烈的相互作用. 传统的基于简化控制方程的非耦合数学模型违背了基本 冲积河流 守恒律, 只能近似地适用于弱冲积河流过程. 建立普遍适用于强、弱冲积过 洪水 程的二维水沙耦合数学模型, 将现有对不可冲刷床面浅水二维流动的、可以 耦合数学模型 泥沙运动 捕捉激波和接触性间断的 WAF TVD 二阶数值方法扩展至可冲刷床面浅水 冰湖溃决洪水 二维水沙运动问题. 应用该耦合模型研究了典型冰湖溃决洪水过程. 近几十年发展了大量的冲积河流数学模型并被广泛应用于研究河流工程、环境、生态与 灾害问题. 但是, 现有模型主要建立在水沙非耦合理论基础之上, 只能近似地适用于输沙强度 小、河床变形很慢的弱冲积过程. 然而, 强冲积过程在自然界广泛存在, 其水流急剧变化, 往 往诱发非常活跃的泥沙运动和快速河床变形. 高含沙洪水经常发生在中国的黄河以及孟加拉 国和巴拉圭等国的一些多沙河流中, 其河床变形速率与水深变化率可能为同一数量级[1], 是一 类典型的强冲积过程. 冲积河流大坝溃决(或拆除)洪水能量大, 必然诱发非常活跃的泥沙运动 与显著河床变形1), 也是典型的强冲积过程(如 1975 年 8 月特大暴雨导致中国河南省板桥、石 漫滩等水库大坝溃决). 冰湖溃决洪水(GLOF, glacier lake outburst flood)通常发生在高原地区陡 峭坡面上[2~6], 水流强度大、侵蚀能力强, 可能急剧冲刷坡面, 诱发泥石流灾害, 伴随着全球气 候变暖, 许多冰湖具有潜在的溃决危险[7], 构成一类典型的强冲积过程. 强冲积过程水流、泥沙与床面之间存在强烈的相互作用, 传统的非耦合数学模拟理论忽略 了其基本力学特征, 其简化的控制方程违背了流体力学守恒定律, 无法适用于强冲积过程. 因 此, 需要建立基于完整守恒律和高性能数值格式的水沙耦合数学模拟理论. 国内有学者建立 1) Zech Y, Spinewine B. Dam-break induced floods and sediment movement-state of the art and need for research. In: First Workshop of EU Project IMPACT. HR Wallingford, 2002 962