魏 來，唐紅梅

(重慶交通大學巖土工程研究所，重慶400074)

中國是世界上滑坡災(zāi)害分布最廣、危害最嚴重的國家之一［1］。每年滑坡災(zāi)害給中國帶來巨大的人員傷亡和經(jīng)濟損失。巖質(zhì)滑坡具有突發(fā)性、發(fā)生規(guī)模大、滑速快等特點，例如意大利瓦伊昂水庫庫岸滑坡、云南昭通頭寨溝滑坡、三峽庫區(qū)千將坪滑坡，都對人們的生命財產(chǎn)造成極大破壞。國內(nèi)外一些學者對降雨誘發(fā)土質(zhì)滑坡的啟動機理進行了深入的研究［2-8］，但是針對巖質(zhì)滑坡的啟動機制研究不多，對發(fā)生大型巖質(zhì)滑坡的機理缺乏清晰認識［9-13］，特別是暴雨誘發(fā)巖質(zhì)滑坡的水力啟動機制研究較少，專門針對軟弱破裂面形成機理鮮有研究。陳喜昌，等［10］認為巖質(zhì)滑坡的產(chǎn)生主要跟雨后地下水壓力升高有關(guān);胡其志，等［12］指出降雨由后緣裂縫滲入順層邊坡結(jié)構(gòu)面，產(chǎn)生靜水壓力和下滲力是滑坡產(chǎn)生的關(guān)鍵條件;吳永，等［13］則提出巖質(zhì)邊坡在震后出現(xiàn)裂縫，降雨導致裂縫的擴展直至與軟弱夾層貫通，水降低了軟弱夾層的物理參數(shù)，從而引起滑坡。

巖質(zhì)滑坡的形成須滿足力學失穩(wěn)條件，在相應(yīng)水力條件下外除了擁有軟弱結(jié)構(gòu)、還必須具有外傾貫通軟弱破裂面［9］，而外傾軟弱破裂面的形成又和巖質(zhì)邊坡物理、裂縫裂隙、水力條件息息相關(guān)。當降雨滲入巖質(zhì)邊坡裂縫時，將造成結(jié)構(gòu)面處水壓力升高，從而促使節(jié)理裂隙的擴展、貫通及泥化。裂隙的擴展又降低了地下水頭的損失，延緩水壓力的衰減，使局部水壓力增加從而又促進軟弱破裂帶的形成，增加地下水揚壓力和下滲力，使滑坡更易發(fā)生。

水力劈裂效應(yīng)是地下水改造巖體結(jié)構(gòu)的一個重要方式，近幾十年隨著大型水利、交通、石油工程的建設(shè)，并在斷裂損傷力學的發(fā)展下取得了很大的研究成果，其成果已經(jīng)在工程上廣泛應(yīng)用［14-17］。筆者以水力學和斷裂力學為基礎(chǔ)，將結(jié)構(gòu)面裂紋的受力模型和巖質(zhì)邊坡水力模型結(jié)合，從力學機制上解釋了巖質(zhì)邊坡結(jié)構(gòu)面水力劈裂的形成過程，推導了臨界水頭公式，為深入研究巖質(zhì)滑坡啟動機理提供科學依據(jù)。

1 巖質(zhì)邊坡結(jié)構(gòu)面水力劈裂分析

巖質(zhì)邊坡結(jié)構(gòu)面水力模型采用Hoek和Bray給出的典型巖石邊坡橫坡面圖，假設(shè)破裂滑動面為一平面且滑動面傾角小于坡角，后緣張裂縫垂直并且走向與坡面走向一致，降雨可沿著后緣張裂縫入滲至滑動面，見圖1。

圖1 巖質(zhì)邊坡結(jié)構(gòu)面水力模型Fig.1 The hydraulic model of structure plane in rock slope

巖質(zhì)邊坡巖體中有大量的節(jié)理裂隙，現(xiàn)先取滑面處據(jù)滑坡頂面豎直距離h處的裂紋進行受力分析。假設(shè)巖體發(fā)生脆性斷裂，應(yīng)力滿足斷裂力學條件時裂紋將會擴展貫通。將巖體內(nèi)的三維裂紋概化為二維的平面閉合單裂紋，取豎直方向剖面研究，受力分析可見圖2。

圖2 含裂紋巖體受力分析Fig.2 Stress analysis of cracked rock mass

巖質(zhì)滑坡一般發(fā)生的深度不大，受到的構(gòu)造應(yīng)力有限，因此只有自重應(yīng)力場的作用。因此，圖2中垂直應(yīng)力σv，水平應(yīng)力σh計算公式分別為:

式中:N0為巖體的側(cè)壓力系數(shù)，其值通常小于1;γ為巖體重度;h為巖體厚度。

裂紋長度為2a，和垂直主應(yīng)力σv作用平面夾角為α，裂紋內(nèi)作用有孔隙水壓力p，并假設(shè)其沿各方向相等。由圖2可知:

由于在斷裂力學中拉為正，所以在式(2)前面加負號。裂紋面上既有正應(yīng)力σn，又有剪應(yīng)力σs，因此可以知道此裂紋擴展方式為Ι-ΙΙ復(fù)合型。當σn＞0時為拉剪復(fù)合，σn＜0為壓剪復(fù)合。

拉剪復(fù)合型時，利用工程近似斷裂準則［18］，失穩(wěn)判別式可表示為:

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

式中:a為裂紋半長。

將式(2)、式(3)分別代入式(5)、式(6)，再將結(jié)果代入式(4)可得:

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

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

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

式中:φ為裂紋面內(nèi)摩擦角，根據(jù)斷裂力學中ΙΙ型裂紋破壞判據(jù)可知:

式中:σ's為有效剪應(yīng)力;KΙΙc為 ΙΙ型斷裂的斷裂韌度。

聯(lián)立式(2)，式(3)，式(10)，式(11)，整理可得拉剪破壞時臨界水壓力pc:

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

若巖體內(nèi)裂紋不是一條，而是一些列周期分布的等長、等間距的平行共線裂紋(圖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:

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

將式(20)分別代入式(18)，式(19)可得:

拉剪破壞時臨界水高度h'1:

當?shù)仁接疫叧两且酝獾闹刀家阎?，則相同臨界水高度h'1下最易發(fā)生水力劈裂效應(yīng)的裂紋角度:

即 α =22.5°。

壓剪破壞時臨界水高度h'1:

當?shù)仁接疫叧两且酝獾闹刀家阎?，則相同臨界水高度h'1下最易發(fā)生水力劈裂效應(yīng)的裂紋角度:

式中:H為坡體總高度;h0為后緣裂縫高度;h1為臨界水高度。若厚度是均布的，則h為常數(shù)。水高度大于h1都會發(fā)生水力劈裂。

若滑坡體為如圖1，則由幾何關(guān)系可知:

2 實例分析

將以上數(shù)據(jù)代入式(21)，并取α=22.5°可得:

由此可見，相同裂紋半長a，連通率η越大發(fā)生拉剪斷裂所需h'1越大。由于h'1＞h0，此邊坡結(jié)構(gòu)面將不會產(chǎn)生拉剪型水力劈裂。

進一步做出不同連通率下裂紋半長和臨界水頭之間的關(guān)系，見圖4。

圖4 不同連通率下裂紋半長與臨界水頭之間的關(guān)系Fig.4 The relationship between half-length crack and critical water head under different connectivity rates

由圖4可見，相同裂紋半長a，發(fā)生壓剪破壞所需水頭h'1隨著η的增加而減少，并隨著η趨近于1而銳減。理論上此邊坡可以發(fā)生壓剪型水力劈裂，當η=0.5時，發(fā)生水力劈裂的最大裂紋半長為0.000 01 m;當η=0.95時，發(fā)生水力劈裂的最大裂紋半長為0.000 073 m。

3 結(jié)論

1)巖質(zhì)邊坡后緣裂縫在暴雨條件下產(chǎn)生的高水頭結(jié)合巖質(zhì)邊坡結(jié)構(gòu)面中的節(jié)理裂隙結(jié)構(gòu)是巖質(zhì)邊坡產(chǎn)生水力劈裂效應(yīng)的主要因素。利用文中推導的臨界水頭公式并結(jié)合實例證明結(jié)構(gòu)面裂隙在一定水頭條件下將會發(fā)生水力劈裂。

2)若發(fā)生拉剪型水力劈裂破壞，與最大主應(yīng)力作用面夾角呈22.5°的裂紋最先發(fā)生破壞，若發(fā)生壓剪型水力劈裂破壞，與最大主應(yīng)力面夾角呈45+φ/2的裂紋最先發(fā)生破壞。隨著水力劈裂現(xiàn)象的發(fā)生，結(jié)構(gòu)面裂隙擴展，貫通。而裂隙的擴展會加強地下水的下滲作用，水頭損失減少，揚壓力進一步加大，將更容易誘發(fā)滑坡。

3)巖質(zhì)邊坡結(jié)構(gòu)面貫通，結(jié)構(gòu)面軟化，巖質(zhì)滑坡發(fā)生之間的機理十分復(fù)雜。筆者推導的臨界水頭公式在解決該類問題方面具有一定的參考價值，但是選取參數(shù)、公式還不夠嚴謹，有待進一步驗證改進。

［1］王思敬，黃鼎成.中國工程地質(zhì)世紀成就［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］李長江，麻土華，李煒，等.滑坡頻度-降雨量的分形關(guān)系［J］.中國地質(zhì)災(zāi)害與防治學報，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］劉禮領(lǐng)，殷申龍.暴雨型滑坡降水入滲機理分析［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］.中國地質(zhì)災(zāi)害與防治學報，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］陳喜昌，谷明成，石勝偉，等.巖質(zhì)滑坡的時間預(yù)報與水力啟動模型［J］.工程地質(zhì)學報，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］胡其志，周輝，肖本林，等.水力作用下順層巖質(zhì)邊坡穩(wěn)定性分析［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］吳永，何思明，羅渝.震后暴雨型巖質(zhì)滑坡啟動機理［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］洪啟超.工程斷裂力學基礎(chǔ)［M］.上海:上海交通大學出版社，1986.

［19］陳洪凱，唐紅梅，王林峰，等.危巖崩塌演化理論及應(yīng)用［M］.北京:科學出版社，2009.