徐文焯,章可易,王子靜,趙迎倩,王田昱,徐思哲,何維予

(南京師范大學(xué) 附屬中學(xué)，江蘇 南京 210003)

相對于現(xiàn)有的流量檢測技術(shù)[1-5]，如：電磁流量計、渦街流量計、渦輪流量計、轉(zhuǎn)子流量計、差壓流量計、質(zhì)量流量計等,超聲波流量計具有與被檢測介質(zhì)的非接觸性、通用性廣、適應(yīng)性強以及安裝與維護(hù)簡便等特點，因此超聲波流量計被廣泛地應(yīng)用于現(xiàn)代工業(yè)的各個領(lǐng)域[6-9]，尤其是石油、天然氣、水利、化工等行業(yè). 其發(fā)展經(jīng)歷了漫長的過程，1931 年德國科學(xué)家Ruttgen[1-2]首次提出了利用時差法測量流量的概念，并進(jìn)行了可行性的理論預(yù)測，但在當(dāng)時的儀器精度下難以實現(xiàn). 以此為基礎(chǔ)，此后出現(xiàn)了大量的理論與實驗驗證，但直到1955 年美國科學(xué)家使用“鳴環(huán)法”[1-2]，即：利用聲音在液體中傳播過程中的多次反射所產(chǎn)生的時間差計算液體的流量，才將超聲波流量計從實驗室應(yīng)用到工業(yè)領(lǐng)域，然而由于此檢測電子系統(tǒng)的過于復(fù)雜，使得該流量計并未被廣泛推廣應(yīng)用. 20 世紀(jì)70 年代隨著電子技術(shù)的快速發(fā)展，超聲波流量計的穩(wěn)定性、可靠性以及測量精度得到了很大提高，同時也促進(jìn)了流量測量原理的迅速發(fā)展,如：時差法、多普勒法、波束偏移法、相關(guān)法、空間濾波法及噪聲法等[4,7]. 目前常用的算法主要是時差法和多普勒法，測量精確度也相對較高[8-10]. 時差法測量原理適用于純凈的或含有顆粒直徑小于1mm的流體材料檢測，而多普勒法檢測原理適合于含有較大顆粒的流體檢測. 目前，超聲波流量計的研究主要集中超聲流量計的設(shè)計以及算法的修正上[10-12]，而對管道內(nèi)的流場分布特性上尚有不足、針對流場測試中的不足之處，本文利用有限元方法，構(gòu)建圓形管道中流場的理論分析模型，并數(shù)值分析管道中的空氣流場在風(fēng)機(jī)驅(qū)動下管道內(nèi)流速的分布特性；實驗上利用時差法測量流量原理，實驗測定相同管道中空氣在風(fēng)機(jī)驅(qū)動下流場的平均流速，并與理論分析結(jié)果相印證.

1 理論與實驗裝置

使用2對中心頻率均為40kHz的超聲波傳感器作為探測器，用于檢測口徑為18cm的有機(jī)玻璃導(dǎo)管中的流體速度. 2對超聲波傳感器的安裝方式如圖1 所示，將4個超聲波傳感器兩兩配對構(gòu)成兩組測試探測器，其中每組中各含1個接收和發(fā)射傳感器，將2組傳感器分別安裝在被測導(dǎo)管的兩側(cè)，且每組中的收/發(fā)超聲探測器和另外一組中的發(fā)/收超聲探測器分別一一相對，并保持聲程一致. 這樣當(dāng)管道中有速度為v0的流體流動時，經(jīng)過相應(yīng)聲程所需要的時間為

(1)

其中,D為管徑，θ為超聲發(fā)射方向與流體流速方向的夾角，k為由于儀器誤差所加的修正值,±分別代表超聲傳播的聲速矢量方向與流體速度矢量方向相同或相反.

圖1 實驗裝置示意圖

測試獲得的兩路超聲信號間的時間差為

(2)

由于v0遠(yuǎn)小于c，化簡(2)式可得

(3)

實驗中采用直徑為18 cm的有機(jī)玻璃管作為測量管道，在管道兩側(cè)成64°夾角放置2組超聲波探頭，實驗中采用函數(shù)信號發(fā)生器驅(qū)動超聲波傳感器，分別將采集到的信號輸入前置放大電路，將經(jīng)前置放大的信號輸入數(shù)字示波器顯示并采集記錄到計算機(jī)，由計算機(jī)進(jìn)行最后的信號處理，以獲取管道中流體的流速.

2 實驗結(jié)果與討論

在如圖1所示的實驗裝置的一端使用小型風(fēng)機(jī)向測試管內(nèi)送風(fēng)，使用DF601熱式智能風(fēng)速變送器標(biāo)定風(fēng)機(jī)的輸出風(fēng)速為2.860 m/s. 使用示波器監(jiān)控2組超聲探頭的輸出信號，其結(jié)果如圖2所示，其中虛線所示為超聲傳播的聲速矢量方向沿管軸的分量與流體速度矢量方向相同的超聲

信號，實線為逆向信號. 實線所示信號滯后于虛線所示信號，其滯后的時間量反映了超聲波信號被管道中流場所調(diào)制而導(dǎo)致的超聲波信號在相同聲程條件下傳輸時間的變化特性.

由互相關(guān)理論可知，互相關(guān)函數(shù)反映了2個相互間存在時延信號的相似程度，而該相似度不僅可以用于判斷信號的存在與否，還可以推斷信號到達(dá)的時間和距離，即可以依據(jù)互相關(guān)函數(shù)最大值時的延遲量，計算出兩接收信號之間的時間差. 依據(jù)此原理，計算逆流和順流時所獲取超聲波信號的互相關(guān)特性，其結(jié)果如圖3所示. 依據(jù)兩通路信號的互相關(guān)系數(shù)得到逆流方向所獲取的聲信號比順流方向所獲取的聲信號滯后約26 個采樣點，由采樣周期可得逆流信號比順流信號共滯后5.2 μs，將該結(jié)果代入式(3)，可得管道中流體的流速v0=0.814 m/s. 與實際測量得風(fēng)機(jī)出口的流速相比，此結(jié)果遠(yuǎn)小于風(fēng)機(jī)出口的風(fēng)速.

圖3 所檢測獲取的兩路信號的互相關(guān)特性

以圖1所示的測量管道為模型，利用有限元分析方法構(gòu)建在標(biāo)準(zhǔn)大氣壓下測量管道中的流體模型，在管道中心處加載口徑為50 mm、流速為2.860 m/s的流場，模擬實驗中所使用的風(fēng)機(jī). 在此條件下數(shù)值分析管道內(nèi)流場，其模擬結(jié)果如圖4 所示. 取管道中與實際測試位置相同處的管道橫切面，該橫切面的流速分布如圖5所示，管道中流速的等速分布圖呈中心對稱，且流場中的速度分布呈中心高、邊沿低的分布特性，其中心最大流速約為2.781 m/s，然后逐漸降低，在邊緣處約為0.080 m/s.

圖4 管道中的數(shù)值模擬流場

圖5 管道中橫截面模擬流場的等速線分布

根據(jù)實驗公式[13]，管內(nèi)流速分布可描述為

(4)

將圖5 中的模擬結(jié)果代入式(4)中，使用最小二乘法擬合管道中的數(shù)值流場速度分布特性，則可得到模擬流場中的平均流速為0.859 m/s. 此結(jié)果與實驗中所得到的流場速度0.814 m/s 接近，偏差為5.2%. 由于在實驗分析中，采用的聲速按照340 m/s的標(biāo)準(zhǔn)值進(jìn)行數(shù)值分析，未對實際測試環(huán)境下的聲速進(jìn)行標(biāo)定[14]；再加上忽略了有機(jī)玻璃管道對流動空氣的黏滯性影響. 綜合上述誤差因素，數(shù)值模擬證明，實驗中所獲取的管道中流場的流速為管道中整個流場速度的平均值，且與管道中的流場均值速度基本相符.

從實驗原理以及測試的過程中的信號處理上看,如在以下幾方面進(jìn)行適當(dāng)改進(jìn)將提高實驗結(jié)果的準(zhǔn)確性. 首先,如在實驗進(jìn)行前使用標(biāo)準(zhǔn)流速進(jìn)行預(yù)先的實驗標(biāo)定將有助于實驗結(jié)果的準(zhǔn)確性；其次,采用高質(zhì)量高頻超聲傳感器也將有效地提高測量結(jié)果,但由于高頻超聲在流體中的衰減較快,因而在選擇傳感器時應(yīng)兼顧測量量程的需求,在測程和準(zhǔn)確性之間尋求最佳的匹配；此外,適當(dāng)?shù)靥岣卟蓸铀俾室部商岣邷y量結(jié)果的準(zhǔn)確性,但其對高頻采樣的硬件電路要求較高,易導(dǎo)致成本的非線性增長.

3 結(jié)束語

采用基于有限元的數(shù)值分析方法分析了管道中空氣流場在外界風(fēng)機(jī)作用下的空氣流速分布特性，并使用時差法建立了針對圓柱形管道的超聲波流場測試計量實驗系統(tǒng)，利用該測試系統(tǒng)可獲取管道內(nèi)流場的平均流速，且該流速與基于有限元的數(shù)值模擬結(jié)果基本一致.

致謝:在本課題的完成過程中，南京大學(xué)聲學(xué)所為本項目提供了相關(guān)的實驗儀器、場地，南京大學(xué)聲學(xué)所張仲寧老師為實驗裝置設(shè)計提出了建設(shè)性的改進(jìn)措施.

[1] Kim M. Ultrasonic flowmeter, Proceedings of the International School of Hydrocarbon Measurement[Z]. 1983:422-424.

[2] 本尼迪克R P. 溫度、壓力、流量測量基礎(chǔ)[M]. 周志烈,譯. 北京：國防工業(yè)版社,1985.

[3] 蔡武昌，孫淮清,紀(jì)綱. 流量測量方法及儀表的選用[M]. 北京：化學(xué)工業(yè)出版社，2001.

[4] 陳玲. 新型流量測量儀表的應(yīng)用發(fā)展[J]. 傳感器與微系統(tǒng)，2007，26(6)：8-11.

[5] 孫淮清. 流量儀表技術(shù)發(fā)展與趨勢[J]. 自動化與儀器儀表，2009(1):1-4.

[6] 姜天仕，吳淑珍. 頻差法超聲波流量計[J]. 應(yīng)用聲學(xué)，1982(1):29-34.

[7] 吳斌，趙剛，夏冰. 淺談超聲波流量計的分類與實際應(yīng)用[J]. 中國計量，2006(11):53-55.

[8] 蘭純純. 時差法超聲波流量計的研究[D]. 重慶：重慶大學(xué)，2006.

[9] 李躍忠. 多聲道超聲波氣體流量測量關(guān)鍵技術(shù)研究[D]. 武漢：華中科技大學(xué)，2009.

[10] 李晶，莫德舉. 基于改進(jìn)時差法的超聲流量計[J]. 中國儀器儀表，2004(12):10-12.

[11] Li Xue, Song Zhengxiang. An ultrasound based liquid pressure measurement method in small diameter pipelines considering the installation and temperature [J]. Sensors，2015,15(4):8253-8265.

[12] Sun Jingyuan, Yan Yong. Non-intrusive measurement and hydrodynamics characterization of gas-solid fluidized beds: a review [J]. Measurement Science and Technology，2016,27(11):112001.1-112001.31.

[13] 楊楊. 基于多普勒效應(yīng)法的超聲波流量計的研究[D]. 南昌:南昌大學(xué)，2007.

[14] 買買提熱夏提·買買提,亞森江·吾普爾,復(fù)爾開提·夏爾丁. 簡談空氣中的聲速與溫度關(guān)系[J]. 物理實驗,2007,27(11):43-47.