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

    基于Labview的不確定環(huán)境下化工廠蒸汽動力系統(tǒng)多目標(biāo)建模與實(shí)時操作優(yōu)化

    2017-07-01 22:37:30石蓬軍曹萃文顧幸生
    關(guān)鍵詞:減壓器發(fā)電機(jī)組間隔

    石蓬軍, 曹萃文, 顧幸生

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

    ?

    基于Labview的不確定環(huán)境下化工廠蒸汽動力系統(tǒng)多目標(biāo)建模與實(shí)時操作優(yōu)化

    石蓬軍, 曹萃文, 顧幸生

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

    針對化工廠蒸汽動力系統(tǒng)運(yùn)行時操作單元蒸汽供需在不確定的時間以不確定的組合形式變化而造成的難以進(jìn)行實(shí)時操作優(yōu)化的問題,建立了多目標(biāo)非線性混合整數(shù)滾動規(guī)劃模型(MO-MINLP),并提出以Labview軟件為集成平臺,利用其數(shù)據(jù)采集與處理模塊實(shí)現(xiàn)對蒸汽動力系統(tǒng)操作運(yùn)行數(shù)據(jù)的實(shí)時采集與處理;按照實(shí)時采集的描述不確定事件開始與結(jié)束的時刻劃分滾動優(yōu)化周期的時間間隔長度。利用上述方法對某真實(shí)化工廠的蒸汽動力系統(tǒng)的數(shù)據(jù)進(jìn)行仿真,并實(shí)時將多目標(biāo)優(yōu)化決策結(jié)果以圖形化的形式顯示在Labview界面上。仿真結(jié)果表明,此方法很好地解決了理論與工業(yè)銜接的問題,具有良好的應(yīng)用價值。

    Labview; 蒸氣動力系統(tǒng); 不確定環(huán)境; 實(shí)時操作優(yōu)化; 多目標(biāo)MINLP

    化工廠蒸汽動力系統(tǒng)的真實(shí)運(yùn)行環(huán)境是動態(tài)和不確定的,如何實(shí)現(xiàn)既能滿足企業(yè)各裝置用汽需求,又能兼顧減少企業(yè)運(yùn)行成本、實(shí)現(xiàn)節(jié)能減排等多種目標(biāo),成為近年來相關(guān)領(lǐng)域的研究熱點(diǎn)。國內(nèi)外相關(guān)領(lǐng)域的研究人員做了大量的研究,其中以數(shù)學(xué)規(guī)劃法最為突出。這種方法首先根據(jù)蒸汽動力系統(tǒng)中的各種設(shè)備建立數(shù)學(xué)模型,然后依據(jù)蒸汽動力系統(tǒng)的拓?fù)浣Y(jié)構(gòu),根據(jù)物料及能量守恒建立以一個或者多個目標(biāo)為優(yōu)化目標(biāo)的整個蒸汽動力系統(tǒng)的數(shù)學(xué)模型。文獻(xiàn)[1-3] 分別建立了LP、NLP數(shù)學(xué)規(guī)劃模型用以解決蒸汽動力系統(tǒng)多目標(biāo)運(yùn)行優(yōu)化問題,但沒有考慮到蒸汽動力系統(tǒng)的運(yùn)行優(yōu)化常伴隨著設(shè)備的選擇及啟停變化。文獻(xiàn)[4-9]分別考慮了在優(yōu)化模型中加入整數(shù)的因素建立混合整數(shù)優(yōu)化模型,將設(shè)備的啟停及選擇采用0-1變量進(jìn)行描述,含有這種變量的數(shù)學(xué)模型稱為混合整數(shù)規(guī)劃模型(MILP、MINLP)。其中Dai等[6]通過改進(jìn)的粒子群算法能夠很快地得到使得目標(biāo)最優(yōu)化的蒸汽動力系統(tǒng)的運(yùn)行策略。Chen等[7]以最小化年總成本為目標(biāo)確定超結(jié)構(gòu)和蒸汽動力系統(tǒng)的最優(yōu)化運(yùn)行配置。Zhang等[8]以相同的方法研究蒸汽動力系統(tǒng)在滿足工藝運(yùn)行需求的前提下最小化能源成本。上述研究均假定蒸汽動力系統(tǒng)處于確定環(huán)境中,然而實(shí)際的蒸汽動力系統(tǒng)運(yùn)行過程中不可避免地伴隨著各種不確定因素,包括沒有規(guī)律的設(shè)備故障、生產(chǎn)環(huán)境的變化以及具有周期性規(guī)律的供需變化。文獻(xiàn)[10-17]考慮了蒸汽動力系統(tǒng)運(yùn)行過程中可能存在的不確定因素,對蒸汽動力系統(tǒng)以不同的優(yōu)化目標(biāo)及策略進(jìn)行了研究。其中張冰劍等[10]建立了以設(shè)備維修時間為不確定因素的多周期運(yùn)行優(yōu)化模型。羅向龍等[11]提出將整個運(yùn)行過程以8 h為1個周期進(jìn)行多周期運(yùn)行優(yōu)化。另外羅向龍等[13]考慮了蒸汽動力系統(tǒng)的汽電需求的不確定因素,將蒸汽動力系統(tǒng)的運(yùn)行過程分為等間隔的多周期運(yùn)行過程。李暉等[15]將蒸汽動力系統(tǒng)優(yōu)化過程中的不確定因素分為汽電需求的確定性變化,并將其分周期處理。池曉等[17]以蒸汽需求量作為不確定參數(shù),在不確定參數(shù)的可行域中將其離散化轉(zhuǎn)換為確定性的多周期規(guī)劃問題。

    在解決蒸汽動力系統(tǒng)運(yùn)行優(yōu)化過程中的不確定問題時,很多學(xué)者將蒸汽動力系統(tǒng)的運(yùn)行過程分成多周期分別進(jìn)行優(yōu)化,這種優(yōu)化需要滿足某些運(yùn)行環(huán)境必須是周期性變化的,且必須事先知道某個周期的運(yùn)行時間長度,而實(shí)際的運(yùn)行過程中多數(shù)的不確定因素是沒有規(guī)律可循的,更不可能準(zhǔn)確地知道在某個狀態(tài)下蒸汽動力系統(tǒng)能夠持續(xù)運(yùn)行的時間長度。針對上述問題,本文提出以Labview軟件為集成平臺,利用其數(shù)據(jù)采集與處理模塊實(shí)現(xiàn)對蒸汽動力系統(tǒng)操作運(yùn)行數(shù)據(jù)的實(shí)時采集與處理;按照實(shí)時采集的描述不確定事件開始與結(jié)束的時刻劃分滾動優(yōu)化周期的時間間隔長度。利用上述方法對某真實(shí)化工廠的蒸汽動力系統(tǒng)的數(shù)據(jù)進(jìn)行仿真,并實(shí)時將多目標(biāo)優(yōu)化決策結(jié)果以圖形化的形式顯示在Labview界面上。仿真結(jié)果表明,本文方法可以很好地解決理論與工業(yè)銜接的問題。

    1 MO-MINLP數(shù)學(xué)模型

    1.1 目標(biāo)函數(shù)

    目標(biāo)1為最小化蒸汽動力系統(tǒng)在t時間間隔、每單位時間的總運(yùn)行成本(CNY/h),其目標(biāo)函數(shù)表達(dá)式如下:

    (3)

    (4)

    (5)

    其中:式(1)中Cpower,t前的負(fù)號表示蒸汽動力系統(tǒng)運(yùn)行過程中產(chǎn)生的電力可以用于自身使用從而減少來自外部電力的成本;式(2)為nb臺燃煤鍋爐在t時間間隔、單位時間燃燒煤的總成本(CNY/h);式(3)為nb臺燃煤鍋爐及nr臺減溫減壓器在t時間間隔、單位時間使用水的總成本(CNY/h);式(4)為蒸汽動力系統(tǒng)中所有設(shè)備在t時間間隔、單位時間的總損耗成本(CNY/h);式(5)為ng組發(fā)電機(jī)組和nt臺汽輪機(jī)在t時間間隔、單位時間產(chǎn)生電力的總收益。

    目標(biāo)2為最小化蒸汽動力系統(tǒng)在t時間間隔、各壓力等級蒸汽管線上單位時間內(nèi)蒸汽的總放空量,其目標(biāo)函數(shù)表達(dá)式如下:

    (6)

    其中:EXTt表示蒸汽動力系統(tǒng)在t時間間隔內(nèi)、單位時間各壓力等級的總蒸汽放空量(t/h)。

    1.2 約束條件

    1.2.1 主要設(shè)備的運(yùn)行模型

    (1) 燃煤鍋爐:鍋爐作為整個蒸汽動力系統(tǒng)的主要產(chǎn)汽單元,通過能量平衡建立其運(yùn)行模型如下:

    (7)

    (8)

    (9)

    式(8)表示在t時間間隔,流入減溫減壓器的減溫水的能量加上流入減溫減壓器的高壓蒸汽的能量等于流出減溫減壓器的低壓蒸汽的能量,SRk,t為減溫減壓器啟停狀態(tài)的0-1變量。式(9)表示在t時間間隔,流入減溫減壓器的減溫水的質(zhì)量流量加上流入減溫減壓器的蒸汽的質(zhì)量流量等于流出減溫減壓器的低壓蒸汽的質(zhì)量流量。

    (2) 汽輪機(jī):汽輪機(jī)采用背壓式汽輪機(jī),它能夠帶動發(fā)電機(jī)組將蒸汽的熱能轉(zhuǎn)換為動能進(jìn)而轉(zhuǎn)換為電能,并且能夠?qū)⒂捎谧龉Χ禍販p壓的蒸汽作為低壓蒸汽供給低壓蒸汽用戶使用,表達(dá)式如下:

    (10)

    (3) 發(fā)電機(jī)組:流入發(fā)電機(jī)組的蒸汽的質(zhì)量流量與發(fā)電功率的關(guān)系與汽輪機(jī)類似,只是發(fā)電機(jī)組出口的蒸汽不再作為汽源回到蒸汽管網(wǎng)中去。SGk,t為發(fā)電機(jī)組的啟停狀態(tài)的0-1變量。

    (11)

    1.2.2 蒸汽供需平衡約束條件 蒸汽動力系統(tǒng)為各生產(chǎn)過程提供蒸汽,在每一壓力等級的蒸汽管網(wǎng)上必須滿足供汽量大于或者等于用汽量,高壓力等級上的用汽單元可以是低壓力等級的供汽單元,各壓力等級的平衡關(guān)系滿足如下一般平衡公式:

    (12)

    1.2.3 設(shè)備負(fù)載及流量上下限 在蒸汽動力系統(tǒng)的運(yùn)行過程中,為了保證蒸汽動力系統(tǒng)安全、連續(xù)地正常運(yùn)行,需要為每一個設(shè)備設(shè)置一定的負(fù)載限制,對某些流量也進(jìn)行合理的限制:

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    式(13)為鍋爐生產(chǎn)蒸汽流量的上下限;式(14)為汽輪機(jī)的蒸汽流量上下限;式(15)為發(fā)電機(jī)組蒸汽流量限制;式(16)為減溫減壓器入口蒸汽流量限制;式(17)為低壓力等級管線用汽單元不確定蒸汽流量的上下限;式(18)為低壓力等級廢熱鍋爐和余熱鍋爐兩類供汽單元不確定蒸汽產(chǎn)量的上下限。

    1.3 模型中的變長時間間隔劃分及滾動優(yōu)化

    蒸汽動力系統(tǒng)運(yùn)行過程中存在操作單元供需變化的事件,這些事件的開始和結(jié)束時刻對應(yīng)某些蒸汽管線的相關(guān)參數(shù)發(fā)生變化。用這些事件的開始與結(jié)束時刻確定當(dāng)前的時間間隔長度,這種時間間隔一定是不等長的。在這些變長時間間隔的開始和結(jié)束時刻,由于原模型中的參數(shù)發(fā)生劇烈變化,需進(jìn)行重新優(yōu)化(稱為“滾動優(yōu)化”)才能正確指導(dǎo)生產(chǎn)的運(yùn)行。優(yōu)化算法采用基于Pareto規(guī)則的多目標(biāo)差分進(jìn)化算法(MODE)。

    2 問題描述及案例研究

    2.1 問題描述

    本文以某真實(shí)化工廠的實(shí)際蒸汽動力系統(tǒng)采集的數(shù)據(jù)作為仿真對象,考慮在操作單元蒸汽供需不確定下對該蒸汽動力系統(tǒng)進(jìn)行多目標(biāo)實(shí)時操作優(yōu)化。圖1所示為某化工廠蒸汽動力系統(tǒng)的工藝流程圖。其中S1壓力等級蒸汽管路上以3臺燃煤鍋爐作為產(chǎn)汽單元,所產(chǎn)生的高壓蒸汽一部分用于汽輪機(jī)組T1和T2發(fā)電,多余的蒸汽通過減溫減壓器D1減溫減壓后進(jìn)入3.82 MPa壓力等級的蒸汽管網(wǎng)中。S2 壓力等級的蒸汽供給發(fā)電機(jī)組G1、G2、蒸汽輪機(jī)組T3、T4發(fā)電,剩余的蒸汽則可以在需要時分別通過減溫減壓器D2降溫降壓后送到S4;或通過減溫減壓器D3降溫降壓后送到S3;或通過減溫減壓器D4降溫降壓后送到S5。S3壓力等級的蒸汽供給汽輪機(jī)組T5、T6發(fā)電,剩余的蒸汽則可以在需要時通過減溫減壓器D5降溫降壓后送到S4。S4壓力等級的蒸汽主要用于各低溫、低壓蒸汽用戶單元所需的蒸汽,S4上還有兩臺余熱鍋爐B4、B5用于產(chǎn)生低溫低壓蒸汽供給S4蒸汽管路,其中多余的蒸汽也可以通過減溫減壓器D6降溫、降壓后輸送到S5。S5壓力等級中的蒸汽主要供給低溫、低壓蒸汽用戶單元使用,也通過兩臺余熱鍋爐B6、B7提供低溫、低壓蒸汽。

    圖1 某化工廠蒸汽動力系統(tǒng)結(jié)構(gòu)圖Fig 1 Structure diagram of steam power system in chemical plant

    在各壓力等級蒸汽管路中,每個壓力等級都有供汽單元和用汽單元,當(dāng)管路中蒸汽過多時則需要進(jìn)行排空處理。實(shí)際生產(chǎn)過程中各用汽單元的用汽量及供汽單元的供汽量往往會產(chǎn)生變化,這使得已有的確定性建模及優(yōu)化方法無法直接應(yīng)用于實(shí)際生產(chǎn)中。本文提出了不確定環(huán)境下化工廠蒸汽動力系統(tǒng)多目標(biāo)建模與實(shí)時滾動操作優(yōu)化方法,有關(guān)前提和假設(shè)如下:

    (1) 蒸汽動力系統(tǒng)中鍋爐、蒸汽輪機(jī)、發(fā)電機(jī)組等設(shè)備的運(yùn)行模型均采用線性化處理后的模型。

    (2) 蒸汽動力系統(tǒng)優(yōu)化的目標(biāo)為最小化蒸汽動力系統(tǒng)的運(yùn)行成本;最小化蒸汽動力系統(tǒng)各壓力等級的蒸汽總放空量;采用基于Pareto準(zhǔn)則的多目標(biāo)差分進(jìn)化算法MODE實(shí)現(xiàn)對MO-MINLP模型的優(yōu)化。所取得的Pareto前沿中,采用文獻(xiàn)[18]中驗(yàn)證過的方法對多個目標(biāo)進(jìn)行歸一化和反歸一化來確定最優(yōu)解。

    (3) 蒸汽動力系統(tǒng)優(yōu)化的決策變量為3臺蒸汽鍋爐的產(chǎn)汽流量;6臺汽輪機(jī)的蒸汽流量;6臺減溫減壓器的入口蒸汽流量;2臺發(fā)電機(jī)組的入口蒸汽流量;描述6臺減溫減壓器啟停狀態(tài)的0-1變量;描述2臺發(fā)電機(jī)組的啟停狀態(tài)的0-1變量。

    (4) 不確定變量為蒸汽壓力等級S4上的2臺余熱鍋爐的供汽流量,由仿真實(shí)時檢測數(shù)據(jù)代入;蒸汽壓力等級S4上的各用汽單元的用汽流量,由仿真實(shí)時檢測數(shù)據(jù)代入;蒸汽壓力等級S5上的兩臺余熱鍋爐的供汽流量,由仿真實(shí)時檢測數(shù)據(jù)代入;蒸汽壓力等級S5上的各用汽單元的用汽流量,由仿真實(shí)時檢測數(shù)據(jù)代入。

    (5) S1~S5壓力等級上的蒸汽壓力及溫度都不變,其熵焓值見表1。

    表1 蒸汽動力系統(tǒng)各壓力等級蒸汽熵焓值

    2.2 案例研究

    如圖1所示,T1~T6為背壓式汽輪機(jī),其背壓輸出蒸汽壓力、溫度及負(fù)載范圍見表2。G1~G2為發(fā)電機(jī)組,負(fù)載范圍見表3,外購電價為0.582 CNY/(kW·h)。B1~B3為燃煤鍋爐,3臺燃煤鍋爐均可產(chǎn)生10 MPa、540 ℃的蒸汽,鍋爐效率均為92.14%,鍋爐的負(fù)載范圍如表4所示,燃煤的價格為650 CNY/t。各不確定參數(shù)負(fù)載范圍見表5。

    采用基于Labview的不確定環(huán)境下滾動優(yōu)化方法對上述案例進(jìn)行優(yōu)化仿真,其中不確定環(huán)境下的變長時間間隔劃分方法見圖2。

    表2 蒸汽輪機(jī)參數(shù)

    表3 發(fā)電機(jī)組負(fù)載范圍及啟停狀態(tài)

    表4 燃煤鍋爐的負(fù)載范圍

    表5 不確定供汽、用汽單元負(fù)載范圍

    3 Labview集成平臺中各個功能模塊的實(shí)現(xiàn)

    Labview是一種圖形化的編程語言,與基于文本的程序開發(fā)語言最大的不同在于其產(chǎn)生的程序是以框圖的形式呈現(xiàn)出來的,其數(shù)據(jù)采集與顯示功能在虛擬儀器方面得到了廣泛的應(yīng)用,在流程工業(yè)中也常被用來作為上位機(jī)程序的開發(fā)。本文設(shè)計的Labview集成平臺主要包括前面板和程序框圖兩個部分,前面板用來顯示數(shù)據(jù),程序框圖用來進(jìn)行數(shù)據(jù)的處理等邏輯控制。Labview集成平臺界面如圖3所示,程序框圖部分如圖4所示,主要包括數(shù)據(jù)采集仿真模塊、算法實(shí)現(xiàn)模塊、結(jié)果顯示模塊。

    圖3 化工廠蒸汽動力系統(tǒng)Labview集成平臺Fig.3 Integrated platform of the steam power system in the chemical plant

    圖4 Labview程序框圖Fig.4 Labview program flow chart

    4 案例仿真結(jié)果及分析

    表6所示為系統(tǒng)優(yōu)化運(yùn)行后的實(shí)時操作數(shù)據(jù),其中tp表示不確定量發(fā)生變化的時間點(diǎn),MOC為當(dāng)前配置下目標(biāo)函數(shù)1的值,EXT表示當(dāng)前配置下目標(biāo)函數(shù)2的值。根據(jù)時間點(diǎn)可以看出每個優(yōu)化周期的時間間隔是不等長的,而且每次優(yōu)化的結(jié)果也是變化的。

    表6 蒸汽動力系統(tǒng)部分實(shí)時操作優(yōu)化數(shù)據(jù)

    表7所示為發(fā)電機(jī)組及減溫減壓器在實(shí)時操作優(yōu)化運(yùn)行過程中的啟停狀態(tài)。圖5所示為單次優(yōu)化結(jié)束后所有種群個體的目標(biāo)值經(jīng)過歸一化后的Pareto分布圖。由Pareto分布圖可以看出每次優(yōu)化兩個相互沖突的目標(biāo)會呈現(xiàn)如圖所示的Pareto前沿。圖6所示為采用Labview集成平臺實(shí)現(xiàn)的某化工廠蒸汽動力系統(tǒng)生產(chǎn)運(yùn)行實(shí)時多目標(biāo)操作優(yōu)化結(jié)果在工藝流程圖中的顯示,其能更直觀地顯示系統(tǒng)的運(yùn)行狀態(tài)。

    圖5 目標(biāo)歸一化后的Pareto前沿Fig.5 Pareto frontier after uniformization

    表7 蒸汽動力系統(tǒng)部分實(shí)時操作優(yōu)化后的設(shè)備啟停狀態(tài)

    圖6 蒸汽動力系統(tǒng)工藝結(jié)構(gòu)及實(shí)時優(yōu)化結(jié)果顯示Fig.6 Process diagram and real-time optimization results display of the steam power system

    5 結(jié) 論

    由仿真結(jié)果及分析可以看出,本文提出的基于Labview的多目標(biāo)混合整數(shù)滾動優(yōu)化方法解決了蒸汽動力系統(tǒng)運(yùn)行過程中蒸汽操作單元供需不確定而難以進(jìn)行優(yōu)化的問題,較之以往通過假定理想情況而進(jìn)行固定時間長度的多周期解決方案更加符合蒸汽動力系統(tǒng)的實(shí)際運(yùn)行環(huán)境。同時采用基于Pareto準(zhǔn)則的混合整數(shù)多目標(biāo)規(guī)劃模型解決了設(shè)備啟停狀態(tài)造成的非線性及多目標(biāo)同時優(yōu)化造成的沖突問題;結(jié)合文獻(xiàn)[18]提出的歸一化方法解決了從Pareto最優(yōu)解集中獲取全局最優(yōu)解的問題。仿真實(shí)驗(yàn)結(jié)果表明,本文提出的基于Labview的化工廠蒸汽動力系統(tǒng)在蒸汽操作單元供需不確定的情況下實(shí)時滾動優(yōu)化方法能夠更好地應(yīng)用于實(shí)際生產(chǎn)環(huán)境,具有一定的實(shí)用價值。

    符號說明:

    Ccoal,t—— 蒸汽動力系統(tǒng)在t時間間隔燃燒煤的成本,CNY/h;

    Cwater,t—— 蒸汽動力系統(tǒng)在t時間間隔用水的成本,CNY/h;

    Cdevice,t—— 蒸汽動力系統(tǒng)在t時間間隔設(shè)備的損耗成本,CNY/h;

    nb—— 蒸汽動力系統(tǒng)中鍋爐的數(shù)量;

    ng—— 蒸汽動力系統(tǒng)中發(fā)電機(jī)組的數(shù)量;

    nl—— 蒸汽動力系統(tǒng)中壓力等級的數(shù)量;

    nr—— 蒸汽動力系統(tǒng)中的減溫減壓器數(shù)量;

    ns—— 某一壓力等級上供汽單元的數(shù)量;

    nt—— 蒸汽動力系統(tǒng)中汽輪機(jī)的數(shù)量;

    nu—— 某一壓力等級上用汽單元的數(shù)量;

    PGk,t—— 發(fā)電機(jī)組k在t時間間隔的發(fā)電功率,kW;

    PTk,t—— 汽輪機(jī)k在t時間間隔的發(fā)電功率,kW;

    SBk,t—— 鍋爐k在t時間間隔的啟停狀態(tài),(0-1);

    SGk,t—— 發(fā)電機(jī)組k在t時間間隔的啟停狀態(tài),(0-1);

    SRk,t—— 減溫減壓器k在t時刻的啟停狀態(tài),(0-1);

    STk,t—— 汽輪機(jī)k在t時間間隔的啟停狀態(tài),(0-1);

    t—— 時間間隔;

    tp—— 不確定參數(shù)發(fā)生變化時的時間點(diǎn);

    Vb—— 鍋爐運(yùn)行過程中的損耗成本,CNY/h;

    Vc—— 煤的價格,CNY/t;

    Ve—— 電價,CNY/t-1;

    Vg—— 發(fā)電機(jī)組運(yùn)行過程中的損耗成本,CNY/h;

    Vr—— 減溫減壓器運(yùn)行過程中的損耗成本,CNY/h;

    Vt—— 汽輪機(jī)運(yùn)行過程中的損耗成本,CNY/h;

    Vw—— 工業(yè)用水的價格,CNY/h;

    ηBk—— 鍋爐k的效率

    [1] PETROULAS T,REKLAITIS G V.Computer-aided synthesis and design of plant utility systems [J].AIChE Journal,1984,30(1):69-78.

    [2] RONG Y,LAHDELMA R.An efficient linear programming model and optimization algorithm for trigeneration [J].Applied Energy,2005,82(1):40-63.

    [3] GUTIERREZ Arriaga C G,SEMA Gonzalez M,PONCE Ortega J M.Multi-objective optimization of steam power plants for sustainable generation of electricity [J].Clean Technologies and Environmental Policy,2013,15(4):551-566.

    [4] ALITVENTAEFF B.Mixed integer non-linear programming and its application to the management of utility networks [J].Engineering Optimization,1991,18(1):183-207.

    [5] OLIVEIRA FRANCISCO A P,MATOS H A.Multi-period synthesis and operational planning of utility systems with environmental concerns [J].Computer and Chemical Engineering,2004,28(5):745-753.

    [6] DAI W Z,YIN H A.A correctional model for multiperiod optimization of steam power systems with consideration of startup/shutdown expenses [J].Applied Mechanics and Materials,2010,40/41:298-303.

    [7] CHEN C L,LIN C Y.Design and optimization of steam distribution systems for steam power plants [J].Industrial and Engineering Chemistry Research,2011,50(13):8097-8109.

    [8] ZHANNG B J,LUO X L,CHEN X Z.Coupling process plants and utility systems for site scale steam integration [J].Industrial and Engineering Chemistry Research,2013,52(41):14627-14636.

    [9] LUO X L,HU J H,ZHAO J.Multi-objective optimization for the design and synthesis of utility systems with emission abatement technology concerns [J].Applied Energy,2014,136:1110-1131.

    [10] 張冰劍,華賁,劉金平,等.石化企業(yè)蒸汽動力系統(tǒng)的多周期優(yōu)化運(yùn)行 [J].華北電力大學(xué)學(xué)報,2005,32(2):90-92.

    [11] 羅向龍,華賁.煉油廠蒸汽動力系統(tǒng)優(yōu)化調(diào)度策略研究 [J].計算機(jī)與應(yīng)用化學(xué),2006,23(11):1057-1060.

    [12] 羅向龍,華賁,張冰劍.石化企業(yè)蒸汽動力系統(tǒng)的多周期運(yùn)行優(yōu)化 [J].計算機(jī)與應(yīng)用化學(xué),2006,23(1):41-45.

    [13] 羅向龍,華賁,尹洪超.蒸汽動力系統(tǒng)多周期設(shè)計與運(yùn)行同步優(yōu)化 [J].陜西科技大學(xué)學(xué)報,2007,25(2):52-57.

    [14] 孫玉梅.生產(chǎn)過程與蒸汽動力系統(tǒng)集成優(yōu)化研究 [D].濟(jì)南:山東大學(xué),2012.

    [15] 李暉,孫力,賀高紅.考慮不確定汽電需求的蒸汽動力系統(tǒng)優(yōu)化設(shè)計 [J].化工學(xué)報,2013,64(1):318-325.

    [16] 蓋麗梅,孫力,劉暢,等.基于帶補(bǔ)償隨機(jī)規(guī)劃的蒸汽動力系統(tǒng)優(yōu)化設(shè)計 [J].化工學(xué)報,2014,65(11):4509-4516.

    [17] 池曉,徐可培,馬學(xué)榮,等.基于LINGO對參數(shù)不確定蒸汽動力系統(tǒng)的優(yōu)化運(yùn)行研究 [J].節(jié)能,2015 (1):21-25.

    [18] 韓庚,曹萃文,顧幸生.某焦化廠供水網(wǎng)絡(luò)多目標(biāo)建模及優(yōu)化[J].化工進(jìn)展,2016,35(4):1033-1041.

    Labview-Based Uncertain Multi-objective Modeling and Real-Time Operational Optimization for Steam Power Systems of Chemical Plants

    SHI Peng-jun, CAO Cui-wen, GU Xing-sheng

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

    In order to overcome the shortcoming of uncertain steam supply and demand of units in the steam power system of chemical plants,this paper proposes a rolling optimization model (MO-MINLP) and develops a novelty.The data acquisition and processing module of integrated platform based on Labview is utilized to collect real-time operational data for the steam power system.The start and end time points of recording uncertain events are used to divide the time interval of every period.Simulation results from the steam power system show that the proposed Labview-based platform can perfectly deal with the integration of the theory and industrial applications and will have a wide application value in other similar systems.

    Labview; steam power system; uncertaint environment; real-time operational optimization; multi-objective MINLP

    1006-3080(2017)03-0375-08

    10.14135/j.cnki.1006-3080.2017.03.013

    2016-10-24

    國家自然科學(xué)基金(61174040,61573144,61673175)

    石蓬軍(1990-),男,碩士生,研究方向?yàn)楣I(yè)過程建模與優(yōu)化。

    曹萃文,E-mail:caocuiwen@ecust.edu.cn

    TP311

    A

    猜你喜歡
    減壓器發(fā)電機(jī)組間隔
    煤氣發(fā)電機(jī)組DEH控制系統(tǒng)的優(yōu)化
    山東冶金(2022年4期)2022-09-14 08:59:30
    間隔問題
    氣體減壓器自動校準(zhǔn)方法的探討
    間隔之謎
    減溫減壓器故障分析與改進(jìn)
    蒸汽系統(tǒng)減溫減壓器質(zhì)量分析及解決措施
    基于PLC控制柴油發(fā)電機(jī)組3D 模型
    八鋼歐冶爐TRT發(fā)電機(jī)組成功并網(wǎng)發(fā)電
    新疆鋼鐵(2015年3期)2015-02-20 14:13:56
    上樓梯的學(xué)問
    德國寶馬閥等量跟蹤快開減溫減壓器的應(yīng)用與改進(jìn)
    91av网一区二区| 国产黄色小视频在线观看| 国产精品免费一区二区三区在线| 国产不卡一卡二| 国产精品嫩草影院av在线观看| 五月伊人婷婷丁香| 麻豆精品久久久久久蜜桃| 国产精品一二三区在线看| 亚洲久久久久久中文字幕| 看黄色毛片网站| 免费电影在线观看免费观看| 亚洲不卡免费看| 人人妻人人澡欧美一区二区| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 美女xxoo啪啪120秒动态图| 精品人妻一区二区三区麻豆| 久久99蜜桃精品久久| 一级毛片我不卡| 在线天堂最新版资源| 久久精品久久久久久噜噜老黄 | 国内精品久久久久精免费| 草草在线视频免费看| 联通29元200g的流量卡| 久久人人爽人人爽人人片va| 国产伦精品一区二区三区四那| 一边亲一边摸免费视频| 丝袜喷水一区| 91在线精品国自产拍蜜月| 国内精品美女久久久久久| 在线观看午夜福利视频| 日韩人妻高清精品专区| 亚洲精品久久国产高清桃花| 少妇高潮的动态图| 高清毛片免费观看视频网站| 男女啪啪激烈高潮av片| 99国产精品一区二区蜜桃av| 国产爱豆传媒在线观看| 成人性生交大片免费视频hd| 亚州av有码| 国产黄色视频一区二区在线观看 | 午夜福利在线观看免费完整高清在 | 91久久精品国产一区二区三区| 韩国av在线不卡| 久久精品国产亚洲av涩爱 | 中文精品一卡2卡3卡4更新| 欧美xxxx性猛交bbbb| 国语自产精品视频在线第100页| 好男人视频免费观看在线| 国产伦精品一区二区三区四那| 亚洲欧美日韩无卡精品| 极品教师在线视频| 寂寞人妻少妇视频99o| 99久久精品国产国产毛片| 寂寞人妻少妇视频99o| a级一级毛片免费在线观看| 男人狂女人下面高潮的视频| 久久九九热精品免费| 亚洲内射少妇av| 久久久久久久久大av| 国产午夜精品论理片| 亚洲欧美精品专区久久| 国产成人一区二区在线| 中出人妻视频一区二区| 最后的刺客免费高清国语| 亚洲熟妇中文字幕五十中出| 久久精品国产自在天天线| 天堂中文最新版在线下载 | 国产黄片美女视频| 亚洲婷婷狠狠爱综合网| 九九在线视频观看精品| 老女人水多毛片| 亚洲一区高清亚洲精品| 亚洲精品国产av成人精品| 国国产精品蜜臀av免费| 亚洲,欧美,日韩| 国产精品伦人一区二区| 亚洲精品国产av成人精品| 噜噜噜噜噜久久久久久91| av免费观看日本| 高清午夜精品一区二区三区 | 天堂网av新在线| 中文字幕精品亚洲无线码一区| 精品免费久久久久久久清纯| 国产精品不卡视频一区二区| 蜜桃亚洲精品一区二区三区| 人人妻人人澡欧美一区二区| 精品久久久噜噜| 狠狠狠狠99中文字幕| 少妇人妻一区二区三区视频| 中文字幕av在线有码专区| 中文字幕制服av| 麻豆国产av国片精品| 我要搜黄色片| 国产黄a三级三级三级人| 久久久久久九九精品二区国产| 网址你懂的国产日韩在线| 亚洲三级黄色毛片| 欧美在线一区亚洲| av福利片在线观看| av天堂中文字幕网| 久久久久久久久久久免费av| 九草在线视频观看| 精品久久久久久久人妻蜜臀av| 亚洲精品乱码久久久v下载方式| 国产视频内射| 成人av在线播放网站| 变态另类丝袜制服| 久久精品国产亚洲av香蕉五月| 亚洲国产精品合色在线| 国产av麻豆久久久久久久| 男女那种视频在线观看| 久久久久久久亚洲中文字幕| 好男人在线观看高清免费视频| 欧美在线一区亚洲| 色综合站精品国产| 国产私拍福利视频在线观看| av在线观看视频网站免费| 国产黄片美女视频| 成人永久免费在线观看视频| 级片在线观看| 极品教师在线视频| 老司机影院成人| 午夜亚洲福利在线播放| 精品午夜福利在线看| 美女黄网站色视频| 亚洲自拍偷在线| 全区人妻精品视频| 国产白丝娇喘喷水9色精品| 老熟妇乱子伦视频在线观看| 日韩欧美在线乱码| 99热精品在线国产| 亚洲欧美日韩高清在线视频| 又爽又黄a免费视频| 干丝袜人妻中文字幕| 国产人妻一区二区三区在| 99国产精品一区二区蜜桃av| 欧美zozozo另类| 在线免费十八禁| 男人狂女人下面高潮的视频| 99热这里只有精品一区| 国产日韩欧美在线精品| av在线亚洲专区| 麻豆久久精品国产亚洲av| 中文欧美无线码| 日本与韩国留学比较| 国产女主播在线喷水免费视频网站 | 身体一侧抽搐| 成人亚洲欧美一区二区av| 亚洲成人av在线免费| 大香蕉久久网| 99热这里只有精品一区| 又爽又黄a免费视频| 天堂影院成人在线观看| 国产精品永久免费网站| 精品人妻一区二区三区麻豆| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一二三区在线看| 国产亚洲精品久久久com| 99riav亚洲国产免费| 在线观看午夜福利视频| 91在线精品国自产拍蜜月| 97在线视频观看| 国产精品人妻久久久久久| 大型黄色视频在线免费观看| 中文亚洲av片在线观看爽| 中文亚洲av片在线观看爽| 我的老师免费观看完整版| 免费av不卡在线播放| 精品熟女少妇av免费看| 1000部很黄的大片| 亚洲精品乱码久久久久久按摩| av在线蜜桃| 国产av在哪里看| 久久国产乱子免费精品| 少妇丰满av| 99久久中文字幕三级久久日本| 精品国产三级普通话版| 26uuu在线亚洲综合色| 国产蜜桃级精品一区二区三区| 亚洲av男天堂| 爱豆传媒免费全集在线观看| www日本黄色视频网| 久久精品人妻少妇| 国产午夜精品一二区理论片| 精品国内亚洲2022精品成人| 欧美高清成人免费视频www| 国产午夜精品一二区理论片| 久久亚洲国产成人精品v| 国产白丝娇喘喷水9色精品| 两个人的视频大全免费| 2022亚洲国产成人精品| 一区福利在线观看| 少妇丰满av| 国产精品一及| 亚洲精品粉嫩美女一区| 国产精品一及| kizo精华| 国产精品一区二区三区四区免费观看| 日韩一本色道免费dvd| 日韩制服骚丝袜av| 一级毛片久久久久久久久女| 国产精品一区www在线观看| 国产综合懂色| 国产成人freesex在线| 久久精品人妻少妇| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品久久久久久成人av| 国内精品久久久久精免费| 免费无遮挡裸体视频| 天堂中文最新版在线下载 | 日本黄大片高清| 国产精品av视频在线免费观看| 欧美性猛交╳xxx乱大交人| 国产精品久久久久久精品电影小说 | 日韩一本色道免费dvd| 欧美日韩精品成人综合77777| 成人漫画全彩无遮挡| 亚洲欧美中文字幕日韩二区| 在线播放国产精品三级| 国产成人freesex在线| 亚洲精品亚洲一区二区| 久久精品国产99精品国产亚洲性色| 亚洲不卡免费看| 少妇熟女aⅴ在线视频| 国产爱豆传媒在线观看| 亚洲一级一片aⅴ在线观看| 午夜亚洲福利在线播放| 精品无人区乱码1区二区| 成人欧美大片| 免费av不卡在线播放| 哪个播放器可以免费观看大片| 色综合站精品国产| 在线免费十八禁| 国产又黄又爽又无遮挡在线| 欧美不卡视频在线免费观看| 亚洲精品乱码久久久久久按摩| 五月玫瑰六月丁香| 国产成人a∨麻豆精品| 女的被弄到高潮叫床怎么办| 97超碰精品成人国产| 国产午夜福利久久久久久| 久久99精品国语久久久| 97超视频在线观看视频| 成人无遮挡网站| 亚洲av电影不卡..在线观看| 日韩国内少妇激情av| 欧美+亚洲+日韩+国产| 床上黄色一级片| 午夜福利在线观看吧| 成人亚洲精品av一区二区| 18+在线观看网站| 精品人妻偷拍中文字幕| 国产男人的电影天堂91| 一级毛片我不卡| 在线播放国产精品三级| 成人高潮视频无遮挡免费网站| 99国产极品粉嫩在线观看| 狂野欧美激情性xxxx在线观看| 日韩亚洲欧美综合| 国产91av在线免费观看| 色5月婷婷丁香| 午夜福利在线在线| 亚洲欧洲日产国产| 午夜老司机福利剧场| 干丝袜人妻中文字幕| 麻豆成人av视频| 自拍偷自拍亚洲精品老妇| 狂野欧美激情性xxxx在线观看| 老司机影院成人| 国产成人a区在线观看| 一级黄色大片毛片| 久久久成人免费电影| 91在线精品国自产拍蜜月| 国产精品蜜桃在线观看 | 亚洲真实伦在线观看| 亚洲av一区综合| 夜夜夜夜夜久久久久| 熟妇人妻久久中文字幕3abv| 深夜a级毛片| 日本爱情动作片www.在线观看| 人妻久久中文字幕网| 亚洲中文字幕日韩| 国产精品免费一区二区三区在线| 免费看美女性在线毛片视频| 婷婷色av中文字幕| 人妻制服诱惑在线中文字幕| 美女大奶头视频| 免费看美女性在线毛片视频| 国产麻豆成人av免费视频| 免费一级毛片在线播放高清视频| 日韩欧美 国产精品| av又黄又爽大尺度在线免费看 | 亚洲欧洲国产日韩| 少妇猛男粗大的猛烈进出视频 | 亚洲国产日韩欧美精品在线观看| 自拍偷自拍亚洲精品老妇| 一个人免费在线观看电影| 国产v大片淫在线免费观看| 白带黄色成豆腐渣| 色综合亚洲欧美另类图片| 亚洲av熟女| 色尼玛亚洲综合影院| 精品熟女少妇av免费看| 国产成人精品婷婷| 精品一区二区三区视频在线| 看黄色毛片网站| 97人妻精品一区二区三区麻豆| 欧美激情在线99| 国产精品久久电影中文字幕| 变态另类成人亚洲欧美熟女| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久末码| 亚洲av成人精品一区久久| 国产女主播在线喷水免费视频网站 | 午夜激情欧美在线| 大香蕉久久网| 国产高清视频在线观看网站| 麻豆乱淫一区二区| 黄色一级大片看看| 久久久久久久久久黄片| 久久精品国产亚洲网站| 日韩一区二区视频免费看| 免费看a级黄色片| 99国产精品一区二区蜜桃av| 国产 一区 欧美 日韩| 亚洲精品自拍成人| 日本爱情动作片www.在线观看| 尤物成人国产欧美一区二区三区| 成年版毛片免费区| 亚洲一级一片aⅴ在线观看| 成人性生交大片免费视频hd| 亚洲在久久综合| 久久99热这里只有精品18| 男人的好看免费观看在线视频| 国产欧美日韩精品一区二区| 日日啪夜夜撸| 国产黄片视频在线免费观看| 国产高清三级在线| 男女啪啪激烈高潮av片| 亚洲av第一区精品v没综合| 亚洲精品乱码久久久v下载方式| 欧美最新免费一区二区三区| 亚洲精品色激情综合| 久99久视频精品免费| 国产精品野战在线观看| 久久热精品热| 春色校园在线视频观看| 免费一级毛片在线播放高清视频| 久久久午夜欧美精品| 免费搜索国产男女视频| 精品久久久久久久久久免费视频| 久久99蜜桃精品久久| 午夜老司机福利剧场| 国产成人午夜福利电影在线观看| 国产极品天堂在线| 亚洲高清免费不卡视频| 小蜜桃在线观看免费完整版高清| 91av网一区二区| 18禁裸乳无遮挡免费网站照片| 国产伦在线观看视频一区| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品女同一区二区软件| 高清毛片免费看| 天天躁夜夜躁狠狠久久av| 丝袜美腿在线中文| 国产 一区 欧美 日韩| 此物有八面人人有两片| 亚洲aⅴ乱码一区二区在线播放| 国产精品电影一区二区三区| 国产 一区精品| 毛片女人毛片| 亚洲成av人片在线播放无| 日韩三级伦理在线观看| 久久99热6这里只有精品| av福利片在线观看| 97在线视频观看| 一本精品99久久精品77| 久久99蜜桃精品久久| 在线观看午夜福利视频| 成年av动漫网址| 在线播放国产精品三级| 麻豆国产av国片精品| 成人亚洲欧美一区二区av| 综合色av麻豆| 熟妇人妻久久中文字幕3abv| 18禁裸乳无遮挡免费网站照片| 亚洲欧美精品综合久久99| 99视频精品全部免费 在线| a级毛片a级免费在线| 中文在线观看免费www的网站| 国产极品天堂在线| 97热精品久久久久久| 欧美日本视频| 日韩欧美三级三区| 免费看av在线观看网站| 亚洲电影在线观看av| 黄色视频,在线免费观看| 国产精品麻豆人妻色哟哟久久 | 久久这里有精品视频免费| 免费观看人在逋| 欧美极品一区二区三区四区| 日韩av不卡免费在线播放| 91在线精品国自产拍蜜月| 可以在线观看毛片的网站| 五月伊人婷婷丁香| 日韩 亚洲 欧美在线| 日韩,欧美,国产一区二区三区 | 综合色av麻豆| 别揉我奶头 嗯啊视频| 日产精品乱码卡一卡2卡三| 亚洲婷婷狠狠爱综合网| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| 日本一二三区视频观看| 亚洲性久久影院| 亚洲18禁久久av| 校园春色视频在线观看| 中文亚洲av片在线观看爽| 国产精品爽爽va在线观看网站| 直男gayav资源| 欧美+亚洲+日韩+国产| 久久午夜福利片| kizo精华| 我的女老师完整版在线观看| 久久99热6这里只有精品| 99久久久亚洲精品蜜臀av| 天天躁日日操中文字幕| 国产成人精品一,二区 | 精品日产1卡2卡| 日本免费一区二区三区高清不卡| 色播亚洲综合网| 亚洲成av人片在线播放无| 五月玫瑰六月丁香| 国产精品久久久久久久久免| 一进一出抽搐gif免费好疼| 亚洲精品影视一区二区三区av| 中文资源天堂在线| 99热只有精品国产| 午夜福利成人在线免费观看| 又粗又爽又猛毛片免费看| 午夜福利视频1000在线观看| 欧美高清成人免费视频www| 亚洲av免费在线观看| 国产v大片淫在线免费观看| 亚洲精品亚洲一区二区| 亚洲精品久久国产高清桃花| 国产精品乱码一区二三区的特点| 啦啦啦啦在线视频资源| 在线观看美女被高潮喷水网站| 国产成人精品一,二区 | 国产单亲对白刺激| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 国产成人a区在线观看| 色噜噜av男人的天堂激情| 男人和女人高潮做爰伦理| 老熟妇乱子伦视频在线观看| 久久久久久久午夜电影| 久久久a久久爽久久v久久| 亚洲精品成人久久久久久| 最近的中文字幕免费完整| 亚洲第一电影网av| 亚洲精品乱码久久久久久按摩| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 久久精品综合一区二区三区| 久久久久久久亚洲中文字幕| 少妇被粗大猛烈的视频| 熟女人妻精品中文字幕| 大香蕉久久网| 我要看日韩黄色一级片| 成人美女网站在线观看视频| 男女做爰动态图高潮gif福利片| 非洲黑人性xxxx精品又粗又长| 插逼视频在线观看| 天美传媒精品一区二区| 久久久久久久久久久丰满| 亚洲精品自拍成人| 午夜精品在线福利| 中国美白少妇内射xxxbb| 天天一区二区日本电影三级| 亚洲内射少妇av| 亚洲精品456在线播放app| 高清午夜精品一区二区三区 | 国产精品久久久久久精品电影小说 | 高清午夜精品一区二区三区 | 老女人水多毛片| 国产不卡一卡二| 最近2019中文字幕mv第一页| 国产精品综合久久久久久久免费| 夜夜爽天天搞| 欧美日本视频| a级毛片免费高清观看在线播放| 国产精品电影一区二区三区| 美女黄网站色视频| 国产在视频线在精品| 精品人妻视频免费看| 热99re8久久精品国产| 中文亚洲av片在线观看爽| av在线播放精品| kizo精华| 麻豆国产97在线/欧美| 日本一本二区三区精品| 91午夜精品亚洲一区二区三区| 女人十人毛片免费观看3o分钟| 尾随美女入室| 色综合站精品国产| 国产亚洲av片在线观看秒播厂 | 亚洲精品国产av成人精品| 国内精品一区二区在线观看| 久久综合国产亚洲精品| 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 亚洲国产精品久久男人天堂| 精品人妻视频免费看| 麻豆乱淫一区二区| 国产高潮美女av| 久久久色成人| 三级男女做爰猛烈吃奶摸视频| 欧美一区二区亚洲| 美女被艹到高潮喷水动态| 久久人人爽人人爽人人片va| 亚洲国产精品国产精品| 深夜a级毛片| 在线免费观看不下载黄p国产| 免费观看在线日韩| 嫩草影院入口| 黄色日韩在线| 欧美zozozo另类| 看十八女毛片水多多多| 国产黄片视频在线免费观看| 中文资源天堂在线| 午夜亚洲福利在线播放| 中文资源天堂在线| 亚洲av一区综合| 色5月婷婷丁香| 乱人视频在线观看| 看片在线看免费视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产亚洲av嫩草精品影院| 国语自产精品视频在线第100页| 日本爱情动作片www.在线观看| 国产激情偷乱视频一区二区| 国产黄a三级三级三级人| 最近最新中文字幕大全电影3| 精品久久久久久久久亚洲| 成人二区视频| 久久99精品国语久久久| 久久久国产成人免费| 黄色一级大片看看| 一级二级三级毛片免费看| 天堂av国产一区二区熟女人妻| 亚洲精品粉嫩美女一区| 亚洲国产精品国产精品| avwww免费| 内地一区二区视频在线| 黄片无遮挡物在线观看| 热99re8久久精品国产| 日韩欧美 国产精品| 精品一区二区三区视频在线| 日本成人三级电影网站| 国产黄片视频在线免费观看| 伊人久久精品亚洲午夜| 免费看美女性在线毛片视频| 国产精品人妻久久久久久| 成人影院久久| 91成人精品电影| 国产精品国产三级国产av玫瑰| 国产亚洲精品第一综合不卡 | 成人国产av品久久久| 亚洲欧美成人综合另类久久久| 久久久久久久国产电影| 亚洲伊人久久精品综合| 97在线人人人人妻| 99久久中文字幕三级久久日本| 涩涩av久久男人的天堂| 色吧在线观看| 新久久久久国产一级毛片| 日本wwww免费看| tube8黄色片| 最近最新中文字幕免费大全7| 亚洲第一av免费看| 大码成人一级视频| 久热这里只有精品99| 97在线视频观看| 熟女av电影| xxxhd国产人妻xxx| 99re6热这里在线精品视频| 啦啦啦视频在线资源免费观看| 国产在线免费精品| av在线播放精品| 亚洲经典国产精华液单| 精品少妇内射三级| 亚洲经典国产精华液单| av国产久精品久网站免费入址| 免费看光身美女| 丝袜美足系列| 国产欧美亚洲国产| 欧美日韩国产mv在线观看视频| 在线天堂最新版资源| 大香蕉97超碰在线| 久久免费观看电影| 99热全是精品| 国产一区二区三区av在线| 香蕉精品网在线| 亚洲欧洲日产国产| 国产又色又爽无遮挡免| 在线观看国产h片| 国产欧美亚洲国产| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 好男人视频免费观看在线| 国产av码专区亚洲av| 午夜福利影视在线免费观看|