孫光林 胡江春 邊亞?wèn)| 陶志剛 梁明星

(1.中原工學(xué)院建筑工程學(xué)院,河南 鄭州 450007；2.深部巖土力學(xué)與地下工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100083；3.煤炭資源與安全開采國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100083)

一般地，對(duì)斷層、空洞等地質(zhì)缺陷的探測(cè)方法主要有：

(1)有損探測(cè)。如鉆探，雖然該方法操作方便，結(jié)果直觀，卻破壞了圍巖的整體穩(wěn)定性，效率低，成本高。

(2)無(wú)損探測(cè)[1-2]。激光探測(cè)法，只能探測(cè)表面現(xiàn)象，不能判斷內(nèi)在情況；電磁法，探測(cè)距離一般在30 m以內(nèi)，精度高，能方便快捷地實(shí)現(xiàn)對(duì)短距離探測(cè)目標(biāo)的預(yù)測(cè)[3]；地震波法，探測(cè)距離一般為50～100 m，精度不如電磁法；地表電法，對(duì)小構(gòu)造探測(cè)效果較差[4]；直流電阻率法，其應(yīng)用仍需進(jìn)一步研究[5]。

地質(zhì)雷達(dá)可確定斷層的位置、走向、傾向、傾角產(chǎn)狀和落差[6]，在以下幾個(gè)方面得到了廣泛應(yīng)用：

(1)隧道。如康富中等[7]在不影響隧道正常運(yùn)轉(zhuǎn)情況下，采用地質(zhì)雷達(dá)對(duì)昆侖山隧道病害狀況進(jìn)行了無(wú)損探測(cè)，為隧道的病害治理提供了可靠依據(jù)。郭亮等[8]將LTD-2100地質(zhì)雷達(dá)用于順層地質(zhì)偏壓隧道圍巖松動(dòng)圈的確定，結(jié)果表明地質(zhì)雷達(dá)在圍巖松動(dòng)圈探測(cè)中是可行的。喻軍等[9]在隧道初內(nèi)襯質(zhì)量檢測(cè)方面做了一定研究工作。楊艷青等[10]應(yīng)用地質(zhì)雷達(dá)對(duì)鐵路隧道整體式襯砌試件進(jìn)行了非接觸檢測(cè)模擬實(shí)驗(yàn)，研究了不同檢測(cè)距離與地質(zhì)雷達(dá)掃描圖像質(zhì)量、精度之間的相互關(guān)系。徐坤等[11]運(yùn)用地質(zhì)雷達(dá)與數(shù)值模擬相結(jié)合的方法對(duì)隧道的安全性進(jìn)行了綜合評(píng)價(jià)。

(2)煤礦井工開采。如宋勁等[12]采用防爆地質(zhì)雷達(dá)對(duì)采煤工作面隱伏鉆桿進(jìn)行了超前探測(cè)定位，在陽(yáng)泉礦區(qū)獲得了較好的應(yīng)用效果。

(3)露天開采。有學(xué)者探索利用地質(zhì)雷達(dá)技術(shù)研究露天礦滑坡體內(nèi)部潛在地質(zhì)缺陷情況，如王永強(qiáng)等[13]運(yùn)用地質(zhì)雷達(dá)對(duì)露天礦山邊坡進(jìn)行了實(shí)地探測(cè)，對(duì)獲得的圖像進(jìn)行了分析解釋，證實(shí)了雷達(dá)探測(cè)結(jié)果與滑坡體地表出露特征一致。另外，地質(zhì)雷達(dá)具有移動(dòng)方便的優(yōu)點(diǎn)，其探測(cè)結(jié)果還能為全站儀、GPS建站選址提供科學(xué)參考。

本文以南芬露天鐵礦334 m平臺(tái)、346 m平臺(tái)和370 m平臺(tái)為研究對(duì)象，運(yùn)用地質(zhì)雷達(dá)無(wú)損探測(cè)技術(shù)方法，探查了該區(qū)域地質(zhì)構(gòu)造及缺陷發(fā)育情況，分析了潛在滑動(dòng)面分布特征，利用數(shù)值模擬方法研究了探測(cè)區(qū)域邊坡穩(wěn)定性情況，為南芬露天礦邊坡災(zāi)害防治提供了可參考的基礎(chǔ)資料。

1 地質(zhì)雷達(dá)探測(cè)機(jī)理

由于不同介質(zhì)介電系數(shù)不同，所以電磁波在不同介質(zhì)中的傳播速度、折射、反射情況也有所不同。

由上式可知，當(dāng)所探測(cè)區(qū)域介質(zhì)界面越多、結(jié)構(gòu)越破碎、組分越復(fù)雜時(shí)，地質(zhì)雷達(dá)所發(fā)射的電磁波被反射和折射的次數(shù)越多，接收天線所接收到的電磁波也越分散，在探測(cè)剖面圖中波形表現(xiàn)得也就越復(fù)雜；反之，當(dāng)所探測(cè)區(qū)域界面越少、結(jié)構(gòu)越單一、組分越簡(jiǎn)單時(shí)，地質(zhì)雷達(dá)所發(fā)射的電磁波被反射和折射的次數(shù)越少，接收天線所接收到的電磁波也越集中，在探測(cè)剖面圖中波形表現(xiàn)得也就越整齊。因此，地質(zhì)雷達(dá)在空洞、斷層、破碎帶等地質(zhì)病態(tài)嚴(yán)重地區(qū)探測(cè)波形尤其復(fù)雜。

2 地質(zhì)雷達(dá)在邊坡探測(cè)中的應(yīng)用

2.1 南芬露天礦探測(cè)區(qū)域情況分析

南芬露天礦礦層頂板圍巖為云母石英片巖和綠泥片巖，上部出露大面積的混合花崗巖。礦層底板為石英綠泥片巖、綠簾角閃片巖，下部為二云母石英片巖和二云母長(zhǎng)石石英片巖。采場(chǎng)下盤334 m平臺(tái)到370 m平臺(tái)區(qū)域深10 m范圍內(nèi)以綠泥角閃巖為主，370 m平臺(tái)屬于老滑體范圍，歷史上曾發(fā)生多次不同規(guī)模的滑坡災(zāi)害，該區(qū)域表面巖體較為松散、破碎，穩(wěn)定性較低。

由于露天礦邊坡各組分介質(zhì)參數(shù)存在差異(表1)，尤其在具有裂隙、滑動(dòng)面特征區(qū)域，地質(zhì)雷達(dá)電磁波反射較強(qiáng)，因此采用地質(zhì)雷達(dá)進(jìn)行探測(cè)是可行的。

表1 露天礦邊坡主要相關(guān)介質(zhì)參數(shù)Table 1 Main medium parameters of open pit mine slope

2.2 地質(zhì)雷達(dá)探測(cè)儀器及探測(cè)方案

根據(jù)現(xiàn)場(chǎng)實(shí)際情況，采用GR-IV型地質(zhì)雷達(dá)，選用200 MHz低頻屏蔽天線，地面耦合方式。時(shí)間窗為220 ns，探測(cè)標(biāo)記點(diǎn)間隔距離為10 m，有效探測(cè)深度為10 m，天線中心通過(guò)標(biāo)記點(diǎn)時(shí)進(jìn)行標(biāo)記定位。露天礦采場(chǎng)下盤邊坡由下至上共布置3條測(cè)線(表2)，平面布置圖如圖1所示，每個(gè)臺(tái)階各輸出1幅剖面圖。

表2 南芬露天礦采場(chǎng)下盤各測(cè)線位置及參數(shù)Table 2 Line location and parametersin the footwall of Nanfen open pit

圖1 測(cè)線布置Fig.1 Line layout diagram

3 探測(cè)結(jié)果及分析

3.1 334 m平臺(tái)探測(cè)結(jié)果及分析

334 m平臺(tái)探測(cè)剖面圖見(jiàn)圖2。圖2表明，地表1 m范圍內(nèi)均出現(xiàn)顯著破碎帶，該破碎帶為修筑簡(jiǎn)易道路時(shí)鋪設(shè)的路基和路面碎石。K0～K5區(qū)間，破碎帶底板埋深從7 m過(guò)渡到6 m，長(zhǎng)約50 m，過(guò)渡界限清晰，呈階梯狀。K4附近出現(xiàn)明顯異常區(qū)，異常區(qū)底板埋深約9 m，寬約3 m，勘查發(fā)現(xiàn)K4標(biāo)記點(diǎn)西側(cè)有明顯裂縫。K5～K14區(qū)間，破碎帶底板界限模糊，出現(xiàn)5處破碎異常區(qū)。K6附近破碎區(qū)高差約5 m，寬約2 m；K8～K9區(qū)間破碎異常區(qū)高差約5 m，寬約3 m；K10～K11區(qū)間，破碎異常區(qū)高差約10 m，寬約3 m；K12～K13區(qū)間，破碎異常區(qū)高差約8 m，寬約4 m?，F(xiàn)場(chǎng)勘查發(fā)現(xiàn)K7標(biāo)記點(diǎn)兩側(cè)有水洼，K10標(biāo)記點(diǎn)附近地表裂縫發(fā)育，K12～K14標(biāo)記點(diǎn)附近發(fā)育多條裂縫。

圖2 334 m平臺(tái)探測(cè)剖面示意Fig.2 334 m platform detection section

因此，K0～K5區(qū)間破碎帶界限較清晰，巖體較完整，自穩(wěn)定性好。K5～K14區(qū)間內(nèi)巖體較為破碎，且地表裂縫發(fā)育，出現(xiàn)明顯裂縫群，基本貫穿整個(gè)平臺(tái)。

3.2 346 m平臺(tái)探測(cè)結(jié)果及分析

346 m平臺(tái)探測(cè)剖面圖見(jiàn)圖3。圖3表明，地表1 m范圍內(nèi)均出現(xiàn)顯著破碎帶，該破碎帶為修筑簡(jiǎn)易道路時(shí)鋪設(shè)的路基和路面碎石。K0～K6區(qū)間，破碎區(qū)呈條帶狀分布，底板埋深9 m，頂板埋深6 m，厚度約3 m，破碎區(qū)長(zhǎng)約60 m，界限清晰，中間無(wú)缺失；K3～K4區(qū)間，出現(xiàn)1個(gè)長(zhǎng)約15 m，厚度約3 m的破碎區(qū)?，F(xiàn)場(chǎng)勘查發(fā)現(xiàn)K4標(biāo)記點(diǎn)西側(cè)有水洼，因此可判斷該區(qū)域受到積水的影響而產(chǎn)生異常。

圖3 346 m平臺(tái)探測(cè)剖面示意Fig.3 346 m platform detection section

因此，346 m平臺(tái)探測(cè)范圍內(nèi)邊坡巖體較為破碎，破碎帶平均厚度約3 m，發(fā)育于地表下6～10 m處，對(duì)邊坡穩(wěn)定性影響較大，一旦該區(qū)域頂板隔水層被雨水或裂隙水貫通，以綠泥角閃巖為主的邊坡巖體強(qiáng)度急劇下降，易引發(fā)局部邊坡滑塌災(zāi)害。

3.3 370 m平臺(tái)探測(cè)結(jié)果及分析

370 m平臺(tái)探測(cè)剖面圖見(jiàn)圖4。圖4表明，地表1.5 m范圍內(nèi)均出現(xiàn)顯著破碎帶，該破碎帶為修筑簡(jiǎn)易道路時(shí)鋪設(shè)的路基和路面碎石。K0～K11區(qū)間，破碎區(qū)呈條帶狀分布，底板埋深9 m，頂板埋深5 m，厚度約4 m，界限清晰，破碎區(qū)長(zhǎng)度約110 m，中間無(wú)缺失。破碎帶上方4 m范圍內(nèi)為完整巖體，巖體內(nèi)基本不存在異常反應(yīng)區(qū)；K11～K13區(qū)間，出現(xiàn)一長(zhǎng)約20 m區(qū)域，該區(qū)域巖體完整，無(wú)任何異常反應(yīng)，該區(qū)域的出現(xiàn)，打斷了破碎帶的延伸，破碎帶在該處尖滅；K13～K16區(qū)間，條帶狀破碎區(qū)繼續(xù)出現(xiàn)，底板埋深9 m，頂板埋深6 m，厚度約3 m，界限清晰，破碎區(qū)長(zhǎng)度約30 m，中間無(wú)缺失，于K16標(biāo)記點(diǎn)處尖滅；K6、K8、K11、K12、K17標(biāo)記點(diǎn)處，出現(xiàn)了垂向條帶狀異常區(qū)。其中，K6處異常區(qū)寬約4 m，高約9 m，現(xiàn)場(chǎng)勘查發(fā)現(xiàn)，K7標(biāo)記點(diǎn)東側(cè)有水洼，由于K6裂隙異常區(qū)的出現(xiàn)，推斷地表水已經(jīng)沿著9 m深的裂隙影響到下部巖體強(qiáng)度；K8處異常區(qū)寬約1 m，高約4 m，地表裂縫不明顯；K11～K12區(qū)間出現(xiàn)多條明顯異常區(qū)，現(xiàn)場(chǎng)勘查可以判斷此異常區(qū)為裂縫異常，多條裂縫組成了裂縫群，裂縫平均深度3 m，寬約1 m，裂縫在地表3 m下尖滅，因此不會(huì)對(duì)深部巖體產(chǎn)生影響；K17標(biāo)記點(diǎn)處出現(xiàn)1條寬約2 m，高約9 m的異常區(qū)，根據(jù)現(xiàn)場(chǎng)記錄可以判斷K16點(diǎn)東側(cè)有明顯寬裂縫，對(duì)照剖面圖可以確定該裂縫已經(jīng)從地表貫穿巖體發(fā)育到9 m深處尖滅。該裂縫對(duì)巖體起到切割破壞作用，地表裂縫兩側(cè)出現(xiàn)明顯落差，有繼續(xù)滑移危險(xiǎn)。

圖4 370 m平臺(tái)探測(cè)剖面示意Fig.4 370 m platform detection section

4 探測(cè)邊坡區(qū)域穩(wěn)定性分析

4.1 邊坡潛在滑面分析

通過(guò)對(duì)南芬露天鐵礦334、346、370 m平臺(tái)的地質(zhì)雷達(dá)物探，獲得了采場(chǎng)中部310～370 m平臺(tái)邊坡破碎帶以條帶狀分布為主。根據(jù)地質(zhì)雷達(dá)物探測(cè)線分布特征和邊坡地形地貌特征，剖分出1條典型剖面線(圖1中1-1′)，按照每個(gè)平臺(tái)相應(yīng)測(cè)點(diǎn)的探測(cè)信息繪制出相應(yīng)剖線的剖面圖，經(jīng)過(guò)剖面圖可以判斷出不同區(qū)域潛在滑動(dòng)面和破碎帶埋深條件和分布特征。剖面1-1′(圖5)按照334 m、346 m和370 m平臺(tái)物探數(shù)據(jù)繪制，其中，334 m平臺(tái)埋深為0～2 m，346 m平臺(tái)埋深為6～10 m，370 m平臺(tái)埋深為6～9 m。

圖5 1-1′ 剖面1-1′滑體及破碎帶示意Fig.5 1-1'section landside mass and fracture zone

同時(shí)，在334、346、370 m平臺(tái)沿剖面1-1′分別布設(shè)B334、B346和B370 3個(gè)位置對(duì)巖土體性質(zhì)進(jìn)行現(xiàn)場(chǎng)鉆探測(cè)試，檢測(cè)結(jié)果見(jiàn)表3。

表3 鉆探檢測(cè)結(jié)果表Table 3 Drilling test results

由表3可知，地質(zhì)雷達(dá)探測(cè)結(jié)果與現(xiàn)場(chǎng)鉆探結(jié)果一致，南芬露天鐵礦采場(chǎng)下盤邊坡潛在滑動(dòng)面埋深較淺，具有淺層滑坡危險(xiǎn)。

4.2 數(shù)值模擬分析

根據(jù)現(xiàn)場(chǎng)工程地質(zhì)條件，采用三維有限差分計(jì)算程序FLAC3D進(jìn)行數(shù)值模擬，邊坡模型尺寸為160 m×200 m×84 m，坡度46°，由4 563個(gè)節(jié)點(diǎn)共16 418個(gè)單元組成。運(yùn)用Mohr-Coulomb屈服準(zhǔn)則，模型底面和側(cè)面邊界位移約束，頂部為自由坡面。模擬得到位移矢量分布見(jiàn)圖6、不平衡力分布見(jiàn)圖7、剪切應(yīng)力分布見(jiàn)圖8。

圖6 位移矢量分布示意Fig.6 Displacement vector distribution map

圖7 不平衡力等值分布示意Fig.7 Unbalanced force equivalent distribution map

圖8 剪切應(yīng)力等值分布示意Fig.8 Equivalent shear stress distribution

由圖6～圖8可得出：

(1)從位移矢量分布圖可以看出，滑體、破碎帶、滑床巖性物理力學(xué)參數(shù)不同，3部分的位移也表現(xiàn)出明顯的差別?；w后緣及前緣位移變形較大，破碎帶后緣位移變形較大，究其原因是由于破碎帶巖體不完整，各力學(xué)參數(shù)降低，造成巖體強(qiáng)度大幅度降低，在自重應(yīng)力作用下，應(yīng)力集中區(qū)首先在破碎帶部位產(chǎn)生，進(jìn)而產(chǎn)生擠壓變形。

(2)從不平衡力分布云圖和剪切應(yīng)力云圖可以看出，剪切破壞區(qū)域主要分布在強(qiáng)度較差的破碎帶后緣，沿坡面向下延伸至破碎帶下緣，表現(xiàn)出一定的非連續(xù)性。不平衡力集中區(qū)域產(chǎn)生于破碎帶與滑體、滑床緊鄰的上下2個(gè)接觸面上。值得注意的是，不平衡力最大值出現(xiàn)在346 m平臺(tái)臨空面處，主要原因是此處巖層受到的約束力較小，加上上部巖體應(yīng)力作用，容易發(fā)生剪切變形破壞。2個(gè)不平衡力集中區(qū)一旦貫通，便構(gòu)成了滑體潛在滑移面。因此，334、346和370 m平臺(tái)區(qū)域穩(wěn)定性受破碎帶作用強(qiáng)烈。

5 結(jié) 論

(1)地質(zhì)雷達(dá)可以探測(cè)淺層滑坡體內(nèi)的地質(zhì)缺陷，本次在南芬露天鐵礦的成功應(yīng)用探明了采場(chǎng)下盤334 m平臺(tái)～370 m平臺(tái)區(qū)域范圍內(nèi)邊坡滑坡體發(fā)育情況，確定了邊坡淺層潛在滑動(dòng)面形態(tài)分布特征。

(2)通過(guò)對(duì)探測(cè)區(qū)域的數(shù)值模擬計(jì)算，分析了破碎帶對(duì)邊坡穩(wěn)定性的不利影響，破碎帶的存在降低了邊坡穩(wěn)定性，且破碎帶極有可能構(gòu)成滑坡體的潛在滑移面。

(3)本文所采用的地質(zhì)雷達(dá)與三維有限差分計(jì)算程序FLAC3D相結(jié)合的綜合分析方法，對(duì)其他類似邊坡工程有一定參考價(jià)值。

[1] 劉盛東，劉 靜，岳建華.中國(guó)礦井物探技術(shù)發(fā)展現(xiàn)狀和關(guān)鍵問(wèn)題[J].煤炭學(xué)報(bào)，2014,39(1):19-20.

Liu Shengdong，Liu Jing，Yue Jianhua.Development status and key problems of Chinese mining geophysical technology[J].Journal of China Coal Society，2014,39(1):19-20.

[2] 程建遠(yuǎn)，石顯新.中國(guó)煤炭物探技術(shù)的現(xiàn)狀與發(fā)展[J].地球物理學(xué)進(jìn)展，2013,28(4):2024-2032.

Cheng Jianyuan，Shi Xianxin.Current status and development of coal geophysical technology in China[J].Progress in Geophys，2013,28(4):2024-2032.

[3] 程久龍，李 飛，彭蘇萍，等.礦井巷道地球物理方法超前探測(cè)研究進(jìn)展與展望[J].煤炭學(xué)報(bào)，2014,39(8):1742-1750.

Cheng Jiulong，Li Fei，Peng Suping，et al.Research progress and development direction on advanced detection in mine roadway working face using geophysical methods[J].Journal of China Coal Society，2014,39(8):1742-1750.

[4] 周建民，金豐年.地下工程塌方的綜合物探研究[J].巖土力學(xué)，2004，25(S2):441-444.

Zhou Jianmin，Jin Fengnian.Study of synthetic geophysical exploration method for collapsed underground engineering[J].Rock and Soil Mechanics，2004，25(S2):441-444.

[5] 張平松,李永盛,胡雄武.巷道掘進(jìn)直流電阻率法超前探測(cè)技術(shù)應(yīng)用探討[J].地下空間與工程學(xué)報(bào)，2013,9(1):135-140.

Zhang Pingsong，Li Yongsheng，Hu Xiongwu.Application and discussion of the advanced detection technology with DC resistivity method in tunnel[J]，Chinese Journal of Underground Space and Engineering，2013,9(1):135-140.

[6] 劉傳孝，蔣金泉，楊永杰.斷層特征參數(shù)地質(zhì)雷達(dá)探測(cè)與識(shí)別的應(yīng)用研究[J].巖土力學(xué)，2004，25(5):832-834.

Liu Chuanxiao，Jiang Jinquan，Yang Yongjie.Research on application of ground penetrating radar to detecting and identifying parameters of faults[J].Rock and Soil Mechanics，2004，25(5):832-834.

[7] 康富中，齊法琳，賀少輝，等.地質(zhì)雷達(dá)在昆侖山隧道病害檢測(cè)中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2010,29(S2):3641-3646.

Kang Fuzhong，Qi Falin，He Shaohui，et al.Application of ground penetrating radar to disease detection of Kunlun Mountain tunnel[J].Chinese Journal of Rock Mechanics and Engineering，2010,29(S2):3641-3646.

[8] 郭 亮，李俊才，張志鋮，等.地質(zhì)雷達(dá)探測(cè)偏壓隧道圍巖松動(dòng)圈的研究與應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2011,30(S1):3009-3015.

Guo Liang，Li Juncai，Zhang Zhicheng，et al.Research on surrounding rock loose zone of tunnel under unsymmetrical loading with ground penetrating radar and its application[J].Chinese Journal of Rock Mechanics and Engineering，2011,30(S1):3009-3015.

[9] 喻 軍，劉松玉，童立元.地質(zhì)雷達(dá)在隧道初襯質(zhì)量檢測(cè)中的應(yīng)用研究[J].巖土力學(xué)，2008，29(11):303-306.

Yu Jun，Liu Songyu，Tong Liyuan.Research on application of GPR to quality detection of first lining of tunnel[J].Rock and Soil Mechanics，2008，29(11):303-306.

[10] 楊艷青，賀少輝，齊法琳，等.鐵路隧道襯砌地質(zhì)雷達(dá)非接觸檢測(cè)模擬試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2011,30(9):1761-1771.

Yang Yanqing，He Shaohui，Qi Falin，et al.Simulation test of GPR non-contact detection on lining of railway tunnel[J].Chinese Journal of Rock Mechanics and Engineering，2011,30(9):1761-1771.

[11] 徐 坤，王志杰，周 藝，等.基于地質(zhì)雷達(dá)技術(shù)隧道安全性評(píng)價(jià)研究[J].地下空間與工程學(xué)報(bào)，2013,9(3):691-696.

Xu Kun，Wang Zhijie，Zhou Yi，et al.Evaluation of tunnel safety based on geological radar technology[J].Chinese Journal of Underground Space and Engineering，2013,9(3):691-696.

[12] 宋 勁，王 磊.探地雷達(dá)探測(cè)采煤工作面隱伏鉆桿研究[J].煤炭學(xué)報(bào)，2014,39(3):537-542.

Song Jin，Wang Lei.Study on the high precision detection of buried drill pipe by ground penetrating radar in the coal face[J].Journal of China Coal Society，2014,39(3):537-542.

[13] 王永強(qiáng)，曹 竹，譚欽文，等.露天礦滑坡體的探地雷達(dá)檢測(cè)技術(shù)[J].煤炭學(xué)報(bào)，2011,36(7):1093-1097.

Wang Yongqiang，Cao Zhu，Tan Qinwen，et al.GPR detection of open-pit landslide[J].Journal of China Coal Society，2011,36(7):1093-1097.