彭洪濤，周俊杰

（1.潞安化工集團(tuán) 余吾煤業(yè)有限公司，山西 長治 046103；2.河北工程大學(xué) 地球科學(xué)與工程學(xué)院，河北 邯鄲 056038）

0 引言

礦井開采中，煤層覆巖運(yùn)動(dòng)規(guī)律及其導(dǎo)水裂隙帶發(fā)育高度的預(yù)測(cè)與探測(cè)，直接影響著礦井回采工作面工作效率及支護(hù)方式。特別是上覆巖層的巖體力學(xué)性質(zhì)較弱或者富含一定的礦井水，隨著開采工作的不斷推進(jìn)，煤層頂板支護(hù)和巖體破壞增強(qiáng)，導(dǎo)水裂隙帶發(fā)育高度和影響距離制約著礦井支護(hù)方式和開采效率[1-2]。如何準(zhǔn)確高效分析和評(píng)價(jià)煤層開采活動(dòng)對(duì)煤層覆巖運(yùn)動(dòng)的影響，是確定支護(hù)方式和手段的關(guān)鍵，也是對(duì)煤層頂板破壞深度確定的方法。針對(duì)覆巖運(yùn)動(dòng)破壞特征分析，錢鳴高院士和宋振騏院士分別提出“砌體梁”和“傳遞巖梁”結(jié)構(gòu)模型，為國內(nèi)外眾多學(xué)者研究采動(dòng)破壞理論[3-4]。目前，導(dǎo)水裂隙帶高度預(yù)測(cè)常采用經(jīng)驗(yàn)公式法、模擬實(shí)驗(yàn)法和現(xiàn)場(chǎng)探測(cè)等多種方法相結(jié)合的研究方法。李民峰對(duì)陜北礦區(qū)利用相似材料模擬分析覆巖中裂隙帶呈現(xiàn)的空間分布特征，認(rèn)為理論分析與相似模擬結(jié)果一致[5]；王沛楠利用理論計(jì)算、數(shù)值模擬和現(xiàn)場(chǎng)實(shí)測(cè)方法，研究礦井特厚煤層頂板破壞深度和控制理論[6]；喬倩等采用回歸分析方法預(yù)測(cè)不同主控因素下導(dǎo)水裂隙帶發(fā)育高度預(yù)測(cè)公式[7]；魏克敏等運(yùn)用有限差分法，對(duì)太平煤礦多煤層采場(chǎng)圍巖及其裂隙帶分布特征進(jìn)行研究[8]；多種煤礦開采工程實(shí)踐表明，礦井需要精細(xì)確定煤層開采中的頂板裂隙帶高度分布特征及預(yù)測(cè)方法。

在余吾煤礦南五采區(qū)的礦井地質(zhì)條件分析與評(píng)價(jià)基礎(chǔ)上，采用經(jīng)驗(yàn)公式計(jì)算、計(jì)算機(jī)數(shù)值模擬計(jì)算、相似材料模擬實(shí)測(cè)計(jì)算方式，得到了不同預(yù)測(cè)方式下南五采區(qū)工作面開采時(shí)頂板覆巖運(yùn)動(dòng)特征及煤層頂板裂隙帶發(fā)育高度數(shù)據(jù)。利用多種方法預(yù)測(cè)的導(dǎo)水裂隙帶發(fā)育高度，結(jié)合礦井其它區(qū)域的導(dǎo)水裂隙帶發(fā)育高度實(shí)測(cè)的數(shù)據(jù)，采用回歸分析方法，修訂南五采區(qū)工作面煤層厚度與導(dǎo)水裂隙帶最大發(fā)育高度的經(jīng)驗(yàn)方程，預(yù)測(cè)了南五采區(qū)導(dǎo)水裂隙帶發(fā)育高度值。

1 概 況

余吾礦井位于屯留北側(cè)，現(xiàn)主采二疊系下統(tǒng)山西組3 號(hào)煤層，平均埋藏深550～600 m，最大埋深大于800 m。煤層厚度5.00～7.25 m，平均5.99 m。煤層穩(wěn)定，頂板一般為泥巖、粉砂質(zhì)泥巖，底板為黑色泥巖、粉砂巖，老底為中細(xì)粒砂巖。南五采區(qū)存在的充水通道主要有斷層或構(gòu)造破碎帶、陷落柱、封閉不良鉆孔以及3 號(hào)煤層開采形成的導(dǎo)水裂隙等。

綜合水文地質(zhì)分析認(rèn)為，礦井正常涌水量為303 m3/h，最大涌水量為340 m3/h。礦井內(nèi)太原組灰?guī)r含水層、山西組及下石盒子組砂巖裂隙含水層是影響3 號(hào)煤開采的主要含水層，上石盒子組砂巖裂隙含水層距離3 煤較遠(yuǎn)，且之間有隔水層，厚度較大，為次要含水層。

2 導(dǎo)水裂隙帶高度預(yù)測(cè)

2.1 導(dǎo)水裂隙帶發(fā)育高度經(jīng)驗(yàn)公式計(jì)算

利用礦井和相鄰礦井相同煤層裂隙帶探查并預(yù)測(cè)成果，建立礦井及其相似礦井地質(zhì)條件的裂隙帶發(fā)育高度理論計(jì)算公式，并利用該計(jì)算模型，進(jìn)行礦井類似地質(zhì)條件區(qū)開采煤層的裂隙帶發(fā)育高度的預(yù)測(cè)。

依據(jù)余吾煤業(yè)關(guān)于3 號(hào)煤層頂板巖石力學(xué)測(cè)試的結(jié)果，覆巖主要為中硬巖。結(jié)合煤層賦存的具體條件和特征以及礦井采煤方式，依據(jù)《煤礦防治水細(xì)則》 中綜放開采導(dǎo)水裂隙帶高度計(jì)算公式，對(duì)3 號(hào)煤采動(dòng)裂隙高度進(jìn)行預(yù)計(jì)：

式中：Hli為導(dǎo)水裂隙帶高度，m；M 為煤層累計(jì)采厚，取6 m。

由公式（1）計(jì)算得，南五采區(qū)工作面以6 m厚度進(jìn)行開采時(shí)，導(dǎo)水裂隙帶發(fā)育最大高度為82.58 m。

2.2 計(jì)算機(jī)模擬預(yù)測(cè)裂隙帶發(fā)育高度

2.2.1 模型及參數(shù)

依據(jù)地質(zhì)資料，采用FLAC3D 軟件，構(gòu)建工作面地質(zhì)模型，工作面的長寬約為1 700 m×296 m，鉆孔1305 揭露3 號(hào)煤層厚度6.4 m。3 號(hào)煤層頂板中直接頂巖性主要為泥巖隔水層；基本頂主要巖性為中、細(xì)砂巖含水層；與基本頂臨近的上覆巖層巖性主要是中、細(xì)砂巖與砂質(zhì)泥巖的互層。

確定模型尺寸參數(shù)X=600 m，Y=500 m，Z=192 m，模型網(wǎng)絡(luò)共剖分單元數(shù)700 000 個(gè)，節(jié)點(diǎn)數(shù)724 271 個(gè)，所有單元類型均為8 節(jié)點(diǎn)六面體單元。模型從下到上分為19 層，各分層的物理參數(shù)見表1，初始模型如圖1 所示。開采參數(shù)中開采寬度設(shè)置為300 m，開采厚度設(shè)置為6 m。

表1 模型各分層物理力學(xué)參數(shù)Table 1 Physical and mechanical parameters of each layer of the model

圖1 計(jì)算機(jī)數(shù)值模擬初始模型Fig.1 Initial model of computer numerical simulation

2.2.2 模擬結(jié)果研究

針對(duì)采區(qū)工作面布置，模擬不同推進(jìn)距離下煤巖層應(yīng)力、應(yīng)變空間分布特征，如圖2 所示。由計(jì)算機(jī)模擬的塑性區(qū)分布圖可以看到，工作面推進(jìn)200 m 時(shí)，采空區(qū)上部主要為剪切破壞，其次為拉伸破壞，巖性為泥巖，塑性區(qū)形狀呈平頂拱型，導(dǎo)水裂隙帶發(fā)育高度約為89 m；工作面推進(jìn)300 m時(shí)，導(dǎo)水裂隙帶發(fā)育高度約為101 m，導(dǎo)水裂隙帶發(fā)育高度增加變緩，巖性為粗砂巖，塑性區(qū)形狀平頂拱型；工作面推進(jìn)400 m 時(shí)，導(dǎo)水裂隙帶發(fā)育高度約為103 m，導(dǎo)水裂隙帶發(fā)育高度基本不再增加，巖性為粗砂巖，塑性區(qū)形狀平頂拱型。

圖2 不同推進(jìn)距離下煤巖層應(yīng)力、應(yīng)變空間分布特征Fig.2 Spatial distribution characteristics of stress and strain of coal strata under different advancing distances

根據(jù)不同采區(qū)工作面推進(jìn)距離下數(shù)值模擬塑性區(qū)發(fā)育高度值，繪制南五采區(qū)導(dǎo)水裂隙帶發(fā)育高度曲線，如圖3 所示。采高6 m 時(shí)，余吾礦南五采區(qū)導(dǎo)水裂隙帶最大發(fā)育高度為103 m，裂采比約17.2。

圖3 南五采區(qū)工作面導(dǎo)水裂隙帶發(fā)育高度曲線Fig.3 The development height curve of water flowing fractured zone in the working face of South No.5 mining area

2.3 相似材料導(dǎo)水裂隙帶發(fā)育高度預(yù)測(cè)

相似材料模型按照1∶200 的比例，推進(jìn)長度為320 m，以砂子為骨料，石膏、碳酸鈣為粘結(jié)材料鋪設(shè)模型，根據(jù)鋪設(shè)巖層的抗壓強(qiáng)度選擇配比方式，再根據(jù)模型的大小及巖層厚度計(jì)算出砂子、石膏、碳酸鈣和水的用量。

模型在開挖320 m 后，巖層的采動(dòng)破壞特征如圖4 所示。開采過程中裂隙富集區(qū)主要集中在前后煤壁，裂隙的發(fā)育程度隨著開采工作面的推進(jìn)，而逐漸發(fā)育，且采空區(qū)越大，裂隙的密度就越大。采厚6 m 的模型，導(dǎo)水裂隙帶發(fā)育的最大高度為105 m 左右；上覆巖層的破斷角，在開切眼一側(cè)為69.1°，在工作面另一端為70.8°；冒落帶高度為38.2 m，冒采比為6.37，導(dǎo)水裂隙帶高度為105 m，裂采比為17.5。

圖4 相似材料模型覆巖采動(dòng)破壞特征Fig.4 Mining failure characteristics of overburden rock in similar material model

3 導(dǎo)水裂隙帶發(fā)育高度綜合分析

在理論公式計(jì)算、數(shù)值模擬及物理模擬的基礎(chǔ)上，結(jié)合礦井其他區(qū)域及相似條件的其他礦井實(shí)測(cè)數(shù)據(jù)，分析認(rèn)為，南五采區(qū)3 號(hào)煤層采高6 m，通過理論公式計(jì)算導(dǎo)水裂隙帶最大高度為82.58 m；計(jì)算機(jī)數(shù)值模擬的最大裂隙高度為103 m；相似材料模擬最大裂隙高度為105 m 左右。結(jié)合礦井在其它工作面（S2106、S2107）進(jìn)行的導(dǎo)水裂隙帶高度現(xiàn)場(chǎng)探查結(jié)果（表2），可以看到，該礦井在采高6 m 左右時(shí)，導(dǎo)水裂隙帶的發(fā)育高度遠(yuǎn)超過經(jīng)驗(yàn)公式的計(jì)算值。

表2 礦井其他采區(qū)裂隙高度現(xiàn)場(chǎng)探測(cè)結(jié)果Table 2 Field detection results of fracture height in other mining areas

導(dǎo)水裂隙帶發(fā)育高度的綜合分析，余吾煤礦礦井裂隙發(fā)育最大高度已不能將經(jīng)驗(yàn)公式預(yù)測(cè)值作為發(fā)育高度預(yù)測(cè)的主要依據(jù)。綜合對(duì)比經(jīng)驗(yàn)公式、數(shù)值模擬、相似材料模擬以及礦井其他采區(qū)裂隙高度現(xiàn)場(chǎng)探查值，基于工程安全考慮，以多種發(fā)育高度預(yù)測(cè)成果與以往礦井現(xiàn)場(chǎng)探查成果進(jìn)行一元線性擬合。煤礦煤層采厚與裂隙高度散點(diǎn)圖如圖5 所示。

圖5 余吾煤礦煤層采厚與裂隙高度散點(diǎn)圖Fig.5 Scatter plot of coal seam mining thickness and fracture height in Yuwu Mine

通過數(shù)據(jù)一元回歸分析，得到礦井煤層采厚（M）與導(dǎo)水裂隙帶最大發(fā)育高度（H）的公式：

綜合分析南五采區(qū)及周圍鉆孔資料，3 號(hào)煤厚度范圍為5.44（鉆孔LA010）～6.4 m（鉆孔1306），則最大導(dǎo)水裂隙帶高度范圍為106.9～118.7 m。當(dāng)南五采區(qū)工作面采高為6 m 時(shí)，最大導(dǎo)水裂隙帶高度為103.1～113.9 m。

4 結(jié)論

（1）依據(jù)相似礦井地質(zhì)條件分析，利用經(jīng)驗(yàn)公式法確定余吾礦南五采區(qū)工作面以6 m 厚度進(jìn)行開采時(shí)，導(dǎo)水裂隙帶發(fā)育最大高度為82.58 m。

（2）不同導(dǎo)水裂隙帶發(fā)育高度預(yù)測(cè)結(jié)果不同，計(jì)算機(jī)數(shù)值模擬采高6 m 時(shí)，南五采區(qū)導(dǎo)水裂隙帶最大發(fā)育高度為103 m，裂采比約17.2；室內(nèi)相似材料模擬預(yù)測(cè)的導(dǎo)水裂隙帶發(fā)育最大高度為105 m左右，裂采比為17.5。

（3）通過多種方法預(yù)測(cè)導(dǎo)水裂隙帶高度數(shù)據(jù)，結(jié)合實(shí)測(cè)數(shù)據(jù)，進(jìn)行數(shù)據(jù)擬合分析，采用綜合數(shù)據(jù)擬合公式預(yù)測(cè)的導(dǎo)水裂隙帶高度為采高為6 m時(shí)，南五采區(qū)最大導(dǎo)水裂隙帶高度為103.1～113.9 m。