基于EEMD-SVM 的混凝土大壩活動(dòng)裂紋電磁輻射信號(hào)識(shí)別方法

摘 要： 為解決混凝土大壩活動(dòng)裂紋電磁輻射信號(hào)與電磁干擾信號(hào)混淆、識(shí)別準(zhǔn)確率低的問(wèn)題，該文提出一種基于集成經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）-支持向量機(jī)（SVM）的混凝土大壩活動(dòng)裂紋電磁輻射信號(hào)識(shí)別方法。該方法利用EEMD 將混凝土大壩周?chē)h(huán)境電磁輻射信號(hào)分解為多個(gè)本征模態(tài)函數(shù)，構(gòu)建漢克爾矩陣和奇異值分解方法，獲取不同本征模態(tài)函數(shù)的平均值、方差和峭度值等信號(hào)特征量，結(jié)合SVM 識(shí)別混凝土大壩活動(dòng)裂紋電磁輻射信號(hào)和周?chē)h(huán)境干擾噪聲。對(duì)400 組混凝土活動(dòng)裂紋電磁輻射信號(hào)與干擾噪聲進(jìn)行比較，所提出的方法識(shí)別成功率高，最終識(shí)別準(zhǔn)確率大于95%，高于其他傳統(tǒng)信號(hào)處理算法的識(shí)別率，可為混凝土大壩活動(dòng)裂紋的識(shí)別提供參考依據(jù)。

關(guān)鍵詞： 混凝土大壩; 活動(dòng)裂紋; 電磁輻射; EEMD-SVM; 信號(hào)識(shí)別

中圖分類(lèi)號(hào)： TB9 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)： 1674–5124（2025）03–0016–06

0 引 言

混凝土大壩作為重要的水利工程結(jié)構(gòu)，承擔(dān)著水資源調(diào)控和防洪任務(wù)。然而，長(zhǎng)期的自然和人為因素的影響，使混凝土大壩存在越來(lái)越多的活動(dòng)裂紋?；顒?dòng)裂紋一旦發(fā)生快速擴(kuò)展，將導(dǎo)致混凝土大壩內(nèi)部結(jié)構(gòu)破壞，易發(fā)生安全事故[1-2]。因此，準(zhǔn)確、及時(shí)地監(jiān)測(cè)和識(shí)別混凝土大壩中的活動(dòng)裂紋是至關(guān)重要的工作。

研究表明，混凝土等脆性材料在開(kāi)裂過(guò)程中，裂紋表面存在電荷的產(chǎn)生和移動(dòng)，從而產(chǎn)生電磁輻射現(xiàn)象[3-4]。然而，混凝土大壩周?chē)姶怒h(huán)境復(fù)雜，活動(dòng)裂紋產(chǎn)生的電磁輻射信號(hào)弱，電磁干擾噪聲時(shí)域波形與活動(dòng)裂紋產(chǎn)生的電磁信號(hào)類(lèi)似。若依靠人工方式識(shí)別活動(dòng)裂紋信號(hào)，則經(jīng)常出現(xiàn)漏檢和誤檢的情況。

由于混凝土大壩活動(dòng)裂紋產(chǎn)生的電磁輻射信號(hào)為非平穩(wěn)信號(hào)，傅里葉分析方法不適用于其電磁輻射信號(hào)的識(shí)別。常用的信號(hào)特征提取方法主要有短時(shí)傅里葉變換[5]（STFT）、小波變換[6]（WT）、希爾伯特黃變換[7]（HHT）、經(jīng)驗(yàn)?zāi)B(tài)分解[8]（EMD）和經(jīng)驗(yàn)小波變換[9]（EWT）。Wei 等[10] 研究了巖石破裂電磁輻射波形的傅里葉譜。結(jié)果表明，在不同應(yīng)力水平下，主頻率隨裂紋的演化而減小。Yin 等[8] 利用希爾伯特黃變換對(duì)煤塊破裂的不同階段進(jìn)行分析，破裂信號(hào)的頻率隨階段的變化而降低。Lou 等[11]利用EWT 對(duì)巖石破裂產(chǎn)生的電磁輻射信號(hào)進(jìn)行分析，與STFT 和HHT 相比，EWT 在時(shí)頻分辨率和克服模態(tài)混疊方面具有顯著優(yōu)勢(shì)。Li 等[12] 用連續(xù)小波變換方法得到不同干擾噪聲和危險(xiǎn)信號(hào)在時(shí)頻上存在顯著差異。然而現(xiàn)有方法主要應(yīng)用于煤層或巖石破裂的研究中，而缺乏混凝土材料開(kāi)裂時(shí)電磁輻射信號(hào)的研究。

因此，本文提出一種基于集成經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）-支持向量機(jī)（SVM）的混凝土大壩活動(dòng)裂紋電磁輻射信號(hào)識(shí)別方法。通過(guò)集成經(jīng)驗(yàn)?zāi)B(tài)分解方法將所需識(shí)別的電磁輻射信號(hào)分解為多個(gè)本征模態(tài)函數(shù)，對(duì)每個(gè)本征模態(tài)函數(shù)構(gòu)建漢克爾矩陣，進(jìn)行奇異值分解，提取其平均值、峭度值和方差等特征量，結(jié)合支持向量機(jī)分類(lèi)方法，構(gòu)建電磁輻射信號(hào)識(shí)別模型，最終實(shí)現(xiàn)混凝土活動(dòng)裂紋電磁輻射信號(hào)和電磁干擾信號(hào)的區(qū)分。

1 集成經(jīng)驗(yàn)?zāi)B(tài)分解

經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）作為一種自適應(yīng)信號(hào)分解方法，原理簡(jiǎn)單，構(gòu)造方便，廣泛應(yīng)用于非線(xiàn)性和非平穩(wěn)信號(hào)分解中。EMD 將非線(xiàn)性和非平穩(wěn)信號(hào)分解成多組本征模態(tài)函數(shù)，實(shí)現(xiàn)信號(hào)在不同頻率范圍內(nèi)的細(xì)節(jié)層次劃分。與小波分析和短時(shí)傅里葉變化等時(shí)頻分析方法相比，其對(duì)信號(hào)局部細(xì)節(jié)的提取能力不受窗函數(shù)和小波基的影響，更具普適性。

但是經(jīng)驗(yàn)?zāi)B(tài)分解方法存在噪聲和模態(tài)混疊問(wèn)題，多頻信號(hào)之間存在相互耦合，無(wú)法提取單一信號(hào)，因此在此基礎(chǔ)上黃鍔學(xué)者提出了一種集成經(jīng)驗(yàn)?zāi)B(tài)分解方法[13]，通過(guò)在原始信號(hào)的多次迭代過(guò)程中加入隨機(jī)噪聲，以減小模態(tài)之間的相互影響和混疊問(wèn)題。混凝土大壩活動(dòng)裂紋產(chǎn)生的電磁輻射信號(hào)為非平穩(wěn)信號(hào)，頻率成分復(fù)雜。因此，本文采用EEMD 方法對(duì)原始信號(hào)進(jìn)行處理，其主要步驟流程如圖1 所示。

S1：x（t） 為混凝土大壩周?chē)h(huán)境電磁輻射信號(hào)，w（t） 為白噪聲序列，取r（t）=x（t）+ w（t），i=0，k=1；設(shè)置終止閾值條件為r（t）lt;δ，取δ=0.2。