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

    集成輕烴回收單元代理模型的氫氣網(wǎng)絡(luò)多目標(biāo)優(yōu)化

    2022-04-26 09:48:48張淑君王詩慧張欣吉旭戴一陽黨亞固周利
    化工學(xué)報(bào) 2022年4期
    關(guān)鍵詞:流股輕烴氫氣

    張淑君,王詩慧,張欣,吉旭,戴一陽,黨亞固,周利

    (四川大學(xué)化學(xué)工程學(xué)院,四川成都 610065)

    引 言

    在過去的十年間,全球變暖趨勢(shì)逐漸加劇,環(huán)境問題不斷突出。如何改善環(huán)境、確保可持續(xù)發(fā)展已成為人類當(dāng)前面臨的嚴(yán)峻問題。2019 年聯(lián)合國(guó)專門召開了氣候行動(dòng)峰會(huì),敦促各國(guó)制定碳減排計(jì)劃[1]。其后,各國(guó)紛紛宣布達(dá)到碳中和目標(biāo)[2]。煉油企業(yè)屬于能源密集型流程工業(yè),是溫室氣體排放大戶[3],雖然國(guó)家和社會(huì)正在努力控制過程工業(yè)的溫室氣體排放量,但當(dāng)前煉油企業(yè)面臨著一系列可能進(jìn)一步導(dǎo)致其溫室氣體排放量增加的因素。受制于全球原油重質(zhì)化加劇及成品油質(zhì)量提升的雙重壓力,煉油企業(yè)不得不加大對(duì)原油的加工比例和處理深度,大量增設(shè)加氫反應(yīng)裝置,氫氣需求量快速增長(zhǎng),氫氣成本已成為煉油企業(yè)僅次于原油成本的第二大成本[4],同時(shí)煉油企業(yè)在生產(chǎn)和使用氫氣的過程中,將導(dǎo)致大量的溫室氣體排放[5]。如何優(yōu)化煉油企業(yè)氫氣網(wǎng)絡(luò),使其達(dá)到氫氣網(wǎng)絡(luò)成本和溫室氣體排放同時(shí)降低,實(shí)現(xiàn)經(jīng)濟(jì)效益和環(huán)境保護(hù)雙贏,這對(duì)大多數(shù)煉油企業(yè)來說是個(gè)巨大的挑戰(zhàn)。

    煉油企業(yè)氫氣網(wǎng)絡(luò)優(yōu)化的主要方法可分為夾點(diǎn)分析法和基于超結(jié)構(gòu)的數(shù)學(xué)規(guī)劃法。Alves等[6]首次提出利用剩余氫量圖來識(shí)別夾點(diǎn)并確定最小氫氣消耗量,之后的研究學(xué)者不斷對(duì)其進(jìn)行改進(jìn),提出氣體級(jí)聯(lián)圖法[7]、源阱負(fù)荷曲線[8]等。夾點(diǎn)法雖可通過圖形形式對(duì)氫氣網(wǎng)絡(luò)進(jìn)行形象的表述,但其在處理實(shí)際生產(chǎn)過程中的約束及大規(guī)模求解等問題上具有一定的限制。Hallale 等[9]提出基于超結(jié)構(gòu)的數(shù)學(xué)規(guī)劃法,通過對(duì)每個(gè)單元及物流節(jié)點(diǎn)建立質(zhì)量和能量衡算方程以及相應(yīng)的操作可行性約束,求解得出滿足約束條件的目標(biāo)函數(shù)最優(yōu)方案。目前數(shù)學(xué)規(guī)劃方法廣泛應(yīng)用于氫氣網(wǎng)絡(luò)優(yōu)化問題中,例如集成提純單元[10?12]以及雜質(zhì)脫除單元[13?15]、不確定性[16?18]、多周期調(diào)度[19]等。雖然已經(jīng)涌現(xiàn)很多的科研成果和成功的應(yīng)用案例,但大多數(shù)研究主要聚焦于氫氣網(wǎng)絡(luò)的耗氫量或費(fèi)用最小化,氫氣網(wǎng)絡(luò)的環(huán)境影響并未受到足夠的關(guān)注。

    在氫氣網(wǎng)絡(luò)中,低分氣中含有豐富的氫氣和輕烴(C1~C5),對(duì)低分氣中輕烴組分進(jìn)行回收,有利于后續(xù)氫氣提純裝置的正常運(yùn)行[20],提高氫氣的利用率;同時(shí),回收的C3+輕烴可作為煉廠原料或產(chǎn)品,一方面可以提高煉油廠整體經(jīng)濟(jì)效益,另一方面燃料系統(tǒng)補(bǔ)充碳排放少的天然氣來代替C3+輕烴[21],有助于減少氫氣網(wǎng)絡(luò)溫室氣體的排放。在氫氣網(wǎng)絡(luò)中回收輕烴組分,已有人做了一定的研究,Deng等[22]提出了將輕烴回收單元與煉廠氫氣網(wǎng)絡(luò)結(jié)合的系統(tǒng)改造方案。使用夾點(diǎn)法與Aspen 過程建模仿真相結(jié)合,以最小氫氣用量為優(yōu)化目標(biāo)。Yang 等[23]將改進(jìn)的夾點(diǎn)法與嚴(yán)格的工藝模擬相結(jié)合,在Matlab和Aspen平臺(tái)實(shí)現(xiàn)了混合方法。

    上述集成輕烴回收單元的氫氣網(wǎng)絡(luò)優(yōu)化均采用的夾點(diǎn)優(yōu)化法,其對(duì)氫氣網(wǎng)絡(luò)中實(shí)際約束考慮可能不足,且研究中對(duì)輕烴回收單元均建立嚴(yán)格的機(jī)理模型,高度復(fù)雜的機(jī)理模型雖可以準(zhǔn)確地捕捉到單元中的物理化學(xué)過程,但其會(huì)使氫氣網(wǎng)絡(luò)優(yōu)化求解計(jì)算成本較大,尤其是將機(jī)理模型嵌入到數(shù)學(xué)規(guī)劃模型中,可能會(huì)導(dǎo)致后續(xù)求解困難。代理模型是一種基于數(shù)據(jù)驅(qū)動(dòng)建立的模型,其具有計(jì)算復(fù)雜度低,且計(jì)算結(jié)果接近嚴(yán)格機(jī)理模型結(jié)果等特點(diǎn),目前在系統(tǒng)優(yōu)化設(shè)計(jì)等研究領(lǐng)域中被廣泛地應(yīng)用[24]。Wang 等[25]和Xia 等[14]分別在氫氣網(wǎng)絡(luò)數(shù)學(xué)規(guī)劃模型中集成了閃蒸單元和脫硫單元的代理模型,以實(shí)際案例證明代理模型可以在氫氣網(wǎng)絡(luò)優(yōu)化模型中發(fā)揮非常好的性能,使得優(yōu)化的計(jì)算成本大幅降低,同時(shí)還能保證優(yōu)化模型的高保真性。

    因此,本研究提出在氫氣網(wǎng)絡(luò)數(shù)學(xué)規(guī)劃模型中集成輕烴回收單元代理模型,在保證模型精度情況下,減少計(jì)算成本,更好地探究氫氣網(wǎng)絡(luò)與輕烴回收單元之間的質(zhì)量集成機(jī)理。同時(shí),將氫氣網(wǎng)絡(luò)的環(huán)境影響納入優(yōu)化目標(biāo),建立一種高效計(jì)算環(huán)境友好型氫氣網(wǎng)絡(luò)多目標(biāo)數(shù)學(xué)規(guī)劃模型,以年度總費(fèi)用為標(biāo)準(zhǔn)評(píng)價(jià)氫氣網(wǎng)絡(luò)的經(jīng)濟(jì)性能,以年度總CO2排放量評(píng)估氫氣網(wǎng)絡(luò)的環(huán)境影響。

    1 基于狀態(tài)空間超結(jié)構(gòu)的煉廠氫氣網(wǎng)絡(luò)

    狀態(tài)空間超結(jié)構(gòu)最早在換熱網(wǎng)絡(luò)設(shè)計(jì)優(yōu)化中提出,其通過對(duì)研究問題中所有設(shè)備單元做出可能的連接,形成網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu),進(jìn)而捕獲豐富的網(wǎng)絡(luò)特征,有助于對(duì)資源網(wǎng)絡(luò)的建模和優(yōu)化[26]。

    本研究采用該模型來表征煉廠氫氣網(wǎng)絡(luò)的結(jié)構(gòu)特征。如圖1 所示,它由6 個(gè)相互連接的部分構(gòu)成:石油精煉單元、氫氣分配單元、脫硫單元、輕烴回收單元、提純單元和壓縮機(jī)單元。氫氣分配單元對(duì)氫源和氫阱之間的氫氣進(jìn)行分配。提純單元主要采用變壓吸附或膜分離的方式,從低于直接回用要求的氫氣流股中回收氫氣資源。脫硫單元將氫氣流股中的有害雜質(zhì)H2S 進(jìn)行脫除,來減少其在整個(gè)系統(tǒng)中的積累。輕烴回收單元從富烴氣體中回收比甲烷或乙烷更重的組分,同時(shí)達(dá)到氫氣富集的作用。壓縮機(jī)單元用來提升流股的壓力水平,使其滿足生產(chǎn)工藝要求。從圖中看出,這五個(gè)單元是一種緊密耦合的關(guān)系。

    圖1 氫氣網(wǎng)絡(luò)的狀態(tài)空間超結(jié)構(gòu)Fig.1 State space superstructure of hydrogen network

    建立氫氣網(wǎng)絡(luò)超結(jié)構(gòu)后,需要將狀態(tài)空間超結(jié)構(gòu)中所有可能的連接情況用數(shù)學(xué)方法表達(dá)出來,進(jìn)而形成該網(wǎng)絡(luò)的數(shù)學(xué)規(guī)劃模型。

    2 輕烴回收單元代理模型的開發(fā)

    2.1 構(gòu)建輕烴回收單元嚴(yán)格機(jī)理模型

    目前回收輕烴的方法主要有吸附分離法、深冷分離法、油吸收法等[27]。考慮到經(jīng)濟(jì)性及可行性,本研究采用油吸收法對(duì)煉廠低分氣中輕烴組分進(jìn)行回收,選用煉廠自產(chǎn)的粗石腦油作為吸收劑。利用Aspen plus 對(duì)輕烴回收單元進(jìn)行建模和模擬,圖2為輕烴回收單元的流程示意圖。

    圖2 輕烴回收單元流程圖Fig.2 Flow chart of light hydrocarbon recovery unit

    如圖所示,該輕烴回收工藝中包括三部分,分別是吸收塔、脫乙烷塔和脫丁烷塔。經(jīng)脫硫處理后的低分氣加壓后進(jìn)入吸收塔。在吸收塔內(nèi),粗石腦油與加壓低分氣逆流接觸,將氣體中大部分C3+輕烴吸收,經(jīng)過吸收后的氣體流股從塔頂排出,此時(shí)流股主要組成為氫氣和甲烷,這股氣體送往變壓吸附進(jìn)行氫氣提純回用。吸收輕烴組分的粗石腦油則從吸收塔塔底被送入脫乙烷塔回收乙烷,之后,來自脫乙烷塔底的粗石腦油送往脫丁烷塔,分離出C3+輕烴組分,這部分輕烴可以被送至乙烯裂解裝置作為原料,填補(bǔ)煉廠輕烴裂解的需求缺口,也可以經(jīng)處理后以LPG形式出售。

    輕烴回收單元是一個(gè)典型的質(zhì)量交換過程,其投資成本和運(yùn)行費(fèi)用隨單元進(jìn)口氫氣流股的流量和輕烴回收程度而變化。回收的輕烴量越多,氫氣和輕烴的再利用率越高,這也就意味著較高的回收運(yùn)行成本,一般要通過更多的吸收劑和更大的分離能耗實(shí)現(xiàn)。

    2.2 輕烴回收單元代理模型的開發(fā)

    考慮到嚴(yán)格的輕烴回收單元機(jī)理模型包含各設(shè)備模型方程、物性計(jì)算方程、流程聯(lián)接方程等,具有變量多、方程維數(shù)大、非線性強(qiáng)等特點(diǎn),直接將其嵌入到氫氣網(wǎng)絡(luò)優(yōu)化模型中進(jìn)行迭代尋優(yōu)將十分耗時(shí),為提高計(jì)算效率,建立代理模型來逼近嚴(yán)格的輕烴回收過程。圖3為構(gòu)建輕烴回收單元代理模型的步驟。

    圖3 輕烴回收單元代理模型的構(gòu)建步驟Fig.3 Steps for constructing the surrogate model of light hydrocarbon recovery unit

    第一步:根據(jù)研究問題及建模需要確定輸入及輸出變量,本研究選取輸入輸出變量具體見圖4。

    圖4 代理模型的輸入與輸出變量Fig.4 Input and output variables of the surrogate model

    第二步:選擇合適的代理模型。在系統(tǒng)優(yōu)化設(shè)計(jì)中,響應(yīng)面模型、克里金和人工神經(jīng)網(wǎng)絡(luò)等代理模型應(yīng)用較為廣泛。針對(duì)輕烴回收過程這種復(fù)雜工業(yè)的高維數(shù)問題,選用響應(yīng)面模型作為輕烴回收單元模型的擬合形式,不僅在優(yōu)化設(shè)計(jì)過程中有助于噪聲的快速收斂,并且得到的模型可以直接用數(shù)學(xué)公式進(jìn)行表達(dá)[28],如式(1)所示,方便將代理模型嵌入氫氣網(wǎng)絡(luò)數(shù)學(xué)規(guī)劃模型中。

    其中,xn表示輸入變量;y表示輸出變量;fnn′(xn,xn′)和f12…N(x1,x2,…,xN)表示二階項(xiàng)和N階項(xiàng)。

    第三步:在模型擬合范圍內(nèi)對(duì)輸入數(shù)據(jù)點(diǎn)進(jìn)行采樣。Garud等[29]研究表明,遵守空間填充標(biāo)準(zhǔn)的采樣技術(shù)對(duì)代理模型的精確度有很大的影響,即代理模型構(gòu)建時(shí)選擇的采樣技術(shù)應(yīng)該具有快速、均勻填充整個(gè)采樣空間的能力。本研究選擇Sobol 隨機(jī)序列采樣法,其在高低維度上均表現(xiàn)出采樣的高效性及空間分布均勻性等優(yōu)點(diǎn)[29]。采樣過程首先在[0,1]N空間內(nèi)進(jìn)行采樣,之后通過反歸一化公式變換為模型輸入范圍的對(duì)應(yīng)的數(shù)值。

    第四步:通過Aspen 對(duì)樣本點(diǎn)進(jìn)行仿真模擬。利用前一節(jié)建立的Aspen 嚴(yán)格輕烴回收單元模擬模型,得到嚴(yán)格機(jī)理模型下輸入值對(duì)應(yīng)的輸出值,并剔除Aspen不收斂的數(shù)據(jù)點(diǎn)。

    第五步:利用輸入值及輸出值來構(gòu)建代理模型,并對(duì)模型進(jìn)行評(píng)價(jià)。對(duì)代理模型的準(zhǔn)確性評(píng)價(jià)指標(biāo)有很多種,本研究選擇判定系數(shù)(R2)、均方根誤差(RMSE)以及殘差圖,來評(píng)估模型的準(zhǔn)確性和可靠性,計(jì)算式見式(2)~式(4)。若模型準(zhǔn)確性滿足要求,則輸出響應(yīng)面模型參數(shù);若準(zhǔn)確性不滿足要求,一方面可以增加新的采樣點(diǎn),另一方面可以改變響應(yīng)面模型的階數(shù),進(jìn)而得到滿足要求的響應(yīng)面模型。

    其中,m是驗(yàn)證模型所用的數(shù)據(jù)點(diǎn)的個(gè)數(shù);yi是第i個(gè)采樣點(diǎn)的輸出值;y?是代理模型在第i個(gè)點(diǎn)處的預(yù)測(cè)值;yˉ是所有采樣點(diǎn)輸出值的平均值;εi代表在第i個(gè)點(diǎn)處的殘差。

    3 基于數(shù)學(xué)規(guī)劃法的氫氣網(wǎng)絡(luò)優(yōu)化模型

    本研究的氫氣網(wǎng)絡(luò)數(shù)學(xué)規(guī)劃模型是依據(jù)第2節(jié)狀態(tài)空間超結(jié)構(gòu)建立的,模型可分為氫氣分配單元、提純單元、脫硫單元、輕烴回收單元、壓縮機(jī)單元和目標(biāo)函數(shù)。

    3.1 氫氣分配單元

    來自氫源的新氫,一般可以直接輸送到氫阱或者經(jīng)壓縮機(jī)壓縮后輸送到氫阱。如式(5)所示,氫源輸送的流股總流量不得超過其產(chǎn)能上限。每個(gè)氫阱入口的質(zhì)量衡算式可以由式(6)和式(7)表示。式(8)和式(9)表示了氫阱出口高分氣和低分氣的衡算方程。式(10)、式(11)分別為氫阱入口的氫氣濃度約束和H2S 濃度約束。通過式(12)確定氫阱入口流股的壓力,當(dāng)多股氫源流股混合時(shí)其壓力取混合流股中的最小壓力。式(13)表示了氫阱入口流股的壓力約束。

    3.2 提純單元

    對(duì)于提純單元而言,其入口的質(zhì)量衡算可由式(14)和式(15)表示。待提純的氫氣流股在提純單元被分為兩部分,如式(16)所示,一部分是含有高濃度氫氣的產(chǎn)品流股,另一部分是殘余氣流股。式(17)為提純過程中的各組分質(zhì)量衡算。高濃度的氫氣產(chǎn)品流股可以直接被送至氫阱利用,也可以送至壓縮機(jī)加壓,如式(18)所示,而殘余氣流股因?yàn)闅錃夂枯^低,一般送至燃料系統(tǒng)燃燒,如式(19)所示。式(20)為提純單元入口的H2S 濃度約束。式(21)用來確定提純單元入口流股的壓力,式(22)給出了提純單元入口流股的壓力約束。式(23)可以根據(jù)提純單元的氫氣回收率計(jì)算出產(chǎn)品流股中的氫氣濃度。產(chǎn)品流股中的其他各組分的濃度變化可以由式(24)和式(25)計(jì)算得出。

    3.3 脫硫單元

    脫硫單元入口的高分氣質(zhì)量衡算由式(26)和式(27)中給出。式(28)和式(29)為高壓脫硫單元進(jìn)口處的壓力約束。脫硫單元低分氣的入口質(zhì)量衡算和壓力約束都與高分氣的類似,在此不列出具體的計(jì)算方程。式(30)為脫硫過程的質(zhì)量衡算式。式(31)和式(32)表示高壓和低壓脫硫單元出口處的質(zhì)量衡算,高壓脫硫單元出口的流股可以直接氫阱回用或送到壓縮機(jī)單元加壓處理;低壓脫硫單元出口的流股含有一定量的氫氣和輕烴組分,具有較高的回收價(jià)值,可送去輕烴回收單元或送至瓦斯燃料系統(tǒng)。

    3.4 輕烴回收單元

    來自低壓脫硫單元的氫氣流股中進(jìn)入輕烴回收單元,回收其中的C3+輕烴組分,同時(shí)完成氫氣的富集。由于輕烴回收單元有入口的壓力限制,所以流股需進(jìn)行加壓處理,如式(33)所示,流股送往壓縮機(jī)單元。輕烴回收單元的入口質(zhì)量衡算由式(34)和式(35)給出。式(36)定義了輕烴回收單元進(jìn)口處壓力的相關(guān)約束。式(37)為輕烴回收單元出口處的質(zhì)量衡算,經(jīng)輕烴回收后氫氣流股送往提純單元,回收的C3+組分送往煉廠其他裝置進(jìn)行處理。需要注意的是,輕烴回收過程中的相關(guān)方程式都是由訓(xùn)練好的代理模型來表示的,因此在這里并沒有列出。

    3.5 壓縮機(jī)單元

    式(38)表示壓縮機(jī)單元入口流股的質(zhì)量衡算,來自氫源、氫阱、提純單元、脫硫單元、輕烴回收單元的流股均可以送到壓縮機(jī)單元進(jìn)行加壓處理。式(39)表示壓縮機(jī)單元出口流股的質(zhì)量衡算,被加壓后的流股可以送至氫阱、提純單元、脫硫單元、輕烴回收單元進(jìn)行進(jìn)一步的利用。如果壓縮機(jī)是現(xiàn)有的,其進(jìn)口處的流量需要進(jìn)行約束,如式(41)所示。壓縮機(jī)入口流股的壓力由式(42)確定。如果現(xiàn)有的壓縮機(jī)無法滿足當(dāng)前的加壓需求,那么需要新增設(shè)壓縮機(jī),新增的壓縮機(jī)的額定功率由進(jìn)出口流股的壓力和流量共同決定[30],如式(43)所示。

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

    3.6.1 經(jīng)濟(jì)性目標(biāo)函數(shù) 選用總年度費(fèi)用(TAC)作為考察系統(tǒng)經(jīng)濟(jì)性能的目標(biāo)函數(shù),如式(44)所示。CH2、Cmdea分別表示氫氣的消耗成本及脫硫單元脫硫劑的消耗成本,計(jì)算式為式(45)和式(46)。Celec為電力成本,主要包括壓縮機(jī)以及輕烴回收單元輸送泵的電耗,如式(47)所示,式(48)為輸送泵的運(yùn)行功率。Cheat為熱交換費(fèi)用,如式(49)所示,由冷/熱公用工程單位費(fèi)用與代理模型計(jì)算出來的換熱量相乘得到。Efuel是燃料系統(tǒng)燃燒廢氣產(chǎn)生熱量的年度效益,這部分可由式(50)計(jì)算得出。ELHR是輕烴回收單元回收C3+輕烴產(chǎn)生的收益,如式(51)所示,注意因?yàn)榛厥樟薈3+組分,燃料系統(tǒng)需額外補(bǔ)充天然氣,故此處C3+輕烴的單價(jià)為市場(chǎng)價(jià)減去相同熱值天然氣的價(jià)格。式(52)和式(53)分別給出了改造管道以及建設(shè)輕烴回收裝置的成本計(jì)算公式[18]。

    3.6.2 環(huán)境影響目標(biāo)函數(shù) 選用總年度CO2排放量(TCE)作為評(píng)估系統(tǒng)環(huán)境影響的目標(biāo)函數(shù),如式(54)所示,煉油企業(yè)氫氣網(wǎng)絡(luò)CO2排放源主要為燃燒排放、工藝排放、間接排放三類[31]。燃燒排放源為燃料系統(tǒng)中廢氫流股和外購(gòu)燃料的燃燒,由于輕烴回收單元回收了流股中C3+輕烴組分,使得送往燃料系統(tǒng)的流股熱值減少,燃料系統(tǒng)需要補(bǔ)充同等熱值的天然氣組分,因此燃料系統(tǒng)的溫室氣體排放量CEfuel等于廢氫流股燃燒產(chǎn)生的CO2排放量加上補(bǔ)充的天然氣燃燒產(chǎn)生的CO2排放量,如式(55)所示。工藝排放源主要是制氫單元,制氫單元CO2排放量與氫氣產(chǎn)量有關(guān),如式(56)所示,制氫單元CO2排放量CEh為制氫單元溫室氣體排放因子CEFh,i與氫氣供應(yīng)量的乘積。間接排放源主要是指煉廠外購(gòu)的電力等能源在生產(chǎn)環(huán)節(jié)中產(chǎn)生的CO2排放,如式(57)所示,來自壓縮機(jī)及泵的運(yùn)行電耗引起的CO2排放CEelec等于電廠的溫室氣體排放因子CEFPower乘耗電量。

    3.6.3 多目標(biāo)求解 對(duì)多目標(biāo)優(yōu)化問題,要使多個(gè)目標(biāo)在給定區(qū)域盡可能最佳,其解通常是一組均衡解,即一組由眾多Pareto 最優(yōu)解組成的最優(yōu)解集合。目前常用的多目標(biāo)優(yōu)化算法包括加權(quán)法[32]、??約束法[33]、遺傳進(jìn)化算法[34]等。本文采用Kim 等[35]提出的自適應(yīng)加權(quán)求和法進(jìn)行目標(biāo)優(yōu)化問題的求解,自適應(yīng)加權(quán)求和法在傳統(tǒng)加權(quán)求和法基礎(chǔ)上進(jìn)行改進(jìn),其能夠產(chǎn)生分布均勻的Pareto 最優(yōu)解并且可以在Pareto 前沿的非凸區(qū)域找到Pareto 最優(yōu)解。自適應(yīng)加權(quán)求和法求解該多目標(biāo)問題步驟如下。

    第一步,將本研究中的多目標(biāo)優(yōu)化問題轉(zhuǎn)化為單目標(biāo)優(yōu)化問題,如式(58)所示,ω1、ω2分別是經(jīng)濟(jì)性能和環(huán)境影響的權(quán)重系數(shù);f1(TAC)、f2(TCE)分別為離差標(biāo)準(zhǔn)化處理后的經(jīng)濟(jì)性能目標(biāo)和環(huán)境影響目標(biāo),式(59)和式(60)所示為離差標(biāo)準(zhǔn)化處理過程。

    第二步,確定初始優(yōu)化點(diǎn)數(shù)N,計(jì)算權(quán)值步長(zhǎng)Δw,如式(61)所示,并由此得出初始的一組權(quán)值,計(jì)算該組權(quán)值下的Pareto 最優(yōu)解。將得到的Pareto 最優(yōu)解繪制于坐標(biāo)系中,計(jì)算各點(diǎn)間距離li,以及各點(diǎn)之間對(duì)應(yīng)的子區(qū)間需要修飾的次數(shù)ni,如式(62)所示,式中l(wèi)avg為所有子區(qū)間的平均距離,C為常數(shù),round為取整函數(shù)。

    第三步,確定子區(qū)間兩端點(diǎn)間的偏差距離。以某一子區(qū)間為例,連接該子區(qū)間兩端點(diǎn)Q1、Q2,構(gòu)建分段線性割線。定義沿割線的偏差距離σJ,并計(jì)算夾角θ,進(jìn)而計(jì)算得到該子區(qū)間各優(yōu)化子目標(biāo)的偏差值σ1和σ2,計(jì)算式如(63)~式(65)所示。其中,和分別是 區(qū)間端 點(diǎn)Q1、Q2處對(duì)應(yīng) 的TAC和TCE值。

    第四步,得到該子區(qū)間各優(yōu)化子目標(biāo)偏差值后,在原優(yōu)化問題中加入新的不等式約束,如式(66)和式(67)所示。子區(qū)間內(nèi)權(quán)重步長(zhǎng)由修飾次數(shù)ni決定,計(jì)算在新約束下該子區(qū)間的不同權(quán)重下的Pareto 最優(yōu)解。對(duì)其他子區(qū)間均重復(fù)該步驟,即可得到所有子區(qū)間內(nèi)Pareto最優(yōu)解。

    4 案例應(yīng)用與結(jié)果分析

    4.1 案例基本情況

    將本文所提出的方法應(yīng)用于中國(guó)西部某煉油企業(yè)的改造案例中。煉油企業(yè)現(xiàn)有的氫氣網(wǎng)絡(luò)由兩個(gè)氫源、五個(gè)氫阱構(gòu)成,氫源為一個(gè)制氫裝置(H2plant)和一個(gè)連續(xù)重整裝置(CCR)。氫阱為兩個(gè)柴油加氫單元(DHT?1 和DHT?2)、一個(gè)汽油加氫單元(GHT)以及兩個(gè)航煤油加氫單元(KHT?1和KHT?2)。工廠現(xiàn)有一套由于經(jīng)常積液而閑置的變壓吸附裝置,裝置運(yùn)行回收率為0.88,產(chǎn)品流股氫氣的純度為99.00%。圖5 為當(dāng)前氫氣網(wǎng)絡(luò)的結(jié)構(gòu)示意圖。表1 列出了當(dāng)前氫氣網(wǎng)絡(luò)中各流股的詳細(xì)信息。表2 為氫氣網(wǎng)絡(luò)中各個(gè)單元之間的管道距離。表3 列出了各氫阱入口的H2濃度下限和H2S 濃度上限。

    表3 氫阱入口流股的濃度約束Table 3 Concentration constraint of the inlet stream of the hydrogen sink

    圖5 當(dāng)前氫氣網(wǎng)絡(luò)的結(jié)構(gòu)示意圖Fig.5 Schematic diagram of the current hydrogen network structure

    表1 氫氣網(wǎng)絡(luò)中的相關(guān)流股的詳細(xì)信息Table 1 Detailed information of related streams in the hydrogen network

    表2 案例中各單元之間的管道距離Table 2 The pipe distance between the units in the case

    4.2 輕烴回收單元代理模型開發(fā)

    在本案例中使用第3節(jié)所描述的代理模型開發(fā)方法,對(duì)輕烴回收單元中的吸收塔、脫乙烷塔及脫丁烷塔進(jìn)行建模。圖4展示了三個(gè)塔代理模型構(gòu)建過程所選的輸入和輸出變量。如圖所示,對(duì)吸收塔選取的輸入變量為入口流股的組分流量(H2、H2S 和C1~C5)以及吸收劑(粗石腦油)的流量。輸出變量為吸收后塔頂氫氣流股中的組分流量(H2、H2S 和C1~C5)。脫乙烷塔輸入變量與吸收塔一致,輸出變量除了脫乙烷后塔底流股中的組分流量外,還包括脫乙烷塔塔頂塔底的換熱負(fù)荷。對(duì)于脫丁烷塔,由于在脫乙烷塔內(nèi)H2和C1 組分幾乎完全被分離,所以脫丁烷塔的輸入變量為入口流股的組分(C2~C5和H2S)流量以及吸收劑的流量。輸出變量為塔頂流股的組分流量以及脫丁烷塔塔頂塔底的換熱負(fù)荷。

    表4列出了吸收塔代理模型的輸入變量取值范圍。這些取值范圍都是根據(jù)實(shí)際工程的正常運(yùn)行條件設(shè)置的。脫乙烷塔代理模型的輸入變量取值范圍則是根據(jù)吸收塔運(yùn)行條件設(shè)置,同樣脫丁烷塔的代理模型的輸入變量取值范圍根據(jù)脫乙烷塔運(yùn)行條件設(shè)置。在確定好代理模型的擬合空間后,采用Sobol 采樣法生成800個(gè)輸入變量數(shù)據(jù)樣本,并通過Aspen 仿真模擬計(jì)算輸出變量的數(shù)值。

    表4 輕烴回收單元吸收塔輸入變量范圍Table 4 Input variable range of absorption tower of light hydrocarbon recovery unit

    將輕烴回收單元的輸入和輸出數(shù)據(jù)進(jìn)行擬合,通過對(duì)模型精度和復(fù)雜性方面綜合考慮進(jìn)行比較,選擇二階響應(yīng)面模型作為嚴(yán)格輕烴回收過程模型的近似。表5 列出了輕烴回收單元模型的驗(yàn)證結(jié)果,對(duì)于得到的二階響應(yīng)面模型,它的R2值達(dá)到0.99,各個(gè)組分輸出變量的RMSE 值均低于0.21。同時(shí),為驗(yàn)證了模型的可靠性,對(duì)輕烴回收單元代理模型做殘差圖,從圖6 可以看出輕烴回收單元代理模型的殘差均在“0”附近隨機(jī)分布,說明殘差值與回歸預(yù)測(cè)值之間無內(nèi)在聯(lián)系。通過對(duì)輕烴回收單元中三個(gè)塔的代理模型驗(yàn)證結(jié)果可以充分證明所選擇的代理模型可以很好地逼近輕烴回收單元吸收過程,并提供可靠的預(yù)測(cè)結(jié)果,同時(shí)實(shí)現(xiàn)了以簡(jiǎn)單的形式捕獲了嚴(yán)格輕烴回收單元的豐富特征。

    表5 輕烴回收單元模型的驗(yàn)證結(jié)果Table 5 Validation results of light hydrocarbon recovery unit model

    圖6 輕烴回收單元代理模型的殘差圖Fig.6 Residual plot of the surrogate model of the light hydrocarbon recovery unit

    4.3 嵌入輕烴回收單元的氫氣分配網(wǎng)絡(luò)優(yōu)化

    將得到的輕烴回收單元代理模型集成到該案例研究中的數(shù)學(xué)規(guī)劃模型中,利用GAMS 軟件求解。最終的數(shù)學(xué)規(guī)劃模型是一個(gè)非線性規(guī)劃模型,為了保證求得的解是全局最優(yōu)解,選擇BARON[36]作為該案例的求解器。計(jì)算機(jī)硬件信息為AMD Ryzen 2.10 GHz/16GB RAM,GAMS 軟件版本為32,求解時(shí)間均小于0.98 CPUs。計(jì)算得到該多目標(biāo)優(yōu)化模型的Pareto 最優(yōu)解集后,繪制Pareto 曲線,每個(gè)點(diǎn)代表了不同權(quán)重系數(shù)下對(duì)應(yīng)的氫氣網(wǎng)絡(luò)結(jié)構(gòu),如圖7 所示,其較好展示出氫氣網(wǎng)絡(luò)總年度費(fèi)用與總年度CO2排放量之間的關(guān)系。曲線最左端表示的是經(jīng)濟(jì)性能權(quán)重為1,即以總年度費(fèi)用最小為優(yōu)化目標(biāo)求解得到的氫氣網(wǎng)絡(luò)結(jié)構(gòu),此時(shí)氫氣網(wǎng)絡(luò)的CO2排放量最大;曲線最右端表示的是環(huán)境影響權(quán)重為1,即以系統(tǒng)年度CO2排放量最小為優(yōu)化目標(biāo)求解得到的氫氣網(wǎng)絡(luò)結(jié)構(gòu),此時(shí)氫氣網(wǎng)絡(luò)的總年度費(fèi)用最大。接下來分別分析這兩點(diǎn)對(duì)應(yīng)的氫氣網(wǎng)絡(luò)優(yōu)化結(jié)構(gòu),來探究集成輕烴回收單元的氫氣網(wǎng)絡(luò)中經(jīng)濟(jì)效益與CO2排放之間的內(nèi)在關(guān)系。

    圖7 多目標(biāo)優(yōu)化后得到的Pareto曲線Fig.7 Pareto curve obtained by multi?objective optimization

    4.3.1 最小年度費(fèi)用為優(yōu)化目標(biāo) 以總年度費(fèi)用最小為優(yōu)化目標(biāo)對(duì)模型進(jìn)行求解,根據(jù)求解結(jié)果可以繪制出圖8 所示的氫氣網(wǎng)絡(luò)優(yōu)化分配方案,該模型的計(jì)算結(jié)果建議安裝輕烴回收單元來回收低分氣中的C3+組分。如圖8 所示,低分氣先送至低壓脫硫單元進(jìn)行脫硫,與原氫氣網(wǎng)絡(luò)不同的是,脫硫后的低分氣送至輕烴回收單元,回收其中的C3+組分,之后送往變壓吸附單元進(jìn)行氫氣提純,提純后氫氣流股中氫氣含量達(dá)99.9%,可直接氫阱回用,氫源新氫供應(yīng)量減少了48.84 mol/s?;厥蛰p烴流股53.60 mol/s,其中C3、C4、C5 組分的摩爾分?jǐn)?shù)分別為43.14%、51.30%、3.09%,煉廠可以對(duì)這部分輕烴進(jìn)行處理,以LPG 形式出售,或送至煉廠的乙烯裂解裝置作為原料使用。對(duì)于系統(tǒng)年度CO2排放,由于新建了輕烴回收單元,使氫氣被更好地回收利用,氫氣網(wǎng)絡(luò)CO2工藝排放較原網(wǎng)絡(luò)減少了17%,同時(shí)輕烴回收單元回收了C3+組分,因此燃料系統(tǒng)補(bǔ)充了更清潔的天然氣燃料,燃料系統(tǒng)CO2排放減少了27%,但輕烴回收單元的壓縮機(jī)及泵的運(yùn)行新增了電耗,導(dǎo)致氫氣網(wǎng)絡(luò)間接溫室氣體排放增加28%。綜合計(jì)算,氫氣網(wǎng)絡(luò)總年度CO2排放較原網(wǎng)絡(luò)系統(tǒng)減少20%,總年度費(fèi)用較原氫氣網(wǎng)絡(luò)系統(tǒng)減少17%,證明了集成輕烴回收單元能有效降低氫氣網(wǎng)絡(luò)的年度費(fèi)用及溫室氣體排放。

    圖8 最小總年度費(fèi)用的氫氣網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.8 Hydrogen network structure diagram with minimum total annual cost

    4.3.2 最小年度CO2排放為優(yōu)化目標(biāo) 以年度CO2排放量最小為優(yōu)化目標(biāo)求解模型,根據(jù)求解結(jié)果繪制圖9 所示的氫氣網(wǎng)絡(luò)優(yōu)化分配方案,它與以年度費(fèi)用最小為優(yōu)化目標(biāo)的氫氣網(wǎng)絡(luò)主要區(qū)別在于:脫硫后的低分氣送至輕烴回收單元,吸收塔內(nèi)吸收劑用量增大,流股中C3+組分幾乎全部被回收,此時(shí)送往變壓吸附單元的氫氣流股濃度為70.3%,回收的輕烴流股增加為64.82 mol/s,C3、C4、C5 組分的摩爾分?jǐn)?shù)分別為51.64%、42.42%、2.55%。因?yàn)榛厥盏腃3+輕烴增多,燃料系統(tǒng)補(bǔ)充了更多清潔的天然氣燃料,該優(yōu)化方案年度CO2排放較原網(wǎng)絡(luò)系統(tǒng)減少21%,但由于吸收劑的用量增加,導(dǎo)致輕烴回收單元的分離能耗增加,使得氫氣網(wǎng)絡(luò)系統(tǒng)年度費(fèi)用增加,總年度費(fèi)用較原氫氣網(wǎng)絡(luò)系統(tǒng)減少16%,如圖10所示。

    圖9 最小總年度CO2排放的氫氣網(wǎng)絡(luò)結(jié)構(gòu)圖Fig.9 Hydrogen network structure diagram with minimum total annual CO2 emissions

    圖10 原氫氣網(wǎng)絡(luò)與兩種優(yōu)化后氫氣網(wǎng)絡(luò)的TAC和TCE對(duì)比Fig.10 Comparison of TAC and TCE between the original hydrogen network and the two optimized hydrogen networks

    4.3.3 年度費(fèi)用與CO2排放之間的權(quán)衡關(guān)系 圖11展示了不同Pareto 最優(yōu)解下氫氣網(wǎng)絡(luò)的CO2排放組成,隨著經(jīng)濟(jì)性能的權(quán)重系數(shù)減小,氫氣網(wǎng)絡(luò)環(huán)境影響越發(fā)重要,通過代理模型計(jì)算出的輕烴回收單元的最佳吸收劑用量在逐漸增加。結(jié)合上述對(duì)年度費(fèi)用最小和年度CO2排放量最小的氫氣網(wǎng)絡(luò)結(jié)構(gòu)的分析,可以得到集成輕烴回收單元的氫氣網(wǎng)絡(luò)年度費(fèi)用與CO2排放之間的權(quán)衡關(guān)系:即隨著輕烴回收單元中吸收劑的增加,吸收的C3+組分隨之增加,直接排放到燃料系統(tǒng)的C3+組分減少,如圖11 所示,燃料系統(tǒng)排放CO2隨之減少;而由于吸收劑用量增加,會(huì)導(dǎo)致輕烴回收單元的分離能耗增加,使得系統(tǒng)間接排放CO2增加,氫氣網(wǎng)絡(luò)的費(fèi)用也同時(shí)增加。但總體而言,隨著輕烴回收單元吸收劑增加,氫氣網(wǎng)絡(luò)總年度CO2排放呈遞減趨勢(shì),而總年度費(fèi)用則呈增加趨勢(shì)。在實(shí)際應(yīng)用中,煉廠可以根據(jù)LPG 市場(chǎng)價(jià)格波動(dòng)、當(dāng)前的碳排放指標(biāo)以及裂解原料供應(yīng)等因素,來調(diào)整氫氣網(wǎng)絡(luò)結(jié)構(gòu),使其達(dá)到期望的經(jīng)濟(jì)性能與環(huán)境影響。

    圖11 不同Pareto最優(yōu)解的CO2排放組成Fig.11 CO2 emission composition of different Pareto optimal solutions

    5 結(jié) 論

    全球變暖是人類當(dāng)前面臨的一個(gè)重要且緊迫的問題,本文針對(duì)目前煉油企業(yè)氫氣需求量持續(xù)增長(zhǎng)導(dǎo)致的煉廠成本及溫室氣體排放增加這一現(xiàn)狀,提出了在基于數(shù)學(xué)規(guī)劃法的氫氣網(wǎng)絡(luò)優(yōu)化模型中集成輕烴回收單元??紤]到直接嵌入輕烴回收單元嚴(yán)格機(jī)理模型可能導(dǎo)致的優(yōu)化模型求解困難的問題,選用代理模型方法對(duì)輕烴回收單元建模,建立的二階響應(yīng)面代理模型的R2達(dá)0.99,殘差均在“0”附近隨機(jī)分布,集成輕烴回收單元代理模型的氫氣網(wǎng)絡(luò)數(shù)學(xué)規(guī)劃模型求解時(shí)間均小于0.98 s,實(shí)現(xiàn)了以較低計(jì)算成本捕獲嚴(yán)格輕烴回收單元的豐富特征。

    本文同時(shí)將氫氣網(wǎng)絡(luò)的環(huán)境影響納入優(yōu)化目標(biāo),分別以年度費(fèi)用和年度CO2排放量評(píng)估氫氣網(wǎng)絡(luò)的經(jīng)濟(jì)性能和環(huán)境影響,建立一種耦合輕烴回收單元代理模型的煉廠氫氣網(wǎng)絡(luò)多目標(biāo)數(shù)學(xué)規(guī)劃模型。將該方法應(yīng)用在某煉廠實(shí)際案例中,計(jì)算表明所提出的方法能夠有效降低氫氣網(wǎng)絡(luò)的年度費(fèi)用及溫室氣體排放,并且揭示集成輕烴回收單元的氫氣網(wǎng)絡(luò)系統(tǒng)經(jīng)濟(jì)性能與環(huán)境影響之間的內(nèi)在關(guān)系。在實(shí)際應(yīng)用中,煉廠可以根據(jù)當(dāng)前的市場(chǎng)波動(dòng)以及碳排放配額指標(biāo)等因素調(diào)整氫氣網(wǎng)絡(luò)結(jié)構(gòu),使其達(dá)到經(jīng)濟(jì)效益與環(huán)境保護(hù)雙贏,來更好地應(yīng)對(duì)碳中和挑戰(zhàn),實(shí)現(xiàn)煉油企業(yè)的可持續(xù)發(fā)展。

    符 號(hào) 說 明

    集合

    CP——壓縮機(jī)集合

    DS——脫硫單元集合(DS=1為高壓脫硫單元,DS=2為低壓脫硫單元)

    I——?dú)湓醇?/p>

    J,K——?dú)溱寮?/p>

    N——流股組分集合

    PF——提純裝置集合

    系數(shù)

    a,b——投資成本系數(shù)

    CV——燃料熱值,kJ/kmol

    ER——人民幣對(duì)美元的匯率

    FR0,CR0——輕烴回收裝置的成本相關(guān)系數(shù)

    R——回收率/脫除率,%

    t——操作時(shí)間,h

    ZJ——年度因子

    η——設(shè)備機(jī)械效率,%

    變量

    C——各項(xiàng)操作費(fèi)用及設(shè)備投資費(fèi)用,CNY

    CE——CO2排放量,t

    CEF——CO2排放因子

    CN——碳原子數(shù)

    D——管道直徑,m

    E——收益,CNY

    F——流股流量,kmol/h

    g——重力加速度,m/s2

    H——輸送泵揚(yáng)程,m

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

    P——壓力,bar(1 bar=105Pa)

    POW——功率,kW·h

    Q——換熱量,k·W

    qv——輸送泵流量,m3/s

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

    TCE——年度總CO2排放,t

    UP——單價(jià),CNY/(kW·a)

    Y——流股濃度,%

    ρ——輸送流體密度,kg/m3

    上角標(biāo)

    d——脫丁烷塔

    HP——高分氣

    in——裝置入口流股

    LP——低分氣

    max——最大值

    min——最小值

    n——生產(chǎn)能力指數(shù)

    out——裝置出口流股

    P——產(chǎn)品流股

    R——?dú)堄鄽饬鞴?/p>

    y——脫乙烷塔

    下角標(biāo)

    C3+——C3及C3以上輕烴組分

    cp——壓縮機(jī)

    cool——冷換熱器

    d——脫丁烷塔

    ds——脫硫裝置

    elec——電

    fuel——燃料

    h——新鮮氫源

    heat——熱能

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

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

    LHR——輕烴回收單元

    mdea——脫硫劑

    n——?dú)淞鞴山M分

    pf——提純裝置

    pipe——管道

    pm——輸送泵

    psa——變壓吸附裝置

    red——熱換熱器

    y——脫乙烷塔

    猜你喜歡
    流股輕烴氫氣
    硫化氫化學(xué)反應(yīng)循環(huán)制氫工藝的流程模擬
    輕烴輔助的生物質(zhì)能源供應(yīng)可行性分析
    云南化工(2020年11期)2021-01-14 00:51:10
    穩(wěn)定輕烴脫硫處理工藝進(jìn)展
    甘肅科技(2020年20期)2020-04-13 00:30:04
    低溫甲醇洗吸收塔產(chǎn)出液再生過程模擬研究
    提高輕烴產(chǎn)量的技術(shù)研究與應(yīng)用
    一種用氫氣還原制備石墨烯的方法
    基于轉(zhuǎn)運(yùn)模型的功交換網(wǎng)絡(luò)綜合
    氫氣對(duì)缺血再灌注損傷保護(hù)的可能機(jī)制
    輕烴回收直接換熱工藝優(yōu)化研究
    小天使·六年級(jí)語數(shù)英綜合(2015年2期)2015-01-22 04:50:59
    国产成年人精品一区二区 | 亚洲精品国产区一区二| 老司机深夜福利视频在线观看| 成人18禁在线播放| 女人被狂操c到高潮| 最新在线观看一区二区三区| 亚洲色图 男人天堂 中文字幕| 久久人妻熟女aⅴ| 国产av精品麻豆| 成在线人永久免费视频| 午夜福利免费观看在线| 琪琪午夜伦伦电影理论片6080| 久久久水蜜桃国产精品网| 久久精品影院6| 琪琪午夜伦伦电影理论片6080| 亚洲狠狠婷婷综合久久图片| 免费少妇av软件| 免费看a级黄色片| 亚洲精品一区av在线观看| 男女午夜视频在线观看| 午夜影院日韩av| 久久国产乱子伦精品免费另类| 久久人人97超碰香蕉20202| 丝袜在线中文字幕| 亚洲全国av大片| 999久久久国产精品视频| 亚洲成人久久性| 亚洲少妇的诱惑av| 国产三级在线视频| 在线观看一区二区三区激情| 欧美av亚洲av综合av国产av| 搡老岳熟女国产| 国产不卡一卡二| 少妇被粗大的猛进出69影院| 久久中文看片网| 国产三级黄色录像| 在线天堂中文资源库| 国产一区二区激情短视频| 国产亚洲欧美在线一区二区| 欧美激情极品国产一区二区三区| 91麻豆精品激情在线观看国产 | 久久久国产成人免费| 免费观看人在逋| 天天添夜夜摸| 精品久久久久久久毛片微露脸| 亚洲人成77777在线视频| 99国产综合亚洲精品| 国产伦人伦偷精品视频| 嫩草影院精品99| 久久久久久久精品吃奶| 亚洲五月婷婷丁香| 后天国语完整版免费观看| 在线观看免费午夜福利视频| 免费不卡黄色视频| 欧美丝袜亚洲另类 | 国内久久婷婷六月综合欲色啪| 国产精品 欧美亚洲| 级片在线观看| 日日摸夜夜添夜夜添小说| 亚洲av成人一区二区三| 午夜91福利影院| 久久久久国产精品人妻aⅴ院| 麻豆成人av在线观看| 不卡一级毛片| 亚洲伊人色综图| avwww免费| av视频免费观看在线观看| 亚洲成人国产一区在线观看| 欧美黑人欧美精品刺激| av有码第一页| 日本免费a在线| 国产成人精品无人区| 国产精品成人在线| 欧美日韩福利视频一区二区| 久久久精品国产亚洲av高清涩受| 精品国产一区二区久久| 少妇的丰满在线观看| 欧美大码av| 美女福利国产在线| 久久久久国产精品人妻aⅴ院| 国产精品免费视频内射| 亚洲精品中文字幕在线视频| 亚洲精品一区av在线观看| 午夜a级毛片| 亚洲熟妇中文字幕五十中出 | 亚洲色图综合在线观看| 看黄色毛片网站| 日本vs欧美在线观看视频| 99精国产麻豆久久婷婷| 国产不卡一卡二| 亚洲国产精品合色在线| 亚洲国产欧美一区二区综合| 国产一区二区激情短视频| av在线天堂中文字幕 | 国产区一区二久久| 欧美另类亚洲清纯唯美| 人人澡人人妻人| 搡老熟女国产l中国老女人| 777久久人妻少妇嫩草av网站| 国产精品综合久久久久久久免费 | 成人永久免费在线观看视频| 又黄又粗又硬又大视频| 水蜜桃什么品种好| 天天躁夜夜躁狠狠躁躁| 国产精品久久久av美女十八| 午夜福利一区二区在线看| 亚洲精品粉嫩美女一区| 一个人观看的视频www高清免费观看 | 色综合站精品国产| 狠狠狠狠99中文字幕| 久久精品亚洲av国产电影网| 日本撒尿小便嘘嘘汇集6| svipshipincom国产片| av国产精品久久久久影院| 香蕉久久夜色| 午夜两性在线视频| 亚洲成人精品中文字幕电影 | www.www免费av| 免费不卡黄色视频| 国产高清国产精品国产三级| 免费一级毛片在线播放高清视频 | 在线观看66精品国产| 神马国产精品三级电影在线观看 | 在线天堂中文资源库| 亚洲av熟女| bbb黄色大片| 成人影院久久| 欧美日韩国产mv在线观看视频| 国产精品久久久久成人av| 黄网站色视频无遮挡免费观看| 91字幕亚洲| 国产又色又爽无遮挡免费看| 国产一区二区三区综合在线观看| 欧美成人性av电影在线观看| 91国产中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 精品久久久久久久久久免费视频 | www.自偷自拍.com| 男女下面进入的视频免费午夜 | 久久久国产成人精品二区 | 亚洲久久久国产精品| 欧美成人性av电影在线观看| 日韩有码中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 在线观看舔阴道视频| 欧美成人午夜精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品久久成人aⅴ小说| 99精国产麻豆久久婷婷| 后天国语完整版免费观看| 成熟少妇高潮喷水视频| 我的亚洲天堂| 手机成人av网站| 男人舔女人的私密视频| 国产精品久久久人人做人人爽| 日韩视频一区二区在线观看| 又大又爽又粗| 久久草成人影院| 黄色毛片三级朝国网站| 亚洲专区国产一区二区| 村上凉子中文字幕在线| 88av欧美| 精品福利永久在线观看| 久久这里只有精品19| 中文字幕最新亚洲高清| www日本在线高清视频| 麻豆久久精品国产亚洲av | 黄色片一级片一级黄色片| 在线播放国产精品三级| 女人爽到高潮嗷嗷叫在线视频| 亚洲熟女毛片儿| 精品一区二区三卡| 女生性感内裤真人,穿戴方法视频| av网站免费在线观看视频| 国产97色在线日韩免费| 亚洲 欧美一区二区三区| 麻豆成人av在线观看| 国产真人三级小视频在线观看| 亚洲熟妇熟女久久| 男人的好看免费观看在线视频 | 精品国产一区二区久久| 国产精品一区二区在线不卡| 超碰成人久久| 自线自在国产av| 美女高潮喷水抽搐中文字幕| 91av网站免费观看| 久久性视频一级片| 日本撒尿小便嘘嘘汇集6| 脱女人内裤的视频| 亚洲五月色婷婷综合| 最好的美女福利视频网| 天堂动漫精品| 欧美在线一区亚洲| 少妇的丰满在线观看| 国产成人精品久久二区二区免费| 欧美亚洲日本最大视频资源| 午夜a级毛片| 俄罗斯特黄特色一大片| tocl精华| 级片在线观看| 老司机亚洲免费影院| www.自偷自拍.com| 男女床上黄色一级片免费看| 精品福利永久在线观看| 欧美丝袜亚洲另类 | 午夜两性在线视频| 在线观看www视频免费| 精品国产一区二区三区四区第35| 国产亚洲精品久久久久5区| 一级作爱视频免费观看| 美女扒开内裤让男人捅视频| 一本综合久久免费| 精品人妻1区二区| av有码第一页| 欧美人与性动交α欧美精品济南到| 多毛熟女@视频| 女人爽到高潮嗷嗷叫在线视频| 丰满饥渴人妻一区二区三| 丁香欧美五月| 久久久久久久午夜电影 | 国产精品1区2区在线观看.| 国产亚洲精品第一综合不卡| 久久香蕉国产精品| 午夜精品在线福利| 正在播放国产对白刺激| 国产精品国产高清国产av| 日本五十路高清| 亚洲精品国产精品久久久不卡| 欧美人与性动交α欧美精品济南到| 成人18禁高潮啪啪吃奶动态图| 又黄又爽又免费观看的视频| 亚洲欧美日韩无卡精品| 淫秽高清视频在线观看| 波多野结衣一区麻豆| 搡老乐熟女国产| av免费在线观看网站| 搡老岳熟女国产| 叶爱在线成人免费视频播放| 无人区码免费观看不卡| 免费观看精品视频网站| 国产一区二区三区视频了| 18禁国产床啪视频网站| 欧美成人性av电影在线观看| 99久久国产精品久久久| 久久青草综合色| 一区二区三区激情视频| 亚洲成人国产一区在线观看| 欧美在线一区亚洲| 精品久久蜜臀av无| 黄色毛片三级朝国网站| 纯流量卡能插随身wifi吗| 日本精品一区二区三区蜜桃| 亚洲五月天丁香| 校园春色视频在线观看| 欧美 亚洲 国产 日韩一| 国产欧美日韩一区二区三区在线| 中文字幕人妻丝袜一区二区| 久久久久九九精品影院| 国产精品永久免费网站| 777久久人妻少妇嫩草av网站| 久久国产精品男人的天堂亚洲| 久久久精品欧美日韩精品| 国产国语露脸激情在线看| 窝窝影院91人妻| 国产精品综合久久久久久久免费 | 欧美日韩黄片免| 国产精品一区二区三区四区久久 | 黄片播放在线免费| 一二三四在线观看免费中文在| 搡老岳熟女国产| 免费观看人在逋| 性欧美人与动物交配| 国产熟女xx| 国产野战对白在线观看| 亚洲性夜色夜夜综合| 亚洲人成伊人成综合网2020| av国产精品久久久久影院| 国产成人欧美| www.999成人在线观看| 亚洲国产欧美一区二区综合| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品合色在线| 日韩欧美一区视频在线观看| 欧美在线一区亚洲| 国产亚洲欧美精品永久| 欧美性长视频在线观看| 亚洲自拍偷在线| 欧美成人性av电影在线观看| 欧美成人免费av一区二区三区| 亚洲自偷自拍图片 自拍| 亚洲精品一区av在线观看| 亚洲九九香蕉| 可以免费在线观看a视频的电影网站| 成人永久免费在线观看视频| 欧美+亚洲+日韩+国产| 黑丝袜美女国产一区| 国产一区二区激情短视频| 超色免费av| 国产高清videossex| 一级a爱视频在线免费观看| 丰满的人妻完整版| 男人的好看免费观看在线视频 | 男男h啪啪无遮挡| 精品电影一区二区在线| 免费高清视频大片| 99香蕉大伊视频| 两个人看的免费小视频| 亚洲精品国产区一区二| 国产黄a三级三级三级人| 精品国产一区二区久久| 精品久久久精品久久久| 久久国产精品男人的天堂亚洲| 夜夜爽天天搞| 国产在线精品亚洲第一网站| 日韩免费av在线播放| 精品国产乱码久久久久久男人| 色综合欧美亚洲国产小说| 老熟妇乱子伦视频在线观看| 一进一出好大好爽视频| 男人操女人黄网站| 国产在线精品亚洲第一网站| 日韩精品青青久久久久久| 一夜夜www| 多毛熟女@视频| 亚洲一区二区三区色噜噜 | 欧美日韩视频精品一区| 亚洲av片天天在线观看| 国产真人三级小视频在线观看| 俄罗斯特黄特色一大片| 久久草成人影院| 在线观看66精品国产| 亚洲熟妇熟女久久| 国产亚洲欧美98| 亚洲专区国产一区二区| 人人妻人人爽人人添夜夜欢视频| 午夜精品久久久久久毛片777| netflix在线观看网站| 午夜福利欧美成人| 久久性视频一级片| 国产1区2区3区精品| 精品久久久久久,| a在线观看视频网站| 91精品国产国语对白视频| 亚洲美女黄片视频| 黄频高清免费视频| 国产av一区在线观看免费| 一进一出抽搐gif免费好疼 | 日日摸夜夜添夜夜添小说| 久久狼人影院| 久久性视频一级片| 亚洲精品av麻豆狂野| 日本欧美视频一区| 18禁裸乳无遮挡免费网站照片 | 亚洲第一青青草原| 国产成人av教育| 国产乱人伦免费视频| 亚洲精品美女久久久久99蜜臀| 久久久久久亚洲精品国产蜜桃av| 九色亚洲精品在线播放| www.999成人在线观看| 久久精品国产99精品国产亚洲性色 | 欧美日韩亚洲综合一区二区三区_| 人妻久久中文字幕网| 一区二区三区激情视频| 露出奶头的视频| 人人妻人人澡人人看| 欧美日韩av久久| 成年人免费黄色播放视频| 1024视频免费在线观看| 在线播放国产精品三级| 老司机靠b影院| 少妇裸体淫交视频免费看高清 | avwww免费| 久久精品91无色码中文字幕| 国产欧美日韩精品一区二区| 美女免费视频网站| 久久午夜亚洲精品久久| 搞女人的毛片| 亚洲精品色激情综合| 麻豆av噜噜一区二区三区| 九九在线视频观看精品| 综合色av麻豆| 欧美高清性xxxxhd video| 每晚都被弄得嗷嗷叫到高潮| 18禁在线播放成人免费| 男女那种视频在线观看| 日韩成人在线观看一区二区三区| 久久久久九九精品影院| 久久久久免费精品人妻一区二区| 九九久久精品国产亚洲av麻豆| 精品人妻熟女av久视频| 网址你懂的国产日韩在线| 哪里可以看免费的av片| 两人在一起打扑克的视频| 麻豆成人av在线观看| 乱人视频在线观看| 乱码一卡2卡4卡精品| 人人妻,人人澡人人爽秒播| 97热精品久久久久久| 波多野结衣巨乳人妻| 男女做爰动态图高潮gif福利片| 简卡轻食公司| .国产精品久久| 色综合站精品国产| 高清在线国产一区| 午夜福利在线在线| 午夜免费激情av| 亚洲,欧美,日韩| 国产成年人精品一区二区| 99久久精品国产亚洲精品| 久久香蕉精品热| 制服丝袜大香蕉在线| 国内精品久久久久精免费| 亚洲片人在线观看| 亚洲人与动物交配视频| 成年版毛片免费区| 琪琪午夜伦伦电影理论片6080| 日韩欧美在线二视频| 欧美又色又爽又黄视频| a在线观看视频网站| 天堂av国产一区二区熟女人妻| www日本黄色视频网| 亚洲成人久久爱视频| 在线a可以看的网站| 舔av片在线| 十八禁国产超污无遮挡网站| 美女免费视频网站| 午夜激情福利司机影院| 九九在线视频观看精品| 亚洲av第一区精品v没综合| 丁香六月欧美| 我的女老师完整版在线观看| 免费看日本二区| 成人毛片a级毛片在线播放| 婷婷六月久久综合丁香| 国产真实乱freesex| 欧美日韩综合久久久久久 | 亚洲美女黄片视频| 麻豆国产av国片精品| 成年版毛片免费区| 久久久精品欧美日韩精品| 夜夜躁狠狠躁天天躁| 亚洲成人中文字幕在线播放| 久久久国产成人免费| 国产黄片美女视频| 欧美性猛交╳xxx乱大交人| 十八禁国产超污无遮挡网站| 五月伊人婷婷丁香| x7x7x7水蜜桃| 亚洲第一区二区三区不卡| 观看美女的网站| 窝窝影院91人妻| 国产免费av片在线观看野外av| 桃色一区二区三区在线观看| 午夜福利在线观看免费完整高清在 | 草草在线视频免费看| 色视频www国产| 亚洲中文字幕一区二区三区有码在线看| 女人十人毛片免费观看3o分钟| 午夜日韩欧美国产| 精品午夜福利在线看| 国内久久婷婷六月综合欲色啪| 午夜福利18| 一进一出抽搐动态| 成人特级黄色片久久久久久久| 免费搜索国产男女视频| 欧美日韩国产亚洲二区| 有码 亚洲区| 十八禁国产超污无遮挡网站| 日韩欧美一区二区三区在线观看| 99久久精品国产亚洲精品| 悠悠久久av| 久久久久国内视频| 国产精品一区二区三区四区免费观看 | xxxwww97欧美| 日本免费一区二区三区高清不卡| 桃红色精品国产亚洲av| 又黄又爽又免费观看的视频| 在线观看美女被高潮喷水网站 | 成人欧美大片| 窝窝影院91人妻| 亚州av有码| 直男gayav资源| 国产不卡一卡二| 久久草成人影院| 国产三级在线视频| 非洲黑人性xxxx精品又粗又长| 永久网站在线| 国产成人影院久久av| h日本视频在线播放| 嫩草影院新地址| 国产精品不卡视频一区二区 | 99热精品在线国产| 88av欧美| 最近最新中文字幕大全电影3| 3wmmmm亚洲av在线观看| 欧美色欧美亚洲另类二区| 久久精品人妻少妇| 午夜精品在线福利| 中文亚洲av片在线观看爽| 性插视频无遮挡在线免费观看| 在线十欧美十亚洲十日本专区| 在线观看美女被高潮喷水网站 | 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 欧美日韩黄片免| 成年人黄色毛片网站| 熟女电影av网| a在线观看视频网站| 亚洲专区国产一区二区| 老鸭窝网址在线观看| 成人高潮视频无遮挡免费网站| 极品教师在线免费播放| 波多野结衣巨乳人妻| 久久久久久久午夜电影| 美女被艹到高潮喷水动态| 亚洲激情在线av| 国产av一区在线观看免费| 99久久精品热视频| 国产精品久久电影中文字幕| 日韩中字成人| 特大巨黑吊av在线直播| 深夜a级毛片| 蜜桃亚洲精品一区二区三区| 亚洲成人精品中文字幕电影| 精品无人区乱码1区二区| 免费在线观看影片大全网站| 日韩高清综合在线| 国产亚洲精品久久久久久毛片| 色播亚洲综合网| 黄色配什么色好看| 嫩草影院新地址| 99国产精品一区二区三区| 少妇裸体淫交视频免费看高清| 在线观看舔阴道视频| 中文字幕熟女人妻在线| 好看av亚洲va欧美ⅴa在| 九九热线精品视视频播放| 日日夜夜操网爽| 成年女人永久免费观看视频| 免费人成在线观看视频色| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添小说| 欧美日韩中文字幕国产精品一区二区三区| 99久久精品一区二区三区| 久久久久久大精品| 9191精品国产免费久久| 久久久精品欧美日韩精品| 国产伦在线观看视频一区| 性色av乱码一区二区三区2| 天天一区二区日本电影三级| av视频在线观看入口| 亚洲av五月六月丁香网| 色播亚洲综合网| av天堂中文字幕网| 日本一二三区视频观看| 动漫黄色视频在线观看| 男人和女人高潮做爰伦理| 国产午夜福利久久久久久| 久9热在线精品视频| 亚洲 国产 在线| 亚洲人成伊人成综合网2020| av在线观看视频网站免费| 欧美黄色淫秽网站| 3wmmmm亚洲av在线观看| 免费黄网站久久成人精品 | 国产精品久久电影中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产亚洲av天美| 国产精品,欧美在线| 精品久久久久久久久久免费视频| 国产欧美日韩精品亚洲av| 国产精品野战在线观看| 国产精品爽爽va在线观看网站| 成人永久免费在线观看视频| 99国产综合亚洲精品| 成人高潮视频无遮挡免费网站| 免费搜索国产男女视频| 日本成人三级电影网站| 久久久久久九九精品二区国产| 1024手机看黄色片| 亚洲精品亚洲一区二区| 久久久久性生活片| 少妇人妻一区二区三区视频| 欧美日韩亚洲国产一区二区在线观看| 色噜噜av男人的天堂激情| 国产久久久一区二区三区| 97超视频在线观看视频| 国产精品野战在线观看| 黄色视频,在线免费观看| 人妻制服诱惑在线中文字幕| 嫩草影院入口| 午夜福利在线观看免费完整高清在 | a级毛片a级免费在线| 内射极品少妇av片p| 亚洲精华国产精华精| 成年女人看的毛片在线观看| 熟女人妻精品中文字幕| 午夜影院日韩av| 亚洲美女搞黄在线观看 | 亚洲一区二区三区不卡视频| 69av精品久久久久久| 全区人妻精品视频| 国产欧美日韩一区二区精品| 人人妻,人人澡人人爽秒播| 色综合欧美亚洲国产小说| 九九久久精品国产亚洲av麻豆| 国产又黄又爽又无遮挡在线| 久久婷婷人人爽人人干人人爱| 老司机午夜福利在线观看视频| 黄色女人牲交| 99热精品在线国产| 亚洲av一区综合| 亚洲精品一卡2卡三卡4卡5卡| 好男人在线观看高清免费视频|