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

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