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

DOC格式论文，方便您的复制修改删减明槽水流的非均匀沙挟沙力研究（作者:_单位: _邮编: _） 摘要：本文从床沙、推移质及悬移质互相交换、衔接的物理图形出发，在力学、随机及紊动扩散分析基础上，提出了颗粒跃移及推移质、悬移质和全沙统一的新的非均匀沙挟沙力公式.文中公布了新近试验取得的57组非均匀沙挟沙力水槽资料，分析提出了反映混合沙中粗细颗粒互相影响的系数的关系，提出了悬移层与推移层交界面高度及浓度预报新的初步关系.应用本文提出的挟沙力系列公式计算验证了大量天然河流及水槽实测分组和总挟沙力资料，计算与实际符合较好. 关键词：非均匀沙 颗粒跃移 悬移质 挟沙力 隐蔽暴露影响 水流挟沙力是泥沙研究中的核心问题，是解决水库、河渠等大量泥沙冲淤问题的关键.天然冲积河流中的泥沙一般都是非均匀沙，悬移质、跃移质及床沙通常都是互相交换、密切相关的，但迄今国内外大多数挟沙力公式从物理图形、原始推导上就与床沙组成无关，实际上是针对且只适用于均匀沙.一些公式经过各种经验性处理以后已用于非均匀沙计算，在生产上得到了广泛应用，但一般都未考虑混合沙中粗细颗粒间的相互影响，分组挟沙力计算精度难以提高，在河床连续冲刷计算中问题暴露往往比较明显.迄今，多数挟沙力公式将悬移质和推移质分割，生产上应用的多是并不匹配的二套公式，对于悬移质和推移质输沙并重的河流，如此作法存在不少问题.水流的非均匀沙挟沙力仍是一个远未解决的重大生产和学术问题. H.A.Einstein1-2早在五十年代即提出了水流的非均匀沙挟沙力理论及公式，半个世纪以来的研究和实践表明，爱氏提出的总的物理图形比较合理.尤其是将床沙、跃移质和悬移质三者联系起来分析而不是将它们看成彼此孤立、分割的事物，这一点是意义深远的.这一物理图形中一个关键的环节是床面与悬移层之间过渡薄层(跃移层或准跃移层)内的水流泥沙运动规律，由于此近底层内规律非常复杂且试验困难，致使爱氏体系的理论发展比较滞缓，爱氏当初推导挟沙力关系时不得不作了很多假定，不少是有问题的.爱氏对颗粒跃移并未进行详细的力学分析和试验，诸多关键的跃移参数相对值均以常数代之，与实际相差甚远，使其挟沙力公式计算成果与实际往往相差很大. R.A.Bagnold, M.S.Yalin及Van Rijn3等很多学者都先后进行了颗粒跃移的力学分析或试验，取得了很多宝贵成果.自1986年以来，王士强、张仁、付晓及胡春宏、惠遇甲等分别对明槽水中颗粒跃移进行了力学分析和试验，取得了很大进展4-6，分别发现颗粒容重对跃移参数的影响规律，与以往所有认识均不相同.在颗粒跃移分析和试验基础上，王士强等提出了计算推移质输沙率的新公式4，进而提出了推移质、悬移质和全沙统一的新的非均匀沙挟沙力公式7，并得到了清华大学活动大钢槽内试验取得的27组实测资料的良好验证8，文献9对此公式进行了进一步的分析和验证，并与国内外其它一些著名的挟沙力公式进行了验证计算精度比较. 近些年来，王光谦、倪晋仁10，邵学军、夏震寰11及胡春宏、惠遇甲6等用不同方法探索研究了近底层内颗粒浓度分布规律，为挟沙力研究的进一步深入提供了宝贵的启示；王士强等考虑近壁层颗粒沉降阻力增大、沉速减小从而得出新的悬移层浓度垂线分布公式12，并进一步总结、验证和应用了新公式.近二年内陈骥完成了新的57组非均匀沙挟沙力试验，在所有以上研究基础上，本文对更多的资料进行了分析，继续研究改进了王士强等新的挟沙力系列公式，进一步提高了该套公式的计算精度. 一、颗粒跃移及推移质、悬移质、全沙挟沙力方程 1.1颗粒跃移方程床面附近的颗粒跃移是悬移质输沙和沙质床面之间过渡层内颗粒的主要运动形式，是研究沙质推移质和悬移质输沙率的基础.作用在床面附近跃移颗粒上的力很多很复杂，文献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*为水流摩阻流速，s、分别为颗粒和水的容量，J为能坡，R为沙粒水力半径，按文献14给出公式求得，平整床面时即为水力半径或水深. (3) (4) (5) 式中 上列颗粒跃移参数公式得到了验证5，与其它著名公式相比9，数值相近但有二点重要区别，上列公式中的相对跃移参数不仅取决于，还另外随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 (6)式中P为FL/W1条件下床面颗粒的起跃概率，A*为待定系数.从物理模式及原始推导看，严格来说，式(6)本只适于跃移推移质，但因以滚动、滑动等形式运动的接触质受到的主要作用力与上述跃移颗粒大体相同，因此控制qb的变量及s/二者应该相同.当式(6)中的待定系数由包括接触质在内的变化范围很大的大量推移质输沙率实测资料适线确定时，该式的适用范围就具有普遍性.式(6)中的变量可由一套解析公式分别求出，相当烦复，为便于应用，由式(6)及其它一系列求其变量的公式概括简化得出公式(7)，两者基本等价，计算值十分相近. (7) 式中：. 图1列出了大量实测点据及公式(7)中不同s/的三条曲线，公式与各种s/的实际资料均相当符合，在高强度输沙部分，与Wilson的二种容重颗粒资料均分别符合相当好.图1还列出了其它著名公式曲线以便互相比较，对于泥沙，公式(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，该式与通常的公式形式相同，仅以修正悬浮指数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 . 设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)中的所有粒径D均为混合沙内各粒径组的平均粒径. 上列公式中的a值可由下述方法求得.因颗粒平均跃移高度ao的计算公式(3)中的系数是根据水槽颗粒跃高试验资料适线确定的，因此在水槽输沙计算中取a=a0,这是合理、自然的，如此计算出的悬移层底部浓度及挟沙力与大量水槽内的实测资料符合相当好. 对于天然河流、水库，无任何实测跃移层高度a的资料，若仍取a=a0,则计算出的悬移层底部浓度一般都远大于实测资料.因天然河流、水库内的水深远大于水槽情况，作为水深方向的物理量a，可以推测和假定随相对水深H/D50而增大，根据一些天然河流、水库实测底部含沙量资料分析反算，得出一般含沙量情况下天然河流内的a可按经验公式(15a)计算，a可称为准跃移层厚度： (15a)对于长江及黄河大部分非高含沙水流情况，按式(15a)计算的悬移质挟沙力与大量实际资料符合较好.但对于黄河少数高浓度情况，若按式(15a)计算又计及沉速修正后得出的悬移层底部含沙量仍明显小于实际情况. 从上述公式(3)可以看出，跃移层厚度随水流容重增大而减小，在高含沙水流情况下，随浓度增加而明显增大.另外，随着浓度增加，对数流速公式中的卡门常数将会减小，即由床面向水面方向的流速变大的梯度增大，床面附近水流垂向紊动强度向水面方向变大的梯度也理应增大，这就意味着相同平均水流强度时含沙浓度愈高，颗粒由跃移转为悬移的高度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反映了床面附近泥沙浓度分布规律，这是一个十分复杂尚待探索的重要问题.作者以往根据推移层内沿深度输沙率分布均匀假定，提出Ks=0.81，据此及上列公式计算验证了很多水槽试验资料，结果较满意，计算与实测挟沙力比值KA在0.52.0之间的资料组次占74%84.7%，比国内外一些著名的水流挟沙力公式计算结果更为接近实际9，深入分析分组输沙率计算与实测值比值看出，粗细颗粒除了在床面上及水中互相影响以外，近壁层泥沙浓度分布规律不同，从文献6也可看出，Ks应随增大而变大，并非常数.文献10和11从试验或分析表明，在上举力大于颗粒水下重量的近壁层范围内，含沙浓度分布会出现上大下小.参考文献6提出的跃移层内浓度分布公式，根据大量分组输沙率资料，作者提出了计算Ks的新的初步经验公式： Ks=0(0.04)(17) (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厚，当放水足够长时间水流输沙达平衡后测量数据.原始床沙分细沙(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.04D35D50D65DmF9.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.45024.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.6692.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.0850.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%，各组床沙级配均取原始级配.表3 非均匀沙挟沙力试验数据(混合沙类)组次NO.流量 Q/s水深 H/m能坡 J/1/1000水温 T/含沙量 C/kg/小于某粒径(mm)沙重% *0.04 全0.06 0.10 沙0.04 床0.06 0.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.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.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.499.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.1821.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.1405.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.0185.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% 由于沙波高度及床沙交换层厚度随机变化，床沙级配随时空变化，因此测取与每一组试验水流输沙相对应的代表性床沙级配i-0是个难点.试验中对每一组输沙平衡后的床沙都进行了测量，取样厚度约为平均沙波高度.经过分析发现，对于粉细沙(F类)，因床沙大多全面起动，床沙级配变化较小，可取原始床沙级配作为各测次代表床沙.M类每测次床沙实测数据则列于表3.分析此次及以往所有试验资料发现，床面形态剪切力不仅对推移质输沙基本无效，而且对以悬移质为主的输沙率影响也不明显. 三、泥沙非均匀性对水流分组挟沙力的影响 床沙和水中泥沙粒径的非均匀性会显著影响各分组挟沙力.在床面上，粗颗粒对细颗粒会产生隐蔽作用，使细颗粒受到的实际水流上举力及拖曳力小于均匀沙情况，反之，粗颗粒相对于细颗粒比较突出暴露，其受到的水流作用力将大于均匀沙情况.当泥沙起动、起跳以后，水中推移层内不同粒径的颗粒因运动速度、轨迹不同而会互相碰撞，悬移泥沙也会因沉速不一而碰撞，在床面附近因含沙浓度大而碰撞愈甚.因细泥沙比粗泥沙跳跃较快，沉降较慢，故互相碰撞、动量传递结果总体上使细颗粒跳跃变慢，沉降加快，而粗颗粒则反之.据上分析看出，水中泥沙与床面泥沙的非均匀性对分组挟沙力的影响总体上是一致的，均使同样水流强度下水流对细颗粒的挟沙能力比均匀沙情况减小而对粗颗粒则增大.由于此种改变作用迄今都只能根据水流分粒径组挟沙力资料反算，难以将起动阶段的隐蔽暴露影响及水中运动阶段的互相碰撞影响加以分开，因此我们在式(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可谓相对粗细度，DK1之粒径组为相对粗颗粒，DK1之粒径组为相对细颗粒.由图可见， KD与DK形成变化趋势很明显的带群，以残差平方和最小为标准求得拟合曲线绘于该图中，即得KD与DK的关系为 (20)图2点群和式(19)、(20)表明，以DK=1即Di=Dp为分界，Dk1时KD1，而Dk1时KD1.分界粒径当量Dp不仅取决于Dm，而且随水流强度增大.随着水流增强，床面粗粒背水面旋涡会增大，在相同Dm时隐蔽的细粒愈多或对细粒的隐蔽作用愈强，同时水中颗粒运动速度增大、含沙浓度增大而碰撞机会增加，这二者都会导致细 粒挟沙力减小，即KD比1减小愈多.但当水流相当强以后，水流继续增强不会使床面粗粒背水面旋涡进一步增大，水中颗粒跃移速度及轨迹变化也愈来愈小.这些反映在式(20)中DK0.4以后DK继续减小(即Di、Dm一定时继续增强)但KD却变化甚微，约为0.4. 根据下面计算验证看出，上列KDDK关系同样适用于天然河流. 分析钱宁和爱氏资料2发现，凡是床面充分扰动情况，床面上粗细颗粒间的隐蔽暴露影响看来基本不再存在，任何DK时KD均取1合适. 应用上列公式计算得出的各试验组次总的及分组挟沙力与实测值对比列于表4，表中KA为挟沙力计算值与实测值之比，PK为KA=0.5-2.0之间的组次占计算总组次的百分比.由表可见，对于总输沙率，PK均超过90%，平均达94%，对分粒径组，PK值平均也超过82%，高者达90%. 表4 输沙率计算与实测对比资料来源类 别资料总数平均PK(%)陈骥细沙总输沙率分组输沙率22 611.19 /95.5 90.2陈骥混合沙总输沙率 分组输沙率35 1050.94 /91.4 79.1爱氏和钱宁总输沙率 分组输沙率26 1361.20 /96.1 81.6总 计总输沙率 分组输沙率83 302/ /94.0 82.5四、计算挟沙力公式验证 应用王士强等公式(3)(20)计算验证王士强等于1989年在清华大学60m长1.2m宽活动循环大钢槽内完成的27组非均匀沙挟沙力试验资料8，验证结果列于表5，并将分粒径组验证结果绘于图3，该资料在上列公式、关系建立中均未利用过. 图3 分组输沙率计算验证表5 王士强等挟沙力公式计算验证(水槽)类 别资料总数平均PK(%)总输沙率271.1692.6分粒径组781.1682.1由表可见，总输沙率及分组输沙率KA=0.52.0之间的组次比例PK分别高达92.6及82.1%，计算精度很高.对于水槽总输沙率资料，国内外一些著名公式计算结果PK值对几十组资料而言超过80%就不多见了.至于非均匀沙分粒径组的挟沙力计算验证，在国内外还很少见. 表6为应用作者公式计算验证大量黄河下游及长江资料结果，由表可见，PK对总输沙率平均约达75%，有些河段超过80%，对分组输沙率均超过60%，对于影响因素十分众多复杂的黄河来说，这一验证结果是比较好的. 表6 王士强挟沙力公式计算验证(天然河流)资料来源类 别资料总数平均PK(%) KA=0.52.0长江荆江总输沙率 分粒径组55 1241.16 /80.0 62.1黄河 利津总输沙率 分粒径组240 3930.94 /71.3 63.1黄河 艾山总输沙率 分粒径组289 5960.86 /77.9 68.5黄河 孙口总输沙率 分粒径组221 4291.06 /79.2 66.2黄河 夹河滩总输沙率 分粒径组187 3191.16 /84.0 66.5天然河流内某组输沙率资料是否为平衡挟沙力主要只能从总量上大概判断，而输沙总量即使平衡时粗颗粒尚可能处于淤积过饱和状态而细颗粒则可能反之，实难判别和选取各粒径组都处于平衡输沙的资料，这也是验证天然河流挟沙力特别是分组挟沙力的难点. 作者曾应用此套公式验证过美国的几条河流，结果也比较满意7.将此套公式与国内外其它5个著名公式同时计算过639组不同容重颗粒的水槽挟沙力资料，结果此套公式与实测资料更为接近，特别是计算物理模型中常用的轻质沙挟沙力的精度更比其它公式要高9.本公式应用于计算黄河一些水库及河道河床冲淤，也取得了比较满意的结果. 五、结语 1.本文提出了同时计算推移质、悬移质及全沙非均匀沙挟沙力公式.该套公式的基本物理图形是床沙、跃移质和悬移质互相交换、衔接. 2.无量纲颗粒跃移参数及推移质输沙率参数不仅随增大，还发现另外随颗粒容重减小，这不仅是理论推导结果，也为试验资料所证实，这就使据此基础得出的全沙挟沙力公式不仅适合于河流泥沙，而且计算轻质沙挟沙力的精度比迄今国内外一些著名公式明显要高.3.沙床近底层含沙浓度分布是决定水流挟沙力的重要因素.根据新近的一些研究成果及实际挟沙力资料分析，本文提出了悬移层底层高度a及浓度Sa的初步新关系.跃移层内浓度分布随增大而愈趋均匀且会出现上大下小现象. 4.非均匀混合沙内的粗细颗粒，在床面上存在相对隐蔽和暴露影响，在水中会互相碰撞，本文提出了非均匀性影响系数KD的新关系. 5.王士强等非均匀沙挟沙力公式已为大量水槽及天然河流资料所验证，计算结果比较符合实际. 6.随着对近底层水沙运动规律研究的发展，作者将继续研究改进这套挟沙力公式. 参考文献1.Einstein, H.A. The Bedload Function for Sediment Transportation in Open Channel Flows. Technical Bulletin. No.1026. U.S. Dept. of Agriculture, Soil Conservation Service, 1950.71.2.Einstein H A, Ning Chien. Transport of Sediment Mixtures with Large Ranges of Grain Sizes. M.R.D., U.S.Corps. Engrs., 1953. 49.3.Van Rijn. L C Sediment Transport, Part : Bed Load Transport, Proc. ASCE, J. Hydr.Engr. 1984，110(10).4.Wang Shiqiang, Zhang Ren. A New Equation of Bed|load Transport. Proc. of 22 Congress IAHR. Lausanne, 1987.5.王士强，张 仁. 颗粒比重对水流挟沙力的影响. 泥沙研究. 1990.6.胡春宏，惠遇甲. 渠挟沙水流运动的力学和统计规律. 北京：科学出版社，1995.7.Wang Shiqiang, Zhang Ren. Sediment Transport Rate for Nonuniform Sand. Proc. Intern. Symp. on Sedim. Trans. Modeling. New Orleans, U.S.A., 1989.8.Wang Shiqiang, Zhang Ren. Experimental St