張麗娟，周 恒，張 濤，游 陽，羅志國，鄒宗樹

(1.東北大學(xué) 材料與冶金學(xué)院，沈陽 110819;2.寶山鋼鐵股份有限公司人才開發(fā)院，上海 200941)

COREX熔融還原煉鐵工藝是奧鋼聯(lián)開發(fā)的一種用煤和球團(tuán)(塊礦)生產(chǎn)鐵水的非高爐煉鐵工藝，具有非焦煤冶煉，環(huán)境友好的特點(diǎn)［1，2］.與傳統(tǒng)高爐相比，COREX預(yù)還原豎爐鼓風(fēng)條件、布料形式、爐型結(jié)構(gòu)等不盡相同，所以爐內(nèi)煤氣停留時(shí)間分布必然呈現(xiàn)出與高爐不同的特征.雖然高爐內(nèi)煤氣流動(dòng)及停留時(shí)間分布研究較多［3～6］，但關(guān)于COREX預(yù)還原豎爐的研究較少.寶鋼引進(jìn)的COREX－3000預(yù)還原豎爐的平均直徑由COREX－2000的6.67 m增加至8.2 m，反應(yīng)器尺寸的加大使?fàn)t內(nèi)煤氣流動(dòng)愈偏離理想流動(dòng)狀態(tài).2#COREX－3000預(yù)還原豎爐在原1#豎爐基礎(chǔ)上安裝了AGD梁，利用爐料下行運(yùn)動(dòng)過程中在梁下方自發(fā)形成的無爐料煤氣通道，將煤氣導(dǎo)入豎爐中心.AGD梁的安裝進(jìn)一步改變了煤氣的初始分布，直接影響了爐內(nèi)煤氣停留時(shí)間分布.因此，有必要對(duì)COREX預(yù)還原豎爐中煤氣停留時(shí)間分布進(jìn)行研究.

本文建立三維、不可壓縮、穩(wěn)態(tài)下還原煤氣在COREX預(yù)還原豎爐內(nèi)停留時(shí)間分布的數(shù)學(xué)模型，考察AGD梁對(duì)爐內(nèi)RTD曲線的影響，并研究不同熔煉率下COREX預(yù)還原豎爐內(nèi)的RTD變化情況，同時(shí)給出不同條件下爐內(nèi)流動(dòng)特性參數(shù)變化，為指導(dǎo)COREX豎爐實(shí)際生產(chǎn)提供理論依據(jù).

1 數(shù)學(xué)模型建立

1.1 控制方程

COREX預(yù)還原豎爐為氣－固逆流反應(yīng)器，相對(duì)于氣體流速而言，固體爐料的下降速度很小，故可將爐內(nèi)爐料近似為固定床，并且忽略氣－固化學(xué)反應(yīng)對(duì)氣體流場的影響.描述煤氣流動(dòng)的基本方程包括連續(xù)性方程、動(dòng)量方程、湍流模型及濃度方程.

1.2 邊界條件及算法

COREX熔化氣化爐產(chǎn)生的還原煤氣經(jīng)冷卻、除塵后通過豎爐圍管導(dǎo)入爐內(nèi)，但部分高溫煤氣通過氣化爐穹頂?shù)腄RI螺旋下降管反竄進(jìn)入豎爐，形成反竄煤氣.因此，將預(yù)還原豎爐的煤氣入口設(shè)為壓力入口，其中圍管煤氣入口表壓為360 kPa，DRI下降管煤氣入口表壓為365 kPa.豎爐頂部煤氣出口為速度出口.COREX預(yù)還原豎爐爐型結(jié)構(gòu)及網(wǎng)格劃分如圖1所示，兩者最大區(qū)別是圍管區(qū)AGD梁的安裝以及頂部煤氣出口的改變.方程求解借助商業(yè)軟件FLUENT，壓力－速度耦合采用SIMPLE算法，控制方程的離散采用一階迎風(fēng)模式，收斂標(biāo)準(zhǔn)殘差小于10－4.

1.3 RTD曲線分析方法

COREX預(yù)還原豎爐內(nèi)還原煤氣的流動(dòng)模式為偏離活塞流及全混流這兩種流動(dòng)狀態(tài)的非理想流動(dòng)狀態(tài)，其中，1#COREX預(yù)還原豎爐頂部煤氣出口為單一出口，因此，COREX預(yù)還原豎爐RTD曲線分析模型可借鑒中間包中流體流動(dòng)特性分析方法［7～10］，本文采用經(jīng)典組合方法.

2#COREX－3000預(yù)還原豎爐在1#爐基礎(chǔ)上進(jìn)行了改裝，其頂部煤氣出口為兩個(gè)，因此2#COREX預(yù)還原豎爐流體流動(dòng)特性描述方法借鑒多流中間包總體RTD曲線分析方法［11，12］.

2 模擬結(jié)果與討論

2.1 AGD梁對(duì)RTD曲線的影響

圖1 網(wǎng)格劃分示意圖Fig.1 Schematic diagram of the grid

圖2 1#及2#COREX預(yù)還原豎爐停留時(shí)間分布密度函數(shù)Fig.2 Density of residence time distribution of 1#and 2#COREX shaft furnace

圖2為COREX－3000預(yù)還原豎爐頂煤氣(標(biāo)準(zhǔn)態(tài))單耗1 050 m3/t，熔煉率150 t/h條件下，安裝AGD前后，還原煤氣停留時(shí)間分布密度函數(shù).由圖2可知，豎爐內(nèi)的煤氣流型為活塞流、全混流和死區(qū)的組合流動(dòng)，且出現(xiàn)兩個(gè)示蹤劑濃度的峰值，分別為圍管煤氣和反竄煤氣內(nèi)的示蹤劑到達(dá)監(jiān)測面所致.隨著AGD梁的安裝，第一峰高度增高，分布變窄，第二峰高度降低且分布變寬，曲線整體寬度變寬.

其中，第二峰為豎爐底部下料管反竄煤氣在檢測面的信息，其主流微元通過流體系統(tǒng)的時(shí)間增長.2#COREX預(yù)還原豎爐AGD梁的安裝能抑制底部下料管煤氣反竄比例［13，14］，在頂部煤氣流量不變條件下，其底部反竄煤氣入口速度將變小;因此，第二峰出峰時(shí)間較1#COREX預(yù)還原豎爐長.2#COREX預(yù)還原豎爐煤氣停留時(shí)間分布描述方法借鑒多流中間包總體RTD曲線分析方法.由圖可知，2#COREX豎爐總體RTD曲線和各出口RTD曲線形態(tài)相差較小，主要原因是2#COREX豎爐頂部兩煤氣出口為對(duì)稱布置，檢測得到的微元停留時(shí)間分布相一致.

表1為COREX－3000預(yù)還原豎爐加裝AGD梁前后爐內(nèi)流動(dòng)特征參數(shù).由表可知，豎爐安裝AGD梁后其煤氣平均停留時(shí)間減小，但減小幅度較小.究其原因是，AGD梁的安裝使?fàn)t內(nèi)圍管煤氣比例增大，圍管水平煤氣平均速度增大，因此第一峰的出峰時(shí)間較早.但整體考慮豎爐底部反竄煤氣影響，其煤氣平均停留時(shí)間減小幅度較小.從表還可以看出，隨著AGD梁的安裝，無量綱方差σ2變大.無量綱方差是用于度量隨機(jī)變量與均值的偏離程度，它與流型一一對(duì)應(yīng)，在活塞流的情況下，方差最小;而在全混流的情況下，方差最大.因此，COREX預(yù)還原豎爐內(nèi)的煤氣流動(dòng)為非理想流動(dòng)，且隨著AGD梁的安裝，無量綱方差變大，與活塞流的偏差越大.主要原因是圍管處還原煤氣可通過AGD梁下方煤氣通道進(jìn)入豎爐中心，增加了豎爐內(nèi)煤氣的反混現(xiàn)象，從而使流型偏離理想活塞流.2#COREX豎爐各出口表征的爐內(nèi)流動(dòng)特征參數(shù)與總體方法得出的流動(dòng)特征參數(shù)相差不大，這主要是由于豎爐頂部兩煤氣出口的對(duì)稱布置所致.

表1 AGD梁對(duì)爐內(nèi)流動(dòng)特征參數(shù)的影響Table 1 Effect of AGD beams on the flow characteristics parameters in furnace

2.2 頂煤氣單耗對(duì)RTD的影響

圖3為不同熔煉率下COREX預(yù)還原豎爐煤氣停留時(shí)間分布密度函數(shù).圖4、圖5分別為平均停留時(shí)間、無量綱方差隨熔煉率的變化圖.從圖中可以看出隨著熔煉率的升高，峰值高度逐漸升高，平均停留時(shí)間逐漸降低，無量綱方差逐漸增大.導(dǎo)致上述結(jié)果原因有:一是在COREX預(yù)還原豎爐內(nèi)還原煤氣主流微元在爐內(nèi)滯留時(shí)間相應(yīng)減小，所以平均停留時(shí)間減小;二是隨著熔煉率的升高，煤氣入口氣量增大，煤氣出口的動(dòng)量增大，爐內(nèi)煤氣的混合程度變大，流型愈偏離理想活塞流，所以無量綱方差逐漸增大.

圖6為COREX預(yù)還原豎爐內(nèi)死區(qū)體積分?jǐn)?shù)隨熔煉率的變化圖.死區(qū)指的是在此區(qū)域內(nèi)的流體運(yùn)動(dòng)得相當(dāng)慢，結(jié)果流體在容器中停留了相當(dāng)長的時(shí)間.為描述不同熔煉率下爐內(nèi)死區(qū)體積分?jǐn)?shù)，本文以150 t/h熔煉率對(duì)應(yīng)的平均停留時(shí)間的兩倍為基準(zhǔn).由圖可知，隨著熔煉率的升高，爐內(nèi)死區(qū)體積分?jǐn)?shù)逐漸降低.其中，1#COREX預(yù)還原豎爐降低幅度較大，當(dāng)熔煉率從150 t/h提升到180 t/h時(shí)，爐內(nèi)死區(qū)體積分?jǐn)?shù)從10.9%降到3.48%，降幅為68.1%.2#COREX預(yù)還原豎爐死區(qū)體積分?jǐn)?shù)降幅較小，僅為15.52%.主要原因是由于2#COREX預(yù)還原豎爐AGD梁的安裝，改變了煤氣初始分布，大部分還原煤氣通過AGD梁下方通道進(jìn)入豎爐.增大熔煉率時(shí)，導(dǎo)氣槽出口還原煤氣速度增量較1#COREX預(yù)還原豎爐低，因此對(duì)豎爐下部死區(qū)體積影響較小，COREX豎爐死區(qū)體積分?jǐn)?shù)降幅較小.

3 結(jié)論

(1)基于中間包流體流動(dòng)特性經(jīng)典組合方法分析得到1#COREX預(yù)還原豎爐的RTD分布，并借鑒多流中間包總體分析法得到2#COREX預(yù)還原豎爐內(nèi)的流動(dòng)特性參數(shù).

(2)隨著COREX－3000預(yù)還原豎爐內(nèi)AGD梁的安裝，爐內(nèi)煤氣平均停留時(shí)間變小，停留時(shí)間分布密度函數(shù)變寬，無量綱方差變大，流型偏離理想活塞流.同時(shí)，AGD梁的安裝，增大了爐內(nèi)死區(qū)體積分?jǐn)?shù).

圖3 不同熔煉率下停留時(shí)間分布密度函數(shù)Fig.3 Density of residence time distribution with different melting rate

圖4 不同熔煉率下平均停留時(shí)間分布Fig.4 Average residence time with different melting rate

(3)隨著COREX預(yù)還原豎爐熔煉率的增大，煤氣平均停留時(shí)間減小，無量綱方差變大.當(dāng)熔煉率從150 t/h增長到180 t/h時(shí)，爐內(nèi)死區(qū)體積分?jǐn)?shù)降低，其中1#COREX豎爐降幅為68.1%，2#COREX豎爐降幅僅為15.52%.

圖5 不同熔煉率下無量綱方差分布Fig.5 Variance of residence time distribution

圖6 不同熔煉率下死區(qū)體積分?jǐn)?shù)Fig.6 Volume fraction of dead zone with different melting rate

［1］ Qu Yingxia，Zou Zongshu，Xiao Yanping.A comprehensive static model for COREX process［J］.ISIJ International，2012，52(12):2186－2193.

［2］周渝生.煤基熔融還原煉鐵新工藝開發(fā)現(xiàn)狀評(píng)述［J］.鋼鐵，2005，40(11):1－8.

(Zhou Y S.Review of current development of coal－based smelting reduction ironmaking process［J］.Iron and Steel，2005，40(11):1 －8.)

［3］ Yagi J，Takeda K，Omori Y.Two－dimensional simulation on the gas flow and heat transfer in the blast furnace［J］.ISIJ International，1982，22(11):884 －892.

［4］ Jimenez J，Mochon J，Ayala J S.Mathmatical model of gas flow distribution in a scale model of a blast furnace shaft［J］.ISIJ Int，2004(3):518 －526.

［5］孫紹杰.用氦氣示蹤法分析高爐爐況［J］.鞍鋼技術(shù)，1982(12):659－660.

(Sun S J.Analysis of blast furnace condition by helium tracing method［J］.Angang Technology，1982(12):659 －660.)

［6］孫紹杰，王喜來.利用85Kr示蹤劑對(duì)高爐內(nèi)煤氣運(yùn)動(dòng)的初步研究［J］.煉鐵，1984(1):53－58.

(Sun S J，Wang X L.Application of 85Kr tracer technique on the study of gas flow in blast furnace［J］.IRONMAKING，1984(1):53 －58.)

［7］ Liu S，Yang X，Du L，et al.Hydrodynamic and mathematical simulationsofflow field and temperature profile n an asymmetrical T－type single－strand continuous casting tundish［J］.ISIJ International，2008，48(12):1712 －1721.

［8］ Sahai Y，Emi T.Melt flow characterization in continuous casting tundishes［J］.ISIJ International，1996，36(6):667－672.

［9］Mazumdar D，Yamanoglu G，Guthrie R I L.Hydrodynamic performance of steelmaking tundish systems:a comparative study of three different tundish designs［J］.Steel Research，1997，68(7):293 －300.

［10］ Kumar A，Koria S C，Mazumdar D.Basis for systematic hydrodynamic analysis of a multi－ strand tundish［J］.ISIJ International，2007，47(11):1618 －1624.

［11］雷洪，趙巖，鮑家琳，等.多流連鑄中間包停留時(shí)間分布曲線總體分析方法［J］.金屬學(xué)報(bào)，2010，46(9):1109－1114.

(Lei H，Zhao Y，Bao J X，et al.Whole analysis approach for residue time distribution curve in multi－strand continuous casting tundish［J］.Acta Metallurgica Sinica，2010，46(9):1109 －1114.)

［12］雷洪，趙巖，邢國成，等.總體分析法在非對(duì)稱兩流中間包水模型中的應(yīng)用［J］.東北大學(xué)學(xué)報(bào) (自然科學(xué)版)，2011，32(4):537－540.

(Lei H，Zhao Y，Xing G C，et al.Use of a comprehensive analytical approach for water modeling of an asymmetrical two－ strand tundish［J］.Journal of Northeastern University(Natural Science)，2011，32(4):537－540.)

［13］徐輝，李曉清，林金嘉，等.COREX豎爐內(nèi)煤氣流分布的三維數(shù)值模擬［J］.世界鋼鐵，2012，1:15－19.

(Xu H，Li X Q，Lin J J，et al.3D numerical simulation of gas flow distribution in COREX shaft furnace［J］.World Iron ＆Steel，2012，1:15 －19.)

［14］ Zhou H，Luo Z G，Zhang T，et al.3D numerical simulation of the influence of AGD beams on gas distribution in COREX－3000 shaft furnace［J］.Advanced Materials Research.2013，712－715:1268－1273.