复式断面河槽流速横向分布及其对滩唇形成的影响.docx

1、水利学报2009年5月SHUILIXUEBAO第40卷第5期Z，文章编号：0559-9350（2009）05-0535-07复式断面河槽流速横向分布及其对滩唇形成的影响吴腾，张红武，张昊（清华大学水沙科学与水利水电工程国家重点实酸室，北京100084）摘要：通过对水流微小控制体进行受力分析，建立了复式断面流速横向分布模型，该模型考虑了侧向副流惯性力的影响c采用实验资料对该模型进行检验，模型计算值与实测值符合较好。运用该模型计算了不同条件下复式断面流速和挟沙力的横向分布，并定量分析了横向分布特性对滩唇形成的影响。分析表明，水流漫滩时.复式断面的横向挟沙力变幅较大，尤其姑在主槽和河漫滩交界处，挟沙

2、力减小迅速，而含沙最减小相对较小，泥沙发生淤积，容易形成滩唇。单从水流挟沙力角度考虑，水流漫滩后水深越小、滩槽的糙率差越大，越容易形成滩唇。关键词：复式断面；流速横向分布；滩唇；河漫滩；侧向副流中图分类号:TV143文献标识码：A滩唇是指河道主槽两侧较高的滩岸边缘部分，它通常高于滩面。滩唇大量存在于具有复式断面的河道中，如黄河下游、淮河干流等，它的形成与河道断面形态和来水来沙条件密切相关。黄河下游河道断面属于典型的复式断面，具有宽阔的河漫滩。滩唇的存在直接影响平滩流量乃至造床流量的确定。复式断面不同于一般的单一断面，当来流量较大时，水流溢出主槽形成漫滩水流，此时过水断面突然展宽，各种水力要素发

3、生突变，其规律与单一断面相异，尤其是在滩和槽的交界面附近，水流结构异常复杂,远非单一断面的水沙运动规律可以描述。复式断面水沙运动特性是河流动力学研究的传统课题之一，目前已取得较多成果，包括复式断面深度平均的水流运动方程，复式断面的动量修正系数特性成，复式断面的水流泥沙特性，复式断面的阻力特性，复式断面流速分布的研究时，复式断面滩槽分流比的研究SE以及通过河工模型试验和数学模型计算对高含沙洪水造床过程的研究等等,，2-，3）o本文在前人研究的基础上，对动量方程进行了简化，同时考虑了侧向副流惯性力，建立了复式断面流速横向分布模型，并基于该模型定量分析断面流速和挟沙力的变化规律。由于流速和挟沙力对断

4、面形态的形成有直接的影响，本文根据断面流速和挟沙力的变化规律进一步探讨滩唇的形成机理，初步解释了黄河下游河道断面普遍存在滩唇的原因，为黄河下游河道的治理提供了参考。1断面流速横向分布模型取水流微小控制体,控制体中心点坐标为（x.y.z）,如图1所示。其中，敞、海、如分别为控制体微段的长度，对*方向进行受力分析，设作用在控制体六个面上的应力分别为九、/“g、f,、代、/;2,其表达式分别为：f,=P*(1)8x机T8xT8?收稿日期:2008-01-03基金项目：国家自然科学基金委员会、水利部黄委会黄河联合研究蓦金项目（50339020）作者简介：矣胸（1972-），男，湖北人，博士生.主要从事

5、水力学及河流动力学的研究。E-mail：wutengQ5mails,23yfy2=2dy8z9r+TdT&a七7d7图1投制体示意图1投制体示意式中，p,为方向上的压应力，j和t“分别为其余两方向的雷诺剪切应力。作用于控制体的表面力等于各面上的应力乘上对应的作用面积，即F,=S-九2）+（%-人2）跄z+S-后）6y6x式中，F,为控制体的综合表面力。将式（1）式（6）代入式（7）可得：匚=（-登+寿+譬设X为单位体受的质傲力在方向的分量，根据牛顿运动定律F=即沿z方向的合力等于该方向的质量与加速度的乘积，有X=pg8x6ydzsn0式中/水流密度，g重力加速度，。为河床底坡角度。将式（8）和

6、式（10）代入式（9）,可得：草d3rdx+dy+dzX=pg8x6ydzsn0式中/水流密度，g重力加速度，。为河床底坡角度。将式（8）和式（10）代入式（9）,可得：草d3rdx+dy+dzpg8x3y$zsin0+=p(10)(11)对于质量力在方向的分量有:对于流速在x方向的全导数有如下关系式：du,du,du,du,8u,(12)37=+蜘在+与有+虬奇式中，x、y、z表示水流方向、横断面方向和垂向，虬、与、虬分别为x、y、z方向的流速。对于恒定流而言箫=0,则式（12）可化为：(13)du,du,dutdu,2T=虬云+%和+匕W将式（13）代入式（11）,有:dydzdu,duM

7、du,由于本文主要探讨流速的横向分布情况，认为不受方向的影响，故可略去/方向的分量，可得到恒定水流横向控制方程：pg8x8y8z%n0+(14)(15)(16).Odrn3r=+奇aT将式（15）沿垂线积分，当断面横向变形不剧烈时，可以得到垂线平均水流运动方程：.8phUMPghsmO+-“=习式中0为断面横向节点的水深，门为边壁剪切应力，取经验公式为阻力系数，乙为垂线平均流速沿方向的分缺，七为垂线平均流速沿y方向的分量。雷诺剪切应力在实际应用中取：g式中，J为节点横向涡黏系数，表达式可写为:“=Mu.（18）式中4为断面节点处水深.为节点处摩阻流速，为无量纲系数，取值可参考相关文献，摩阻流速

8、又可表示为:(19)将式（19）代入式（18）得到:(20)将式（17）、式（20）代入式（16）可得：(21)皿sin。,彖（p叫#）5仇翌）-如0=迦曾打式中，啰为侧向副流惯性力，与水体的重力和边壁的剪切力相关，所以可以假设:(22)斗（pgAsin。-死，k为无量纲系数，在应用中可通过实测资料反求。最后将式（22）代入式（21）可得漫滩水流横向控制方程：(23)（1-A）（溯sin6-gpU，+彖（p（f）翌卜本文采用差分法离散式（23）,直接进行求解。2模型的验证2418E12u602418E12u600102030405060708090100图2实验水槽示意2.1流速的验证选用中国

9、水利水电科学研究院的实验资料进行验证分析，断面形状如图2。实验的基本数据为：主槽宽0.3m,滩地宽0.7m,滩槽高差0.06m,坡度为0.001,水槽主槽糙率为0.011.滩地糙率0.013o在计算中，取水的运动黏度取为1.139x10_6m2/so图3为不同水位下的断面流速验证图，由于试验水槽成中心对称，对称轴左右两侧的流速分布相同，文中仅给出了半槽宽的流速分布图。由图可知，通过数值计算直接求解式（23）得到的结果与实测值变化趋势相同，且数值较为接近。图中主槽的流速明显大于滩地的流速，并且在滩槽交界处，流速变化较快；同时也可以看出，小水深时滩槽流速差异较大，随着水深的增加，这种差异逐渐减小。

10、理论而言，河漫滩糙率较主槽大，尤其是水深不大时糙率对流速的影响尤为明显，当水深增大时，糙率对流速的影响减弱，滩槽流速的差异减小，数值计算结果与理论分析一致。因此，直接求解式（23）的流速横向分布基本可以反映实际情况，可用于复式断面的研究。2.2挟沙力的验证本文第一作者曾在曹志先教授指导下采用文献16中的实验资料对挟沙力变化进行验证。实验采用一滩一槽的复式断面：主槽宽0.6m,滩宽1.4m,滩槽高差0.1m,床面纵向坡降0.1%,试验沙为天然沙，中值粒径0.105mm。图4为采用流速横向分布模型计算的不同水位时断面的挟沙力，其中，水流挟沙力公式可采用Gu。或张红武给出的公式计算。例如，水沙河流条

11、件下，张红武公式简化为S=。皿3不In（酒）（24）(b)水深8.25cm流速巍证图距跚cm605040302010(S.E0帝*2040(_*/实测值计算值试验槽!距离/cm4020(c)水深9.08cm流速救证图(e)水深10.75cm速验证图(d)水深10.03cm流速验证图(0水深12.33cm流速验证图图4挟沙力我证式中，S.为水流挟沙力，3为泥沙沉速，0#为床沙中径，本文的水流挟沙力均采用该公式进行计算。可以看出计算值与实测值的变化规律相同，在漫滩时断面挟沙力迅速减小，然后随水深逐渐增大，式(23)计算的流速可以用于水流挟沙力的计算。3滩唇的成因分析3.1单一断面与复式断面的比较分

12、析图5和图6分别为水深0.05m时矩形断面的流速和挟沙力横向分布图，断面宽0.8m,糙率0.011,比降0.001：图7和图8为水深0.2m时复式断面的流速和挟沙力横向分布图，对应的断面形态与图2相似，主槽宽0.8m,滩槽总宽度为4m,主槽糙率为0.011,滩地糙率为0.025,河槽比降O.OOlo对于单一河槽，全断面挟沙力大小较为均匀，仅在河槽的两岸挟沙力迅速减小。对于复式断面，断面横向水流挟沙力变化幅度较大，尤其是在主槽和河漫滩交界处，挟沙力减小迅速。对于复式断面，当水流漫滩时，水流由主槽流向滩地，含沙最变化不大，而滩地挟沙力降低较快，导致泥沙快速淤积，水体的含沙量也快速减小，沿横向的淤积

13、量将逐渐减小，故而在滩槽交界附近淤积量最大,远离主槽处淤积量逐渐减小，久而久之，在主槽边缘的滩地上逐渐形成高出附近滩面的自然沙埴，成为(-sEbw图6单一河槽挟沙力横向分布图5单一河槽流速横向分布图7复式河槽流速横向分布图8夏式河情挟沙力横向分布3.2不同水深对滩唇形成的影响图9和图10为不同水深条件下流速和挟沙力沿断面横向分布图,计算的断面形态如图2,具体参数为，半槽宽0.4m,半滩宽2m,滩槽高差0.06m,河槽的糙率为0.011,河漫滩糙率为0.025,计算的水深分别为0.1m、0.2m、0.3m、0.4m。由图可知，不同水深条件下，断面的横向流速分布较为类似，主槽流速大，河漫滩流速小，

14、滩和槽的交界附近流速迅速减小。挟沙力沿断面的横向变化规律与流速变化规律相似，不过在滩槽交界面附近挟沙力变化较为复杂，存在一定的波动，主要是因为滩槽交界附近虽然流速快速减小，但是对应的水深也减小较多，由水流挟沙力的计算式可知，水流挟沙力与流速的三次方成正比，与水深成反比，挟沙力的大小处决于何者占优，故该处的挟沙力会出现一定的波动。当水深较小时，滩和槽的流速差很大，挟沙力也相差较大；当水深增大时，滩槽间的流速差逐渐减小，水流挟沙力的差值也逐渐减小。水流漫滩过程中，水流由主槽流向滩地，滩槽交界附近的水体含沙量与主槽的含沙fit相差不大，水深较小时，滩上的挟沙力远小于主槽，说明滩地水体的挟沙能力严重不

15、足，泥沙淤积程度大，当水深增大时，滩和槽的挟沙力差逐渐减小，滩地水体的挟沙能力逐渐与含沙量相适应，泥沙的淤积逐渐减小，由此可以看出漫滩后水深越小越容易发生淤积，越容易形成滩唇。3.3不同糙率对滩唇形成的影响图11和图12分别为不同糙率条件下流速和挟沙力沿断面横向的分布图，该断面的具体参数为，半槽宽0.4m,半滩宽2m,滩槽高差0.06m,河槽的水深0.Im,河槽的糙率为0.011,滩地的糙率分别选取0.011、0.015、0.02、0.03。由图可知糙率对滩地流速影响明显，滩地的糙率越小，流速越大，挟沙力也越大，主槽和滩地上的挟沙力差别越小;反之，滩地的糙率越大，流速越小，挟沙力也越小，主槽和

16、滩地挟沙力的差异也越大。无论是糙率如何变化，从主槽到河漫滩的过渡过程中，皆出现挟沙力迅速减小的区域，且在主槽和滩地交界面附近。水流刚漫滩时，主槽的水流进入滩地，水体的含沙量与主槽的含沙量相差不大，此时滩地上水体的含沙量处于超饱和状态，水流为达到平衡进行自动调整，泥沙在滩地发生淤积，尤其是水流挟沙力梯度最大的地方淤积最为明显，易于形成滩唇。图9不同水深时流速9向同I.20.80.47s提图10不同水深时挟沙力横向分布图11不同糙率时流速横向分布2.533.5半槽鬣mn=0.01ln=0.0l5n=0.02n=0.03半槽竟m图13不同滩宽时流速横向分布图14不同滩宽时挟沙力横向分布不同糙率时挟沙

17、力横向分布3.4不同滩槽宽度对滩唇形成的影响图13和图14为不同宽度的河漫滩条件下流速和挟沙力横向分布图，具体参数为，半槽宽0.4m,滩槽高差0.06m,河槽的水深0.1m,河槽的糙率为0.011,河漫滩糙率为0.025,河漫滩宽度分别为lm、2m、3m、4m。由图可知，当河漫滩宽度发生变化时，主槽和滩地的流速大小基本不发生变化，对应的挟沙力也变化不大。在主槽过流一定而滩地过流随河漫滩宽度增加而增加的条件下，河漫滩的宽度变化对滩唇形成的影响不明显。4结论对复式断面的水流微小控制体进行受力分析，建立了断面流速横向分布模型，并采用实验资料对该模型进行检验，模型的计算结果与实验资料符合较好,该模型可

18、以用于滩唇成因的分析。采用该模型分别探讨不同水深、不同糙率和不同滩地宽度时复式断面的流速分布和挟沙力分布规律。结果表明，复式断面主槽流速大，河漫滩流速小，滩和槽的交界附近流速迅速减小，挟沙力沿断面的横向变化规律与流速变化规律相近;在水流漫滩后，水深越小滩槽的挟沙力差值越大，滩槽的糙率差越大，滩槽的挟沙力差值也会越大。水流漫滩过程中，水流由主槽流向滩地,滩槽交界附近的水体含沙量与主槽的含沙量相差540不大，若滩和槽的挟沙力相差较大，则泥沙淤积最增加，越容易形成滩唇，即水流漫滩后水深越小、滩槽糙率差越大，滩槽间的挟沙力的差异越大，越容易形成滩唇。参考文献：(1张红武，张清，江恩惠.黄河下游河道造床

19、流量的计算方法J.泥沙研究，1994(4):30-38.2ShionoK,KnightDW.TurbulentopenchannelflowswithvariabledepthacrossthechannelJ.Fluid.Meeh.,1991,222:617-646.3YangK,CaoS,LiuX.TransversekineticenergycorrectioncoefficientinhydrauliccomputationofcompoundchannelsCProceedingsoftheInternationalSymposiumonRiverSedimentation,Octo

20、ber,Yichang,China,2004：1247-1252.4王双明，孙西欢，潘光在复式断面河道泥沙数值模拟J.太原理工大学学报.2000,31(1):60-63.5杨克君，曹叔尤，刘兴年.复式河槽综合糙率计算方法比较与分析J.水利学报,2005,36(7):780-786.6李彪.胡旭跃.徐立军.复式断面滩槽流速分布研究综述J.水道港口，2005,26(4):228-232.7童汉毅，槐文信，张礼卫.复式断面明渠的流速比和流量比J.武汉大学学报(工学版),2003,36(5):9-13.8许唯临.复式河道漫滩水流计算方法研究口.水利学报,2002(6):21-26.(9KnightD,

21、DemetriouJ.FloodplainandmainchannelflowinteractionJ.JournalofHydraulicEngineering,(ASCE),1983,109(8):1073-1092.10 王树东.漫滩水流二维流速分布及水力学计算J.水利学报,1986(11):51-58.11 刘星，贺晓春.天然河道断面流速计算J.水运工程,1998(11):29-31.12 张红武，江恩惠，等.黄河高含沙洪水模型的相似律M.郑州：河南科学技术出版社,1994.13 张红武，江恩惠，吕昕.黄河高含沙洪水异常现象成因分析J.人民黄河,1993(9):8-11.14 周宜林，

22、王明甫.滩槽挟沙水流流速横向分布研究J.泥沙研究,1996,9(3):56-63.15 胡春宏，吉祖稳，等.复式断面非均匀沙滩槽交换规律研究R.北京：水利水电科学研究院，1997.16 丁君松，王树东.漫滩水流的水流结构及悬沙运动J.泥沙研究，1989(1):82-87.17 GuoJ.LogarithmicmatchinganditsapplicationincomputationalhydraulicsandsedimenttransportJ.J.Hydr.Res.,IAHR,2002,40(5)：555-566.TransversedistributionofflowvelocityincompoundchannelanditseffectsonfloodplainsleveesWUTeng,ZHANGHong-wu,ZHANGHao(TtinghuaUniversity,Beijing100084,China)Abstract:Basedontheanalysisofforcesactingoncompoundchannelflows,amodeloftransverseflowvelocitywasestablishedwithinertiaforceoflateralsecondaryflowtakenintoaccou