煤礦開采地表沉陷盆地邊界的再認識*

李春意，車宇航，王石巖

(1. 河南理工大學 測繪與國土信息工程學院，河南 焦作 454000；2.塔斯馬尼亞大學 科學與工程學院，塔斯馬尼亞 霍巴特 7005； 3.河南省基礎地理信息中心, 河南 鄭州 450003)

0 引言

煤炭是我國主要的能源和工業(yè)原料，在為我國經(jīng)濟和社會帶來飛速發(fā)展的同時，大規(guī)模井工煤礦的開采破壞了圍巖的應力平衡狀態(tài)，誘發(fā)巖層發(fā)生移動、變形和破壞，覆巖移動傳播至地表會對建(構(gòu))筑物、水體和鐵路產(chǎn)生不同程度的影響[1]。根據(jù)一般地質(zhì)采礦條件下沉陷盆地的特征，可把全盆地劃分為3個邊界：最外邊界、危險移動邊界和裂縫邊界，分別以邊界角、移動角和裂縫角3個角量參數(shù)來界定[2]。在建(構(gòu))筑物下、鐵路下壓煤開采設計和保護煤柱留設時，需要考慮邊界角和移動角這2個關(guān)鍵參數(shù)，即準確確定沉陷盆地的邊界和最外側(cè)臨界變形值點的位置，尤其當受護對象級別較高時更是如此。

在保護煤柱設計時，我國《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規(guī)范》(安監(jiān)總煤裝[2017]66號)(以下簡稱“三下”開采規(guī)范)中依據(jù)地表建(構(gòu))筑物的重要性、用途以及受采動影響后可能產(chǎn)生的后果，把受護對象劃分為Ⅰ～Ⅴ5個等級。對于一般的建(構(gòu))筑物、鐵路保護煤柱采用移動角進行留設，但特級受護對象保護煤柱應采用邊界角留設，如高速鐵路、核電站以及大型水庫等[3]。從建筑物采動損害評價的角度分析，移動角可以確定建(構(gòu))筑物的臨界變形值，邊界角可以確定沉陷盆地的范圍，在處理工農(nóng)關(guān)系和計算采動損害建筑物補償費用時，較小的邊界角偏差可能會使沉陷盆地范圍顯著擴大，增加建筑物補償費用。此外，我國以10 mm下沉點作為沉陷盆地邊界還有待進一步商榷，因為在厚含水松散層礦區(qū)，沉陷盆地邊緣雖然下沉量很小，但水平移動量卻相對較大[4-6]。因此，從工程的實際需要出發(fā)，科學合理地確定沉陷盆地邊界和邊界角具有研究價值和現(xiàn)實意義。

1 地表沉陷盆地邊界的界定

我國“三下”開采規(guī)范中明確給出了邊界角的概念，即在充分或接近充分采動條件下，沉陷盆地主斷面上下沉10 mm的點與采空區(qū)邊界之間的連線和水平線在煤柱一側(cè)的夾角，如圖1和圖2所示，其中：A，B，C，D為沉陷盆地邊界；δ0，β0，γ0分別為走向、下山和上山方向的邊界角，(°)；φ為松散層移動角，(°)。

圖1 走向方向邊界角示意Fig.1 Schematic shape oflimit angle alone strike

圖2 傾向方向邊界角示意Fig.2 Schematic shape of limit angle alone dip

為了求出邊界角，首先需要確定地表沉陷盆地的邊界。關(guān)于煤礦開采后地表沉陷盆地的邊界國內(nèi)外認識并不統(tǒng)一，如表1所示。

表1 不同國家沉陷盆地邊界的界定Table 1 Definition of surface trough margin for different countries

注：Wm為地表最大下沉值。

我國 “三下”開采規(guī)范雖然把沉陷盆地主斷面上10 mm的下沉點作為沉陷盆地的邊界，但并未給出理論依據(jù)。實際上，在地表沉陷日常觀測中，觀測站各點的高程按四等水準施測[13]，從技術(shù)層面考慮，四等水準要求每千米往返測高差中數(shù)的全中誤差不超過±10 mm。因此，我國以10 mm下沉點作為沉陷盆地的邊界是基于施測等級和施測精度要求而確定的。但在特定條件下，如礦區(qū)覆蓋有超過100 m的含水松散層，若仍以10 mm下沉點作為沉陷盆地的邊界欠妥。因為淮南、淮北和邢東煤礦實測結(jié)果顯示：當?shù)V區(qū)范圍內(nèi)有厚含水松散層覆蓋時，下沉系數(shù)偏大，沉陷盆地收斂非常緩慢，盆地邊緣實測水平移動值遠大于下沉值[4-5,14-15]，如圖3和圖4所示。

圖3 邢東煤礦沉陷盆地邊緣下沉和水平移對比Fig.3 Comparison between settlement and horizontal movement values at the edge of subsidence basin in Xingdong pit

圖4 淮北和淮南礦區(qū)沉陷盆地邊緣下沉值和水平移動對比Fig.4 Comparison between settlement and horizontal movement values at the edge of subsidence basin in Huaibei and Huainan coal mining area

從建(構(gòu))筑物采動損害與評價的角度出發(fā)，上述厚含水松散層開采條件下若仍以下沉10 mm的點作為沉陷盆地的邊界是不合適的。此外，對于下沉或傾斜比較敏感的受特級保護的帶狀構(gòu)筑物(如高速鐵路)，筆者認為也不能一概按10 mm下沉點作為沉陷盆地的邊界，而應以二等水準為依據(jù)實施精密水準測量，即要求每千米往返測高差中數(shù)的全中誤差不超過±2 mm。因此，建議此時以2 mm下沉點作為沉陷盆地的邊界。

2 厚松散層下開采地表沉陷盆地邊界角的確定

通過上述分析可知，厚含水松散層下開采時，沉陷盆地邊緣由于水平移動范圍大于下沉范圍，應該以水平移動量來確定沉陷盆地的邊界。根據(jù)《煤礦測量規(guī)程》(能源煤總[1989]25號)和《礦山開采沉陷學》，對地表移動觀測站開展全面觀測時平面坐標應按照5″導線(一級導線)標準進行施測[2,13]，導線點的點位中誤差可用極坐標法來確定，方法如圖5所示，A、B分別為地表移動觀測站控制點；C為待定點；其大地坐標分別假定為A(xA,yA)，B(xB,yB)，C(xC,yC)，導線AC的長度為L，mm；其與AB的夾角為β，(°)；α為已知導線邊AB的坐標方位角，(°)；γ為導線AC的坐標方位角，(°)；其中β、L可通過實測獲取。

圖5 觀測站點位中誤差的確定Fig.5 Determination of root mean square error for surface monitoring landmark

導線AB的坐標方位角α可依據(jù)A、B點的位置通過坐標反算獲得，即：

(1)

導線AC的方位角γ為：

γ=α-β

(2)

采用坐標正算，C點的坐標可表達為：

(3)

根據(jù)誤差傳播定律，C點的橫坐標xC的點位中誤差mxC為：

(4)

式中：mL為導線邊L的測距中誤差，mm；mβ為β的測角中誤差，mm；ρ=206 265″為常數(shù)。

為了求取式(4)的極值，該式可化為：

(5)

根據(jù)《工程測量規(guī)范》(GB 50026—2007)[16]，一級導線測角中誤差為5″，最大導線邊不超過0.5 km，測距中誤差不超過15 mm，由式(5)可得，待定點C點位中誤差mxC的極大值mxCmax為：

mxCmax=mL=15 mm

(6)

同理可得待定點C縱坐標的點位中誤差極大值myCmax=mxCmax=15 mm。綜上，厚含水松散層開采條件下沉陷盆地邊緣15 mm的水平移動點可作為沉陷盆地的邊界，進而再計算邊界角。

3 采動程度與邊界角的協(xié)同作用關(guān)系

3.1 工程背景

在開采沉陷學中，采動程度n是指工作面傾向方向尺寸D(假定工作面為走向長壁開采)與平均開采深度H0的比值，即n=D/H0，無量綱。n≤1/4～1/3時為極不充分采動；1/4～1/31.2～1.4時為超充分采動[17]。為了探求采動程度與邊界角尤其是與上山和下山邊界角之間的協(xié)同作用關(guān)系，采用UDEC數(shù)值分析方法來進行研究。作者曾以霍州煤電辛置礦五采區(qū)為工程背景，結(jié)合實測數(shù)據(jù)研究了間歇開采地表沉陷的規(guī)律及巖移機理[16]，本研究仍以該區(qū)域為工程背景，工作面布置如圖6所示。模擬區(qū)域涉及9個回采工作面，煤層開采厚度3 m，煤層傾角7°，平均開采深度455 m。

圖6 工作面布置Fig.6 Arrangement diagram of longwall coal panel

根據(jù)工作面地質(zhì)采礦條件，并考慮邊界角等因素，決定沿地層傾向主斷面建立尺寸為1 800 m×580 m的平面應變數(shù)值分析模型，模擬從極不充分開采至充分開采的全過程，采出空間從70 m(首采工作面603B，采動程度0.15)至550 m(回采全部9個工作面，采動程度1.2)，開采順序、覆巖地層結(jié)構(gòu)及工作面幾何參數(shù)如圖7和表2所示。

圖7 模擬地層結(jié)構(gòu)及工作面回采順序Fig.7 Simulated overlying strata and extraction sequences of longwall coal panels

考慮到煤層開采后巖石在不同圍壓條件下具有明顯的彈塑性變形特征，并發(fā)生塑性破壞和拉剪破壞，本數(shù)值模型中塊體采用摩爾-庫倫本構(gòu)模型，節(jié)理采用面接觸的庫倫滑移模型，巖體的破壞準則采用Mohr-Coulomb屈服準則，具體節(jié)理參數(shù)、材料塊體參數(shù)和實際開采情況見文獻[18]。采用給出的參數(shù)按實際開采情況進行模擬，全采后地表沉陷數(shù)值模擬結(jié)果與實測結(jié)果對比，如圖8所示，模擬相對中誤差12.8%，表明二者在宏觀上是吻合的，能夠開展工作面采動程度與邊界角協(xié)同作用關(guān)系研究。

表2 回采工作面幾何參數(shù)及邊界角計算Table 2 Geometric parameters of excavated panels and calculated limit angles

3.2 采動程度與邊界角之間的協(xié)同作用分析

為了直觀反映采動程度與地表沉陷盆地邊界角之間的協(xié)同作用過程，針對每次開挖，當應力迭代計算達到平衡后，提取地表下沉曲線并繪制對應的邊界角，如圖9所示。

圖8 下沉值實測結(jié)果與數(shù)值模擬結(jié)果對比Fig.8 Comparison of in-situ and numerical simulation subsidence values

由圖9可知，當工作面開采70 m時(采動程度僅為0.15)，屬于極不充分開采，下山邊界角為46°，上山邊界角為47.1°。當開挖空間擴大至170 m(采動程度n=0.37)時，為非充分開采，與極不充分開采相比，下山方向和上山方向沉陷盆地邊界角均有所減小，分別為45.4°和44.7°。隨著采出空間的繼續(xù)擴大，無論是下山方向還是上山方向沉陷盆地邊界角均呈遞增趨勢，當開采區(qū)間為550 m時，下山邊界角為55.2°，上山邊界角為49.5°。上述現(xiàn)象也說明，極不充分開采條件下工作面采動影響范圍(工作面邊界至同側(cè)沉陷盆地邊界之間的水平距離)并沒有明顯減小，但地表最大下沉值卻很小(模擬地表最大下沉138 mm)，地表沉陷盆地十分平緩，這對于建(構(gòu))筑物下壓煤開采是非常有利的，與文獻[19]中實測地表沉陷特征相一致。

圖9 采動程度與邊界角的協(xié)同作用Fig.9 Synergistic interaction between mining degree and limit angle

根據(jù)上述離散元數(shù)值分析結(jié)果，采動程度與上、下山邊界角之間的協(xié)同作用關(guān)系如圖10所示。

圖10 采動程度與上下山邊界角之間的關(guān)系Fig.10 Interrelation between mining degree and limit angle at rise-side and dip-side

由圖10可知，二者均呈現(xiàn)先緩慢減小再逐漸增加的趨勢，且通過研究發(fā)現(xiàn)二者之間均服從Lorentz函數(shù)分布規(guī)律，采動程度與上山、下山邊界角之間的回歸函數(shù)模型如式(7)和式(8)所示：

(7)

(8)

式中各參數(shù)含義同前，式(7)和式(8)的相關(guān)度分別為0.88和0.99，說明Lorentz函數(shù)模型能夠表達采動程度與上、下山邊界角之間的內(nèi)在聯(lián)系。

基于式(7)和式(8)，在礦井后續(xù)開采過程中，可依據(jù)采動程度n或開采尺寸對上、下山邊界角進行預測，確定沉陷范圍，以便提前采取有效措施對沉陷區(qū)范圍內(nèi)的受護對象進行加固或維修，排除險情，確保地面建(構(gòu))筑物的安全、正常使用。

4 結(jié)論

1)不同國家對邊界角的認識并不統(tǒng)一，沉陷盆地邊界角隨受護對象的工程屬性和地質(zhì)采礦條件的差異而變化。我國一般以10 mm下沉點作為沉陷盆地的邊界，對于特級保護對象，建議以二等水準為依據(jù)實施精密水準測量，以2 mm下沉點作為沉陷盆地邊界。

2)建立了厚含水松散層條件下煤礦開采地表沉陷監(jiān)測點點位中誤差函數(shù)模型，沉陷盆地邊緣15 mm的水平移動點可作為沉陷盆地的邊界，并由此計算邊界角。

3)隨著采動程度的增加，沉陷盆地邊界角呈現(xiàn)先緩慢減小再逐漸增加的趨勢，二者服從Lorentz函數(shù)分布規(guī)律，相關(guān)度分別為0.88和0.99，該函數(shù)模型能夠表達采動程度與研究區(qū)域上、下山邊界角之間的內(nèi)在聯(lián)系。