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

1、冲积河道冲刷过程中横向展宽的初步模拟摘要：本文首先描绘了冲积河道横向展宽现象，分析了河岸横向冲刷、崩塌的机理及展宽原因。然后建立细沙河流一维非均匀悬移质泥沙数学模型，同时结合河道横向展宽形式，预测冲积河道冲刷过程中，粘性土组成的河岸的横向展宽过程。此外，还分析了不同来水来沙条件及河岸土体特性对河道横向展宽和河床冲刷的影响，指出了影响冲积河道横向展宽的主要因素。关键词：冲积河道河槽横向展宽粘性河岸1引言在冲积河道中，河床包括河岸、滩地和主槽总是处在不断变化开展之中，其纵向变形和横向变形的程度取决于床沙质来量和水流挟沙力之间，以及河岸抗冲性能和水流冲刷才能之间的比照关系。在河流上修建水利工程后，河

2、床变形将更加明显。如水库在蓄水拦沙期间或非汛期，下泄水流的含沙量很小，下游河道往往发生长间隔 的冲刷。河道在纵向冲刷下切的同时，通常伴有横向展宽现象。据统计，三门峡水库蓄水运用期间，黄河下游由于受长时间的清水冲刷，造成滩地大量崩塌，其中最严重的是花园口至高村河段，约有200k2的滩地崩塌，滩地的大量崩塌使二滩之间的河槽宽度增加，如花园口至东坝头的二滩宽度由2563增加到36331。此外当水库淤满后，为恢复库容，通常降低坝前水位，采用泄空冲刷的运行方式来排出库区内淤沙。在泄空排沙过程中，库区内的主槽不断冲刷下切，滩槽高差不断加大，导致失稳崩塌，引起河槽宽度增加。由于岸滩的横向冲刷、崩塌，不仅严重

3、影响岸边土地的规划利用、居民的正常生活，而且还危及堤防的平安稳定。因此，模拟河道纵向冲刷过程时，同时模拟冲积河道的横向展宽过程，对水库和护岸工程的规划设计相当重要。但是现有的很多泥沙数学模型，很少模拟河道的横向变形。本文首先分析了冲积河道横向展宽的机理。在建立细沙河流的一维恒定非均匀泥沙数学模型的根底上，引入sanandThrne2,3提出的粘性河岸的横向展宽形式，模拟冲积河道在冲刷过程中的横向展宽过程。根据一概化河段的模拟结果，预测不同水沙条件下的库区河道横断面的变化过程，同时分析不同的河岸土体特性对横向展宽和冲刷过程的影响。2横向展宽机理2.1横向展宽的主要原因一般来讲，冲积河道的横向展宽

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

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

6、的河岸而言，确定要比无粘性组成的河岸复杂得多。粘性土的起动不同于较粗的散体泥沙，一般呈团粒、块状冲起或一块一块的剥落，即粘性土以团粒形式起动，破坏时呈块状。因此粘性土抗冲的临界切应力应为细颗粒泥沙成团起动的切应力，而不是单个细泥沙颗粒的起动切应力4。3一维泥沙数学模型的建立3.1根本方程水流连续方程(1-1)水流运动方程(1-2)悬移质不平衡输移方程(1-3)河床变形方程(1-4)式中Q为流量；A、B为过水面积和河宽；Z、h为断面平均水位、水深；Ask为第k粒径组泥沙的冲淤变形面积；Sk,S*k为断面悬移质分组含沙量和分组水流挟沙力；k为第k粒径组泥沙沉速；为泥沙干密度，g为重力加速度；k为第

7、k粒径组泥沙恢复饱和系数，n为糙率；x、t为间隔 及时间。为简化计算，采用非耦合解法，即先用二分法解水流方程(1-1)、(1-2)，求出有关水力要素后，再求解(1-3)式，得各断面含沙量。最后用(1-4)式计算河床冲淤面积。3.2有关参数的处理3.2.1水流挟沙力计算本文采用在黄河上应用较广、考虑因素较全面的张红武挟沙力公式5。该公式通过对二维水流单位水体的能量平衡方程式沿垂线积分，经分析整理得出包括全部悬移质泥沙在内的水流挟沙力公式，即S*=K(0.0022+SV)u3/s-/ghln(h/6D50)(2)上式单位均为kgs制。s、为泥沙和浑水容重；为Karan常数，与含沙量有关；SV为体积

8、比含沙量；为非均匀沙的代表沉速；D50为床沙的中值粒径，K、为参数，一般分别取2.5、0.62。上式不仅适用于一般挟沙水流，而且更适用高含沙紊流，充分考虑含沙量对挟沙才能的影响。经其别人分别开展的验证结果说明，该式计算精度明显优于其它公式。3.2.2挟沙力级配计算悬移质分组挟沙力级配是数学模型计算中的关键问题，然而由于天然河道中水沙条件的复杂性，至今仍未得到较好的解决。本文在计算中采用李义天提出的方法6，这种做法是在输沙平衡时，认为第k粒径组泥沙在单位时间内沉降在床面上的总沙量等于冲起的总沙量，然后根据垂线平均含沙量和河底含沙量之间的关系，确定悬沙级配和床沙级配的关系。其特点是同时考虑了水流条

9、件和床沙组成对挟沙力级配的影响。3.2.3床沙级配计算为模拟冲刷过程中的床沙粗化现象，本模型将床沙分为两大层，最上层的床沙活动层及其该层以下的分层记忆层。床沙活动层的厚度为H，相应的级配为Pbk。分层记忆层可根据实际情况共分n层，各层的厚度及相应的级配分别为Hn、Pnk。在计算中，当各粒径组泥沙发生淤积时，那么记忆层数相应增加，即为n+1层，且该层的级配为t时刻的床沙活动层级配Ptbk。当各粒径组泥沙均发生冲刷时，根据冲刷量的大小，记忆层数相应减少假设干层，且级配作相应的调整。某断面的各粒径组的冲淤厚度Hsk，及总的冲淤厚度Hs，那么床沙级配的调整计算，通常可分为以下2种情况计算：(1)总的冲

10、淤厚度Hs0的情况，那么床沙活动层的级配可用下式计算Pbkt+t=Hsk+Ptbk(Ht-Hs)/Ht+t(3-1)式中Ptbk、Pbkt+t分别为第t时刻、t+t时刻的床沙活动层的级配。Ht、Ht+t分别为第t时刻、t+t时刻的床沙活动层的厚度。(2)总的冲淤厚度Hs0的情况，那么床沙活动层的级配可用下式计算Pbkt+t=Hsk+PtbkHt+Hs.Prek/Ht+t(3-2)式中Ptbk、Pbkt+t、Ht、Ht+t同前面。Prek为假设干个记忆层内的床沙平均级配。3.2.4恢复饱和系数的取值恢复饱和系数，反映了悬移质不平衡输沙时，含沙量向水流挟沙力靠近的恢复速度。对于的取值，既与来水来沙

11、条件有关，也与河床断面形态有关，是一个非常复杂的参数，它在河床变形计算中起着很大的作用。但对如何详细取值，目前还没有统一的规定。尽管不少研究者8认为的取值应大于1，即认为为近底含沙量与垂线平均含沙量的比值。当前，许多模型所使用的恢复饱和系数的取值都是经历性的，如韩其为9等人通过实测资料率定计算，南方河流的在淤积时为0.25，冲刷时为1.0。象黄河这样的细沙多沙河流，由于泥沙的有效沉速比较小，在一样的条件下，淤积时含沙量的沿程衰减过程或冲刷时的含沙量沿程增加过程，要更慢一些。因此黄河上的应该比南方河流校如对不同的粒径组采用一样的值，那么河床冲淤强度与泥沙粒径或沉速成正比，当河床处于冲刷状态时，泥

12、沙粒径越粗，河床冲刷量越大；泥沙粒径越细，河床冲刷量越校结果使河床发生细化现象，这显然与实际情况不符。为此本文在计算中采用韦直林10在黄河下游泥沙数学模型中用的计算方法，即对不同的粒径组泥沙采用不同的值(4)式中k为第k粒径组的泥沙沉速，单位为/s。转贴于论文联盟.ll.论文联盟.LL.编辑。4河道横向展宽的计算形式为预测冲积河道演变过程中的河宽变化，过去常用的方法多为冲淤平衡方法11，包括许多经历公式、理论和假设。如通过分析、整理和统计天然河流实测资料得出的河床宽深关系式，如张红武的河相关系式5；如各种极值假说，常见的有水流功率最小12、水流能耗率最小13、输沙率最大14等。由于这些方法都是

13、用来预测河道由冲淤不平衡状态到冲淤平衡状态时河宽调整的大小，无法确定河道处于冲淤不平衡情况下的河宽变化。因此本文采用文献2,3提出的横向展宽计算形式，该形式为水动力学土力学的分析方法，即采用水动力学模型计算河床冲深，用土力学方法分析河岸是否会发生崩塌，主要用于预测粘性土体组成的河岸的横向展宽过程。详细计算过程如下。4.1横向冲刷间隔 在tse时间内，粘性河岸被水流横向冲刷后退的间隔 为B=t/602.185410-2(-)/se-1.3(5)式中s河岸土体的容重(KN/3)；B为t时间内河岸因水流侧向冲刷而后退的间隔 ()；为作用在河岸上的水流切应力(N/2)；为河岸土体的起动切应力(N/2)

14、。4.2河岸稳定性分析当河床冲深Z，河槽冲宽B后，相应的河岸高度增加，岸坡变陡，河岸稳定性降低，因此本文根据土力学中的边坡稳定性关系，用来计算河岸初次崩塌时临界条件下的河岸高度、破坏面和程度面的夹角等参数。在河岸稳定性分析中采用如下假定：如河岸土体由粘性沙组成且垂向分布均匀；河岸土体崩塌时的破坏面为斜面，且通过坡脚；在稳定性分析中不考虑其他因素的影响，尽管有些因素在特定情况下对河岸崩塌很重要；此外本文仅考虑河岸坡度大于60度的陡坡稳定情况。4.2.1河岸初次崩塌图1为河岸发生初次崩塌时的河岸形态，在初始河岸高度H0，初始河岸坡度i0的情况下，根据水动力学模型计算得到的床面冲刷深度Z，由式(5)

15、计算横向冲刷宽度B，确定冲刷后的河岸高度H1以及转折点以上的河岸高度H2，即可计算相对河岸高度的实测值(H1/H2)。当河岸发生初次崩塌时，破坏面与程度面的夹角为，即=0.5tg-1(H1/H2)(1.0-k2).tg(i0)+(6)式中k为河岸上部拉伸裂缝的深度Ht与河岸高度H1之比，一般取0.5；为河岸土体的内摩擦角。由6式求出后，便可采用土力学中的边坡稳定性分析，计算将要发生崩塌时相对河岸高度的分析解(H1/H2)，即图1河岸初次崩塌示意图Skethfinitialbankfailure(7)式中1=(1-k2)(sins-s2tg),2=2(1-k)/(sH2),3=(sinstg-s

16、in2)/tg(i0)，其中为河岸土体的凝聚力。4.2.2河岸二次崩塌(8)式中1=sins-s2tg，2=2(1-k)/(sH2)。在(H1/H2)、(H1/H2)的情况下，那么河岸的边坡稳定性分析，可采用类似河岸发生初始崩塌时的方法进展判断。此外，本文在计算中还采用以下假定：认为直接从河岸冲刷下来的和河岸崩塌后产生的泥沙都均匀的铺在床面上。也就是认为河床的实际冲刷量等于由悬移质不平衡输沙引起的冲刷量减去河岸的泥沙冲刷量(包括河岸土体崩塌的部分)。转贴于论文联盟.ll.5算例及结果分析由于目前缺少冲积河道横向展宽的实测资料，故本文采用一概化的模拟河段，研究冲积河道河槽展宽对冲刷过程的影响，同

17、时在考虑河岸可冲刷的情况下，分析了不同的来水来沙条件及河岸土体特性对河道展宽和河床冲刷过程的影响。5.1河道概况假定有一概化的顺直河段，长100k，河床初始纵比降为=0.0002，床沙中值粒径为0.038，河道初始糙率为0.02。河道断面形态为梯形，初始河岸坡度为60，初始河岸高度为8.0，河槽宽度为400。在河段出口段修建有一大坝。由于泥沙大量淤积在库区内，决定降低坝前水位，产生溯源冲刷，以便排出库内淤沙。在排沙过程中，坝前水位保持不变为105.0，上游来流量为2000.03/s，且水中根本不挟带泥沙，即含沙量根本为0。粘性河岸的土体特性：容重s为18KN/3、内摩擦角为14、凝聚力为20K

18、N/2、河岸土体的起动切应力为3.0N/2。在计算中，把模拟河段划分为100小段，计算时段恒为2天。假定计算初始时刻的全河段的床沙级配一样。5.2计算结果分析为研究冲积河道横向展宽对冲刷过程的影响，本文同时计算了河岸不可冲刷和河岸可冲刷两种情况。从左岸的坡脚位置变化图图3中可以看出：刚开始时，河道横向冲刷相当迅速，尤其在库尾流速较大的河段；而在坝前段，由于水深大，流速小，因此对河岸的冲刷相当弱。随着时间增加，这种冲刷趋势逐渐减缓。原因在于河槽的横向冲刷速率取决于水流作用在河岸上的切应力，当河道在冲刷过程中横向展宽后，过水面积增大，使相对减小，最后导致横向冲刷速率减校图3左岸坡脚位置变化Left

19、banktelatinvariatinindegradatin表1不同河岸类型的河床冲刷过程Degradatinpressesfdifferentbanktypes河岸类型冲刷天数河床冲深()河槽宽度()Vbed(3/)Vbank(3/)VTT(3/)Vbed/VTT(%)Vbank/VTT(%)(I)301.06400.0424.00.0424.0100.00.0(II)0.80407.7322.9114.1437.073.926.1(I)602.11400.0844.00.0844.0100.00.0(II)1.66414.9679.1199.4878.577.322.7(I)903.15

20、400.01260.00.01260.0100.00.0(II)2.50421.91029.2290.31319.578.022.0(I)1204.16400.01664.00.01664.0100.00.0(II)3.30428.61370.5385.81756.378.022.0河岸类型(I)：河岸稳定且不可冲刷；河岸类型(II)：河岸可冲刷Vbed：从河床上冲起的泥沙体积；Vbank：从河岸上冲走及崩塌的体积；VTT：总的冲刷体积图4表示进口断面的断面形态随时间的变化过程，从图中可知，由于粘性河岸高度相对较高，坡度较陡，故河岸在横向冲刷及河床纵向下切后，又会在重力作用下发生崩塌，河岸上部

21、分土体由于拉伸裂缝的形成，保持直立状态，河岸下部分沿一斜面滑动，这是因为当河岸发生在初次崩塌后，在随后的失稳破坏中，河岸以平行后退的方式崩塌。此外，从表1中还可以看出，在4个月中，当河岸不可冲刷时，河床总共冲深4.16，而当河岸可冲刷时，河床仅冲深3.30，两者相差近1.0。从冲刷量来看，在4个月中，从河岸冲刷下的泥沙数量为385.83/，为整个断面冲刷量的22%；从河床上冲起的泥沙量为1370.53/，低于河岸不可冲刷的情况。图4河岸可冲刷时断面(S=1)形态变化nfiguratinvariatinfrsssetinN.1图5床沙平均粒径沿程变化Variatinfeansizefbedate

22、rialindegradatin从上述分析来看，在冲积河道可横向展宽的情况下，由两个重要原因可以减缓河道冲刷：一是河槽宽度增大，导致过水面积增大，流速较小，水流挟沙力减小，从而减少对床面的冲刷。二是从河岸冲刷和崩塌下来的泥沙数量部分满足了悬移质水流挟沙力的要求，因此可减少从床面冲起的泥沙数量。此外，还有河床冲刷过程中产生的床沙粗化现象，可引起河道阻力增加，也可减缓对河道的冲刷。如图5表示床沙平均粒径的沿程变化情况。在水库泄空冲刷过程中，经过120天后，床沙有明显粗化现象，尤其在河床冲深较多的库尾段。5.3不同来水来沙条件对河道冲刷和横向展宽的影响众说周知，冲积河道的河床演变与上游的来水来沙条件

23、密不可分。上游来水来沙条件的变化，尤其在河岸可冲刷且不稳定的情况下，必然会对下游河床变形产生影响。图6(a)、6(b)分别表示不同来水条件下的河槽展宽间隔 和河床冲刷深度在120天后的变化情况。当来流流量从20003/s增加到30003/s时，河道在库尾段展宽间隔 增加10多，而坝前段展宽间隔 增加相对较少。而从河床冲刷深度的沿程变化情况来看，在库尾段冲刷深度增加较少，但在坝前段，河床冲深增加较多。这是因为当流量增加的过程中，在库尾段河岸冲刷量增加较多，部分满足了水流挟沙力的要求，故河床冲深较少；而在坝前段，由于河岸冲刷量增加较少，但流速增加多，即水流挟沙力有较大的增加，所以河床冲刷深度有较大

24、的增加。图6(a)不同来水条件的河槽展宽间隔 Variatinfidenningdistanefdifferentinflnditins图6(b)不同来水条件的河床冲刷深度Variatinfsuringdepthfdifferentinflnditins图7(a)、7(b)分别表示不同来沙条件下的河槽展宽间隔 和河床冲刷深度在120天后的变化情况。从图中可知，来流含沙量从0kg/3变化到2kg/3，沿程各断面的河槽展宽间隔 根本未变，而河床冲刷深度那么有较大的减少，尤其在库尾段。在坝前段，由于流速较小，从上游冲刷下来的含沙量较大，水流挟沙力根本到达饱和，故河床冲刷深度很小，有个别断面出现了淤积

25、情况。因此当上游来沙条件有较小的变化时，对下游河道横向展宽影响较小，而对河床纵向冲刷影响较大，这是因为上游来沙条件变化对下游河道的水流挟沙力的影响最显著，对其它条件影响较校从图6和图7比较来看，来水条件比来沙条件的变化对河道横向展宽的影响更显著。图7(a)不同来水条件的河槽展宽间隔 Variatinfidenningdistanefdifferentinflnditins图7(b)不同来水条件的河床冲刷深度Variatinfsuringdepthfdifferentinflnditins5.4不同类型的河岸土体对河道冲刷和横向展宽的影响从前面的分析可知，对粘性土组成的河岸而言，不同特性的河岸土

26、体对河道的冲刷过程有明显的影响。一般可用以下参数表示粘性河岸的土体特性：、s、和。这些参数取值不同，对河岸冲刷量、河槽展宽间隔 和河床冲刷深度的计算结果都有不同程度的影响。从图8(a)、8(b)来看，当河岸土体的抗冲性能减弱时，即从3.0N/2减小到2.5N/2后，河槽展宽间隔 在库尾段有明显增加，而在坝前段由于水深大，水流作用在河岸上的切应力小，因此两者增加较少。从河床的冲刷深度来看，河岸抗冲才能减弱，容易受水流横向冲刷，失稳崩塌的泥沙数量增多，因此部分可以弥补从床面冲起的泥沙，故河床冲刷深度反而减校图8(a)不同河岸抗冲性的河槽展宽间隔 Variatinfidenningdistanefr

27、differentsurresistanefbanksil图8(b)不同河岸抗冲性的河床冲刷深度Variatinfsuringdepthfrdifferentsurresistanefbanksil河岸土体的内摩擦角越大，表示河岸土体的抗剪强度大，稳定性越好，河岸越不易发生失稳崩塌。从图9(a)、9(b)来看，当从14增加到20后，河槽展宽间隔 根本不变，而河岸冲刷深度有一定程度的增加，这是由于从河岸冲刷或崩塌下来的泥沙数量减少的缘故。通过对其它参数的比照分析计算说明，如河岸土体的s、，在一定范围内变化时，对河岸的冲刷量、河道展宽间隔 及河床的冲刷深度都没有明显的影响。图9(a)不同河岸土体的

28、河槽展宽间隔 Variatinfidenningdistanefrdifferentbanksilfritinangle图9(b)不同河岸土体的河床冲刷深度Variatinfsuringdepthfrdifferentbanksilfritinangle6结语本文把一维泥沙数学模型同河道的横向展宽形式结合，模拟冲积河道冲刷过程中的横向展宽过程，由于影响粘性河岸冲刷、失稳崩塌的因素比较复杂，在采用假设干假定的根底上，得到以下结论。1.在河岸可冲刷的条件下，冲积河道的横向展宽可以大大减缓对河床的纵向下切速率。2.不同的来水来沙条件，对河道展宽和河床的冲刷下切影响不同。来水条件对河道的展宽影响较大，

29、而来沙条件对河床冲刷下切的影响较大。3.不同土体特性的粘性河岸对河道横向展宽和河床纵向冲刷的影响不同，尤其是河岸抗冲的临界切应力对两者的影响最大，其它参数如河岸土体容重、内摩擦角及凝聚力等对两者的影响较校尽管这些参数之间互相联络，但由于对计算结果影响大，因此必须正确确实定参数。4.本文假定河道断面形态为梯形，因此河岸稳定性分析比较简单。而对于天然河道的实际断面形态，不但稳定性分析复杂，而且解决从河岸冲刷下来的泥沙在横向如何分布问题，以及转化为床沙和悬移质的比例问题，都是相当困难的，这些问题都有待于进一步的解决。转贴于论文联盟.ll.参考文献1赵业安等。黄河下游河道演变根本规律。河南：黄河水利出版社，1998年。2sanandThrne,Riverban