朱以斌

(中國(guó)石化上海石油化工股份有限公司芳烴部，上海200540)

中國(guó)石化上海石油化工股份有限公司(以下簡(jiǎn)稱上海石化)1#芳烴抽提裝置是根據(jù)德國(guó)克虜伯公司提供的N-甲酰嗎啉(NFM)抽提蒸餾專利技術(shù)進(jìn)行設(shè)計(jì)的國(guó)產(chǎn)化芳烴抽提裝置。裝置處理加氫汽油420 kt/a，生產(chǎn)純苯191.1 kt/a、甲苯98 kt/a；副產(chǎn)抽余油58.8 kt/a、C8+餾分72.8 kt/a。乙烯副產(chǎn)物裂解汽油經(jīng)過(guò)加氫后主要含有芳烴和非芳烴兩類物質(zhì)，兩者的蒸汽壓比較接近，用一般的蒸餾方法難以將它們分離，通過(guò)在原料中加入選擇性溶劑NFM，被NFM溶解的芳烴不易揮發(fā)，進(jìn)而用蒸餾的方法將芳烴與非芳烴分離。

影響抽余油中苯質(zhì)量分?jǐn)?shù)的因素有溶劑比、溶劑進(jìn)料溫度、塔壓、塔底溫度、塔真空度等。抽余油中苯的質(zhì)量分?jǐn)?shù)是影響芳烴收率的主要因素，由于影響抽余油中苯質(zhì)量分?jǐn)?shù)的因素多，變量與變量之間有很強(qiáng)的相關(guān)性，生產(chǎn)機(jī)理比較復(fù)雜，很難從復(fù)雜數(shù)據(jù)中提取有價(jià)值的信息，制約了收率的提高。裝置的工藝機(jī)理復(fù)雜，難以用機(jī)理模型解釋，而進(jìn)行數(shù)據(jù)分析時(shí)，受變量多、且非線性，變量之間存在較強(qiáng)耦合性等因素制約，很難用一般的數(shù)理統(tǒng)計(jì)方法處理。對(duì)于小樣本集，運(yùn)用神經(jīng)網(wǎng)絡(luò)等較新興的機(jī)器學(xué)習(xí)方法訓(xùn)練出來(lái)的模型結(jié)構(gòu)，存在過(guò)擬合與欠擬合、局部極小點(diǎn)問(wèn)題。

由思華數(shù)據(jù)有限公司研發(fā)的數(shù)據(jù)挖掘技術(shù)(DMOS)優(yōu)化軟件為工程師提供了一個(gè)有效分析生產(chǎn)數(shù)據(jù)的工具。該軟件能處理多變量、非線性、變量之間強(qiáng)關(guān)聯(lián)、數(shù)據(jù)又被“噪聲”污染的復(fù)雜數(shù)據(jù)。DMOS采用數(shù)據(jù)挖掘技術(shù)，綜合應(yīng)用數(shù)理統(tǒng)計(jì)、模式識(shí)別、機(jī)器學(xué)習(xí)、計(jì)算機(jī)可視化等技術(shù)，從大量的數(shù)據(jù)集合中發(fā)現(xiàn)隱藏的、有價(jià)值的規(guī)律進(jìn)行優(yōu)化生產(chǎn)。采用DMOS優(yōu)化軟件可以進(jìn)行裝置操作參數(shù)的優(yōu)化，提供優(yōu)化操作方案，達(dá)到降低抽余油中苯質(zhì)量分?jǐn)?shù)的目的。

20世紀(jì)90年代中期，在統(tǒng)計(jì)學(xué)習(xí)理論基礎(chǔ)上發(fā)展起來(lái)的支持向量機(jī)算法在解決小樣本、非線性及高維模式識(shí)別問(wèn)題的過(guò)程中所取得了核函數(shù)應(yīng)用等方面的突出進(jìn)展，參照這一優(yōu)化算法，DMOS優(yōu)化軟件不斷改進(jìn)，以期在處理復(fù)雜數(shù)據(jù)方面明顯的優(yōu)勢(shì)成為工程化的優(yōu)化軟件。

1 DMOS優(yōu)化軟件在降低抽余油中苯質(zhì)量分?jǐn)?shù)的運(yùn)用

從C6、C7餾分中分離芳烴的生產(chǎn)過(guò)程涉及到兩個(gè)塔：苯抽提蒸餾塔(DA-4503塔)和苯產(chǎn)品汽提塔(DA-4504塔)。以抽余油中的苯質(zhì)量分?jǐn)?shù)作為目標(biāo)變量(y)，選取11個(gè)工藝參數(shù)作為相關(guān)變量(xi)：苯餾分進(jìn)料量(x1)、溶劑比(x2)、進(jìn)料苯質(zhì)量分?jǐn)?shù)(x3)、DA-4503塔回流比(x4)、DA-4503塔塔頂壓力(x5)、DA-4503塔塔底溫度(x6)、DA-4503塔靈敏溫度(x7)、苯餾分溫度(x8)、溶劑溫度(x9)、DA-4504塔塔底溫度(x10)、DA-4503塔真空度(x11)。收集上述變量數(shù)據(jù)，我們得到DMOS軟件數(shù)據(jù)處理報(bào)告，并提出優(yōu)化操作調(diào)整建議，實(shí)施后取得了較好的效果。

1.1 優(yōu)化前的生產(chǎn)狀況

為了對(duì)優(yōu)化前變量之間的相關(guān)性等有所了解，以確定采用數(shù)據(jù)分析的方法，首先用DMOS優(yōu)化軟件的數(shù)理統(tǒng)計(jì)模塊對(duì)生產(chǎn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析(見表1)。

從表1中可以看出：此生產(chǎn)狀況下抽余油中苯質(zhì)量分?jǐn)?shù)尚不夠穩(wěn)定，平均值為2.3%，標(biāo)準(zhǔn)差為1.1%。抽余油中苯質(zhì)量分?jǐn)?shù)在2%以上的樣本點(diǎn)占全體樣本的56%，抽余油中苯質(zhì)量分?jǐn)?shù)在2%以下的樣本點(diǎn)占全體樣本的44%，反映出生產(chǎn)狀況不是很好，分布離散。優(yōu)化的目標(biāo)是降低抽余油中的苯質(zhì)量分?jǐn)?shù)的均值，并減少標(biāo)準(zhǔn)差，使生產(chǎn)過(guò)程更平穩(wěn)。

表1 優(yōu)化前抽余油中苯質(zhì)量分?jǐn)?shù)

為了進(jìn)一步研究變量之間、以及變量與目標(biāo)之間的關(guān)系，用DMOS軟件進(jìn)行了相關(guān)性分析，并列出變量之間的相關(guān)因子(見表2)。

表2 相關(guān)因子

表2中目標(biāo)變量y與x5、x8、x10的相關(guān)因子分別為：0.662，0.804，-0.708，這些變量與產(chǎn)品質(zhì)量有著較大的線性相關(guān)性。另外，許多變量之間的線性相關(guān)系數(shù)較大，如：x5與x6、x5與x8、x5與x10、x6與x8、x6與x9、x6與x10、x8與x10。數(shù)據(jù)中存在強(qiáng)共線性問(wèn)題，用DMOS優(yōu)化軟件中模式識(shí)別方法進(jìn)行數(shù)據(jù)分析。

1.2 生產(chǎn)數(shù)據(jù)的模式識(shí)別分析

運(yùn)用數(shù)理統(tǒng)計(jì)方法對(duì)數(shù)據(jù)進(jìn)行分析，可以深層次幫助人們更好地了解生產(chǎn)情況，獲得一些有意義的信息。但要進(jìn)行深入地分析，并獲得生產(chǎn)優(yōu)化的規(guī)律，需要用到DMOS的模式識(shí)別方法。

模式識(shí)別方法的原理可以用高維空間的聚類進(jìn)行分析和理解。既然優(yōu)化問(wèn)題常常是一個(gè)多因素問(wèn)題，就可以將所有相關(guān)變量構(gòu)建一個(gè)高維空間，并把收集到的數(shù)據(jù)在這個(gè)高維空間中表達(dá)。根據(jù)“物以類聚”的道理，優(yōu)類樣本點(diǎn)必然會(huì)形成優(yōu)化區(qū)，通過(guò)合適的算法，將高維數(shù)據(jù)空間投影到兩維平面(稱為特征面)，就可以直觀地在平面圖形上看到樣本點(diǎn)的分布和優(yōu)化區(qū)形狀，進(jìn)而研究?jī)?yōu)化操作的規(guī)律。

(1)建訓(xùn)練集

將優(yōu)化前3個(gè)月數(shù)據(jù)構(gòu)成訓(xùn)練集，定義抽余油中苯質(zhì)量分?jǐn)?shù)<2%的為優(yōu)類樣本，抽余油中苯質(zhì)量分?jǐn)?shù)>2%的為差類樣本。訓(xùn)練集中優(yōu)類樣本數(shù)33個(gè)，占41%；差類樣本數(shù)41個(gè)，占59%。

(2)建模

圖1 特征圖

從圖1中可以看出：優(yōu)類點(diǎn)明顯聚集在一起(有32個(gè)樣本點(diǎn)，占84%，落在優(yōu)化區(qū)內(nèi))，形成了一個(gè)優(yōu)化區(qū)，將生產(chǎn)操作參數(shù)控制在特征平面的優(yōu)化區(qū)內(nèi)，則有84%的概率使抽余油中苯質(zhì)量分?jǐn)?shù)小于2%，與此前的生產(chǎn)相比，優(yōu)類點(diǎn)概率提高40%。

(3)操作參數(shù)對(duì)目標(biāo)變量的作用和影響

特征圖只表示了樣本點(diǎn)在高維數(shù)學(xué)空間中的分布規(guī)律，它可以告訴我們優(yōu)化規(guī)律的存在，但不能直接提示操作的方向。為了與工藝參數(shù)聯(lián)系起來(lái)，使用軟件的權(quán)重圖工具分析，通過(guò)這個(gè)圖形工具可以將特征圖和工藝參數(shù)聯(lián)系起來(lái)，分析各個(gè)工藝參數(shù)在優(yōu)化生產(chǎn)中的作用和影響。

權(quán)重圖在橫坐標(biāo)上標(biāo)出了10個(gè)變量，根據(jù)優(yōu)化軟件的研判方法，方塊在橫坐標(biāo)上方表示該變量若增加有利于接近目標(biāo)值(苯質(zhì)量分?jǐn)?shù)降低)，反之，方塊在橫坐標(biāo)下方表示該變量若減小有利于接近目標(biāo)值(苯質(zhì)量分?jǐn)?shù)降低)，并且根據(jù)方塊的高度判斷各個(gè)變量對(duì)目標(biāo)影響的力度(見圖2)。

圖2 權(quán)重圖

通過(guò)對(duì)特征圖和權(quán)重圖的分析，對(duì)各個(gè)主要變量的作用得到以下初步結(jié)論：

x10、x6、x8、x5、x4、x9等變量對(duì)目標(biāo)值的影響有較大作用。x1、x2、x7、x11等變量對(duì)目標(biāo)值的影響作用較小，但可作參考。

(4)優(yōu)化方案

根據(jù)實(shí)際生產(chǎn)情況，參照裝置設(shè)計(jì)工藝條件，通過(guò)運(yùn)用DMOS軟件的特征圖和權(quán)重圖對(duì)數(shù)據(jù)進(jìn)行分析，運(yùn)用研判法提出優(yōu)化操作建議(見表3)。

表3 優(yōu)化操作建議

注：上下箭頭表示該變量繼續(xù)增大或減少，對(duì)優(yōu)化生產(chǎn)有利，可視情況繼續(xù)調(diào)。

2 應(yīng)用效果

上述操作建議在隨后的生產(chǎn)中進(jìn)行實(shí)施，運(yùn)行一段時(shí)間后，收集優(yōu)化實(shí)施后的數(shù)據(jù)，并對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析(見表4)。

表4 優(yōu)化實(shí)施后的抽余油中苯質(zhì)量分?jǐn)?shù)

如表4所示：優(yōu)化后的生產(chǎn)狀況下產(chǎn)品質(zhì)量穩(wěn)定，抽余油中苯質(zhì)量分?jǐn)?shù)的平均值為0.97%，標(biāo)準(zhǔn)差為0.69%。通過(guò)與優(yōu)化前產(chǎn)品質(zhì)量作比較，抽余油中苯質(zhì)量分?jǐn)?shù)的平均值由2.3%降至0.97%，降低了1.33個(gè)百分點(diǎn)；標(biāo)準(zhǔn)差由1.1%降至0.69%，下降了0.42個(gè)百分點(diǎn)。抽余油的設(shè)計(jì)產(chǎn)能為58.8 kt/a，由此推算，優(yōu)化后石油苯產(chǎn)量可增產(chǎn)782 t/a。

3 結(jié)論

(1)應(yīng)用DMOS優(yōu)化軟件對(duì)芳烴抽提裝置抽余油中苯質(zhì)量分?jǐn)?shù)實(shí)施優(yōu)化后，抽余油中苯質(zhì)量分?jǐn)?shù)明顯下降，增效明顯。

(2)運(yùn)用DMOS優(yōu)化軟件指導(dǎo)優(yōu)化生產(chǎn)，投入少、見效快、無(wú)風(fēng)險(xiǎn)，是數(shù)據(jù)處理實(shí)現(xiàn)優(yōu)化生產(chǎn)的有效工具之一，具有實(shí)用價(jià)值。

(3)DMOS優(yōu)化軟件適用于內(nèi)部機(jī)理復(fù)雜、影響因素多的生產(chǎn)裝置，可以有效地處理復(fù)雜的生產(chǎn)數(shù)據(jù)，對(duì)優(yōu)化生產(chǎn)操作具有指導(dǎo)意義；利用現(xiàn)有的數(shù)據(jù)資源，通過(guò)優(yōu)化生產(chǎn)達(dá)到增效，同樣對(duì)各類煉油化工裝置的優(yōu)化生產(chǎn)具有推廣價(jià)值。