水泥生料調(diào)配運行工況識別

1、水泥生料調(diào)配運行工況識別1、數(shù)據(jù)處理1.1數(shù)據(jù)濾波水泥生料配料過程主要是依據(jù)在線分析儀與熒光分析儀兩種分析儀的數(shù)據(jù)來改變配比使得生料質(zhì)量穩(wěn)定。由于水泥生料的制備過程中由于在線分析儀檢測模型問題引起檢測到的水泥生料成分?jǐn)?shù)據(jù)會有一定的噪聲存在，并且這些噪聲何時發(fā)生并不能提前預(yù)知，因此產(chǎn)生了隨機誤差，僅根據(jù)原始數(shù)據(jù)進行控制會生產(chǎn)出大量的質(zhì)量不合格的水泥生料。這些隨機誤差是不可避免且數(shù)值不正確的，因此要采用合適的數(shù)字濾波的方式來消除。經(jīng)過一系列的仿真實驗，證明限幅濾波可以把錯誤信息剔除保留有用信息，因此該方法對數(shù)據(jù)的處理能夠滿足工況識別要求。以對水泥生料成分中的硅為例，在線分析儀對硅的成分檢測會在一定

2、范圍內(nèi)波動，其干擾一般為隨機干擾。因此采用限幅濾波法進行數(shù)據(jù)預(yù)處理。限幅濾波是對連續(xù)采樣值進行判斷，如若采樣值在允許最大偏差范圍內(nèi)則保留，若大于最大允許最大范圍則采用上一個采樣值進行代替，最后再對數(shù)據(jù)滾動求平均。由于實際生產(chǎn)過程中si的波動較小一般在11.3%-12.3%之間波動，僅少數(shù)情況下較低或較高，大于13%或者小于10.5%的數(shù)則不符合實際情況為錯誤值，在多次嘗試后濾波上下限設(shè)在10.5%,13%，并且si的擾動較大滯后時間較長根據(jù)實驗和現(xiàn)場應(yīng)用得出濾波均值平滑時間采用十五分鐘滾動求平均較好，平滑公式如下：由圖2可以看出Ca的擾動較小波動較為穩(wěn)定一般在40%上下波動經(jīng)過實驗以及現(xiàn)場應(yīng)用

3、得出Ca含量的在線濾波上下限為42%,44%，均值濾波時間為5分鐘，滾動進行。平滑公式為：由于Al,Fe的成分含量很低并且在斷料時會出現(xiàn)為的零況發(fā)生，且?guī)缀鯖]有噪聲存在，因此不需要對這兩種成分進行限幅，僅做了平滑濾波，時間為五分鐘樣本數(shù)圖1在線分析儀S數(shù)據(jù)濾波曲j2111111fB-L.111-樣即UJVU刪沖習(xí)木圖3在線分析儀Fe數(shù)據(jù)濾波圖4在線分析儀Al數(shù)據(jù)濾波由以上濾波效果圖以及實際應(yīng)用結(jié)果來看在水泥生料配料過程中采用限幅濾波的方法處理現(xiàn)場數(shù)據(jù)效果較好。1.2在線分析儀數(shù)據(jù)修正本文所研究所依托的項目主要是依靠在線分析儀所測的的數(shù)據(jù)進行研究。由于熒光分析儀和在線分析儀兩種儀器采樣點相距很近

4、可以看作是在固定的時間內(nèi)檢測的水泥生料是同一個樣本。圖5所示的曲線為在線分析儀和熒光分析儀監(jiān)測樣本變化曲線，紅色的為在線分析儀檢測的成分變化，藍(lán)色曲線為熒光分析儀檢測的成分變化，從中可以看出盡管兩種儀器所測量的數(shù)值上有所差異，但是所反映水泥生料成分變化的趨勢是一致的。圖5熒光分析儀在線分析儀對水泥生料四種成分檢測對比圖化驗室在對水泥生料化驗時使用熒光分析儀進行化驗，給出的生料控制質(zhì)量指標(biāo)也是根據(jù)熒光分析儀化驗出的結(jié)果計算得出，因此要找出一種補償方法使在線分析儀的數(shù)據(jù)盡可能的和熒光分析儀一致。這樣才可以保證水泥生料的正常生產(chǎn)。經(jīng)過對兩種數(shù)據(jù)的對比研究最后得出了一種補償方法，控制軟件根據(jù)補償后的數(shù)

5、據(jù)來對進行控制。由于熒光分析儀為一個小時一個樣本所以在補償時是一個小時補償一次，補償公式如下：X=0.4*(xrt1(i)sf(i)-i-0,6*(xrf(j)sf(j)+sf其中x為補償后的數(shù)，為上一個小時熒光分析儀（連續(xù)/瞬時樣）結(jié)果，為上一個小時在線分析儀對應(yīng)的（連續(xù)/瞬時樣）結(jié)果；為當(dāng)前小時的時熒光分析儀（連續(xù)/瞬時樣）結(jié)果，sf（j）為當(dāng)前小時在線分析儀對應(yīng)的（連續(xù)/瞬時樣）結(jié)果，sf為在線分析儀實時檢測到的數(shù)?？梢詮南聢D圖6所示三條曲線綠色曲線（即最上面的曲線）為在線分析儀經(jīng)過濾波處理后的石灰石飽和系數(shù)原始數(shù)據(jù)，藍(lán)色曲線（中間曲線）為熒光分析儀檢測的石灰石飽和系數(shù)數(shù)據(jù)，紅色曲線（最

6、下面）為在線分析儀經(jīng)過補償修正后的石灰石飽和系數(shù)數(shù)圖6石灰石飽和系數(shù)熒光值在線分析儀值補償后值對比2調(diào)配運行工況2.1趨勢分析技術(shù)及其相關(guān)定義尺度的唯一性三種特性。Janusz禾口Venkatasubramanian在1991年7所示，具體包括趨勢分析指的趨勢是屬于過程的趨勢，反映特征值的變化過程即出發(fā)時刻，發(fā)展速度。趨勢應(yīng)有(1)存在性；(2)特征的穩(wěn)定性；(3)A,B,C,D,E,F,基于解釋過程趨勢的定性趨勢分析法是由提出的，并給出了趨勢的等級化的描述，如圖G七種基元。圖7趨勢變化七種基元圖生料配料過程趨勢分析對于水泥生料配料過程，要想實現(xiàn)水泥生料的自動化的生產(chǎn)，僅憑幾個主要的工藝參數(shù)進

7、行數(shù)字分析，或者對某些過程參數(shù)做簡單的判斷，并不能滿足水泥生料配料自動化生產(chǎn)的要求。對于生料配料這樣具有大滯后特性的生產(chǎn)環(huán)節(jié)，水泥生料三率值在絕大多數(shù)情況下在一定時間內(nèi)不會有急劇變化的情況發(fā)生。而在一定的時間內(nèi)，可以觀察石灰石飽和系數(shù)、硅率、鋁氧率的趨勢變化，這些趨勢特征可以為操作人員或水泥生料自動配料程序提供正確的參考信息。在水泥生料配料過程參數(shù)進行定量分析的基礎(chǔ)上，對大量的過程數(shù)據(jù)中提取有益于工況識別的信息，融入了趨勢分析法，以便生產(chǎn)高質(zhì)量的水泥生料。實現(xiàn)水泥生料配料生料三率值的趨勢變化的識別和提取還要考慮下面兩個問題：1. 噪聲影響在工業(yè)實際生產(chǎn)必然存在噪聲，近紅外線在線分析儀安裝在生料

8、運輸斜槽上，內(nèi)部干擾因素較多，也不容易判別，因此分析儀檢測到的數(shù)據(jù)存在較多的“毛刺”。噪聲的存在會使趨勢的特征提取造成一定影響，可能會使提取到的趨勢特征不能反映實際情況，因此對噪聲的處理要在使趨勢不失真的前提下進行。結(jié)合現(xiàn)場實際情況和仿真得出采用限幅濾波法然后進行均值濾波處理的方法較為合理。2. 時間窗口寬度時間窗的寬度是指趨勢分析中需要用到的每個變量的時間窗寬度。比較大的時間窗寬度能夠反映整體的變化趨勢，但時間窗過大數(shù)據(jù)量就會多使得計算復(fù)雜，而且可能一些具有意義的趨勢會被忽略。水泥生料配料反應(yīng)過程較長，生料三率值在短時間內(nèi)變化較為緩慢，因此時間窗的長度選取至為重要。根據(jù)水泥生料配料過程的工藝

9、特性和總結(jié)人工經(jīng)驗得出由于石灰石飽和系數(shù)的趨勢變化較為頻繁，而硅率和鋁氧率則變化緩慢長期較為穩(wěn)定，本文所依托的水泥生料生產(chǎn)線采用輥壓磨經(jīng)過對數(shù)據(jù)分析和現(xiàn)場總結(jié)得知在配比改變五分鐘左右開始變化十多分鐘后反應(yīng)結(jié)束，因此對水泥生料三率值系數(shù)趨勢識別采用的時間窗寬度為五分鐘。而且五分鐘能夠表征整體變化趨勢，還能夠?qū)ν蛔兊那闆r進行表征。本文主要是對水泥生料三率值的趨勢提取，即將三率值的定性趨勢從含有噪聲的過程數(shù)據(jù)中準(zhǔn)確無誤的提取出來，并用定義出趨勢。趨勢提取的方法包括基于多尺度趨勢提取方法，多項式擬合趨勢提取方法，知識庫法等。定性趨勢包含3個方面的特性：(1) 存在性，即這個趨勢是存在的，對于實際水泥生

10、料制備過程中，趨勢是在線分析儀檢測出的水泥生料成份變化的描述。(2) 特征的穩(wěn)定性，在較小的擾動干擾下，趨勢應(yīng)該是相對穩(wěn)定的狀態(tài)。如若沒有穩(wěn)定性，那么這個過程只能被認(rèn)為是一個隨機的過程。(3) 尺度唯一性，趨勢分析的結(jié)果由觀測尺度及數(shù)據(jù)空間組織的不同而不同。但如果觀測尺度固定，時間間隔也固定那么在固定的時間內(nèi)變量的趨勢是唯一的。圖8為現(xiàn)場實際運行數(shù)據(jù)，可以看出實際運行中硅率和鐵率的波動情況很小基本趨于穩(wěn)定易于調(diào)節(jié)；主要是石灰石飽和系數(shù)波動較大，其原因是因為石灰石飽和系數(shù)允許波動誤差為士0.02所以波動較大。實際生產(chǎn)中要結(jié)合其趨勢和指標(biāo)的相對位臵進行調(diào)節(jié)，經(jīng)驗證建立模型的方法無法滿足生產(chǎn)要求，在

11、水泥生料生產(chǎn)過程中結(jié)合操作員的人工經(jīng)驗以及實際情況，得出水泥生料配料過程適合采用特征知識庫這一方法對趨勢進行識別文獻針對生料控制系統(tǒng)建立了專家控制系統(tǒng)，但是僅是根據(jù)三率值與指標(biāo)的差來做出調(diào)節(jié)方案。但是僅根據(jù)這一點進行調(diào)節(jié)并不能使用任何情況，例如：石灰石飽和系數(shù)偏大但其趨勢是下降的，如若減少石灰石配比極有可能會造成石灰石飽和系數(shù)偏低。因此其專家系統(tǒng)不太完善，還應(yīng)把當(dāng)時的三率值趨勢與其值結(jié)合起來考慮。2.3趨勢識別及應(yīng)用結(jié)合現(xiàn)場實際情況以及水泥生料三率值的允許波動范圍，針對三率值趨勢情況和其與指標(biāo)的相對位臵做了水泥生料三率值趨勢情況劃分，見表1。表1水泥生料三率值趨勢情況劃分工也帖iiE工況罪別旳

12、時,J=L劃輕'嚴(yán)宦伉KH”5M-FE的狛j勺加?？傂钥?吃召卜豈旳土龍SMFIZ旳仇対KH旳Ml好狂“IL腔站均溝平雜SMkm旳tn旳叫址-、FE旳tfl刈“ft".11.吃始旳刈屮務(wù)莒HFE的航為"優(yōu)S5-M的戰(zhàn)內(nèi)”視”.H站弟旳血中甩工說SIMFE鞘訕旳”優(yōu)"-KH的的対”良"”U血坍為"”j't,P-'PSMKH的譏旳”冊、FEmtri沿“良11蛙到旳“+”迪K卜!ferrjffL/j-ft-.SM的ItL為“R-.11.也松劉“十”或、”KH-SM-rE的儻旳IL匕側(cè)平E?KH的<Jtt為刊優(yōu)Sh4h.F

13、F的ifL為"EtL，J=L©第均屮檜的蝕為帀優(yōu)",KMFE的g沁"Hr上LZ対垃屮檜FEgffl羽-優(yōu)SKM,SMfTJat為"耳"”上L©勢垃屮檜ICHSMferrj<ii旳“釀".J1L也第為f口瓷SM的口:旳帀七亡碎卜KHfcn<ifn上L=sm鈣劉”令“我FE帕恆劉n.KH甲smwjinjaj”elr.匸L出韜劉“瞥稈我“嚴(yán)SM,FE的廿L為“憂"旦趟齡十悽;KH的位対4趨診珀SM,FW的也.為"老"上L蘭給嚇豐皂；的左L*J屮“趨豹為“SM上也軸下程1；FE.K

14、H的恤為屮、洛魚給為"+-1-嚴(yán)£21R1(ft沖”良-型”憂“.上L趨豹平年總；FF.KH邯忸為+、g曲九”+曲P屮FFJ±趨將、嚴(yán)fS；32tKH的世1為+“上L超輕為"*fefrjfn-咸”良”.±1rts；sM序“伍為電”uhrmn況>粽”一h.超丹決j“*“承Z_"SMFE的位劉"彳：1£".bKta；KHfTJfltL>Jw_,iSS9-HEPE(fliE5."M.iiWKfeiSMWs<nIP*fl-.KH擴Mil.加k.ij故曲Kh-Sm-FE的債駕人*

15、5;"+J-LZ*W*_M沖況KH-SIM,FE的偵邸劉"¥"+JLQ御亠盧色SM的佰，三“龍”屮實“比“l(fā)l.A?竭型.KH-FE的伯Ayr龍-l-l.7y*KH-SIVI-FE的俏丄勺IL匕嗎平錘KH-5M-FE的脩.1勺妁"中".丨L匕妁対”+S世生料三率值的波動是圍繞指標(biāo)進行的因此根據(jù)三率值與指標(biāo)的位臵進行了等級劃分，三率值的“優(yōu)”，“良”，“中”，“差”的評價等級劃分。因為石灰石飽和系數(shù)允許波動范圍為士0.02，硅率鐵率允許波動范圍為士0.1，經(jīng)過對大量數(shù)據(jù)分析研究以及詢問現(xiàn)場操作人員，經(jīng)過仿真測試得出劃分評價標(biāo)準(zhǔn)如表2所示

16、。表2參數(shù)范圍對應(yīng)的評價等級優(yōu)良中差KH實吋勲與指林：0-0.010.02-0.03或>0.03或-0.01-0<-0.030.05”0.1臥丨(1.2>0.2咸或或川血-01“-0.05I-E實時蠡與指揃0”(1.0S0.05-0.10.1-仏2>0.2戒或或-0.05-0.2-n,|<-0.2針對趨勢的識別采用了特征知識庫識別的方法，首先判斷時間窗內(nèi)三率值和指標(biāo)的相對位，以判斷石灰石飽和系數(shù)為例一個集合規(guī)則的例子如下：當(dāng)前時間窗內(nèi)每分鐘的石灰石飽和系數(shù)值和指標(biāo)的差在一灰石飽和系數(shù)值為優(yōu)在指標(biāo)附近EndifIf當(dāng)前時間窗內(nèi)每分鐘的石灰石飽和系數(shù)值和指標(biāo)的差在輸出

17、當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)值為良，在指標(biāo)上方EndifIf當(dāng)前時間窗內(nèi)每分鐘的石灰石飽和系數(shù)值和指標(biāo)的差在一輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)值為良，在指標(biāo)下方0.01和0.01內(nèi)Then輸出當(dāng)前時間窗內(nèi)石0.01和0.02內(nèi)Then0.01和一0.02內(nèi)ThenEndif以此類推分別判斷三率值和指標(biāo)的相對位臵評價等級。然后再判斷時間窗內(nèi)的趨勢，其趨勢可分為:平穩(wěn)，上升，下降，先上升后下降，先下降后上升，不連續(xù)上升，不連續(xù)下降其中的一個規(guī)則集合如下：If時間窗內(nèi)的結(jié)束值一初始值在一0.01和0.01內(nèi)Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢平穩(wěn)EndifIf時間窗內(nèi)的結(jié)束值一初始值大于0.01內(nèi)

18、ThenIf時間窗內(nèi)的值連續(xù)遞增Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢連續(xù)上升Elseif時間窗的極大值時間窗初始值A(chǔ)nd時間窗的極大值時間窗結(jié)束值Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢先上升后下降Elseif時間窗的極大值時間窗初始值A(chǔ)nd時間窗的極大值=時間窗結(jié)束值Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢不連續(xù)上升EndifEndifIf時間窗內(nèi)的結(jié)束值一初始值小于一0.01內(nèi)ThenIf時間窗內(nèi)的值連續(xù)遞減Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢連續(xù)下降Elseif時間窗的極小值時間窗初始值A(chǔ)nd時間窗的極小值時間窗結(jié)束值。Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢先下降后上

19、升Elseif時間窗的極小值時間窗初始值A(chǔ)nd時間窗的極小值=時間窗結(jié)束值Then輸出當(dāng)前時間窗內(nèi)石灰石飽和系數(shù)趨勢不連續(xù)下降EndifEndif對硅率和鐵率判斷平穩(wěn)上升標(biāo)準(zhǔn)時經(jīng)過仿真測試得出波動范圍在0.03內(nèi)可認(rèn)為波動趨勢為穩(wěn)定，大于0.03則認(rèn)為上升或者下降。把三種率值的情況綜合起來就是表1劃分出來的趨勢模版。但是當(dāng)識別出時間窗內(nèi)三率值的趨勢和值后并不能為水泥生料配料自動配料過程給出較好的調(diào)節(jié)指導(dǎo)，因而在趨勢識別的基礎(chǔ)上進一步的對趨勢和率值進行分析計算以及結(jié)合人工經(jīng)驗給出調(diào)節(jié)初步策，判斷出是調(diào)節(jié)兩種配比還是三種或者四種配比。把各種情況綜合起來建立規(guī)則集如表3列舉了部分規(guī)則，其中0為不需要

20、對相應(yīng)的原材料繼續(xù)配比進行調(diào)節(jié)，當(dāng)為1時應(yīng)對相應(yīng)的原材料進行配比調(diào)節(jié)（完整的調(diào)節(jié)策略見附錄）。表3識別結(jié)果調(diào)節(jié)策略三*粘孫以屋其世的評倚野拔怙況比砂樸配比比怏應(yīng)潰冏L匕KH,SM,Fc值為優(yōu)趙菇不必葦世000OKU.SM00OOKH,SM,FE丨SM.KU偏周I杯00O位KH.SM,FEfflZiE.KSM-F託平dKH鞭近U標(biāo)00OfflKH,SM,FFISM.Kill平相”Fe偏離口標(biāo)000OifiKHSM,FE血為山JIKEiT3，SM萍亠打杯tn.000OKH,SM,FE血為陸且KHSM偏滯口標(biāo)肌.n0noFe靠近目標(biāo)傭KH*SX1.FEin>j曲且SMKII.Fe僞馬口標(biāo)n0o

21、o他KH.SM,FE值為住ii.SMKH，F(xiàn)e菲近日標(biāo)n0oofl1!KH,SM,FE<&?JHL1LSMKU偏禺碼標(biāo)m.n00oI<'i'lH柿、J|'|KH,SM,FE11SMV,Kri車砸切標(biāo).Fc000ofa離目標(biāo)愷KHSM,FE譏為巾耳平能1I11KH,SMFEfiJ*中FLEM.FETftKI1惱肉H*尸1111fitKli.Shi.FEWL為中.H.SM.FE平IBKI1岸近H札1111ffLtnKII，SM，1-1liJ'l1IISM，KIII"橙心桝近冃標(biāo)flK.IISM；口悄.為1冃一JCIIm.re仇為優(yōu)芯孫不

22、必：號曲：KH,SM:,'i-JJL11栩値it伯:対仇枱松十必獰垃；it伯:対仇枱松十必獰垃；SMft誠圣Jl.fiHe仙為優(yōu)血列平必曲諂；KJ2X1.帝虹曰標(biāo)iiLi,sM臨離曰標(biāo)iilf.1址伯:旳優(yōu)越始不必芳也：SM伯.旳IP疋徒SM俯離口】上值対優(yōu)趨必號也；KU,SMHi>J|I'LjKKU怖肉E枸;伉“菲近II標(biāo)垃偵為優(yōu)畫釗不，姑號龍；KIL.住伯珀1迤塔II郴輜I近冃標(biāo)fiTSM優(yōu)體第不必考宓、KHFcft為中或帶冃柳偏L為B標(biāo)伯：擁譏根捌芯約和配比時節(jié)帥糞劃分汨廠II況樓版z34水泥丨料調(diào)配応行T況創(chuàng)竹1:洵11祝2KiS31:況4需訥卩閱種扯比馬調(diào)疔二種

23、嘲比需調(diào)沾貝1種"比如表所示將工況類別分為四種分別為調(diào)節(jié)兩種，三種四種配比以及不調(diào)節(jié)配比。根據(jù)趨勢識別得出配比調(diào)節(jié)種類后，控制算法可以選取不同的模型進行控制。配比調(diào)節(jié)的主要依據(jù)公式為其中C.S.A.F介別表小-水泥!一料中傳Ca、Si.Al.I-e的禽r：+匚I一匚I、舟一町、叫一丹.|、J|一./.|介別表小:四種原材料中CaSiARFc的含量*兀一兀分別表示四種原材料的配比+以上的單位均為r|分比”求出水泥£士料成分后可以根據(jù)公式(2.1)(2.2)(2.3)計算出水泥生料ZZ率值“判斷原材料配比改變策略是結(jié)合人工經(jīng)驗以及根據(jù)此i-f算步驟逆推岀來的結(jié)呆Rd時"

24、;2.4仿真驗證選取50組歷史數(shù)據(jù)三率值三個指標(biāo)分別設(shè)為0.95,2.55,2.1進行識別仿真數(shù)據(jù)和識別結(jié)果如表4表5仿真運行結(jié)果圖KHSHFeKII位置0.9612.552AS0.9532,542J40.9552.5斗2J4優(yōu)指標(biāo)附近0.9422542J50»412.542.15KH趨勞別付鬼闕超嗚F<?位置卜k趨勢下降優(yōu)抬標(biāo)附近平莊0.9412.542150.9512,542,150.9522.542.140.9552.552.130.9532.562.13優(yōu)指標(biāo)附近平穩(wěn)優(yōu)指標(biāo)附近平葩黑標(biāo)附平穩(wěn)近0.9552.562.130.9522.572.120.9562.5623優(yōu)捋標(biāo)附近0.9632.562.130.9652552A2不連續(xù)上升憂指標(biāo)附近平穩(wěn)捫標(biāo)附平穩(wěn)近0.9632562.120r9742.572J00/)7?2.572.10中指標(biāo)上