黄河下游高含沙洪水演进及河床冲淤过程的数值模拟.pdfVIP

泥 沙 研 究 年 月 第 期 黄河下游高含沙洪水演进及河床冲淤过程的数值模拟 贺 莉 夏军强 王光谦 清华大学水沙科学与水利水电工程国家重点实验室 北京 摘要 高含沙洪水通过黄河下游河道时往往发生严重淤积 因此 数值模拟研究高含沙洪水的演进过程 对理 论研究及实际工程需求均具有重要意义 以往的泥沙数学模型通常不考虑高含沙量对水流控制方程的影响 而这将会影响高含沙洪水演进过程的模拟 本文建立了黄河下游高含沙洪水演进过程的一维非恒定非均匀 沙数学模型 该模型不仅在水流控制方程组中采用浑水的连续方程与动量方程 而且通过引入滩槽划分及 二 级悬河 处理等技术考虑黄河下游复杂地形下的水沙演进 以 年的典型高含沙洪水过程为例 分析高含 沙量对洪水演进的影响 并率定出模型中的恢复饱和系数 不饱和系数等关键参数 将率定后的模型参数应 用于 年下游实测洪水的模拟 验证模型的适用性 计算结果表明 考虑高含沙的影响后 洪峰流量 水 量 传播时间的模拟误差减小 说明在数值模拟高含沙洪水的演进中有必要采用浑水控制方程 关键词 黄河下游 高含沙洪水 浑水连续方程 浑水运动方程 洪水演进 河床冲淤 中图分类号 文献标识码 文章编号 概述 黄河最突出的问题是水少沙多 导致下游河道不断淤积 抬高 防洪形势日趋严重 而发生在黄河 中下游的高含沙洪水 是造成下游河道严重淤积的重要原因之一 据统计 至 年期间 三门峡 站洪水中共有 次洪水的瞬时最大含沙量超过 ? 而这 次高含沙洪水在下游共淤积泥沙 亿 占来沙量的 占 至 年下游河道总淤积量的 因此 数值模拟高含沙洪水 的演进过程及河道冲淤特性 具有重要的理论意义 一维泥沙数学模型因其所需资料相对较少 计算速度快 且计算精度能满足实际工程要求 在模拟 长河段 长时段的河床变形过程时仍发挥着重要作用 且理论和实践上都比较成熟 在实际工程应用中 较多 目前 国内已开发了很多适用于低含沙水流的泥沙数学模型 适用于高含沙水流的数学模型 也得到了较大发展 如张红武模型采用修正的泥沙连续方程及河床变形方程 引入其提出的水流挟 沙力和河床阻力公式 但是目前大多数模型都没有考虑高含沙量对洪水演进的影响 即取浑水密 度 近似等于清水密度 该假定一般适用于低含沙水流的情况 若设 则浑水运动方程可化 为定床中的运动方程 与清水运动方程完全相同 而这种假设在高含沙量条件下不成立 因此在模拟 高含沙洪水的演进时需要考虑采用浑水控制方程组 本文在此建立一维非恒定非均匀沙数学模型 不仅在水流控制方程中采用了浑水连续方程与动量 方程 在泥沙输移方程中则选取考虑了高含沙影响的挟沙力等公式 以充分考虑高含沙量对泥沙的沉降 特性 粘性等的作用 此外 为适应黄河下游复式断面中高含沙洪水的演进过程 还引入滩槽划分及 二 级悬河 处理等技术 首先以 年的典型高含沙洪水过程为例 分析高含沙量对洪水演进过程中峰 收稿日期 基金项目 创新群体研究基金 国家自然科学基金 作者简介 贺莉 女 博士研究生 主要研究方向为河流动力学 值 传播时间等的影响 并率定出模型中恢复饱和系数 不饱和系数等关键参数 然后将率定的参数应用 于 年发生的洪水 对模型进行进一步的验证 计算结果表明 高含沙量对水沙演进过程中的传播 时间 洪峰流量 水量及沙峰均有一定影响 总的来说 考虑高含沙影响后 洪峰流量 水量 传播时间的 模拟误差减小 模拟结果与实测更为符合 模型更能用于黄河下游高洪水演进过程的模拟 因此 在模 拟高含沙洪水演进时有必要采用浑水控制方程 模型控制方程和数值解法 控制方程 一维非恒定