丁 湘，周新河，閆 鑫，吳 訪，蒲治國，孟廣山

(1.中煤能源研究院有限責任公司，陜西省西安市，710054；2.中煤沖擊地壓與水害防治研究中心，內蒙古自治區(qū)鄂爾多斯市, 017000；3.中天合創(chuàng)能源有限責任公司，內蒙古自治區(qū)鄂爾多斯市, 017300)

煤層開采擾動導致原巖應力失去平衡，底板巖層在應力作用下向采空區(qū)自由空間運移，導致底板巖體發(fā)生剪切膨脹破壞，破裂的巖層將局部或全部喪失隔水能力[1]。針對煤層采動導致底板破壞的問題，相關專家學者進行了大量研究。孟祥瑞等[2]研究了煤層采動后的應力分布規(guī)律和底板破壞深度；張金才[3]從彈性力學的角度出發(fā)，基于半無限體理論研究了底板破壞機理；張風達等[4]結合塑性力學滑移線場理論研究了影響底板破壞深度的因素；馬文著等[5]指出影響底板巖體破壞的主要因素包括開采深度、開采范圍和煤層產狀等，底板巖層結構完整的條件下，塑性破壞區(qū)的形態(tài)接近“勺形”；宋文成等[6]基于彈塑性理論構建了底板破壞判據(jù)模型，并用FLAC3D軟件構建數(shù)值模擬模型分析了底板破壞形態(tài)特征。

多煤層礦井上部煤層開采會引起底板巖層的破壞和移動變形，導致圍巖應力平衡狀態(tài)發(fā)生變化和底板巖層破壞，從而加重上部采空區(qū)積水對其下部煤層采掘作業(yè)的威脅程度。底板破壞帶深度的確定是評價下部煤層采掘活動受上覆采空區(qū)積水威脅程度的重要指標，也是進行采空區(qū)積水探放鉆孔參數(shù)設計的關鍵數(shù)據(jù)，因此，需要確定底板采動破壞帶深度，為下部煤層采掘時的采空區(qū)積水防治提供基礎數(shù)據(jù)[7-9]。

1 工作面概況

根據(jù)井下實際采掘資料分析，2-1號煤層一盤區(qū)21102工作面和21103工作面采空區(qū)內地勢低洼處存在積水，積水面積約67.66萬m2，積水量約35.09萬m3，上部采空區(qū)積水對下部22103工作面的掘進及回采構成較大威脅。

圖1 各工作面位置關系

2 2-1號煤層底板巖層破壞帶發(fā)育深度

2.1 2-1號煤層底板巖層特征

2.2 2-1號煤層底板破壞帶深度計算

2.2.1 “三下”規(guī)范經驗公式

《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規(guī)范》(以下簡稱“三下”規(guī)范)中提供了多種底板破壞深度計算方法。在工作面斜長30～205 m、采深100～1 000 m條件下，有以下3個公式均可用來計算煤層采后底板破壞帶發(fā)育深度。

h=0.303L0.8

(1)

h=0.700 7+0.107 9L

(2)

(3)

式中：h——底板采動破壞深度，m；

H——開采深度，m；

L——工作面斜長，m；

α——煤層傾角，(°)。

2.2.2 多元線性回歸模型

由于“三下”規(guī)范經驗公式未考慮煤層開采厚度、底板抗破壞能力及構造的影響，這些因素在深埋條件下，對圍巖的應力變化具有較大影響。因此需要對深埋條件下煤層采后底板破壞深度進行重新統(tǒng)計分析。

收集了25組埋深在300 m以上的工作面回采后底板破壞深度實測值[10]，以煤層采深(H)、煤層傾角(α)、煤層采厚(M)、工作面斜長(L)、底板抗破壞能力(R)和是否有斷層構造(F)等6個因素為影響指標，構建了煤層采后底板破壞深度預測回歸模型。部分礦井工作面回采后實測底板破壞深度見表1?；貧w模型為：

表1 部分礦井工作面回采后實測底板破壞深度

(4)

式中：M——煤層采厚，m；

R——底板抗破壞能力；

F——是否有斷層構造，分別取1或0。

依據(jù)F檢驗結果，顯著性P值為0.000，水平上呈現(xiàn)顯著性，回歸方程相關性系數(shù)R為0.890，擬合程度較高，回歸模型擬合效果如圖2所示。

圖2 回歸模型擬合效果

2.2.3 塑性力學理論計算公式

根據(jù)塑性力學理論[11-14]，底板破壞帶深度可采用下述公式計算：

(5)

式中：Ls——煤層塑性區(qū)寬度，估算為0.015倍采深，m；

φ0——底板巖體權重平均內摩擦角，(°)。

將研究區(qū)相關數(shù)值代入式(1)～(5)，得到底板采動破壞帶深度計算值，見表2。

表2 研究區(qū)底板采動破壞帶深度h計算值 m

通過計算可以發(fā)現(xiàn)，依據(jù)“三下”規(guī)范經驗公式計算底板破壞帶深度為30.59～35.94 m，多元線性回歸模型計算結果為42.75～42.90 m，塑性力學理論計算公式計算結果為18.27 m。

2.3 2-1號煤層底板破壞帶數(shù)值模擬

2.3.1 建立數(shù)值模型

綜合分析21102和21103工作面內鉆孔柱狀數(shù)據(jù)，對地層結構進行整理分析，考慮到相鄰巖層的巖性及物理力學參數(shù)具有一定的相似性，因此，對鄰近薄巖層進行概化合并，利用FLAC3D有限差分軟件建立數(shù)值模擬模型。數(shù)值計算模型如圖3所示。采用分段開挖的方式模擬開挖過程，工作面沿走向(X軸正向)推進，每步沿走向開挖20 m、傾向開挖320 m的平衡狀態(tài)。

圖3 數(shù)值計算模型

2.3.2 模型結果分析

煤層的開挖使得圍巖的應力場遭到擾動，在應力重新分布的過程中，底板巖層形成一定的塑性破壞區(qū)[15-16]。煤層采動底板破壞過程數(shù)值模擬如圖4所示。工作面回采前期，隨著工作面的不斷推進，采空區(qū)范圍逐漸擴大，底板巖層在應力再平衡過程中不斷發(fā)生破壞(推進20～120 m)；當工作面回采到一定階段，煤層頂板巖層垮落壓實采空區(qū)，底板巖體卸荷到一定程度后逐漸減弱，煤層底板巖層的破壞程度逐漸減弱(推進120～140 m)，直至達到相對平衡狀態(tài)，底板破壞深度不再增大，此時底板破壞深度即為煤層采后底板破壞帶發(fā)育深度(推進160 m)。根據(jù)數(shù)值模擬結果分析，2-1號煤層采動底板破壞過程達到穩(wěn)定狀態(tài)后，底板破壞帶最大深度在16 m左右。

圖4 煤層采動底板破壞過程數(shù)值模擬

2.4 2-1號煤層底板破壞帶鉆孔窺視

2.4.1 鉆孔窺視探測方法

采用鉆孔窺視法對葫蘆素2-1號煤層底板采動破壞帶深度進行實測，考慮到需要施工采空區(qū)積水探放鉆孔，窺視工作利用2-1號煤層采空區(qū)探放水鉆孔進行。仰孔窺視底板采動破壞帶示意如圖5所示。根據(jù)底板巖層結構和破壞帶深度預估情況，考慮到煤層間距的不確定性及預估結果可能存在誤差，鉆孔可能會提前進入底板破壞帶。因此，在鉆進至15 m(垂高8.6 m，距2-1號煤層底板17.4 m)時即開始窺視，窺視時重點觀察圍巖裂隙發(fā)育情況，直至鉆進2-1號煤層采空區(qū)。

圖5 仰孔窺視底板采動破壞帶示意

2.4.2 探測過程

在22103工作面輔運巷F1鉆場施工F1-1K、F1-2K窺視鉆孔，綜合對比分析獲得2-1號煤層一盤區(qū)底板采動破壞帶深度實測結果，見表3。

表3 F1鉆場窺視鉆孔參數(shù)

2.4.3 探測結果

F1-1K窺視鉆孔鉆進過程中共窺視4次，F(xiàn)1-2K窺視鉆孔鉆進過程中共窺視1次。窺視鉆孔采動裂隙發(fā)育情況如圖6所示。由圖6可以看出，F(xiàn)1-1K孔采動裂隙發(fā)育下限在孔深11.5 m處，套管外至孔深11.5 m范圍內未發(fā)現(xiàn)采動裂隙，此處垂高6.3 m，即F1-1K孔窺視位置2-1號煤層底板采動破壞帶深度為19.6 m；F1-2K孔采動裂隙發(fā)育下限在孔深11.9 m處，套管外至孔深11.9 m范圍內未發(fā)現(xiàn)采動裂隙，此處垂高6.7 m，即，F(xiàn)1-2K孔窺視位置2-1號煤層底板采動破壞帶深度為19.0 m。綜合對比2個窺視孔，通過視覺識別出的采動裂隙下限距2-1號煤層底板最大垂距為19.6 m，即窺視位置2-1號煤層底板采動破壞深度為19.6 m。

圖6 窺視鉆孔采動裂隙發(fā)育情況

2.5 2-1號煤層底板破壞帶深度確定

綜合分析理論計算、數(shù)值模擬和鉆孔窺視探測結果，“三下”規(guī)范經驗公式、多元線性回歸模型計算結果與實際探測窺視結果相差較大，公式適用性較差。塑性力學理論公式計算結果與實際探測窺視結果基本一致，且計算方便，可適用于近距離深埋煤層底板破壞深度的計算。數(shù)值模擬結果與實際探測窺視結果較為接近，但該方法流程復雜，應用難度較大。鉆孔探測窺視結果為實測值，可信度最高，但施工及觀測困難，不宜作為常用手段。因此，葫蘆素煤礦2-1號煤層底板采動破壞帶深度取實測值19.6 m。類似條件下簡易計算可采用塑性力學理論公式。

3 結論

(1)利用“三下”規(guī)范經驗公式、多元線性回歸模型、塑性力學理論公式計算、數(shù)值模擬和鉆孔窺視探測方法，研究了深埋條件下煤層開采后底板破壞深度。

(2)根據(jù)鉆孔窺視探測、數(shù)值模擬和理論計算，綜合確定2-1號煤層開采后底板最大破壞帶深度為19.6 m。

(3)該區(qū)域近距離煤層深埋條件下“三下”規(guī)范經驗公式適用性較差。數(shù)值模擬和鉆孔窺視方法復雜，應用不便。塑性力學理論計算公式可用于該地區(qū)底板破壞帶深度計算。