王延貴

(中國(guó)水利水電科學(xué)研究院國(guó)際泥沙研究培訓(xùn)中心,北京 100048)

河流岸灘挫落崩塌機(jī)理及其分析模式

王延貴

(中國(guó)水利水電科學(xué)研究院國(guó)際泥沙研究培訓(xùn)中心,北京 100048)

為探討岸灘挫落崩塌發(fā)生的機(jī)理和分析模式,通過(guò)分析河流岸灘土壤特性和地質(zhì)結(jié)構(gòu),結(jié)合岸灘土體應(yīng)力分布與變化特點(diǎn),深入分析了岸灘挫落崩塌的力學(xué)機(jī)制和崩塌過(guò)程,指出在河床沖刷下切過(guò)程中,黏性岸灘頂部出現(xiàn)的縱向裂隙會(huì)促使岸灘挫落崩塌的發(fā)生,發(fā)生過(guò)程相對(duì)簡(jiǎn)單和短暫,破壞面為平面;根據(jù)岸灘崩塌體力學(xué)平衡分析,探討了岸灘挫落崩塌的分析模式,導(dǎo)出了挫落崩塌的穩(wěn)定分析方程,該方程反映了河岸土壤特性、岸坡形態(tài)、河床沖刷、岸灘滲流等因素的共同作用關(guān)系,將其應(yīng)用于岸灘臨界崩塌高度公式的推導(dǎo),并利用實(shí)測(cè)資料對(duì)挫落崩塌分析模式的合理性進(jìn)行了檢驗(yàn),結(jié)果表明該分析模式是合理的,可以在實(shí)踐中應(yīng)用。

挫落崩塌;崩塌機(jī)理;臨界崩塌高度;岸灘;分析模式

河流岸灘崩塌是河道側(cè)向演變的重要形式,在沖積河流中普遍存在。岸灘崩塌是河道水流動(dòng)力與河岸土質(zhì)邊界結(jié)構(gòu)作用的結(jié)果,其中河岸土質(zhì)邊界結(jié)構(gòu)條件是岸灘崩塌的內(nèi)因,河道水流動(dòng)力條件是岸灘崩塌的外因。岸灘挫落崩塌是重要的崩岸類型,除與河岸形態(tài)、土壤性質(zhì)等有關(guān)外,還與河道水沙條件、河床沖刷狀態(tài)等有關(guān)[1-2]。對(duì)于具有縱向裂隙的黏性土岸灘,當(dāng)河床沖刷下切到一定程度,岸坡土體下部失去支撐,岸坡發(fā)生挫落崩塌。河流上修建水庫(kù)后,攔截了大量泥沙,改變了進(jìn)入下游河道的水沙過(guò)程,導(dǎo)致下游河道嚴(yán)重沖刷和下切,促使岸灘挫落崩塌發(fā)生。長(zhǎng)江三峽蓄水以來(lái),造成下游荊江河道大幅度沖刷,沙市三八灘、文村夾、石首彎道進(jìn)口、調(diào)關(guān)彎道、監(jiān)利烏龜洲等河段深泓高程降低幅度較大,如文村夾附近1998年10月至2009年10月間最大沖深11.2 m,致使三峽水庫(kù)下游河道崩岸狀況發(fā)生變化,三峽水庫(kù)試驗(yàn)性蓄水期2010年和2011年長(zhǎng)江中下游發(fā)生崩岸險(xiǎn)情分別為64處和65處,崩岸長(zhǎng)度分別為4.8km和44.8km,在荊江河段新沙洲、荊江門、張家墩、界牌等發(fā)生較大的崩岸[3-4],其中挫落崩塌占有相當(dāng)?shù)谋戎?給兩岸防洪及建筑物的安全產(chǎn)生重要影響。因此,深入研究河道沖刷下切對(duì)河流岸灘挫落崩塌的影響機(jī)理及崩塌分析模式是非常重要的。

1 岸灘挫落崩塌機(jī)理

1.1 挫落崩塌的力學(xué)機(jī)制

岸灘崩塌發(fā)生與否主要是由岸灘崩體下滑力(矩)與阻滑力(矩)的平衡關(guān)系決定的,當(dāng)下滑力(矩)大于阻滑力(矩)時(shí),岸灘崩體失穩(wěn)發(fā)生崩塌。岸灘崩體受力狀況與岸灘土體內(nèi)的應(yīng)力變化有很大的關(guān)系,崩體應(yīng)力場(chǎng)受到外部條件的影響而發(fā)生變化。河床沖刷下切后,岸灘沖刷部分失去支撐作用,岸灘內(nèi)應(yīng)力場(chǎng)也將發(fā)生變化。王延貴等[1,5-7]曾利用土力學(xué)擋土墻的破壞原理揭示河床沖刷下切后的岸灘挫落崩塌機(jī)理,結(jié)合土力學(xué)理論給出岸灘表面以下土體的應(yīng)力分布,如圖1所示。河床沖刷過(guò)程中,岸灘崩體內(nèi)應(yīng)力變化經(jīng)歷了彈性變形狀態(tài)到塑性破壞的極限平衡狀態(tài),河床沖刷初期,由于河流岸灘土體一側(cè)的泥沙逐漸被沖走,其坡度變陡,岸灘土體水平壓力減小,此階段處于彈性變形狀態(tài)(圖1(a))。若河床進(jìn)一步?jīng)_刷,岸灘土體將發(fā)展為塑性極限平衡狀態(tài)(圖1(b)),此時(shí)岸灘頂部會(huì)出現(xiàn)拉應(yīng)力,并出現(xiàn)與岸線平行的張性裂隙,張性裂隙深度H′可由極限平衡狀態(tài)下水平側(cè)向壓力分布求得[1,5-7],即

圖1 河床沖刷下切后岸灘挫落崩塌力學(xué)機(jī)制示意圖

式中和圖1中:σ為岸灘土體應(yīng)力;σz為岸灘土體自重應(yīng)力;σx為岸灘土體水平側(cè)向壓力;τ為土體切應(yīng)力;τf為土體抗剪強(qiáng)度;c為土體黏聚力;θ為土體內(nèi)摩擦角;ρ為土體密度。據(jù)文獻(xiàn)[1,5]的研究,岸灘頂面裂隙的產(chǎn)生會(huì)促進(jìn)岸灘挫落崩塌的發(fā)生,裂隙越深,岸灘崩體越不穩(wěn)定,挫落崩塌越易發(fā)生。

1.2 挫落崩塌的發(fā)生過(guò)程

對(duì)于黏性土壤岸灘而言,一方面,由于岸灘頂部土壤的收縮及張拉應(yīng)力的作用而產(chǎn)生裂縫[6-7];另一方面,河床沖刷下切會(huì)促使岸灘崩體內(nèi)的應(yīng)力發(fā)生變化,使岸灘頂部產(chǎn)生張性裂隙。表面張性裂隙將會(huì)使得岸灘的穩(wěn)定性降低,當(dāng)河床沖刷達(dá)到一定深度時(shí),岸灘土體下部將逐漸失去支撐,使得崩體沿著破壞面滑落而形成挫落崩塌。岸灘挫落崩塌一般沿河岸呈條形,橫向崩塌寬度較小,發(fā)生挫落崩塌的過(guò)程較簡(jiǎn)單和短暫,破壞面為平面。挫落崩塌發(fā)生后,新出露的岸壁直立,仍可發(fā)生新的挫落崩塌。岸灘崩體沿破壞面進(jìn)入河道的方式主要有旋轉(zhuǎn)倒入河內(nèi)的倒塌(圖2(a))和平滑入河的滑塌(圖2(b))兩種模式[8-9]。由于這兩種崩塌入河模式的崩塌機(jī)理并無(wú)不同,因此其崩體穩(wěn)定評(píng)價(jià)分析方法應(yīng)是相同的。挫落崩塌崩體穩(wěn)定評(píng)價(jià)可采用崩體受力平衡法進(jìn)行分析,通過(guò)計(jì)算崩體的安全系數(shù),判斷岸灘崩體的穩(wěn)定性。

圖2 挫落崩塌模式與發(fā)生過(guò)程示意圖

2 岸灘挫落崩塌的分析模式

2.1 挫落崩塌模式分析

對(duì)于黏性岸灘頂面出現(xiàn)縱向裂隙的情況,崩體最可能沿著滑裂面崩塌,如圖3所示。岸灘崩體作為一個(gè)整體,且假定為剛性,崩體所受的力主要包括有效重力W、破壞面處的支撐力N和阻滑力Pτ、滲透力Pd及外力P0。設(shè)整個(gè)岸灘崩體斷面面積為A,崩塌斷面滲透線以上和以下的崩體斷面面積分別為Au和Ad,且

式中:H和α分別為河岸高度和破壞面的傾角;α0為沖刷后岸腳段的坡角;Hi和αi分別為河岸斷面第i(i= 1,2)個(gè)折點(diǎn)和折線的深度與坡角,具體如圖3所示。

圖3 挫落崩塌穩(wěn)定分析示意圖

崩體的有效重力為

式中:ρ和ρsat分別為土體的密度和飽和密度;ρw為水的密度`;g為重力加速度。

崩體所受的滲透水動(dòng)力為

式中:J為土體內(nèi)滲流梯度。

聯(lián)解x和y軸向力的平衡方程:

式中:DF為崩體的下滑力。

崩體破壞面上的阻滑力為

岸灘崩體的穩(wěn)定系數(shù)定義為崩體阻滑力與下滑力的比值,相應(yīng)的穩(wěn)定系數(shù)為

式中:β為滲透力與水平方向的夾角。將式(2)及

式中:Sαθ為主要反映岸灘土體內(nèi)摩擦角和破壞面傾角影響的參數(shù);Sβα和nJγ為主要反映滲透力影響的參數(shù);Sα為主要反映河岸形態(tài)和破壞面角度影響的參數(shù);St為河岸強(qiáng)度系數(shù),反映河岸土壤黏性與重力的對(duì)比關(guān)系;SA為主要反映河岸崩體大小和形態(tài)影響的參數(shù)。方程(10)為岸灘崩體穩(wěn)定分析模式方程,主要反映了河岸土壤特性、岸坡形態(tài)、河床沖刷、岸灘滲流等因素的共同作用關(guān)系。折線岸灘破壞面的傾角主要取決于河岸形態(tài)、裂隙深度和土壤內(nèi)摩擦角,可由式(11)確定[10]:

2.2 典型岸灘形態(tài)挫落崩塌的分析模式

Thorne等[9]深入研究了洪水期的河岸崩塌問(wèn)題,指出枯水期的河流岸灘崩塌影響因素與洪水期具有明顯的不同。在枯水期間,岸灘水流浸泡和滲流對(duì)河流岸灘崩塌的作用較小,而河流岸灘高度、河床沖刷等則是影響岸灘崩塌的主要因素[11-12]。因此,在探討枯水期河流岸灘崩塌的分析模式時(shí),不需要考慮岸灘水流浸泡和滲流的作用,即f=0。枯水期岸灘崩體穩(wěn)定分析模式方程變?yōu)?/p>

式(12)就是枯水期岸灘崩體穩(wěn)定分析方程。在枯水期,對(duì)于不同岸灘形態(tài)的挫落崩塌,崩體穩(wěn)定分析方程不變,皆為式(12),所變化的是方程中的一些系數(shù)。最典型的岸灘為簡(jiǎn)單邊坡(簡(jiǎn)單邊坡是指河岸邊坡形態(tài)為直線的邊坡),簡(jiǎn)單邊坡岸灘的首次崩塌和二次崩塌有很大的差異。根據(jù)簡(jiǎn)單邊坡岸灘崩塌的特點(diǎn),求得幾種典型岸灘形態(tài)挫落崩塌穩(wěn)定分析方程的主要變化參數(shù),如表1所示。

3 挫落崩塌分析模式的應(yīng)用與檢驗(yàn)

在分析河流岸灘崩塌過(guò)程中,文獻(xiàn)[1,5]提出了岸灘臨界崩塌高度的概念,指出河床沖刷下切過(guò)程中,能夠維持岸灘穩(wěn)定的最大高度稱為岸灘臨界崩塌高度,對(duì)應(yīng)的岸灘崩體穩(wěn)定系數(shù)K=1。當(dāng)河流岸灘高度大于岸灘臨界崩塌高度時(shí),岸灘將發(fā)生挫落崩塌。在黏性岸灘頂面出現(xiàn)張性裂縫的情況下,取穩(wěn)定系數(shù)K=1,求解方程(10)便得岸灘臨界崩塌高度[1,5]:

表1 簡(jiǎn)單邊坡穩(wěn)定分析方程的主要參數(shù)

式(13)表明,岸灘臨界崩塌高度的主要影響因素包括河岸土壤特性、岸坡形態(tài)、河床沖刷、岸灘滲流等。作為算例,文獻(xiàn)[1,5]對(duì)洪水期的臨界崩塌高度的變化過(guò)程進(jìn)行了估計(jì)和分析,結(jié)果表明洪水期和洪水退落末期是河岸容易崩塌的重要時(shí)期。

在枯水期的低水位狀態(tài),根據(jù)方程(12)可求得相應(yīng)的岸灘臨界崩塌高度[1,5]:

枯水期岸灘臨界崩塌高度與河岸強(qiáng)度系數(shù)成正比,與河岸裂隙深度和岸灘坡度成反比,這進(jìn)一步說(shuō)明河流岸灘縱向裂隙會(huì)促進(jìn)岸灘崩塌。

岸灘崩塌機(jī)理十分復(fù)雜,影響因素很多,而開展岸灘崩塌、邊界條件和水動(dòng)力條件的觀測(cè)工作難度較大,因而岸灘崩塌分析模式的驗(yàn)證工作受到很大的限制。鑒于岸灘臨界崩塌高度公式是建立在崩塌模式基礎(chǔ)之上的,且相比之下岸灘臨界崩塌高度在實(shí)驗(yàn)室或天然河流中較容易獲取,因此,可以利用岸灘臨界崩塌高度公式的合理性來(lái)間接檢驗(yàn)河岸崩塌分析模式的適用性。

在天然河道中,針對(duì)河岸土質(zhì)特性、河岸形態(tài)及水流條件等同時(shí)開展觀測(cè)的情況是比較少的,且一般觀測(cè)的河岸高度和坡度都是穩(wěn)定狀態(tài)的,比臨界值小20％以上[5,13]。根據(jù)文獻(xiàn)[14]提供的天然河流崩塌資料和文獻(xiàn)[1,5]的試驗(yàn)資料,進(jìn)行了岸灘崩塌高度的對(duì)比分析。在此需要說(shuō)明的是,文獻(xiàn)[14]提供的Niger河崩塌資料為穩(wěn)定河岸高度與坡度的關(guān)系,沒(méi)有區(qū)分河岸崩塌類型,實(shí)際上也很難區(qū)分,因此該資料僅用于檢驗(yàn);模型沙崩塌高度試驗(yàn)在水槽內(nèi)進(jìn)行,采用清水下切的試驗(yàn)方式,岸灘崩塌多為挫落崩塌。圖4(a)和(b)分別為Niger河穩(wěn)定河岸高度和實(shí)驗(yàn)室模型沙崩塌高度的檢驗(yàn)情況。由圖4可見(jiàn),無(wú)論是Niger河的野外實(shí)測(cè)資料,還是模型沙的崩塌試驗(yàn)研究,岸灘實(shí)測(cè)臨界崩塌高度與計(jì)算值基本一致,表明河流岸灘崩塌分析模式是合理的,可以在實(shí)踐中應(yīng)用。

圖4 岸灘挫落崩塌分析模式的檢驗(yàn)

4 結(jié) 論

a.從河岸土體應(yīng)力分析出發(fā),闡述了黏性河岸挫落崩塌的機(jī)制,指出在河床沖刷過(guò)程中,岸灘頂部會(huì)出現(xiàn)縱向裂隙,縱向裂隙會(huì)加速岸灘挫落崩塌的發(fā)生。

b.岸灘挫落崩塌沿河岸一般呈較小寬度的條形崩塌,崩塌過(guò)程簡(jiǎn)單和短暫,破壞面為平面,崩體沿破壞面進(jìn)入河道的方式主要有倒塌和滑塌兩種。

c.通過(guò)岸灘崩體受力和穩(wěn)定分析,建立了挫落崩塌的分析模式,給出了挫落崩塌的穩(wěn)定分析方程;結(jié)合岸灘臨界崩塌高度的分析,野外和室內(nèi)實(shí)測(cè)資料檢驗(yàn)表明該分析模式是合理的。

[1]王延貴.沖積河流岸灘崩塌機(jī)理的理論分析與試驗(yàn)研究[D].北京:中國(guó)水利水電科學(xué)研究院,2003.

[2]張幸農(nóng),蔣傳豐,陳長(zhǎng)英,等.江河崩岸的影響因素分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2009,39(1):36-40. (ZHANG Xingnong,JIANG Chuanfeng,CHEN Changying, et al.Influencing factors for bank collapse in fluvial rivers [J].Journal of Hohai University:Natural Sciences, 2009,39(1):36-40.(in Chinese))

[3]中華人民共和國(guó)水利部.中國(guó)河流泥沙公報(bào):2010 [M].北京:中國(guó)水利水電出版社,2011.

[4]中華人民共和國(guó)水利部.中國(guó)河流泥沙公報(bào):2011 [M].北京:中國(guó)水利水電出版社,2012.

[5]王延貴,匡尚富.岸灘臨界崩塌高度的研究[J].水利學(xué)報(bào),2007,38(11):1158-1165.(WANG Yangui,KUANG Shangfu.Critical height of bank collapse[J].Journal of Hydraulic Engineering,2007,38(10):1158-1165.(in Chinese))

[6]蔡偉銘,胡中雄.土力學(xué)與基礎(chǔ)工程[M].北京:中國(guó)建筑工業(yè)出版社,1991.

[7]陳希哲.土力學(xué)地基基礎(chǔ)[M].北京:清華大學(xué)出版社, 1984.

[8]ASCE Task Committee on Hydraulic,Bank Mechanics, and Modeling of Riverbank Width Adjustment.River width adjustmentⅠ:processes and mechanisms[J].Journal of Hydraulic Engineering,1998,124(9):881-902.

[9]THORNE C R,TOVEY N K.Stability of composite riverbanks[J].Earth Surface Processes and Landforms, 1981,6:469-484.

[10]JUMIKIS A R.Soil mechanics[M].Princeton:D.Van Nostrand Company,Inc.1962.

[11]王延貴,匡尚富,黃永健,等.洪水期岸灘崩塌有關(guān)問(wèn)題的研究[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào),2003,1 (2):90-97.(WANG Yangui,KUANG Shangfu,HUANG Yongjian.Study on bank collapse during flooding periods [J].Journal of China Institute of Water Resources and Hydropower Research,2003,1(2):90-97.(in Chinese))

[12]王延貴,匡尚富.河岸淘刷及其對(duì)河岸崩塌的影響[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào),2005,3(4):251-257. (WANG Yangui,KUANG Shangfu.Study on bank erosion and its influence on bank collapse[J].Journal of China Institute of Water Resources and Hydropower Research, 2005,3(4):251-257.(in Chinese))

[13]官?gòu)d水庫(kù)坍岸研究小組.水庫(kù)灘岸研究[M].北京:水利電力出版社,1958.

[14]ABAM T K S.Factors affecting distribution of instability of river banks in the Niger Delta[J].Engineering Geology, 1993,35:123-133.

Mechanism and analysis mode of sinking failures of river banks

//WANG Yangui(International Research＆Training Center on Erosion and Sedimentation,China Institute of Water Resources and Hydropower Research,Beijing100048,China)

In order to study the mechanism and analysis mode of sinking failures,the paper analyzes the mechanical mechanism and process of sinking failures by analyzing the soil properties and geologic structures of the river bank and combining the stress distribution and variation characteristics of the river bank soil.It reveals that longitudinal cracks appeared at the top of the viscous bank accelerate the sinking failures during the process of bed degradation,and the process of the sinking failures is relatively simple and short with a plane failure surface.Based on the analysis on the mechanical equilibrium of bank collapses,the analysis mode of sinking failures was studied and a stable analytical equation of sinking failures was derived.The equation shows the interaction among soil properties,slope form,bed scour and bank seepage,and it can be employed to derive the critical height formulas of bank sinking failures.The rationality of the analysis mode was verified by the field measurements.The results suggest that it is rational and can be used in engineering applications.

sinking failure;mechanism of sinking failure;critical failure height;river bank;analysis mode

10.3880/j.issn.10067647.2013.05.005

TV143+.3

A

10067647(2013)05002105

20121129 編輯:熊水斌)

國(guó)家自然科學(xué)基金(51179208);中國(guó)水利水電科學(xué)研究院科研專項(xiàng)(沙集1334)

王延貴(1963—),男,山東鄆城人,教授級(jí)高級(jí)工程師,博士,主要從事水力學(xué)及河流動(dòng)力學(xué)研究。E-mail:wangyg@iwhr.com