程鐵棟,王運(yùn)來,張志釗,易其文,尹寶勇,袁海平

(1.江西理工大學(xué) 電氣工程與自動(dòng)化學(xué)院，江西 贛州 341000；2.江西理工大學(xué) 理學(xué)院，江西 贛州 341000；3.合肥工業(yè)大學(xué) 土木與水利工程學(xué)院，安徽 合肥 230009)

巖石聲發(fā)射(acoustic emission，AE)是指巖石在受到外部載荷或內(nèi)力作用下，其內(nèi)部含有的大量微孔洞、裂隙發(fā)生了聚合、擴(kuò)展及相互貫穿，儲(chǔ)存的能量以彈性波形式釋放的過程[1]。聲發(fā)射技術(shù)作為一種動(dòng)態(tài)檢測(cè)方法，被廣泛應(yīng)用于巖石力學(xué)工程領(lǐng)域的災(zāi)害監(jiān)測(cè)、預(yù)警方面[2]。通過對(duì)聲發(fā)射傳感器采集的巖石聲發(fā)射信號(hào)進(jìn)行分析，可以獲取巖石破壞的時(shí)間、地點(diǎn)、機(jī)制等信息。由于背景噪聲、機(jī)械噪聲等因素的干擾，傳感器所采集的聲發(fā)射信號(hào)波形極為復(fù)雜，導(dǎo)致聲發(fā)射信號(hào)的識(shí)別異常困難。因此，巖石聲發(fā)射信號(hào)的準(zhǔn)確識(shí)別對(duì)巖體失穩(wěn)災(zāi)變的預(yù)警研究具有重要意義。

隨著聲發(fā)射技術(shù)在各個(gè)領(lǐng)域的廣泛應(yīng)用，對(duì)信號(hào)去噪能力的要求逐漸提高。劉東瀛等[3]提出一種基于相關(guān)系數(shù)原理的經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition，EMD)去噪方法；戴聰聰?shù)萚4]提出一種將總體經(jīng)驗(yàn)?zāi)B(tài)分解法(ensemble empirical mode decomposition，EEMD)與云相似度理論相結(jié)合的去噪方法；姜長(zhǎng)泓等[5]通過對(duì)小波閾值法的改進(jìn)，提出基于平移不變小波的聲發(fā)射信號(hào)去噪方法。上述去噪方法雖然取得了不錯(cuò)的效果，但也存在不足之處。例如，EMD在分解信號(hào)的過程中易出現(xiàn)模態(tài)混疊及端點(diǎn)效應(yīng)等問題[6]，EEMD在分解信號(hào)的過程中會(huì)存在殘余的白噪聲[7]，小波閾值去噪在分解層數(shù)的確定和閾值的選擇等方面存在困難[8]。

特征提取是聲發(fā)射信號(hào)識(shí)別的關(guān)鍵，常用的聲發(fā)射信號(hào)特征提取方法主要包括基于參數(shù)的分析方法和基于波形的分析方法[9]。參數(shù)分析法通過對(duì)信號(hào)的統(tǒng)計(jì)特征參數(shù)進(jìn)行分析，是目前聲發(fā)射信號(hào)領(lǐng)域應(yīng)用最廣泛的方法。姚旭龍等[10]利用能量貢獻(xiàn)率構(gòu)建了巖石破裂關(guān)鍵事件聲發(fā)射信號(hào)的優(yōu)選方法；司莉等[11]通過對(duì)聲發(fā)射信號(hào)的遞歸定量分析，提出基于遞歸定量分析與支持向量機(jī)相結(jié)合的聲發(fā)射信號(hào)識(shí)別方法。這些方法為聲發(fā)射信號(hào)的識(shí)別提供了多個(gè)途徑并取得不錯(cuò)的效果，但也存在一些不足：能量貢獻(xiàn)率作為單一特征，提供的聲發(fā)射源信息有限；遞歸定量分析存在計(jì)算量大的問題。

變分模態(tài)分解(variational mode decomposition，VMD)是Dragomiretskiy等[12]提出的一種新的非遞歸式模態(tài)分解算法，該方法可有效避免EEMD在分解信號(hào)過程中產(chǎn)生殘余的白噪聲。然而，VMD對(duì)信號(hào)進(jìn)行分解時(shí)，需要提前設(shè)置本征模態(tài)分量(intrinsic mode function，IMF)的個(gè)數(shù)K，且K的取值影響算法的分解效果[13]。目前時(shí)頻域特征參數(shù)已被廣泛應(yīng)用于故障診斷領(lǐng)域，但多個(gè)特征融合參數(shù)在聲發(fā)射領(lǐng)域應(yīng)用較少。

鑒于此，本研究利用改進(jìn)的VMD算法對(duì)巖石聲發(fā)射信號(hào)進(jìn)行分解，結(jié)合排列熵篩選出最優(yōu)IMF分量并進(jìn)行重構(gòu)，然后提取重構(gòu)信號(hào)的多個(gè)特征參數(shù)，將多個(gè)特征參數(shù)融合并輸入到多元宇宙優(yōu)化算法(multi-verse optimizer，MVO)優(yōu)化的支持向量機(jī)(support vector machine，SVM)，對(duì)巖石聲發(fā)射信號(hào)進(jìn)行識(shí)別。

1 改進(jìn)的VMD算法

1.1 VMD算法

VMD算法的實(shí)質(zhì)是通過構(gòu)建并求解約束變分模型，得到多個(gè)具有特定稀疏性的IMF分量。VMD算法構(gòu)建的約束變分模型為：

(1)

式中：uk為VMD分解得到的各個(gè)IMF分量；ωk為各個(gè)IMF分量對(duì)應(yīng)的中心頻率；δ(t)為單位脈沖函數(shù)；j為虛數(shù)單位；f(t)為原始信號(hào)。

為了求解該約束變分模型，通過引入二次懲罰因子和拉格朗日乘子項(xiàng)將有約束變分問題轉(zhuǎn)換為如下無約束變分問題：

(2)

1.2 VMD算法的改進(jìn)

VMD算法需要提前設(shè)置IMF分量個(gè)數(shù)，且通過主觀經(jīng)驗(yàn)難以合理設(shè)置該參數(shù)。本研究利用奇異值最佳有效秩階次自動(dòng)獲得VMD分解的最優(yōu)K值。排列熵是一種檢測(cè)動(dòng)力學(xué)突變和時(shí)間序列隨機(jī)性的方法[14]，其大小反映時(shí)間序列的隨機(jī)程度，可以根據(jù)排列熵值評(píng)定IMF分量的含噪程度[15]。改進(jìn)VMD去噪算法的主要步驟如下。

1) 信號(hào)的奇異值分解

假設(shè)信號(hào)序列為X={x1,x2,…,xN}，將X構(gòu)造為p×q階Hankel矩陣：

(3)

對(duì)Hankel矩陣H進(jìn)行奇異值分解，得到：

(4)

式中：U∈Rp×p稱為左奇異矩陣，V∈Rq×q稱為右奇異矩陣，且均為正交矩陣；Σ=diag(σ1,σ2,…,σr),為對(duì)角矩陣。對(duì)角矩陣Σ中的元素σi(i=1,2,…,r)稱作H的奇異值，且奇異值為降序排列，即σ1≥σ2≥…≥σr≥0，r=min(p，q)為矩陣的H秩。

2) 求相鄰兩個(gè)奇異值之間的差值

若k階奇異值與k-1階奇異值之差最大，且與k+1階奇異值之差最小，即max(σk-σk-1), min(σk-σk+1)時(shí)，則將此時(shí)的k值確定為奇異值最佳有效秩階次。

3) 信號(hào)的VMD分解

將k設(shè)定為VMD分解的最優(yōu)K值并對(duì)信號(hào)進(jìn)行分解，得到K個(gè)IMF分量。

4) 計(jì)算K個(gè)IMF分量的排列熵值，對(duì)一個(gè)長(zhǎng)度為N的IMF分量進(jìn)行相空間重構(gòu)，得到矩陣

(5)

式中：m為嵌入維數(shù)，t為延遲時(shí)間，M=N-(m-1)t。

將矩陣Y中的每行按升序重新排列，得到向量中各元素位置的列索引構(gòu)成的一組符號(hào)序列:

S(l)={j1,j2,…,jm}，l=1,2,…,d。

(6)

式中:m維相空間映射不同的符號(hào)序列總共有m!種，d≤m!。

每一種符號(hào)序列出現(xiàn)的次數(shù)t與m!種符號(hào)序列的總次數(shù)M的比值作為該符號(hào)序列出現(xiàn)的概率，即：

(7)

進(jìn)而得到排列熵的計(jì)算式為：

(8)

將排列熵值進(jìn)行歸一化處理，即：

(9)

對(duì)每一個(gè)IMF分量都進(jìn)行排列熵的計(jì)算，可以得到K個(gè)IMF分量的排列熵值。

5) 信號(hào)重構(gòu)

(10)

2 多特征融合提取方法

2.1 梅爾頻率倒譜系數(shù)

梅爾頻率倒譜系數(shù)(Mel frequency cepstral coefficients，MFCC)是一種廣泛應(yīng)用于音頻場(chǎng)景和語音識(shí)別的特征參數(shù)，其計(jì)算簡(jiǎn)單、區(qū)分能力較為突出[16]。MFCC提取信號(hào)特征的過程如下。

1) 語音信號(hào)一般需要經(jīng)過預(yù)加重、分幀以及加窗處理，變成單幀的短時(shí)信號(hào)，而本研究的聲發(fā)射信號(hào)在進(jìn)行特征提取之前已經(jīng)是單幀的信號(hào)，可將每一幀信號(hào)視為一組信號(hào)進(jìn)行處理。

2) 對(duì)聲發(fā)射信號(hào)作快速傅里葉變換，得到頻域數(shù)據(jù)

(11)

3) 對(duì)頻域數(shù)據(jù)取模后再平方，得到信號(hào)的功率譜

E(h)=|X(h)|2。

(12)

4) 將信號(hào)的功率譜通過Q個(gè)Mel濾波器組進(jìn)行計(jì)算，得到：

(13)

式中，Hi(h)為第i個(gè)濾波器的傳遞函數(shù)，其計(jì)算式為：

(14)

式中，f(i)為中心頻率。

5) 將S(i)取對(duì)數(shù)并進(jìn)行離散余弦變換，可以得到：

(15)

式中，a為Mel頻率倒譜系數(shù)的個(gè)數(shù)，即特征向量的維數(shù)。a的取值范圍通常為12～16，增大該值會(huì)使得計(jì)算量相應(yīng)增加[17]，因而本研究中MFCC作為特征向量時(shí)的維數(shù)確定為12維。

2.2 短時(shí)能量

短時(shí)能量(short-time energy，SE)是對(duì)聲音信號(hào)在時(shí)域上的分析。由于聲發(fā)射信號(hào)是一種非線性、非平穩(wěn)信號(hào)，能量會(huì)隨時(shí)間變化，因此引入短時(shí)能量作為聲發(fā)射信號(hào)在時(shí)域方面的特征。短時(shí)能量的計(jì)算式為：

(16)

式中：xn(m)為輸入的第n幀信號(hào)；N為一幀信號(hào)的長(zhǎng)度；SEn為第n幀信號(hào)的能量，即短時(shí)能量。

2.3 譜質(zhì)心

譜質(zhì)心(spectral centroid，SC)是描述音色屬性的重要物理參數(shù)之一，是頻率成分的重心[18]。對(duì)信號(hào)在一定頻率范圍內(nèi)的能量進(jìn)行加權(quán)平均，便可得到該信號(hào)的譜質(zhì)心，計(jì)算式為：

(17)

式中：fi為信號(hào)的頻率，E(fi)為信號(hào)的能量譜，SC為信號(hào)做離散傅里葉變換的譜質(zhì)心。

2.4 特征融合

將上述3個(gè)特征向量以列向量的形式進(jìn)行組合，最終得到一個(gè)融合的特征向量

F=[F1,F2, …,F13,F14]。

(18)

式中：F1,F2, …,F12為MFCC，表征信號(hào)的頻域特征；F13為SE，表征信號(hào)的時(shí)域特征，與MFCC在時(shí)頻域起互補(bǔ)作用；F14為SC，反映信號(hào)的頻率分布和能量分布情況。

3 建立MVO-SVM模型

SVM是Vapnik[19]在統(tǒng)計(jì)學(xué)習(xí)理論和結(jié)構(gòu)風(fēng)險(xiǎn)最小化的基礎(chǔ)上提出的一種通用學(xué)習(xí)方法，該方法適合解決小樣本和非線性問題。MVO是Mirjalili等[20]提出的一種基于物理機(jī)理的啟發(fā)式優(yōu)化算法，具有待調(diào)節(jié)的參數(shù)少、搜索效率高及優(yōu)化能力強(qiáng)等優(yōu)點(diǎn)，已經(jīng)被成功應(yīng)用于許多工程問題中。SVM在應(yīng)用過程中，需要提前設(shè)置合適的懲罰參數(shù)γ。本研究利用MVO對(duì)SVM的參數(shù)進(jìn)行尋優(yōu)，可以避免因人工選取參數(shù)不合適導(dǎo)致模型分類精度降低問題。利用訓(xùn)練集樣本作為SVM的輸入數(shù)據(jù)，并引入MVO算法優(yōu)化SVM，從而建立MVO-SVM模型。如圖1所示，MVO-SVM模型的建立過程包括5個(gè)步驟。

圖1 MVO-SVM模型的建立過程

2) 定義SVM對(duì)訓(xùn)練數(shù)據(jù)的分類誤差作為MVO優(yōu)化SVM的適應(yīng)度函數(shù)，當(dāng)適應(yīng)度函數(shù)達(dá)到最小值時(shí)，對(duì)應(yīng)的宇宙分量作為當(dāng)前最優(yōu)分量。

3) 依據(jù)式(20)更新宇宙。

(20)

4) 當(dāng)算法運(yùn)行達(dá)到最大迭代次數(shù)時(shí)，停止更新宇宙，此時(shí)輸出的全局最優(yōu)宇宙即為懲罰參數(shù)C和核函數(shù)參數(shù)γ的最優(yōu)值。

5) 以懲罰參數(shù)C和核函數(shù)參數(shù)γ的最優(yōu)值重構(gòu)SVM模型。

4 實(shí)驗(yàn)分析

4.1 實(shí)驗(yàn)方案

實(shí)驗(yàn)室模擬巖石聲發(fā)射信號(hào)，聲發(fā)射信號(hào)采集設(shè)備的工作原理如圖2所示。

圖2 聲發(fā)射信號(hào)采集設(shè)備的工作原理

聲發(fā)射信號(hào)的采集使用江西理工大學(xué)自主研發(fā)的PIE-3-L型傳感器如圖3所示，傳感器的采樣頻率為145 kHz，聲發(fā)射信號(hào)的采樣點(diǎn)數(shù)為4 000。實(shí)驗(yàn)時(shí)，在傳感器與巖石的接觸面之間涂抹少量黃油，然后將傳感器放置于巖石的固定位置。

圖3 聲發(fā)射信號(hào)采集裝置

在距離聲發(fā)射傳感器20 cm處進(jìn)行4種不同信號(hào)的采集測(cè)試：①爆破信號(hào)。將鞭炮置于巖石底下，用鞭炮爆炸產(chǎn)生的信號(hào)模擬礦山現(xiàn)場(chǎng)巖石爆破產(chǎn)生的信號(hào)；②斷鉛信號(hào)。按國(guó)際無損檢測(cè)界規(guī)定，選擇0.5 mmHB鉛筆芯進(jìn)行斷鉛實(shí)驗(yàn)，以模擬巖石破斷聲發(fā)射信號(hào)；③電鉆信號(hào)。利用電鉆對(duì)巖石進(jìn)行鉆孔，模擬礦山現(xiàn)場(chǎng)的工況聲源；④敲擊信號(hào)。用金屬棒保持一定的力度敲擊巖石，模擬礦山現(xiàn)場(chǎng)的另一工況聲源。實(shí)驗(yàn)分別采集不同類型的信號(hào)各105組，總計(jì)采集到數(shù)據(jù)420組。

4.2 改進(jìn)的VMD去噪

從采集的聲發(fā)射信號(hào)中取一幀長(zhǎng)度為4 000點(diǎn)的信號(hào)，通過計(jì)算信號(hào)的奇異值最佳有效秩階次，得到VMD分解預(yù)設(shè)的分量個(gè)數(shù)K=3，再用VMD算法對(duì)信號(hào)進(jìn)行分解得到IMF1、IMF2、IMF3，進(jìn)而計(jì)算3個(gè)IMF分量的排列熵值分別為0.442 1、0.504 1、0.845 5，最后確定IMF1和IMF2為最優(yōu)分量并重構(gòu)，重構(gòu)后的信號(hào)如圖4所示。由圖4可見，改進(jìn)的VMD去噪后聲發(fā)射信號(hào)波形較為光滑，在剔除大部分噪聲的同時(shí)保留了聲發(fā)射信號(hào)的主要特征。對(duì)采集的多組聲發(fā)射信號(hào)分別運(yùn)用該方法去噪，均能達(dá)到較好的效果。

圖4 改進(jìn)的VMD去噪結(jié)果

原始信號(hào)分別經(jīng)EEMD、小波閾值和改進(jìn)的VWD去噪處理后的頻譜如圖5所示。圖5(a)為原始聲發(fā)射信號(hào)的頻譜，信號(hào)包含了部分高頻噪聲，頻譜分布范圍寬廣。圖5(b)為EEMD去噪后的信號(hào)頻譜，可以看出EEMD去噪效果一般，無法抑制大部分噪聲。圖5(c)為小波閾值去噪后的信號(hào)頻譜圖，表明小波閾值去噪能抑制高頻噪聲，對(duì)低頻噪聲也能起到較好的抑制作用，但該方法需要人為設(shè)定閾值，在信號(hào)重構(gòu)時(shí)會(huì)丟失部分有用信息。圖5(d)為改進(jìn)的VMD去噪頻譜圖，可以看出，該方法對(duì)低頻和高頻噪聲均能達(dá)到良好的抑制作用，且去噪后的頻譜集中分布于某一頻率，與EEMD去噪和小波閾值去噪算法相比，改進(jìn)的VMD去噪算法效果更好。

圖5 原始聲發(fā)射信號(hào)及各方法去噪后的頻譜

為定量評(píng)價(jià)改進(jìn)的VMD去噪算法的性能，用信噪比(signal-to-noise ratio，SNR)對(duì)3種算法的去噪效果進(jìn)行比較，結(jié)果見表1。

表1 3種去噪算法參數(shù)對(duì)比

由表1可見，相比于EEMD、小波閾值去噪，改進(jìn)的VMD去噪算法的SNR值分別提高了5.23、0.84 dB，說明改進(jìn)的VMD降噪方法相對(duì)較好。

4.3 融合特征向量提取

利用融合的特征向量對(duì)采集的爆破、斷鉛、電鉆和敲擊信號(hào)進(jìn)行特征提取，再對(duì)所有信號(hào)樣本的特征值進(jìn)行歸一化處理。為展示不同類別信號(hào)樣本的差異，表2列出了每一類信號(hào)中4個(gè)樣本的14維特征值。其中，F(xiàn)1～F12為MFCC的12維特征參數(shù)，F(xiàn)13為SE的一維特征參數(shù)，F(xiàn)14為SC的一維特征參數(shù)。由表2可見，不同類別信號(hào)的特征參數(shù)存在較大差異。

表2 4類信號(hào)的特征參數(shù)對(duì)比

4.4 聲發(fā)射信號(hào)識(shí)別

1) 信號(hào)預(yù)處理。利用改進(jìn)的VMD方法對(duì)采集的4類信號(hào)共計(jì)420組數(shù)據(jù)進(jìn)行去噪處理。

2) 構(gòu)建融合特征向量，將多個(gè)特征參數(shù)融合作為分類識(shí)別的特征向量。劃分訓(xùn)練數(shù)據(jù)和測(cè)試數(shù)據(jù)，從已有的4類信號(hào)共420組數(shù)據(jù)中，按照訓(xùn)練數(shù)據(jù)與測(cè)試數(shù)據(jù)為4∶1的劃分方式對(duì)信號(hào)進(jìn)行隨機(jī)分層采樣，以保證訓(xùn)練數(shù)據(jù)和測(cè)試數(shù)據(jù)中的類別數(shù)目均衡。設(shè)定爆破信號(hào)的標(biāo)識(shí)類別為1，斷鉛信號(hào)的標(biāo)識(shí)類別為2，電鉆信號(hào)的標(biāo)識(shí)類別為3，敲擊信號(hào)的標(biāo)識(shí)類別為4。

3) 將訓(xùn)練集和測(cè)試集分別作歸一化處理。目的是將所有的特征值大小限定在[0，1]，防止某個(gè)特征值過大或者過小。

4) 建立MVO-SVM模型。把MVO算法尋優(yōu)得到的模型參數(shù)作為懲罰參數(shù)、核函數(shù)參數(shù)。采用RBF核函數(shù)訓(xùn)練SVM模型，利用訓(xùn)練好的SVM模型對(duì)測(cè)試集分類，結(jié)果如圖6所示。

圖6 MVO-SVM分類識(shí)別結(jié)果

圖6中，測(cè)試集樣本數(shù)量為84個(gè)，類別標(biāo)簽為1、2、3、4，樣本實(shí)際標(biāo)簽與樣本預(yù)測(cè)標(biāo)簽重合表示該樣本被正確識(shí)別。由圖6可見，類別1樣本中有1個(gè)樣本被誤識(shí)別為類別3、類別2樣本中有1個(gè)樣本被誤識(shí)別為類別4，類別3樣本中有1個(gè)樣本被誤識(shí)別為類別2，類別4樣本中有1個(gè)樣本被誤識(shí)別為類別1。綜上，MVO-SVM對(duì)爆破、斷鉛、電鉆及敲擊信號(hào)的識(shí)別準(zhǔn)確率為95.24%。

為了測(cè)試去噪方法的有效性與融合特征向量的優(yōu)越性，采用改進(jìn)的VMD去噪、EEMD去噪、小波閾值去噪等方法，對(duì)比不同特征提取方法在MVO-SVM模型上的識(shí)別準(zhǔn)確率，如表3所示。

表3 不同去噪方法下各種特征提取方法的識(shí)別準(zhǔn)確率

由表3可知：①在不同的去噪方法下，采用融合特征向量提取聲發(fā)射信號(hào)的特征，可以獲得較高的識(shí)別準(zhǔn)確率，而單一的特征向量則無法獲得較好的效果，且僅SE和SC作為特征向量時(shí)，聲發(fā)射信號(hào)的識(shí)別準(zhǔn)確率均較低；②在不同的特征提取方法下，利用改進(jìn)的VMD去噪后識(shí)別效果整體上優(yōu)于EEMD和小波閾值去噪方法；③在EEMD和小波閾值去噪方法上，單一特征向量的識(shí)別準(zhǔn)確率相同，而融合特征向量的識(shí)別準(zhǔn)確率顯著提高，且在兩種去噪方法上的識(shí)別效果存在差異。因此，基于改進(jìn)的VMD與融合特征向量構(gòu)建的聲發(fā)射信號(hào)識(shí)別方法準(zhǔn)確率高，對(duì)4種不同類別信號(hào)均能進(jìn)行有效識(shí)別。

5 結(jié)論

1) 提出一種改進(jìn)的VMD算法，利用奇異值最佳有效秩階次確定VMD分解的IMF分量個(gè)數(shù)K，并利用排列熵篩選出最優(yōu)IMF分量并重構(gòu)。將改進(jìn)的VMD算法應(yīng)用于聲發(fā)射信號(hào)去噪，提高了信號(hào)的信噪比。

2) 利用MFCC、SE和SC等3種特征構(gòu)建融合特征向量F，以此計(jì)算信號(hào)的特征值，使各類信號(hào)樣本的特征值差異更加明顯，為MVO-SVM模型識(shí)別聲發(fā)射信號(hào)提供了新的方法。

3) 構(gòu)建融合特征向量的特征提取方法，在EEMD去噪和小波閾值去噪的基礎(chǔ)上均可以取得較好的識(shí)別效果，而改進(jìn)VMD去噪的融合特征向量的識(shí)別準(zhǔn)確率最高，可以達(dá)到95.24%。表明基于改進(jìn)的VMD去噪與融合特征向量構(gòu)建的巖石聲發(fā)射信號(hào)識(shí)別方法是可行的，能有效識(shí)別出聲發(fā)射信號(hào)。

本研究提出的基于改進(jìn)VMD去噪與多特征融合的聲發(fā)射信號(hào)識(shí)別方法對(duì)室內(nèi)模擬巖石聲發(fā)射信號(hào)可以達(dá)到較高的識(shí)別準(zhǔn)確率，但對(duì)于實(shí)際礦山巖石聲發(fā)射信號(hào)的識(shí)別還需進(jìn)一步驗(yàn)證，后續(xù)將針對(duì)現(xiàn)場(chǎng)巖石聲發(fā)射信號(hào)進(jìn)行研究，建立泛化性能更好的分類模型，以期得到較好的識(shí)別效果，為實(shí)際工程應(yīng)用提供參考。