孫志強(qiáng) ，陳延平

(1. 中南大學(xué) 能源科學(xué)與工程學(xué)院，湖南 長沙，410083；2. 中南大學(xué) 流程工業(yè)節(jié)能湖南省重點實驗室，湖南 長沙，410083)

在能源、動力、石油、化工等流程工業(yè)中，流量是與生產(chǎn)效率和操作安全密切相關(guān)的最為重要的過程參數(shù)之一。作為一類精度較高的通用流量儀表[1]，渦街流量計近年來得到了迅速發(fā)展[2]，已被廣泛應(yīng)用于液體、氣體、蒸汽和部分混相流的測量[3-4]。渦街流量計的原理是利用垂直插入流體中的非流線型旋渦發(fā)生體產(chǎn)生旋渦脫落形成卡門渦街，當(dāng)滿足一定雷諾數(shù)條件時，旋渦脫落頻率與流速成正比，此時采用適當(dāng)?shù)臋z測技術(shù)獲取頻率，即可實現(xiàn)流量測量。在實際使用過程中，渦街流量計不可避免地會受到與之相連接的管道和設(shè)備振動的影響，致使其輸出信號中含有大量噪聲，嚴(yán)重時甚至導(dǎo)致信號完全失真[5]。如何從含有噪聲的信號中準(zhǔn)確提取出旋渦脫落頻率，日益成為優(yōu)化和改進(jìn)渦街流量計測量性能的核心問題之一[6]。綜合實際可行性與應(yīng)用成本，近年來數(shù)字信號處理方法被廣泛應(yīng)用于渦街流量計信號的去噪和旋渦脫落頻率的估計[7]。Amadi等[8]在對原始測量信號進(jìn)行譜估計的基礎(chǔ)上，建立了渦街流量計信號表達(dá)式。Hondoh等[9]利用頻譜估計技術(shù)對渦街流量計輸出信號進(jìn)行濾噪，以提高其測量精度。Xu等[10]通過嵌入功率譜分析技術(shù)，開發(fā)了具有實用價值的渦街流量計信號處理系統(tǒng)。孫志強(qiáng)等通過對管道內(nèi)渦街流動穩(wěn)定性[11]和流場結(jié)構(gòu)[12]的分析，利用Hilbert-Huang變換對渦街流量計信號進(jìn)行去噪[13]，探討了渦街流量計信號的高階統(tǒng)計量特征[14]、功率譜式表征[15]與測量介質(zhì)類型和流量的關(guān)系。這些研究為渦街流量計旋渦脫落頻率的準(zhǔn)確提取奠定了基礎(chǔ)。然而，各種信號處理方法由于受到本身理論的限制，只有在特定參數(shù)條件下才能得到滿意的頻率提取結(jié)果，其通用性尚待驗證與提高。因此，尋找適合于渦街流量計信號分析的方法，研究其性能參數(shù)的優(yōu)化設(shè)置，對旋渦脫落頻率的準(zhǔn)確提取與渦街流量計測量精度的提高都具有重要意義。

1 基于Burg算法的AR模型譜估計

現(xiàn)代譜估計以隨機(jī)過程的參數(shù)模型為基礎(chǔ)，能夠基于有限數(shù)據(jù)從頻域內(nèi)提取出淹沒在噪聲中的有用信號。AR模型(自回歸模型)譜估計是最為常用的一種用于信號特征分析的譜分析方法，其模型參數(shù)的精確估計只需求解一組線性方程即可獲得。Burg算法是建立在數(shù)據(jù)基礎(chǔ)上的AR模型系數(shù)求解的有效算法，已在通信濾波、圖像處理等領(lǐng)域得到成功應(yīng)用[16]。AR模型的Burg算法遞推過程建立在數(shù)據(jù)序列基礎(chǔ)上，避開了序列的自相關(guān)函數(shù)估計，所以，與自相關(guān)法相比，具有較好的頻率分辨率[17]。

利用格型濾波器結(jié)構(gòu)，計算各階前向和后向預(yù)測誤差為：

式中：λk為格型濾波器的反射系數(shù)。

前向和后向預(yù)測誤差的平均功率為：

為使Pk最小，令?Pk/?λk=0，得到反射系數(shù)為：

然后利用Levinson遞推公式求出AR模型的參數(shù)為：

在Burg算法中，階次直接關(guān)系到AR模型譜估計的質(zhì)量。階次過低，譜峰難以出現(xiàn)；階次過高，可能產(chǎn)生虛假峰值。因此，階次的合理選擇是Burg算法能否準(zhǔn)確提取渦街流量計旋渦脫落頻率的關(guān)鍵。

2 實驗

實驗在自行搭建的實驗臺上進(jìn)行。實驗系統(tǒng)通過管道將蓄水池、水泵、穩(wěn)壓水罐、電磁流量計、渦街流量計連接成閉合回路，如圖1所示。測量介質(zhì)為室溫下的水，首先由水泵將其從蓄水池中抽入穩(wěn)壓水罐中，經(jīng)消除流動振蕩后，依次流過提供參考流量讀數(shù)的電磁流量計(精度為0.5級)和實驗用渦街流量計，最后又流回蓄水池循環(huán)使用。實驗的水流量范圍為3～35 m3/h，其大小通過管路中的手動閥進(jìn)行調(diào)節(jié)。

圖1 實驗系統(tǒng)Fig.1 Experimental system

實驗用渦街流量計的儀表系數(shù)為9 217.08 m-3，測量管內(nèi)徑為50 mm，旋渦發(fā)生體為三角形柱體，迎流面寬度為14 mm。渦街流量計的旋渦脫落信號通過軸向管壁差壓法[20]獲取，上游和下游的取壓孔分別位于旋渦發(fā)生體迎流面之前50 mm和之后10 mm。采用的差壓傳感器的動態(tài)響應(yīng)時間為1 ms，其輸出的電壓信號經(jīng)放大后由 Tektronix TDS3054B示波器采集和存儲，采樣頻率設(shè)置為1 kHz，每次測量的采樣時間為10 s。

總之，在概念教學(xué)中，教學(xué)方法有許多，合情推理作為促進(jìn)學(xué)生概念學(xué)習(xí)的有效途徑之一，教師應(yīng)把握好課堂教學(xué)的時機(jī)，找準(zhǔn)合情推理的方法，從而使概念學(xué)習(xí)過程更容易被學(xué)生理解與接受，進(jìn)而全面提升概念教學(xué)質(zhì)量。

圖2所示為流量為11 m3/h時原始的渦街流量計輸出信號片段。圖 2(a)所示為流動平穩(wěn)且無管道振動時的管壁差壓信號，可以看到管壁差壓做頻率恒定的正弦波動，信號幅度變化較小，反映了渦街流量計內(nèi)形成了規(guī)則的卡門渦街。圖2(b)所示為來流振蕩時的管壁差壓信號，此時管壁差壓整體上出現(xiàn)低頻波動，信號幅度差異顯著，表明渦街流量計中雖然存在旋渦脫落，但其穩(wěn)定性受到影響，必須對原始信號進(jìn)行適當(dāng)處理方能有效提取旋渦脫落頻率。

圖2 渦街流量計輸出信號Fig.2 Output signals of vortex flowmeter

3 結(jié)果與分析

對水流量為3～35 m3/h范圍內(nèi)的19組渦街流量計輸出信號，采用基于Burg算法的100階AR模型進(jìn)行譜估計，提取的旋渦脫落頻率結(jié)果如表1所示。其中，參考頻率通過渦街流量計儀表系數(shù)計算得到，相對誤差的計算以參考頻率作為基準(zhǔn)。從表1可知：基于Burg算法的AR模型譜估計對于不同流量下的渦街流量計輸出信號都能準(zhǔn)確提取出其中旋渦脫落頻率。除了在個別流量(12.94 m3/h)下相對誤差稍大外，其余流量下旋渦脫落估計頻率的相對誤差均在±3%以內(nèi)。

表1 旋渦脫落頻率估計性能Table 1 Estimation results of vortex shedding frequency

由于AR模型的階次直接關(guān)系到Burg算法的譜估計精度和計算效率，因此，討論階次對渦街流量計旋渦脫落頻率提取性能的影響具有重要價值。為此，從1至100逐漸增加AR模型的階次，得到Burg算法在不同階次時的旋渦脫落頻率估計值，如圖3所示，其中虛線代表旋渦脫落參考頻率。

圖3 AR模型階次與旋渦脫落頻率的關(guān)系Fig.3 Relationship between AR model’s order and vortex shedding frequency

圖3表明AR模型階次對基于Burg算法的旋渦脫落頻率估計性能具有顯著影響。當(dāng)階次較小時，頻率估計值誤差較大，特別是在旋渦脫落頻率較低時尤為突出。在流量為3.36 m3/h、旋渦脫落頻率為8.60 Hz時，AR模型階次須大于32，估計頻率的相對誤差才能控制在±5%以內(nèi)。進(jìn)一步的試算結(jié)果表明，當(dāng)旋渦脫落頻率小于6 Hz時，基于Burg算法的AR模型譜估計對旋渦脫落頻率的提取將失效，此時譜估計曲線出現(xiàn)劇烈波動。隨著階次的提高，Burg算法的頻率估計性能逐漸穩(wěn)定，與參考頻率趨于一致，旋渦脫落頻率的譜估計誤差最終恒定在一個較低水平。隨著頻率的增大，譜估計所需的最低AR模型階次不斷降低。

利用 Burg算法實現(xiàn)高效高精度提取旋渦脫落頻率的關(guān)鍵在于選擇合適的AR模型階次。對于不同頻率的渦街流量計輸出信號，應(yīng)選用匹配的AR模型階次。由圖 3可知：Burg算法達(dá)到較高精度時的最低AR模型階次隨旋渦脫落頻率的降低而呈現(xiàn)遞增趨勢。

在圖3計算結(jié)果的基礎(chǔ)上，定義AR模型最小階次，即在該階次下，基于Burg算法的旋渦脫落估計頻率與參考頻率的相對誤差在±3%以內(nèi)。圖 4所示為AR模型最小階次與旋渦脫落頻率之間的關(guān)系?？梢姡弘S著旋渦脫落頻率增加，AR模型最小階次降低。

圖4 AR模型最小階次與旋渦脫落頻率的關(guān)系Fig.4 Relationship between minimum AR model’s order and vortex shedding frequency

以正向取整ceil函數(shù)來描述AR模型最小階次P與旋渦脫落頻率f的關(guān)系，即

經(jīng)對圖4中的數(shù)據(jù)進(jìn)行擬合后，得到

以上擬合的殘差平方和為0.93，擬合系數(shù)為0.99，表明該曲線具有較高的擬合度。在采用基于Burg算法的AR模型譜估計時，根據(jù)式(8)可以先設(shè)定一個AR模型階次的參考值，以保證旋渦脫落頻率提取的精度。不建議采用低于參考值的階次，否則可能導(dǎo)致譜估計誤差過大。

選取相同實驗條件下不同流量的渦街流量計輸出信號驗證式(8)的有效性。先按照式(8)計算出階次的參考值，然后利用Burg算法計算該階次下AR模型參數(shù)，估計旋渦脫落頻率，結(jié)果列于表2。其中，參考頻率為通過渦街流量計儀表系數(shù)計算得到的頻率，相對誤差的計算均以參考頻率作為基準(zhǔn)。

表2 最小AR模型階次Burg估計與FFT結(jié)果Table 2 Results of Burg estimation under minimum AR model’s order and FFT

表 2同時還給出了采用快速傅里葉變換(fast Fourier transform，F(xiàn)FT)計算得到旋渦脫落頻率及其相對誤差?？梢姡趨⒖茧A次下Burg算法頻率估計的相對誤差小于 6%，每個流量下的頻率估計性能均優(yōu)于FFT的計算結(jié)果，表明該方法對于渦街流量計旋渦脫落頻率的提取精度較高，在工業(yè)上具有較好的實用性。

4 結(jié)論

(1) 基于Burg算法的AR模型譜估計適合于渦街流量計旋渦脫落頻率的提取，頻率估計值具有較高的精度，滿足工程應(yīng)用的要求。

(2) AR模型的階次對Burg算法的譜估計精度和計算效率具有重要的影響，對于不同頻率的渦街流量計輸出信號，應(yīng)選用匹配的AR模型階次。

(3) 建立了頻率相對誤差小于3%的AR模型最小階次與旋渦脫落頻率的擬合關(guān)系式，AR模型最小階次隨旋渦脫落頻率的增大而減小。

[1] Pankanin G L. The vortex flowmeter: various methods of investigating phenomena[J]. Measurement Science and Technology, 2005, 16(3): R1-R16.

[2] SUN Zhiqiang. Design and performance of the converging-diverging vortex flowmeter[J]. Metrology and Measurement Systems, 2011, 18(1): 129-136.

[3] SUN Zhiqiang. Mass flow measurement of gas-liquid bubble flow with the combined use of a Venturi tube and a vortex flowmeter[J]. Measurement Science and Technology, 2010,21(5): 055403.

[4] SUN Zhiqiang, ZHANG Hongjian. Measurement of the flow rate and volume void fraction of gas-liquid bubble flow using a vortex flow meter[J]. Chemical Engineering Communications,2010, 197(2): 145-157.

[5] Miau J J, Hu C C, and Chou J H. Response of a vortex flowmeter to impulsive vibrations[J]. Flow Measurement and Instrumentation, 2000, 11(1): 41-49.

[6] 黃詠梅, 張宏建, 孫志強(qiáng). 渦街流量計的研究[J]. 傳感技術(shù)學(xué)報, 2006, 19(3): 776-782.HUANG Yongmei, ZHANG Hongjian, SUN Zhiqiang. The development of vortex flowmeter[J]. Chinese Journal of Sensors and Actuators, 2006, 19(3): 776-782.

[7] Zhang T, Sun H, Wu P. Wavelet denoising applied to vortex flowmeters[J]. Flow Measurement and Instrumentation, 2004,15(5/6): 325-329.

[8] Amadi E, Zhu H, Higham E H. Analysis of signals from vortex flowmeter[J]. Flow Measurement and Instrumentation, 1993,4(4): 225-231.

[9] Hondoh M, Wada M, Andoh T, et al. A vortex flowmeter with spectral analysis signal processing[C]//Proceedings of Sicon'01 Sensors for Industry Conference. Rosemont: IEEE Press, 2001:35-40.

[10] Xu K, Huang Y, Lü X. Power-spectrum-analysis-based signal processing system of vortex flowmeters[J]. IEEE Transactions on Instrumentation and Measurement, 2006, 55(3): 1006-1011.

[11] 孫志強(qiáng), 周宏亮, 植曉琴, 等. 管內(nèi)氣液兩相流渦街穩(wěn)定性分析[J]. 中南大學(xué)學(xué)報: 自然科學(xué)版, 2008, 39(6): 1164-1169.SUN Zhiqiang, ZHOU Hongliang, ZHI Xiaoqin, et al. Analysis of the stability of gas–liquid two-phase vortex street in conduits[J]. Journal of Central South University: Science and Technology, 2008, 39(6): 1164-1169.

[12] 桑文慧, 孫志強(qiáng), 周孑民. 有限流道內(nèi)低雷諾數(shù)二維圓柱繞流數(shù)值模擬[J]. 中南大學(xué)學(xué)報: 自然科學(xué)版, 2012, 43(3):1166-1170.SANG Wenhui, SUN Zhiqiang, ZHOU Jiemin. Numerical simulation of two-dimensional flow around a circular cylinder at low Reynolds numbers in finite channel[J]. Journal of Central South University: Science and Technology, 2012, 43(3):1166-1170.

[13] SUN Zhiqiang, ZHOU Jiemin, ZHOU Ping. Application of Hilbert-Huang transform to denoising in vortex flowmeter[J].Journal of Central South University of Technology, 2006, 13(5):501-505.

[14] 孫志強(qiáng), 張宏建. 渦街流量計信號的高階統(tǒng)計量特征[J]. 傳感技術(shù)學(xué)報, 2007, 20(7): 1542-1545.SUN Zhiqiang, ZHANG Hongjian. Higher-order statistics characteristics of vortex flowmeter signals[J]. Chinese Journal of Sensors and Actuators, 2007, 20(7): 1542-1545.

[15] 孫志強(qiáng), 張宏建. 渦街流量計信號能量的功率譜式表征與應(yīng)用[J]. 傳感技術(shù)學(xué)報, 2007, 20(8): 1800-1804.SUN Zhiqiang, ZHANG Hongjian. Power spectrum based representation of vortex flowmeter signal energy and its application[J]. Chinese Journal of Sensors and Actuators, 2007,20(8): 1800-1804.

[16] 韓春林, 何子述. 基于 Burg疊代的多徑信號DOA估計算法[J]. 通信學(xué)報, 2003, 24(7): 24-29.HAN Chunlin, HE zishu. A DOA estimation algorithm for multipath signals based on Burg iterations[J]. Journal of China Institute of Communications, 2003, 24(7): 24-29.

[17] 姚文俊. 自相關(guān)法和Burg法在AR模型功率譜估計中的仿真研究[J]. 計算機(jī)與數(shù)字工程, 2007, 35(10): 32-34.YAO Wenjun. Research on AR model power spectrum estimation based on the algorithm and Burg algorithm[J].Computer and Digital Engineering, 2007, 35(10): 32-34.

[18] ZHANG Hongjian, HUANG Yongmei, SUN Zhiqiang. A study of mass flow rate measurement based on the vortex shedding principle[J]. Flow Measurement and Instrumentation, 2006,17(1): 29-38.