張永興,盧 黎,張四平,胡岱文
(重慶大學(xué)土木工程學(xué)院,重慶 400045)
危巖(Overhanging rock)是典型的地質(zhì)災(zāi)害,危巖崩塌具有突發(fā)性強(qiáng),爆發(fā)能量高,沖擊力巨大的特點(diǎn),常造成嚴(yán)重的災(zāi)害。2009年6月重慶武隆雞尾山山體崩塌,造成9人死亡,64人失蹤的慘重事故;2004年貴州納雍山山體崩塌,39人死亡,4人失蹤;2004年9月,萬洲太白巖崩塌80 m3,造成 1人重傷。1948年重慶市洪崖洞危巖崩塌,毀壞房屋幾十間,傷亡300多人。
差異風(fēng)化型危巖是指在軟硬巖互層的沉積巖地區(qū),陡傾的巖質(zhì)邊坡表面風(fēng)化剝落速度不一致,下伏的軟巖退進(jìn)形成巖腔,上部的硬巖向外懸挑,構(gòu)成危巖。西南地區(qū)崩塌災(zāi)害頻發(fā),與差異風(fēng)化型危巖的廣泛分布有直接關(guān)系。以重慶為例,侏羅系中統(tǒng)上沙溪廟組(J2S)沉積巖廣泛分布,泥巖、砂巖互層出露,砂巖抗風(fēng)化能力強(qiáng),泥巖風(fēng)化速度快,泥巖基座風(fēng)化內(nèi)縮,形成巖腔,其上部的砂巖懸挑日益嚴(yán)重,當(dāng)達(dá)到一定的破壞準(zhǔn)則時(shí),砂巖上形成主控裂隙,從而發(fā)育為差異風(fēng)化型危巖。
危巖的分析與治理具有較高的難度,相關(guān)規(guī)范對(duì)危巖的規(guī)定多從治理原則出發(fā),缺乏具體規(guī)定。國內(nèi)的學(xué)者對(duì)危巖的研究也多從工程出發(fā),具體問題具體分析。陳洪凱以萬洲太白崖危巖為例,對(duì)三峽庫區(qū)危巖進(jìn)行了分類[1],胥良以川西108國道高位崩塌為例,研究了崩塌落石的運(yùn)動(dòng)規(guī)律[2],文獻(xiàn)[3]和文獻(xiàn)[4]也是國內(nèi)分析與治理危巖的典型代表。在這些分析中,對(duì)差異風(fēng)化型危巖的形成多用定性描述的方式進(jìn)行,缺乏科學(xué)定量的方法。
差異風(fēng)化型危巖主控裂隙的形成與諸多因素相關(guān)[5-6],邊坡中張性地應(yīng)力的釋放是重要原因之一。邊坡張性地應(yīng)力的存在已被工程地質(zhì)理論分析證實(shí)[7],而張性地應(yīng)力與差異風(fēng)化型危巖形成關(guān)系的研究尚不多見。隨著巖腔的風(fēng)化發(fā)展,根據(jù)危巖體上的受力情況,判斷臨界巖腔深度可以預(yù)測(cè)危巖破壞模式和破壞時(shí)機(jī)。對(duì)張性地應(yīng)力的分析可以判斷危巖的發(fā)育范圍,從而為加固工程提供科學(xué)的依據(jù)。
在巖腔形成過程中,可將風(fēng)化與剝落視為兩個(gè)相互影響的過程。邊坡上的軟巖先風(fēng)化,再因風(fēng)雨重力等原因崩解剝落。地形越陡,降雨沖刷越直接,剝落越快。當(dāng)巖腔形成后,由于懸挑硬巖的遮蔽,退進(jìn)的軟巖不直接受雨水沖刷,剝落速度減慢,風(fēng)化進(jìn)程減緩,而且?guī)r腔越深,這種風(fēng)化進(jìn)程越慢。這種非線性進(jìn)程是差異風(fēng)化型危巖巖腔發(fā)育的重要特點(diǎn)。
根據(jù)實(shí)地觀測(cè)和資料調(diào)研,表1列出了重慶及周邊地區(qū)代表性泥巖邊坡的風(fēng)化速率[8-10]。
表1 泥巖風(fēng)化剝落進(jìn)度
根據(jù)對(duì)巖腔中原住民的調(diào)查和觀測(cè)結(jié)果,以及表1的資料,建議用如下經(jīng)驗(yàn)公式估算巖腔發(fā)育進(jìn)度:
式中:xi為時(shí)間,年;b1為巖腔發(fā)育深度,cm;n為計(jì)算的總時(shí)間周期。
邊坡中的張性地應(yīng)力對(duì)危巖的形成具有重要作用。另外,由于巖腔不斷深入,硬巖日益懸挑,在自重作用下的彎曲拉剪應(yīng)力不斷增高,再與張拉地應(yīng)力疊加,最后引發(fā)破壞,在某一危險(xiǎn)截面上發(fā)育主控裂隙面,從而形成危巖。
圖1 危巖形成力學(xué)模型
按平面問題考慮,主控裂隙形成前的力學(xué)分析如圖1所示,假定主控裂隙沿m-n面形成,則m-n面上的彎矩為:
式中:γ是砂巖重度;β是巖層傾角。E、ks別是水平地震力和水平地震系數(shù);F、frc、kc分別是泥巖支承反力、泥巖抗壓強(qiáng)度,泥巖強(qiáng)度發(fā)揮系數(shù)。
將硬巖視為深梁,考慮到懸臂深梁的應(yīng)力集中,在一般危巖條件下,簡(jiǎn)化算法的彎矩引起的最大彎曲應(yīng)力宜放大兩倍[11],則有:
將彎曲應(yīng)力σfmax與地應(yīng)力引起的拉應(yīng)力疊加,得到最小主應(yīng)力(剪應(yīng)力作用較小,可忽略),最大主應(yīng)力仍為地應(yīng)力引起。顯然b1越大,主應(yīng)力差越大,這種應(yīng)力差達(dá)到一定的破壞準(zhǔn)則后,引發(fā)m-n面破壞,形成主控裂隙。因此,在一定的破壞準(zhǔn)則下可研究b1與其它變量的關(guān)系,見算例。
對(duì)于破壞準(zhǔn)則,結(jié)合應(yīng)用較廣泛的 Hoek-Brown準(zhǔn)則[12]進(jìn)行分析。
式中:σ1、σ3分別是大小主應(yīng)力;σc為巖石單軸抗壓強(qiáng)度;m、s為材料參數(shù)。
在拉剪條件下,存在如下關(guān)系[13]:
式中:τf為巖石純剪切強(qiáng)度。Rt為巖石單軸抗拉強(qiáng)度。
應(yīng)當(dāng)指出,當(dāng)主控裂隙形成后,地應(yīng)力會(huì)充分釋放,主應(yīng)力差迅速減小,系統(tǒng)重新恢復(fù)平衡[14]。因此,張性地應(yīng)力在危巖的形成過程中具有階段性的重要作用。
主控裂隙形成后,在水,風(fēng)化,植物根劈作用下不斷發(fā)展,形成貫通裂隙,此時(shí),危巖上的力主要為重力,軟巖反力,地震力和水壓力,受力模型如圖2所示。
圖2 危巖破壞的力學(xué)模型
對(duì)于貫通式裂隙,由于兩側(cè)裂隙的排水作用,一般降雨條件下,裂隙水高度不會(huì)超過裂隙深度的一半[15]。同時(shí),裂隙尖端處僅會(huì)存在毛細(xì)水,不考慮毛細(xì)壓力的有利作用,按裂隙高度一半計(jì)算水壓力Ew。則危巖塊的破壞模式和相應(yīng)的判據(jù)如下:
(1)軟巖基座破壞
式中,pmax為泥巖基座最大壓力,按式(9)、式(10)計(jì)算。
當(dāng)pmin≥0時(shí)
當(dāng)時(shí)pmin<0時(shí)
(2)平動(dòng)破壞
(3)轉(zhuǎn)動(dòng)破壞
危巖繞o點(diǎn)轉(zhuǎn)動(dòng),臨界狀態(tài)時(shí)危巖體與泥巖基座脫開。則有:
從3種破壞模式可以看出,隨著b1的增大,危巖趨向于破壞,可根據(jù)(8)到(12)式計(jì)算分析破壞的控制模式。
按上述理論對(duì)危巖的形成與破壞進(jìn)行分析需要測(cè)試邊坡地應(yīng)力情況。下面結(jié)合地應(yīng)力測(cè)試進(jìn)行分析和計(jì)算示例。
重慶市渝中區(qū)洪崖洞危巖邊坡位于嘉陵江南側(cè)沖刷岸,解放碑向斜軸部。解放碑向斜為平緩開闊式褶皺,其軸向NE5°~15°,巖層產(chǎn)狀平緩 ,一般為3°~ 10°。
岸坡上陡下緩,中上部為陡坡,坡度 70°~80°,其間分布有多種形式的危巖。下部為30°~50°中緩坡,呈臺(tái)階分布直至嘉陵江[16]。如圖3所示。
圖3 洪崖洞危巖邊坡典型地質(zhì)剖面圖
洪崖洞地區(qū)主要以走向 70°~90°及320°~330°兩組裂隙發(fā)育,裂隙傾角較陡,集中在60°~80°之間。裂隙寬度通常在 0.5~1 cm,最寬的達(dá)到20 cm。
陡崖帶坡體主要由侏羅系中統(tǒng)上沙溪廟組(J2S)砂巖體和其下泥巖基座組成。坡體上部的長石砂巖呈灰白色,中粒結(jié)構(gòu),出露厚度 20.00~25.00 m。坡體中下部的紫紅色泥巖,層厚50~60 m,構(gòu)成危巖軟弱基座。
采用套孔應(yīng)力解除法測(cè)試了洪崖洞危巖邊坡內(nèi)的地應(yīng)力,測(cè)試選定WK7與WK8號(hào)危巖之間的一面完整砂巖巖壁作為測(cè)試地點(diǎn)。從鉆孔取芯看,測(cè)孔僅在4m左右有一道裂隙,更深處的巖石完整,地應(yīng)釋放不充分,是一個(gè)良好的測(cè)試場(chǎng)所。測(cè)試地點(diǎn)距地表深度h=21.53m,測(cè)試鉆孔深度13m,方位角170°,傾角4°,采用KX81型空心包體全應(yīng)力計(jì)和專門設(shè)計(jì)的定向安裝棒定位安裝,見圖2,測(cè)試示意圖見圖3。測(cè)試成果列于表2。
圖4 定向安裝棒示意圖
圖5 地應(yīng)力測(cè)試示意圖
表2 主應(yīng)力大小及方位
砂巖實(shí)測(cè)容重24.2 kN/m3,對(duì)于埋深較淺的地層,地應(yīng)力中間主應(yīng)力往往由重力引起。重力應(yīng)力σz=γ H=24.2×21.53=0.52 MPa。理論計(jì)算結(jié)果與實(shí)測(cè)結(jié)果中間主應(yīng)力非常接近,而中間主應(yīng)力值傾角接近垂直,從而驗(yàn)證了測(cè)試結(jié)果。
將邊坡走向與實(shí)測(cè)的地應(yīng)力大小和方向繪于圖4中。可見主拉應(yīng)力作用面與邊坡走向接近垂直,而已有危巖優(yōu)勢(shì)卸荷裂隙面方向約90°,與主拉應(yīng)力近于垂直。
圖6 測(cè)試結(jié)果與邊坡走向關(guān)系圖
根據(jù)格里菲斯理論,巖石的開裂始于內(nèi)部隨機(jī)存在的細(xì)微裂隙,在一定應(yīng)力條件下,這些微裂隙端部產(chǎn)生明顯的應(yīng)力集中,當(dāng)應(yīng)力集中超過巖石的局部抗拉強(qiáng)度時(shí),新的宏觀裂隙開始出現(xiàn)。格里菲斯強(qiáng)度準(zhǔn)則可作為巖石中新裂隙始發(fā)的判別條件,如公式(13)、(14)所示。
洪崖洞長石砂巖抗拉強(qiáng)度為4.56 MPa,將判據(jù)與地應(yīng)力測(cè)試結(jié)果一起繪于圖8中,可以直觀地看出,測(cè)點(diǎn)應(yīng)力處于非破裂區(qū)中。在測(cè)點(diǎn)部位,單純的地應(yīng)力作用下,尚不具備產(chǎn)生新裂隙的條件。另據(jù)調(diào)查顯示[16],洪崖洞砂巖卸荷裂隙在坡頂向后5~10 m內(nèi)密集發(fā)育。理論分析和現(xiàn)場(chǎng)調(diào)查結(jié)合,可以證明在此范圍以內(nèi)為巖石卸荷松動(dòng)區(qū),是當(dāng)前危巖和潛在次級(jí)危巖的發(fā)育區(qū),而測(cè)點(diǎn)部位表現(xiàn)出相對(duì)穩(wěn)定的狀態(tài),可見危巖發(fā)育范圍約10~12 m左右。
圖7 Griffith破裂判別示意圖
以洪崖洞WK10-5號(hào)危巖為例分析,目前該危巖巖腔深度1.35 m。計(jì)算參數(shù)為:a=4.7 m,b=3.8 m,β=4.5°,砂巖抗壓強(qiáng)度 36.52 MPa,抗拉強(qiáng)度4.56 MPa,粘聚力1.25 MPa,內(nèi)摩擦角36.16°,泥巖抗壓強(qiáng)度7.06 MPa。地應(yīng)力采用表2的實(shí)測(cè)數(shù)據(jù)。當(dāng)b1=0時(shí),m-n上沒有彎矩,可據(jù)此條件反算kc標(biāo)準(zhǔn)值為0.015。圖8列出的計(jì)算結(jié)果表明了破壞時(shí)的巖腔深度與kc的相關(guān)關(guān)系。圖中斜率較大的曲線簇是 σ1-σ3,斜率較小的曲線簇是(下同)。
可見,按標(biāo)準(zhǔn)kc分析,當(dāng)巖腔深度1.07 m時(shí),主控裂隙形成。事實(shí)上,由于泥巖表層風(fēng)化和向臨空面的塑性流變或由于地下水作用,在巖層分界面處出現(xiàn)泥化現(xiàn)象等原因,kc會(huì)逐漸變小,而kc的減小對(duì)主控裂隙形成時(shí)機(jī)影響很大,當(dāng) kc減小到0.01,主控裂隙形成時(shí),巖腔發(fā)育深度僅0.46 m。實(shí)際工程中可根據(jù)具體情況,合理確定kc取值,以預(yù)測(cè)主控裂隙開展時(shí)機(jī)。將巖腔臨界發(fā)育深度代入式(1),即可估算主控裂隙形成時(shí)間。
圖9 ks對(duì)主控裂隙形成的影響
圖9的結(jié)果顯示了水平地震系數(shù)ks對(duì)主控裂隙形成的影響?;鶞?zhǔn)水平地震系數(shù)取為0.05,當(dāng)其上下變化時(shí),曲線簇交點(diǎn)變化不是很大??梢?這個(gè)參數(shù)的變化對(duì)巖腔發(fā)育深度的影響不如kc明顯。
如果考慮到巖體本身存在的裂隙,可將巖石強(qiáng)度參數(shù)降低進(jìn)行類似分析,同時(shí)應(yīng)注意到,裂隙巖體會(huì)引起地應(yīng)力大量釋放,裂隙巖體中殘余地應(yīng)力與完整巖石中存在的地應(yīng)力之比為釋放系數(shù)kr,圖10列出了巖石強(qiáng)度參數(shù)折半和不同kr條件下的計(jì)算結(jié)果??梢?kr的相關(guān)性也是很明顯的。
圖10 kr對(duì)主控裂隙形成的影響
如前所述,當(dāng)主控裂隙形成后,地應(yīng)力會(huì)充分釋放,系統(tǒng)重新恢復(fù)平衡。這一點(diǎn)從張性地應(yīng)力在應(yīng)力差中所占的比例可以明顯看出。因此,張性地應(yīng)力在危巖的發(fā)育過程中只具有階段性作用。
按(8)~(12)式對(duì)WK10-5可能破壞模式進(jìn)行預(yù)測(cè),結(jié)果繪于圖11中,為了直觀比較,將縱軸單位進(jìn)行了統(tǒng)一化處理,不影響結(jié)果。當(dāng)危巖巖腔發(fā)育深度達(dá)到1.83m時(shí),基底壓力線與基底抗力線相交,表明泥巖基座壓壞;巖腔深度接近1.9 m時(shí),發(fā)生轉(zhuǎn)動(dòng)破壞[17];直到巖腔深度達(dá)到2.0時(shí),平動(dòng)推力與抗力仍未取得交點(diǎn),事實(shí)上已超出定義域,這種破壞模式不會(huì)發(fā)生。將控制性破壞模式的巖腔深度代入式(1),即可預(yù)估發(fā)生破壞的時(shí)間。
圖11 典型危巖破壞模式預(yù)測(cè)
差異風(fēng)化型危巖的形成與破壞模式表明,試圖對(duì)危巖用大規(guī)模清除的辦法來治理應(yīng)當(dāng)慎重[18],清除相當(dāng)于進(jìn)行了人工卸荷,去除當(dāng)前危巖后,后面的巖石繼續(xù)發(fā)展為新危巖。
差異風(fēng)化型危巖的形成與破壞是一個(gè)漸進(jìn)的過程,它與地形、地貌,地質(zhì)構(gòu)造,危巖上的各種作用力和巖腔發(fā)育時(shí)間歷程都有關(guān)系。通過對(duì)差異風(fēng)化型危巖體地應(yīng)力測(cè)試和形成破壞機(jī)理的研究,可以得出下面幾個(gè)結(jié)論:
1)分析表明,差異風(fēng)化型危巖的形成與張性地應(yīng)力和巖腔發(fā)育深度相關(guān),而危巖的破壞模式主要受巖腔發(fā)育深度影響。
2)采用空心包體全應(yīng)力計(jì)測(cè)試了危巖邊坡的地應(yīng)力,并驗(yàn)證了測(cè)試結(jié)果的可靠性。將地應(yīng)力測(cè)試結(jié)果與地勘資料綜合分析,可以判定危巖體發(fā)育范圍,確定工程設(shè)計(jì)的穩(wěn)定母巖位置。
3)用地應(yīng)力測(cè)試結(jié)果對(duì)差異風(fēng)化型危巖進(jìn)行計(jì)算分析,結(jié)果表明主控裂隙的形成與地應(yīng)力大小、巖腔發(fā)育深度和泥巖強(qiáng)度發(fā)揮系數(shù)的相關(guān)性較強(qiáng),最后提出了控制性破壞模式的預(yù)測(cè)方法。
[1]陳洪凱,唐紅梅.長江三峽水庫區(qū)危巖分類及宏觀判據(jù)研究[J].中國地質(zhì)災(zāi)害與防治學(xué)報(bào),2005,16(4):53-57.CHEN HONG-KAI,TANG HONG-M EI.Classification and identify of perilous rock in the area of the Three Gorges Reservoir[J].The Chinese Journal of Geological Hazard and Control,2005,16(4):53-57.
[2]胥良,李云貴,劉艷梅.川西108國道高位崩塌成因與運(yùn)動(dòng)特征[J].水文地質(zhì)工程地質(zhì),2008(3):28-31.XU LIANG,LI YUN-GUI,LIU YAN-MEI.Formation mechanism and movementcharacteristic of ahigh altitude collapse on the No108 National Highway in west Sichuan Province[J].Hydrogeology and Engineering Geology,2008(3):28-31.
[3]張立乾,郭富民,楊國興,等.媧皇宮危巖穩(wěn)定性評(píng)價(jià)及其加固對(duì)策[J].建筑結(jié)構(gòu),2007,37(4):68-71.ZHANGLI-QIAN,GUOFU-MIN,YANGGUOXING.Stability evaluation and reinforcement precept to crags of regina's replace[J].Building Structures,2007,37(4):68-71.
[4]黃達(dá),黃潤秋,裴向軍,等.溪洛渡水電站某危巖體穩(wěn)定性及加固措施研究[J].巖土力學(xué),2008,29(5):1425-1429.HUANG DA,HUANG RUN-QIU,PEI XIANG-JUN.Study on stability and reinforcement measures of an unstable rock masses of Xiluodu Hydropower Station[J].Rock and Soil Mechanics,2008,29(5):1425-1429.
[5]張倬元,王士天,王蘭生,等.工程地質(zhì)分析原理[M].2版.北京:地質(zhì)出版社,2002.
[6]李建林.三峽工程巖石拉剪斷裂特性的試驗(yàn)研究[J].地下空間,2002,22(2):149-152.LI JIAN-LIN.A Study on tension-shear crack property of rock related to the three gorgesproject[J].Underground Space,2002,22(2):149-152.
[7]重慶渝碚地質(zhì)工程勘察院.重慶市渝中區(qū)洪崖洞危巖崩塌綜合整治工程勘察報(bào)告[R].2001.
[8]陳洪凱.三峽庫區(qū)危巖鏈?zhǔn)揭?guī)律的地貌學(xué)解譯[J].重慶交通大學(xué)學(xué)報(bào):自然科學(xué)版,2008,27(1):91-95.CHEN HONG-KAI.Geomorphology research on chained regularity of perilous rock in the area of the Three Gorges Reservoir[J].Journal Of Chongqing Jiaotong University:Nature Science Edition,2008,27(1):91-95.
[9]楊宗才,張俊云,周德培.紅層泥巖邊坡快速風(fēng)化特性研究[J].巖石力學(xué)與工程學(xué)報(bào),2006,25(2):275-283.YANG ZONG-CAI,ZHANG JUN-YUN,ZHOU DEPEI.Study on fast weathering characteristics of red bed mudstone slope[J].Chinese Journal of Rock Mechanics and Engineering 2006,25(2):275-283.
[10]胡杰剛,俞敏,全洪波,等.桂柳高速公路邊坡巖石風(fēng)化速度的研究[J].水文地質(zhì)工程地質(zhì),2003,(4):67-71.HUJIE-GANG,YUM IN,QUANHONG-BO.Weathering velocity of slope rock in Guilin-Liuzhou expressway[J].Hydrogeology and Engineering Geology,2003,(4):67-71.
[11]盧黎.重慶洪崖洞危巖邊坡破壞和加固機(jī)理研究[D].重慶:重慶大學(xué),2003.
[12]王桂芳.懸臂深梁的應(yīng)力分析[J].四川聯(lián)合大學(xué)學(xué)報(bào):工程科學(xué)版,1999,3(6):31-38.WANG GUI-FANG.Stress analysis for cantilever deep beam[J].JournalofSichuan Union University:Engineering Science Edition,1999,3(6):31-38.
[13]CHEN HONG-KAI,TANG HONG-MEI.Method to calculate fatigue fracture life ofcontrol fissure in perilous rock[J].Applied Mathematics and Mechanics,2007,28(5):643-649.
[14]OMER M,ABDEL K U.Fracture mechanisms of offset rock joints-A laboratory investigation[J].Geotechnical and Geological Engineering,2004,22:545-562.
[15]GOODMAN R E.Introduction to rock mechanics[M].New York:John Wiley&Sons Inc,1980.
[16]CRAVERO M,IABICHINO G.Analysis of the flexural failure of an overhanging rock slab[J].International Journal of Rock Mechanics and Mining Sciences,2004,41(5):2B 24:1-6
[17]PARONUZZI,PAOLO SERAFINI.Walter stress state analysis of a collapsed overhanging rock slab:A case study[J].Engineering Geology,2009,108(2):65-75.
[18]GIANI GP,GIACOMINI A,M IGLIAZZA M.Segalini a experimental and theoretical studies to improve rock fall analysis and protection work design[J].Rock Mechanics and Rock Engineering,2004,37(5):369-389.