• <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)
    啦啦啦观看免费观看视频高清| av免费观看日本| 久久精品国产99精品国产亚洲性色| 99国产精品一区二区蜜桃av| 亚洲综合色惰| 午夜福利在线观看免费完整高清在| 成人一区二区视频在线观看| 色视频www国产| 午夜精品在线福利| 国产精品一区二区在线观看99 | 欧美性猛交黑人性爽| 精品国产一区二区三区久久久樱花 | 日韩精品有码人妻一区| 91狼人影院| 国产91av在线免费观看| 亚洲成人久久爱视频| 一本久久精品| 一夜夜www| 成年女人永久免费观看视频| 精品国产三级普通话版| 婷婷色麻豆天堂久久 | 成年女人看的毛片在线观看| 尾随美女入室| 久久久精品大字幕| 特大巨黑吊av在线直播| 一区二区三区免费毛片| 男女边吃奶边做爰视频| 伦精品一区二区三区| 水蜜桃什么品种好| 深爱激情五月婷婷| 五月玫瑰六月丁香| 婷婷色综合大香蕉| 日韩av在线免费看完整版不卡| 网址你懂的国产日韩在线| 中文字幕av在线有码专区| 午夜久久久久精精品| 干丝袜人妻中文字幕| 日韩av在线大香蕉| 午夜亚洲福利在线播放| 久久精品影院6| 亚洲成人精品中文字幕电影| 岛国毛片在线播放| 午夜视频国产福利| 亚洲av成人精品一区久久| 色尼玛亚洲综合影院| 99热这里只有是精品在线观看| 精品久久久久久久末码| 综合色av麻豆| 日韩大片免费观看网站 | 亚洲精品亚洲一区二区| 日本免费在线观看一区| 不卡视频在线观看欧美| 久久久久久久久中文| 日韩欧美国产在线观看| 久久精品国产亚洲av天美| 韩国av在线不卡| 日韩高清综合在线| 一个人免费在线观看电影| 日韩强制内射视频| 久久久久久久久大av| 日韩欧美精品v在线| 两个人的视频大全免费| 亚洲中文字幕一区二区三区有码在线看| 国产精品99久久久久久久久| 免费一级毛片在线播放高清视频| 国产单亲对白刺激| 天美传媒精品一区二区| 精品午夜福利在线看| 国产成人免费观看mmmm| 少妇丰满av| 日本与韩国留学比较| 亚洲精品乱码久久久久久按摩| 亚洲av中文字字幕乱码综合| 国产色婷婷99| 午夜精品国产一区二区电影 | .国产精品久久| 精品酒店卫生间| 看免费成人av毛片| 精品酒店卫生间| 老司机影院毛片| 亚洲,欧美,日韩| 亚洲成av人片在线播放无| 欧美潮喷喷水| 免费观看精品视频网站| 18禁动态无遮挡网站| 视频中文字幕在线观看| 久久综合国产亚洲精品| 一边亲一边摸免费视频| 国产精品一区二区三区四区久久| 亚洲av成人精品一二三区| 岛国毛片在线播放| 伊人久久精品亚洲午夜| 特大巨黑吊av在线直播| 精品人妻熟女av久视频| 亚洲五月天丁香| 亚洲精品乱久久久久久| 只有这里有精品99| 日本午夜av视频| 国产一区二区三区av在线| 欧美精品一区二区大全| 成人亚洲欧美一区二区av| 日韩国内少妇激情av| 日本一本二区三区精品| 国产精品麻豆人妻色哟哟久久 | 最近中文字幕高清免费大全6| 日日撸夜夜添| 亚洲av免费在线观看| 国产精品福利在线免费观看| 国产精品综合久久久久久久免费| 在线a可以看的网站| 最近中文字幕2019免费版| 免费无遮挡裸体视频| 人妻制服诱惑在线中文字幕| 午夜免费男女啪啪视频观看| 免费av观看视频| 午夜福利成人在线免费观看| av线在线观看网站| 级片在线观看| 国产成人精品一,二区| 国产精品野战在线观看| 国产免费男女视频| 精品熟女少妇av免费看| 蜜臀久久99精品久久宅男| 99久久精品热视频| 99在线人妻在线中文字幕| 免费看美女性在线毛片视频| 欧美xxxx黑人xx丫x性爽| 国产中年淑女户外野战色| 中文字幕制服av| 天堂av国产一区二区熟女人妻| 欧美日韩国产亚洲二区| 成人亚洲欧美一区二区av| 婷婷色麻豆天堂久久 | 亚洲精品,欧美精品| 欧美又色又爽又黄视频| av天堂中文字幕网| 免费av不卡在线播放| 伊人久久精品亚洲午夜| 黄色一级大片看看| 国产黄片美女视频| 精品国产露脸久久av麻豆 | 国模一区二区三区四区视频| 好男人视频免费观看在线| 午夜免费激情av| 国产淫片久久久久久久久| 99热这里只有是精品在线观看| 成人二区视频| 国产精品人妻久久久影院| 床上黄色一级片| 欧美激情在线99| 亚洲国产精品专区欧美| 午夜福利网站1000一区二区三区| 亚洲激情五月婷婷啪啪| 日本黄色视频三级网站网址| 小说图片视频综合网站| 亚洲中文字幕日韩| 男的添女的下面高潮视频| 男的添女的下面高潮视频| 国产三级在线视频| 久久这里只有精品中国| 午夜a级毛片| 亚洲av福利一区| 免费看美女性在线毛片视频| 卡戴珊不雅视频在线播放| 一本久久精品| 成人无遮挡网站| 亚洲av熟女| 亚洲最大成人手机在线| 麻豆乱淫一区二区| 老司机福利观看| 乱码一卡2卡4卡精品| 日韩欧美在线乱码| 亚洲精华国产精华液的使用体验| 人人妻人人澡欧美一区二区| 亚洲中文字幕一区二区三区有码在线看| 国产乱人视频| 美女大奶头视频| 高清av免费在线| 日韩成人av中文字幕在线观看| 午夜福利在线观看吧| 午夜激情福利司机影院| 亚洲人成网站在线播| 一级爰片在线观看| 国产成人午夜福利电影在线观看| 蜜桃亚洲精品一区二区三区| 九九爱精品视频在线观看| videos熟女内射| 日本黄色片子视频| 好男人在线观看高清免费视频| 欧美成人a在线观看| 国产黄片视频在线免费观看| 麻豆国产97在线/欧美| 成人漫画全彩无遮挡| videossex国产| 亚洲成人久久爱视频| 级片在线观看| 高清视频免费观看一区二区 | av专区在线播放| 国产老妇女一区| 久久精品久久久久久噜噜老黄 | 国产精品电影一区二区三区| .国产精品久久| 91aial.com中文字幕在线观看| 免费电影在线观看免费观看| 视频中文字幕在线观看| 六月丁香七月| 91精品国产九色| 在线播放无遮挡| 一区二区三区高清视频在线| 欧美日韩精品成人综合77777| 色综合站精品国产| 深爱激情五月婷婷| 国产精品熟女久久久久浪| 日本黄色片子视频| 十八禁国产超污无遮挡网站| 精品少妇黑人巨大在线播放 | 色哟哟·www| 国内揄拍国产精品人妻在线| 久久婷婷人人爽人人干人人爱| 嫩草影院精品99| 欧美色视频一区免费| 国产 一区 欧美 日韩| 国产精品.久久久| 国产乱人视频| 亚洲精品色激情综合| 亚洲国产精品成人综合色| 晚上一个人看的免费电影| 国产高清三级在线| 日韩欧美在线乱码| 国产伦在线观看视频一区| 国产精品久久久久久久电影| 日日啪夜夜撸| 亚洲精品乱码久久久久久按摩| 天天躁夜夜躁狠狠久久av| 欧美日韩精品成人综合77777| 菩萨蛮人人尽说江南好唐韦庄 | 听说在线观看完整版免费高清| 少妇猛男粗大的猛烈进出视频 | 在线播放无遮挡| 国产成人freesex在线| 午夜福利在线观看免费完整高清在| 国产精品国产高清国产av| 成人高潮视频无遮挡免费网站| 美女被艹到高潮喷水动态| 国产亚洲精品av在线| 天堂√8在线中文| 久久久精品94久久精品| 七月丁香在线播放| 最近的中文字幕免费完整| 国产精品99久久久久久久久| 亚洲av中文av极速乱| 岛国在线免费视频观看| 亚洲人成网站高清观看| 91aial.com中文字幕在线观看| 丝袜美腿在线中文| 欧美日本视频| 亚洲欧洲日产国产| av黄色大香蕉| 国产av在哪里看| 春色校园在线视频观看| 亚洲精品一区蜜桃| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| 国产高清有码在线观看视频| 超碰97精品在线观看| 国产精品野战在线观看| 国产熟女欧美一区二区| 国产三级在线视频| 国产成年人精品一区二区| 综合色丁香网| 一本一本综合久久| 老女人水多毛片| 日韩国内少妇激情av| 五月玫瑰六月丁香| 国产在视频线在精品| 久久99蜜桃精品久久| 精品一区二区三区视频在线| 精品国产露脸久久av麻豆 | 国产一区有黄有色的免费视频 | 精品一区二区免费观看| 成人二区视频| 婷婷色麻豆天堂久久 | 91久久精品电影网| 春色校园在线视频观看| 99视频精品全部免费 在线| 我要搜黄色片| 天堂av国产一区二区熟女人妻| 久久久成人免费电影| 成年免费大片在线观看| 少妇高潮的动态图| 午夜福利在线在线| 淫秽高清视频在线观看| 丰满少妇做爰视频| 久久精品国产亚洲av天美| 午夜视频国产福利| 你懂的网址亚洲精品在线观看 | 中文乱码字字幕精品一区二区三区 | 久久99蜜桃精品久久| 亚洲最大成人手机在线| 欧美日韩综合久久久久久| 美女xxoo啪啪120秒动态图| 99热这里只有精品一区| 欧美潮喷喷水| 久久久久网色| 国产一级毛片在线| 国产午夜精品论理片| 国产午夜精品久久久久久一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 国产老妇女一区| 国产高清三级在线| 日本午夜av视频| 午夜爱爱视频在线播放| 免费观看在线日韩| 一区二区三区乱码不卡18| 国产成人freesex在线| 色综合站精品国产| 99热精品在线国产| 色噜噜av男人的天堂激情| 亚洲精品国产av成人精品| 日本猛色少妇xxxxx猛交久久| 久久韩国三级中文字幕| 亚洲精品一区蜜桃| 夫妻性生交免费视频一级片| 免费观看a级毛片全部| 91久久精品国产一区二区三区| 国产亚洲av嫩草精品影院| 久久精品国产亚洲av天美| 亚洲色图av天堂| 日韩三级伦理在线观看| 特大巨黑吊av在线直播| 久久久午夜欧美精品| 欧美色视频一区免费| 精品国产露脸久久av麻豆 | 在线免费观看的www视频| 中文字幕亚洲精品专区| 男女下面进入的视频免费午夜| 国产69精品久久久久777片| 日韩av在线免费看完整版不卡| 毛片女人毛片| 免费观看在线日韩| 久久久久九九精品影院| 国产熟女欧美一区二区| 大香蕉97超碰在线| 国产一级毛片七仙女欲春2| 国产精品久久视频播放| 日韩在线高清观看一区二区三区| 色综合站精品国产| 小蜜桃在线观看免费完整版高清| 噜噜噜噜噜久久久久久91| 日韩一区二区视频免费看| 亚洲国产欧美在线一区| 欧美97在线视频| 国产精品,欧美在线| 精品人妻一区二区三区麻豆| 国产成人免费观看mmmm| 日韩欧美国产在线观看| 国产不卡一卡二| 国产精品熟女久久久久浪| 免费看光身美女| 成人漫画全彩无遮挡| 中文字幕av在线有码专区| 少妇的逼好多水| 国产久久久一区二区三区| 麻豆久久精品国产亚洲av| 日韩av在线大香蕉| 看十八女毛片水多多多| 国产黄a三级三级三级人| 国产一级毛片在线| 日韩中字成人| 亚洲av电影不卡..在线观看| 午夜a级毛片| 亚洲欧美精品综合久久99| 村上凉子中文字幕在线| 又爽又黄无遮挡网站| 亚洲av中文av极速乱| 中文字幕精品亚洲无线码一区| 亚洲精品aⅴ在线观看| 好男人在线观看高清免费视频| 欧美三级亚洲精品| 日韩亚洲欧美综合| 亚洲一级一片aⅴ在线观看| 一级毛片我不卡| 亚洲av成人av| 视频中文字幕在线观看| 欧美成人免费av一区二区三区| 久久精品91蜜桃| 日韩av在线免费看完整版不卡| 日本与韩国留学比较| 免费观看人在逋| 99热这里只有是精品在线观看| 国产伦在线观看视频一区| 美女大奶头视频| 国产片特级美女逼逼视频| 免费一级毛片在线播放高清视频| www日本黄色视频网| 在现免费观看毛片| 亚洲欧美成人精品一区二区| 成人三级黄色视频| 99久久精品国产国产毛片| 亚洲成av人片在线播放无| 我的女老师完整版在线观看| 国产乱来视频区| 乱系列少妇在线播放| 91狼人影院| 国产乱人视频| 国产成人a∨麻豆精品| 大又大粗又爽又黄少妇毛片口| 国产乱人视频| 亚洲精品影视一区二区三区av| 欧美一级a爱片免费观看看| 自拍偷自拍亚洲精品老妇| 校园人妻丝袜中文字幕| 免费观看精品视频网站| 午夜福利网站1000一区二区三区| 搡老妇女老女人老熟妇| 国产又黄又爽又无遮挡在线| 亚洲精品久久久久久婷婷小说 | 中文字幕久久专区| 国产视频首页在线观看| 久久久久久久久久成人| 一个人观看的视频www高清免费观看| 国产白丝娇喘喷水9色精品| 国产精品久久久久久av不卡| 熟女电影av网| 亚洲欧美日韩高清专用| 久久久午夜欧美精品| 女人被狂操c到高潮| 哪个播放器可以免费观看大片| 中文字幕制服av| 亚洲av成人精品一二三区| 中文在线观看免费www的网站| 国产欧美另类精品又又久久亚洲欧美| 国产激情偷乱视频一区二区| 国产高清国产精品国产三级 | 六月丁香七月| 亚洲国产高清在线一区二区三| 国产老妇伦熟女老妇高清| 男女啪啪激烈高潮av片| 久久精品久久久久久久性| 少妇被粗大猛烈的视频| 国产乱人偷精品视频| 日韩欧美三级三区| 可以在线观看毛片的网站| a级毛片免费高清观看在线播放| 日日摸夜夜添夜夜添av毛片| 日日干狠狠操夜夜爽| 亚洲av二区三区四区| 禁无遮挡网站| 中文字幕av成人在线电影| 少妇人妻精品综合一区二区| 久久国产乱子免费精品| 天堂影院成人在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲精华国产精华液的使用体验| av在线播放精品| 一个人看的www免费观看视频| 一级毛片久久久久久久久女| 成人午夜精彩视频在线观看| 黄色一级大片看看| 全区人妻精品视频| 国产不卡一卡二| 一区二区三区四区激情视频| 天堂网av新在线| 欧美人与善性xxx| 久久精品夜夜夜夜夜久久蜜豆| 色尼玛亚洲综合影院| 国产成年人精品一区二区| 国产在线男女| 久久综合国产亚洲精品| 国产精品国产高清国产av| 日日摸夜夜添夜夜爱| 欧美区成人在线视频| 亚洲欧美日韩高清专用| 欧美日本亚洲视频在线播放| 男人狂女人下面高潮的视频| 久久久成人免费电影| 亚洲三级黄色毛片| 久久综合国产亚洲精品| 一级黄色大片毛片| 国产精品美女特级片免费视频播放器| 91aial.com中文字幕在线观看| 欧美成人午夜免费资源| 黑人高潮一二区| 成人午夜精彩视频在线观看| 综合色丁香网| 99视频精品全部免费 在线| 国产真实乱freesex| 国产精品福利在线免费观看| 99久久人妻综合| 国产精品国产三级国产av玫瑰| 99热精品在线国产| 日本av手机在线免费观看| 国产黄片美女视频| 老司机影院成人| 国产亚洲5aaaaa淫片| 一卡2卡三卡四卡精品乱码亚洲| 色综合亚洲欧美另类图片| 久久久精品大字幕| 一级毛片久久久久久久久女| 精品国产一区二区三区久久久樱花 | 国产成人a∨麻豆精品| 寂寞人妻少妇视频99o| 国产片特级美女逼逼视频| 亚洲精品,欧美精品| 国产久久久一区二区三区| 搡女人真爽免费视频火全软件| 国产日韩欧美在线精品| 免费看美女性在线毛片视频| 精品国产露脸久久av麻豆 | 亚洲欧美精品自产自拍| 99在线人妻在线中文字幕| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 亚洲国产色片| 国产成人精品久久久久久| 亚洲综合精品二区| 午夜久久久久精精品| 欧美+日韩+精品| 日日撸夜夜添| 91精品国产九色| 2021天堂中文幕一二区在线观| 欧美色视频一区免费| h日本视频在线播放| 99九九线精品视频在线观看视频| 久久精品熟女亚洲av麻豆精品 | 国产伦理片在线播放av一区| 精品午夜福利在线看| 成人美女网站在线观看视频| 国产男人的电影天堂91| 国产大屁股一区二区在线视频| av免费在线看不卡| 日本三级黄在线观看| 特级一级黄色大片| 麻豆一二三区av精品| 长腿黑丝高跟| 秋霞伦理黄片| 日韩强制内射视频| 99久国产av精品| 2022亚洲国产成人精品| 热99re8久久精品国产| 天天躁夜夜躁狠狠久久av| a级毛色黄片| 免费看a级黄色片| 欧美成人精品欧美一级黄| 精品久久久久久久久亚洲| 18禁动态无遮挡网站| 在线观看一区二区三区| 日本三级黄在线观看| 亚洲av免费在线观看| 久久久久久久久中文| 噜噜噜噜噜久久久久久91| 两性午夜刺激爽爽歪歪视频在线观看| 综合色丁香网| 国产激情偷乱视频一区二区| 99久久精品热视频| 麻豆av噜噜一区二区三区| av天堂中文字幕网| 青春草亚洲视频在线观看| 国产久久久一区二区三区| 亚洲激情五月婷婷啪啪| 91精品伊人久久大香线蕉| 亚洲av中文字字幕乱码综合| 久久久成人免费电影| 啦啦啦啦在线视频资源| 亚洲第一区二区三区不卡| 91久久精品国产一区二区成人| 尤物成人国产欧美一区二区三区| 国产精品国产三级国产av玫瑰| 国产探花在线观看一区二区| 丰满乱子伦码专区| 亚洲在久久综合| 国产精品久久久久久av不卡| 天堂中文最新版在线下载 | 高清视频免费观看一区二区 | 日本欧美国产在线视频| 精品人妻偷拍中文字幕| 韩国av在线不卡| 日本色播在线视频| 尤物成人国产欧美一区二区三区| 天堂√8在线中文| 中文资源天堂在线| 一夜夜www| 夫妻性生交免费视频一级片| 日韩亚洲欧美综合| av在线观看视频网站免费| 天堂影院成人在线观看| 日韩成人av中文字幕在线观看| 日韩欧美 国产精品| 久久久久久久久久久丰满| 偷拍熟女少妇极品色| 日本-黄色视频高清免费观看| 嘟嘟电影网在线观看| 国产黄a三级三级三级人| 成人亚洲欧美一区二区av| 国产精品国产三级国产专区5o | 日韩在线高清观看一区二区三区| 男插女下体视频免费在线播放| 亚洲乱码一区二区免费版| 99热精品在线国产| 一边摸一边抽搐一进一小说| 欧美日本亚洲视频在线播放| 亚洲三级黄色毛片| av国产久精品久网站免费入址| 成人毛片60女人毛片免费| 久久精品国产亚洲网站| 国产高清三级在线| 精品久久久久久久人妻蜜臀av| 久久综合国产亚洲精品| 少妇高潮的动态图| 国内精品宾馆在线| 国内少妇人妻偷人精品xxx网站| 国产精品伦人一区二区| 精品久久久久久电影网 | av在线观看视频网站免费| 久久久久久久国产电影|