基于圖像分析的底吹攪拌反應(yīng)器混合特性

肖清泰，王仕博，李鵬，高勤，徐建新，王華

（昆明理工大學(xué)省部共建復(fù)雜有色金屬資源清潔利用國家重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650093）

基于圖像分析的底吹攪拌反應(yīng)器混合特性

肖清泰，王仕博，李鵬，高勤，徐建新，王華

（昆明理工大學(xué)省部共建復(fù)雜有色金屬資源清潔利用國家重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650093）

以富氧底吹熔池熔煉過程為研究對象，基于氣泡擾動形成流場圖像測算灰度級強(qiáng)度，測算與比較底吹攪拌反應(yīng)器內(nèi)的局部混合特性。參照某企業(yè)底吹煉銅爐實(shí)際尺寸，搭建氣體噴吹水模型試驗(yàn)平臺；利用高速攝像機(jī)拍攝不同設(shè)計(jì)參數(shù)下的試驗(yàn)過程，捕獲并存儲水模型中噴口上方流場的彩色圖像；截取未受裝置干擾的研究區(qū)域，測算其綠色成分的灰度級強(qiáng)度均值與標(biāo)準(zhǔn)差。研究發(fā)現(xiàn)：氣泡擾動流場RGB圖像的綠色成分在本試驗(yàn)條件下對于真實(shí)圖像氣泡輪廓特征提取更為準(zhǔn)確；綠色成分圖像的水平剖線針對單排噴口與雙排噴口分別呈現(xiàn)單峰與雙峰，驗(yàn)證了以灰度級強(qiáng)度表征混合特性的準(zhǔn)確性；測算研究區(qū)域的灰度級強(qiáng)度均值時(shí)間序列和標(biāo)準(zhǔn)差時(shí)間序列，最終發(fā)現(xiàn)：在熔池動態(tài)平衡階段，大部分時(shí)間序列在統(tǒng)計(jì)學(xué)意義上服從正態(tài)分布。本研究為推進(jìn)圖像分析技術(shù)在底吹攪拌反應(yīng)器流動與混合特性研究方面提供了參考依據(jù)。

氣液兩相流；分布；攪拌容器；混合；成像

冶金爐窯中金屬熔體、煙氣等高溫流體的多相流動參數(shù)測試關(guān)系到冶金過程反應(yīng)的效率與冶金爐窯的安全使用壽命[1]。特別地，富氧底吹銅熔池熔煉過程即呈現(xiàn)流體混合、傳熱傳質(zhì)、動量傳遞及化學(xué)反應(yīng)等復(fù)雜的多相流動現(xiàn)象[2]。金屬熔體混合質(zhì)量的優(yōu)劣直接影響傳熱傳質(zhì)過程，因此研究銅熔體在底吹銅熔煉爐內(nèi)的混合過程及影響參數(shù)十分必要。

水模型實(shí)驗(yàn)法對熔池熔煉爐內(nèi)各類參數(shù)測定提供了便利手段[3]。在研究方法方面，眾多國內(nèi)外專家學(xué)者采用水模型實(shí)驗(yàn)方法或結(jié)合數(shù)值模擬方法，對底吹銅熔煉爐水模型展開系列研究。例如，閆紅杰等[4]采用數(shù)值模擬方法指出冶金過程高溫熔體多相流流動特性模擬效果最佳的數(shù)學(xué)模型是Realizablek-ε模型，并基于此探討了氧槍傾角及直徑對于底吹熔池?cái)嚢栊Ч挠绊憴C(jī)制；張振揚(yáng)等[5]采用數(shù)值模擬方法得到了氧槍直徑、氧槍間距及氧槍傾角等氧槍結(jié)構(gòu)參數(shù)的最優(yōu)組合；邵品等[6]采用Eulerian-Eulerian模型計(jì)算與比較了雙噴嘴在不同噴氣角度下底吹冰銅吹煉爐內(nèi)的氣-液兩相流行為及氣體含量；王東興等[7]基于底吹煉銅轉(zhuǎn)爐的水模型實(shí)驗(yàn)，揭示了底吹造锍轉(zhuǎn)爐內(nèi)噴嘴數(shù)量、噴嘴角度、噴嘴直徑、氣流速度等對轉(zhuǎn)爐熔池內(nèi)氣泡大小分布、氣含率（gas hold-up）及液面噴濺面積的影響規(guī)律；余躍等[8]利用流體體積（volume of fluid，VOF）模型描述底吹爐水力模型實(shí)驗(yàn)，對比研究了4種噴口結(jié)構(gòu)對爐內(nèi)流動過程的影響。在富氧底吹煉銅工藝特性與優(yōu)化方面，郭學(xué)益研究組開展了許多深入系列的基礎(chǔ)研究，進(jìn)行了底吹銅熔池熔煉機(jī)理分析[9]、揭示了底吹爐內(nèi)氧勢及硫勢的梯度波動機(jī)制[10]、研究了渣-銅锍間多組元含量的映射關(guān)系及熔煉過程優(yōu)化[11]。上述研究對優(yōu)化底吹熔煉爐的結(jié)構(gòu)設(shè)計(jì)與生產(chǎn)過程控制具有參考價(jià)值。

在常見的氣液兩相混合系統(tǒng)中，局部氣泡大小分布（bubble size distribution，BSD）是刻畫氣泡形態(tài)學(xué)特征的重要參數(shù)[8]。LI等[12]與YANG等[13]利用數(shù)值模擬方法對高壓環(huán)境下氣液兩相流中單個氣泡的形成、運(yùn)動及形狀等特征參數(shù)進(jìn)行了研究。近年來，數(shù)字圖像處理技術(shù)作為一種常見的現(xiàn)代化信息處理技術(shù)，依其非接觸、可視化的優(yōu)勢在多相流參數(shù)檢測方面的應(yīng)用引起了中外學(xué)者的廣泛關(guān)注。LAU等[14]利用非接觸式的圖像分析技術(shù)研究了泡狀流中BSD，提出了一種表征稠密泡狀流的圖像測量方法。KARN等[15]提出了一種多水平圖像分析技術(shù)來檢測較大尺寸的氣泡，并解決了通風(fēng)式水翼中湍流尾流的氣泡聚類問題。XIAO等[16-17]結(jié)合圖像分析技術(shù)與統(tǒng)計(jì)學(xué)線性回歸模型、假設(shè)檢驗(yàn)理論對有機(jī)工質(zhì)直接接觸式沸騰換熱過程混合特性進(jìn)行了細(xì)致研究，提出了量化混合均勻性與不同混合瞬態(tài)有效識別的技術(shù)手段。上述研究所采用的高速攝影法結(jié)合圖像分析手段對于冶金反應(yīng)過程多相攪拌混合體系混合流型快速準(zhǔn)確地非接觸式評定具有重要參考價(jià)值。

基于上述，數(shù)字圖像處理技術(shù)作為富氧底吹銅熔池熔煉過程水模型實(shí)驗(yàn)的研究手段尚未充分報(bào)道，尤其是在融合數(shù)理統(tǒng)計(jì)理論“強(qiáng)大數(shù)據(jù)建模能力”的技術(shù)研發(fā)方面仍具有更深層次的研究空間，亟待豐富與發(fā)展。鑒于此，本研究通過參考實(shí)際底吹煉銅爐結(jié)構(gòu)參數(shù)制作水模型試驗(yàn)臺，并基于圖像識別技術(shù)深入挖掘底吹攪拌反應(yīng)器混合特性。

1 實(shí)驗(yàn)部分

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

基于水模型原理，考慮本實(shí)驗(yàn)原型與模型主要尺寸的幾何相似開展實(shí)驗(yàn)研究。參照某煉銅企業(yè)底吹熔煉爐實(shí)際尺寸，制作底吹爐的水模型，如圖1所示。其中，冷態(tài)攪拌系統(tǒng)忽略了實(shí)際底吹熔池熔煉過程的傳熱單元，攪拌反應(yīng)器主體采用透明材質(zhì)制作，因此有利于利用可視化手段對氣液兩相在其中的混合特性開展針對性研究。

圖1 底吹熔池熔煉過程水模型試驗(yàn)臺

本試驗(yàn)系統(tǒng)主要由底吹攪拌單元（透明有機(jī)玻璃，60L）、圖像采集單元（高速攝像機(jī)，AVI格式，50幀/秒）以及數(shù)據(jù)處理單元（PC，Intel(R) Core(TM)i5-6300HQ CPU @ 2.30GHz ）等三大單元組成，如圖1所示。其中，玻璃轉(zhuǎn)子流量計(jì)規(guī)格為16～160L/h，噴口直徑為4mm，攪拌槽內(nèi)液位為180mm，噴口氣體流量保持45L/h不變。三腳架的主要作用是穩(wěn)定高速攝像機(jī)，使得高速攝像機(jī)與底吹攪拌反應(yīng)器間距離以及高速攝像機(jī)高度保持不變，保證試驗(yàn)拍攝的視頻基本一致?？諝鈮嚎s機(jī)（型號Z-0.12/8，功率2.0KW）的主要作用是向爐內(nèi)輸送氣體。在試驗(yàn)過程中，利用高速攝像機(jī)記錄試驗(yàn)的全過程，隨后借助于視頻處理軟件保存為bmp格式圖片。氣泡圖像處理過程均在MATLAB環(huán)境下進(jìn)行，通過內(nèi)置函數(shù)、自行設(shè)計(jì)程序完成本文計(jì)算。本試驗(yàn)工況設(shè)計(jì)具體參數(shù)如表1所示。其中，噴口傾角具有0°、7°與14°三種情形，除C4為雙排噴口外，其余噴口均為單排設(shè)計(jì)。

表1 水模型試驗(yàn)設(shè)計(jì)中的因素與水平

1.2 圖像預(yù)處理

氣體對熔池的攪拌存在一個動態(tài)平衡過程，單純分析攪拌初期不夠準(zhǔn)確，下文考慮熔池動態(tài)平衡階段氣泡擾動流場的混合瞬態(tài)情形，圖2即給出了工況C1第250幀時(shí)動態(tài)平衡流場圖像的處理過程。其中：圖2(a)為實(shí)驗(yàn)獲取的原始RGB（彩色）氣泡圖像，噴口上方的氣柱震蕩區(qū)域和水面上方的氣泡破碎區(qū)域可明顯識別。圖2(b)是依據(jù)如下加權(quán)平均公式獲取的氣泡灰度級圖像，見式(1)。

式中，C(i,j)為RGB圖像中第i行、第j列的灰度級強(qiáng)度（grayscale intensity）；R(i,j)、G(i,j)及B(i,j) 分別為RGB圖像中第i行、第j列的紅色、綠色及藍(lán)色分量。圖2(c)為依據(jù)Otsu法（Otsu，1979年）進(jìn)行灰度圖像二值化（binarization）操作之后得到的氣泡邊界輪廓（boundary contour）圖，由此可知噴口上方附近的微小氣泡被誤消除，因此氣泡圖像灰度化操作并未有效提取原始RGB圖像中的明顯特征。圖2(d)、圖2(e)、圖2(f)分別展示了圖2(a)的紅色、藍(lán)色及綠色成分，圖2(g)、圖2(h)及圖2(i)分別對應(yīng)上圖的氣泡邊界輪廓圖。由此可知，底吹攪拌反應(yīng)過程中氣泡RGB圖像的紅色成分較好地提取了氣泡擾動流場的上半部分特征，其藍(lán)色成分較好地提取了氣泡擾動流場的下半部分特征，而只有其綠色成分較為完整地提取了氣泡擾動的整個局部混合區(qū)域的流場特征。

1.3 實(shí)驗(yàn)方法

富氧氣體由底吹熔煉爐底部的噴口高速噴吹進(jìn)入熔池內(nèi)部產(chǎn)生氣泡并形成氣柱，不僅起到氧化劑作用促進(jìn)反應(yīng)進(jìn)行，而且提供了金屬熔體流動的動力來源。為首先從直觀上直接觀察氣泡形成并長大的過程，特選取圖像幀數(shù)分別為f=9，f=11，f=13，f=15，f=17，f=19，f=21，f=23及f=25等時(shí)刻的流場瞬態(tài)彩色圖像，如圖3(a)～圖3(i)所示。由此可知，在氣體噴吹的初始階段，氣泡逐漸形成并長大；當(dāng)氣泡長大到一定程度后，會與噴口端面發(fā)生些許分離；隨著氣量持續(xù)地鼓入，大氣泡脫離于后續(xù)氣泡，后續(xù)氣泡將繼續(xù)生長為大氣泡，因此氣泡攪拌形成的動態(tài)流場實(shí)現(xiàn)了爐內(nèi)多相流體的接觸混合。

人們普遍采用化學(xué)脫色法進(jìn)行流體混合質(zhì)量評價(jià)，但每次試驗(yàn)需依賴觀察者主觀判斷，且只能獲取混合時(shí)間（mixing time）參數(shù)，無法知曉整個混合過程中的瞬態(tài)混合狀況。此外，本試驗(yàn)研究重點(diǎn)在于如何改進(jìn)底吹爐噴口結(jié)構(gòu)設(shè)計(jì)以便金屬熔體更快更好混合，流體性質(zhì)并非主要研究因素。因此，依據(jù)現(xiàn)存試驗(yàn)條件，本文提出利用氣泡擾動的流場圖像信息定義富氧底吹銅熔池熔煉爐內(nèi)氣液兩相流直接接觸混合系統(tǒng)的局部混合度，其計(jì)算公式如式(2)、式(3)所示。

圖2 單噴口工況下動態(tài)平衡流場的氣泡圖像處理過程

圖3 底吹爐噴口上方氣柱的形成過程

式中，ci為局部混合區(qū)域的灰度級強(qiáng)度；n為局部混合區(qū)域的像素個數(shù)。越大，說明圖像元素更多地趨于255，即局部混合區(qū)域白色部分占優(yōu)，氣體攪拌較好；反之，說明局部混合區(qū)域的混合質(zhì)量較差。S越小，說明圖像元素波動較小，但只是局部區(qū)域的混合均勻性較好的充分條件。

2 結(jié)果與討論

2.1 剖線分析

為驗(yàn)證本試驗(yàn)條件下所獲混合瞬態(tài)圖像能夠真實(shí)反映氣液混合效果，即圖像元素差異反映了氣液接觸系統(tǒng)中目標(biāo)（氣相）與背景（液相）的分離，以取自熔池動態(tài)平衡階段的圖2(e)為研究對象進(jìn)行剖線分析。圖4展示了圖2(e)中兩條水平線上的灰度級強(qiáng)度變化大小，其中橫軸為圖像水平像素?cái)?shù)（1280）、縱軸為對應(yīng)水平剖線上的灰度值。A→B提取的是氣泡擾動流場上表面的圖像特征，可以發(fā)現(xiàn)其像素灰度值基本穩(wěn)定在200附近，存在一些空隙的波動，此是由于底吹大氣泡在液體與空氣接觸表面破碎造成的；C→D提取的是氣泡擾動流場的中間區(qū)域特征，可以發(fā)現(xiàn)在C1穩(wěn)定時(shí)刻，局部混合區(qū)域的中間位置會形成較為穩(wěn)定的豎直氣柱，因此區(qū)域中間位置的灰度級強(qiáng)度值與兩端位置的灰度級強(qiáng)度值呈現(xiàn)明顯差異，整體序列出現(xiàn)單峰。C→D段的中間位置呈現(xiàn)噴口鼓出氣體形成的大氣泡，A→B段整體呈現(xiàn)了大氣泡破碎形成小氣泡或微氣泡的流型特征，因此從圖像分析的角度，兩者呈現(xiàn)相同的灰度級強(qiáng)度，說明氣泡的圖像元素值較大。

圖4 氣泡圖像[圖3(e)]頂部和中部的剖線灰度級強(qiáng)度

為了論證從灰度級強(qiáng)度判斷不同噴口設(shè)計(jì)優(yōu)劣的科學(xué)性與合理性，表2給出了本試驗(yàn)設(shè)計(jì)中6組工況剖線截取的位置參數(shù)，用以說明圖像元素強(qiáng)度差異能夠有效提取氣液兩相的圖像特征。其中A→B段全部截取的是氣泡擾動流場氣液分界面附近的水平剖線，C→D段和E→F段為氣泡擾動流場局部混合區(qū)域中間位置附近的水平剖線；由于在單噴口0°情形下的氣柱在圖像中豎直向上，因此不考慮E→F段（以Null表示）。為從BSD角度準(zhǔn)確分析和提取氣泡圖像特征，不同實(shí)驗(yàn)案例C→D段和E→F段的選擇依據(jù)實(shí)際特征明顯程度而定。

表2 不同工況條件下的氣泡圖像剖線位置參數(shù)

圖5(a)與圖5(b)取自熔池動態(tài)平衡階段，分別顯示了工況C2與C3中氣體攪拌流場混合瞬態(tài)的綠色成分圖像；圖6(a)與圖6(b)亦取自熔池動態(tài)平衡階段，分別顯示了實(shí)驗(yàn)案例C5及C6中氣體攪拌流場混合瞬態(tài)彩色圖像的綠色成分。而圖5(c)與圖5(d)、圖6(c)與圖6(d)則分別反映了圖5(a)與圖5(b)、圖6(a)與圖6(b)中水平剖線自A至B、自C至D及自E至F的灰度級強(qiáng)度值，用以觀測氣液兩相混合圖像的像素差異。鑒于實(shí)驗(yàn)案例C2與C3為單噴口結(jié)構(gòu)情形，發(fā)現(xiàn)其對應(yīng)的灰度級強(qiáng)度時(shí)間序列圖呈現(xiàn)單峰形狀特征；鑒于工況C5與C6為單排雙噴口結(jié)構(gòu)情形，發(fā)現(xiàn)其對應(yīng)的灰度級強(qiáng)度時(shí)間序列圖呈現(xiàn)雙峰的形狀特征；說明氣相的圖像元素對應(yīng)較高的強(qiáng)度值，而液相的圖形元素對應(yīng)較低的強(qiáng)度值，因此符合實(shí)際圖像特征提取，可從灰度級強(qiáng)度的角度判斷不同噴口設(shè)計(jì)的優(yōu)劣。

圖4、圖5(c)及圖5(d)分別對應(yīng)0°、7°及14°共3種單噴口傾角的灰度變化情形，圖5(c)中灰度級分布更為集中，氣液接觸表面氣泡破碎程度處于另外兩者之間，但局部混合區(qū)域灰度級整體較高，說明氣柱形態(tài)更加堅(jiān)挺，因此7°時(shí)的流場內(nèi)呈現(xiàn)出了較為強(qiáng)烈而連續(xù)的氣泡攪拌效果。比較圖6(c)與圖6(d)，可以發(fā)現(xiàn)后者雙峰間的距離要長于前者，說明前者的局部混合攪拌區(qū)域較大，即7°案例比14°工況表現(xiàn)更劇烈的攪拌效果。另外，局部混合區(qū)域的灰度級較大的范圍更大，說明雙排雙噴口設(shè)計(jì)創(chuàng)造了更加劇烈的攪拌效果，而雙排的雙噴口結(jié)構(gòu)參數(shù)明顯優(yōu)于單排單噴口的結(jié)構(gòu)參數(shù)，圖像分析結(jié)果符合實(shí)際。

除此之外，噴口排布方式對熔池的氣液兩相流的攪拌效果還主要體現(xiàn)在熔池氣含率、氣液兩相流混合時(shí)間、氣體在熔體內(nèi)停留時(shí)間等多個方面。但是鑒于實(shí)驗(yàn)條件所限，氣體在熔體內(nèi)停留時(shí)間的測定較繁瑣，且測定氣液兩相流混合時(shí)間缺乏實(shí)時(shí)性、干擾流場狀態(tài)的可能性較大（如電導(dǎo)率法）。對于混合瞬態(tài)時(shí)刻的情形，考慮基于常用的Otsu法計(jì)算熔池內(nèi)的氣含率，以此對本文所提的研究方法加以論證。

觀察發(fā)現(xiàn)：綠色通道圖像中的白色區(qū)域?yàn)闅馀莘植紖^(qū)域，深色區(qū)域?yàn)橐合啻媪魠^(qū)域，但是透明的容器壁亦呈現(xiàn)較高的像素值，因此選取圖2(e)、圖5(a)與圖5(b)、圖6(a)與圖6(b)中的上半?yún)^(qū)域（第1～360行，避免攪拌容器的壁面干擾熔池內(nèi)氣液兩相圖像元素）為研究區(qū)域，以研究區(qū)域內(nèi)氣相部分與氣液混合部分的占比作為熔池氣含率，對比研究區(qū)域內(nèi)灰度級強(qiáng)度值的均值和標(biāo)準(zhǔn)差變化。表3給出了6種情況下熔池氣含率與灰度級強(qiáng)度均值、標(biāo)準(zhǔn)差的比較結(jié)果（工況C1與C4中氣柱的位置相同，故未展現(xiàn)工況C4的對應(yīng)圖像），說明本文所提方法不僅與生產(chǎn)實(shí)際，而且與現(xiàn)存手段的研究結(jié)果基本一致，即雙噴口設(shè)計(jì)要優(yōu)于單噴口設(shè)計(jì)，但規(guī)避了灰度圖像的二值化，沒有閾值選擇的困擾，更加精確地體現(xiàn)圖像元素的差異。另外，工況C3的灰度標(biāo)準(zhǔn)差最小，此是由氣體攪拌區(qū)域太小、液相仍占據(jù)大部分范圍而引起的，體現(xiàn)了灰度信息的敏感性。

2.2 波動演化特性

鑒于高速攝像機(jī)在捕獲底吹熔池熔煉水模型中氣泡擾動流場圖像方面存在的局限性，截取原始圖像中不包含底部壁面的部分（圖像第1～400行）作為氣泡擾動流場波動演化特性表征的研究區(qū)域?；谇拔乃觯?.1節(jié)部分），選取原始?xì)馀軷GB圖像中的綠色成分，利用MATLAB軟件中的內(nèi)置函數(shù)（in-house functions）獲取所選區(qū)域灰度級像素的均值與標(biāo)準(zhǔn)差（即、S），其隨時(shí)間演化趨勢分別如圖7與圖8所示（對應(yīng)C1）。圖示均反映了案例C1從實(shí)驗(yàn)開始時(shí)刻到結(jié)束時(shí)刻氣泡擾動流場圖像特征的波動特性，灰度級像素均值序列顯示在初始階段的研究區(qū)域呈現(xiàn)由小變大的趨勢，在達(dá)到穩(wěn)定階段以后基本在130上下浮動；灰度級像素標(biāo)準(zhǔn)差序列同樣顯示出從小變大的初始特征，并在穩(wěn)定階段波動于55上下。

圖5 單噴口工況條件下實(shí)驗(yàn)案例C2與C3的氣泡攪拌圖像及其剖線灰度級強(qiáng)度

圖6 雙噴口工況條件下實(shí)驗(yàn)案例C5及C6的氣泡攪拌圖像及其剖線灰度級強(qiáng)度

表4給出了實(shí)驗(yàn)案例研究區(qū)域的氣含率時(shí)間序列、灰度級像素均值與標(biāo)準(zhǔn)差時(shí)間序列的統(tǒng)計(jì)參數(shù)值：序列均值和序列標(biāo)準(zhǔn)差。由圖示可知，雙噴口的結(jié)構(gòu)設(shè)計(jì)要優(yōu)于單噴口設(shè)計(jì)（C4、C5與C6灰度均值的數(shù)學(xué)期望均大于C1、C2與C3，與氣含率的變化趨勢相一致，但規(guī)避了圖像的深度處理，有利于避免信息損失），符合實(shí)際，亦說明本方法的有效性；各評價(jià)指標(biāo)隨時(shí)間變化值的標(biāo)準(zhǔn)差說明了實(shí)驗(yàn)過程的穩(wěn)定程度，標(biāo)準(zhǔn)差時(shí)間序列給出了評價(jià)混合質(zhì)量優(yōu)劣的充分條件；在單噴口情形下，0°設(shè)計(jì)擾動的流場強(qiáng)度更大；在單排雙噴口情形下，14°設(shè)計(jì)要優(yōu)于7°設(shè)計(jì)，此是因?yàn)槊總€噴口氣量相同，所選取的區(qū)域僅僅包含了靠近空氣的液面部分，14°設(shè)計(jì)對于該區(qū)域的攪拌作用更為明顯。

表3 熔池氣含率與灰度級強(qiáng)度均值、標(biāo)準(zhǔn)差比較結(jié)果

圖7 實(shí)驗(yàn)案例C1研究區(qū)域內(nèi)灰度級強(qiáng)度隨時(shí)間變化的均值序列

圖8 實(shí)驗(yàn)案例C1研究區(qū)域內(nèi)灰度級強(qiáng)度隨時(shí)間變化的標(biāo)準(zhǔn)差序列

表4 不同實(shí)驗(yàn)案例研究區(qū)域內(nèi)氣含率與灰度級強(qiáng)度時(shí)間序列統(tǒng)計(jì)參數(shù)

另外，為測試試驗(yàn)工況的實(shí)際數(shù)據(jù)是否服從正態(tài)分布，采用Lillietest測試（鑒于樣本容量很小，不采用Jarque-Bera測試）。對于均值序列，除C4外的結(jié)果為0，說明大部分穩(wěn)定階段的均值試驗(yàn)數(shù)據(jù)均服從正態(tài)分布；對于標(biāo)準(zhǔn)差序列，C1、C2及C6計(jì)算結(jié)果為0，C3、C4及C5計(jì)算結(jié)果為1，說明穩(wěn)定階段的標(biāo)準(zhǔn)差值試驗(yàn)數(shù)據(jù)并非全部服從正態(tài)分布。一方面，試驗(yàn)案例間的差異由于不同試驗(yàn)設(shè)計(jì)參數(shù)造成的，即氣流在上升過程中會發(fā)生搖擺震蕩；另一方面，試驗(yàn)過程中實(shí)驗(yàn)儀器的震動帶來一定的誤差。

3 結(jié)論

（1）以富氧底吹熔池熔煉過程為研究對象，搭建底吹煉銅爐水模型試驗(yàn)平臺，捕獲攪拌反應(yīng)器內(nèi)動態(tài)平衡階段氣泡擾動流場的RGB圖像。發(fā)現(xiàn)在本試驗(yàn)條件下，加權(quán)平均灰度化不能準(zhǔn)確提取流場圖樣特征；利用邊界輪廓圖，發(fā)現(xiàn)綠色成分對圖像特征識別與提取更為準(zhǔn)確。

（2）對6種試驗(yàn)案例中取自熔池動態(tài)平衡階段的混合瞬態(tài)進(jìn)行水平剖線分析，發(fā)現(xiàn)對于0°、7°及14°等單排單噴口案例與0°雙排的雙噴口案例，水平剖線的灰度級強(qiáng)度值呈現(xiàn)單峰；對于7°與14°單排雙噴口案例，水平剖線的灰度級強(qiáng)度值呈現(xiàn)雙峰；表明G成分能夠用于測度與比較氣泡擾動流場的攪拌混合質(zhì)量。

（3）測算研究區(qū)域的灰度級強(qiáng)度均值時(shí)間序列和標(biāo)準(zhǔn)差時(shí)間序列，以熔池氣含率作為論證參數(shù)，說明了灰度方法的有效性；同時(shí)發(fā)現(xiàn)大部分熔池動態(tài)平衡階段的時(shí)間序列在統(tǒng)計(jì)學(xué)意義上服從正態(tài)分布；采用數(shù)字圖像處理技術(shù)研究底吹攪拌反應(yīng)器的混合特性，為推進(jìn)透明或半透明攪拌反應(yīng)器氣液兩相混合特性研究提供了參考依據(jù)，一定程度上豐富和發(fā)展了冶金反應(yīng)過程多相流特性的研究手段。

[1] SRIDHAR R，TOGURI J M，SIMEONOV S. Copper losses and thermodynamic considerations in copper smelting[J]. Metallurgical and Materials Transactions B，1997，28（2）：191-200.

[2] SERGEI A D，ARTHUR D P. A thermodynamic database for copper smelting and converting[J]. Metallurgical and Materials Transactions B，1999，30（4）：661-669.

[3] 肖清泰，王華，王仕博，等. 基于改進(jìn)0-1測試的多相宏觀混合效果混沌表征[J]. 化學(xué)工程，2016，44（2）：46-51.XIAO Q T，WANG H，WANG S B，et al. Chao characterization of multiphase macromixing efficiency based on modified 0-1 test[J].Chemical Engineering，2016，44（2）：46-51.

[4] 閆紅杰，劉方侃，張振揚(yáng)，等. 氧槍布置方式對底吹熔池熔煉過程的影響[J]. 中國有色金屬學(xué)報(bào)，2012，22（8）：2393-2400.YAN H J，LIU F K，ZHANG Z Y，et al. Influence of lance arrangement on bottom-blowing bath smelting process[J]. The Chinese Journal of Nonferrous Metals，2012，22（8）：2393-2400.

[5] 張振揚(yáng)，陳卓，閆紅杰，等. 富氧底吹熔煉爐內(nèi)氣液兩相流動的數(shù)值模擬[J]. 中國有色金屬學(xué)報(bào)，2012，22（6）：1826-1834.ZHANG Z Y，CHEN Z，YAN H J，et al. Numerical simulation of gas-liquid multi-phase flows in oxygen enriched bottom-blowing furnace[J]. The Chinese Journal of Nonferrous Metals，2012，22（6）：1826-1834.

[6] 邵品，張廷安，劉燕，等. 底吹冰銅吹煉爐中氣-液流動狀況的數(shù)學(xué)模擬[J]. 東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，33（9）：1303-1306，1318.SHAO P，ZHANG T A，LIU Y，et al. Numerical simulation of gas-liquid flow behavior in a copper matte bottom-blowing converter[J]. Journal of Northeastern University（Natural Science），2012，33（9）：1303-1306，1318.

[7] 王東興，張廷安，劉燕，等. 氧氣底吹造锍過程中氣泡行為的水模型實(shí)驗(yàn)[J]. 東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，34（12）：1755-1758.WANG D X，ZHANG T A，LIU Y，et al. Water model study of bubble behavior in matte smelting process with oxygen bottom blowing[J].Journal of Northeastern University（Natural Science），2013，34（12）：1755-1758.

[8] 余躍，溫治，劉訓(xùn)良，等. 噴槍結(jié)構(gòu)對底吹煉銅爐流場影響的模擬及實(shí)驗(yàn)研究[J]. 中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，45（12）：4129-4137.YU Y，WEN Z，LIU X L，et al. Simulation and experimental of influence of nozzle structure on bottom-blowing furnace flowing process[J]. Journal of Central South University（Science and Technology），2014，45（12）：4129-4137.

[9] 郭學(xué)益，王親猛，廖立樂，等. 銅富氧底吹熔池熔煉過程機(jī)理及多相界面行為[J]. 有色金屬科學(xué)與工程，2014，5（5）：28-34.GUO X Y，WANG Q M，LIAO L L，et al. Mechanism and multiphase interface behavior of copper sulfide smelting in oxygen-enriched bottom blowing furnace[J]. Nonferrous Metals Science and Engineering，2014，5（5）：28-34.

[10] 郭學(xué)益，王親猛，田慶華，等. 基于區(qū)位氧勢硫勢梯度變化下銅富氧底吹熔池熔煉非穩(wěn)態(tài)多相平衡過程[J]. 中國有色金屬學(xué)報(bào)，2015，25（4）：1072-1079.GUO X Y，WANG Q M，TIAN Q H，et al. Non-steady multiphase equilibrium process of copper oxygen-enriched bottom blowing bath smelting with gradual change of oxygen and sulfur potential of different positions in furnace[J]. The Chinese Journal of Nonferrous Metals，2015，25（4）：1072-1079.

[11] 郭學(xué)益，王親猛，田慶華，等. 氧氣底吹熔煉工藝分析及過程優(yōu)化[J]. 中國有色金屬學(xué)報(bào)，2016，26（3）：689-698.GUO X Y，WANG Q M，TIAN Q H，et al. Analysis and optimization of oxygen bottom blowing copper smelting process[J]. The Chinese Journal of Nonferrous Metals，2016，26（3）：689-698.

[12] LI Y，ZHANG J P，F(xiàn)AN L S. Discrete-phase simulation of single bubble rise behavior at elevated pressures in a bubble column[J].Chemistry Engineering Science，2000，55（20）：4597-4609.

[13] YANG G，DU B，F(xiàn)AN L. Bubble formation and dynamics in gas-liquid-solid fluidization—a review[J]. Chemistry Engineering Science，2007，62（1）：2-27.

[14] LAU Y M，DEEN N G，KUIPERS J A M. Development of an image measurement technique for size distribution in dense bubbly flows[J].Chemical Engineering Science，2013，94：20-29.

[15] KARN A，ELLIS C，ARNDT R，et al. An integrative image measurement technique for dense bubbly flows with a wide size distribution[J]. Chemical Engineering Science，2015，122：240-249.[16] XIAO Q，PAN J，XU J，et al. Hypothesis-testing combined with image analysis to quantify evolution of bubble swarms in a direct-contact boiling heat transfer process [J]. Applied Thermal Engineering，2016，113：851-857.

[17] XIAO Q，XU J，WANG H. Quantifying the evolution of flow boiling bubbles by statistical testing and image analysis: toward a general model[J]. Scientific Reports，2016，6：31548.

[18] XIAO Q，WANG S，ZHANG Z，et al. Analysis of sunspot time series（1749—2014）by means of 0-1 test for chaos detection[C]//The 11th International Conference on Computational Intelligence and Security，2015：215-218.

Study on mixing characteristics in bottom-blowing stirred reactor using image analysis

XIAO Qingtai，WANG Shibo，LI Peng，GAO Qin，XU Jianxin，WANG Hua
（State Key Laboratory of Complex Nonferrous Metal Resources Clean Utilization，Kunming University of Science and Technology，Kunming 650093，Yunnan，China）

The oxygen-enriched bottom-blowing smelting process is taken as research object. In order to measure and compare the local mixing characteristics in the bottom-blowing stirred reactor，the grayscale intensity was obtained based on field images of bubble perturbation. According to the actual size of bottom-blowing copper melting furnace，the test platform for water model of gas blowing was built. By watching high-speed videos of experimental process at different design parameters，the field RGB images of water model above the nozzles were captured and saved. Only the area undisturbed by apparatus was extracted as the study area. Finally，the grayscale intensity average and standard deviation of green component can be obtained. It was found that the green component of bubble RGB images extracted the contour features of actual bubble images more than other grayed images at our experimental conditions. The horizontal profiles of green component present obvious single apex and double peak，corresponding to the operating conditions of single-line nozzle and double-line nozzles，respectively. These results proved the validity of grayscale intensity for characterizing the mixing process. The mean and standard deviation time series of grayscale intensity in study area were calculated.Results showed that most times series follow Gaussian distributions statistically at the dynamic equilibrium stage of the bath. This work provides some academic bases and references on pushing the research about adopting image analysis techniques to investigate flow and mixing in the bottom-blowing stirred reactors.

gas-liquid flow；distributions；stirred vessel；mixing；imaging

TF811

：A

：1000-6613（2017）09-3195-08

10.16085/j.issn.1000-6613.2017-0071

2017-01-13；修改稿日期：2016-03-02。

國家自然科學(xué)基金（51666006，51406071，U150220046）、張文海院士工作站（2015IC005）及云南省科技領(lǐng)軍人才項(xiàng)目（2015HA019）。

肖清泰（1989—），男，博士研究生，主要從事工業(yè)圖像處理研究。E-mail：qingtaixiao2016@kmust.edu.cn。聯(lián)系人：徐建新，副教授，碩士生導(dǎo)師，主要從事多相流可視化研究。E-mail：xujianxina@163.com。