多節(jié)理石灰?guī)r質(zhì)高陡邊坡崩塌落石分析

馬 俊，郭曉麗，田 野，楊宇，

（1.山西晉焦高速公司，山西 晉城 048000；2.山西省交通科技研發(fā)有限公司 巖土與地下工程重點實驗室，山西 太原 030032）

0 引言

高陡巖質(zhì)邊坡崩塌落石是一種分布廣泛且影響嚴(yán)重的地質(zhì)災(zāi)害，主要通過現(xiàn)場調(diào)查方法進(jìn)行評價與預(yù)測，但是在很多山區(qū)，受地形地貌影響，難以進(jìn)行直接測量，增大了落石災(zāi)害評價的難度。

隨著技術(shù)的不斷發(fā)展與完善，各類非接觸測量技術(shù)被廣泛應(yīng)用于地理信息識別中，三維激光掃描儀、無人機(jī)、衛(wèi)星遙感、Lidar、InSAR、熱成像等，通過這些手段，可以準(zhǔn)確且快速地獲取地面信息，大大降低人為測量的成本和危險。以上測量手段獲取的地理信息多以數(shù)字?jǐn)?shù)據(jù)表達(dá)，最常見的是點云數(shù)據(jù)，點云數(shù)據(jù)是諸多地物點的集合，根據(jù)測量精度不同，點云的密度也有所差別。點云數(shù)據(jù)中，每個點包含X、Y、Z絕對坐標(biāo)，可以確定空間位置。此外還含有其他關(guān)鍵信息，例如地物顏色、地理坐標(biāo)等。根據(jù)點云數(shù)據(jù)，可以在室內(nèi)實現(xiàn)對測量地區(qū)地形地貌的重構(gòu)，并建立高程圖、地形圖、地貌圖等，構(gòu)建三維地質(zhì)模型，為后續(xù)數(shù)值計算提供基礎(chǔ)。

激光掃描結(jié)果有很多種后期處理方法，學(xué)者們基于激光掃描進(jìn)行了深入研究與探討。Fanti[1]利用地面激光掃描技術(shù)提取了巖體結(jié)構(gòu)面幾何結(jié)構(gòu)，將落石類型分為彎曲傾倒和楔塊破壞。Ravanel[2]利用陸地激光掃描識別了地貌，通過長期掃描，可以為環(huán)境變化提供基礎(chǔ)數(shù)據(jù)。Marsella[3]采用地面激光掃描重建了邊坡的幾何結(jié)構(gòu),描述巖石不連續(xù)面的分布和方向，有助于落石危險和風(fēng)險評估。Van Veen[4]以激光掃描為手段，對落石信息進(jìn)行了檢測與識別，構(gòu)建了邊坡變形頻率與震級的關(guān)系。Fanos[5]利用激光雷達(dá)獲取了高精度數(shù)字地形模型，從落石源、誘發(fā)因素、運(yùn)動形式等方面進(jìn)行了全面分析，并對落石進(jìn)行了風(fēng)險評估與預(yù)測。Obanawa[6]利用無人機(jī)運(yùn)動和地面激光掃描測量對復(fù)雜地貌進(jìn)行測量，構(gòu)建了多期三維點云，定量評價了巖體體積變化特征。Vanneschi[7]利用地面激光掃描獲得邊坡表面幾何結(jié)構(gòu)的詳細(xì)三維幾何信息，進(jìn)一步利用三維概念數(shù)值模擬，對摩擦角、不連續(xù)剪切和法向剛度、不連續(xù)間距和方向等關(guān)鍵參數(shù)進(jìn)行了一系列敏感性分析。Dong[8]利用地面激光掃描方法，描述了從原始點云數(shù)據(jù)重建塊狀巖體三維表示的處理方法。塊體離散元是一種以非連續(xù)性介質(zhì)力學(xué)為基礎(chǔ)的數(shù)值模擬方法，適用于塊裂介質(zhì)和板裂介質(zhì)的穩(wěn)定性分析，尤其適用于崩塌落石的計算。

本文利用三維激光掃描技術(shù)，對太行山南麓的高速公路高陡邊坡進(jìn)行了調(diào)查與測量，選取了其中一個典型邊坡進(jìn)行穩(wěn)定性分析，通過建立三維地質(zhì)數(shù)字模型，提出了一套適用于高陡巖質(zhì)邊坡危巖落石的穩(wěn)定性評價方法。從崩塌落石誘發(fā)機(jī)制和運(yùn)動機(jī)理出發(fā)，將巖質(zhì)崩塌分為三大類：摩擦-滑動型、鎖固-斷裂型和碎裂型。研究結(jié)果可以為其他地區(qū)崩塌落石研究提供參考，為崩塌地質(zhì)災(zāi)害防治提供理論基礎(chǔ)。

1 研究區(qū)域與方法

1.1 研究區(qū)域

研究區(qū)域位于山西省太行山南麓，上覆第四系全新統(tǒng)洪積塊石（Q4），下部以奧陶系中下統(tǒng)馬家溝組灰?guī)r（O2）為主，基巖產(chǎn)狀154°～161°∠19°～23°，受太行山構(gòu)造作用影響，巖體節(jié)理裂隙發(fā)育，節(jié)理面無填充物，平均間距約2～4 cm，巖體內(nèi)部有小型溶洞分布。邊坡長度約90 m，高度69 m，整體傾向6°，坡面平均傾角52°，其分布區(qū)域及工程地質(zhì)剖面如圖1所示。

圖1 研究地點（單位：m）

巖塊夾雜有很多細(xì)小石英礦脈，硬度較高，錘擊有清脆響聲，難溶解于水，平均密度為2.42 g/cm3。干燥巖塊平均縱波波速為5 500 m/s，干燥巖塊單軸抗壓強(qiáng)度為50 MPa，飽和單軸抗壓強(qiáng)度為47 MPa，表面回彈值最高8.9 MPa，最低5.1 MPa，主要礦物成分見表1。

表1 礦物成分

1.2 儀器設(shè)備

三維激光掃描技術(shù)是一種新興的非接觸式空間信息獲取技術(shù)，可以對目標(biāo)實現(xiàn)快速、精確、無接觸測量的作用，該次測試采用的是德國Z+F 5010型三維激光掃描儀，可實現(xiàn)水平360°垂直320°全方位點云數(shù)據(jù)掃描與捕捉，最高精度可達(dá)到0.3 mm，最遠(yuǎn)測距187.3 m，可自動檢查數(shù)據(jù)質(zhì)量，檢查標(biāo)志質(zhì)量，查找補(bǔ)測盲區(qū)，數(shù)據(jù)實時拼接。

圖2 Z+F 5010三維激光掃描儀

2 三維激光結(jié)果

三維激光掃描儀掃描結(jié)果如圖3所示，圖中可以清晰地分辨出坡體形狀，由于表層覆蓋有植物，需要進(jìn)行進(jìn)一步處理。

圖3 三維掃描全景

首先對點云進(jìn)行去噪處理（圖4a），去除道路、車輛、植物等干擾因素，再剔除邊坡邊緣的體外孤立點，獲得了較為清晰準(zhǔn)確的點云集合。在此基礎(chǔ)上，將點云通過算法連接，形成面，并通過手動填充孔洞，刪除凸出部分，調(diào)整坡形，使其與實際坡形相同（圖4b）。隨后利用平面工具，對邊坡表面的若干結(jié)構(gòu)面進(jìn)行標(biāo)注，得到平面并編號（圖4c），通過調(diào)取平面信息，可獲取平面的法向向量，進(jìn)而換算得到結(jié)構(gòu)面的產(chǎn)狀。

圖4 三維地質(zhì)結(jié)構(gòu)面建立

以處理后的點云數(shù)據(jù)為基礎(chǔ)，生成等高線，通過插值法進(jìn)行填充，構(gòu)建邊坡三維地質(zhì)模型（圖5），實現(xiàn)模型可視化重構(gòu)?？梢詮膱D中看到明顯的3級臺階，以及白云巖和石灰?guī)r的界面。

圖5 三維地質(zhì)模型建立

3 工程巖體穩(wěn)定性分析

3.1 赤平極射投影分析

依據(jù)三維激光點云數(shù)據(jù)計算得到了12組具有代表性的巖體結(jié)構(gòu)面產(chǎn)狀，列于表2。

表2 巖體主要結(jié)構(gòu)面產(chǎn)狀

依據(jù)節(jié)理特征，繪制出坡面與兩組節(jié)理面的赤平極射投影圖。兩組節(jié)理面的投影交點處于邊坡面外側(cè)，邊坡整體處于穩(wěn)定狀態(tài)，發(fā)生楔形體破壞的可能性較小，但是根據(jù)現(xiàn)場觀測，巖體表層遭受強(qiáng)烈風(fēng)化作用，巖體破碎，容易發(fā)生小塊體崩塌。

圖6 赤平投影圖

3.2 落石能量分布特征

根據(jù)三維激光掃描點云數(shù)據(jù)得到的地形結(jié)果，可得到邊坡剖面圖，按照工程地質(zhì)情況，邊坡分為3種材料類型，頂部為碎石土，中下部為堅硬石灰?guī)r，坡腳為瀝青路面。綜上可知，危巖體以小塊居多，因此在進(jìn)行能量計算時，按照平均重量10 kg為準(zhǔn)，坡體轉(zhuǎn)折部位定位初始點，初始水平速度和豎直速度均為0 m/s，在重力作用下發(fā)生運(yùn)動。

落石整體上發(fā)生了3次運(yùn)動形式的變化，邊坡下部發(fā)生一次，路面上發(fā)生兩次，其中在路面上第一次以彈跳為主，彈跳高度較高，跳躍后能量急劇降低，使得路面上的第二次運(yùn)動形式轉(zhuǎn)為滾動。以邊坡不同位置為基礎(chǔ)，分為7個關(guān)鍵點，進(jìn)行落石能量分析。整體呈現(xiàn)正態(tài)變化，從A到D點，能量逐漸增大，渡過D點后，能量發(fā)生大幅度衰減，在經(jīng)過E-F段時，能量發(fā)生二次衰減，最終散落到G點，能量消耗殆盡。

由此可知，邊坡下部D點位置承受最大落石能量，而落石集中掉落于坡腳的E-F段內(nèi)，在路面發(fā)生跳躍，極易對邊坡下部和坡腳造成破壞。

圖7 落石運(yùn)動軌跡及能量分布

4 危巖落石巖體風(fēng)險分級

基于劉洪亮等[9]提出的可拓學(xué)理論，依據(jù)三維激光掃描獲取的信息，進(jìn)行落石風(fēng)險評價方法，結(jié)合工程巖體質(zhì)量分級標(biāo)準(zhǔn)，充分考慮多種影響因素的相互作用，對危巖落石風(fēng)險評價方法進(jìn)行了改進(jìn)和完善。

根據(jù)研究發(fā)現(xiàn)，影響危巖落石的主要因素可分為三大類：落石運(yùn)動特征、工程地質(zhì)條件以及自然環(huán)境條件，其中具體影響分類列于表3。其中按照對工程影響程度大小分為5級：極低風(fēng)險、低風(fēng)險、中等風(fēng)險、高風(fēng)險以及極高風(fēng)險，每級對應(yīng)的評分等級為1、2、3、4、5分。將落石運(yùn)動特征細(xì)分為7種次級影響因素，基于落石運(yùn)動學(xué)特征、落石能量以及落石影響范圍對影響程度進(jìn)行分級。工程地質(zhì)條件分為8類，依據(jù)工程巖體質(zhì)量進(jìn)行分級。自然環(huán)境條件分為6類。

表3 危巖落石影響因素

由于篇幅限制，計算原理及過程在此暫不列出，可參見劉洪亮等人的計算方法。本文研究的邊坡具體參數(shù)列于表4。最終計算得到的變量特征值為4.16，危巖巖體整體質(zhì)量為Ⅳ級，屬于高風(fēng)險巖體。

表4 太行山南麓某高陡巖質(zhì)邊坡指標(biāo)參數(shù)

對于單個邊坡，不同位置的巖體風(fēng)化程度有所差異，落石運(yùn)動特征也不相同，在進(jìn)行落石防護(hù)時，若選取保守防護(hù)措施，會導(dǎo)致工程造價過高。因此細(xì)化斜坡巖體風(fēng)險分級，可以為防護(hù)措施提供依據(jù)，便于針對性防護(hù)，可以有效節(jié)省成本，達(dá)到最優(yōu)支護(hù)效果。將上述計算方法嵌入程序，基于CAD進(jìn)行二次開發(fā)，針對單體邊坡不同位置的落石特征、巖性差異、風(fēng)化程度進(jìn)一步分級。

圖8 邊坡巖體分級

依據(jù)本節(jié)研究結(jié)果，可以構(gòu)建高陡巖質(zhì)邊坡崩塌落石的力學(xué)模型。針對多節(jié)理硬質(zhì)巖體，根據(jù)節(jié)理發(fā)育程度和接觸類型不同，建立了3種變形破壞模式：摩擦-滑動型、鎖固-斷裂型、碎裂型。具體特征列于表5。

表5 崩塌力學(xué)模式

5 結(jié)語

對太行山南麓灰?guī)r質(zhì)高陡邊坡進(jìn)行調(diào)查，以其中具有代表性的邊坡為對象，通過三維激光掃描，構(gòu)建了邊坡三維地質(zhì)模型，并以此為基礎(chǔ)，利用塊體離散元進(jìn)行了穩(wěn)定性分析，同時建立了落石能量與運(yùn)動路徑的關(guān)系，依據(jù)巖體受力特征和節(jié)理發(fā)育情況對高陡巖質(zhì)邊坡崩塌模式進(jìn)行了分類。得到了以下認(rèn)識：

a）三維激光掃描是一種有效可靠的非接觸測量方法，可以準(zhǔn)確地識別地物信息，并進(jìn)行一系列后續(xù)處理與分析，可極大縮減勞動量，提高工作效率。

b）高陡邊坡崩塌主要受到層理和節(jié)理控制，層理、節(jié)理與坡形三者的空間關(guān)系共同決定了邊坡的穩(wěn)定性。

c）節(jié)理發(fā)育的碳酸鹽巖邊坡在發(fā)生崩塌落石后，還需注意其引起的次生災(zāi)害，重點探測坡體內(nèi)部的溶洞。

d）落石運(yùn)動路徑受坡形影響較大，落石能量與運(yùn)動路徑呈正態(tài)分布的關(guān)系，當(dāng)落石運(yùn)動至邊坡中下部時，能量達(dá)到最大，在經(jīng)歷一系列彈跳、滾動等行為后，能量迅速衰減，大量落石集中在坡腳，形成塊石堆積體。

e）按照巖體結(jié)構(gòu)面發(fā)育特征和受力特征，建立了3種高陡巖質(zhì)邊坡崩塌力學(xué)模式：摩擦-滑動型、鎖固-斷裂型以及碎裂型，不同模式具有不同的運(yùn)動形式和破壞程度。