數字式時差法超聲流量計的設計與實現

王飛 官龍騰 張福春 黃躍峰 李勇平

(中國科學院上海應用物理研究所1,上海 201800;中國科學院大學2,北京 100049;中國科學院核輻射與核能技術重點實驗室3,上海 201800)

數字式時差法超聲流量計的設計與實現

王飛1,2官龍騰1,2張福春1,3黃躍峰1,3李勇平1,3

(中國科學院上海應用物理研究所1,上海 201800;中國科學院大學2,北京 100049;中國科學院核輻射與核能技術重點實驗室3,上海 201800)

基于時差法原理,采用數字化方案設計實現了一種超聲波流量計。系統以TMS320F28335數字信號控制器(DSC)為核心,擴展ADS807流水線ADC實現高速數據采集,模擬前端實現超聲發射、接收信號調理。DSC軟件算法通過互相關及擬合算法求解順逆流傳播時間差,希爾伯特變換與小波變換結合確定渡越時間,在此基礎上計算流量。經過流量實測證明,設計的超聲流量計充分發揮了數字化方案的優勢,較好地完成了流量測量工作。

時差法 超聲流量 數字信號控制器 互相關 渡越時間

0 引言

目前,超聲波流量計根據測量原理的不同,大致可分為傳播速度差法(時間差法、相位差法、頻差法)、多普勒法、波束偏移法、噪聲法等類型。其中,時間差法超聲波流量計具有準確度高、性能穩定的特點,被廣泛應用于工業環境下大多數清潔均勻流體流量的測量。

時間差法的基本原理是:超聲換能器順逆流交替發送、接收超聲信號,由于流體的流動超聲信號順流傳播時間會小于其逆流傳播時間,通過順逆流傳播時間及其時間差異可以獲取流速、流量等信息[1-2]。

本文基于時間差法原理,采用數字化方案設計實現了一種超聲波流量計。具體介紹了設計所采用的流量計算模型以及信號處理算法,并以TMS320F28335為核心實現了整個測量過程的控制和信號處理算法。

1 原理介紹

測量原理說明如圖1所示。

圖1 測量原理圖Fig.1 Illustration of the measuring principle

順逆流換能器以Z法對射式外夾在管壁上,路徑1為流體靜止情況超聲傳播路徑,路徑2為順流超聲傳播路徑。超聲傳播路徑上的流體線平均流速為v,超聲在流體中的傳播時間為tfl,于是有波束偏移:

路徑1、2傳播時間差:

式中:cs為聲楔聲速;θ為聲楔角度。

從而得到順逆流傳播時間差:

進而得到超聲傳播路徑上的流體線平均流速:

將v乘以流體力學校正因子Kf,得到沿管道截面的平均流速。最終,流量計算公式為:

式中:A為管道橫截面積。

2 硬件設計

超聲波流量計由超聲波換能器、電子學系統組成,硬件框圖如圖2所示。

圖2 硬件框圖Fig.2 Block diagram of the hardware

超聲波換能器是實現電能和聲能相互轉換的器件,設計采用壓電超聲換能器,中心頻率為2 MHz。直接面向超聲換能器的模擬前端實現順逆流換能器的激勵、接收信號調理。由于實際工況下傳播介質對超聲的衰減不同,為了將接收信號穩定在一定的幅度范圍里,采用閉環增益控制,通過采樣計算調節放大器增益。模擬前端對目標信號應有足夠的放大能力,其對干擾的抑制能力也是一個重要考量,可采用平衡差分發射、平衡差分接收電路。

數字部分核心器件采用TI公司的TMS320F28335 DSC,其具有豐富的片上資源和軟件資源支持,既能完成對測量過程的控制,又能勝任復雜的信號處理運算。由于數字信號控制器(DSC)兼具微控制器(MCU)的功能和數字信號處理器(DSP)的處理能力,將其作為控制和計算的核心,是一種高度集成、高性能、高可靠性的解決方案。

高速數據采集部分采用ADS807,這是一款12 bit、最高采樣率53 MS/s的高速流水線ADC,將其接在DSC的XINTF上,實際工作采樣率25 MS/s,采樣間隔40 ns。數據采集流程:DSC的ePWM1產生25 MHz占空比50%的方波作為ADC的輸入時鐘,產生PWM波的同時產生EPWM1SOCA事件觸發DMA數據傳輸,通過XINTF讀取ADC轉換數據到DSC片內RAM[3-4]。

人機接口部分,通過XINTF擴展TFT LCD,用于設置參數交互、測量參數顯示、接收波形顯示。鍵盤掃描芯片CH452負責鍵盤事務,處理用戶輸入。

3 算法及實現

3.1 窗口技術

將發射時刻、順逆流信號到達時刻在同一根時間軸上用實心圓標示,其示意圖如圖3所示。虛線圓表示不受流體流動影響即流體靜止情況下的超聲信號到達時刻,聲路計算過程中稱之為預計到時tE,Δt為順逆流傳播時間差。

圖3 窗口技術示意圖Fig.3 Schematic diagram of the window technology

超聲信號傳播遵循斯涅耳定律,儀器根據管道參數換能器參數、流體參數、聲程數推導出聲路,給出換能器的安裝距離,計算出預計到時tE。再根據最大流速計算Δt的最大值,即超聲信號可能的最早、最晚到時,進而確定出信號窗口位置,并留有一定裕量,只在窗口位置進行接收信號的采集。窗口技術減小了數據采集處理開銷,避免了窗外信號對測量的影響,還可以檢驗參數設置、換能器安裝是否正確。

3.2 時差算法

將采集的順逆流波形疊加顯示如圖4所示。

圖4 順逆流波形Fig.4 Downstream and upstream waveforms

從圖4可以看出,兩路波形在時間上有先后差異,這種差異和流速成正比。時差算法本質上是為了求解具有相似性的兩路波形在時間上的差異,這里采用互相關法,將順逆流波形作互相關,互相關函數峰值位置所對應的時間位移即順逆流傳播時間差Δt。

3.2.1 互相關

數字域中信號x(n)和y(n)的互相關函數定義為:

當m為正、負值時,互相關函數值分別對應x(n)保持不動,將y(n)左移、右移m個時間間隔后兩個序列對應相乘再相加的結果,反映了不同時延下兩信號的相似程度。實際中數據長度N常為有限值,則求rxy(m)估計值的一種方法是:

直接從時域計算互相關,計算量較大,特別是當N、m比較大時。根據維納-辛欽定理,互相關函數與互功率譜是一對傅里葉變換對,利用FFT從頻域實現互相關快速計算[5],計算步驟如下。

①xN(n)補N個零,yN(n)補N個零,通過FFT得X2N(k)、Y2N(k);

圖5 互相關計算結果示意圖Fig.5 Result of cross-correlation calculation

3.2.2 互相關擬合

數字域中直接通過互相關峰確定時間差的準確度與采樣率有關,采樣率越高,采樣間隔越小,準確度越高,誤差越小。然而過高的硬采樣率會增加系統的開銷,卻不能帶來性能上質的提升。相比之下,可以通過軟件算法上數據插值、數據擬合的方法來逼近[6-7]。本文的做法是,首先對順逆流原始數據插值以提高采樣率,進而計算互相關得到離散互相關函數峰值位置作為粗定位結果。然后以離散互相關函數峰值位置為中心,在其左右各取5個點作最小二乘法二次曲線擬合,計算出二次曲線對稱軸位置,得到其與中心位置的偏差作為細定位結果。最后將粗定位結果與細定位結果相加,達到對互相關函數峰值位置的一個較好逼近,再乘以采樣間隔得到最終的時差計算結果。上述擬合過程如圖6所示。

圖6 互相關擬合過程Fig.6 Cross-correlation and fitting

3.3 渡越時間

順流到時tD,逆流到時tU,順逆流到時取平均得平均傳播時間,對應于流體靜止情況下的傳播路徑。該路徑下的預計到時與實測傳播時間的比值還可以作為一個診斷參數判斷流量計是否工作正常。超聲信號平均傳播時間減去超聲在聲楔中的傳播時間、管壁中的傳播時間,即得到超聲在流體中的傳播時間tfl。

運用小波的奇異點定位功能,采用希爾伯特變換與小波變換結合確定渡越時間。由于超聲接收信號的近似光滑性,直接利用小波定位效果不好,通過希爾伯特變換將信號轉換成可用于小波奇異點定位的信號,再使用小波變換進行定位[8-10]。

對原始信號進行希爾伯特變換得到解析信號,解析信號的模乘以相角得到轉換信號,如圖7所示。轉換信號的突變周期與原始信號周期相同,突變點恰是原始信號的局部極小值點。用小波檢測出突變點,進而確定信號到達時刻。

運用雙正交樣條小波對轉換信號進行小波分解,將尺度為4的小波系數與原始信號疊加,如圖8所示。由圖8可見,小波系數明顯體現出轉變信號突變點的位置,即原始信號局部極小值點位置。找到滿足信號周期約束第一個極小值點,向前移3/4個信號周期在小范圍內通過擬合方式確定出信號過零點作為到達時刻,計算得出渡越時間,進而求得tfl。

圖7 希爾伯特變換轉換信號Fig.7 Hilbert transform conversion signal

圖8 信號到達時刻確定Fig.8 Determination of the arrival time

除了本文介紹的方法,還可以采用互相關、功率統計等多種信號處理方法、在多個域中對信號進行處理確定超聲信號到達時刻[10]。

3.4 實現優化

①實現上使用了來自TMS320C28x FPU Library的支持。該算法庫是TI公司針對C28x+FPU開發的高度優化的函數庫,效率極高。它主要包含FFT等信號處理常用運算,使得程序員能最大限度地發揮C28x+FPU的計算處理能力。

②TMS320F28335具有片內Flash,然而150 MHz主頻時片內Flash指令執行速度僅為90～95 MI/s,相比之下片內RAM的指令執行速度為150 MI/s。于是,上電將關鍵代碼拷貝到DSC片內RAM中,加快執行速度。

4 試驗結果

流量檢定在水流量標定平臺上進行,管道內徑50 mm,試驗介質為清水。檢定過程設置5個流量點,每個流量點重復檢定6次,瞬時體積流量結果如表1所示。

表1 檢定結果Tab.1 Test results

從檢定結果可以看出,方案設計較好地完成了流量測量工作,且測量結果具有較高的準確性、重復性。

5 結束語

本文研究了時間差法超聲流量測量,并采用TMS320F28335實現了一種數字化超聲流量測量方案。根據時差法流量測量的基本原理,推導了所采用的流量計算模型,研究并實現了時間差算法、渡越時間算法等關鍵算法。根據流量測量的需求,對系統的總體架構、軟硬件進行了設計,給出了相應的實現方案。經過系統測試和流量實測證明,方案具有較高可行性,為后續的研究、開發工作奠定了基礎。

[1] 馮若.超聲手冊[M].南京:南京大學出版社,1999.

[2] 同濟大學聲學研究室.超聲工業測量技術[M].上海:上海人民出版社,1977.

[3] 李文.基于DSP的超聲波流量計的設計[D].濟南:山東科技大學,2011.

[4] 張文磊.基于TMS320F28335的超聲波流量計[D].濟南:山東科技大學,2011.

[5] 胡廣書.數字信號處理[M].2版.北京:清華大學出版社,2003. [6]龔艷.基于MSP430F1611超聲互相關流量計方案優化[D].北京:北京大學,2009.

[7] 桂永芳.相關法超聲波流量計二次儀表的研究[D].杭州:浙江大學,2009.

[8] 付光懷,楊惠中.基于雙正交樣條小波的超聲回波信號檢測[J].儀器儀表學報,2008,29(4):309-310.

[9] 呂清華,唐惠強.基于小波的信號到達時間定位方法[J].武漢理工大學學報,2006,30(5):874-875.

[10] 張濤,蒲誠,趙宇洋.傳播時間法超聲流量計信號處理技術進展述評[J].化工自動化及儀表,2009,36(4):2-5.

Design and Implementation of the Digital Transit-time Ultrasonic Flowmeter

Based on principle of propagation time difference,the digital transit time ultrasonic flowmeter is implemented by adopting digitized design strategy.The system is composed of the digital signal controller(DSC)TMS320F28335 as the kernel,and ADS807 pipelined A/D converter as the high speed data acquisition device,the ultrasonic signal transmitting and conditioning of receiving signal are performed by analog front end(AFE).The difference of propagation time of the ultrasonic signal upon upstream and downstream in the fluid is solving through cross correlation and fitting algorithms in the software,the combination of Hilbert transform and wavelet transform is determining the propagation time,then the flow rate is calculated on this basis.The practical tests verify that this ultrasonic flowmeter gives full play of the superiority of digitized strategy,and well accomplishes the job for flow measurement.

Transit-time method Ultrasonic flow Digital signal controller Cross-correlation Propagation time

TH814

A

中國科學院戰略性先導科技專項基金資助項目(編號:XDA02010300)。

修改稿收到日期:2014-03-04。

王飛(1989-),男,現為中國科學院上海應用物理研究所信號與信息處理專業在讀碩士研究生;主要從事超聲流量測量的研究。