沈 坤,劉有勢(shì),劉錄光,南晨晨

(湖南師范大學(xué)信息科學(xué)與工程學(xué)院，中國 長(zhǎng)沙 410081)

在實(shí)際煉焦生產(chǎn)過程中，為提高配煤的準(zhǔn)確性，配煤專家根據(jù)實(shí)際煉焦工況及配合煤質(zhì)量指標(biāo)預(yù)測(cè)焦炭質(zhì)量，并基于預(yù)測(cè)結(jié)果為配煤方案的優(yōu)化提供數(shù)據(jù)支撐。焦炭質(zhì)量預(yù)測(cè)的理論依據(jù)是配合煤的成焦理論，而在實(shí)際生產(chǎn)過程中，普遍基于煉焦生產(chǎn)數(shù)據(jù)，采用焦炭質(zhì)量預(yù)測(cè)模型進(jìn)行計(jì)算。

用于焦炭質(zhì)量預(yù)測(cè)的基礎(chǔ)數(shù)據(jù)主要包括小焦?fàn)t試驗(yàn)數(shù)據(jù)和實(shí)際煉焦生產(chǎn)數(shù)據(jù)[1-2]，其區(qū)別在于試驗(yàn)數(shù)據(jù)可以提供更精細(xì)的煤、焦質(zhì)量指標(biāo)信息，但成本高、耗時(shí)長(zhǎng)，難以在日常煉焦生產(chǎn)過程中使用，因此主要用于焦化及配煤理論研究。實(shí)際煉焦生產(chǎn)數(shù)據(jù)受焦化企業(yè)化驗(yàn)室技術(shù)水平、化驗(yàn)成本及質(zhì)量管理的影響僅包含有限的煤、焦質(zhì)量指標(biāo)信息，但這些指標(biāo)可以按班組檢測(cè)，形成指導(dǎo)生產(chǎn)的質(zhì)量報(bào)表，且滿足焦炭質(zhì)量國家標(biāo)準(zhǔn)的檢測(cè)要求，能真實(shí)地反映煉焦生產(chǎn)過程，可直接用于焦炭質(zhì)量預(yù)測(cè)。

用于焦炭質(zhì)量預(yù)測(cè)的方法主要包括以回歸分析為代表的統(tǒng)計(jì)學(xué)方法[2-3]和以神經(jīng)網(wǎng)絡(luò)為代表的智能方法[4-7]。當(dāng)前，具有逼近任意非線性函數(shù)的神經(jīng)網(wǎng)絡(luò)成為了焦炭質(zhì)量預(yù)測(cè)的主要技術(shù)手段，其中采用BP神經(jīng)網(wǎng)絡(luò)[4]、RBF神經(jīng)網(wǎng)絡(luò)[5]以及自適應(yīng)神經(jīng)模糊系統(tǒng)[6]構(gòu)建的焦炭質(zhì)量預(yù)測(cè)方法得到了廣泛的應(yīng)用。但上述文獻(xiàn)主要關(guān)注于神經(jīng)網(wǎng)絡(luò)的設(shè)計(jì)以及輸入變量的選取，并偏重于從成焦理論上分析配合煤質(zhì)量指標(biāo)與焦炭質(zhì)量指標(biāo)之間的關(guān)系，從煉焦工藝上分析煉焦工況對(duì)焦炭質(zhì)量的影響，而對(duì)用于訓(xùn)練神經(jīng)網(wǎng)絡(luò)的實(shí)際煉焦生產(chǎn)數(shù)據(jù)及其關(guān)系缺乏深入的分析。雖然文獻(xiàn)[7]提到了用于焦炭質(zhì)量預(yù)測(cè)的模型數(shù)據(jù)時(shí)序匹配困難，但并沒有給出具體可行的數(shù)據(jù)匹配方法。在實(shí)際煉焦生產(chǎn)過程中，煤、焦質(zhì)量及煉焦工況數(shù)據(jù)以質(zhì)量報(bào)表和生產(chǎn)工藝報(bào)表的形式呈現(xiàn)，這些報(bào)表數(shù)據(jù)是配煤專家進(jìn)行配煤的依據(jù)，因此從報(bào)表中分析和提取有效信息并完成時(shí)序匹配是進(jìn)行焦炭質(zhì)量預(yù)測(cè)的必要工作。

由以上分析可知，基礎(chǔ)數(shù)據(jù)和預(yù)測(cè)方法是實(shí)現(xiàn)焦炭質(zhì)量預(yù)測(cè)的關(guān)鍵要素，已有的研究主要關(guān)注于預(yù)測(cè)方法的設(shè)計(jì)，而很少涉及對(duì)實(shí)際煉焦生產(chǎn)數(shù)據(jù)的分析與匹配。本文針對(duì)實(shí)際焦化企業(yè)的報(bào)表數(shù)據(jù)，分析、提取煉焦生產(chǎn)過程中煤、焦、工況數(shù)據(jù)及其匹配關(guān)系，并以此為基礎(chǔ)采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建符合實(shí)際煉焦生產(chǎn)過程的焦炭質(zhì)量預(yù)測(cè)方法。

1 實(shí)際煉焦生產(chǎn)數(shù)據(jù)分析

焦化企業(yè)在煉焦生產(chǎn)過程中，由化驗(yàn)室按班組采集單煤、配合煤、焦炭樣本進(jìn)行化驗(yàn)，并以質(zhì)量報(bào)表的形式呈現(xiàn)各樣本的質(zhì)量指標(biāo)。本文對(duì)湖南某焦化企業(yè)5.5 m搗固焦質(zhì)量日?qǐng)?bào)表進(jìn)行分析，選取某日配合煤與焦炭的質(zhì)量指標(biāo)數(shù)據(jù)如表1所示。由表1可知，化驗(yàn)室3班均檢測(cè)配合煤的水分(Mt)、灰分(Ad)、揮發(fā)分(Vdaf)、硫分(St,d)、G值、細(xì)度及細(xì)度分布。對(duì)于焦炭質(zhì)量，化驗(yàn)室3班均檢測(cè)焦炭的水分(Mt)、灰分(Ad)、揮發(fā)分(Vdaf)、硫分(St,d)、平均粒度及粒度分布，而只有每天的白班才檢測(cè)焦炭的M40與M10，并由白班在干熄焦工段取樣檢測(cè)焦炭的CRI與CSR指標(biāo),分別為20.7%和69.9%。

表1 配合煤與焦炭質(zhì)量日?qǐng)?bào)表

考慮到配合煤細(xì)度分布可以由細(xì)度表示，因此本文采用配合煤的Mt，Ad，Vdaf，St,d和G值及細(xì)度指標(biāo)作為焦炭質(zhì)量預(yù)測(cè)模型的輸入變量。焦炭質(zhì)量指標(biāo)中的平均粒度是粒度分布的函數(shù)，且焦炭粒度受篩焦工藝影響，而焦炭的水分主要由熄焦工藝決定，因此選取焦炭的Ad，Vdaf，St,d，M40，M10，CRI及CSR指標(biāo)作為焦炭質(zhì)量預(yù)測(cè)模型的輸出變量。

基于上述質(zhì)量日?qǐng)?bào)表中配合煤與焦炭質(zhì)量指標(biāo)的選取，本文收集了湖南某焦化企業(yè)2018年全年的配合煤與焦炭質(zhì)量數(shù)據(jù)，并對(duì)數(shù)據(jù)進(jìn)行圖表分析，部分指標(biāo)的分析結(jié)果如圖1和圖2所示。由圖1可知，配合煤各質(zhì)量指標(biāo)的平均值均在要求范圍內(nèi)，其中水分在年初出現(xiàn)了超過上限值的情況，這可能與南方雨季空氣濕度高有關(guān)，揮發(fā)分穩(wěn)定在限值范圍內(nèi)，G值和細(xì)度出現(xiàn)了較大范圍的波動(dòng)，且部分?jǐn)?shù)據(jù)超出限值。在圖2中，焦炭的冷態(tài)強(qiáng)度(M40，M10)及熱態(tài)強(qiáng)度(CSR，CRI)指標(biāo)均在限值范圍內(nèi)。

圖1 配合煤質(zhì)量指標(biāo)統(tǒng)計(jì)圖

圖2 焦炭質(zhì)量指標(biāo)統(tǒng)計(jì)圖

煉焦工況數(shù)據(jù)由工況日?qǐng)?bào)表提供，本文選取對(duì)焦炭質(zhì)量有影響的最短結(jié)焦時(shí)間(Ctime)、機(jī)側(cè)直行溫度(CMT)、焦側(cè)直行溫度(CCT)、直行溫度平均系數(shù)K均(KM)及直行溫度安定系數(shù)K安(KA)等工況數(shù)據(jù)進(jìn)行分析，如表2所示。最短結(jié)焦時(shí)間、機(jī)側(cè)直行溫度、焦側(cè)直行溫度按班組統(tǒng)計(jì)，K均與K安按天統(tǒng)計(jì)。

根據(jù)上述工況日?qǐng)?bào)表，經(jīng)剔除無效、未記錄數(shù)據(jù)后的湖南某焦化企業(yè)2018年全年工況數(shù)據(jù)統(tǒng)計(jì)圖如圖3所示。由圖3可知，大部分工況數(shù)據(jù)滿足生產(chǎn)要求，但仍存在超出規(guī)范值的情況，其中最短結(jié)焦時(shí)間在一段時(shí)間內(nèi)稍微低于要求值，企業(yè)在該時(shí)段內(nèi)可能處于強(qiáng)化生產(chǎn)過程，與此同時(shí)機(jī)、焦側(cè)直行溫度先超出上限值再回到平均值附近，這說明現(xiàn)場(chǎng)操作人員針對(duì)強(qiáng)化生產(chǎn)的影響修改了加熱制度。從K均曲線可以看出，雖然其平均值可以滿足要求，但仍存在低于下限值的情況，它反映了焦?fàn)t測(cè)溫火道平均溫度的均勻性。

表2 工況日?qǐng)?bào)表

圖3 工況數(shù)據(jù)統(tǒng)計(jì)圖

2 數(shù)據(jù)匹配關(guān)系分析

為了分析配合煤、焦炭及工況數(shù)據(jù)的匹配關(guān)系，根據(jù)企業(yè)生產(chǎn)過程與管理要求，在圖4中給出了一個(gè)煉焦周轉(zhuǎn)周期所涉及的班組、物流位置、操作及數(shù)據(jù)等信息。配合煤經(jīng)粉碎后由皮帶送入煤塔，制樣人員在皮帶上取配合煤煤樣，化驗(yàn)室人員給出該煤樣的質(zhì)量數(shù)據(jù)；為保證生產(chǎn)的連續(xù)性，配合煤將在煤塔中儲(chǔ)存一個(gè)班組時(shí)間(8 h)，再經(jīng)搗固后裝煤；配合煤在炭化室內(nèi)經(jīng)過26 h的干餾煉焦，再經(jīng)過出焦、熄焦后送至篩焦樓；經(jīng)篩焦工藝處理的焦炭由皮帶送出，制樣人員在皮帶上取焦樣，化驗(yàn)室人員給出該焦樣的質(zhì)量數(shù)據(jù)。對(duì)于工況數(shù)據(jù)，在每個(gè)班組交班的時(shí)候給出當(dāng)前班組的工況信息，對(duì)于KM與KA等按日統(tǒng)計(jì)的數(shù)據(jù)，則由中班最后統(tǒng)計(jì)。

圖4 數(shù)據(jù)時(shí)序關(guān)系圖

2.1 配合煤與焦炭數(shù)據(jù)匹配關(guān)系

由圖4可知，晚班制樣的配合煤經(jīng)過一系列煉焦工序后，在次日中班可以得到其對(duì)應(yīng)的焦炭質(zhì)量指標(biāo)數(shù)據(jù)，由此可以認(rèn)為當(dāng)日晚班取樣的配合煤與次日中班取樣的焦炭相匹配。即焦炭質(zhì)量指標(biāo)數(shù)據(jù)與配合煤質(zhì)量指標(biāo)數(shù)據(jù)在時(shí)序上相差5個(gè)班組。在配煤專家的配煤經(jīng)驗(yàn)中，配合煤到焦炭的時(shí)間間隔一般取1.5個(gè)結(jié)焦時(shí)間。在本文中，結(jié)焦時(shí)間為26 h，因此分析得到的5個(gè)班組間隔與經(jīng)驗(yàn)值基本一致。在上述實(shí)際焦炭質(zhì)量指標(biāo)分析中，每班組均可得到焦炭的Ad，Vdaf，St和d，因此對(duì)于這3個(gè)數(shù)據(jù)可按間隔5個(gè)班組的時(shí)序找到相匹配的配合煤班組數(shù)據(jù)。焦炭的M40和M10僅由白班給出，因此也可按間隔5個(gè)班組的時(shí)序找到與其匹配的配合煤數(shù)據(jù)。在表1中，在干熄焦工段取樣檢測(cè)焦炭的CSR及CRI指標(biāo)，結(jié)合圖4可知焦炭的CSR及CRI指標(biāo)需按間隔4個(gè)班組的時(shí)序找到與之匹配的配合煤班組數(shù)據(jù)。

2.2 工況與焦炭數(shù)據(jù)匹配關(guān)系

在圖4中，晚班制樣的配合煤在炭化室內(nèi)受到當(dāng)日白班、中班及次日晚班、白班工況的影響，即中班獲取的焦炭質(zhì)量數(shù)據(jù)與當(dāng)日白班、晚班及前一日中班、白班的工況數(shù)據(jù)相匹配。在表2給出的工況數(shù)據(jù)中，機(jī)、焦側(cè)直行溫度每日3班均可獲得，且直行溫度影響整個(gè)煉焦過程，因此可按上述匹配關(guān)系進(jìn)行數(shù)據(jù)匹配。最短結(jié)焦時(shí)間由每個(gè)班組統(tǒng)計(jì)，由圖4可知，當(dāng)日中班獲取的焦炭質(zhì)量數(shù)據(jù)與當(dāng)日白班統(tǒng)計(jì)的最短結(jié)焦時(shí)間相匹配。K均和K安系數(shù)按日統(tǒng)計(jì)，因此在圖4中，當(dāng)日中班獲取的焦炭質(zhì)量數(shù)據(jù)與當(dāng)日及前一日的K均和K安系數(shù)匹配，而由當(dāng)日白班或晚班獲取的焦炭質(zhì)量數(shù)據(jù)僅與前一日的K均和K安系數(shù)匹配。

2.3 匹配數(shù)據(jù)組及PCA降維

基于上述數(shù)據(jù)匹配關(guān)系，本文將圖2和圖3所示的某焦化企業(yè)2018年實(shí)際煉焦生產(chǎn)數(shù)據(jù)分為3個(gè)匹配數(shù)據(jù)組，如表3所示。由于存在圖4所示的時(shí)序關(guān)系，并剔除了部分記錄不完整的數(shù)據(jù)組，因此表3給出的數(shù)據(jù)規(guī)模稍少于圖2與圖3。此外，對(duì)于工況指標(biāo)數(shù)據(jù)，考慮到焦炭質(zhì)量受多個(gè)班組工況數(shù)據(jù)的共同影響，且圖3所示的實(shí)際煉焦生產(chǎn)過程工況相對(duì)穩(wěn)定，因此取多個(gè)班組工況數(shù)據(jù)的平均值具有實(shí)際意義。

由表3可知，在匹配數(shù)據(jù)組中配合煤指標(biāo)包含6項(xiàng)原始數(shù)據(jù)指標(biāo)，工況指標(biāo)包含5項(xiàng)數(shù)據(jù)指標(biāo)。為降低數(shù)據(jù)指標(biāo)的維度，對(duì)原始指標(biāo)數(shù)據(jù)進(jìn)行主成分分析(PCA)，分析結(jié)果如表4所示?；谥鞒煞挚傌暙I(xiàn)率大于90%的原則，選取經(jīng)PCA分析后的配合煤指標(biāo)3個(gè)主成分，分別命名為M1，M2和M3，各主成分的貢獻(xiàn)率如表4所示。對(duì)于工況數(shù)據(jù)，經(jīng)PCA分析后，其兩個(gè)主成分的總貢獻(xiàn)率達(dá)到99%，因此選取工況指標(biāo)主成分C1與C2。

表3 匹配數(shù)據(jù)組

表4 數(shù)據(jù)組的PCA分析

3 焦炭質(zhì)量預(yù)測(cè)

依據(jù)上述實(shí)際煉焦生產(chǎn)數(shù)據(jù)的分組與匹配關(guān)系，本文采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建焦炭質(zhì)量預(yù)測(cè)模型。BP神經(jīng)網(wǎng)絡(luò)具有逼近非線性函數(shù)的能力，并已廣泛應(yīng)用于焦炭質(zhì)量預(yù)測(cè)算法的設(shè)計(jì)。針對(duì)表4中數(shù)據(jù)組C_Base，C_M40M10和C_CRICSR分別構(gòu)建BP神經(jīng)網(wǎng)絡(luò)。在MATLAB軟件中，數(shù)據(jù)經(jīng)標(biāo)幺化操作后，采用標(biāo)準(zhǔn)函數(shù)newff構(gòu)建3層神經(jīng)網(wǎng)絡(luò)，其中輸入層神經(jīng)元個(gè)數(shù)為5個(gè)，分別對(duì)應(yīng)表4中配合煤主成分(M1，M2和M3)及工況主成分(C1和C2)，中間層神經(jīng)元個(gè)數(shù)為9個(gè)，輸出層神經(jīng)元與焦炭質(zhì)量指標(biāo)相對(duì)應(yīng)。利用train函數(shù)訓(xùn)練網(wǎng)絡(luò)，最后由sim函數(shù)進(jìn)行預(yù)測(cè)計(jì)算。選取數(shù)據(jù)組的前2/3部分作為訓(xùn)練樣本，后1/3部分作為測(cè)試樣本。為了對(duì)比分析實(shí)際煉焦工況數(shù)據(jù)對(duì)焦炭質(zhì)量預(yù)測(cè)精度的影響，對(duì)各數(shù)據(jù)組分別構(gòu)建了不考慮工況的神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

基于上述實(shí)際煉焦生產(chǎn)數(shù)據(jù)組及BP神經(jīng)網(wǎng)絡(luò)，本文進(jìn)行了考慮工況與不考慮工況的焦炭質(zhì)量預(yù)測(cè)分析，各焦炭質(zhì)量指標(biāo)的預(yù)測(cè)誤差如圖5所示。由圖5可知，考慮工況與不考慮工況對(duì)應(yīng)的預(yù)測(cè)誤差曲線幾乎一致，這說明穩(wěn)定、順行、合規(guī)的工況對(duì)焦炭質(zhì)量預(yù)測(cè)的影響較小。而無論是否考慮工況，焦炭灰分的預(yù)測(cè)誤差可以控制在±4%以內(nèi)，揮發(fā)分的預(yù)測(cè)誤差集中在±5%以內(nèi)，硫分的預(yù)測(cè)誤差大部分在±5%以內(nèi)，M40的預(yù)測(cè)誤差控制在±3%以內(nèi)，M10的預(yù)測(cè)誤差集中在±10%以內(nèi)，CRI的預(yù)測(cè)誤差大部分在±10%以內(nèi)，CSR的預(yù)測(cè)誤差在±5%以內(nèi)。

圖5 焦炭質(zhì)量預(yù)測(cè)誤差分析

表5給出了上述各焦炭質(zhì)量指標(biāo)預(yù)測(cè)誤差絕對(duì)值的平均值。由表5可知，除了CRI指標(biāo)，其他指標(biāo)的預(yù)測(cè)誤差絕對(duì)值的平均值均小于5%。工況參與預(yù)測(cè)后，焦炭質(zhì)量指標(biāo)CRI與CSR的預(yù)測(cè)精度得到了提高。

表5 絕對(duì)誤差平均值

進(jìn)一步在表6中給出了上述預(yù)測(cè)計(jì)算的誤差均方根，可以看出在預(yù)測(cè)模型中考慮實(shí)際煉焦工況可以提高焦炭質(zhì)量的預(yù)測(cè)精度，且部分指標(biāo)優(yōu)于文獻(xiàn)[5]的計(jì)算結(jié)果，這說明精確的數(shù)據(jù)匹配可以有效地提高焦炭質(zhì)量的預(yù)測(cè)精度，并彌補(bǔ)預(yù)測(cè)算法的性能。

表6 誤差均方根

4 結(jié)論

在實(shí)際煉焦生產(chǎn)過程中，焦炭質(zhì)量的預(yù)測(cè)精度取決于對(duì)煉焦生產(chǎn)數(shù)據(jù)的處理以及預(yù)測(cè)方法的設(shè)計(jì)，本文分析了實(shí)際焦化企業(yè)的報(bào)表數(shù)據(jù)，完成了數(shù)據(jù)的分析、提取、匹配及PCA降維，并結(jié)合BP神經(jīng)網(wǎng)絡(luò)，構(gòu)建了符合實(shí)際煉焦生產(chǎn)過程的焦炭質(zhì)量預(yù)測(cè)模型，計(jì)算結(jié)果表明穩(wěn)定、順行、合規(guī)的煉焦工況對(duì)焦炭質(zhì)量預(yù)測(cè)精度的影響較小，而精確的數(shù)據(jù)匹配有利于提高預(yù)測(cè)精度，并可降低對(duì)預(yù)測(cè)方法性能的要求。此外，本文對(duì)實(shí)際煉焦生產(chǎn)數(shù)據(jù)的分析以及焦炭質(zhì)量預(yù)測(cè)方法的設(shè)計(jì)是工業(yè)大數(shù)據(jù)分析與數(shù)據(jù)挖掘的一種嘗試，并可為其他流程化工業(yè)生產(chǎn)過程的數(shù)據(jù)分析與挖掘提供借鑒。