竇林名，曹晉榮，曹安業(yè)，柴彥江，白金正，闞吉亮

(1.中國礦業(yè)大學(xué) 深部煤炭資源開采教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221116；2.中國礦業(yè)大學(xué) 礦業(yè)工程學(xué)院，江蘇 徐州 221116)

0 引 言

采礦誘發(fā)地震是指地面或幾百米淺層和上千米深層的礦山開采引起的地震活動(dòng)，簡稱礦震[1]。在煤礦中，礦震與沖擊地壓不能一概而論，礦震不一定會(huì)導(dǎo)致沖擊地壓的發(fā)生[2]。但少數(shù)強(qiáng)礦震發(fā)生后，可能誘發(fā)煤礦井下沖擊地壓[3-5]和煤與瓦斯突出[6-7]等災(zāi)害，有時(shí)甚至導(dǎo)致地面晃動(dòng)、地表塌陷、建筑物損壞等嚴(yán)重后果[8]，在造成人員傷亡和設(shè)備損壞的同時(shí)，容易引發(fā)社會(huì)問題。

目前，針對煤礦礦震類型與特征、發(fā)生機(jī)制、傳播規(guī)律及監(jiān)測預(yù)警技術(shù)，國內(nèi)外學(xué)者開展了諸多研究。朱佩武[9]區(qū)分了礦震、天然地震和爆破的震相特征，并據(jù)此將礦震分為沖擊型與重力型2類。李鐵等[10]對海孜煤礦9個(gè)強(qiáng)礦震進(jìn)行震源機(jī)制分析，結(jié)合區(qū)域地應(yīng)力測量結(jié)果，提出了海孜礦強(qiáng)礦震發(fā)生機(jī)理及防治思路。曹安業(yè)[11]通過矩張量反演方法研究了采動(dòng)煤巖震源破裂機(jī)理，結(jié)果表明該方法可以較好地揭示礦震震源破裂特征。魏東等[12]基于煤礦現(xiàn)場微震監(jiān)測，結(jié)合覆巖空間結(jié)構(gòu)理論，分析得到了相鄰采空區(qū)下工作面回采關(guān)鍵層失穩(wěn)誘發(fā)礦震機(jī)理及防治方案。王樹立等[13]研究了超厚高位紅層砂巖破斷礦震活動(dòng)規(guī)律，強(qiáng)礦震的發(fā)生和紅層破斷密切相關(guān)。陳學(xué)華等[14]研究了地壘構(gòu)造區(qū)域內(nèi)工作面礦震分布特征，認(rèn)為斷層活化導(dǎo)致礦震能量及頻次升高。苗小虎等[15]對一例異常礦震信號進(jìn)行深入挖掘，提出礦震震動(dòng)破壞誘發(fā)高應(yīng)力區(qū)發(fā)生沖擊地壓機(jī)理。曹安業(yè)等[16]分析了礦震震動(dòng)波能量的傳播模式并給出了礦震誘沖的條件。為進(jìn)一步精確描述礦震震動(dòng)波傳播過程，潘一山等[17]基于均勻各向同性介質(zhì)假設(shè)，通過直接求解三維波動(dòng)微分方程得到了描述震動(dòng)波在三維空間中傳播的三維激震模型。朱廣安等[18]通過數(shù)值模擬發(fā)現(xiàn)斷層對礦震應(yīng)力波起著明顯的阻隔和衰減作用。高明仕等[19]通過在不同介質(zhì)中進(jìn)行震動(dòng)波傳播試驗(yàn)，得出巖土介質(zhì)中能量與傳播距離呈乘冪關(guān)系衰減。為了準(zhǔn)確預(yù)測礦震沖擊災(zāi)害，竇林名等[20]結(jié)合微震監(jiān)測和互聯(lián)網(wǎng)等技術(shù)，搭建了礦震沖擊災(zāi)害遠(yuǎn)程在線預(yù)警平臺，并成功應(yīng)用于現(xiàn)場。

已有的礦震分類多是定性描述，缺乏量化指標(biāo)，同時(shí)礦震危險(xiǎn)性的判別值得關(guān)注。煤巖體中原生節(jié)理、采動(dòng)裂縫以及地質(zhì)構(gòu)造等不連續(xù)面對礦震震動(dòng)波傳播影響顯著，理論分析及有限元數(shù)值模擬研究方法存在一定局限性。筆者分析了多個(gè)礦區(qū)不同地質(zhì)和開采條件下礦震震源信息、波形和井下礦壓顯現(xiàn)情況，根據(jù)礦震震動(dòng)特征及危險(xiǎn)程度對其進(jìn)行分類。采用UDEC數(shù)值模擬研究了節(jié)理巖體震動(dòng)波傳播規(guī)律，揭示了不連續(xù)面對震動(dòng)波應(yīng)力傳播的作用機(jī)制。提出了礦震震動(dòng)波能量衰減擬合分析和危險(xiǎn)性礦震判別方法，以某礦為工程背景，分別分析了實(shí)體煤區(qū)域和采空區(qū)內(nèi)礦震震動(dòng)波衰減特性。結(jié)果表明，現(xiàn)場實(shí)測與模擬結(jié)果具有較好的一致性。

1 煤礦礦震類型

礦震是煤礦采掘過程中煤巖體對區(qū)域或局部應(yīng)力調(diào)整的一種響應(yīng)，通常伴隨著能量釋放與震動(dòng)，其本質(zhì)原因是煤巖體中存在高應(yīng)力或高應(yīng)力差。如圖1所示，在煤礦開采中，煤體宏觀破裂及失穩(wěn)、巖層破斷及運(yùn)動(dòng)、斷層活化等均會(huì)引發(fā)礦震。因此，采深、斷層、褶曲、煤柱、堅(jiān)硬頂板以及回采速度等是影響礦震頻次和能量的重要因素。

圖1 礦震成因示意Fig.1 Schematic of causes of coal mine tremor

基于SOS微震監(jiān)測系統(tǒng)及礦井礦震遠(yuǎn)程在線監(jiān)測預(yù)警平臺，捕獲并積累了大量不同地質(zhì)、開采技術(shù)條件下煤礦開采過程中的礦震震動(dòng)信號。在海量數(shù)據(jù)的基礎(chǔ)上，通過對礦震能量、波形特征、震源位置、震動(dòng)時(shí)煤壁震動(dòng)速度以及煤礦井下礦壓顯現(xiàn)特征等的綜合分析，提出將礦震分為采動(dòng)破裂型、巨厚覆巖型和高能震動(dòng)型3種類型，見表1。其中采動(dòng)破裂型礦震是指采掘過程中由煤層及附近頂?shù)装鍘r層破裂產(chǎn)生的礦震，其能量一般小于104J，屬于采掘狀態(tài)下采場周圍煤巖體破裂有序釋放能量的正?，F(xiàn)象；巨厚覆巖型礦震是指距煤層100 m以上、厚度大于100 m、巖石強(qiáng)度相對不大的巨厚巖層在采空區(qū)上方破斷、滑移產(chǎn)生的礦震，能量大于105J，這類礦震大部分能被地震臺網(wǎng)記錄；高能震動(dòng)型礦震是指能量104J以上，且震源位于采掘工作面附近實(shí)體煤及其頂?shù)装鍘r層之中，此類礦震也可能被地震臺網(wǎng)記錄，根據(jù)發(fā)生主體不同，可將高能震動(dòng)型礦震進(jìn)一步細(xì)分為煤體內(nèi)爆型、頂板失穩(wěn)型和斷層活化型。

表1 煤礦礦震分類Table 1 Classification of coal mine tremor

礦震是否對井下產(chǎn)生危害與其能量大小及相對采場的位置有關(guān)，并不是所有的高能量礦震均會(huì)誘發(fā)沖擊地壓等礦井災(zāi)害。采動(dòng)破裂型礦震雖然分布于采場周圍的煤巖體中，但震動(dòng)能量小，屬于安全性礦震；巨厚覆巖型礦震產(chǎn)生的能量大，但距離采掘工作面較遠(yuǎn)且位于采空區(qū)中，震動(dòng)波衰減較快，因此一般情況下對井下工作人員及設(shè)備的威脅不大，但在特殊地質(zhì)和開采因素下，巨厚覆巖型礦震也具有危險(xiǎn)性；高能震動(dòng)型礦震震源處于或鄰近采場空間，同時(shí)釋放能量較大，較強(qiáng)的震動(dòng)波突然作用于采場周圍煤巖體中的應(yīng)力場和滲流場時(shí)，極易誘發(fā)沖擊地壓、煤與瓦斯突出、突水等礦井災(zāi)害。結(jié)合礦震發(fā)生時(shí)煤礦井下礦壓顯現(xiàn)情況，從礦震對井下人員或設(shè)備危害性大小的角度出發(fā)，可以將礦震分為正常礦震和危險(xiǎn)性礦震2類，如圖2所示。

圖2 礦震危險(xiǎn)程度分類Fig.2 Classification of degree of danger of coal mine tremor

2 礦震震動(dòng)波傳播數(shù)值模擬

礦震震動(dòng)波傳播規(guī)律是揭示采動(dòng)動(dòng)載誘發(fā)沖擊地壓機(jī)理的重要基礎(chǔ)[21]。采用UDEC數(shù)值模擬方法研究巖體中震動(dòng)波傳播引起的應(yīng)力響應(yīng)和不連續(xù)面對震動(dòng)波傳播作用機(jī)制。由于考慮了不連續(xù)面的作用，該數(shù)值模擬方法被廣泛應(yīng)用于采礦工程相關(guān)問題研究中[22,23]。CAI等[24]采用UDEC研究了多組平行裂縫對巖體彈性波衰減的影響規(guī)律。DENG等[25]研究了節(jié)理巖體中地下結(jié)構(gòu)在爆炸誘發(fā)應(yīng)力波作用下的損傷機(jī)理。

2.1 數(shù)值模型與方案

在UDEC數(shù)值模型中，煤巖體由一系列塊體和塊體間的節(jié)理構(gòu)成，節(jié)理按照力學(xué)準(zhǔn)則允許發(fā)生剪切滑動(dòng)或張拉位移，節(jié)理失效與否取決于其強(qiáng)度與應(yīng)力狀態(tài)，當(dāng)微觀損傷累積到一定程度后可觀察到煤巖體的宏觀破壞。數(shù)值模型如圖3所示，長60 m、高50 m，節(jié)理間距設(shè)置為10 m，模型頂部的箭頭代表所施加應(yīng)力的方向，為研究震動(dòng)波傳播過程中質(zhì)點(diǎn)應(yīng)力及速度變化規(guī)律，共布設(shè)5個(gè)監(jiān)測點(diǎn)，高度H分別為5、15、25、35、45 m。

圖3 UDEC數(shù)值模型Fig.3 UDEC numerical model

巖層塊體采用各向同性的線彈性模型，節(jié)理本構(gòu)關(guān)系為Mohr-Coulomb滑移模型，模型物理力學(xué)參數(shù)采用試錯(cuò)法進(jìn)行驗(yàn)證，巖石塊體和節(jié)理性質(zhì)見表2，表中E為巖石彈性模量，kn、ks分別指節(jié)理面的法向、切向剛度。將礦震震動(dòng)波簡化為半正弦拉伸應(yīng)力波，幅值為7 MPa、頻率為10 Hz，如圖4所示。巖石的波速一般為4 000 ～7 000 m/s，可知該模型下震動(dòng)波波長為400～700 m，模型最小尺寸為0.5 m，遠(yuǎn)小于波長的1/8～1/10，震動(dòng)波參數(shù)滿足在模型中傳播的要求。

圖4 礦震震動(dòng)波示意Fig.4 Schematic of coal mine tremor waves

表2 震動(dòng)波傳播數(shù)值模擬參數(shù)Table 2 Numerical simulation parameters of shock wave propagation under impact load

為模擬礦震震動(dòng)波在節(jié)理巖體中傳播，將數(shù)值模擬計(jì)算分為靜態(tài)和動(dòng)態(tài)2個(gè)階段。在靜態(tài)計(jì)算模式下，分別向模型垂直、水平方向上施加5.0、7.5 MPa的原始應(yīng)力，固定模型四周邊界，將數(shù)值模型計(jì)算至平衡。待模型靜力平衡后，改變邊界條件，輸入震動(dòng)參數(shù)，打開動(dòng)力分析模式進(jìn)行計(jì)算。

動(dòng)力分析需選用合理的邊界條件和阻尼。為減少震動(dòng)波在模型邊界反射的影響，改變平衡后模型的邊界條件，將兩側(cè)及下部邊界設(shè)置為黏滯邊界，模型頂部變?yōu)閼?yīng)力邊界。阻尼采用Rayleigh阻尼，Rayleigh阻尼可以表示為

C=λM+ηK

(1)

式中：C為黏性阻尼矩陣；M為質(zhì)量矩陣；K為剛度矩陣；λ、η分別為質(zhì)量阻尼比例系數(shù)與剛度阻尼比例系數(shù)。

對于Rayleigh阻尼，需要確定2個(gè)參數(shù)：臨界阻尼比ξmin和模型的固有頻率fmin。其中臨界阻尼比ξmin一般為2%～5%，當(dāng)模型可能發(fā)生大變形或者節(jié)理可能產(chǎn)生大的位移時(shí)，應(yīng)選用較小的臨界阻尼比，固有頻率fmin通常為輸入波的頻率[26]。因此，選用的臨界阻尼比為2%，固有頻率為10 Hz。

為了研究不連續(xù)面對震動(dòng)波傳播的作用機(jī)制，根據(jù)節(jié)理是否發(fā)生破壞共設(shè)置了2組數(shù)值試驗(yàn)，由于應(yīng)力波為拉伸應(yīng)力，模型更易出現(xiàn)拉伸破壞，因此試驗(yàn)中僅改變節(jié)理抗拉強(qiáng)度，其余參數(shù)相同，具體如下：①模型中設(shè)置節(jié)理面，節(jié)理面抗拉強(qiáng)度為10 MPa；②模型中設(shè)置節(jié)理面，節(jié)理面抗拉強(qiáng)度為3 MPa。

2.2 震動(dòng)波傳播應(yīng)力響應(yīng)過程

數(shù)值模型中的應(yīng)力響應(yīng)可以真實(shí)地反映震動(dòng)波傳播過程，選取震動(dòng)波作用1個(gè)完整周期內(nèi)不同時(shí)刻的應(yīng)力進(jìn)行分析，節(jié)理面抗拉強(qiáng)度為10 MPa時(shí)的垂直應(yīng)力如圖5所示。

圖5 節(jié)理面抗拉強(qiáng)度10 MPa下垂直應(yīng)力Fig.5 Vertical stress of joint suface under tensile strength 10 MPa

在原巖梯度應(yīng)力作用下，初始模型整體處于壓應(yīng)力狀態(tài)，應(yīng)力均為負(fù)值且絕對值隨著監(jiān)測點(diǎn)高度H的增加而減小。震動(dòng)波施加于模型頂部后，模型上部的應(yīng)力狀態(tài)首先發(fā)生改變，應(yīng)力降低，表明初始垂直應(yīng)力與震動(dòng)波所傳遞的拉伸應(yīng)力出現(xiàn)疊加，對原巖應(yīng)力有較大影響。隨著作用時(shí)間增加，震動(dòng)波自上而下進(jìn)行傳播，模型中、下部的應(yīng)力依次出現(xiàn)波動(dòng)，震動(dòng)波自施加邊界向外部擴(kuò)散。半正弦震動(dòng)波作用1個(gè)周期內(nèi)，模型中的垂直應(yīng)力呈現(xiàn)“負(fù)-正-負(fù)”變化特征，應(yīng)力狀態(tài)也隨之顯著改變，由初始受壓變?yōu)槭芾詈笥只貧w到受壓狀態(tài)。

震動(dòng)波傳播過程中測點(diǎn)處的應(yīng)力、速度變化曲線如圖6所示。模型各點(diǎn)處的應(yīng)力、速度均與時(shí)間呈半正弦函數(shù)關(guān)系，與震動(dòng)波呈相同的波動(dòng)特征。但各測點(diǎn)的應(yīng)力和速度響應(yīng)也存在差異性，在震動(dòng)波作用下，各測點(diǎn)應(yīng)力及速度響應(yīng)過程不同，距震動(dòng)波作用處越遠(yuǎn)，應(yīng)力或速度發(fā)生改變的時(shí)刻越晚，具有明顯的傳遞特征。當(dāng)節(jié)理面強(qiáng)度較高時(shí)，巖體并沒有發(fā)生破壞，震動(dòng)波在全模型中傳播，作用1個(gè)周期后，各點(diǎn)處的速度重新恢復(fù)為0，應(yīng)力狀態(tài)也基本復(fù)原。

圖6 節(jié)理面抗拉強(qiáng)度10 MPa下測點(diǎn)應(yīng)力及速度變化Fig.6 Simulated stress and velocity of measuring point with tensile strength 10 MPa joint suface

2.3 不連續(xù)面的應(yīng)力阻斷機(jī)制

節(jié)理面抗拉強(qiáng)度為3 MPa時(shí)的垂直應(yīng)力如圖7所示。對比圖5可知，震動(dòng)波作用時(shí)間為0.01～0.04 s時(shí)，不同抗拉強(qiáng)度巖體中垂直應(yīng)力分布情況相同，主要是由于震動(dòng)初期震動(dòng)波應(yīng)力較低，巖體的應(yīng)力狀態(tài)由初始應(yīng)力與震動(dòng)波所傳遞的應(yīng)力共同決定，此階段節(jié)理面上分布的拉應(yīng)力較小，節(jié)理面未發(fā)生張拉破壞。當(dāng)作用時(shí)間為0.05 s時(shí)，震動(dòng)波應(yīng)力達(dá)到7 MPa，應(yīng)力疊加使得模型上部應(yīng)力狀態(tài)由初始4 MPa左右的壓應(yīng)力轉(zhuǎn)變?yōu)榧s3 MPa的拉應(yīng)力，上分層節(jié)理面的拉應(yīng)力超過其抗拉強(qiáng)度時(shí)，節(jié)理發(fā)生張拉破壞，垂直方向上缺乏來自頂部的約束，中、下部巖體中的應(yīng)力迅速得到釋放，巖體中部應(yīng)力小于1 MPa。此時(shí)上分層已經(jīng)脫離下部模型，震動(dòng)波傳遞的拉應(yīng)力只作用于上分層，導(dǎo)致上分層測點(diǎn)的質(zhì)點(diǎn)震動(dòng)速度在拉應(yīng)力持續(xù)作用下出現(xiàn)突然升高，如圖8所示。隨著震動(dòng)波作用時(shí)間進(jìn)一步增加，震動(dòng)波應(yīng)力逐漸下降，逐漸恢復(fù)至壓應(yīng)力狀態(tài)，但節(jié)理破壞具有不可逆性，因此僅模型上分層中的應(yīng)力狀態(tài)發(fā)生變化，其下方各點(diǎn)不再受震動(dòng)波的影響，不連續(xù)面阻斷了震動(dòng)波的應(yīng)力傳遞過程。

圖7 節(jié)理面抗拉強(qiáng)度3 MPa下垂直應(yīng)力Fig.7 Vertical stress of joint surface under tensile strength 3 MPa

圖8 節(jié)理面抗拉強(qiáng)度3 MPa測點(diǎn)應(yīng)力及速度Fig.8 Simulated stress and velocity of measuring point under tensile strength 3 MPa of joint surface

基于上述規(guī)律，可以得出礦震震動(dòng)波傳播過程中不連續(xù)面的應(yīng)力阻斷機(jī)制：不連續(xù)面為巖體中的弱面，通常最容易發(fā)生破壞。當(dāng)震動(dòng)波傳遞的應(yīng)力與不連續(xù)面處的原巖應(yīng)力相互疊加超過不連續(xù)面強(qiáng)度時(shí)，不連續(xù)面發(fā)生破壞，進(jìn)而產(chǎn)生位移，介質(zhì)不連續(xù)性阻斷了應(yīng)力波的傳播。

3 礦震傳播能量衰減規(guī)律與實(shí)測

為驗(yàn)證礦震震動(dòng)波傳播數(shù)值模擬結(jié)果，討論工程現(xiàn)場中危險(xiǎn)性礦震的判別方法，采用理論分析和現(xiàn)場實(shí)測的方法，從能量角度出發(fā)，進(jìn)一步研究了礦區(qū)內(nèi)震動(dòng)波傳播過程中的衰減規(guī)律。

3.1 震動(dòng)波能量衰減規(guī)律

煤礦礦震對井下工作面或巷道破壞程度與震動(dòng)波在煤巖體中的衰減特性密切相關(guān)。礦震震動(dòng)波傳播受傳播距離、煤巖體介質(zhì)、不連續(xù)面作用以及震源破裂特征等眾多因素影響，采動(dòng)作用下震源破裂模式和尺度等震源特征異常復(fù)雜，具有很大的不確定性。為便于統(tǒng)計(jì)分析，將礦震簡化為點(diǎn)震源，震動(dòng)波從震源到達(dá)接收點(diǎn)時(shí)，其震動(dòng)速度的幅值會(huì)發(fā)生變化，其變化主要來源于波前的擴(kuò)散和介質(zhì)的吸收，前者與距離r成反比，后者與距離r呈負(fù)指數(shù)的關(guān)系，震動(dòng)波傳播過程中微震監(jiān)測系統(tǒng)臺站所記錄的質(zhì)點(diǎn)震動(dòng)速度衰減規(guī)律為

Ai=C1e-αri/ri

(2)

式中：C1為震源處質(zhì)點(diǎn)震動(dòng)速度；α為速度吸收系數(shù)；ri為震源至臺站i的距離；Ai為臺站i記錄的幅值。

對任意彈性媒介，波在彈性媒介中傳播時(shí)，各質(zhì)點(diǎn)在其平衡位置附近振動(dòng)，從而使彈性媒介具有動(dòng)能；又因?yàn)檎駝?dòng)造成了各質(zhì)點(diǎn)之間的距離發(fā)生改變，媒介發(fā)生形變，從而使之具有形變勢能。由于煤巖為脆性材料，震動(dòng)波在其中傳播時(shí)所能引起的煤巖體形變勢能改變量較小，可認(rèn)為動(dòng)能是礦震震動(dòng)波能量傳遞的主要部分。速度與動(dòng)能Ek,i的關(guān)系見式(3)，結(jié)合震動(dòng)速度衰減規(guī)律，可以得出震動(dòng)波傳播過程中能量衰減規(guī)律為

(3)

式中：Ek,i為臺站i所記錄波形的動(dòng)能；m為震動(dòng)巖體的質(zhì)量；Ai為臺站i記錄的幅值。

聯(lián)立式(2)和式(3)，得：

Ek,i=C2ri-2e-βri，β=2α

(4)

式中：C2為震源震動(dòng)能量；ri為震源至臺站i的距離；β為能量吸收系數(shù)。

基于微震監(jiān)測數(shù)據(jù)，通過擬合可以確定礦井某個(gè)區(qū)域一定時(shí)期內(nèi)質(zhì)點(diǎn)震動(dòng)速度衰減規(guī)律，根據(jù)速度吸收系數(shù)α計(jì)算出能量吸收系數(shù)β，結(jié)合微震監(jiān)測系統(tǒng)中所記錄的震源能量，即可得到礦震震動(dòng)波傳播過程能量衰減規(guī)律。由于不連續(xù)面對震動(dòng)波傳播有重要作用，擬合應(yīng)分為2種條件：未開采的實(shí)體煤區(qū)域和裂縫場較發(fā)育的采空區(qū)區(qū)域。

3.2 現(xiàn)場實(shí)測

某礦主采煤層為3煤層，平均厚度8.4 m，開采深度為530～780 m。采區(qū)內(nèi)已連續(xù)回采3個(gè)相鄰工作面，形成較大的采空區(qū)，其余為未開采的實(shí)體煤區(qū)域，為監(jiān)測工作面回采過程中礦震發(fā)生情況，沿采空區(qū)周圍均勻布置了微震監(jiān)測臺站，6306工作面位置及微震布設(shè)如圖9所示。

圖9 工作面及微震臺站布置Fig.9 Layout of working face and dispersion of micro-seismic stations

基于上述分析，選取6306工作面不同開采時(shí)期的微震數(shù)據(jù)，分別對實(shí)體煤區(qū)域和采空區(qū)中的礦震震動(dòng)速度衰減規(guī)律進(jìn)行擬合分析，為減小差異性，微震事件震源位置與能量大小應(yīng)相近，震源位置及選用的臺站分布如圖10所示。

圖10 震源位置及臺站選取Fig.10 Location of seismic source and selection of stations

圖11、圖12分別為實(shí)體煤區(qū)域、采空區(qū)礦震震動(dòng)速度衰減擬合結(jié)果，隨著傳播距離增加，震動(dòng)速度呈負(fù)指數(shù)型衰減。礦震在實(shí)體煤區(qū)域傳播時(shí)，震動(dòng)速度吸收系數(shù)平均值為0.003 75；當(dāng)震動(dòng)波經(jīng)過采空區(qū)時(shí)，震動(dòng)速度吸收系數(shù)升高，平均值為0.007 6，表明不連續(xù)面對震動(dòng)波的衰減作用明顯，驗(yàn)證了數(shù)值模擬結(jié)果的正確性。

圖11 實(shí)體煤區(qū)域震動(dòng)速度衰減Fig.11 Attenuation of seismic velocity in solid coal area

圖12 采空區(qū)區(qū)域震動(dòng)速度衰減Fig.12 Attenuation of vibration velocity in gob

根據(jù)震動(dòng)速度吸收系數(shù)與能量吸收系數(shù)之間的關(guān)系，可以分別得到實(shí)體煤、采空區(qū)條件下礦震震動(dòng)波能量吸收系數(shù)分別為0.007 5、0.015 2。如果已知礦震震源能量、震動(dòng)波能量吸收系數(shù)和巷道抗沖擊能力Emax，根據(jù)式(5)即可求解得出安全距離rs，將實(shí)際震源至巷道的距離與安全距離進(jìn)行比較，如果實(shí)際距離小于安全距離，表明對井下而言，此次礦震具有危險(xiǎn)性。

Emax=C2rs-2e-βrs

(5)

式中:Emax為巷道抗沖擊能力;rs為安全距離。

4 結(jié) 論

1)綜合分析礦震震源信息、波形特征和井下礦壓顯現(xiàn)情況，提出將礦震分為采動(dòng)破裂型、巨厚覆巖型和高能震動(dòng)型3種類型。根據(jù)礦震對采場及巷道是否具有危險(xiǎn)性將其分為正常礦震和危險(xiǎn)性礦震。

2)采用UDEC數(shù)值模擬再現(xiàn)了節(jié)理巖體中礦震震動(dòng)波傳播過程，揭示了煤巖體中不連續(xù)面對礦震震動(dòng)波傳播的應(yīng)力阻斷機(jī)制：當(dāng)震動(dòng)波傳遞的應(yīng)力與原巖應(yīng)力疊加超過不連續(xù)面強(qiáng)度時(shí)，不連續(xù)面失效產(chǎn)生位移，介質(zhì)不連續(xù)性阻斷了震動(dòng)波應(yīng)力傳播。

3)以某礦為背景，實(shí)測得到礦區(qū)實(shí)體煤區(qū)域和采空區(qū)區(qū)域的質(zhì)點(diǎn)震動(dòng)速度吸收系數(shù)分別為0.0037 5和0.007 6，能量吸收系數(shù)分別為0.007 5、0.015 2，2種區(qū)域內(nèi)震動(dòng)波衰減差異明顯。

4)提出了基于微震監(jiān)測數(shù)據(jù)進(jìn)行礦震震動(dòng)波能量衰減規(guī)律分析的方法，綜合分析礦震位置及能量、巷道抗沖擊能力和震動(dòng)波能量衰減規(guī)律，求解得出安全距離，將實(shí)際震源至巷道的距離與安全距離進(jìn)行比較，如果實(shí)際距離小于安全距離，表明對井下而言，此次礦震具有危險(xiǎn)性。提出的礦震類型、震動(dòng)波傳播規(guī)律以及危險(xiǎn)性礦震判別準(zhǔn)則，對煤礦礦震的認(rèn)知、防控及相關(guān)災(zāi)害治理具有一定的參考價(jià)值。