冲积河道冲刷过程中横向展宽的初步模拟论文

《冲积河道冲刷过程中横向展宽的初步模拟论文》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、冲积河道冲刷过程中横向展宽的初步模拟论文摘要：本文首先描述了冲积河道横向展宽现象，分析了河岸横向冲刷、崩塌的机理及展宽原因。然后建立细沙河流一维非均匀悬移质泥沙数学模型，同时结合河道横向展宽模式，预测冲积河道冲刷过程中，粘性土组成的河岸的横向展宽过程。此外，还分析了不同来水来沙条件及河岸土体特性对河道横向展宽和河床冲刷的影响，指出了影响冲积河道横向展宽的主要因素。关键词：冲积河道河槽横向展宽粘性河岸1引言在冲积河道中，河床（包括河岸、滩地和主槽）总是处在不断变化发展之中.freel2的滩地崩塌，

2、滩地的大量崩塌使二滩之间的河槽宽度增加，如花园口至东坝头的二滩宽度由2563m增加到3633m1。此外当水库淤满后，为恢复库容，通常降低坝前水位，采用泄空冲刷的运行方式来排出库区内淤沙。在泄空排沙过程中，库区内的主槽不断冲刷下切，滩槽高差不断加大，导致失稳崩塌，引起河槽宽度增加。由于岸滩的横向冲刷、崩塌，不仅严重影响岸边土地的规划利用、居民的正常生活，而且还危及堤防的安全稳定。因此，模拟河道纵向冲刷过程时，同时模拟冲积河道的横向展宽过程，对水库和护岸工程的规划设计相当重要。但是现有的很多泥沙数学

3、模型，很少模拟河道的横向变形。本文首先分析了冲积河道横向展宽的机理。在建立细沙河流的一维恒定非均匀泥沙数学模型的基础上，引入OsmanandThorne2,3提出的粘性河岸的横向展宽模式，模拟冲积河道在冲刷过程中的横向展宽过程。根据一概化河段的模拟结果，预测不同水沙条件下的库区河道横断面的变化过程，同时分析不同的河岸土体特性对横向展宽和冲刷过程的影响。2横向展宽机理2.1横向展宽的主要原因一般来讲，冲积河道的横向展宽过程是河岸土体和近岸水流相互作用的结果，即河岸土体的抗冲性和近岸水流的冲刷侵蚀力

4、之间相互作用的结果。除了河岸上植被生长情况、河道内水位变化、渗流、管涌等因素影响河道横向展宽外，但主要是以下两种情况，是导致河道横向展宽的主要原因：(1)通过水流直接横向冲刷导致展宽。当近岸水流切应力大于土体抗冲的切应力时，岸坡上的泥沙颗粒(或团粒)将被水流带走，从而导致河道展宽。(2)通过河岸崩塌导致展宽。对粘性土组成的河岸而言，主要以这种方式展宽，即床面冲刷导致岸高增加，或者水流淘刷河岸坡脚，使岸坡变陡，都会降低稳定性，当土体内部的滑动力大于其抗滑力时，河岸部分土体发生失稳、崩塌导致河宽增加

5、。对于非粘性土组成的岸滩而言，当岸坡坡度超过天然情况下泥沙的水下休止角时，边坡失去稳定，发生滑动。2.2影响横向展宽的主要因素影响冲积河道横向展宽的主要因素有：作用于近岸土体上的水流切应力、近岸土体的几何、物理化学特性（包括抗冲的临界切应力、河岸高度等）。对作用于近岸上的水流切应力τ而言，在某一特定的过水断面下，通常为断面平均水流切应力的0.75～0.77倍，但是由于河床冲深，过水断面形态改变。因此上述关系式就不可能存在，因此如果不能较好的估计河岸切应力的分布，可用断面的平均水流切应力代替。对于

6、无粘性土组成的河岸而言，土体以颗粒的状态被冲走，通常河岸土体抗冲的临界切应力τc可用泥沙颗粒在斜坡上的起动拖曳力表示。对于粘性土组成的河岸而言，确定τc要比无粘性组成的河岸复杂得多。粘性土的起动不同于较粗的散体泥沙，一般呈团粒、块状冲起或一块一块的剥落，即粘性土以团粒形式起动，破坏时呈块状。因此粘性土抗冲的临界切应力应为细颗粒泥沙成团起动的切应力，而不是单个细泥沙颗粒的起动切应力4。3一维泥沙数学模型的建立3.1基本方程水流连续方程(1-1)水流运动方程(1-2)悬移质不平衡输移方程(1-3)河

7、床变形方程(1-4)式中Q为流量；A、B为过水面积和河宽；Z、h为断面平均水位、水深；Ask为第k粒径组泥沙的冲淤变形面积；Sk,S*k为断面悬移质分组含沙量和分组水流挟沙力；ωk为第k粒径组泥沙沉速；ρ′为泥沙干密度，g为重力加速度；αk为第k粒径组泥沙恢复饱和系数，n为糙率；x、t为距离及时间。为简化计算，采用非耦合解法，即先用二分法解水流方程(1-1)、(1-2)，求出有关水力要素后，再求解(1-3)式，得各断面含沙量。最后用(1-4)式计算河床冲淤面积。3.2有关参数的处理3.2.1水流

8、挟沙力计算本文采用在黄河上应用较广、考虑因素较全面的张红武挟沙力公式5。该公式通过对二维水流单位水体的能量平衡方程式沿垂线积分，经分析整理得出包括全部悬移质泥沙在内的水流挟沙力公式，即S*=K(0.0022+SV)u3/κγs-γm/γmghωmln(h/6D50)m(2)上式单位均为kgms制。γs、γm为泥沙和浑水容重；κ为Karman常数，与含沙量有关；SV为体积比含沙量；ωm为非均匀沙的代表沉速；D50为床沙的中值粒径，K、m为参数，一般分别取2.5、0.62。上式不仅适用于一般挟沙水流