亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        甲苯歧化及烷基轉移單元用能優(yōu)化與改進

        2016-10-20 02:51:56陳清林
        石油學報(石油加工) 2016年5期
        關鍵詞:優(yōu)化

        陳 婷,陳清林

        (1.江漢大學 武漢研究院,湖北 武漢 430056; 2.江漢大學 化學與環(huán)境工程學院,湖北 武漢 430056; 3.中山大學 化學與化學工程學院 低碳化學與過程節(jié)能廣東省重點實驗室,廣東 廣州 510275)

        ?

        甲苯歧化及烷基轉移單元用能優(yōu)化與改進

        陳婷1,2,陳清林3

        (1.江漢大學 武漢研究院,湖北 武漢 430056; 2.江漢大學 化學與環(huán)境工程學院,湖北 武漢 430056; 3.中山大學 化學與化學工程學院 低碳化學與過程節(jié)能廣東省重點實驗室,廣東 廣州 510275)

        在流程模擬軟件Aspen Plus對甲苯歧化及烷基轉移單元進行流程模擬的基礎上,采用三環(huán)節(jié)能量結構模型對整個單元進行用能分析及評價,并運用該分析所揭示的用能瓶頸合理設計過程的換熱網(wǎng)絡(HEN),以減少或消除過程重復加熱或冷卻,降低不合理損,有效實現(xiàn)過程節(jié)能。基于改造后的換熱網(wǎng)絡,對分餾塔的操作參數(shù),包括進料溫度和位置進行優(yōu)化;基于過程能量分析,提出有效節(jié)能改進措施,提高整個甲苯歧化及烷基轉移單元的能源利用效率約為6%。

        三環(huán)節(jié)能量結構模型;甲苯歧化及烷基轉移單元;流程模擬;節(jié)能

        芳烴是石化工業(yè)中規(guī)模和產(chǎn)量僅次于乙烯和丙烯的重要有機化工原料,其衍生物在橡膠、塑料和化纖等化工產(chǎn)品生產(chǎn)中得到廣泛應用。隨著紡織工業(yè)及石油化工行業(yè)的發(fā)展,芳烴的需求量持續(xù)增加[1-3]。芳烴裝置中的歧化及烷基轉移單元,具有反應器部分溫位高、壓力大,同時精餾塔數(shù)量多、負荷大等特點。因此,對甲苯歧化及烷基轉移單元進行節(jié)能降耗對于整個芳烴裝置的節(jié)能有著重要的意義。

        目前,已經(jīng)有部分過程系統(tǒng)優(yōu)化方法和技術應用在甲苯歧化及烷基轉移單元的工藝節(jié)能改造中。Chen等[4]基于熱力學分析方法,首先運用塔總組合曲線(CGCC)首先對塔系中各單塔的操作參數(shù)進行優(yōu)化,再運用夾點技術分析背景工藝過程物流關系,進而分析了塔系間及其與背景工藝間的耦合關聯(lián)特性,提出了一種新型的基于CGCC的塔系熱量集成與優(yōu)化方法,并且成功運用于甲苯歧化單元和二甲苯單元,經(jīng)濟分析評價表明,改造方案的投資回收期僅0.26 a。Zhang等[5]同步優(yōu)化甲苯歧化及烷基轉移單元分餾塔的進料溫度及其換熱網(wǎng)絡,不僅能降低整個單元的公用工程消耗,還打破了分餾塔與物流間的換熱網(wǎng)絡的邊界問題。

        華賁、陳清林等[6-9]以過程能量結構“三環(huán)節(jié)”模型為基礎,建立了過程系統(tǒng)能量分析與綜合優(yōu)化系統(tǒng)方法,并將其發(fā)展為過程能量結構理論;提出了過程系統(tǒng)能量結構“三環(huán)節(jié)”模型(即按能量功能可將過程系統(tǒng)劃分為能量的工藝利用、能量回收和能量轉換3個環(huán)節(jié)),按能量流認識過程系統(tǒng)的結構,建立了嚴格的定量化的能量分析、分析和經(jīng)濟分析模型,進而將其發(fā)展成為以子系統(tǒng)經(jīng)濟優(yōu)化、分解協(xié)調優(yōu)化和復合措施等為特色的過程能量綜合優(yōu)化的普遍化策略方法,形成了過程系統(tǒng)的能量結構理論[10-13];在過去的幾十年,已開展了一系列將過程能量綜合優(yōu)化方法應用于操作優(yōu)化的研究。Chen等[14]基于“三環(huán)節(jié)”能量流結構模型和“三環(huán)節(jié)”經(jīng)濟模型,對延遲焦化裝置進行了用能分析與改進,與現(xiàn)有裝置相比,實施改進措施后,裝置能耗下降約37.2%。華賁等[15]運用經(jīng)濟學相關理論對換熱過程的熱損失進行了計價研究,對換熱器投資估算方法進行深入分析,發(fā)展了基于換熱器最優(yōu)傳熱溫差的數(shù)學模型,開發(fā)了可應用于工程設計過程的換熱器最優(yōu)傳熱溫差、計算機輔助優(yōu)化計算軟件;基于該方法進行換熱網(wǎng)絡的設計與綜合,可顯著降低過程操作和投資費用,并且證明用原來考慮物料-能量-設備投資三維權衡的思路來指導設備不動的操作調優(yōu)完全可行,并為多家石化企業(yè)進行了操作調優(yōu)工程研究,取得了較好的經(jīng)濟效益。

        過程模擬的目的是通過建立數(shù)學模型并對數(shù)學模型進行求解,研究過程系統(tǒng)的整體穩(wěn)態(tài)特性,從而探索系統(tǒng)和設備改進的可行性,并確定化工系統(tǒng)在已知條件下的預期效果[16]。美國Aspen Tech公司推出的Aspen Plus軟件在過程分析、設計以及參數(shù)優(yōu)化中應用廣泛,已被用于設計與改進工廠和工藝流程,獲得了經(jīng)過驗證的巨大經(jīng)濟效益[17]。筆者首先借助Aspen Plus軟件對甲苯歧化及烷基轉移單元進行流程模擬,然后運用三環(huán)節(jié)能量結構模型,對整個單元進行用能分析及評價;在所揭示的用能瓶頸基礎上,合理設計過程的換熱網(wǎng)絡(HEN),減少或消除過程重復加熱或冷卻,降低過程不合理損,有效實現(xiàn)過程節(jié)能;并基于改造后的換熱網(wǎng)絡,對分餾塔的操作參數(shù),包括進料溫度和位置進行優(yōu)化,從而達到提高整個甲苯歧化及烷基轉移單元的能源利用效率的目的。

        1 芳烴聯(lián)合裝置中歧化及烷基轉移單元工藝流程

        歧化及烷基轉移單元是芳烴聯(lián)合裝置中的重要生產(chǎn)單元。它以芳烴聯(lián)合裝置中較為廉價的甲苯和C9芳烴及少量C10芳烴為主要原料,在歧化催化劑和氫氣存在的條件下,發(fā)生烷基轉移反應,生成苯和主要用于生產(chǎn)對二甲苯的C8芳烴。圖1為簡化的歧化及烷基轉移單元的分離部分的原則工藝流程。本單元包括3座分餾塔、5臺換熱器和1臺加熱爐。反應器出料與循環(huán)氫分離之后,進入換熱器E-11、E-13換熱,然后進入汽提塔C-11。C-11的底部流出物流與芳烴抽提裝置來的芳烴混合,通過換熱器E-14之后送入白土塔D-11,除去少量的烯烴。D-11的流出物進入苯塔C-12,且從第3塊塔板得到產(chǎn)品苯。從苯塔C-12塔底出來的物流進入甲苯塔C-13,甲苯塔C-13的塔頂產(chǎn)品即為甲苯,C-13的底部流出物為二甲苯分餾單元的原料。

        2 歧化及烷基轉移單元三環(huán)節(jié)用能分析

        流程模擬軟件Aspen Plus可以對不同流體的氣-液、液-液平衡常數(shù)以及熵、焓、密度等參數(shù)進行模擬計算。歧化及烷基轉移單元的物系是含石油餾分、芳烴、H2及水的混合物,組分多,物性復雜。

        圖1 芳烴聯(lián)合裝置中甲苯歧化及烷基轉移單元流程

        對于該臨氫物系,可采用Soave-Redich-Kwong (SRK)狀態(tài)方程法和Peng-Robinson(PR)方法進行模擬計算。筆者對汽提塔、分餾塔以及換熱網(wǎng)絡的模擬熱力學方法均采用SRK狀態(tài)方程法,能級ε和物理計算的基準壓力設為0.1 MPa,基準溫度為288 K。

        各分餾塔模擬過程中,根據(jù)進、出料位置劃分不同塔段,采用分段選擇塔板效率的方法進行模擬計算;依據(jù)各分餾塔的實際作用與特點,按其實際分離能力的不同,根據(jù)經(jīng)驗值進行適當選擇。通常,在模擬過程中需要在一定范圍內適當調整塔板效率,直到塔板上各點的溫度、物料流量以及產(chǎn)品餾程等數(shù)據(jù)與采集的實際操作數(shù)據(jù)大致吻合[18]。對歧化反應產(chǎn)物等難以直接測定具體組分的物流,可依據(jù)物料平衡的原則,利用產(chǎn)品物流反推反應物組成的方法,并綜合考慮產(chǎn)品質量、產(chǎn)品分布等實際情況對其進行微調,以得到較準確的模擬所需的物流組成數(shù)據(jù)。

        2.1能量轉換環(huán)節(jié)

        從能量分析結果看出,整個裝置的外供入能達到63.61 MW,能量轉換環(huán)節(jié)用能效率為82.56%;從分析結果可知,能量轉換環(huán)節(jié)的效率為36.76%,能量轉換環(huán)節(jié)的大部分過程損由加熱爐燃料燃燒造成。燃燒過程將高能級的燃料化學能轉化為能級較低的熱能,該過程的不可逆性導致過程損大。

        2.2能量利用環(huán)節(jié)

        2.3能量回收環(huán)節(jié)

        3 基于用能分析結果對歧化及烷基轉移單元進行用能優(yōu)化改進

        根據(jù)對歧化單元原有流程的用能瓶頸分析,提出相應的改進措施,包括:(1)結合歧化原料油換熱流程的優(yōu)化,抽提混芳溫度提升至162℃,避免大溫差物流直接混合導致混合過程損;(2)結合汽提塔進料換熱流程優(yōu)化,控制汽提塔進料與汽提塔底油換熱器的換熱量,使汽提塔底油與抽提混芳混合后的溫度仍能達到177℃,取消白土塔進料加熱;(3)歧化反應產(chǎn)物采用熱高分流程,從源頭減少裝置的低溫余熱量;(4)取消白土塔進料加熱器及白土塔進出料換熱器,使苯塔進料溫度由原來的136℃提高至177℃,降低塔底再沸器加熱負荷。

        3.1能量回收環(huán)節(jié)用能改進

        3.1.1歧化單元原料換熱流程用能優(yōu)化改進

        從圖1可以看到,原流程中將66℃的成品塔頂氣、192℃的重芳烴塔頂甲苯和154℃的循環(huán)甲苯進行混合,混合時溫差較大,造成較大的損。在原流程中增加換熱器,逐步降低反應流出物溫度,可極大減小大溫差換熱過程損,改進后流程如圖2所示。從圖2可以看出,根據(jù)進料和補充進料換熱流程進行局部改進,在循環(huán)甲苯進料處新增換熱器EN-02,與來自二甲苯塔頂產(chǎn)品換熱,將循環(huán)甲苯溫度提高至164.0℃。結合汽提塔進料換熱流程優(yōu)化,通過控制汽提塔進料與汽提塔底油換熱器的換熱量,使汽提塔底油與抽提混芳混合后的溫度仍能達到177℃,取消白土塔進料加熱,在該過程中,苯塔進料溫度由原來的136℃提高至177℃。

        圖2 歧化及烷基轉移單元用能改進后流程簡圖

        采用模擬軟件Aspen Plus對改進后流程進行模擬計算,并基于模擬計算結果,作出改進前后的能級(ε)-熱量(Q)圖,結果示于圖3。圖3中灰色陰影部分的面積代表原有預加氫產(chǎn)物換熱流程的損,紅色陰影部分的面積表示改造后換熱過程損。由圖3可知,流程改進后減少了3.80 MW的蒸汽用量,苯塔熱負荷降低了2.56 MW。大溫差換熱損由現(xiàn)有的7.64 MW降至為5.73 MW,降低了約25.0%。

        圖3 歧化及烷基轉移單元進料預熱流程用能改進前后的ε-Q圖

        3.1.2歧化單元反應產(chǎn)物分離用能優(yōu)化改進

        在一定溫度、壓力下,歧化及烷基轉移反應才能進行。根據(jù)相應的工藝要求,反應后的物料需要進行分離。反應產(chǎn)物約361℃,直接與103℃的混合進料換熱,反應產(chǎn)物的溫度降至107℃,進一步空冷水冷至37℃,而后又需要升溫至一定的溫度進入后續(xù)的汽提塔中。在此換熱網(wǎng)路中,歧化及烷基轉移的反應產(chǎn)物存在重復冷卻、升溫的過程,造成能量的浪費以及傳熱損增大。采用類熱高分工藝來代替現(xiàn)有的冷高分流程,并且以最小過程損為優(yōu)化目標,裝置內合理配匹熱源和熱阱,同時換熱過程中遵循能量梯級利用,降低傳熱損,盡可能多地回收反應產(chǎn)物的熱量,以提高后續(xù)苯塔的進料溫度,從而達到降低塔底再沸負荷的目的,減少裝置能耗。由圖2可知,361℃的反應產(chǎn)物與131℃的反應進料換熱,反應產(chǎn)物的溫度降至137℃,然后直接進入熱分離罐進行分離。對比改進前后的流程可知,改進后的流程不僅避免物流的重復冷卻、升溫過程,同時也將傳熱溫差從24.5℃降至15.5℃,基于模擬計算結果得到的改進前后的ε-Q圖示于圖4。圖4中灰色陰影部分面積代表原有預加氫產(chǎn)物換熱流程的損,紅色陰影部分面積表示改造后換熱過程損。從圖4可知,此過程不僅減小了損約0.95 MW,同時空冷器A-01負荷相應減少至23.14 MW,較之前負荷30.56 MW減小了7.42 MW,下降幅度為24.3%。

        圖4 甲苯歧化及烷基轉移單元反應產(chǎn)物換熱流程用能改進前后的ε-Q圖

        3.2能量利用環(huán)節(jié)用能優(yōu)化與改進

        3.2.1歧化及烷基轉移單元塔系的操作參數(shù)優(yōu)化

        采用流程模擬軟件Aspen Plus對汽提塔、苯塔和甲苯塔分別進行操作參數(shù)的優(yōu)化,得到的汽提塔C-11、苯塔C-12和甲苯塔C-13的塔底熱負荷與進料位置的關系如圖5所示,改造前后的能耗列于表1。由圖5(a)和表1可知,優(yōu)化改進方案中,汽提塔C-11的進料板數(shù)從第18塊塔板下移至第9塊塔板,回流比從3.02降低至2.41;冷凝器的負荷由改造前的6.74 MW降至改造后的5.51 MW,再沸器負荷從改造前的11.33 MW的降至改造后的10.09 MW,塔板數(shù)由改造前的33塊增加至改造后的35塊。圖5(b)表明,苯塔的進料板位置也不合理。優(yōu)化方案將進料板從第15塊調整為第24塊,冷凝器的負荷由改造前的17.72 MW降至改造后的16.59 MW,再沸器負荷從16.85 MW下降至15.58 MW。甲苯塔C-13有49塊塔板,進料板為第20塊塔板,進料溫度為140℃,其回流比為2.50。圖5(c)表明,甲苯塔C-13進料板位置不合理。

        表1 甲苯歧化及烷基轉移單元汽提塔、苯塔和甲苯塔改造前后操作參數(shù)

        圖5 甲苯歧化及烷基轉移單元汽提塔、苯塔和甲苯塔熱負荷(H)與進料位置(N)的關系

        汽提塔C-11和甲苯塔C-13的塔頂溫度分別為135℃ 和180℃,直接冷卻排棄;白土塔D-11的底部流出物預熱D-11的進料物流,之后送至C-12(見圖1)。因此,可以考慮將C-13的塔頂氣相物流首先預熱汽提塔C-11進料物流,而汽提塔C-11的塔頂氣相物流可以用來預熱白土塔D-11的進料。這樣,白土塔的底部流出物流(190℃)能有足夠的熱量將甲苯塔進料從136℃預熱至177℃,相應地進料位置由第20塊塔板上移至第24塊塔板。由表1可知,冷凝器的負荷由改造前的26.85 MW降至改造后的26.43 MW,再沸器負荷則從30.97 MW下降至28.68 MW。

        4 甲苯歧化及烷基轉移單元改進后的裝置用能評價

        基于優(yōu)化改進后流程的模擬結果,對改進后的流程進行三環(huán)節(jié)能量與分析,結果列于表2~4。

        表2 甲苯歧化及烷基轉移單元能量回收環(huán)節(jié)改造前后用能分析結果

        表3 甲苯歧化及烷基轉移單元能量利用環(huán)節(jié)改造前后用能分析結果

        表4 甲苯歧化及烷基轉移單元能量轉換環(huán)節(jié)改造前后用能分析結果

        5 結 論

        運用Aspen Plus模擬甲苯歧化及烷基轉移單元流程,并基于三環(huán)節(jié)能量結構模型,對整個單元進行用能分析及評價。根據(jù)三環(huán)節(jié)用能分析所揭示的用能瓶頸,合理優(yōu)化設計過程的換熱網(wǎng)絡(HEN),減少或消除過程重復加熱或冷卻,降低過程不合理損,有效實現(xiàn)過程節(jié)能,同時對分餾塔的操作參數(shù),包括進料溫度和位置進行優(yōu)化。對能量利用和能量回收環(huán)節(jié)的優(yōu)化和調整后,總供給能量從63.61 MW下降至60.05 MW,整個甲苯歧化及烷基轉移單元的能源利用效率提高約為6%,表明該優(yōu)化改進策略合理。

        [1]邱江.芳徑生產(chǎn)技術現(xiàn)狀及研究進展[J].當代化工,2006,35(5):313-317.(QIU Jiang.Present situation and research progress of aromatic-hydrocarbon production technology[J].Contemporary Chemical Industry,2006,35(5):313-317.)

        [2]韓鳳山,林克芝.世界芳經(jīng)生產(chǎn)技術的發(fā)展趨勢[J].當代石油石化,2006,14(5):30-35.(HAN Fengshan,LIN Kezhi.Development trend of aromatic production and technology in the world[J].Petroleum & Petrochemical Today,2006,14(5):30-35.)

        [3]米多.國內外芳烴生產(chǎn)及消費預測[J].江蘇化工,2004,32(6):59-62.(MI Duo.The forecasting for consumption of domestic and aromatics production[J].Jiangsu Chemical Industry,2004,32(6):59-62.)

        [4]CHEN T,ZHANG B J,CHEN Q L.Heat integration for fractionating systems in para-xylene plant based on column optimization[J].Energy,2014,72(1):311-321.

        [5]ZHANG B,LUO X,LIU K,et al.Simultaneous target of HEN and columns with variable feed temperatures for a toluene disproportionation plant[J].Industrial & Engineering Chemistry Research,2014,53(25):10429-10438.

        [6]華賁.煉油工藝過程用能的理論基礎與實用分析[J].石油煉制,1982,(2):25-31.(HUA Ben.Theoretical basis and practical analysis of energy use in refining processes[J].Petroleum Refining,1982,(2):25-31.)

        [7]HUA B,YIN Q,WU G.Energy optimization through exergy-economic evaluation[J].Trans of the ASME Journal of Energy Resources Technology,1989,111(3):148-153.

        [8]華賁.過程系統(tǒng)的能量綜合和優(yōu)化[J].化工進展,1994,(5):6-14.(HUA Ben.Energy integration and optimization of process systems[J].Chemical Industry and Engineering Progress,1994,(5):6-14.)

        [9]陳清林,尹清華,王松平,等.過程系統(tǒng)能量流結構模型及其應用[J].化工進展,2003,22(3):239-240.(CHEN Qinglin,YIN Qinghua,WANG Songping,et al.Energy structural model of process systems and its application[J].Chemical Industry and Engineering Progress,2003,22(3):239-240.)

        [10]華賁.過程能量綜合的研究進展與展望[J].石油化工,1996,25(1):62-70.(HUA Ben.Integration of energy process progress and prospects petrochemical technology[J].Petrochemical Technology 1996,25(1):62-70.)

        [12]華賁.工藝過程用能分析與綜合[M].北京:烴加工出版社,1989.

        [13]HUA B,WU G,CHEN Z,et al.Exergy-economic optimization for energy systems in chemical process[J].Chemical Engineering (China),1991,19(2):64-75.

        [14]CHEN Q L,WANG S P,YIN Q H,et al.Energy-use analysis and improvement of delayed coking units[J].Energy,2004,29(12):2225-2237.

        [16]張程.基于ASPEN PLUS的常減壓裝置過程模擬與換熱網(wǎng)絡優(yōu)化[D].北京:中國石油大學,2011.

        [17]吳瑤.基于Aspen Plus對耦合化學鏈燃燒的焦爐煤氣重整制氫系統(tǒng)性能的研究[D].重慶:重慶大學,2014.

        [18]王洪江.流程模擬計算中關于塔板效率選擇的技巧[J].石化技術,2003,10(1):31-33.(WANG Hongjiang.The choice of column stage efficiencies in simulation[J].Petrochemical Industry Technology,2003,10(1):31-33.)

        Energy-Use Analysis and Optimization of a Toluene Disproportionation andTransalkylation Unit

        CHEN Ting1,2,CHEN Qinglin3

        (1.Wuhan Research Institute,Jianghan University,Wuhan 430056,China; 2.School of Chemistry and Environmental Engineering,Jianghan University,Wuhan 430056,China; 3.Key Laboratory of Low-Carbon Chemistry & Energy Conservation of Guangdong Province, School of Chemistry and Chemical Engineering,Sun Yat-Sen University,Guangzhou 510275,China)

        Based on the simulation of the process of a disproportionation unit by the process simulation software Aspen Plus,energy-use analysis and optimization of the toluene disproportionation and transalkylation unit were performed by a three-link energy structural model,from which it is found poor design of the heat exchanger network (HEN)and inappropriate fractionating towers operation were the bottlenecks of energy-use.To relieve the drawbacks,new heat exchange networks were proposed,and the feed temperature and the feed location of the fractionating towers were optimized.In the optimized process the energy requirement was decreased by about 6% as compared to that of the original unit,which suggested that the proposed strategies were beneficial to the design and operation of the continuous catalytic reforming unit and disproportionation unit.

        three-link energy structural model; the toluene disproportionation and transalkylation unit; process simulation; energy-use optimization

        2015-10-19

        江漢大學武漢研究院開放性課題(jhunwyy2015324)資助

        陳婷,女,講師,博士研究生,從事過程系統(tǒng)工程方面的研究

        陳清林,男,教授,從事過程系統(tǒng)工程方面的研究;E-mail:chqlin@mail.sysu.edu.cn

        1001-8719(2016)05-0959-08

        TE624

        Adoi:10.3969/j.issn.1001-8719.2016.05.013

        猜你喜歡
        優(yōu)化
        超限高層建筑結構設計與優(yōu)化思考
        PEMFC流道的多目標優(yōu)化
        能源工程(2022年1期)2022-03-29 01:06:28
        民用建筑防煙排煙設計優(yōu)化探討
        關于優(yōu)化消防安全告知承諾的一些思考
        一道優(yōu)化題的幾何解法
        由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
        圍繞“地、業(yè)、人”優(yōu)化產(chǎn)業(yè)扶貧
        事業(yè)單位中固定資產(chǎn)會計處理的優(yōu)化
        消費導刊(2018年8期)2018-05-25 13:20:08
        4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
        幾種常見的負載均衡算法的優(yōu)化
        電子制作(2017年20期)2017-04-26 06:57:45
        精品亚亚洲成av人片在线观看| 久久免费的精品国产v∧| 无码专区中文字幕DVD| 日本高清一区二区三区视频| 国产一区二区精品人妖系列在线 | 久久综合给合久久97色| 久久伊人精品中文字幕有尤物 | 亚洲日韩国产一区二区三区在线| 正在播放国产多p交换视频| av中文字幕在线资源网| 亚洲国产性夜夜综合另类| 久久精品黄aa片一区二区三区| 99国产精品视频无码免费| 国产大片在线观看三级| 伊人久久精品亚洲午夜| 中文字字幕在线精品乱码| 亚洲中文欧美日韩在线人| 亚洲av中文字字幕乱码软件 | 亚洲美女主播内射在线| 亚洲国产精品无码久久久| 内射交换多p国产| 日韩精品一区二区亚洲av性色| 久久国产成人午夜av免费影院| 日本黄网站三级三级三级| 国产中文制服丝袜另类| 激情视频国产在线观看| 亚洲av综合色区无码一区| 亚洲精华国产精华液的福利 | 国产av无码专区亚洲av男同 | 久久精品国产亚洲vr| 国产一区二区三区porn| 日本丰满少妇xxxx| 最近中文字幕视频高清| 成在线人免费无码高潮喷水| 国产精品福利高清在线| 幻女bbwxxxx在线视频| 精品人妻丰满久久久a| 国产精品亚洲一区二区三区在线 | 一区二区精品国产亚洲 | 国模精品一区二区三区| 国产午夜无码视频免费网站|