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

    廠際氫氣網(wǎng)絡(luò)多周期集成的分步優(yōu)化方法

    2021-10-04 15:11:06蔣迎花韓儒松康麗霞劉永忠
    化工學(xué)報(bào) 2021年9期
    關(guān)鍵詞:流股公用網(wǎng)絡(luò)結(jié)構(gòu)

    蔣迎花,韓儒松,康麗霞,2,劉永忠,2,3

    (1 西安交通大學(xué)化工系,陜西西安 710049;2 陜西省能源化工過程強(qiáng)化重點(diǎn)實(shí)驗(yàn)室,陜西西安 710049;3 熱流科學(xué)與工程教育部重點(diǎn)實(shí)驗(yàn)室,陜西西安 710049)

    引言

    在能源化工產(chǎn)業(yè)中,氫氣既是煤化工和石油化工的副產(chǎn)品,也是現(xiàn)代煤化工廠和煉廠生產(chǎn)潔凈燃料的重要原料。在一定區(qū)域或化工園區(qū)內(nèi),能源化工企業(yè)內(nèi)部一般均建有獨(dú)立的制氫裝置和氫氣網(wǎng)絡(luò)系統(tǒng)。對(duì)于不同類型的能源化工企業(yè),由于生產(chǎn)工藝、生產(chǎn)產(chǎn)品和生產(chǎn)周期不同,對(duì)氫氣的需求也具有不同的特點(diǎn),這將為區(qū)域或園區(qū)內(nèi)能源化工企業(yè)之間共享和協(xié)調(diào)氫氣資源提供可能。采用系統(tǒng)集成方法在區(qū)域或園區(qū)內(nèi)協(xié)調(diào)優(yōu)化氫氣的供給和需求對(duì)于氫氣資源的合理配置和生產(chǎn)成本的整體降低具有重要意義[1-5]。在企業(yè)的實(shí)際生產(chǎn)中,由于氫氣網(wǎng)絡(luò)操作參數(shù)(如流率、壓力和純度等)受原料性質(zhì)、工藝變化、催化劑活性衰減、季節(jié)更替、產(chǎn)品市場(chǎng)價(jià)格等內(nèi)外部因素的影響,考慮變化的操作參數(shù)對(duì)氫氣網(wǎng)絡(luò)的優(yōu)化設(shè)計(jì)非常重要[6-10]。在工業(yè)園區(qū)中,各企業(yè)氫氣網(wǎng)絡(luò)的操作參數(shù)在時(shí)間上的變化特征并不相同。因此,為了有效降低區(qū)域內(nèi)氫氣消耗和制氫能耗,有必要研究和開發(fā)系統(tǒng)優(yōu)化方法和工具,對(duì)區(qū)域內(nèi)各企業(yè)的氫氣網(wǎng)絡(luò)系統(tǒng)從時(shí)間和空間的維度上進(jìn)行整體協(xié)調(diào)優(yōu)化,以達(dá)到區(qū)域內(nèi)氫氣網(wǎng)絡(luò)系統(tǒng)的協(xié)同最優(yōu)和安全操作[11-14]。

    在時(shí)間維度上,對(duì)于氫氣網(wǎng)絡(luò)的優(yōu)化設(shè)計(jì)可采用多周期方法考慮變化操作參數(shù)對(duì)系統(tǒng)優(yōu)化設(shè)計(jì)的影響[15-18]。對(duì)于多周期氫氣網(wǎng)絡(luò)系統(tǒng)的優(yōu)化設(shè)計(jì),通??刹捎猛椒ê头植椒▽?duì)氫氣網(wǎng)絡(luò)系統(tǒng)進(jìn)行優(yōu)化設(shè)計(jì)[19-20]。從設(shè)計(jì)模型全局出發(fā),同步法同時(shí)權(quán)衡了系統(tǒng)中各設(shè)計(jì)變量的影響,可得到氫氣網(wǎng)絡(luò)系統(tǒng)的最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)和氫氣調(diào)度策略[21]。然而,隨著氫氣網(wǎng)絡(luò)系統(tǒng)的設(shè)計(jì)模型規(guī)模擴(kuò)大和工況數(shù)量增加,優(yōu)化設(shè)計(jì)模型可能難于獲得可行解。對(duì)于分步法,雖然難以獲得氫氣網(wǎng)絡(luò)系統(tǒng)優(yōu)化設(shè)計(jì)的最優(yōu)解,但往往可以簡(jiǎn)化求解模型并提高求解效率[22]。因此,在實(shí)際的問題中,可采用分步法對(duì)大規(guī)模的多周期氫氣網(wǎng)絡(luò)優(yōu)化設(shè)計(jì)模型進(jìn)行求解。

    在空間維度上,廠際氫氣網(wǎng)絡(luò)系統(tǒng)優(yōu)化的分步集成方法得到廣泛應(yīng)用[23]。Deng 等[24]采用改進(jìn)的問題表法確定廠際氫氣網(wǎng)絡(luò)集成中氫氣公用工程用量目標(biāo)。該方法先確定每個(gè)廠的氫氣公用工程用量,然后通過廠際間氫氣回用進(jìn)一步減少氫氣公用工程用量。此外,Deng 等[25]還提出了改進(jìn)問題表法構(gòu)建廠際氫氣網(wǎng)絡(luò)集成的兩種策略,一是先進(jìn)行單廠集成,然后通過廢物流股將不同的廠連接起來,另一策略是將所有廠的流股當(dāng)作一個(gè)整體進(jìn)行集成。針對(duì)廠際氫氣網(wǎng)絡(luò)的優(yōu)化設(shè)計(jì),Kang 等[26]提出了優(yōu)化中間管道中氫氣的純度等級(jí)和壓力等級(jí)的分步優(yōu)化方法。Lou等[23]提出了兩步優(yōu)化法,先使用結(jié)合夾點(diǎn)的轉(zhuǎn)運(yùn)模型確定各廠的氫氣消耗,然后再分別優(yōu)化各廠的氫氣網(wǎng)絡(luò)結(jié)構(gòu)。因此,分步集成法大多是首先優(yōu)化氫氣網(wǎng)絡(luò)中的氫氣公用工程消耗量,然后再優(yōu)化氫氣系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)。雖然上述方法可有效求解廠際氫氣網(wǎng)絡(luò)的優(yōu)化設(shè)計(jì)問題,但大多數(shù)分步集成方法并未考慮操作參數(shù)多周期變化下的廠際氫氣網(wǎng)絡(luò)系統(tǒng)集成,也缺乏針對(duì)大規(guī)模多周期的廠際氫氣網(wǎng)絡(luò)系統(tǒng)優(yōu)化設(shè)計(jì)的高效方法。

    本文提出了一種多周期廠際氫氣網(wǎng)絡(luò)優(yōu)化設(shè)計(jì)的分步求解方法。與現(xiàn)有的廠際氫氣系統(tǒng)集成方法不同,本文方法分三步實(shí)現(xiàn)多周期廠際氫氣系統(tǒng)的優(yōu)化設(shè)計(jì),即是首先通過求解單周期廠際氫氣網(wǎng)絡(luò)優(yōu)化模型確定各子周期中廠際氫氣傳輸量,然后采用多周期單廠優(yōu)化設(shè)計(jì)模型確定各廠內(nèi)的氫氣系統(tǒng)結(jié)構(gòu),再采用多周期廠際氫氣系統(tǒng)優(yōu)化模型確定廠際間的氫氣流股匹配結(jié)構(gòu),最終獲取整個(gè)氫氣系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)和氫氣調(diào)度方法。

    1 問題描述

    圖1給出工業(yè)園區(qū)中廠際氫氣網(wǎng)絡(luò)系統(tǒng)的多周期集成設(shè)計(jì)概念圖,其中圖1(a)表示工業(yè)園區(qū)中多廠氫氣流率的動(dòng)態(tài)變化特性,圖1(b)表示廠際氫氣網(wǎng)絡(luò)多周期集成的超結(jié)構(gòu)。在工業(yè)園區(qū)有N個(gè)涉及氫氣供需的工廠,每個(gè)廠均包含氫氣公用工程(HU)、氫源(HS)、氫阱(HK)和燃?xì)庀到y(tǒng)(HF)。這些廠通過共享氫氣公用工程實(shí)現(xiàn)氫氣的集中供應(yīng),廠內(nèi)則采用氫氣管道和提純器實(shí)現(xiàn)連接,系統(tǒng)中的壓縮機(jī)用于提升流股的壓力,而棄用的過程氫氣流股排放至燃?xì)庀到y(tǒng)。

    圖1 廠際氫氣網(wǎng)絡(luò)多周期集成設(shè)計(jì)Fig.1 Integration of an inter-plant hydrogen network under multi-period operation

    考慮變化操作參數(shù)的影響,可將廠際氫氣網(wǎng)絡(luò)集成問題轉(zhuǎn)化為一個(gè)廠際氫氣網(wǎng)絡(luò)多周期優(yōu)化設(shè)計(jì)問題。本文針對(duì)該廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)問題可描述為:(1)已知各廠的氫氣流股操作參數(shù)和(2)氫氣公用工程、電力和設(shè)備等費(fèi)用參數(shù)。目標(biāo)是通過廠間和廠內(nèi)氫氣系統(tǒng)的集成,獲取多周期操作下廠際氫氣網(wǎng)絡(luò)的最優(yōu)系統(tǒng)結(jié)構(gòu)和氫氣調(diào)度策略,使得整個(gè)工業(yè)園區(qū)氫氣系統(tǒng)的年總成本費(fèi)用最低。

    為了描述廠際多周期氫氣系統(tǒng)的結(jié)構(gòu),定義集合如下:

    N={n|n=1,2,…,N}表示工業(yè)園區(qū)的用氫工廠集合;HK={k|k=1,2,…,K}表示氫阱集合;HU表示氫氣公用工程集合;HS 表示氫源集合;US=HU ∪HS={i|i=1,2,…,I}表示整個(gè)系統(tǒng)中的氫源集合;PU={m|m=1,2,…,M}表示系統(tǒng)中提純器集合;P={p|p=1,2,…,P}表示周期數(shù)集合。

    除集合N 外,以上集合均可以寫成由N個(gè)子集的并集,例如,集合US=US1∪US2∪…∪USN。

    2 廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)模型

    對(duì)于廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)問題,其特點(diǎn)是氫氣系統(tǒng)涉及多廠協(xié)調(diào)和多周期共享等,相應(yīng)的數(shù)學(xué)規(guī)劃模型規(guī)模隨廠數(shù)和周期數(shù)增大;在氫氣網(wǎng)絡(luò)系統(tǒng)的優(yōu)化中,操作費(fèi)用占總費(fèi)用的比例較顯著?;谶@些特點(diǎn),本文提出一種分步方法求解廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)問題,該方法的實(shí)施流程如圖2 所示,圖中F 表示氫氣公用工程傳輸量,S表示氫氣系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)。

    圖2 廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)模型分步求解流程Fig.2 A schematic diagram of the three-step method to integrate an inter-plant hydrogen network under multi-period operations

    該方法由三步構(gòu)成:(1)Step 1:求解單周期廠際氫氣網(wǎng)絡(luò)模型,獲取每個(gè)公用工程跨廠傳輸?shù)募儦淞?;?)Step 2:采用多周期單廠氫氣網(wǎng)絡(luò)優(yōu)化模型,在(1)的基礎(chǔ)上確定廠內(nèi)氫氣系統(tǒng)的結(jié)構(gòu)、流量及設(shè)備參數(shù);(3)Step 3:求解廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)模型確定廠際氫氣網(wǎng)絡(luò)的結(jié)構(gòu),最終獲得工業(yè)園區(qū)整個(gè)氫氣網(wǎng)絡(luò)的系統(tǒng)結(jié)構(gòu)和氫氣調(diào)度策略。在以上三步中均以年度化費(fèi)用最低為優(yōu)化目標(biāo),分別用三組優(yōu)化模型(Model T、Model P 和Model Q)表示,各步具體實(shí)施過程分述如下。

    2.1 Step 1:廠際氫氣傳輸量的優(yōu)化

    廠際氫氣傳輸量可通過求解單周期廠際氫氣網(wǎng)絡(luò)優(yōu)化模型獲取。在本文中,單周期廠際氫氣網(wǎng)絡(luò)優(yōu)化模型表示為Model T。

    2.1.1 目標(biāo)函數(shù) 該模型以氫氣網(wǎng)絡(luò)的年度化費(fèi)用最小為優(yōu)化目標(biāo),表示為

    其中,t為子周期的時(shí)間長(zhǎng)度;ei為氫氣公用工程的單價(jià)。

    燃?xì)馐找婵杀硎緸?/p>

    氫氣管道投資成本可表示為

    其中,di,fuel為過程氫源i和燃?xì)庀到y(tǒng)之間的距離。二元變量zp,i,fuel用于判斷過程氫源到燃?xì)庀到y(tǒng)的連接是否存在。

    壓縮機(jī)的投資成本可表示為

    壓縮機(jī)的投資費(fèi)用表示為

    其 中,zp,a,b∈{zp,i,k,zp,i,m,zp,m,k}?i∈US,m∈PU,k∈HK,p∈P;acom和bcom為壓縮機(jī)的費(fèi)用系數(shù);為升壓二元變量,用來判斷某個(gè)匹配連接是否需要壓縮機(jī)來提升壓力。

    提純器的投資成本可表示為

    式中,zp,m用于判斷提純器是否存在的二元變量;apur和bpur為提純器的費(fèi)用系數(shù)。

    2.1.2 約束條件 該模型的約束條件包括:(1)氫源和氫阱流股約束;(2)提純器約束;(3)匹配連接約束;(4)壓縮機(jī)約束和(5)廠際連接約束。

    (1)氫源和氫阱流股約束 氫源流量平衡可表示為

    其中,F(xiàn)p,i、Fp,i,k和Fp,i,m分別表示氫源流率、氫源輸送至氫阱的流率和氫源輸送至提純器的流率。對(duì)于過程氫源(i∈HS)而言,F(xiàn)p,i,fuel是輸送到燃?xì)庀到y(tǒng)的流量;而對(duì)于公用工程氫源(i∈HU)而言,F(xiàn)p,i,fuel是過剩制氫能力。

    氫阱流量物料平衡可表示為

    其中,F(xiàn)p,k和Fp,m,k分別表示氫阱流率和提純器輸送到氫阱的流率。

    氫阱的最小流量約束為

    式中,是氫阱要求的最小流率。

    氫阱的氫氣純度約束為

    (4)壓力約束 如果一個(gè)氫阱或提純器的壓力等級(jí)高于供給氫源,需要通過壓縮機(jī)來提高氫源壓力。本文設(shè)置升壓參數(shù)來判斷某個(gè)匹配連接是否需要壓縮機(jī)來提升壓力的取值由式(26)確定。

    (5)廠際連接約束 氫氣公用工程氫氣禁止輸送至提純器,即

    過程氫源不允許輸送至其他廠的提純器和氫阱,即滿足

    Model T 考慮了廠際集成對(duì)每個(gè)單廠氫氣系統(tǒng)的影響,為一個(gè)混合整數(shù)線性規(guī)劃(MILP)問題。按周期數(shù)求解P 次Model T 即可獲得各子周期總費(fèi)用最低時(shí)的廠際氫氣公用工程的傳輸量。

    2.2 Step 2:廠內(nèi)網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化

    基于第一步所獲取的廠際氫氣公用工程傳輸量,第二步將以單廠氫氣網(wǎng)絡(luò)為對(duì)象,考慮變化操作參數(shù)對(duì)單廠氫氣網(wǎng)絡(luò)的影響,構(gòu)建多周期單廠氫氣網(wǎng)絡(luò)的數(shù)學(xué)優(yōu)化模型,以確定各個(gè)廠的氫氣網(wǎng)絡(luò)結(jié)構(gòu)。在本文中,多周期單廠氫氣網(wǎng)絡(luò)數(shù)學(xué)優(yōu)化模型表示為Model P。與第一步獲取廠際氫氣公用工程傳輸量的Model T 不同,在第二步中,廠內(nèi)網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化的Model P有以下變化:

    (1)Model T的研究對(duì)象是工業(yè)園區(qū)整個(gè)廠際氫氣網(wǎng)絡(luò)系統(tǒng),而Model P 的研究對(duì)象則是工業(yè)園區(qū)內(nèi)單廠的氫氣網(wǎng)絡(luò),即Model P 的對(duì)象是Model T 的一個(gè)子系統(tǒng);

    (2)將工業(yè)園區(qū)多廠認(rèn)為是一個(gè)氫氣網(wǎng)絡(luò)系統(tǒng)時(shí),Model P與Model T的差別在于周期數(shù)的不同,即Model T中子周期之間是相互獨(dú)立的,子周期之間無相互關(guān)系;而在Model P 中,所有子周期共享一套氫氣網(wǎng)絡(luò)結(jié)構(gòu),是相互影響的。

    基于以上差別,本文在Model T 的基礎(chǔ)上,構(gòu)建廠內(nèi)網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化模型Model P。

    2.2.1 目標(biāo)函數(shù) 廠內(nèi)網(wǎng)絡(luò)結(jié)構(gòu)優(yōu)化的目標(biāo)為單廠氫氣網(wǎng)絡(luò)的年度化費(fèi)用最低,可表示為

    操作費(fèi)用的各項(xiàng)差別主要體現(xiàn)在周期數(shù)的不同,具體如下。

    氫氣網(wǎng)絡(luò)的投資費(fèi)用由設(shè)備決定,所有的周期共享一套設(shè)備。因此,單廠氫氣網(wǎng)絡(luò)的各項(xiàng)投資費(fèi)用有以下變化。

    氫氣管道投資成本可表示為

    管道容量需要滿足各子周期的要求,即

    每個(gè)壓縮機(jī)的投資費(fèi)用滿足

    壓縮機(jī)的最大功率不低于各子周期的功率

    提純器的成本為

    提純器容量需要滿足各子周期的要求,即

    其中,zm為判斷第m個(gè)提純器是否存在的二元變量。

    2.2.2 約束條件 過程氫源與燃?xì)庀到y(tǒng)的連接約束為

    提純器的最大流量約束為

    此外,還應(yīng)滿足第一步獲取的廠際氫氣公用工程量約束。

    對(duì)于任意廠,接收的氫氣公用工程量應(yīng)滿足

    式中,HKn表示廠n中的氫阱的集合;代表第n個(gè)廠外部的公用工程的集合。

    對(duì)于任意廠,為了滿足氫氣公用工程向外廠輸出氫氣的要求,其供應(yīng)給廠內(nèi)氫阱的消耗量不得高于第一步確定的結(jié)果,即

    該模型的(1)氫源和氫阱流股約束,(2)提純器約束,(3)匹配連接約束,(4)壓力約束和(5)廠際連接約束等相關(guān)約束與上一步的Model T類似。

    Model P 考慮了變化的操作參數(shù)對(duì)每個(gè)單廠氫氣網(wǎng)絡(luò)的影響,為一個(gè)混合整數(shù)線性規(guī)劃(MILP)問題。依次求解N次Model P 可獲得各廠在年度化費(fèi)用最低時(shí)的氫氣網(wǎng)絡(luò)結(jié)構(gòu)。綜合每個(gè)廠的網(wǎng)絡(luò)結(jié)構(gòu),可以得到該方案下整個(gè)網(wǎng)絡(luò)設(shè)計(jì)方案和費(fèi)用。

    2.3 Step 3:廠際網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化

    為了降低廠際連接投資,可以對(duì)廠際結(jié)構(gòu)再進(jìn)行優(yōu)化。由于廠間連接只涉及氫氣公用工程流股和氫阱流股,與過程氫源無關(guān),因此目標(biāo)函數(shù)和約束條件得到簡(jiǎn)化。將過程氫源和提純器輸送到氫阱的流量固定為第二步得到的值,求解該模型以改善氫氣公用工程氫氣在廠間的分配情況,得到最終的氫氣網(wǎng)絡(luò)設(shè)計(jì)方案。本文定義第三步中的多周期廠際氫氣網(wǎng)絡(luò)數(shù)學(xué)優(yōu)化模型為Model Q。

    2.3.1 目標(biāo)函數(shù) 該模型的目標(biāo)函數(shù)為

    壓縮機(jī)的電費(fèi)

    壓縮機(jī)的功率

    氫氣管道投資成本

    壓縮機(jī)投資成本

    管道容量需要滿足各子周期的要求,即

    壓縮機(jī)的最大功率不低于各子周期的功率

    2.3.2 約束條件 公用工程的消耗量不得高于其最大流量。

    氫阱的輸入流率可表示為

    氫阱的氫純度約束為

    式(56)~式(58)中只有Fp,i,k(i∈HU)是變量。

    公用工程與氫阱之間的流股的流量上下限約束表示為

    Model Q 的約束條件還包括了壓力約束,與Model T類似。

    Model Q 考慮了氫氣公用工程選取對(duì)廠際氫氣系統(tǒng)的影響,也是一個(gè)混合整數(shù)線性規(guī)劃(MILP)問題?;诠潭ǖ膹S內(nèi)網(wǎng)絡(luò)結(jié)構(gòu),Model Q 可獲得廠際氫氣系統(tǒng)在年度化費(fèi)用最低時(shí)的氫氣網(wǎng)絡(luò)結(jié)構(gòu)。

    3 案例分析

    3.1 基礎(chǔ)數(shù)據(jù)

    本節(jié)采用一個(gè)工業(yè)園區(qū)包含3個(gè)廠的氫氣網(wǎng)絡(luò)系統(tǒng)為例闡述所提分步方法在解決廠際氫氣網(wǎng)絡(luò)多周期優(yōu)化問題的有效性和優(yōu)勢(shì)。

    該工業(yè)園區(qū)含三個(gè)化工廠,分別為廠A、廠B 和廠C。園區(qū)中共有22 個(gè)氫源流股(其中S1、S8、S9、S10 和S19 為氫氣公用工程)和16 個(gè)氫阱流股。在本文中,假設(shè)氫氣網(wǎng)絡(luò)除檢修和停工時(shí)間后其年運(yùn)行時(shí)間為8000 h,并將運(yùn)行時(shí)間分為7 個(gè)子周期。各子周期下含氫流股的操作參數(shù)如表1 所示,氫氣公用工程的單價(jià)見表2。需要說明的是,在表1中氫氣公用工程的流率表示的是氫氣公用工程最大的生產(chǎn)能力,而非實(shí)際的生產(chǎn)量。各廠之間的空間距離和其他費(fèi)用參數(shù)可見文獻(xiàn)[20,26-27]。假設(shè)每個(gè)廠內(nèi)僅設(shè)置一臺(tái)提純器。

    本文中的全部結(jié)果均是在GAMS 24.1軟件平臺(tái)上,以CPLEX 作為混合整數(shù)線性規(guī)劃問題的求解器[31],在主頻2.93 GHz/ 6 GB RAM 的計(jì)算機(jī)上進(jìn)行計(jì)算。

    3.2 廠際氫氣網(wǎng)絡(luò)多周期優(yōu)化設(shè)計(jì)三步法的計(jì)算結(jié)果

    本節(jié)將以上述基礎(chǔ)數(shù)據(jù)為案例闡述本文所提方法的實(shí)施步驟和有效性。

    (1)Step 1:廠際氫氣傳輸量的優(yōu)化

    采用Step 1 的Model T 按子周期數(shù)進(jìn)行7 次求解,可以得到該工業(yè)園區(qū)廠際氫氣傳輸量,如表3所示。

    表3 給出了7 個(gè)子周期下3 個(gè)廠之間的氫氣公用工程傳輸情況。由表可知,廠際氫氣網(wǎng)絡(luò)經(jīng)Step 1 優(yōu)化后,廠A 的氫氣公用工程流股S1 主要供給本廠,也有少部分供給廠C;而廠B 中的氫氣公用工程流股S8 則是向A、B、C 三廠供給。在不同的子周期下,廠際間的氫氣公用工程傳輸量是不同的,這主要是由于在各子周期內(nèi)氫源氫阱流股基礎(chǔ)數(shù)據(jù)差異導(dǎo)致的。由表1可知,如果按最大公用工程計(jì)算,廠C 的氫源是過剩的,但在氫氣公用工程流股可選擇的前提下,廠C 并不選擇廠內(nèi)的氫氣公用工程流股,而是選擇了廠A 中的S1 和廠B 中的S8。這是權(quán)衡投資費(fèi)用和操作費(fèi)用后的結(jié)果,也說明了采用氫氣公用工程實(shí)現(xiàn)廠際集成的有效性。

    表1 各子周期下廠際氫氣網(wǎng)絡(luò)的流股數(shù)據(jù)Table 1 Flowrate data of inter-plant hydrogen network in each subperiod

    表2 各廠氫氣公用工程的價(jià)格Table 2 Prices of hydrogen utilities in each plant

    表3 各廠在子周期中消耗的公用工程流量Table 3 Hydrogen utility flowrates consumed by each plant in each subperiod

    (2)Step 2:廠內(nèi)氫氣網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化

    在已知廠際氫氣傳輸量的前提下,采用Step 2中的Model P,按廠的個(gè)數(shù),經(jīng)過3 次求解可以得到各廠在年度化費(fèi)用最低時(shí)的氫氣網(wǎng)絡(luò)結(jié)構(gòu)。綜合每個(gè)廠的網(wǎng)絡(luò)結(jié)構(gòu),可以得到該方案下整個(gè)工業(yè)園區(qū)氫氣網(wǎng)絡(luò)的設(shè)計(jì)方案,如圖3 所示。需要說明的是,此處Step 2所得到的結(jié)果是基于Step 1的優(yōu)化結(jié)果,表示為Step 1-2,即由Step 1和Step 2兩步組成。

    圖3 廠際氫氣網(wǎng)絡(luò)的結(jié)構(gòu)(Step 1-2和Step 1-2-3)Fig.3 The structure of inter-plant hydrogen network(Step 1-2 and Step 1-2-3)

    (3)Step 3:廠際網(wǎng)絡(luò)結(jié)構(gòu)的優(yōu)化

    根據(jù)Step 2 所得各廠氫氣網(wǎng)絡(luò)結(jié)構(gòu),固定各廠廠內(nèi)氫氣網(wǎng)絡(luò)結(jié)構(gòu),通過求解Step 3中的Model Q 可優(yōu)化廠際氫氣公用工程流股和氫阱流股之間匹配關(guān)系和氫氣調(diào)度策略,繼而可獲得廠際氫氣系統(tǒng)在年度化費(fèi)用最低時(shí)的氫氣網(wǎng)絡(luò)結(jié)構(gòu),如圖3 所示。經(jīng)Step 1、Step 2 和Step 3 順次優(yōu)化的方法表示為Step1-2-3。需要說明的是,由Step 1-2和Step 1-2-3 獲得的廠際網(wǎng)絡(luò)結(jié)構(gòu)是相同的,但其廠際的調(diào)度策略和年度化費(fèi)用是有差別的,如表4所示。在表4中Step 2-3 則為由Step 2 和Step 3 組合的優(yōu)化方法,Step 1-2-3 為本文所提出的分步求解廠際氫氣網(wǎng)絡(luò)多周期集成問題的優(yōu)化設(shè)計(jì)方法。

    表4 多廠集成與分廠集成的氫氣系統(tǒng)費(fèi)用對(duì)比Table 4 Comparison of cost for the inter-plant integration and the individual plant integration

    由表4 可知,由前兩步所獲得的廠際氫氣網(wǎng)絡(luò)的費(fèi)用為5.432×108CNY·a-1,而經(jīng)三步優(yōu)化后系統(tǒng)的費(fèi)用為5.422×108CNY·a-1。盡管兩種方法所得到的廠際氫氣網(wǎng)絡(luò)的匹配結(jié)構(gòu)是一樣的,但三步法所獲得的系統(tǒng)費(fèi)用略低,說明經(jīng)Step 3 優(yōu)化后廠際氫氣公用工程流股和氫阱流股匹配關(guān)系和氫氣調(diào)度策略應(yīng)有所調(diào)整。附表A1和附表A2給出了兩種方法所得廠際匹配流股的流率,這也是導(dǎo)致兩種方法所得到的系統(tǒng)總費(fèi)用產(chǎn)生差異的原因。

    同時(shí),本文采用Step 2和Step 3也可以求解該廠際氫氣網(wǎng)絡(luò)的多周期集成問題,與本文所構(gòu)建三步法的總費(fèi)用相比,其優(yōu)化后的系統(tǒng)總費(fèi)用較高,并且也高于由前兩步優(yōu)化后系統(tǒng)的總費(fèi)用。這說明在針對(duì)廠際氫氣網(wǎng)絡(luò)進(jìn)行多周期優(yōu)化設(shè)計(jì)時(shí),優(yōu)化廠際氫氣公用工程傳輸量(Step 1)是非常必要的,因?yàn)樵跉錃饩W(wǎng)絡(luò)優(yōu)化中,操作費(fèi)用占到總費(fèi)用的比例很大。

    圖4 給出經(jīng)Step 2 和Step 3 優(yōu)化后的廠際氫氣網(wǎng)絡(luò)結(jié)構(gòu)。Step 2-3 所得到的氫氣網(wǎng)絡(luò)的匹配數(shù)最少,結(jié)構(gòu)最為簡(jiǎn)單。這也說明了氫氣網(wǎng)絡(luò)的優(yōu)化路徑為通過增加較小的投資費(fèi)用而實(shí)現(xiàn)較大操作費(fèi)用的降低。對(duì)比圖3 和圖4 可見,兩種方法得到氫氣網(wǎng)絡(luò)結(jié)構(gòu)中廠際氫氣傳輸方向和具體的廠際流股匹配均不完全相同。這是由于在Step 2-3 方法中,首先優(yōu)化的是多周期單廠的氫氣網(wǎng)絡(luò),繼而再考慮廠際的集成,縮小了廠際氫氣網(wǎng)絡(luò)集成的尋優(yōu)空間,導(dǎo)致兩種方法所得到的廠際結(jié)構(gòu)不同,進(jìn)而使得Step 2-3 方法的系統(tǒng)費(fèi)用比Step 1-2 和Step 1-2-3 所獲得氫氣網(wǎng)絡(luò)的費(fèi)用高。但以上三種方法所得氫氣網(wǎng)絡(luò)的費(fèi)用是相近的,也說明了本文所提方法的有效性。

    圖4 二步法獲得的氫氣網(wǎng)絡(luò)結(jié)構(gòu)(Step 2-3)Fig.4 The hydrogen network structure obtained by the 2-step method(Step 2-3)

    3.3 關(guān)于不同廠際氫氣網(wǎng)絡(luò)多周期集成方法的分析與討論

    為了進(jìn)一步闡明本文方法用于廠際氫氣網(wǎng)絡(luò)多周期集成的有效性,本節(jié)將從氫氣系統(tǒng)費(fèi)用、網(wǎng)絡(luò)結(jié)構(gòu)和數(shù)學(xué)模型求解時(shí)間等方面對(duì)比分析本文所提分步方法與同步方法以及其他分步方法的廠際氫氣網(wǎng)絡(luò)多周期集成的效果。表5給出不同方法在年度化費(fèi)用最低情況下廠際氫氣網(wǎng)絡(luò)的費(fèi)用和流股匹配數(shù)以及數(shù)學(xué)優(yōu)化模型的求解時(shí)間。

    表5 中,結(jié)構(gòu)合并法(structure-merged method)是針對(duì)各子周期采用單周期廠際氫氣網(wǎng)絡(luò)優(yōu)化模型(即Model T)進(jìn)行優(yōu)化,然后將各子周期的氫氣網(wǎng)絡(luò)結(jié)構(gòu)進(jìn)行合并,進(jìn)而得到整個(gè)廠際氫氣網(wǎng)絡(luò)的設(shè)計(jì)和操作方案;固定結(jié)構(gòu)法(structure-fixed method)是首先對(duì)各子周期采用單周期廠際氫氣網(wǎng)絡(luò)優(yōu)化模型(即Model T)進(jìn)行優(yōu)化,然后以此網(wǎng)絡(luò)結(jié)構(gòu)作為初始結(jié)構(gòu)對(duì)其他子周期進(jìn)行優(yōu)化,進(jìn)而得到廠際優(yōu)化結(jié)果。而同步法(simultaneous method)則是直接對(duì)廠際多周期氫氣網(wǎng)絡(luò)進(jìn)行優(yōu)化設(shè)計(jì)。

    表5 不同廠際氫氣網(wǎng)絡(luò)多周期集成方法對(duì)比Table 5 Comparison of the integration of inter-plant hydrogen network by different methods

    由表5可見,就氫氣網(wǎng)絡(luò)的經(jīng)濟(jì)性而言,由于同步法所獲得的是最優(yōu)解,相比于其他分步法其年度化費(fèi)用最低。結(jié)構(gòu)合并法和兩步法(Step 2-3)所得到的氫氣網(wǎng)絡(luò)年度化費(fèi)用相對(duì)較高。這是因?yàn)檫@兩種方法并沒有對(duì)多周期和多廠這兩個(gè)方面進(jìn)行優(yōu)化,結(jié)構(gòu)合并法重點(diǎn)優(yōu)化的是多廠結(jié)構(gòu),而兩步法(Step 2-3)則考慮了多周期的影響。而本文方法所得到氫氣網(wǎng)絡(luò)的年度化費(fèi)用與同步法優(yōu)化后的費(fèi)用相差不大,這也是因?yàn)楸疚乃岬娜椒ňC合考慮了多周期和多廠兩個(gè)方面的集成潛能,因而本文所提方法所得到的解更接近于同步法。

    由表5可見,就氫氣網(wǎng)絡(luò)的結(jié)構(gòu)復(fù)雜性而言,對(duì)于優(yōu)化后的氫氣網(wǎng)絡(luò)結(jié)構(gòu),分步方法所獲得的結(jié)構(gòu)比同步法所得的結(jié)構(gòu)簡(jiǎn)單。兩步法(Step 2-3)和固定結(jié)構(gòu)法(structure-fixed method)兩種方法所獲得氫氣網(wǎng)絡(luò)的匹配數(shù)最少,但從費(fèi)用上說固定結(jié)構(gòu)法有一定優(yōu)勢(shì)。而對(duì)于本文所提方法,與大多數(shù)分步方法相比,雖然其獲得的氫氣網(wǎng)絡(luò)結(jié)構(gòu)相對(duì)復(fù)雜,但是與同步方法相比仍具有一定優(yōu)勢(shì)。同時(shí),從經(jīng)濟(jì)性和網(wǎng)絡(luò)結(jié)構(gòu)上分析,與同步方法相比,分步方法犧牲了經(jīng)濟(jì)性從而獲得了較為簡(jiǎn)單的廠際氫氣系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)。

    由表5可見,就優(yōu)化設(shè)計(jì)的計(jì)算效率而言,分步方法所需要的計(jì)算時(shí)間遠(yuǎn)低于同步方法的計(jì)算時(shí)間。這是因?yàn)?,分步方法是根?jù)求解問題的特征進(jìn)行模型分解,將大問題的模型拆分為若干個(gè)小問題分別進(jìn)行求解,雖然所需求解的模型數(shù)目增多,但是大大降低了模型的規(guī)模和求解難度,因此減少了整個(gè)問題的求解時(shí)間。在本文案例中,本文所提方法與其他分步法所需要的計(jì)算時(shí)間差別不大。與同步法相比,本文所提方法在計(jì)算效率上具有較大優(yōu)勢(shì),在不改變廠際多周期氫氣網(wǎng)絡(luò)結(jié)構(gòu)復(fù)雜度前提下,本文三步法的計(jì)算效率可提高約30倍。

    與其他分步法相比,本文方法綜合考慮了多周期和多廠對(duì)氫氣網(wǎng)絡(luò)系統(tǒng)集成的影響。在權(quán)衡氫氣網(wǎng)絡(luò)的費(fèi)用和氫氣網(wǎng)絡(luò)匹配數(shù)以及模型的求解時(shí)間方面,本文方法不僅降低了廠際氫氣網(wǎng)絡(luò)多周期優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型規(guī)模,提高了計(jì)算效率,還可獲得接近全局最優(yōu)解的計(jì)算結(jié)果,并降低整個(gè)氫氣系統(tǒng)網(wǎng)絡(luò)系統(tǒng)結(jié)構(gòu)的復(fù)雜度。因此,相對(duì)于同步法和其他分步法,本文所提出的三步法均具有一定優(yōu)勢(shì)。

    4 結(jié)論

    針對(duì)廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)問題,本文采用廠內(nèi)直接回用和提純回用與廠際氫氣公用工程共享相結(jié)合的集成策略,并考慮操作參數(shù)變化對(duì)氫氣系統(tǒng)網(wǎng)絡(luò)集成的影響,提出了一種用于求解廠際氫氣系統(tǒng)多周期優(yōu)化設(shè)計(jì)的三步求解方法。該方法以費(fèi)用最小化為優(yōu)化目標(biāo),首先采用廠際氫氣網(wǎng)絡(luò)的單周期優(yōu)化設(shè)計(jì)模型獲取各子周期下的氫氣公用工程的傳輸量,然后采用單廠氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)模型獲得各廠內(nèi)的氫氣系統(tǒng)結(jié)構(gòu),再采用廠際氫氣系統(tǒng)的多周期優(yōu)化設(shè)計(jì)模型確定廠際氫氣網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu),并確定化工園區(qū)中廠際氫氣系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)和氫氣調(diào)度方案。在分步求解中均為混合整數(shù)線性規(guī)劃問題。

    本文采用一個(gè)三廠7周期的廠際氫氣網(wǎng)絡(luò)優(yōu)化設(shè)計(jì)問題為案例,闡述了本文方法的實(shí)施步驟和有效性,并與固定結(jié)構(gòu)法、結(jié)構(gòu)合并法、同步集成法等方法的計(jì)算結(jié)果進(jìn)行了對(duì)比和分析。研究表明,與其他分步方法相比,本文方法在經(jīng)濟(jì)性上是接近于同步法的最優(yōu)解;與同步法相比,本文在計(jì)算效率上具有顯著優(yōu)勢(shì),可在不改變廠際氫氣網(wǎng)絡(luò)多周期設(shè)計(jì)網(wǎng)絡(luò)結(jié)構(gòu)復(fù)雜度的前提下,將計(jì)算效率提高約30 倍。因此,本文所構(gòu)建的三步求解策略為廠際氫氣網(wǎng)絡(luò)的多周期優(yōu)化設(shè)計(jì)提供了一種有效方法。

    符號(hào)說明

    參數(shù)

    Af——年度化因子

    a——固定投資因子

    b——可變投資因子

    c——比定壓熱容,J·(kg·K)-1

    d——管道長(zhǎng)度,m

    e——能源價(jià)格,CNY·mol-1或CNY·kJ-1

    P——壓力,MPa

    r——提純器回收率

    T——壓縮機(jī)的輸入氣體溫度,K

    t——子周期操作時(shí)長(zhǎng),h

    y——?dú)錃饧兌龋?(mol)

    δ——用于判斷管道連接是否跨廠

    γ——?dú)錃獾谋榷▔簾崛菖c比定容熱容之比

    η——壓縮機(jī)的效率

    連續(xù)變量

    C——費(fèi)用,CNY

    F——?dú)錃饬髀?,mol·s-1

    TAC——年度化總費(fèi)用,CNY

    W——壓縮機(jī)功率,kW

    二元變量

    z——表示是否存在的二元變量

    z′——輔助變量

    集合

    HS——過程氫源

    HU——公用工程

    HK——?dú)溱?/p>

    N——廠

    P——子周期

    PU——提純器

    US——?dú)湓?/p>

    上角標(biāo)

    com——壓縮機(jī)

    ele——電力

    fuel——燃?xì)?/p>

    H2——?dú)錃?/p>

    heat——熱量

    in——廠內(nèi)

    inv——投資

    L——下界

    oper——操作

    out——廠外

    pipe——管道

    prod——提純器產(chǎn)品流股的純度

    pur——提純器

    resd——提純器尾氣

    U——上界

    下角標(biāo)

    fuel——燃?xì)庀到y(tǒng)

    i——?dú)湓?/p>

    k——?dú)溱?/p>

    m——提純器

    n——廠的編號(hào)

    p——子周期

    附錄

    附表A1 廠際氫氣網(wǎng)絡(luò)氫氣公用工程流股與氫阱流股匹配的流量(Step 1-2)Table A1 The flowrate transported from hydrogen utilities to hydrogen sinks in the inter-plant hydrogen network(Step 1-2)

    附表A2 廠際氫氣網(wǎng)絡(luò)氫氣公用工程流股與氫阱流股匹配的流量(Step 1-2-3)Table A2 The flowrate transported from hydrogen utilities to hydrogen sinks in the inter-plant hydrogen network(Step 1-2-3)

    猜你喜歡
    流股公用網(wǎng)絡(luò)結(jié)構(gòu)
    硫化氫化學(xué)反應(yīng)循環(huán)制氫工藝的流程模擬
    低溫甲醇洗吸收塔產(chǎn)出液再生過程模擬研究
    一個(gè)公用品牌的養(yǎng)成——橫山羊肉是咋樣“吃香”的
    公用電梯自動(dòng)取消??垦b置初步設(shè)計(jì)
    電子制作(2017年17期)2017-12-18 06:40:45
    醫(yī)生私車公用撞傷人 醫(yī)院擔(dān)責(zé)
    基于互信息的貝葉斯網(wǎng)絡(luò)結(jié)構(gòu)學(xué)習(xí)
    知識(shí)網(wǎng)絡(luò)結(jié)構(gòu)維對(duì)于創(chuàng)新績(jī)效的作用機(jī)制——遠(yuǎn)程創(chuàng)新搜尋的中介作用
    滬港通下A+ H股票網(wǎng)絡(luò)結(jié)構(gòu)演化的實(shí)證分析
    復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)比對(duì)算法研究進(jìn)展
    基于轉(zhuǎn)運(yùn)模型的功交換網(wǎng)絡(luò)綜合
    亚洲国产高清在线一区二区三| 久久久色成人| 国产 一区 欧美 日韩| 九草在线视频观看| 亚洲熟女精品中文字幕| 搡老乐熟女国产| 国产精品嫩草影院av在线观看| 免费黄网站久久成人精品| 内射极品少妇av片p| 男女啪啪激烈高潮av片| 亚洲最大成人av| 在现免费观看毛片| 一级毛片aaaaaa免费看小| 黄色日韩在线| 99久久中文字幕三级久久日本| 久久99热这里只频精品6学生| 国产精品一区www在线观看| av在线老鸭窝| 国产有黄有色有爽视频| 国产欧美另类精品又又久久亚洲欧美| 春色校园在线视频观看| 观看免费一级毛片| 永久网站在线| 午夜激情久久久久久久| 国产淫片久久久久久久久| 丝袜喷水一区| 国产av码专区亚洲av| 九九久久精品国产亚洲av麻豆| 一区二区三区四区激情视频| 丝瓜视频免费看黄片| 伊人久久精品亚洲午夜| 亚洲成人精品中文字幕电影| 91在线精品国自产拍蜜月| 热re99久久精品国产66热6| 在线播放无遮挡| 久久精品人妻少妇| 欧美日韩视频高清一区二区三区二| 少妇熟女欧美另类| 欧美人与善性xxx| 欧美亚洲 丝袜 人妻 在线| 欧美另类一区| 国产成人aa在线观看| 亚洲丝袜综合中文字幕| 国内少妇人妻偷人精品xxx网站| 又粗又硬又长又爽又黄的视频| 婷婷色综合www| 精品人妻一区二区三区麻豆| 亚洲欧美精品自产自拍| 精品久久国产蜜桃| 乱系列少妇在线播放| 日韩电影二区| 国产精品福利在线免费观看| 91在线精品国自产拍蜜月| 欧美日韩亚洲高清精品| 欧美区成人在线视频| 久久久成人免费电影| 亚洲婷婷狠狠爱综合网| 高清在线视频一区二区三区| 欧美三级亚洲精品| 国国产精品蜜臀av免费| 国产成人aa在线观看| 午夜视频国产福利| 亚洲欧美日韩无卡精品| 精品少妇久久久久久888优播| 美女高潮的动态| 日日摸夜夜添夜夜添av毛片| 精品久久久久久久久亚洲| 91在线精品国自产拍蜜月| 亚洲av不卡在线观看| 岛国毛片在线播放| 精品国产乱码久久久久久小说| 91在线精品国自产拍蜜月| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 丰满乱子伦码专区| 又爽又黄无遮挡网站| 内射极品少妇av片p| 美女主播在线视频| 国产成人免费无遮挡视频| 熟妇人妻不卡中文字幕| 亚洲欧美一区二区三区黑人 | 亚洲最大成人中文| av免费在线看不卡| 99久国产av精品国产电影| 国产成人a区在线观看| 亚洲欧美精品自产自拍| 亚洲国产欧美人成| 欧美xxxx性猛交bbbb| 亚洲人成网站在线观看播放| 日韩中字成人| 色播亚洲综合网| 极品少妇高潮喷水抽搐| 国产av码专区亚洲av| 大码成人一级视频| 99热全是精品| 亚洲在久久综合| 国产一区二区在线观看日韩| 国语对白做爰xxxⅹ性视频网站| 国产色爽女视频免费观看| 国产午夜福利久久久久久| 麻豆国产97在线/欧美| 免费黄色在线免费观看| 国产成人一区二区在线| 成人国产麻豆网| 国语对白做爰xxxⅹ性视频网站| 在线观看一区二区三区激情| 男女啪啪激烈高潮av片| 久久韩国三级中文字幕| 可以在线观看毛片的网站| 日韩av不卡免费在线播放| av在线天堂中文字幕| 欧美成人精品欧美一级黄| 国产精品一区二区三区四区免费观看| av又黄又爽大尺度在线免费看| 成年人午夜在线观看视频| 另类亚洲欧美激情| 精品国产三级普通话版| 亚洲精品国产av成人精品| 中文字幕制服av| 看十八女毛片水多多多| 一级片'在线观看视频| 国产精品99久久99久久久不卡 | 日日摸夜夜添夜夜爱| 免费黄网站久久成人精品| 白带黄色成豆腐渣| 大片电影免费在线观看免费| 女人被狂操c到高潮| 国产毛片a区久久久久| 国产成人91sexporn| 色婷婷久久久亚洲欧美| 人人妻人人爽人人添夜夜欢视频 | h日本视频在线播放| 草草在线视频免费看| 国产免费一区二区三区四区乱码| a级一级毛片免费在线观看| 亚洲欧美精品专区久久| 日韩欧美精品免费久久| 欧美日本视频| 日韩伦理黄色片| 一本久久精品| 亚洲欧美成人精品一区二区| 七月丁香在线播放| 人妻系列 视频| 成人高潮视频无遮挡免费网站| 菩萨蛮人人尽说江南好唐韦庄| 高清欧美精品videossex| 美女被艹到高潮喷水动态| 国产黄频视频在线观看| 男插女下体视频免费在线播放| 亚洲三级黄色毛片| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 精品视频人人做人人爽| 日韩伦理黄色片| 毛片女人毛片| 免费看光身美女| 国产成人免费观看mmmm| 美女高潮的动态| 黄色日韩在线| 亚洲av电影在线观看一区二区三区 | 久久久久久久大尺度免费视频| 亚洲成人av在线免费| 国产国拍精品亚洲av在线观看| 人妻系列 视频| 网址你懂的国产日韩在线| 97超视频在线观看视频| 国产中年淑女户外野战色| 成人高潮视频无遮挡免费网站| 97超视频在线观看视频| 特大巨黑吊av在线直播| 久久久久久九九精品二区国产| 成人漫画全彩无遮挡| 99热这里只有精品一区| 日韩av在线免费看完整版不卡| 国产av国产精品国产| 中文在线观看免费www的网站| 国产精品一区二区性色av| 一级二级三级毛片免费看| 熟妇人妻不卡中文字幕| 日日啪夜夜撸| 一本久久精品| 国产伦理片在线播放av一区| 黑人高潮一二区| 青春草国产在线视频| 人人妻人人澡人人爽人人夜夜| 又粗又硬又长又爽又黄的视频| 女人被狂操c到高潮| 精品国产一区二区三区久久久樱花 | 久久99热这里只有精品18| 国产人妻一区二区三区在| 十八禁网站网址无遮挡 | 3wmmmm亚洲av在线观看| 亚洲欧美清纯卡通| 国产探花极品一区二区| 五月玫瑰六月丁香| 免费看av在线观看网站| 精品熟女少妇av免费看| 少妇猛男粗大的猛烈进出视频 | 欧美日韩视频高清一区二区三区二| 极品少妇高潮喷水抽搐| 亚洲精品日韩在线中文字幕| 最近手机中文字幕大全| 日本熟妇午夜| 国产高清国产精品国产三级 | 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 在线免费十八禁| 赤兔流量卡办理| 亚洲av中文av极速乱| 色哟哟·www| 亚洲av不卡在线观看| 国产女主播在线喷水免费视频网站| 少妇裸体淫交视频免费看高清| 观看免费一级毛片| 视频区图区小说| 亚洲欧美中文字幕日韩二区| 岛国毛片在线播放| 亚洲最大成人av| av在线app专区| 亚洲精品亚洲一区二区| 高清欧美精品videossex| 亚洲,欧美,日韩| xxx大片免费视频| 亚洲成人一二三区av| 欧美日韩一区二区视频在线观看视频在线 | 国产成人aa在线观看| h日本视频在线播放| 男女国产视频网站| 成人综合一区亚洲| 美女国产视频在线观看| 18禁动态无遮挡网站| 午夜福利网站1000一区二区三区| av在线播放精品| 成人一区二区视频在线观看| 色吧在线观看| 亚洲av在线观看美女高潮| 在线观看av片永久免费下载| 日本一二三区视频观看| 国产精品.久久久| 欧美少妇被猛烈插入视频| 伊人久久国产一区二区| 2021天堂中文幕一二区在线观| 大香蕉97超碰在线| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 午夜亚洲福利在线播放| 亚洲欧美日韩无卡精品| 青春草亚洲视频在线观看| 久久久久性生活片| 亚洲av中文字字幕乱码综合| 神马国产精品三级电影在线观看| 狂野欧美激情性xxxx在线观看| 人妻系列 视频| 麻豆成人av视频| 久久热精品热| 国产视频首页在线观看| 成年人午夜在线观看视频| 国产淫语在线视频| 永久网站在线| 91aial.com中文字幕在线观看| 久久久久久久精品精品| 国产精品久久久久久精品古装| 亚洲国产精品999| 汤姆久久久久久久影院中文字幕| 欧美日韩视频高清一区二区三区二| 欧美3d第一页| 白带黄色成豆腐渣| 最新中文字幕久久久久| 一级毛片我不卡| 国产探花极品一区二区| 久久久久久国产a免费观看| 欧美成人午夜免费资源| 精品人妻视频免费看| 人妻少妇偷人精品九色| 国产男女超爽视频在线观看| 国产亚洲午夜精品一区二区久久 | 少妇人妻一区二区三区视频| 国产精品久久久久久精品电影小说 | 高清午夜精品一区二区三区| 国产黄片视频在线免费观看| 欧美人与善性xxx| 国产精品三级大全| 乱码一卡2卡4卡精品| 日本色播在线视频| 久久精品久久久久久久性| 国内少妇人妻偷人精品xxx网站| 亚洲va在线va天堂va国产| 亚洲欧美日韩东京热| 亚洲最大成人av| 国产探花在线观看一区二区| 亚洲av男天堂| 伊人久久国产一区二区| 九九久久精品国产亚洲av麻豆| 男人狂女人下面高潮的视频| 中文欧美无线码| 最近的中文字幕免费完整| 永久网站在线| 热99国产精品久久久久久7| 国产精品.久久久| 在线观看美女被高潮喷水网站| 日韩欧美一区视频在线观看 | 3wmmmm亚洲av在线观看| 极品教师在线视频| 亚洲欧美日韩无卡精品| 亚洲成人一二三区av| 婷婷色麻豆天堂久久| 成年人午夜在线观看视频| 精品国产露脸久久av麻豆| 女人十人毛片免费观看3o分钟| av在线亚洲专区| 少妇的逼好多水| 日本一本二区三区精品| 亚州av有码| 精品亚洲乱码少妇综合久久| 草草在线视频免费看| 如何舔出高潮| 欧美bdsm另类| 熟妇人妻不卡中文字幕| 亚洲天堂av无毛| 国产精品国产三级国产专区5o| 国产永久视频网站| 国产人妻一区二区三区在| 欧美xxⅹ黑人| 亚洲精品日本国产第一区| 亚洲电影在线观看av| 嫩草影院入口| av天堂中文字幕网| 久久午夜福利片| 国产成人精品福利久久| 毛片一级片免费看久久久久| 久久久久久久久久久免费av| 美女内射精品一级片tv| 老司机影院毛片| 久久6这里有精品| 精品少妇黑人巨大在线播放| 校园人妻丝袜中文字幕| 久久精品国产亚洲av涩爱| 青春草亚洲视频在线观看| 国产精品国产av在线观看| 日韩av免费高清视频| 午夜精品一区二区三区免费看| 久久久久久久久大av| 久久精品国产a三级三级三级| 日韩伦理黄色片| 国产精品久久久久久久久免| 五月伊人婷婷丁香| 免费黄网站久久成人精品| 欧美bdsm另类| 久久精品国产a三级三级三级| 久久久久精品久久久久真实原创| 亚洲美女搞黄在线观看| 亚洲国产精品专区欧美| 国产 一区 欧美 日韩| 高清日韩中文字幕在线| 国产色爽女视频免费观看| 我的女老师完整版在线观看| 国产亚洲av片在线观看秒播厂| 国产成人91sexporn| 国产在线一区二区三区精| 少妇人妻 视频| 18禁动态无遮挡网站| av国产免费在线观看| 久久人人爽人人片av| 亚洲精品日韩在线中文字幕| 国产亚洲一区二区精品| av免费观看日本| 国产精品久久久久久精品古装| 国产午夜精品久久久久久一区二区三区| 亚洲精品成人av观看孕妇| 国内少妇人妻偷人精品xxx网站| 简卡轻食公司| 久久久午夜欧美精品| 久久人人爽人人片av| 2022亚洲国产成人精品| 韩国高清视频一区二区三区| 亚洲精品日韩av片在线观看| 欧美极品一区二区三区四区| 国产精品av视频在线免费观看| 欧美高清性xxxxhd video| 日韩伦理黄色片| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 99热这里只有精品一区| 高清日韩中文字幕在线| 男人添女人高潮全过程视频| 亚洲欧洲日产国产| 中国美白少妇内射xxxbb| 欧美一区二区亚洲| 黄色一级大片看看| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区在线观看99| 国产成人a∨麻豆精品| 黄色日韩在线| 亚洲最大成人手机在线| av网站免费在线观看视频| 日韩 亚洲 欧美在线| 国产黄色免费在线视频| 又爽又黄无遮挡网站| 精华霜和精华液先用哪个| 中文字幕免费在线视频6| 一区二区av电影网| 国产一区有黄有色的免费视频| 精品午夜福利在线看| 成年女人在线观看亚洲视频 | 女人久久www免费人成看片| 国产成人a∨麻豆精品| 大陆偷拍与自拍| 免费观看的影片在线观看| 永久网站在线| 少妇人妻精品综合一区二区| 日本与韩国留学比较| 啦啦啦啦在线视频资源| 亚洲欧美成人综合另类久久久| 高清在线视频一区二区三区| 又粗又硬又长又爽又黄的视频| 久久这里有精品视频免费| 狂野欧美激情性xxxx在线观看| 色5月婷婷丁香| 少妇裸体淫交视频免费看高清| 看免费成人av毛片| 能在线免费看毛片的网站| 亚洲欧美精品自产自拍| 国产中年淑女户外野战色| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 久久久久久伊人网av| 99热这里只有是精品50| 丰满乱子伦码专区| 日本熟妇午夜| 一区二区三区乱码不卡18| 日本免费在线观看一区| 亚洲人与动物交配视频| 97人妻精品一区二区三区麻豆| 丝瓜视频免费看黄片| 亚洲精品乱久久久久久| 看黄色毛片网站| 观看美女的网站| 白带黄色成豆腐渣| 热99国产精品久久久久久7| 天天躁夜夜躁狠狠久久av| 色视频www国产| 五月伊人婷婷丁香| 最近最新中文字幕大全电影3| 一边亲一边摸免费视频| 中文资源天堂在线| av卡一久久| 亚洲电影在线观看av| 熟女人妻精品中文字幕| 国产午夜精品久久久久久一区二区三区| 国产欧美日韩一区二区三区在线 | 丰满人妻一区二区三区视频av| 国产淫语在线视频| 在线观看美女被高潮喷水网站| 白带黄色成豆腐渣| 久久久色成人| 国产精品偷伦视频观看了| 国产爽快片一区二区三区| 七月丁香在线播放| 亚洲人成网站在线观看播放| 亚洲欧美日韩卡通动漫| 制服丝袜香蕉在线| 久久99热6这里只有精品| 成人国产麻豆网| 亚洲av成人精品一区久久| 人人妻人人澡人人爽人人夜夜| 亚洲高清免费不卡视频| 亚洲欧美一区二区三区黑人 | 春色校园在线视频观看| 亚洲精品视频女| 超碰97精品在线观看| 亚洲国产av新网站| av专区在线播放| 99九九线精品视频在线观看视频| 波野结衣二区三区在线| 国产精品一区www在线观看| 欧美日韩在线观看h| 日本三级黄在线观看| 久久人人爽人人爽人人片va| 在线精品无人区一区二区三 | 国产日韩欧美亚洲二区| 97在线视频观看| 日韩国内少妇激情av| 午夜老司机福利剧场| 天天躁夜夜躁狠狠久久av| 蜜桃亚洲精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 欧美变态另类bdsm刘玥| 久久97久久精品| 精品熟女少妇av免费看| 韩国高清视频一区二区三区| 国产精品久久久久久精品电影| 秋霞伦理黄片| 嫩草影院新地址| 夜夜看夜夜爽夜夜摸| 久久人人爽人人片av| 视频区图区小说| 亚洲三级黄色毛片| 校园人妻丝袜中文字幕| 永久网站在线| 国产精品久久久久久精品古装| 国产精品一及| 日韩在线高清观看一区二区三区| 精品久久久久久久末码| 男人舔奶头视频| 99热这里只有精品一区| 精品久久国产蜜桃| 伦理电影大哥的女人| 亚洲国产色片| 哪个播放器可以免费观看大片| 老师上课跳d突然被开到最大视频| 综合色av麻豆| 久久女婷五月综合色啪小说 | 欧美三级亚洲精品| 日本黄色片子视频| 国产精品国产三级国产av玫瑰| 欧美老熟妇乱子伦牲交| 一区二区三区免费毛片| 亚洲精品亚洲一区二区| 亚洲精品一区蜜桃| 欧美另类一区| 国产精品一区二区在线观看99| 国产成人aa在线观看| 男女下面进入的视频免费午夜| 中文字幕久久专区| 亚洲欧美日韩卡通动漫| 成人毛片a级毛片在线播放| 18禁动态无遮挡网站| 午夜激情久久久久久久| 成人高潮视频无遮挡免费网站| 久久久色成人| eeuss影院久久| 久久精品熟女亚洲av麻豆精品| 亚洲欧美日韩卡通动漫| 国内少妇人妻偷人精品xxx网站| 亚洲精品国产av蜜桃| 免费电影在线观看免费观看| 欧美日本视频| 美女cb高潮喷水在线观看| 国产午夜精品久久久久久一区二区三区| 国产综合懂色| 超碰av人人做人人爽久久| 一级毛片黄色毛片免费观看视频| 美女高潮的动态| 国产欧美另类精品又又久久亚洲欧美| 亚洲,欧美,日韩| 2021少妇久久久久久久久久久| 欧美日韩视频精品一区| 亚洲三级黄色毛片| 卡戴珊不雅视频在线播放| 久久久亚洲精品成人影院| 蜜桃久久精品国产亚洲av| 一级毛片电影观看| 99精国产麻豆久久婷婷| 欧美丝袜亚洲另类| 99热网站在线观看| 麻豆成人av视频| 日本欧美国产在线视频| 亚洲精品日韩av片在线观看| 国产淫语在线视频| 嫩草影院入口| av在线观看视频网站免费| 又爽又黄a免费视频| 久久精品综合一区二区三区| 亚洲色图av天堂| 精品一区二区三区视频在线| 成人黄色视频免费在线看| videos熟女内射| 狂野欧美激情性xxxx在线观看| 网址你懂的国产日韩在线| 丝袜脚勾引网站| 自拍偷自拍亚洲精品老妇| 亚洲av不卡在线观看| 日韩av免费高清视频| 日本熟妇午夜| 一本一本综合久久| 日韩欧美精品v在线| 国产亚洲av片在线观看秒播厂| 夫妻午夜视频| 少妇高潮的动态图| 在现免费观看毛片| 日本与韩国留学比较| 中文字幕人妻熟人妻熟丝袜美| 男女边吃奶边做爰视频| 免费大片黄手机在线观看| 一本一本综合久久| 日本色播在线视频| 最近最新中文字幕大全电影3| 中国国产av一级| 欧美区成人在线视频| 日韩一区二区三区影片| 国产成人91sexporn| 免费看不卡的av| 欧美成人一区二区免费高清观看| 久久热精品热| 亚洲av不卡在线观看| 一边亲一边摸免费视频| 欧美人与善性xxx| 成人无遮挡网站| 久久久午夜欧美精品| 日日啪夜夜爽| 噜噜噜噜噜久久久久久91| 黄色日韩在线| 国产 精品1| 欧美日韩精品成人综合77777| 男人添女人高潮全过程视频| 亚洲精品久久午夜乱码| 制服丝袜香蕉在线| 蜜桃亚洲精品一区二区三区| 国产成人精品婷婷| www.av在线官网国产| 51国产日韩欧美| 国内揄拍国产精品人妻在线| 青春草视频在线免费观看| 日韩成人av中文字幕在线观看| 亚洲电影在线观看av|