唐國(guó)維,趙 雪

(東北石油大學(xué)計(jì)算機(jī)與信息技術(shù)學(xué)院,黑龍江大慶163318)

在煉油生產(chǎn)過(guò)程中,油品調(diào)和是一項(xiàng)重要工序[1]。近些年,國(guó)內(nèi)對(duì)于高標(biāo)號(hào)汽油的各種環(huán)保油品的需求較多,而大多數(shù)煉油廠用現(xiàn)行技術(shù)生產(chǎn)的成品油一次合格率不高,不能滿足用戶的需求。曹萃文等[2]引入油品非線性調(diào)和模型,提出了集成控制原油線性調(diào)和、基于訂單的短期生產(chǎn)計(jì)劃、成品油非線性調(diào)和及庫(kù)存管理的算法;王繼東等[3]采用遺傳算法解決油品調(diào)和調(diào)度問(wèn)題;趙小強(qiáng)[4]針對(duì)該調(diào)和問(wèn)題給出了一種基于基本微粒群算法的解決方法,但粒子群算法收斂速度較慢且容易陷入局部極值點(diǎn)。下面,筆者引入一種動(dòng)態(tài)改變權(quán)重的粒子群優(yōu)化算法,對(duì)煉油企業(yè)生產(chǎn)調(diào)度的油品調(diào)和決策進(jìn)行了研究。

1 成品油調(diào)和優(yōu)化模型的建立

煉油生產(chǎn)裝置規(guī)模龐大,工藝復(fù)雜,影響煉油生產(chǎn)的參數(shù)和變量繁多。筆者結(jié)合煉油工藝和某煉油廠實(shí)際生產(chǎn)情況,選取對(duì)煉油生產(chǎn)有決定意義裝置的變量和參數(shù)作為決策量。一般油品調(diào)和模型以費(fèi)用最低為目標(biāo)。油品調(diào)和模型[5]可歸結(jié)為數(shù)學(xué)規(guī)劃問(wèn)題。目標(biāo)函數(shù)由調(diào)和產(chǎn)品的價(jià)值、調(diào)和組分的價(jià)值及調(diào)和操作費(fèi)用組成。其目標(biāo)函數(shù)為:

式中,Pi(i=1,…,I)表示調(diào)和產(chǎn)品的種類;Cj(j=1,…,J)表示調(diào)和組分的種類;t(t=1,…,T)表示時(shí)間間隔長(zhǎng)度;BCj,t、B Pi,t分別為Cj和Pi的生產(chǎn)量;D Pi,Cj為Cj調(diào)和成Pi的成本;F為調(diào)和產(chǎn)品利潤(rùn);GPj為Pi的單位市場(chǎng)價(jià);GCj為Cj的單位市場(chǎng)價(jià);λPi,Cj=0或 1,表示是否有參與調(diào)和。

約束條件包括如下內(nèi)容:①產(chǎn)品的生產(chǎn)量約束(min B Pi,t≤BPi,t≤max BPi,t),min BPi,t、max BPi,t分別為BPi,t的最小和最大生產(chǎn)量;②產(chǎn)品和調(diào)組份平衡關(guān)系;③產(chǎn)品庫(kù)存量平衡關(guān)系(VPi,t=V Pi,i-1+BPi,t-XPi,t),XPi,t為P i的需求量;VCj,t、VPi,t分別為Cj和P i的庫(kù)存量;④產(chǎn)品的庫(kù)存量約束分別為VPi,t最小和最大庫(kù)存量;⑤產(chǎn)品的質(zhì)量滿足(Qm,Pi,t=f q(BPi,Cj,t,Qm,Cj,t,NPi,Cj,t)),NPi,Cj,t為Cj用于P i的比例,BPi,Cj,t為Pi所用的Cj的生產(chǎn)量,Qm(m=1,…,M)表示調(diào)和產(chǎn)品屬性種類;Qm,Cj,t、Qm,Pi,t分別為Cj和Pi的屬性值;⑥產(chǎn)品質(zhì)量約束(min Qm,Pi≤Qm,Pi,t≤max Qm,Pi),min Qm,Pi、max Qm,Pi分別為Qm,Pi,t的屬性的最小值和最大值;;⑧組份庫(kù)存量約束(min VCj,t≤VCj,t≤max VCj,t),min VCj,t、max VCj,t分別為VCj,t的最小和最大庫(kù)存量;⑨組份質(zhì)量滿足(min Qm,Cj≤Qm,Cj,t≤max Qm,Cj),min Qm,Cj和⑦組分庫(kù)存量滿足max Qm,Cj為Qm,Cj,t屬性的最小和最大值;⑩調(diào)和組份滿足

在調(diào)和模型中,對(duì)約束條件的處理很重要。為此,筆者采用罰函數(shù)方法,把優(yōu)化計(jì)算中約束優(yōu)化問(wèn)題轉(zhuǎn)化為無(wú)約束優(yōu)化求解。對(duì)于非線性油品調(diào)和優(yōu)化問(wèn)題采用了如下罰函數(shù):

式中,qi(x)=max{0,gi(x)};Q(qi(x))為多段賦值函數(shù),Q(qi(x))的選取規(guī)則為Q(qi(x))=10(qi(x)＜1)或Q(qi(x))=100(qi(x)≥1);r(qi(x))為罰函數(shù)的冪函數(shù),r(qi(x))的選取規(guī)則為r(qi(x))=1(qi(x)＜1)或r(qi(x))=2(qi(x)≥1);h(x)為Rn上的連續(xù)函數(shù)。

2 粒子群算法描述

粒子群優(yōu)化算法 (PSO)是一種基于群體智能的優(yōu)化技術(shù)[6],其從隨機(jī)解出發(fā),通過(guò)迭代尋找最優(yōu)解,通過(guò)適應(yīng)度來(lái)評(píng)價(jià)解的品質(zhì),通過(guò)追隨當(dāng)前搜索到的最優(yōu)值來(lái)尋找全局最優(yōu)。

油品調(diào)和是直接影響煉油廠經(jīng)濟(jì)效益的一個(gè)關(guān)鍵環(huán)節(jié)。在油品調(diào)和過(guò)程中,成品油性質(zhì)指標(biāo)通常是組分性質(zhì)的非線性函數(shù)。筆者在基本粒子群算法的基礎(chǔ)上,動(dòng)態(tài)改變慣性權(quán)重[7],使迭代中的慣性權(quán)重ω隨粒子的位置和目標(biāo)函數(shù)的變化而變化。

在粒子群算法中,假設(shè)Xi=(x i1,xi2,…,xiD)T表示第i個(gè)粒子,其中D是粒子的維數(shù)。其經(jīng)歷的最好位置表示為Pi=(pi1,pi2,…,piD)T,而整個(gè)群體經(jīng)歷的最好位置表示為Pg=(pg1,pg2,…,pgD)T,粒子i的速度為Vi=(vi1,v i2,…,v iD)T。按追隨當(dāng)前最優(yōu)粒子的原理,粒子i將按式(3)改變速度,式(4)改變位置[8]:

式中,n為當(dāng)前的進(jìn)化代數(shù);c1、c2為學(xué)習(xí)因子;rand()為分布于(0,1)的隨機(jī)數(shù)。

采用粒子平均維間距來(lái)描述種群的維多樣性。設(shè)S為種群粒子數(shù)目,xid為第i個(gè)粒子在第d維上的值,ˉx d表示所有粒子第d維上坐標(biāo)值均值,則定義第d維的種群多樣性為:

式(5)用平均維間距M(d)描述了種群某一維粒子相互之間分布的離散程度,設(shè)定ω=0.8。在多樣性最差、某一維所有點(diǎn)聚集在一起時(shí),可設(shè)定ω=0.1進(jìn)行精細(xì)的開(kāi)發(fā)。因此,得到ω的計(jì)算公式:

式中,K為壓縮因子,在 (0,1)間取值。具體算法的步驟如下:①設(shè)定算法相關(guān)參數(shù)值;②初始化種群中每個(gè)粒子的位置和速度;③評(píng)價(jià)每一個(gè)粒子,更新粒子的最優(yōu)位置與群體的最優(yōu)位置;④采用上述方法動(dòng)態(tài)調(diào)整慣性權(quán)重值,并按式 (3)和式 (4)更新粒子的位置和速度;⑤判斷算法是否滿足中止條件,若滿足則算法結(jié)束;否則,轉(zhuǎn)回式 (6)。

3 成品油調(diào)和優(yōu)化過(guò)程的實(shí)現(xiàn)

以某煉油廠汽油調(diào)和為例,一重催、二重催生產(chǎn)的組份汽油,經(jīng)過(guò)加入少量的重整汽油、乙烯芳烴和MTBE,就可以調(diào)和成各種牌號(hào)的乙醇汽油組份油。組分油和成品油的性能指標(biāo)如表1和表2所示。

表1 組分油性能指標(biāo)

表2 成品油性能指標(biāo)

調(diào)和乙醇汽油組分油90#、93#、97#的產(chǎn)量分別為X 1、X 2、X3;調(diào)和乙醇汽油組分油90#使用一重催生產(chǎn)的組分汽油、使用二重催生產(chǎn)的組分汽油、重整汽油、乙烯芳烴、MTBE的量分別為X 11、X 12、X 13、X14、X 15;93#使用一重催生產(chǎn)的組分汽油、使用二重催生產(chǎn)的組分汽油、重整汽油、乙烯芳烴、MTBE的量分別為X 21、X 22、X23、X 24、X 25;97#使用一重催生產(chǎn)的組分汽油、使用二重催生產(chǎn)的組分汽油、重整汽油、乙烯芳烴、MTBE的量為 X31、X32、X33、X34、X35。

則:

物料平衡約束為:

收入量約束為:

產(chǎn)量約束為:

辛烷值約束為:

為了方便比較各算法的性能,將相關(guān)參數(shù)設(shè)置如下:①基本粒子群算法中,c1=c2=2;慣性權(quán)重因子ω采用線性遞減的方式從0.9減少到0.4。②在筆者引入的改進(jìn)算法中,c1=c2=1.62;慣性權(quán)重按式(6)進(jìn)行動(dòng)態(tài)調(diào)整。2種算法分別進(jìn)行最大迭代次數(shù)為G=1500和2000,種群規(guī)模P=40和80。為了減少偶然性的影響,在迭代次數(shù)相同的條件下,各算法測(cè)試均運(yùn)行50次,然后取平均值。其結(jié)果如表3所示。

從表3可以看出,與基本算法相比,通過(guò)改進(jìn)算法得到的調(diào)和成品油的利潤(rùn)提高了10400元。這是因?yàn)楦倪M(jìn)算法略微加快了算法的收斂速度,同時(shí)有利于跳出局部極值,具有較好的全局搜索能力,從而保證了乙醇汽油的質(zhì)量,使煉油廠的經(jīng)濟(jì)效益得到提高。

表3 算法計(jì)算結(jié)果比較

4 結(jié) 語(yǔ)

結(jié)合煉油廠實(shí)際生產(chǎn)情況,選取有決定意義的裝置的變量和參數(shù),建立油品調(diào)和模型,以煉油廠最大利潤(rùn)為目標(biāo),引入改變權(quán)重的粒子群優(yōu)化算法來(lái)解決,既保證了尋優(yōu)過(guò)程中粒子群體的多樣性,又加快了算法的全局收斂能力,保證了產(chǎn)品的質(zhì)量,使煉油廠效益得到提高。因此,粒子群優(yōu)化算法對(duì)煉油廠改進(jìn)油品調(diào)和調(diào)度方法、實(shí)施調(diào)和調(diào)度優(yōu)化有現(xiàn)實(shí)意義。

[1]蔡智.油品調(diào)和技術(shù)-煉油工業(yè)技術(shù) [M].北京:中國(guó)石化出版社,2006.

[2]曹萃文,顧幸生.非線性油品優(yōu)化調(diào)和在煉油廠中的應(yīng)用 [J].化工自動(dòng)化及儀表,2005,32(6):12-15.

[3]王繼東,王萬(wàn)良.基于遺傳算法的汽油調(diào)和生產(chǎn)優(yōu)化研究 [J].化工自動(dòng)化及儀表,2005,32(1):6-9.

[4]趙小強(qiáng),榮岡.基于微粒群優(yōu)化算法的不確定性調(diào)和調(diào)度 [J].中國(guó)化學(xué)工程學(xué)報(bào),2005,13(4):535-541.

[5]薛美盛,李祖奎,吳剛,等.成品油調(diào)和調(diào)度優(yōu)化模型及其應(yīng)用 [J].石油煉制與化工,2005(11):64-48.

[6]Kennedy J,Eberhart R C.Partical Swarm Optimization[A],Proceeding of the 1995 IEEE International Conference on Neural Network[C].Perth,Australia,1995:1942-1948.

[7]王啟付,王戰(zhàn)江,王書亭.一種動(dòng)態(tài)改變慣性權(quán)重的粒子群優(yōu)化算法 [J].中國(guó)機(jī)械工程,2005,16(11):945-948.

[8]龍文,梁昔明,董淑華,等.動(dòng)態(tài)調(diào)整慣性權(quán)重的粒子群算法 [J].計(jì)算機(jī)應(yīng)用,2009,29(8):2240-2242.