地下轉露天礦山復雜采空區探測工程實踐

譚卓英 鐘 文 胡天壽 馮 蕭

(1.北京科技大學土木與環境工程學院，北京100083;2.金屬礦山高效開采與安全教育部重點實驗室，北京100083; 3.江西銅業股份有限公司永平銅礦，江西上饒334506;4.江西理工大學工程研究院，江西贛州341000)

露天開采相比于地下開采具有許多明顯的優勢。例如，礦石開采能力大、礦石回收率高、工作安全系數高、施工作業環境好等［1］。近年來，隨著采礦生產力水平的不斷提高，國內外大部分具備相應條件的礦山都實施了地下轉露天開采。盡管該類礦山總體上數量不多，但在礦石儲量、開采條件等因素允許的情況下盡早轉變為露天開采對于促進當今地下開采礦山企業的高能效、低能耗、可持續性發展來說具有積極的戰略意義［2-4］。然而在地下轉露天開采礦山的區域內必然存在著大量形態各異、成因復雜的采空區，這些采空區由于若干年的荒廢變得更加難以標定和探測，時刻威脅露天礦山的正常生產，給露天開采人員及設備帶來極大的安全隱患。因此，在眾多采空區探測方法中找到一種更加適用于地下轉露天礦山復雜采空區探測的技術方法具有重要的現實和長遠意義。本文主要介紹利用探地雷達(GPR)對地下轉露天開采境界內的淺地表復雜采空區進行探測的技術方法及其地質效果。

1 采空區的分類及其探測技術

1.1 采空區的分類

在礦山領域，采空區分類形式有很多，從采空區的成因和演化的角度出發可以將其分為4類［5-6］。

(1)規劃開采形成的采空區。該類采空區主要是指新中國成立以后，正規礦山企業合理規劃開采所形成的采空區，其特點是有較為詳細的采場、巷道及相關輔助設施的設計圖紙資料，采空區的位置和形態可以很容易地在圖紙上標定出來。

(2)古代開采和非法盜采形成的采空區。從已發掘的古代采礦遺址來看，我國的采礦歷史可以追溯到宋代甚至更早，大量的古代開采活動在推動生產力水平前進的同時也遺留下許多未知的采空區。非法盜采則是從古至今都無法杜絕的民間采礦形式，粗放式的采礦工藝不僅對有限的礦產資源造成了極大的浪費，其形成的采空區也為今后的正規開采埋下了安全隱患。這2種采空區俗稱老窿，其特點是基本沒有標定采空區所需參考資料，有必要對其進行采空區探測與處理。

(3)充填不充分或充填體沉降形成的采空區。在對采空區進行充填處理時，無論是何種充填材料，都難以做到充分充填采空區，再加上充填體長時間的沉降，頂板和圍巖在外部人工或自然因素的干擾下發生2次破碎，使充填不充分區域擴大化、復雜化，形成采空區。該類采空區的特點是規模一般較小，空區內往往夾雜2次破碎充填物，通過現有資料能夠判斷大概的位置和形態。

(4)2次垮塌形成的采空區。在早期地下開采過程中存在這一些未經處理、規模較大的采空區，這些采空區原本可以標定，但由于長期受到降雨入滲、爆破震動、人工擾動和應力釋放等外部因素的影響發生了2次垮塌，原采空區周邊的巖體產生了大規模的破壞，相當于形成一個體積更大、形態更為復雜的采空區。該類采空區雖然不多見，但其破壞力極強，應當進行適當的排查。

通過對以上4類采空區類型和特點進行歸納分析可知，規劃開采形成的采空區易于標定和控制，不會對于地下轉露天礦山的日常生產造成太大影響;古代開采和非法盜采形成的采空區、充填不充分或充填體沉降形成的采空區以及2次垮塌形成的采空區，由于形態復雜且位置不明，對安全生產構成威脅，需要進行合理有效的采空區探測工作。

1.2 采空區探測方法

目前，用于探測地下采空區的技術方法日新月異，但歸納起來還是隸屬于鉆探方法和地球物理方法這兩個大類。鉆探方法的特點是形象直觀，但容易造成以孔代區、以偏概全，尤其是對于古代開采、民間非法盜采且廢棄多年的復雜采空區探測精度難以掌控。地球物理方法最大的特點就是一次性探測區域大，所以在空區探測領域發展迅速。其中，具有代表性的地球物理方法包括高密度電法、地震類方法、GPR等。高密度電法對水較為敏感，所以在判斷采空區含水量時優勢明顯，但總體探測精度不高;地震類方法主要應用于埋深較大的采空區，對于淺部及中深部采空區探測有一定的局限性;GPR法則適用于埋深較淺的采空區，其具有操作方便簡單、圖像清晰直觀、解譯通俗易懂、成本經濟高效等優點，非常適合地下轉露天礦山復雜采空區的探測普查工程［3］。

2 GPR法的基本原理

GPR的探測方式是采用單個固定形式的反射/接收天線組合，通過在地表移動該天線組合實現對地下區域的反射或透射探測。其工作原理可以簡單概括為發射器將高頻電磁波轉化為一系列寬頻帶短脈沖，通過發射天線T穿透地表射入地下，在透射過地下各種巖土體介質或采空區后再反射回地表被接收天線R接收，形成閉合回路。其工作原理如圖1所示。探地雷達的脈沖波雙程走時可以表示為

式中，t為高頻電磁波的雙程走時;z為反射體與探測地表之間的距離;v為地下巖土體介質中的波速;x為發射與接收天線之間的距離［4］。

圖1 GPR基本原理示意Fig.1 GPR detection p rincip les

如果已知地下巖土體介質中的波速v，則可以將探測過程記錄的精確時間t代入(1)式求出反射體與探測地表之間的距離［5］。高頻電磁波在地下巖土體介質中傳播時，其強度和波形將會因為透射介質不同的電性質及幾何形態而呈現出特有的形態［7］。因此，根據雙程走時t、波幅以及波形等數據資料，就可以推斷出探測區域地下巖土體介質的大致結構(包括采空區的基本位置、形態)。通常輸出的雷達數據圖是以脈沖反射波波形的方式保存，以黑、白2種顏色分別表示波形的正、負峰，也可以通過灰階或彩色來展現，最終在軸向剖面以波形圖等灰度圖或等色線圖的形式即可較為直觀地表現出地下各種不同的巖土體介質反射面［8-9］。GPR探測采空區基本原理如圖2所示。

圖2 GPR探測采空區基本原理示意Fig.2 The p rincip les of goaf detection GPR

3 工程實例

銀山九區銅硫金露天礦是國內較為有代表性的地下轉露天礦山。在露天開采境界內的邊坡下賦存一定數量的未充填或廢石充填不充分的規劃性采空區，包括有標高+72m中段采空區、+22m中段采空區、－33 m中段采空區及－88 m中段采空區等，其中可能有個別未充填空區已經演化發展為大型2次垮塌采空區。另外在露天采場北東邊坡還存在少量位置不詳的古代開采和民間非法盜采所形成的采空區(老窿)。以上這些采空區經過長時間的演化，具有復雜的形態特征，對安全生產構成了極大威脅，需要對其開展適當的普查和詳查工作。

3.1 GPR現場探測技術方法

本次復雜采空區探測工程所采用的儀器為IDS公司生產的RIS－K2型探地雷達，其性能指標如表1所示。根據露天礦區已掌握的歷史資料分析了可能賦存的復雜采空區域，著重選取了5個控制測點，其中，測點1位于北邊坡+120 m平臺;測點2位于東邊坡+72 m平臺;測點3位于北部運巖公路附近的+96 m平臺;測點4與測點5分別位于北部運巖公路的南段和北段。測點分布如圖3所示，測點3測線布置如圖4所示。

表1 RIS－K 2型探地雷達性能指標一覽Table1 Ground penetrating radar(GPR) performance index of RIS-K 2

圖3 采空區探測GPR測點分布Fig.3 Surveying points of goaf detection GPR

圖4 測點3的測線布置Fig.4 Line layout of point 3

3.2 數據采集處理與圖像解譯

本次數據采集處理采用RIS－K2專用軟件，該軟件能夠對所接收高頻電磁波進行自動修復、過濾以及F－K濾波和反褶積處理，從而輸出高信噪比的雷達走時剖面圖，大大提高了有用信號識別的效率，為后續的圖像解譯工作奠定基礎［5］。

GPR圖像資料解譯的基本原則是真實、全面、準確。首先是將軟件處理后得到雷達走時剖面圖轉換為探地雷達剖面圖(波形圖和灰度圖)，再根據反射波組在圖像中呈現出的波形與灰度特性，加以同相軸定位追蹤，從而推斷出異常反射波組的巖土介質信息，并依據多組剖面解譯信息建立起整個探測區域的空間探測結果［5-6］。以下以測點3為例進行解釋說明。

測點3位于北部運巖公路附近的+96 m平臺。圖4為測點3的測線布置圖，其中橫向測線在邊坡平臺上平行分布，共3條，間距為3 m;而縱向測線則垂直橫向測線布置，共8條，間距同樣為3 m。

橫向測線探測資料如圖5、圖6和圖7所示，通過分析可以看出，電磁波在正常類連續介質中有著能量逐漸衰減、波幅較小等規律，而+96 m平臺下4 m到18 m的范圍內電磁波波形雜亂，波幅變化異常，由此推測該處可能存在破碎帶或老窿充填區;平臺下28 m以下部分突然出現了雙曲線形態的強反射回波信號特征，這是穹形頂空洞在GPR波形圖上出現的典型特征［7］，據此可推測該處可能存在采空區，寬度至少有4.5 m，另外從空區的分界線分析判斷有可能是巷道空區或2次垮塌形成的采空區。縱向測線探測資料圖8和圖9所示，通過分析可以看出，+96 m平臺下4 m到18 m的范圍內存在大量的破碎充填物;平臺下28 m以下存在巷道空區或2次垮塌形成的采空區，其位置與橫向測線所測得的位置高度吻合。通過綜合分析橫向測線和縱向測線的灰度圖及波形圖，就能夠基本確定采空區的位置及其形態。

圖5 Z1橫向測線GPR灰度圖和波形圖Fig.5 GPR grayscale and waveform of horizontal line Z1

圖6 Z2橫向測線GPR灰度圖和波形圖Fig.6 GPR grayscale and waveform of horizontal line Z2

圖7 Z3橫向測線GPR灰度圖和波形圖Fig.7 GPR grayscale and waveform of horizontal line Z3

圖8 H1～H3橫向測線GPR灰度圖和波形圖Fig.8 GPR grayscale and waveform of vertical line H 1～H 3

圖9 H4～H8橫向測線GPR灰度圖和波形模擬圖Fig.9 GPR grayscale and waveform of vertical line H 4～H 8

在探測現場能看到一些意外揭露且已經延伸至地表的采空區，如圖10所示。這些采空區大多夾雜著碎石，雖然位置較為隱蔽，寬度也不大，但從其空間位置、揭露標高以及延展方向來看，基本上能和本次探測的結果聯系起來，從而進一步證明了本次采空區探測工程的可靠性。

本次現場探測的5個測點結果具體歸納如下:測點1，探測地點位于北邊坡+120 m平臺，采空區距探測地表高度約為17 m，且在距探測地表3 m處有1條巖(土)體分層，采空區內含有大量的充填物，且上部較破碎。測點2，探測地點位于東邊坡的+72 m平臺，探測深度范圍內未發現明顯空區，距探測地表2 m左右有少量破碎巖(土)體，基本不構成安全隱患。測點3，探測地點位于靠近北部運巖公路北邊坡的+96 m平臺，空區距探測地表高度約為28 m，+96 m平臺下4m到18 m的范圍內存在破碎層或破碎帶，空區拱頂位于+96 m平臺與+108 m平臺之間，并且空區寬度至少有4.5 m。測點4，探測地點位于靠近北部運巖公路南段，距探測地表高度20～27 m的范圍內存在3個巷道空區或2次應力釋放空區，空區寬度為3 m左右。測點5，探測地點位于北部運巖公路北段，探測深度范圍內未發現明顯空區，距探測地表20 m左右以下區域可能存在破碎帶或古河道的充填區，寬度約為8 m。

圖10 揭露空區及其破碎充填物Fig.1 0 Exposed goaf and its fracture fillings

4 結語

在地下轉露天開采過程中，可能遇到的復雜采空區包括規劃開采形成的采空區、古代開采和非法盜采形成的采空區、充填不充分或充填體沉降形成的采空區以及2次垮塌形成的采空區。由于后3種采空區的形成區域未知或需要預測，對礦山的正常生產構成巨大的安全隱患。工程實踐表明，GPR法對于探測地表淺部各種采空區、充填破碎帶等不良地質結構不僅結果準確可靠，而且經濟高效，在同類探測方法中優勢明顯，尤其是在探測地下轉露天礦山中賦存的復雜采空區效果顯著，值得在相關領域應用推廣。

［1］ 王海君，李克民，陳樹召，等.地下轉露天開采的采空區危害與探測技術探討［J］.煤炭工程，2010，23(10):69-71.

Wang Haijun，Li Kemin，Chen Shuzhao，et al.Discussion on goaf dangers and detection technology for mining converted from underground to surface［J］.Coal Engineering，2010，23(10):69-71.

［2］ 門業凱，范廣海，邵安林，等.露天鐵礦采空區綜合探測方法研究與應用［J］.金屬礦山，2010(5):124-127.

Men Yekai，Fan Guanghai，Shao Anlin，et al.The study of integrated detection method to themined-out areas in open pit ironsmine and its application［J］.Metal Mine，2010(5):124-127.

［3］ 李 彬，庹先國，張 賡，等.高密度電法在礦山地質災害勘察中的應用［J］.金屬礦山，2013(6):87-89.

Li Bin，Tuo Xianguo，Zhang Geng，et al.Application of high density resistivitymethod to the survey of geologic disasters in mine［J］.Metal Mine，2013(6):87-89.

［4］ 劉 圓，朱建新，彭 林.隱覆采空區的地質雷達探測［J］.現代礦業，2011(7):66-67.

Liu Yuan，Zhu Jianxin，Peng Lin.Hidden geological radar detection ofmined-out area［J］.Modern Mine，2011(7):66-67.

［5］ 董 杰，孟慶生，尹明泉.綜合物探方法在采空塌陷區調查中的應用［J］.物探與化探，2013，37(3):557-560.

Dong Jie，Meng Qingsheng，Ying Mingquan.The application of integrated geophysicalmethods to investigatingmine-out collapse geological disasters in a goldmine［J］.Geophysical＆Geochemical Exploration，2013，37(3):557-560.

［6］ 丁新啟，張 華，李 遠.采空區影響下的高邊坡穩定性綜合分析［J］.有色金屬:礦山部分，2009，61(3):39-43.

Ding Xinqi，Zhang Hua，Li Yuan.Analysis on the stability of high slope effected by mined-out area［J］.Nonferrous Metals:Mining Section，2009(3):39-43.

［7］ 程久龍，胡克峰，王玉和，等.探地雷達探測地下采空區的研究［J］.巖土力學，2004，25(S1):79-82.

Cheng Jiulong，Hu Kefeng，Wang Yuhe，et al.Research on detecting of underground mined-out areas by using GPR［J］.Rock and Soil Mechanics，2004，25(S1):79-82.

［8］ 劉敦文.地下巖體工程災害隱患雷達探測與控制研究［D］.長沙:中南大學，2001.

Liu Dunwen.Research on GPR Detection and Hazard Control of Underground Rock Mass Engineering［D］.Changsha:Central South U-niversity，2001.

［9］ 金解放，趙 奎.基于探地雷達原位試驗的空區探測研究［J］.礦冶工程，2011，31(5):4-8.

Jin Jiefang，Zhao Kui.Study on mined-out area detection based on in-situ testwith ground penetrating radar［J］.Mining and Metallurgical Engineering，2011，31(5):4-8.