殷興景,熊玉亭,蘇中地

(1.浙江蒼南儀表廠,浙江蒼南 325800;2.中國計量學(xué)院計量測試工程學(xué)院,杭州 310018)

0 引 言

近年來,很多大管道氣體流量檢測儀表采用的原理是速度面積法,即僅測取管道中一點(diǎn)或多點(diǎn)的流速來推算管道流量,安裝形式為插入式。這類儀表的共同特點(diǎn):結(jié)構(gòu)簡單、安裝維修方便、價格低廉和重復(fù)性好,是工控系統(tǒng)中檢測大管道氣體流量性價比較高的儀表[1]。在管道口徑(大于300 mm)時,均速管流量計壓損小、節(jié)能、結(jié)構(gòu)簡單和拆裝方便等優(yōu)勢更加明顯。

均速管流量計是基于早期皮托管(Pitot)測速原理發(fā)展起來的,它是以測管道中直徑(圓管)或長與寬(矩形管)上幾點(diǎn)的流速來推算流量的一種插入式流量儀表。

在國內(nèi)倡導(dǎo)建設(shè)節(jié)約型社會的前提下,均速管流量計被廣泛運(yùn)用于大口徑、大流量氣體介質(zhì)管道的測量中。西氣東輸工程中,在內(nèi)徑為1 000 mm的干線管道上,選用了50臺均速管流量計,占總量的52%[2]。此外,均速管流量計在國內(nèi)電力、冶金、石化等行業(yè)中,也得到廣泛應(yīng)用。

針對均速管流量計輸出差壓小、精確度低,忽視管內(nèi)徑對精確度的影響等缺點(diǎn)推出的一項專利產(chǎn)品——均速環(huán)流量計。它由雙文丘里管測低壓,提高了輸出差壓,用多根均速管充分反映了管內(nèi)的流速分布等一系列措施,改善均速管的技術(shù)特性,已引起國內(nèi)外廠商及用戶的關(guān)注[3]。

1 方法分析

筆者建議,在可能的情況下,針對不同管道的特定截面使用精度較高的流速計測量其速度分布曲線及一定流速范圍內(nèi)的速度分布曲線的變化趨勢,就可以根據(jù)速度分布公式推出該面的速度分布情況,最后通過速度面積法得出流量。

對于圓管內(nèi)速度分布式的研究,最經(jīng)典且最常用的是Nikaradse于1932年提出的指數(shù)分布模型(如圖1所示)。v是管道截面上一點(diǎn)的流速;r代表管道內(nèi)相對于中心軸線的位移;R代表管道半徑,則v和r在理論上滿足公式:

式中 vmax——管道中心處的最大速度,vmax=——管內(nèi)平均流速,n——指數(shù),隨雷諾數(shù) ReD變化的常數(shù)[4]。

圖1 流速分布圖

該模型所描述的流速分布在管壁附近 (r= R)和管道中心 (r=0～ 0.2R)兩處,無論雷諾數(shù)為多少,均與實際情況有較大偏離。因此,用該模型來擬合實際管道中的流速分布吻合度不夠高。

從邊界層的半經(jīng)驗理論出發(fā),考察圓管中充分發(fā)展的湍流問題。整個邊界層可分為粘性底層、緩沖層和湍流核心層三部分[5],以下三位學(xué)者分別提出了按邊界層分段的速度分布公式:

a)Von Karman表達(dá)式:

b)Taylor表達(dá)式:

c)Spalding表達(dá)式[6]:

式中 v+——無因次流速;y+——無因次距離;常數(shù)B=5.2,K=0.41。

這些經(jīng)典的分段式模型與實驗數(shù)據(jù)吻合得較好,但必須考慮自變量范圍分段計算的過程,增加了實際工程計量的難度。為解決這一問題,國內(nèi)有學(xué)者于20世紀(jì)80年代對Spalding式進(jìn)行了改進(jìn),得到了統(tǒng)一描述整個管道截面上的速度分布模型如下[7-8]:

把式(11)和(12)代入式(10),可得到在不同ReD下之間的關(guān)系,它就是以ReD為參量的的一族曲線,當(dāng)預(yù)先給定一ReD值時,可借助計算機(jī)而得到不同下的值。

文獻(xiàn)[13]中通過對比光滑圓管中指數(shù)公式(1)和統(tǒng)一分布式(10)的解,證明了統(tǒng)一分布式與實際情況更吻合。

在圖1中,選擇一個具有穩(wěn)定流速分布的面,在該面上選擇多條直徑上的測點(diǎn),用皮托管進(jìn)行流速測量,當(dāng)測得足夠多的點(diǎn)后,可獲得該面的速度分布曲線。

然而,由于工業(yè)現(xiàn)場各阻力件,管道粗糙度以及前直管段長度不足等因素對流動的影響,管道截面的流速分布可能有畸變,不再是對稱的。因此,必須有針對性地對需要獲得流量值的管道的某一截面進(jìn)行流速分布測量,獲得流速分布方程,再根據(jù)該流速方程以及管道中某點(diǎn)的坐標(biāo)推算出流速。

2 實驗數(shù)據(jù)與結(jié)果討論

實驗在300 mm管道氣體流量裝置上進(jìn)行,實驗裝置框架如圖2所示,整套裝置長50D,實驗段前后各有20 Hz的直管段長度。此套裝置采用吸氣式,動力源為變頻器(6～50 Hz)控制的風(fēng)機(jī),正常工作時管內(nèi)流速范圍為1.30～11.74 m/s。標(biāo)準(zhǔn)表采用精確度為±1%的TBQ系列氣體渦輪流量計,流量范圍為320～6 500 m3/h,工作溫度為-30～60℃,大氣壓力為86～106 kPa,相對濕度為5%～95%RH。

圖2 實驗裝置框架

由于變頻器的頻率在30 Hz以上時,管道受風(fēng)機(jī)擾動,產(chǎn)生一定振動,影響測量精確度,實驗中選用10 Hz的頻率,實驗溫度為20℃,使用皮托管測量,皮托管半徑為2 mm。實驗中,通過調(diào)節(jié)變頻器頻率,測量管道內(nèi)不同流速下,在圖2所示的測量面上的流速分布。10 Hz頻率下測量面的流速分布圖,實驗數(shù)據(jù)如圖3所示,從圖3中可以看出,該面上是充分發(fā)展的湍流速度分布。分別取直徑上坐標(biāo)為-140,-100,-50,0,50, 100,140七個點(diǎn),利用式(1)和式(10)計算點(diǎn)流速,如圖4所示。

圖3 測量面上的流速分布曲線

圖4 實驗值與擬合值的對比曲線

為準(zhǔn)確計量平均速度,如圖5所示,在測量面的一條直徑上均勻地取10個點(diǎn),r值分別為±15 mm, ±45 mm,±75 mm,±105 mm和±135 mm。采用速度面積法獲得該面的流量值,與式(1)和式(10)的流量對比,并與渦輪流量計測得的數(shù)值進(jìn)行比較,由于渦輪流量計的示值的最后一位會發(fā)生上下跳動,取一定時間內(nèi)的平均值,見表1所列。

表1 各種方法的流量值對比 m3/h

通過計算發(fā)現(xiàn),通過公式擬合的流速與實際流速在大部分區(qū)域基本吻合,且統(tǒng)一分布式(10)較式(1)在管中心部分更接近實際值。

圖5 圓管速度面積法測點(diǎn)位置

3 結(jié)束語

該文通過測量圓管內(nèi)的某個面上的流速分布,獲得流速分布曲線后,用指數(shù)分布方程和統(tǒng)一分布式方程擬合該曲線,發(fā)現(xiàn)擬合曲線上的點(diǎn)流速與實際流速較接近。通過這種方法,可以在之后的流速測量中直接利用該面的流速分布曲線及統(tǒng)一分布式方程推算各點(diǎn)的流速值。

該方法必須對不同的管道分別測取某個面的流速分布,獲得對應(yīng)的流速分布曲線,然后用統(tǒng)一分布式方程進(jìn)行擬合,并對方程做必要的修正。最后達(dá)到通過方程和流速分布曲線直接推算流速的目的。

[1] 毛新業(yè),李偉健.大管道氣體流量檢測儀表[J].電力設(shè)備, 2008,9(12):8-11.

[2] 張東飛,蘇中地,史曉妍,等.不同雷諾數(shù)下均速管流量計流量系數(shù)的確定[J].中國計量學(xué)院學(xué)報,2007,18(1): 18-21.

[3] 毛新業(yè).大管道氣體流量檢測儀表[J].傳感器與儀器儀表, 2006,10(10):56-58.

[4] 孔 瓏.工程流體力學(xué)[M].北京:中國電力出版社,2007: 106-107.

[5] 蘇中地.非牛頓冪律流體圓管湍流的新速度分布式[J].中國計量學(xué)院學(xué)報,1997,(1):111-115.

[6] SPALDING D B.A Single Formula for the Law of the Wall [J].Journal of Applied Mechanics,1961.

[7] 張維一,龔家彪.圓管中充分發(fā)展湍流的速度分布與平均流速位置的研究[J].計量學(xué)報.1983,4(1):1-5.

[8] 薛恨春,龔家彪.充分發(fā)展圓管流速分布統(tǒng)一模型及其平均流速位置的研究[J].計量學(xué)報,1985,6(4):274-278.

[9] 張東飛,蘇中地,何馨雨.一種適用于均速管測量的管內(nèi)統(tǒng)一速度分布式[J].中國計量學(xué)院學(xué)報,2007,3: 200-203.