氬氧精煉低碳鉻鐵終點(diǎn)碳含量預(yù)測(cè)模型

江 虹， 王宏志

（1.長(zhǎng)春工業(yè)大學(xué) 電氣與電子工程學(xué)院，吉林 長(zhǎng)春 130012；2.長(zhǎng)春工業(yè)大學(xué) 計(jì)算機(jī)科學(xué)與工程學(xué)院，吉林 長(zhǎng)春 130012）

0 引 言

冶金行業(yè)是我國(guó)實(shí)現(xiàn)節(jié)能減排目標(biāo)的重點(diǎn)，也是一個(gè)高能耗行業(yè)。氬氧煉鋼（Argon Oxygen Decarburization，AOD）是通過(guò)高壓的混合氧氣和氬氣吹入AOD爐內(nèi)，從而在較低溫度下也能實(shí)現(xiàn)脫碳保鉻的目的［1］。終點(diǎn)檢測(cè)的主要對(duì)象是終點(diǎn)溫度和碳含量，并決定是否能出鋼。如果這兩個(gè)指標(biāo)不滿足出鋼要求，控制模型將全啟動(dòng)進(jìn)行調(diào)整［2］，如果和控制終點(diǎn)要求相差較大，會(huì)對(duì)鋼包的壽命產(chǎn)生影響，而且各種調(diào)整材料的單耗增加［3］。由于冶煉過(guò)程復(fù)雜，無(wú)法確定各個(gè)計(jì)算參數(shù)的具體數(shù)值，對(duì)傳統(tǒng)的建模方法有局限。而神經(jīng)網(wǎng)絡(luò)有著明顯的優(yōu)越性，可以不受非線性模型的限制，廣泛應(yīng)用于非線性系統(tǒng)的建模和辨識(shí)［4］。

1 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的建立

1.1 神經(jīng)網(wǎng)絡(luò)技術(shù)

神經(jīng)網(wǎng)絡(luò)為解決那些復(fù)雜的、不確定的非線性控制問(wèn)題提供了有效的解決途徑［5］。由于BP（Back Propagation）網(wǎng)絡(luò)用的是負(fù)梯度下降法調(diào)節(jié)，會(huì)出現(xiàn)局部極小及收斂速度慢等問(wèn)題，因此用于函數(shù)逼近時(shí)，權(quán)值的調(diào)整方法存在局限性［6］。和BP網(wǎng)絡(luò)在分析能力、逼近能力、學(xué)習(xí)速度等方面進(jìn)行對(duì)比，明顯RBF（Radial Basis Function）網(wǎng)絡(luò)都要優(yōu)于BP網(wǎng)絡(luò)。這是因?yàn)镽BF網(wǎng)絡(luò)的輸出為隱單元輸出的線性加權(quán)和，學(xué)習(xí)速度加快［7］，而且RBF網(wǎng)絡(luò)使用徑向基函數(shù)，并且激活函數(shù)用的是高斯函數(shù)，神經(jīng)元里的輸入空間區(qū)域非常小，隱含層空間到輸出空間之間為線性映射，從而在實(shí)時(shí)控制要求較高的控制場(chǎng)合比較實(shí)用［8］。

RBF網(wǎng)絡(luò)輸入層可以將實(shí)驗(yàn)數(shù)據(jù)導(dǎo)入到隱含層，而輸出層通常由簡(jiǎn)單的線性函數(shù)構(gòu)成［9］，其結(jié)構(gòu)如圖1所示。

圖1 RBF網(wǎng)絡(luò)結(jié)構(gòu)

1.2 預(yù)測(cè)模型的建立

RBF網(wǎng)絡(luò)的學(xué)習(xí)算法需要確定徑向基函數(shù)中心、方差、隱含層到輸出層的權(quán)值。

1.2.1 RBF網(wǎng)絡(luò)中心的選取

文中采用的是近鄰法，以這種方法來(lái)選取中心，具體步驟如下：

1）選擇高斯函數(shù)寬度r；

2）從第一對(duì)實(shí)驗(yàn)數(shù)據(jù)對(duì)（X1，Y1）開(kāi)始，將它作為學(xué)習(xí)樣本，在X1上建立一個(gè)聚類中心，使得C1＝X1；

3）從第二對(duì)實(shí)驗(yàn)數(shù)據(jù)對(duì)（X2，Y2），求出X2到C1這個(gè)聚類中心的距離｜X2－C1｜。若｜X2－C1｜＜r，則C1為X2的鄰近聚類；若｜X2－C1｜＞r，則將X2作為一個(gè)新聚類中心，并C2＝X2，這就需要再添加一個(gè)隱單元C2；

4）進(jìn)行到第k個(gè)實(shí)驗(yàn)數(shù)據(jù)對(duì)（Xk，Yk）時(shí)，已經(jīng)找到M個(gè)聚類中心，中心點(diǎn)C1，C2，…，CM，分別求出到這M個(gè)聚類中心的距離，｜Xk－Ci｜，i＝1，2，…，M，設(shè)｜Xk－Cj｜為最小距離，即Cj為Xk最近鄰聚類。

1.2.2 RBF網(wǎng)絡(luò)權(quán)值調(diào)整

網(wǎng)絡(luò)權(quán)值調(diào)整是神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的一部分。為了能夠進(jìn)行有效的調(diào)整，在研究分析調(diào)整方法和大量閱讀文獻(xiàn)之后，尋找選擇了一種最小二乘法，該方法具有加權(quán)遺忘因子，可以得到比較符合期望的權(quán)值，具體步驟如下：

1）賦予權(quán)ωi（i＝1，2，…，M）初值；

2）令循環(huán)變量k＝1；

3）計(jì)算隱層節(jié)點(diǎn)輸出

式中：Q（0）——m×m的正定對(duì)角陣，m為隱含層節(jié)點(diǎn)數(shù)；

ρ——遺忘因子，0＜ρ＜1。

4）計(jì)算目標(biāo)累積誤差。鐵合金模型預(yù)測(cè)值與冶煉實(shí)際終點(diǎn)碳含量之間的累積誤差E（k），它對(duì)應(yīng)第k個(gè)樣本。

φk——隱節(jié)點(diǎn)輸出向量；

yd（k）——期望輸出。

5）判斷是否E（k）＜E，如果E（k）＜E，則結(jié)束訓(xùn)練；如果E（k）＞E，判斷是否k＜N。如果k＜N，則k自加并計(jì)算隱層節(jié)點(diǎn)輸出；否則重新賦初值。其中N為輸入實(shí)驗(yàn)數(shù)據(jù)數(shù)量，E為建模前預(yù)先設(shè)定的誤差的目標(biāo)。

1.2.3 數(shù)據(jù)的標(biāo)準(zhǔn)化處理

根據(jù)AOD爐冶煉鐵合金生產(chǎn)低碳鉻鐵的實(shí)際入爐情況，冶煉過(guò)程中加入的入爐鐵水量、氧氣量、氬氣量，以及補(bǔ)吹階段加入的鉻礦、硅鐵等物質(zhì)的數(shù)量級(jí)和單位相差較大，要想將這些作為模型的輸入量，就必須將這些變量的單位等作出標(biāo)準(zhǔn)化處理。

其中

式中，i＝1，2，…，7，對(duì)應(yīng)于輸入變量；k＝1，2，…，p，對(duì)應(yīng)于樣本。

對(duì)輸入量做了標(biāo)準(zhǔn)化處理，輸出量yc也要進(jìn)行標(biāo)準(zhǔn)化處理。這樣才能保證模型的兩端與實(shí)際情況是對(duì)等的。下面開(kāi)始處理輸出量，神經(jīng)元訓(xùn)練學(xué)習(xí)后，就會(huì)得到模型預(yù)測(cè)的碳含量y，然而它并不能代表實(shí)際的訓(xùn)練結(jié)果，對(duì)此可以按下式進(jìn)行處理：

y——訓(xùn)練后網(wǎng)絡(luò)的實(shí)際輸出值；

y∧——通過(guò)標(biāo)準(zhǔn)化數(shù)據(jù)訓(xùn)練網(wǎng)絡(luò)后的輸出值。

2 實(shí)驗(yàn)與分析

根據(jù)前面已介紹過(guò)的預(yù)測(cè)模型，可以確定輸入層節(jié)點(diǎn)個(gè)數(shù)為7，輸出層節(jié)點(diǎn)為1。對(duì)隱含層的節(jié)點(diǎn)數(shù)選取10，再對(duì)網(wǎng)絡(luò)性能的狀況進(jìn)行分析，然后從兩邊各選一個(gè)數(shù)，比如4和14，對(duì)上面這幾種情況得到的預(yù)測(cè)性能進(jìn)行對(duì)比分析，就可以選擇出一個(gè)比較合理的隱含層節(jié)點(diǎn)個(gè)數(shù)。具體預(yù)測(cè)結(jié)果見(jiàn)表1。

表1 不同隱含層節(jié)點(diǎn)數(shù)對(duì)應(yīng)的預(yù)測(cè)結(jié)果對(duì)比

分析表1可知，在誤差準(zhǔn)則和精度都相同的情況下，隱含層節(jié)點(diǎn)數(shù)為14時(shí)的預(yù)測(cè)效果較好，所以合理的隱含層節(jié)點(diǎn)數(shù)確定為14。

在誤差準(zhǔn)則E＝0.002條件下對(duì)表1中的數(shù)據(jù)做標(biāo)準(zhǔn)化處理，當(dāng)終點(diǎn)碳含量模型中的精度rc＝0.1，隱含層節(jié)點(diǎn)數(shù)為14時(shí)，建立碳含量預(yù)測(cè)模型，并對(duì)碳含量作出5次預(yù)報(bào)，預(yù)報(bào)結(jié)果如圖2所示。

圖2 RFB預(yù)報(bào)模型結(jié)果

RBF神經(jīng)網(wǎng)絡(luò)終點(diǎn)預(yù)測(cè)與實(shí)際碳含量對(duì)比見(jiàn)表2。

表2 RBF神經(jīng)網(wǎng)絡(luò)終點(diǎn)預(yù)測(cè)與實(shí)際碳含量對(duì)比

觀察表2可知偏差均在0.02以內(nèi)，該模型的準(zhǔn)確性比較高。

3 結(jié) 語(yǔ)

對(duì)冶金終點(diǎn)碳含量的預(yù)測(cè)模型使用RBF神經(jīng)網(wǎng)絡(luò)。首先用實(shí)際得到的冶煉數(shù)據(jù)作為對(duì)象，并且對(duì)這些數(shù)據(jù)作標(biāo)準(zhǔn)化處理，減少預(yù)測(cè)誤差和對(duì)訓(xùn)練樣本數(shù)據(jù)造成的影響。然后確定RBF網(wǎng)絡(luò)中心，調(diào)整網(wǎng)絡(luò)權(quán)值，進(jìn)而得出終點(diǎn)預(yù)測(cè)模型。通過(guò)預(yù)測(cè)模型調(diào)整控制量，當(dāng)預(yù)測(cè)偏差較大時(shí)，通過(guò)控制模型實(shí)時(shí)計(jì)算出所需加入的冷卻劑或補(bǔ)吹的氧氣量，并加以調(diào)整，從而達(dá)到出鋼要求。最后通過(guò)實(shí)驗(yàn)證明，用RBF神經(jīng)網(wǎng)絡(luò)構(gòu)建的氬氧精煉低碳鉻鐵終點(diǎn)碳含量的預(yù)測(cè)模型，預(yù)測(cè)精度比較準(zhǔn)確，可以為控制模型提供有價(jià)值的參考信息。

［1］ 李彪，曹吉花.AOD精煉鐵合金終點(diǎn)預(yù)測(cè)模型的建立［J］.宿州學(xué)院學(xué)報(bào)，2012，8（27）：49－52.

［2］ 陳國(guó)翠.我國(guó)鉻鐵合金生產(chǎn)與發(fā)展［J］.鐵合金，2000，31（5）：40－42.

［3］ 于歡歡，許紅巖，張洋，等.模糊控制在AOD爐冶煉鉻鐵合金噴濺系統(tǒng)中的應(yīng)用［J］.長(zhǎng)春工業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，32（1）：1－6.

［4］ 劉文林，錢曉龍，張雪松.ControlLogix系統(tǒng)在AOD爐控制系統(tǒng)改造中的應(yīng)用［J］.沈陽(yáng)工程學(xué)院學(xué)報(bào)，2006，2（1）：54－57.

［5］ 王海，陳宏偉，張麗榮.氧氣頂吹轉(zhuǎn)爐自動(dòng)控制系統(tǒng)的開(kāi)發(fā)與應(yīng)用［J］.電工技術(shù)雜志，2004（6）：57－60.

［6］ 舒杰輝，魏季和，池和冰，等.不銹鋼AOD精煉過(guò)程數(shù)學(xué)模擬進(jìn)展［J］.煉鋼，2004，20（3）：53－56.

［7］ 池和冰，魏季和，舒杰輝，等.不銹鋼AOD轉(zhuǎn)爐精煉過(guò)程數(shù)學(xué)模擬初探［J］.上海金屬，2007，29（4）49－58.

［8］ 許紅巖，于歡歡.AOD爐鐵合金冶煉終點(diǎn)預(yù)報(bào)模型［J］.長(zhǎng)春工業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，32（4）：370－374.

［9］ 王登峰，倪紅衛(wèi).人工神經(jīng)網(wǎng)絡(luò)在轉(zhuǎn)爐煉鋼終點(diǎn)預(yù)報(bào)中的應(yīng)用研究［J］.鋼鐵研究，2005，4（2）：27－31.