于東偉 湯 偉 游彥卿 王孟效 張 秦

(1.陜西科技大學(xué)輕工與能源學(xué)院，陜西西安，710021;2.陜西科技大學(xué)自動(dòng)化研究所，陜西咸陽，712000;3.四川銀鴿竹漿紙業(yè)有限公司，四川瀘州，646300)

置換蒸煮系統(tǒng) (Displacement Digester Systems，簡稱DDS)是由快速置換加熱蒸煮技術(shù) (RDH)制漿技術(shù)發(fā)展而來的，可通過多個(gè)置換階段實(shí)現(xiàn)化學(xué)品與熱能的置換回用，是一項(xiàng)高效節(jié)能、綠色環(huán)保的間歇式制漿技術(shù)［1］。由于置換蒸煮過程取樣困難，缺乏穩(wěn)定、可靠、價(jià)廉、精度高的卡伯值在線傳感器和測量儀表。長期以來，蒸煮終點(diǎn)的判定是置換蒸煮系統(tǒng)一直未能很好解決的難題。與其他傳統(tǒng)的間歇式蒸煮過程不同，置換蒸煮是一個(gè)連續(xù)的生產(chǎn)過程，因而，傳統(tǒng)的離線式實(shí)驗(yàn)判定方法因存在較大滯后性而不再適用。工藝技術(shù)人員雖然制定了標(biāo)準(zhǔn)升溫曲線，但由于實(shí)際升溫曲線很難與標(biāo)準(zhǔn)升溫曲線擬合，所以由升溫曲線指導(dǎo)的人工經(jīng)驗(yàn)判別法存在很大誤差。針對這一難題，本研究在與傳統(tǒng)判別方法分析比較的基礎(chǔ)上，采用分階段計(jì)算H-因子的方法，減少了H-因子的計(jì)算誤差。提出了蒸煮有效溫度的處理機(jī)制，避免了因儀表故障引起的錯(cuò)誤動(dòng)作。對于置換蒸煮終點(diǎn)的判定，本研究結(jié)合蒸煮有效溫度處理機(jī)制，提出了分階段計(jì)算H-因子的在線計(jì)算方法。以STEP7和WinCC為軟件開發(fā)平臺(tái)，開發(fā)了在線計(jì)算的西門子PLC軟件，并完成了硬件配置和程序編寫，實(shí)現(xiàn)了置換蒸煮終點(diǎn)的有效預(yù)測。

1 H-因子的定義及在DDS中的復(fù)雜性分析

蒸煮是一個(gè)十分復(fù)雜的物理化學(xué)過程，是置換蒸煮的重要環(huán)節(jié)?？ú凳羌垵{重要的質(zhì)量指標(biāo)，是評價(jià)過程控制的一個(gè)重要參數(shù)。如果蒸煮不夠，紙漿卡伯值過大，會(huì)產(chǎn)生質(zhì)量不合格的生漿;如果蒸煮過頭，紙漿卡伯值過小、得率減少，且纖維遭到破壞、強(qiáng)度降低，會(huì)產(chǎn)生質(zhì)量不合格的爛漿［2］。研究表明，蒸煮過程中紙漿卡伯值與H-因子的對數(shù)有良好的線性關(guān)系。通過H-因子可以反應(yīng)蒸煮的程度，因此，H-因子的計(jì)算非常重要。

1.1 H-因子的定義

在堿法蒸煮過程中的影響因素很多，而且這些影響因素彼此之間又是相互作用的，如溫度和時(shí)間是兩個(gè)彼此關(guān)聯(lián)的影響因素，蒸煮溫度高了則蒸煮時(shí)間相應(yīng)需要縮短，蒸煮溫度低了則蒸煮時(shí)間相應(yīng)需要延長。因此，1957年加拿大造紙研究所K.E.Vroom提出將蒸煮溫度與蒸煮時(shí)間兩個(gè)因素結(jié)合成一個(gè)變量即H-因子，用以控制蒸煮過程，使成漿的質(zhì)量控制在要求的范圍之內(nèi)。H-因子并不是蒸煮溫度和蒸煮時(shí)間的簡單結(jié)合，而是與蒸煮溫度有關(guān)的相對反應(yīng)速率常數(shù)和蒸煮時(shí)間的結(jié)合為H-因子［3］。再具體地說，H-因子就是在不同蒸煮溫度下的相對反應(yīng)速率常數(shù)對其相應(yīng)的蒸煮時(shí)間所做出的曲線下求的面積。

1.2 H-因子的計(jì)算

H-因子是關(guān)于溫度的泛函數(shù)，其定義見公式(1)。

結(jié)合阿倫尼烏斯 (Arrhenius)公式，可得公式(2)。

式中，K為相對反應(yīng)速率常數(shù);Ea為活化能，三個(gè)脫木素階段的活化能值并不相同;則H-因子計(jì)算公式也并不相同;R為氣體常數(shù) (8.3 J/mol)，lnK為反應(yīng)常數(shù)，T為蒸煮過程中的絕對溫度。因?yàn)檎糁箦伾现邢聹囟炔⒉灰恢?，H-因子計(jì)算中的溫度為蒸煮鍋上、中、下溫度的平均值。

設(shè)100℃時(shí)的反應(yīng)速率常數(shù)K為1時(shí)，則公式(2)可變?yōu)楣?(3)。

由公式 (3)求得lnK0=43.2，進(jìn)而可計(jì)算出硫酸鹽木漿蒸煮時(shí)任意溫度下的相對反應(yīng)速率常數(shù)KRT，見公式 (4)。

由公式 (1)和公式 (4)可計(jì)算H-因子，見公式(5)。在DCS系統(tǒng)中計(jì)算H-因子需經(jīng)過離散化處理，其計(jì)算見公式 (6)。

2 有效溫度處理機(jī)制及分段計(jì)算H-因子的方法

2.1 有效溫度的定義和計(jì)算

由于紙漿的實(shí)際蒸煮溫度滿足加和性，因此實(shí)際蒸煮過程的平均溫度見公式 (7)。

式中，V為蒸煮鍋的體積;T是蒸煮鍋中各點(diǎn)的溫度。

置換蒸煮立式蒸煮鍋為圓柱型，在鍋內(nèi)有上、中、下3個(gè)循環(huán)溫度測量點(diǎn)。以中部為坐標(biāo)原點(diǎn)，建立直角坐標(biāo)系如圖1［4］所示，平均溫度計(jì)算見公式 (8)。

因?yàn)橐哉糁箦佒胁繛樽鴺?biāo)原點(diǎn)建模，所以蒸煮鍋鍋體的總高度H=2h，S為圓柱體的截面積。由公式(8)可以看出，當(dāng)蒸煮液能良好循環(huán)，溫度分布對于相同的豎坐標(biāo)z是等溫時(shí)，溫度只是豎坐標(biāo)z的函數(shù)。利用熱工學(xué)知識(shí)可以建立溫度T隨豎坐標(biāo)z變化的數(shù)學(xué)模型，該模型與藥液流速有關(guān)，而流速又是變動(dòng)的，因此，該數(shù)學(xué)模型非常復(fù)雜，為了計(jì)算方便，實(shí)際工程上可以簡化成簡單的蒸煮有效溫度。

圖1 置換蒸煮立式蒸鍋?zhàn)鴺?biāo)原理圖

在循環(huán)加熱階段，蒸煮鍋黑液經(jīng)中部篦子從蒸煮鍋出來，經(jīng)中壓蒸汽加熱后，由循環(huán)泵送至蒸煮鍋頂部和底部［5］。因此 H-因子計(jì)算基于蒸煮鍋 TI-X03(中部)抽出來的黑液、進(jìn)入蒸煮鍋TI-X01(頂部)和TI-X02(底部)的黑液平均溫度。在實(shí)際應(yīng)用中，CPL、GLV等國外公司采用蒸煮鍋中間、頂部和底部溫度的加權(quán)平均值。因而，與H-因子相關(guān)的蒸煮有效溫度THF=((TI－X01)+(TI－X02)+2×(TIX03))/4)，其目標(biāo)值設(shè)置在160～170℃之間。當(dāng)這三個(gè)溫度中的任意一個(gè)溫度出現(xiàn)問題 (如在循環(huán)加熱過程中的管道震動(dòng)、斷線、接線松動(dòng)等)，此時(shí)控制器讀取到的溫度數(shù)據(jù)將產(chǎn)生錯(cuò)誤，大于蒸煮目標(biāo)溫度的設(shè)定值，系統(tǒng)將錯(cuò)誤地切換到保溫階段，使循環(huán)蒸煮無法正常進(jìn)行。蒸煮鍋有效溫度處理 (一)如表1所示，當(dāng)任何一個(gè)溫度計(jì)故障，控制器將產(chǎn)生錯(cuò)誤動(dòng)作。

針對上述情況，溫度計(jì)算過程中應(yīng)增加有效的處理機(jī)制。因?yàn)樵谘h(huán)加熱階段蒸煮液由中部抽出經(jīng)加熱器加熱后，從頂部、底部集流管道送入蒸煮鍋。因此，正常情況下，中部循環(huán)溫度比頂部、底部循環(huán)溫度低。同時(shí)頂部和底部蒸煮溫差不能超過5℃，否則溫差會(huì)干擾黑液循環(huán)。因此，頂部、底部的蒸煮溫度與中部的溫差不能超過2.5℃。通過上面的分析，本研究對溫度的處理如下。

第一步設(shè)置故障等級。三個(gè)溫度變送器同時(shí)故障的優(yōu)先級最高 (Ⅰ)，兩個(gè)溫度變送器同時(shí)故障的優(yōu)先級其次 (Ⅱ)，一個(gè)溫度變送器故障的優(yōu)先級較低(Ⅲ)，沒有溫度變送器故障的優(yōu)先級最低 (Ⅳ)。

第二步針對不同優(yōu)先級設(shè)計(jì)有效算法，如圖2所示。Ⅰ故障的概率很低，如果正常運(yùn)行時(shí)一旦出現(xiàn)，系統(tǒng)將直接提示最高級別的錯(cuò)誤報(bào)警，停止蒸煮并需要及時(shí)修理。Ⅱ故障出現(xiàn)后，當(dāng)中部循環(huán)溫度正常時(shí)，THF=TI-X03;當(dāng)頂部循環(huán)溫度正常時(shí)，THF=TI-X01;當(dāng)?shù)撞垦h(huán)溫度正常時(shí)，THF=TI-X02。Ⅲ故障出現(xiàn)后，當(dāng)?shù)撞繙囟裙收蠒r(shí)，THF=(TI-X01+2×TI-X03)/3;當(dāng)中部溫度故障時(shí)，THF=(TI-X01+TI-X02)/2;當(dāng)頂部溫度故障時(shí)，THF=(TI-X02+2×TI-X03)/3;Ⅳ表示沒有溫度變送器的故障，按THF=(TI-X01)+(TI-X02)+2×(TI-X03)/4)處理［6］。

表1 蒸煮鍋有效溫度處理 (一)

圖2 蒸煮溫度的處理機(jī)制

經(jīng)過處理后，可以得到蒸煮鍋有效溫度處理(二)見表2，避免了因?yàn)閮x表故障引起的錯(cuò)誤操作。

表2 蒸煮鍋有效溫度處理 (二)

制漿的蒸煮過程是一個(gè)十分復(fù)雜的物理化學(xué)過程，蒸煮過程中紙漿的卡伯值與H-因子的對數(shù)有良好的線性關(guān)系，因此H-因子可以反應(yīng)蒸煮程度。而在實(shí)際生產(chǎn)中，H-因子又無法在線直接測量。本研究通過蒸煮溫度T、活化能Ea等二次變量來獲得對主導(dǎo)變量H-因子的最佳估計(jì)。雖然處理后的有效溫度不是蒸煮鍋的真實(shí)有效溫度，但可以通過有效的溫度處理機(jī)制得到蒸煮溫度的最優(yōu)測量，避免了因控制器讀取到錯(cuò)誤溫度數(shù)據(jù)而產(chǎn)生誤操作，減少了H-因子的計(jì)算誤差。雖然計(jì)算值與實(shí)際值仍存在偏差，但也足以保證在實(shí)際生產(chǎn)允許的誤差范圍內(nèi)。

2.2 分段計(jì)算H-因子

理論上脫木素過程分為三個(gè)階段進(jìn)行，各個(gè)階段的活化能不同，這直接涉及到H-因子的計(jì)算值問題，而一般數(shù)學(xué)模型都沒有過多考慮這個(gè)問題，只是將活化能設(shè)置為大量脫木素階段的固定值134 kJ/mol，模型以大量脫木素階段為起點(diǎn)，并未考慮殘余脫木素階段與大量脫木素階段的不同，因此無法真實(shí)跟蹤實(shí)際的變化過程。本研究分階段計(jì)算H-因子，綜合考慮脫木素三個(gè)階段的活化能對于H-因子計(jì)算的影響。

由于各個(gè)脫木素階段的轉(zhuǎn)折點(diǎn)是一個(gè)模糊的區(qū)域，因此如何判斷蒸煮過程各個(gè)脫木素階段的起點(diǎn)成為急需解決的問題。文獻(xiàn) ［7］中提出了一種簡單實(shí)用的起點(diǎn)判斷法:首先在實(shí)驗(yàn)室通過多次蒸煮實(shí)驗(yàn)分析，以有效堿濃度對木素的得率作圖，確定該原料各個(gè)脫木素階段的起點(diǎn)溫度范圍，然后在該溫度范圍內(nèi)選定一個(gè)溫度作為該脫木素階段的起點(diǎn)溫度。根據(jù)不同階段的起點(diǎn)溫度確定各自的活化能。查閱相關(guān)資料［8］，表明不同纖維原料、不同制漿方法的活化能不同，以蒸煮針葉木硫酸鹽漿為例，初始脫木素階段的活化能很低約為40 kJ/mol，大量脫木素階段活化能約為134 kJ/mol，殘余脫木素階段活化能約為大量脫木素段的2/3，即90 kJ/mol。分別計(jì)算各階段的H-因子，再求其累加和，為最終的H-因子。

3 在線計(jì)算方法實(shí)現(xiàn)

3.1 硬件配置

西門子S7-400系列PLC作為現(xiàn)場控制站點(diǎn)主控制器，以工業(yè)以太網(wǎng)和PROFIBUS-DP為通信基礎(chǔ)，以DELL工控機(jī)為操作員站和工程師站。主控制柜和遠(yuǎn)程I/O柜分散布置，由站點(diǎn)模塊ET-200M及相應(yīng)的I/O模塊組成從站，可置于任何生產(chǎn)現(xiàn)場，現(xiàn)場的儀表、閥門或低壓開關(guān)柜的測控點(diǎn)都可就近與其聯(lián)接。

(1)控制站點(diǎn)配置

考慮到整個(gè)置換蒸煮控制系統(tǒng)的實(shí)時(shí)性和精確性，這里選用處理速度較高、內(nèi)存資源較大的西門子S7 412-3為控制器主體，滿足了處理大量數(shù)據(jù)和復(fù)雜算法的要求。

·數(shù)字量輸入輸出模塊 (DIO):選用SM323，1塊 (16點(diǎn)DI、16點(diǎn)DO/塊)?？商幚?6點(diǎn)數(shù)字量輸入信號和16點(diǎn)數(shù)字量輸出信號。

·模擬量輸入模塊 (AI):選用SM331，1塊 (8點(diǎn)/塊)。處理8路模擬量輸入信號。

·模擬量輸出模塊 (AO):選用SM332，1塊(4點(diǎn)/塊)。處理4路模擬量輸出信號。

(2)儀表部分配置

蒸煮溫度是H-因子計(jì)算的關(guān)鍵參數(shù)，因此溫度的測量極其重要。溫度計(jì)選擇鎧裝熱電阻，它的外保護(hù)管套采用不銹鋼，內(nèi)充高密度氧化物絕緣體，具有很強(qiáng)的抗污染性能和優(yōu)良的機(jī)械強(qiáng)度。而且直徑小、易彎曲、抗震性好、熱響應(yīng)時(shí)間快，非常適合蒸煮階段鍋溫的測量。

(3)操作員站和工程師站

操作員站具體配置為:Intel雙核CPU/1G內(nèi)存/160G硬盤/19’LCD/CP5611通訊卡/HP1020打印機(jī)/上海山特2000VA帶穩(wěn)壓功能UPS;便攜式工程師站為DELL商用筆記本電腦，平時(shí)可辦公用，只有在系統(tǒng)維護(hù)時(shí)才接入現(xiàn)場的工業(yè)以太網(wǎng)。各個(gè)操作站、工程師站可通過以太網(wǎng)共享此打印機(jī)，方便地記錄生產(chǎn)狀況及系統(tǒng)運(yùn)行情況。

3.2 軟件設(shè)計(jì)

本控制系統(tǒng)軟件分為上位機(jī)和下位機(jī)兩部分。下位機(jī)采用SIEMENS STEP7 V5.4軟件平臺(tái)進(jìn)行軟件編程;用于監(jiān)控的上位機(jī)軟件采用 SIEMENS WinCC V6.0 SP3平臺(tái)。

(1)溫度的設(shè)計(jì)

本研究中蒸煮溫度的計(jì)算是通過比較選擇不同的公式計(jì)算的，當(dāng)中部循環(huán)溫度正常時(shí)，THF=TIX03;當(dāng)頂部循環(huán)溫度正常時(shí)，THF=TI-X01;當(dāng)?shù)撞垦h(huán)溫度正常時(shí)，THF=TI-X02。當(dāng)?shù)撞繙囟裙收蠒r(shí)，THF=(TI-X01+2×TI-X03)/3;當(dāng)中部溫度故障時(shí)，THF=(TI-X01+TI-X02)/2;當(dāng)頂部溫度故障時(shí)，THF=(TI-X02+2×TI-X03)/3;沒有溫度變送器的故障，按THF=(TI-X01)+(TI-X02)+2×(TIX03))/4處理。溫度的計(jì)算通過下位機(jī)STEP 7軟件來實(shí)現(xiàn)，利用語句表的方法來編寫溫度計(jì)算程序。在編寫的溫度計(jì)算程序中，包含了多個(gè)變量，表3為各個(gè)變量的地址及含義。

表3 溫度計(jì)算程序相關(guān)變量

因在下位機(jī)軟件STEP 7中編寫的溫度計(jì)算程序相對簡單，文章就不再詳細(xì)贅述。

(2)H-因子的設(shè)計(jì)

軟件中H-因子的設(shè)計(jì)計(jì)算也分階段進(jìn)行，即分別計(jì)算初始脫木素階段、大量脫木素階段和殘余脫木素階段的H-因子，再求其累加和，為最終的H-因子。H-因子的計(jì)算通過下位機(jī)STEP7軟件來實(shí)現(xiàn)，利用語句表 (STL)的方法來編寫H-因子計(jì)算程序。在編寫的H-因子計(jì)算程序中，包含了多個(gè)變量，表4為各個(gè)變量的地址及含義［9］。

表4 H-因子程序相關(guān)變量

在分步計(jì)算H-因子時(shí)，活化能Ea的賦值十分關(guān)鍵，直接決定H-因子的計(jì)算是否精確。因此，STEP7中活化能Ea的賦值分為三種情況。當(dāng)采樣溫度小于T1(大量脫木素階段的起始溫度)時(shí)，將初始脫木素階段的Ea值放入MD44中，即給此階段的活化能賦值。同時(shí)將 H-因子的初始值置零;當(dāng)采樣溫度大于T1，小于T2(殘余脫木素階段的起始溫度)時(shí)，將大量脫木素階段的Ea值放入MD44中;同理，采樣溫度大于T2時(shí)，將殘余脫木素階段活化能賦值給Ea。

然后根據(jù)H-因子累積計(jì)算公式:Ht+Δt= Ht+ 0.5 ［elnK0－Ea/RT+elnK0－Ea/RTt+Δt］Δt，在 STEP 7 中進(jìn)行H-因子計(jì)算程序的設(shè)計(jì)。其中，T為上一時(shí)刻的采樣溫度;Δt為采樣時(shí)間;T+Δt為此次采樣溫度;Ht為上一時(shí)刻的 H-因子計(jì)算值;Ht+Δt為此次H-因子計(jì)算值。

3.3 運(yùn)行結(jié)果

此系統(tǒng)目前已經(jīng)在四川銀鴿竹漿紙業(yè)的置換蒸煮車間成功運(yùn)行，圖3和圖4給出了此系統(tǒng)在實(shí)際運(yùn)行中H-因子和蒸煮鍋溫度的控制效果畫面。對蒸煮有效溫度和H-因子進(jìn)行了監(jiān)控，有效地降低了蒸煮粗漿卡伯值的波動(dòng)范圍，為后續(xù)洗滌、篩選、漂白工序以及成紙的質(zhì)量提供了良好保證，同時(shí)也取得了明顯的節(jié)能增效效果。

4 結(jié)語

在置換蒸煮工藝機(jī)理的基礎(chǔ)上，對置換蒸煮H-因子計(jì)算方法進(jìn)行了研究。對H-因子采用分段計(jì)算，減少了H-因子的計(jì)算誤差。同時(shí)提出了蒸煮有效溫度的處理機(jī)制，避免了因?yàn)閮x表故障引起的錯(cuò)誤動(dòng)作。設(shè)計(jì)選用西門子 S7-400PLC，并通過 STEP7、WinCC軟件完成程序編寫和畫面繪制，對置換蒸煮終點(diǎn)進(jìn)行了有效控制。

圖3 實(shí)際運(yùn)行中DDS蒸煮鍋溫度監(jiān)控畫面

圖4 H-因子的累計(jì)值

［1］DING Shi-huo，ZHANG Ming-feng，WANG Wu-xiong.The Latest Development of Bach Cooking Technique of RDH in Displacement Digester System［J］.China Pulp ＆ Paper，2005，24(6):62.

丁仕火，張銘峰，王武雄.DDSTM置換蒸煮系統(tǒng)RDH間歇蒸煮技術(shù)新進(jìn)展［J］.中國造紙，2005，24(6):62.

［2］LV Ding-yun，HAN Xiao-juan，TANG Wei.The Latest Development of Displacement Batch Cooking［J］.China Pulp ＆ Paper，2007，26(12):63.呂定云，韓小娟，湯 偉.置換間歇式蒸煮的新進(jìn)展［J］.中國造紙，2007，26(12):63.

［3］Guo Chang-yu，Sun Jia-lin.Calculation and Controlling of H-factor in Continuous Digestion［J］.Paper and Paper Making，2005，24(6):68.

郭長裕，孫佳林.連續(xù)蒸煮H-因子的計(jì)算和控制［J］.紙和造紙，2005，24(6):68.

［4］LI Xiang-yang，ZHU Xue-feng，HUANG Dao-ping.Comprehensive H-Factor Based on Mean Kappa Number in Batch Pulping Process［J］.China Pulp ＆ Paper，2001，20(3):24.

李向陽，朱學(xué)峰，黃道平.間歇蒸煮過程中基于平均卡伯值的綜合H-因子的研究［J］.中國造紙，2001，20(3):24.

［5］Huang Jun-mei，Tang Wei，Xu Bao-hua.Development and Application of Displacement Digester System［J］.Control and Instruments in Chemical Industry，2010，37(8):1.

黃俊梅，湯 偉，許保華.置換蒸煮系統(tǒng)(DDS)的發(fā)展及應(yīng)用［J］.化工自動(dòng)化及儀表，2010，37(8):1.

［6］Zhang Qin.The Design and Application of Dissolving Pulp DDS Control System Based on PCS7［D］.Xi'an:Shaanxi University of Science and Technology，2013.

張 秦.基于PCS7的溶解漿DDS控制系統(tǒng)［D］.西安:陜西科技大學(xué)，2013.

［7］Li Yan，Zhang Jian，Zhu Xue-feng.Application of DCS and Soft Sensor Technique during Batch Cooking Process［J］.Control and Instruments in Chemical Industry，2003，30(4):26.

李 艷，張 健，朱學(xué)峰.DCS控制系統(tǒng)及軟測量技術(shù)在間歇紙漿蒸煮過程的應(yīng)用［J］.化工自動(dòng)化及儀表，2003，30(4):26.

［8］Li Jun-hu.Shallowly Discussed the H-Factor and the H-Factor's Application in the Alkaline Cooking［J］.Jiangsu Paper Making，2008，24(1):18.

李軍虎.淺談H-因子及H-因子在堿法蒸煮中的應(yīng)用［J］.江蘇造紙，2008，24(1):18.

［9］Huang Jun-mei.Study on Pulp Cooking Endpoint Soft Measurement in Displacement Digester System［D］.Xi'an:Shaanxi University of Science and Technology，2011.

黃俊梅.置換蒸煮(DDS)過程蒸煮終點(diǎn)軟測量的研究［D］.西安:陜西科技大學(xué)，2011. CPP