程 明，馮紅亮，肖定國，韓德龍

一種便攜式高精度超聲流量儀

程 明1，馮紅亮2，肖定國1，韓德龍2

（1 北京理工大學(xué)機械與車輛學(xué)院 北京 100081 2 北京遙測技術(shù)研究所 北京 100094）

研究一種超聲波傳播方向平行于液體流動方向的流量管和以此為基礎(chǔ)的超聲流量儀。針對流量系統(tǒng)的新型流量管結(jié)構(gòu)推導(dǎo)兩種流量算法并進行誤差分析，其中一種算法可以減小由聲速變化引起的誤差，從而提高系統(tǒng)的測量穩(wěn)定性，也為流量管的改進方向提供了思路。使用測顯比曲線修正管道中流體流速分布不均勻的影響以提升測量準確性，試驗證明，在測量范圍內(nèi)，系統(tǒng)測量的平均示值誤差不超過0.2%。

超聲流量計；流量算法；管道流體；流量計校準

引 言

超聲波流量計作為一種流量測量儀表不僅應(yīng)用于工業(yè)流體計量，在醫(yī)療、河流、海洋監(jiān)測等計量領(lǐng)域也有廣泛的應(yīng)用[1]。超聲波流量計可用于多種介質(zhì)的流量測量，尤其對具有腐蝕性、高粘度、非導(dǎo)電等特性的流體流量測量有顯著的優(yōu)勢[2]。

當超聲波在流體中傳播時，由于超聲波速與流體流速在速度場上的疊加，導(dǎo)致超聲波束在相同的距離內(nèi)順著流體流動方向傳播與逆著流體流動方向傳播所用的時間不一樣，順流時所需時間短，逆流時所需時間長[3]。時差法超聲流量測量正是通過測量超聲波在順流與逆流中的時間差來計算流體的流動速度。為了降低順逆流時間的測量誤差，時差法需要超聲波有較長的超聲傳播路程，才能在較大流速下也可以獲得較高的測量精度[4]。目前，成熟的時差法和多普勒超聲流量計多采用外夾式，時差法超聲流量計在標定后測量精確度可以達到0.33%~0.5%[5]，多普勒超聲流量計測量精確度能達到0.4%左右[6]。但是，外夾式超聲流量計會遇到超聲波經(jīng)歷多次界面反射和透射導(dǎo)致的相對信號較弱的問題以及管道狀態(tài)較為復(fù)雜（管道厚度、腐蝕等情況各不相同）導(dǎo)致的應(yīng)用均一性很難得到保障的問題。另外，實際使用的管道直徑普遍較小，超聲波傳播路程很短，這些都會導(dǎo)致外夾式超聲流量計測量誤差較大，測量精度較低。

對于超聲流量測量來說，準確獲得流體流動速度是關(guān)鍵，但是管道流體的流速分布并不均勻，需要使用修正或者標定的方法獲得不同速度下的截面平均流速，才能準確地得到流體的體積流量[7]。流量管中液體的流速分布與液體的流動狀態(tài)有關(guān)，關(guān)于管道流體的截面流速分布也有不少實驗、理論研究和仿真研究[8-10]，甚至根據(jù)經(jīng)驗公式給出了截面平均流速與最大流速之比隨雷諾數(shù)變化的計算公式[11]。目前，關(guān)于截面平均流速的公式大多只能針對特定流動狀態(tài)，且本身存在較大的誤差，當流速變化范圍大、跨越多種流動狀態(tài)時，利用理論模型，根據(jù)測出的流速推算管截面流體的平均流速誤差較大。

本文提出了一種基于超聲波平行于液體流動方向的流量管（以下簡稱聲束平行流量管）的超聲流量儀，推導(dǎo)了該流量管結(jié)構(gòu)下的兩種流速流量計算方法。對比分析了兩種算法的誤差，并提出改進的方向，實現(xiàn)了測量時流體外部環(huán)境的統(tǒng)一和超聲波在流體中傳播時間的大幅提升，這有利于提高測量精度。使用標定曲線擬合公式的方式修正管道內(nèi)徑測量誤差等系統(tǒng)誤差和管道流體流速分布不均勻的影響，實現(xiàn)了較高的測量精度。

1 超聲流量儀的系統(tǒng)組成

本文研究的超聲波流量儀的特點是采用超聲波平行于液體流動方向的流量管進行測量，測量系統(tǒng)的構(gòu)成如圖1所示。超聲波探頭A和探頭B分別安裝在流量管兩端，兩探頭的聲束軸與流體的流動方向平行。在流量管長度已知的情況下，探頭A發(fā)射的超聲波被探頭B接收，可以計算獲得超聲波順流傳播時間；探頭B發(fā)射的超聲波被探頭A接收，可以計算獲得超聲波逆流傳播時間，得到的超聲波的順逆流傳播時間可以用于計算液體流速流量值。

圖1 超聲流量儀的模塊

超聲流量儀主要組成部分包括：超聲波探頭（探頭A、探頭B）、超聲波收/發(fā)通道切換模塊、單通道超聲檢測卡、便攜式計算機、流量管及相應(yīng)的計算顯示軟件等。進行流量測量時，首先是探頭A與超聲檢測卡的脈沖發(fā)射端接通，探頭B與超聲檢測卡的脈沖接收端接通，探頭A發(fā)射的超聲波通過流量管中的流體被探頭B接收。接收端口接收到超聲波信號后，由集成在超聲檢測卡中的相應(yīng)模塊對接收到的脈沖信號進行處理（包括放大、濾波等）、A/D轉(zhuǎn)換、以及數(shù)據(jù)采集等操作，計算出該方向超聲波的傳播時間。經(jīng)過一次脈沖發(fā)射與接收后，通道切換模塊在計算機和超聲檢測卡的控制下自動切換發(fā)射/接收通道，更改為探頭B發(fā)射、探頭A接收，重復(fù)進行上述動作進行第二次數(shù)據(jù)采集和計算。計算機通過以太網(wǎng)與超聲檢測卡進行通訊，對順逆流條件下超聲檢測卡兩次采集存儲的數(shù)據(jù)進行運算處理得出流量值并顯示，使用包括拉格朗日插值和互相關(guān)算法提高超聲波傳播時間的測量精度。該系統(tǒng)支持根據(jù)選用流量管型號和被測液體參數(shù)進行參數(shù)配置。

基于此超聲流量測量系統(tǒng)進行了流速流量相關(guān)算法的研究和實驗。該系統(tǒng)的工控系統(tǒng)部分進行了便攜式集成設(shè)計，工控系統(tǒng)的儀器機箱尺寸為360mm×276mm×103mm，重量小于10kg，具有體積小、便攜性操控性好、精度高的特點；流量管進行了模塊化設(shè)計，兩端安裝超聲探頭的基座部分可以重復(fù)使用，流量管道部分可以選用預(yù)制的不同管徑和長度的直管模塊以滿足不同的測量測試需求。

本文超聲流量測量系統(tǒng)使用直管型號為DN20的流量管，直管長度為400mm。超聲波探頭選用了奧林巴斯公司的C540-SM型水浸探頭，中心頻率為1MHz。經(jīng)過實驗驗證，對于水介質(zhì)流體，在常溫下、聲程為800mm時，無論流量管中流體處于順流還是逆流狀態(tài)，超聲流量測量系統(tǒng)依然可以獲得穩(wěn)定性良好、強度滿足要求的超聲信號。

2 流量管結(jié)構(gòu)及流量計算原理

2.1 聲束平行流量管結(jié)構(gòu)

聲束平行流量管的兩端是對稱結(jié)構(gòu)，左端的結(jié)構(gòu)如圖2所示。該流量管的結(jié)構(gòu)能夠避免流體徑向流速對測量的干擾，同時可以通過增大超聲波傳播距離來提高超聲波順逆流傳播時差測量精度。流量管的整體長度隨中間直管的長度不同而不同，本文使用的管道內(nèi)徑為20mm，直管長度可以在400mm~600mm范圍內(nèi)變化。

相較于外夾式“Z”型超聲探頭布置結(jié)構(gòu)，新型流量管由于使用了平行于流動方向的超聲探頭布置方式，能在較小管徑的直管中使用，同時也獲得了較大的聲程，因此在小管徑小流量的測量應(yīng)用中優(yōu)勢明顯。

1-中間直管，2-連接器，4-45°三通，5-連接螺釘，6-超聲探頭，9-探頭座，10-連接件，3,7,8-密封件

這兩種流量管在具有相同的管道直徑、靜流聲速、流速、只有布置結(jié)構(gòu)不相同的理想情況下，二者關(guān)鍵參數(shù)的計算對比結(jié)果如表1所示。

由表1可知，在相同管徑下，新型結(jié)構(gòu)具有較大的聲程，同時新流量管可以通過加長中間直流管的方式持續(xù)增大聲程，進而增大順逆流聲時差，減小流量測量誤差。

表1 同等條件下兩種結(jié)構(gòu)理論聲時差計算

2.2 聲束平行流量管的流量測量模型

聲束平行流量管流速測量原理如圖3所示。根據(jù)流量管的結(jié)構(gòu)和探頭安裝方式，推導(dǎo)流量計算公式如下。

設(shè)在測量開始時流體充滿整個流量管，流體從靠近探頭A的進水口流入，經(jīng)過流量管從靠近探頭B的出水口流出。

圖3 新型流量管流速測量示意圖

先對測量直管段內(nèi)的流體流動做簡化處理：在靠近A、B探頭兩端的Aa段和bB段，近似認為流體處于靜止狀態(tài)；測量流體在直管ab段從左向右的軸線流動速度（進出水管道的軸心線與測量直管的軸心線相交與a、b兩點）?；谝陨虾喕幚?，探頭A發(fā)射的超聲波被探頭B接收時，得到順流超聲波傳播時間為

式中，表示軸向流動距離ab，1表示兩端流體近似靜止段長度，c表示靜流中聲波傳播速度，m表示流體的流速。

探頭B發(fā)射、探頭A接收時，得到逆流傳播時間為

聯(lián)立式（1）和式（2）得

由于流量管內(nèi)流體的流速分布并不均勻，在不同流動速度下其流速沿管截面徑向是變化的，因此，由式（4）計算得到的m在一般情況下不等于平均流速，而體積流量計算需要使用管截面內(nèi)的平均流速，因此還需要對所測的流速m進行修正，才能較準確的計算出管道流量。管道的體積流量可由下式計算：

式中，為管道體積流量，為管道內(nèi)徑，v為修正系數(shù)。

靜流聲速c可利用流量管在靜流條件下根據(jù)聲傳播時間和傳播距離現(xiàn)場測出，測量計算得到靜流聲速為

式中，0為靜流條件下超聲波從探頭A到探頭B的或從探頭B到探頭A的傳播時間。

此外，液體流速和流量還可用另一種改進的算式計算：

聯(lián)立式（7）和式（8）得

2.3 誤差分析

分析流量的計算式（5）可知，管道內(nèi)徑的測量誤差和流速m的測量誤差都會造成流量值的誤差。管道內(nèi)徑的測量誤差為固定值，可以使用標定的方法包含到修正系數(shù)v中；液體流動速度m的測量誤差是流量測量的主要誤差來源，因此主要分析速度m的測量誤差。

對式（4）進行不確定度分析，可得相對不確定為

式中，為測得流速的不確定度，c為聲速的不確定度，σAB和σBA分別為測得順逆流聲傳播時間的不確定度。

管內(nèi)液體流速為1m/s左右時聲時差大約為300ns，系統(tǒng)對時間的分辨力為1ns，假定實驗中流體溫度為25℃，實驗中溫度升高5℃，將會導(dǎo)致聲速增大12.4m/s，此時流速m的相對誤差大約為1.2%。不考慮系統(tǒng)對時間的分辨力，只考慮溫度變化時m的相對誤差大約為0.8%。

由此可見，按式（4）計算流速時，由溫度變化引起的靜流聲速c的變化對流速測量準確性的影響較大。

對式（9）進行不確定度分析，可得

由式（12）可知，當1趨近于0時，由溫度變化引起的靜流聲速c的變化對流速測量準確性的影響很弱，式（12）轉(zhuǎn)化為

分析式（11）、式（12）和式（13）可知，在1足夠小的情況下，式（12）表示的相對誤差明顯小于式（11），說明用式（9）計算流速比用式（4）計算流速更加穩(wěn)定準確，更有利于消除由溫度變化引起的靜流聲速c的變化對流速測量準確性的影響。本文系統(tǒng)采用式（10）計算流量。

3 試驗和校準

3.1 系統(tǒng)標定

如上文所述，為了獲得準確的平均流速，本文沒有使用精度不高且有局限的仿真模型和經(jīng)驗公式，而是采用了標準流量裝置標定的方法來確定修正系數(shù)v。

圖4 標定系統(tǒng)模塊和校準現(xiàn)場

標定的方法如圖4所示，超聲流量儀的流量管通過連接管接入流量標準裝置的回路中，通過對比流量標準裝置的流量和超聲流量儀的顯示流量來標定超聲流量儀。用于標定的標準流量系統(tǒng)是采用計時稱重法的流量測量系統(tǒng)。標定時使用的液體為純水，溫度為室溫。由于計時稱重法流量測量值準確性高，視為真實值。標定時記錄不同流速下超聲流量儀顯示的流量值s和計時稱重法測得的流量值，單位均為g/s計算它們的比值，該比值即為修正系數(shù)v，將v的倒數(shù)記為。對不同流量下標定得到的值與超聲流量儀顯示值s進行5次多項式擬合，并將擬合的方程用于修正測量結(jié)果，計算對應(yīng)所測流量準確值。

標定得到的s曲線如圖5所示。由于流量管直徑的測量誤差是固定值，所以該曲線的變化趨勢也反映了流量管中截面等效平均流速a和最大流速m的關(guān)系。因此，超聲流量儀配合高精度流量標準裝置還可以進行流速分布的相關(guān)研究。

3.2 系統(tǒng)校準

為了評價超聲流量儀測量液體流量的準確性，使用了中國計量科學(xué)研究院的0.01m3/h~200m3/h水流量基準裝置進行評價，該裝置的證書編號為國基證[2017]第015號，測量不確定度為0.05%，采用計時稱重法計算流量。實驗采用的液體為25.5℃的水，校準時的系統(tǒng)連接的表前直管長度大于40倍直徑，呈V型安裝，連接形式與標定時相同。

圖5 Qs-k曲線

校準數(shù)據(jù)如表2所示，計時稱重流量標準裝置測得的流量值稱為流量標準值，超聲流量儀顯示的流量值稱為流量示值，示值誤差為流量示值與流量標準值的差值和流量標準值的比值百分數(shù)。對覆蓋測量范圍的3個不同流量的測量準確性進行了評價，每個流量測量3次，記錄每次的示值誤差，3次示值誤差的平均值為該流量測量準確性的評價結(jié)果。從表中可以看出，在整個測量范圍內(nèi)，3個流量測點的最大平均示值誤差為0.17%，表明該超聲流量儀在測量范圍內(nèi)有不高于0.2%的平均測量誤差。

表2 校準評價數(shù)據(jù)

4 結(jié) 論

本文提出了基于聲束平行流量管的超聲流量儀，研究了在該流量管結(jié)構(gòu)下流體流速的計算方法和流量計算修正方法，并使用高精度計時稱重流量測量標準裝置進行了校準和評價，得到結(jié)論如下：

①聲束平行流量管超聲探頭平行液體流動方向的布置方式經(jīng)過簡化計算，能在較小管徑上獲得更長的聲程和更大的聲時差。這降低了測量誤差，本文所述系統(tǒng)測量平均示值誤差不高于0.2%。

②根據(jù)聲束平行流量管的結(jié)構(gòu)推導(dǎo)了時差法流量測量的流速計算方法，改進了流速計算算式，削弱了由溫度變化引起的靜流聲速c的變化對流速測量準確性的影響，降低了測量誤差。

③提出了標定得到s曲線，擬合s關(guān)系方程修正流速分布誤差的方法，并分析指出該曲線不僅可以修正管道流體徑向截面流速分布不均而造成的流速測量誤差，還可用于管道內(nèi)流體流速分布的反演。

[1] 薛蕾, 彭亮, 張超. 超聲流量計在液體實流標定中的難點分析與解決方案[J]. 科技經(jīng)濟導(dǎo)刊, 2019, 27(10): 179. XUE Lei, PENG Liang, ZHANG Chao. Analysis and solution of difficulties in ultrasonic flowmeter calibration of liquid real flow[J]. Technology and Economic Guide, 2019, 27(10): 179.

[2] 李定川. 超聲波流量計工程應(yīng)用概述[J]. 智慧工廠, 2016(6): 84–88. LI Dingchuan. The composition and basic principle of ultrasonic flow meter and its application in production[J]. Smart factory, 2016(6): 84–88.

[3] SVILAINIS L. The time-of-flight estimation accuracy versus digitization parameters[J]. Ultragarsas, 2016, 63(1): 12–17.

[4] ORVATINIA M. Improvement of accuracy in multi-path ultrasonic flow meters[J]. Sensors & Transducers, 2019, 231(3): 1–9.

[5] 熊光德. 新型天然氣超聲波流量計量技術(shù)[J]. 天然氣與石油, 2002(2): 57–61. XIONG Guangde. New technology of natural gas ultrasonic flow rate measurement[J]. Natural Gas and Oil, 2002(2): 57–61

[6] 王興國, 張路鑫, 堯亮富, 等. 管道內(nèi)流體流量的超聲多普勒測量[J]. 電子測量與儀器學(xué)報, 2019, 33(10): 96–103. WANG Xingguo, ZHANG Luxin, YAO Liangfu, et al. Measurement of pipeline flow using ultrasonic Doppler method[J]. Journal of Electronic Measurement and Instrumentation, 2019, 33(10): 96–103.

[7] AGOLOM M O, LUCAS G, WEBILOR R O. Measurement of velocity profiles in transient single and multiphase flows using inductive flow tomography[J]. Flow Measurement and Instrumentation, 2018, 62: 246–254.

[8] 劉敦利, 胡鶴鳴. 圓形管道截面的流速分布實驗與分析[J]. 儀器儀表標準化與計量, 2019(6): 43–45. LIU Dunli, HU Heming. Experiment and analysis on velocity distribution of circular pipeline cross-section[J]. Instrument Standardization & Metrology, 2019(6): 43–45.

[9] 宋振華. 超聲流量計流場分析及導(dǎo)波流量檢測方法研究[D]. 北京: 北京工業(yè)大學(xué), 2017. SONG Zhenhua. Study on flow field analysis of ultrasonic flowmeter and flow measurement based on guided waves[D]. Beijing: Beijing University of Technology, 2017.

[10] ESCUE A, CUI J. Comparison of turbulence models in simulating swirling pipe flows[J]. Applied Mathematical Modelling, 2010, 34(10): 2840–2849.

[11] 何龍德. 管道截面的流速分布參數(shù)及其在渦街流量計標定中的應(yīng)用[J]. 氣動實驗與測量控制, 1995(2): 67–73. HE Longde. Flow velocity parameters of pipeline section and its application in vortex flowmeter calibration[J]. Journal of Experiments in Fluid Mechanics, 1995(2): 67–73.

Portable high-precision ultrasonic flow measurement system

CHENG Ming1, FENG Hongliang2, XIAO Dingguo1, HAN Delong2

（1. School of Mechanical Engineering, Beijing Institute of Technology, Beijing, 100081, China; 2. Beijing Research Institute of Telemetry, Beijing 100094, China）

This paper presents a portable ultrasonic flow measurement system based on a new flow tube with the ultrasonic sensor parallel to the direction of liquid flow. Based on the new flow tube structure of the flow measurement system, two flow calculating algorithms are deduced and error analysis is performed. It is found that one of the algorithms can reduce the error caused by the change in sound velocity, which improves the measurement stability of the system and can also suggest ideas for improvement. This system uses the Qs-k ratio curve to correct the influence of uneven flow velocity distribution of the moving fluid. Within the measurement range, experiment shows that the error of the measured value of the system does not exceed 0.2%.

Ultrasonic flowmeter; Flow algorithm; Pipeline fluid; Flow meter calibration

TH814

A

CN11-1780(2020)02-0020-07

2020-03-20

Email:ycyk704@163.com

TEL:010-68382327 010-68382557

程 明 1994年生，碩士研究生，主要研究方向為超聲流量測量和超聲顯微相關(guān)技術(shù)。

馮紅亮 1976年生，高級工程師，研究方向為傳感器及測試技術(shù)。

肖定國 1959年生，副教授，碩士生導(dǎo)師，主要研究方向為超聲掃查檢測與材料特性測試、超聲顯微檢測與測量技術(shù)、超聲檢測儀器與換能器校準技術(shù)。

韓德龍 1994年生，碩士研究生，主要研究方向為測量及傳感技術(shù)。