利用流體中溫度噪聲信號的相關分析進行流速測量.docx

1、第14卷第2期J993年6月資訊htip:/石。哈爾濱船舶工程學院學報JournalofHarbinShipbuildingEngineeringInstituteV6.I4N2Jun.,1993利用流體中L度噪聲信號的相關分析進行流速測量*宿成基防志明”夏虹（船殖動力工程系）（摘要】通過時流體中存在的溫度噪聲信號的相關分析，進行相關技術用于測量氣體流速的實驗研究與理論分析.首先,利用粕關測量系統(tǒng)的敷學模型在PC機上進行了仿真研充，為分析各參斂對相關函數的影響、改進測量方法以及研霓誤差產生的原因，從理論上提供了依據.其次，筒要介紹了所用的實驗裝實洌結果以及進一步提高制量精度的途徑.關鍵詞：相蕓

2、些維;里把述遂，混度噪聲信號，竺Mil儀、中國圖書疝疝27亟蘭一。引言在工業(yè)生產過程中，有些難測流體例如低溫流體、反應堆一回路載熱劑、兩相流體等速度與流量的測量，長時間以來一直是一個亟待解決的問題.近年來，國內外測雖工作者將相關技術用于測蜃流體的速度.根據流體的4$點，采用不同原理的傳感器從流體中獲取流動噪聲信號，經過相關分析，對互相關函數“尋峰求得流體的渡越時間，從而得到速度與流量.相關測延技術不但能夠實現(xiàn)非接觸測量，而且可以測低非穩(wěn)態(tài)流體的瞬時速度與流量，顯示出其獨特的優(yōu)越性.本文以鳳洞中的氣流作為研究對象，采用熱絲作為敏感元件，對從氣流中獲取的溫度噪聲信號進行相關處理，測出氣流的速度與流

3、量.1系統(tǒng)模型與仿真如果將介于兩個傳感器之間的測量段視為線性時不變系統(tǒng)，該系統(tǒng)可以用單位脈沖響應函數g,Q）描述（圖1）系統(tǒng)的輸入信號泌），輸出信號兇）是由流體中的溫度噪聲所引起的.兩個傳感器與相應的濾波.放大器的傳遞函數分別用幻。）思2。）表示.通過它們的作用得到信號x（0,；（0.輸入與輜出信號通過卷積相聯(lián)系.即，r-+叩）«（/）（1）收稿日期:I992TM-I7L兄Q)=J：工(，一T)gI(七)頃+氣(。=x(/)*g(。+氣Q)P(O=J：V(，一T)g2(rt)dr1+七Q)=xQ)g2O)+七Q)利用上式進行互相關函數運算.則可得到相關測量系統(tǒng)數學模型的基本結構=(4

4、)考虛到流動截面上速度分布不均勻，測量的互相關函數由不同速度下的互相關函數R,匕K)與速度分布密度*(K)加權得到：&行6)=U_8侄，(丁)-"(/力訂*W(5)其中.速度分布密度(外以正態(tài)'分布函數計算.由式中可見，互相關函數0r,(t)的計算需要經過兩次卷積.分別由下列參數參加運算：(T)為輸入信號雙f)的自相關函數；為測量段流體的單位脈沖響應函數；為兩個信號傳遞通道(傳感器與濾波放大部分)的窗函數.隨機輸入過程叮以用自相關函數A(t)描述.它可作為觀察流體中噪聲信號時間與空間結構的一個尺度，通過試驗確定出與速】度分布因素有關的表達式：*=t-t2T產試_7式中

5、，君，弓時間常數；1平均速度；K比例系數.考慮到流體流動過程中的速度波動、溫度與濃度波動，可以用一組微分方程表示氣流的流動過程.在一定的條件下，可以推導出帶有通用擴散系數的一般型式：(7)芝+產一。±4=0式中a被傳遞員；D用擴散系數.根據一定的初始條件與邊界條件，可以得到流體的傳遞函數G,(s)與單位脈沖響應函數&Q)：如=斜叩式中L感器間的距離；To均渡越時間.傳感器與濾波放大部分可視為具有帶通泄波器的傳遞特性.它的窗函數可以通過對濾波器獗譜進行富立葉反變換得到K1JiL1i«if=就布定C°0）式中T"l高、低通渡波器的時間常數；K,放大系

6、數以上導出了互相關函數夫“y）的基本結構以及模型函數夫。?。?8白）占以*）.為了討論各參數對互相關函數形狀、極值位置的影響，將該模型在IBMPC機上進行數字商真.由于在時域作卷積計算時間比較長，而且誤差較大，因此式（5）的計算是在頻域完成的，然后通過富立葉反變換，得到一紐時域的互相關函數曲線.圖2示出了不同流速下的互相關函數曲線,顯然，在速度較大時，互相關函數曲線比較尖銳、峰值明顯，這對提高測量精度有利，隨著速度的減小，互相關函數逐漸變得平坦.圖3表示了傳感器間距離改變時互相關函數曲線的變化,由圖3可見，隨著傳感器距離的增大，互相關函數曲線變寬、最大值減小，而且導致測得的渡越時間偏小.因此，

7、在相關測量中，兩個傳感器的距離應盡量小.但是，另一方面應防止距離過小時測距相對誤差的影響，要兼顧上述兩個因素，選取最佳的傳感器距離.圖4示出了當Peclet數變化時，流體的單位脈沖響應函數&0）的變化情況.當Peclct數祓小時，曲線峰值明顯左移，不對稱度增加.可見擴散過程的存在，增大了渡越時間測量中的系統(tǒng)誤差.為此，當擴散系數較大時，應當對流動模型進行修正.圖4不同Pcclct時的互相關畫教曲線0200250300350T（ms）圖6布通濾波頻率對互相關函斂的影響濾波影響的增加，互相關函數的峰值位置向偏小方向移動,同時,曲線出現(xiàn)明顯的不對圖5示出了低通波波器流波頻率對互相關函數的影響

8、.當酒波器截止頻率減小時.互相關函數曲線變寬,峰值右移，這是由于低通濾波的作用.使得示蹤粒子中的高頻分量被削弱.而低頻分最得到加強所致.因而使互相關函數極值向增大渡越時間方向移動.當4>150Hz后，相關函數峰值位置基本保持不變.圖6示出了高通淀波頻率對互相關函數的影響.與低通泌波的作用相反，隨著高通稱.高通濾波對互相關函數的影響，一般不會成為測量的主要何購.它主要起抑制信號直流成分的作用.通?？扇≥^小值.2試驗裝置與測量結果試驗采用兩個熱絲探頭（d=5“m）作為傳感器，將它們固定在標準風洞內,兩個探頭間的距離為L（圖1）當氣流流過熱絲時.由于熱絲很細，氣流溫度的微小波動可在熱絲上感應出

9、隨機流動噪聲信號,熱絲質量小，響度時間為微秒級.可以增大信號帶寬.使得兩信號互相關函數的峰值變得尖銳.有利于提高測量精度.圖8示出了由熱線測得信號的自功率譜密度函數.由圖可見，隨機信號帶寬達2kHz左右.將熱絲采集的隨機信號經過泌波、放大后.送入由高速信號處理器TMS320-C25與IBM-PC微機組成的信號處理系統(tǒng).TMS320-C25是一個32位的單片機，它采用哈佛結構流水線作業(yè)以及半溢出檢測技術，因此完成1024點復數FFT運算僅需14ms.大大加快了隨機信號處理的速度,由IBM-PC微機與高速信號處理器構成主從機系統(tǒng)進行相關函數的運算、尋峰求得渡越時間.試驗中采用丹麥DANTEC公司生

10、產的DISA熱線風速儀作為標準測速儀器，可以用來對相關法測出的流速進行標定利用該試驗裝置及研制的相關計算軟件，設計了一組試驗，分別檢驗當流速變化、改變兩個傳感器間距以及波波頻率等參敷對信號互相關函數的影響，實驗結果證明上述仿真計算的各項結論是正確的.164H1(8、e<圖7試毆裝JL筒圖圖8信號的自功率譜密度函教曲線圖9空氣流堆標定的線在進行流速標定試驗中，調整標準風洞的壓力，可以對各種不同的氣流速度進行測量.通過將相關法測出的塑匕與熱線風速儀測得的流速V進行比較,得出圖9所示的速度標定曲線當氣流速度r>3m/s時湘關測重的相對誤差低于5%,當y<3m/$時，測魚的相對誤差較

11、大.這是因為當流速過低時，兩傳感器的距離£固定，流體流經兩個傳感器之間的渡越時間7。較大，因而流動模型的變化亦比較大，使得兩個信號的相關程度變差，互相關函數峰值明顯降低,相關測量的誤差隨之增大，這一結論與圖2所示仿真計算結果是一致的.按照互相關函數的定義：&)=limj；xQ)yQ+t)dt(11)F-*«積分時間應為無窮大，實際上，只能在有限長的時間內進行運算.因而只能得到互相關函數的估值夫刈().另外對于各種傳感器而言，只能采集到有限帶寬的隨機信號，加之信號中所蔻加無關噪聲的影響，必定使實際上計算得到的互相關西數峰值位置與理論值產生偏差.渡越時間的方差為/(c)

12、=*/(礦7J?w(0)/夫“(0)(12)式中B信號帶寬；T互相關函數積分時間；&(0),&(0)分別為=0時信號xQ)與噪聲信號心)的自相關函數；K常數.由上式可見，渡越時間估值的方差與信號帶寬瓦積分時間以及流動過程的信噪比入(0)/&(0)成反比.必須注意到，方差。2(t)與帶寬3的三次方成反比，這說明信號帶寬的微小變動將會引起方差急劇變化.因此，在選擇傳感器時必須認真考慮其濾波特性.以便盡可能取得寬帶隨機信號.測量系統(tǒng)上、下游信號通道的特性匹配，也會對測量精度造成很大影響.如果兩通道特性不相同，會引起較大的測最誤差，嚴重者其至會使兩個信號變成彼此不相關的信號，無

13、法完成相關測最的任務.由式(10)可知，如果兩個信號通道傳遞特性相同，則可得到一個對于*軸(t=0)對稱的窗函數Fm2(C反之，當兩信號通道特性不同,其窗函數的對稱特性會被破壞，則會影響互相關函數的峰值位置，給測量帶來系統(tǒng)誤差.因此，在相關測鼠中，應當盡松保證協(xié)個信號通道特性的良好匹配，提高測最精度.3結論通過對相關測最系統(tǒng)的仿真研究與對氣體流速的實際測量，可以得出以下結論：(1) 通過對單相液體中的溫度噪聲信號進行相關處理，濯定流體的速度與流總是可行的.(2) 計算機模擬與試驗結果一致證明，系統(tǒng)模型是正確的.數字仿真為在實際測信中確定合理的傳感器距離、帶通濾波頻率等提供了理論依據.(3) 實

14、際測楫結果表明，通過增大信號帶寬，減小信號采樣間隔等措施，可以使相關測最達到較高的精度.命考文獻MeschF.SpeedandFlowMeasurementbyanIntclfigentCorrelationSystem.ISA90,NewOrleans,USA,1990,109KoppermannC.GeschwindigkeitsundDurcflussmessungincinphasiRenStroemunzenmitICreuz-Korrelationsverfahren,FortschrBer.VDI-I.Rdbe8Nr.63,1983:24-27XuChengji.Korrelat

15、iveGcschwindigkeitsmessunganKraftstofTtniteinemInfrarot-Seiuor.ForschungsberichtUni.Kaisruhc1982:18-19BeckMS,DrineJ.ParticleVelocnyandMassMeamrementinPneumaticConveryors.PowderTechnolosy1968(2):335陸志明湘關測*技術應用于單相淹體淹*測的研究：學位論文.哈爾濱船姻工程學院,】990FluidVelocityMeasurementbyUsingtheCorrelationAnalysisoftheTem

16、peratureNoiseSignalsinFluidXuChengjiLuZbisiagXUHoag(Dept,ofMarinePowerEng.)AbstractInthispaperbymeansofcorrelationanalysisofthetemperaturenoiseinfluid,experimentandtheoreticalanalysisareperformedforapplyingthistechnologytothemeasurementofgasflowvelocity.First,themathematicalmodelofcorrelationmeasurementsys>temissimulatedonanIBMPCcomputer.Itprovidestheoreticalevidenceforanalysingtheparameterinfluenceoncrosscorrelationfunction,improv