一般形式的挟沙力公式及其分析.pdfVIP

泥 沙 研 究 2014年 2月 JournalofSedimentResearch 第 1期 一 般形式的挟沙力公式及其分析 李瑞杰 ，郑俊 ，张文虎 ，曹宏生 (1．河海大学 海岸灾害及 防护教育部重点实验室 ，江苏 南京 210098；2．国家海洋局 海域管理技术重点 实验室，辽宁 大连 116023；3．河海大学 环境海洋实验室，江苏 南京 210098) 摘要 ：将浑水潜流动量方程延拓到整个含沙水流水深 的范围，在稳定条件下，推 出平衡含沙量 的表达式，在此 基础上采用水流最小能耗率原理，得到一般形式的挟沙力公式，挟沙力系数与阻力系数成正 比。论证 了在不 同条件下，公式可以与刘家驹、韩其为、维利卡诺夫、窦国仁和张瑞瑾等公式的形式一致 ，并指出各公式表达形 式的不同是由于阻力系数与水流雷诺数的不同依赖关系导致的。从理论上说明了各公式的适用范围，其中张 瑞瑾公式中的指数范围应为0．67—1，韩其为公式指数 11／12即0．92对应为紊流过渡 区水流流态 ；刘家驹公 式适用于淤泥质底床 ，且其系数不是常数 ；维利卡诺夫公式和窦 国仁公式适用于紊流粗糙 区水流 ，维利卡诺夫 公式的系数是相对糙率 的函数 ，与雷诺数无关。通过河流及近岸实测资料对文 中导出的公式进行 了验证 ，结 果表明计算值与实测值吻合很好 。 关键词 ：挟沙力；阻力系数 ；雷诺数 ；紊流 中图分类号 ：TV143．4 文献标识码 ：A 文章编号 ：0468—155X(2014)01—0001 7 1 前言 挟沙力问题是泥沙运动力学的关键问题之一 ，其确定对于挟沙水流的泥沙输运 、底床冲淤以及水底 边界条件都具有至关重要的意义和价值，也关系到河流、河 口和海洋泥沙运动及底床冲淤演变预测的精 度和成败。为此很多研究人员对挟沙力问题从不 同角度进行 了探讨和研究 ，取得了大量的研究成 果 ，促进了对挟沙力问题认识的提高。 关于挟沙力与各物理量的关系已经进行过大量的探讨 ，基本方法包括经验分析、因次分析和能量平 衡等方法。由于影响挟沙力的各种要素的多样性和相互关系的复杂性 ，以及认知方式、研究方法 、水域 特征、资料来源的差异 ，导致现有的挟沙力公式种类繁多、结构各异。仅从动力影响因子流速来看 ，指数 从 0．0～4．0范围存在着大量不同的用于不同水域的量化公式 ，这还不包含泥沙本身特性 的相关因子 的 影响。目前被人们普遍认可的挟沙能力统一公式并不多见 ，从理论上严格导出挟沙力公式并解释挟 沙力公式形式多样的原因，还没有见到报道。 本文从浑水潜流动量方程 出发，根据最小能耗率原理 对挟沙力 的一般形式进行 了推导和分 析 ，并根据推导所得挟沙力公式，解释现有挟沙力公式 的多样性。最后根据 Grass¨ 的方法 ，将所得一 般关系推广应用到海洋挟沙力的计算公式，并用实测资料对推得的公式进行了验证和 比较 ，其结果与河 流和海洋中的实测数据吻合很好。 2 一般形式的挟沙力公式 Parker等 “根据流体质量守恒和动量守恒定律导出了浑水潜流方程 中的动量平衡方程为 收稿 日期 ：2013-01-21 基金项 目：国家 自然科学基金项 目 ；国家海洋局海洋公益性行业科研专项经费(201205005，201105006) 作者简介 ：李瑞杰(1963一)，男 ，山东青岛人 ，博士 ，博士生导师 ，主要从事物理海洋学 、河 口海岸物质输运研究 。 E—mail：qli@ hhu．edu．en OcD2 + 一 g +RgcDS u (1) 式中u为垂向平均流速 ；D为浑水潜流层厚度 ；R=p／o一1，p为泥沙颗粒密度 ，p为水体密度 ；g为重 力加速度；c为垂向平均体积含沙量 ；S为底坡 ，S= ，对于大多数河流情况 S远小于 1，叼为床面坐 标 ；C。为阻力系数； 为沿