張麗芬,魏貴春,曾夏生,廖武林
清江水布埡庫岸木竹坪滑坡失穩(wěn)探討
張麗芬1,魏貴春1,曾夏生2,廖武林1
(1.中國地震局 地震研究所,武漢 430071;2.湖北中南勘察基礎工程有限公司,武漢 430081)
木竹坪滑坡發(fā)生于2007年5月10日,總方量約800萬m3,是水布埡庫岸的重大滑坡災害之一。根據(jù)實地調(diào)查資料,闡述了滑坡區(qū)基本地質(zhì)條件及滑坡基本地質(zhì)地貌特征。結(jié)合滑坡穩(wěn)定性計算分析了滑坡的成因機制,認為木竹坪滑坡的發(fā)生是內(nèi)外因共同作用的結(jié)果。坡體前緣受清江支流的弧形切割,為滑坡提供了臨空條件,而薄層灰?guī)r及第四系坡積物又為易發(fā)生滑坡的地層巖性條件。在內(nèi)部條件具備的前提下,滑坡發(fā)生前的連續(xù)降雨及水庫蓄水均使得坡體物質(zhì)的力學強度降低,故引發(fā)了大規(guī)模的滑坡失穩(wěn)。
滑坡;地質(zhì)特征;滑坡成因;水布埡水庫蓄水
據(jù)統(tǒng)計,我國已發(fā)現(xiàn)新老滑坡近30萬處,其中災害性滑坡1.5萬處。每年因滑坡災害造成的經(jīng)濟損失在10億元以上[1]。
2007年5月10日凌晨2點,巴東縣清太平鎮(zhèn)水布埡庫岸發(fā)生了大型山體滑坡,滑坡體急劇下滑,將總量約800萬m3的山體和木竹坪村二組16棟民房推入清江支流中,致使清江一級支流磨刀河出口地段近400 m河道被滑坡堆積物堵塞,形成堰塞湖,并激起10 m多高的涌浪,還造成清金線交通中斷,大小船舶停航,直接經(jīng)濟損失達4 000多萬元,所幸未造成人員傷亡。此次滑坡影響范圍較大,給人民的財產(chǎn)安全帶來了較大的危害,引起了社會各界的高度關(guān)注。本文在深入調(diào)查滑坡賦存的地質(zhì)環(huán)境條件、地貌特征和變形破壞跡象的基礎上,對滑坡的成因機制進行了初步的分析與探討,能為以后深入研究類似滑坡成因等提供一定的參考。
木竹坪滑坡體位于湖北省巴東縣清太坪鎮(zhèn)木竹坪村轄區(qū)內(nèi),位于清江北岸,近東西流向清江與南北流向磨刀河交匯處西側(cè)的山岬處,距水布埡大壩約25 km,距清太坪鎮(zhèn)約20 km(如圖1所示),地理坐標東經(jīng)110°12.948′,北緯30°23.993′。滑坡區(qū)內(nèi)夏季濕熱多雨,冬季寒冷而多霧。年平均氣溫14~15℃,相對濕度70%~80%,降水豐沛。
圖1 木竹坪滑坡地理位置示意圖Fig.1 Location of the Muzhuping landslide
木竹坪滑坡體所在的區(qū)域地貌為二、三級剝夷面的交界部位。地貌形態(tài)主要為峰叢槽谷及溶丘洼地。二級剝夷面高程1 250~1 550 m,主要分布在野三關(guān)以南磨刀河兩側(cè)山體部分。三級剝夷面高程950~1 200 m,以溶丘洼地為其特征。區(qū)內(nèi)碳酸鹽巖廣為分布,天坑、溶蝕洼地、溶洞及石芽、溶溝、溶槽巖溶地貌較為發(fā)育。
滑坡區(qū)內(nèi)及附近出露地層主要有二疊系、三疊系及第四系,其中三疊系下統(tǒng)大冶組碳酸鹽巖出露面積最大,巖溶較為強烈;巖性主要為灰白色-深灰色薄層至中厚層石灰?guī)r,層理清晰,夾黃綠、土黃色頁巖,下部為二者互層;二疊系分布于清太坪鎮(zhèn)雙社坪村一帶,主要為龍?zhí)督M灰色頁巖、泥灰?guī)r互層,頂部為炭質(zhì)頁巖,含少量燧石結(jié)核,夾含煤層;第四系零星發(fā)育。
木竹坪滑坡體位于清太坪復向斜內(nèi),該復向斜位于長陽背斜和野三關(guān)背斜之間,走向北東15°,由于長陽背斜的抑制作用,其北段在靠近弧形轉(zhuǎn)彎部位,產(chǎn)生了同步的弧形拐曲。由三疊系地層組成的清太坪向斜、雙社坪背斜、大路坡向斜等,在空間上呈右行雁列分布,兩翼地層傾角30~60°,于轉(zhuǎn)折部位,地層平緩,傾角10°左右,或近于水平。由于常順軸向波狀起伏,在雙社坪背斜的高點上有二疊系出露。雙社坪背斜和清太坪向斜因受到同期縱斷裂的切錯而形態(tài)不甚完整。區(qū)域附近斷裂較少,多為北北東向,如位于滑坡區(qū)西側(cè)的白果園斷裂、桃李溪斷裂以及位于滑坡區(qū)東側(cè)的橋河斷裂等。這些斷裂規(guī)模均較小,且第四紀以來無活動記錄。滑坡體區(qū)域地質(zhì)構(gòu)造見圖2。
圖2 滑坡體區(qū)域地質(zhì)構(gòu)造圖Fig.2 Geological map around Muzhuping landslide
滑坡區(qū)補給區(qū)為河間地塊的各級剝夷面,接受大氣降水補給,由于剝夷面上的峰叢槽谷、溶丘洼地在大面積內(nèi)沒有與地表水系相連的溪溝,故降水連同剝夷面上的溶隙泉,除蒸發(fā)、植物吸收、人畜飲用外,全部轉(zhuǎn)入地下補給地下水,其中降水沿溶隙下滲為本區(qū)的主要形式。
木竹坪滑坡體地表總體形態(tài)呈圈椅狀,上窄下寬。后緣呈弧形,側(cè)緣呈扇形展開,前緣呈弧形,高程230~280 m。前緣臨近清江支流磨刀河河床最寬約600 m,后緣高程800~860 m,寬約200 m?;麦w順主滑方向縱向長約400 m,按平均厚度20 m計算,面積約36萬m2,體積約700萬m3。主滑方向SE100°。木竹坪滑坡全貌如圖3所示。
按照邊坡的物質(zhì)種類劃分,木竹坪滑坡為一土質(zhì)邊坡,上部為黃褐色碎石土偶夾大塊石,土體為黃褐色粉質(zhì)粘土,粘性不好,碎石粒徑1~4 cm,呈次棱角狀-次圓狀。下部基巖為三疊系下統(tǒng)大冶組(T1d)灰色、灰白色灰?guī)r。按滑體各個部位的不同堆積特征,大致可分為2個部分:①后緣,為滑體向下滑動后遺留在后壁的滑帶物質(zhì),以泥土堆積為主,夾小塊石沿斜坡分布約數(shù)米;②中部及前緣有灰?guī)r出露,前緣仍有部分粘土和碎石土。
圖3 木竹坪滑坡全貌Fig.3 General view of Muzhuping landslide
木竹坪滑坡剖面如圖4所示。滑坡周界明顯,滑坡體后緣稍低于山頂,高程約800~860 m;前緣直抵清江支流磨刀河的東岸,高程在230 m左右,滑坡體前后緣高差在600 m左右。滑坡的后壁形成一上陡下緩的滑床光面,坡角約50°,滑床光面上可見多處滑痕,清晰且擦痕方向基本一致,約100°。光面上主要是黃褐色的覆蓋層碎石土,有些部位出露灰?guī)r?;聜?cè)壁與后壁連在一起,呈弧形?;律现饕獨埩粲邢蚪婢弮A的一級臺面,坡度5~15°。
滑坡除了邊界特征明顯外,滑動中往往在滑體內(nèi)形成許多特殊地貌。木竹坪滑坡體上鼻狀鼓丘、馬刀樹和地表裂縫等滑坡地貌十分顯著。由于滑坡體上下各段各塊的滑動時間、滑動速度常常不一致,在滑坡體表面還形成了一些錯臺、陡壁,即滑坡臺階。
圖4 木竹坪滑坡剖面示意圖(滑動前、滑動后)Fig.4 Profile map of Muzhuping landslide(before sliding and after sliding)
該滑坡于2007年5月5日地面開始出現(xiàn)局部裂縫現(xiàn)象,8日裂縫不斷加劇,9日局部開始出現(xiàn)崩塌,10日凌晨便出現(xiàn)了第一次大規(guī)模滑坡,滑坡體變形破壞現(xiàn)象非常明顯。
滑坡體后部變形較強烈,中后部坡體發(fā)生垮塌,滑體后緣產(chǎn)生較大拉裂縫,呈弧形展布。在滑坡體中部亦有鼓脹變形現(xiàn)象。
滑坡體后緣右側(cè),地面明顯下沉形成面積較大的塌陷,深度約30 cm??拷麦w右緣邊界的土坯房已經(jīng)垮塌,附近的磚混房屋破壞嚴重但未垮塌,后墻出現(xiàn)明顯的空心磚錯位約5 cm,墻體及地面出現(xiàn)多處較大的拉裂縫,房屋外墻基處裂縫長約3 m,最寬處約2 cm。正房地面裂縫長約3 m,最寬處約5 cm。土坯房墻體脫落。
滑坡體右側(cè)外部,房屋多處開裂產(chǎn)生破壞,且產(chǎn)生多處較大地裂縫。屋內(nèi)地基下陷,最大沉降處約7 cm。墻體傾斜,墻面脫落,門無法開合。房屋后墻窗口出現(xiàn)明顯裂縫,最大寬約3 cm。在某戶村民房屋的外墻基出現(xiàn)一條長約10 m、寬約5 cm的地裂縫,走向190°。據(jù)居民反映,房屋是在滑坡后5~6 d后才出現(xiàn)的裂縫。墻體原來有輕微裂縫,但也是在滑坡發(fā)生后,房屋裂縫明顯拉大。
在邊坡穩(wěn)定性計算中應用的評價方法較多,針對滑裂面的形態(tài)特征以及滿足的力學條件選用不同的極限平衡方法。針對木竹坪滑坡圓弧狀滑裂面的形態(tài)特征,又考慮到瑞典圓弧法不滿足每一土條力的平衡條件,一般計算出的安全系數(shù)偏低[2],故本文選用簡化畢肖普方法進行穩(wěn)定性計算。
基于剛性極限平衡的簡化Bishop法[2,3]滿足整體力矩平衡條件,滿足各條塊力的多邊形閉合條件,近似考慮了土條間相互作用力的影響,其計算簡圖如圖5所示。
圖5 木竹坪滑坡條分模型及簡化畢肖普法受力示意圖Fig.5 Slice model of Muzhuping landslide and the stressing diagram of Bishop's simplified method
圖中,Ei和Xi分別表示土條間的法向力和切向力,畢肖普假設Xi=Xi+1,即略去土條間的切向力;Wi為土條自重,在浸潤線上、下分別按濕重度和飽和重度計算;Qi為水平力;Ni和Ti分別為土條底部的總法向力(包括有效法向力及孔隙水壓力)和總切向力;βi為第i條塊滑面處傾角;bi為第i土條寬度;xi為第i土條對圓心的力臂長度;R為圓弧半徑;ei為土條水平作用力到圓心的力矩。
選取主滑方向剖面進行穩(wěn)定性計算,對滑坡體中取任意條塊i進行垂直方向的靜力分析,得
根據(jù)滿足土坡穩(wěn)定安全系數(shù)Fs的極限平衡條件,有
式中ci,φi為有效抗剪強度指標。
考慮整個滑動坡體的整體力矩平衡條件,各個土條的作用力對圓心的力矩之和為零,滑動面的正壓力Ni通過圓心,也不產(chǎn)生力矩,只有重力Wi和切向力Ti對圓心產(chǎn)生力矩,即
整理后得
其中,
由于其他參數(shù)均已知,但Fs出現(xiàn)在等式兩邊,故采用試算-迭代法求算即可。
滑坡區(qū)附近的水布埡水庫2007年4月21日二期蓄水,水位從259 m高程上升到290 m,5月上旬該區(qū)持續(xù)降雨,9 d總降雨量超過100 mm,而整個5月份降雨量僅約130 mm??紤]到影響木竹坪滑坡穩(wěn)定性的主要因素,擬采用天然狀況、降雨、蓄水以及組合條件下的4種工況進行滑坡穩(wěn)定性分析。根據(jù)現(xiàn)場調(diào)查,采用地質(zhì)類比和反演分析,選取相應的物理參數(shù)(表1),利用簡化Bishop法進行計算。
表1 木竹坪滑坡穩(wěn)定性計算參數(shù)Table 1 Main parameters in Muzhuping landslide stabilitycalculation
(1)天然狀態(tài),僅考慮自重;
(2)自重+降雨;
(3)自重+290 m水庫蓄水;
(4)自重+降雨+290 m水庫蓄水;
計算結(jié)果顯示,木竹坪滑坡在水布埡水庫未蓄水、滑體前緣未浸泡的天然狀態(tài)下,滑坡體的穩(wěn)定系數(shù)F=1.53,整體穩(wěn)定性很好;在連續(xù)降雨和自重共同影響的情況下,F(xiàn)=1.39,與蓄水前天然狀態(tài)相比,降低了9.2%,但整體穩(wěn)定性仍較好,主要原因是滑體自重遠大于降雨的作用。水布埡水庫從4月21日正式進行二期蓄水,按計劃壩前水位將達347 m,實際上5月10日滑坡滑動時,該處的實際水位為290 m。僅考慮水庫蓄水未降雨時,滑坡的穩(wěn)定系數(shù)F=1.06,比蓄水前降低30.7%;第4種工況則與滑坡滑動時的實際情況吻合,由于庫區(qū)從未發(fā)生過水位從259 m上升到290 m的情況,因此,滑坡對這種工況的抗擾動可能較低,當疊加上降雨時,抗擾動能力更低。計算得出F=0.99,與蓄水前天然狀態(tài)下比較,明顯降低35.3%。
邊坡失穩(wěn)是多種因子耦合的變形,最終被外界誘發(fā)因子激發(fā)的一種地質(zhì)現(xiàn)象。內(nèi)因是變化的基礎,外因是變化的條件,邊坡失穩(wěn)的根本原因取決于組成斜坡的巖土體性質(zhì)和結(jié)構(gòu),但外界因子如水的作用可改變斜坡巖土體的性質(zhì)、強度、結(jié)構(gòu)和狀態(tài),有時還可改變斜坡的外形,成為誘發(fā)邊坡失穩(wěn)的重要因子[4]。本文立足于對滑坡的實地考察,通過對掌握的相關(guān)資料綜合分析,來進行該滑坡成因的初步探討。
4.1.1 地貌條件
有一定高差和坡度的地形地貌是形成滑坡、崩塌所必需的臨空條件[5]。木竹坪滑坡體規(guī)模較大,為一上陡下緩的坡體,后緣陡傾角約50°左右,前緣平臺坡度較緩,約5~15°,前后緣高差約600 m。前緣直抵清江支流磨刀河河床,河流的長期侵蝕與弧形切割,為滑坡提供了臨空條件。
4.1.2 地層巖性組合
實踐證明[6],在我國,凡是下述地層分布的地區(qū),都是容易發(fā)生滑坡的地區(qū):黏性土、黃土、類黃土和各種成因的松散、松軟沉積物;砂巖、頁巖和泥巖的互層地層;煤系地層;灰?guī)r、頁巖或泥灰?guī)r的互層地層等。木竹坪滑坡體表層主要是較厚層的黏土和碎石土等松散坡積物,容易在水的作用下發(fā)生飽水變形等復雜的物理和化學變化,使得土體抗剪強度下降,且下伏滑床基巖為三疊系下統(tǒng)大冶組薄層灰?guī)r,這種特征的巖土結(jié)構(gòu),在外界因素的觸發(fā)下,堆積物易沿下伏灰?guī)r接觸帶失穩(wěn)下滑。由此可見,該滑坡體具備了產(chǎn)生滑坡的基本條件。
影響木竹坪滑坡穩(wěn)定性的因素中,除地層巖性、地形地貌外,水文地質(zhì)條件對滑坡穩(wěn)定性影響也很大。統(tǒng)計發(fā)現(xiàn)[7],90%以上的邊坡失穩(wěn)與水有關(guān),在邊坡失穩(wěn)的發(fā)生過程中,水的作用是一個至關(guān)重要的外在因素。無論是庫水還是降水等其它各種水,其引起斜坡體內(nèi)部的結(jié)構(gòu)及物理化學等因素均是向不利于穩(wěn)定的方向變化的。
4.2.1 水布埡水庫蓄水對坡體的影響
庫岸邊坡的失穩(wěn)與庫水位的變化有著密切的關(guān)系,主要表現(xiàn)在巖土體及結(jié)構(gòu)面的遇水軟化、庫水對巖土體靜水壓力作用、巖土體含水量和容重的增加、水-巖物理化學作用等幾個方面[8]。
據(jù)了解,水布埡水庫于2007年4月21日21時開始下閘進行二期蓄水,截至5月10日,水布埡壩前水位海拔高程由259 m升至290 m。木竹坪滑坡發(fā)生在水布埡水庫二期蓄水開始后20 d,其水位上升與滑坡發(fā)生有密切關(guān)系。庫水對滑坡的作用首先使得涉水區(qū)的巖土體出現(xiàn)局部破壞,木竹坪滑坡體涉水區(qū)處于滑坡的前緣。前緣高程230~280 m,水布埡蓄水至290 m后,滑坡體剪出口及前緣10~60 m的部分浸沒在水位線以下。在水位上升過程中,滑坡前緣土體的力學狀態(tài)發(fā)生改變,在庫水作用下,土體中細顆粒被帶走,土體結(jié)構(gòu)重新調(diào)整。由于庫水的浸泡作用使得土體軟化造成內(nèi)摩擦角降低,同時c值由于土體中膠結(jié)物質(zhì)的減少而降低,從而使得滑坡前緣土體在庫水作用下抗剪強度降低,滑體前緣出現(xiàn)破壞現(xiàn)象。
此外,庫水位的逐漸抬升,使得滑坡地下水位上升,對前緣阻滑段滑體產(chǎn)生浮托作用[9],導致巖土體自重減小,引起邊坡前緣“阻滑段”正壓力的“損失”及抗滑力的下降。破壞了坡體的平衡,當達到臨界值時,滑坡開始大規(guī)?;瑒?。
4.2.2 集中降雨對坡體的影響
近年來,由于氣候條件的變化,集中降雨誘發(fā)了大量的滑坡。盡管降雨誘發(fā)滑坡的力學作用非常復雜,從滑坡誘發(fā)機制上還是可以概括表現(xiàn)為促進滑移面剪應力增大及促使抗剪強度降低。降雨誘發(fā)滑坡,是通過促進滑移面剪應力增大、促使滑坡體抗剪強度降低來實現(xiàn)的[10]。
據(jù)當?shù)鼐用穹从?,在木竹坪滑坡發(fā)生前,該地區(qū)5月初連續(xù)降雨,9 d總降雨量超過100 mm。在降雨入滲前,邊坡土體尤其是上部土體處于非飽和狀態(tài)。在降雨過程中,雨水從2個方面誘發(fā)滑坡:①對邊坡土體,降雨入滲起加載作用,由于河水水位的上升降低了滑坡體內(nèi)的水力梯度,減小了滑坡體內(nèi)地下水的排泄量,使大量降雨滯留于滑坡體中,從而使土體的含水量增大;由于水的作用,滑動面上的有效正應力降低,抗剪強度減小,致使滑坡體的穩(wěn)定性降低。②降雨改變邊坡土體的力學性能,由于該滑坡體表層為黏土、碎石土等坡積物,集中降雨將大量的降水經(jīng)地表滲透轉(zhuǎn)入滑坡體中。雨水滲入滑體,使巖土軟化,導致其內(nèi)聚力下降,基質(zhì)吸力減小,抗剪強度降低。邊坡土體的自重增加和強度降低這兩個不利因素在降雨過程中同時影響邊坡的穩(wěn)定性,達到一定程度就引發(fā)了滑坡。
綜上所述,木竹坪滑坡體具備了滑動的基本物質(zhì)條件和臨空條件,滑坡發(fā)生前的連續(xù)降雨和庫區(qū)蓄水均是滑坡體運動的誘發(fā)因素。
通過對滑坡的現(xiàn)場調(diào)查以及成因的初步討論分析,認為該滑坡體是在降雨與水庫蓄水耦合作用下,前緣失穩(wěn)導致滑動的牽引式滑坡。在水布埡水庫繼續(xù)蓄水至400 m高程后將繼續(xù)淹沒滑坡體,加上該區(qū)降雨充沛,滑面在水的浸泡下軟化,力學性質(zhì)降低,黏聚力及摩擦系數(shù)減小,滑坡體將進一步滑移。建議盡快進行滑坡的治理,安排好當?shù)鼐用竦陌徇w工作。
致謝:感謝姚運生老師和甘家思研究員給論文提出的寶貴修改意見和建議!
[1] 王曉峰,劉光焰,劉均利,等.降雨入滲與滑坡關(guān)系研究綜述[J].人民黃河,2006,28(8):25-27.
[2] 方 云,林 彤,譚松林.土力學[M].武漢:中國地質(zhì)大學出版社,2003.
[3] SL274-2001.碾壓式土石壩設計規(guī)范[S].
[4] 成永剛.邊坡失穩(wěn)的水巖相互作用機理分析[J].探礦工程,2005,(8):38-42.
[5] 楊 林,余 勤,陳福春.南京市滑坡形成機理及誘發(fā)條件研究[J].江蘇地質(zhì),2005,29(4):215-218.
[6] 肖學沛,李天斌,張志龍.三峽水庫區(qū)嚴家坡滑坡成因分析及穩(wěn)定性評價[J].中國地質(zhì)災害防治學報,2005,16(1):39-43.
[7] 喬 娟,羅先啟,張立仁,等.三峽庫區(qū)水環(huán)境對庫岸斜坡失穩(wěn)的作用機理探討[J].水科學與工程技術(shù),2005,(6):22-27.
[8] 朱冬林,任光明,聶德新,等.庫水位變化對水庫滑坡穩(wěn)定性影響的預測[J].水文地質(zhì)工程地質(zhì),2002,(3):6-9.
[9] 張年學.長江三峽工程庫區(qū)順層岸坡研究[M].北京:地震出版社,1993.
[10]黃 濤,羅喜元.地表水入滲環(huán)境下邊坡穩(wěn)定性的模型試驗研究[J].巖石力學與工程學報,2004,23(16):2671-2675.
Failure Mechanism of Muzhuping Landslide in Shuibuya Reservoir Area,Qingjiang River
ZHANG Li-fen1,WEI Gui-chun1,ZENG Xia-sheng2,LIAO Wu-lin1
(1.Institute of Seismology,China Earthquake Administration,Wuhan 430071,China;2.Investigating and Foundation Engineering Co.Ltd in Middle South,Wuhan 430081,China)
Muzhuping landslide,700 million m3in volume,which occurred in the Shuibuya Reservoir area on May 10,2007,is one of the serious landslide disasters in this area after second stage impounding.And it caused a serious loss of property.According to the on-the-spot survey on this landslide,the paper introduces the geological background in landslide area and the geologic and geomorphic characteristics.By means of the landslide stability calculation,it is considered that the occurrence of Muzhuping landslide was the union action result of inner and outer causes.The front area of the landslide was cut by the branch of Qingjiang River,which offered a free face.And the thin-layer limestone and quaternary alluvial purple soil easily inclined to slide.The result shows that the vulnerably geological structure of landslide body,concentrated rainfall and reservoir impoundment caused this landslide.
landslide;geological characteristic;landslide cause;Shuibuya reservoir impoundment
TU457
A
1001-5485(2009)04-0013-05
2008-07-23;
2008-11-10
國家科技支撐項目(2006BAC01B04-02-02);中國地震局地震研究所基金項目(IS200756052)
張麗芬(1981-),女,山東泰安人,研實員,從事水庫誘發(fā)地震研究工作,(電話)027-87166247(電子信箱)zhanglf112@163.com。
(編輯:劉運飛)