文章編號1000-5269（2024）06-0049-08

DOI：10.15958/j.cnki.gdxbzrb.2024.06.08

摘要：厘清采動(dòng)覆巖導(dǎo)水裂隙演化規(guī)律對頂板防治水研究具有重要意義。以龍鳳煤礦9#煤層120910綜采工作面為工程背景，通過相似模擬和數(shù)值模擬方法，對緩傾斜淺埋煤層覆巖裂隙演化規(guī)律進(jìn)行了研究。結(jié)果表明：煤層采動(dòng)頂板導(dǎo)水裂隙帶發(fā)育高度將貫穿長興組含水層底界，工作面推進(jìn)過程中存在突水危險(xiǎn)性，并將覆巖裂隙演化過程分為4個(gè)階段，即裂隙發(fā)育階段、頂板垮落階段、裂隙壓實(shí)階段、整體穩(wěn)定階段；同時(shí)引入分形理論，定量描述了采動(dòng)覆巖裂隙演化過程，裂隙演化4個(gè)階段的分形損傷度ω分別為0.676、0.702、0.638、0.639，當(dāng)ω≥0.68時(shí)，巖體破裂現(xiàn)象越明顯。研究結(jié)果可為相似地質(zhì)條件下上覆巖層頂板防治水研究提供參考。

關(guān)鍵詞：緩傾斜煤層；裂隙演化；分形損傷度；相似模擬；數(shù)值模擬

中圖分類號：TD82

文獻(xiàn)標(biāo)志碼：A

黔北煤田以淺埋煤層為主，受復(fù)雜地質(zhì)條件影響，礦區(qū)內(nèi)頂板冒落和突水突泥等礦山災(zāi)害頻發(fā)，致災(zāi)機(jī)理復(fù)雜［1-3］。地下開采活動(dòng)致使上覆巖層應(yīng)力重新分布，引起上覆巖層發(fā)生變形破壞并形成導(dǎo)水裂隙帶，從而引發(fā)地下水運(yùn)移和地表非連續(xù)變形等一系列工程問題，因此研究緩傾斜淺埋煤層上覆煤巖裂隙發(fā)育規(guī)律對煤礦安全高效開采具有重要作用。

目前，從采動(dòng)作用對巖體的損傷程度來看，采動(dòng)影響對巖體的裂隙擴(kuò)展程度隨著推進(jìn)距離、推進(jìn)速度不斷變化。因此，如何建立合理的采動(dòng)損傷評價(jià)指標(biāo)，客觀準(zhǔn)確衡量采動(dòng)損傷程度和采動(dòng)裂隙之間的對應(yīng)關(guān)系已成為采礦領(lǐng)域亟待解決的難題。近年來，隨著數(shù)字圖像處理技術(shù)的發(fā)展，通過得到工作面周圍裂隙的分布圖像，采用分形損傷度研究采場巖體裂隙損傷分形特征，為煤礦開采過程中巖體損傷及破壞程度定量表征提供新思路。國內(nèi)外學(xué)者對采動(dòng)裂隙拓展發(fā)育過程和形態(tài)特征進(jìn)行了大量研究，王創(chuàng)業(yè)等［4］利用分形理論將淺埋煤層工作面推進(jìn)過程概括為3個(gè)特征階段。ZHANG等［5］研究了淺埋煤層演化特征并推導(dǎo)了相應(yīng)結(jié)構(gòu)的力學(xué)計(jì)算模型，將開采工作面的裂隙演化和巖石破裂特征分為3個(gè)區(qū)域。文獻(xiàn)［6-7］揭示了淺埋煤層覆巖裂隙發(fā)育及演化規(guī)律。趙毅鑫等［8］采用通用離散單元法程序（universal discrete element code， UDEC）建立了采動(dòng)覆巖裂隙與耗散能量的關(guān)系，進(jìn)而得到導(dǎo)水裂隙帶高度與采高和工作面推進(jìn)距離的擬合公式。文獻(xiàn)［9-10］通過相似模擬試驗(yàn)對近距離煤層群重復(fù)采動(dòng)下覆巖裂隙發(fā)育規(guī)律進(jìn)行研究。文獻(xiàn)［11-13］通過相似模擬試驗(yàn)分析不同開采條件下的覆巖垮落形態(tài)及裂隙演化分布特點(diǎn)。文獻(xiàn)［14-16］通過分形維數(shù)對裂隙網(wǎng)絡(luò)滲透特性進(jìn)行研究，表明不同區(qū)域下裂隙網(wǎng)絡(luò)滲透率與分維值呈正相關(guān)性。郝鑫宇［17］通過三維電法分析不同開采狀態(tài)下覆巖電阻率變化特征，并對厚煤層覆巖裂隙演化特征進(jìn)行研究。GAO等［18］對放頂煤條件下裂隙演化規(guī)律與支架支承壓力變化趨勢進(jìn)行研究，建立了比較完整的定性（半定量）的分析模型。

綜上所述，眾多學(xué)者在覆巖運(yùn)移和裂隙發(fā)育規(guī)律等方面取得了有益成果，然而，針對黔北煤田緩傾斜淺埋煤層覆巖裂隙演化過程大多處于定性描述階段，未能定量化表征上覆巖層裂隙的分布特征和空間形態(tài)。因此，本文以黔北煤田龍鳳煤礦9#煤層120910為工程背景，通過相似模擬試驗(yàn)、UDEC數(shù)值模擬、分形理論等研究手段，對采動(dòng)影響下的覆巖裂隙分布形態(tài)和動(dòng)態(tài)演化過程進(jìn)行研究，并通過分形損傷度模型來定量表征開采擾動(dòng)下巖體破壞過程中的裂隙孕育及拓展規(guī)律，為黔北煤田煤礦安全高效開采提供指導(dǎo)及依據(jù)。

1工程概括

龍鳳煤礦位于金沙縣新華鄉(xiāng)，為年生產(chǎn)規(guī)模90萬t生產(chǎn)礦井，采用長臂走向后退式采煤法，當(dāng)前主采煤層為二疊系上統(tǒng)龍?zhí)督M中上部9#煤層，煤厚2.00～4.17 m，平均厚2.98 m。其中120910采煤工作面埋深介于180～200 m，煤層傾角7°～11°，屬于緩傾斜淺埋煤層。根據(jù)龍鳳煤礦綜合柱狀圖（圖1），結(jié)合關(guān)鍵層理論，判別出區(qū)內(nèi)主關(guān)鍵層為長興組石灰?guī)r，亞關(guān)鍵層為龍?zhí)督M泥質(zhì)粉砂巖。

2傾斜淺埋煤層覆巖裂隙演化相似模擬試驗(yàn)

2.1相似模擬試驗(yàn)?zāi)Ｐ?/p>

根據(jù)龍鳳煤礦120910工作面地質(zhì)資料，構(gòu)建二維物理相似模擬試驗(yàn)?zāi)Ｐ?。本次試?yàn)平臺(tái)幾何尺寸長3.0 m、寬0.25 m和高1.5 m，選取河沙作為骨料，石膏和石灰作為膠結(jié)物，采用云母粉在不同巖層間進(jìn)行分層，將龍鳳煤礦實(shí)際地質(zhì)條件按幾何相似比1∶100進(jìn)行開挖（下文描述距離單位已通過幾何相似比換算），相似模擬試驗(yàn)煤巖層材料配比參數(shù)見表1。

模型完成21 d養(yǎng)護(hù)后進(jìn)行開挖（圖2）。9#煤層左右兩側(cè)各留60 m保護(hù)煤柱，從左到右每10 m開挖一次，每隔5 s高速攝像機(jī)拍照采樣，記錄相似模型變形情況，直至工作面開挖完成后上覆巖層趨于穩(wěn)定階段時(shí)停止采集。

2.2試驗(yàn)結(jié)果及分析

由圖3可知，工作面覆巖垮落過程主要分為4個(gè)階段：

當(dāng)工作面推進(jìn)至30 m時(shí)，直接頂出現(xiàn)離層裂隙，有向下發(fā)生垮落的趨勢，處于裂隙發(fā)育階段。當(dāng)工作面推進(jìn)至90 m時(shí)，在采動(dòng)影響和上覆巖層彈性能釋放共同作用下，直接頂和基本頂同時(shí)垮落，導(dǎo)致兩側(cè)橫向破斷裂隙和縱向裂隙大量增加，

導(dǎo)水裂隙帶高度為14.6 m；與此同時(shí)，5#煤層頂板也開始出現(xiàn)離層裂隙，進(jìn)而誘導(dǎo)其在推進(jìn)過程中發(fā)生垮落，處于頂板垮落階段。當(dāng)工作面推進(jìn)至150 m時(shí)，頂板懸露面積不斷增加，當(dāng)基本頂再次達(dá)到極限破壞強(qiáng)度時(shí)，基本頂發(fā)生周期來壓，此時(shí)基本頂上覆巖層橫向裂隙進(jìn)一步發(fā)育；由于兩側(cè)所受的支承壓力不一致，采空區(qū)右側(cè)裂隙明顯發(fā)育，但由于此時(shí)巖體垮落現(xiàn)象明顯，使得采空區(qū)大量裂隙發(fā)生壓實(shí)，導(dǎo)致裂隙總量大幅下降，此時(shí)導(dǎo)水裂隙帶高度達(dá)到26.7 m，處于裂隙壓實(shí)階段。當(dāng)工作面推進(jìn)至180 m時(shí)，采空區(qū)兩側(cè)裂隙縱向發(fā)育，與之前

所形成的縱向裂隙貫通，使得采空區(qū)上部裂隙帶呈梯形狀動(dòng)態(tài)彎曲下沉，此時(shí)，由于新生裂隙生成的速率和閉合裂隙閉合速率達(dá)到動(dòng)態(tài)平衡，處于整體穩(wěn)定階段，確定導(dǎo)水裂隙帶發(fā)育高度為48.8 m；根據(jù)地質(zhì)資料，9#煤層頂板距長興組含水層底界30～52 m，平均39 m，工作面回采過程中存在突水危險(xiǎn)性，威脅工作面安全開采。

根據(jù)試驗(yàn)結(jié)果可得：隨著工作面不斷推進(jìn)，裂隙主要分布在上覆巖層和工作面兩側(cè)，由于在開采擾動(dòng)下，采場巖體不斷發(fā)生采動(dòng)損傷，采動(dòng)覆巖裂隙逐步向前向上發(fā)育，覆巖裂隙發(fā)育總體表現(xiàn)為開切眼和停采線兩側(cè)高度發(fā)育，中部壓實(shí)的“馬鞍形”分布規(guī)律；采空區(qū)兩側(cè)的豎向破斷裂隙區(qū)明顯大于其他區(qū)域，是采動(dòng)裂隙充分發(fā)育區(qū)和含水層卸壓優(yōu)勢通道。

3緩傾斜淺埋煤層覆巖裂隙演化數(shù)值模擬

3.1模型建立

根據(jù)龍鳳煤礦120910工作面地質(zhì)背景，建立UDEC數(shù)值計(jì)算模型，塊體本構(gòu)模型選用Mohr-Coulomb模型，節(jié)理本構(gòu)模型采用Joint contact-Coulomb slip criterion模型。設(shè)計(jì)模型尺寸300 m×160 m，如圖4所示。

為消除數(shù)值模型的邊界效應(yīng)，在模型左、右兩側(cè)分別設(shè)立60 m的邊界煤柱。在數(shù)值模型中，模型內(nèi)各個(gè)單元均受到重力影響，模型上部為自由界面，模型左右兩側(cè)僅限制水平位移，底部為固定邊界，各煤巖層物理力學(xué)參數(shù)見表2。

如上所述，按照黔北煤田龍鳳煤礦綜采工作面9#煤層實(shí)際情況賦予物理力學(xué)參數(shù)、初始條件和邊界條件后，計(jì)算至初始地應(yīng)力平衡狀態(tài)；對9#煤層實(shí)行自左向右開采，從開切眼的位置，每次推進(jìn)10 m，直至工作面回采結(jié)束，待采空區(qū)自然垮落至平衡狀態(tài)。

3.2試驗(yàn)結(jié)果及分析

通過UDEC導(dǎo)出工作面不同推進(jìn)距離下覆巖垮落各階段裂隙孕育特征示意圖，如圖5所示。

通過圖5可得：當(dāng)工作面推進(jìn)至30 m時(shí)，采場新生裂隙開始發(fā)育，并向工作面推進(jìn)方向和覆巖上方拓展，處于裂隙發(fā)育階段；當(dāng)工作面推進(jìn)到90 m時(shí)，基本頂發(fā)生冒落，導(dǎo)致兩側(cè)橫向破斷裂隙和縱向裂隙大量增加，裂隙網(wǎng)絡(luò)發(fā)育復(fù)雜，處于頂板垮落階段；當(dāng)工作面推進(jìn)到150 m時(shí)，覆巖發(fā)生垮落，采空區(qū)兩側(cè)的裂隙始終保持較發(fā)育，而中部的大量裂隙逐漸趨于壓實(shí)，導(dǎo)致裂隙總量大幅下降，處于裂隙壓實(shí)階段；當(dāng)工作面推進(jìn)到180 m時(shí)，采空區(qū)裂隙逐漸壓實(shí)，出現(xiàn)大量裂隙閉合現(xiàn)象，而采空區(qū)兩側(cè)會(huì)出現(xiàn)部分新生裂隙，直至閉合裂隙和新生裂隙處于動(dòng)態(tài)平衡階段，此時(shí)整體處于穩(wěn)定階段。

通過物理試驗(yàn)和數(shù)值模擬結(jié)果進(jìn)行對比分析，發(fā)現(xiàn)淺埋煤層覆巖裂隙演化數(shù)值模擬結(jié)果與相似模擬試驗(yàn)結(jié)果基本吻合，故可用數(shù)值模擬試驗(yàn)裂隙圖來對覆巖裂隙演化規(guī)律進(jìn)行研究。

4覆巖裂隙演化過程分形表征

采動(dòng)裂隙分形維數(shù)是對覆巖裂隙發(fā)育程度的數(shù)字化表征。裂隙發(fā)育程度與分形維數(shù)呈正相關(guān)，分形維數(shù)值越大，其表征的裂隙網(wǎng)絡(luò)結(jié)構(gòu)越復(fù)雜。因此通過數(shù)值模擬試驗(yàn)裂隙圖，結(jié)合分形理論對淺埋煤層覆巖裂隙演化過程進(jìn)行定量化表征。

4.1覆巖裂隙分形維數(shù)計(jì)算

PCAS 是用于孔隙（顆粒）及裂隙圖像識別與分析的專業(yè)軟件。為進(jìn)一步表明淺埋煤層覆巖裂隙分布特征，運(yùn)用PCAS圖像處理技術(shù)將圖5轉(zhuǎn)化為二值圖（圖6），然后導(dǎo)入Fraclab求解分形維數(shù)D，計(jì)算式如下［19］：

D=－limr→0lg N（r）lg r （1）

式中：D為分形維數(shù)值；r為正方形方格大?。籒（r）為含裂隙格子數(shù)。

4.2覆巖裂隙演化過程分形維數(shù)變化規(guī)律

由圖6可見：在采動(dòng)影響下覆巖裂隙網(wǎng)絡(luò)的分形相關(guān)系數(shù)均在98%以上，具有顯著的分形特征。隨著工作面的推進(jìn)，采動(dòng)覆巖裂隙分形維數(shù)呈現(xiàn)動(dòng)態(tài)變化，分形維數(shù)將覆巖裂隙發(fā)育4個(gè)階段分為上升區(qū)，突降區(qū)和穩(wěn)定區(qū)3個(gè)區(qū)（圖7）。

第一階段：Ⅰ—Ⅱ階段（上升區(qū)），在此階段裂隙分形值不斷增加，分維值從1.351上升到1404，工作面推進(jìn)至30～90 m過程中，煤層頂板破斷后不斷垮落，在90 m處分形維數(shù)達(dá)到峰值；第二階段：Ⅱ—Ⅲ階段（突降區(qū)），裂隙分維值從1404突降至1.276，工作面推進(jìn)至90～150 m過程中，關(guān)鍵層破斷發(fā)生大面積垮落，9#煤層上覆煤巖垮落壓實(shí)后，壓實(shí)裂隙發(fā)育，且原有裂隙網(wǎng)絡(luò)的壓實(shí)速率遠(yuǎn)大于壓實(shí)裂隙再發(fā)育的速率；第三階段：Ⅲ—Ⅳ階段（穩(wěn)定區(qū)），裂隙分維值在1.276～1291區(qū)間穩(wěn)定波動(dòng)，工作面推進(jìn)至150～180 m過程中，上覆煤巖處于整體垮落壓實(shí)階段，其中新裂隙生成速率與舊裂隙閉合速率趨于動(dòng)態(tài)平衡。

4.3采動(dòng)損傷分形評價(jià)

近年來，部分學(xué)者將裂縫、孔洞和節(jié)理等視為“損傷場”的連續(xù)介質(zhì)，研究其在巖體中的分布狀態(tài)，并采用分形維數(shù)作為特征參數(shù)表征巖體損傷程度，構(gòu)建采動(dòng)巖體破壞程度分維數(shù)值Dt與試件內(nèi)部“裂隙場”ω之間的關(guān)系［20］，即：

ω=Dt－D0Dt，max－D0（2）

式中：ω為分形損傷度；D0為采動(dòng)前巖體內(nèi)部初始損傷面積的分形維數(shù)；Dt為采動(dòng)后巖體內(nèi)部損傷面積的分形維數(shù)；Dt，max為巖體達(dá)到最大損傷面積時(shí)的分形維數(shù)，對于二維平面問題，Dt，max=2，對于三維立體問題，Dt，max=3。

通過分形損傷度評價(jià)開挖擾動(dòng)下巖體的破壞程度，損傷度越大，巖體破壞也越嚴(yán)重。不同推進(jìn)距離下的分形損傷度如圖8所示。從圖8可得：當(dāng)分形損傷度ωlt;0.60時(shí)，巖體中沒有產(chǎn)生宏觀貫通裂隙，巖體較完整，有宏觀裂隙產(chǎn)生但巖體未產(chǎn)生宏觀損傷，即I區(qū)；當(dāng)0.60≤ωlt;0.64時(shí)，巖體中裂隙發(fā)育較少，為裂隙稀疏區(qū)，即II區(qū)；當(dāng)0.64≤ωlt;0.68時(shí)，巖體被裂隙分成若干碎塊，巖體破碎程度較高，為裂隙過渡區(qū)，即III區(qū)；當(dāng)ω≥0.68時(shí)，巖體幾乎完全破壞，形成致密的宏觀裂隙，為裂隙密集區(qū)，即IV區(qū)。

當(dāng)煤層頂板受到開采擾動(dòng)后，巖體的分形損傷度急劇上升，隨著工作面推進(jìn)至30 m時(shí)，ω達(dá)到0.676，當(dāng)分形損傷度首次處于III區(qū)時(shí)，巖體開始產(chǎn)生大量離層裂隙，誘發(fā)巖體發(fā)生垮落，此時(shí)處于裂隙發(fā)育階段。隨著采動(dòng)不斷進(jìn)行，當(dāng)工作面推進(jìn)至90 m時(shí)，巖體的分形損傷度近似呈線性增加，ω達(dá)到最大值0.702，增加了3.8%，此時(shí)損傷度處于IV區(qū)，巖體表面宏觀裂隙開始大面積擴(kuò)散發(fā)育，不斷向上向前拓展，并且發(fā)生大面積垮塌現(xiàn)象，此時(shí)處于頂板垮落階段。當(dāng)工作面推進(jìn)至150 m時(shí)，ω達(dá)到0.638，分形損傷度也從IV區(qū)—III區(qū)—II區(qū)逐漸遞減，在推進(jìn)過程中，數(shù)值模擬產(chǎn)生的覆巖裂隙會(huì)被大量壓實(shí)，巖體裂隙壓實(shí)的速率遠(yuǎn)大于巖體裂隙產(chǎn)生的速率，故而導(dǎo)致分形損傷度降低，表明巖體破壞的速度逐漸減小，此時(shí)處于裂隙壓實(shí)階段。當(dāng)工作面推進(jìn)至180 m時(shí)，分形損傷度處于II區(qū)，ω達(dá)到0.639，此時(shí)新生裂隙發(fā)育速率與裂隙壓實(shí)速率處于動(dòng)態(tài)平衡，表明該工作面巖層結(jié)構(gòu)達(dá)到新的穩(wěn)定狀態(tài)，此時(shí)處于整體穩(wěn)定階段。

當(dāng)處于分形損傷IV區(qū)時(shí)，為防止在采動(dòng)階段造成頂板大變形或因水分滲透導(dǎo)致強(qiáng)度降低等不利情況，從而降低巷道圍巖的整體承載性能，應(yīng)當(dāng)采取及時(shí)有效的支護(hù)措施，具體包括安裝錨桿錨索進(jìn)行加固，并輔以工字鋼進(jìn)行強(qiáng)化支護(hù)。

5結(jié)論

本文以貴州黔北煤田緩傾斜淺埋煤層為研究對象，綜合運(yùn)用相似模擬試驗(yàn)、數(shù)值模擬和分形理論等研究方法，對采動(dòng)影響下緩傾斜淺埋煤層采動(dòng)覆巖運(yùn)移及裂隙演化規(guī)律進(jìn)行研究，在定性分析的基礎(chǔ)上，定量化表征覆巖裂隙動(dòng)態(tài)演化過程，研究結(jié)果可為類似地質(zhì)條件下頂板巖層控制及水害防治提供參考。研究得到的主要結(jié)論如下：

1）緩傾斜淺埋煤層覆巖裂隙演化過程可分為4個(gè)階段：裂隙發(fā)育階段、頂板垮落階段、裂隙壓實(shí)階段和整體穩(wěn)定階段；覆巖裂隙發(fā)育總體表現(xiàn)為采空區(qū)兩側(cè)高度發(fā)育，中部壓實(shí)的“馬鞍形”分布規(guī)律，最終確定導(dǎo)水裂隙帶發(fā)育高度為48.8 m，回采過程中將貫穿長興組含水層底界，存在頂板突水危險(xiǎn)性。

2）覆巖裂隙演化過程與裂隙網(wǎng)絡(luò)分維值變化整體相似，裂隙網(wǎng)絡(luò)分維值呈上升—突降—穩(wěn)定（1.351～1.404～1.276），其中：新生裂隙發(fā)育到頂板垮落壓實(shí)階段對應(yīng)上升區(qū)，頂板垮落壓實(shí)到壓實(shí)裂隙發(fā)育階段對應(yīng)突降區(qū)，壓實(shí)裂隙發(fā)育到整體穩(wěn)定階段對應(yīng)穩(wěn)定區(qū)。

3）隨著開采不斷進(jìn)行，分形維數(shù)不斷增大，表明裂隙網(wǎng)絡(luò)發(fā)育越復(fù)雜，同時(shí)分形損傷度ω也增加，巖體宏觀破壞程度越高。裂隙演化4個(gè)階段的分形損傷度ω分別為0.676、0.702、0.638、0.639，當(dāng)ω處于裂隙密集區(qū)（IV區(qū)）時(shí)，巖體破裂現(xiàn)象越明顯。

參考文獻(xiàn)：

[1]陳洪松， 岳躍民， 王克林. 西南喀斯特地區(qū)石漠化綜合治理： 成效、問題與對策［J］. 中國巖溶， 2018， 37（1）： 37-42.

［2］陳相標(biāo)， 丁文榮. 石林喀斯特巖溶區(qū)生態(tài)系統(tǒng)服務(wù)時(shí)空演變及權(quán)衡協(xié)同關(guān)系［J］. 水土保持研究， 2023， 30（2）： 285-293.

［3］唐雨生， 蘇培東， 喻洪平， 等.中國巖溶區(qū)鐵路工程主要地質(zhì)問題研究［J］. 地下空間與工程學(xué)報(bào)， 2022， 18（增刊1）： 296-304.

［4］王創(chuàng)業(yè)， 王敬， 劉赟. 淺埋煤層覆巖裂隙分形特征實(shí)驗(yàn)研究［J］. 中國礦業(yè)， 2020， 29（增刊1）： 486-489.

［5］ZHANG J， HE Y F， YANG T， et al. Coevolution mechanism and branch of pillar-overburden fissures in shallow coal seam mining［J］. Energy Science amp; Engineering， 2023， 11（5）： 1630-1642.

［6］黃慶享， 韓金博. 淺埋近距離煤層開采裂隙演化機(jī)理研究［J］. 采礦與安全工程學(xué)報(bào)， 2019， 36（4）： 706-711.

［7］劉國建， 楊富強(qiáng)， 高學(xué)豐， 等. 淺埋薄基巖采場覆巖破斷特征及其導(dǎo)水裂隙帶發(fā)育規(guī)律研究［J］. 煤炭工程， 2023， 55（8）： 108-113.

［8］趙毅鑫， 令春偉， 劉斌， 等. 淺埋超大采高工作面覆巖裂隙演化及能量耗散規(guī)律研究［J］. 采礦與安全工程學(xué)報(bào)， 2021， 38（1）： 9-18， 30.

［9］梁勇， 薛彥平， 胡勝勇. 近距離煤層群重復(fù)采動(dòng)覆巖破壞規(guī)律分析研究［J］. 礦業(yè)研究與開發(fā)， 2021， 41（4）： 72-75.

［10］李友偉， 張玉軍， 肖杰. 多煤層重復(fù)采動(dòng)覆巖破壞高度發(fā)育規(guī)律研究［J］. 煤炭工程， 2022， 54（7）： 97-103.

［11］LIU H， DENG K Z， ZHU X J， et al. Effects of mining speed on the developmental features of mining-induced ground fissures［J］. Bulletin of Engineering Geology and the Environment， 2019， 78（8）： 6297-6309.

［12］李建華， 王帥， 賀艷軍， 等. 煤層群疊加開采采空區(qū)覆巖裂隙演化規(guī)律研究［J］. 煤炭工程， 2021， 53（12）： 92-96.

［13］郭瑞瑞， 鄭凱歌， 席杰， 等. 煤層群開采覆巖運(yùn)移規(guī)律及“三帶”高度確定［J］. 煤礦安全， 2023， 54（8）： 106-117.

［14］張俊， 盛建龍， 黃濤， 等. 巖體裂隙網(wǎng)絡(luò)滲透特性與分形特征研究［J］. 化工礦物與加工， 2020， 49（11）： 29-34.

［15］梁濤， 劉曉麗， 王思敬. 采動(dòng)裂隙拓展規(guī)律及滲透特性分形研究［J］. 煤炭學(xué)報(bào)， 2019， 44（12）： 3729-3739.

［16］鄒恩隆， 李志華， 馬軍， 等. 采動(dòng)裂隙動(dòng)態(tài)演化定量分析及其應(yīng)用研究［J］. 礦業(yè)研究與開發(fā)， 2023， 43（7）： 101-107.

［17］郝鑫宇. 覆巖采動(dòng)裂隙發(fā)育規(guī)律的電法無損探測實(shí)驗(yàn)研究［D］. 徐州：中國礦業(yè)大學(xué)， 2019.

［18］GAO M Z， ZHANG R， XIE J， et al. Field experiments on fracture evolution and correlations between connectivity and abutment pressure under top coal caving conditions［J］. International Journal of Rock Mechanics and Mining Sciences， 2018， 111： 84-93.

［19］LI J， DU Q， SUN C X. An improved box-counting method for image fractal dimension estimation［J］. Pattern Recognition， 2009， 42（11）： 2460-2469.

［20］楊仁樹， 許鵬. 爆炸作用下介質(zhì)損傷破壞的分形研究［J］. 煤炭學(xué)報(bào)， 2017， 42（12）： 3065-3071.

（責(zé)任編輯：周曉南）

Abstract：

It is of great significance to clarify the evolution law of water-conducting fissure in mining overlying rock for the study of water control in roof. Based on 120910 fully mechanized mining face of No.9 coal seam in Longfeng Coal Mine as engineering background， the evolution law of overlying rock fractures in gently inclined shallow coal seam is studied by means of similarity simulation and numerical simulation. The results show that the development height of water-conducting fracture zone in coal seam mining roof will run through the bottom boundary of the Changxing Formation aquifer， and there is a risk of water inrush during the advancing process of the working face. The evolution process of overburden rock fracture can be divided into four stages， namely， fracture development stage， roof caving stage， fracture compaction stage and overall stability stage. At the same time， fractal theory is introduced to quantitatively describe the fracture evolution process of mining overburden rock. The fractal damage degree ω in four stages of fracture evolution is 0.676， 0.702， 0.638 and 0.639， respectively. When ω≥0.68， the rock mass fracture phenomenon is more obvious. The results can provide reference for the study of water control in overlying strata roof under similar geological conditions.

Key words：

gently inclined coal seam; crack evolution; fractal damage degree; analogy simulation; numerical simulation

收稿日期：2024-05-16

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目（51964006，52164007）；貴州省科技支撐計(jì)劃資助項(xiàng)目（黔科合支撐[2022]一般247，黔科合支撐[2022]一般248）；貴州省高層次創(chuàng)新型人才培養(yǎng)資助項(xiàng)目（黔科合人才[2016]4011號）；貴州省礦山動(dòng)力災(zāi)害預(yù)警與控制技術(shù)科技創(chuàng)新人才團(tuán)隊(duì)資助項(xiàng)目（黔科合平臺(tái)人才[2019]5619）

作者簡介：蒯?。?997—），男，在讀碩士，研究方向：礦山壓力與巖層控制，E-mail：kuaichen97@163.com.

*通訊作者：左宇軍，E-mail：zuo_yujun@163.com.