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

    甲醇模塊化生產(chǎn)中分時(shí)儲(chǔ)熱系統(tǒng)的優(yōu)化設(shè)計(jì)

    2022-03-03 06:01:22孔昕山黃仁星康麗霞劉永忠
    化工學(xué)報(bào) 2022年2期
    關(guān)鍵詞:流股儲(chǔ)熱冷器

    孔昕山,黃仁星,康麗霞,2,3,劉永忠,2,3

    (1 西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院,陜西 西安 710049; 2 新能源系統(tǒng)工程與裝備陜西省高校工程研究中心,陜西 西安 710049;3陜西省能源化工過(guò)程強(qiáng)化重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710049)

    引 言

    隨著風(fēng)電和光伏發(fā)電裝機(jī)量的迅猛擴(kuò)增,波動(dòng)性和間歇性的可再生能源發(fā)電使得傳統(tǒng)電力系統(tǒng)面臨巨大壓力和挑戰(zhàn)[1-2]。近年來(lái),“電轉(zhuǎn)液”(power to liquid,PtL)技術(shù)為可再生能源的消納提供了新途徑,通過(guò)可再生能源發(fā)電電解水生產(chǎn)氫氣,通過(guò)二氧化碳加氫制備甲醇,不僅可以消納波動(dòng)的可再生能源,而且可以提供易于存貯和運(yùn)輸?shù)幕瘜W(xué)品[3-4]。利用可再生能源生產(chǎn)大宗化學(xué)品需要考慮可再生能源供給的間歇性和產(chǎn)品市場(chǎng)需求的波動(dòng)性對(duì)生產(chǎn)系統(tǒng)的影響[5]。

    可再生能源的波動(dòng)性和產(chǎn)品市場(chǎng)需求變化都將影響甲醇生產(chǎn)過(guò)程優(yōu)化設(shè)計(jì)和操作。傳統(tǒng)的規(guī)?;B續(xù)甲醇生產(chǎn)系統(tǒng)通常需要頻繁調(diào)整生產(chǎn)工況,操作難度大且生產(chǎn)裝置難以在最優(yōu)工況運(yùn)行[6-7]。而甲醇的模塊化生產(chǎn)過(guò)程則易于通過(guò)“數(shù)量放大”方式調(diào)整產(chǎn)品數(shù)量和生產(chǎn)地點(diǎn),可快速而靈活地適應(yīng)生產(chǎn)工況和產(chǎn)品市場(chǎng)的變化[8-9]。模塊化生產(chǎn)系統(tǒng)通過(guò)標(biāo)準(zhǔn)化生產(chǎn)設(shè)備、生產(chǎn)線或生產(chǎn)系統(tǒng),使生產(chǎn)系統(tǒng)便捷而快速地適應(yīng)工況變化,并使生產(chǎn)裝置發(fā)揮最佳效能[10-11]。因而小規(guī)模的模塊化生產(chǎn)在消納可再生能源的生產(chǎn)過(guò)程中較傳統(tǒng)大規(guī)模生產(chǎn)更具經(jīng)濟(jì)性和操作靈活性[12-13]。然而,在模塊化的甲醇生產(chǎn)系統(tǒng)中,生產(chǎn)模塊的頻繁啟停將導(dǎo)致生產(chǎn)系統(tǒng)中流股的預(yù)熱、冷卻和反應(yīng)熱移除等均隨時(shí)間波動(dòng)。為了提高生產(chǎn)系統(tǒng)的經(jīng)濟(jì)性需通過(guò)能量?jī)?chǔ)存和調(diào)度提升波動(dòng)工況條件下生產(chǎn)系統(tǒng)能量利用效率。

    對(duì)于波動(dòng)的供給與需求之間的熱量匹配問(wèn)題,設(shè)立儲(chǔ)熱系統(tǒng)是平抑波動(dòng)的主要途徑之一[14-15]。儲(chǔ)熱系統(tǒng)主要通過(guò)顯熱、潛熱和化學(xué)反應(yīng)熱儲(chǔ)存熱量[16]。Farzan等[17]研究了由吸熱能力強(qiáng)、通過(guò)顯熱儲(chǔ)熱的瀝青材料構(gòu)成的瀝青太陽(yáng)能空氣加熱器的熱動(dòng)力學(xué)和熱效率。Kumar 等[18]建立了立式圓柱形填充床潛熱儲(chǔ)熱系統(tǒng)的數(shù)值模型,分析了傳熱流體和封裝相變材料的瞬態(tài)變化,著重探究封裝形狀、孔隙率等對(duì)儲(chǔ)熱性能的影響。Liu 等[19]通過(guò)熱力學(xué)和動(dòng)力學(xué)評(píng)估高溫下金屬氧化物Mn2O3/Mn3O4和Co3O4/CoO 氧化還原反應(yīng)的化學(xué)反應(yīng)儲(chǔ)熱平衡溫度、反應(yīng)熱和反應(yīng)速率。當(dāng)供需熱量在時(shí)間和空間上不匹配時(shí),引入儲(chǔ)能介質(zhì)和儲(chǔ)罐可提高能量利用率,實(shí)現(xiàn)熱量的分時(shí)調(diào)度,即將時(shí)間和空間再分配[20-21]。例如,對(duì)于間歇性的太陽(yáng)能熱供應(yīng)和需求的匹配問(wèn)題,Wan 等[22]引入了雙罐熔融鹽儲(chǔ)能系統(tǒng),解決日與季節(jié)在時(shí)間尺度上不匹配的問(wèn)題。Lou等[23]考慮了冷熱流體的強(qiáng)混合干擾溫度分層,導(dǎo)致儲(chǔ)罐的熱性能降低的問(wèn)題,對(duì)單罐儲(chǔ)熱進(jìn)行了模擬研究。Mao 等[24]考慮了儲(chǔ)罐中相變材料PCM 的動(dòng)態(tài)儲(chǔ)能特性及其對(duì)系統(tǒng)性能的影響。對(duì)于利用太陽(yáng)能的儲(chǔ)熱系統(tǒng),由于其熱量的供需品位固定,一般僅需設(shè)立單個(gè)儲(chǔ)罐或高溫和低溫兩個(gè)溫位的儲(chǔ)罐進(jìn)行儲(chǔ)熱[25]。而對(duì)于以生產(chǎn)線為模塊的模塊化甲醇生產(chǎn)系統(tǒng),在生產(chǎn)線啟停過(guò)程中熱量的供給和需求為不同溫位和數(shù)量的冷熱流股,則需要采用多種品位的儲(chǔ)熱介質(zhì)和儲(chǔ)罐實(shí)現(xiàn)能量的分時(shí)儲(chǔ)存和調(diào)度,這就需要對(duì)儲(chǔ)熱儲(chǔ)罐的數(shù)量、儲(chǔ)熱儲(chǔ)罐的溫位和儲(chǔ)罐容量進(jìn)行優(yōu)化設(shè)計(jì)。

    針對(duì)上述問(wèn)題,本文以甲醇模塊化生產(chǎn)為例,提出了模塊化生產(chǎn)系統(tǒng)中分時(shí)儲(chǔ)熱系統(tǒng)的三步優(yōu)化設(shè)計(jì)方法,利用夾點(diǎn)法獲得的無(wú)分流換熱網(wǎng)絡(luò)為依據(jù)設(shè)立儲(chǔ)罐儲(chǔ)熱,通過(guò)儲(chǔ)熱介質(zhì)與工藝流股的換熱匹配關(guān)系確定儲(chǔ)罐的儲(chǔ)熱溫區(qū),建立基于分時(shí)儲(chǔ)熱的線性規(guī)劃模型,以最少的儲(chǔ)罐數(shù)量和最優(yōu)儲(chǔ)罐容量配置為目標(biāo),在滿足生產(chǎn)模塊啟停的能量需求條件下進(jìn)行分時(shí)儲(chǔ)熱系統(tǒng)的優(yōu)化設(shè)計(jì),同時(shí)最小化儲(chǔ)熱系統(tǒng)的投資費(fèi)用和操作費(fèi)用。通過(guò)熱量的分時(shí)調(diào)度,將前階段熱量存儲(chǔ)后供后續(xù)階段模塊啟動(dòng)使用,提高生產(chǎn)系統(tǒng)的能量利用率,為分時(shí)儲(chǔ)熱系統(tǒng)的優(yōu)化設(shè)計(jì)提供分析工具。

    1 模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng)

    模塊化化工生產(chǎn)過(guò)程與傳統(tǒng)的間歇化工過(guò)程既有聯(lián)系也有區(qū)別,兩者在能量利用方式上具有不同的特點(diǎn)。間歇化工過(guò)程中的能量供需在生產(chǎn)過(guò)程排產(chǎn)確定后是確定的,即生產(chǎn)周期一般是確定的;而模塊化化工生產(chǎn)過(guò)程需根據(jù)原料供給或產(chǎn)品需求的不確定性變化調(diào)整生產(chǎn)模塊的啟動(dòng)或停止,生產(chǎn)模塊的啟停與原料供應(yīng)或產(chǎn)品需求的周期規(guī)律相適應(yīng),具有一定的不確定性。同時(shí)間歇化工生產(chǎn)過(guò)程中流股的初始和目標(biāo)溫度、起止時(shí)間以及熱負(fù)荷是給定的;而在模塊化化工生產(chǎn)過(guò)程中,雖然單個(gè)生產(chǎn)模塊的排產(chǎn)是既定的,但生產(chǎn)模塊調(diào)度的數(shù)量將隨原料供給或產(chǎn)品需求的不確定性變化,即模塊啟停過(guò)程中可進(jìn)行熱集成流股的熱負(fù)荷會(huì)發(fā)生變化。對(duì)于間歇化工生產(chǎn)過(guò)程,出現(xiàn)在同一時(shí)間段的流股可采用直接換熱方式進(jìn)行匹配,對(duì)跨時(shí)段流股的換熱需采用儲(chǔ)熱介質(zhì)進(jìn)行間接換熱;而在模塊化化工生產(chǎn)中,生產(chǎn)模塊實(shí)質(zhì)是小規(guī)模連續(xù)化生產(chǎn),生產(chǎn)模塊的啟停過(guò)程為非穩(wěn)態(tài),所有流股均出現(xiàn)在同一時(shí)間段,難以動(dòng)態(tài)地進(jìn)行直接換熱,需采用儲(chǔ)熱介質(zhì)進(jìn)行間接換熱。

    化學(xué)品模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng)如圖1 所示。在化學(xué)品模塊化生產(chǎn)系統(tǒng)中,隨著原料供應(yīng)和產(chǎn)品需求的變化,在不同的時(shí)間區(qū)段,生產(chǎn)系統(tǒng)由生產(chǎn)能力相同而數(shù)量變化的生產(chǎn)線為模塊構(gòu)成,其中反應(yīng)熱的移除采用分時(shí)儲(chǔ)熱的策略。根據(jù)生產(chǎn)需要,整個(gè)生產(chǎn)時(shí)間劃分為K 個(gè)生產(chǎn)周期,縱坐標(biāo)NU為當(dāng)前周期與上一周期相比生產(chǎn)線變化數(shù)量的絕對(duì)值,即生產(chǎn)線啟停的數(shù)量。系統(tǒng)中熱量的分時(shí)分配由儲(chǔ)熱介質(zhì)承擔(dān)。由圖可見(jiàn),在生產(chǎn)線的啟停過(guò)程中生產(chǎn)處于非穩(wěn)態(tài),工藝流股間無(wú)法直接換熱,需通過(guò)儲(chǔ)熱介質(zhì)將生產(chǎn)線啟停過(guò)程中的熱量?jī)?chǔ)存在相應(yīng)儲(chǔ)罐中,同時(shí)儲(chǔ)罐也將通過(guò)儲(chǔ)熱介質(zhì)向生產(chǎn)線供給所需的熱量,從而實(shí)現(xiàn)模塊化生產(chǎn)的分時(shí)儲(chǔ)熱和熱量調(diào)度。此外,在該儲(chǔ)熱系統(tǒng)中,每個(gè)儲(chǔ)罐配備換熱器,儲(chǔ)罐中儲(chǔ)熱介質(zhì)的熱量可根據(jù)需要廢棄,也可通過(guò)換熱器轉(zhuǎn)運(yùn)到其他溫區(qū)的儲(chǔ)罐中。

    圖1 化學(xué)品模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng)Fig.1 Time-sharing heat storage system for modular production of chemicals

    針對(duì)上述模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng),給定:(1)生產(chǎn)線的工藝參數(shù);(2)模塊化生產(chǎn)中模塊數(shù)量隨周期的調(diào)度情況;(3)系統(tǒng)中儲(chǔ)罐、換熱設(shè)備等組件的操作參數(shù)和儲(chǔ)能介質(zhì)的物性參數(shù)。本文將通過(guò)儲(chǔ)罐設(shè)置及流股匹配、儲(chǔ)罐儲(chǔ)熱溫區(qū)確定和儲(chǔ)罐容量配置及調(diào)度三個(gè)步驟,確定分時(shí)儲(chǔ)熱系統(tǒng)中儲(chǔ)罐的優(yōu)化設(shè)計(jì)和運(yùn)行方案,以解決波動(dòng)條件下模塊化生產(chǎn)過(guò)程中能量的儲(chǔ)存和調(diào)度問(wèn)題,提升系統(tǒng)的能量利用效率。

    2 模塊化生產(chǎn)中分時(shí)儲(chǔ)熱系統(tǒng)的優(yōu)化設(shè)計(jì)

    2.1 儲(chǔ)熱系統(tǒng)的優(yōu)化設(shè)計(jì)策略

    對(duì)以生產(chǎn)線為模塊的化學(xué)品模塊化生產(chǎn)系統(tǒng),本文提出的分時(shí)儲(chǔ)熱系統(tǒng)優(yōu)化設(shè)計(jì)策略分為以下三個(gè)步驟。

    步驟1:儲(chǔ)罐設(shè)置及流股匹配。確定分時(shí)儲(chǔ)熱系統(tǒng)中儲(chǔ)罐的數(shù)量和儲(chǔ)熱介質(zhì)與工藝流股的換熱匹配關(guān)系。

    步驟2:儲(chǔ)罐儲(chǔ)熱溫區(qū)確定。根據(jù)換熱匹配關(guān)系確定相應(yīng)的儲(chǔ)罐儲(chǔ)熱溫度約束關(guān)系,進(jìn)而確定儲(chǔ)罐的儲(chǔ)熱溫區(qū)。

    步驟3:儲(chǔ)罐容量配置及調(diào)度。建立以儲(chǔ)熱系統(tǒng)總費(fèi)用最小為目標(biāo)的儲(chǔ)罐容量配置優(yōu)化數(shù)學(xué)模型,求解模型并獲得儲(chǔ)罐容量的最優(yōu)配置以及熱量調(diào)度方案。

    以下將詳述這三個(gè)步驟。

    2.2 步驟1:儲(chǔ)罐設(shè)置及流股匹配

    在確定分時(shí)儲(chǔ)熱系統(tǒng)中儲(chǔ)罐的數(shù)量以及儲(chǔ)熱介質(zhì)與流股的換熱匹配關(guān)系時(shí),本文基于連續(xù)生產(chǎn)中單條生產(chǎn)線的無(wú)分流換熱網(wǎng)絡(luò)。為表述方便,定義以下集合:h∈[1,H]表示換熱網(wǎng)絡(luò)中的熱流股;c∈[1,C]表示換熱網(wǎng)絡(luò)中的冷流股;i∈[1,I]表示劃分流股的組別;j∈[1,Ji],?i∈[1,I]表示第i組中的流股;其中H、C、I、Ji分別表示熱流股、冷流股、劃分流股組別和第i組中的流股的數(shù)量。定義以下參數(shù):Sij表示第i組中第j條流股;STij表示第i組中第j個(gè)儲(chǔ)罐。

    圖2給出儲(chǔ)熱儲(chǔ)罐的設(shè)置流程。圖中熱流股尋找方向?yàn)榻禍胤较?,冷流股尋找方向?yàn)樯郎胤较?,流股目?biāo)溫度即為換熱網(wǎng)絡(luò)中換熱器的流股出口溫度,終溫即為換熱網(wǎng)絡(luò)中流股的最終溫度。該策略的儲(chǔ)熱方法以換熱網(wǎng)絡(luò)為依據(jù),目的在于盡可能合并換熱以減少儲(chǔ)罐的數(shù)量。按照?qǐng)D2 中流程使用的分時(shí)儲(chǔ)熱方法,每進(jìn)行一次換熱需要設(shè)立相應(yīng)的儲(chǔ)熱儲(chǔ)罐,將流股劃分為I 組,每組包括Ji條流股,則

    圖2 儲(chǔ)熱儲(chǔ)罐設(shè)置流程Fig.2 Heat storage tank setup flowchart

    需要注意的是,同一條流股可以劃分在不同的組中,同一組中只能使用同一種介質(zhì),種類應(yīng)根據(jù)工藝條件選擇。在得到的每組流股中儲(chǔ)熱介質(zhì)的流向?yàn)橐粋€(gè)循環(huán),從上一級(jí)儲(chǔ)罐中流出,與對(duì)應(yīng)的工藝流股換熱后流入下一級(jí)儲(chǔ)罐。最末一級(jí)儲(chǔ)罐中的介質(zhì)與對(duì)應(yīng)工藝流股換熱后流入第一級(jí)儲(chǔ)罐。

    2.3 步驟2:儲(chǔ)罐儲(chǔ)熱溫區(qū)確定

    通過(guò)步驟1確定了工藝流程中所需儲(chǔ)罐的數(shù)量以及每個(gè)儲(chǔ)罐中儲(chǔ)存介質(zhì)與工藝流股換熱的匹配關(guān)系。為保持分時(shí)儲(chǔ)熱過(guò)程的穩(wěn)定性,令換熱過(guò)程中儲(chǔ)熱介質(zhì)的進(jìn)出口溫度和儲(chǔ)罐儲(chǔ)熱溫度不隨周期發(fā)生變化。對(duì)任意的第i組流股,若其主流股Si1為熱流股,則其余流股為冷流股,在換熱過(guò)程中傳熱溫差需不小于最小傳熱溫差,儲(chǔ)熱溫度滿足以下約束

    2.4 步驟3:儲(chǔ)罐容量配置及調(diào)度

    2.4.1 目標(biāo)函數(shù) 為了確定儲(chǔ)熱系統(tǒng)中儲(chǔ)罐的容量和操作中熱量的最優(yōu)調(diào)度,本節(jié)中以儲(chǔ)熱系統(tǒng)總費(fèi)用最小為目標(biāo)對(duì)儲(chǔ)熱系統(tǒng)進(jìn)行優(yōu)化。該儲(chǔ)熱系統(tǒng)的成本主要包括儲(chǔ)罐、換熱器的投資和維護(hù)費(fèi)用,以及換熱器的操作費(fèi)用??傎M(fèi)用可表示為

    式中,k∈[1, K] 表示周期;Δk為周期的時(shí)間長(zhǎng)度;Ctotal為系統(tǒng)的總費(fèi)用;Ccap為各設(shè)備的投資費(fèi)用;Cmai為各設(shè)備的維護(hù)費(fèi)用;HEij為第i組中第j個(gè)換熱器;CRF為各設(shè)備的投資回收因子;OPHEij為換熱器HEij的操作費(fèi)用;PHEijk為第k周期內(nèi)換熱設(shè)備HEij的平均功率。

    式中,r為利率;n為各設(shè)備的使用壽命。

    2.4.2 約束條件 對(duì)任意第i組流股,以主流股Si1為熱流股為例,其余流股為冷流股。則儲(chǔ)罐廢棄熱量需通過(guò)冷卻器冷卻至最后一級(jí)儲(chǔ)罐儲(chǔ)存溫度,所需配置的換熱器為冷卻器,需滿足以下約束。

    (1)儲(chǔ)熱介質(zhì)的能量守恒 任意周期內(nèi),儲(chǔ)熱介質(zhì)與對(duì)應(yīng)流股換熱過(guò)程滿足能量守恒,即

    式中,F(xiàn)ij,k表示在第k周期內(nèi)與第i組中第j條流股換熱的儲(chǔ)熱介質(zhì)流率;cpi表示在第i組中使用的儲(chǔ)熱介質(zhì)的比熱容;NUk表示在第k周期內(nèi)生產(chǎn)線的啟停數(shù)量;CPij,k表示第i組中第j條流股在第k周期內(nèi)的平均熱容流率。

    (2)儲(chǔ)熱儲(chǔ)罐的質(zhì)量守恒 儲(chǔ)罐中儲(chǔ)熱介質(zhì)為下一周期供能,并且多余熱量可以廢棄,由冷卻器冷卻后輸入最后一個(gè)儲(chǔ)罐,其輸入和輸出流率滿足

    除了為后續(xù)階段生產(chǎn)線的啟停供能外,儲(chǔ)罐中的熱量還可以用于其他適用途徑,如產(chǎn)生蒸汽或作為熱泵和熱機(jī)的熱源等,若無(wú)適用途徑,則通過(guò)冷卻器冷卻廢棄。其負(fù)荷可以表示為

    對(duì)任意的第i組流股,若主流股Si1為冷流股,其余流股為熱流股,則儲(chǔ)罐廢棄能量需通過(guò)加熱器加熱至最后一級(jí)儲(chǔ)罐儲(chǔ)存溫度,所需配置的換熱器為加熱器。

    式(7)~式(19)即構(gòu)成了一個(gè)線性規(guī)劃模型,適用于模塊化化工生產(chǎn)中分時(shí)儲(chǔ)熱系統(tǒng)優(yōu)化。

    3 甲醇模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng)分析

    3.1 甲醇模塊化生產(chǎn)基礎(chǔ)參數(shù)分析

    本文將上述分時(shí)儲(chǔ)熱系統(tǒng)的優(yōu)化設(shè)計(jì)策略應(yīng)用于可再生能源驅(qū)動(dòng)的甲醇模塊化生產(chǎn)系統(tǒng)。該甲醇模塊化生產(chǎn)系統(tǒng)以生產(chǎn)線為模塊,根據(jù)可再生能源供給和產(chǎn)品需求的變化,投入的生產(chǎn)模塊數(shù)量進(jìn)行相應(yīng)變化。圖3給出甲醇模塊化生產(chǎn)中單條生產(chǎn)線的簡(jiǎn)化工藝流程[26]。CO2和H2按一定比例混合后升溫加壓至反應(yīng)條件,在反應(yīng)器中反應(yīng)后出料閃蒸分離,氣相回流,液相精制后得出產(chǎn)物甲醇。

    圖3 甲醇模塊化生產(chǎn)中單生產(chǎn)模塊的簡(jiǎn)化工藝流程Fig.3 Simplified process flow of single production module in methanol modular production

    本文計(jì)算的基準(zhǔn)時(shí)間取30 天,圖4 給出其中一個(gè)典型天中24 h 內(nèi)生產(chǎn)線數(shù)量的調(diào)度情況。在30天內(nèi)由于氣候條件和市場(chǎng)情況相似,可認(rèn)為可再生能源與下游市場(chǎng)需求的波動(dòng)特性相同,即生產(chǎn)線的調(diào)度情況在每天內(nèi)的變化特性相同。由于在生產(chǎn)線穩(wěn)態(tài)運(yùn)行時(shí)冷熱流股可以互相匹配換熱,因此在分時(shí)儲(chǔ)熱系統(tǒng)中,針對(duì)熱量的儲(chǔ)存與調(diào)度僅發(fā)生在生產(chǎn)線的啟停過(guò)程。每發(fā)生一次生產(chǎn)線數(shù)量的變化記為一個(gè)周期,取生產(chǎn)線由未工作狀態(tài)啟動(dòng)至穩(wěn)態(tài)運(yùn)行狀態(tài)的時(shí)間,以及由穩(wěn)態(tài)運(yùn)行狀態(tài)關(guān)停至未工作狀態(tài)的時(shí)間,即生產(chǎn)線啟停時(shí)間為1 h,周期時(shí)間長(zhǎng)度Δk=1 h。由圖4可知,在一天內(nèi)生產(chǎn)線調(diào)度8次,即1天為8個(gè)周期,30天為240個(gè)周期,K=240。

    圖4 一個(gè)典型天24 h內(nèi)生產(chǎn)線數(shù)量的調(diào)度情況Fig.4 Scheduling situation of the number of production lines for 24 h in a typical day

    圖5給出了該工藝流程中公用工程換熱器的熱負(fù)荷。式(9)和式(10)右邊的流股熱量變化項(xiàng)即是相應(yīng)換熱器的熱負(fù)荷,即

    圖5 一個(gè)典型天24 h內(nèi)生產(chǎn)線換熱器熱負(fù)荷變化情況Fig.5 Heat exchanger heat duty variation of the production line for 24 h in a typical day

    3.2 甲醇模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng)優(yōu)化設(shè)計(jì)過(guò)程

    3.2.1 儲(chǔ)罐的數(shù)量以及儲(chǔ)熱介質(zhì)與流股的換熱匹

    配關(guān)系 針對(duì)甲醇的模塊化生產(chǎn)系統(tǒng),本文采用2.2節(jié)的方法確定分時(shí)儲(chǔ)熱系統(tǒng)的儲(chǔ)罐數(shù)量。圖6給出單一模塊甲醇生產(chǎn)流程的無(wú)分流換熱網(wǎng)絡(luò)。由圖可見(jiàn),除一股冷卻公用工程(冷卻水)外,該換熱網(wǎng)絡(luò)包括1 股熱流股和3 股冷流股。所需換熱最多的流股為熱流股,確定為主流股,與其他3股冷流股可劃分為1個(gè)組,共包含4股流股,需設(shè)立4個(gè)儲(chǔ)罐。

    圖6 甲醇生產(chǎn)工藝流程的無(wú)分流換熱網(wǎng)絡(luò)Fig.6 Heat exchanger network without stream splits for methanol production process

    為了使儲(chǔ)熱介質(zhì)用量盡可能少,可盡可能取較大的進(jìn)出口溫差。同時(shí),為使儲(chǔ)熱介質(zhì)在4 個(gè)儲(chǔ)罐間的流量接近,避免儲(chǔ)熱介質(zhì)的累積,ST11, ST12,ST13, ST14的儲(chǔ)熱溫度分別取為296.85℃, 291.85℃,154.75℃,105℃。

    3.2.3 儲(chǔ)罐的最優(yōu)容量配置 本節(jié)采用2.4 節(jié)中所建立的線性規(guī)劃模型對(duì)儲(chǔ)罐的容量進(jìn)行優(yōu)化。對(duì)于儲(chǔ)罐廢棄熱量所配備換熱器選擇,由于該案例中主流股為熱流股,儲(chǔ)罐所需配備換熱器為冷卻器,而儲(chǔ)罐的儲(chǔ)熱溫度遠(yuǎn)高于環(huán)境溫度,且廢棄熱量換熱溫差大,因此可采用空冷器,避免使用冷卻公用工程。儲(chǔ)熱介質(zhì)工作溫度在30~296.85℃之間,可采用Dynalene 公司的導(dǎo)熱油Dynalene MT,該導(dǎo)熱油可在0~350℃間工作并保持熱穩(wěn)定性。本文中年利率取為5%。表1 給出儲(chǔ)罐和空冷器的基本參數(shù)。

    表1 儲(chǔ)罐和空冷器的費(fèi)用計(jì)算參數(shù)Table 1 Cost calculation data for storage tanks and air coolers

    本文在PYOMO平臺(tái)上求解該模型,采用Gurobi求解器,相對(duì)誤差為1%,計(jì)算機(jī)CPU 為Intel(R)Core(TM)i5-10400 CPU@2.90GHz。

    3.3 計(jì)算結(jié)果的分析與討論

    在本文構(gòu)建的分時(shí)儲(chǔ)熱系統(tǒng)中,為了討論儲(chǔ)熱熱量的廢棄對(duì)分時(shí)儲(chǔ)熱系統(tǒng)總費(fèi)用的影響,本文分為以下兩種情形進(jìn)行分析和討論。

    情形S1:儲(chǔ)熱熱量完全利用,無(wú)熱量廢棄。

    情形S2:儲(chǔ)熱熱量盡可能利用,但允許熱量廢棄。

    3.3.1 設(shè)計(jì)特性分析 表2給出了兩種情形下甲醇模塊化生產(chǎn)的分時(shí)儲(chǔ)熱系統(tǒng)中儲(chǔ)罐容量配置結(jié)果。對(duì)比情形S1,情形S2 中儲(chǔ)罐ST12容量下降87.68%,儲(chǔ)罐ST13容量下降99.04%,儲(chǔ)罐ST14容量下降99.06%,這說(shuō)明當(dāng)允許系統(tǒng)廢棄熱量時(shí),儲(chǔ)罐存儲(chǔ)的熱量只需滿足后續(xù)階段生產(chǎn)的熱量需求,不需要過(guò)剩熱量的儲(chǔ)存,因而情形S2中儲(chǔ)罐的容量配置顯著下降。然而,在兩種情形下,儲(chǔ)罐ST11容量的配置并未發(fā)生明顯變化。這是由于熱量的廢棄將引起額外的空冷器投資費(fèi)用和操作費(fèi)用,而儲(chǔ)罐ST11儲(chǔ)存介質(zhì)的溫度較高,換熱溫差大,對(duì)應(yīng)空冷器的操作負(fù)荷大,其操作費(fèi)用高于其余儲(chǔ)罐,因此,系統(tǒng)會(huì)優(yōu)先選擇存儲(chǔ)熱量而不是廢棄熱量。

    表2 情形S1和情形S2中儲(chǔ)罐容量配置Table 2 The capacity configuration of storage tanks in Scenario 1 and Scenario 2

    此外,在兩種情形下,儲(chǔ)罐ST11的容量都要遠(yuǎn)大于儲(chǔ)罐ST12、ST13和ST14。這是由于在每組流股中,儲(chǔ)熱介質(zhì)從主流股中存儲(chǔ)的熱量要供給其余流股,其對(duì)應(yīng)的儲(chǔ)罐儲(chǔ)熱需求量大,導(dǎo)致其容量配置較大。

    表3 給出兩種情形下的空冷器容量配置結(jié)果。在情形S1 中,系統(tǒng)僅對(duì)儲(chǔ)罐ST14配置了空冷器。這是由于為了使儲(chǔ)熱介質(zhì)與每股流股換熱后恢復(fù)到初始溫度,需要為最后一個(gè)儲(chǔ)罐配置空冷器,使介質(zhì)的溫度滿足換熱時(shí)的溫差約束。在情形S2中,系統(tǒng)對(duì)每個(gè)儲(chǔ)罐配置了空冷器(即是下文中的S2(1234)情形),允許儲(chǔ)罐中的熱量廢棄。

    表3 情形S1和情形S2中空冷器熱負(fù)荷配置Table 3 Heat exchanger heat duty configuration in Scenario 1 and Scenario 2

    表4給出了兩種情形下系統(tǒng)的設(shè)備費(fèi)用與操作費(fèi)用計(jì)算結(jié)果。在情形S1中,儲(chǔ)罐投資費(fèi)用占比總費(fèi)用94.66%,空冷器投資費(fèi)用占比0.68%,操作費(fèi)用占比4.66%。而在情形S2 中,儲(chǔ)罐投資費(fèi)用占比下降到69.77%,空冷器投資費(fèi)用占比2.41%,操作費(fèi)用占比27.82%。對(duì)比情形S1,情形S2 中系統(tǒng)儲(chǔ)熱系統(tǒng)總費(fèi)用下降了40.38%,儲(chǔ)罐投資費(fèi)用下降了56.06%,空冷器投資費(fèi)用增加了111.27%,操作費(fèi)用增加了255.68%。

    表4 情形S1和情形S2中設(shè)備投資費(fèi)用與操作費(fèi)用Table 4 Equipment costs and operating costs in Scenario 1 and Scenario 2

    為權(quán)衡允許廢棄熱量所造成的空冷器投資和操作費(fèi)用與系統(tǒng)節(jié)省儲(chǔ)罐費(fèi)用的關(guān)系,本文依據(jù)設(shè)立空冷器的不同組合,可將情形S2 分為S2(14)、S2(24)、S2(34)、S2(124)、S2(134)、S2(234)和S2(1234)7 種子情形,S2(1234)即是上文中的情形S2。例如,S2(24)表示僅設(shè)立空冷器HE12和HE14,而S2(234)則表示設(shè)置空冷器HE12、HE13和HE14,依此類推。在這7種子情形中,空冷器HE14必須設(shè)置。本文將討論空冷器組合設(shè)置對(duì)系統(tǒng)費(fèi)用的影響。圖7給出了包括情形S1 內(nèi)共8 種情形下系統(tǒng)設(shè)備費(fèi)用與操作費(fèi)用對(duì)比情況。

    圖7 不同情形下系統(tǒng)費(fèi)用的對(duì)比Fig.7 Comparison of the cost of the system in different scenarios

    由圖可見(jiàn),在系統(tǒng)總費(fèi)用中儲(chǔ)罐的投資費(fèi)用占比最大,其次是空冷器操作費(fèi)用,空冷器投資費(fèi)用最小。因此,儲(chǔ)罐容量差異引起的儲(chǔ)罐投資費(fèi)用的變化對(duì)系統(tǒng)總費(fèi)用的影響最大。在上述各情形中,所有儲(chǔ)罐均設(shè)立空冷器時(shí)(即情形S2(1234))的系統(tǒng)總費(fèi)用最小,空冷器投資和操作費(fèi)用最大。而僅設(shè)立空冷器HE14時(shí)(即情形S1)系統(tǒng)總費(fèi)用最大,空冷器投資和操作費(fèi)用最小,其余情形均處于這兩種情形之間。上述變化趨勢(shì)說(shuō)明,當(dāng)系統(tǒng)允許儲(chǔ)存的熱量廢棄,勢(shì)必會(huì)加大空冷器的負(fù)荷,從而增加額外的空冷器投資費(fèi)用和操作費(fèi)用,同時(shí)也會(huì)降低儲(chǔ)罐的存儲(chǔ)容量,從而降低儲(chǔ)罐投資費(fèi)用。

    對(duì)比情形S2(24)、S2(34)和S2(14)下的費(fèi)用組成可以看出,S2(24)情形下系統(tǒng)總費(fèi)用最大,儲(chǔ)罐投資費(fèi)用最大,即設(shè)立空冷器HE12對(duì)系統(tǒng)總費(fèi)用的影響最小,其原因在于:在四個(gè)儲(chǔ)罐中,ST12容量配置最小,可節(jié)省的費(fèi)用最少。設(shè)立空冷器HE11對(duì)系統(tǒng)總費(fèi)用影響略小于HE13,這是由于雖然設(shè)立HE11可大幅降低儲(chǔ)罐ST11的容量配置,節(jié)省較多的儲(chǔ)罐成本,但儲(chǔ)罐ST11較高的儲(chǔ)存溫度同時(shí)也造成了更大的空冷器操作負(fù)荷,引起額外的操作費(fèi)用,綜合兩者,設(shè)立HE13對(duì)總費(fèi)用的影響最大。這與表3中空冷器熱負(fù)荷配置結(jié)果相吻合。

    上述結(jié)果表明,在儲(chǔ)罐后設(shè)置空冷器廢棄熱量可以降低相應(yīng)儲(chǔ)罐的容量配置,降低儲(chǔ)罐的投資費(fèi)用,但同時(shí)也會(huì)造成額外的空冷器投資和操作費(fèi)用。對(duì)于儲(chǔ)熱系統(tǒng)而言,需要在熱量的廢棄與存儲(chǔ)間進(jìn)行權(quán)衡,這取決于儲(chǔ)罐可減少的容量配置、儲(chǔ)罐的存儲(chǔ)溫度以及空冷器的熱負(fù)荷配置等因素。

    3.3.2 儲(chǔ)罐運(yùn)行特性分析 圖8 給出在S1 和S2 兩種情形下儲(chǔ)熱系統(tǒng)中儲(chǔ)罐的儲(chǔ)能過(guò)程特性。由圖中的局部放大圖可見(jiàn),在情形S1 和情形S2 中,所有儲(chǔ)罐的負(fù)荷因子都呈現(xiàn)出了明顯的周期性變化規(guī)律,每8 個(gè)小周期的變化形成了一個(gè)大周期。這樣的變化規(guī)律與3.1 節(jié)中一個(gè)典型天為8 個(gè)周期的設(shè)定相吻合。

    圖8 情形S1和情形S2儲(chǔ)罐的儲(chǔ)能過(guò)程特性Fig.8 Energy storage process of storage tanks in Scenario 1 and Scenario 2

    在情形S1 中,在整體上4 個(gè)儲(chǔ)罐的負(fù)荷因子都接近線性變化,其中儲(chǔ)罐ST11、ST12和ST13的負(fù)荷因子隨周期延長(zhǎng)逐漸增大,而儲(chǔ)罐ST14的負(fù)荷因子隨周期延長(zhǎng)逐漸減小。這是由于在情形S1 中儲(chǔ)罐中的熱量不允許廢棄,熱量的調(diào)度僅取決于生產(chǎn)中工藝流股熱量的供給量和需求量,儲(chǔ)罐負(fù)荷因子的變化規(guī)律在每個(gè)大周期內(nèi)是固定的,在整體上的變化趨勢(shì)取決于在每個(gè)大周期內(nèi)的變化。其中儲(chǔ)罐的容量配置取決于整個(gè)過(guò)程中負(fù)荷最大的時(shí)刻,這就導(dǎo)致儲(chǔ)罐的利用率相對(duì)較低。

    在情形S2 中,儲(chǔ)罐ST12、ST13和ST14的負(fù)荷因子在整體上隨周期呈現(xiàn)波動(dòng),在最小和最大負(fù)荷狀態(tài)下波動(dòng)。這是由于在情形S2 中儲(chǔ)罐中允許廢棄熱量,系統(tǒng)在熱量滿足需求后優(yōu)先選擇廢棄熱量以減少存量,儲(chǔ)罐中存儲(chǔ)的熱量被頻繁調(diào)度。與情形S1相比,儲(chǔ)罐的利用率明顯提升。儲(chǔ)罐ST11的負(fù)荷因子變化規(guī)律與情形S1 中相同,這與3.3.1 節(jié)中的分析結(jié)論相吻合。

    上述研究表明,允許熱量的廢棄可以在滿足需求的前提下減少儲(chǔ)罐中儲(chǔ)熱量,提升儲(chǔ)罐的利用率,減小儲(chǔ)罐的容量配置。從情形S1中儲(chǔ)罐的負(fù)荷因子變化趨勢(shì)也可以推斷,當(dāng)計(jì)算的基準(zhǔn)時(shí)間尺度更大時(shí),其線性變化趨勢(shì)也會(huì)隨之延長(zhǎng),存儲(chǔ)的熱量將不斷累積,導(dǎo)致儲(chǔ)罐的容量配置需求不斷增大。為了使儲(chǔ)熱系統(tǒng)恢復(fù)到初始狀態(tài),必須設(shè)置空冷器排出一個(gè)大周期中積累的過(guò)剩熱量。

    4 結(jié) 論

    在模塊化化學(xué)品生產(chǎn)系統(tǒng)中,生產(chǎn)模塊的頻繁啟停將導(dǎo)致生產(chǎn)系統(tǒng)中流股的預(yù)熱、冷卻和反應(yīng)熱移除等均隨時(shí)間變化,換熱流股之間無(wú)法直接進(jìn)行熱回收,降低了生產(chǎn)系統(tǒng)中能量利用效率。為解決該問(wèn)題,本文提出了模塊化生產(chǎn)系統(tǒng)中分時(shí)儲(chǔ)熱設(shè)計(jì)儲(chǔ)熱系統(tǒng)的三步優(yōu)化設(shè)計(jì)方法。通過(guò)儲(chǔ)罐設(shè)置及流股匹配、儲(chǔ)罐儲(chǔ)熱溫區(qū)確定和儲(chǔ)罐容量配置及調(diào)度三個(gè)步驟,確定分時(shí)儲(chǔ)熱系統(tǒng)中儲(chǔ)罐的優(yōu)化設(shè)計(jì)和運(yùn)行方案,以解決波動(dòng)條件下的模塊化生產(chǎn)過(guò)程中能量的儲(chǔ)存和調(diào)度問(wèn)題,提升系統(tǒng)的能量利用效率。

    本文以可再生能源驅(qū)動(dòng)的甲醇模塊化生產(chǎn)系統(tǒng)為例,在給定生產(chǎn)工藝流程的工藝參數(shù)和生產(chǎn)線模塊的調(diào)度情況下采用分時(shí)儲(chǔ)熱策略優(yōu)化設(shè)計(jì)了儲(chǔ)熱系統(tǒng),并闡述和分析了所提出方法的實(shí)施過(guò)程。為了研究分時(shí)儲(chǔ)熱系統(tǒng)總費(fèi)用的影響因素,本文分為兩種情形分析和討論了儲(chǔ)熱熱量的廢棄對(duì)儲(chǔ)熱系統(tǒng)經(jīng)濟(jì)性的影響。研究表明,允許熱量廢棄使儲(chǔ)罐僅需滿足后續(xù)階段的換熱需求可大幅降低儲(chǔ)罐的容量配置,并提高設(shè)備的使用率;雖然熱量廢棄將額外增加換熱器設(shè)備費(fèi)用和操作費(fèi)用,但可使儲(chǔ)罐的費(fèi)用下降56.06%,總費(fèi)用下降

    40.38%。

    符 號(hào) 說(shuō) 明

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

    CP——流股熱容流率,MW/℃

    CRF——投資回收因子

    cp——儲(chǔ)熱介質(zhì)比熱容,J/(g·℃)

    F——流率,t/h

    LF——負(fù)荷因子

    NU——設(shè)備啟停數(shù)量

    OP——換熱器操作費(fèi)用,USD/(MWh)

    P——換熱器功率,MW

    Q——儲(chǔ)罐容量,t

    r——年利率

    S——?jiǎng)澐纸M別中的流股

    T——熱流股和儲(chǔ)熱介質(zhì)進(jìn)出口溫度,℃

    t——冷流股進(jìn)出口溫度,℃

    上角標(biāo)

    HE——換熱器

    in——輸入

    n——設(shè)備壽命

    out——輸出

    ST——儲(chǔ)熱儲(chǔ)罐

    waste——廢棄

    下角標(biāo)

    c——冷流股

    cap——投資費(fèi)用

    h——熱流股

    i——組別

    j——組中流股

    k——第k個(gè)時(shí)間段

    mai——維護(hù)費(fèi)用

    max——最大值

    min——最小值

    rated——額定值

    total——總費(fèi)用

    猜你喜歡
    流股儲(chǔ)熱冷器
    滑油空冷器壓降控制研究
    堿金屬熔鹽修飾MgO對(duì)其儲(chǔ)熱性能的影響研究
    硫化氫化學(xué)反應(yīng)循環(huán)制氫工藝的流程模擬
    空調(diào)表冷器的研究現(xiàn)狀與發(fā)展趨勢(shì)
    制冷(2019年2期)2019-12-09 08:10:42
    低溫甲醇洗吸收塔產(chǎn)出液再生過(guò)程模擬研究
    基于轉(zhuǎn)運(yùn)模型的功交換網(wǎng)絡(luò)綜合
    太陽(yáng)能熱發(fā)電儲(chǔ)熱系統(tǒng)綜述
    鋁合金空-空中冷器符合翹片管材料選擇與結(jié)構(gòu)匹配
    大容量熱化學(xué)吸附儲(chǔ)熱原理及性能分析
    儲(chǔ)熱技術(shù)基礎(chǔ)(I)——儲(chǔ)熱的基本原理及研究新動(dòng)向
    亚洲专区中文字幕在线| 亚洲精品色激情综合| 久久精品国产亚洲av涩爱 | 国产久久久一区二区三区| 一夜夜www| 99在线人妻在线中文字幕| 欧美又色又爽又黄视频| 国产伦在线观看视频一区| 亚洲avbb在线观看| 91久久精品国产一区二区成人| 亚洲美女视频黄频| 1000部很黄的大片| 看免费成人av毛片| 天堂影院成人在线观看| 久久久国产成人精品二区| 两个人视频免费观看高清| 日日撸夜夜添| 久久欧美精品欧美久久欧美| 亚洲色图av天堂| a级一级毛片免费在线观看| 欧美日韩国产亚洲二区| 91麻豆av在线| 国产淫片久久久久久久久| 超碰av人人做人人爽久久| 日韩欧美三级三区| 国产男靠女视频免费网站| 欧美区成人在线视频| 亚洲,欧美,日韩| 在线播放国产精品三级| 国产视频内射| 性色avwww在线观看| 久久精品综合一区二区三区| 麻豆久久精品国产亚洲av| 欧美+日韩+精品| 联通29元200g的流量卡| 蜜桃久久精品国产亚洲av| 一进一出抽搐动态| 18禁裸乳无遮挡免费网站照片| 女生性感内裤真人,穿戴方法视频| 成人国产麻豆网| 亚洲自拍偷在线| 在线观看av片永久免费下载| 亚洲av五月六月丁香网| 欧美日本视频| 欧美潮喷喷水| 国产 一区 欧美 日韩| 久久久国产成人精品二区| 人妻少妇偷人精品九色| 一本久久中文字幕| 天天躁日日操中文字幕| 综合色av麻豆| 欧美成人a在线观看| 日本欧美国产在线视频| 精品日产1卡2卡| 国产探花在线观看一区二区| 黄色一级大片看看| 悠悠久久av| 久久国产乱子免费精品| 日日啪夜夜撸| 精品国产三级普通话版| 美女大奶头视频| 亚洲真实伦在线观看| 午夜福利在线在线| 国产极品精品免费视频能看的| 欧美绝顶高潮抽搐喷水| 免费观看精品视频网站| 亚洲国产精品久久男人天堂| 日本在线视频免费播放| 亚洲欧美精品综合久久99| 免费看日本二区| bbb黄色大片| 美女cb高潮喷水在线观看| 一个人免费在线观看电影| 精品久久久久久久久av| 他把我摸到了高潮在线观看| 91av网一区二区| 精品人妻偷拍中文字幕| 日韩中文字幕欧美一区二区| 国产精品乱码一区二三区的特点| 成人美女网站在线观看视频| 成人二区视频| 一级黄片播放器| 午夜精品在线福利| 中亚洲国语对白在线视频| 制服丝袜大香蕉在线| 琪琪午夜伦伦电影理论片6080| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕高清在线视频| 黄色丝袜av网址大全| 成人鲁丝片一二三区免费| 国产精品嫩草影院av在线观看 | 国产成人aa在线观看| 日本一二三区视频观看| 成人二区视频| av黄色大香蕉| 久久人妻av系列| 国产在视频线在精品| 又爽又黄无遮挡网站| 国产久久久一区二区三区| 国产大屁股一区二区在线视频| 中文字幕精品亚洲无线码一区| 日本一二三区视频观看| 久久午夜亚洲精品久久| 日韩 亚洲 欧美在线| 嫁个100分男人电影在线观看| 夜夜爽天天搞| 舔av片在线| 免费电影在线观看免费观看| av中文乱码字幕在线| 午夜免费成人在线视频| 亚洲真实伦在线观看| 久久久久久九九精品二区国产| 免费观看精品视频网站| 一本精品99久久精品77| 国产亚洲欧美98| 99在线视频只有这里精品首页| 日日干狠狠操夜夜爽| 狠狠狠狠99中文字幕| 国产伦人伦偷精品视频| 九九热线精品视视频播放| 久久久久久久久久黄片| 91麻豆精品激情在线观看国产| 最好的美女福利视频网| 亚洲成人中文字幕在线播放| 久久香蕉精品热| 亚洲国产精品合色在线| 国产午夜福利久久久久久| 伊人久久精品亚洲午夜| 欧美日本视频| 婷婷亚洲欧美| 又黄又爽又免费观看的视频| 中国美白少妇内射xxxbb| 色综合婷婷激情| 韩国av一区二区三区四区| 中国美白少妇内射xxxbb| 日韩国内少妇激情av| 国产精品99久久久久久久久| 久久久久免费精品人妻一区二区| 91在线观看av| 欧美日韩国产亚洲二区| 免费看光身美女| 日韩国内少妇激情av| 一区二区三区激情视频| 少妇熟女aⅴ在线视频| 久久久久免费精品人妻一区二区| 久久久精品大字幕| 国产在视频线在精品| 成人美女网站在线观看视频| 亚洲av中文字字幕乱码综合| 亚洲午夜理论影院| 搡老熟女国产l中国老女人| av专区在线播放| 又爽又黄a免费视频| 别揉我奶头 嗯啊视频| 精品久久久噜噜| 欧美日本视频| 日韩,欧美,国产一区二区三区 | 人妻丰满熟妇av一区二区三区| 亚洲午夜理论影院| 日本一本二区三区精品| 欧美色欧美亚洲另类二区| 一级黄片播放器| 欧美色欧美亚洲另类二区| 亚洲人与动物交配视频| 麻豆精品久久久久久蜜桃| 欧美3d第一页| 亚洲欧美精品综合久久99| 蜜桃久久精品国产亚洲av| 亚洲图色成人| 国产综合懂色| 可以在线观看毛片的网站| 国产精品人妻久久久久久| 中国美白少妇内射xxxbb| 内射极品少妇av片p| 欧美精品啪啪一区二区三区| 三级男女做爰猛烈吃奶摸视频| 日韩中文字幕欧美一区二区| 成人av一区二区三区在线看| 欧美日韩黄片免| 成人毛片a级毛片在线播放| 又粗又爽又猛毛片免费看| 2021天堂中文幕一二区在线观| 亚洲精品日韩av片在线观看| 国产不卡一卡二| 国产精品久久久久久精品电影| 日日撸夜夜添| 久久精品影院6| 成人特级av手机在线观看| 噜噜噜噜噜久久久久久91| 超碰av人人做人人爽久久| 成人午夜高清在线视频| 嫩草影院入口| 美女cb高潮喷水在线观看| 成熟少妇高潮喷水视频| 成人欧美大片| 国产精品电影一区二区三区| 亚洲国产日韩欧美精品在线观看| 伊人久久精品亚洲午夜| 小说图片视频综合网站| 一进一出抽搐gif免费好疼| 人人妻人人澡欧美一区二区| 国产一区二区在线观看日韩| 99久久精品热视频| 国产男人的电影天堂91| 国产一区二区亚洲精品在线观看| 18禁在线播放成人免费| 99久久无色码亚洲精品果冻| 国产亚洲av嫩草精品影院| 午夜激情福利司机影院| 极品教师在线视频| 午夜久久久久精精品| 成人二区视频| 午夜爱爱视频在线播放| 久久精品国产99精品国产亚洲性色| 欧美日韩亚洲国产一区二区在线观看| 亚洲图色成人| 赤兔流量卡办理| 欧美不卡视频在线免费观看| 日本-黄色视频高清免费观看| 欧美一区二区精品小视频在线| 久久人人精品亚洲av| 久久亚洲精品不卡| 国产精品日韩av在线免费观看| 亚洲av第一区精品v没综合| 中国美女看黄片| 久99久视频精品免费| 久久久国产成人精品二区| 黄色丝袜av网址大全| 99热精品在线国产| 日韩国内少妇激情av| 国产精品98久久久久久宅男小说| 永久网站在线| 日本五十路高清| 精品国产三级普通话版| 天堂√8在线中文| 国产精品女同一区二区软件 | 精品人妻一区二区三区麻豆 | 日韩亚洲欧美综合| 久久久国产成人精品二区| 色尼玛亚洲综合影院| 看片在线看免费视频| 熟妇人妻久久中文字幕3abv| 国产亚洲精品久久久久久毛片| 午夜福利18| 五月玫瑰六月丁香| 舔av片在线| 国产精品嫩草影院av在线观看 | 国产色爽女视频免费观看| 日韩欧美在线乱码| 免费观看在线日韩| 国产三级在线视频| 亚洲av.av天堂| 欧美一区二区亚洲| 国产淫片久久久久久久久| 久久草成人影院| 色综合婷婷激情| 变态另类成人亚洲欧美熟女| 欧美一区二区精品小视频在线| 少妇被粗大猛烈的视频| 免费高清视频大片| 免费观看人在逋| 麻豆成人午夜福利视频| 欧美另类亚洲清纯唯美| 免费在线观看日本一区| 黄色日韩在线| 国产午夜精品久久久久久一区二区三区 | 国产黄色小视频在线观看| a在线观看视频网站| 亚洲人成网站在线播| 国产精品野战在线观看| 99久久九九国产精品国产免费| 色综合色国产| 91久久精品电影网| av天堂在线播放| 欧美黑人欧美精品刺激| 女人十人毛片免费观看3o分钟| 一进一出抽搐动态| 精品久久久噜噜| 黄色一级大片看看| 亚洲 国产 在线| 久久这里只有精品中国| 丰满乱子伦码专区| 嫁个100分男人电影在线观看| 最新在线观看一区二区三区| 乱码一卡2卡4卡精品| 不卡视频在线观看欧美| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 校园人妻丝袜中文字幕| 国产淫片久久久久久久久| 国产精品一区二区免费欧美| av福利片在线观看| 国产真实伦视频高清在线观看 | 国产久久久一区二区三区| 亚洲欧美日韩东京热| 人人妻,人人澡人人爽秒播| 日韩欧美三级三区| 天天躁日日操中文字幕| 成人二区视频| 搡老熟女国产l中国老女人| 亚洲国产高清在线一区二区三| 色综合亚洲欧美另类图片| 日日干狠狠操夜夜爽| 亚洲人成伊人成综合网2020| 国产久久久一区二区三区| 在线播放无遮挡| 91午夜精品亚洲一区二区三区 | 亚洲一区高清亚洲精品| 午夜免费男女啪啪视频观看 | 看免费成人av毛片| 国产女主播在线喷水免费视频网站 | 国产一区二区激情短视频| 精品久久久久久久末码| 亚洲久久久久久中文字幕| 看黄色毛片网站| 午夜福利在线观看免费完整高清在 | 国产精品爽爽va在线观看网站| 亚洲最大成人手机在线| 老女人水多毛片| 色哟哟哟哟哟哟| 色哟哟·www| 嫩草影院新地址| 国产美女午夜福利| 成人精品一区二区免费| 婷婷精品国产亚洲av| 日韩欧美国产一区二区入口| 免费不卡的大黄色大毛片视频在线观看 | 国产成人影院久久av| 精品99又大又爽又粗少妇毛片 | 久久人人爽人人爽人人片va| 色精品久久人妻99蜜桃| 国产69精品久久久久777片| 在线免费观看不下载黄p国产 | 亚洲中文字幕一区二区三区有码在线看| 日本-黄色视频高清免费观看| 婷婷六月久久综合丁香| 有码 亚洲区| av女优亚洲男人天堂| 淫妇啪啪啪对白视频| 精品一区二区三区人妻视频| 88av欧美| 日韩高清综合在线| 国产精品久久电影中文字幕| 老司机深夜福利视频在线观看| 人人妻人人澡欧美一区二区| a级毛片免费高清观看在线播放| 精品无人区乱码1区二区| 久久久久久久亚洲中文字幕| 夜夜看夜夜爽夜夜摸| 免费电影在线观看免费观看| 久久久久国内视频| 中文字幕人妻熟人妻熟丝袜美| 国产老妇女一区| 亚洲av.av天堂| 性色avwww在线观看| 精品久久久久久久久久久久久| 精品福利观看| av国产免费在线观看| 综合色av麻豆| 久99久视频精品免费| a级一级毛片免费在线观看| 精品人妻偷拍中文字幕| 男女边吃奶边做爰视频| 日韩 亚洲 欧美在线| 波多野结衣巨乳人妻| 欧美+亚洲+日韩+国产| 国产一区二区亚洲精品在线观看| 乱系列少妇在线播放| 日日摸夜夜添夜夜添小说| ponron亚洲| 亚洲无线在线观看| 国内精品久久久久久久电影| 日韩在线高清观看一区二区三区 | 天堂√8在线中文| 国产真实乱freesex| 在线国产一区二区在线| 人人妻,人人澡人人爽秒播| 亚洲av一区综合| 久久草成人影院| 亚洲经典国产精华液单| 97人妻精品一区二区三区麻豆| 一区二区三区四区激情视频 | 亚洲第一区二区三区不卡| АⅤ资源中文在线天堂| 国产麻豆成人av免费视频| 赤兔流量卡办理| 亚洲一区二区三区色噜噜| 十八禁网站免费在线| 国产淫片久久久久久久久| 在线播放无遮挡| 一a级毛片在线观看| 欧美xxxx黑人xx丫x性爽| 级片在线观看| 色哟哟哟哟哟哟| 久久久精品欧美日韩精品| 国内精品久久久久久久电影| 国产伦一二天堂av在线观看| avwww免费| 亚洲最大成人手机在线| 午夜视频国产福利| 国产 一区精品| 久久午夜福利片| 人妻丰满熟妇av一区二区三区| 国产爱豆传媒在线观看| 欧美一区二区亚洲| 中文字幕免费在线视频6| 在线观看66精品国产| 美女黄网站色视频| 国产精品一区www在线观看 | 久久精品国产清高在天天线| 精品一区二区三区视频在线| 18+在线观看网站| 一级a爱片免费观看的视频| 欧美成人a在线观看| 精品久久久久久久人妻蜜臀av| 精品一区二区免费观看| 免费在线观看成人毛片| 亚洲精华国产精华液的使用体验 | 国产日本99.免费观看| 欧美日韩瑟瑟在线播放| 亚州av有码| 免费看av在线观看网站| 国产熟女欧美一区二区| 噜噜噜噜噜久久久久久91| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区激情视频| 国内毛片毛片毛片毛片毛片| 亚洲在线观看片| 精品人妻视频免费看| 18禁黄网站禁片午夜丰满| 熟女电影av网| 少妇裸体淫交视频免费看高清| 国产 一区 欧美 日韩| 97碰自拍视频| videossex国产| 波多野结衣巨乳人妻| 日韩中字成人| 成年人黄色毛片网站| 18+在线观看网站| 久久精品国产亚洲av涩爱 | 热99re8久久精品国产| 国产精品乱码一区二三区的特点| 久久久久久伊人网av| 狂野欧美白嫩少妇大欣赏| 少妇熟女aⅴ在线视频| 免费av毛片视频| 亚洲成人久久爱视频| 一个人看视频在线观看www免费| 一进一出好大好爽视频| 乱系列少妇在线播放| 欧美日韩国产亚洲二区| 他把我摸到了高潮在线观看| 在现免费观看毛片| 真人做人爱边吃奶动态| 午夜福利欧美成人| 久久久国产成人免费| 天天一区二区日本电影三级| 亚洲一级一片aⅴ在线观看| 国内精品久久久久久久电影| 白带黄色成豆腐渣| АⅤ资源中文在线天堂| 国产v大片淫在线免费观看| 动漫黄色视频在线观看| 我要搜黄色片| 亚洲在线观看片| 18禁黄网站禁片免费观看直播| 国产精华一区二区三区| 尾随美女入室| 国产精品免费一区二区三区在线| 一级黄片播放器| 国产综合懂色| 免费观看在线日韩| 在线观看美女被高潮喷水网站| 中文字幕人妻熟人妻熟丝袜美| 亚洲av中文av极速乱 | 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 成人av在线播放网站| 欧美不卡视频在线免费观看| 亚洲内射少妇av| 国产精品久久电影中文字幕| 日韩欧美 国产精品| 99国产精品一区二区蜜桃av| 免费人成在线观看视频色| 欧美bdsm另类| 18禁黄网站禁片午夜丰满| 欧美不卡视频在线免费观看| 少妇的逼好多水| 九九爱精品视频在线观看| 久久精品人妻少妇| 日韩高清综合在线| 熟妇人妻久久中文字幕3abv| 欧美精品国产亚洲| 亚洲不卡免费看| 欧美又色又爽又黄视频| 国产成人影院久久av| 国产免费男女视频| 亚洲熟妇熟女久久| 99久久成人亚洲精品观看| av国产免费在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲经典国产精华液单| 亚洲狠狠婷婷综合久久图片| 国产高清三级在线| 精华霜和精华液先用哪个| 日本撒尿小便嘘嘘汇集6| 亚洲一区二区三区色噜噜| 性插视频无遮挡在线免费观看| 欧美日韩国产亚洲二区| АⅤ资源中文在线天堂| 亚洲人成网站在线播| 免费看a级黄色片| 高清毛片免费观看视频网站| 欧美日韩国产亚洲二区| 波多野结衣高清作品| 国产单亲对白刺激| 国产精品综合久久久久久久免费| 91久久精品国产一区二区成人| 国产精品一区二区免费欧美| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 搞女人的毛片| 女生性感内裤真人,穿戴方法视频| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久久久亚洲 | 中文亚洲av片在线观看爽| 女人十人毛片免费观看3o分钟| 看黄色毛片网站| 亚洲精品色激情综合| 欧美国产日韩亚洲一区| 亚洲欧美日韩无卡精品| 能在线免费观看的黄片| 国产熟女欧美一区二区| 夜夜爽天天搞| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产高清三级在线| 黄色女人牲交| 亚洲人成网站在线播| 午夜影院日韩av| 三级毛片av免费| 悠悠久久av| 成年人黄色毛片网站| 三级男女做爰猛烈吃奶摸视频| 一区二区三区四区激情视频 | 精品久久久久久久末码| .国产精品久久| 亚洲va在线va天堂va国产| av福利片在线观看| 高清毛片免费观看视频网站| 日本黄色视频三级网站网址| 日韩高清综合在线| 久久香蕉精品热| 精品久久久久久久末码| 亚洲av二区三区四区| 欧美激情在线99| 国内精品一区二区在线观看| 男女之事视频高清在线观看| 亚洲性久久影院| 99九九线精品视频在线观看视频| 尾随美女入室| 偷拍熟女少妇极品色| 老司机深夜福利视频在线观看| 91在线精品国自产拍蜜月| 网址你懂的国产日韩在线| 成人午夜高清在线视频| 少妇熟女aⅴ在线视频| av在线天堂中文字幕| 中国美白少妇内射xxxbb| 久久久久久大精品| 在线播放无遮挡| 99九九线精品视频在线观看视频| 看黄色毛片网站| 免费大片18禁| 国产成人一区二区在线| 我要看日韩黄色一级片| 免费一级毛片在线播放高清视频| 色哟哟·www| 午夜日韩欧美国产| or卡值多少钱| 国产伦一二天堂av在线观看| 亚洲七黄色美女视频| 黄色一级大片看看| 男女视频在线观看网站免费| 国产在线精品亚洲第一网站| АⅤ资源中文在线天堂| 人妻丰满熟妇av一区二区三区| 国产成人福利小说| 亚洲中文日韩欧美视频| 亚洲最大成人av| 日日摸夜夜添夜夜添av毛片 | 精品久久久久久成人av| 久久这里只有精品中国| 级片在线观看| 97超视频在线观看视频| 国产 一区 欧美 日韩| 精品乱码久久久久久99久播| 国产在线精品亚洲第一网站| 麻豆国产97在线/欧美| 午夜日韩欧美国产| 在线观看美女被高潮喷水网站| 成人高潮视频无遮挡免费网站| 一个人观看的视频www高清免费观看| av福利片在线观看| 亚洲最大成人手机在线| 国产高潮美女av| 日本黄色片子视频| 欧美bdsm另类| 久久精品国产鲁丝片午夜精品 | 精品日产1卡2卡| 国产精品av视频在线免费观看| 久久久久久久久久久丰满 | 久久久久久久精品吃奶| 麻豆成人av在线观看| 国产欧美日韩一区二区精品|