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

    基于模擬的單/雙壓蒸發(fā)有機(jī)朗肯循環(huán)結(jié)構(gòu)選擇

    2021-11-12 07:21:26唐建平康麗霞劉永忠
    關(guān)鍵詞:夾點(diǎn)臨界溫度工質(zhì)

    唐建平, 康麗霞, 2, 劉永忠, 2

    基于模擬的單/雙壓蒸發(fā)有機(jī)朗肯循環(huán)結(jié)構(gòu)選擇

    唐建平1, 康麗霞1, 2, 劉永忠1, 2

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

    針對(duì)單/雙壓蒸發(fā)有機(jī)朗肯循環(huán)(ORC)系統(tǒng)的工質(zhì)和結(jié)構(gòu)選擇問題,結(jié)合Aspen Plus模擬和夾點(diǎn)理論分析,以120~200 ℃的低溫?zé)嵩礊槔?,?duì)單雙壓蒸發(fā)ORC系統(tǒng)在不同工質(zhì)條件下的熱力學(xué)性能和換熱過程進(jìn)行了研究。結(jié)果表明:熱源和工質(zhì)換熱形成的夾點(diǎn)類型是單/雙壓蒸發(fā)ORC系統(tǒng)結(jié)構(gòu)選擇的關(guān)鍵。對(duì)于單壓蒸發(fā)ORC,隨著工質(zhì)臨界溫度的降低,工質(zhì)在最優(yōu)條件下與熱源形成的夾點(diǎn)類型由泡點(diǎn)夾點(diǎn)(VPP)逐漸變?yōu)榕蔹c(diǎn)和進(jìn)口雙夾點(diǎn)(VPP&PPP),并最終形成單獨(dú)的進(jìn)口夾點(diǎn)(PPP);在形成VPP&PPP雙夾點(diǎn)時(shí),單壓蒸發(fā)系統(tǒng)效率達(dá)到最大;研究還指出VPP&PPP雙夾點(diǎn)是單/雙蒸發(fā)結(jié)構(gòu)選擇的臨界點(diǎn),并據(jù)此獲得了單/雙壓蒸發(fā)系統(tǒng)結(jié)構(gòu)選擇的定量關(guān)聯(lián)式,為ORC系統(tǒng)工質(zhì)和結(jié)構(gòu)的選擇提供了定量分析基礎(chǔ)。

    低溫余熱;有機(jī)朗肯循環(huán);雙壓蒸發(fā)結(jié)構(gòu);過程模擬;夾點(diǎn)

    1 前 言

    有機(jī)朗肯循環(huán)(organic Rankine cycle,ORC)已被證實(shí)是實(shí)際回收利用低溫余熱最有效的方式之一[1]。近年來,為了最大限度地提升ORC系統(tǒng)性能,各種不同類型的ORC結(jié)構(gòu)相繼被提出和應(yīng)用,如同流回?zé)峤Y(jié)構(gòu)[2]、側(cè)線回?zé)峤Y(jié)構(gòu)[3]、雙(多)壓蒸發(fā)結(jié)構(gòu)[4]、有機(jī)閃蒸結(jié)構(gòu)[5]、超臨界循環(huán)結(jié)構(gòu)[6]等,極大地推動(dòng)了ORC在余熱回收領(lǐng)域的進(jìn)一步發(fā)展和應(yīng)用[7]。針對(duì)ORC系統(tǒng)結(jié)構(gòu)的選擇和優(yōu)化也成為研究者們關(guān)注的重點(diǎn)之一。

    基礎(chǔ)ORC結(jié)構(gòu)主要包含工質(zhì)泵、蒸發(fā)器、透平機(jī)和冷凝器4個(gè)單元模塊。雙壓蒸發(fā)循環(huán)結(jié)構(gòu)是在此基礎(chǔ)上,將整個(gè)熱回收過程分成2個(gè)壓力級(jí),分別產(chǎn)生不同品質(zhì)的蒸氣,不但可以提升工質(zhì)與余熱流股的換熱效率,提高余熱回收量,還可進(jìn)一步減小系統(tǒng)的不可逆損失[8],是目前ORC系統(tǒng)中最具潛力的結(jié)構(gòu)之一。Du等[4]構(gòu)建了雙壓蒸發(fā)ORC循環(huán)中各單元模塊的數(shù)學(xué)模型,并以系統(tǒng)的發(fā)電輸出為目標(biāo)對(duì)雙壓力系統(tǒng)進(jìn)行分析。結(jié)果表明,雙壓蒸發(fā)ORC系統(tǒng)的輸出功比基礎(chǔ)結(jié)構(gòu)高18.9%,傳熱面積增加42.14%,且隨著熱源流量的上升,高壓蒸發(fā)級(jí)與低壓蒸發(fā)級(jí)的輸出功比值降低。Li等[9]以總輸出和總導(dǎo)熱之比為目標(biāo),建立了雙壓蒸發(fā)ORC系統(tǒng)的火積理論模型,對(duì)比和分析了串聯(lián)和并聯(lián)2種不同的雙壓力蒸發(fā)循環(huán)結(jié)構(gòu)。結(jié)果表明,串聯(lián)結(jié)構(gòu)相對(duì)于并聯(lián)結(jié)構(gòu)更優(yōu),也更適合工業(yè)推廣,且串聯(lián)的雙壓蒸發(fā)ORC結(jié)構(gòu)相較于單壓蒸發(fā)結(jié)構(gòu)有更好的熱力學(xué)性能,但單壓蒸發(fā)ORC結(jié)構(gòu)的經(jīng)濟(jì)性優(yōu)于雙壓蒸發(fā)ORC結(jié)構(gòu)[10]。Li等[11]通過對(duì)小規(guī)模串聯(lián)雙壓蒸發(fā)結(jié)構(gòu)的熱力學(xué)和經(jīng)濟(jì)性分析指出,熱源溫度的升高和流量的增加可有效降低雙壓蒸發(fā)結(jié)構(gòu)的總費(fèi)用,使特定工況下雙壓蒸發(fā)結(jié)構(gòu)的經(jīng)濟(jì)性優(yōu)于單壓蒸發(fā)結(jié)構(gòu)。

    Manente等[12]在熱源進(jìn)口溫度為100~200 ℃時(shí),對(duì)單壓和串聯(lián)雙壓蒸發(fā)結(jié)構(gòu)熱力學(xué)性能進(jìn)行了對(duì)比研究。結(jié)果表明,熱源進(jìn)口溫度在100~125 ℃時(shí),雙壓蒸發(fā)ORC結(jié)構(gòu)下系統(tǒng)輸出最高可提升29%,但在熱源進(jìn)口溫度位于150~200 ℃時(shí),雙壓蒸發(fā)ORC系統(tǒng)的優(yōu)勢(shì)減弱。特別地,當(dāng)工質(zhì)的臨界溫度低于熱源進(jìn)口溫度40 ℃以上時(shí),雙壓蒸發(fā)結(jié)構(gòu)將不再適用。在此基礎(chǔ)上,Li等[13]提出了熱源溫度與工質(zhì)臨界溫度的定量表達(dá)式來確定單/雙壓蒸發(fā)結(jié)構(gòu)的適用性。除此之外,他們還將串聯(lián)雙壓蒸發(fā)ORC結(jié)構(gòu)拓展到高壓蒸發(fā)為超臨界狀態(tài),低壓蒸發(fā)為亞臨界狀態(tài)的結(jié)構(gòu)[14]。該結(jié)構(gòu)的發(fā)電輸出較單一的超臨界ORC循環(huán)和基礎(chǔ)ORC循環(huán)提升19.9% 和48.9%。Wang等[15]以系統(tǒng)輸出為目標(biāo),構(gòu)建了雙壓蒸發(fā)結(jié)構(gòu)和超臨界結(jié)構(gòu)下ORC系統(tǒng)的數(shù)學(xué)模型。計(jì)算結(jié)果表明,高壓力級(jí)對(duì)整個(gè)ORC系統(tǒng)的性能影響更大,且隨著高壓力級(jí)蒸發(fā)壓力的不斷提升,雙壓力結(jié)構(gòu)下熱回收部分的熱量交換更加充分。此外,還定義了邊界溫度用于判斷雙壓力結(jié)構(gòu)和超臨界結(jié)構(gòu)的適用性,高于邊界溫度,超臨界結(jié)構(gòu)更適用,反之,則雙壓力結(jié)構(gòu)更適用。

    綜上所述,目前針對(duì)雙壓蒸發(fā)ORC循環(huán)系統(tǒng)的研究多關(guān)注系統(tǒng)的熱力學(xué)性能或經(jīng)濟(jì)性能分析,缺少統(tǒng)一的結(jié)構(gòu)選擇依據(jù)。盡管部分研究者[12-13]提出了定性規(guī)則或定量的關(guān)系式,用于判斷單/雙壓蒸發(fā)ORC結(jié)構(gòu)在不同熱源溫度下的適用性,但研究方法以熱力學(xué)建模分析為主,缺少對(duì)系統(tǒng)內(nèi)部機(jī)理的深入探究和分析,加之研究中選擇的候選流體有限,使得確定的規(guī)則準(zhǔn)確性和普適性不足。為此,本研究針對(duì)實(shí)際工業(yè)生產(chǎn)中最常見的低溫?zé)嵩礈囟确秶?,挑選了盡可能多的候選流體,通過對(duì)單雙壓蒸發(fā)ORC結(jié)構(gòu)在不同熱源進(jìn)口溫度下適用性的模擬研究和分析,闡明了單/雙壓蒸發(fā)結(jié)構(gòu)選擇的內(nèi)部機(jī)理和系統(tǒng)結(jié)構(gòu)對(duì)系統(tǒng)性能的影響,并總結(jié)了單/雙壓蒸發(fā)ORC結(jié)構(gòu)的選擇原則。

    2 單/雙壓蒸發(fā)ORC系統(tǒng)模擬

    2.1 單/雙壓蒸發(fā)ORC系統(tǒng)結(jié)構(gòu)

    本研究中所探究的單壓蒸發(fā)循環(huán)(single pressure evaporation cycle,SPEC)和雙壓蒸發(fā)循環(huán)(dual pressure evaporation cycle,DPEC)結(jié)構(gòu)及其對(duì)應(yīng)的-(溫-熵)圖如圖1和2所示,圖中為溫度。在SPEC中,工質(zhì)泵(1→2)將冷凝器出口的液體加壓并輸送到熱回收裝置,熱回收裝置包括預(yù)熱器(2→3)和蒸發(fā)器(3→4) 2部分,在熱回收裝置中,工質(zhì)回收余熱中的熱量,由過冷液體變成高溫高壓蒸氣狀態(tài)。熱回收裝置出口的高溫高壓蒸氣進(jìn)入透平機(jī)(4→5)膨脹后變成低溫低壓乏氣,同時(shí)產(chǎn)生發(fā)電輸出,透平機(jī)出口的乏氣進(jìn)入冷凝器(4→5→1),被公用工程冷凝成飽和液體,如圖1所示。

    圖1 單壓蒸發(fā)循環(huán)結(jié)構(gòu)圖和t-S圖

    研究采用的雙壓蒸發(fā)循環(huán)結(jié)構(gòu)為串聯(lián)的雙壓力蒸發(fā)結(jié)構(gòu),該結(jié)構(gòu)在單壓蒸發(fā)ORC結(jié)構(gòu)的基礎(chǔ)上增加了一個(gè)壓力級(jí)(3→5→6→7→8),將熱回收裝置分為高壓(high pressure,HP)和低壓(low pressure,LP) 2個(gè)部分,產(chǎn)生2種不同品質(zhì)的工質(zhì)蒸氣,并分別進(jìn)入高壓透平機(jī)(7→8)和低壓透平機(jī)(9→10),高壓透平機(jī)出口的蒸氣壓力與低壓熱回收裝置的蒸發(fā)壓力相等,低壓部分預(yù)熱器出口的飽和液體分流一部分通過高壓工質(zhì)泵進(jìn)入高壓熱回收部分。需要注意的是,該串聯(lián)蒸發(fā)結(jié)構(gòu)較并聯(lián)蒸發(fā)結(jié)構(gòu)具有更好的熱回收性能[9],且其中透平機(jī)的布置方式也選擇了熱力學(xué)性能和經(jīng)濟(jì)性更優(yōu)的交互透平機(jī)布置方式[16],如圖2所示。

    圖2 雙壓蒸發(fā)循環(huán)結(jié)構(gòu)和t-S圖

    2.2 候選工質(zhì)

    工質(zhì)的選擇直接影響ORC系統(tǒng)的熱回收性能。大量研究表明干、等熵工質(zhì)具有更好的熱回收性能,而且能夠有效避免透平機(jī)出口的液相流體對(duì)透平機(jī)造成損害[18]。在挑選工質(zhì)時(shí),有機(jī)工質(zhì)的臨界溫度與熱源溫度也存在一定的匹配關(guān)系,即低臨界溫度的工質(zhì)適用于溫度較低的熱源,而高溫?zé)嵩锤m合選擇臨界溫度較高的工質(zhì)[19],最優(yōu)的工質(zhì)的臨界溫度低于熱源的進(jìn)口溫度約30~50 ℃,臨界溫度相近的工質(zhì),其飽和蒸氣曲線的斜率越小,也就是臨界壓力越高的工質(zhì)有更高的輸出[20],當(dāng)工質(zhì)的臨界溫度低于熱源進(jìn)口溫度55 ℃以上時(shí),熱回收效率會(huì)變得很低[21]。

    表1 候選工質(zhì)的物性

    在工業(yè)領(lǐng)域,一般的低溫余熱為生產(chǎn)過程產(chǎn)生的低于200 ℃的余熱氣、冷凝水、熱水等。為此,依據(jù)以上的工質(zhì)選擇原則,研究選取了如表1所示的17種干、等熵工質(zhì)作為候選工質(zhì),表1還給出了候選工質(zhì)臨界溫度t和臨界壓力c。17種工質(zhì)的臨界溫度為204.35~85.8 ℃,可適用于絕大多數(shù)200 ℃以下的低溫?zé)嵩础?/p>

    2.3 ORC系統(tǒng)模擬

    2.3.1 基本假設(shè)

    鑒于其詳盡的物性數(shù)據(jù)庫(kù)、精確的物性計(jì)算方法、成熟完備的單元模型以及廣泛的工業(yè)應(yīng)用,研究采用Aspen軟件對(duì)單/雙壓力蒸發(fā)ORC系統(tǒng)進(jìn)行模擬。根據(jù)表1中候選流體的性質(zhì),選用Peng-Robinson(P-R)方程作為物性方程[22],因?yàn)镻-R方程在臨界區(qū)域內(nèi)有較高的準(zhǔn)確性,同時(shí)也能對(duì)物質(zhì)的液相密度等物性能進(jìn)行較好的預(yù)測(cè)[23],在有機(jī)流體物性模擬預(yù)測(cè)方面有廣泛的應(yīng)用[24-25]。

    模擬中用到的假設(shè)條件歸納如下:

    1) 選取入口溫度為120~200 ℃水作為熱源流股,壓力為20′105Pa;模擬中以10 ℃為步長(zhǎng)提高熱源入口溫度;熱源流量為1 kg×s-1,不限制熱源的出口溫度。

    2) 除了透平機(jī)和工質(zhì)泵的壓力變化外,只考慮蒸發(fā)器和冷凝器的壓降,且設(shè)定為10 kPa[26],忽略系統(tǒng)中其他附屬設(shè)備(如管道等)的壓降。

    3) 考慮到干、等熵工質(zhì)的過熱對(duì)系統(tǒng)的熱回收性能影響較小[3],忽略工質(zhì)的過熱,假設(shè)蒸發(fā)器出口為飽和蒸氣。

    4) 冷凝器出口的工質(zhì)狀態(tài)為40 ℃的飽和液體,且為了適應(yīng)實(shí)際工業(yè)場(chǎng)景,整個(gè)裝置的壓力應(yīng)大于等于大氣壓。

    5) 系統(tǒng)的最小傳熱溫差為5 ℃,透平機(jī)和工質(zhì)泵的等熵效率為常數(shù),取值為0.8,透平機(jī)出口氣相分率為1。

    2.3.2 評(píng)價(jià)指標(biāo)

    1) 系統(tǒng)效率

    ORC系統(tǒng)效率,sys定義為系統(tǒng)的總輸出(tur-pum)與可回收余熱總量av的比值,如式(1)所示:

    因?yàn)楦哂诠べ|(zhì)冷凝溫度5 ℃(最小傳熱溫差)的熱源皆具有被回收利用的潛力,因此文中高于45 ℃的余熱均是可回收余熱。特別地,對(duì)于雙壓蒸發(fā)結(jié)構(gòu)而言,系統(tǒng)的凈輸出功應(yīng)是所有透平機(jī)輸出功和所有工質(zhì)泵消耗功率的差值。

    2) 熱回收效率

    系統(tǒng)的熱回收效率,thm定義為系統(tǒng)的總輸出(tur-pum)與實(shí)際回收余熱re的比值,如式(2)所示:

    其中,系統(tǒng)實(shí)際回收的余熱量可表示為蒸發(fā)器負(fù)荷和預(yù)熱器負(fù)荷的總和。

    3) 其他約束

    本文將蒸發(fā)器出口溫度定義為決策變量,且蒸發(fā)溫度應(yīng)介于工藝限制的最高溫度和冷凝溫度之間,如式(3)所示。

    對(duì)雙壓蒸發(fā)循環(huán)而言,低壓蒸發(fā)級(jí)的蒸發(fā)溫度應(yīng)介于高壓力級(jí)蒸發(fā)溫度和冷凝溫度之間,如式(4)所示。

    式(3)和(4)中,con和eva分別表示冷凝和蒸發(fā)溫度;WH,in為熱源進(jìn)口溫度;c為工質(zhì)臨界溫度;Δpp和Δd分別為最小傳熱溫差、工質(zhì)的最高蒸發(fā)溫度與其臨界溫度的最小差值;上標(biāo)LP和HP分別為低壓和高壓蒸發(fā)過程。需要注意的是,為了獲得最優(yōu)的操作條件,在模擬過程中,將在以上給定的變化范圍內(nèi),以2 ℃為步長(zhǎng)對(duì)蒸發(fā)溫度進(jìn)行離散。

    3 單壓蒸發(fā)ORC系統(tǒng)分析

    3.1 固定熱源的單壓蒸發(fā)ORC系統(tǒng)分析

    在固定熱源入口溫度為200 ℃時(shí),對(duì)單壓蒸發(fā)的基礎(chǔ)ORC結(jié)構(gòu)進(jìn)行模擬,得到不同候選工質(zhì)下對(duì)應(yīng)的最優(yōu)蒸發(fā)溫度和最大的系統(tǒng)效率,如圖3所示,圖中橫坐標(biāo)的工質(zhì)依照其臨界溫度由高到低排列。

    由圖3可以看出當(dāng)熱源進(jìn)口溫度為200 ℃時(shí),系統(tǒng)效率最大時(shí)對(duì)應(yīng)的工質(zhì)為R245ca,最大的系統(tǒng)效率為0.153,當(dāng)熱源進(jìn)口溫度為200 ℃時(shí),系統(tǒng)對(duì)應(yīng)的可用熱源的熱量為730 kW,系統(tǒng)最大的輸出為111.89 kW。同時(shí)從圖3可以看出,R245ca的最優(yōu)蒸發(fā)溫度也是最高的,為169 ℃,意味著使用工質(zhì)R245ca可以產(chǎn)生最高品質(zhì)的蒸氣。由圖3還可看出,ORC的系統(tǒng)效率隨工質(zhì)臨界溫度降低呈現(xiàn)先增大后減小的趨勢(shì),在工質(zhì)為R245ca時(shí)出現(xiàn)一個(gè)極大值,對(duì)應(yīng)最大的系統(tǒng)效率為0.153。除R245ca外,R245fa2,R603,R245fa等工質(zhì)下對(duì)應(yīng)的系統(tǒng)效率也較高,其臨界溫度與熱源進(jìn)口溫度的溫差在26~45 ℃。溫差超出該范圍時(shí)系統(tǒng)效率都較低。特別地,工質(zhì)R603的臨界溫度R245fa高約6 ℃,但兩者對(duì)應(yīng)的系統(tǒng)效率和發(fā)電輸出卻基本相同,這是因?yàn)镽603的臨界壓力高于R245ca,說明工質(zhì)的臨界壓力較高也有利于提升系統(tǒng)效率。因此,在候選工質(zhì)的臨界溫度相近時(shí),應(yīng)優(yōu)先選擇臨界壓力較大的工質(zhì)。圖4給出了部分候選工質(zhì)下系統(tǒng)熱回收的圖。其中,工質(zhì)的溫度變化都由一個(gè)線性增長(zhǎng)的預(yù)熱階段和一個(gè)溫度恒定的等溫蒸發(fā)階段構(gòu)成。

    圖3 熱源進(jìn)口溫度200 ℃下不同工質(zhì)對(duì)應(yīng)的最優(yōu)蒸發(fā)溫度和系統(tǒng)效率

    圖4 熱源進(jìn)口溫度200 ℃下部分工質(zhì)的系統(tǒng)熱回收t-Q圖

    Fig 4 Diagram of system heat recovery of some working fluids at waste heat inlet temperature of 200 ℃

    由圖4可以看出,在最優(yōu)條件下,臨界溫度較高的工質(zhì)與熱源流股換熱時(shí)形成的夾點(diǎn)類型為泡點(diǎn)夾點(diǎn)(vaporization pinch point,VPP)[27],工質(zhì)預(yù)熱段的負(fù)荷曲線偏離熱源負(fù)荷曲線較遠(yuǎn),余熱出口溫度較高;隨著工質(zhì)臨界溫度的下降,工質(zhì)的蒸發(fā)溫度升高,工質(zhì)等溫蒸發(fā)部分的負(fù)荷曲線上移,預(yù)熱部分的負(fù)荷曲線逐漸靠近熱源負(fù)荷曲線,當(dāng)工質(zhì)為R245ca時(shí)形成VPP和進(jìn)口夾點(diǎn)(preheating pinch point,PPP)[27]雙夾點(diǎn),隨著工質(zhì)臨界溫度進(jìn)一步降低,工質(zhì)等溫蒸發(fā)部分的負(fù)荷曲線逐漸下移,蒸發(fā)溫度降低,工質(zhì)與熱源僅在進(jìn)口處形成PPP。因此,對(duì)于給定進(jìn)口溫度的熱源,在單壓蒸發(fā)系統(tǒng)的最優(yōu)條件下,工質(zhì)與熱源形成的夾點(diǎn)類型隨著工質(zhì)的臨界溫度降低的變化趨勢(shì)為:VPP→VPP&PPP→PPP,其系統(tǒng)效率隨著工質(zhì)臨界溫度的降低先增大,后減少,并在形成VPP&PPP雙夾點(diǎn)時(shí)達(dá)到最大,此時(shí)回收的熱量最多,產(chǎn)生的蒸氣品質(zhì)也最高。當(dāng)工質(zhì)與熱源形成VPP夾點(diǎn)時(shí),工質(zhì)預(yù)熱部分負(fù)荷曲線與熱源負(fù)荷曲線相距較遠(yuǎn),熱源出口溫度與工質(zhì)進(jìn)口溫度間溫差較大,可通過增加一級(jí)低溫蒸發(fā),進(jìn)一步回收余熱,提升系統(tǒng)輸出和系統(tǒng)效率。然而,當(dāng)工質(zhì)與熱源形成PPP夾點(diǎn)時(shí),僅使用單壓力級(jí)就可將余熱回收完全,對(duì)應(yīng)的熱源出口溫度較低,此時(shí)添加第2級(jí)蒸發(fā)將影響第1級(jí)蒸發(fā)產(chǎn)生的蒸氣品質(zhì)和數(shù)量,從而降低整個(gè)系統(tǒng)的發(fā)電輸出和系統(tǒng)效率。因此,當(dāng)工質(zhì)與熱源形成VPP夾點(diǎn)時(shí),系統(tǒng)具備多級(jí)蒸發(fā)潛力,而形成PPP夾點(diǎn)時(shí),則不具備該潛力。換言之,VPP&PPP雙夾點(diǎn)狀態(tài)既是單壓蒸發(fā)結(jié)構(gòu)下系統(tǒng)所能達(dá)到的最優(yōu)狀態(tài),也是判斷系統(tǒng)是否具有雙壓蒸發(fā)潛力的臨界狀態(tài)。因此,對(duì)于給定進(jìn)口溫度的熱源,若其最優(yōu)工質(zhì)已知,則當(dāng)候選工質(zhì)的臨界溫度高于單壓蒸發(fā)的最優(yōu)工質(zhì)的臨界溫度時(shí),系統(tǒng)可利用多級(jí)蒸發(fā)提升系統(tǒng)的效率和輸出,反之則不可。

    3.2 熱源溫度變化對(duì)單壓蒸發(fā)ORC系統(tǒng)的影響

    以表1中的候選工質(zhì)數(shù)據(jù)為基礎(chǔ),討論熱源進(jìn)口溫度變化對(duì)單壓蒸發(fā)ORC系統(tǒng)工質(zhì)選擇和性能指標(biāo)的影響。設(shè)熱源的進(jìn)口溫度從200 ℃開始,以10 ℃為步長(zhǎng)依次降低直至120 ℃,最終的優(yōu)化結(jié)果分別匯總在表2和圖5中。

    表2 不同熱源進(jìn)口溫度下的最優(yōu)的操作參數(shù)

    圖5 不同熱源進(jìn)口溫度下最大的系統(tǒng)效率和熱回收效率

    由表2的結(jié)果可知,固定熱源進(jìn)口溫度下得到的工質(zhì)選擇原則和夾點(diǎn)變化規(guī)律在熱源進(jìn)口溫度變化的條件下仍然是適用的。即不同熱源溫度下選擇的最優(yōu)工質(zhì)仍然滿足熱源進(jìn)口與工質(zhì)臨界溫度的溫差(26~45 ℃),且熱源流股與工質(zhì)形成VPP&PPP雙夾點(diǎn)時(shí),單壓蒸發(fā)系統(tǒng)效率達(dá)到最大,如200 ℃和160 ℃下的結(jié)果所示。當(dāng)然該狀態(tài)仍可作為判斷系統(tǒng)是否具有雙壓蒸發(fā)潛力的臨界狀態(tài)。

    由圖5也可看出,熱源的入口溫度越高,其系統(tǒng)效率和熱回收效率越高,且系統(tǒng)效率曲線和熱回收效率曲線非常接近。這是因?yàn)?,在最?yōu)操作點(diǎn)下,工質(zhì)的負(fù)荷曲線與熱源流股的負(fù)荷曲線無限接近而形成VPP&PPP雙夾點(diǎn),使得系統(tǒng)效率和熱回收效率也無限接近,且偏離雙夾點(diǎn)狀態(tài)越遠(yuǎn),系統(tǒng)效率和熱回收效率的差距越大。同時(shí),從工質(zhì)的臨界溫度和蒸發(fā)溫度的溫差(5 ℃)也可看出,在最優(yōu)工作狀態(tài)下,工質(zhì)被加熱到臨界溫度附近以獲取盡可能高品質(zhì)的蒸氣,使系統(tǒng)效率最大化。而工質(zhì)的臨界溫度和蒸發(fā)溫度的溫差是為了保證整個(gè)過程的絕對(duì)亞臨界狀態(tài)而給定的工藝限制,見式(3)。

    需要特別指出的是,表2和圖5中都只有200 ℃和160 ℃時(shí)形成的VPP&PPP雙夾點(diǎn),其他熱源溫度下系統(tǒng)效率達(dá)到最大時(shí)并未形成雙夾點(diǎn),而是無限接近雙夾點(diǎn)。這是因?yàn)楹蜻x工質(zhì)的臨界溫度是不連續(xù)的,而表2中熱源流股的進(jìn)口溫度也都是以一定步長(zhǎng)離散的,從而無法保證工質(zhì)的臨界溫度與熱源溫度能夠完全按照“雙夾點(diǎn)規(guī)則”實(shí)現(xiàn)匹配和對(duì)應(yīng),而是盡可能地接近。為了驗(yàn)證這一點(diǎn),對(duì)熱源溫度進(jìn)行微調(diào),以使其與對(duì)應(yīng)的工質(zhì)達(dá)到VPP&PPP雙夾點(diǎn)狀態(tài),如表3所示。

    由表3可知,對(duì)熱源進(jìn)口溫度進(jìn)行微調(diào)后,所有熱源溫度下都可與其對(duì)應(yīng)的工質(zhì)達(dá)到雙夾點(diǎn)最優(yōu)狀態(tài)。同時(shí),在最優(yōu)工況下,工質(zhì)的臨界溫度隨著熱源入口溫度的升高而升高,兩者存在一定的數(shù)學(xué)關(guān)系,擬合結(jié)果如圖6所示。可以看出,最優(yōu)工質(zhì)的臨界溫度與熱源的進(jìn)口溫度基本呈線性關(guān)系,兩者的擬合關(guān)系式為

    由上一節(jié)的分析可知,如果工質(zhì)的臨界溫度高于單壓蒸發(fā)最優(yōu)工質(zhì)的臨界溫度,則系統(tǒng)具有雙壓蒸發(fā)的潛力,反之則不具備雙壓蒸發(fā)的潛力。那么,對(duì)于變熱源進(jìn)口溫度的情況,如果工質(zhì)的臨界溫度位于圖6中曲線下方,則系統(tǒng)不具備雙壓蒸發(fā)的潛力,反之,則系統(tǒng)具備雙壓蒸發(fā)潛力。因此,該線性關(guān)系式可用于任意熱源進(jìn)口溫度下的單壓蒸發(fā)和雙壓蒸發(fā)結(jié)構(gòu)的選擇。

    表3 調(diào)整后不同熱源進(jìn)口溫度下的最優(yōu)的操作參數(shù)

    圖6 熱源進(jìn)口溫度與其最優(yōu)工質(zhì)臨界溫度的關(guān)系

    4 雙壓蒸發(fā)ORC系統(tǒng)分析

    4.1 雙壓蒸發(fā)工藝討論

    雙壓蒸發(fā)ORC結(jié)構(gòu)是在基礎(chǔ)ORC結(jié)構(gòu)的基礎(chǔ)上增加了一個(gè)壓力級(jí),將熱回收分為高壓和低壓兩部分,同時(shí)產(chǎn)生兩種不同品質(zhì)的蒸氣,其結(jié)構(gòu)和-圖如圖2所示。雙壓力結(jié)構(gòu)中換熱部分的-(溫-焓)圖如圖7所示,包含2段預(yù)熱和2段等溫蒸發(fā)部分。

    圖7 雙壓蒸發(fā)ORC熱回收部分的t-H圖

    圖8 不同低壓蒸發(fā)級(jí)蒸發(fā)溫度下最大的系統(tǒng)效率和熱回收效率

    由圖7中工質(zhì)負(fù)荷曲線的拐點(diǎn)可知,雙壓力過程可能產(chǎn)生3類夾點(diǎn),即高壓泡點(diǎn)夾點(diǎn) (HVPP),低壓泡點(diǎn)夾點(diǎn) (LVPP)和PPP。不難看出,在HVPP&LVPP&PPP 3夾點(diǎn)情況下,兩級(jí)蒸發(fā)都能夠產(chǎn)生盡可能高品質(zhì)的蒸氣,系統(tǒng)的熱回收也最完全,是雙壓蒸發(fā)循環(huán)的最理想的工況點(diǎn)。然而,由單壓力級(jí)情況的分析可知,在僅形成VPP時(shí)系統(tǒng)才具備2級(jí)蒸發(fā)的潛力。換言之,選擇雙壓力結(jié)構(gòu)的前提是單壓蒸發(fā)系統(tǒng)形成VPP,且無法形成PPP。因此,對(duì)于適用于雙壓蒸發(fā)系統(tǒng)的流體來說,其最優(yōu)工況只可能是同時(shí)形成LVPP&HVPP雙夾點(diǎn)。那么,雙壓力系統(tǒng)中低壓力回路的功能就是在保證蒸氣品質(zhì)的基礎(chǔ)上,盡可能增加熱回收,以提升總的系統(tǒng)效率。

    4.2 低壓級(jí)蒸發(fā)溫度對(duì)系統(tǒng)的影響

    在固定高壓力級(jí)蒸發(fā)溫度的基礎(chǔ)上,分析低壓蒸發(fā)級(jí)蒸發(fā)溫度對(duì)系統(tǒng)性能的影響。本節(jié)選取熱源進(jìn)口溫度為200 ℃,工質(zhì)為R141b,高壓力級(jí)的蒸發(fā)溫度為140 ℃。分析結(jié)果如圖8所示。換熱部分的圖見圖9。

    圖9 雙壓蒸發(fā)結(jié)構(gòu)不同低壓蒸發(fā)溫度對(duì)應(yīng)的換熱回收部分的t-Q圖

    由圖8可以看出,雙壓蒸發(fā)ORC結(jié)構(gòu)下的系統(tǒng)效率隨著低壓蒸發(fā)級(jí)的蒸發(fā)溫度的增加先增大后減小,存在最大的系統(tǒng)效率,為0.124 4,較單壓蒸發(fā)循環(huán)下的系統(tǒng)效率提升20.78%。需要注意的是,系統(tǒng)效率變化曲線兩個(gè)端點(diǎn)處的低壓蒸發(fā)溫度分別為冷凝溫度(40 ℃)和高壓蒸發(fā)溫度(140 ℃),對(duì)應(yīng)的ORC系統(tǒng)為單壓蒸發(fā)循環(huán)結(jié)構(gòu)。因此,該變化曲線的含義是:當(dāng)?shù)蛪赫舭l(fā)溫度高于40 ℃時(shí),系統(tǒng)可考慮添加低壓蒸發(fā)級(jí),且隨著低壓蒸發(fā)溫度升高,雙壓蒸發(fā)ORC系統(tǒng)效率逐漸提升,在低壓蒸發(fā)溫度為90 ℃時(shí)達(dá)到最大;此時(shí)若繼續(xù)升高低壓蒸發(fā)溫度,雙壓蒸發(fā)ORC的系統(tǒng)效率將開始降低,直至低壓蒸發(fā)溫度達(dá)到140 ℃時(shí),系統(tǒng)重新變?yōu)閱螇毫?jí)結(jié)構(gòu)。除此之外,雙壓蒸發(fā)ORC結(jié)構(gòu)下的系統(tǒng)效率隨低壓蒸發(fā)級(jí)蒸發(fā)溫度的變化幅度不大,不超過13.7%。換言之,在雙壓蒸發(fā)ORC結(jié)構(gòu)中,低壓蒸發(fā)級(jí)對(duì)系統(tǒng)性能的影響較小,高壓力級(jí)才是決定雙壓蒸發(fā)循環(huán)熱回收性能的關(guān)鍵。

    由圖8中還可看出,系統(tǒng)的熱回收效率隨著低壓蒸發(fā)級(jí)的蒸發(fā)溫度的提升而增大。這是因?yàn)殡S著低壓蒸發(fā)級(jí)蒸發(fā)溫度的升高, 系統(tǒng)的熱回收量減少,使得系統(tǒng)的熱回收效率得以提升,如圖9所示。圖9中藍(lán)色曲線為高壓部分的工質(zhì)負(fù)荷曲線,黑色曲線為低壓部分的工質(zhì)負(fù)荷曲線。由圖9可知,雙壓蒸發(fā)循環(huán)相比單壓蒸發(fā)循環(huán)可回收更多的余熱,降低熱源出口溫度,提升熱回收效率。低壓蒸發(fā)級(jí)蒸發(fā)溫度在40~140 ℃內(nèi)變化時(shí),雙壓蒸發(fā)循環(huán)可比單壓蒸發(fā)循環(huán)多回收204 kW余熱熱量,占總的可用熱量的27.9%,且熱回收量隨著低壓蒸發(fā)級(jí)蒸發(fā)溫度的升高而減少,低壓蒸發(fā)級(jí)蒸發(fā)溫度越高,熱源的出口溫度越高,回收的熱量越少,使得系統(tǒng)的熱回收效率越高。

    5 結(jié) 論

    針對(duì)單壓和雙壓蒸發(fā)ORC系統(tǒng)的工質(zhì)和結(jié)構(gòu)選擇問題,在熱源進(jìn)口溫度為200~120 ℃,選取了17種候選工質(zhì),對(duì)單壓和雙壓蒸發(fā)ORC結(jié)構(gòu)適用條件進(jìn)行模擬研究,得到如下結(jié)論:

    (1) 對(duì)于單壓蒸發(fā)ORC結(jié)構(gòu),隨著工質(zhì)臨界溫度的降低,工質(zhì)與熱源流股在最優(yōu)情況下形成的夾點(diǎn)類型由VPP逐漸變?yōu)閂PP&PPP雙夾點(diǎn),并最終形成單獨(dú)的PPP;在形成VPP&PPP雙夾點(diǎn)時(shí),系統(tǒng)效率達(dá)到最大。

    (2) 熱源和工質(zhì)形成的夾點(diǎn)類型是ORC系統(tǒng)結(jié)構(gòu)選擇的關(guān)鍵指標(biāo)。在單壓蒸發(fā)最優(yōu)條件下,夾點(diǎn)類型為VPP時(shí),可選擇雙壓力結(jié)構(gòu)進(jìn)一步提升系統(tǒng)效率,而夾點(diǎn)類型為PPP時(shí),系統(tǒng)并不具備通過雙壓蒸發(fā)結(jié)構(gòu)提升效率的潛能,采用單壓蒸發(fā)結(jié)構(gòu)更優(yōu)。

    (3) 不同熱源條件下,熱源進(jìn)口溫度與其單壓蒸發(fā)結(jié)構(gòu)下的最優(yōu)工質(zhì)的臨界溫度呈線性關(guān)系,該關(guān)系式也可作為任意熱源進(jìn)口溫度下ORC系統(tǒng)結(jié)構(gòu)選擇的依據(jù)。

    (4) 雙壓蒸發(fā)循環(huán)可有效地提升系統(tǒng)的熱回收量,降低熱源出口溫度,提升系統(tǒng)效率和發(fā)電輸出。對(duì)于雙壓蒸發(fā)結(jié)構(gòu)的系統(tǒng),其系統(tǒng)效率最大時(shí)對(duì)應(yīng)的換熱夾點(diǎn)類型為HVPP&LVPP雙夾點(diǎn)。

    (5) 在雙壓蒸發(fā)ORC結(jié)構(gòu)中,高壓蒸發(fā)級(jí)是決定熱回收性能的關(guān)鍵,低壓蒸發(fā)級(jí)對(duì)系統(tǒng)性能的影響較小。

    [1] ZIVIANI D, BEYENE A, VENTURINI M. Advances and challenges in ORC systems modeling for low grade thermal energy recovery [J]. Applied Energy, 2014, 121: 79-95.

    [2] SALEH B, KOGLBAUER G, WENDLAND M,. Working fluids for low-temperature organic Rankine cycles [J]. Energy, 2007, 32(7): 1210-1221.

    [3] MAGO P J, CHAMRA L M, SRINIVASAN K K,. An examination of regenerative organic Rankine cycles using dry fluids [J]. Applied Thermal Engineering, 2008, 28(8): 998-1007.

    [4] DU Y, YANG Y, HU D,. Off-design performance comparative analysis between basic and parallel dual-pressure organic Rankine cycles using radial inflow turbines [J]. Applied Thermal Engineering, 2018, 138: 18-34.

    [5] HO T, MAO S S, GREIF R. Increased power production through enhancements to the organic flash cycle (OFC) [J]. Energy, 2012, 45(1): 686-695.

    [6] LE V L, FEIDT M, KHEIRI A,. Performance optimization of low-temperature power generation by supercritical ORCs (organic Rankine cycles) using low GWP (global warming potential) working fluids [J]. Energy, 2014; 67: 513-526.

    [7] LECOMPTE S, HUISSEUNE H, BROEK M,. Review of organic Rankine cycle (ORC) architectures for waste heat recovery [J]. Renewable & Sustainable Energy Reviews, 2015; 47: 448-461.

    [8] LI T L, WANG Q, ZHU J,. Thermodynamic optimization of organic Rankine cycle using two-stage evaporation [J]. Renewable Energy, 2015, 75(C): 654-664.

    [9] LI T L, ZHANG Z, LU J,. Two-stage evaporation strategy to improve system performance for organic Rankine cycle [J]. Applied Energy, 2015, 150: 323-334.

    [10] LI T L, MENG N, LIU J,. Thermodynamic and economic evaluation of the organic Rankine cycle (ORC) and two-stage series organic Rankine cycle (TSORC) for flue gas heat recovery [J]. Energy Conversion and Management, 2019, 183: 816-829.

    [11] LI J, YANG Z, HU S,. Thermo-economic analyses and evaluations of small-scale dual-pressure evaporation organic Rankine cycle system using pure fluids [J]. Energy, 2020, 206: 413-431.

    [12] MANENTE G, LAZZARETTO A, BONAMICO E. Design guidelines for the choice between single and dual pressure layouts in organic Rankine cycle (ORC) systems [J]. Energy, 2017, 123: 413-431.

    [13] LI J, GE Z, DUAN Y,. Parametric optimization and thermodynamic performance comparison of single-pressure and dual-pressure evaporation organic Rankine cycles [J]. Applied Energy, 2018, 217: 409-421.

    [14] LI J, GE Z, DUAN Y,. Design and performance analyses for a novel organic Rankine cycle with supercritical-subcritical heat absorption process coupling [J]. Applied Energy, 2019, 235: 1400-1414.

    [15] WANG M, ZHANG J, ZHAO S,. Performance investigation of transcritical and dual-pressure organic Rankine cycles from the aspect of thermal match [J]. Energy Conversion and Management, 2019, 197: 111850.

    [16] LI J, GE Z, LIU Q,. Thermo-economic performance analyses and comparison of two turbine layouts for organic Rankine cycles with dual-pressure evaporation [J]. Energy Conversion and Management, 2018, 164: 603-614.

    [17] CHEN H J, GOSWAMI D Y, STEFANAKOS E K. A review of thermodynamic cycles and working fluids for the conversion of low-grade heat [J]. Renewable & Sustainable Energy Reviews, 2010, 14(9): 3059-3067.

    [18] BROWN J S, BRIGNOLI R, QUINE T. Parametric investigation of working fluids for organic Rankine cycle applications [J]. Applied Thermal Engineering, 2015, 90: 64-74.

    [19] LONG R, BAO Y J, HUANG X M,. Exergy analysis and working fluid selection of organic Rankine cycle for low grade waste heat recovery [J]. Energy, 2014, 73: 475-483.

    [20] HAERVIG J, SORENSEN K, CONDRA T J. Guidelines for optimal selection of working fluid for an organic Rankine cycle in relation to waste heat recovery [J]. Energy, 2016, 96: 592-602.

    [21] VIVIAN J, MANENTE G, LAZZARETTO A. A general framework to select working fluid and configuration of ORCs for low-to-medium temperature heat sources [J]. Applied Energy, 2015, 156: 727-746.

    [22] PENG D, ROBINSON D B. A new two-constant equation of state [J]. Industrial & Engineering Chemistry Fundamentals, 1976, 15(1): 59-64.

    [23] KAMATH R S, BIEGLER L T, GROSSMANN I E. An equation-oriented approach for handling thermodynamics based on cubic equation of state in process optimization [J]. Computers & Chemical Engineering, 2010, 34(12): 2085-2096.

    [24] CHEN C L, LI P Y, LE S N T. Organic rankine cycle for waste heat recovery in a refinery [J]. Industrial & Engineering Chemistry Research, 2016, 55(12): 3262–3275.

    [25] YU H S, EASON J, BIEGLER L T,. Process integration and superstructure optimization of organic Rankine cycles (ORCs) with heat exchanger network synthesis [J]. Computers & Chemical Engineering, 2017, 107(S1): 257-270.

    [26] TCHANCHE B F, PAPADAKIS G, LAMBRINOS G,. Fluid selection for a low-temperature solar organic Rankine cycle [J]. Applied Thermal Engineering, 2009, 29(11/12): 2468-2476.

    [27] YU H S, FENG X, WANG Y F. A new pinch based method for simultaneous selection of working fluid and operating conditions in an ORC (organic Rankine cycle) recovering waste heat [J]. Energy, 2015, 90: 36-46.

    Simulation-based structure selection of single and dual pressure evaporation organic Rankine cycles

    TANG Jian-ping1, KANG Li-xia1,2, LIU Yong-zhong1,2

    (1. Department of Chemical Engineering, Xi’an Jiaotong University, Xi’an 710049, China;2. Shaanxi Key Laboratory of Energy Chemical Process Intensification, Xi’an 710049, China)

    For the selection of the working fluid and structure of single and dual evaporation organic Rankine cycle (ORC), a simulation analysis combined with pinch theory by Aspen Plus was conducted to explore the thermodynamic performance and heat transfer process of single and dual pressure evaporation ORCs with different working fluids under the inlet temperature of waste heat at 120-200 ℃. Results show that the type of pinch point formed by the waste heat and the working fluid was the key to the structure selection of ORC. For single pressure evaporation cycle, with the decrease of critical temperature of working fluids, the pinch point of working fluid and waste heat under the optimal conditions gradually changed from vaporization pinch point (VPP) to vaporization and preheaing pinch points (VPP&PPP), and finally formed a preheating pinch point (PPP). The system efficiency of single pressure evaporation ORC reached the maximum when VPP&PPP was formed. Besides, VPP&PPP were also the critical point for the selection of single and dual evaporation structures, and a quantitative relationship was then obtained for the efficient selection of the working fluid and ORC structure in actual industrial process.

    low-temperature waste heat; organic Rankine cycle; dual pressure evaporation structure; process simulation; pinch point

    1003-9015(2021)05-0896-09

    TQ 021.8

    A

    10.3969/j.issn.1003-9015.2021.05.017

    2020-11-30;

    2021-01-23。

    國(guó)家自然科學(xué)基金(21808179,21878240)。

    唐建平(1996-),男,湖北葉堰人,西安交通大學(xué)碩士生。通信聯(lián)系人:康麗霞,E-mail:lx_kang@mail.xjtu.edu.cn

    猜你喜歡
    夾點(diǎn)臨界溫度工質(zhì)
    海洋溫差能發(fā)電熱力循環(huán)系統(tǒng)的工質(zhì)優(yōu)選
    ORC系統(tǒng)蒸發(fā)器夾點(diǎn)溫差的數(shù)值分析和實(shí)驗(yàn)研究
    Bogoliubov-Tolmachev-Shirkov模型臨界溫度和能隙解的數(shù)值方法
    采用R1234ze(E)/R245fa的非共沸混合工質(zhì)有機(jī)朗肯循環(huán)系統(tǒng)實(shí)驗(yàn)研究
    采用二元非共沸工質(zhì)的有機(jī)朗肯循環(huán)熱力學(xué)分析
    若干低GWP 純工質(zhì)在空調(diào)系統(tǒng)上的應(yīng)用分析
    機(jī)械A(chǔ)uto CAD教學(xué)中夾點(diǎn)操作的探究
    RDX基炸藥熱起爆臨界溫度的測(cè)試及數(shù)值計(jì)算
    高于臨界溫度的頁(yè)巖吸附甲烷數(shù)據(jù)預(yù)測(cè)
    以色列科學(xué)家開發(fā)出光控超導(dǎo)材料可通過不同類型的光照改變臨界溫度
    亚洲午夜精品一区,二区,三区| 人妻久久中文字幕网| 精品久久久久久久久久久久久 | 午夜福利免费观看在线| 天堂动漫精品| 国产欧美日韩精品亚洲av| 嫩草影院精品99| 国产高清有码在线观看视频 | 国产精品 国内视频| 亚洲国产精品成人综合色| 母亲3免费完整高清在线观看| 日本成人三级电影网站| 国产精品免费视频内射| 色婷婷久久久亚洲欧美| 人人妻人人看人人澡| 老司机福利观看| 丰满人妻熟妇乱又伦精品不卡| 黑人欧美特级aaaaaa片| 成人一区二区视频在线观看| 成人永久免费在线观看视频| 成年版毛片免费区| 夜夜看夜夜爽夜夜摸| 又黄又爽又免费观看的视频| 国产亚洲av高清不卡| 老司机靠b影院| 亚洲精品一区av在线观看| 中文字幕人成人乱码亚洲影| 久久中文看片网| 精品久久久久久久毛片微露脸| 啪啪无遮挡十八禁网站| 母亲3免费完整高清在线观看| 成年女人毛片免费观看观看9| 99久久99久久久精品蜜桃| 欧美一区二区精品小视频在线| 亚洲男人天堂网一区| 国产成人影院久久av| 欧美国产日韩亚洲一区| 亚洲 欧美 日韩 在线 免费| 欧美激情 高清一区二区三区| 久久精品成人免费网站| 国产精品一区二区免费欧美| 不卡一级毛片| 日本精品一区二区三区蜜桃| 国产麻豆成人av免费视频| 精品国产乱子伦一区二区三区| 午夜两性在线视频| 国产精品,欧美在线| 精品久久久久久久末码| 亚洲在线自拍视频| 精品一区二区三区视频在线观看免费| 国产一区二区在线av高清观看| 好男人在线观看高清免费视频 | 一区二区日韩欧美中文字幕| 色精品久久人妻99蜜桃| 国产人伦9x9x在线观看| 欧美三级亚洲精品| 国产三级黄色录像| 国产97色在线日韩免费| 少妇熟女aⅴ在线视频| 久久久久亚洲av毛片大全| 成人免费观看视频高清| x7x7x7水蜜桃| 中文字幕精品免费在线观看视频| 国产精品香港三级国产av潘金莲| 国产精品久久视频播放| 精品久久久久久久久久免费视频| 久久久久久九九精品二区国产 | 99精品久久久久人妻精品| 欧美最黄视频在线播放免费| 亚洲真实伦在线观看| 国产一区二区三区在线臀色熟女| 中出人妻视频一区二区| 国产亚洲精品久久久久久毛片| 国产91精品成人一区二区三区| 国内毛片毛片毛片毛片毛片| 熟女少妇亚洲综合色aaa.| 熟妇人妻久久中文字幕3abv| 久久婷婷人人爽人人干人人爱| 人成视频在线观看免费观看| 国产高清videossex| 一级黄色大片毛片| 免费av毛片视频| 亚洲精品久久国产高清桃花| 看免费av毛片| 18禁观看日本| 国产成人精品无人区| 一卡2卡三卡四卡精品乱码亚洲| 不卡av一区二区三区| 久久狼人影院| 青草久久国产| netflix在线观看网站| 日本五十路高清| 91av网站免费观看| 国产精品永久免费网站| 国产成人欧美| 欧美日韩中文字幕国产精品一区二区三区| www国产在线视频色| 亚洲午夜精品一区,二区,三区| 少妇被粗大的猛进出69影院| 成年人黄色毛片网站| 怎么达到女性高潮| 亚洲成人精品中文字幕电影| a在线观看视频网站| 国产在线精品亚洲第一网站| 男人舔奶头视频| 亚洲性夜色夜夜综合| 国产97色在线日韩免费| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| a级毛片在线看网站| 中文在线观看免费www的网站 | 亚洲国产欧洲综合997久久, | 美女午夜性视频免费| 一级片免费观看大全| 美女扒开内裤让男人捅视频| 身体一侧抽搐| 一夜夜www| 特大巨黑吊av在线直播 | 狠狠狠狠99中文字幕| 国产极品粉嫩免费观看在线| 50天的宝宝边吃奶边哭怎么回事| 久久精品亚洲精品国产色婷小说| 亚洲五月色婷婷综合| 日韩中文字幕欧美一区二区| av电影中文网址| 亚洲黑人精品在线| 99精品欧美一区二区三区四区| 欧美日韩乱码在线| 欧美日韩福利视频一区二区| 中文字幕最新亚洲高清| 熟妇人妻久久中文字幕3abv| 91成人精品电影| 精华霜和精华液先用哪个| 亚洲熟妇熟女久久| 久久精品91无色码中文字幕| 欧美中文综合在线视频| 亚洲色图av天堂| av欧美777| 91av网站免费观看| 欧美中文日本在线观看视频| 成人亚洲精品av一区二区| 国产精品亚洲美女久久久| 18禁观看日本| 听说在线观看完整版免费高清| 午夜亚洲福利在线播放| 视频在线观看一区二区三区| 黄色片一级片一级黄色片| 亚洲第一欧美日韩一区二区三区| 国产久久久一区二区三区| 夜夜躁狠狠躁天天躁| 欧美性长视频在线观看| av天堂在线播放| 国产精品一区二区三区四区久久 | 久久久久国产一级毛片高清牌| 国产免费男女视频| 在线视频色国产色| 亚洲av第一区精品v没综合| 女人被狂操c到高潮| 欧美乱色亚洲激情| 丁香欧美五月| 亚洲精品中文字幕一二三四区| 狂野欧美激情性xxxx| 超碰成人久久| 亚洲精品一区av在线观看| 叶爱在线成人免费视频播放| 国产成人一区二区三区免费视频网站| 波多野结衣巨乳人妻| 午夜免费成人在线视频| 日本免费a在线| 亚洲第一欧美日韩一区二区三区| 欧美激情久久久久久爽电影| 欧美在线黄色| а√天堂www在线а√下载| www国产在线视频色| 国产伦在线观看视频一区| 亚洲真实伦在线观看| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区四区五区乱码| 精品国产国语对白av| 精品人妻1区二区| 香蕉丝袜av| 国产成人欧美在线观看| 欧美性长视频在线观看| 99re在线观看精品视频| 欧美成人性av电影在线观看| 国产激情久久老熟女| 黑人巨大精品欧美一区二区mp4| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 成人国语在线视频| 岛国在线观看网站| 中文字幕精品免费在线观看视频| 男人舔女人的私密视频| 在线永久观看黄色视频| 一区二区三区精品91| 亚洲欧洲精品一区二区精品久久久| 日本成人三级电影网站| 在线免费观看的www视频| 欧美成狂野欧美在线观看| 免费搜索国产男女视频| 成年女人毛片免费观看观看9| 欧美绝顶高潮抽搐喷水| 日本五十路高清| 国产一区二区三区在线臀色熟女| 99热6这里只有精品| 免费在线观看日本一区| 久热爱精品视频在线9| 亚洲欧美日韩无卡精品| 日韩视频一区二区在线观看| 青草久久国产| 女人被狂操c到高潮| 欧美日本视频| 久久精品人妻少妇| 老汉色av国产亚洲站长工具| 十八禁人妻一区二区| 欧美国产日韩亚洲一区| 欧美日韩精品网址| 国产真人三级小视频在线观看| 69av精品久久久久久| 特大巨黑吊av在线直播 | 国产一区二区三区在线臀色熟女| 亚洲人成电影免费在线| 最新美女视频免费是黄的| 精品国产亚洲在线| 国产成年人精品一区二区| 成年人黄色毛片网站| 亚洲av中文字字幕乱码综合 | 在线观看午夜福利视频| 亚洲国产中文字幕在线视频| 日韩精品中文字幕看吧| 久久99热这里只有精品18| 1024视频免费在线观看| 国内久久婷婷六月综合欲色啪| 欧美日韩瑟瑟在线播放| 观看免费一级毛片| 亚洲精品在线美女| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 老司机午夜福利在线观看视频| 中文字幕精品亚洲无线码一区 | 老熟妇乱子伦视频在线观看| 国产又爽黄色视频| 男女那种视频在线观看| 国产精品1区2区在线观看.| 久热这里只有精品99| 亚洲黑人精品在线| 亚洲精品国产一区二区精华液| 免费观看精品视频网站| 精品久久久久久成人av| 国产视频一区二区在线看| 亚洲无线在线观看| 欧美国产日韩亚洲一区| 曰老女人黄片| 一二三四社区在线视频社区8| 欧美黑人巨大hd| 老司机午夜十八禁免费视频| 日韩国内少妇激情av| 国产一区在线观看成人免费| 黄片大片在线免费观看| 成人精品一区二区免费| 亚洲专区国产一区二区| 亚洲国产欧洲综合997久久, | 国产精品亚洲美女久久久| 99久久久亚洲精品蜜臀av| 叶爱在线成人免费视频播放| 精品无人区乱码1区二区| 日日爽夜夜爽网站| 亚洲自偷自拍图片 自拍| 亚洲色图av天堂| 国产国语露脸激情在线看| 50天的宝宝边吃奶边哭怎么回事| 一级黄色大片毛片| 亚洲男人天堂网一区| 精品久久蜜臀av无| 亚洲三区欧美一区| 欧美黑人欧美精品刺激| 欧美丝袜亚洲另类 | 男女做爰动态图高潮gif福利片| 国产一区二区激情短视频| 最新在线观看一区二区三区| 亚洲人成网站在线播放欧美日韩| 国产精品一区二区三区四区久久 | 欧美日韩黄片免| 久久久久久大精品| 女人被狂操c到高潮| 国产免费男女视频| 色综合亚洲欧美另类图片| 丝袜美腿诱惑在线| 国产97色在线日韩免费| 亚洲色图av天堂| 国产精品久久久人人做人人爽| 成人午夜高清在线视频 | 国产成人系列免费观看| 日韩视频一区二区在线观看| aaaaa片日本免费| 国产伦一二天堂av在线观看| 波多野结衣av一区二区av| 一级片免费观看大全| 亚洲在线自拍视频| 久久久久久九九精品二区国产 | 欧美日韩亚洲综合一区二区三区_| 91麻豆av在线| 国产精品99久久99久久久不卡| 国产91精品成人一区二区三区| 麻豆成人午夜福利视频| 国产熟女xx| 欧美黑人精品巨大| 欧美色视频一区免费| 黑人操中国人逼视频| 久久精品aⅴ一区二区三区四区| 高清毛片免费观看视频网站| 国产99久久九九免费精品| 久久人人精品亚洲av| 一级片免费观看大全| 夜夜躁狠狠躁天天躁| 国产精品影院久久| 老鸭窝网址在线观看| 国产精品一区二区免费欧美| 国产精品影院久久| 亚洲片人在线观看| 免费搜索国产男女视频| 啦啦啦韩国在线观看视频| 日本五十路高清| 国产成人精品无人区| 久久午夜综合久久蜜桃| 12—13女人毛片做爰片一| 成在线人永久免费视频| 欧美国产日韩亚洲一区| 亚洲全国av大片| 国产成人精品久久二区二区91| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 国产一区二区三区在线臀色熟女| 亚洲精品国产一区二区精华液| 999精品在线视频| 一夜夜www| 黑人欧美特级aaaaaa片| 最好的美女福利视频网| 日韩av在线大香蕉| 人妻丰满熟妇av一区二区三区| 欧美av亚洲av综合av国产av| 丝袜人妻中文字幕| 窝窝影院91人妻| 精品国产一区二区三区四区第35| 午夜福利视频1000在线观看| 91在线观看av| 国产av一区二区精品久久| x7x7x7水蜜桃| 亚洲美女黄片视频| 精品久久蜜臀av无| 国产精品香港三级国产av潘金莲| 午夜成年电影在线免费观看| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 丁香欧美五月| 亚洲男人的天堂狠狠| 国产又爽黄色视频| 午夜久久久在线观看| 日韩精品免费视频一区二区三区| 亚洲国产精品久久男人天堂| 久久99热这里只有精品18| 久久伊人香网站| 日本 av在线| 成人亚洲精品av一区二区| 99re在线观看精品视频| 日韩欧美在线二视频| 国产高清激情床上av| 校园春色视频在线观看| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 宅男免费午夜| 啦啦啦观看免费观看视频高清| 亚洲国产欧洲综合997久久, | 久久人妻福利社区极品人妻图片| 欧美性长视频在线观看| 精品免费久久久久久久清纯| 正在播放国产对白刺激| 国产精品99久久99久久久不卡| 看黄色毛片网站| 色综合婷婷激情| 18禁黄网站禁片免费观看直播| 人成视频在线观看免费观看| 啦啦啦韩国在线观看视频| 男人舔女人下体高潮全视频| 美女 人体艺术 gogo| 在线十欧美十亚洲十日本专区| 黄色视频不卡| 久久伊人香网站| 99国产极品粉嫩在线观看| 中国美女看黄片| 又黄又粗又硬又大视频| 亚洲熟妇中文字幕五十中出| 日韩欧美国产在线观看| 丁香欧美五月| 97超级碰碰碰精品色视频在线观看| 亚洲三区欧美一区| 久热这里只有精品99| 成在线人永久免费视频| 国产av又大| 精品国内亚洲2022精品成人| 久久人人精品亚洲av| 亚洲九九香蕉| 天堂影院成人在线观看| 免费看a级黄色片| 人成视频在线观看免费观看| 在线国产一区二区在线| 精品乱码久久久久久99久播| 精品欧美一区二区三区在线| 欧美日韩黄片免| 69av精品久久久久久| 成人精品一区二区免费| 在线免费观看的www视频| www.精华液| 美女午夜性视频免费| 精品欧美一区二区三区在线| 成人三级黄色视频| 免费在线观看成人毛片| 精品不卡国产一区二区三区| 露出奶头的视频| 亚洲人成网站高清观看| 亚洲国产毛片av蜜桃av| 成人手机av| 亚洲色图av天堂| 天天添夜夜摸| 50天的宝宝边吃奶边哭怎么回事| 日韩精品免费视频一区二区三区| x7x7x7水蜜桃| 妹子高潮喷水视频| 欧美黑人巨大hd| 一本综合久久免费| 国产一区二区激情短视频| 久久国产精品人妻蜜桃| 日韩欧美一区视频在线观看| 中文字幕精品亚洲无线码一区 | 午夜福利欧美成人| 久久久久国产一级毛片高清牌| 国产男靠女视频免费网站| 无人区码免费观看不卡| 麻豆久久精品国产亚洲av| 精品欧美国产一区二区三| 国产区一区二久久| 久久性视频一级片| 天堂影院成人在线观看| 国产又爽黄色视频| 精品一区二区三区四区五区乱码| 老司机靠b影院| 嫁个100分男人电影在线观看| 日本成人三级电影网站| 免费人成视频x8x8入口观看| 国产免费男女视频| 成在线人永久免费视频| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看| 国产久久久一区二区三区| 美女免费视频网站| 精品一区二区三区四区五区乱码| 一进一出抽搐gif免费好疼| 久久九九热精品免费| 女同久久另类99精品国产91| 免费人成视频x8x8入口观看| 日本熟妇午夜| 免费看美女性在线毛片视频| 热re99久久国产66热| 在线天堂中文资源库| 最近最新中文字幕大全电影3 | 色播亚洲综合网| 一二三四在线观看免费中文在| 熟女电影av网| 亚洲国产精品成人综合色| netflix在线观看网站| 国产午夜精品久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品在线观看二区| 18美女黄网站色大片免费观看| 国产精品99久久99久久久不卡| 听说在线观看完整版免费高清| 精品高清国产在线一区| 久久久久久国产a免费观看| 亚洲在线自拍视频| 亚洲精品色激情综合| 成人18禁高潮啪啪吃奶动态图| 在线永久观看黄色视频| 午夜成年电影在线免费观看| 午夜a级毛片| 人妻久久中文字幕网| 精华霜和精华液先用哪个| 桃色一区二区三区在线观看| 男人舔女人的私密视频| 国产麻豆成人av免费视频| 18禁观看日本| bbb黄色大片| 一进一出抽搐gif免费好疼| 国产真实乱freesex| 亚洲精品中文字幕在线视频| 成人手机av| 黄色 视频免费看| 麻豆国产av国片精品| 国产精品综合久久久久久久免费| 国产欧美日韩一区二区精品| 国产成人欧美| 99在线视频只有这里精品首页| 亚洲国产欧美日韩在线播放| 在线观看日韩欧美| 久久久久久九九精品二区国产 | 高清在线国产一区| 久久中文看片网| 怎么达到女性高潮| 亚洲精品粉嫩美女一区| 黄色毛片三级朝国网站| 亚洲欧美一区二区三区黑人| 久久精品91蜜桃| 很黄的视频免费| 久久婷婷人人爽人人干人人爱| 日日爽夜夜爽网站| 免费搜索国产男女视频| 亚洲专区中文字幕在线| 精品日产1卡2卡| 中文字幕人妻丝袜一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 免费电影在线观看免费观看| 脱女人内裤的视频| av视频在线观看入口| 在线观看www视频免费| 亚洲欧美激情综合另类| 国产一卡二卡三卡精品| 色播在线永久视频| 丝袜人妻中文字幕| 国产三级在线视频| 久久午夜综合久久蜜桃| 国内毛片毛片毛片毛片毛片| 亚洲精品在线美女| 老司机午夜福利在线观看视频| 午夜免费观看网址| 制服人妻中文乱码| 老汉色av国产亚洲站长工具| 欧美国产日韩亚洲一区| 亚洲天堂国产精品一区在线| 亚洲片人在线观看| 成人亚洲精品av一区二区| 亚洲国产精品999在线| 99久久国产精品久久久| 韩国精品一区二区三区| 免费在线观看亚洲国产| 久久国产乱子伦精品免费另类| 女人被狂操c到高潮| 啦啦啦免费观看视频1| 国产成人啪精品午夜网站| 黄色女人牲交| 制服人妻中文乱码| 亚洲熟妇熟女久久| 亚洲av片天天在线观看| 听说在线观看完整版免费高清| 国产精品乱码一区二三区的特点| 欧美黑人巨大hd| 国产99白浆流出| 国产麻豆成人av免费视频| 国产男靠女视频免费网站| 国产v大片淫在线免费观看| 久久精品aⅴ一区二区三区四区| 两个人看的免费小视频| 精品国产美女av久久久久小说| 两个人看的免费小视频| 国产av不卡久久| 亚洲午夜理论影院| 满18在线观看网站| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 黑人巨大精品欧美一区二区mp4| 久久久久九九精品影院| 欧美一区二区精品小视频在线| 国内精品久久久久久久电影| 美女免费视频网站| 91大片在线观看| 国产三级黄色录像| 成人永久免费在线观看视频| 大型av网站在线播放| 国语自产精品视频在线第100页| 中亚洲国语对白在线视频| 国产精品一区二区精品视频观看| 日韩高清综合在线| 黄色 视频免费看| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品色激情综合| 精品国产国语对白av| 亚洲一码二码三码区别大吗| 一二三四在线观看免费中文在| 最近在线观看免费完整版| 伊人久久大香线蕉亚洲五| 亚洲第一欧美日韩一区二区三区| 久久精品国产亚洲av香蕉五月| 怎么达到女性高潮| av在线天堂中文字幕| 亚洲一区中文字幕在线| 日韩大码丰满熟妇| 日韩欧美 国产精品| 十分钟在线观看高清视频www| 在线观看一区二区三区| 男女那种视频在线观看| 视频区欧美日本亚洲| 老司机深夜福利视频在线观看| 免费看十八禁软件| 黄色 视频免费看| 首页视频小说图片口味搜索| 国产精品久久久人人做人人爽| 国产99久久九九免费精品| 少妇裸体淫交视频免费看高清 | 国产高清videossex| 男女之事视频高清在线观看| 精品电影一区二区在线| 久久人妻av系列| 日韩av在线大香蕉| 免费看a级黄色片| 国产精品电影一区二区三区|