李秉芮， 劉娜， 井上雅弘

(山東科技大學(xué) 安全與環(huán)境工程學(xué)院， 山東 青島 266500)

0 引言

通風(fēng)智能化是指通過動態(tài)采集通風(fēng)信息實(shí)現(xiàn)通風(fēng)狀態(tài)分析、故障診斷和風(fēng)量調(diào)控，是礦井通風(fēng)的主流發(fā)展方向[1-2]。信息采集是通風(fēng)智能化的基礎(chǔ)，然而目前井下風(fēng)速測量常用的機(jī)械式風(fēng)表、壓差式和超聲波渦街式風(fēng)速監(jiān)測傳感器的測量范圍大多為0.3～15 m/s，誤差為±0.2 m/s[3]，啟動風(fēng)速幾乎都大于0.3 m/s[4-8]。《煤礦安全規(guī)程》第一百三十六條規(guī)定：輸送機(jī)巷，采區(qū)進(jìn)、回風(fēng)巷，采煤工作面，掘進(jìn)中的煤巷和半煤巖巷的最低允許風(fēng)速為0.25 m/s；掘進(jìn)中的巖巷和其他通風(fēng)人行巷道的最低允許風(fēng)速為0.15 m/s。現(xiàn)有風(fēng)速傳感器難以滿足《煤礦安全規(guī)程》的風(fēng)速測量要求，同時(shí)也難以滿足通風(fēng)系統(tǒng)異常診斷及智能調(diào)控的需要，是通風(fēng)智能化需要解決的技術(shù)難題之一[9-10]。

超聲波風(fēng)速測量儀可以適應(yīng)惡劣環(huán)境，具有較高的精度和較寬的測量范圍，因而被廣泛應(yīng)用于氣象、軍事、礦山等領(lǐng)域[6,11-13]。利用超聲波測風(fēng)速有渦街法和時(shí)差法2種技術(shù)途徑。之前煤礦普遍使用的超聲波風(fēng)速傳感器多采用渦街法，如GFW15型礦用風(fēng)速傳感器[5]，但渦街法對風(fēng)流穩(wěn)定性及安裝精度要求較高。時(shí)差法不破壞風(fēng)流流場，無壓力損失，幾乎不受風(fēng)流溫度、壓力等因素影響，是未來礦井風(fēng)速測量技術(shù)的發(fā)展方向。國外生產(chǎn)的風(fēng)流監(jiān)測設(shè)備如Accutron MAQS和Zephyr AQSTM均采用超聲波時(shí)差法。國內(nèi)早期開發(fā)的基于超聲波時(shí)差法的礦用風(fēng)速測量儀使用2對超聲波換能器，通過交替發(fā)射與接收超聲波信號實(shí)現(xiàn)風(fēng)速測量[14]。文獻(xiàn)[15]在此基礎(chǔ)上進(jìn)行了改進(jìn)，使用1對超聲波換能器完成風(fēng)速測量，最低測量風(fēng)速達(dá)0.2 m/s，性能優(yōu)于傳統(tǒng)的風(fēng)速測量儀，但由于風(fēng)速越低，測量誤差越大，該風(fēng)速測量儀仍不能滿足礦井低風(fēng)速條件下精準(zhǔn)測風(fēng)需要。針對該問題，本文基于超聲波時(shí)差法風(fēng)速測量原理，設(shè)計(jì)了一種具有低啟動風(fēng)速的高精度礦用超聲波風(fēng)速測量儀。

1 超聲波時(shí)差法風(fēng)速測量原理

超聲波在空氣中的傳播速度與風(fēng)速疊加，在順風(fēng)方向上傳播速度加快，在逆風(fēng)方向上傳播速度變慢，在一定檢測條件下與風(fēng)速函數(shù)對應(yīng)，可通過計(jì)算求出精確風(fēng)速。超聲波在空氣中的傳播速度受溫度的影響，但若采用2個(gè)通道分別檢測2個(gè)相反方向的超聲波傳播速度，則溫度對超聲波傳播速度產(chǎn)生的影響可以忽略不計(jì)。超聲波時(shí)差法風(fēng)速測量原理如圖1所示。在巷道相對的兩側(cè)巷壁分別布置A1，A2和B1，B22組超聲波探頭。其中A1，B2為超聲波發(fā)射探頭，A2，B1為超聲波接收探頭。發(fā)射探頭與接收探頭的距離為L，其連線與風(fēng)流方向的夾角為α。

圖1 超聲波時(shí)差法風(fēng)速測量原理Fig.1 Principle of wind speed measurement of ultrasonic time-difference method

設(shè)巷道的風(fēng)流方向?yàn)閄方向(巷道軸方向)，與風(fēng)流方向垂直的方向?yàn)閅方向，沿Y方向的風(fēng)速分布為v(Y),當(dāng)A1，B2同時(shí)發(fā)射超聲波時(shí)，順風(fēng)方向超聲波由A1到B1的傳播時(shí)間為

(1)

式中:C為巷道空氣中超聲波的傳播速度，m/s，C=331.45+0.607t,t為巷道空氣溫度，℃;l為測距的積分變量。

逆風(fēng)方向由B2到A2的傳播時(shí)間為

(2)

用平均風(fēng)速vm表示沿Y方向的風(fēng)速分布v(Y)，并代入式(1)和式(2)得

(3)

(4)

由式(3)和式(4)可得

(5)

(6)

當(dāng)L已知時(shí)，測量順風(fēng)與逆風(fēng)狀態(tài)超聲波的傳播時(shí)間，根據(jù)式(5)即可求出沿Y方向的平均風(fēng)速，且與風(fēng)流溫度無關(guān)。當(dāng)超聲波探頭受外界影響而發(fā)生位移時(shí)，可根據(jù)式(6)計(jì)算位移量并校準(zhǔn)探頭間距。

2 超聲波風(fēng)速測量儀設(shè)計(jì)和信號調(diào)理

2.1 超聲波風(fēng)速測量儀設(shè)計(jì)

超聲波風(fēng)速測量儀工作原理如圖2所示。超聲波發(fā)射探頭A1，B2和接收探頭A2，B1交錯(cuò)設(shè)置，形成2個(gè)超聲波傳播通路，分別用于測算超聲波在特定風(fēng)速場中的順風(fēng)和逆風(fēng)傳播時(shí)間。

圖2 超聲波風(fēng)速測量儀工作原理Fig.2 Working principle of ultrasonic anemometer

風(fēng)速測量儀由單片機(jī)、超聲波收發(fā)模塊及紅外定向裝置等組成。單片機(jī)集放大電路、計(jì)時(shí)器、振蕩器、觸發(fā)電路、濾波器、積分電路、信號轉(zhuǎn)換器于一體，完成時(shí)差法風(fēng)速測量。紅外定向裝置由紅外線發(fā)射器E1，E2及標(biāo)志區(qū)F1，F(xiàn)2組成，位于超聲波探頭A2，B2上方。E1，E2分別發(fā)射紅外線至對面標(biāo)志區(qū)F1，F(xiàn)2，進(jìn)行風(fēng)速測量儀安裝方向校準(zhǔn)，以減小由安裝方向偏離引起的測量誤差。

風(fēng)速測量儀工作電壓為5 V，通過觸發(fā)電路觸發(fā)振蕩器，產(chǎn)生具有一定頻率的脈沖信號，激勵超聲波發(fā)射探頭發(fā)射規(guī)定頻率的超聲波，即將電信號轉(zhuǎn)換為超聲波信號。在發(fā)射超聲波的同時(shí)，單片機(jī)啟動計(jì)時(shí)器。2個(gè)超聲波接收探頭接收超聲波信號，并將其轉(zhuǎn)換為電信號輸入單片機(jī)，同時(shí)中斷計(jì)時(shí)器。信號經(jīng)過放大、濾波處理后，再經(jīng)過積分電路處理成方波。低電平持續(xù)時(shí)間就是超聲波在空氣中的傳播時(shí)間，在順風(fēng)與逆風(fēng)條件下低電平持續(xù)時(shí)間不同，由此可計(jì)算出超聲波傳播時(shí)差，根據(jù)式(5) 可計(jì)算出風(fēng)速并判斷風(fēng)向，然后實(shí)時(shí)輸出風(fēng)速。

2.2 超聲波信號調(diào)理

超聲波信號調(diào)理是指對超聲波信號進(jìn)行與或非邏輯運(yùn)算。對振蕩器產(chǎn)生的電信號進(jìn)行信號處理及門電路運(yùn)算，得到傳播時(shí)間差造成的脈沖寬度，再通過濾波器及積分電路將電信號轉(zhuǎn)換為穩(wěn)定的輸出電壓。

超聲波風(fēng)速測量儀信號調(diào)理如圖3所示。超聲波風(fēng)速測量儀由IO口輸出一個(gè)10 μs的高電平信號至TRIG(控制端)，激勵超聲波發(fā)射探頭A1，B2發(fā)射8個(gè)40 kHz方波(波形O)，同時(shí)A1和B2的TRIG輸出一個(gè)低電平。B1接收到A1發(fā)射的超聲波后，換能器產(chǎn)生脈沖信號a，同時(shí)發(fā)生電平跳變，從B1的ECHO(接收端)輸出一個(gè)高電平，經(jīng)積分電路整形為方波(波形P)并存儲，低電平寬度為超聲波在順風(fēng)條件下的傳播時(shí)間T1。A2接收到B2發(fā)射的超聲波后，換能器產(chǎn)生脈沖信號b，同時(shí)發(fā)生電平跳變，從A2的ECHO輸出一個(gè)高電平，經(jīng)積分電路整形為方波(波形Q)并存儲, 低電平寬度為超聲波在逆風(fēng)條件下的傳播時(shí)間T2。圖3中，Tr為超聲波風(fēng)速測量儀的1個(gè)發(fā)射周期，即從超聲波發(fā)射探頭A1，B2第1次發(fā)射8個(gè)40 kHz方波到第2次發(fā)射方波的時(shí)間；T0為待機(jī)時(shí)間，即超聲波發(fā)射探頭發(fā)射8個(gè)40 kHz方波后至下一次發(fā)射的等待時(shí)間。

圖3 超聲波風(fēng)速測量儀信號調(diào)理Fig.3 Signal conditioning of ultrasonic anemometer

3 性能參數(shù)設(shè)定及其影響因素分析

超聲波風(fēng)速測量儀的測量精度除了受頻率、待機(jī)時(shí)間等性能參數(shù)影響外，還受探頭方向、聲噪聲等因素干擾。

3.1 超聲波頻率

超聲波在空氣中傳播的過程中，強(qiáng)度會隨距離增大而衰減，頻率越高，衰減率越大，長距離測量越困難；但低頻超聲波的工作噪聲較大，易對測量精度產(chǎn)生影響。按照圖1所示的超聲波風(fēng)速測量儀布置方式，10 m以內(nèi)的測距基本可以滿足井下巷道測風(fēng)要求。根據(jù)理論和實(shí)踐經(jīng)驗(yàn)，10 m以內(nèi)測距最適宜的超聲波頻率為40 kHz[16]，因此，本文選用特性頻率為40 kHz的超聲波換能器。

3.2 待機(jī)時(shí)間

若超聲波風(fēng)速測量儀待機(jī)時(shí)間過短，則測量空間內(nèi)發(fā)射探頭的余振、回響信號等聲噪聲會降低測量精度；若待機(jī)時(shí)間過長，則測量響應(yīng)速度慢。因此，需設(shè)置合理的待機(jī)時(shí)間。設(shè)探頭間距L=10 m，α=0，空氣溫度t=0 ℃，巷道平均風(fēng)速vm=10 m/s，考慮反射波的影響，超聲波從發(fā)射到返回所需的傳播時(shí)間為

(7)

由式(7)可知，待機(jī)時(shí)間大于60 ms時(shí)，可避免反射波的影響。若探頭間距L=1 m，則待機(jī)時(shí)間應(yīng)大于6 ms。

設(shè)定探頭間距L=1 m，待機(jī)時(shí)間分別為50，100，200 ms，連續(xù)進(jìn)行300次測試，風(fēng)速標(biāo)準(zhǔn)差及風(fēng)速波動范圍見表1?？梢钥闯?，當(dāng)待機(jī)時(shí)間大于超聲波折返所需傳播時(shí)間時(shí)，風(fēng)速標(biāo)準(zhǔn)差很小，實(shí)測值波動小，表明測量穩(wěn)定性好，幾乎不受待機(jī)時(shí)間的影響?？紤]風(fēng)速測量儀的響應(yīng)速度，設(shè)定待機(jī)時(shí)間為100 ms。

表1 不同待機(jī)時(shí)間的風(fēng)速標(biāo)準(zhǔn)差及波動范圍Table 1 Standard deviation and fluctuation range of wind speed under different stand-by time

3.3 超聲波探頭方向性

超聲波具有良好的方向性，超聲波場中至波源充分遠(yuǎn)處同一橫截面上各點(diǎn)的聲壓不同，軸線上的聲壓最高[17]。一般實(shí)用性的超聲波探頭半擴(kuò)散角為4～15°，即測距為1 m時(shí)，超聲波探頭方向的最大允許錯(cuò)位為0.268 m。因此，設(shè)超聲波探頭方向錯(cuò)位0.25 m，采用4種布置方式進(jìn)行室內(nèi)實(shí)驗(yàn)，如圖4所示。超聲波連續(xù)發(fā)射100次。

(a) 布置方式1

(b) 布置方式2

(c) 布置方式3

(d) 布置方式4

不同超聲波探頭布置方式的實(shí)測風(fēng)速如圖5所示，風(fēng)速標(biāo)準(zhǔn)差及波動范圍見表2。

(a) 布置方式1實(shí)測風(fēng)速

(b) 布置方式2實(shí)測風(fēng)速

(c) 布置方式3實(shí)測風(fēng)速

(d) 布置方式4實(shí)測風(fēng)速

從圖5和表2可看出，超聲波探頭方向錯(cuò)位越嚴(yán)重，風(fēng)速標(biāo)準(zhǔn)差和波動范圍越大，表明探頭方向?qū)y量的準(zhǔn)確性與可靠性有一定影響。為減小探頭方向錯(cuò)位造成的誤差，本超聲波風(fēng)速測量儀設(shè)置了紅外定向裝置。

表2 不同超聲波探頭布置方式的風(fēng)速標(biāo)準(zhǔn)差及波動范圍Table 2 Standard deviation and fluctuation range of wind speed under different ultrasonic probe arrangements m/s

3.4 數(shù)據(jù)平滑處理

對原始風(fēng)速波形曲線進(jìn)行平滑處理的目的是消除隨機(jī)噪聲，提高輸出風(fēng)速的穩(wěn)定性。通過比較分析，本文選用滑動平均法，其原理是對n個(gè)非平穩(wěn)數(shù)據(jù){yi}(i=0,1,…,n-1)，逐個(gè)滑動取m個(gè)相鄰數(shù)據(jù)計(jì)算算術(shù)平均值：

(8)

式中：fk為第k個(gè)算術(shù)平均值；yk+i為原始風(fēng)速測量值。

m的取值直接影響平滑效果:m值較大有利于抑制頻繁隨機(jī)起伏的隨機(jī)誤差，但可能導(dǎo)致高頻變化的確定性成分被平均；m值較小不利于抑制隨機(jī)誤差。數(shù)據(jù)處理中多設(shè)置m值為5～11。為確認(rèn)m值對風(fēng)速波形曲線平滑效果的影響，分別取m值為4，6，8，10，12，對測距1.0 m的原始風(fēng)速波形進(jìn)行數(shù)據(jù)平滑處理，結(jié)果如圖6所示。

(a) 原始風(fēng)速波形

(b) 平滑處理后波形

m值與風(fēng)速標(biāo)準(zhǔn)差的關(guān)系如圖7所示。可以看出，m值越大，風(fēng)速標(biāo)準(zhǔn)差越小，風(fēng)速變化范圍越小，風(fēng)速越穩(wěn)定，但風(fēng)速測量儀響應(yīng)速度變慢。當(dāng)m≥8時(shí)，風(fēng)速標(biāo)準(zhǔn)差變化小，可以認(rèn)為m對風(fēng)速平均值的影響很小，因此，本文取m=8進(jìn)行數(shù)據(jù)平滑處理。

圖7 m值與風(fēng)速標(biāo)準(zhǔn)差的關(guān)系Fig.7 Relation between m and the standard deviation of wind speed

4 性能指標(biāo)分析

(1) 分辨率。設(shè)α=0，對式(5)求T1的導(dǎo)數(shù)，得

(9)

(10)

已知超聲波風(fēng)速測量儀計(jì)時(shí)器的分辨率為1 μs，設(shè)L=5 m，C=340 m/s，vm=5 m/s，代入式(10)得dvm=0.011 9 m/s，即超聲波風(fēng)速測量儀的測量分辨率為0.01 m/s。

(2) 測量誤差。已知超聲波探頭的測距誤差|dl|≤3 mm，換算為超聲波傳播時(shí)間誤差得

(11)

將dT1代入式(10)得dvm=0.103 m/s，則超聲波風(fēng)速測量儀誤差為0.1 m/s。

(3) 啟動風(fēng)速。因?yàn)椴捎梅菣C(jī)械式傳動方式，理論上超聲波風(fēng)速測量儀的啟動風(fēng)速應(yīng)該等于其分辨率。根據(jù)實(shí)驗(yàn)室內(nèi)無風(fēng)、微風(fēng)狀態(tài)的實(shí)測，在風(fēng)速小于0.1 m/s的情況下，超聲波風(fēng)速測量儀均獲得了穩(wěn)定的實(shí)測值，證明其啟動風(fēng)速低于0.1 m/s。

5 結(jié)論

(1) 基于超聲波時(shí)差法風(fēng)速測量原理，設(shè)計(jì)了具有低啟動風(fēng)速的高精度礦用超聲波風(fēng)速測量儀。通過2個(gè)發(fā)射探頭同時(shí)發(fā)射超聲波，形成2個(gè)超聲波傳播通路，溫度對風(fēng)速測量的影響可以忽略不計(jì)。

(2) 采用紅外定向裝置，可減小超聲波探頭方向錯(cuò)位引起的測量誤差。選用滑動平均法對原始風(fēng)速波形曲線進(jìn)行平滑處理，以消除隨機(jī)噪聲，提高輸出風(fēng)速的穩(wěn)定性。

(3) 理論計(jì)算和實(shí)驗(yàn)結(jié)果表明：該風(fēng)速測量儀的分辨率為0.01 m/s，測量誤差為±0.1 m/s，啟動風(fēng)速低于0.1 m/s，可以滿足礦井低風(fēng)速條件下的精準(zhǔn)測風(fēng)要求，今后將通過實(shí)踐進(jìn)行進(jìn)一步驗(yàn)證。