糙率减小与洪峰增值关系的数值模拟研究.pdfVIP

泥 沙 研 究 2014年 12月 JournalofSedimentResearch 第 6期 糙率减小与洪峰增值关系的数值模拟研究 张小峰 ，元媛 ，‘赵梓涵 (1．武汉大学 水资源与水电工程科学国家重点实验室 ，湖北 武汉 430072： 2．山东黄河勘测设计研究院，山东 济南 250013) 摘要 ：以矩形断面河槽为研究对象 ，采用一维非恒定水流数学模型，模拟计算 了糙率变化对下游洪峰增值过程 的影响。结果表明：在河床断面不发生冲淤 ，断面形状沿程不变 ，且只考虑清水的情况下，当糙率不变时，随着 洪水逐渐 向下游传播 ，受河道槽蓄作用 的影响，洪峰流量沿程逐渐衰减 ，峰型趋于矮胖 ；当糙率随时间减小且 上游变化率大于下游河道 时 ，发生流量沿程异常增大的现象 ，且变化率越大 ，增值越明显 。通过理想情况下 的 单因素分析 ，认为糙率随时间沿程减小时，下游洪峰会发生增值现象。此外 ，进 口洪峰形态对洪水演进有一定 的影 响。峰型越矮胖 ，下游 断面的洪峰流量增幅越大 ；当入 口为连续洪峰时 ，因两峰叠加 ，洪峰异常增值有所 加剧 。 关键词 ：洪峰增值 ；糙率变化 ；洪水演进 ；数值模拟 中图分类号 ：TV143 文献标识码 ：A 文章编号 ：0468—155X(2Ol4)06-0036-08 根据水量平衡原理 ，洪水在宽浅河道演进时 ，受到河槽蓄滞洪作用的影响，洪峰逐渐坦化。近年来 ， 黄河下游洪水演进的过程中频繁出现洪峰流量异常增值的现象，增大了洪水分析和预测的难度 ，也对小 浪底水库调度运用造成 了诸多困难。关于洪峰异常增值 的成因，科研工作者进行 了多方面的探索和研 究。张红武…将 “92．8洪水 ”洪峰增值归因于漫滩洪水归槽。弓增喜 指出花 园口站 “05·7洪水”发生 洪水增值现象是刚刚结束的调水调沙使河道普遍刷深 、河床形态变优造成 的。崔 占峰 ，张小峰等 发 现在含沙量较大的情况下 ，水流连续方程里的浑水与床面的交换项对流量影响是明显 的，解释 了在象黄 河这样的多沙河流 中存在的 “水变沙”、“沙变水”的现象。Cao 认为河床强烈冲刷会导致洪水峰值 的 迅速增大。丁赞等 等根据浅水流动方程的基本特性 ，研究 了水流运动与河床变形耦合时洪水传播 的 特性 ，分析了河床变形对黄河高含沙水流增值的影响。在河道阻力方面，江恩惠等 ’探讨了洪峰增值 的水动力学机理 ，认为河道糙率变小是造成洪峰异常增值的根本原因，并采用能够合理反映含沙量及泥 沙级配变化对河道糙率影响的黄河河道准二维水沙数学模型进行 了数值模拟 ，验证 了上述观点。李国 英 通过研究发现，在高含沙洪水演进过程 中，由于含沙量 、悬沙级配 、床沙级配、床面形态、河道断面 形态及河岸周界情况等因素使得河道糙率大幅度减小 ，是产生洪峰增值的主要原因。 上述研究从不同角度对黄河洪水异常增值 的发生过程或者发生机理进行了探索。其中由于水流与 泥沙的共同作用导致综合糙率不断减小 ，引发洪峰异常增值的观点研究较多。但 目前对糙率引发增峰 观点多采用实测大断面资料进行河床冲淤变化和洪水复演 ，未排除河道冲淤变化等因素。本文基于一 维洪水演进数学模型，采用沿程不变的矩形断面进行概化计算 ，研究清水条件下，糙率单一因素以不 同 方式变化对洪水演进过程的影响 。 1 一维非恒定流数学模型介绍 描述河道一维非恒定水流运动的基本方程为圣维南方程组 ，包括水流连续方程和水流运动方程。 收稿 日期：2014-05429 基金项 目：国家 自然科学基金项 t：l 作者简介 ：张小峰(1962一)，男 ，浙江省嵊州人，教授 ，主要从事水力学及河流动 力学方面研究 。Em·ail：zhangxfwuhee@263．net 36 水流 连续 方程 一 oQ+AO ：0 (1) Ox Ot 水流运动方程 Qo