• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于偏最小二乘的Kriging代理模型在加氫裂化建模中的應(yīng)用

    2017-07-01 22:37:30鐘偉民
    關(guān)鍵詞:航煤加氫裂化柴油

    喬 成, 鐘偉民, 范 琛

    (華東理工大學(xué)化工過(guò)程先進(jìn)控制和優(yōu)化技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200237)

    ?

    基于偏最小二乘的Kriging代理模型在加氫裂化建模中的應(yīng)用

    喬 成, 鐘偉民, 范 琛

    (華東理工大學(xué)化工過(guò)程先進(jìn)控制和優(yōu)化技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,上海 200237)

    提出了一種改進(jìn)的代理模型方法 (Kriging with Partial Least Squares,KPLS)。該方法在Kriging模型的基礎(chǔ)上引入偏最小二乘的思想,利用偏最小二乘方法構(gòu)建新的Kriging模型的高斯核函數(shù)。將該模型應(yīng)用于加氫裂化過(guò)程建模,有效地提高了航煤、柴油質(zhì)量收率的預(yù)測(cè)精度。采用 GLAMP(Global and local search strategy)優(yōu)化算法對(duì)建立的KPLS模型進(jìn)行優(yōu)化,仿真結(jié)果顯示航煤、柴油質(zhì)量收率得到了顯著提升。

    加氫裂化; Kriging代理模型; 偏最小二乘; 收率預(yù)測(cè); GLAMP優(yōu)化算法

    加氫裂化是一種轉(zhuǎn)換率很高的煉油二次加工過(guò)程,在加熱、氫氣和催化劑存在的條件下,使重質(zhì)油發(fā)生加氫、裂化和異構(gòu)化反應(yīng)。其中加氫反應(yīng)用來(lái)去除碳?xì)浠衔镏械牧蚝偷入s質(zhì),并在這個(gè)過(guò)程中產(chǎn)生飽和的碳?xì)浠衔?。在高溫、高壓以及催化劑的條件下,裂化和異構(gòu)化反應(yīng)將重質(zhì)油生成高十六烷值的柴油和航煤等輕質(zhì)油[1]。由于加氫裂化生產(chǎn)工藝靈活多變,產(chǎn)品油收率高、質(zhì)量好,且硫、氮含量低,已成為重要的原油二次加工手段[2]。

    加氫裂化模型一般采用機(jī)理建模法,但是由于未知的化學(xué)反應(yīng)、非線性關(guān)系以及大量的變量使得建模過(guò)程變得十分復(fù)雜[3]。繁重和耗時(shí)的計(jì)算以及難以在機(jī)理模型的基礎(chǔ)上進(jìn)行操作參數(shù)的優(yōu)化等問(wèn)題都是機(jī)理模型存在的缺點(diǎn)[4],但是Kriging代理模型卻能有效替代機(jī)理模型去處理這一系列復(fù)雜難題,它只需要輸入輸出的數(shù)據(jù)而不用了解詳細(xì)的系統(tǒng)知識(shí)。由于Kriging模型運(yùn)行時(shí)間短,對(duì)多峰值、非線性函數(shù)具有較好的預(yù)測(cè)能力,因此在眾多學(xué)科中得到了越來(lái)越多的應(yīng)用。Han等[5]為了解決傳統(tǒng)的CFD模型對(duì)機(jī)翼氣動(dòng)系數(shù)、極限阻力等數(shù)值求取困難、耗時(shí)長(zhǎng)的問(wèn)題,提出了一種梯度增強(qiáng)的Kriging(GEK)代理模型,該模型在傳統(tǒng)的Kriging模型中引入了梯度信息。實(shí)際結(jié)果表明,該模型能夠在短時(shí)間內(nèi)對(duì)各種工況環(huán)境下的機(jī)翼氣動(dòng)系數(shù)、極限阻力做出準(zhǔn)確預(yù)測(cè)。在化工領(lǐng)域中,如何通過(guò)對(duì)精餾塔改造降低生產(chǎn)成本一直是困擾眾多學(xué)者的問(wèn)題,基于傳統(tǒng)的Hysys機(jī)理模型難以對(duì)這種混合整數(shù)非線性問(wèn)題進(jìn)行優(yōu)化求解,故Quirante等[6]通過(guò)建立精餾塔的Kriging代理模型,并在代理模型基礎(chǔ)上采用GAMS-BARON優(yōu)化方法對(duì)精餾塔的塔高、塔板數(shù)、回流比、再沸比等參數(shù)進(jìn)行優(yōu)化,最終達(dá)到降低設(shè)備成本和操作成本的目的。祁榮賓等[7]為了解決PX氧化過(guò)程優(yōu)化時(shí)間長(zhǎng)、效率低的問(wèn)題,提出了建立Kriging代理模型代替流程仿真模型。結(jié)果表明,PX氧化過(guò)程代理模型能夠準(zhǔn)確地預(yù)測(cè)醋酸燃燒損失、PX燃燒損失、4-CBA含量,并在代理模型的基礎(chǔ)上利用粒子群算法進(jìn)行優(yōu)化,不僅能夠獲得全局最優(yōu)解,還能極大地縮短搜索時(shí)間。

    由于加氫裂化反應(yīng)過(guò)程建模中輸入端變量較多,Han等[5]采用的梯度增強(qiáng)的Kriging方法只適合在低維度建模優(yōu)化中使用,隨著維度的不斷增加,其高斯相關(guān)矩陣容易陷入病態(tài)化,影響預(yù)測(cè)精度。為了提高Kriging模型在高維建模中的預(yù)測(cè)精度,結(jié)合加氫裂化反應(yīng)器特點(diǎn),本文提出了一種Kriging模型與偏最小二乘法(PLS)相結(jié)合的建模方法,能夠有效地提高石油產(chǎn)品收率的預(yù)測(cè)精度。同時(shí)結(jié)合GLAMP優(yōu)化算法,對(duì)操作參數(shù)和原料性質(zhì)進(jìn)行優(yōu)化,提高航煤、柴油的質(zhì)量收率。

    1 KPLS代理模型的結(jié)構(gòu)

    1.1 Kriging模型

    Kriging也被稱(chēng)為數(shù)學(xué)回歸過(guò)程[8],是一種基于統(tǒng)計(jì)學(xué)的插值技術(shù)。該技術(shù)利用最優(yōu)線性無(wú)偏估計(jì)方法,通過(guò)對(duì)先驗(yàn)函數(shù)值的分布對(duì)未知點(diǎn)的函數(shù)值進(jìn)行預(yù)測(cè)。

    為了實(shí)現(xiàn)對(duì)函數(shù)f(x)的近似,F(x)被定義為其響應(yīng)值,公式如下:

    F(x)=u+Z(x)

    (1)

    其中:x是d維的向量;u為一個(gè)常數(shù);Z(x)為一個(gè)隨機(jī)過(guò)程,其期望為0,方差是σ2,協(xié)方差為σ2R(θ,x,x′),R(θ,x,x′)是相關(guān)函數(shù),一般采用高斯函數(shù),定義如下:

    R(θi,xi,x′i)=exp(-θi|xi-x′i|2)

    (2)

    將上述公式整合,獲得Kriging模型的高斯協(xié)方差核函數(shù)如下:

    ?θi∈R+

    (3)

    (4)

    預(yù)測(cè)方差為

    (5)

    (6)

    1.2 偏最小二乘法

    偏最小二乘(PLS)是一種將輸入變量映射到新的空間來(lái)建立輸入輸出變量之間聯(lián)系的多元統(tǒng)計(jì)方法,通過(guò)搜索新的多維坐標(biāo)方向來(lái)更好地解釋輸出向量y的特性[9]。為了最好地獲取新的多維坐標(biāo)方向,對(duì)采樣矩陣X(n×d)及其響應(yīng)向量y求取第一主成分t(1),其計(jì)算公式如下:

    t(1)=Xw(1)

    (7)

    (8)

    其中w(1)是矩陣XTyyTX最大特征值對(duì)應(yīng)的特征向量,包含矩陣X在第一主成分中的權(quán)重。在求得第一主成分t(1)后,利用式(9)計(jì)算矩陣X和向量y的殘差:

    (9)

    這里的p(1)(l×d向量)、c1分別是X和y在第一主元t(1)上的自回歸系數(shù),X(0)=X,y(0)=y。

    在計(jì)算第二主成分t(2)過(guò)程中,只需將式(9)中的X和y分別替代成其殘差X(1)、y(1),其他的主成分可以依次通過(guò)相同的方法迭代求得,因此主成分t(l)可表示為

    (10)

    W*=W(PTW)-1

    (11)

    1.3 構(gòu)建KPLS模型

    (12)

    (13)

    在構(gòu)建KPLS模型時(shí),將式(13)得到的新高斯協(xié)方差核函數(shù)序列代替原Kriging模型中公式(3)的核函數(shù),步驟如下:

    (1) 首先獲取足夠的初始化實(shí)驗(yàn)數(shù)據(jù)(X,y)。

    (2) 推導(dǎo)出預(yù)測(cè)輸出值的公式,確定核函數(shù)k的結(jié)構(gòu)(這里定義為高斯指數(shù)函數(shù)結(jié)構(gòu))。

    (3) 利用PLS方法構(gòu)建新的高斯協(xié)方差核函數(shù)序列kkpls1:h:

    (a) 初始化PLS算法,初始迭代值l=1;

    (b) 如果l≠1,通過(guò)公式(9) 計(jì)算殘差Xl-1和yl-1;

    (d) 通過(guò)式(13)定義新的核函數(shù)kkpls1:h;

    (e) 如果迭代參數(shù)l達(dá)到上限,返回步驟(3),否則繼續(xù)循環(huán);

    (f) 更新參數(shù)l=l+1。

    (4) 利用極大似然估計(jì)的方法計(jì)算出參數(shù)θ。

    圖1所示為建立KPLS模型的流程圖。

    圖1 構(gòu)建KPLS模型的流程圖

    2 建立加氫裂化反應(yīng)過(guò)程KPLS代理模型

    2.1 加氫裂化反應(yīng)過(guò)程描述

    圖2示出了典型的單段串聯(lián)加氫裂化的工藝流程。原料經(jīng)泵(P1001A/B)升壓,與氫氣混合后進(jìn)入加熱爐(F1001)加熱,隨后進(jìn)入加氫精制反應(yīng)器(R1001)進(jìn)行加氫精制反應(yīng)。加氫精制反應(yīng)器中設(shè)有多個(gè)催化劑床層,床層間設(shè)急冷氫注入設(shè)施。反應(yīng)流出物進(jìn)入加氫裂化反應(yīng)器(R1002)進(jìn)行加氫裂化反應(yīng),在兩個(gè)反應(yīng)器之間設(shè)急冷氫注入點(diǎn),同樣加氫裂化反應(yīng)器中也擁有多個(gè)催化劑床層,床層間設(shè)有急冷氫注入設(shè)施。由加氫反應(yīng)器出來(lái)的反應(yīng)流出物經(jīng)換熱器(E1001)溫度降至260 ℃左右,進(jìn)入熱高壓分離器(V1002)。熱高分氣體再經(jīng)熱高分氣空冷器冷卻至49 ℃左右進(jìn)入冷高壓分離器(V1003)。熱高分氣冷卻后在冷高壓分離器中進(jìn)行油、水、氣三相分離。頂部分離出的氫氣進(jìn)入循環(huán)氫脫硫塔后,經(jīng)壓縮機(jī)(C1001)升壓后返回系統(tǒng)使用。底部流出的生成油減壓后進(jìn)入低壓分離器(V1006),脫除水,并釋放出部分溶解氣體(燃料氣)。加熱后的生成油被送入主汽提塔(T1001),在1.0~1.2 MPa蒸出液化氣,塔底液體加熱至320 ℃左右后進(jìn)入分餾塔,得到石腦油、航空煤油、低凝柴油和塔底油(尾油)。

    圖2 單段串聯(lián)加氫裂化簡(jiǎn)化工藝流程圖Fig.2 Process flow diagram of single series hydrocracking

    2.2 KPLS代理模型變量的選取

    加氫裂化在石油產(chǎn)品加工過(guò)程中,產(chǎn)品的收率主要受到R1001和R1002反應(yīng)器床層溫度、壓力、氫油體積比以及原料性質(zhì)的影響。根據(jù)實(shí)際情況分析,本文選取了10個(gè)對(duì)產(chǎn)品質(zhì)量收率有顯著影響的輸入變量,主要包括:原料密度、氫油體積比、原料氮含量、原料硫含量、R1001床層1溫度、R1001床層2溫度、R1001入口壓力、R1002床層1溫度、R1002床層2溫度、R1002床層3溫度。由于加氫裂化的主要產(chǎn)品是柴油和航煤,因此輸出端的變量選擇航煤、柴油的質(zhì)量收率。其他產(chǎn)品如石腦油、液化氣等由于篇幅的限制不一一介紹。

    2.3 KPLS代理模型的建立與驗(yàn)證

    在確定輸入輸出變量后,通過(guò)拉丁立方采樣的方法(LHS)獲取200個(gè)采樣點(diǎn),其中每一個(gè)輸入變量的取值范圍都是根據(jù)某石化企業(yè)實(shí)際工況確定的。在建模過(guò)程中,以Kriging模型為基礎(chǔ),利用PLS方法分別對(duì)柴油、航煤質(zhì)量收率構(gòu)建3個(gè)包含不同主成分的KPLS模型,并用留一法交叉驗(yàn)證模型的有效性。預(yù)測(cè)誤差結(jié)果見(jiàn)表1。

    表1 加氫裂化反應(yīng)過(guò)程KPLS代理模型對(duì)航煤、柴油質(zhì)量收率的預(yù)測(cè)誤差

    表1中的統(tǒng)計(jì)結(jié)果表明,加氫裂化反應(yīng)過(guò)程KPLS代理模型對(duì)航煤、柴油收率的預(yù)測(cè)精度比傳統(tǒng)的Kriging模型有了大幅的提高,其中以航煤質(zhì)量收率為預(yù)測(cè)目標(biāo)的KPLS代理模型在包含3個(gè)主成分時(shí)預(yù)測(cè)誤差最小,誤差為2.4%。而以柴油質(zhì)量收率為預(yù)測(cè)目標(biāo)的KPLS代理模型在包含2個(gè)主成分時(shí)預(yù)測(cè)誤差最小,誤差為2.2%。因此,選擇KPLS(3 Components)作為航煤收率預(yù)測(cè)的加氫裂化反應(yīng)過(guò)程代理模型,選擇KPLS(2 Components)作為柴油收率預(yù)測(cè)的加氫裂化反應(yīng)過(guò)程代理模型。表1中誤差計(jì)算公式如下:

    (14)

    表2顯示了Hysys機(jī)理模型與各Kriging代理模型在預(yù)測(cè)航煤、柴油質(zhì)量收率時(shí)所花費(fèi)CPU的運(yùn)行時(shí)間。表中的數(shù)據(jù)表明,在進(jìn)行100組數(shù)據(jù)預(yù)測(cè)過(guò)程中,代理模型能夠極大地提高航煤、柴油收率的預(yù)測(cè)效率,并且引入偏最小二乘方法的KPLS模型與傳統(tǒng)Kriging模型在預(yù)測(cè)時(shí)間上基本一致,并沒(méi)有出現(xiàn)明顯的差異。

    表2 Hysys機(jī)理模型與各代理模型預(yù)測(cè)100組數(shù)據(jù)的平均時(shí)間(50組航煤,50組柴油)

    3 加氫裂化反應(yīng)過(guò)程KPLS代理模型優(yōu)化

    GLAMP算法是一種平衡全局搜索和局部搜索的代理模型優(yōu)化算法,該算法結(jié)合KPLS模型曲線光滑、梯度易求的特點(diǎn),采用Dennis[13]多起點(diǎn)序列二次規(guī)劃法(SQP)來(lái)找出局部最優(yōu)點(diǎn),并利用自適應(yīng)搜索步長(zhǎng)策略(SAMP),對(duì)局部最優(yōu)點(diǎn)附近區(qū)域進(jìn)行深入搜索。而全局搜索則利用KPLS模型對(duì)未開(kāi)發(fā)區(qū)域具有出色預(yù)測(cè)性的特點(diǎn),保證了算法的穩(wěn)定性,其算法流程如圖3所示。

    圖3 GLAMP算法的流程圖Fig.3 Flow diagram of GLAMP algorithm

    在KPLS代理模型的基礎(chǔ)上,通過(guò)GLAMP優(yōu)化算法搜索最優(yōu)的操作參數(shù)和原料性質(zhì),本文將最大搜索迭代次數(shù)設(shè)置為100。為了比較GLAMP算法的優(yōu)越性,本文引入了DDS、CAND、EGO 3個(gè)常用的Kriging模型傳統(tǒng)優(yōu)化算法進(jìn)行對(duì)比。

    圖4、圖5分別示出了基于KPLS代理模型的航煤、柴油的質(zhì)量收率圖。分別采用GLAMP、EGO、CAND、DDS 4種算法對(duì)加氫裂化過(guò)程的KPLS代理模型進(jìn)行優(yōu)化,在對(duì)航煤質(zhì)量收率的優(yōu)化中,GLAMP算法的優(yōu)化效果最為有效,其次是CAND和EGO算法,表現(xiàn)最差的是DDS算法。4條曲線中,GLAMP算法無(wú)論是相同迭代次數(shù)的優(yōu)化值,還是最終的優(yōu)化值都比其他3個(gè)算法效果好。最終GLAMP算法將航煤的質(zhì)量收率由48.32%提高到51.12%,提高了近2.8%。在對(duì)柴油質(zhì)量收率優(yōu)化過(guò)程中,GLAMP算法在20代前的優(yōu)化效果并沒(méi)有CAND和DDS算法好,但是憑借著強(qiáng)大的深度搜索能力,在后期的優(yōu)化過(guò)程中取得了優(yōu)異的效果。DDS算法和CAND算法的優(yōu)化結(jié)果基本一致,而EGO算法的優(yōu)化效果最差。最終GLAMP算法將柴油的質(zhì)量收率由35.31%提高到36.50%,提高了近1.2%。表3列出了GLAMP優(yōu)化的最優(yōu)操作參數(shù)和原料性質(zhì)。

    圖4 基于KPLS模型的航煤質(zhì)量收率優(yōu)化圖Fig.4 Optimization graph of fuel mass yield based on KPLS model

    圖5 基于KPLS模型的柴油質(zhì)量收率優(yōu)化圖Fig.5 Optimization graph of diesel mass yield based on KPLS model

    表3 GLAMP優(yōu)化的最優(yōu)操作參數(shù)和原料性質(zhì)

    4 結(jié) 論

    為了解決高維狀態(tài)下模型預(yù)測(cè)精度低的問(wèn)題,在Kriging模型的基礎(chǔ)上,引入PLS思想,提出了一種改進(jìn)的代理模型——KPLS代理模型。將該代理模型用于加氫裂化反應(yīng)過(guò)程的建模中,并對(duì)航煤和柴油的收率進(jìn)行預(yù)測(cè),結(jié)果表明KPLS代理模型與傳統(tǒng)Kriging模型相比,能夠有效地減少預(yù)測(cè)誤差,提高模型預(yù)測(cè)精度,同時(shí)兼具代理模型在預(yù)測(cè)時(shí)間上的高效性,具有廣闊的應(yīng)用前景。在KPLS代理模型的基礎(chǔ)上,利用GLAMP優(yōu)化算法對(duì)航煤和柴油的質(zhì)量收率進(jìn)行優(yōu)化,獲取最優(yōu)的操作參數(shù)和原料性質(zhì),極大地提高了企業(yè)的經(jīng)濟(jì)效益。

    [1] MAPLES R E.Petroleum Refinery Process Economics[M].USA:PennWell,2000.

    [2] LYONS W C,PLISGA G J.Standard Handbook of Petroleum and Natural Gas Engineering[M].USA:Gulf Professional Publishing,2011.

    [3] BHUTANI N,RANGAIAH G P,RAY A K.First-principles,data-based,and hybrid modeling and optimization of an industrial hydrocracking unit[J].Industrial and Engineering Chemistry Research,2006,45(23):7807-7816.

    [4] 盧建翔,周華,師佳,等.工業(yè)加氫裂化反應(yīng)器模型的建立[J].石油學(xué)報(bào)(石油加工),2010,26(6):966-971.

    [5] HAN Z H,G?RTZ S,ZIMMERMANN R.Improving variable-fidelity surrogate modeling via gradient-enhanced Kriging and a generalized hybrid bridge function[J].Aerospace Science and Technology,2013,25(1):177-189.

    [6] QUIRANTE N,JAVALOYES J,CABALLERO J A.Rigorous design of distillation columns using surrogate models based on Kriging interpolation[J].AIChE Journal,2015,61(7):2169-2187.

    [7] 祁榮賓,何鎵梁,靳文浩,等.基于Kriging代理模型的PX氧化過(guò)程優(yōu)化[J].華東理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,38(6):712-717.

    [8] CRESSIE N.The origins of Kriging[J].Mathematical Geology,1990,22(3):239-252.

    [9] ALBERTO P,GONZLEZ F.Partialleast squares regression on symmetric positive-definite matrices[J].Revista Colombiana de Estadistica,2012,36(1):177-192.

    [10] MANNE R.Analysis of two partial-least-squares algorithms for multivariate calibration[J].Chemometrics & Intelligent Laboratory Systems,1987,2(1/3):187-197.

    [11] 段星辰,杜文莉.基于 bootstrap GEI 算法的碳二加氫反應(yīng)器代理模型[J].化工學(xué)報(bào),2015,66(12):4904-4909.

    [12] ALHAJREE I,ZAHEDI G,MANAN Z A,etal.Modeling and optimization of an industrial hydrocracker plant[J].Journal of Petroleum Science & Engineering,2011,78(3):627-636.

    [13] DENNIS Jr J E,SCHNABEL R B.Numerical Methods for Unconstrained Optimization and Nonlinear Equations[M].USA:Society for Industrial & Applied,1996.

    [14] 王國(guó)榮.神經(jīng)網(wǎng)絡(luò)技術(shù)在加氫裂化產(chǎn)品收率預(yù)測(cè)模型中的應(yīng)用[J].石化技術(shù)與應(yīng)用,2015,33(4):349-353.

    [15] YE K Q,LI W,SUDJIANTO A.Algorithmic construction of optimal symmetric latin hypercube designs[J].Journal of Statistical Planning & Inference,2000,90(1):145-159.

    [16] DONG H,SONG B,WANG P,etal.A kind of balance between exploitation and exploration on Kriging for global optimization of expensive functions[J].Journal of Mechanical Science & Technology,2015,29(5):2121-2133.

    Kriging Agent Model Based on Partial Least Squares in the Application of the Hydrocracking Modeling

    QIAO Cheng, ZHONG Wei-min, FAN Chen

    (Key Laboratory of Advanced Control and Optimization for Chemical Process,Ministry of Education,East China University of Science and Technology,Shanghai 200237,China)

    This paper proposes a modified agent modeling method,Kriging with partial least squares (KPLS).By means of Kriging model,we use the partial least squares method to establish a new Gaussian kernel function.Compared with the traditional Kriging model,the proposed KPLS model can effectively improve the accuracy of the fuel and diesel yield prediction.Besides,the GLAMP (global and local search strategy) search algorithm is used to optimize the KPLS model.The simulation results show that the yield of diesel and fuel is significantly improved.

    hydrocracking; Kriging surrogate model; partial least squares; yield prediction; GLAMP optimization algorithm

    1006-3080(2017)03-0383-06

    10.14135/j.cnki.1006-3080.2017.03.014

    2016-09-26

    國(guó)家自然科學(xué)基金(61422303,21376077);上海市人才發(fā)展資金;中央高?;緲I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金

    喬 成(1991-),男,碩士生,主要研究方向?yàn)楣I(yè)過(guò)程建模與優(yōu)化。

    鐘偉民,E-mail:wmzhong@ecust.edu.cn

    TQ221.242;TQ052

    A

    猜你喜歡
    航煤加氫裂化柴油
    航煤增長(zhǎng)130%
    煉油廠增加航空煤油產(chǎn)量技術(shù)措施分析
    云南化工(2021年12期)2022-01-17 08:42:10
    加氫裂化裝置脫丁烷塔頂換熱器失效原因分析及預(yù)防措施
    用履帶收割機(jī)送柴油
    上海石化首次實(shí)現(xiàn)成品航煤直接調(diào)和
    加氫裂化裝置循環(huán)氫壓縮機(jī)干氣密封失效處理
    柴油發(fā)電機(jī)負(fù)荷計(jì)算
    金陵石化新年首船航煤出口韓國(guó)
    加氫裂化工藝技術(shù)研究
    HC裝置摻煉催化柴油改善柴油品質(zhì)技術(shù)應(yīng)用
    河南科技(2014年8期)2014-02-27 14:07:59
    国产成人精品久久二区二区免费| 91麻豆av在线| 黄色丝袜av网址大全| 很黄的视频免费| 久久精品91蜜桃| 最近在线观看免费完整版| 91字幕亚洲| 啦啦啦观看免费观看视频高清| 亚洲欧美一区二区三区黑人| 最新中文字幕久久久久 | 老司机午夜十八禁免费视频| 日本a在线网址| 极品教师在线免费播放| 成年女人永久免费观看视频| 久久天堂一区二区三区四区| 少妇丰满av| 免费看a级黄色片| 亚洲国产精品久久男人天堂| 制服人妻中文乱码| 99久久精品国产亚洲精品| 中出人妻视频一区二区| 久久久久久久久中文| 免费在线观看影片大全网站| av欧美777| 日本五十路高清| 国产欧美日韩精品一区二区| 一级黄色大片毛片| av欧美777| 综合色av麻豆| 后天国语完整版免费观看| 欧美成人免费av一区二区三区| 高潮久久久久久久久久久不卡| 一个人看视频在线观看www免费 | 男女床上黄色一级片免费看| 久久亚洲精品不卡| 97人妻精品一区二区三区麻豆| 老汉色∧v一级毛片| 黄色日韩在线| 99国产精品一区二区三区| 国产成人精品久久二区二区免费| 人妻丰满熟妇av一区二区三区| 国产又黄又爽又无遮挡在线| 夜夜夜夜夜久久久久| 在线观看免费午夜福利视频| 美女cb高潮喷水在线观看 | 国产精品野战在线观看| 亚洲av中文字字幕乱码综合| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人不卡在线观看播放网| 好男人电影高清在线观看| 亚洲欧美日韩无卡精品| 亚洲一区二区三区色噜噜| 成年女人永久免费观看视频| 18禁观看日本| 变态另类成人亚洲欧美熟女| 波多野结衣巨乳人妻| 最近最新中文字幕大全电影3| 亚洲av第一区精品v没综合| 国产高潮美女av| 免费在线观看影片大全网站| 亚洲 欧美 日韩 在线 免费| 国产精品1区2区在线观看.| 欧美绝顶高潮抽搐喷水| www国产在线视频色| 免费在线观看亚洲国产| 别揉我奶头~嗯~啊~动态视频| 一本一本综合久久| 1024手机看黄色片| 亚洲欧美激情综合另类| 午夜激情欧美在线| 热99re8久久精品国产| 哪里可以看免费的av片| 在线观看免费视频日本深夜| 1000部很黄的大片| 99国产极品粉嫩在线观看| 一级a爱片免费观看的视频| 亚洲真实伦在线观看| 亚洲成人久久性| 老司机在亚洲福利影院| 每晚都被弄得嗷嗷叫到高潮| av国产免费在线观看| 午夜a级毛片| 久久中文看片网| or卡值多少钱| 69av精品久久久久久| 亚洲无线在线观看| 日韩三级视频一区二区三区| 免费看美女性在线毛片视频| 亚洲国产欧美一区二区综合| 熟女人妻精品中文字幕| 色综合欧美亚洲国产小说| 99久久成人亚洲精品观看| 国产精品一区二区免费欧美| 亚洲欧美精品综合久久99| 最近最新中文字幕大全电影3| 国产精品一区二区免费欧美| 欧美高清成人免费视频www| 亚洲国产精品合色在线| 亚洲中文字幕一区二区三区有码在线看 | 这个男人来自地球电影免费观看| 国产精品美女特级片免费视频播放器 | 精品免费久久久久久久清纯| 日日干狠狠操夜夜爽| 欧美性猛交黑人性爽| 国产精品爽爽va在线观看网站| 午夜视频精品福利| 老熟妇仑乱视频hdxx| 亚洲精品在线美女| 日韩精品青青久久久久久| 日韩成人在线观看一区二区三区| 国产精品av视频在线免费观看| 亚洲国产看品久久| 精华霜和精华液先用哪个| 在线永久观看黄色视频| 亚洲熟妇中文字幕五十中出| 哪里可以看免费的av片| 偷拍熟女少妇极品色| 一个人看视频在线观看www免费 | 午夜久久久久精精品| 黑人欧美特级aaaaaa片| 国产欧美日韩精品一区二区| 老司机深夜福利视频在线观看| 久久久久国产精品人妻aⅴ院| 欧美不卡视频在线免费观看| 国产高潮美女av| 国产欧美日韩一区二区三| 黄频高清免费视频| a级毛片在线看网站| 在线看三级毛片| 黄色丝袜av网址大全| 桃色一区二区三区在线观看| 婷婷精品国产亚洲av在线| 日本熟妇午夜| 手机成人av网站| 三级男女做爰猛烈吃奶摸视频| 国产成人精品久久二区二区免费| 欧美日本视频| 757午夜福利合集在线观看| 成人一区二区视频在线观看| 非洲黑人性xxxx精品又粗又长| 久久久久九九精品影院| 国产亚洲精品久久久com| 男人和女人高潮做爰伦理| 欧美黑人巨大hd| 日本 欧美在线| 午夜福利视频1000在线观看| 亚洲成av人片免费观看| 久久午夜亚洲精品久久| 午夜激情欧美在线| 成人亚洲精品av一区二区| 国内久久婷婷六月综合欲色啪| 精品国内亚洲2022精品成人| 国产亚洲精品av在线| 国产人伦9x9x在线观看| 天堂影院成人在线观看| 小蜜桃在线观看免费完整版高清| 757午夜福利合集在线观看| 国产免费男女视频| 亚洲精品在线美女| 国产亚洲精品久久久久久毛片| 黄色 视频免费看| 国产av一区在线观看免费| 在线观看美女被高潮喷水网站 | 亚洲熟妇熟女久久| 夜夜爽天天搞| 免费看日本二区| 看黄色毛片网站| 免费在线观看影片大全网站| 51午夜福利影视在线观看| 久久伊人香网站| 黄色成人免费大全| 国产精品av视频在线免费观看| 亚洲国产精品成人综合色| 中出人妻视频一区二区| 亚洲欧美精品综合一区二区三区| 国产高清激情床上av| 久久亚洲真实| www.熟女人妻精品国产| 欧美日韩黄片免| 国产91精品成人一区二区三区| 国产精品 国内视频| 日本在线视频免费播放| 麻豆成人av在线观看| 香蕉av资源在线| 精品久久久久久久末码| 制服人妻中文乱码| 老汉色∧v一级毛片| 亚洲国产精品sss在线观看| 国产一区二区三区在线臀色熟女| 亚洲五月天丁香| 色av中文字幕| 亚洲国产精品合色在线| 男女做爰动态图高潮gif福利片| 老鸭窝网址在线观看| www.精华液| 99国产精品一区二区三区| 久久国产乱子伦精品免费另类| 日韩欧美免费精品| 久久久国产成人精品二区| 欧美乱色亚洲激情| 日韩欧美一区二区三区在线观看| 高清毛片免费观看视频网站| 精品一区二区三区av网在线观看| 一进一出抽搐gif免费好疼| 床上黄色一级片| 中国美女看黄片| 首页视频小说图片口味搜索| 成年免费大片在线观看| 天天添夜夜摸| 黄色视频,在线免费观看| 久久亚洲真实| a级毛片在线看网站| 操出白浆在线播放| av女优亚洲男人天堂 | 曰老女人黄片| 精品国产超薄肉色丝袜足j| 在线免费观看的www视频| 好男人在线观看高清免费视频| 91老司机精品| 99久久久亚洲精品蜜臀av| 日韩欧美精品v在线| 热99在线观看视频| 国产亚洲av嫩草精品影院| 啦啦啦观看免费观看视频高清| 欧美日韩综合久久久久久 | 熟女电影av网| 白带黄色成豆腐渣| www.自偷自拍.com| 久久这里只有精品中国| 亚洲精品国产精品久久久不卡| 两个人视频免费观看高清| 久99久视频精品免费| 亚洲精品乱码久久久v下载方式 | 亚洲人成伊人成综合网2020| 午夜日韩欧美国产| 99re在线观看精品视频| 99riav亚洲国产免费| 亚洲 欧美一区二区三区| 久久久久亚洲av毛片大全| 噜噜噜噜噜久久久久久91| 久99久视频精品免费| 国产aⅴ精品一区二区三区波| 国产精品一区二区免费欧美| 欧美绝顶高潮抽搐喷水| 99国产精品一区二区三区| 免费在线观看影片大全网站| 巨乳人妻的诱惑在线观看| 婷婷六月久久综合丁香| 91av网站免费观看| 啪啪无遮挡十八禁网站| 夜夜躁狠狠躁天天躁| 国产三级中文精品| 俄罗斯特黄特色一大片| 久久久色成人| 夜夜躁狠狠躁天天躁| 日本五十路高清| 狠狠狠狠99中文字幕| 国产高清videossex| 悠悠久久av| 亚洲九九香蕉| av黄色大香蕉| 男人舔奶头视频| 国内毛片毛片毛片毛片毛片| 成人18禁在线播放| 精品久久久久久久毛片微露脸| 啪啪无遮挡十八禁网站| 中国美女看黄片| 亚洲熟妇熟女久久| 女警被强在线播放| 波多野结衣巨乳人妻| 波多野结衣高清无吗| 综合色av麻豆| 狂野欧美激情性xxxx| 欧美日韩综合久久久久久 | 国产精品一区二区精品视频观看| 欧美日本亚洲视频在线播放| 叶爱在线成人免费视频播放| 老司机福利观看| 在线a可以看的网站| 精品乱码久久久久久99久播| 欧美日韩福利视频一区二区| 国产午夜福利久久久久久| 亚洲七黄色美女视频| 在线国产一区二区在线| 日本精品一区二区三区蜜桃| 九色国产91popny在线| 午夜a级毛片| 国产69精品久久久久777片 | 美女扒开内裤让男人捅视频| 国产男靠女视频免费网站| 最近在线观看免费完整版| av天堂中文字幕网| 一个人观看的视频www高清免费观看 | 国产精品久久久久久亚洲av鲁大| 麻豆国产av国片精品| 日本成人三级电影网站| 亚洲欧美日韩卡通动漫| 成人永久免费在线观看视频| 亚洲性夜色夜夜综合| 窝窝影院91人妻| 免费电影在线观看免费观看| 久久午夜综合久久蜜桃| 国产主播在线观看一区二区| 成人亚洲精品av一区二区| 最好的美女福利视频网| 久久久久性生活片| 国产精品久久久久久亚洲av鲁大| 热99re8久久精品国产| 人妻夜夜爽99麻豆av| 成人无遮挡网站| 欧美不卡视频在线免费观看| 91在线精品国自产拍蜜月 | 亚洲成人精品中文字幕电影| e午夜精品久久久久久久| 国产爱豆传媒在线观看| 在线观看免费午夜福利视频| 久久久久久九九精品二区国产| 青草久久国产| 一区二区三区激情视频| 男人的好看免费观看在线视频| 网址你懂的国产日韩在线| 亚洲七黄色美女视频| 好看av亚洲va欧美ⅴa在| 色综合站精品国产| 小蜜桃在线观看免费完整版高清| 亚洲无线在线观看| 精品国内亚洲2022精品成人| 九色成人免费人妻av| 亚洲av成人精品一区久久| 久久精品国产亚洲av香蕉五月| 老司机午夜福利在线观看视频| 免费看a级黄色片| 无人区码免费观看不卡| 午夜福利在线在线| 亚洲欧美激情综合另类| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 精品人妻1区二区| 三级毛片av免费| ponron亚洲| 久9热在线精品视频| 999久久久精品免费观看国产| 性色avwww在线观看| 一个人看的www免费观看视频| 俄罗斯特黄特色一大片| 91老司机精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲色图av天堂| 亚洲国产精品999在线| 99riav亚洲国产免费| 成人18禁在线播放| 一区二区三区激情视频| or卡值多少钱| 日本在线视频免费播放| av天堂中文字幕网| 美女高潮喷水抽搐中文字幕| 国内揄拍国产精品人妻在线| 极品教师在线免费播放| 国产探花在线观看一区二区| 91麻豆精品激情在线观看国产| 欧美中文日本在线观看视频| 老汉色∧v一级毛片| 国产成人影院久久av| 一级毛片高清免费大全| 亚洲,欧美精品.| 中文字幕精品亚洲无线码一区| 黑人欧美特级aaaaaa片| 99国产极品粉嫩在线观看| 中亚洲国语对白在线视频| 人妻丰满熟妇av一区二区三区| 村上凉子中文字幕在线| 欧美在线一区亚洲| 成人欧美大片| 亚洲18禁久久av| 国产高清三级在线| 久久久水蜜桃国产精品网| 窝窝影院91人妻| 国产三级黄色录像| 中文字幕av在线有码专区| 亚洲国产欧美人成| 搡老熟女国产l中国老女人| 叶爱在线成人免费视频播放| 特级一级黄色大片| а√天堂www在线а√下载| 757午夜福利合集在线观看| 国产成人精品无人区| 白带黄色成豆腐渣| a在线观看视频网站| 在线十欧美十亚洲十日本专区| 19禁男女啪啪无遮挡网站| 中文在线观看免费www的网站| 视频区欧美日本亚洲| 精品久久蜜臀av无| 老熟妇仑乱视频hdxx| 国产精品女同一区二区软件 | 亚洲欧美日韩无卡精品| 禁无遮挡网站| 99国产精品一区二区三区| 波多野结衣高清作品| 国产av一区在线观看免费| 久久欧美精品欧美久久欧美| 亚洲色图 男人天堂 中文字幕| 国语自产精品视频在线第100页| av在线天堂中文字幕| 国产av麻豆久久久久久久| 亚洲激情在线av| 亚洲欧美精品综合久久99| 亚洲精品在线观看二区| 在线观看日韩欧美| 淫妇啪啪啪对白视频| 男女做爰动态图高潮gif福利片| 午夜精品久久久久久毛片777| 国产乱人伦免费视频| 国产私拍福利视频在线观看| 亚洲欧美激情综合另类| 在线观看66精品国产| 噜噜噜噜噜久久久久久91| 精品日产1卡2卡| 黄色 视频免费看| 日本 欧美在线| 99国产精品99久久久久| www.999成人在线观看| 久久久国产成人精品二区| 精品一区二区三区视频在线观看免费| 91av网站免费观看| 成人欧美大片| 日本与韩国留学比较| 久久性视频一级片| 亚洲自偷自拍图片 自拍| 十八禁网站免费在线| 亚洲国产中文字幕在线视频| 日本与韩国留学比较| 欧美乱妇无乱码| 一个人看的www免费观看视频| 国产精品女同一区二区软件 | 亚洲 欧美一区二区三区| 国产主播在线观看一区二区| 真人一进一出gif抽搐免费| 免费在线观看视频国产中文字幕亚洲| 国产真人三级小视频在线观看| 午夜成年电影在线免费观看| 在线十欧美十亚洲十日本专区| 久久久久久国产a免费观看| 精品日产1卡2卡| 国产精品综合久久久久久久免费| а√天堂www在线а√下载| 啦啦啦韩国在线观看视频| 巨乳人妻的诱惑在线观看| 男人和女人高潮做爰伦理| 超碰成人久久| 亚洲国产看品久久| aaaaa片日本免费| 国产精品久久久久久久电影 | a级毛片在线看网站| 久久久久久久午夜电影| 国产69精品久久久久777片 | 久久午夜亚洲精品久久| 婷婷亚洲欧美| 国产伦人伦偷精品视频| 少妇的逼水好多| 波多野结衣高清作品| 国产精品 欧美亚洲| 亚洲精品一区av在线观看| 精品熟女少妇八av免费久了| 国产精品久久久久久久电影 | 一夜夜www| 亚洲午夜理论影院| 国产成年人精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 神马国产精品三级电影在线观看| 1024手机看黄色片| 99视频精品全部免费 在线 | 欧美日韩瑟瑟在线播放| 亚洲国产欧洲综合997久久,| 亚洲国产精品久久男人天堂| 亚洲电影在线观看av| 一二三四社区在线视频社区8| 亚洲精品乱码久久久v下载方式 | 免费看光身美女| 男女床上黄色一级片免费看| svipshipincom国产片| 国产三级中文精品| 麻豆成人午夜福利视频| 精品国产亚洲在线| 夜夜爽天天搞| 757午夜福利合集在线观看| 制服丝袜大香蕉在线| 欧美日韩国产亚洲二区| 国产黄a三级三级三级人| 美女扒开内裤让男人捅视频| 亚洲美女黄片视频| 别揉我奶头~嗯~啊~动态视频| 免费av不卡在线播放| 国产熟女xx| 国产亚洲欧美在线一区二区| 免费观看精品视频网站| 黄色女人牲交| 老司机午夜十八禁免费视频| 不卡一级毛片| 母亲3免费完整高清在线观看| 亚洲成av人片在线播放无| 99精品欧美一区二区三区四区| 亚洲精品456在线播放app | 一进一出抽搐gif免费好疼| 欧美在线一区亚洲| 午夜日韩欧美国产| 十八禁网站免费在线| 中文字幕最新亚洲高清| 国产成人系列免费观看| 可以在线观看的亚洲视频| 国产99白浆流出| 国产精华一区二区三区| 婷婷亚洲欧美| 观看美女的网站| 成人欧美大片| 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 亚洲美女黄片视频| 非洲黑人性xxxx精品又粗又长| 三级国产精品欧美在线观看 | 三级国产精品欧美在线观看 | 黑人巨大精品欧美一区二区mp4| 免费大片18禁| 成人精品一区二区免费| 久久中文看片网| 国产精品永久免费网站| 午夜影院日韩av| 在线观看免费午夜福利视频| 久久香蕉精品热| 99国产精品一区二区蜜桃av| 欧美乱码精品一区二区三区| 亚洲18禁久久av| 国内久久婷婷六月综合欲色啪| 少妇人妻一区二区三区视频| 操出白浆在线播放| 熟女电影av网| 桃红色精品国产亚洲av| 国产精品1区2区在线观看.| 国内毛片毛片毛片毛片毛片| 亚洲国产欧美网| 久久久久久国产a免费观看| 99热精品在线国产| 男女床上黄色一级片免费看| 成人18禁在线播放| 亚洲美女视频黄频| 欧美日韩亚洲国产一区二区在线观看| 无遮挡黄片免费观看| 精品一区二区三区四区五区乱码| 国产精品综合久久久久久久免费| 亚洲成av人片免费观看| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 日韩精品中文字幕看吧| 国产高清有码在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 中文字幕av在线有码专区| 久久香蕉国产精品| 精品久久久久久久人妻蜜臀av| 成人特级av手机在线观看| 日韩欧美在线乱码| 噜噜噜噜噜久久久久久91| 看片在线看免费视频| 九色成人免费人妻av| 人妻丰满熟妇av一区二区三区| 日本精品一区二区三区蜜桃| 国产综合懂色| 婷婷六月久久综合丁香| 欧美另类亚洲清纯唯美| 美女免费视频网站| 一卡2卡三卡四卡精品乱码亚洲| 97人妻精品一区二区三区麻豆| 色综合欧美亚洲国产小说| 一区福利在线观看| 精品午夜福利视频在线观看一区| 村上凉子中文字幕在线| 99久国产av精品| 丁香欧美五月| 欧美另类亚洲清纯唯美| 亚洲片人在线观看| 日韩欧美在线二视频| 久久婷婷人人爽人人干人人爱| 天堂影院成人在线观看| 宅男免费午夜| 男女那种视频在线观看| 亚洲av成人av| 欧美最黄视频在线播放免费| 久久香蕉国产精品| 国产激情偷乱视频一区二区| 亚洲熟妇熟女久久| 99久久精品一区二区三区| 久久久久久国产a免费观看| 亚洲av中文字字幕乱码综合| 淫秽高清视频在线观看| 9191精品国产免费久久| 亚洲性夜色夜夜综合| 亚洲午夜理论影院| 观看美女的网站| 欧美激情久久久久久爽电影| 欧美激情在线99| 免费av不卡在线播放| 丝袜人妻中文字幕| 久久精品91蜜桃| 午夜精品久久久久久毛片777| 国产高清视频在线播放一区| 可以在线观看的亚洲视频| 欧美+亚洲+日韩+国产| 两人在一起打扑克的视频| 国产成人系列免费观看| 欧美中文综合在线视频| 久久精品国产清高在天天线| 国产激情欧美一区二区| 亚洲av第一区精品v没综合| 亚洲中文字幕日韩| 国产成人aa在线观看|