高鐵緊鄰深厚采空區的綜合勘察技術及應用★

岳建剛

(中鐵第四勘察設計院集團有限公司，湖北 武漢 430000)

鐵路通過煤礦采空區時，一般以繞避為主，同時，通過保留一定寬度的安全礦柱。近年來，隨著高速鐵路跨越式發展，選線時不得不穿越或緊鄰采空區。采空區本身具有隱蔽性、復雜性，尤其對采深超過1 000 m的深厚采空區，如何準確地查明其空間范圍及邊界，是鐵路工程地質勘察的一大難題。

對于深厚采空區，深孔鉆探實施難度大、周期長、代價大，常規勘察以電法、磁法及地震勘探等手段為主，但受制于設備探測深度的局限性，很難準確查明其采空區的空間分布特征。

本文以魯南高鐵古城煤礦為例，針對采深大于1 000 m的深厚采空區，創新勘察手段，通過InSAR技術、深源物探技術，圈定了采空區的范圍及規模，并結合現場調查進行驗證分析，確定采空區的周界，為鐵路選線奠定基礎[1]。

1 工程概況

1.1 采空區概況

華北某設計時速為350 km/h的高速鐵路，沿線分布有多處采空區，煤層埋深1 000 m～1 200 m，厚約8 m～10 m，因沿線設站條件受限等原因，線路需穿行于煤礦采空區之間的夾心地，其中，DK276～DK280段左右兩側均分布有采空區(見圖1)，且近年來煤礦一直有采掘活動，為確保鐵路工程建設及運營安全，需準確查明煤礦采空區的邊界及范圍[2]。

以線路右側古城煤礦為例，目前開采的煤層為3煤層，主要集中在南部的31,32采區，采深1 200 m，開采方式為條帶開采，全部采用綜采放頂煤工藝開采。

1.2 開采歷史

古城煤礦3107,3209工作面(如圖2所示)距離本線最近，其采空區對鐵路安全影響最大。根據竣工資料顯示，3107工作面距離線路最近約920 m，地面標高約55 m，工作面標高-1 080 m～-1 165 m，煤層厚度為8.6 m，傾向東南，傾角7°。該工作面自2017年1月開始回采，2018年3月回采結束，回采長度701 m。采用保護煤柱條帶式開采，綜采放頂煤一次采全高生產工藝。

3209工作面距離線路最近約835 m，南北走向，位于3107工作面東側300 m，地面標高約55 m，工作面標高-1 189 m～-1 196 m，寬度174 m，煤層厚度8.7 m，走向近南北向，偏北西向，傾向北東東。采用保護性煤柱條帶式開采，綜采放頂煤一次采全高生產工藝。該工作面2017年9月開始回采，2018年9月停采，退采長度470 m，靠線路側的工作面停采時間已近36個月。

2 InSAR解譯

2.1 解譯方法

合成孔徑雷達干涉測量技術(Interferometric Synthetic Aperture Radar，InSAR)是利用雷達反射，對地表典型點的反射與固定點的相位差，分析兩者之間的變形差，可宏觀上反映一定規模的地表變形特征和趨勢[3]。采空區形成后，地表必然會出現不同程度的變形，在以“山西”為代表的煤礦采空區研究中曾采用InSAR技術開展過采空區地表沉降監測工作，并且取得了良好的效果，是當前研究地表變形的有效技術手段。

2.2 數據源

針對古城煤礦采空區，根據《地面沉降干涉雷達數據處理技術規程》的有關規定，充分考慮大范圍、精細化監測的應用需求，結合煤礦開采歷史，選擇SARadarSat-2高分數據及Sentinel-1號衛星數據作為主要數據源。共獲取了2016年—2018年共15期SARadarSat-2高分影像以及2019年—2020年共32期Sentinel-1號衛星影像。

SARadarSat-2數據及Sentinel-1數據均無模糊、散焦、條帶、增益過度等現象存在；所有數據幅度影像圖覆蓋范圍保持一致，無跳變現象存在；且重復軌道數據影像相干性好，無數量缺失，數據公共區域滿足工作區監測需求[4-5]。

2.3 解譯分析

古城煤礦位于DK276+350～DK279+926段右側，根據古城煤礦地表變形監測報告，條帶開采下沉系數0.11～0.15，計算得到預估沉降量約0.8 m。古城煤礦3107工作面從2017年1月開采至2018年3月停采，3209工作面從2017年9月開采至2018年9月停采。

根據2017年—2018年度的年累積地面沉降等值線圖(如圖3(a)，圖3(b)所示)進行分析，可得出以下結論：

1)2017年—2018年期間，古城煤礦3107及3209工作面采空區沉降漏斗中心向線路方向偏移。

2)地表沉降漏斗中心位于3107工作面的北面，3209工作面的西側，移動盆地呈東西向展布，距離線路最近為1.06 km。沉降漏斗中心年累計沉降量持續保持在200 mm左右。

3)InSAR解譯分析成果與3107及3209工作面開采時間、開采時序、形態等資料相符，InSAR解譯成果可很好地反映采空區的時序特征。

根據2019年—2020年古城煤礦段(DK276+350～DK279+926)半月累計地面沉降等值線圖(如圖4(a),圖4(b)所示)進行分析，可得出以下結論：

1)2019年—2020年3107及3209工作面已開采停止，此處沉降中心位置未變化，下沉量未明顯增大。

2)新的地表沉降漏斗中心已偏離線路，并向西移動，逐漸與新開采區相對應。

3)舊采掘面所在區域地面年沉降速率逐漸減緩，新采掘面所在區域年地面沉降速率加快。

2.4 解譯成果

1)根據2017年—2020年InSAR解譯成果，DK276+350～DK279+926段累計最大沉降量為0.5 m，與煤礦監測資料中最大沉降0.546 m基本一致，說明古城煤礦在該時間段內持續存在采掘活動。

2)結合既有礦區資料，對比礦區開采面的分布和InSAR解譯的沉降等值線中心位置，每年的沉降漏斗基本和當年主要開采面重合，且沉降漏斗基本與工作面中心一致，說明地表變形沉降是由采空引起。

3)2017年—2020年InSAR解譯采空區變形趨勢特征表明，InSAR變形中心區位于3107工作面偏北側、3209西側，3107與3209采空區引起的地面變形已于2019年初穩沉，2020年的InSAR解譯成果已無變形發展。

3 深源物探

3.1 探測方法

3.1.1 大噸位可控源寬線地震反射波法

高分辨率地震勘探技術可以滿足“數千米”勘探深度，煤層被采出后原有狀態與特征發生變化，采空區附近圍巖發生塌落、破碎，使得地層的波阻抗發生變化，從而導致其地球物理特征的變化。考慮勘探深度以及精度的要求，選用了當前國際上最先進的美國Mertz公司生產的大噸位28 t大能量可控震源車(如圖5所示)作為本次縱波地震勘探的激發源，實際有效勘探深度超過2 000 m。

3.1.2 可控源音頻大地電磁法(CSAMT)

可控源大地電磁法(CSAMT)是工程勘察中適用1 000 m深度以上最有效的探測手段，是通過人工控制場源的一種電磁測深法，它通過改變發射源的發射頻率達到測深目的，用測量相互正交的電場和磁場分量計算卡尼亞視電阻率和阻抗相位，根據反演電阻率和深度關系圖進行地質解釋[6]。

一般條件下，煤礦采空區形成后，會伴隨產生導水裂隙帶，地下水會將巖體中的鈣離子、鐵離子等溶離出來，采空區在視電阻率圖中表現為明顯低阻異常，而圍巖則表現為高阻狀態。本次CSAMT使用加拿大鳳凰地球物理公司V8多功能大地電磁儀(見圖6)，實際有效勘探深度大于1 500 m。

3.2 測線布置及解譯分析

根據周邊地形、地貌以及建構筑物的分布特征，結合收集的采空區資料，寬線地震反射、可控源大地電磁勘探布置原則上保證覆蓋古城煤礦3107,3209工作面，以查明采空區邊界，同時兼顧探明采空區邊界以南近線路一側是否還存在其他的采空區、構造及其發育情況。

3.2.1 地震反射波法

寬線地震反射共布設反射采集測線10條。其中南北向測線7條，覆蓋3107,3209工作面南北邊界，東西向測線3條，覆蓋工作面東西邊界。10條測線按四個寬線地震勘探區以及一個二維地震勘探測線組合，共計生成地震寬線剖面12條(見圖7)，測線長度22.7 km。

根據煤系地層各反射波的相對關系及波組特征來解釋煤層的賦存形態、斷裂構造、采空區等地質特征。解譯成果如下：

（4）基于GIS的選線工作輔助系統，其可靠性與計算結果準確性及現有數據信息有直接關系，同時也會受到用戶自身選擇與判斷的影響。只有正確的選擇與判斷，才能充分發揮工具應有的作用和效果。

1)本次地震勘探12條測線成果中僅NS505線、GCL4線未見明顯異常，其余10條測線均存在煤層反射波下錯、波組缺失以及頻率變低等現象，反映為典型的煤礦采空區特征。

2)煤礦采空區南側邊界以北整體表現出明顯的采空區特征，本次勘探南北4條測線均在平面上穿越3107工作面并有部分測線進入其北側的3201工作面。

3)南北向測線GCL302存在多處斷層構造，與古城煤礦3209采面揭示斷層平面位置一致，綜合東側GCL4線勘探成果，表明GCL302以東無采掘作業面(見圖7)。

4)東西向測線NS505無明顯煤礦采空區異常反應，推斷古城煤礦采掘工作面向南未穿越NS505測線(見圖8)。

3.2.2 可控源大地電磁法

可控源大地電磁(CSAMT)共布設測線6條(見圖9)，東西向布置1條，南北向布置5條。配合地震反射波法綜合確定采空區異常，測線長度11.7 km。

根據測線的視電阻率剖面來分析區域地層背景、基巖完整性等地質特征。解譯成果如下：

1)探測范圍內整體電阻率呈層狀分布，與區域地層地質背景一致，該段主要是以二疊系、石炭系泥巖、砂巖為主，呈層狀構造。

2)從電阻率分布規律來看，中淺層(500以淺)地層電阻率較低并且縱橫向變化不大，反映500 m以上基巖較完整。

3)CSAMT5-9～CSAMT5-17段中深部整體電阻率極低，反映為典型的煤礦采空區特征，推斷為3209采空區異常。

4)CSAMT5-22～CSAMT5-28段在中深部仍存在明顯的低電阻率異常特征，推斷為基巖破碎異常區(見圖10)。

3.3 解譯成果

通過對比以上兩種深源物探成果，對古城煤礦3107及3209工作面采空區范圍進行綜合分析(見圖11)。

1)通過寬線地震反射及可控源大地電磁探測，查明了探測區域內的采空區異常范圍，并圈定了3109及3207采空區的范圍及邊界；且兩種方法對采空區解譯的范圍基本一致，說明采空區的勘探范圍準確有效。

2)由于采空區電性在空間分布上具有明顯的體積效應，受電磁方法理論限制，CSAMT圈定采空區邊界較地震反射圈定的邊界偏大。

3)寬線地震反射勘探采用了大能量(28 t可控源)、小道距(5 m)、全數字檢波技術進行數據采集，勘探精度高于常規煤田地震反射勘探精度，理論上水平分辨率可達2.5 m、縱向分辨率5.0 m。

4)對比采空區綜合物探解釋成果與古城煤礦采掘工程平面圖揭示采空區平面位置，地震反射圈定3107,3209采面采空區邊界控制點較采掘工程平面采掘邊界大，這與實際煤礦采空區“冒落帶”發育范圍有直接關系，總體來看在采空區長軸方向“冒落帶”變形范圍更大，這也很好的證明了地震反射勘探結果的準確性與可行性。

4 調查監測

采空區形成后往往引起地表一定程度的變形，主要表現為房屋、道路等建構筑物開裂。根據現場3107工作面附近地表調查，采空區地表范圍內道路可見變形主要為裂縫，裂縫寬度2 mm～30 mm，裂縫長度從幾米到幾十米不等；房屋主要為豎向裂縫，裂縫寬度1 mm～5 mm[7]。

1)2019年7月～8月，通過現場走訪、問詢，對古城煤礦3107及3209工作面周邊的地表構筑物進行了調查。通過調查發現，現場地表變形以沉降變形為主，采空區中心位于基本農田內，未見塌陷坑，不影響小麥與玉米的種植，小沂河未發生斷流現象。根據現場在距離3209工作面南側200 m～500 m水平范圍內的楊莊村房屋及小路裂縫的調查情況，對差異沉降十分敏感的砌體結構、條形基礎民用建筑，局部建筑房屋存在張開裂縫(紋)，裂縫寬2 mm～5 mm(見圖12)。現場調查的房屋裂縫最遠點至3209工作面邊界水平距離約400 m。地表變形量較小，村莊居民正常生活、道路交通和農業生產未受影響。

2)2020年5月、8月到現場再次進行地表變形調查，采空區中心種植玉米地，整齊規劃，村民安居樂業，裂縫(紋)無加劇，范圍也未擴展，也未發現新的變形點，與采空區的開采時序吻合對應。

3)《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規范》(安監總煤裝[2017]66號)中對于長度或者變形縫區段內長度不大于20 m的磚混結構建筑物損壞等級由輕微至嚴重分為Ⅰ,Ⅱ,Ⅲ,Ⅳ四個等級，其中Ⅰ級損壞等級地表變形值的限值為：水平變形ε≤±2 mm/m，傾斜i≤±3 mm/m，曲率K≤±0.2×10-3/m。

4)3107及3209工作面南側地表建筑物變形值均在Ⅰ級損壞等級范圍內，Ⅰ級損壞區的影響范圍遠離線路，距離Ⅰ級損壞區影響邊界420 m(見圖13)。

5 結語

地質選線對于鐵路走向具有重要的意義，本文以古城煤礦3107和3209工作面為例，通過“空-源-地”三位一體的手段，即從InSAR的高位“空”中影像分析地表區域變形特征，利用可控源音頻大地電磁法和大噸位可控源地震縱波反射法的深層物理“源”地球物理勘探方法，查明了古城煤礦深厚采空區的范圍及邊界，同時通過“地”表建構筑物的調查研究，對采空區的范圍邊界進行驗證。勘探獲取的煤礦采空區邊界基本與竣工資料一致，確保了采空區的范圍合理，為下一步鐵路選線及采空區評價提供可靠的地質依據[8-9]。

1)3107工作面竣工資料顯示的工作面邊界距線路920 m，物探揭示煤層及異常區距線路900 m～920 m；3209工作面竣工資料顯示邊界距線路835 m，物探揭示煤層及異常區距線路830 m～835 m；物探揭示的采空區邊界與竣工資料顯示一致，說明深源物探方法能夠較好地適應超深采空區的勘察。

2)InSAR解譯中心沉降區距離線路約1 000 m，2018年底 的InSAR解譯的中心沉降區距離線路約1 100 m，InSAR解譯的影響范圍及形態，以及物探的解譯范圍與3107及3209工作面較為吻合。

3)地面調查房屋裂縫、裂紋(線路位于Ⅰ級損壞區影響邊界外420 m)與采空區的范圍也有較好的對應性，表明3107及3209采空區位置是比較準確、可靠的。