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