本發(fā)明涉及鋼鐵燒結過程生產(chǎn)節(jié)能降耗領域，具體涉及一種建立綜合焦比預測模型和預測綜合焦比的方法。
背景技術：
：鋼鐵工業(yè)是國民經(jīng)濟支柱產(chǎn)業(yè)之一，鋼鐵工業(yè)的發(fā)展也將決定我國國民經(jīng)濟的發(fā)展。隨著鋼鐵行業(yè)的快速發(fā)展，其所帶來的能源消耗以及環(huán)境問題也日益突顯。燒結過程又是煉鐵過程的一個重要環(huán)節(jié)，其過程生產(chǎn)的燒結礦是高爐煉鐵的主要原料，并且燒結過程也是鋼鐵生產(chǎn)過程除高爐煉鐵過程外最大的耗能工序，其能耗約占鋼鐵冶金總能耗的10％～15％，在能源消耗中，主要有焦粉燃燒、燃氣點火和電能消耗，其中焦粉燃料消耗約為80％，燃氣消耗約為6％，電能消耗約為13.5％，其他約為0.5％。隨著近幾年來全球生態(tài)環(huán)境的持續(xù)惡化，各國也越來越重視對環(huán)境的保護，“可持續(xù)發(fā)展”、“綠色制造”和“低碳經(jīng)濟”等經(jīng)濟發(fā)展觀念被越來越多的國家接受，并融入到各國的經(jīng)濟發(fā)展之中。由此可見，通過提高燒結過程碳能源利用率(即提高碳效)，將是實現(xiàn)我國鋼鐵行業(yè)生產(chǎn)過程節(jié)能減排的重要途徑之一。目前，國際上著名鋼鐵企業(yè)所采用的鋼鐵生產(chǎn)方式是帶式抽風燒結方式，其生產(chǎn)過程主要是先把原料充分混合得到混合料，然后將混合料平鋪在臺車上，混合料在點火爐處點燃料層表面，料層隨臺車移動而移動，此時臺車下方的風箱開始進行負壓抽風，料層將自上而下地進行燃燒，直到在燒結終點處燒穿料層，最終完成燒結造塊過程。燒結過程的工藝流程見附圖1。但帶式抽風燒結方式在燒結生產(chǎn)中，充分混合的混合料在料層中燃燒會產(chǎn)生1300℃左右的高溫，使得燒結混合料會在這樣的高溫環(huán)境下會發(fā)生一系列的化學物理變化，這樣會使混合料層出現(xiàn)分層現(xiàn)象，其中料層自上而下可分為燒結礦層、燃燒層、預熱干燥層、過濕層和生料層等，其中燒結料層的分層現(xiàn)象見附圖2。另外，燒結生產(chǎn)過程中，焦粉作為燒結過程主要的能量來源，無法在保證燒結礦質量和產(chǎn)量的前提下提高焦粉利用率，而且，燒結過程涉及到配料、制粒、布料、燒結點火和燒結終點等等工序，并且該過程具有非線性、強耦合、機理復雜和工藝流程長等特征。因此，為了保證燒結礦的質量和產(chǎn)量，需要對燒結過程碳效指標進行準確有效預測。技術實現(xiàn)要素：有鑒于此，本發(fā)明的實施例提供了一種能對燒結過程碳效指標進行準確預測的建立綜合焦比預測模型和預測綜合焦比的方法。本發(fā)明的實施例提供一種建立綜合焦比預測模型和預測綜合焦比的方法，包括以下步驟：(1)選定綜合焦比作為衡量燒結過程碳效的指標，并計算綜合焦比；(2)確定影響綜合焦比的燒結參數(shù)，對燒結參數(shù)的生產(chǎn)歷史數(shù)據(jù)進行時序配準和平均值濾波處理，得到樣本數(shù)據(jù)庫；(3)對步驟(2)得到的樣本數(shù)據(jù)庫采用主成分分析法對燒結參數(shù)進行數(shù)據(jù)降維，并重新組合，得到主成分變量；(4)以步驟(3)得的主成分變量作為輸入變量，以步驟(1)得的綜合焦比作為輸出變量進行重復計算和驗證，建立綜合焦比預測模型；(5)分析待燒結礦的燒結參數(shù)，對待燒結礦的燒結參數(shù)進行數(shù)據(jù)降維，并重新組合，得到待燒結礦的主成分變量；(6)將步驟(5)得到的待燒結礦的主成分變量輸入步驟(4)建立的綜合焦比預測模型，綜合焦比預測模型的輸出變量即為待燒結礦燒結過程的綜合焦比。進一步，所述步驟(1)中，綜合焦比是生產(chǎn)每噸鐵水消耗的燃料所折合成的焦炭數(shù)量，綜合焦比的的計算式如下：式中：η表示綜合焦比(Kg/t)，wc表示焦粉配比(％)，p表示成品率(％)，q表示燒損率(％)，QD表示燒結的大成礦產(chǎn)量(Kg/h)，QX表示燒結的小成礦產(chǎn)量(Kg/h)，QF表示燒結的返礦量(Kg/h)，QP表示燒結的鋪底料量(Kg/h)。進一步，所述步驟(2)中，燒結參數(shù)通過分析燒結過程機理得到，燒結參數(shù)包括垂直燃燒速度、上升點、上升點溫度、燒結終點、燒結終點溫度、返礦、風箱負壓、料層厚度、臺車速度和焦粉配比。進一步，所述步驟(2)中，燒結過程存在時滯，為了確保燒結參數(shù)在時序上一致，對燒結參數(shù)的生產(chǎn)歷史數(shù)據(jù)進行時序配準。進一步，所述步驟(2)中，采用燒結過程中大成礦的波動周期作為采樣周期對燒結參數(shù)的生產(chǎn)歷史數(shù)據(jù)進行平均值濾波處理，得到樣本數(shù)據(jù)庫。進一步，所述燒結過程中大成礦的波動周期為45min。進一步，所述步驟(3)中，數(shù)據(jù)降維包括以下步驟：(3.1)對樣本數(shù)據(jù)庫進行標準化處理；設樣本數(shù)據(jù)庫有n組樣本數(shù)據(jù)，每組數(shù)據(jù)中有u個判別指標，樣本數(shù)據(jù)庫用矩陣X＝(xij)n×u表示，xij為第j個判別指標的第i組樣本數(shù)據(jù)，標準化公式如下：式中：Xscalar是樣本數(shù)據(jù)標準化后得到的數(shù)據(jù)，Xmin是樣本數(shù)據(jù)中的最小值，Xmax是樣本數(shù)據(jù)中的最大值，Xactual是樣本數(shù)據(jù)中的實際值；(3.2)計算樣本數(shù)據(jù)的相關系數(shù)rij和相關系數(shù)矩陣C；計算式如下：C＝(rij)u×u(i＝1,2,…,u；j＝1,2,…,u)式中：xai為第i個判別指標第a組樣本數(shù)據(jù)；為第i個判別指標所有樣本數(shù)據(jù)的平均值；xaj為第j個判別指標第a組樣本數(shù)據(jù)；為第j個判別指標所有樣本數(shù)據(jù)的平均值；(3.3)計算相關系數(shù)矩陣C的u個特征值，記為λ1≥λ2≥…≥λu≥0；(3.4)計算每個燒結參數(shù)的貢獻率確定主成分變量的個數(shù)；設第d個燒結參數(shù)的貢獻率為以作為前m個燒結參數(shù)的累積貢獻率，以累積貢獻率達到85％及以上時，燒結參數(shù)的個數(shù)確定為主成分變量的個數(shù)；(3.5)設m個燒結參數(shù)的累積貢獻率達到85％及以上，則主成分變量個數(shù)為m個，可得到主成分變量為：式中：yd表示第d個主成分變量，ld表示m個特征值所對應的單位特征向量，(d＝1,2,…,m)。進一步，所述步驟(4)中，根據(jù)最小二乘支持向量機法建立綜合焦比預測模型。進一步，所述最小二乘支持向量機法建立綜合焦比預測模型包括以下步驟：(4.1)設主成分變量和綜合焦比組成的樣本數(shù)據(jù)集為yk∈Rm為m維輸入變量，Yk∈R為對應的目標輸出，即綜合焦比；(4.2)采用非線性映射將輸入變量映射到高維特征空間，得到：式中：ω是權函數(shù)，表示將輸入變量非線性映射到高維空間，b是偏置量；(4.3)基于結構風險最小化原理，將上述回歸問題轉化為約束優(yōu)化問題：式中：γ是懲罰因子，其是衡量參數(shù)風險和擬合誤差相對重要程度的參數(shù)，ek是擬合誤差；(4.4)對(4.3)中的約束優(yōu)化問題采用拉格朗日乘子法、正定矩陣和卡羅需-庫恩-塔克條件進行求解，得到綜合焦比預測模型為：式中：αk是拉格朗日乘子，K(y,yk)＝exp(-||y-yk||2/2σ2)為高斯徑向基核函數(shù)，σ是核函數(shù)寬度。與現(xiàn)有技術相比，本發(fā)明具有以下有益效果：(1)本發(fā)明確定綜合焦比作為衡量燒結過程碳效的指標，并給出其相關計算，為實現(xiàn)燒結過程節(jié)能降耗奠定基礎；(2)本發(fā)明通過燒結過程機理分析，確定影響綜合焦比的燒結參數(shù)，采用主成分分析法對這些燒結參數(shù)進行降維重組獲得主成分變量，可有效克服燒結參數(shù)之間的耦合問題；(3)本發(fā)明基于主成分變量，采用最小二乘支持向量機法，建立綜合焦比預測模型，能夠反映燒結過程的非線性特性，有效保證預測模型的精確和合理；(4)本發(fā)明基于燒結過程的生產(chǎn)歷史數(shù)據(jù)，進行模型的仿真實驗，能夠在實際生產(chǎn)過程中廣泛應用。附圖說明圖1是帶式抽風燒結方式的燒結工藝流程圖。圖2是帶式抽風燒結方式的燒結料層分層現(xiàn)象的示意圖。圖3是本發(fā)明一實施例的流程圖。圖4是本發(fā)明實施例的綜合焦比預測結果與實際值的對比圖。圖5是本發(fā)明實施例的綜合焦比預測相對誤差的結果圖。具體實施方式下面結合附圖和具體實施例對本發(fā)明做進一步說明。實施例請參照附圖3，本實施例包括以下步驟：(1)綜合焦比是生產(chǎn)每噸鐵水消耗的燃料所折合成的焦炭數(shù)量，提高焦炭利用率是保證燒結礦產(chǎn)量和質量的前提，因此，選定綜合焦比作為衡量燒結過程碳效的指標，并計算綜合焦比，綜合焦比的的計算式如下：式中：η表示綜合焦比(Kg/t)，wc表示焦粉配比(％)，p表示成品率(％)，q表示燒損率(％)，QD表示燒結的大成礦產(chǎn)量(Kg/h)，QX表示燒結的小成礦產(chǎn)量(Kg/h)，QF表示燒結的返礦量(Kg/h)，QP表示燒結的鋪底料量(Kg/h)。(2)通過分析燒結過程機理得到影響綜合焦比的燒結參數(shù)，燒結參數(shù)包括垂直燃燒速度、上升點(BRP)、BRP溫度、燒結終點(BTP)、BTP溫度、返礦、風箱負壓、料層厚度、臺車速度和焦粉配比，垂直燃燒速度、BRP、BRP溫度、BTP、BTP溫度、返礦、風箱負壓、料層厚度、臺車速度和焦粉配比的生產(chǎn)歷史數(shù)據(jù)以日報表的形式保存在操作室工控機的本地數(shù)據(jù)庫中，根據(jù)日報表的數(shù)據(jù)，收集影響綜合焦比的垂直燃燒速度、BRP、BRP溫度、BTP、BTP溫度、返礦、風箱負壓、料層厚度、臺車速度和焦粉配比一個月的生產(chǎn)歷史數(shù)據(jù)，組成原始樣本數(shù)據(jù)；燒結過程中的各個工序之間相互影響，從配料生產(chǎn)到燒結礦的破碎完成大約需要2個小時，對燒結礦化驗需要2個小時，從此可見燒結生產(chǎn)過程存在時滯現(xiàn)象，為了確保燒結生產(chǎn)過程燒結參數(shù)在時序上保持一致，需對燒結參數(shù)的生產(chǎn)歷史數(shù)據(jù)進行時序配準處理；根據(jù)燒結工藝規(guī)定，大成礦每隔一段時間會分流到鋪底料倉，大成礦呈現(xiàn)周期性變化，周期大約為45min，由于大成礦的周期性波動會造成綜合焦比的波動，故采用大成礦的波動周期作為采樣周期，選定采樣周期為45min對燒結參數(shù)的生產(chǎn)歷史數(shù)據(jù)進行平均值濾波處理，得到燒結參數(shù)垂直燃燒速度、BRP、BRP溫度、BTP、BTP溫度、返礦、風箱負壓、料層厚度、臺車速度和焦粉配比的生產(chǎn)歷史數(shù)據(jù)共900組，以此，建立樣本數(shù)據(jù)庫；(3)對步驟(2)得到的樣本數(shù)據(jù)庫采用主成分分析法對燒結參數(shù)進行數(shù)據(jù)降維，并重新組合，得到主成分變量；數(shù)據(jù)降維包括以下步驟：(3.1)對樣本數(shù)據(jù)庫進行標準化處理；設樣本數(shù)據(jù)庫有n組樣本數(shù)據(jù)，每組數(shù)據(jù)中有u個判別指標，樣本數(shù)據(jù)庫用矩陣X＝(xij)n×u表示，xij為第j個判別指標的第i組樣本數(shù)據(jù)，標準化公式如下：式中：Xscalar是樣本數(shù)據(jù)標準化后得到的數(shù)據(jù)，Xmin是樣本數(shù)據(jù)中的最小值，Xmax是樣本數(shù)據(jù)中的最大值，Xactual是樣本數(shù)據(jù)中的實際值；(3.2)計算樣本數(shù)據(jù)的相關系數(shù)rij和相關系數(shù)矩陣C；計算式如下：C＝(rij)u×u(i＝1,2,…,u；j＝1,2,…,u)式中：xai為第i個判別指標第a組樣本數(shù)據(jù)；為第i個判別指標所有樣本數(shù)據(jù)的平均值；xaj為第j個判別指標第a組樣本數(shù)據(jù)；為第j個判別指標所有樣本數(shù)據(jù)的平均值；(3.3)計算相關系數(shù)矩陣C的u個特征值，記為λ1≥λ2≥…≥λu≥0；(3.4)計算每個燒結參數(shù)的貢獻率確定主成分變量的個數(shù)；設第d個燒結參數(shù)的貢獻率為以作為前m個燒結參數(shù)的累積貢獻率，以累積貢獻率達到85％及以上時，燒結參數(shù)的個數(shù)確定為主成分變量的個數(shù)；(3.5)設m個燒結參數(shù)的累積貢獻率達到85％及以上，則主成分變量個數(shù)為m個，可得到主成分變量為：式中：yd表示第d個主成分變量，ld表示m個特征值所對應的單位特征向量，(d＝1,2,…,m)。據(jù)此，從數(shù)據(jù)樣本庫的900組數(shù)據(jù)中隨機選取800組數(shù)據(jù)采用主成分分析法對燒結參數(shù)進行數(shù)據(jù)降維，并重新組合，得到數(shù)個互不相關的主成分變量，分析結果如表1所示；表1樣本數(shù)據(jù)庫燒結參數(shù)的主成分變量分析結果成分λ(特征值)貢獻率/％累積貢獻率/％垂直燃燒速度3.526434.4834.48BRP1.844018.0352.51BRP溫度1.226311.9964.5BTP1.093210.6975.19BTP溫度0.89648.7783.96返礦0.75477.3891.34風箱負壓0.43804.2895.62料層厚度0.36703.5999.21臺車速度0.07250.7199.92焦粉配比0.00760.08100由表1可以得出，樣本數(shù)據(jù)庫的前6個燒結參數(shù)累積貢獻概率已達到91.34％，則可用前6個燒結參數(shù)作為主成分變量來表示原樣本數(shù)據(jù)信息，并且這些主成分變量能夠解決影響綜合焦比的燒結參數(shù)之間的耦合問題，再求得燒結參數(shù)所對應的6個主成分變量所對應的特征向量為：l1＝[0.70-0.39-0.26-0.02-0.180.370.260.15-0.08-0.12]l2＝[0.380.020.040.100.21-0.20-0.18-0.12-0.340.77]l3＝[0.300.090.420.490.34-0.04-0.300.19-0.16-0.46]l4＝[0.03-0.30-0.25-0.20-0.34-0.24-0.77-0.03-0.13-0.18]l5＝[-0.020.37-0.36-0.290.420.45-0.310.410.010.06]l6＝[-0.37-0.360.240.34-0.260.38-0.130.47-0.030.33]則可計算影響綜合焦比的燒結參數(shù)的主成分變量為：其中，tc(c＝1,2,…,10)表示樣本數(shù)據(jù)庫中垂直燃燒速度、BRP、BRP溫度、BTP、BTP溫度、返礦、風箱負壓、料層厚度、臺車速度和焦粉配比的數(shù)據(jù)；(4)以步驟(3)得的6個主成分變量作為輸入變量，以步驟(1)得的綜合焦比作為輸出變量進行重復計算，根據(jù)最小二乘支持向量機法建立綜合焦比預測模型，并用樣本數(shù)據(jù)庫中剩余的100組數(shù)據(jù)進行綜合焦比預測模型驗證；最小二乘支持向量機法是一種統(tǒng)計學的建模方法，其原理是通過將低維空間中線性不可分問題，采用核函數(shù)映射到高維空間，從而使其線性可分，采用最小二乘對高維空間的最優(yōu)分類面進行求解，把不等式約束變?yōu)榈仁郊s束，這樣就可以使用該方法進行綜合焦比預測；包括以下步驟：(4.1)設主成分變量和綜合焦比組成的樣本數(shù)據(jù)集為yk∈Rm為m維輸入變量，Yk∈R為對應的目標輸出，即綜合焦比；(4.2)采用非線性映射將輸入變量映射到高維特征空間，得到：式中：ω是權函數(shù)，表示將輸入變量非線性映射到高維空間，b是偏置量；(4.3)基于結構風險最小化原理，將上述回歸問題轉化為約束優(yōu)化問題：式中：γ是懲罰因子，其是衡量參數(shù)風險和擬合誤差相對重要程度的參數(shù)，ek是擬合誤差；(4.4)對(4.3)中的約束優(yōu)化問題采用拉格朗日乘子法、正定矩陣和卡羅需-庫恩-塔克條件進行求解，得到綜合焦比預測模型為：式中：αk是拉格朗日乘子，K(y,yk)＝exp(-||y-yk||2/2σ2)為高斯徑向基核函數(shù)，σ是核函數(shù)寬度。將預測結果和實際結果進行對比，對比圖見附圖4，預測的相對誤差見附圖5，由圖可知，綜合焦比預測結果的相對誤差在[-3％,2％]之內(nèi)，因此，建立的綜合焦比預測模型具有可行性；(5)分析待燒結礦的燒結參數(shù)，對待燒結礦的燒結參數(shù)進行數(shù)據(jù)降維，并重新組合，得到待燒結礦的主成分變量；(6)將步驟(5)得到的待燒結礦的主成分變量輸入步驟(4)建立的綜合焦比預測模型，綜合焦比預測模型的輸出變量即為待燒結礦燒結過程的綜合焦比。本方法預測精度高，能夠滿足實際燒結過程生產(chǎn)要求，可為燒結過程的碳效優(yōu)化奠定基礎。在不沖突的情況下，本文中上述實施例及實施例中的特征可以相互結合。以上所述僅為本發(fā)明的較佳實施例，并不用以限制本發(fā)明，凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進等，均應包含在本發(fā)明的保護范圍之內(nèi)。當前第1頁1 2 3