武中鵬，劉 宏，董秀群，鄧凱倫

(1.貴州大學(xué)自然資源部喀斯特環(huán)境與地質(zhì)災(zāi)害重點(diǎn)實(shí)驗(yàn)室，貴州 貴陽(yáng) 550025；2.天津水保工程咨詢(xún)有限公司，天津 300022)

0 引言

據(jù)前人的研究可知，引起崩塌的因素有很多，且它們之間相互作用、相互影響。其中，當(dāng)臨空面和貫通的節(jié)理裂隙存在時(shí)，重力是引起崩塌的首要因素。崩塌危巖體作為不規(guī)則的塊體，在運(yùn)動(dòng)過(guò)程中極易裹挾巖屑，且移動(dòng)距離會(huì)超過(guò)其原始巖坡高度。隨著城鎮(zhèn)化進(jìn)度的不斷加快，地質(zhì)災(zāi)害發(fā)生的頻率也越來(lái)越高，崩塌地質(zhì)災(zāi)害具有突發(fā)性，崩落速度快，極易在短時(shí)間內(nèi)造成巨大的損失，因此，研究單體危巖體崩塌災(zāi)害的危險(xiǎn)性評(píng)價(jià)具有重要的意義。

目前，對(duì)于崩塌的研究還主要局限于崩塌的形成及失穩(wěn)破壞機(jī)理，而對(duì)小尺度的單體崩塌危巖體的危險(xiǎn)性評(píng)價(jià)研究較少。在早期，胡厚田[1]曾對(duì)各種崩塌落石的運(yùn)動(dòng)速度、運(yùn)動(dòng)距離和落石軌跡進(jìn)行了研究；唐紅梅等[2]就崩塌的啟動(dòng)機(jī)制、碰撞過(guò)程中的能量損耗及滾動(dòng)過(guò)程做了深入的研究，總結(jié)出了運(yùn)動(dòng)軌跡方程；張路青等[3]研究了落石對(duì)流動(dòng)人員和車(chē)輛的風(fēng)險(xiǎn)評(píng)價(jià)；鄭光等[4]研究了采礦區(qū)對(duì)崩塌危巖體失穩(wěn)破壞的影響；黨亞倩[5]采用蒙特卡洛法分析風(fēng)頭山崩塌隱患體的穩(wěn)定性狀態(tài)；葉萬(wàn)軍等[6]根據(jù)黃土崩塌常見(jiàn)的運(yùn)動(dòng)形式，提出崩塌體自由飛落及滾動(dòng)情況下影響范圍的理論公式。在前人對(duì)崩塌災(zāi)害危險(xiǎn)性評(píng)價(jià)研究的基礎(chǔ)上，本文以新發(fā)鄉(xiāng)托海小流域樊家?guī)r崩塌體為例，采用失穩(wěn)概率和落實(shí)軌跡模擬，圈畫(huà)出落石影響范圍，初步建立單體危巖體危險(xiǎn)性評(píng)價(jià)方法。

1 樊家?guī)r崩塌基本特征

崩塌前，樊家?guī)r危巖體上部主要為中風(fēng)化灰?guī)r，風(fēng)化裂隙較為發(fā)育，節(jié)理間距0.15～1 m，結(jié)構(gòu)面結(jié)合程度差。危巖體坡向221°，坡角74°，巖層產(chǎn)狀21°∠10°，節(jié)理主要有兩組，分別為J1：185°∠81°、J2：288°∠86°。受差異性風(fēng)化作用的影響，在危巖體下部有凹巖腔發(fā)育，發(fā)育深度為2～3 m，沿節(jié)理面有黏性土、巖屑、巖塊局部充填。巖體在重力及裂隙水或地震作用下局部失穩(wěn)，使巖體發(fā)生墜落。下部巖體(高程1 530～1 820 m)為粉砂巖，坡度在40°～50°。巖體整體較堅(jiān)硬，巖體較破碎，屬Ⅳ類(lèi)巖體。坡腳崩塌堆積體以灰?guī)r為主，沿著危巖體下部斜坡呈扇形堆積，平均坡度為30°，結(jié)構(gòu)松散。區(qū)內(nèi)堆積體塊體普遍較小，一般在0.6 m以?xún)?nèi)。此外，有少量直徑為1.2～2.4 m的塊石在該區(qū)中下部散落。根據(jù)危巖體坡面、巖層、節(jié)理J1和節(jié)理J2的傾向、傾角(表1)繪制赤平投影(圖1)，可判斷危巖體目前處于穩(wěn)定狀態(tài)。危巖體剖面如圖2所示。

表1 赤平投影參數(shù)Table 1 Parameters of stereographic projection

圖1 赤平投影分析圖Fig.1 The analysis diagram of stereographic projection

圖2 樊家?guī)r危巖體剖面Fig.2 The dangerous rock body profile of Fanjiayan

2 危巖體失穩(wěn)概率計(jì)算

根據(jù)所查資料得知，用蒙特卡洛法計(jì)算危巖體的失穩(wěn)概率比較精確。蒙特卡洛法又稱(chēng)隨機(jī)抽樣法、統(tǒng)計(jì)試驗(yàn)法等，其原理是通過(guò)隨機(jī)抽樣和大量的數(shù)據(jù)統(tǒng)計(jì)得到失穩(wěn)概率。本研究根據(jù)蒙特卡洛法原理編寫(xiě)了失穩(wěn)概率分析程序，危巖體的穩(wěn)定性計(jì)算采用重慶市地方標(biāo)準(zhǔn)《地質(zhì)災(zāi)害防治工程勘察規(guī)范》(DB 50/143—2003)中的推薦公式(式(1))。不僅將密度、黏聚力和內(nèi)摩擦角等幾個(gè)常規(guī)的參數(shù)考慮在內(nèi)，還增加了裂隙傾角、裂隙占比、充水占比等因素，以提高結(jié)果的準(zhǔn)確性。在程序中輸入各種參數(shù)的均值和方差，就可以得到危巖體的失穩(wěn)概率值。在具體計(jì)算時(shí)，每次抽樣計(jì)算得到穩(wěn)定性系數(shù)F，對(duì)穩(wěn)定性系數(shù)進(jìn)行判別，抽樣總次數(shù)為1 500次，每次抽樣結(jié)束時(shí)統(tǒng)計(jì)F小于1的次數(shù)，其值與總次數(shù)的比值可以近似為失穩(wěn)概率。

墜落式后緣有陡傾裂隙的懸挑式危巖體穩(wěn)定性計(jì)算公式：

(1)

式中：ζ——危巖抗彎力矩計(jì)算系數(shù)，依據(jù)潛在破壞面形態(tài)取值，一般可取1/12～1/6，當(dāng)潛在破壞面矩形時(shí)取1/6；

flk——危巖體抗拉強(qiáng)度標(biāo)準(zhǔn)值/kPa，根據(jù)巖石抗拉強(qiáng)度標(biāo)準(zhǔn)值乘以0.2的折減系數(shù)確定；

a0——危巖體重心到潛在破壞面水平距離/m；

b0——危巖體重心到過(guò)潛在破壞面形心的鉛錘距離/m；

c——危巖體黏聚力標(biāo)準(zhǔn)值/kPa；

φ——危巖體內(nèi)摩擦角標(biāo)準(zhǔn)值/(°)。

根據(jù)大量的工程實(shí)例和結(jié)合規(guī)范，本文主要考慮了4種不同的工況：天然工況、天然+暴雨(強(qiáng)度重現(xiàn)期按20a考慮)、天然+地震、天然+暴雨+地震(分級(jí)用工況1、工況2、工況3、工況4表示)。由于C和φ值(危巖體的巖石力學(xué)參數(shù)取值如表2所示)受野外采樣、室內(nèi)試驗(yàn)等影響較大，所以作為隨機(jī)變量，得出不同工況下的危巖體失穩(wěn)概率(表3)，同時(shí)得到影響因子的正態(tài)分布圖(圖3、圖4)。

表2 危巖體巖石力學(xué)參數(shù)Table 2 The rock mechanics parameters of dangerous rock mass

表3 危巖體的失穩(wěn)概率Table 3 The instability probability of dangerous rock mass

圖3 自重和地震工況下參數(shù)正態(tài)分布Fig.3 Normal distribution of parameters under gravity and seismic conditions

圖4 暴雨工況下參數(shù)正態(tài)分布Fig.4 Normal distribution of parameters under rainstorm conditions

3 危險(xiǎn)范圍確定

不同形狀的崩塌落石在運(yùn)動(dòng)過(guò)程中具有共性，均以滾動(dòng)和彈跳兩種形態(tài)為主[7]。運(yùn)動(dòng)軌跡受到地形起伏、植被覆蓋度、地層巖性等因素的影響。落石運(yùn)動(dòng)軌跡的模擬在Rockfall原理上進(jìn)行改進(jìn)，將地形起伏、地層巖性、植被覆蓋度和法向、切向恢復(fù)系數(shù)(恢復(fù)系數(shù)取值如表4所示)作為評(píng)價(jià)因子，模擬出落石軌跡(圖5)。將研究區(qū)地形等高線(xiàn)作為模擬落石軌跡的基礎(chǔ)，即避免了人為選取剖面產(chǎn)生的誤差，又可以對(duì)整個(gè)研究區(qū)危巖體進(jìn)行落石軌跡模擬。根據(jù)落石軌跡模擬結(jié)果、地形等高線(xiàn)和居民房屋分布特征，將落石經(jīng)過(guò)的區(qū)域劃分為高危險(xiǎn)區(qū)，臨近落石軌跡的范圍劃分為中危險(xiǎn)區(qū)(圖6)。

表4 法向恢復(fù)系數(shù)和切向恢復(fù)系數(shù)取值Table 4 Values of normal recovery coefficient and tangential recovery coefficient

圖5 樊家?guī)r落石軌跡模擬圖Fig.5 The simulation diagram of Fanjiayan rockfall trajectory

圖6 樊家?guī)r崩塌落石影響范圍劃分Fig.6 The influence scope division of collapse and rockfall in Fanjiayan

4 防治措施建設(shè)

崩塌落石的防治主要包括工程措施、生物措施、規(guī)避、監(jiān)測(cè)預(yù)警和加強(qiáng)公眾教育等。由于樊家?guī)r崩塌危巖體范圍大，處于落石影響區(qū)的居民較多，故提出了有針對(duì)性的防治措施(圖7)，例如InSAR監(jiān)測(cè)、裂隙變形監(jiān)測(cè)、雨量監(jiān)測(cè)、規(guī)劃避難路線(xiàn)等。

圖7 樊家?guī)r崩塌防治措施布設(shè)圖Fig.7 Layout of prevention and control measures for fanjiayan collapse

5 結(jié)論

(1)采用蒙特卡洛法對(duì)危巖體的失穩(wěn)概率進(jìn)行分析，主要考慮了4種工況，結(jié)果表明，危巖體在暴雨和地震的工況下最易失穩(wěn)破壞，其失穩(wěn)概率為0.124。

(2)對(duì)于崩塌落石路徑的模擬，本文在Rockfall軟件基礎(chǔ)上進(jìn)行了改進(jìn)，根據(jù)模擬出來(lái)的落石軌跡確定了落石的影響范圍并對(duì)影響區(qū)進(jìn)行危險(xiǎn)度的劃分，分為高危險(xiǎn)區(qū)、中危險(xiǎn)區(qū)和低危險(xiǎn)區(qū)。

(3)根據(jù)崩塌危巖體的特征、不同程度的危險(xiǎn)區(qū)和地形，對(duì)崩塌危巖體的防治提出了InSAR監(jiān)測(cè)、裂縫變形監(jiān)測(cè)、設(shè)置雨量監(jiān)測(cè)站、放置警示牌和逃避路線(xiàn)指示牌等針對(duì)性的措施。