王心哲， 韓 敏

（大連理工大學(xué)電子與信息工程學(xué)院，遼寧大連 116024）

0 引 言

轉(zhuǎn)爐煉鋼是一個(gè)降碳升溫的過程，主要目的是在吹煉終點(diǎn)時(shí)得到目標(biāo)碳含量和溫度滿足工藝要求的合格鋼水.通常，對(duì)轉(zhuǎn)爐煉鋼生產(chǎn)的控制包括靜態(tài)控制和動(dòng)態(tài)控制.靜態(tài)控制模型是動(dòng)態(tài)控制模型的基礎(chǔ)，其精度直接影響到動(dòng)態(tài)控制的效果，如果靜態(tài)模型不精確，僅依靠吹煉末期的動(dòng)態(tài)控制來(lái)調(diào)整碳溫，不利于終點(diǎn)碳含量和溫度的命中［1～3］.因此，靜態(tài)控制模型的研究至關(guān)重要.

目前已有的靜態(tài)控制模型主要有4種：理論模型、統(tǒng)計(jì)模型、增量模型和神經(jīng)網(wǎng)絡(luò)模型，但由于煉鋼過程反應(yīng)的復(fù)雜性和隨機(jī)性，許多影響因素很難用準(zhǔn)確的數(shù)學(xué)方程或統(tǒng)計(jì)方法進(jìn)行描述，且神經(jīng)網(wǎng)絡(luò)模型受限于網(wǎng)絡(luò)結(jié)構(gòu)和參數(shù)的調(diào)整，已有模型的控制精度受到很大影響［4～6］.

案例推理（case－based reasoning，CBR）方法是人工智能發(fā)展較為成熟的一個(gè)分支，它基于歷史經(jīng)驗(yàn)和知識(shí)進(jìn)行推理，已在故障診斷、自動(dòng)糾錯(cuò)系統(tǒng)和決策支持等領(lǐng)域取得成功應(yīng)用［7～11］，但在復(fù)雜工業(yè)過程的建模與控制方面仍處于探索階段.其主要思想是根據(jù)當(dāng)前問題的描述，在案例庫(kù)中搜索與當(dāng)前問題相似的案例，再利用案例庫(kù)中搜索到的相似案例的解，構(gòu)造當(dāng)前問題的解.

本文使用案例推理方法建立轉(zhuǎn)爐煉鋼靜態(tài)控制模型，并將因果關(guān)系分析引入案例推理系統(tǒng)，分析各變量對(duì)解空間的影響程度，選擇對(duì)解空間影響較大的變量作為案例的條件屬性.

1 模型的建立

案例推理是一種類比推理方法，主要有以下步驟：案例描述、案例檢索、案例重用、案例修正和案例保存.首先建立對(duì)問題（案例）的描述，根據(jù)案例描述，在案例庫(kù)中進(jìn)行檢索，找到最相似的若干案例；再利用搜索到的相似問題的解構(gòu)造當(dāng)前問題的解；最后，將得到的解輸出給用戶，如果成功解決問題，則將此案例存入案例庫(kù).基于案例推理的轉(zhuǎn)爐煉鋼靜態(tài)控制模型結(jié)構(gòu)如圖1所示.

1.1 案例描述

案例描述是案例推理過程的基礎(chǔ)，是案例檢索和案例重用有效性的保證.案例描述就是確定用于描述案例的各屬性和結(jié)構(gòu)的過程，通常根據(jù)背景知識(shí)或經(jīng)驗(yàn)確定.本文將因果關(guān)系分析引入案例推理系統(tǒng)，用于對(duì)案例屬性的選擇.

圖1 轉(zhuǎn)爐煉鋼案例推理靜態(tài)控制模型結(jié)構(gòu)Fig.1 Structure of CBR model for static control of converter steelmaking

1.1.1 因果關(guān)系 因果關(guān)系的思想最早由Wiener提出，并由Granger將其公式化［12］，Granger因果關(guān)系被廣泛應(yīng)用于各領(lǐng)域.因果關(guān)系分析的主要思想是，對(duì)于一個(gè)變量的當(dāng)前值，如果用該變量的過去值和第2個(gè)變量的過去值估計(jì)與僅由該變量自身的過去值估計(jì)相比，得到的誤差項(xiàng)的方差減小，則認(rèn)為第2個(gè)變量與該變量存在因果關(guān)系，且第2個(gè)變量為該變量的因.

對(duì)于一個(gè)固定的變量x（t），首先考慮基于m個(gè)歷史測(cè)量值的自回歸預(yù)測(cè)：

其中εx（t）為預(yù)測(cè)誤差，其大小可表示為方差形式var（εx（t））.在自回歸的基礎(chǔ)上，引入另一組固定變量y（t）的歷史測(cè)量值，則可得到下面的預(yù)測(cè)：

如果得到的預(yù)測(cè)誤差的方差小于自回歸預(yù)測(cè)誤差的方差，即var（εx｜y（t））＜var（εx（t）），則變量y（t）對(duì)變量x（t）具有因果關(guān)系.

二變量間的因果關(guān)系分析雖然可以得到是否存在因果關(guān)系，但無(wú)法確定是直接關(guān)系還是間接關(guān)系，即存在如圖2所示的兩種情況：

其中，圖2（a）表示變量z對(duì)x的因果關(guān)系為間接關(guān)系，是通過變量y實(shí)現(xiàn)的.圖2（b）則表示變量z與x間存在直接的因果關(guān)系.為此，使用條件因果關(guān)系分析來(lái)確定變量間的因果關(guān)系屬于哪種模式.

圖2 兩種因果關(guān)系模式Fig.2 Two patterns of causal interactions

與二變量間因果關(guān)系的分析相似，已知變量z與x之間存在因果關(guān)系，在式（2）的基礎(chǔ)上引入變量z得到的預(yù)測(cè)表達(dá)式為

如果得到的預(yù)測(cè)誤差的方差小于式（2）的回歸預(yù)測(cè)誤差的方差， 即 var（εx｜yz（t））＜var（εx｜y（t）），則變量z（t）對(duì)變量x（t）具有直接的因果關(guān)系，否則說明具有間接的因果關(guān)系.

1.1.2 案例屬性的確定 根據(jù)轉(zhuǎn)爐煉鋼生產(chǎn)的工況信息、工藝目標(biāo)和過程數(shù)據(jù)構(gòu)建案例的問題屬性和解屬性.影響靜態(tài)控制的因素很多，如工藝目標(biāo)要求、吹煉過程操作以及原料成分等.轉(zhuǎn)爐煉鋼靜態(tài)控制模型的案例描述包括兩種屬性：

（1）離散屬性

離散屬性主要考慮吹煉終點(diǎn)鋼水的目標(biāo)溫度和碳含量.根據(jù)操作規(guī)程，冶煉的鋼種不同時(shí)，鋼水終點(diǎn)目標(biāo)溫度和碳含量的高低有所不同.因此，案例的離散屬性包括鋼水的目標(biāo)溫度和目標(biāo)碳含量.由此，離散屬性集合可表示為｛目標(biāo)溫度，目標(biāo)碳含量｝.

（2）連續(xù)屬性

連續(xù)屬性的構(gòu)建主要考慮設(shè)備信息和原料信息，包括氧槍的槍齡、轉(zhuǎn)爐的爐齡、鐵水成分和溫度、廢鋼裝入量和輔原料加入量.由此，連續(xù)屬性集合可表示為｛設(shè)備信息，鐵水信息，原料加入量信息｝.可用于描述連續(xù)屬性的變量包括表示設(shè)備信息的氧槍的槍齡和轉(zhuǎn)爐的爐齡；表示鐵水信息的鐵水的碳含量、硅含量、錳含量、磷含量、硫含量和鐵水溫度；表示原料加入量信息的廢鋼裝入量、石灰加入量和白云石加入量.

如果選擇所有的變量作為描述案例的屬性，不僅會(huì)增加計(jì)算案例間相似度時(shí)的計(jì)算量，而且那些不重要的屬性作為噪聲會(huì)影響相似案例的選擇，從而降低模型計(jì)算的精度.同時(shí)，案例屬性過多還會(huì)使案例庫(kù)的規(guī)模增長(zhǎng)過快，不利于案例庫(kù)的維護(hù).由此，本文利用Granger因果關(guān)系和條件因果關(guān)系對(duì)案例的連續(xù)屬性進(jìn)行確定，選擇重要的變量并剔除冗余變量，以提高案例推理系統(tǒng)的運(yùn)行效率和有效性.

使用因果關(guān)系分析確定案例屬性后，案例的構(gòu)建描述為如圖3所示的4層結(jié)構(gòu).

圖3 案例屬性結(jié)構(gòu)Fig.3 Structure of case attributes

1.2 案例檢索

案例檢索是案例推理過程的關(guān)鍵，是構(gòu)造當(dāng)前問題精確解的基礎(chǔ)，其核心步驟是歷史案例與當(dāng)前案例間相似度的計(jì)算.目前較常用的案例檢索方法有最鄰近法、歸納索引法和知識(shí)引導(dǎo)法等［13］.通過案例檢索，獲得與當(dāng)前問題最相似的若干案例，檢索得到的相似案例的解用于案例重用.

根據(jù)案例描述的結(jié)構(gòu)，分別計(jì)算離散屬性和連續(xù)屬性相似度，整體相似度按照案例描述的4層結(jié)構(gòu)計(jì)算.

1.2.1 離散屬性相似度的計(jì)算 案例的離散屬性與生產(chǎn)的工藝目標(biāo)相關(guān)，主要包括鋼水目標(biāo)溫度和目標(biāo)碳含量，這也是衡量鋼水是否合格的兩個(gè)指標(biāo).離散屬性相似度的計(jì)算式如下：

式中：simd為離散屬性的相似度；simaim＿C為目標(biāo)碳含量的相似度，simaim＿t為目標(biāo)溫度的相似度.目標(biāo)溫度和目標(biāo)碳含量相似度的計(jì)算式為

其中fc為案例庫(kù)中某個(gè)案例的離散屬性值，f為當(dāng)前案例的離散屬性值.

1.2.2 連續(xù)屬性相似度的計(jì)算 與離散屬性相似度的計(jì)算方法相同，連續(xù)屬性相似度采用最鄰近法計(jì)算，并且按照連續(xù)屬性本身的3層結(jié)構(gòu)進(jìn)行計(jì)算.

頂層為連續(xù)屬性相似度，相似度的計(jì)算表達(dá)式如下式所示：

式中：simc為案例連續(xù)屬性的整體相似度；sim i為一個(gè)屬性集合的相似度，分別表示設(shè)備信息、鐵水信息、原料加入量信息這3個(gè)屬性集合的相似度；w i為對(duì)應(yīng)每個(gè)屬性集合的相似度權(quán)值，由經(jīng)驗(yàn)確定.

各屬性集合中的屬性元素均為連續(xù)型變量，采用如下的公式計(jì)算各屬性集合的相似度：

式中：m為第i個(gè)屬性集合中屬性的個(gè)數(shù)；f ij為當(dāng)前爐次第i個(gè)屬性集合中第j個(gè)屬性值；fc，ij為案例庫(kù)中第i個(gè)屬性集合中第j個(gè)屬性值；sim ij（f ij，fc，ij）為第i個(gè)屬性集合中第j個(gè)屬性的相似度，即底層相似度，計(jì)算公式如式（5）所示；w ij為第i個(gè)屬性集合中第j個(gè)屬性相似度的權(quán)值，其值的確定過程如下.

首先，計(jì)算每個(gè)屬性相似度的標(biāo)準(zhǔn)差：

其中n為案例庫(kù)中案例的個(gè)數(shù)；sim k為當(dāng)前案例與案例庫(kù)中第k個(gè)案例間第j個(gè)屬性的相似度；為當(dāng)前案例與案例庫(kù)中案例間第j個(gè)屬性相似度的均值.

再利用式（8）的計(jì)算結(jié)果確定各屬性權(quán)值的大小，如下式所示：

對(duì)于相同的屬性，各個(gè)sim kj的值越接近，則σj越小，即不同的案例在該屬性上的差異越小，也表示該屬性在區(qū)分各個(gè)案例的相似程度上作用不大，則被賦予相對(duì)較小的權(quán)值.相反的情況被賦予較大的權(quán)值.

1.2.3 案例整體相似度的表達(dá) 案例整體相似度simt可由離散屬性相似度和連續(xù)屬性相似度的集合形式表達(dá)，即整體相似度包括離散屬性相似度和連續(xù)屬性相似度兩個(gè)元素，表示為simt＝｛simd，simc｝.

對(duì)于當(dāng)前案例，首先計(jì)算當(dāng)前案例與案例庫(kù)中案例間離散屬性的相似度，選擇滿足一定離散屬性相似度要求的案例；在此基礎(chǔ)上，計(jì)算當(dāng)前案例與選擇案例間的連續(xù)屬性相似度，并將計(jì)算結(jié)果按相似度的大小降序排序，選擇與當(dāng)前案例連續(xù)屬性相似度值最大的k個(gè)歷史案例作為最終的相似案例，相似案例的解在下一步的案例重用中使用.

根據(jù)案例描述的結(jié)構(gòu)，將案例檢索過程按照離散屬性和連續(xù)屬性分兩步進(jìn)行，基于離散屬性的案例檢索是連續(xù)屬性案例檢索的基礎(chǔ)，案例檢索的最終結(jié)果由連續(xù)屬性的案例檢索結(jié)果直接決定.

1.3 案例重用和修正

案例重用為當(dāng)前問題提出建議解或者最終解，案例重用的方式直接影響到對(duì)當(dāng)前問題解的精確性.

通過案例檢索得到與當(dāng)前問題最相似的k個(gè)案例的解，使用相似案例的解構(gòu)造當(dāng)前問題的解，計(jì)算表達(dá)式為

式中：S為構(gòu)造的當(dāng)前案例解；si為第i個(gè)相似案例的解；k為相似案例的個(gè)數(shù)；w i為權(quán)值，其大小為simc.

1.4 案例存儲(chǔ)和維護(hù)

案例的存儲(chǔ)和維護(hù)對(duì)案例推理系統(tǒng)的運(yùn)行十分必要.案例存儲(chǔ)，可以豐富案例庫(kù)中的案例.但隨著生產(chǎn)的進(jìn)行會(huì)造成案例庫(kù)的過度膨脹，降低系統(tǒng)的運(yùn)行效率，需要對(duì)案例庫(kù)進(jìn)行維護(hù)，以保證運(yùn)行效率.如果當(dāng)前爐次的吹煉結(jié)果滿足工藝要求，且與歷史案例的相似度小于一個(gè)閾值，則本爐次作為一個(gè)成功案例存入案例庫(kù).

2 仿真研究

使用某鋼廠一座150 t轉(zhuǎn)爐現(xiàn)場(chǎng)采集的實(shí)際生產(chǎn)數(shù)據(jù)進(jìn)行仿真研究，選擇630爐不使用復(fù)吹且終點(diǎn)碳含量和溫度均命中的爐次.選擇前500爐數(shù)據(jù)作為案例庫(kù)中的案例，對(duì)其余130組進(jìn)行靜態(tài)控制計(jì)算.

使用1.2的案例檢索方法在案例庫(kù)中搜索與第501爐相似的爐次案例，選擇與當(dāng)前爐次最相似的前13個(gè)案例的解計(jì)算當(dāng)前爐次的吹氧量和冷卻劑加入量.計(jì)算結(jié)束后，將第501爐數(shù)據(jù)存儲(chǔ)在案例庫(kù)中，作為成功案例用于對(duì)下一爐的計(jì)算，依此類推，完成所有130爐的靜態(tài)控制計(jì)算，并保證案例的更新.

2.1 模型計(jì)算結(jié)果

按照上述方法對(duì)實(shí)際數(shù)據(jù)進(jìn)行仿真，得到靜態(tài)吹氧量的模型計(jì)算值與實(shí)際值之間的均方根誤差為302.36 m3，誤差限度為500 m3時(shí)的準(zhǔn)確率為90.00%；冷卻劑加入量的模型計(jì)算值與實(shí)際值之間的均方根誤差為0.65 t，誤差限度為1 t時(shí)的準(zhǔn)確率為86.15%.使用案例推理計(jì)算的吹氧量Qo和冷卻劑加入量Qc與實(shí)際值間的曲線圖分別如圖4和5所示.

2.2 與選擇全部變量作為案例屬性模型的比較

將因果關(guān)系分析用于對(duì)案例屬性的確定，由此減小了案例描述的復(fù)雜度，且能控制案例庫(kù)的規(guī)模，有利于案例庫(kù)的維護(hù).為了進(jìn)一步驗(yàn)證該方法的有效性，將應(yīng)用因果關(guān)系分析確定案例屬性的模型（C－CBR）計(jì)算結(jié)果與選擇全部變量作為案例屬性的模型（CBR）計(jì)算結(jié)果進(jìn)行比較，結(jié)果如表1所示.

圖4 吹氧量計(jì)算值與實(shí)際值間曲線Fig.4 Curves of calculation values and real values of oxygen blown amount

圖5 冷卻劑加入量計(jì)算值與實(shí)際值間曲線Fig.5 Curves of calculation values and real values of coolant addition amount

表1 不同案例描述時(shí)的模型計(jì)算結(jié)果比較Tab.1 Comparison of the results with different case description

表1中，吹氧量的準(zhǔn)確率的誤差范圍為小于500 m3，冷卻劑量的準(zhǔn)確率的誤差范圍為小于1 t.

由表1可以看出，引入因果關(guān)系分析對(duì)案例屬性確定后，不會(huì)影響模型的計(jì)算精度.

2.3 與現(xiàn)有方法的比較

將使用案例推理方法建立的靜態(tài)控制模型的計(jì)算結(jié)果與文獻(xiàn)［4、5］中的靜態(tài)控制方法進(jìn)行比較.文獻(xiàn)［4］中對(duì)120 t轉(zhuǎn)爐使用增量回歸方法建模；文獻(xiàn)［5］中對(duì)250 t轉(zhuǎn)爐使用遺傳算法結(jié)合徑向基神經(jīng)網(wǎng)絡(luò)方法（GA－RBF）計(jì)算，比較結(jié)果分別見表2和3.

通過比較可以看出，使用案例推理方法建立的轉(zhuǎn)爐煉鋼靜態(tài)控制模型無(wú)論在精度還是在準(zhǔn)確率上均優(yōu)于已有方法.

表2 吹氧量的比較Tab.2 Comparison of oxygen blown amount

表3 冷卻劑加入量的比較Tab.3 Comparison of coolant addition amount

3 結(jié) 論

使用因果分析方法對(duì)案例屬性進(jìn)行確定，能夠在保證模型精度的前提下減小案例描述的復(fù)雜度，降低案例庫(kù)的存儲(chǔ)空間，提高案例庫(kù)的維護(hù)效率.仿真結(jié)果顯示，使用該方法建立的靜態(tài)控制模型的效果優(yōu)于已有方法，對(duì)實(shí)際生產(chǎn)具有一定的指導(dǎo)意義.

［1］黃金俠，金寧德.轉(zhuǎn)爐冶煉終點(diǎn)靜態(tài)控制預(yù)測(cè)模型［J］.煉鋼，2006，22（1）：45－48

［2］陳忠偉，袁守謙.LD轉(zhuǎn)爐冶煉的靜態(tài)數(shù)學(xué)模型及實(shí)現(xiàn)［J］.煉鋼，2000，16（5）：31－34

［3］朱光俊，梁本川.轉(zhuǎn)爐煉鋼靜態(tài)控制優(yōu)化模型［J］.煉鋼，1999，15（4）：25－28

［4］孟祥寧，張海鷹，朱苗勇.轉(zhuǎn)爐煉鋼過程靜態(tài)控制模型的改進(jìn)［J］.材料與冶金學(xué)報(bào)，2004，3（4）：246－249

［5］王建輝，徐 林，方曉柯，等.基于一類GA－RBF神經(jīng)網(wǎng)絡(luò)的轉(zhuǎn)爐煉鋼靜態(tài)模型控制［J］.東南大學(xué)學(xué)報(bào)，2005，35（2）：90－94

［6］丁 容，劉 瀏 轉(zhuǎn)爐煉鋼過程人工智能靜態(tài)控制模型［J］.鋼鐵，1997，32（1）：22－26

［7］CHANG Pei－chann，LIU Chen－h(huán)ao，LAI R K.A fuzzy case－based reasoning model for sales forecasting in print circuit board industries［J］.Expert Systems with Applications，2008，34（3）：2049－2058

［8］HUANG Mu－jung，CHEN Mu－yen，LEE Showchin.Integrating data mining with case－based reasoning for chronic diseases prognosis and diagnosis［J］.Expert Systems with Applications，2007，32（3）：856－867

［9］NEGNY S，LE LANN J－M.Case－based reasoning for chemical engineering design［J］.Chemical Engineering Research and Design，2008，8（6）：648－658

［10］GUILLERMO C R，NEGNY S，LE LANN J－M. Case－based reasoning and TRIZ：A coupling for innovative conception in chemical engineering［J］.Chemical Engineering and Processing，2009，4（1）：239－249

［11］TUNG Yuan－h(huán)sin，TSENG Shian－shyong，WENG Jui－feng，etal.A rule－based CBR approach for expert finding and problem diagnosis［J］.Expert Systems with Applications，2010，37（3）：2427－2438

［12］GRANGER C W J.Investigation causal relations by econometric models and cross－spectral methods［J］.Econometrica，1969，37（3）：424－438

［13］李偉明，穆志純.基于案例推理的冷軋負(fù)荷分配模型研究［J］.系統(tǒng)仿真學(xué)報(bào)，2006，18（8）：2162－2165