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

1、明槽水流的非均匀沙挟沙力研究    摘要：本文从床沙、推移质及悬移质互相交换、衔接的物理图形出发，在力学、随机及紊动扩散分析基础上，提出了颗粒跃移及推移质、悬移质和全沙统一的新的非均匀沙挟沙力公式.文中公布了新近试验取得的57组非均匀沙挟沙力水槽资料，分析提出了反映混合沙中粗细颗粒互相影响的系数的关系，提出了悬移层与推移层交界面高度及浓度预报新的初步关系.应用本文提出的挟沙力系列公式计算验证了大量天然河流及水槽实测分组和总挟沙力资料，计算与实际符合较好. 关键词：非均匀沙 颗粒跃移 悬移质 挟沙力 隐蔽暴露影响 水流挟沙力是泥沙研究中的核心问题，是解决水

2、库、河渠等大量泥沙冲淤问题的关键.天然冲积河流中的泥沙一般都是非均匀沙，悬移质、跃移质及床沙通常都是互相交换、密切相关的，但迄今国内外大多数挟沙力公式从物理图形、原始推导上就与床沙组成无关，实际上是针对且只适用于均匀沙.一些公式经过各种经验性处理以后已用于非均匀沙计算，在生产上得到了广泛应用，但一般都未考虑混合沙中粗细颗粒间的相互影响，分组挟沙力计算精度难以提高，在河床连续冲刷计算中问题暴露往往比较明显.迄今，多数挟沙力公式将悬移质和推移质分割，生产上应用的多是并不匹配的二套公式，对于悬移质和推移质输沙并重的河流，如此作法存在不少问题.水流的非均匀沙挟沙力仍是一个远未解决的重大生产和学术问题.

3、 H.A.Einstein1-2早在五十年代即提出了水流的非均匀沙挟沙力理论及公式，半个世纪以来的研究和实践表明，爱氏提出的总的物理图形比较合理.尤其是将床沙、跃移质和悬移质三者联系起来分析而不是将它们看成彼此孤立、分割的事物，这一点是意义深远的.这一物理图形中一个关键的环节是床面与悬移层之间过渡薄层(跃移层或准跃移层)内的水流泥沙运动规律，由于此近底层内规律非常复杂且试验困难，致使爱氏体系的理论发展比较滞缓，爱氏当初推导挟沙力关系时不得不作了很多假定，不少是有问题的.爱氏对颗粒跃移并未进行详细的力学分析和试验，诸多关键的跃移参数相对值均以常数代之，与实际相差甚远，使其挟沙力公式计算成果与实际

4、往往相差很大. R.A.Bagnold, M.S.Yalin及Van Rijn3等很多学者都先后进行了颗粒跃移的力学分析或试验，取得了很多宝贵成果.自1986年以来，王士强、张仁、付晓及胡春宏、惠遇甲等分别对明槽水中颗粒跃移进行了力学分析和试验，取得了很大进展4-6，分别发现颗粒容重对跃移参数的影响规律，与以往所有认识均不相同.在颗粒跃移分析和试验基础上，王士强等提出了计算推移质输沙率的新公式4，进而提出了推移质、悬移质和全沙统一的新的非均匀沙挟沙力公式7，并得到了清华大学活动大钢槽内试验取得的27组实测资料的良好验证8，文献9对此公式进行了进一步的分析和验证，并与国内外其它一些著名的挟沙力公

5、式进行了验证计算精度比较. 近些年来，王光谦、倪晋仁10，邵学军、夏震寰11及胡春宏、惠遇甲6等用不同方法探索研究了近底层内颗粒浓度分布规律，为挟沙力研究的进一步深入提供了宝贵的启示；王士强等考虑近壁层颗粒沉降阻力增大、沉速减小从而得出新的悬移层浓度垂线分布公式12，并进一步总结、验证和应用了新公式.近二年内陈骥完成了新的57组非均匀沙挟沙力试验，在所有以上研究基础上，本文对更多的资料进行了分析，继续研究改进了王士强等新的挟沙力系列公式，进一步提高了该套公式的计算精度.一、颗粒跃移及推移质、悬移质、全沙挟沙力方程1.1颗粒跃移方程床面附近的颗粒跃移是悬移质输沙和沙质床面之间过渡层内颗粒的主要运

6、动形式，是研究沙质推移质和悬移质输沙率的基础.作用在床面附近跃移颗粒上的力很多很复杂，文献4主要考虑了占优势影响的颗粒重力(W)、水的浮力(W-W,W为水下重力)、水流上举力(FL)及水流拖曳力(FDx)或阻力(FDy)，其它作用力因较为次要或不很明确而未直接考虑.垂直床面(或水流)的y向及沿水流x向的颗粒跃移运动的微分方程分别为 FL-W-FDy=(Wd2y)/(gdt2)(1) FDx=(Wd2x)/(gdt2).        (2)式中g为重力加速度，t为时间.通过适当简化对上列方程积分，得出了各跃移参数解析式，由跃移试验资料

7、适线确定了这些解析式中的系数4.下列简化公式是从较复杂的解析式概括得出的，二者基本等价，计算结果十分相近.式中D为颗粒粒径，a0为颗粒平均跃移高度，Lm为平均跃移长度，ub为平均跃移速度，u*为水流摩阻流速，s、分别为颗粒和水的容量，J为能坡，R为沙粒水力半径，按文献14给出公式求得，平整床面时即为水力半径或水深.         (3)      (4) (5)式中 上列颗粒跃移参数公式得到了验证5，与其它著名公式相比9，数值相近但有二点重要区别，上列公式中的相对跃移参数不仅取决于，

8、还另外随s/而变，同时的幂指数并非常数.对于泥沙，当从0.06增至2.5时，ao/D由1增至11，Lm/D由6增至800，ub/u*由4增至12，尤以跃长变化范围最大，若如Einstein假设为常数100，则对于其它大部分水流强度情况必然产生很大误差.1.2推移质输沙率公式颗粒以跃移形式运动的单宽推移质输沙率qb应等于单位时间、单位宽度内跳跃通过某断面的各种跃长的所有颗粒的总重量，可由下式4确定： qb=A*spubLm               

9、60;               (6)式中P为FL/W1条件下床面颗粒的起跃概率，A*为待定系数.从物理模式及原始推导看，严格来说，式(6)本只适于跃移推移质，但因以滚动、滑动等形式运动的接触质受到的主要作用力与上述跃移颗粒大体相同，因此控制qb的变量及s/二者应该相同.当式(6)中的待定系数由包括接触质在内的变化范围很大的大量推移质输沙率实测资料适线确定时，该式的适用范围就具有普遍性.式(6)中的变量可由一套解析公式分别求出，相当烦复，为便于应用，由

10、式(6)及其它一系列求其变量的公式概括简化得出公式(7)，两者基本等价，计算值十分相近.                         (7)式中：. 图1列出了大量实测点据及公式(7)中不同s/的三条曲线，公式与各种s/的实际资料均相当符合，在高强度输沙部分，与Wilson的二种容重颗粒资料均分别符合相当好.图1还列出了其它著名公式曲线以便互相比较，对于泥沙，公式

11、(7)与Yalin公式十分接近，幂指数m3随增大而减小，高强度时接近1.5. 公式(7)得到了其它资料的验证9. 图1 和、s/的关系1.3悬移质及全沙挟沙力公式水流的悬移质输沙一方面依赖于紊动扩散，同时取决于床沙扬起及跃移泥沙浓度.水流的悬移质单宽输沙率qs由式(8)确定：   (8)式中sy及uy分别为离床面y高度的含沙浓度及水流速，H为水深，a为悬移层底层离床面高度.式中流速垂线分布按Einstein总结的公式计算13： uy=5.75u*log(30.2x1y/D65) .(9) 含沙浓度分布sy由紊动扩散方程并考虑床面附近沉速减小的因素确定，可由作者提出的式(10)计算12

12、，该式与通常的公式形式相同，仅以修正悬浮指数z1代替通常的悬浮指数z，z1随z增加及粒径D的减小而比z小得愈多.    z=0/u* .式中为卡门常数，0为泥沙沉速，Sa为y=a处的含沙浓度.    将式(9)、(10)代入式(8)积分可得式(11)为：qs=11.6su*Saa(P1I1+I2) ,(11)式中： P1=2.303lg(30.2x1H)/(D65) , I1=0.216(Az1-1)/(1+A)z1)A(1-y)/(y)z1dy , I2=0.216(Az1-1)/(1+A)z1)A(1-y)/(y)z1lnydy

13、. 设Sa=KsSm, sm=qb/(aubs), sm为推移层内的平均含沙浓度，Ks为近壁层浓度分布特性系数.则式(11)可表达为 qs=11.6Ksqbu*/ub(P1I1+I2), (12)全沙输沙率qT为： qT=qb+qs=qb(1+11.6Ks(P1I1+I2)u*/ub) .(13)对于非均匀混合沙，每个粒径组的水流挟沙力与床沙中该粒径组所占比重io相对应，当不考虑粗细颗粒互相影响时，则 iTqT=ibqb+isqs=ibqb(1+11.6Ks(P1I1+I2)u*/ub), (14)式中： ibqb=i0Dg, io,iT,ib,is分别为该粒径组在床沙及水中全沙、推移质、悬移

14、质输沙率中所占重量百分比.式(14)中的所有粒径D均为混合沙内各粒径组的平均粒径. 上列公式中的a值可由下述方法求得.因颗粒平均跃移高度ao的计算公式(3)中的系数是根据水槽颗粒跃高试验资料适线确定的，因此在水槽输沙计算中取a=a0,这是合理、自然的，如此计算出的悬移层底部浓度及挟沙力与大量水槽内的实测资料符合相当好. 对于天然河流、水库，无任何实测跃移层高度a的资料，若仍取a=a0,则计算出的悬移层底部浓度一般都远大于实测资料.因天然河流、水库内的水深远大于水槽情况，作为水深方向的物理量a，可以推测和假定随相对水深H/D50而增大，根据一些天然河流、水库实测底部含沙量资料分析反算，得出一般含

15、沙量情况下天然河流内的a可按经验公式(15a)计算，a可称为准跃移层厚度： (15a)对于长江及黄河大部分非高含沙水流情况，按式(15a)计算的悬移质挟沙力与大量实际资料符合较好.但对于黄河少数高浓度情况，若按式(15a)计算又计及沉速修正后得出的悬移层底部含沙量仍明显小于实际情况. 从上述公式(3)可以看出，跃移层厚度随水流容重增大而减小，在高含沙水流情况下，随浓度增加而明显增大.另外，随着浓度增加，对数流速公式中的卡门常数将会减小，即由床面向水面方向的流速变大的梯度增大，床面附近水流垂向紊动强度向水面方向变大的梯度也理应增大，这就意味着相同平均水流强度时含沙浓度愈高，颗粒由跃移转为悬移的高

16、度a离床面愈近.根据以上初步分析及实际资料反算，在公式(15a)中增加卡门常数的影响因子，得出了包括高浓度水流情况的计算a的初步经验公式(15)  . (15)式中卡门常数据一些实测资料13可按式(16)计算.=0.03421-0.24153E-0.04142E2,(E-3)(16)=0.365-0.007E,(E=-3-5) =0.4(E-5) E=lg(s-)/(）(isSi)/(VJ) 式中S为重量计含沙量，i为分组沉速，V为流速. 式(12)或式(13)中比值Ks反映了床面附近泥沙浓度分布规律，这是一个十分复杂尚待探索的重要问题.作者以往根据推移层内沿深度输沙率分布

17、均匀假定，提出Ks=0.81，据此及上列公式计算验证了很多水槽试验资料，结果较满意，计算与实测挟沙力比值KA在0.52.0之间的资料组次占74%84.7%，比国内外一些著名的水流挟沙力公式计算结果更为接近实际9，深入分析分组输沙率计算与实测值比值看出，粗细颗粒除了在床面上及水中互相影响以外，近壁层泥沙浓度分布规律不同，从文献6也可看出，Ks应随增大而变大，并非常数.文献10和11从试验或分析表明，在上举力大于颗粒水下重量的近壁层范围内，含沙浓度分布会出现上大下小.参考文献6提出的跃移层内浓度分布公式，根据大量分组输沙率资料，作者提出了计算Ks的新的初步经验公式： Ks=0(0.04)(17)(

18、0.040.2) (0.2) 式(17)中分界大体上与动床床面形态突变相应14，0.04时qs=0, Ks当然为0；=0.2时大体上沙垄发育最高，相应Ks=0.81，即此前设定Ks为0.81常数大体上是一种平均情况；2时式中Ks1，即近壁层内含沙浓度出现上大下小分布，此时冲积床面一般为高能态逆行沙波.应用公式(17)计算Ks使挟沙力计算精度进一步提高.从定性上看，式(17)是比较合理的，但定量上尚待今后继续研究.    二、非均匀沙挟沙力水槽试验成果 本次试验在长16m宽0.5m的循环可调坡水槽内进行，槽内预先铺沙15cm厚，当放水足够长时间水流输沙达平

19、衡后测量数据.原始床沙分细沙(F)及混合沙(M)两大类，各类原始床沙级配及特征粒径列于表1，表中颗粒分散度=(D84/D50+D50/D16)/2. 试验床沙分散度虽不很大，但与黄河、长江等中下游沙质河流床沙分散度相近.试验共得57组资料，流量Q为0.0210.126m3/s，水深H为0.0940.231m，能坡J为(1.7-12.9)×10-3，含沙量C为13.7229kg/m3.表2、3分别为F及M类沙试验数据.所谓全沙指包括悬移质和推移质的水中全部输沙. 表1 原始床沙级配及特征值试验类别原始床沙级配(mm)，(%)特征粒径(mm)分散度0.020.060.010.020.04

20、D35D50D65DmF9.3744.3392.751001000.0500.0650.0780.0641.873a M1.4423.2043.9788.261000.0800.1070.1250.1121.797b0.8921.2845.9184.421000.0800.1050.1260.1231.920c0.3213.2833.141001000.1020.1210.1440.1512.104表2 非均匀沙挟沙力试验数据(细沙类)组次NO.流量Q/s水深H/m能坡J/1/1000水温T/含沙量C/kg/小于某粒径沙重%(全沙)0.020.060.10 *10.1260.1632.4502

21、4.082.5125.8964.8790.0320.0570.2091.72020.826.5232.6472.2694.2430.1260.2234.74022.396.4023.4357.9690.0540.1030.2193.84023.050.9022.3268.1298.7450.0860.2313.13023.050.6218.0165.4197.1060.0710.1993.87021.554.4827.7870.1297.3670.0530.1892.59022.034.5627.1273.1699.1380.1180.20612.90023.5136.9419.1751.66

22、92.9690.1030.1916.88023.0114.7925.2963.1696.23100.0930.1735.10023.2133.1528.1662.3496.07110.0670.1615.80023.0105.8633.2072.2997.03120.0570.1405.03023.091.8121.3359.3795.80130.0310.1222.16022.839.6730.9376.3997.21140.0260.1161.94022.825.6835.4976.7698.49150.0410.1422.59023.552.2333.2474.9498.31160.08

23、50.1516.25025.0141.3929.2565.3596.77170.1120.1875.60025.5215.7927.9367.8198.06180.0700.1286.47024.8164.3519.6249.5194.73190.0810.1528.10025.0219.8323.5959.0994.84200.0560.1289.00025.0229.3224.8362.4696.84210.0500.1403.02023.084.5327.4668.7096.51220.0430.1651.72023.228.9620.1868.3997.33* D0.2mm沙重占100

24、%，各组床沙级配均取原始级配.表3 非均匀沙挟沙力试验数据(混合沙类)组次NO.流量Q/s水深H/m能坡J/1/1000水温T/含沙量C/kg/小于某粒径(mm)沙重% *0.04 全0.060.10沙0.04床0.060.10沙10.1080.2035.59025.5101.3427.3252.3578.627.9623.2043.9720.1090.2006.45027.076.2838.0761.2992.627.9623.2043.9730.0650.1883.01027.049.3821.3046.9495.307.8622.2046.1340.0910.1994.70026.563

25、.8633.5455.4891.798.2026.9255.3050.0450.1971.72020.513.7542.5260.2691.6810.3332.6464.6560.0650.1703.87020.844.5138.4761.7695.814.7919.0749.0970.0380.1582.15022.022.3340.4161.7495.085.4321.7046.3180.0880.1576.88024.074.4322.5344.6591.527.9826.8851.7890.0310.1293.00024.031.5734.0356.0187.236.5623.6248

26、.50100.0580.1284.30020.753.7727.9248.3993.118.0327.5352.28110.0790.1457.00022.8118.1835.7157.2789.544.1517.5140.66120.0210.0973.01021.934.9236.3759.2787.124.4518.4243.57130.0480.0976.46022.485.0040.0561.9995.634.8518.9542.64140.0400.1013.44022.255.5035.7859.9587.419.9233.8268.13150.0730.1313.01022.4

27、99.3431.6951.7583.515.5722.4351.01160.1090.1926.45023.578.4934.3351.4387.246.0821.2845.91170.0680.1703.66023.045.1332.5753.5186.2211.4635.4667.94180.0960.1895.16021.064.8730.0548.2984.775.3323.0653.49190.0780.2023.44021.537.3628.4048.9481.6011.3433.0663.33200.0560.2032.37020.029.0832.2953.2687.526.1

28、821.7547.41210.0710.1644.09020.549.1334.7454.0186.6413.2336.8865.08220.0390.1702.15018.525.1647.8969.3197.875.7922.2749.52230.0850.1675.38019.778.2929.7353.0384.1010.4031.5861.09240.0440.1253.44018.147.9831.2452.4988.205.5022.8651.33250.0550.1293.87019.052.4535.4858.4989.5812.6838.9671.12260.0730.14

29、05.80018.079.3524.7545.1384.928.9827.1054.00270.0300.1113.01019.041.5437.7456.5189.544.9819.9744.80280.0320.0943.44018.551.8733.4552.8089.8410.6332.0060.52290.0410.0943.87019.063.6138.5661.2586.677.7425.7552.65300.0630.1323.66019.749.7032.5655.2786.923.6014.1830.71310.0600.1055.16021.074.0526.8247.0

30、185.134.8418.6941.92320.0710.1712.79019.534.9030.5050.4882.425.4918.7736.30330.0500.1592.59019.527.7238.3663.2591.473.8414.5131.95340.0680.1992.58019.423.8131.3353.5491.867.3526.6556.39350.0480.1083.87020.043.8231.7756.6190.784.2416.5535.05*  全沙与床沙中D0.4mm沙重均占总重的100% 由于沙波高度及床沙交换层厚度随机变化，床沙级配

31、随时空变化，因此测取与每一组试验水流输沙相对应的代表性床沙级配i/-0是个难点.试验中对每一组输沙平衡后的床沙都进行了测量，取样厚度约为平均沙波高度.经过分析发现，对于粉细沙(F类)，因床沙大多全面起动，床沙级配变化较小，可取原始床沙级配作为各测次代表床沙.M类每测次床沙实测数据则列于表3.分析此次及以往所有试验资料发现，床面形态剪切力不仅对推移质输沙基本无效，而且对以悬移质为主的输沙率影响也不明显.    三、泥沙非均匀性对水流分组挟沙力的影响床沙和水中泥沙粒径的非均匀性会显著影响各分组挟沙力.在床面上，粗颗粒对细颗粒会产生隐蔽作用，使细颗粒受到的实际

32、水流上举力及拖曳力小于均匀沙情况，反之，粗颗粒相对于细颗粒比较突出暴露，其受到的水流作用力将大于均匀沙情况.当泥沙起动、起跳以后，水中推移层内不同粒径的颗粒因运动速度、轨迹不同而会互相碰撞，悬移泥沙也会因沉速不一而碰撞，在床面附近因含沙浓度大而碰撞愈甚.因细泥沙比粗泥沙跳跃较快，沉降较慢，故互相碰撞、动量传递结果总体上使细颗粒跳跃变慢，沉降加快，而粗颗粒则反之.据上分析看出，水中泥沙与床面泥沙的非均匀性对分组挟沙力的影响总体上是一致的，均使同样水流强度下水流对细颗粒的挟沙能力比均匀沙情况减小而对粗颗粒则增大.由于此种改变作用迄今都只能根据水流分粒径组挟沙力资料反算，难以将起动阶段的隐蔽暴露影响

33、及水中运动阶段的互相碰撞影响加以分开，因此我们在式(14)中统一加上泥沙非均匀性影响系数KD，则式(14)变为式(18)：iTqT=KD(ibqb+isqs)=KDi0Dq(1+11.6Ks(P1I1+I2)u*/ub). (18)图2 KD和DK关系在分析KD变化规律时，将实测非均匀沙各粒径组输沙率iTqT与据式(14)不考虑非均匀沙内相互影响计算求得的iTqT之比作为KD，将本次试验及文献2的资料求得的KD和DK关系点绘于图2.纵坐标DK按下式计算DK=(Di)/(Dp) , (19) 式中Di为混合沙中各粒径组平均粒径，Dp为粗细颗粒分界粒径当量，Dm则为床沙平均粒径.DK可谓相对粗细度

34、，DK1之粒径组为相对粗颗粒，DK1之粒径组为相对细颗粒.由图可见， KD与DK形成变化趋势很明显的带群，以残差平方和最小为标准求得拟合曲线绘于该图中，即得KD与DK的关系为                  (20)图2点群和式(19)、(20)表明，以DK=1即Di=Dp为分界，Dk1时KD1，而Dk1时KD1.分界粒径当量Dp不仅取决于Dm，而且随水流强度增大.随着水流增强，床面粗粒背水面旋涡会增大，在相同Dm时隐蔽的细粒愈多或对细粒的隐蔽作用愈强，同时水中颗粒运动速度增大、含沙浓度增大而碰撞机会增加，这二者都会导致细粒挟沙力减小，即KD比1减