沙波与刚性非淹没植物复合床面推移质泥沙输移特性研究_叶志恒.pdf
《沙波与刚性非淹没植物复合床面推移质泥沙输移特性研究_叶志恒.pdf》由会员分享,可在线阅读,更多相关《沙波与刚性非淹没植物复合床面推移质泥沙输移特性研究_叶志恒.pdf(5页珍藏版)》请在咨信网上搜索。
1、第4 1卷第5期2023年5月水 电 能 源 科 学W a t e r R e s o u r c e s a n d P o w e rV o l.4 1 N o.5M a y 2 0 2 3D O I:1 0.2 0 0 4 0/j.c n k i.1 0 0 0-7 7 0 9.2 0 2 3.2 0 2 2 1 2 8 9沙波与刚性非淹没植物复合床面推移质泥沙输移特性研究叶志恒1,唐立模2,倪培桐1,苗 青1(1.广东省水利水电科学研究院,广东 广州5 1 0 6 1 0;2.河海大学水利水电学院,江苏 南京 2 1 0 0 9 8)摘要:鉴于复合床面条件下刚性非淹没植物改变推移质输沙
2、特性的相关研究成果较少,为此采用水槽试验和理论分析方法,研究了平行排列的刚性非淹没植物对泥沙起动谢尔兹数和推移质输沙强度的影响。结果表明,平行排列的刚性非淹没植物使得泥沙起动谢尔兹数增大了9 7.2%,推移质输沙率降低了8 0%以上。通过理论分析提出了适用于刚性非淹没植物平行排列条件下恩格隆推移质输沙率公式,经验证该公式具有较高的精确度。关键词:沙波;刚性非淹没植物;复合床面;推移质泥沙;输移特性中图分类号:TV 1 4 3 文献标志码:A 文章编号:1 0 0 0-7 7 0 9(2 0 2 3)0 5-0 0 4 3-0 5收稿日期:2 0 2 2-0 6-2 2,修回日期:2 0 2 2
3、-0 8-0 9基金项目:广东省水利科技创新项目(2 0 2 1-0 7)作者简介:叶志恒(1 9 9 3-),男,硕士、工程师,研究方向为河流动力学E-m a i l:y e z h i h e n g 4 2 51 6 3.c o m1 引言河流中每颗植水生植物均可视为圆柱绕流,当其成片出现后形成了复杂的“树群”效应,改变了流态,进而改变了泥沙输移特性,而泥沙输移特性的改变又影响着流态,三者相互作用形成了一个复杂的相互反馈模式。早在1 8 7 0年代L I R M等1研究了植物的不同排列方式对输沙的影响,指出采用梅花形交错排列阻水效果和拦沙能力明显优于平行排列;YAG E R E M等 2
4、采用刚性圆柱体模拟挺水植物,研究了挺水植物梅花状排列对推移质输沙率的影响,指出含植物条件下推移质输沙率增加或减小取决于平均流速和边界剪切应力;孙志林等3通过水槽试验表明水生植物的存有利于泥沙淤积。2 0世纪末,毕慈芬等4将沙棘种植在河沟中形成“柔性坝”,通过水槽试验研究了“柔性坝”促使泥沙沉积的原因和规律;邱秀云等5采用塑料草模拟沙棘等柔性植物进行试验,表明底坡越大,植物坝上游泥沙运动速度的衰减率就越小;杨庆楠等6发现植物根系密度、直径、表面积密度等特征与产沙量存在负相关线性关系。上述研究表明植物的存在改变了泥沙的输移特性,植物和沙波的复合床面条件下泥沙输移特性的研究可为湿地保护、江心洲行洪安
5、全、滩地冲刷等提供科学的理论依据。为此,本文通过水槽试验和理论分析,探求了复合床面下刚性非淹没植物平行排列对泥沙起动及推移质输沙率的影响程度和机理。2 试验布置与方法2.1 试验布置和测控系统试验在长1 2 m、宽0.6 m、高0.6 m的可变坡矩形水槽中进行,末端设置长1.3 m、宽0.6 m的沉砂池,并安装0.1 g的电子天平,用以实时称量推移质泥沙水下质量。采用直径0.6 c m、高度2 1 c m的圆柱形刚棒模拟刚性植物,以间隔5 c m平行排列于厚度1 c m的塑料板中,植物密度采用植株比法7:=N d2/(4S)(1)其中S=B L式中,N为植物株数;d为植株直径;S为床面面积;B
6、为宽度;L为长度,本文试验植物密度为0.0 0 2 8 2 7。试验设置2种床面形态(图1):工况1。沙波和刚性非淹没植物复合床面;工况2。单一沙波床面。床沙区始于水槽入口2 m处,铺设厚度9 c m,总长度8 m,于床沙区最前铺设长2 m、坡比为15的过渡区。试验区距入口68 m处。在工况1中设置植物区,植物区始于水槽入口4 m处,铺设总长度6 m。(b)2工况(a)1工况2 m2 m2 m6 m12 mV刚性非淹没植物9 cm厚床沙沉砂池沉砂池9 cm厚床沙2 m2 m2 m6 m12 m图1 工况1、2试验布置示意图F i g.1 S c h e m a t i c d i a g r
7、a m o f e x p e r i m e n t a l l a y o u t o f w o r k i n g c o n d i t i o n s 1 a n d 2试验前将泥沙铺设为初始工况,通过调节流量、坡降和尾门保证试验区水深(1 0 c m)延程不变,此时认为试验段已达到均匀流状态。采用电磁流量计、水槽刻度尺和电子天平读取试验数据。2.2 床面泥沙试验采用中值粒径D5 0=0.3 8 1 mm、不均匀系数1.5、比重为2.6 5的石英砂。根据沙玉清8过渡区公式计算得到泥沙沉速为5.5 0 c m/s。一般条件下植物的存在将减小泥沙沉速,当植物面积密度在0.0 7 1及以
8、下且粒径小于0.2 c m的颗粒可不考虑植物的影响9。因此可认为复合床面泥沙沉速为5.5 0 c m/s。根据泥沙特性分区图1 0可判断本文试验水流强度泥沙的运动以在床面层内的滑动、滚动、跳跃为主。3 试验结果3.1 泥沙起动流速试验采用K R AME R H1 1提出的弱动作为泥沙起动的判别标准,通过测量泥沙起动时的流速及公式验证,得出2种工况下泥沙起动流速值。(1)工况1。当整个床面区域泥沙均开始弱动,测量得到泥沙起动流速为1 8.6 5 c m/s。采用T ANG HONGWU等1 2刚性非淹没植物条件下泥沙起动流速公式:Uv c=0.3 6s-g DhD()1/6dh D/4-()0.
9、3 1 9(2)式中,Uv c为刚性非淹没植物条件下泥沙起动流速;s为泥沙比重,2 6.5 k N/m3;为水比重,9.8 k N/m3;g为重力加速度,9.8 1 m/s2;D为泥沙颗粒直径,0.3 8 11 0-3m;h为水深,0.1 m。计算得到复合床面泥沙起动流速为1 7.4 0 c m/s,与实测值误差为7.2%,该误差值与泥沙密实度、起动的判别标准等有关。本文将采用实测值1 8.6 5 c m/s作为复合床面泥沙起动流速。(2)工况2。缓慢增大流速至床面泥沙起动,测量得到泥沙起动流速为2 3.9 0 c m/s。采用沙莫夫、冈恰洛夫和张瑞瑾公式1 0计算该工况下泥沙起动流速范围为2
10、 2.7 02 8.0 0 c m/s,实测值处于计算范围内,因此采用实测值作为单一沙波床面泥沙起动流速。3.2 推移质输沙率(1)工况1。王浩1 3认为含刚性非淹没植物条件下当水流流速达到起动流速后,泥沙经过5 h冲刷可达到平衡。试验过程中发现当冲刷时长达到5 h后,电子天平中泥沙质量增长率趋于稳定,以56 h之内的推移质输移泥沙水下质量计算单宽输沙率:q b=W/(B T)(3)式中,qb为以水下重量计算的单宽推移质输沙率;W 为第56 h之内实测泥沙水下重量;B为水槽宽度,0.6 m;T为推移质输沙时间,3 6 0 0 s,计算得到推移质输沙率见表1。表1 工况1推移质输沙率T a b.
11、1 B e d l o a d t r a n s p o r t r a t e i n w o r k i n g c o n d i t i o n 1组次植物区平均流速Uv比降J/输沙率qb组次植物区平均流速Uv比降J/输沙率qb11 8.6 41.4 2062 3.8 72.1 4 0.0 0 0 9 8 121 9.6 61.5 50.0 0 0 1 8 072 4.2 92.2 1 0.0 0 1 1 3 032 0.4 91.6 60.0 0 0 3 2 382 6.1 02.5 1 0.0 0 1 7 5 042 1.7 61.8 30.0 0 0 5 9 592 7.9 8
12、2.9 1 0.0 0 3 4 7 052 2.6 71.9 60.0 0 0 6 7 6注:植物区平均流速、输沙率单位分别为c m/s、k g/(sm)。(2)工况2。采用电子天平测量第56 h内推移质水下重量,得到无植物工况下的推移质输沙率,试验结果见表2。表2 工况2推移质输沙率T a b.2 B e d l o a d t r a n s p o r t r a t e i n w o r k i n g c o n d i t i o n 2组次流速/(ms-1)比降J/输沙率qb/(k g(sm-1)1 0.2 5 00.3 1 61.2 31 0-52 0.2 6 00.3 5
13、00.0 0 0 1 3 63 0.2 7 90.4 2 50.0 0 0 9 8 24 0.2 8 60.4 5 80.0 0 1 4 5 05 0.3 0 90.5 5 80.0 0 2 7 5 06 0.3 3 30.6 6 60.0 0 3 9 0 07 0.3 5 70.7 8 30.0 0 6 2 5 04 试验结果分析4.1 含刚性非淹没植物剪切力线性叠加分析方法含植物条件下床面总剪切力表达式1 4为:=g h J(1-h*)(4)44水 电 能 源 科 学 2 0 2 3年 第4 1卷第5期叶志恒等:沙波与刚性非淹没植物复合床面推移质泥沙输移特性研究其中h*=l/h式中,为单位
14、床面总剪切应力;为水的密度;J为水力坡降;h为水深;l为植物被淹没高度,非淹没植物h*=1。GA R D E R J等1 5指出复合床面工况下总剪切力可分解为植物剪切力和床面剪切力。其中刚性非淹没植物剪切应力表达式为:v=CD vN d l U2v/2(5)式中,CD v为植物拖曳力系数,无量纲数;N为单位面积上植株个数;d为刚性非淹没植物直径;Uv为植物区平均流速。S TON E B M等 1 6指出植物区水流平均速度Uv可用含植物明渠水流中平均流速U替代,用以求解植物拖曳力。植物区流速可表示为:Uv=U=q/(1-)(6)式中,q为单宽流量。对同一水力条件下植物剪切力仅与被淹没高度有关,植
15、物剪切力沿水深呈线性分布。同时B a g n o l d假设挟沙水流的总剪切力沿水深亦呈线性分布1 2,因此本文将总剪切应力沿水深进行分解为植物剪切力和床面剪切力线性叠加而成,作用于床面的剪切力影响着泥沙运动,称之为有效剪切力。4.2 植物谢尔兹数计算方法将剪切力无量纲化,得到总谢尔兹数与植物谢尔兹数v表达式为:=g h J(1-h*)/(s-)D(7)v=CD vN d l V2v/2(s-)D(8)4.3 有效谢尔兹数计算方法有效谢尔兹数是指有效剪切力无量纲化后的数值。对含刚性非淹没植物条件下,有效谢尔兹数e v为总谢尔兹数除去植物谢尔兹数v:e v=-v(9)对工况2,有效谢尔兹数e即为
16、总谢尔兹数:e=(1 0)4.4 输沙强度计算将推移质输沙率换算成无量纲推移质输沙强度,采用下式:=q b/sD(s-)g D/(1 1)计算得到工况1、2输沙强度和有效谢尔兹数(表3、4)。计算得到复合床面条件下临界谢尔兹数为c v=0.0 7 1,单一沙波条件下临界谢尔兹数为c v=0.0 3 8。4.5 刚性非淹没植物对泥沙输移特性的影响根据试验数据可知,复合床面泥沙起动谢尔表3 复合床面谢尔兹数与输沙强度计算T a b.3 C a l c u l a t i o n o f s h i e l d s n u m b e r a n d s a n d t r a n s p o r
17、t i n t e n s i t y o f t h e c o m p o u n d b e d f o r m组次植物区平均流速Uv相对流速Uv/Uv c总谢尔兹数植物谢尔兹数v有效谢尔兹数e vv输沙强度e11 8.6 41.0 00.2 2 50.1 5 40.0 7 10.6 8 4021 9.6 61.0 50.2 4 60.1 7 10.0 7 50.6 9 6 0.0 0 2 2 732 0.4 91.1 00.2 6 30.1 8 50.0 7 80.7 0 5 0.0 0 4 0 842 1.7 61.1 70.2 9 10.2 0 80.0 8 30.7 1 7 0.
18、0 0 7 5 152 2.6 71.2 20.3 1 10.2 2 50.0 8 60.7 2 6 0.0 0 8 5 362 3.8 71.2 80.3 4 00.2 4 90.0 9 10.7 3 5 0.0 1 2 4 072 4.2 91.3 00.3 5 00.2 5 80.0 9 20.7 3 7 0.0 1 4 3 082 6.1 01.4 00.3 9 80.2 9 70.1 0 10.7 4 6 0.0 2 2 1 092 7.9 81.5 00.4 6 10.3 4 00.1 2 10.7 3 8 0.0 4 3 7 0注:植物区平均流速Uv单位为c m/s。表4 单一沙
19、波床面谢尔兹数与输沙强度计算T a b.4 C a l c u l a t i o n o f s h i e l d s n u m b e r a n d s a n d t r a n s p o r t i n t e n s i t y o f s i n g l e s a n d d u n e b e d f o r m组次流速/(c ms-1)相对流速U/Uc总谢尔兹数输沙强度(干)1 2 5.0 01.0 50.0 3 8 20.0 0 0 1 5 52 2 6.0 01.1 00.0 4 2 20.0 0 1 7 2 03 2 7.9 01.1 70.0 5 1 30.0
20、 1 2 4 0 04 2 8.6 01.2 00.0 5 5 30.0 1 8 3 0 05 3 0.9 01.3 00.0 6 7 40.0 3 4 6 0 06 3 3.3 01.4 00.0 8 0 50.0 4 9 3 0 07 3 5.7 01.5 00.0 9 4 50.0 7 8 9 0 0兹数为单一床面的1.9 7倍,刚性非淹没植物存在使得泥沙更难以起动。试验范围内,复合床面与单一沙波床面输沙强度()均随谢尔兹数()的增大而增大(图2)。相同有效谢尔兹数条件下,复合 床面推移质 输沙率仅为 单一沙波床 面的101010102101-=1/)(10-410-310-210-11
21、00工况工况12图2 输沙强度与水流参数散点F i g.2 S c a t t e r p o i n t o f s e d i m e n t t r a n s p o r t i n t e n s i t y a n d w a t e r f l o w p a r a m e t e r s1 0%2 0%。由于复合床面下总谢尔兹数中植物谢尔兹数占比高达6 8%以上,水流总剪切力中大部分作用于刚性非淹没植物,减小了作用于床面泥沙的剪切力,使得推移质输沙强度减小。此外水流的紊动猝发对床面泥沙推移层厚度有着较大影响1 7,复合床面条件下,植物根茎成为拟序结构发展的障碍物,部分猝发过程
22、并非完整的喷射-清扫周期1 8,因此本文认为刚性非淹没植物的存在减小了推移层厚度,导致相近有效剪切力条54件下含刚性非淹没植物推移质输沙强度小于单一沙波床面推移质输沙强度。采用剪切力为特征指标分析刚性非淹没植物对泥沙输移特性影响,这与恩格隆推移质输沙率公式采用的特征指标一致1 0,因此借鉴恩格隆,通过本次试验数据,得出适用于刚性非淹没植物条件下恩格隆推移质输沙率公式:=1 1.6(-c)-0.7c()(1 2)式中,为输沙强度;c为泥沙起动临界谢尔兹数;为谢尔兹数。经典恩格隆公式适用于单一沙波床面条件下推移质输沙率,式(1 2)中谢尔兹数即为复合床面中有效谢尔兹数。D E S UT T E R
23、 R等1 9指出,推移质公式在使用前均需应用当地的推移质实测资料重新标定公式中的参数,可将经典恩格隆公式改写为:v=K e v-c v()e v-c v()(1 3)式中,K、为复合床面恩格隆推移质输沙率公式参数;v为复合床面条件下为输沙强度。公式中K、与推移质运动所处的流速、床面动摩擦系数等特征要素有关。由于刚性非淹没植物的存在增大了河道糙率2 0改变了泥沙运动特性,为得到适用于本次复合床面推移质泥沙输沙率公式需采用本次试验结果重新率定两个参数。采用最小二乘法率定得到复合床面恩格隆推移质输沙率公式:v=4.2e v-c v()e v-0.6 3c v()(1 4)式(1 4)与试验成果拟合度
24、较高(图3)。实测值0.050.040.030.020.0100 01.0 02.0 03.0 04.0 05.理论值实测值yx=图3 输沙率实测值与理论值对比F i g.3 C o m p a r i s o n o f m e a s u r e d a n d t h e o r e t i c a l v a l u e s o f b e d l o a d t r a n s p o r t采用UME S H C KOTHYA R I等1 4实测数据对式(1 4)精确度进行检验(图4),检验结果表明式(1 4)对K o t h y a r i数据总体上拟合度较高,实测数据较为均分布
25、在曲线两侧,仅部分实测点与曲线偏差较大,这与不同试验工况下泥沙临界谢尔兹数、泥沙密实程度和泥沙特性有关。推移质输沙公式精度量化指标计算公式2 1为:10-410-310-210-1100101102101010101-=1/)0V0V(工况公式1Kothyari14v图4 公式(1 4)与室内试验成果对比F i g.4 C o m p a r i s o n o f f o r m u l a(1 4)w i t h l a b o r a t o r y e x p e r i m e n t a l r e s u l t sRRMS E R=1Nni=1Xi-YiYi()2(1 5)式中
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 刚性 淹没 植物 复合 推移 泥沙 特性 研究 叶志恒
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。