湯 偉 于東偉 張怡真 黨世紅 王 鋒

(1.陜西科技大學(xué)輕工與能源學(xué)院,陜西西安,710021；2.浙江力諾流體控制科技股份有限公司,浙江瑞安,325200)

?

基于遞推PLS的置換蒸煮終點(diǎn)軟測(cè)量建模研究

湯 偉1于東偉1張怡真1黨世紅1王 鋒2

(1.陜西科技大學(xué)輕工與能源學(xué)院,陜西西安,710021；2.浙江力諾流體控制科技股份有限公司,浙江瑞安,325200)

基于離線數(shù)據(jù)得到的卡伯值軟測(cè)量模型無法隨蒸煮過程變化而更新,針對(duì)這一問題,應(yīng)用基于移動(dòng)窗口遞推的偏最小二乘法(PLS)建立了新的卡伯值軟測(cè)量模型。通過對(duì)制漿過程和蒸煮終點(diǎn)影響因素的分析,確定了建模所需的輔助變量；利用特定大小的移動(dòng)窗口中的最新數(shù)據(jù)來更新該測(cè)量模型,能夠確保模型對(duì)蒸煮過程的跟蹤能力。仿真結(jié)果表明,與傳統(tǒng)的PLS模型相比,該模型可更好地預(yù)測(cè)置換蒸煮的蒸煮終點(diǎn)。

置換蒸煮;軟測(cè)量建模;蒸煮終點(diǎn);遞推PLS

置換蒸煮系統(tǒng)(Displacement Digester Systems,DDS)是以間歇式反應(yīng)器為基礎(chǔ)的紙漿蒸煮系統(tǒng),可通過多個(gè)置換工段實(shí)現(xiàn)化學(xué)品與熱能的回收利用,是一項(xiàng)高效節(jié)能、綠色環(huán)保的深度脫木素制漿技術(shù)。因置換蒸煮技術(shù)節(jié)能減排效果顯著,其已被廣泛用于制漿領(lǐng)域。以溶解漿生產(chǎn)為例,置換蒸煮過程主要包括裝鍋、通氣和水解、中和、加熱蒸煮、置換、卸料等過程[1]。其中,蒸煮時(shí)間直接影響成漿質(zhì)量、得率以及能耗,是蒸煮過程中最關(guān)鍵的控制因素。

置換蒸煮過程中紙漿取樣困難,且缺乏穩(wěn)定、可靠、價(jià)廉、精度高的卡伯值在線傳感器和測(cè)量?jī)x表。因此,蒸煮終點(diǎn)判定是置換蒸煮系統(tǒng)一直未能解決的問題。雖然有些采用離線數(shù)據(jù)得到的卡伯值軟測(cè)量模型解決了卡伯值難以直接測(cè)量的難題,但基于批量數(shù)據(jù)的模型一旦建立便不能隨蒸煮過程的變化實(shí)時(shí)進(jìn)行模型更新。針對(duì)這一問題,本課題應(yīng)用遞推PLS(偏最小二乘法)建立了蒸煮終點(diǎn)軟測(cè)量模型。通過分析工藝流程和蒸煮終點(diǎn)影響因素,確定建模所需的輔助變量；利用特定大小的移動(dòng)窗口中的最新數(shù)據(jù)更新模型,從而確保模型對(duì)蒸煮過程的跟蹤能力。仿真結(jié)果表明,與傳統(tǒng)PLS模型相比,該模型能更好地預(yù)測(cè)置換蒸煮終點(diǎn)。

1 卡伯值軟測(cè)量模型建模變量的選擇

卡伯值和得率是蒸煮過程的2個(gè)重要指標(biāo),是評(píng)價(jià)蒸煮工段的重要參數(shù)。如果蒸煮程度不足,雖然紙漿得率較高,但卡伯值過大,紙漿可能成為“生漿”；如果蒸煮過頭,卡伯值過小,得率降低,且纖維遭到破壞,紙漿強(qiáng)度降低,影響使用。因此,實(shí)際生產(chǎn)過程中會(huì)根據(jù)原料和蒸煮工藝不同,在綜合考慮成漿質(zhì)量和得率的情況下,確定合適的紙漿卡伯值設(shè)定值。由于卡伯值的影響因素十分復(fù)雜,蒸煮條件又時(shí)時(shí)變化,很難建立準(zhǔn)確的機(jī)理模型，但是,蒸煮過程中有大量的過程數(shù)據(jù)和快捷準(zhǔn)確的檢測(cè)方法,適合采用統(tǒng)計(jì)方法建模。因此,本實(shí)驗(yàn)采用遞推PLS建立卡伯值軟測(cè)量模型。

1.1 卡伯值影響因素分析

為了準(zhǔn)確建立卡伯值的預(yù)測(cè)模型,首先需對(duì)卡伯值影響因素進(jìn)行分析,以避免建模中過程變量信息不全或存在冗余,降低模型精度。生產(chǎn)實(shí)踐表明,影響蒸煮終點(diǎn)卡伯值的主要因素是蒸煮溫度、蒸煮時(shí)間、有效堿濃度、H因子等,這些因素錯(cuò)綜復(fù)雜地對(duì)卡伯值產(chǎn)生影響[2]。在蒸煮過程中,蒸煮溫度和時(shí)間是2個(gè)彼此關(guān)聯(lián)的影響因素,一般由H因子來表示它們的綜合影響,其值為不同蒸煮溫度下的反應(yīng)速度對(duì)蒸煮時(shí)間所得曲線下包圍的面積。雖然H因子中已包含溫度和時(shí)間,但如果去掉其中任何一個(gè)變量,可能導(dǎo)致模型包含的有用信息不足、魯棒性減弱。故本研究將其全部作為二次變量[3]。加熱循環(huán)工段的流量與蒸煮溫度相互耦合,會(huì)導(dǎo)致升溫、保溫工段蒸煮鍋內(nèi)部溫差增大,蒸煮不均勻。另外，各個(gè)工段白液用量也會(huì)對(duì)蒸煮終點(diǎn)卡伯值產(chǎn)生一定影響。此外,當(dāng)蒸煮用堿量一定時(shí),合理控制液比可以保證適當(dāng)?shù)恼糁笏幰簼舛?有利于藥液在原料中的浸透,有利于蒸煮鍋內(nèi)藥液循環(huán)。與傳統(tǒng)蒸球蒸煮和立鍋蒸煮不同,置換蒸煮系統(tǒng)在水解、中和、熱充和置換4個(gè)過程中都會(huì)進(jìn)行蒸煮液的填充和回流。在整個(gè)蒸煮過程中,蒸煮藥液體積一直在變化且無法測(cè)量,無法計(jì)算液比,故忽略液比[4]。

因此,將H因子、有效堿濃度等13個(gè)可測(cè)量變量作為輸入變量(見表1)用于建立卡伯值軟測(cè)量模型。

表1 卡伯值建模變量表

1.2 潛變量數(shù)量確定

在PLS回歸中,確定潛變量數(shù)量非常關(guān)鍵。潛變量數(shù)量不足會(huì)使模型包含的有用信息不足,不僅模型對(duì)建模樣本的擬合精度不夠,模型也不能擁有對(duì)測(cè)試樣本良好的預(yù)測(cè)能力。相反,潛變量數(shù)量過多會(huì)使系統(tǒng)噪聲混入模型,造成模型對(duì)建模樣本的過擬合,而對(duì)測(cè)試樣本的預(yù)測(cè)精度不高。因此,適宜的潛變量數(shù)量對(duì)于提高模型的預(yù)測(cè)精度至關(guān)重要。在多數(shù)情況下,PLS回歸并不需要將全部的潛變量應(yīng)用于回歸建模中[5]。

(j=1,2,…,p)

(1)

Y=(y1,…,yp)T的預(yù)測(cè)誤差平方和為：

(2)

將樣本集分為若干組,剔除1組,利用余下的樣本來建立模型,再用剔除的那組樣本作為檢驗(yàn)樣本,計(jì)算模型在檢驗(yàn)樣本上的預(yù)測(cè)誤差,然后重復(fù)上述步驟,直至將每組數(shù)據(jù)都剔除過1次。將每組數(shù)據(jù)的模型預(yù)測(cè)誤差求和得到預(yù)測(cè)殘差累積平方和。分別計(jì)算取不同數(shù)量的特征向量時(shí)所對(duì)應(yīng)的殘差累計(jì)平方和；當(dāng)殘差累計(jì)平方和取最小值時(shí),潛變量數(shù)量作為模型最后保留的潛變量數(shù)量。建模時(shí)對(duì)應(yīng)潛變量數(shù)量的訓(xùn)練集誤差曲線如圖1所示。由圖1可確定潛變量數(shù)量為6個(gè)。

圖1 潛變量數(shù)量訓(xùn)練集誤差曲線

2 基于移動(dòng)窗口的遞推PLS算法

2.1 基本原理

作為一種基于數(shù)據(jù)回歸模型的軟測(cè)量建模方法,PLS在蒸煮終點(diǎn)卡伯值預(yù)測(cè)的建模與控制中得到廣泛重視。但隨著時(shí)間的推移,蒸煮原料的特性、工藝流程和蒸煮工段都可能發(fā)生變化,因此，必須對(duì)模型進(jìn)行在線更新,才能具有較強(qiáng)的工況變化跟蹤能力。傳統(tǒng)的PLS方法是基于批量數(shù)據(jù)的模型,即使用整批離線數(shù)據(jù)來確定模型參數(shù)。當(dāng)新數(shù)據(jù)增加到一定程度時(shí),把新、舊數(shù)據(jù)結(jié)合起來重新建立模型。這樣,基于批量數(shù)據(jù)的模型一旦建立便不能及時(shí)適應(yīng)過程變化。此外,舊數(shù)據(jù)被重復(fù)利用,而且隨著時(shí)間的推移,需要存儲(chǔ)的舊數(shù)據(jù)信息越來越多,在線訓(xùn)練模型時(shí)所涉及的數(shù)據(jù)量也越來越大。

本研究將移動(dòng)窗口的概念和PLS有效結(jié)合在一起,利用移動(dòng)窗口內(nèi)的數(shù)據(jù)進(jìn)行卡伯值建模。窗口按固定步長(zhǎng)向前移動(dòng),每向前移動(dòng)1次,卡伯值測(cè)量模型將在線更新1次,解決了固定模型對(duì)過程變化跟蹤能力差的問題。其中,移動(dòng)窗口長(zhǎng)度不能太小,否則卡伯值預(yù)測(cè)模型會(huì)有較大誤差,太大又不利于在線計(jì)算；步長(zhǎng)越小,新數(shù)據(jù)起的作用越大,卡伯值測(cè)量模型更新的頻率越高[6-10]。

本實(shí)驗(yàn)以四川銀鴿竹漿有限公司置換蒸煮車間2014年11月的84組蒸煮數(shù)據(jù)為實(shí)驗(yàn)數(shù)據(jù)樣本,將這些數(shù)據(jù)分成兩部分,第一部分為前40組數(shù)據(jù),作為訓(xùn)練樣本用來建立卡伯值預(yù)測(cè)模型；第二部分為后44組數(shù)據(jù),作為測(cè)試樣本,檢驗(yàn)?zāi)Ｐ蛯?duì)卡伯值的預(yù)測(cè)精度。移動(dòng)窗口長(zhǎng)度為24,步長(zhǎng)為1。

2.2 理論推導(dǎo)

給定2個(gè)經(jīng)標(biāo)準(zhǔn)化的矩陣x∈Rn·m和y∈Rn·p,假設(shè)它們線性相關(guān)：

y=xcpls+v

(3)

式中,cpls為回歸系數(shù)矩陣,v為殘差矩陣。

用傳統(tǒng)的PLS回歸建模時(shí),要對(duì)x和y進(jìn)行如下分解：

(4)

(5)

式中,t1和u1為第1對(duì)輸入和輸出的分向量,p1和q1為相應(yīng)的負(fù)荷向量,E1和F1為殘差矩陣。

在構(gòu)造了式(4)和式(5)的外部模型后,特征向量間通過一元線性回歸模型建立聯(lián)系：

u1=b1t1+r1

(6)

式中,r1為殘差向量,b1為回歸因子：

(7)

式(7)確定了PLS的內(nèi)部模型。此時(shí),再計(jì)算回歸系數(shù)矩陣：

y=uqT=tbqT=xpbqT

(8)

cpls=pbqT

(9)

3 MATLAB仿真及結(jié)果分析

將表1所列的13個(gè)過程變量作為建模輸入,建立了基于遞推PLS算法的置換蒸煮卡伯值預(yù)測(cè)模型。在實(shí)際生產(chǎn)過程中還有加熱循環(huán)工段的流速、各組分濃度等其他變量會(huì)對(duì)卡伯值產(chǎn)生影響,但因其影響較小和模型運(yùn)算速度問題,本課題建立的卡伯值預(yù)測(cè)模型將其忽略。

圖2為PLS預(yù)測(cè)模型與遞推PLS預(yù)測(cè)模型的預(yù)測(cè)效果比較，圖3為這2種模型的預(yù)測(cè)誤差對(duì)比。從圖2和圖3可以看出,在曲線前段,遞推PLS模型和傳統(tǒng)PLS模型的預(yù)測(cè)值很接近,但隨時(shí)間推移,過程發(fā)生變化,PLS模型預(yù)測(cè)值出現(xiàn)漂移,而遞推PLS模型顯示出更強(qiáng)的跟蹤過程變化的能力,具有更高的預(yù)測(cè)精度；同時(shí),遞推PLS模型更穩(wěn)定,不會(huì)出現(xiàn)如PLS模型的預(yù)測(cè)值偏離較大的問題。但與卡伯值實(shí)測(cè)值相比,這2個(gè)模型都存在一定的預(yù)測(cè)誤差。因此,為進(jìn)一步改進(jìn)軟測(cè)量效果,可以考慮與人工神經(jīng)網(wǎng)絡(luò)相結(jié)合構(gòu)建非線性遞推PLS的自適應(yīng)模型。

為了衡量模型的預(yù)測(cè)精度,本課題采用模型預(yù)測(cè)的平均預(yù)測(cè)誤差(MPE)、均方誤差(MSPE)、均方根誤差(RMSE)3個(gè)指標(biāo)來衡量。2種模型的誤差對(duì)比見表2。

(10)

(11)

(12)

紙漿實(shí)測(cè)卡伯值與2種模型預(yù)測(cè)的卡伯值對(duì)比見表3。由表3可知,與傳統(tǒng)PLS模型相比,基于遞推PLS建立的卡伯值預(yù)測(cè)模型所得預(yù)測(cè)卡伯值更接近卡伯值實(shí)測(cè)值,其誤差更小。雖然該方法測(cè)量結(jié)果仍有一定誤差(小于±1),但該誤差結(jié)果已處于實(shí)際生產(chǎn)允許的誤差范圍內(nèi)。

圖2 2種模型預(yù)測(cè)效果比較

表2 2種模型的誤差對(duì)比

表3 卡伯值實(shí)測(cè)值與模型預(yù)測(cè)值的對(duì)比

4 結(jié) 語

針對(duì)置換蒸煮終點(diǎn)難以判定這一難題,以四川銀鴿竹漿紙業(yè)有限公司置換蒸煮車間的蒸煮數(shù)據(jù)為數(shù)據(jù)樣本,采用遞推PLS(偏最小二乘法)建立了蒸煮終點(diǎn)軟測(cè)量模型。利用特定大小的移動(dòng)窗口中的最新數(shù)據(jù)更新模型,從而確保模型對(duì)蒸煮過程的跟蹤能力；該廠實(shí)際蒸煮數(shù)據(jù)的驗(yàn)證結(jié)果后證明了該建模方法的有效性。雖然該模型預(yù)測(cè)的卡伯值與卡伯值實(shí)測(cè)值仍存在±1的誤差,但是該誤差結(jié)果已處于實(shí)際生產(chǎn)允許的誤差范圍內(nèi)。仿真結(jié)果表明,與傳統(tǒng)的基于PLS的預(yù)測(cè)模型相比,本實(shí)驗(yàn)建立的基于遞推PLS的預(yù)測(cè)模型能夠更好地預(yù)測(cè)置換蒸煮過程的蒸煮終點(diǎn)。

[1] 張 秦. 基于PCS7的溶解漿DDS控制系統(tǒng)[D]. 西安： 陜西科技大學(xué), 2013.

[2] 呂定云, 韓小娟, 湯 偉. 置換間歇式蒸煮的新進(jìn)展[J]. 中國(guó)造紙, 2007, 26(12)：68.

[3] 于東偉, 湯 偉, 游彥卿. 置換蒸煮過程H因子在線計(jì)算方法研究與實(shí)現(xiàn)[J]. 中國(guó)造紙, 2014, 33(4)： 56.

[4] 湯 偉, 王 震, 甘文濤. 蒸煮鍋內(nèi)溫差DMC-PID串級(jí)解耦控制[J]. 中國(guó)造紙, 2014, 33(12)： 47.

[5] 袁 勇. 基于多向偏最小二乘的間歇過程軟測(cè)量建模研究[D]. 沈陽： 東北大學(xué), 2008.

[6] 汪小勇, 梁 均, 劉育明, 等. 基于遞推PLS的自適應(yīng)軟測(cè)量模型及其應(yīng)用 [J]. 浙江大學(xué)學(xué)報(bào), 2005, 39(5)： 676.

[7] Hector J Galicia, He Q Peter, Wang Jin. A Reduced Order Soft Sensor Approach and Its Application to a Continuous Digester[J]. Journal of Process Control, 2011, 21(4): 489.

[8] 李春夫, 葉 昊, 王桂增. 基于多向PLS方法的間歇過程質(zhì)量預(yù)測(cè) [J]. 系統(tǒng)仿真學(xué)報(bào), 2004, 39(6)： 1168.

[9] Shen Wenhao, Chen Xiaoquan. Measuring and Controlling Model of Pulp Kappa Number with Spectroscopy during Batch Sulfite Pulping Process [J] . Ind.Eng.Chem.Res., 2009, 48(19)： 8980.

[10] 李向陽, 朱學(xué)峰, 劉煥彬. 間歇制漿蒸煮過程的混合建模方法研究[J]. 中國(guó)造紙學(xué)報(bào), 2001, 16(2)： 24.

(責(zé)任編輯：王 巖)

Research on Cooking End Point Modeling Based on Recursive PLS for Displacement Cooking

TANG wei1,*YU Dong-wei1ZHANG Yi-zhen1DANG Shi-hong1WANG Feng2

(1.CollegeofLightIndustryandEnergy,ShaanxiUniversityofScienceandTechnology,Xi’an,ShaanxiProvince, 710021；2.ZhejiangLinuoFluidControlTechnologyCo.,Ltd,Ruian,ZhejiangProvince, 325200)

(*E-mail： wtang906@163.com)

Kappa number model based on offline data can not be updated with the pulping process changes. By combining rolling window with recursive partial least squares (PLS), a new adaptive Kappa number model was proposed in the paper. In order to determine the supplementary variables for modeling instrumental, the pulping process and the factors influencing the cooking end point were analyzed. The model could be updated by using the latest data in a rolling window with specific dimension, if could ensure the ability of the model to follow the track of cooking process. Simulation results with real industrial data showed that the adaptive Kappa number model had higher precision and better prediction of cooking end point of displacement cooking than the traditional PLS model.

displacement digester； optical fiber on-line measurement； cooking end point； multivariable binomial regression

2015- 03-13

國(guó)家國(guó)際科技合作項(xiàng)目(2010DFB43660)；陜西省重點(diǎn)科技創(chuàng)新團(tuán)隊(duì)計(jì)劃項(xiàng)目(2014KCT-15)。

湯 偉,男,1971年生；教授；主要研究方向：工業(yè)智能控制及工業(yè)高級(jí)過程控制。 E-mail:wtang906@163.com

TP273

A

1000- 6842(2015)03- 0047- 04