韓志軍，曹 威，孫曉云，程 恩，王明明，吳世星

（1.石家莊鐵道大學(xué)中鐵建安工程設(shè)計(jì)院，河北石家莊050043；2.石家莊鐵道大學(xué)電氣與電子工程學(xué)院，河北石家莊050043；）

非周期變尺度隨機(jī)共振在錨桿無損檢測(cè)中的應(yīng)用

韓志軍1，曹 威2，孫曉云2，程 恩2，王明明2，吳世星2

（1.石家莊鐵道大學(xué)中鐵建安工程設(shè)計(jì)院，河北石家莊050043；2.石家莊鐵道大學(xué)電氣與電子工程學(xué)院，河北石家莊050043；）

在強(qiáng)噪聲下檢測(cè)錨桿微弱信號(hào)一直是實(shí)際工程中的難題。隨機(jī)共振與以往去噪方法不同，其利用噪聲能量來增強(qiáng)被噪聲埋沒的微弱信號(hào)。經(jīng)典的隨機(jī)共振只適合于小參數(shù)信號(hào)（信號(hào)頻率、信號(hào)幅值、噪聲強(qiáng)度均遠(yuǎn)小于1）的檢測(cè)，而應(yīng)力波反射法采集到的錨桿信號(hào)屬于衰減振蕩的暫態(tài)信號(hào)，具有非周期性、大參數(shù)等特點(diǎn)，因此采用非周期變尺度隨機(jī)共振方法處理錨桿信號(hào)。針對(duì)非周期隨機(jī)共振，綜合峭度指標(biāo)與互相關(guān)系數(shù)指標(biāo)的各自優(yōu)缺點(diǎn)，采用加權(quán)峭度作為其度量指標(biāo)。變尺度隨機(jī)共振法對(duì)幅值進(jìn)行歸一化、對(duì)頻率進(jìn)行線性壓縮，使錨桿信號(hào)滿足經(jīng)典隨機(jī)共振的要求。通過處理非周期大參數(shù)的仿真信號(hào)和實(shí)際工程中的錨桿振動(dòng)信號(hào)，驗(yàn)證了非周期變尺度隨機(jī)共振的有效性。

錨桿；隨機(jī)共振；加權(quán)峭度；變尺度

錨桿作為巖土體加固、支護(hù)圍巖的重要桿件結(jié)構(gòu)，廣泛應(yīng)用于邊坡、隧道、礦山、壩體中［1-2］。在錨

桿損傷檢測(cè)中，由于惡劣的檢測(cè)環(huán)境、傳感器的數(shù)據(jù)漂移以及檢測(cè)儀器自身所帶來的干擾，導(dǎo)致采集到的數(shù)據(jù)常常淹沒在強(qiáng)噪聲下，給損傷識(shí)別帶來了極大的困難。在強(qiáng)噪聲背景下檢測(cè)錨桿微弱信號(hào)一直是錨桿無損檢測(cè)領(lǐng)域備受關(guān)注的熱點(diǎn)問題。常用的降噪方法有自適應(yīng)濾波法、小波變換法［3］、短時(shí)傅里葉變換［4］、高階譜分析、經(jīng)驗(yàn)?zāi)B(tài)分解［5］等。這些方法在降噪的同時(shí)也削弱和扭曲了有用信號(hào)，導(dǎo)致有用信號(hào)缺失，使得檢測(cè)結(jié)果不準(zhǔn)確。而隨機(jī)共振（stochastic resonance，SR）與以往的信號(hào)處理方法不同，隨機(jī)共振利用噪聲能量來增強(qiáng)微弱信號(hào)，使淹沒在強(qiáng)噪聲中的微弱信號(hào)達(dá)到共振，極大地提高了輸出信噪比［6］。

由于受絕熱近似理論限制，經(jīng)典的隨機(jī)共振只能分析周期性的小參數(shù)信號(hào)（信號(hào)頻率、信號(hào)幅值、噪聲強(qiáng)度均遠(yuǎn)小于1）［7］。應(yīng)力波反射法采集到的錨桿信號(hào)屬于衰減振蕩的暫態(tài)信號(hào)，具有非周期性、大參數(shù)等特點(diǎn)，并不適用于經(jīng)典的隨機(jī)共振?；诖?，本文提出了一種非周期變步長(zhǎng)檢測(cè)大參數(shù)的隨機(jī)共振方法，對(duì)于錨桿信號(hào)檢測(cè)提供了新的思路。

1 雙穩(wěn)態(tài)隨機(jī)共振模型

雙穩(wěn)態(tài)系統(tǒng)是研究隨機(jī)共振現(xiàn)象的基礎(chǔ)，其中該系統(tǒng)的典型模型為具有雙勢(shì)阱性質(zhì)的郎之萬方程［8］，具體公式見式（1）。

式中：U′表示一個(gè)四階的帶有兩個(gè)勢(shì)阱和一個(gè)勢(shì)壘的對(duì)稱雙穩(wěn)態(tài)勢(shì)函數(shù)；s（t）表示待檢測(cè)信號(hào)；z（t）表示均值為零、方差為D的高斯白噪聲。

雙穩(wěn)態(tài)勢(shì)函數(shù)的表達(dá)式見式（2）。則勢(shì)函數(shù)對(duì)應(yīng)的圖形見圖1。

式中：a表示系統(tǒng)參數(shù)，該參數(shù)是一個(gè)影響系統(tǒng)穩(wěn)態(tài)點(diǎn)位置的正數(shù)；b表示影響系統(tǒng)穩(wěn)態(tài)點(diǎn)位置的另一個(gè)參數(shù)，也是正數(shù)。

圖1 對(duì)稱雙穩(wěn)態(tài)勢(shì)函數(shù)示意圖

圖1中，兩個(gè)勢(shì)阱的位置為x=±xm，xm=勢(shì)壘的高ΔU=a2/4b。

2 非周期隨機(jī)共振特性

在早期隨機(jī)共振的研究中，大多處理的都是周期性的信號(hào)，這大大縮小了隨機(jī)共振的應(yīng)用范圍。Collins于1995年提出了非周期隨機(jī)共振（Aperiodic Stochastic Resonance，ASR）理論，指出當(dāng)輸入為非周期信號(hào)和噪聲組成時(shí)，非周期信號(hào)也會(huì)與噪聲產(chǎn)生共振，噪聲能量轉(zhuǎn)移到非周期信號(hào)，非周期信號(hào)得到放大與增強(qiáng)［9］。ASR理論的提出擴(kuò)寬了隨機(jī)共振的應(yīng)用范圍。

ASR與SR最大的區(qū)別在于評(píng)價(jià)指標(biāo)的不同。非周期信號(hào)頻譜與周期信號(hào)頻譜不一樣，非周期信號(hào)不具有確定頻率，頻譜不再僅分布于單個(gè)頻率點(diǎn)附近，因此ASR不能采用經(jīng)典SR的測(cè)量指標(biāo)—信噪比。檢測(cè)ASR常用的度量指標(biāo)是峭度或互相關(guān)系數(shù)。峭度K不受信號(hào)絕對(duì)水平影響，對(duì)信號(hào)中的沖擊成分十分敏感，因此可以作為評(píng)價(jià)檢測(cè)衰減振動(dòng)信號(hào)的指標(biāo)［10］。K定義為信號(hào)四階中心矩與二階中心距的平分的比值，見式（3）。

式中：x=｛x1，x2，…，xN｝表示離散信號(hào)序列；N表示信號(hào)長(zhǎng)度表示離散信號(hào)序列x的均值。

當(dāng)信號(hào)中存在多個(gè)沖擊時(shí)，由于參數(shù)選擇不當(dāng)，會(huì)導(dǎo)致某些沖擊不明顯或被破壞，此時(shí)K值反而變大。因此，選取峭度K作為ASR的檢測(cè)指標(biāo)存在不足。

互相關(guān)系數(shù)C表示的是兩個(gè)信號(hào)之間的相關(guān)程度。C定義見如式（4）。

式中：x（n）和y（n）表示兩個(gè)能量有限的確定性信號(hào)；T表示數(shù)據(jù)的長(zhǎng)度。C越大，表征兩個(gè)信號(hào)越相似，|C|≤1。

如果已知無噪聲的輸入信號(hào)，C可以作為檢測(cè)指標(biāo)。但在實(shí)際工程中，往往無法預(yù)知無噪聲的輸入信號(hào)，限制了C作為ASR檢測(cè)指標(biāo)。

加權(quán)峭度指標(biāo)整合了峭度和互相關(guān)系數(shù)的優(yōu)點(diǎn)，既保證了對(duì)沖擊成分的敏感性又兼顧了輸出信號(hào)與原始信號(hào)的相似性，定義見式（5）。

式中：sgn（）表示符號(hào)函數(shù)，保證輸出信號(hào)和原始信號(hào)的相位盡量一致；K表示峭度指標(biāo)；C表示互相關(guān)系數(shù)；r表示C的指數(shù)，取正實(shí)數(shù)，調(diào)整輸入信號(hào)與輸出的指標(biāo)。

本文將加權(quán)峭度Kw作為非周期信號(hào)的測(cè)量指標(biāo)，加權(quán)峭度Kw越大，SR效果越好。

3 變尺度隨機(jī)共振原理

Mc Namara等學(xué)者于1989年提出絕熱近似理論［7］，該理論指出系統(tǒng)若產(chǎn)生SR，輸入信號(hào)頻率、輸入信號(hào)幅值和噪聲均應(yīng)遠(yuǎn)遠(yuǎn)小于1，即w＜＜1，A＜＜1，D＜＜1，這對(duì)于錨桿損傷檢測(cè)信號(hào)受到了很大限制。

本文采用變尺度的方法來實(shí)現(xiàn)工程中大參數(shù)信號(hào)的SR。首先對(duì)實(shí)測(cè)信號(hào)的幅值進(jìn)行歸一化處理、對(duì)頻率按一定比例進(jìn)行壓縮（設(shè)壓縮尺度比為R），使實(shí)測(cè)信號(hào)幅值與頻率均滿足小參數(shù)隨機(jī)共振條件。其次將壓縮的小參數(shù)信號(hào)輸入到雙穩(wěn)態(tài)系統(tǒng)得到輸出響應(yīng)譜，分析響應(yīng)譜的頻譜特性。最后，按壓縮尺度比還原實(shí)測(cè)信號(hào)。該方法的實(shí)質(zhì)就是對(duì)大參數(shù)的實(shí)測(cè)信號(hào)進(jìn)行線性壓縮，以符合SR的要求，因此并沒有改變?cè)瓬y(cè)量信號(hào)的任何性質(zhì)。

變尺度壓縮的具體步驟如下所述：①對(duì)信號(hào)幅值線性歸一化，即x'=x/max（abs（x））；②確定一個(gè)頻率壓縮尺度比R，采樣頻率壓縮后變?yōu)閒sr= fs/R；③時(shí)間軸被重新定義為t'=Rt，步長(zhǎng)變?yōu)閔' =1/fsr；④根據(jù)四階龍格庫塔法數(shù)值求解雙穩(wěn)態(tài)系統(tǒng)的輸出響應(yīng)；⑤按壓縮尺度還原輸出信號(hào)。

移頻變尺度SR流程圖見圖2。

圖2 移頻變尺度隨機(jī)共振流程圖

4 仿真驗(yàn)證

為了驗(yàn)證本文的非周期變尺度SR算法，構(gòu)造一個(gè)非周期大參數(shù)衰減的振動(dòng)仿真信號(hào)，見式（6）。

式中：n=3；A1=22；A2=30；A3=28；f1=150 Hz；f2=310Hz；f3=440Hz；N（t）為高斯白噪聲，其噪聲強(qiáng)度D=90，采樣頻率fs=20k Hz，采樣點(diǎn)數(shù)N= 4096。加噪信號(hào)的時(shí)域圖和頻譜圖如圖3所示。從加噪信號(hào)的時(shí)域圖和頻譜圖上可以看出信號(hào)完全被噪聲淹沒，無法識(shí)別原始信號(hào)。

對(duì)信號(hào)進(jìn)行變尺度處理，首先對(duì)信號(hào)幅值線性歸一化。其次，根據(jù)圖3（b）頻譜圖，確定頻率壓縮尺度比R=10000，壓縮采樣頻率為fsr=2Hz，步長(zhǎng)h=0.5。歸一化后的信號(hào)s（t）見圖4。

圖3 加噪信號(hào)的時(shí)域圖和頻譜圖

圖4 變尺度信號(hào)的時(shí)域圖和頻譜圖

對(duì)歸一化信號(hào)進(jìn)行SR處理。根據(jù)仿真信號(hào)特點(diǎn)，當(dāng)SR模型的參數(shù)a=0.1、b=4.8時(shí)，加權(quán)峭度Kw為最大值2.6059。經(jīng)過隨機(jī)共振后如圖5所

示，頻譜上可看到譜峰fs1=0.01465 Hz，fs2= 0.03076Hz，fs3=0.04395 Hz。按照頻率壓縮比例換算后頻率f1′=fs1R=146.5 H，f2′=fs2R= 307.6 Hz f3′=fs3R=439.5 Hz，f1′與f1、f2′與f2、f3′與f3基本相同，因此SR能夠很好地檢測(cè)出衰減的振蕩信號(hào)。

5 實(shí)驗(yàn)驗(yàn)證

5.1 錨桿測(cè)試實(shí)驗(yàn)

實(shí)驗(yàn)對(duì)象采用工程中長(zhǎng)度為5m、直徑為22mm的螺紋鋼錨桿。實(shí)驗(yàn)裝置采用chirp信號(hào)發(fā)射探頭、IEPE加速度傳感器、東華DH5923動(dòng)態(tài)信號(hào)測(cè)試儀器，實(shí)驗(yàn)采集過程見圖6。其中采樣頻率fs=100k Hz。

圖5 輸出信號(hào)的時(shí)域圖和頻譜圖

圖6 實(shí)驗(yàn)裝置圖

東華DH5923動(dòng)態(tài)信號(hào)測(cè)試儀器包含動(dòng)態(tài)信號(hào)測(cè)試系統(tǒng)所需的信號(hào)調(diào)理器（應(yīng)變、振動(dòng)等調(diào)理器）、直流電壓放大器、低通濾波器、抗混濾波器、16位A/D轉(zhuǎn)換器、以及采樣控制和計(jì)算機(jī)通訊的全部硬件，而且提供了充分考慮用戶方便操作本系統(tǒng)所需的控制軟件及分析軟件，是以計(jì)算機(jī)為基礎(chǔ)、智能化的動(dòng)態(tài)信號(hào)測(cè)試分析系統(tǒng)。該測(cè)試分析系統(tǒng)除了滿足多通道、高精度、高速動(dòng)態(tài)信號(hào)的測(cè)量需求，還以DMA方式傳送數(shù)據(jù)，保證了數(shù)據(jù)傳輸?shù)母咚?、穩(wěn)定、不漏碼。除此之外，該系統(tǒng)的便攜式特點(diǎn)使得它不僅可以應(yīng)用在實(shí)驗(yàn)室中，而且可以很方便地應(yīng)用在野外現(xiàn)場(chǎng)。

采集到強(qiáng)噪聲下的錨桿振動(dòng)信號(hào)，如圖7（a）所示，經(jīng)過FFT變換得頻譜圖，如圖7（b）所示。由圖7可知，錨桿振動(dòng)信號(hào)淹沒在強(qiáng)噪聲下。

由于錨桿振動(dòng)信號(hào)屬于非周期大參數(shù)信號(hào)，需要采用變尺度法對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行預(yù)處理。首先對(duì)信號(hào)幅值線性歸一化；其次，根據(jù)圖7（b）頻譜圖，確定頻率壓縮比R=50000。則壓縮采樣頻率為fsr= 2 Hz，步長(zhǎng)h=0.5。經(jīng)處理后的信號(hào)時(shí)域圖和頻譜圖如圖8所示。

圖7 加噪信號(hào)的時(shí)域圖和頻譜圖

圖8 變尺度信號(hào)的時(shí)域圖和頻譜圖

根據(jù)信號(hào)特點(diǎn)，當(dāng)SR模型的參數(shù)a=0.01、b= 0.2時(shí)，加權(quán)峭度Kw為最大值2.6011。SR輸出時(shí)域圖和頻譜圖見圖9。

從圖9可知，錨桿信號(hào)主頻被凸顯，壓縮主頻為f=0.0625 Hz，還原后的原始主頻為ff=f×R= 3125Hz。經(jīng)過SR之后噪聲能量從高頻轉(zhuǎn)移到低頻

區(qū)域，高頻噪聲成分被削減，低頻原始信號(hào)被放大，可以檢測(cè)出淹沒在噪聲中的原始信號(hào)。

圖9 隨機(jī)共振信號(hào)的時(shí)域圖和頻譜圖

6 結(jié) 論

針對(duì)錨桿信號(hào)具有非周期性、大參數(shù)的特點(diǎn)，本文以經(jīng)典小參數(shù)隨機(jī)共振為基礎(chǔ)，采用非周期變尺度隨機(jī)共振方法來處理錨桿信號(hào)。對(duì)于非周期信號(hào)，加權(quán)峭度指標(biāo)兼顧了峭度與互相關(guān)系數(shù)的優(yōu)缺點(diǎn)，因此，采用加權(quán)峭度作為度量指標(biāo)。對(duì)于大參數(shù)信號(hào)，運(yùn)用變尺度的方法對(duì)幅值和頻率進(jìn)行了線性壓縮，既保留了原始信號(hào)的特征又滿足小參數(shù)隨機(jī)共振的要求。通過仿真與實(shí)驗(yàn)驗(yàn)證了算法的有效性與優(yōu)越性，為工程中檢測(cè)被噪聲埋沒的錨桿信號(hào)提供了良好的解決途徑。

［1］ 孫曉云，王志遠(yuǎn)，程久龍，等.基于譜峭度的錨桿參數(shù)檢測(cè)方法［J］.中國礦業(yè)，2015，24（12）：147-151.

［2］ Xiaoyun Sun，Jiulong Cheng，Donghui Liu，et al.Application of Information Processing Technology on Nondestructive Detecting of Rock Bolts，Int.J.Modelling，Identification and Control，2012，16（3）：241-245

［3］ 孫曉云，程久龍.現(xiàn)代信息處理技術(shù)在錨桿無損檢測(cè)中的應(yīng)用探討［J］.中國礦業(yè)，2008，17（11）：87-89，92.

［4］ 徐健，張語勍，李彥斌，等.短時(shí)傅里葉變換和S變換用于檢測(cè)電壓暫降的對(duì)比研究［J］.電力系統(tǒng)保護(hù)與控制，2014，42（24）：44-48.

［5］ 張麗華，劉海波；郭金鑫.基于EMD和LS-SVM的復(fù)合地基沉降［J］.中國礦業(yè)，2014，23（11）：141-144

［6］ BRUCE M，KURT W.Theory of stochastic resonance［J］. Physical Review A，1989，39（9）：4854-4869.

［7］ Mc Namara B，Wiesenfeld K.Theory of stochastic resonance［J］.Rhys.Rev.A，1989，39（9）：4854-4869.

［8］ 韓東穎，丁雪娟，時(shí)培明.基于自適應(yīng)變尺度頻移帶通隨機(jī)共振降噪的EMD多頻微弱信號(hào)檢測(cè)［J］.機(jī)械工程學(xué)報(bào)，2013，49（8）：10-18.

［9］ J.J.Collins，Carson C.Chow，Thomas T.Imhoff.Aperiodic stochastic resonance in excitable systems［J］.Physical review.E，1995，52（4）：51-54.

［10］ Antoni J，RandallR B.The spectral kurtosis：application to the vibratory surveillance and diagnostics of rotating machines［J］.Mechanical Systems and Signal Processing，2006，20（2）：308-331.

內(nèi)蒙古石墨產(chǎn)業(yè)發(fā)展聯(lián)盟正式成立

內(nèi)蒙古石墨產(chǎn)業(yè)發(fā)展聯(lián)盟近日在呼和浩特內(nèi)蒙古礦業(yè)（集團(tuán)）有限責(zé)任公司正式成立。

該聯(lián)盟由內(nèi)蒙古自治區(qū)人民政府支持，經(jīng)內(nèi)蒙古自治區(qū)經(jīng)信委批準(zhǔn)，內(nèi)蒙古礦業(yè)（集團(tuán)）有限責(zé)任公司組織牽頭，聯(lián)合廈門大學(xué)、蘇州中材非金屬礦工業(yè)設(shè)計(jì)研究院有限公司、內(nèi)蒙古金彩實(shí)業(yè)集團(tuán)、內(nèi)蒙古元亨石墨礦業(yè)有限公司等8家單位共同發(fā)起籌建。經(jīng)大會(huì)選舉，內(nèi)蒙古礦業(yè)（集團(tuán)）有限責(zé)任公司當(dāng)選為聯(lián)盟理事長(zhǎng)單位。目前，聯(lián)盟已有會(huì)員單位41家，涵蓋國內(nèi)石墨采選業(yè)、石墨產(chǎn)品研究、開發(fā)、營銷、管理、科研和制造等相關(guān)領(lǐng)域，并已與因成功從石墨中剝離出石墨烯而獲得諾貝爾獎(jiǎng)的英國曼徹斯特大學(xué)教授康斯坦丁·諾沃肖洛夫建立聯(lián)系。

據(jù)悉，內(nèi)蒙古石墨資源豐富，尤其是優(yōu)質(zhì)大鱗片石墨儲(chǔ)量占比高達(dá)80%，位居全國首位。其中查汗木胡魯大磷片石墨礦資源儲(chǔ)量在國內(nèi)乃至世界罕見，又因其獨(dú)特的地理區(qū)位優(yōu)勢(shì)、成熟的產(chǎn)業(yè)政策優(yōu)勢(shì)、低廉的成本競(jìng)爭(zhēng)優(yōu)勢(shì)，為發(fā)展石墨烯、石墨新材料產(chǎn)業(yè)提供了得天獨(dú)厚的便利條件。

相關(guān)人士表示，內(nèi)蒙古石墨聯(lián)盟的成立有利于促進(jìn)各方優(yōu)勢(shì)互補(bǔ)，推動(dòng)協(xié)同創(chuàng)新，促進(jìn)重大產(chǎn)出，實(shí)現(xiàn)共贏發(fā)展，將加快推進(jìn)內(nèi)蒙古區(qū)域創(chuàng)新體系的深度融合，對(duì)內(nèi)蒙古石墨產(chǎn)業(yè)發(fā)展建設(shè)具有里程碑意義。

來自全國78家企業(yè)的代表及各省市區(qū)行業(yè)管理部門共計(jì)200余人參加了成立大會(huì)。大會(huì)期間，還組織了三場(chǎng)高端講座，業(yè)內(nèi)專家及知名人士就石墨深加工、微波石墨烯、石墨烯未來發(fā)展趨勢(shì)等方面做了講解。

Application on aperiodic scale transformation stochastic resonance in nondestructive detection of bolt

HAN Zhi-jun1，CAO Wei2，SUN Xiao-yun2，CHENG En2，WANG Ming-ming2，WU Shi-xing2
（1.China Railway Jian'an Engineering Design Institute，University of Shijiazhuang Railway，Shijiazhuang 050043，China；（2.School of Electrical and Electronic Engineering，University of Shijiazhuang Railway，Shijiazhuang 050043，China）

The weak signal detection of bolt in strong noise has been a real engineering problem. Stochastic resonance is different from conventional denoising method，which uses the noise energy to enhance the weak signal buried in noise.The classical stochastic resonance is only suitable for detecting small parameter signal（signal frequency、signal amplitude、noise intensity is much less than 1），but bolt signal by stress wave reflection method belongs to the damped oscillation transient signal，which has the aperiodicity and large parameter characteristics.Therefore，aperiodic scale transformation stochastic resonance method is adopt to process bolt signal.Based on the advantages and disadvantages of aperiodic stochastic resonance with kurtosis index and correlation coefficient respectively，weighted kurtosis index is used as measurement index.Amplitude is normalization processing and frequency is linear compression in scale transformation stochastic resonance，which meets the requirements of classical stochastic resonance for bolt signal.By processing aperiodic large parameter simulation signal and vibration signal of bolts in practical project，the validity of aperiodic scale transformation stochastic resonance is verified.

anchor bolt；stochastic resonance；weighted kurtosis；scale transformation

TD353

A

1004-4051（2016）09-0154-05

2016-04-30

國家自然科學(xué)基金項(xiàng)目“錨桿受荷條件下聲學(xué)特征及承載力預(yù)測(cè)方法研究”資助（編號(hào)：51274144）；河北省自然科學(xué)基金項(xiàng)目“基于導(dǎo)波的錨桿無損檢測(cè)方法研究”資助（編號(hào)：13205002D）；河北省自然科學(xué)基金項(xiàng)目“基于導(dǎo)納頻率響應(yīng)函數(shù)的錨桿損傷檢測(cè)與定位方法研究”資助（編號(hào)：E2014210075）；河北省省級(jí)研究生創(chuàng)新資助項(xiàng)目“基于RBF神經(jīng)網(wǎng)絡(luò)錨桿錨固質(zhì)量檢測(cè)方法的研究”資助（編號(hào)：yc2016014）

韓志軍（1971-），男，河北趙縣人，高級(jí)工程師，主要研究電氣設(shè)計(jì)或無損檢測(cè)。E-mail：13703112153@163.com。

孫曉云（1971-），女，教授，博士生導(dǎo)師，主要從事礦山安全監(jiān)測(cè)的研究工作。E-mail：sunxyheb@stdu.edu.cn。