盧紹文，丁進(jìn)良，劉琳琳

（東北大學(xué)流程工業(yè)綜合自動化國家重點(diǎn)實(shí)驗室，遼寧 沈陽 110819）

引 言

我國鐵礦資源主要以貧礦為主，品位低，選別困難。一般需要通過豎爐對鐵礦石進(jìn)行高溫還原磁化焙燒，使弱磁性礦物變成強(qiáng)磁性礦物，然后經(jīng)過磁選得到鐵精礦[1]。豎爐過程是包含復(fù)雜的物理、化學(xué)變化的工業(yè)過程，由于給礦、給氣物料運(yùn)動方式的特殊性，無法實(shí)時地檢測關(guān)鍵過程質(zhì)量數(shù)據(jù)。這些難測過程變量難以控制，對生產(chǎn)過程運(yùn)行品質(zhì)的影響難以定量和定性評估。而且豎爐過程是高耗能過程，面臨節(jié)能降耗的重要挑戰(zhàn)。因此，如何使豎爐過程在優(yōu)化狀態(tài)下穩(wěn)定運(yùn)行，是相關(guān)企業(yè)關(guān)心的熱點(diǎn)問題[2]。

豎爐還原焙燒過程的主要工藝指標(biāo)是磁選管回收率，它反映了礦石焙燒的質(zhì)量與效率。磁選管回收率需要滿足一定的工藝要求才能使得后續(xù)的磁選操作選出品位合格的精礦[3]。磁選管回收率無法在線連續(xù)檢測，且與豎爐主要操作變量，如燃燒室溫度、還原煤氣流量、礦石搬出時間等之間的動態(tài)特性難以用精確模型描述，表現(xiàn)出強(qiáng)非線性、強(qiáng)耦合、隨工況運(yùn)行條件變化而變化的復(fù)雜特性。許多企業(yè)仍然依靠操作員按照生產(chǎn)邊界條件，憑經(jīng)驗給出豎爐底層控制回路設(shè)定值，來控制磁選管回收率，使其進(jìn)入工藝目標(biāo)區(qū)間。但當(dāng)邊界條件發(fā)生變化時，操作員常常不能及時準(zhǔn)確地調(diào)整設(shè)定值，使得磁選管回收率不能滿足工藝要求。

豎爐焙燒的運(yùn)行優(yōu)化控制技術(shù)可以很好地解決上述問題[2,4-5]。該技術(shù)在回路控制上層增加了優(yōu)化設(shè)定控制層，通過控制回路預(yù)設(shè)定模型、前饋補(bǔ)償器、反饋補(bǔ)償器、磁選管回收率預(yù)報模型，動態(tài)地給出豎爐過程底層控制回路設(shè)定值。確保邊界條件變化下產(chǎn)品質(zhì)量和生產(chǎn)運(yùn)行穩(wěn)定，同時使得磁選管回收率滿足工藝要求。

豎爐焙燒過程運(yùn)行優(yōu)化控制方法設(shè)計的關(guān)鍵是建立磁選管回收率相對于控制操作參數(shù)、邊界條件的動態(tài)響應(yīng)關(guān)系。由于精確數(shù)學(xué)模型難建，需要進(jìn)行大量工業(yè)現(xiàn)場實(shí)驗，獲得磁選管回收率的動態(tài)響應(yīng)數(shù)據(jù)。工業(yè)實(shí)驗可能導(dǎo)致生產(chǎn)偏離正常工藝范圍，進(jìn)而產(chǎn)生不合格產(chǎn)品；而且采樣和化驗分析的頻率比正常運(yùn)行更頻繁。采樣動態(tài)模擬技術(shù)，進(jìn)行虛擬實(shí)驗，能夠極大節(jié)約運(yùn)行優(yōu)化控制方法設(shè)計的實(shí)驗成本，降低風(fēng)險。

目前常用的流程模擬軟件，如Aspen Plus、Pro/Ⅱ、ChemCAD及HYSYS大多用于化工過程，但仍缺乏用于冶金過程的模型單元。METSIM是冶金領(lǐng)域中較成熟的軟件包[6-8]。由許多模塊組成。這些模塊含有設(shè)計規(guī)范的方程組子集和每一步工藝執(zhí)行特性，允許對流程中的每一個操作單元進(jìn)行單獨(dú)分析。METSIM能夠模擬冶金和選礦過程中常用單元操作，內(nèi)嵌物性數(shù)據(jù)庫，采用序貫?zāi)K法進(jìn)行物料和能量平衡。但目前報道的應(yīng)用主要是冶金工藝的穩(wěn)態(tài)模擬。實(shí)際上，METSIM具備動態(tài)模擬的功能，完全能夠用于先進(jìn)控制、優(yōu)化控制方法的仿真實(shí)驗，但這方面的應(yīng)用仍未見報道。為此，本文采用METSIM軟件建立豎爐焙燒過程的動態(tài)模擬對象，進(jìn)而將其與真實(shí)的豎爐焙燒過程控制系統(tǒng)和運(yùn)行優(yōu)化控制系統(tǒng)連接，建立半實(shí)物仿真系統(tǒng)，為運(yùn)行優(yōu)化控制方法提供檢驗、測試工具。

1 豎爐焙燒還原過程工藝描述

豎爐還原焙燒過程將原礦石在豎爐加熱帶加熱到 700～800℃以后，礦石自動落入還原帶，在550～600℃溫度下與煤氣中的還原氣體發(fā)生化學(xué)反應(yīng)，從而將弱磁性的赤鐵礦還原為強(qiáng)磁性的磁鐵礦[9]，其工藝流程如圖1所示。

豎爐過程生產(chǎn)系統(tǒng)從上至下可以分為給礦系統(tǒng)、燃燒系統(tǒng)、高溫還原系統(tǒng)、冷卻搬出系統(tǒng)4個子系統(tǒng)。其中給礦系統(tǒng)主要包括礦石的輸送和預(yù)熱兩個過程，通過皮帶輸送機(jī)將原礦石輸送到爐頂存儲槽，礦石先到達(dá)爐膛上部的預(yù)熱帶，與加熱帶出氣口的氣體發(fā)生熱交換，使其平均溫度達(dá)到 150～200℃；燃燒系統(tǒng)包括燃燒和加熱兩個過程，預(yù)熱的礦石進(jìn)入加熱帶，通過助燃空氣和燃燒煤氣發(fā)生化學(xué)反應(yīng)所釋放的熱量，對礦石進(jìn)行加熱，同時為了保證礦石在進(jìn)入還原帶時有足夠高的溫度進(jìn)行還原反應(yīng)，必須使燃燒控制室內(nèi)的溫度控制在 1050～1150℃，這就要求煤氣與助燃空氣配比適當(dāng)，使其充分發(fā)生反應(yīng)，發(fā)揮最大的燃燒率；高溫還原系統(tǒng)是整個系統(tǒng)中最重要的環(huán)節(jié)，加熱的礦石進(jìn)入還原帶與還原煤氣相接觸，使其溫度降低在 570℃左右時與還原煤氣中的CO和H2進(jìn)行還原反應(yīng)；冷卻搬出系統(tǒng)將還原帶出口的礦石在水封池進(jìn)行冷卻至400℃，然后由搬出系統(tǒng)運(yùn)出至布料系統(tǒng)以備磨礦過程使用。

2 豎爐過程METSIM動態(tài)模擬

2.1 物料衡算

燃燒加熱和還原氣體物流走向如圖2所示。燃燒單元即燃燒煤氣中CO和H2等與空氣中的O2在一定的溫度、壓力下發(fā)生劇烈的氧化反應(yīng)并放熱。其中燃燒煤氣為高爐煤氣，其主要成分如表1所示。助燃空氣為普通空氣，其成分為 O2（21%）和 N2（79%）。還原單元主要是低溫的還原煤氣被高溫的礦石加熱后與礦石中的氧化鐵發(fā)生化學(xué)反應(yīng)，從而將弱磁性的Fe2O3還原為強(qiáng)磁性的Fe3O4。礦石磁性的提高有利于后續(xù)選別工序的選別，提高金屬回收率。

圖1 豎爐還原焙燒工藝流程及儀表Fig. 1 Process and instrument diagram of shaft furnace roasting process

圖2 豎爐焙燒中的物流走向Fig. 2 Material flows in shaft furnace roasting process

豎爐燃燒室內(nèi)的氣體進(jìn)行復(fù)雜的燃燒、流動及能量轉(zhuǎn)換過程，影響燃燒的因素不僅是化學(xué)反應(yīng)動力因素，還包括其他流動、傳熱、傳質(zhì)等因素[10]，若全面考慮這些因素，模型將會變得異常復(fù)雜，以至于不能適用于系統(tǒng)動態(tài)特性的研究，因此，需要在不失過程特性且具有一定精度的前提下對模型進(jìn)行簡化。在燃燒環(huán)節(jié)發(fā)生的燃燒反應(yīng)為

表1 高爐煤氣成分Table 1 Components of blast furnace gas/%

表2 焦?fàn)t煤氣成分Table 2 Components of coke oven gas/%

還原帶單元中還原煤氣為焦?fàn)t煤氣，其成分如表2所示，從表中可以看出，在反應(yīng)時組分CnHm、CO、H2、CH4等都將不同程度地參與化學(xué)反應(yīng)，但結(jié)合各組分的比例，只有甲烷和氫氣是有可能影響還原反應(yīng)，同時又考慮到甲烷在高溫狀態(tài)下裂解生成的氫氣[11]比例較小，所以本文只描述氫氣還原氧化鐵的過程，暫不考慮裂解部分，依據(jù)某赤鐵礦選礦廠提供的數(shù)據(jù)，鐵礦石的成分如表3所示。其中發(fā)生的化學(xué)反應(yīng)為[9]

表3 某選礦廠赤鐵礦鐵礦石成分Table 3 Component percentages of feeding hematite ore to shaft furnace

2.2 能量衡算

燃燒室內(nèi)的能量來源有兩部分：一部分是化學(xué)反應(yīng)放熱；另一部分是氣體自身帶入的熱量。燃燒室內(nèi)的能量流失也包括兩部分：一部分是高溫?zé)煔鈳ё叩臒崃?；另一部分是向爐壁傳熱照成的損耗。燃燒室內(nèi)的穩(wěn)態(tài)能量平衡方程為[12]

式中，為燃料全部燃燒時的發(fā)熱量，Hr為燃料在輸入系統(tǒng)時帶入的熱焓，Hk為燃燒時所用空氣在輸入系統(tǒng)時帶入的熱焓，Qlj為燃燒產(chǎn)物由于熱離解而吸收的熱量，Qcr為熱損失，Qwr為未燃燒產(chǎn)物的熱量，Hqy為煙氣在燃燒溫度時具有的熱焓，W為煙氣對外做功耗散的能量。

加熱帶內(nèi)主要進(jìn)行的是燃燒過程釋放的熱量，通過高溫氣體向低溫礦石的傳熱過程，其中主要的傳熱過程有高溫氣體對原礦石的對流傳熱，氣體對爐壁的對流傳熱，以及原礦石對爐壁的輻射傳熱。加熱帶的能量平衡方程為[13]

式中，ΣV為爐膛內(nèi)氣體總體積，m3；為爐膛內(nèi)氣體集總溫度，℃；Gi為加熱氣體在溫度下的水當(dāng)量，kJ·℃?1；為爐膛內(nèi)礦石的集總溫度，℃；t0為爐壁溫度，℃；Ui為爐膛內(nèi)礦石的總質(zhì)量，kg；Wi為爐膛內(nèi)礦石的水當(dāng)量，kJ·℃?1；hA為有效面積傳熱系數(shù)；δig為黑體輻射系數(shù)；為黑體輻射傳熱面積，m2。

還原帶物料的能量是以熱力學(xué)能的形式存儲的，它的變化率等于流入的原礦石帶入的熱量減去原礦石與還原氣體發(fā)生的對流換熱，高溫礦石與還原氣體的輻射換熱，礦石與爐壁的輻射換熱，發(fā)生的化學(xué)反應(yīng)的吸熱以及流出礦石所帶出的熱量，所以其能量平衡方程為

式中，Ms為礦石的總質(zhì)量，Cs為礦石比熱容，ws為固體質(zhì)量流量，Ts為還原帶礦石的溫度，為入口礦石的溫度，hsg為氣固對流傳熱系數(shù)，為對氣固對流換熱面積，Tg為還原帶氣體溫度，δsg為氣固輻射傳熱系數(shù)，為氣固輻射換熱面積，δso為固體與爐壁之間的輻射傳熱系數(shù)，為礦石與爐壁之間的輻射換熱面積，ΔH為化學(xué)反應(yīng)吸熱速率。

2.3 METSIM模型

采用METSIM對豎爐焙燒過程中的燃燒和還原單元建立如圖3所示模型。運(yùn)行豎爐METSIM模型，物料衡算和能量衡算結(jié)果分別如表4、表5所示。

圖3 豎爐焙燒METSIM模型結(jié)構(gòu)Fig. 3 Structure of shaft furnace roasting process in METSIM

表4 固相物料衡算結(jié)果Table 4 Results of mass balance (solid phase)

表5 能量衡算結(jié)果Table 5 Results of heat balance/ kJ·h?1

建模過程中選用METSIM模型庫中的SPP容器用來模擬礦石的給礦系統(tǒng)，在其輸入物流中配比鐵礦石的成分（成分見表3），同時考慮到原料礦石在爐膛入料口被加熱帶的出口氣體先進(jìn)行預(yù)熱達(dá)到平均溫度為150～200℃，將進(jìn)入加熱帶的原料礦石的初始溫度設(shè)置為 200℃，近而模擬了原礦石的預(yù)熱過程。用BRN容器來模擬燃燒室，在其Reactions屬性框中根據(jù)式（1）和式（2）配置 CO和 H2與O2的化學(xué)反應(yīng)方程式，這里設(shè)置為完全反應(yīng)（反應(yīng)率=100%）。為了模擬放熱氣體對進(jìn)入爐膛下落的礦石進(jìn)行加熱過程，選用BIN儲料罐模擬加熱帶，在其內(nèi)部高溫氣體與礦石進(jìn)行匯合加熱。

燃燒和還原環(huán)節(jié)是豎爐還原焙燒過程中的關(guān)鍵環(huán)節(jié)，采用METSIM進(jìn)行動態(tài)模擬的主要難點(diǎn)在于確定不同溫度下氫氣還原氧化鐵的反應(yīng)速率。根據(jù)低溫下氫氣還原氧化鐵的動力學(xué)研究，在正常工藝操作溫度范圍內(nèi)氫氣還原氧化鐵的速率R（%）與反應(yīng)溫度T（℃）近似為線性關(guān)系[14]

根據(jù)實(shí)驗數(shù)據(jù)進(jìn)行一元線性回歸，獲得系統(tǒng)分別為：a= ?0.6，b= ?0.002。選取 METSIM 中的 TNK 反應(yīng)器模型來模擬還原帶。根據(jù)之前的理論分析在Reactions屬性框中配置Fe2O3和H2的還原反應(yīng)方程式(3)。依據(jù)之前分析得到的一元線性回歸方程，將化學(xué)反應(yīng)率的計算模式設(shè)置為Enter Expression，輸入反應(yīng)率方程式(8)的APL方程0.3+0.002×(VTEC s7?450)，其中，VTEC s7為燃燒室進(jìn)入還原帶的氣流的溫度。

3 豎爐過程半實(shí)物仿真系統(tǒng)

豎爐焙燒過程的控制主要包括燃燒室溫度控制回路和還原煤氣流量控制回路。其中燃燒室溫度控制采用串級控制，外環(huán)為燃燒室溫度，內(nèi)環(huán)為燃燒煤氣流量和助燃空氣流量控制。還原煤氣流量設(shè)定值由操作員根據(jù)經(jīng)驗或優(yōu)化設(shè)定系統(tǒng)給出。在實(shí)際生產(chǎn)過程中，礦石的大小、種類、成分頻繁地變化，人工給定設(shè)定值通常不能及時地根據(jù)工況變化進(jìn)行調(diào)整，難以將磁選管回收率控制在目標(biāo)值范圍內(nèi)。為此引入豎爐過程運(yùn)行優(yōu)化控制。

豎爐過程運(yùn)行優(yōu)化控制是在運(yùn)行過程中根據(jù)期望的磁選管回收率指標(biāo)及工況變化自動調(diào)整底層控制回路設(shè)定值。通過底層控制系統(tǒng)的輸出跟蹤給定的設(shè)定值，從而將磁選管回收率控制在目標(biāo)值范圍內(nèi)。豎爐過程運(yùn)行優(yōu)化控制主要由采用案例推理、神經(jīng)網(wǎng)絡(luò)、專家系統(tǒng)、模糊控制等技術(shù)建立的控制回路優(yōu)化設(shè)定模型、爐況診斷模型、智能預(yù)報模型、前饋與反饋補(bǔ)償模型、多變量智能解耦控制系統(tǒng)組成，能夠在豎爐焙燒運(yùn)行過程中智能地給出燃燒室溫度、還原煤氣流量等控制回路的優(yōu)化設(shè)定值，使得各回路跟蹤優(yōu)化設(shè)定值，從而保證磁選管回收率進(jìn)入最優(yōu)目標(biāo)值范圍內(nèi)[2]。

圖4 豎爐焙燒過程運(yùn)行優(yōu)化控制半實(shí)物仿真系統(tǒng)Fig. 4 Hardware-in-loop simulation system of integrated control of shaft furnace roasting process

為了實(shí)現(xiàn)對豎爐焙燒過程運(yùn)行優(yōu)化控制系統(tǒng)接近工況實(shí)際的驗證研究，建立了半實(shí)物仿真實(shí)驗平臺，如圖4所示。該仿真系統(tǒng)主要由METSIM虛擬對象仿真系統(tǒng)、過程監(jiān)控系統(tǒng)、豎爐焙燒過程運(yùn)行優(yōu)化控制系統(tǒng)、虛擬儀表與執(zhí)行機(jī)構(gòu)裝置系統(tǒng)、羅克韋爾的RSLogix5000可編程邏輯控制器（PLC）組成。其中，組成部分除了虛擬對象、虛擬儀表與執(zhí)行機(jī)構(gòu)裝置，其余部分均與工業(yè)現(xiàn)場一致。

模擬豎爐焙燒過程的動態(tài)過程的 METSIM 動態(tài)模型建立過程如下：首先設(shè)置運(yùn)行周期為每毫秒執(zhí)行50次動態(tài)輸出，共采集1000個時間點(diǎn)數(shù)據(jù)。設(shè)置燃燒煤氣流量的控制變量 USSP1，并在METSIM的 DDE通道中配置 EXCEL標(biāo)簽名USSP1，設(shè)置數(shù)據(jù)的傳輸方向為INPUT。同時需要將過程輸出值也通過接口軟件輸出到EXCEL中，用于外部PLC控制器的PID計算，即燃燒室溫度TEMP3，設(shè)置數(shù)據(jù)的傳輸方向為OUTPUT；在外部PLC中設(shè)置PID控制器，并在標(biāo)簽庫中統(tǒng)一命名，如燃燒煤氣控制器為 rsmqPID，燃燒室溫度temp_out，燃燒煤氣流量USSTR1，主定時器timer；最后，在數(shù)據(jù)通信接口軟件中，將仿真過程所需要的PLC標(biāo)簽名預(yù)先寫入單元格中。

4 實(shí)驗驗證研究

以豎爐焙燒過程多變量智能優(yōu)化設(shè)定方法為驗證對象，依托本文介紹的METSIM模型和半實(shí)物仿真系統(tǒng)進(jìn)行虛擬實(shí)驗。

考慮到豎爐底層控制中只有加熱煤氣與燃燒室溫度形成閉環(huán)（通過加熱煤氣和形成空氣定比例控制），而還原煤氣流量和礦石搬出時間均為開環(huán)，通過運(yùn)行優(yōu)化控制系統(tǒng)給定設(shè)定值。因此，本節(jié)的實(shí)驗假定還原煤氣和搬出時間不變情況下，考察調(diào)節(jié)燃燒室溫度設(shè)定值對關(guān)鍵過程變量燃燒室溫度、還原帶溫度以及質(zhì)量指標(biāo)磁選管回收率的動態(tài)影響。這里參考文獻(xiàn)[15]，設(shè)定還原煤氣流量為2950 m3·h?1，設(shè)定初始燃燒室溫度為 1000℃，冷卻搬出時間為6.5 min。

首先，驗證磁選管回收率隨溫度變化趨勢。將燃燒室的溫度設(shè)定值從的1000℃，每次提高10℃，緩慢調(diào)節(jié)到1130℃。記錄磁選管回收率隨溫度變化的趨勢。實(shí)驗結(jié)果如圖5所示，可以看出：磁選管回收率和還原帶溫度隨燃燒室溫度的變化基本呈線性關(guān)系，這與之前設(shè)定的焙燒還原線性反應(yīng)率關(guān)系一致。

其次，驗證還原帶溫度和還原帶磁選管回收率相對于燃燒室溫度的階躍響應(yīng)。赤鐵礦豎爐焙燒工藝要求燃燒室溫度范圍在1000～1100℃之間。為此，將燃燒室溫度從1000℃直接調(diào)節(jié)到1100℃，記錄并觀察還原帶溫度、磁選管回收率對于溫度設(shè)定值階躍的動態(tài)響應(yīng)。等待磁選管回收率穩(wěn)定以后，進(jìn)一步將燃燒室溫度設(shè)定值由 1100℃調(diào)節(jié)到 1150℃，記錄并觀察還原帶溫度、磁選管回收率對于溫度設(shè)定值階躍的動態(tài)響應(yīng)。實(shí)驗結(jié)果如圖6所示。從圖中可以看出，在燃燒室初始溫度1000℃時，磁選管回收率穩(wěn)定在 81.0%左右，還原帶溫度穩(wěn)定在554℃。當(dāng)將燃燒室溫度設(shè)定值調(diào)整為 1100℃后，曲線經(jīng)過小范圍的波動很快收斂在1100℃，同時磁選管回收率穩(wěn)定在 82.3%、還原帶溫度穩(wěn)定到 571℃。這與文獻(xiàn)[9]中的實(shí)驗數(shù)據(jù)吻合。

圖5 燃燒室設(shè)定為1000～1130℃范圍內(nèi)磁選管回收率Fig. 5 Simulated MTRR under temperature

圖6 磁選管回收率階躍響應(yīng)曲線Fig. 6 Simulated MTRR response curve to combustion chamber temperature

5 結(jié) 論

（1）本文介紹的基于 METSIM 的豎爐還原焙燒仿真軟件，能準(zhǔn)確快速地模擬豎爐焙燒動態(tài)過程。在豎爐優(yōu)化設(shè)定方法設(shè)計中，以作為虛擬實(shí)驗對象，獲得在實(shí)際中難以在線連續(xù)檢測的質(zhì)量指標(biāo)，對于關(guān)鍵過程變量的動態(tài)響應(yīng)關(guān)系。這不僅可以用于檢驗過程控制策略，對于以節(jié)能降耗和提高生產(chǎn)質(zhì)量的優(yōu)化控制方法的設(shè)計，也提供了具有實(shí)用價值實(shí)驗手段。

（2）由于豎爐過程是一個復(fù)雜的工業(yè)過程，邊界條件變化頻繁，多變量之間的強(qiáng)非線性，在建模的過程中一些復(fù)雜的機(jī)理模型都進(jìn)行了簡化。所開發(fā)的METSIM模型中，還原反應(yīng)速率與燃燒室溫度與還原煤氣流量間的參數(shù)需要進(jìn)一步進(jìn)行現(xiàn)場實(shí)驗進(jìn)行調(diào)試得到準(zhǔn)確的參數(shù)值。目前開發(fā)的模型中，涉及豎爐過程中搬出時間的部分是根據(jù)豎爐的臺時產(chǎn)量進(jìn)行了理論推導(dǎo)，進(jìn)而在METSIM模型中實(shí)時地顯示了搬出時間，但還不能進(jìn)行實(shí)時的控制。

[1]Chai Tianyou(柴天佑), Ding Jinliang(丁進(jìn)良), Su Chunyi(蘇春翌).Hybrid intelligent optimal control method for operation of complex industrial processes [J].ACTA Automatica Sinica(自動化學(xué)報), 2008,34(5): 505-515

[2]Chai Tianyou, Ding Jinliang, Wu Fenghua. Hybrid intelligent control for optimal operation of shaft furnace roasting process [J].Control Engineering Practice, 2011, 19(3): 264-275

[3]Wu Fenghua(吳峰華), Yue Heng(岳恒), Chai Tianyou(柴天佑).Product Quality Monitoring System for Roasting Process of Shaft Furnace [J].Journal of Northeastern University:Natural Science(東北大學(xué)學(xué)報:自然科學(xué)版), 2007, 28(7): 913-916

[4]Zhou Ping(周平), Dai Wei(代偉), Chai Tianyou(柴天佑).Exploration of operational optimization control system for shaft furnace roasting process and its experiment study [J].Control Theory& Applications(控制理論與應(yīng)用), 2012, 29(12): 1565-1572

[5]Yan Aijun(嚴(yán)愛軍), Chai Tianyou(柴天佑), Yue Heng(岳恒).Multivariable intelligent optimizing control approach for shaft furnance roasting process [J].Acta Automatica Sinica(自動化學(xué)報),2006, 32(4): 636-640

[6]Hernández P C, Taboada M E, Graber T A, Galleguillos H R.Crystallization of hydrated ferric arsenate. process design using metsim [J].Industrial & Engineering Chemistry Research, 2009,48(23): 10522-10531

[7]Lü Yeqing(呂頁清), Zheng Shili(鄭詩禮), Wang Shaona(王少娜), Du Hao(杜浩), Zhang Yi(張懿). Application of METSIM in process design of the cleaner chromate production using sub-molten salt [J].The Chinese Journal of Process Engineering(過程工程學(xué)報), 2013,13(1):99-104

[8]Wang Ming(王明), Zan Qiming(諶啟明). Application of METSIM software in metallurgical engineering design [J].Rare Metals and Cemented Carbides(稀有金屬與硬質(zhì)合金), 2009, 37(4): 64-67

[9]Xu Wei(徐偉). Design and development of the software system for optimization control of shaft furnace roasting process[D]. Shenyang:Northeastern University, 2008

[10]Cen Kefa(岑可法), Yao Qiang(姚強(qiáng)), Luo Zhongyang(駱仲泱).Advanced Combustion(高等燃燒學(xué))[M]. Hangzhou: Zhejiang University Press, 2002

[11]Li Jianxiong(李建雄), Wang Zhenyan(王振艷). Study on hydrogen production by catalytic decomposition of methane [J].Journal of Henan Mechanical and Electrical Engineering College(河南機(jī)電高等專科學(xué)校學(xué)報), 2008, 16(1): 18-19

[12]Zhang Songshou(張松壽), Tong Zhengming(童正明), Zhou Wenzhu(周文鑄). Engineering Combustion(工程燃燒學(xué))[M].Shanghai: Shanghai Jiao Tong University Press, 1987

[13]Kitaev B I, Yaroshenko Y G, Suchkov V D. Heat Exchange in Shaft Furnaces[M]. London: Pergamon Press,1967

[14]Pang Jianming(龐建明), Guo Peimin(郭培民), Zhao Pei(趙沛), Cao Chaozhen(曹朝真), Zhang Dianwei(張殿偉). Study on kinetics of hematite reduction by hydrogen at low temperature [J].Iron and Steel(鋼鐵), 2008, 43(7):1-6

[15]Yan Aijun(嚴(yán)愛軍), Chai Tianyou(柴天佑), Wang Pu(王普).Multivariable intelligent optimizing setting method for shaft furnace roasting process [J].Journal of System Simulation(系統(tǒng)仿真學(xué)報),2008, 20(8): 2044-2047