靳海波，連以誠，劉鑫，林金芝，楊索和，何廣湘

?

基于電阻層析成像技術(shù)的氣升式內(nèi)循環(huán)鼓泡反應(yīng)器流體力學(xué)研究

靳海波，連以誠，劉鑫，林金芝，楊索和，何廣湘

(北京石油化工學(xué)院化學(xué)工程學(xué)院，北京，102617)

采用電阻層析成像(ERT)技術(shù)，對高為1 024 mm、直徑為165 mm，導(dǎo)流筒高為590 mm、直徑為80 mm的氣升式內(nèi)循環(huán)反應(yīng)器進(jìn)行流體力學(xué)參數(shù)的測量?？疾毂碛^氣速對上升區(qū)和下降區(qū)內(nèi)等3個(gè)截面處氣含率的影響以及反應(yīng)器內(nèi)流型的變化情況。研究結(jié)果表明：在上升區(qū)，氣含率隨表觀氣速的變化與普通鼓泡塔情況一樣，而在下降區(qū)，一部分氣泡由于氣液湍動被帶入下降區(qū)；隨著表觀氣速的增加，液體的循環(huán)速度逐漸增大，更多的氣泡被拖曳到下降區(qū)，下降區(qū)內(nèi)氣含率增加明顯。通過對ERT圖像時(shí)間序列的疊加，可以清晰得出在上升區(qū)低氣速下氣泡離散式上升，隨著表觀氣速的增加，以氣泡群的方式流動且氣泡聚并明顯，并出現(xiàn)氣泡群的擺動；在下降區(qū)，隨著液體循環(huán)速度增加，下降區(qū)氣液層消失，氣含率逐漸增大。

氣升式內(nèi)循環(huán)反應(yīng)器；電阻層析成像；氣含率；時(shí)間序列流動結(jié)構(gòu)

氣升式內(nèi)循環(huán)反應(yīng)器是氣?液、氣?液?固多相之間進(jìn)行質(zhì)量、動量、熱量傳遞以及化學(xué)反應(yīng)的重要反應(yīng)器之一。內(nèi)循環(huán)流動模式改進(jìn)了反應(yīng)器內(nèi)流體流動形式，增強(qiáng)了流體的流動，強(qiáng)化了相間的傳質(zhì)、傳熱效率。因具有結(jié)構(gòu)簡單、造價(jià)低、運(yùn)行成本低、無運(yùn)動部件、傳熱和傳質(zhì)特性好、流體流動均勻、混合效率高等特點(diǎn)，在石油工業(yè)、生物化工、化學(xué)工業(yè)、濕法冶金、煤液化及環(huán)境保護(hù)等方面有著廣泛的應(yīng)用，如費(fèi)托合成、一步法合成甲醇和二甲醚、汽油脫硫、重油加氫、生物廢水處理以及發(fā)酵工程等[1?2]。電層析成像技術(shù)(ERT)是20世紀(jì)90年代發(fā)展起來的一種以兩相流或多相流為主要對象的過程參數(shù)二維或三維分布狀況的在線實(shí)時(shí)監(jiān)測技術(shù)。它將傳統(tǒng)的單點(diǎn)、局部的測量，發(fā)展為過程參數(shù)的多點(diǎn)、截面分布式的測量，在不破壞、不干擾流體流動的情況下，獲得設(shè)備內(nèi)部兩相/多相流體的分布信息，因此，在多相流測量過程中顯示出明顯的優(yōu)勢[3?5]。ERT技術(shù)的測量原理是基于不同的媒質(zhì)具有不同的電導(dǎo)率，判斷出處于敏感場的物體電導(dǎo)率分布，便可知物場的媒質(zhì)分布狀況。當(dāng)場內(nèi)的電導(dǎo)率分布變化時(shí)，場的分布也會隨之變化，導(dǎo)致場內(nèi)電勢分布變化，從場域邊界上的測量電壓也要發(fā)生變化。邊界測量電壓變化包含了場域內(nèi)電導(dǎo)率變化的信息。通過實(shí)際對象的邊界測量電壓與連續(xù)相為滿管狀態(tài)下對應(yīng)位置測量電壓的對比，運(yùn)用相應(yīng)的成像算法，便可重建實(shí)際對象的電導(dǎo)率分布，實(shí)現(xiàn)可視化測量，同時(shí)也可以提取流型、相含率等特征參數(shù)。由于電阻層析成像(ERT)技術(shù)具有可視化、無輻射、非侵入、成本低等優(yōu)點(diǎn)，被廣泛用于測量多相流動力學(xué)參數(shù)[6?15]。JIN等[13?15]利用電阻層析成像技術(shù)對鼓泡塔內(nèi)氣液相行為參數(shù)進(jìn)行了測量研究，并得到很好的結(jié)果。WANG等[16]利用ERT技術(shù)對鼓泡塔內(nèi)氣液兩相流進(jìn)行測量，并發(fā)展了ERT的可視化技術(shù)。韓玉環(huán)等[17]利用ERT對氣液固三相外環(huán)流反應(yīng)器相含率和氣液流動結(jié)構(gòu)進(jìn)行研究，取得了良好的結(jié)果。GUMERY等[18]則利用ERT對氣升式內(nèi)循環(huán)反應(yīng)器的混合動力學(xué)參數(shù)進(jìn)行了測量，研究了流體黏度、氣體分布器結(jié)構(gòu)對反應(yīng)器混合效率的影響。本文作者主要利用電阻層析成像技術(shù)對氣升式內(nèi)循環(huán)反應(yīng)器進(jìn)行測量，考察表觀氣速對反應(yīng)器上升區(qū)和下降區(qū)的影響，并研究反應(yīng)器內(nèi)流體流型在表觀氣速影響下的變化，為進(jìn)一步實(shí)驗(yàn)研究奠定基礎(chǔ)。

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

圖1所示為實(shí)驗(yàn)裝置。實(shí)驗(yàn)系統(tǒng)分為3部分：進(jìn)氣系統(tǒng)，反應(yīng)器和測量系統(tǒng)(見圖1)。進(jìn)氣系統(tǒng)由空氣壓縮機(jī)、流量計(jì)等組成。整個(gè)反應(yīng)器由有機(jī)玻璃制成，高1 024 mm，外徑165 mm；導(dǎo)流筒高590 mm，外徑80 mm；導(dǎo)流筒底隙高40 mm(導(dǎo)流筒底部到反應(yīng)器底部的距離)。進(jìn)氣口開孔18 mm。電阻層析成像測量系統(tǒng)，有傳感電極和數(shù)據(jù)采集、成像系統(tǒng)組成。其中電極片長20 mm，寬6 mm，厚2 mm。從上到下，ERT電極相隔距離依次是30和250 mm。導(dǎo)流筒內(nèi)部電極和外部鏈接的部分套有絕緣密封套?？諝庥蛇M(jìn)氣孔進(jìn)入反應(yīng)器導(dǎo)流筒內(nèi)，導(dǎo)致反應(yīng)器內(nèi)電導(dǎo)率的變化，通過ERT測量、圖像處理系統(tǒng)可以實(shí)時(shí)地反映出反應(yīng)器內(nèi)混合流體流動的情況。

本實(shí)驗(yàn)采用ITS電阻層析成像系統(tǒng)(P2+ v7.0?2009)，實(shí)驗(yàn)前先校正ERT數(shù)據(jù)采集系統(tǒng)(DAS)。實(shí)驗(yàn)過程中，在常溫、常壓下，表觀氣速從0.043 m/s到0.433 m/s，每次分別對Plane1，Plane2和Plane3進(jìn)行數(shù)據(jù)采集，每次采集400幀圖像。由于導(dǎo)流筒采用有機(jī)玻璃制成，可以看出絕緣體，通過提高激勵電流來消除導(dǎo)流筒對下降區(qū)測量結(jié)果的影響。

圖1 實(shí)驗(yàn)裝置

2 實(shí)驗(yàn)結(jié)果與討論

2.1 局部氣含率

氣含率是氣液反應(yīng)器中重要的動力學(xué)參數(shù)之一，它直接影響到相接觸面積，最終影響到相間傳質(zhì)、傳熱以及反應(yīng)速率等。因此，對氣含率及其分布特性的考察對于研究反應(yīng)器放大規(guī)律有著非常重要的意義。

電阻層析成像系統(tǒng)每幀可測量316個(gè)電導(dǎo)數(shù)據(jù)，然后根據(jù)Maxwell方程得到下式，可計(jì)算得出局部氣含率為

(2)

圖2所示為在Plane1，Plane2和Plane3處氣含率隨著表觀氣速的變化。從圖2可以看出：上升區(qū)(Plane1)、下降區(qū)(Plane2, Plane3)均隨表觀氣速的增加而增大；在上升區(qū)，氣含率隨表觀氣速的變化和常規(guī)鼓泡塔內(nèi)情況一樣，在低氣速(小于0.13 m/s)時(shí)，上升區(qū)內(nèi)氣泡呈離散型運(yùn)動，表觀氣速對氣含率的影響為增加氣泡的數(shù)量，此時(shí)的流型為鼓泡流；當(dāng)氣速繼續(xù)增大時(shí)，氣泡在上升區(qū)內(nèi)聚并開始加劇，使得大氣泡數(shù)量增多，進(jìn)而上升區(qū)內(nèi)氣含率增長速率變大，但由于大氣泡在上升區(qū)內(nèi)上升速度快，停留時(shí)間短，使得在上升區(qū)內(nèi)氣含率變化又趨于平衡增長，并且小于低氣速下氣含率的增長速率；而在下降區(qū)，由于Plane2處靠近導(dǎo)流筒頂端，而在導(dǎo)流筒頂端由于氣液分離湍動劇烈，一部分氣泡由于氣液湍動流動被帶入下降區(qū)頂端，因此，低氣速時(shí)，Plane2處氣含率和上升區(qū)氣含率相差不大。但此時(shí)由于表觀氣速小，下降區(qū)液體流動比較緩慢，氣泡只是在下降區(qū)頂端部分，氣泡受到的升力遠(yuǎn)遠(yuǎn)大于液體對氣泡的曳力，此時(shí)氣泡逃逸出下降區(qū)，因此低氣速時(shí)Plane3處的氣含率幾乎為0，并且在下降區(qū)有明顯的氣液層，在氣液層以上氣泡數(shù)量很多，氣含率很大，在氣液層以下，氣含率幾乎為0。隨著表觀氣速的增加，液體在整個(gè)反應(yīng)器內(nèi)開始循環(huán)流動，并且循環(huán)速度逐漸增大，更多的氣泡被拖曳到下降區(qū)，下降區(qū)內(nèi)的氣液層也隨之逐漸消失，此時(shí)Plane3處的氣含率增加速率突然加快。當(dāng)增大表觀氣速時(shí)，Plane3處的氣泡流動趨于穩(wěn)定，氣含率也基本趨于穩(wěn)定。而在Plane2處，由于下降區(qū)頂端湍動比較劇烈，低氣速時(shí)氣含率受到液體循環(huán)速度的影響較小，之后，液體循環(huán)速度增加，越來越多的氣泡被拖曳到下降區(qū)，氣含率隨表觀氣速的增加而增大，但增長速率趨于恒定，此處氣泡在每個(gè)階段都處于穩(wěn)定階段，氣泡停留時(shí)間長。因此，Plane2是氣液傳質(zhì)、傳熱效率比較高的地方。

1—Plane1；2—Plane2；3—Plane3。

2.2 徑向氣含率

氣含率的徑向分布也是氣升式反應(yīng)器的重要參數(shù)之一，分析氣含率徑向分布情況也有助于理解流型轉(zhuǎn)變、液體混合、傳質(zhì)以及傳熱性能。氣泡在反應(yīng)器內(nèi)所受的徑向力對氣含率的徑向分布意義很大，氣泡所受的徑向力包括升力、壁面力、湍動擴(kuò)散力等。由于氣泡所受的徑向力不同，導(dǎo)致氣含率在不同的流型下有不同的徑向分布[19]。OHNUKI等[20]通過對氣含率徑向分布的研究，認(rèn)為氣液兩相流中氣含率的徑向分布情況有2種：1) 在均勻鼓泡區(qū)或較低氣速下，氣含率為邊壁峰或均勻分布；2) 在非均勻鼓泡區(qū)或較高氣速下，氣含率呈現(xiàn)明顯的中心峰分布。圖3~5分別為Plane1~3處不同表觀氣速下氣含率的徑向分布。在較低氣速下，氣含率的幾乎呈均勻分布，而在高氣速下，氣含率呈峰中心分布較明顯，這和OHNUKI等[20]的結(jié)論是一致的。

圖3 Plane1處氣含率徑向分布

圖4 Plane2處氣含率徑向分布

圖5 Plane3處氣含率徑向分布

2.3 ERT可視化

電阻層析成像最重要的優(yōu)點(diǎn)之一是能夠?qū)崟r(shí)、可視化地顯示反應(yīng)器內(nèi)流動狀態(tài)，從 ERT所采集的圖像疊加得到的時(shí)間序列圖中可以反映出上升區(qū)、下降區(qū)流體流動形態(tài)。通過對每一橫截面時(shí)間序列圖進(jìn)行T3D處理發(fā)現(xiàn)：在上升區(qū)，低氣速下氣泡離散式上升，隨著表觀氣速的增加，液體開始在反應(yīng)器內(nèi)循環(huán)流動，隨著液體循環(huán)速度增大，上升區(qū)內(nèi)氣泡聚并明顯，此時(shí)上升區(qū)內(nèi)以氣泡群的方式流動，并且出現(xiàn)了氣泡群的擺動，而在氣液分離時(shí)，被拖曳到下降區(qū)內(nèi)的氣泡數(shù)量以及氣含率也隨著表觀氣速的增加而增加。因此，在下降區(qū)，隨著液體循環(huán)速度增加，下降區(qū)氣液層消失，氣含率逐漸增大。這和實(shí)驗(yàn)中所觀察到現(xiàn)象是一致的。

(a) Plane 1；(b) Plane 2；(c) Plane 3

3 結(jié)論

1) 電阻層析成像技術(shù)對含有導(dǎo)流筒的兩相流反應(yīng)器測量是可行的，并對電阻層析成像采集圖像進(jìn)行處理，給出了氣相在反應(yīng)器整個(gè)截面內(nèi)軸向分布的情況。在上升區(qū)，氣含率隨表觀氣速的變化和常規(guī)鼓泡塔內(nèi)情況一樣；而在下降區(qū)，由于在導(dǎo)流筒頂端氣液分離湍動劇烈，一部分氣泡由于氣液湍動流動被帶入下降區(qū)頂端；隨著表觀氣速的增加，液體的循環(huán)速度逐漸增大，更多的氣泡被拖曳到下降區(qū)，下降區(qū)內(nèi)氣含率增加速率加快。

2) 氣含率隨表觀氣速的增加而增大，在較低氣速下，氣含率分布均勻，而在高氣速下，氣含率呈峰中心分布較明顯。

3) 在上升區(qū)，低氣速下氣泡離散式上升；隨著表觀氣速的增加，以氣泡群的方式流動且上升區(qū)內(nèi)氣泡聚并明顯，并且出現(xiàn)了氣泡群的擺動；在下降區(qū)，隨著液體循環(huán)速度增加，下降區(qū)氣液層消失，氣含率逐漸增大。

參考文獻(xiàn)：

[1] 何廣湘, 楊索和, 靳海波. 氣升式環(huán)流反應(yīng)器的研究進(jìn)展[J]. 化學(xué)工業(yè)與工程. 2008, 25(1): 65?71. HE Guangxiang, YANG Suohe, JIN Haibo. Development of air loop reactors[J]. Chemical Industrial and Engineering, 2008, 25(1):65-71.

[2] NIE Dashi, CUI Yingying, ZHANG Qiang. Characteristics and application of airlift loop reactor[J]. Chemical Industry & Engineering, 2004, 25(6): 6?10.

[3] WANG F, MARASHDEH Q, FAN L S, et al. Electrical capacitance, electrical resistance, and positron emission tomography techniques and their applications in multi-phase flow systems[J]. Adv in Chem Eng, 2009, 37: 179?222.

[4] ALBION K J, BRIENS L, BRIENS C, et al. Multiphase flow measurement techniques for slurry transport[J]. Inter J of Chem Reac Eng, 2011, 9(1): 1?50.

[5] MENG Z, HUANG Z, WANG B, et al. Air-water two-phase flow measurement using a venture meter and an electrical resistance tomography sensor[J]. Flow Meas Instrum, 2010, 21(3): 268?176.

[6] WANG M, MANN R, DICKIN F J. Electrical resistance tomography sensing system for industrial applications[J]. Chem Eng Commun, 1999, 175(1): 49?70.

[7] JIN Haibo, WANG M, WILLIAMS R A. The effect of sparger geometry on gas bubble ?ow behaviors using electrical resistance tomography[J]. Chin J Chem Eng, 2006, 14(1): 127?131.

[8] WILLIAMS R A, WANG M. Dynamic imaging of process plant reactors and separators using electrical process tomography[J]. Oil Gas Sci Technol, 2000, 55(2): 185?186.

[9] VLAEV D S, WANG M, DYAKOWSKI T, et al. Detecting ?lter-cake pathologies in solid-liquid ?ltration: semitech scale demonstrations using electrical resistance tomography[J]. Chem Eng J, 2000, 77(1/2): 87?92.

[10] WANG M. Impedance mapping of particulate multiphase ?ows[J]. Flow Meas Instrum, 2005, 16(2/3): 183?189.

[11] 董峰, 劉小平, 鄧湘, 等. 電阻層析成像技術(shù)在兩相管流測量中的應(yīng)用[J]. 化工自動化及儀表, 2001, 28(6): 50?54. DONG Feng, LIU Xiaoping, DENG Xiang, et al. Application of electrical resistance tomography in two-phase flow measurement [J]. Control and Instruments in Chemical Industry, 2001, 28(6): 50?54

[12] JIN H B, YANG S H, WANG M, et al. Measurement of gas holdup profiles in a gas liquid cocurrent bubble column using electrical resistance tomography[J]. Flow Meas Instrum, 2007, 18(5/6): 191?196.

[13] JIN H, WANG M, WILLIAMS R A. Analysis of bubble behaviors in bubble column using electrical resistance tomography[J]. Chem Eng J, 2007, 130(2/3): 179?185.

[14] JIN H, LIAN Y, YANG S, et al. The Parameters Measurement of air-water two phase flow using electrical resistance tomography (ERT) technique in a bubble column[J]. Flow Meas Instrum, 2013, 31(6): 55?60.

[15] JIN H, LIAN Y, QIN Y, et al. Distribution characteristics of holdups in a multi-stage bubble column using electrical resistance tomography[J]. Particuology, 2013, 11(2): 225?231.

[16] WANG M, JIA X, BENNET M, et al. Bubble column measurement and control using electrical resistance tomography[C]// XIE H, WANG Y, JIANG Y, eds. Computer Application in the Minerals Industries. A A Balkema/Rotterdam, Netherlands, 2001: 459?464.

[17] 韓玉環(huán), 楊索和, 靳海波, 等. 氣液固三相外環(huán)流反應(yīng)器相含率分布與氣液流動結(jié)構(gòu)[J]. 過程工程學(xué)報(bào), 2010, 10(5): 862?867. HAN Yuhuan, YANG Suohe, JIN Haibo, et al. Distribution characteristics of phase holdups and flow structure of gas?liquid in a three-phase external loop airlift reactor[J]. The Chinese Journal of Process Engineering, 2010, 10(5): 862?867.

[18] GUMERY F, EIN-MOZAFFARI F, DAHMAN Y. Macromixing hydrodynamic study in draft-tube airlift reactors using electrical resistance tomography[J]. Bioprocess Biosyst Eng, 2011, 34(2): 135?144.

[19] LUCAS D, KREPPER E, PRASSER H M. Prediction of radial gas profiles in vertical pipe flow on the basis of bubble size distribution[J]. Int J Therm Sci, 2001, 40(3): 17?225.

[20] OHNUKI A, AKIMOTO H. Experimental study on transition of flow pattern and phase distribution in upward air-water two-phase flow along a large vertical pipe[J]. Intl J Multi Flow, 2000, 6(3): 67?386.

(編輯 陳愛華)

Hydrodynamic parameters in an airlift inner-loop column using electrical resistance tomography

JIN Haibo, LIAN Yicheng, LIU Xin, LIN Jinzhi, YANG Suohe, HE Guangxiang

(Department of Chemical Engineering, Beijing Institute of Petrochemical Technology, Beijing 102617, China)

The hydrodynamic parameters are studied in an airlift inner-loop reactor (height 1 024 mm, diameter 165 mm) with draft tube (height 590 mm, diameter 80 mm) using electrical resistance tomography (ERT) technology. The influence of gas superficial velocity on gas holdup of three cross sections was presented and the transformation of flow structure was also discussed. The results show that gas holdups in the rise area increased with the increase of the superficial gas velocity. In addition, liquid circulation velocity increases gradually and a part of bubbles drops into down-comer zone due to the gas-liquid turbulent with the increase of superficial gas velocity, and the gas holdup in down-comer zone increases significantly. Moreover, the transformation of flow structure in the reactor can be seen clearly with the result of composition of ERT time series images. Bubbles in the rising area increase with the discrete type under the low gas velocity, and the bubble swarms rise with the spiral escalation due to the rising bubble coalescence with the increase of superficial gas velocity. Meanwhile, this technique offers the possibility for noninvasive and nonintrusive visualization of ?ow ?elds in the reactor and provides a foundation for the future research.

airlift inner-loop reactor; electrical resistance tomography; gas holdup; flow structure with time series

10.11817/j.issn.1672-7207.2016.11.042

TQ021

A

1672?7207(2016)11?3935?05

2016?01?03；

2016?04?15

國家自然科學(xué)基金資助項(xiàng)目(21073020)；北京市屬高等學(xué)校人才強(qiáng)教深度計(jì)劃項(xiàng)目(PHR200906139)；北京市屬高等學(xué)校高層次人才引進(jìn)與培養(yǎng)計(jì)劃項(xiàng)目(CIT&TCD20130325) (Project(20776018) supported by the National Natural Science Foundation of China；Project(PHR200906139) supported by Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality; Project(CIT&TCD20130325) supported by the Importation and Development of High-Caliber Talents Project of Beijing Municipal Institutions)

靳海波，博士，教授，從事多相流反應(yīng)工程研究；E-mail: jinhaibo@bipt.edu.cn