趙奎，　胡源，　王曉軍，2，　馮蕭，　卓毓龍

（1.江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州　341000；2.西部礦業(yè)博士后科研工作站，西寧810006）

采場(chǎng)頂板穩(wěn)定性位移與聲發(fā)射組合識(shí)別研究

趙奎1，胡源1，王曉軍1，2，馮蕭1，卓毓龍1

（1.江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州341000；2.西部礦業(yè)博士后科研工作站，西寧810006）

地下工程頂板變形破壞是一個(gè)復(fù)雜的非線性過(guò)程，單一監(jiān)測(cè)結(jié)果難以對(duì)其穩(wěn)定性狀態(tài)準(zhǔn)確判別.基于數(shù)理統(tǒng)計(jì)的方法推導(dǎo)了頂板變形過(guò)程中位移監(jiān)測(cè)樣本的標(biāo)準(zhǔn)差，根據(jù)結(jié)果分別表征了無(wú)變形、穩(wěn)定變形和加速變形3個(gè)頂板變形階段.結(jié)合當(dāng)前巖體破壞聲發(fā)射特征參數(shù)分析試驗(yàn)結(jié)果，對(duì)頂板巖體變形過(guò)程的位移與聲發(fā)射監(jiān)測(cè)結(jié)果進(jìn)行組合分析，提出了頂板穩(wěn)定性識(shí)別的4種模式.結(jié)合某礦工程實(shí)例，對(duì)回采過(guò)程頂板破壞進(jìn)行了準(zhǔn)確預(yù)測(cè).研究結(jié)果為地下工程頂板穩(wěn)定性監(jiān)測(cè)結(jié)果分析提供了數(shù)據(jù)處理新方法，加強(qiáng)了頂板穩(wěn)定性識(shí)別的準(zhǔn)確性.

標(biāo)準(zhǔn)差；加速變形；聲發(fā)射；頂板穩(wěn)定性

在上覆巖層壓力及開(kāi)挖力學(xué)效應(yīng)影響下，地下采場(chǎng)及硐室頂板是整個(gè)工程中相對(duì)薄弱的部分，當(dāng)開(kāi)挖空間、應(yīng)力分布等發(fā)生變化時(shí)，極有可能誘發(fā)其失穩(wěn)破壞，對(duì)正常生產(chǎn)及生命安全形成極大威脅[1-3].據(jù)統(tǒng)計(jì)數(shù)據(jù)分析，地下工程作業(yè)過(guò)程中，頂板事故是各種安全事故中最為突出的一類(lèi).因此，頂板穩(wěn)定性的識(shí)別預(yù)測(cè)對(duì)地下工程安全作業(yè)有著極為重要的意義.

目前對(duì)頂板穩(wěn)定性測(cè)試分析的方法較多，監(jiān)測(cè)方法主要包括變形觀測(cè)法、聲發(fā)射法、位移監(jiān)測(cè)分析法等[4-6]，數(shù)據(jù)分析方法多采用經(jīng)驗(yàn)、解析和半定量的方法[7-10].但頂板變形破壞是一個(gè)復(fù)雜的非線性過(guò)程，從承載到失穩(wěn)整個(gè)過(guò)程中所釋放的信息較多，僅通過(guò)某種監(jiān)測(cè)技術(shù)和分析難以全面反映頂板的應(yīng)力、應(yīng)變分布和破壞情況，故對(duì)頂板穩(wěn)定性評(píng)價(jià)也存在較大的局限性和片面性.為此，針對(duì)頂板巖體變形過(guò)程所釋放的最為直接的多方面信息，利用有效的數(shù)學(xué)方法分別對(duì)其進(jìn)行量化提取并組合分析，以實(shí)現(xiàn)地下工程頂板穩(wěn)定性的準(zhǔn)確評(píng)價(jià)與預(yù)測(cè).

1　頂板變形位移識(shí)別

地下工程頂板一旦被揭露，即發(fā)生相應(yīng)的變形，該時(shí)段的變形稱(chēng)為開(kāi)挖變形，隨著開(kāi)挖變形結(jié)束，對(duì)不同的巖體而言，頂板變形增量逐漸減小或消失，其變形處于穩(wěn)定階段，當(dāng)巖體處于非穩(wěn)定狀態(tài)時(shí)，其變形增量又逐步增大，該時(shí)段稱(chēng)為加速變形階段，直至頂板破壞失穩(wěn).因此，頂板巖體在開(kāi)挖完成后整個(gè)承載過(guò)程中可能出現(xiàn)的位移變化階段主要包括無(wú)變形或穩(wěn)定變形階段、加速變形階段.

1.1無(wú)變形階段

該階段說(shuō)明開(kāi)挖變形完成后不產(chǎn)生任何時(shí)效變形，如果位移值發(fā)生細(xì)微變化說(shuō)明是由于測(cè)量?jī)x器的誤差而引起的，因此，實(shí)際監(jiān)測(cè)的位移值的變化就是儀器測(cè)量誤差值.即：

其中di為每次監(jiān)測(cè)的位移值，δi為測(cè)量誤差，測(cè)量誤差δi由儀器的精密程度引起，屬于隨機(jī)誤差，通常是由多種相互獨(dú)立的因素造成的許多微小誤差的總和.其變化服從正態(tài)分布N（μ，σ）.共監(jiān)測(cè)n組數(shù)據(jù)后，測(cè)量值的標(biāo)準(zhǔn)差和平均值為：

進(jìn)一步推得：

式（5）中：

此處δi服從正態(tài)分布N （μ，σ），所以δ1+δ2+…+δn= nμ，ˉ=μ，由此，對(duì)式（6）進(jìn)行整理：

由于δi服從正態(tài)分布N（μ，σ），對(duì)其進(jìn)行標(biāo)準(zhǔn)化變換，得到服從標(biāo)準(zhǔn)正態(tài)分布N（0，1），進(jìn)行推導(dǎo)可得：

將式（8）代入式（7）得：

由于Ti服從標(biāo)準(zhǔn)正態(tài)分布N（0，1），故T1+T2+…+ Ti=0，而服從自由度為n的卡方分布，即：

其均值E（χ2）=n，所以將以上結(jié)果代入式（2），整理后可得測(cè)量值的標(biāo)準(zhǔn)差：

式（11）中σ為測(cè)量誤差的標(biāo)準(zhǔn)差，即測(cè)量?jī)x器的精度，n為測(cè)量次數(shù)，可見(jiàn)隨著測(cè)量次數(shù)增多，測(cè)量值的標(biāo)準(zhǔn)差逐漸增大，所以應(yīng)該盡量提高測(cè)量精度，以降低測(cè)量中的隨機(jī)誤差.同時(shí)，由于式（11）推導(dǎo)過(guò)程中測(cè)量值的標(biāo)準(zhǔn)差計(jì)算中只考慮了測(cè)量?jī)x器誤差造成的變化，而將該階段頂板的實(shí)際變形值定為零.因此，如果頂板處于無(wú)變形階段，則：

1.2穩(wěn)定變形階段

該階段表明頂板產(chǎn)生了變形，因此該階段實(shí)際監(jiān)測(cè)的位移值包含2個(gè)部分，一部分為頂板變形產(chǎn)生的位移，另一部分為測(cè)量?jī)x器引起的誤差值.即：

式（13）中wi為頂板實(shí)際產(chǎn)生的位移值，有關(guān)測(cè)量?jī)x器引起的誤差值δi前面已經(jīng)論述，此處僅討論頂板實(shí)際產(chǎn)生的位移值.

其均值為：

其標(biāo)準(zhǔn)差為：

如果頂板巖體處于穩(wěn)定變形階段，則可以近似認(rèn)為在等時(shí)間間隔內(nèi)頂板實(shí)際產(chǎn)生的位移值的變化率相等，即等時(shí)間間隔內(nèi)頂板位移值的增量一致.

w0為開(kāi)始監(jiān)測(cè)之前的位移值，取值為0，由此可知：

將式（17）、式（18）、式（19）代入式（16）整理得：

將式（20）代入式（15）得頂板實(shí)際位移值的標(biāo)準(zhǔn)差：

由于t為穩(wěn)定變形階段位移的增量，監(jiān)測(cè)過(guò)程很難準(zhǔn)確獲得其實(shí)際值，故將其轉(zhuǎn)化為最終監(jiān)測(cè)的位移值，即：S=nt，S為頂板實(shí)際變化位移值的總和.故式（21）轉(zhuǎn)化為：

根據(jù)式（13）將測(cè)量?jī)x器引起的誤差值考慮在內(nèi)，得到穩(wěn)定變形階段實(shí)際監(jiān)測(cè)頂板位移值的標(biāo)準(zhǔn)差為：

由此可以得到，在穩(wěn)定變形階段：

1.3加速變形階段

當(dāng)頂板巖體處于加速變形階段時(shí)，可近似認(rèn)為在等時(shí)間間隔內(nèi)頂板實(shí)際產(chǎn)生的位移值的變化率逐漸增大，即在相等時(shí)間間隔內(nèi)，頂板位移增量逐漸變大.

根據(jù)1.2節(jié)推導(dǎo)過(guò)程，可最終得到在加速變形階段：

通過(guò)對(duì)頂板位移監(jiān)測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，可以判別頂板變形階段，但整個(gè)地下工程施工過(guò)程中，受復(fù)雜因素的影響，頂板應(yīng)力狀態(tài)不斷變化，頂板3種變形階段也可以相互轉(zhuǎn)變，換言之，即使進(jìn)入加速變形階段，頂板也不一定失穩(wěn).僅憑位移監(jiān)測(cè)難以分析頂板的穩(wěn)定性狀態(tài).

2　頂板失穩(wěn)聲發(fā)射指標(biāo)識(shí)別

巖體變形破壞過(guò)程伴隨裂紋產(chǎn)生和擴(kuò)展，裂紋產(chǎn)生與擴(kuò)展又伴隨著應(yīng)變能的釋放，即巖體破壞聲發(fā)射現(xiàn)象.針對(duì)巖體破壞的聲發(fā)射特征，國(guó)內(nèi)外學(xué)者進(jìn)行了大量的研究，通過(guò)室內(nèi)試驗(yàn)及現(xiàn)場(chǎng)測(cè)試，揭示了巖石破壞過(guò)程與聲發(fā)射特征參數(shù)的內(nèi)在關(guān)系.

2.1巖石承載過(guò)程聲發(fā)射關(guān)聯(lián)分維數(shù)變化規(guī)律

筆者及相關(guān)研究發(fā)現(xiàn)巖石破壞過(guò)程聲發(fā)射特征參數(shù)具有明顯分形特征[11]，針對(duì)這一研究結(jié)果，筆者在聲發(fā)射測(cè)量地應(yīng)力的試驗(yàn)中發(fā)現(xiàn)[12]：不同巖性的多組巖石試樣在單軸加載全過(guò)程中，其聲發(fā)射關(guān)聯(lián)分維數(shù)總體在開(kāi)始階段和中間階段，聲發(fā)射關(guān)聯(lián)分維數(shù)總體呈上升趨勢(shì)，但在上升過(guò)程中也出現(xiàn)了不同程度的多次起伏，即在加載過(guò)程中多次出現(xiàn)關(guān)聯(lián)分維數(shù)的上升和下降.當(dāng)加載至巖石峰值應(yīng)力前夕（即將破壞），聲發(fā)射關(guān)聯(lián)分維數(shù)呈現(xiàn)降低的特征，圖1是其試驗(yàn)結(jié)果曲線圖.

2.2巖石承載過(guò)程聲發(fā)射長(zhǎng)程相關(guān)指數(shù)變化規(guī)律

許福樂(lè)等[13]提出了用聲發(fā)射序列的長(zhǎng)程波動(dòng)來(lái)描述巖石破壞過(guò)程，通過(guò)不同煤巖試樣單軸加載全過(guò)程聲發(fā)射試驗(yàn)，分析了整個(gè)過(guò)程聲發(fā)射強(qiáng)度序列的長(zhǎng)程相關(guān)性，研究了長(zhǎng)程相關(guān)指數(shù)H的變化規(guī)律，試驗(yàn)結(jié)果表明，在加載初期，H值增長(zhǎng)緩慢，隨著加載繼續(xù)進(jìn)行，H值不斷增大，并達(dá)到最大值，隨后又出現(xiàn)降低.但在整個(gè)過(guò)程中多次出現(xiàn)升高-降低的模式，直至加載至巖石達(dá)到加速非線性變形階段 （瀕臨破壞前夕），聲發(fā)射強(qiáng)度長(zhǎng)程相關(guān)指數(shù)H值出現(xiàn)了“最大-減小”模式，由此得出結(jié)論可以利用巖石破壞聲發(fā)射強(qiáng)度長(zhǎng)程相關(guān)指數(shù)H值的“最大-減小”模式判別巖體的失穩(wěn)，見(jiàn)圖2.

2.3巖體破壞聲發(fā)射相對(duì)強(qiáng)弱指標(biāo)變化規(guī)律

王寧等[14]研究了評(píng)價(jià)采場(chǎng)巖體破壞的聲發(fā)射相對(duì)強(qiáng)弱量化指標(biāo)，對(duì)現(xiàn)場(chǎng)監(jiān)測(cè)的聲發(fā)射能率、總事件、大事件進(jìn)行標(biāo)準(zhǔn)化處理，得到聲發(fā)射事件率的相對(duì)強(qiáng)弱指標(biāo)CR和聲發(fā)射能率相對(duì)強(qiáng)弱指標(biāo)ER，從圖3可以看出，在整個(gè)巖體承載聲發(fā)射監(jiān)測(cè)過(guò)程中，CR和ER總體呈現(xiàn)波動(dòng)趨勢(shì)，但在失穩(wěn)前兆期，相對(duì)強(qiáng)弱指標(biāo)出現(xiàn)迅速回落.

圖1　不同巖樣單軸加載過(guò)程聲發(fā)射關(guān)聯(lián)分維數(shù)變化曲線Fig.1　Relevant dim ensions change curve of different rock specim ens during loading p rocess

圖2　煤巖試樣單軸加載過(guò)程聲發(fā)射強(qiáng)度長(zhǎng)程相關(guān)指數(shù)H值變化曲線Fig.2　Variation curves of long-term correlation exponent H of acoustic em ission of coal-rock during loading p rocess

圖3　巖體承載過(guò)程聲發(fā)射CR/ER變化曲線Fig.3　The change curve of acoustic em ission CR/ER of rock m ass during loading p rocess

從筆者[15]及相關(guān)人員對(duì)巖體承載聲發(fā)射特性參數(shù)研究結(jié)果來(lái)看[16]，都說(shuō)明了同樣一點(diǎn)，巖體加速破壞時(shí)，聲發(fā)射監(jiān)測(cè)指標(biāo)頻度與量度激增，巖體在失穩(wěn)前夕，聲發(fā)射監(jiān)測(cè)指標(biāo)頻度與量度出現(xiàn)了迅速回落，對(duì)應(yīng)于聲發(fā)射b值的升降模式，聲發(fā)射長(zhǎng)程相關(guān)指數(shù)H值的最大-減小模式，以及聲發(fā)射相對(duì)強(qiáng)弱指標(biāo)CR、ER的起伏模式.因此巖體失穩(wěn)前，其聲發(fā)射監(jiān)測(cè)指標(biāo)頻度（量度）出現(xiàn)激增-回落現(xiàn)象.但反之并不成立，從多次研究結(jié)果的變化曲線上來(lái)分析，在整個(gè)監(jiān)測(cè)過(guò)程中，這種激增-回落現(xiàn)象在變化曲線中多次出現(xiàn)，也就是說(shuō)這種激增-回落現(xiàn)象并不能說(shuō)明巖體一定失穩(wěn)，所以單純利用聲發(fā)射監(jiān)測(cè)指標(biāo)的這一性質(zhì)并不能準(zhǔn)確判斷其穩(wěn)定性狀態(tài).

3　位移與聲發(fā)射組合模式

位移與聲發(fā)射是巖體穩(wěn)定性監(jiān)測(cè)中最有效、最直觀的伴生信息，上述分析可知單一方法無(wú)法準(zhǔn)確判定頂板的穩(wěn)定性狀態(tài)，對(duì)其進(jìn)行組合分析可以增加分析結(jié)果的準(zhǔn)確性.對(duì)頂板不同的變形階段和聲發(fā)射信息指標(biāo)進(jìn)行組合分析，可以確定以下幾種組合模式.

1）無(wú)變形穩(wěn)定模式.該模式下監(jiān)測(cè)周期內(nèi)，頂板位移監(jiān)測(cè)結(jié)果分析滿(mǎn)足式 （12），頂板處于無(wú)變形階段，因此頂板不會(huì)破壞，也不會(huì)產(chǎn)生相應(yīng)大量激增的聲發(fā)射現(xiàn)象，即使有也是局部巖體受結(jié)構(gòu)控制產(chǎn)生的小規(guī)模破裂.據(jù)此，該位移與聲發(fā)射的組合模式可以分解為以下2種，見(jiàn)式（27）.均說(shuō)明頂板處于穩(wěn)定狀態(tài).

2）平穩(wěn)變形穩(wěn)定模式.該模式下整個(gè)監(jiān)測(cè)周期內(nèi)頂板位移監(jiān)測(cè)結(jié)果分析滿(mǎn)足式（24），說(shuō)明頂板處于平穩(wěn)變形階段，該階段內(nèi)由于巖體產(chǎn)生了實(shí)際變形，變形結(jié)果可能形成2種情況，其一為巖體只發(fā)生變形并未產(chǎn)生內(nèi)部結(jié)構(gòu)的破裂，因此聲發(fā)射相對(duì)“平靜”；其二為巖體變形引發(fā)內(nèi)部結(jié)構(gòu)破壞，產(chǎn)生一定量的聲發(fā)射現(xiàn)象.據(jù)此，該位移與聲發(fā)射的組合模式也可分為2種情況，見(jiàn)式（28）.該組合分析模式也表明頂板處于穩(wěn)定狀態(tài).

3）加速變形穩(wěn)定模式.該模式下整個(gè)監(jiān)測(cè)周期內(nèi)頂板位移監(jiān)測(cè)結(jié)果分析滿(mǎn)足式（26），說(shuō)明頂板實(shí)際已由平穩(wěn)變形階段過(guò)渡到加速變形階段，由于頂板巖體的加速變形，導(dǎo)致頂板巖體內(nèi)部產(chǎn)生一定數(shù)量的新裂紋，并隨著變形進(jìn)一步增加，巖體內(nèi)部裂紋逐步擴(kuò)展，聲發(fā)射實(shí)際監(jiān)測(cè)結(jié)果也體現(xiàn)為其能量和頻度等特征參量激增，但此時(shí)只是加速變形階段巖體內(nèi)部結(jié)構(gòu)在破壞，而頂板依然可以繼續(xù)承載，尚能保持穩(wěn)定，組合模式見(jiàn)式（29）.

4）加速變形失穩(wěn)模式.該模式下整個(gè)監(jiān)測(cè)周期內(nèi)頂板位移監(jiān)測(cè)結(jié)果分析滿(mǎn)足式（26），頂板處于加速變形階段，同時(shí)聲發(fā)射特征參量分析結(jié)果出現(xiàn)激增-回落現(xiàn)象.說(shuō)明此時(shí)頂板巖體大部分已經(jīng)破壞，內(nèi)部裂紋已經(jīng)出現(xiàn)貫通融合，致使聲發(fā)射特征參量頻度（量度）回落.該組合分析模式表明頂板巖體已進(jìn)入整體失穩(wěn)階段，見(jiàn)式（30），頂板處于失穩(wěn)前兆期.

4　應(yīng)用實(shí)例

某銅礦對(duì)階段水平礦柱（頂?shù)字?shí)施回采[15]，水平礦柱上方為上部中段松散塊石充填體，回采過(guò)程主要采用分層進(jìn)路式上采方式，在上部散體壓力作用下，進(jìn)路頂板的穩(wěn)定性是安全回采的關(guān)鍵.為此，對(duì)進(jìn)路頂板實(shí)施位移與聲發(fā)射組合監(jiān)測(cè)，柱17#采場(chǎng)7#進(jìn)路回采過(guò)后在進(jìn)路中布置頂板沉降儀監(jiān)測(cè)頂板位移，同時(shí)利用水平礦柱上部中段坑道在頂板上布置鉆孔，安裝聲發(fā)射連續(xù)監(jiān)測(cè)儀，測(cè)點(diǎn)布置見(jiàn)圖4.

圖4　監(jiān)測(cè)儀器布置Fig.4　Layout ofmonitoring equipm ent

從儀器安裝至進(jìn)路充填共持續(xù)監(jiān)測(cè)41 d，取得位移監(jiān)測(cè)數(shù)據(jù)15組，聲發(fā)射監(jiān)測(cè)處理數(shù)據(jù)40組.見(jiàn)

表1　頂板位移監(jiān)測(cè)結(jié)果Table1　Results of roof of disp lacem entm onitoring

表2　頂板聲發(fā)射監(jiān)測(cè)處理結(jié)果Table 2　Resu lts of roof disp lacem ent acoustic em ission m onitoring

利用第1節(jié)中的處理方法對(duì)位移監(jiān)測(cè)值進(jìn)行分析，由于D2未能取得完整監(jiān)測(cè)結(jié)果，分別計(jì)算D1、D3測(cè)點(diǎn)位移監(jiān)測(cè)值的標(biāo)準(zhǔn)差.

分別計(jì)算出D1、D3測(cè)點(diǎn)無(wú)變形階段和穩(wěn)定變形階段的標(biāo)準(zhǔn)差上限值，測(cè)量?jī)x器的精度σ為0.01.

根據(jù)計(jì)算結(jié)果可知：

由此可知D1測(cè)點(diǎn)所在區(qū)域頂板已經(jīng)處于加速變形階段，而D3測(cè)點(diǎn)所在區(qū)域頂板尚處于穩(wěn)定變形階段.頂板是否出現(xiàn)失穩(wěn)前兆應(yīng)配合聲發(fā)射測(cè)點(diǎn)的監(jiān)測(cè)結(jié)果進(jìn)行組合分析.

對(duì)聲發(fā)射能率監(jiān)測(cè)結(jié)果繪制變化曲線，見(jiàn)圖5.

圖5　聲發(fā)射特征參數(shù)變化曲線Fig.5　Change curve of acoustic em ission parameters

從圖5可知，在整個(gè)監(jiān)測(cè)周期內(nèi)，聲發(fā)射能率先后出現(xiàn)了3次較大的起伏，其中尤以第3次出現(xiàn)了聲發(fā)射的激增-回落，配合該監(jiān)測(cè)期內(nèi)頂板位移監(jiān)測(cè)結(jié)果，發(fā)現(xiàn)截止至12月23日，頂板D1測(cè)點(diǎn)區(qū)域已經(jīng)處于加速變形階段，符合位移與聲發(fā)射組合分析的加速變形不穩(wěn)定模式，由此判斷D1測(cè)點(diǎn)附近區(qū)域頂板已經(jīng)處于失穩(wěn)前兆期，遂撤出該進(jìn)路的設(shè)備及監(jiān)測(cè)儀器，12月27日早班發(fā)現(xiàn)該進(jìn)路端頭與上盤(pán)巖體接觸帶附近頂板巖體小規(guī)?？逅?7日小夜班對(duì)該進(jìn)路實(shí)施充填.

5　結(jié)　論

1）頂板位移監(jiān)測(cè)數(shù)據(jù)數(shù)理統(tǒng)計(jì)分析結(jié)果表明：地下工程開(kāi)挖后，頂板位移主要經(jīng)歷3種變形過(guò)程，即無(wú)變形過(guò)程、穩(wěn)定變形過(guò)程與加速變形過(guò)程.在整個(gè)地下工程的復(fù)雜施工過(guò)程中，3個(gè)過(guò)程可以相互過(guò)渡轉(zhuǎn)換.

2）巖石承載過(guò)程多次聲發(fā)射特性試驗(yàn)結(jié)果說(shuō)明，巖體失穩(wěn)前，其聲發(fā)射監(jiān)測(cè)指標(biāo)頻度（量度）出現(xiàn)激增-回落現(xiàn)象.但試驗(yàn)全過(guò)程顯示該現(xiàn)象出現(xiàn)并不能預(yù)測(cè)巖體必然失穩(wěn).

3）僅憑位移或聲發(fā)射單一監(jiān)測(cè)技術(shù)并不能對(duì)頂板巖體穩(wěn)定性做出準(zhǔn)確預(yù)判，對(duì)2種測(cè)試結(jié)果進(jìn)行組合分析得到4種組合模式，其中加速變形失穩(wěn)模式是頂板巖體失穩(wěn)的真正前兆.

[1]楊宇江，高晗，李元輝.附加載荷作用下空區(qū)頂板穩(wěn)定性分析[J].東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，32（9）：1332-1335.

[2]Swift G M，Reddish D J.St ability problems associated with an abandoned ironstone mine[J].Bulletin of Engineering Geology and the Environment，2002，61（3）：227-239.

[3]朱維申，何滿(mǎn)潮.復(fù)雜條件下圍巖穩(wěn)定性與巖體動(dòng)態(tài)施工力學(xué)[M].北京：科學(xué)出版社，1996.

[4]王恩元，賈慧霖，李忠輝，等.用電磁輻射法監(jiān)測(cè)預(yù)報(bào)礦山采空區(qū)頂板穩(wěn)定性[J].煤炭學(xué)報(bào)，2006，31（1）：16-19.

[5]文光才，李建功，鄒銀輝，等.礦井煤巖動(dòng)力災(zāi)害聲發(fā)射監(jiān)測(cè)適用條件初探[J].煤炭學(xué)報(bào)，2011，36（2）：278-282.

[6]李庶林.試論微震監(jiān)測(cè)技術(shù)在地下工程中的應(yīng)用[J].地下空間與工程學(xué)報(bào)，2009，5（1）：122-128.

[7]姜忻良，崔奕，李圓，等.天津地鐵盾構(gòu)施工地層變形實(shí)測(cè)及動(dòng)態(tài)模擬[J].巖土力學(xué)，2005，26（10）：1 612-616.

[8]Peck R B.Deep excavations and tunneling in soft ground[C]// Proceedingsofthe7thInternational ConferenceonSoil Mechanics and Foundation Engineering.Mexico City：[s.n.]，1969：225-290.

[9]Nomikos P P，Sofianos A I，Tsoutrelis C E.Structural response of verticallymulti-jointed roof rock beams[J].International Journal of Rock Mechanics and Mining Sciences，2002，39（1）：79-94.

[10]王金安，李大鐘，尚新春.采空區(qū)堅(jiān)硬頂板流變破斷力學(xué)分析[J].北京科技大學(xué)學(xué)報(bào)，2011，33（2）：142-148．

[11]趙奎，王更峰，王曉軍，等.巖石聲發(fā)射Kaiser點(diǎn)信號(hào)頻帶能量分布和分形特征研究[J].巖土力學(xué)，2008，29（11）：3082-3088.

[12]趙奎，閆道全，王曉軍，等.聲發(fā)射測(cè)量地應(yīng)力綜合分析方法與實(shí)驗(yàn)驗(yàn)證[J].巖土工程學(xué)報(bào)，2012，34（8）：1403-1411.

[13]許福樂(lè)，王恩元，宋大釗，等．煤巖破壞聲發(fā)射強(qiáng)度長(zhǎng)程相關(guān)性和多重分形分布研究[J]．巖土力學(xué)，2011，32（7）：2111-2116

[14]王寧，韓志型，王月明，等．評(píng)價(jià)巖體穩(wěn)定性的聲發(fā)射相對(duì)強(qiáng)弱指標(biāo)[J]．巖土工程學(xué)報(bào)，2005（2）：190-192.

[15]王曉軍，付振濤，程秋亭，等．合理布筋方式提高頂板充填體穩(wěn)定性試驗(yàn)研究[J]．有色金屬科學(xué)與工程，2011，2（1）：73-76.

[16]胡京濤，趙奎，胡慧明，等．尾砂膠結(jié)充填體聲發(fā)射特征的分形分析[J]．有色金屬科學(xué)與工程，2011，2（2）：78-82.

Stope roofing stability disp lacem ent and its identification by acoustic em ission features

ZHAO Kui1,HU Yuan1,WANG Xiaojun1，2,FENG Xiao1,ZHUO Yulong1
（1.School of Resources and Environmental Engineering,JiangxiUniversity of Science and Technology,Ganzhou 341000,China；2.Post-doctoral Scientific Research Workstation，Western Mining Co.Ltd.，Xining 810006，China）

The roof deformation of underground projects is a complex nonlinear process.The singlemonitoring means is insufficient to determine its roofing stability accurately.Standard deviation of displacementmonitoring sample for roofing deformation process is deduced based onmathematical statistics.Three deformation stages for stope roofing,non-deformation,stable deformation and accelerated deformation are characterized according to the fore-mentioned results.Acoustic emission characteristics are analyzed during rock deformation and failure according to the combined analysis results of displacementand acoustic emission.The stability criterion of roof is divided into four kinds of modes.Combined with the engineering example of a mine,the roof failure is accurately predicted during mining process.The result provides a new data processing method for stability analysis of roof of underground projects by improving the accuracy of roof stability criterion.

standard deviation;accelerated deformation;acoustic emission;roofing stability

TD326

A

10.13264/j.cnki.ysjskx.2015.03.014

2014-12-30

國(guó)家自然科學(xué)基金資助項(xiàng)目（51304083）；江西省科技支撐計(jì)劃項(xiàng)目（20141BBE50005）

趙奎（1969-），男，博士，教授，主要從事巖石力學(xué)測(cè)試與分析技術(shù)、巖體穩(wěn)定性方面研究，E-mail：yglmf_zk@163.com.