吳 超，王壽松

基于軌跡圖像的霧滴顆粒運(yùn)動(dòng)速度光電測量方法

吳 超1，王壽松2

（1.江蘇中煙工業(yè)有限責(zé)任公司南京卷煙廠綜合管理處，南京210094；2.南京執(zhí)銳科技實(shí)業(yè)有限公司，南京210019）

為了對雙介質(zhì)噴嘴霧滴顆粒的運(yùn)動(dòng)速度分布進(jìn)行測量，建立了基于波長為532nm連續(xù)激光器和工業(yè)相機(jī)組成的在線測量系統(tǒng)，并獲得了霧滴的軌跡圖像。采用Steger方法對獲得的粒子軌跡圖像進(jìn)行了圖像處理，取得了霧滴顆粒的速度分布的數(shù)據(jù)。結(jié)果表明，所開發(fā)的基于粒子軌跡圖像霧滴速度分布光電測試方法可以滿足在線測量的要求。

激光技術(shù)；軌跡圖像；速度分布；雙介質(zhì)噴嘴；圖像處理；煙草加香

引 言

煙草制絲加工過程中，往往需要使用雙介質(zhì)噴嘴對煙草進(jìn)行加料、增濕和加香處理，因而對雙介質(zhì)的噴嘴的工作特性進(jìn)行研究，可以提高煙草的加料、增濕和加香效果，對節(jié)能增效具有重要的意義。然而，目前該方面還基本停留在工藝參量優(yōu)化研究方面，對直接影響煙草均勻吸收的噴嘴工作狀態(tài)判斷也只停留在經(jīng)驗(yàn)判斷方面。所以，有必要對噴嘴的工作狀態(tài)及效果進(jìn)行定量測試研究，而其中，霧滴顆粒的速度分布是給料效果的重要影響因素之一。

對霧滴顆粒速度的測量有多種方法，如激光多普勒測速法（laser Doppler velocimetry，LDV）［1］、高速攝影［2］、粒子成像測速（particle image velocimetry，PIV）［3-4］、數(shù)字全息技術(shù)［5］等等。然而這些方法都有各自的優(yōu)缺點(diǎn)，如LDV只能進(jìn)行單點(diǎn)測量，而且對測量環(huán)境要求高；PIV設(shè)備和數(shù)字全息技術(shù)設(shè)備復(fù)雜且昂貴；高速攝影的圖像分辨率低、軟件處理速度慢、價(jià)格也比較昂貴［6］。而最近發(fā)展起來的一種基于粒子軌跡圖像的測試方法，具有結(jié)構(gòu)簡單、成本低、可視化測量的特點(diǎn)。目前國內(nèi)外不少學(xué)者嘗試將其用于在線測量。如2007年，ZHANG等人用單幀圖像測量汽輪機(jī)低壓缸內(nèi)二次水滴的速度和流動(dòng)角［7］；2009年，CHINNAYYA在4MW煤粉燃燒器實(shí)驗(yàn)機(jī)組的水平輸粉管道上安裝了圖像探針，并利用該方法測量了顆粒物的速度［8］；2011年，WU等人對煤粉顆粒的粒徑和速度進(jìn)行了測量［9］；2011年，ZHANG等人利用該方法對固相二相流速度場進(jìn)行了測量，重點(diǎn)研究了圖像處理算法［4］。目前該方法應(yīng)用難點(diǎn)在于圖像處理算法尚不夠完善。常見的光條中心線提取方法有灰度閾值法、極值法和梯度閾值法等，這些方法實(shí)現(xiàn)簡單，但精度不高。STEGER利用Hessian矩陣得到圖像中的道路和血管等類似條紋的法線方向，然后求法線方向上的極值點(diǎn)得到條紋中心線的亞像素位置［10-11］。Steger方法具有精度高、魯棒性好等優(yōu)點(diǎn)，廣泛應(yīng)用于航拍圖像和醫(yī)學(xué)圖像的自動(dòng)分割。

本文中首先建立了基于激光片光源的在線測試系統(tǒng)，對煙草行業(yè)所用的雙介質(zhì)噴嘴的速度分布進(jìn)行了測量，獲得了霧滴的軌跡圖像，并采用Steger方法進(jìn)行圖像處理，得到了霧滴顆粒軌的速度分布。

1 測量原理

針對運(yùn)動(dòng)顆粒的測試中，為了防止顆粒圖像的模糊，需要盡量縮短顆粒圖像的曝光時(shí)間。通?？梢圆捎脙煞N方式，一種是通過加快CCD的快門速度；另外一種方式是采用瞬態(tài)探測光源照明。由于現(xiàn)有的工業(yè)相機(jī)曝光時(shí)間通常都在數(shù)微秒量級，當(dāng)顆粒具有較高的速度，或者測量系統(tǒng)具有較大的放大率的時(shí)候，單純利用工業(yè)相機(jī)往往無法避免顆粒圖像的模糊，因而測試高速運(yùn)動(dòng)物體速度常需要瞬態(tài)光源。如PIV方法測速，通常采用調(diào)Q的脈沖激光，即持續(xù)時(shí)間為納秒量級的光脈沖作為瞬態(tài)照明光源。該光源可以克服拍攝時(shí)因顆粒運(yùn)動(dòng)引起圖像的模糊，然而這樣做顯著地增加了設(shè)備成本。

基于顆粒軌跡圖像的方法充分利用相機(jī)長曝光的特點(diǎn)，該方法還要適當(dāng)減緩相機(jī)快門速度，使顆粒獲得較長的曝光時(shí)間。這樣獲得的顆粒圖像實(shí)際上是顆粒運(yùn)動(dòng)的軌跡。對顆粒軌跡圖像進(jìn)行處理，即可獲得顆粒的運(yùn)動(dòng)速度。相應(yīng)的測試系統(tǒng)中采用連續(xù)激光器和工業(yè)相機(jī)各一個(gè)，相對于傳統(tǒng)的測試方法顯著降低了成本。

設(shè)相機(jī)的曝光時(shí)間為te，該曝光時(shí)間可以通過程序設(shè)置后自動(dòng)執(zhí)行。設(shè)顆粒的軌跡長度分別為L，設(shè)鏡頭的成像倍率為β，則顆粒的粒徑d可表示為：

當(dāng)L?d時(shí)，則由上式可得到：

當(dāng)運(yùn)動(dòng)顆粒在光源照明下的軌跡足夠長時(shí)，可以忽略粒徑的影響，并由（2）式得到顆粒的運(yùn)動(dòng)速度。這是本文中涉及的測試系統(tǒng)的測量原理。相機(jī)的曝光時(shí)間te可以通過軟件設(shè)置得到，而鏡頭的成像倍率β也能通過實(shí)驗(yàn)測量得到，于是如何通過圖像處理的方法獲得顆粒軌跡的長度成為該測量方法的關(guān)鍵。

由于存在畸變，因而鏡頭的成像倍率β并不是固定值。因而為了獲得更加準(zhǔn)確的結(jié)果，本文中采用國際象棋盤圖案的標(biāo)定板和開源軟件Open CV中的標(biāo)定算法對該測試系統(tǒng)的成像鏡頭進(jìn)行標(biāo)定，標(biāo)定板尺寸為20mm×20mm。由于標(biāo)定算法比較成熟，這里不再贅述。

另外，相機(jī)的真實(shí)門寬也與標(biāo)稱值有所偏差，因而也需要對相機(jī)的快門的真實(shí)曝光時(shí)間進(jìn)行標(biāo)定。作者采用一個(gè)調(diào)Q的脈沖光源進(jìn)行標(biāo)定，測試原理如圖1所示。

Fig.1 Schematic diagram of camera gate width calibration setup

通過電腦設(shè)置工業(yè)相機(jī)的門寬tw，并通過數(shù)字延時(shí)器DG535改變調(diào)Q脈沖激光和工業(yè)相機(jī)之間的延時(shí)時(shí)間td。測試系統(tǒng)中，程序自動(dòng)以td×1%的步長掃描－2tw＜td＜2tw的延時(shí)區(qū)間（如果td×1%大于相機(jī)最小可設(shè)置的步長，以相機(jī)最小步長掃描），在每步延時(shí)獲取相機(jī)采集到的圖像，并判斷采集圖像的曝光亮度。如果曝光亮度超過設(shè)定的閾值（系統(tǒng)圖像的亮度閾值設(shè)為50，其中相機(jī)的最高亮度為255），則認(rèn)為激光輸出的脈沖光在門寬之內(nèi)，如果不滿足，則說明相機(jī)沒有采集到脈沖激光。通過程序可以得到一組td，求得最大值td，max和最小值td，min，兩者之差即為工業(yè)相機(jī)在標(biāo)稱門寬tw的真實(shí)門寬te。

2 實(shí)驗(yàn)裝置

霧滴顆粒速度測量裝置如圖2所示。該裝置有待測的雙介質(zhì)噴嘴、連續(xù)激光器、工業(yè)相機(jī)和電腦組成。其中連續(xù)激光波長為532nm，輸出功率1.2W，激光片光源厚度約為1mm。此處的片光源不宜過厚，這樣測試區(qū)域的霧滴顆粒密度就不會(huì)過高，否則嚴(yán)重的重疊會(huì)給圖像處理帶來困難；而片光源太薄的話，顆粒易穿出測試區(qū)域，加大了速度測量的誤差。綜合考慮，本文中采用的厚度為1mm，而且使片光源穿過了噴射區(qū)域的中心軸位置，這樣可以基本保證霧滴顆粒能夠被片光源照亮。噴嘴采用空氣引射水，引射的氣壓可以調(diào)節(jié)。

Fig.2 Schematic diagram of test system

測試區(qū)域霧化顆粒的軌跡圖像由數(shù)字相機(jī)實(shí)時(shí)采集，并傳輸給電腦進(jìn)行實(shí)時(shí)處理。相機(jī)的曝光時(shí)間最低為10μs，分辨率為1024×768。在測量中，需根據(jù)霧滴顆粒的軌跡圖像改變相機(jī)的門寬。這是因?yàn)槿绻鄼C(jī)的門寬太寬，得到的霧滴顆粒的軌跡圖像重疊嚴(yán)重，將給圖像處理帶來困難；如果相機(jī)的門寬太窄，一方面會(huì)引起霧滴顆粒圖像的曝光不足，另外一方面，粒子軌跡會(huì)比較短，給測量結(jié)果帶來誤差。

3 測量結(jié)果和軌跡圖像處理

圖3所示為壓縮空氣壓力約為0.25MPa時(shí)的典型圖像，測試區(qū)域距離噴嘴約1m的距離。

Fig.3 Typical image obtained when compressed air pressure is0.25MPa

為了獲得高精度的測試結(jié)果，作者采用Steger方法［10-11］。相對于傳統(tǒng)的灰度閾值法、極值法和梯度閾值法等算法，Steger方法能給出更高精度的測量結(jié)果。該方法利用Hesslan矩陣得到方向，然后求得法線方向上的極值點(diǎn)得到條紋中心線的亞像素位置。Steger方法具有精度高、魯棒性好等優(yōu)點(diǎn)。作者采用與ZHOU等人對結(jié)構(gòu)光相類似的處理方法［12］。具體的操作如下。

首先利用模板高斯卷積計(jì)算霧滴顆粒軌跡圖像上各點(diǎn)的Hessian矩陣確定其法線方向，在法線方向進(jìn)行泰勒級數(shù)展開求得亞像素坐標(biāo)，并將相鄰的中心點(diǎn)連接。如果兩個(gè)霧滴顆粒過近，得到的中心點(diǎn)可能位于兩條直線上，此時(shí)可利用直線的交叉點(diǎn)對直線進(jìn)行拆分。最后對共線中心點(diǎn)進(jìn)行融合，構(gòu)成新的直線，通過長度的約束，得到霧滴軌跡對應(yīng)的中心點(diǎn)，然后進(jìn)行直線擬合，即可得到霧滴顆粒的軌跡圖像。

提取的軌跡如圖4所示。由圖4可以看出，該方法基本能夠準(zhǔn)確獲取圖中大部分的霧滴顆粒軌跡。結(jié)合之前圖像定標(biāo)的結(jié)果，即可得到軌跡在世界坐標(biāo)系中的長度。已知相機(jī)的曝光時(shí)間，通過簡單的計(jì)算就可以得到霧滴顆粒的速度。

Fig.4 Trajectory map extracted with Steger method

基于該算法，本文中開發(fā)了實(shí)時(shí)在線測試的軟件。程序運(yùn)行時(shí)，軟件自動(dòng)采集圖像，并自動(dòng)提取顆粒圖像進(jìn)行分析和統(tǒng)計(jì)，程序設(shè)定統(tǒng)計(jì)數(shù)量為5000個(gè)，當(dāng)統(tǒng)計(jì)數(shù)量超過5000個(gè)時(shí)，程序就會(huì)停止采集圖像，并輸出結(jié)果。

圖5中給出了霧滴顆粒速率分布，從圖中可以看出，其速率主要分布在1m／s～3m／s之間，個(gè)別速率接近5m／s。作者認(rèn)為這是由于壓縮空氣氣壓不穩(wěn)定造成的，也可能是由于霧化不充分，存在的較大的顆粒具有較高的速度，這需要進(jìn)一步的研究。另外，仔細(xì)觀察圖4可以發(fā)現(xiàn)，圖中有兩個(gè)的霧滴顆粒具有較高的速度，這與測試的結(jié)果相符。

Fig.5 Droplet particle velocity distribution

系統(tǒng)的測試誤差主要來源于兩個(gè)方面：一個(gè)方面是成像鏡頭不可避免地存在一些失真，系統(tǒng)采用的鏡頭標(biāo)稱失真率約為0.1%，該失真也會(huì)造成一些誤差；另一個(gè)方面，系統(tǒng)忽略了霧滴粒徑的影響，當(dāng)霧滴顆粒的運(yùn)動(dòng)軌跡足夠長時(shí)，該誤差可忽略不計(jì)。

另外需要指出的是，當(dāng)霧滴的濃度過高時(shí)，將有大量的霧滴軌跡重疊，雖然本文中的算法能夠區(qū)別一定的重疊，但是若發(fā)生大量的重疊，特別是霧滴軌跡首尾相連的情況發(fā)生概率提高時(shí)，就會(huì)給系統(tǒng)帶來顯著的誤差，此時(shí)可通過提高測試系統(tǒng)的放大率，同時(shí)縮短相機(jī)的快門來避免霧滴軌跡圖像的大量重疊，進(jìn)而提高測量誤差。

4 結(jié) 論

針對煙草行業(yè)所用的雙介質(zhì)噴頭噴霧顆粒速度分布測量的問題，本文中開發(fā)了基于軌跡圖像處理的在線測試系統(tǒng)，其測量原理簡單、價(jià)格低廉。對系統(tǒng)獲得的圖像以Steger方法為基礎(chǔ)，實(shí)現(xiàn)了其亞像素精度的提取相關(guān)信息。由于Steger方法具有較好的魯棒性，以該方法為基礎(chǔ)開發(fā)的程序，實(shí)現(xiàn)了軌跡圖像的在線實(shí)時(shí)測量，并獲得了理想的結(jié)果，為煙草行業(yè)所用的雙介質(zhì)噴嘴的霧化狀態(tài)的表征提供了技術(shù)手段。

［1］ TSUJIY，MORIKAWA Y.LDV measurements of an air-solid twophase flow in a horizontal pipe［J］.Journal of Fluid Mechanics，1982，120（1）：385-409.

［2］ WANG D Zh，HUANG Zh，TONG Ch J，et al.Characteristics of diesel sprays by high speed video camera［J］.Journal of Shanghai Jiaotong University，2000，34（4）：453-457（in Chinese）.

［3］ ANGARITA-JAIMESD，TOWERS C E，TOWERS D P.Three component multi-phase velocimetry measurements on a GDI spray using optically efficient fluorescent tracers［J］.Experiments in Fluids，2012，52（4）：949-962.

［4］ ZHANG Y Y，ZHANG GM，XU M.Laser diagnostics for spray of spark ignition direct injection（SIDI）combustion system［J］.Journal of Automotive Safety and Energy，2011，2（4）：294-307（in Chinese）.

［5］ XU Q，CAON，CAO L，et al.Experimental research of electrostatic spray with digital holography［J］.Laser Technology，2013，37（2）：143-146（in Chinese）.

［6］ ZHANG J J，F(xiàn)AN X L，CAIX Sh.Measurement for the velocity vector field of gas-solid two phase flow using single frame and single exposure imaging［J］.Journal of Engineering Thermophysics，2012，33（1）：79-82（in Chinese）.

［7］ ZHANG H，CAIX Sh，SHANG Zh T，et al.The measurement of size，velocity and flow angle of coarse water with xingle image method［J］.Thermal Turbine，2008，37（1）：26-29（in Chinese）.

［8］ CHINNAYYA A，CHTAB A，SHAO JQ，et al.Characterization of pneumatic transportation of pulverised coal in a horizontal pipeline through measurement and computational modelling［J］.Fuel，2009，88（12）：2348-2356.

［9］ WU X Ch，WANG H，HU Q，et al.Velocity and size measurement of coal particle by trajectory imaging［J］.Journal of Zhejiang University（Engineering Science Edition），2011，45（8）：1458-1462（in Chinese）.

［10］ STEGER C.Extraction of curved lines from images［C］／／In 13th International Conference on Pattern Recognition.Vienna，Austria：IEEE，1996：251-255.

［11］ STEGER C.An unbiased detector of curvilinear structures［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，1998，20（2）：113-125.

［12］ ZHOU FQ，WANG F，ZHANG G J.Three-step extraction method for line-stripe in structured light vision［J］.Chinese Journal of Mechanical Engineering，2008，44（11）：215-219（in Chinese）.

Measurement method of the drop let particle velocity based on the trajectory images

WUChao1，WANGShousong2
（1.Integrated Management Office，Nanjing Cigarette Factory，China Tobacco Jiangsu Industrial Co.Ltd.，Nanjing210094，China；2.Nanjing Zhirui Technology Industrial Co.Ltd.，Nanjing 210019，China）

In order to measure the velocity distribution of dual-media nozzle droplet particles，an on-line measurement system was established based on continuous laser at 532nm wavelength and industrial cameras.The droplet trajectory images were obtained.The Steger method was used to process the particle trajectory images.The velocity distribution of droplet particles was extracted.The results show that the test method for measuring the speed distribution of the droplet particles based on particle trajectory images，meets the requirements of the online measurement.

laser technique；track image；velocity distribution；dual media nozzle；image processing；tobacco flavoring

TP391.41

A

10.7510／jgjs.issn.1001-3806.2014.02.015

1001-3806（2014）02-0214-04

吳 超（1972-），男，工程碩士，工程師，主要研究方向?yàn)闄C(jī)械設(shè)計(jì)、電氣控制。

E-mail：wuchao_ny＠126.com

2013-05-20；

2013-07-10