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

    基于多尺度建模的煉油化工過程報警根源分析

    2018-04-12 02:39:34胡瑾秋張來斌
    石油學(xué)報(石油加工) 2018年2期
    關(guān)鍵詞:分餾塔根源煉油

    胡瑾秋, 蔡 爽, 張來斌

    (中國石油大學(xué) 機(jī)械與儲運工程學(xué)院, 北京 102249)

    煉油化工生產(chǎn)過程從原料到產(chǎn)品的生產(chǎn)需要經(jīng)過許多工序和復(fù)雜的加工單元,通過多次反應(yīng)或分離完成。其流程復(fù)雜,工藝條件苛刻,常伴有高溫、高壓、低溫、真空、大流量、高轉(zhuǎn)速等極端條件,在煉油廠區(qū),若某單一設(shè)備或工藝過程出現(xiàn)故障,極易借助系統(tǒng)單元之間的相互依存、相互制約關(guān)系,觸發(fā)連鎖效應(yīng),由一種擾動或故障經(jīng)過一段時間傳播,從一個地域空間傳播至另一個更廣闊的地域空間,引起衍生事故,帶來難以估量的社會與環(huán)境風(fēng)險,引發(fā)一系列故障鏈直至引起事故或災(zāi)難。為了及時找出故障原因以有效抑制故障傳播,保證生產(chǎn)過程安全及產(chǎn)品質(zhì)量,國內(nèi)外學(xué)者對煉油化工過程報警根源分析方法進(jìn)行了一系列的研究。

    報警泛濫是工業(yè)過程報警系統(tǒng)中存在的主要問題之一。對過程變量報警根源進(jìn)行分析,可及時抑制異常擾動的進(jìn)一步傳播以減少報警泛濫現(xiàn)象的發(fā)生。為了辨識煉油化工過程的擾動傳播根源,基于過程知識建立過程拓?fù)淠P陀兄谶M(jìn)一步分析報警根源,如多層流模型(Multilevel flow model,MFM)[1],因果模型(Cause-effect model,CE)[2],符號有向圖(Signed directed graph,SDG)[3]等。Dahlstrand[1]采用多層流方法表達(dá)過程變量間的因果關(guān)系。Wan等[2]將CE模型引入到一個邏輯和系統(tǒng)化的報警系統(tǒng)設(shè)計方法中。符號有向圖是因果圖的一種,對過程變量間的因果關(guān)系和故障傳播路徑提供了簡單和圖形化描述,并且可以在缺乏詳細(xì)過程知識的條件下根據(jù)經(jīng)驗或基本定律來建立符號有向圖模型[4-7]。知識建模方法簡單實用,易于理解,但嚴(yán)重依賴于過程知識不包括定量信息;此外,模型容易受人的主觀判斷影響[8]。數(shù)據(jù)驅(qū)動方法如時滯關(guān)聯(lián)分析[9],信息理論方法(傳遞熵,直接傳遞熵)[10-11],貝葉斯網(wǎng)絡(luò)模型(Bayesian network,BN)[12]等,也可用于實現(xiàn)報警根源的辨識,以便從根本上抑制滋擾報警和報警泛濫現(xiàn)象的出現(xiàn)。關(guān)聯(lián)性并不一定代表因果性,但如果引入時間延遲,可找出兩時間序列間關(guān)聯(lián)系數(shù)的最大值,其對應(yīng)的時間關(guān)系可以近似地反映因果關(guān)系?;ハ嚓P(guān)函數(shù)(Cross-correlation function,CCF)[9]可用于估計兩報警序列間的時間延遲,通過兩序列間的時滯和關(guān)聯(lián)系數(shù),生成對應(yīng)的因果矩陣和關(guān)聯(lián)系數(shù)矩陣,結(jié)合一致性檢驗和拓?fù)浼僭O(shè),可以構(gòu)建因果網(wǎng)絡(luò)[13]。但互相關(guān)函數(shù)會產(chǎn)生大量冗余時滯干擾分析結(jié)果,有時由該數(shù)據(jù)驅(qū)動方法建立的擾動傳播路徑可能與實際過程不符。因此有必要結(jié)合過程知識檢驗傳播路徑的合理性。Bauer等[11]和Gao等[14]采用信息理論計算兩變量間的傳遞熵,從本質(zhì)上反映了信息傳遞所帶來的因果性。Duan等[10]提出直接傳遞熵概念來檢測兩變量間是否有一個直接的物質(zhì)流或信息流,以提高因果圖的準(zhǔn)確性。Gao等[15]提出基于解析結(jié)構(gòu)模型的信息融合方法,通過構(gòu)建連接矩陣和可達(dá)矩陣構(gòu)建因果網(wǎng)絡(luò),實現(xiàn)報警根源的辨識。概率符號有向圖[16]是一種基于實現(xiàn)最大系統(tǒng)可信度的報警信號挑選方法,該符號有向圖的建立使過程變量間的因果關(guān)系可視化。數(shù)據(jù)驅(qū)動方法可定量分析過程變量間存在的因果關(guān)系,但對于分析變量的挑選依然常常根據(jù)主觀經(jīng)驗,對于一些復(fù)雜的大型系統(tǒng),常涉及大量監(jiān)控變量,主觀挑選模型變量可能遺漏某些重要原因變量,以致無法推繹出引發(fā)報警的擾動根源。

    因此筆者采用多尺度建模的方法,通過對煉油化工生產(chǎn)系統(tǒng)的空間尺度進(jìn)行劃分,針對不同的風(fēng)險過程,分析不同尺度下的過程風(fēng)險表征指標(biāo),建立多尺度下的煉油化工過程報警根源分析模型,并對不同尺度下的模型進(jìn)行應(yīng)用對比分析,從而選擇合適的尺度進(jìn)行建模,保證煉油化工過程報警根源分析的準(zhǔn)確性。

    1 時滯分析方法基本原理

    時滯分析方法可分析不同過程變量間的時滯和關(guān)聯(lián)關(guān)系,從而推斷他們固有的因果關(guān)系。

    時滯分析方法采用互相關(guān)函數(shù)[9]計算過程變量間考慮時延性的關(guān)聯(lián)函數(shù)值。對于兩過程變量xi和xj的離散時間序列數(shù)據(jù),定義互相關(guān)函數(shù)φxi,xj(k)為

    (1)

    n為過程變量個數(shù),過程變量xi和xj的時間序列長度為N,其均值為μxi和μxj,標(biāo)準(zhǔn)差為σxi和σxj,xig表示標(biāo)準(zhǔn)化后的變量xi在第g個時刻的變量值,φxi,xj(k)為具有時間意義的k的單值函數(shù)。

    定義φxi,xj(k)的最大值φmax和最小值φmin為

    (2)

    記φmax對應(yīng)的k為kmax,φmin對應(yīng)的k為kmin,xi、xj間的時滯λi,j表示為

    (3)

    時滯λi,j用于確定傳播方向。若λi,j>0,實際傳播方向為xi→xj;反之,傳播方向為xj→xi。

    兩時間序列關(guān)聯(lián)系數(shù)ρi,j記為

    (4)

    因關(guān)聯(lián)系數(shù)由統(tǒng)計方法計算得到,每2個時間序列可得一確定值,但若關(guān)聯(lián)系數(shù)過小,考慮兩變量間的關(guān)聯(lián)性將沒有意義。因此有必要設(shè)置合適的閾值對關(guān)聯(lián)系數(shù)的顯著性水平進(jìn)行檢驗。筆者采用如下的假設(shè)檢驗方法,通過式(5)計算閾值(ρth(N))[9],若|ρi,j|<ρth(N),表明關(guān)聯(lián)系數(shù)未通過顯著性檢驗,此時考慮兩變量間的關(guān)聯(lián)性將沒有意義。ρi,j的正負(fù)反映了變量間的關(guān)聯(lián)關(guān)系。若ρi,j>0,兩變量具有正相關(guān)關(guān)系,即xi的增加可能引起xj的增加。反之,若ρi,j<0, 兩變量具有負(fù)相關(guān)關(guān)系。

    ρth(N)=1.85N-0.41+2.37N-0.53

    (5)

    然而,若φmax和φmin的大小接近,兩變量的因果關(guān)系將變得模糊,方向性難以確定。為了辨識φmax和φmin的差異性,定義方向性指數(shù)(ψ)如式(6):

    (6)

    通過式(7)計算方向性指數(shù)閾值(σψ(N))[9],若ψ≥σψ(N),表明兩變量具有顯著的因果關(guān)系。

    σψ(N)=0.46N-0.16

    (7)

    2 煉油化工過程報警根源分析方法

    報警系統(tǒng)作為復(fù)雜煉油化工過程的重要組成部分,可有效監(jiān)控過程的運行狀態(tài)。報警系統(tǒng)設(shè)置是否合理直接影響生產(chǎn)過程的運行。當(dāng)過程設(shè)備異常運行時,將以聲光形式觸發(fā)報警。收到報警后,操作者將結(jié)合過程知識調(diào)查故障原因并采取必要措施阻止異常的進(jìn)一步惡化,使設(shè)備恢復(fù)到正常運行范圍。如今,工業(yè)過程報警系統(tǒng)的實際問題是報警數(shù)過多。實際上,操作者每1 h可能面臨數(shù)以百計的報警,大量的報警干擾操作者的判斷,使他們浪費大量時間處理次要報警而忽視關(guān)鍵報警。因此,需要一個合理的策略進(jìn)行報警根源分析,及時辨識擾動根源,避免異常風(fēng)險的進(jìn)一步發(fā)展。

    提出一種煉油化工過程報警根源分析方法,通過對煉油化工生產(chǎn)系統(tǒng)的空間尺度進(jìn)行劃分,針對不同的風(fēng)險過程,分析不同尺度下的過程風(fēng)險表征指標(biāo),從而通過時滯分析并結(jié)合過程知識構(gòu)建擾動傳播因果圖,最后提出一種報警根源搜索方法辨識最可能的擾動傳播路徑以找出報警的根本原因。

    2.1 煉油化工生產(chǎn)系統(tǒng)的空間尺度劃分

    煉油廠區(qū)包含多套生產(chǎn)裝置,如煉油廠主要由煉油過程和輔助設(shè)備兩大部分組成,經(jīng)過多個物理及化學(xué)的煉油過程將原油煉制成各種石油產(chǎn)品。

    煉油生產(chǎn)裝置按生產(chǎn)目的可分為原油分離裝置、重油輕質(zhì)化裝置、油品改質(zhì)及油品精制裝置、油品調(diào)合裝置、氣體加工裝置、制氫裝置等。由于生產(chǎn)方案不同,煉油廠中所包含的煉油過程的種類和規(guī)模大小各不相同,或者說復(fù)雜程度都有所不同。一般來說,規(guī)模大的煉油廠其復(fù)雜程度更高。此外,輔助設(shè)施是維持煉油廠正常運轉(zhuǎn)和生產(chǎn)所必需的,主要的輔助設(shè)施有:供電系統(tǒng)、供水系統(tǒng)、供水蒸氣系統(tǒng)、原油和產(chǎn)品儲運系統(tǒng)、三廢處理系統(tǒng)等。這里的研究不考慮輔助設(shè)施。

    煉油化工生產(chǎn)系統(tǒng)的空間尺度可定義為煉油化工生產(chǎn)系統(tǒng)空間范圍的大小。這里的空間范圍是指在一定尺度的系統(tǒng)中進(jìn)行風(fēng)險分析的范圍。根據(jù)煉油化工生產(chǎn)系統(tǒng)的規(guī)模大小,將空間尺度分為大尺度、中尺度和小尺度3種類型。

    其中,大尺度主要包括生產(chǎn)裝置級及比其空間范圍更大的煉油化工生產(chǎn)系統(tǒng)。其研究區(qū)域的范圍小至單個生產(chǎn)裝置(如常減壓裝置、催化裂化裝置等),大至整個煉油廠區(qū)。

    中尺度則主要指單元級。其研究區(qū)域的范圍為生產(chǎn)裝置中的某一單元。

    傳統(tǒng)分析中,單元的定義為整體中自為一組或自成系統(tǒng)的獨立單位。在安全生產(chǎn)風(fēng)險管理中,為了保證風(fēng)險辨識的科學(xué)性、系統(tǒng)性、合理性和可操作性而將研究對象整體按照某種規(guī)則劃分為相對獨立的部分。因此,這里的單元是指將某個生產(chǎn)裝置按照其生產(chǎn)功能劃分的各個相對獨立的部分。如催化裂化裝置可劃分為反應(yīng)再生單元、分餾單元、吸收穩(wěn)定單元等。

    小尺度主要是指設(shè)備級,其研究區(qū)域的范圍為某一過程設(shè)備,這里的設(shè)備包括各種塔器、反應(yīng)器、換熱器、各類儲罐等,如提升管反應(yīng)器、分餾塔、回?zé)捰凸薜取?/p>

    具體的尺度劃分如表1所示。

    表1 煉油化工生產(chǎn)系統(tǒng)空間尺度劃分Table 1 Spatial scale division of refinery production system

    2.2 多尺度風(fēng)險表征指標(biāo)辨識

    隨著煉油化工過程工藝、技術(shù)的提高和裝置大型化發(fā)展,保證過程的安全、平穩(wěn)、長周期運行成為煉油化工企業(yè)日常管理的重點,也是實現(xiàn)煉油企業(yè)效益最大化的有效途徑。因此,為避免發(fā)生各類非計劃停車、泄漏、火災(zāi)、爆炸等事故,需對煉油化工過程風(fēng)險及其風(fēng)險表征指標(biāo)進(jìn)行辨識研究,以便準(zhǔn)確、及時的對煉油化工過程風(fēng)險進(jìn)行狀態(tài)監(jiān)測及報警根源診斷,有效抑制異常擾動的進(jìn)一步傳播。

    通過統(tǒng)計現(xiàn)場煉油化工過程采集的監(jiān)控變量數(shù)據(jù),可根據(jù)某一風(fēng)險發(fā)展過程中一段時間內(nèi)各變量的報警數(shù)量及其擾動變化趨勢,辨識不同尺度下與風(fēng)險發(fā)展密切相關(guān)的監(jiān)控變量作為過程風(fēng)險表征指標(biāo)。這里的多尺度風(fēng)險表征指標(biāo)辨識方法主要包括如下兩個步驟:

    (1)統(tǒng)計某煉油化工企業(yè)近兩年事故記錄中每次某一風(fēng)險發(fā)展過程前后的一段時間內(nèi)不同尺度下各監(jiān)控變量的報警數(shù)量,因報警數(shù)的多少反映了風(fēng)險發(fā)生時對該變量的影響大小,因此,這里根據(jù)專家經(jīng)驗選取多次統(tǒng)計數(shù)據(jù)的平均報警數(shù)多于5個的變量作為待選的風(fēng)險表征指標(biāo)。

    (2)因受風(fēng)險影響較大的變量必然會產(chǎn)生明顯的擾動,可根據(jù)其擾動變化趨勢的大小進(jìn)一步確定過程風(fēng)險表征指標(biāo)。這里采用最小二乘線性擬合方法對上述不同尺度下的待選風(fēng)險表征指標(biāo)在每次風(fēng)險發(fā)展過程中一段時間內(nèi)的數(shù)據(jù)進(jìn)行擬合,并以所求斜率的平均值作為該變量的擾動變化率,其絕對值的大小可近似反映變量的擾動變化趨勢,若擾動變化率的絕對值過小,說明該變量并未產(chǎn)生大的波動,即無法較好的表征風(fēng)險發(fā)展過程。因此若擾動變化率的絕對值小于某一閾值(根據(jù)專家經(jīng)驗設(shè)置)時,不使用相應(yīng)的變量作為沖塔風(fēng)險表征指標(biāo)。

    其中,擾動變化率定義如下:

    定義1:對于某一過程變量xi,考慮以時刻κ為中心,選擇時間間隔為[κ-m,κ+m] (時間序列長度為2m+1)的變量xi的時間序列進(jìn)行最小二乘線性擬合如式(8),所求斜率ai的大小作為變量xi的擾動變化率。

    令xig=aitg+bi,g=κ-m,…,κ,…,κ+m,i= 1,…,n,n為過程變量個數(shù),xig為標(biāo)準(zhǔn)化后的變量xi在第g個時刻的變量值,最小二乘線性擬合公式如式(8):

    (8)

    對變量xi的標(biāo)準(zhǔn)化公式如式(9):

    (9)

    通過上述步驟,即可確定不同尺度下的過程風(fēng)險表征指標(biāo)。

    2.3 因果圖的建立

    采用時滯分析方法分析不同過程變量間的時滯和關(guān)聯(lián)關(guān)系,可推斷他們固有的因果關(guān)系,該方法適用于運行相對平穩(wěn)的化工過程[9]。時滯分析的基本假設(shè)是最先引起擾動產(chǎn)生的變量是將風(fēng)險引入整個過程的根原因。擾動變量不同于異常報警變量。最早出現(xiàn)的報警變量是在故障影響下最先超出預(yù)設(shè)操作范圍的變量,但它不一定是最早的擾動,也不一定是引發(fā)報警的根本原因。因此,當(dāng)報警出現(xiàn)時,有必要通過分析擾動傳播路徑辨識報警根源從而消除報警。

    變量間的時滯信息和因果關(guān)系可以通過時滯分析方法來獲得,從而構(gòu)造出兩個上三角矩陣。時滯矩陣反映了受擾動影響的兩變量間在時間發(fā)展上的先后關(guān)系。其中元素的正負(fù)反映變量間的因果關(guān)系。最初受到擾動的變量是后續(xù)變量發(fā)生擾動的原因。關(guān)聯(lián)系數(shù)矩陣中元素的正負(fù)反映了變量間的正負(fù)相關(guān)關(guān)系。時滯矩陣L和關(guān)聯(lián)系數(shù)矩陣P有下面的形式:

    (10)

    (11)

    兩過程變量間時滯λi,j和關(guān)聯(lián)系數(shù)記ρi,j分別由式(3)及式(4)求得,并根據(jù)式(6)計算兩變量的方向性指數(shù)ψi,j,根據(jù)式(5)計算關(guān)聯(lián)顯著性閾值ρth(N),通過式(7)計算方向性指數(shù)閾值σψ(N),若|ρi,j|<ρth(N)或ψi,j<σψ(N),ρi,j和λi,j在式(10)和式(11)中將設(shè)為0。

    基于這些特征,通過時滯分析建立過程擾動傳播因果圖,如圖1所示。因果圖由有向弧和代表過程變量xi(i=1,2,…6) 的節(jié)點組成。若λi,j>0,有向弧由xi指向xj,即由上級原因變量指向下級結(jié)果變量,反之,則傳播方向相反。虛線弧表示兩變量負(fù)相關(guān),即ρi,j<0,實線弧表示兩變量正相關(guān),即ρi,j>0。

    圖1 擾動傳播因果圖Fig.1 Disturbance propagation causality diagramxi(i=1,2,…6) represents the ith process variable node;The dotted line indicates a negative correlation between two variables;The solid line indicates a positive correlation between two variables

    但互相關(guān)函數(shù)會產(chǎn)生大量冗余時滯干擾分析結(jié)果。有時由該數(shù)據(jù)驅(qū)動方法建立的擾動傳播路徑可能與實際過程不符。因此有必要結(jié)合過程知識檢驗傳播路徑的合理性,對因果圖進(jìn)行修正。

    2.4 報警根源搜索方法

    當(dāng)一個報警出現(xiàn)時,報警可能不僅受到相鄰變量的影響,也受到傳播路徑中其他變量的影響。若變量值超出預(yù)設(shè)控制限,將觸發(fā)相應(yīng)報警。但是該報警不一定是最早產(chǎn)生的擾動,也不一定是系統(tǒng)異常的根本原因。在過程中的其它變量受到擾動前發(fā)生擾動的變量,受到故障的影響最為直接,故障通過該變量在過程中不斷傳播進(jìn)而影響其他變量導(dǎo)致報警的產(chǎn)生。因此找出最初擾動、確定報警的根本原因,操作人員可迅速采取有效措施處理報警。這里提出一種報警根源搜索方法,包括如下4個步驟:

    步驟1:相同的報警可能由不同傳播路徑中產(chǎn)生的擾動引起。為了辨識報警的根本原因,當(dāng)一個變量xj發(fā)生報警時,將其作為下級變量,根據(jù)2.3節(jié)所提方法建立因果圖,據(jù)此搜索與該變量直接相連的所有上級變量xi(i=1,2,…,I,I為上級變量個數(shù)),例如,圖1中的x3報警,可搜索到其2個上級變量x2和x5。

    步驟2:若變量xj在κ時刻發(fā)生報警,通過式(8)和式(9)計算xj各相關(guān)上級變量xi的擾動變化率ai,其值大小可近似反映下級變量受各相關(guān)上級變量擾動的影響大小。

    步驟3:比較各相關(guān)上級變量的擾動變化率ai的絕對值大小,將ai絕對值最大值對應(yīng)的上級變量作為其上級原因變量。上級變量的擾動變化率ai的絕對值越大,其對下級變量的影響越大,因此,通過ai絕對值的大小可以比較下級變量受各相關(guān)上級變量擾動的影響大小,將ai絕對值的最大值對應(yīng)的上級變量作為其上級原因變量。若所求擾動變化率的絕對值過小,說明該變量的時間序列數(shù)據(jù)趨于平穩(wěn),并未產(chǎn)生大的擾動傳播至下級變量,因此,根據(jù)專家經(jīng)驗及歷史數(shù)據(jù)統(tǒng)計,若ai的絕對值小于閾值C,將不考慮該上級變量。

    步驟4:重復(fù)步驟2和步驟3,將搜索到的最后一個上級變量作為報警的根本原因。這里有必要確定ai所對應(yīng)的上級變量變化趨勢是否與變量間的正負(fù)相關(guān)關(guān)系一致。也就是說,若兩變量負(fù)相關(guān),一個變量的變化率為正,另一個應(yīng)該為負(fù);若兩變量正相關(guān),一個變量的變化率為正,另一個也應(yīng)該為正。因此,若對應(yīng)于ai的上級變量變化趨勢與他們之間的正負(fù)相關(guān)關(guān)系一致,擾動可能由該上級變量向下級變量傳播;若不一致,根據(jù)ai值的大小依次考慮其他上級變量是否滿足要求,并重復(fù)步驟2和步驟3,將搜索到的最后一個上級變量作為報警的根本原因。

    3 煉油化工過程報警根源分析方法步驟

    煉油化工過程報警根源分析方法流程圖如圖2所示,具體方法步驟如下:

    步驟1:對煉油化工生產(chǎn)系統(tǒng)的空間尺度進(jìn)行劃分,劃分標(biāo)準(zhǔn)見2.1節(jié);

    步驟2:采用2.2節(jié)所提方法辨識不同尺度下與風(fēng)險發(fā)展密切相關(guān)的監(jiān)控變量作為過程風(fēng)險表征指標(biāo);

    步驟3:采用時滯分析方法建立過程風(fēng)險擾動傳播因果圖,并結(jié)合過程知識檢驗傳播路徑的合理性,對因果圖進(jìn)行修正,具體方法見2.3節(jié);

    步驟4:當(dāng)某一過程變量發(fā)生報警時,采用所提報警根源搜索方法辨識報警的根本原因,具體方法見2.4節(jié)。

    圖2 煉油化工過程報警根源分析方法流程圖Fig.2 Diagram of alarm root cause analysis in refinery process

    4 案例分析

    催化裂化裝置是將重質(zhì)油轉(zhuǎn)化為輕質(zhì)油的關(guān)鍵裝置,其運行狀態(tài)不僅關(guān)系到整個廠區(qū)的安全,同時也決定了產(chǎn)品的收率。由于催化裂化裝置是煉油廠最主要和最關(guān)鍵的裝置,保持其安全平穩(wěn)長周期運行是提高催化裂化裝置的經(jīng)濟(jì)效益、避免各類非計劃停車、泄漏、火災(zāi)、爆炸等事故的關(guān)鍵所在。從催化裂化裝置運行情況來看,隨著摻渣比的提高,原料重質(zhì)化,劣質(zhì)化,反應(yīng)再生系統(tǒng)、分餾系統(tǒng)結(jié)焦、分餾塔沖塔、催化劑跑損、膨脹節(jié)開裂、油漿泵磨損、儀表故障以及人為操作事故等都是導(dǎo)致催化裂化裝置無法安全平穩(wěn)運行的主要因素。分餾單元是催化裂化生產(chǎn)裝置的一部分,分餾過程運行相對平穩(wěn),是一個非常復(fù)雜的物理變化過程。自反應(yīng)器頂部來的高溫油氣進(jìn)入分餾塔人字擋板底部,與人字擋板頂部返回的的循環(huán)油漿逆流接觸,油氣自下而上被冷卻洗滌,經(jīng)分餾后得到的富氣、粗汽油、輕循環(huán)油、重循環(huán)油及油漿。沖塔故障是分餾過程中經(jīng)常發(fā)生的現(xiàn)象之一。當(dāng)汽液相負(fù)荷過大時,氣體通過塔板的壓降增大,會使降液管中液面高度增加,液相負(fù)荷增加時,出口堰上液面高度增加,當(dāng)液體充滿整個降液管時,上下塔板液體連成一片,分餾完全破壞,導(dǎo)致沖塔的發(fā)生,嚴(yán)重影響產(chǎn)品質(zhì)量。造成分餾塔沖塔的原因有很多,如塔盤掀翻或損壞;塔盤有臟物堵塞浮閥或升氣孔;有臟物沉淀在塔板上,使溢流堰過低,液層過薄;塔盤結(jié)鹽;換熱器故障,取熱不能滿足,導(dǎo)致回流溫度高引起沖塔;各種機(jī)泵故障,如一中回流泵、回?zé)捰捅?、頂循泵等故障,機(jī)泵故障由帶水會不上量(即泵抽空)及各種設(shè)備自身問題引起,包括腐蝕、磨損、機(jī)封泄漏、入口有雜物損壞葉輪、電機(jī)故障等。為了保證產(chǎn)品質(zhì)量,避免沖塔故障的發(fā)生,這里對分餾塔沖塔風(fēng)險過程報警根源分析進(jìn)行多尺度建模研究,并對不同尺度下的模型進(jìn)行應(yīng)用對比分析,從而選擇合適的尺度進(jìn)行建模以保證過程風(fēng)險報警根源分析的準(zhǔn)確性。

    4.1 尺度劃分及多尺度沖塔風(fēng)險表征指標(biāo)辨識

    就分餾塔沖塔風(fēng)險而言,可主要從分餾塔本身及整個分餾單元兩個尺度進(jìn)行考慮,即設(shè)備級小尺度和單元級中尺度兩個尺度。

    分餾單元本身屬中尺度研究范圍,其研究區(qū)域為催化裂化生產(chǎn)裝置中的分餾單元如圖3所示,主要包括分餾塔、分餾塔頂油氣分離器、回?zé)捰凸蕖⒃嫌途彌_罐、輕柴油汽提塔及油漿外甩器。對分餾單元的沖塔風(fēng)險表征指標(biāo)進(jìn)行辨識,某石油化工企業(yè)催化裂化分餾單元的主要監(jiān)控變量如表2所示。

    表2 分餾單元主要監(jiān)控變量Table 2 Main monitoring variables of fractional distillation unit

    通過統(tǒng)計某石油化工企業(yè)現(xiàn)場催化裂化過程采集的監(jiān)控變量數(shù)據(jù)(采樣間隔為5 s),可根據(jù)分餾塔沖塔風(fēng)險發(fā)展過程中一段時間內(nèi)各變量的報警數(shù)量及其擾動變化趨勢,辨識不同尺度下與風(fēng)險發(fā)展密切相關(guān)的監(jiān)控變量,作為過程風(fēng)險表征指標(biāo)。

    統(tǒng)計該企業(yè)近2年事故記錄中每次分餾塔風(fēng)險發(fā)展過程前、后的10 min內(nèi)不同尺度下各監(jiān)控變量的報警數(shù)量,如表3所示,因報警數(shù)的多少反映了風(fēng)險發(fā)生時對該變量的影響大小,因此,這里根據(jù)專家經(jīng)驗選取多次統(tǒng)計數(shù)據(jù)的平均報警數(shù)多于5個的變量作為待選的風(fēng)險表征指標(biāo),如表3所示。

    圖3 催化裂化分餾單元Fig.3 Fractional distillation unit

    表3 待選風(fēng)險表征指標(biāo)Table 3 Possible risk indicators

    因受風(fēng)險影響較大的變量必然會產(chǎn)生明顯的擾動,可根據(jù)其擾動變化趨勢的大小進(jìn)一步確定過程風(fēng)險表征指標(biāo)。這里采用最小二乘線性擬合方法對上述不同尺度下的待選風(fēng)險表征指標(biāo)在每次風(fēng)險發(fā)展過程中前、后10 min內(nèi)(采樣間隔為5 s)的數(shù)據(jù)進(jìn)行擬合,并以所求斜率絕對值的平均值作為該變量的擾動變化率(計算如式(8)),如表3所示,其大小可近似反映變量的擾動變化趨勢,若擾動變化率過小,說明該變量并未產(chǎn)生大的波動,即無法較好的表征風(fēng)險發(fā)展過程。因此若擾動變化率小于某一閾值(這里根據(jù)專家經(jīng)驗設(shè)為1.5×10-3)時,不使用相應(yīng)的變量作為沖塔風(fēng)險表征指標(biāo)。

    通過上述步驟,即可確定單元級中尺度下的沖塔風(fēng)險表征指標(biāo),如表4所示。

    表4 小尺度及中尺度沖塔風(fēng)險表征指標(biāo)Table 4 Risk indicators of small scale and mesoscale

    分餾塔本身屬小尺度研究范圍,若僅從小尺度考慮時,沖塔風(fēng)險表征指標(biāo)包括表4中的分餾塔頂壓力、分餾塔塔底液位、分餾塔塔頂溫度、分餾塔塔底攪拌蒸汽流量、分餾塔底溫度、人字擋板上溫度、重柴油側(cè)線餾出口溫度及輕柴油側(cè)線餾出口溫度。

    4.2 分餾塔沖塔風(fēng)險因果圖的建立

    表4列出了小尺度及中尺度的沖塔風(fēng)險表征指標(biāo),基于互相關(guān)函數(shù)計算并結(jié)合過程知識建立沖塔風(fēng)險因果圖。

    4.2.1中尺度沖塔風(fēng)險因果圖的建立

    首先,通過式(1)計算過程風(fēng)險表征變量間考慮時延性的關(guān)聯(lián)函數(shù)值,并通過式(2)和式(3)得到過程變量間時滯以確定風(fēng)險傳播方向。通過式(4)可確定兩變量間的時延關(guān)聯(lián)系數(shù),如果關(guān)聯(lián)系數(shù)沒有通過顯著性檢驗式(5),兩變量間的關(guān)聯(lián)性將沒有意義。此外,通過式(6)式(7)進(jìn)一步確定兩變量因果關(guān)系的顯著性,從而構(gòu)造出時滯矩陣L和關(guān)聯(lián)系數(shù)矩陣P:

    過程變量個數(shù)n為18,時間序列長度N為300?;跁r滯矩陣L和關(guān)聯(lián)系數(shù)矩陣P,建立沖塔風(fēng)險過程擾動傳播因果圖,并結(jié)合過程知識檢驗傳播路徑的合理性,對因果圖進(jìn)行修正,結(jié)果如圖4 所示。

    圖4 中尺度沖塔風(fēng)險過程擾動傳播因果圖Fig.4 Cause and effect diagram of risk disturbancepropagation in mesoscalexi(i=1,2,…18) represents the ith process variable node corresponding to the risk indicators of mesoscale in Table 4 respectively The dotted lines indicate a negative correlation between two variables;The solid lines indicate a positive correlation between two variables

    4.2.2小尺度沖塔風(fēng)險因果圖的建立

    根據(jù)設(shè)備級小尺度的沖塔風(fēng)險表征指標(biāo),如表4 所示,同理,基于互相關(guān)函數(shù)計算并結(jié)合過程知識建立小尺度沖塔風(fēng)險過程擾動傳播因果圖,如圖5所示。

    圖5 小尺度沖塔風(fēng)險過程擾動傳播因果圖Fig.5 Cause and effect diagram of riskdisturbance propagation in small scalexi represents the ith process variable node corresponding to the risk indicators of small scale in Table 4 respectively The dotted lines indicate a negative correlation between two variables;The solid lines indicate a positive correlation between two variables

    4.3 分餾塔沖塔風(fēng)險報警根源搜索

    以發(fā)生在某石油化工企業(yè)的分餾塔沖塔事件為例,采用所提報警根源搜索方法進(jìn)行分析。2016年7月24日,某石油化工企業(yè)催化裂化分餾塔發(fā)生沖塔,經(jīng)現(xiàn)場人員調(diào)查分析,沖塔發(fā)生原因為一中回流泵故障,引起一中回流量過低,使得分餾塔塔頂溫度升高,重柴油側(cè)線餾出口溫度升高,從而引起塔底液位發(fā)生低報警。

    4.3.1中尺度沖塔風(fēng)險報警根源搜索

    下面首先基于單元級中尺度的沖塔風(fēng)險過程擾動傳播因果圖,進(jìn)行報警根源搜索分析,如圖6所示。

    圖6 中尺度沖塔風(fēng)險過程報警根源搜索圖Fig.6 The searching map for alarm root cause ofthe risk process in mesoscalexi(i=1,2,…18) represents the ith process variable node corresponding to the risk indicators of mesoscale in Table 4 respectively The dotted lines indicate a negative correlation between two variables;The solid lines indicate a positive correlation between two variables

    步驟1:因分餾塔塔底液位(x7)發(fā)生報警,所以將其作為下級變量,根據(jù)因果圖搜索與其直接相連的所有上級變量包括油漿下返塔流量(x6)、油漿上返塔流量(x3)、重柴油側(cè)線餾出口溫度(x5)和分餾塔塔頂溫度(x16)。

    步驟2:根據(jù)式(8)可得這些上級變量的擾動變化率如表5所示。

    步驟3:通過比較各相關(guān)上級變量的擾動變化率

    的絕對值,其最大值為3.5×10-3,對應(yīng)x16,大于閾值C(這里根據(jù)專家經(jīng)驗設(shè)為1.5×10-3),將上級變量x16作為其可能的上級原因變量。

    步驟4:進(jìn)一步確定該上級變量x16的變化趨勢是否與兩變量x7和x16間的正負(fù)相關(guān)關(guān)系一致。從圖中可辨識出x7與x16負(fù)相關(guān),因分餾塔塔底液位(x7)降低發(fā)生低報警,且所求x16的變化率為正,即分餾塔塔頂溫度升高,與兩變量間的正負(fù)相關(guān)關(guān)系一致,因此將上級變量分餾塔塔頂溫度(x16)作為上級原因變量。繼續(xù)搜索x16的上級變量包括一中回流量(x9)、分餾塔頂壓力(x14)、分餾塔頂循環(huán)油流量(x12)和分餾塔塔底攪拌蒸汽流量(x4),根據(jù)式(8)可得這些上級變量的擾動變化率如表5所示,通過比較各相關(guān)上級變量的擾動變化率的絕對值,其最大值為3.3×10-3,對應(yīng)x9,大于閾值C=1.5×10-3,因x16與x9負(fù)相關(guān),所求x9的變化率為負(fù),x16的變化率為正,與變量間的正負(fù)相關(guān)關(guān)系一致,因此將上級變量一中回流量(x9)作為上級原因變量,因x9無上級變量,所以一中回流量過低即為報警的根本原因,而一中回流量過低是由一中回流泵故障引起,該結(jié)論與現(xiàn)場人員調(diào)查分析所得的沖塔發(fā)生原因一中回流泵故障一致。

    表5 各上級變量的擾動變化率Table 5 The disturbance rate of each parent variable

    4.3.2小尺度沖塔風(fēng)險報警根源搜索

    下面基于設(shè)備級小尺度的沖塔風(fēng)險過程擾動傳播因果圖,進(jìn)行報警根源搜索分析,如圖7所示。

    步驟1:因分餾塔塔底液位(x7)發(fā)生報警,所以將其作為下級變量,根據(jù)因果圖搜索與其直接相連的所有上級變量包括重柴油側(cè)線餾出口溫度(x5)和分餾塔塔頂溫度(x16)。

    步驟2:根據(jù)式(8)可得這些上級變量的擾動變化率如表6所示:

    步驟3:通過比較各相關(guān)上級變量的擾動變化率的絕對值,其最大值為3.5×10-3,對應(yīng)x16,大于閾值C=1.5×10-3,將上級變量x16作為其可能的上級原因變量。

    圖7 小尺度沖塔風(fēng)險過程擾動傳播因果圖Fig.7 The searching map for alarm root cause ofthe risk process in small scalexi represents the ith process variable node corresponding to the risk indicators of small scale in Table 4 respectively The dotted lines indicate a negative correlation between two variables;The solid lines indicate a positive correlation between two variables

    表6 各上級變量的擾動變化率Table 6 The disturbance rate of each parent variable

    步驟4:從圖中可辨識出x7與x16負(fù)相關(guān),因分餾塔塔底液位(x7)降低發(fā)生低報警,且所求x16的變化率為正,即分餾塔塔頂溫度升高,與變量間的正負(fù)相關(guān)關(guān)系一致,因此將上級變量分餾塔塔頂溫度(x16)作為上級原因變量,繼續(xù)搜索x16的上級變量包括分餾塔頂壓力(x14)和分餾塔塔底攪拌蒸汽流量(x4),根據(jù)式(8)可得這些上級變量的擾動變化率如表6所示,通過比較各相關(guān)上級變量的擾動變化率的絕對值,其最大值為1.4,小于閾值C=1.5×10-3,因此該最大值對應(yīng)的上級變量分餾塔塔底攪拌蒸汽流量(x4)不再被考慮,搜索終止,所以分餾塔塔頂溫度(x16)過高即為報警的根本原因,該結(jié)論與現(xiàn)場人員調(diào)查分析所得的沖塔發(fā)生原因一中回流泵故障并不一致,即從設(shè)備級小尺度上進(jìn)行分析并未辨識出引發(fā)報警的根本原因。

    5 結(jié) 論

    通過對煉油化工生產(chǎn)系統(tǒng)的空間尺度進(jìn)行劃分,采用多尺度建模的方法,針對分餾塔沖塔風(fēng)險過程,分析其不同尺度下的過程風(fēng)險表征指標(biāo),建立多尺度下的沖塔風(fēng)險報警根源分析模型,并對不同尺度下的模型進(jìn)行應(yīng)用對比分析,從而選擇合適的尺度進(jìn)行建模,保證煉油化工過程報警根源分析的準(zhǔn)確性。

    對于過程風(fēng)險報警根源分析而言,需對整個風(fēng)險單元進(jìn)行分析,僅考慮與風(fēng)險直接相關(guān)的設(shè)備可能遺漏某些重要原因變量,以致無法推繹出引發(fā)報警的擾動根源。

    對于煉油化工過程風(fēng)險狀態(tài)監(jiān)測及故障診斷而言,選擇合適的尺度進(jìn)行建模分析,同樣有助于保證過程風(fēng)險狀態(tài)監(jiān)測及故障診斷的準(zhǔn)確性。

    對于存在反饋控制的復(fù)雜過程,可能會出現(xiàn)閉環(huán)問題,未來將進(jìn)一步研究存在閉環(huán)問題的煉油化工過程報警根源分析。

    [1] DAHLSTRAND F. Consequence analysis theory for alarm analysis[J]. Knowledge-Based Systems, 2002, 15(1): 27-36.

    [2] WAN Y, YANG F, Lü N, et al. Statistical root cause analysis of novel faults based on digraph models[J]. Chemical Engineering Research and Design, 2013, 91(1): 87-99.

    [3] YANG F, SIRISH L S, XIAO D. 2010 Conference on Control and Fault-Tolerant Systems (SysTol)[C]//France:IEEE, 2010: 387-392.

    [4] MAURYA M R, RENGASWAMY R, VENKATASUBRAMANIAN V. A systematic framework for the development and analysis of signed digraphs for chemical processes 1 Algorithms and analysis[J]. Industrial & Engineering Chemistry Research, 2003, 42(20): 4789-4810.

    [5] MAURYA M R, RENGASWAMY R, VENKATASUBRAMANIAN V. A systematic framework for the development and analysis of signed digraphs for chemical processes 2 Control loops and flow sheet analysis[J]. Industrial & Engineering Chemistry Research, 2003, 42(20): 4811-4827.

    [6] 李安峰, 夏濤, 張貝克, 等. 化工過程SDG建模方法[J].系統(tǒng)仿真學(xué)報, 2003, 15(10): 1364-1368. (LI Anfeng, XIA Tao, ZHANG Beike, et al. SDG modeling approach for chemical engineering process[J]. Journal of System Simulation, 2003, 15(10): 1364-1368.)

    [7] 王偉, 于達(dá)仁, 趙輝, 等. 基于符號有向圖模型的故障診斷方法[J].動力工程學(xué)報, 2007, 27(5): 736-741. (WANG Wei, YU Daren, ZHAO Hui, et al. Fault diagnosis method based on sign directed graph models[J]. Journal of Power Engineering, 2007, 27(5): 736-741.)

    [8] 朱群雄, 高慧慧, 劉菲菲, 等. 過程工業(yè)報警系統(tǒng)研究進(jìn)展[J].計算機(jī)與應(yīng)用化學(xué), 2014, 31(2): 129-134. (ZHU Qunxiong, GAO Huihui, LIU Feifei, et al. Research progress of alarm systems in the process industries[J]. Computers and Applied Chemistry, 2014, 31(2): 129-134.)

    [9] BAUER M, THORNHILL N F. A practical method for identifying the propagation path of plant-wide disturbances[J].Journal of Process Control, 2008, 18(7): 707-719.

    [10] DUAN P, YANG F, CHEN T, et al. American Control Conference (ACC)[C]//Canada: IEEE, 2012: 3522-3527.

    [11] BAUER M, COX J W, CAVENESS M H, et al. Finding the direction of disturbance propagation in a chemical process using transfer entropy[J].IEEE Transactions on Control Systems Technology, 2007, 15(1): 12-21.

    [12] ABELE L, ANIC M, GUTMANN T, et al. Combining knowledge modeling and machine learning for alarm root cause analysis[J].IFAC Proceedings Volumes, 2013, 46(9): 1843-1848.

    [13] HAN L, GAO H, XU Y, et al. Combining FAP, MAP and correlation analysis for multivariate alarm thresholds optimization in industrial process[J].Journal of Loss Prevention in the Process Industries, 2016, 40: 471-478.

    [14] GAO J, TULSYAN A, YANG F, et al. A transfer entropy method to quantify causality in stochastic nonlinear systems[J].IFAC Proceedings Volumes, 2016, 49(7): 454-459.

    [15] GAO H, XU Y, GU X, et al. Systematic rationalization approach for multivariate correlated alarms based on interpretive structural modeling and Likert scale[J].Chinese Journal of Chemical Engineering, 2015, 23(12): 1987-1996.

    [16] PENG D, GU X, XU Y, et al. Integrating probabilistic signed digraph and reliability analysis for alarm signal optimization in chemical plant[J].Journal of Loss Prevention in the Process Industries, 2015, 33: 279-288.

    猜你喜歡
    分餾塔根源煉油
    粉煤熱解分餾系統(tǒng)主要工藝和控制
    廣州化工(2022年13期)2022-08-01 01:11:14
    探究煉油設(shè)備腐蝕與防護(hù)技術(shù)
    中國石化煉油銷售有限公司
    石油瀝青(2020年1期)2020-05-25 06:53:56
    一道電偏轉(zhuǎn)創(chuàng)新題的命題根源剖析
    幫孩子找出問題的根源
    從小到大,由大走強(qiáng)—40載煉油節(jié)節(jié)高
    “?!睕r空前 煉油之踵
    湊合是離婚的根源
    焦化分餾塔頂結(jié)鹽問題的分析及處理
    化工科技(2016年1期)2016-06-04 23:15:00
    傳統(tǒng)媒體商業(yè)模式坍塌的根源
    聲屏世界(2015年4期)2015-02-28 15:19:33
    国产男人的电影天堂91| 国产精品99久久久久久久久| ponron亚洲| 色噜噜av男人的天堂激情| 欧美一区二区国产精品久久精品| 亚洲欧美一区二区三区国产| 亚洲三级黄色毛片| 全区人妻精品视频| 我的女老师完整版在线观看| 插阴视频在线观看视频| 十八禁国产超污无遮挡网站| 亚洲欧美精品自产自拍| 菩萨蛮人人尽说江南好唐韦庄 | 欧美成人午夜免费资源| 国产精品乱码一区二三区的特点| 2021天堂中文幕一二区在线观| 国产亚洲午夜精品一区二区久久 | 99国产精品一区二区蜜桃av| 午夜福利成人在线免费观看| 亚洲成人av在线免费| 精品久久久久久久久亚洲| 18禁在线播放成人免费| 久久久国产成人精品二区| 久久久久久久久中文| 欧美一区二区精品小视频在线| 日韩三级伦理在线观看| 欧美xxxx黑人xx丫x性爽| 欧美另类亚洲清纯唯美| 在线a可以看的网站| 久久精品久久久久久噜噜老黄 | 亚洲av二区三区四区| 亚洲精品成人久久久久久| 欧美一区二区亚洲| 中文字幕制服av| 日韩国内少妇激情av| 只有这里有精品99| 最近中文字幕高清免费大全6| 国产精品国产三级国产专区5o | 嫩草影院精品99| 亚洲欧美一区二区三区国产| 欧美人与善性xxx| 日本与韩国留学比较| 噜噜噜噜噜久久久久久91| 国产麻豆成人av免费视频| 久久这里有精品视频免费| av在线蜜桃| 欧美精品国产亚洲| 成人性生交大片免费视频hd| 久久久久免费精品人妻一区二区| 五月伊人婷婷丁香| av在线亚洲专区| 熟女电影av网| 色哟哟·www| 国产片特级美女逼逼视频| 又爽又黄a免费视频| 成年版毛片免费区| 午夜福利在线在线| 国产国拍精品亚洲av在线观看| 两个人视频免费观看高清| 能在线免费看毛片的网站| 九九爱精品视频在线观看| 国语自产精品视频在线第100页| 精华霜和精华液先用哪个| 国产黄色小视频在线观看| 超碰97精品在线观看| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 久久午夜福利片| 久久精品国产鲁丝片午夜精品| 一个人看的www免费观看视频| 国内揄拍国产精品人妻在线| 亚洲四区av| 亚洲在线观看片| 身体一侧抽搐| 日本猛色少妇xxxxx猛交久久| 麻豆成人午夜福利视频| 久久亚洲国产成人精品v| 自拍偷自拍亚洲精品老妇| 国产精品女同一区二区软件| av线在线观看网站| 永久免费av网站大全| 精品久久久久久久久av| 亚洲一区高清亚洲精品| 免费无遮挡裸体视频| 亚洲国产精品专区欧美| 成人av在线播放网站| 国产在线男女| 蜜臀久久99精品久久宅男| 青春草亚洲视频在线观看| 国产亚洲一区二区精品| 久久久久久久久久成人| 免费av毛片视频| 久久久久九九精品影院| 联通29元200g的流量卡| 91午夜精品亚洲一区二区三区| 91久久精品国产一区二区三区| 韩国av在线不卡| 久久久久免费精品人妻一区二区| 亚洲精品乱码久久久久久按摩| 日韩av不卡免费在线播放| 国产精品美女特级片免费视频播放器| 最近中文字幕高清免费大全6| 久久韩国三级中文字幕| 国产精品精品国产色婷婷| 欧美bdsm另类| 国产一区二区三区av在线| 久久亚洲精品不卡| 亚洲欧美精品专区久久| 亚洲自拍偷在线| 亚洲自拍偷在线| 最新中文字幕久久久久| 3wmmmm亚洲av在线观看| 亚洲欧美成人综合另类久久久 | 三级国产精品片| 天堂网av新在线| 久久久久国产网址| 亚洲性久久影院| 免费av观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品99久久久久久久久| 国产精品人妻久久久久久| 国产色婷婷99| 国产成人一区二区在线| 国产精品不卡视频一区二区| 亚洲激情五月婷婷啪啪| 久久精品夜夜夜夜夜久久蜜豆| 成人特级av手机在线观看| 国产黄色小视频在线观看| 青青草视频在线视频观看| 久久久久久大精品| 69av精品久久久久久| 午夜爱爱视频在线播放| 国产精品永久免费网站| 亚洲国产日韩欧美精品在线观看| 欧美激情久久久久久爽电影| 国产女主播在线喷水免费视频网站 | 超碰97精品在线观看| 18+在线观看网站| 熟妇人妻久久中文字幕3abv| 亚洲av福利一区| 秋霞在线观看毛片| 一个人看的www免费观看视频| 亚洲av免费在线观看| 水蜜桃什么品种好| 亚洲色图av天堂| 久久久久久久久久黄片| a级毛色黄片| 一个人看的www免费观看视频| av又黄又爽大尺度在线免费看 | 成年女人看的毛片在线观看| av在线天堂中文字幕| 成人亚洲精品av一区二区| 美女脱内裤让男人舔精品视频| 国产精品三级大全| 久久精品国产鲁丝片午夜精品| 成人毛片60女人毛片免费| 我要看日韩黄色一级片| 国产亚洲av片在线观看秒播厂 | 久久久久精品久久久久真实原创| 成人性生交大片免费视频hd| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 18禁动态无遮挡网站| 色网站视频免费| 99在线人妻在线中文字幕| 啦啦啦啦在线视频资源| av在线蜜桃| 一级毛片久久久久久久久女| a级毛色黄片| 黄色欧美视频在线观看| 国产精品人妻久久久久久| 高清视频免费观看一区二区 | 亚洲天堂国产精品一区在线| 乱码一卡2卡4卡精品| 老女人水多毛片| 18禁动态无遮挡网站| 熟女人妻精品中文字幕| 六月丁香七月| 亚洲综合精品二区| 人妻少妇偷人精品九色| 国产成人aa在线观看| 国产成人a区在线观看| 人人妻人人澡人人爽人人夜夜 | av专区在线播放| 精品酒店卫生间| 高清午夜精品一区二区三区| 欧美日韩在线观看h| 国产成人福利小说| 国产一区二区在线av高清观看| 欧美激情在线99| 六月丁香七月| 非洲黑人性xxxx精品又粗又长| 欧美97在线视频| 精品一区二区三区人妻视频| 黑人高潮一二区| 久久国内精品自在自线图片| 高清在线视频一区二区三区 | 久久精品久久久久久噜噜老黄 | 精品酒店卫生间| 看非洲黑人一级黄片| 校园人妻丝袜中文字幕| 久久精品国产亚洲网站| 日本黄色片子视频| 搡女人真爽免费视频火全软件| 久久久久久大精品| 51国产日韩欧美| 一级av片app| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在| 国产一区二区亚洲精品在线观看| 又爽又黄a免费视频| 久久99精品国语久久久| 亚洲国产色片| 国产精品久久久久久久电影| 一区二区三区四区激情视频| 97超碰精品成人国产| 日日撸夜夜添| 国产精品久久久久久久久免| 久久久午夜欧美精品| 久久午夜福利片| 91精品一卡2卡3卡4卡| www.色视频.com| 麻豆乱淫一区二区| 一区二区三区高清视频在线| 一边摸一边抽搐一进一小说| 国产精品乱码一区二三区的特点| 97在线视频观看| 别揉我奶头 嗯啊视频| 久久精品人妻少妇| 禁无遮挡网站| 亚洲成人av在线免费| 一级二级三级毛片免费看| 国产大屁股一区二区在线视频| 黄片无遮挡物在线观看| 久99久视频精品免费| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 色尼玛亚洲综合影院| 只有这里有精品99| 欧美xxxx性猛交bbbb| 在线a可以看的网站| 国产精品女同一区二区软件| 国产午夜精品一二区理论片| 国产av不卡久久| 亚洲美女搞黄在线观看| 免费观看a级毛片全部| 国产淫片久久久久久久久| 大香蕉久久网| 久久99热这里只频精品6学生 | 日产精品乱码卡一卡2卡三| 性色avwww在线观看| 色5月婷婷丁香| 国产成年人精品一区二区| 晚上一个人看的免费电影| 久久久久久九九精品二区国产| av国产久精品久网站免费入址| 中国美白少妇内射xxxbb| 国产精品人妻久久久影院| 精品欧美国产一区二区三| 1024手机看黄色片| 欧美精品国产亚洲| 99久久无色码亚洲精品果冻| 成人亚洲精品av一区二区| 国产午夜精品一二区理论片| 精品一区二区三区人妻视频| 国产女主播在线喷水免费视频网站 | 好男人在线观看高清免费视频| 国产精品人妻久久久久久| 亚洲精品色激情综合| 欧美成人免费av一区二区三区| 亚洲精品乱码久久久久久按摩| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 91久久精品国产一区二区成人| 久久亚洲国产成人精品v| 国产男人的电影天堂91| 一级爰片在线观看| 亚洲一级一片aⅴ在线观看| 女的被弄到高潮叫床怎么办| av播播在线观看一区| av女优亚洲男人天堂| 亚洲成色77777| 国产黄色小视频在线观看| 小蜜桃在线观看免费完整版高清| 国产一区二区在线av高清观看| 国产综合懂色| 精品久久久久久电影网 | 亚洲怡红院男人天堂| 99热这里只有是精品在线观看| 国产精品一区二区性色av| 91在线精品国自产拍蜜月| 国产在线男女| 欧美性猛交黑人性爽| 国产一区有黄有色的免费视频 | 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 少妇被粗大猛烈的视频| 中文字幕制服av| 亚洲自拍偷在线| 最近2019中文字幕mv第一页| ponron亚洲| 大香蕉97超碰在线| 插逼视频在线观看| 性色avwww在线观看| 女人久久www免费人成看片 | 纵有疾风起免费观看全集完整版 | 最近视频中文字幕2019在线8| 欧美97在线视频| 国产私拍福利视频在线观看| 99久久成人亚洲精品观看| 最新中文字幕久久久久| 国产精品久久电影中文字幕| 精品国产三级普通话版| 午夜福利成人在线免费观看| 1000部很黄的大片| 视频中文字幕在线观看| 韩国高清视频一区二区三区| 黄色欧美视频在线观看| 两个人的视频大全免费| 国产精品一区二区三区四区久久| 亚洲综合精品二区| 天堂网av新在线| .国产精品久久| 久久人人爽人人爽人人片va| av卡一久久| 内地一区二区视频在线| 亚洲第一区二区三区不卡| 亚洲av免费在线观看| 午夜爱爱视频在线播放| 麻豆一二三区av精品| 中文天堂在线官网| 亚洲精品国产成人久久av| 久久久久性生活片| 两性午夜刺激爽爽歪歪视频在线观看| 精品人妻视频免费看| 18禁在线播放成人免费| 麻豆精品久久久久久蜜桃| 大香蕉久久网| 国国产精品蜜臀av免费| 中文字幕免费在线视频6| 亚洲中文字幕一区二区三区有码在线看| 乱人视频在线观看| 一级黄片播放器| 联通29元200g的流量卡| ponron亚洲| 99久久无色码亚洲精品果冻| 亚洲国产日韩欧美精品在线观看| 色综合亚洲欧美另类图片| 三级国产精品片| 黄色配什么色好看| 久久精品综合一区二区三区| 久久久a久久爽久久v久久| 精华霜和精华液先用哪个| 久久精品国产亚洲网站| 成人毛片60女人毛片免费| 99久久无色码亚洲精品果冻| 亚洲性久久影院| 精品久久国产蜜桃| 天堂影院成人在线观看| 国产 一区精品| 国产真实乱freesex| 久久鲁丝午夜福利片| 亚洲18禁久久av| 老司机福利观看| 亚洲精品乱码久久久久久按摩| 亚洲高清免费不卡视频| 欧美人与善性xxx| 成年女人永久免费观看视频| 欧美极品一区二区三区四区| 亚洲av日韩在线播放| 好男人视频免费观看在线| 国产亚洲精品av在线| 国产精品人妻久久久影院| 国产v大片淫在线免费观看| 寂寞人妻少妇视频99o| 最近视频中文字幕2019在线8| 亚洲av电影不卡..在线观看| 两个人的视频大全免费| 日本爱情动作片www.在线观看| 国内精品宾馆在线| 国产成人一区二区在线| 七月丁香在线播放| 久久久国产成人精品二区| 22中文网久久字幕| 少妇猛男粗大的猛烈进出视频 | 亚洲精品影视一区二区三区av| 久久久久免费精品人妻一区二区| 美女黄网站色视频| 国产精品国产三级国产专区5o | 一边摸一边抽搐一进一小说| 国产亚洲精品av在线| 久久久久久久久久黄片| 国产精品久久久久久久久免| 精品酒店卫生间| 91久久精品国产一区二区三区| 99久国产av精品| 欧美日本亚洲视频在线播放| 亚洲av不卡在线观看| 3wmmmm亚洲av在线观看| av黄色大香蕉| 亚洲av福利一区| av专区在线播放| 日本猛色少妇xxxxx猛交久久| 麻豆一二三区av精品| 亚洲最大成人av| 伦精品一区二区三区| 天美传媒精品一区二区| 人人妻人人澡人人爽人人夜夜 | 国产人妻一区二区三区在| 久久热精品热| 成人鲁丝片一二三区免费| 午夜日本视频在线| 中文乱码字字幕精品一区二区三区 | 免费黄色在线免费观看| 日韩,欧美,国产一区二区三区 | 免费人成在线观看视频色| 久久精品人妻少妇| 亚洲高清免费不卡视频| 99热精品在线国产| 男人舔奶头视频| 视频中文字幕在线观看| 内地一区二区视频在线| 麻豆久久精品国产亚洲av| 成人午夜高清在线视频| 少妇的逼水好多| 国产精品一区二区三区四区久久| 一卡2卡三卡四卡精品乱码亚洲| 婷婷色麻豆天堂久久 | 你懂的网址亚洲精品在线观看 | 国产午夜精品久久久久久一区二区三区| 免费看av在线观看网站| 亚洲怡红院男人天堂| 又黄又爽又刺激的免费视频.| 少妇的逼好多水| 亚洲成av人片在线播放无| 一夜夜www| 夫妻性生交免费视频一级片| 欧美色视频一区免费| 男人舔奶头视频| 国产精品伦人一区二区| 高清日韩中文字幕在线| 天堂av国产一区二区熟女人妻| 日韩强制内射视频| 精品酒店卫生间| 亚洲欧美清纯卡通| 日韩精品青青久久久久久| 女人被狂操c到高潮| 色噜噜av男人的天堂激情| videos熟女内射| 国产成人福利小说| 精品少妇黑人巨大在线播放 | 我的老师免费观看完整版| 欧美xxxx性猛交bbbb| 国产 一区精品| 在线播放国产精品三级| 久久人人爽人人爽人人片va| 午夜精品一区二区三区免费看| 久久久久精品久久久久真实原创| 日本一本二区三区精品| 久久久精品欧美日韩精品| 国产精品乱码一区二三区的特点| 国产免费视频播放在线视频 | 免费av毛片视频| 国产探花在线观看一区二区| 亚洲av电影在线观看一区二区三区 | 嘟嘟电影网在线观看| 久久婷婷人人爽人人干人人爱| av在线老鸭窝| 最近视频中文字幕2019在线8| 国产伦在线观看视频一区| 舔av片在线| 中文字幕精品亚洲无线码一区| 国产 一区 欧美 日韩| 亚洲欧美成人综合另类久久久 | 久久99热6这里只有精品| 亚洲av中文av极速乱| 99久久无色码亚洲精品果冻| 人人妻人人澡人人爽人人夜夜 | 亚洲内射少妇av| av在线老鸭窝| 菩萨蛮人人尽说江南好唐韦庄 | 午夜福利在线观看免费完整高清在| 国语对白做爰xxxⅹ性视频网站| 一级黄色大片毛片| 亚洲高清免费不卡视频| 三级经典国产精品| 毛片一级片免费看久久久久| 女人久久www免费人成看片 | 老司机影院毛片| 久久午夜福利片| 小说图片视频综合网站| 97热精品久久久久久| 色网站视频免费| 亚洲va在线va天堂va国产| 在线天堂最新版资源| 日本一二三区视频观看| 亚洲国产高清在线一区二区三| 亚洲欧美清纯卡通| 免费观看在线日韩| 男人的好看免费观看在线视频| 99视频精品全部免费 在线| 能在线免费看毛片的网站| 国产老妇女一区| 午夜免费激情av| 免费观看在线日韩| 久久久国产成人免费| 国产女主播在线喷水免费视频网站 | 国产国拍精品亚洲av在线观看| 久久久a久久爽久久v久久| 99久久中文字幕三级久久日本| 亚洲天堂国产精品一区在线| 国产精品久久久久久精品电影| 老司机福利观看| 久久精品久久久久久噜噜老黄 | 国产一区二区在线观看日韩| 18禁动态无遮挡网站| 亚洲精品日韩av片在线观看| 国产精品伦人一区二区| 亚洲精品色激情综合| 国产伦一二天堂av在线观看| 国语对白做爰xxxⅹ性视频网站| 久久久午夜欧美精品| 色尼玛亚洲综合影院| 一级毛片电影观看 | 日本免费a在线| 插逼视频在线观看| 成年免费大片在线观看| 亚洲av免费在线观看| 久久久久久久久久成人| 午夜精品在线福利| 国产 一区精品| 夜夜爽夜夜爽视频| 亚洲av熟女| 中文天堂在线官网| 亚洲人成网站在线播| 亚洲欧洲日产国产| 国产精品一区二区在线观看99 | 九九爱精品视频在线观看| 男女边吃奶边做爰视频| 女人久久www免费人成看片 | 在线观看一区二区三区| 99热精品在线国产| 久久韩国三级中文字幕| 日韩精品青青久久久久久| 国产乱人偷精品视频| 激情 狠狠 欧美| 国产熟女欧美一区二区| 成人鲁丝片一二三区免费| 亚洲av免费在线观看| 免费看光身美女| av女优亚洲男人天堂| 在线观看一区二区三区| 少妇的逼水好多| 成人无遮挡网站| 日日啪夜夜撸| 中文字幕人妻熟人妻熟丝袜美| 国产精品精品国产色婷婷| 欧美日韩一区二区视频在线观看视频在线 | 黄片wwwwww| 成人一区二区视频在线观看| 色播亚洲综合网| 亚洲精品日韩在线中文字幕| 精品久久久久久成人av| 狂野欧美激情性xxxx在线观看| 九九在线视频观看精品| 久久精品国产99精品国产亚洲性色| 日韩三级伦理在线观看| www.av在线官网国产| 简卡轻食公司| 国产伦精品一区二区三区视频9| 免费av观看视频| 免费观看的影片在线观看| 久久国产乱子免费精品| 秋霞伦理黄片| 97热精品久久久久久| 亚洲自拍偷在线| 啦啦啦韩国在线观看视频| 性插视频无遮挡在线免费观看| 日日啪夜夜撸| 亚洲在线自拍视频| 久久99热这里只频精品6学生 | 禁无遮挡网站| 精品久久久久久久久亚洲| 成人午夜精彩视频在线观看| 一区二区三区乱码不卡18| 中文欧美无线码| 99热精品在线国产| 我要搜黄色片| 国产av码专区亚洲av| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 黄色一级大片看看| 91久久精品国产一区二区成人| 高清午夜精品一区二区三区| 最近手机中文字幕大全| 亚洲精品乱码久久久久久按摩| 建设人人有责人人尽责人人享有的 | 精品久久久久久久久久久久久| 亚洲伊人久久精品综合 | 一区二区三区高清视频在线| 日韩一区二区三区影片| 亚洲精品成人久久久久久| 亚洲美女搞黄在线观看| 国产欧美另类精品又又久久亚洲欧美| av专区在线播放| 欧美一区二区亚洲| 亚洲国产高清在线一区二区三| 午夜福利成人在线免费观看| 噜噜噜噜噜久久久久久91| 寂寞人妻少妇视频99o| 亚洲图色成人| 免费观看人在逋| 国产乱人视频| 久久久久久久亚洲中文字幕| av卡一久久|