明槽水流的非均匀沙挟沙力研究.docxVIP

明槽水流的非均匀沙挟沙力研究 摘要 本文从床沙、推移质及悬移质 互相交换、衔接的物理图形出发，在力学、 随机及紊动扩散分析基础上， 提出了颗粒跃 移及推移质、悬移质和全沙统一的新的非均 匀沙挟沙力公式.文中公布了新近试验取得 的57组非均匀沙挟沙力水槽资料，分析提 出了反映混合沙中粗细颗粒互相影响的系 数的关系，提出了悬移层与推移层交界面高 度及浓度预报新的初步关系.应用本文提出 的挟沙力系列公式计算验证了大量天然河 流及水槽实测分组和总挟沙力资料，计算与 实际符合较好. 关键词非均匀沙颗粒跃移悬移质挟沙力 隐蔽暴露影响 水流挟沙力是泥沙研究中的核心问题， 是解决水库、河渠等大量泥沙冲淤问题的关 键.天然冲积河流中的泥沙一般都是非均匀 沙，悬移质、跃移质及床沙通常都是互相交 换、密切相关的，但迄今国内外大多数挟沙 力公式从物理图形、原始推导上就与床沙组 成无关，实际上是针对且只适用于均匀沙 . 一些公式经过各种经验性处理以后已用于 非均匀沙计算，在生产上得到了广泛应用， 但一般都未考虑混合沙中粗细颗粒间的相 互影响，分组挟沙力计算精度难以提高，在 河床连续冲刷计算中问题暴露往往比较明 显.迄今，多数挟沙力公式将悬移质和推移 质分割，生产上应用的多是并不匹配的二套 公式，对于悬移质和推移质输沙并重的河流， 如此作法存在不少问题.水流的非均匀沙挟 沙力仍是一个远未解决的重大生产和学术 问题. 等很多学者都先后进行了颗粒跃移的 力学分析或试验，取得了很多宝贵成果.自 1986年以来，王士强、张仁、付晓及胡春宏、 惠遇甲等分别对明槽水中颗粒跃移进行了 力学分析和试验，取得了很大进展:4-6 : 分别发现颗粒容重对跃移参数的影响规律， 与以往所有认识均不相同.在颗粒跃移分析 和试验基础上，王士强等提出了计算推移质 输沙率的新公式:4］,进而提出了推移质、 悬移质和全沙统一的新的非均匀沙挟沙力 公式［7 ］,并得到了清华大学活动大钢槽内 试验取得的27组实测资料的良好验证［8］, 文献:9］对此公式进行了进一步的分析和 验证，并与国内外其它一些着名的挟沙力公 式进行了验证计算精度比较. 近些年来，王光谦、倪晋仁:10］，邵 学军、夏震寰［11］及胡春宏、惠遇甲［6］ 等用不同方法探索研究了近底层内颗粒浓 度分布规律，为挟沙力研究的进一步深入提 供了宝贵的启示；王士强等考虑近壁层颗粒 沉降阻力增大、沉速减小从而得出新的悬移 层浓度垂线分布公式［12］,并进一步总结、 验证和应用了新公式.近二年内陈骥完成了 新的57组非均匀沙挟沙力试验，在所有以 上研究基础上，本文对更多的资料进行了分 析，继续研究改进了王士强等新的挟沙力系 列公式，进一步提高了该套公式的计算精度 一、颗粒跃移及推移质、悬移质、全沙挟沙 力方程 颗粒跃移方程 床面附近的颗粒跃移 是悬移质输沙和沙质床面之间过渡层内颗 粒的主要运动形式，是研究沙质推移质和悬 移质输沙率的基础.作用在床面附近跃移颗 粒上的力很多很复杂，文献:4]主要考虑 了占优势影响的颗粒重力(W)、水的浮力 (W‘ -W,W为水下重力)、水流上举力(FL) 及水流拖曳力(FDx)或 阻力(FDy),其它作用力因较为次要或不很 明确而未直接考虑.垂直床面(或水流)的y 向及沿水流x向的颗粒跃移运动的微分方程 分别为 FL-W -FDy=(Wd2y)/(gdt2) (1) FDx=(Wd2x)/(gdt2). (2) 式中g为重力加速度，t为时间.通过适当简 化对上列方程积分，得出了各跃移参数解析 式，由跃移试验资料适线确定了这些解析式 中的系数:4].下列简化公式是从较复杂的 解析式概括得出的，二者基本等价，计算结 果十分相近.式中D为颗粒粒径，a0为颗粒 平均跃移高度，Lm为平均跃移长度，ub为 平均跃移速度，u*为水流摩阻流速，Y S、 丫分别为颗粒和水的容量，J为能坡，R 为沙粒水力半径，按文献:14］给出公式求 得，平整床面时即为水力半径或水深 . (3) (4) (5) 式中 上列颗粒跃移参数公式得到了验证:5］ 与其它着名公式相比:9］,数值相近但有二 点重要区别，上列公式中的相对跃移参数不 仅取决于6),还另外随y S/ y而变，同时 ?的籍指数并非常数.对于泥沙，当?从增 至时，ao/D由1增至11,Lm/D由6增至800, ub/u*由4增至12,尤以跃长变化范围最大， 若如Einstein假设为常数100,则对于其它 大部分水流强度情况必然产生很大误差 . 推移质输沙率公式 颗粒以跃移形式 运动的单宽推移质输沙率 qb应等于单位时 间、单位宽度内跳跃通过某断面的各种跃长 的所有颗粒的总重量，可由下式［4］确定： qb=A* 丫 spubLm (6) 式中P为FL/W 1条件下床面颗粒的起跃