魏 來,唐紅梅
(重慶交通大學巖土工程研究所,重慶400074)
中國是世界上滑坡災害分布最廣、危害最嚴重的國家之一[1]。每年滑坡災害給中國帶來巨大的人員傷亡和經濟損失。巖質滑坡具有突發性、發生規模大、滑速快等特點,例如意大利瓦伊昂水庫庫岸滑坡、云南昭通頭寨溝滑坡、三峽庫區千將坪滑坡,都對人們的生命財產造成極大破壞。國內外一些學者對降雨誘發土質滑坡的啟動機理進行了深入的研究[2-8],但是針對巖質滑坡的啟動機制研究不多,對發生大型巖質滑坡的機理缺乏清晰認識[9-13],特別是暴雨誘發巖質滑坡的水力啟動機制研究較少,專門針對軟弱破裂面形成機理鮮有研究。陳喜昌,等[10]認為巖質滑坡的產生主要跟雨后地下水壓力升高有關;胡其志,等[12]指出降雨由后緣裂縫滲入順層邊坡結構面,產生靜水壓力和下滲力是滑坡產生的關鍵條件;吳永,等[13]則提出巖質邊坡在震后出現裂縫,降雨導致裂縫的擴展直至與軟弱夾層貫通,水降低了軟弱夾層的物理參數,從而引起滑坡。
巖質滑坡的形成須滿足力學失穩條件,在相應水力條件下外除了擁有軟弱結構、還必須具有外傾貫通軟弱破裂面[9],而外傾軟弱破裂面的形成又和巖質邊坡物理、裂縫裂隙、水力條件息息相關。當降雨滲入巖質邊坡裂縫時,將造成結構面處水壓力升高,從而促使節理裂隙的擴展、貫通及泥化。裂隙的擴展又降低了地下水頭的損失,延緩水壓力的衰減,使局部水壓力增加從而又促進軟弱破裂帶的形成,增加地下水揚壓力和下滲力,使滑坡更易發生。
水力劈裂效應是地下水改造巖體結構的一個重要方式,近幾十年隨著大型水利、交通、石油工程的建設,并在斷裂損傷力學的發展下取得了很大的研究成果,其成果已經在工程上廣泛應用[14-17]。筆者以水力學和斷裂力學為基礎,將結構面裂紋的受力模型和巖質邊坡水力模型結合,從力學機制上解釋了巖質邊坡結構面水力劈裂的形成過程,推導了臨界水頭公式,為深入研究巖質滑坡啟動機理提供科學依據。
巖質邊坡結構面水力模型采用Hoek和Bray給出的典型巖石邊坡橫坡面圖,假設破裂滑動面為一平面且滑動面傾角小于坡角,后緣張裂縫垂直并且走向與坡面走向一致,降雨可沿著后緣張裂縫入滲至滑動面,見圖1。

圖1 巖質邊坡結構面水力模型Fig.1 The hydraulic model of structure plane in rock slope
巖質邊坡巖體中有大量的節理裂隙,現先取滑面處據滑坡頂面豎直距離h處的裂紋進行受力分析。假設巖體發生脆性斷裂,應力滿足斷裂力學條件時裂紋將會擴展貫通。將巖體內的三維裂紋概化為二維的平面閉合單裂紋,取豎直方向剖面研究,受力分析可見圖2。

圖2 含裂紋巖體受力分析Fig.2 Stress analysis of cracked rock mass
巖質滑坡一般發生的深度不大,受到的構造應力有限,因此只有自重應力場的作用。因此,圖2中垂直應力σv,水平應力σh計算公式分別為:

式中:N0為巖體的側壓力系數,其值通常小于1;γ為巖體重度;h為巖體厚度。
裂紋長度為2a,和垂直主應力σv作用平面夾角為α,裂紋內作用有孔隙水壓力p,并假設其沿各方向相等。由圖2可知:

由于在斷裂力學中拉為正,所以在式(2)前面加負號。裂紋面上既有正應力σn,又有剪應力σs,因此可以知道此裂紋擴展方式為Ι-ΙΙ復合型。當σn>0時為拉剪復合,σn<0為壓剪復合。
拉剪復合型時,利用工程近似斷裂準則[18],失穩判別式可表示為:

式中:KΙc為裂紋Ι型斷裂的斷裂韌度,KΙ和KΙΙ可以分別用式(5)、式(6)表示:

式中:a為裂紋半長。
將式(2)、式(3)分別代入式(5)、式(6),再將結果代入式(4)可得:

從式(7)中可推導出壓剪裂紋失穩臨界水壓力pc為:

再將式(1)代入式(8)可得:

當σn<0時為壓剪復合型,此時在工程上可以近似認為屬于純ΙΙ型裂紋斷裂問題,KΙ=0。因為裂紋面受壓力閉合,產生摩擦力σntanφ和凝聚力c,因此有效剪應力為:

式中:φ為裂紋面內摩擦角,根據斷裂力學中ΙΙ型裂紋破壞判據可知:

式中:σ's為有效剪應力;KΙΙc為 ΙΙ型斷裂的斷裂韌度。
聯立式(2),式(3),式(10),式(11),整理可得拉剪破壞時臨界水壓力pc:

再將式(1)代入式(12)可得:

若巖體內裂紋不是一條,而是一些列周期分布的等長、等間距的平行共線裂紋(圖3),則可由文獻[18]知斷裂韌度 K'Ιc、K'ΙΙc計算表達式如下:

式中:η為裂紋連通率,η=a/b。

圖3 周期性平行共線裂紋示意Fig.3 Schematic diagram of periodic parallel collinear cracks

分別將式(16),式(17)代入式(9),式(13)可得:
拉剪破壞時臨界水壓力pc':

壓剪破壞時臨界水壓力p'c:

當降雨沿著后緣張裂縫入滲至滑動面,會對周圍產生孔隙水壓力p,假設在某處(見圖1中的h1點)恰好達到臨界水壓力pc,裂縫開裂:

將式(20)分別代入式(18),式(19)可得:
拉剪破壞時臨界水高度h'1:

當等式右邊除α角以外的值都已知,則相同臨界水高度h'1下最易發生水力劈裂效應的裂紋角度:

即 α =22.5°。
壓剪破壞時臨界水高度h'1:

當等式右邊除α角以外的值都已知,則相同臨界水高度h'1下最易發生水力劈裂效應的裂紋角度:


式中:H為坡體總高度;h0為后緣裂縫高度;h1為臨界水高度。若厚度是均布的,則h為常數。水高度大于h1都會發生水力劈裂。
若滑坡體為如圖1,則由幾何關系可知:

將以上數據代入式(21),并取α=22.5°可得:
由此可見,相同裂紋半長a,連通率η越大發生拉剪斷裂所需h'1越大。由于h'1>h0,此邊坡結構面將不會產生拉剪型水力劈裂。
進一步做出不同連通率下裂紋半長和臨界水頭之間的關系,見圖4。

圖4 不同連通率下裂紋半長與臨界水頭之間的關系Fig.4 The relationship between half-length crack and critical water head under different connectivity rates
由圖4可見,相同裂紋半長a,發生壓剪破壞所需水頭h'1隨著η的增加而減少,并隨著η趨近于1而銳減。理論上此邊坡可以發生壓剪型水力劈裂,當η=0.5時,發生水力劈裂的最大裂紋半長為0.000 01 m;當η=0.95時,發生水力劈裂的最大裂紋半長為0.000 073 m。
1)巖質邊坡后緣裂縫在暴雨條件下產生的高水頭結合巖質邊坡結構面中的節理裂隙結構是巖質邊坡產生水力劈裂效應的主要因素。利用文中推導的臨界水頭公式并結合實例證明結構面裂隙在一定水頭條件下將會發生水力劈裂。
2)若發生拉剪型水力劈裂破壞,與最大主應力作用面夾角呈22.5°的裂紋最先發生破壞,若發生壓剪型水力劈裂破壞,與最大主應力面夾角呈45+φ/2的裂紋最先發生破壞。隨著水力劈裂現象的發生,結構面裂隙擴展,貫通。而裂隙的擴展會加強地下水的下滲作用,水頭損失減少,揚壓力進一步加大,將更容易誘發滑坡。
3)巖質邊坡結構面貫通,結構面軟化,巖質滑坡發生之間的機理十分復雜。筆者推導的臨界水頭公式在解決該類問題方面具有一定的參考價值,但是選取參數、公式還不夠嚴謹,有待進一步驗證改進。
[1]王思敬,黃鼎成.中國工程地質世紀成就[M].北京:科學出版社,2004:302.
[2]Beniamino S,Pasquale V.A real time forecasting model for landslides triggered by rainfall[J].Meccanica,1996,31:73-85.
[3]Polemio M,Sdao F.The role of rainfall in the landslide hazard:the case of the Avigliano urban area[J].Engineering Geology,1999,53(3/4):297-309.
[4]Lorella M,Roberto V.Experimental analysis and modeling of shallow landslides[J].Landslides,2007(4):291-296.
[5]Leonardo C,Sabatino C,Manuel P,et al.Modeling of rainfall-induced shallow landslides of the flow-type[J].Geotechnical and Geoenvironmental Engineering,2010,136(1):85-98.
[6]Zonghu L,Yang H,Dalia K,et al.Assessment of shallow landslides from Hurricane Mitch in central America using a physically based model[J/OL].Environment Earth Science,2011[2011-3-9].http://www.springerlink.com/content/6735761014145537/fulltext.html.
[7]李長江,麻土華,李煒,等.滑坡頻度-降雨量的分形關系[J].中國地質災害與防治學報,2010,21(1):87-93.Li Changjiang,Ma Tuhua,Li Wei,et al.Fractal relation of landslide frequency and rainfall[J].The Chinese Journal of Geological Hazard and Control,2010,21(1):87-93.
[8]劉禮領,殷申龍.暴雨型滑坡降水入滲機理分析[J].巖土力學,2008,29(4):1061-1066.Liu Liling,Yin Shenlong.Analysis of rainfall infiltration mechanism of rainstorm landslide[J].Rock and Soil Mechanics,2008,29(4):1061-1066.
[9]王蘭生.意大利瓦依昂水庫滑坡考察[J].中國地質災害與防治學報,2007,18(3):145-148.Wang Lansheng.Vaiont landslide visit in Italy[J].The Chinese Journal of Geological Hazard and Control,2007,18(3):145-148.
[10]陳喜昌,谷明成,石勝偉,等.巖質滑坡的時間預報與水力啟動模型[J].工程地質學報,2009,17(6):750-760.Chen Xichang,Gu Mingcheng,Shi Shengwei,et al.Model of hydrodynamic start-up for prediction of rock landslide[J].Journal of Engineering Geology,2009,17(6):750-760.
[11]Abbas T,Kazuo T.Assessment of the stability of rock slopes by the slope stability rating classification system[J/OL].Rock Mechanics and Rock Engineering,2009[2009-5-5].http://www.springerlink.com/content/w42714tt52515862/fulltext.html.
[12]胡其志,周輝,肖本林,等.水力作用下順層巖質邊坡穩定性分析[J].巖土力學,2010,31(11):3594-3598.Hu Qizhi,Zhou Hui,Xiao Benlin,et al.Analysis of stability of rock bedded slope under hydraulic pressure[J].Rock and Soil Mechanics,2010,31(11):3594-3598.
[13]吳永,何思明,羅渝.震后暴雨型巖質滑坡啟動機理[J].巖土工程學報,2011,33(1):135-139.Wu Yong,He Siming,Luo Yu.Failure mechanism of post-earthquake rock landslide under landslide[J].Chinese Journal of Geotechnical Engineering,2011,33(1):135-139.
[14]黃潤秋,王賢能,陳龍生.深埋隧道涌水過程的水力劈裂作用分析[J].巖石力學與工程學報,2000,19(5):573-576.Huang Runqiu,Wang Xianneng,Chen Longsheng.Hydro-splitting off analysis on underground water in deep-lying tunnels and its effect on water gushing out[J].Chinese Journal of Rock Mechanics and Engineering,2000,19(5):573-576.
[15]朱珍德,胡定.裂隙水壓對巖體強度的影響[J].巖土力學,2000,21(1):64-67.Zhu Dezhen,Hu Ding.The effect of intestitial water pressure on rock mass strength[J].Rock and Soil Mechanics,2000,21(1):64-67.
[16]李宗利,張宏朝,任青文,等.巖石裂紋水力劈裂分析與臨界水壓計算[J].巖土力學,2005,26(8):1216-1220.Li Zhongli,Zhang Hongchao,Ren Qingwen,et al.Analysis of hydraulic fracturing and calculation of critical internal water pressure of rock fracture[J].Rock and Soil Mechanics,2005,26(8):1216-1220.
[17]陳五一,趙顏輝.土石壩心墻水力劈裂計算方法研究[J].巖石力學與工程學報,2008,27(7):1380-1386.Chen Wuyi,Zhao Yanhui.Study of calculation method of hydraulic fracturing for core of earth-rockfill dam[J].Chinese Journal of Rock Mechanics and Engineering,2008,27(7):1380-1386.
[18]洪啟超.工程斷裂力學基礎[M].上海:上海交通大學出版社,1986.
[19]陳洪凱,唐紅梅,王林峰,等.危巖崩塌演化理論及應用[M].北京:科學出版社,2009.