加焦方式對氣化爐內氣流分布的影響

李海峰，游　洋，韓立浩，羅志國，蔡九菊，鄒宗樹

(東北大學材料與冶金學院，沈陽110819)

?

加焦方式對氣化爐內氣流分布的影響

李海峰，游洋，韓立浩，羅志國，蔡九菊，鄒宗樹

(東北大學材料與冶金學院，沈陽110819)

摘要:布料模式決定了料床的空隙度，而料床的空隙度分布決定了煤氣流的二次分布.本文建立了研究三維氣化爐爐料結構和煤氣流分布的物理模型和數學模型，物理模型采用熱電偶測溫的方法，從爐內氣體溫度分布信息考察了氣體的流動情況，由于物理實驗無法獲得內部空隙度分布信息，故基于離散單元法模型，以Fluent軟件為載體，利用多孔介質模型并加入用戶自定義函數，通過數學模型進一步研究了不同加焦方式下氣化爐內煤氣流分布的影響機理，獲得了氣化爐內煤氣的速度場和流線.物理模擬與數值模擬結果相吻合.通過模擬計算獲得的非均勻床層內氣體流動規律的認識對COREX氣化爐加焦工藝有借鑒意義.

關鍵詞:COREX熔化氣化爐;煤氣流分布;加焦模式;物理模擬;數值模擬

COREX熔融還原法是奧鋼聯開發的非焦煉鐵技術，也是第一個實現工業化的熔融還原技術.通過COREX熔融還原煉鐵將作為鋼鐵企業的前道工序，為煉鋼工序生產鐵水［1，2］.COREX熔融還原裝置主體分為上下兩部分，上部的預還原豎爐和下部的熔化氣化爐.其中下部的熔化氣化爐內的煤氣流分布從形成到排出爐外經過了兩次分布.在風口處回旋區內，由鼓入的氧氣與半焦燃燒生成高溫氣體，形成了煤氣流的初始分布;煤氣流經料床至排出爐外為其第二次分布，料床內爐料構成的空隙度分布決定著煤氣流的第二次分布，而料床的空隙度由爐料的布料模式及爐料運動行為決定.

本文根據氣化爐布料模式［3，4］以及借鑒前人關于高爐布料研究方法［5～7］的文獻，對不同加焦方式下可能形成的料床結構進行研究分析.加焦方式分為三種，即均勻混合、中心區域加焦、中間區域加焦.其中，均勻混合是目前COREX熔化氣化爐采用的一種布料方式，爐料混合均勻是與高爐層狀爐料結構差異最大的地方，故其爐內煤氣流有別于高爐煤氣發展方式.在實際的高爐生產中，時常會采用中心加焦或中間加焦的方式來調節煤氣流分布，這種方式可促使煤氣合理分布，以達到提高煤氣利用率的目標.然而此布料方式(中心加焦、中間加焦的布料模式)在熔化氣化爐并未采用，其對煤氣流分布影響的研究文獻尚少.故本文首次提出在熔化氣化爐內使用中心或中間加焦的方式調節爐內煤氣流分布狀況，通過物理模型和數學模型相結合的方法，考察了不同加焦方式對氣化爐內煤氣流分布的影響程度.

1　物理模型建立

首先在實驗室搭建了熔化氣化爐物理模擬實驗裝置，如圖1所示.該裝置以寶鋼COREX 3000原型，按照30∶1進行縮小的三維扁平模型，其中模型本體為不銹鋼結構，模型厚度為100 mm;主體兩側各有兩個風口由側壁插入，向下傾角4(°)，深入爐內2 mm;實驗物料由模型頂部加入.采用3 mm的玉米粒子模擬塊煤(或焦炭)和3 mm的石蠟粒子模擬DRI.由風口處鼓入熱風并通過調節螺旋排料器排料，轉速控制排料量，模擬焦炭燃燒過程.為了觀察物料的運動及熔化過程，正面采用可視化的鋼化玻璃面板，同時在背面開有196個溫度采集孔，以便測量爐內氣體溫度分布信息.經前人文獻［8－11］調研，在進行常溫實驗時，通常采用壓力傳感器測定壓力場.在高溫實驗時，因模擬物料石蠟熔化會堵塞壓力傳感器的測量孔，導致無法測量，因此，本實驗通過測定床層溫度場來間接反映煤氣流動的信息.另外，由于氣體與溫度探頭是全接觸，而顆粒與溫度探頭為點接觸，因此本文認為所測得的溫度場為氣體溫度場.為保證爐體的保溫效果，減輕壁面帶來的影響，觀察面板采用兩層鋼化玻璃結構，除觀察面板外，爐體周圍覆蓋保溫材料.

當氣體由風口鼓入后，在爐內進行二次分布，由于裝置背面的熱電偶在模型內部與顆粒為點接觸，周圍主要被鼓入的氣體圍繞，故熱電偶測得的數據可看作為是顆粒周圍氣體的溫度.考慮到爐內氣體溫度變化大致可反映出爐內氣流分布狀況，所以采用熱電偶測得的溫度值來間接分析爐內氣流分布狀況.

所謂加焦就是借助從爐頂向爐內另外添加少量焦炭以形成焦柱，用于減小爐內狹小范圍內的礦焦比，使其料柱透氣性改善，從而通過更多氣流.本文模擬的焦炭加入量參考高爐中心加焦量，中心加焦寬度選擇為3 cm，經過爐型體積計算，得到中心加焦量占加焦總量的8.7%;因高爐中暫無中間區域加入焦炭的研究，故本文中間區域加焦量與中心加焦量一致.

圖1　實驗裝置示意圖Fig.1 Schematic diagram of experimental apparatus1—熔化氣化爐模型; 2—溫度變送器; 3—工控機; 4—風口; 5—爐缸底部機械傳動機構; 6—爐底擋板; 7—熱電偶; 8—氣體分配器; 9—熱風機; 10—流量計; 11—鼓風機

2　數學模型建立

2.1離散單元模型及模型參數選取

離散單元法［12］適用于模擬散體物料的運動，通過對每一個顆粒進行受力分析，獲得顆粒的運動行為.模型中顆粒與顆粒發生碰撞時的作用力如圖2所示，其中顆粒i、顆粒j均受到兩種力和兩種力矩的作用.作用力包括顆粒自身重力、顆粒－顆?；蝾w粒－壁面之間的接觸力.力矩包括切向力矩和滾動摩擦力矩.根據牛頓第二定律，可獲得顆粒的運動方程，公式如下:

式中涉及的符號意義及顆粒間或顆粒與壁面間的接觸力及力矩的計算公式見文獻［13］.

實驗室具備IBM P55A工作站，其計算能力上限為50萬顆粒，并不能完成工業尺度數量級計算.因此，為簡化計算，并根據B.Wright等［14］的研究，本文將氣化爐實際幾何尺寸按1∶30縮小，并取其3D扁平模型進行計算，風口回旋區假設為球形，計算區域示意圖如圖3所示.表1給出了COREX熔化氣化爐物料運動DEM模型的模擬條件.

圖2　顆粒間的受力分析圖Fig.2 Schematic diagram of forces in two particles

圖3　計算幾何區域Fig.3 Computation geometry

表1　爐料物理性質和模擬計算條件Table 1 Physical properties of burden and simulated conditions

2.2煤氣流分布的流體力學模型及邊界條件設置

在建立流體力學模型過程中，為減少模型復雜性，作了如下簡化和假設:

( 1)爐內化學反應狀態已達到穩態，風口回旋區內只有氣體流動;

( 2)填充床、軟熔區及半焦床等各區域間的分界線明確，數據均來自物理實驗結果［15］.連續性方程

動量守恒方程

湍流方程

本模擬研究使用Ergun方程［16］來推導多孔介質中氣流的黏性和慣性阻力損失系數.

每一方向上的黏性阻力損失系數和慣性阻力損失系數為

式( 3)～( 9)中涉及的符號意義見文獻［17］.

通過以上控制方程和假設，建立了氣化爐內煤氣流分布的三維數學模型，方程求解借助商業軟件FLUENT，采用SIMPLE算法求解，方程的離散采用一階迎風格式，收斂準則殘差小于1× 10－4.采用ANSYS ICEM對氣化爐模型進行網格劃分，網格類型為四面體結構性網格，采用Map方法劃分，為了使模擬結果更加準確，對氣體入口—風口處的網格進行了加密處理.模型計算區域主要有自由空間、填充床、軟熔區域、半焦床、回旋區、死料柱等多個區域，均按多孔介質處理.其中半焦床、軟熔區域、回旋區大小和高度以本實驗室所進行的COREX氣化爐的熱態試驗為依據［15，18］.出口壓力設置為0.45 MPa，模擬未考慮爐內反應，只考慮入口氣體速度、成分等條件.

3　結果分析與討論

3.1不同加焦方式下的熱態物理模擬結果

采用58#石蠟、風量14m3·h－1、風溫100℃、料層高度370 mm、排料速度為0.5 r·min－1、爐內礦/(塊煤+焦炭)體積比0.5∶1、加焦粒度為5 mm條件下，對比爐料混裝、中心加焦8.7%、中間加焦8.7%三種情況下，模型內部各區域溫度的變化情況，分析加焦方式對熔化氣化爐內部氣流分布的影響，物理實驗結果如圖4所示.

對比分析圖4( a)與圖4( b)的物理實驗結果，可知采用中心加焦方式，爐內中心溫度即平行于縱向的x =0豎直線上溫度分布較爐內無中心焦柱時溫度分布稍高.說明爐內中心加焦使得中心料柱透氣性較好，氣流順行，故氣流溫度較高.同時中心加焦方式使得爐子邊緣處等溫線傾斜度較大，說明中心加焦對邊緣氣流有一定的抑制作用.

對比分析圖4( a)與圖4( c)可知，中間區域加焦時，與混裝條件下相比而言，在平行于縱軸的x =8 cm豎直線上，同一高度處，中間區域加焦方式下的熱電偶溫度值高，這是由于中間區域料柱空隙度增大，透氣性得到改善，煤氣流往中間發展，從而使中間區域溫度升高.中間焦柱的加入，使爐體中心區域的溫度也降低.這說明中間區域加焦氣流在爐體中間得到發展，邊緣及中心氣流得到抑制，氣流在爐內呈兩道氣流向上發展.

圖4　加焦方式對爐內溫度分布的影響Fig.4 Influence of different coke charging on temperature distribution in the furnace( a)—無中心焦柱; ( b)—中心加焦8.7%( 5 mm) ; ( c)—中間加焦8.7%( 5 mm)

3.2加焦方式對空隙度的影響

因物理實驗無法獲得顆粒堆積的內部信息，故采用DEM方法獲得了中心加焦、中間區域加焦時的氣化爐內的空隙度分布的結果，如圖5、圖6所示.由圖5可以看出，中心加焦使得氣化爐中心空隙度數值較大，而采用混裝布料的氣化爐兩側部分空隙度較小，空隙度數值的大小從側面可以說明煤氣流的走向，即氣流會沿著空隙度大的位置流動，空隙度小的位置會對煤氣流產生一定的抑制作用.由圖6可以看出，由于中間區域焦炭的加入使得熔化氣化爐兩側中部空隙度數值較大，而采用混裝布料的熔化氣化爐其他位置空隙度數值較小，從側面說明其中心部分透氣透液性將會變差，而氣化爐兩側中間區域透氣性將會得到改善.

對比分析中心加焦與中間加焦形成的料柱空隙度模擬結果，可獲得不同的加焦方式會使得熔化氣化爐內物料空隙度分布存在不同程度的差異，這勢必會造成爐內煤氣流分布的不一致，因此，針對熔化氣化爐采用的混裝布料模式以及其日臻成熟的研究手段，下文將對加焦方式下熔化氣化爐內煤氣流分布進行研究，詳細分析不同加焦方式下爐內煤氣流分布趨勢，以期提出較為合理的加焦方式.

圖5　中心加焦方式下空隙度分布圖Fig.5 Voidage distribution under central coke charing

3.3加焦方式對煤氣流分布的影響

根據上述離散單元模型的模擬結果，將各區域內空隙度分布以UDF形式導入到流體力學計算軟件的計算區域網格內，并對每一個網格內的空隙度進行UDF人工編程設定，不同加焦方式下其各區域的空隙度分布值列于表2中.

圖6　中間區域加焦方式下空隙度分布圖Fig.6 Voidage distribution under intermediate coke charing

表2　不同加焦方式下各區域空隙度分布Table 2 Voidage distribution of each zone under different coke charging patterns

根據流體力學多孔介質模型建立氣化爐煤氣流分布數學模型，模型獲得了不同加焦方式下氣化爐內氣體速度矢量、流線分布信息，如圖7所示.圖中可以看出，氣體從風口處鼓入氣化爐內，爐料的存在使得其對氣流有一定的阻礙作用，從而使氣流速度減小，并且其流向由于軟熔區域的存在而發生變化，究其原因是爐內各區域空隙度分布不同所致.焦炭在爐中心位置和爐中間位置的加入，會使得爐內空隙度分布存在明顯不同，中心加焦和中間區域加焦時，在爐中心位置和中間區域位置的空隙度較大，與未加入焦柱相比，氣體在此兩位置的流動不會出現明顯的方向變化，如圖7中氣體流線所示.

從圖7氣體速度矢量圖和流線圖還可以看出，同一加焦方式下，加焦位置處氣流速度要高于未加焦炭時的氣流速度，即加焦位置處會形成一個阻力較小的煤氣流發展通路，從側面說明加入焦柱后會對氣化爐透氣性有一定的改善，這與物理實驗得出的結論類似.

模型還可進一步獲得不同加焦方式下爐內氣體在徑向上及z = 0截面處的軸向上的速度分布信息，模擬結果如圖8所示.圖中可以看出在在不同加焦方式下，氣體流速在不同高度上的徑向分布以及z =0截面處不同徑向爐內速度在軸向上的分布趨勢大體是一致的，在氣化爐內加入焦柱后，煤氣流速整體得到增加，在氣化爐內部煤氣流速平均增加了約0.4 m·s－1，加焦位置處煤氣流速增加較多，約為1.1～1.6 m·s－1，如圖8中黑色橢圓線框所示.雖然焦炭的引入使得煤氣流速得到了提高，爐內透氣性得到明顯改善，但是加焦范圍的增加難免會影響煤氣利用率，因此，應嚴格控制加焦量.

4　結論

本文建立了研究三維氣化爐爐料結構和煤氣流分布的物理模型和數學模型.物理模型采用熱電偶測溫的方式，從爐內氣體溫度分布信息考察了氣體的流動情況，但物理實驗無法獲得內部信息，故基于離散單元法模型，以Fluent軟件為載體，利用多孔介質模型并加入用戶自定義函數，通過數學模型進一步研究了不同加焦方式下的氣化爐內煤氣流的分布的影響機理，獲得了氣化爐內煤氣的速度場和流線.其模擬結果對COREX氣化爐采用加焦措施指導調節煤氣流分布具有借鑒意義，其主要結論如下.

圖8　不同加焦方式爐內煤氣速度的徑向(左列)及軸向(右列)分布圖Fig.8 The gas velocity distribution in radial ( left) and axial ( right) directions under different coke charging patterns( a)、( d)—混裝(未加焦柱) ; ( b)、( e)—中心加焦; ( c)、( f)—中間區域加焦

( 1)利用DEM方法計算得到的不同加焦方式下爐內空隙度分布，對氣化爐內空隙度分布進行了設置，為今后計算氣化爐氣流分布空隙度設置提供了依據.

( 2)在氣化爐內加入焦柱后，煤氣流速整體得到增加，在氣化爐內部煤氣流速平均增加了約0.4 m·s－1，加焦位置處煤氣流速增加較多，約為1.1～1.6 m·s－1.

( 3)氣化爐內焦柱的加入會使得加焦位置形成煤氣流發展通路，進而改善氣化爐透氣性，但也應控制焦炭的加入量，避免氣流過度發展，影響煤氣流利用率.

參考文獻:

［1］周渝生，錢暉.COREX熔融還原煉鐵新工藝［J］.世界鋼鐵，2005，5( 1) : 22－26.( Zhou Yusheng，Qian Hui.COREX smelting reduction iron －making process［J］.World Iron＆Steel.2005，5( 1) : 22 －26.)

［2］方覺.非高爐煉鐵工藝與理論［M］.北京:冶金工業出版社，2010: 8－9.( Fang Jue.The technology and theory of non－blast furnace ironmaking［M］.Beijing: Metallurgy Industry Press，2010: 8 －9.)

［3］周渝生.煤基熔融還原煉鐵新工藝發展現狀評述［C］/ /中國非高爐煉鐵年會.沈陽，2006: 42－57.( Zhou Yusheng.The review of current development of doal －based smelting reduction Iron－making process［C］/ /China Non－Blast Furnace Iron－making Conference.Shenyang，2006: 42－57.)

［4］李維國.COREX－3000生產現狀和存在問題的分析［J］.寶鋼技術，2008( 6) : 11－18.( Li Weiguo.Operation status quo and technical problems of COREX－3000［J］.Baosteel Technology，2008 ( 6) : 11 －18.)

［5］畢學工.高爐過程數學模型及計算機控制［M］.北京:冶金工業出版社，1996.( Bi Xuegong.The mathematical model and computer control during the blast furnace process［M］.Beijing: Metallurgy Industry Press，1996.)

［6］儲滿生，郭憲臻，沈峰滿，等.高爐數學模型的進展［J］.中國冶金，2007，17( 4) : 10－14.( Chu Manshen，Guo Xianzhen，Shen Fengman，et al.Development of blast furnace mathematical model［J］.China Metallurgy，2007，17( 4) : 10－14.)

［7］于海彬，陳義勝.軟熔帶對煤氣流動影響的數值模擬［J］.包頭鋼鐵學院學報，2006，25( 1) : 1－4.( Yu Haibin，Chen Yisheng.Numerical simulation of influence on gas flow caused by cohesive zone［J］.Journal of Baotou University of Iron and Steel Technology，2006，25 ( 1) : 1 －4.)

［8］陳世超，李文忠，杜榮山，等.高爐軟熔帶氣體力學模型試驗研究［J］.鋼鐵，1981，16( 07) : 27－38.( Chen Shichao Li Wenzhong，Du Rongshan，et al.Gas pneumatic model study at cohesive zone of the blast furnace ［J］.Iron and Steel，1981，16 ( 07) : 27－38.)

［9］武文華，宋建成，賈陰隆.高爐軟熔帶氣體力學冷模擬研究［J］.鋼鐵研究，1993，( 06) : 38－42.( Wu Wenhua，Song Jiancheng，Jia Yinlong.Aerodynamic cold model study on fusion zone in BF［J］.Research on Iron ＆Steel，1993，( 06) : 38－42.)

［10］杜鶴桂，車傳仁.高爐軟熔帶對氣流分布影響的研究［J］.東北工學院學報，1980，( 01) : 72－84.( Du Hegui，Che Chuanren.On the influence of softeningmelting zone on gas distribution in blast furnace［J］.Journal of Northeastern University of Technology，1980，( 01 ) : 72 －84.)

［11］湯清華，安云沛，成蘭伯，等.高爐軟熔帶的形成及其對煤氣流分布的影響［J］.鋼鐵，1987，22( 06) : 4－10.( Tang Qinghua，AnYunPei，Cheng Lanbo，et al.Formation of cohesive zone in blast furnace and its influence on gas flow ［J］.Iron and Steel，1987，22( 06) : 4－10.)

［12］Cundall P A，Strack O D L.A discrete numerical model for granular assemblies［J］.Geotechnique，1979，29 ( 1) : 47 －65.

［13］Zhu H P，Zhou Z Y，Yang R Y，et al.Discrete particle simulation of particulate systems: theoretical developments ［J］.Chemical Engineering Science，2007，62: 3378－3396.

［14］Wright B，Zulli P，Zhou Z Y，et al.Gas－solid flow in an ironmaking blast furnace—I: physical modeling［J］.Powder Technology，2011，208( 1) : 86－97.

［15］周曉雷，韓立浩，王祺，等.利用溫度梯度分析COREX熔化氣化爐內區域［J］.東北大學學報(自然科學版)，2011，32( 8) : 1137－1140.( Zhou Xiaolei，Han Lihao，Wang Qi，et al.Analysis of zones in COREX melter gasifier by temperature gradient［J］.Journal of Northeastern University ( Natural Science)，2011，32( 8) : 1137－1140.)

［16］王成善，常青，冉登高，等.填充床中氣體流動特性的數值研究［J］.金屬材料與冶金工程，2010，38( 2) : 28－33.( Wang Chengshan，Chang Qing，Rab Denggao，et al.Numerical research on the characters of gas flow in packed bed ［J］.Metal Materials and Metallurgy Engineering，2010，38 ( 2) : 28－33.)

［17］李曉清，李海峰，羅志國.COREX熔化氣化爐布料模式對煤氣流分布的影響［J］.材料與冶金學報，2012，11( 4) : 235－241.( Li Xiaoqing，Li Haifeng，Luo Zhiguo.Influence of burden distribution mode on the gas flow distribution in COREX melter gasifier［J］.Journal of Materials and Metallurgy，2012，11( 4) : 235－241.)

［18］孫俊杰，狄瞻霞，韓立浩，等.熔化氣化爐風口回旋區的形成研究［C］/ /第十六屆全國冶金反應工程學術會.唐山，2012: 123－128.( Sun Junjie，Di Zhanxia，Han Lihao，et al.Study on the existence of raceway in the COREX melter gasifier［C］/ / Metallurty Reaction Engineering Conference.Tangshan，2012: 123－128.)

Influence of coke charging patterns on the gas flow distribution in COREX melter gasifier

Li Haifeng，You Yang，Han Lihao，Luo Zhiguo，Cai Jiuju，Zou Zongshu
( School of Materials＆Metallurgy，Northeastern University，Shenyang 110819，China)

Abstract:Charging mode determines the burden voidage which directly affects the second distribution of the gas flow.In this paper，a 3D mathematical model and physical model for gas flow was established to simulate the gas flow distribution under different coke charging patterns in COREX melter gasifier.An experiment for measuring the gas temperature with a thermocouple was carried out，which can be used to analyse the distribution information of the gas flow.The physical model was unable to gain the internal porosity.Based on the discrete element method ( DEM)，Fluent software and porous media method，Gas velocity vectors，streamlines were obtained.The results showed that the voidage distribution caused by coke charging patterns in the charging process has a significant effect on the gas flow distribution.The simulation results are in good agreement with those obtained by experiment.The authors believe that the gas flow paths above are available for coke charging patterns in COREX melter gasifier.

Key words:COREX melter gasifier; gas flow distribution; coke charging patterns; physical simulation; numerical simulation

作者簡介:李海峰( 1982—)，男，博士，E－mail: lihf@ smm.neu.edu.cn.

基金項目:國家科技支撐計劃項目( 2011BAE04B01，2011BAE04B02) ;國家自然科學基金資助項目( 51104037).

收稿日期:2015-05-11.

doi:10.14186/j.cnki.1671－6620.2016.01.003

中圖分類號:TF 557

文獻標識碼:A

文章編號:1671-6620( 2016) 01-0012-08