馮錦艷，劉旭杭，于志全

(北京航空航天大學(xué) 交通科學(xué)與工程學(xué)院，北京 100191)

大傾角煤層采動(dòng)裂隙演化規(guī)律

馮錦艷，劉旭杭，于志全

(北京航空航天大學(xué) 交通科學(xué)與工程學(xué)院，北京 100191)

通過離散元UDEC數(shù)值模擬計(jì)算，分析了大傾角煤層開采過程中，擾動(dòng)巖體裂隙的演化規(guī)律，同時(shí)采用相似材料模擬實(shí)驗(yàn)結(jié)果驗(yàn)證了數(shù)值計(jì)算結(jié)果的可靠性。大傾角煤層采動(dòng)巖體裂隙發(fā)展初期以平行于工作面傾向的離層裂隙為主，后期離層裂隙的高度基本恒定，豎向裂隙擴(kuò)展加速，直至延伸到地表，形成錯(cuò)斷式裂縫和塌陷坑。在此基礎(chǔ)上，建立了大傾角煤層采動(dòng)裂隙分形維數(shù)D和采深的關(guān)系式，以及逾滲概率p和采深的關(guān)系式,實(shí)現(xiàn)了定量描述裂隙的演化規(guī)律。通過擬合分形維數(shù)D和逾滲概率p的關(guān)系式，可求得分形維數(shù)的峰值，進(jìn)而可求得對(duì)應(yīng)的采深，此采深即為破碎覆巖壓密沿傾斜方向重新形成穩(wěn)定結(jié)構(gòu)所需的采深。計(jì)算結(jié)果可用于預(yù)先估算大傾角煤層開采埋深范圍內(nèi)是否可形成基本頂穩(wěn)定結(jié)構(gòu)，進(jìn)而可采取不同的措施指導(dǎo)工作面的安全開采以及巷道的合理支護(hù)。

大傾角煤層;裂隙;分形維數(shù);逾滲概率;穩(wěn)定結(jié)構(gòu)

煤層傾角是確定采煤方法的重要因素之一，隨著綜合機(jī)械化采煤技術(shù)的不斷進(jìn)步，當(dāng)綜采技術(shù)逐漸向傾斜煤層和急傾斜煤層開采拓展時(shí)，開始出現(xiàn)“大傾角煤層”稱謂，中國(guó)按傾角的三類劃分方法已不能滿足開采實(shí)踐與學(xué)術(shù)發(fā)展的要求。目前，多數(shù)文獻(xiàn)認(rèn)為“大傾角煤層”是指傾角為35°～55°的煤層，是傾斜和急斜煤層的一部分[1-3]。國(guó)際主要采煤國(guó)家(前蘇聯(lián)、美國(guó)、德國(guó)、波蘭等)將傾角介于35°～55°的煤層稱為大斜或陡斜煤層，這與中國(guó)的大傾角煤層稱謂一致。

大傾角煤層約占中國(guó)煤炭探明儲(chǔ)量的20%以及產(chǎn)量的10%，受特殊成煤環(huán)境的控制，50%以上的大傾角煤層為無(wú)煙煤和優(yōu)質(zhì)焦煤，是中國(guó)保護(hù)性開采的稀缺煤種。在經(jīng)濟(jì)發(fā)展相對(duì)滯后的西部地區(qū)，煤炭資源開發(fā)已成為西部地區(qū)區(qū)域經(jīng)濟(jì)發(fā)展的重要支柱，而西部礦區(qū)50%以上的礦井存在大傾角煤層，大傾角煤層占煤炭總儲(chǔ)量的30%左右[3]，如主要產(chǎn)煤省(區(qū))的云南、貴州、四川、重慶、新疆、寧夏、甘肅等，要提高煤炭產(chǎn)業(yè)技術(shù)的經(jīng)濟(jì)效益，實(shí)現(xiàn)安全開采，就必須解決大傾角煤層開采中遇到的各種難題。

在經(jīng)濟(jì)快速發(fā)展的中國(guó)東部地區(qū)，大強(qiáng)度開采方式已經(jīng)使許多礦區(qū)的淺部資源(儲(chǔ)量)接近枯竭，如山東充州礦區(qū)、河北邢臺(tái)礦區(qū)和開灤礦區(qū)、安徽淮南和淮北礦區(qū)、江蘇徐州礦區(qū)、上海大屯礦區(qū)等，礦區(qū)必須轉(zhuǎn)向復(fù)雜的大傾角或急傾斜煤層開采，要保持這些礦區(qū)高產(chǎn)高效和可持續(xù)發(fā)展同樣需要解決大傾角煤層開采中遇到的各種問題。在老工業(yè)基地東北地區(qū)，鶴崗礦區(qū)一直致力于大傾角(急傾斜)煤層開采，雙鴨山礦區(qū)的主力礦井也面臨著條件異常復(fù)雜的大傾角煤層開采問題，要振興這些區(qū)域的工業(yè)經(jīng)濟(jì)，就必須振興其支柱產(chǎn)業(yè)—煤炭企業(yè)的經(jīng)濟(jì)發(fā)展。由此可見，大傾角煤層的安全開采問題是資源開發(fā)中的重大工程問題。

大傾角煤層是國(guó)際采礦界公認(rèn)的難采煤層，是世界范圍內(nèi)綜合機(jī)械化開采的“技術(shù)禁區(qū)”。早在20世紀(jì)七八十年代，前蘇聯(lián)、德國(guó)、法國(guó)、英國(guó)和波蘭開展了相關(guān)研究，西班牙和捷克等也做了一些工作，但其綜合機(jī)械化開采基礎(chǔ)理論、核心技術(shù)和關(guān)鍵裝備研究等一直未取得突破。近20多年來，中國(guó)歷經(jīng)幾代科研和技術(shù)工作者的努力，如石平五[4]、章黎明[5]、楊印朝[6]、解盤石[7]、王紅偉[8]、王金安[9-10]等許多科研人員的努力解決了大傾角煤層走向長(zhǎng)壁綜合機(jī)械化開采巖層控制基礎(chǔ)理論問題，攻克了綜合機(jī)械化開采關(guān)鍵技術(shù)，自主研發(fā)、制造了大傾角煤層綜合機(jī)械化開采裝備，建立了完善的工作面安全防護(hù)體系，形成了大傾角煤層走向長(zhǎng)壁綜合機(jī)械化開采理論和技術(shù)體系，這對(duì)中國(guó)乃至全世界的大傾角煤層開采具有重要意義。

大傾角煤層實(shí)現(xiàn)綜合機(jī)械化開采的同時(shí)，受到許多安全問題的困擾，如瓦斯泄漏、底板突水、地表水入滲、地面塌陷和地表裂縫等，這些安全問題均和擾動(dòng)巖體的裂隙演化規(guī)律密切相關(guān)。煤巖體中的裂隙場(chǎng)是流體滲透和運(yùn)移的主要通道和富集區(qū)，如何定量或定性描述煤巖體裂隙網(wǎng)絡(luò)的演化規(guī)律，一直是巖體力學(xué)領(lǐng)域關(guān)切的重大問題之一。

20世紀(jì)七八十年代，MNADELBROT[11]創(chuàng)立了分形幾何學(xué)，利用分形維數(shù)描述自然界不規(guī)則及無(wú)序的現(xiàn)象和行為。AVILES[12]，XIE[13-14]，BABADAGLI[15]等的研究表明，巖體中的采動(dòng)裂隙分布具有分形特征。中國(guó)的研究人員謝和平[16]、YU[17]、周宏偉[18]、王志國(guó)[19]等同樣揭示了采動(dòng)巖體裂隙網(wǎng)絡(luò)具有分形特性，并給出了采動(dòng)覆巖裂隙網(wǎng)絡(luò)分形維數(shù)與采寬間的統(tǒng)計(jì)關(guān)系。除此之外周福軍等[20]建立了二維、三維計(jì)算機(jī)裂隙網(wǎng)絡(luò)模型，并且研究了巖體不連續(xù)面分布的分形特征;王國(guó)艷等[21]利用RFPA 軟件揭示了采動(dòng)巖體的裂隙網(wǎng)絡(luò)分形特征。

除了分形方法外，逾滲理論(Percolation)也被引入用于描述覆巖裂隙的演化規(guī)律，它是處理強(qiáng)無(wú)序隨機(jī)幾何結(jié)構(gòu)的重要方法之一，從本質(zhì)上講屬于概率論的一個(gè)分支，概念最早是由 BROADBENT 等[22]在 1957 年提出，用于研究流體在隨機(jī)介質(zhì)中的滲流問題。馮增朝等[23-24]建立了煤體孔隙-裂隙的逾滲研究方法，并且研究了低滲透煤層的瓦斯?jié)B透性;栗東平等[25]分析了水平煤層采動(dòng)煤巖體裂隙網(wǎng)絡(luò)的演化特征，得到了逾滲概率和分形維數(shù)的變化關(guān)系。

可見，國(guó)內(nèi)外學(xué)者對(duì)裂隙網(wǎng)絡(luò)的分形及逾滲特征做了大量研究，揭示了采用逾滲方法和分形理論研究開采巖層移動(dòng)和裂隙演化規(guī)律的可行性和科學(xué)性，但對(duì)大傾角煤層開采條件下的采動(dòng)巖體裂隙網(wǎng)絡(luò)逾滲和分形特征以及兩者間的關(guān)系研究還不多見。鑒于此，本文結(jié)合某45°大傾角煤層工作面，開展數(shù)值模擬研究，對(duì)開采擾動(dòng)煤巖體的裂隙網(wǎng)絡(luò)逾滲與分形特征進(jìn)行研究，以期能夠定量描述采動(dòng)巖體裂隙的發(fā)展規(guī)律。

1 工程概況

以甘肅省某礦井煤層工作面為例，煤層傾角45°，地表為12 m厚的黃土，有1煤和2煤兩層主采煤層，主采煤層水平間距60 m，其間為粗砂巖細(xì)砂巖互層，1煤直接頂為8 m厚的粉砂巖，下有0.5 m厚的高炭質(zhì)泥巖偽頂，極易垮落，基本頂為中粗細(xì)砂巖，厚度大于25 m，直接底為0.8 m厚的泥巖，基本底為18 m厚的粗砂巖。

礦井采用走向長(zhǎng)壁綜采放頂煤采煤法從地下埋深-150 m水平開采，一次垂直采高50 m，每個(gè)水平留有垂直高15 m的煤柱，工作面走向長(zhǎng)400 m。

2 數(shù)值計(jì)算模型及參數(shù)

根據(jù)所選區(qū)域垂向剖面的地形特征、地質(zhì)條件以及巖體結(jié)構(gòu)特點(diǎn)，采用UDEC2D進(jìn)行數(shù)值模擬計(jì)算，地質(zhì)模型如圖1所示，模型尺寸長(zhǎng)×寬=1 300 m×495 m。模型兩側(cè)設(shè)置水平約束，模型底側(cè)固定、頂面自由?；夭砷_挖過程為：在階段水平內(nèi)采用自左向右的扒皮式回采，階段水平間采用自上而下的下行式一次采全高垮落式放頂煤開采模式，主采煤層的開采順序如圖2所示，為了更好地探索采動(dòng)巖體的裂隙。數(shù)值計(jì)算采用的參數(shù)見表1，2。

圖1 垂向剖面地質(zhì)模型示意Fig.1 Geological model diagram of vertical section

圖2 主采煤層開采順序示意Fig.2 Sketch map of step-by-step excavation process in main mineable coal seam

巖土層密度d/(kg·m-3)體積模量K/GPa剪切模量G/GPa抗拉強(qiáng)度σT/MPa黏聚力c/MPa內(nèi)摩擦角φ/(°)地表黃土層140003020001002231煤層直接頂以上巖層240346281500300341煤層直接頂巖層240040201000200301與2煤層間巖層250060401500250321與2煤層150065201100180302煤層底板巖層240810075350085043風(fēng)化帶巖層26006515120022032

表2斷層與節(jié)理力學(xué)參數(shù)
Table2Mechanicsparametersoffaultsandjoints

巖土層體積模量K/GPa剪切模量G/GPa抗拉強(qiáng)度σT/MPa黏聚力c/MPa內(nèi)摩擦角φ/(°)地表黃土層03020001002201煤層直接頂以上巖層60400100120151煤層直接頂巖層40200020101與2煤層間巖層40300600140281與2煤層60401500250322煤層底板巖層100100250065040風(fēng)化帶巖層7040100016030

3 采動(dòng)巖體裂隙計(jì)算結(jié)果分析

圖3給出了由數(shù)值模擬計(jì)算得到的開采擾動(dòng)巖體裂隙演化圖。從圖3可以看出，當(dāng)開采埋深較小時(shí)(-200 m)，上覆巖層主要產(chǎn)生平行于工作面傾向的離層裂隙，裂隙帶高度為90 m左右。隨著大傾角煤層開采深度的增加，上覆巖層裂隙帶高度逐漸增加，當(dāng)開采埋深為-320 m時(shí)，離層裂隙高度為170 m左右;當(dāng)開采埋深為-400 m時(shí)，離層裂隙帶高度為220 m左右，此時(shí)離層裂隙貫穿至地表。由此可見，大傾角煤層采動(dòng)巖體裂隙包含平行于工作面的離層裂隙和垂直于工作面的豎向裂隙兩種，當(dāng)離層裂隙垂直于工作面發(fā)展到一定高度后，可與豎向裂隙結(jié)合貫穿至地表。

圖3 采動(dòng)巖體裂隙Fig.3 Fractures evolution of mining rocks

計(jì)算結(jié)果同時(shí)揭示出，多個(gè)主采煤層之間的裂隙擴(kuò)展速度小于工作面上覆擾動(dòng)巖體的裂隙擴(kuò)展速度。

為了驗(yàn)證數(shù)值計(jì)算結(jié)果的準(zhǔn)確性，同時(shí)進(jìn)行了相似材料模擬實(shí)驗(yàn)。試驗(yàn)結(jié)果顯示，隨著開采深度的增加，裂隙逐漸擴(kuò)張，當(dāng)開采到埋深-400 m時(shí)，采動(dòng)巖體離層裂隙快速發(fā)展，直至貫通地表，誘發(fā)大范圍的錯(cuò)斷和塌陷(圖4)，這與數(shù)值計(jì)算結(jié)果基本一致，說明了數(shù)值計(jì)算結(jié)果的可靠性，可在此基礎(chǔ)上進(jìn)行進(jìn)一步分析。

圖4 覆巖裂隙貫通至地表Fig.4 Fractures were through the earth’s surface

4 采動(dòng)裂隙演化的分形計(jì)算

分形維數(shù)是定量表征事件分形性質(zhì)的重要參數(shù)，這對(duì)于定量描述采動(dòng)裂隙的演化規(guī)律具有重要意義。本文采用改變粗視化程度的方法求取分形維數(shù)，該方法適用于具有統(tǒng)計(jì)意義的自相似事件求取分形維數(shù)，Mandelbrot提出的方格統(tǒng)計(jì)法就是該法的一例。用不同尺度r的方格網(wǎng)覆蓋所研究的覆巖裂隙，計(jì)數(shù)每一次覆蓋時(shí)圍巖裂隙所占有的方格數(shù)N(r)，兩者之間的關(guān)系為

將這種關(guān)系表示在雙對(duì)數(shù)坐標(biāo)系中，可以得到lgN(r)～lgr關(guān)系曲線，在雙對(duì)數(shù)圖中呈簡(jiǎn)單的直線，其中斜率正好為-D。

圖5給出了分形研究中的覆巖裂隙網(wǎng)格劃分示意圖，根據(jù)相似材料模擬實(shí)驗(yàn)結(jié)果，覆巖和底板的破壞影響范圍為長(zhǎng)×寬=500 m×495 m，在此基礎(chǔ)上，r采取了5個(gè)數(shù)量，分別為250，125，62.5，31.25和15.625 m。根據(jù)不同采深情況下的采動(dòng)巖體裂隙lgN(r)～lgr曲線，可以得出分形維數(shù)D，斷層和節(jié)理均按照裂隙處理。

圖5 采動(dòng)巖體分形幾何格子覆蓋法示意Fig.5 Lattice layout of mining rocks by fractal theory

將裂隙的分形維數(shù)與采深的關(guān)系繪制如圖6所示，可以看出當(dāng)采深小于400 m時(shí)，采動(dòng)巖體裂隙的分形維數(shù)隨開采深度的增加而增加，且呈線性增加趨勢(shì)，擬合得到擾動(dòng)巖體裂隙分形維數(shù)與埋深H的關(guān)系滿足式(2):

圖6 采動(dòng)巖體裂隙分形維數(shù)與采深關(guān)系曲線Fig.6 Relationship between the fractal dimension and mining depth of mining rocks

5 采動(dòng)裂隙演化的逾滲概率計(jì)算

為了更好地定量描述采動(dòng)巖體裂隙的演化規(guī)律，為安全開采提供理論依據(jù)，同時(shí)采用逾滲方法計(jì)算了采動(dòng)巖體的裂隙發(fā)展規(guī)律。

座逾滲也稱為格點(diǎn)逾滲，是逾滲的兩種基本模型之一，主要考慮的是網(wǎng)格中各點(diǎn)是否被占據(jù)的問題，每個(gè)格點(diǎn)被占據(jù)和不被占據(jù)的概率分別為p和1-p，同樣每個(gè)格點(diǎn)占據(jù)概率與周圍相鄰格點(diǎn)的狀態(tài)無(wú)關(guān)。對(duì)于煤巖體而言，整體可看作被裂隙和固體顆粒完全充滿，由煤層開采引起的平行于煤層傾向的離層裂隙和豎向破斷裂隙貫通，構(gòu)成了許多連通的裂隙團(tuán)，簡(jiǎn)稱為團(tuán)。隨著開采深度的增加，裂隙逐步擴(kuò)展且貫通，相鄰的團(tuán)形成更大的團(tuán)，其中占格點(diǎn)最大的團(tuán)稱為最大團(tuán)。隨著最大團(tuán)所包含的裂隙數(shù)逐漸增加，當(dāng)其連通平面區(qū)域的兩個(gè)對(duì)稱邊界時(shí)稱為逾滲團(tuán)，逾滲概率可定義為

其中，p為逾滲概率;M(S)為最大團(tuán)中被占格點(diǎn)的數(shù)量;Sn為網(wǎng)格中格點(diǎn)的總數(shù)量。當(dāng)為二維正方形網(wǎng)絡(luò)時(shí)，S為每行(列)中所包含的網(wǎng)格數(shù)，此時(shí)n=2。

通過編制計(jì)算機(jī)程序，對(duì)數(shù)值模擬計(jì)算得到的裂隙網(wǎng)絡(luò)二值圖確定各個(gè)逾滲團(tuán)，進(jìn)而可確定最大團(tuán)，屏蔽除最大團(tuán)以外的逾滲團(tuán)及孤立的裂隙，將僅含最大團(tuán)的裂隙網(wǎng)絡(luò)圖轉(zhuǎn)化為[0，1]矩陣，含有裂隙的格點(diǎn)轉(zhuǎn)化為1，不含裂隙的格點(diǎn)轉(zhuǎn)化為0;0的位置被固體顆粒占據(jù)，為不滲透格點(diǎn)。統(tǒng)計(jì)矩陣中1的個(gè)數(shù)，即為最大團(tuán)被占格點(diǎn)的數(shù)量M(S)，將其除以裂隙網(wǎng)絡(luò)中格點(diǎn)總數(shù)即可得到逾滲概率p。工作面不同開采深度的裂隙圖對(duì)應(yīng)的逾滲概率p與開采埋深H之間的關(guān)系如圖7所示，擬合公式為

從計(jì)算結(jié)果可以看出，隨著大傾角煤層開采深度的增加，采動(dòng)巖體裂隙的逾滲概率呈拋物線式增加趨勢(shì)，當(dāng)埋深小于330 m時(shí)，逾滲概率增加趨勢(shì)稍緩，基本呈直線關(guān)系;當(dāng)埋深大于330 m時(shí)，逾滲概率增加幅度加劇，說明采動(dòng)裂隙加劇了連通速度。

圖7 采動(dòng)巖體裂隙逾滲概率與采深關(guān)系曲線Fig.7 Relationship between the percolation probability and mining depth of mining rocks

6 采動(dòng)裂隙的分形維數(shù)與逾滲概率的關(guān)系

由于各礦區(qū)埋深條件不同，如果能消除埋深的影響，建立采動(dòng)巖體裂隙的分形維數(shù)與逾滲概率之間的關(guān)系，對(duì)于推廣研究煤巖體的裂隙定量描述將具有重要意義。

由式(2)和(4)可以看出，采動(dòng)巖體裂隙的分形維數(shù)與逾滲概率之間存在關(guān)系，通過聯(lián)立公式可以消除開采深度的影響，有利于定量描述開采埋深小于400 m的煤巖體裂隙演化規(guī)律。經(jīng)過推導(dǎo)，由式(2)和(4)推導(dǎo)出采動(dòng)巖體裂隙分形維數(shù)D與逾滲概率p之間的關(guān)系為

式(5)的圖形表示如圖8所示，采動(dòng)巖體分形維數(shù)與逾滲概率之間存在拋物線的關(guān)系式。根據(jù)公式推導(dǎo)預(yù)測(cè)，當(dāng)采深為421 m時(shí)，分形維數(shù)達(dá)到峰值1.76;當(dāng)采深小于421 m時(shí)，分形維數(shù)D隨逾滲概率的增加而增加，即本工作面開采深度范圍內(nèi)，裂隙的分形維數(shù)一直呈增加的趨勢(shì);當(dāng)采深大于421 m時(shí)，分形維數(shù)D隨逾滲概率的增加而減小。計(jì)算結(jié)果充分說明了隨著開采深度的增加，上部覆巖重新壓密，形成穩(wěn)定結(jié)構(gòu)。

圖8 采動(dòng)巖體裂隙逾滲概率與分形維數(shù)關(guān)系曲線Fig.8 Relationship between the percolation probability and fractal dimension of mining rocks

圖9給出了采動(dòng)巖體裂隙閉合率隨采深的變化曲線，可以看出，裂隙閉合率隨著開采深度的增加也在增加，當(dāng)采深在400 m左右時(shí)，裂隙閉合率呈現(xiàn)恒定值趨勢(shì)，即閉合裂隙和總體裂隙都在增加，但其比率保持恒定，驗(yàn)證了隨著開采深度的增加，大傾角煤層采動(dòng)巖體重新壓密，沿傾斜方向可形成新的穩(wěn)定結(jié)構(gòu)[26](圖10)。

圖9 采動(dòng)巖體裂隙閉合率與埋深的關(guān)系曲線Fig.9 Relationship between the crack closure rate and mining depth of mining rocks

圖10 采動(dòng)巖體愈合的力學(xué)機(jī)制Fig.10 Healing mechanics mechanism of mining rocks

通過以上分析可知，當(dāng)單一采用分形維數(shù)D或者單一采用逾滲概率p來描述采動(dòng)巖體裂隙演化規(guī)律時(shí)，在開采范圍內(nèi)均呈增加的變化趨勢(shì)，較難尋找到峰值點(diǎn)。為了更好地定量描述采動(dòng)巖體裂隙的發(fā)展規(guī)律，尋找到分形維數(shù)的峰值點(diǎn)，可通過擬合逾滲概率p和分形維數(shù)D的關(guān)系式獲得，依二者的關(guān)系可較容易地確定分形維數(shù)D峰值對(duì)應(yīng)的開采埋深，即為基本頂沿傾斜方向重新形成穩(wěn)定結(jié)構(gòu)需要的采深。

根據(jù)數(shù)值模擬計(jì)算或者相似材料模擬試驗(yàn)結(jié)果可知，當(dāng)在采深范圍內(nèi)基本頂可以形成承載結(jié)構(gòu)時(shí)，可有效利用其穩(wěn)定性實(shí)現(xiàn)安全開采，當(dāng)基本頂在采深范圍內(nèi)尚未形成承載結(jié)構(gòu)時(shí)，則需加強(qiáng)工作面和巷道的支護(hù)措施，確保安全生產(chǎn)。這對(duì)大傾角煤層采動(dòng)巖體裂隙的演化規(guī)律研究以及安全開采提供了理論依據(jù)。

7 閉合裂隙的導(dǎo)水性探討

為了充分探討大傾角煤層采動(dòng)巖體裂隙的愈合機(jī)制，對(duì)擾動(dòng)區(qū)域進(jìn)行了滲流計(jì)算，滲流計(jì)算參數(shù)見表3。

表3水文計(jì)算參數(shù)
Table3Hydrologicalcalculationparameters

巖土層零法向應(yīng)力張開度/mm殘余應(yīng)力張開度/mm節(jié)理滲透系數(shù)/(Pa·s)-1地表黃土層25073571煤層直接頂以上巖層1507251煤層直接頂巖層20073001與2煤層間巖層1007101與2煤層18072002煤層底板巖層1007150風(fēng)化帶巖層1507100

計(jì)算結(jié)果如圖11,12所示，根據(jù)數(shù)值計(jì)算得到的工作面滲流矢量場(chǎng)，可計(jì)算工作面的總流量。計(jì)算結(jié)果顯示，工作面內(nèi)的流量隨開采深度的增加先呈增加的趨勢(shì)，當(dāng)埋深在400 m時(shí)，流量出現(xiàn)了略微減小的趨勢(shì)，驗(yàn)證了閉合裂隙具有了一定程度的阻水作用。

圖11 開采到地下-320 m的巖體滲流矢量場(chǎng)Fig.11 Total water inflow when mining to the burial depth of -320 m

圖12 工作面涌水量與采深的關(guān)系Fig.12 Relationship between total water inflow of the working faces and mining depth

8 結(jié) 論

(1)大傾角煤層綜放開采，上部覆巖以離層裂隙為主，隨著開采深度的增加，離層裂隙的高度逐漸增加，直至恒定;后期開采會(huì)進(jìn)一步導(dǎo)致豎向裂隙加速擴(kuò)展，直至連通地表，出現(xiàn)錯(cuò)斷式裂縫和塌陷坑。

(2)在一定采深范圍內(nèi)，大傾角煤層采動(dòng)巖體裂隙的分形維數(shù)D以及逾滲概率p均隨采深的增加逐漸增加，說明了裂隙的不斷擴(kuò)展趨勢(shì)。

(3)采動(dòng)巖體裂隙的分形維數(shù)D與逾滲概率p之間的關(guān)系是描述裂隙演化規(guī)律的有效方法。通過擬合逾滲概率p和分形維數(shù)D的關(guān)系式獲得，依二者的關(guān)系可較容易地確定分形維數(shù)D峰值對(duì)應(yīng)的開采埋深，即為基本頂沿傾斜方向重新形成穩(wěn)定結(jié)構(gòu)需要的采深，當(dāng)小于此埋深時(shí)，需加強(qiáng)工作面的支護(hù)，當(dāng)大于此埋深時(shí)，可充分利用基本頂沿傾斜方向的穩(wěn)定性結(jié)構(gòu)，適當(dāng)減弱支護(hù)結(jié)構(gòu)，可為現(xiàn)場(chǎng)工作面安全開采以及巷道合理支護(hù)提供理論依據(jù)。

[1] 伍永平，劉孔智，贠東風(fēng)，等.大傾角煤層安全高效開采技術(shù)研究進(jìn)展[J].煤炭學(xué)報(bào)，2014，39(8)：1611-1618. WU Yongping，LIU Kongzhi，YUN Dongfeng，et al.Research progress on the safe and efficient mining technology of steeply dipping seam[J].Journal of China Coal Society，2014，39(8)：1611-1618.

[2] 伍永平，贠東風(fēng)，張淼豐.大傾角煤層綜采基本問題研究[J].煤炭學(xué)報(bào)，2000，25(5)：465-468. WU Yongping，YUN Dongfeng，ZHANG Miaofeng.Study on the elementary problems of full-mechanized coal mining in greater pitching seam[J].Journal of China Coal Society，2000，25(5)：465-468.

[3] 王紅偉.大傾角煤層開采覆巖結(jié)構(gòu)特征分析[D].西安：西安科技大學(xué)，2010. WANG Hongwei.Analysis of overburden structure feature in steeply dipping seam mining[D].Xi’an：Xi’an University of Science and Technology,2010.

[4] 石平五，高召寧.急斜特厚煤層開采圍巖與覆蓋層破壞規(guī)律[J].煤炭學(xué)報(bào)，2003，28(1)：13-16. SHI Pingwu，GAO Zhaoning.The failure laws of surrounding rocks and overlying bed in the steep special thickness seam mining[J].Journal of China Coal Society，2003，28(1)：13-16.

[5] 章黎明.王家山煤礦大傾角厚煤層綜放采場(chǎng)礦壓規(guī)律研究[J].煤炭科學(xué)技術(shù)，2007，35(12)：22-26. ZHANG Liming.Research on ground behavior law of fully mechanized caving mining face in deep inclined seam of Wangjiashan Mine[J].Coal Science and Technology.2007，35(12)：22-26.

[6] 楊印朝，徐文彬，潘衛(wèi)東，等.大傾角厚煤層綜采工作面設(shè)備改造與生產(chǎn)管理[J].現(xiàn)代礦業(yè)，2015(11)：234-237. YANG Yinzhao，XU Wenbin，PAN Weidong,et al.The equipment modification and production management of full-mechanized coal mining in greater pitching seam[J].Modern Mining，2015(11)：234-237.

[7] 解盤石，伍永平.大傾角煤層長(zhǎng)壁大采高開采煤壁片幫機(jī)理及防控技術(shù)[J].煤炭工程，2015，47(1)：74-77. XIE Panshi，WU Yongping.Mechanism and control methods of rib spalling in steeply dipping thick seam in fully-mechanized longwall mining with large mining height[J].Coal Engineering，2015，47(1)：74-77.

[8] 王紅偉，伍永平，曹沛沛，等.大傾角煤層開采大型三維可加載相似模擬試驗(yàn)[J].煤炭學(xué)報(bào)，2015，40(7)：1505-1511. WANG Hongwei，WU Yongping，CAO Peipei，et al.Large scale loadable 3D-simulation tests on mining steeply dipping seam[J].Journal of China Coal Society，2015，40(7)：1505-1511.

[9] 王金安，張基偉，高小明，等.大傾角厚煤層長(zhǎng)壁綜放開采基本頂破斷模式及演化過程(I)——初次破斷[J].煤炭學(xué)報(bào)，2015，40(6)：1353-1360. WANG Jin’an，ZHANG Jiwei，GAO Xiaoming，et al.Fracture mode and evolution of main roof stratum above longwall fully mechanized top coal caving in steeply inclined thick coal seam (I)：initial fracture[J].Journal of China Coal Society，2015，40(6)：1353-1360.

[10] 王金安，張基偉，高小明，等.大傾角厚煤層長(zhǎng)壁綜放開采基本頂破斷模式及演化過程(II)——周期破斷[J].煤炭學(xué)報(bào)，2015，40(8)：1737-1745. WANG Jin’an，ZHANG Jiwei，GAO Xiaoming，et al.Fracture mode and evolution of main roof stratum above longwall fully mechanized top coal caving in steeply inclined thick coal seam (II):Periodie fracture[J].Journal of China Coal Society，2015，40(8)：1737-1745.

[11] MANDELBROT B B.The fractal geometry of nature[M].New York:W.H.Freeman,1982.

[12] AVILES C A,SCHOLZ C H,BOATWRIGHT J.Fractal analysis applied to characteristic segments of the San Andreas fault[J].Journal of Geophysical Research,1987,92(B1):331-334.

[13] XIE H.Fractals in rock mechanics[M].Netherlands:Balkema A A Publishers,1993.

[14] XIE H,SANDERSON D J,PEACOCK D C P.A fractal model and energy dissipation for an echelon fractures[J].Engineering Fracture Mechanics,1994,48(5):655-662.

[15] BABADAGLI Tayfun.Fractal analysis of 2D fracture networks of geothermal reservoirs in south-western Turkey[J].Journal of Volcanology and Geothermal Research,2001,11(2):83-103.

[16] 謝和平，于廣明，楊倫，等.采動(dòng)巖體分形裂隙網(wǎng)絡(luò)研究[J].巖石力學(xué)與工程學(xué)報(bào)，1999，18(2)：147-151. XIE Heping,YU Guangming,YANG Lun,et al.Research on the fractal effects of crack network in overburden rock stratum[J].Chinese Journal of Rock Mechanics and Engineering,1999,18(2):147-151.

[17] YU G M,XIE H,ZHAO J F,et al.Fractal evolution of a crack network in overburden rock stratum[J].Discrete Dynamics in Nature and Society,2000,5(1):47-52.

[18] 周宏偉，張濤，薛東杰，等.長(zhǎng)壁工作面覆巖采動(dòng)裂隙網(wǎng)絡(luò)演化特征[J].煤炭學(xué)報(bào),2011,36(12):1957-1962. ZHOU Hongwei,ZHANG Tao,XUE Dongjie,et al.Evolution of mining-induced crack network in overburden strata of longwall face[J].Journal of China Coal Society,2011,36(12):1957-1962.

[19] 王志國(guó)，周宏偉，謝和平.深部開采上覆巖層采動(dòng)裂隙網(wǎng)絡(luò)演化的分形特征研究[J].巖土力學(xué),2009,30(8):2403-2408. WANG Zhiguo,ZHOU Hongwei,XIE Heping.Research on fractal characterization of mined crack network evolution in overburden rock stratum under deep mining[J].Rock and Soil Mechanics，2009，30(8):2403-2408.

[20] 周福軍，陳劍平，牛岑岑.裂隙化巖體不連續(xù)面密度的分形研究[J].巖石力學(xué)與工程學(xué)報(bào)，2013，32(S1)：2624-2631. ZHOU Fujun,CHEN Jianping,NIU Cencen.Study of discontinuity density of fractured rock masses based on fractal theory[J].Chinese Journal of Rock Mechanics and Engineering，2013，32(S1):2624-2631.

[21] 王國(guó)艷，于廣明，于永江，等.采動(dòng)巖體裂隙分維演化規(guī)律分析[J].采礦與安全工程學(xué)報(bào)，2012，29(6)：859-863. WANG Guoyan,YU Guangming,YU Yongjiang,et al.Study on cracks fractal evolution laws of mining rock mass[J].Journal of Mining & Safety Engineering,2012,29(6):859-863.

[22] BROADBENT S K,HAMMERSLEY J M.Percolation Processes[J].Mathematical Proceedings of the Cambridge Philosophical Society,1957，53(3):629-641.

[23] 馮增朝，趙陽(yáng)升，文再明.煤巖體孔隙裂隙雙重介質(zhì)逾滲機(jī)理研究[J].巖石力學(xué)與工程學(xué)報(bào),2005，24(2)：236-240. FENG Zengchao,ZHAO Yangsheng,WEN Zaiming.Percolation mechanism of fractured coal rocks as dual-continua[J].Chinese Journal of Rock Mechanics and Engineering，2005，24(2)：236-240.

[24] 馮增朝,趙陽(yáng)升,呂兆興.強(qiáng)隨機(jī)分布裂隙介質(zhì)的二維逾滲規(guī)律研究[J].巖石力學(xué)與工程學(xué)報(bào),2006,25(S2):3904-3908. FENG Zengchao,ZHAO Yangsheng,Lü Zhaoxing.Research on laws of 2D percolation of fully random distribution fracture media[J].Chinese Journal of Rock Mechanics and Engineering，2006，25(S2)：3904-3908.

[25] 栗東平，周宏偉，薛東杰，等.煤巖體采動(dòng)裂隙網(wǎng)絡(luò)的逾滲與分形特征關(guān)系研究[J].巖土力學(xué)，2015，36(4)：1135-1140. LI Dongping,ZHOU Hongwei,XUE Dongjie,et al.Relationship between percolation and fractal properties of mining-induced crack network in coal and rock masses[J].Rock and Soil Mechanics，2015，36(4)：1135-1140.

[26] 馮錦艷，王金安，朱建明，等.急傾斜長(zhǎng)壁工作面開采老頂平衡結(jié)構(gòu)阻水性能[J].煤炭學(xué)報(bào)，2009,34(2):156-158. FENG Jinyan,WANG Jin’an,ZHU Jianming,et al.The effect of water-block for the bearing structure in main roof rocks due to mining in steep long-wall seam[J].Journal of China Coal Society，2009，34(2)：156-158.

Numericalsimulationstudyonthemining-inducedfractureevolutionofsteepcoalseam

FENG Jinyan，LIU Xuhang,YU Zhiquan

(SchoolofTransportationScienceandEngineering,BeihangUniversity,Beijing100191,China)

Using the discrete element numeric simulation program UDEC,the evolution law of mining-induced fractures in the mining process in steep coal seam was studied.Combined with the study of similar material simulation experiment results,the result by numerical calculation was verified.The results indicate that bed separation fissures are the main type at the initial stage,the height of the bed separation fissures remain constant and the vertical fissures speed up to extend in the later period.When the vertical fissures extend to the surface,ladder cracks and collapse pits may be formed.The relational expression between fractal dimension (D) of the mining-induced fracture and mining depth (H) was established,and the relational expression between percolation probability (p) of the mining-induced fracture and mining depth (H) was also established based on the numerical simulation,so that the evolution law of mining-induced fracture can be quantified.By fitting the relation between the fractal dimension (D) and penetration probability (p),the peak value of fractal dimension can be obtained and the corresponding mining depth,which is needed for re-forming the stable structure of broken rock masses along the dip direction,can also be obtained.The study can be used to estimate whether steep coal seam can form a stable main roof structure in the range of mining depth in advance,and then guide the rational support measures of roadway and the safe mining of working face.

steep coal seam;fracture;fractal dimension;percolation probability;stable structure

10.13225/j.cnki.jccs.2016.1859

TD823.21;TD311

:A

:0253-9993(2017)08-1971-08

國(guó)家自然科學(xué)基金資助項(xiàng)目(41302273)

馮錦艷(1978—)，女，天津人，副教授，碩士生導(dǎo)師。E-mail:fengjinyan226@163.com

馮錦艷，劉旭杭，于志全.大傾角煤層采動(dòng)裂隙演化規(guī)律[J].煤炭學(xué)報(bào),2017,42(8):1971-1978.

FENG Jinyan，LIU Xuhang,YU Zhiquan.Numerical simulation study on the mining-induced fracture evolution of steep coal seam[J].Journal of China Coal Society,2017,42(8):1971-1978.doi:10.13225/j.cnki.jccs.2016.1859