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

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

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

1 理論與實驗裝置

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

(1)

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

圖1 實驗裝置示意圖

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

(2)

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

(3)

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

2 實驗結果與討論

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

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

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

圖3 所檢測獲取的兩路信號的互相關特性

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

圖4 管道中的數值模擬流場

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

根據實驗公式[13]，管內流速分布可描述為

(4)

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

從實驗原理以及測試的過程中的信號處理上看,如在以下幾方面進行適當改進將提高實驗結果的準確性. 首先,如在實驗進行前使用標準流速進行預先的實驗標定將有助于實驗結果的準確性；其次,采用高質量高頻超聲傳感器也將有效地提高測量結果,但由于高頻超聲在流體中的衰減較快,因而在選擇傳感器時應兼顧測量量程的需求,在測程和準確性之間尋求最佳的匹配；此外,適當地提高采樣速率也可提高測量結果的準確性,但其對高頻采樣的硬件電路要求較高,易導致成本的非線性增長.

3 結束語

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

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

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

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

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

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

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

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

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

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

[9] 李躍忠. 多聲道超聲波氣體流量測量關鍵技術研究[D]. 武漢：華中科技大學，2009.

[10] 李晶，莫德舉. 基于改進時差法的超聲流量計[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] 楊楊. 基于多普勒效應法的超聲波流量計的研究[D]. 南昌:南昌大學，2007.

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