李允生, 徐德生, 馬志鋒, 周海軍, 郭文豪

(1.山東東山古城煤礦有限公司， 山東 濟(jì)寧 272000；2.徐州弘毅科技發(fā)展有限公司， 江蘇 徐州 221018；3.中國(guó)礦業(yè)大學(xué) 礦業(yè)工程學(xué)院， 江蘇 徐州 221116)

0 引言

近幾年，隨著煤礦采深不斷加大，沖擊地壓顯現(xiàn)越來(lái)越突出[1]。沖擊地壓的發(fā)生與地質(zhì)構(gòu)造、煤柱留設(shè)等密切相關(guān)，尤其在復(fù)雜地質(zhì)構(gòu)造區(qū)域，煤柱受多個(gè)構(gòu)造、支承壓力疊加影響，應(yīng)力集中度高，復(fù)雜構(gòu)造區(qū)域的覆巖破斷或垮落可能誘發(fā)強(qiáng)礦震、沖擊地壓等現(xiàn)象[2-4]。針對(duì)復(fù)雜地質(zhì)構(gòu)造區(qū)域的沖擊地壓危險(xiǎn)監(jiān)測(cè),研究者提出了眾多方法，如微震法[5]、鉆屑法[6]、地音法[7]和電磁輻射法[8]等。但每種方法僅反映了煤巖某一方面的特征，且大多方法僅能監(jiān)測(cè)而不能預(yù)警。如微震法只能對(duì)區(qū)域和局部微震進(jìn)行監(jiān)測(cè)而不能預(yù)警，微震法中的微震監(jiān)測(cè)也多為事后記錄，無(wú)法進(jìn)行預(yù)測(cè)；鉆屑法為點(diǎn)監(jiān)測(cè)，可用于檢驗(yàn)和預(yù)警；地音法和電磁輻射法也為局部監(jiān)測(cè)，不能預(yù)警。上述方法的監(jiān)測(cè)指標(biāo)單一，針對(duì)性差，不能在復(fù)雜地質(zhì)構(gòu)造區(qū)起到良好的監(jiān)測(cè)預(yù)警效果。

礦震震動(dòng)波波速層析成像(Computerized Tomography，CT)技術(shù)作為復(fù)雜地質(zhì)條件下新興的沖擊危險(xiǎn)監(jiān)測(cè)預(yù)警技術(shù)，部分學(xué)者對(duì)此展開(kāi)了相關(guān)研究。鞏思園[9]圍繞實(shí)驗(yàn)室研究得到的震動(dòng)波波速與應(yīng)力間的試驗(yàn)關(guān)系模型，利用震動(dòng)波信號(hào)建立了礦震震動(dòng)波波速層析成像模型，建立了波速異常、波速梯度變化等預(yù)警指標(biāo)，用于預(yù)測(cè)預(yù)報(bào)沖擊危險(xiǎn)。竇林名等[10-11]基于縱波波速開(kāi)采區(qū)域沖擊危險(xiǎn)的理論基礎(chǔ)，建立了采掘區(qū)域沖擊危險(xiǎn)性動(dòng)態(tài)震動(dòng)波CT 探測(cè)技術(shù)，并給出了沖擊危險(xiǎn)判別準(zhǔn)則。程剛等[12]以某礦3下113回采工作面為工程背景，研究了震動(dòng)波CT反演技術(shù)在工作面受采空、斷層、煤柱與采動(dòng)影響的多重作用下的應(yīng)用效果。曹安業(yè)等[13]基于某礦7192回采工作面的覆巖空間結(jié)構(gòu)及其運(yùn)移規(guī)律研究，揭示了由斷層錯(cuò)動(dòng)導(dǎo)致的工作面動(dòng)力顯現(xiàn)及大能量礦震的頻發(fā)機(jī)理，同時(shí)利用震動(dòng)波CT技術(shù)對(duì)工作面沖擊危險(xiǎn)性進(jìn)行動(dòng)態(tài)預(yù)警。以上研究表明，CT反演能較好反映現(xiàn)場(chǎng)工作面的危險(xiǎn)性，但是在實(shí)際操作過(guò)程中，CT反演由于其反演周期長(zhǎng)，反演結(jié)果往往不能及時(shí)反映工作面的當(dāng)前危險(xiǎn)程度。

綜上所述，單獨(dú)使用微震監(jiān)測(cè)技術(shù)，其微震監(jiān)測(cè)范圍僅能反映微震所在區(qū)域的危險(xiǎn)性，給現(xiàn)場(chǎng)監(jiān)測(cè)帶來(lái)了困難；單獨(dú)使用CT反演監(jiān)測(cè)在煤礦現(xiàn)場(chǎng)可操作性有限，需要復(fù)雜的計(jì)算和判斷，對(duì)數(shù)據(jù)分析人員的能力要求和專(zhuān)業(yè)要求很高，普適性較低。因此需要找到一種煤礦分析人員普遍好理解且現(xiàn)場(chǎng)好操作、執(zhí)行明確的監(jiān)測(cè)預(yù)警方法。鑒此，本文以山東濟(jì)寧東山古城煤礦3105工作面為工程背景，提出了一種微震-應(yīng)力多維信息監(jiān)測(cè)預(yù)警方法，將短期的微震監(jiān)測(cè)和中長(zhǎng)期震動(dòng)波CT反演監(jiān)測(cè)技術(shù)相結(jié)合對(duì)工作面沖擊危險(xiǎn)區(qū)域進(jìn)行判識(shí)，并結(jié)合基于中長(zhǎng)期CT反演的大范圍監(jiān)測(cè)結(jié)果，利用工作面周邊的微震數(shù)據(jù)對(duì)其修正，實(shí)現(xiàn)對(duì)工作面中長(zhǎng)期應(yīng)力演變及短期微震集聚的實(shí)時(shí)、動(dòng)態(tài)監(jiān)測(cè)預(yù)警。

1 工程背景

古城煤礦位于山東省濟(jì)寧市兗州區(qū)東郊，本區(qū)地勢(shì)平坦，為沖積平原，地勢(shì)東高西低，一般高程為+51～+56 m。其中3105工作面位于31采區(qū)的中部，工作面下至31采區(qū)集中軌道，上部隔離煤柱緊鄰2107采空區(qū)，工作面東部為3106工作面，西南為31采區(qū)未開(kāi)拓區(qū)域。地面標(biāo)高為+52.7～+53.6 m，工作面標(biāo)高為-1 060～-1 140 m。

在回采期間，工作面前方進(jìn)入到臨3106采空區(qū)不規(guī)則煤柱及斷層影響區(qū)域，礦震頻發(fā)，且礦震頻次及能量呈升高趨勢(shì)。3105工作面局部復(fù)雜地質(zhì)及采掘區(qū)域平面圖如圖1所示。

圖1 3105工作面局部復(fù)雜地質(zhì)及采掘區(qū)域平面圖Fig.1 Plan of local complex geology and mining area of 3105 working face

工作面所采煤層為3號(hào)煤層，平均厚度為8.6 m，煤層傾角為0～15°，堅(jiān)固性系數(shù)f為2～3，煤層結(jié)構(gòu)簡(jiǎn)單。工作面回采范圍內(nèi)兩巷道及切眼在掘進(jìn)過(guò)程中共揭露落差H為0～8 m的正斷層12條。3105工作面綜合柱狀圖如圖2所示。

圖2 3105工作面綜合柱狀圖Fig.2 Comprehensive histogram of 3105 working face

2 3105工作面微震時(shí)空演化分析

古城煤礦安裝有SOS微震監(jiān)測(cè)系統(tǒng)，通過(guò)3105工作面布置的6個(gè)微震探頭對(duì)工作面采掘過(guò)程中的煤巖體震動(dòng)信號(hào)進(jìn)行采集、記錄和分析，通過(guò)對(duì)多組波形處理，確定微震事件三維定位和能量，微震事件定位原理如圖3所示，圖中x0，y0，z0和x，y，z分別為震源和臺(tái)站的空間坐標(biāo)，t0和ti(i=1,2,…,6,為微震探頭數(shù)量)分別為震源的發(fā)生時(shí)間和臺(tái)站的接收時(shí)間。

圖3 微震事件定位原理Fig.3 Microseismic positioning principle

從2019年8月開(kāi)始，3105工作面微震頻次能量逐漸升高，如圖4所示。7月底之前微震能量處于低水平狀態(tài)，且時(shí)序發(fā)展平穩(wěn)，但從07-28開(kāi)始，連續(xù)5 d微震頻次能量上升，進(jìn)入8月，微震能量及頻次處于高水平狀態(tài)，時(shí)序發(fā)展極其不穩(wěn)定，微震平面分布如圖5所示。由圖5可知，工作面前方煤柱區(qū)內(nèi)微震能量、頻次明顯增加，表明工作面超前應(yīng)力范圍已開(kāi)始進(jìn)入到臨3106采空區(qū)不規(guī)則工作面區(qū)域，且工作面前方開(kāi)采揭露斷層，進(jìn)一步造成微震事件積聚。在2019-08-01—15，微震大能量事件(能量大于104J)分布明顯比7月多，微震分布范圍相對(duì)7月較廣，且工作面前方逐漸變窄的區(qū)段煤柱區(qū)域出現(xiàn)了較多微震事件。

圖4 微震能量及頻次時(shí)序活動(dòng)曲線(xiàn)Fig.4 Time series activity curves of microseismic energy and frequency

(a) 2019-07-15—31的微震平面分布

(b) 2019-08-01—15的微震平面分布圖5 微震平面分布Fig.5 Plane distribution of microseisms

通過(guò)近1個(gè)月微震事件的時(shí)空分布異常情況判斷，可以斷定工作面前方變窄的區(qū)段煤柱區(qū)域?yàn)閼?yīng)力異常區(qū)域，但僅通過(guò)微震事件的分布不足以精確判斷高應(yīng)力區(qū)和準(zhǔn)確實(shí)施卸壓措施的范圍，需要通過(guò)更加精準(zhǔn)的CT反演技術(shù)進(jìn)行危險(xiǎn)區(qū)域劃分和預(yù)警，為卸壓措施提供可靠的依據(jù)。

3 復(fù)雜應(yīng)力區(qū)CT反演及其應(yīng)用結(jié)果

3.1 震動(dòng)波CT反演原理

震動(dòng)波CT反演技術(shù)就是利用傳感器測(cè)定的巖體內(nèi)初始平均波速與監(jiān)測(cè)到的微破裂位置為初始數(shù)據(jù)，利用人工標(biāo)定波形的P波到時(shí)，確定初始傳播時(shí)間，通過(guò)反演計(jì)算不斷調(diào)整微破裂位置與速度模型，最終重構(gòu)和反演巖樣內(nèi)部在不同加載階段的速度場(chǎng)分布規(guī)律，探測(cè)原理如圖6所示。研究結(jié)果表明，煤巖體在應(yīng)力作用下會(huì)引起震動(dòng)波波速的改變，震動(dòng)波傳播通過(guò)工作面煤巖體時(shí)，煤巖體上所受的應(yīng)力越高，震動(dòng)傳播的速度就越快[9]。

圖6 震動(dòng)波CT反演技術(shù)探測(cè)原理Fig.6 Exploration principle of seismic wave CT inversion technology

3.2 工作面復(fù)雜應(yīng)力區(qū)CT反演結(jié)果分析

基于微震異常情況，為更準(zhǔn)確地預(yù)測(cè)沖擊危險(xiǎn)區(qū)域，使卸壓方案有的放矢，古城煤礦從2019年8月開(kāi)始，以15 d的微震數(shù)據(jù)為基礎(chǔ)，進(jìn)行被動(dòng)CT反演，選擇測(cè)站為6、12、13、15、18、19(圖1)。選取2019-08-01—15、08-16—31及09-01—15三個(gè)階段進(jìn)行對(duì)比分析，分別選取212、192、162 個(gè)微震有效波形，選取標(biāo)高為-1 160 m的波速反演結(jié)果切片。3個(gè)階段的反演結(jié)果如圖7所示。

(a) 08-01—15

(b) 08-16—31

(c) 09-01—15圖7 3105工作面復(fù)雜應(yīng)力區(qū)CT反演結(jié)果Fig.7 CT inversion results of complex stress area in 3105 working face

沖擊礦壓事故多發(fā)生于煤層巷道內(nèi)[14-15]。通過(guò)圖1可以發(fā)現(xiàn)，3105工作面材料巷和運(yùn)輸巷分別受到斷層和煤柱的影響，沖擊風(fēng)險(xiǎn)較高。因此結(jié)合古城煤礦歷史CT反演數(shù)據(jù)，以?xún)上锔卟ㄋ?≥4.2 km/s)反演區(qū)域作為研究對(duì)象對(duì)兩巷微震和沖擊危險(xiǎn)演化進(jìn)行分析：第1階段(08-01—15)的反演結(jié)果中高波速區(qū)域?yàn)椴牧舷?50～880 m(230 m)及運(yùn)輸巷590～780 m(190 m)、930～960 m(30 m)區(qū)域；第2階段(08-16—31)的反演結(jié)果中高波速區(qū)域?yàn)椴牧舷?00～690 m(90 m)及運(yùn)輸巷690～720 m(30 m)區(qū)域；第3個(gè)階段(09-01—15)的反演結(jié)果的高波速區(qū)域?yàn)椴牧舷?50～910 m(60 m)及運(yùn)輸巷950～1 030 m(80 m)區(qū)域內(nèi)。對(duì)于3105工作面，其波速及應(yīng)力異常區(qū)域主要集中于材料巷和運(yùn)輸巷超前100～200 m范圍，分布范圍最大可達(dá)230 m。因此，沖擊危險(xiǎn)卸壓區(qū)域應(yīng)位于工作面超前兩巷位置。

根據(jù)各階段的CT反演結(jié)果，針對(duì)性地對(duì)高波速區(qū)域加強(qiáng)監(jiān)測(cè)，并及時(shí)采取有效的卸壓措施。實(shí)施卸壓措施前后材料巷和運(yùn)輸巷3個(gè)階段的高波速區(qū)域最大范圍分別為230,90,60 m和190,120,30 m，即高波速區(qū)域范圍呈現(xiàn)變小的趨勢(shì)，并且由于斷層X(jué)F17影響,材料巷相對(duì)比運(yùn)輸巷沖擊危險(xiǎn)程度高。綜合反演結(jié)果可得出，08-16—31的應(yīng)力顯現(xiàn)低于08-01—15的應(yīng)力顯現(xiàn)，09-01—15的應(yīng)力顯現(xiàn)低于08-16—31的應(yīng)力顯現(xiàn)。

4 解危措施及效果檢驗(yàn)

通過(guò)每15 d 的工作面前方CT反演結(jié)果，給出確定的強(qiáng)沖擊和中等沖擊危險(xiǎn)區(qū)域范圍，依據(jù)分區(qū)管理的原則，對(duì)高波速?zèng)_擊危險(xiǎn)區(qū)域采取了及時(shí)的煤層大直徑鉆孔[16]補(bǔ)強(qiáng)卸壓和煤體爆破[17]解危措施，同時(shí)對(duì)該區(qū)域頂板進(jìn)行深孔爆破卸壓[18]，爆破層位為煤層上方厚度為21 m的中砂巖層。頂板爆破原理如圖8所示，通過(guò)采空區(qū)頂板爆破工程，主要實(shí)現(xiàn)2個(gè)卸壓目的：① 通過(guò)對(duì)采空區(qū)懸頂爆破，促進(jìn)采空區(qū)頂板的彎曲下沉乃至垮斷，縮短采空區(qū)頂板懸頂長(zhǎng)度，實(shí)現(xiàn)對(duì)煤柱靜載的卸壓及頂板周期垮斷時(shí)動(dòng)載強(qiáng)度的降低。② 通過(guò)爆破，降低頂板巖體的完整程度，促進(jìn)其裂隙發(fā)展乃至破壞，降低巖體的應(yīng)力集中程度，同時(shí)能夠提高其抵抗動(dòng)載沖擊的能力。

圖8 頂板爆破原理Fig.8 Principle of roof blasting

卸壓具體參數(shù)如下。

大直徑鉆孔卸壓參數(shù)：孔徑為110～150 mm，孔深為20 m，間距為1 m。

煤層爆破卸壓參數(shù)：孔深為15 m，間距為4 m，孔徑為110 mm，裝藥量為11 kg，裝藥長(zhǎng)度為4 m，封孔長(zhǎng)度為4 m，起爆間隔時(shí)間為30 min，采用煤礦二級(jí)乳化炸藥，一次起爆一孔。

切頂爆破卸壓參數(shù)：致裂層為煤層上覆厚度為21 m的中砂巖層，孔深為23 m，孔徑為94 mm，間距為6～10 m，裝藥量為30 kg，裝藥長(zhǎng)度為 10 m，封孔長(zhǎng)度為13 m，起爆間隔時(shí)間為30 min，采用煤礦二級(jí)乳化炸藥，一次起爆一孔。頂板爆破卸壓工程如圖9所示。

圖9 頂板爆破卸壓工程Fig.9 Schematic diagram of pressure relief blasting project of roof

對(duì)圖7(a)CT反演結(jié)果中的工作面前方200 m范圍外危險(xiǎn)區(qū)域?qū)嵤╉敯灞?、煤體爆破卸壓工程，實(shí)施日期為2019-08-16以后。從圖7(b)和圖7(c)CT反演結(jié)果可看出，實(shí)施爆破卸壓后，工作面前方的沖擊危險(xiǎn)性明顯降低，表現(xiàn)為工作面前方高沖擊危險(xiǎn)區(qū)域減小及危險(xiǎn)區(qū)域前移(遠(yuǎn)離工作面)，卸壓效果明顯，高波速區(qū)域明顯遠(yuǎn)離工作面及巷道，且范圍減小。

每15 d的微震統(tǒng)計(jì)數(shù)據(jù)如圖10所示，從圖10可看出，工作面回采至縮窄煤柱區(qū)域，微震頻次和總能量呈逐漸上升趨勢(shì)，微震頻次由420上升至515，上升了22.6%，微震事件總能量由10.3×105J上升至12×105J，上升了16.5%。在采取了頂板爆破卸壓措施后，微震頻次及總能量下降明顯，15 d和30 d后的微震頻次由515降低至387、338，分別下降了24.9%、34.4%；微震能量由12×105J降低至8.98×105、5.98×105J，分別降低了25.2%、50.2%，卸壓效果明顯。

圖10 每15 d微震頻次及能量曲線(xiàn)Fig.10 Microseismic frequency and energy curves every 15 days

5 結(jié)論

(1) 每天微震活動(dòng)頻次、能量升高及異常變化可作為微震活動(dòng)的異常情況，并為實(shí)施CT監(jiān)測(cè)預(yù)警提供依據(jù)。

(2) 利用震動(dòng)波CT反演技術(shù)能夠有效地預(yù)測(cè)沖擊危險(xiǎn)區(qū)域，而且能夠區(qū)分危險(xiǎn)等級(jí)。3105工作面波速及應(yīng)力異常區(qū)域主要分布在工作面超前100～200 m的巷道區(qū)域。震動(dòng)波CT反演技術(shù)有效提高了預(yù)警準(zhǔn)確性，同時(shí)可指導(dǎo)實(shí)施針對(duì)性卸壓措施，降低了卸壓勞動(dòng)強(qiáng)度。

(3) 通過(guò)每15 d一次的CT反演給出的巷道沖擊危險(xiǎn)結(jié)果，分別實(shí)施大直徑鉆孔卸壓和頂板深孔爆破方案。微震監(jiān)測(cè)和CT反演結(jié)果表明，采取爆破卸壓措施后，微震頻次及總能量下降明顯，高波速區(qū)域的微震頻次及能量分別由515、12×105J降低至338、5.98×105J，分別下降了34.4%、50.2%，取得了很好的沖擊防治效果。