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

    基于改進(jìn)恒熱傳輸模型的精餾模擬初始化

    2023-10-14 07:53:10張帆陶少輝陳玉石項曙光
    化工進(jìn)展 2023年9期
    關(guān)鍵詞:塔板沸物精餾塔

    張帆,陶少輝,陳玉石,項曙光

    (1 青島科技大學(xué)化工學(xué)院過程系統(tǒng)工程研究所,山東 青島 266042;2 石化盈科信息技術(shù)有限責(zé)任公司,上海 200040)

    精餾是化工分離過程中的重要操作單元之一,在化工生產(chǎn)及應(yīng)用方面占有舉足輕重的地位[1-4]。精餾過程模擬是通過給定的進(jìn)料條件以及塔模擬設(shè)置迭代求解非線性MESH 方程組,從而得出塔板上溫度、壓力、流量以及組成分布的過程[5-10]。而精餾模擬的初始化是指為MESH 方程組求解提供一組可行、高效的初始解,從而加速精餾模擬的收斂,或者促使本來難收斂的精餾模擬收斂的過程。

    精餾模擬計算中應(yīng)用的初始化算法中,最常用的是Boston 和Sullivan[11]基于內(nèi)外層計算所提出的方法,雖然該方法給出的初值精度并不是很高,但是對于大多數(shù)分離體系仍能收斂;針對涉及共沸物的強(qiáng)非理想體系精餾過程,Venkataraman和Lucia[12]提出了考慮共沸物的初始化方法;另外,Ivakpou和Kasiri[13]提出了基于改進(jìn)Kremser 法[14]的初始化方法,并加入Broyden 阻尼因子[15-16]用以加速迭代計算的收斂;而針對原油精餾過程,Rabeau等[17]提出了一種減少組分?jǐn)?shù)和不考慮組分間交互作用的特殊初始化方法。此外,F(xiàn)letcher 和Morton[18]基于精餾塔全回流模擬、Barttfeld 與Aguirre[19]針對基于可逆精餾模型分別提出了不同的初始化方法。

    現(xiàn)有的初始化方法中,多采用恒摩爾流(constant molar overflow,CMO)假設(shè)[20-22]來進(jìn)行流量初值的估算。然而恒摩爾流假設(shè)成立的前提是,精餾過程中各組分間的蒸發(fā)焓相等,且其他熱效應(yīng)可忽略不計或抵消,但是忽略組分間的蒸發(fā)焓的差異必然導(dǎo)致塔板間流量估值的較大誤差。本文針對Rev[23]提出的恒熱傳輸(constant heat transport,CHT)模型進(jìn)行改進(jìn),在研究組分間蒸發(fā)焓差異對塔板流量的影響的基礎(chǔ)上,考慮了溫度對組分蒸發(fā)焓的影響,并針對具有多進(jìn)料或多側(cè)線采出復(fù)雜精餾過程,推導(dǎo)出了塔板汽液相流量分布的簡捷計算公式,從而提出了一種新的精餾模擬初始化方法。并通過烴類和強(qiáng)非理想體系的兩個分離實例研究了各種初始解的準(zhǔn)確程度,以及其代入模擬計算后的收斂情況。

    1 恒熱傳輸模型及其改進(jìn)

    1.1 Rev恒熱傳輸模型

    Rev認(rèn)為,由于塔板間物流顯熱的變化遠(yuǎn)小于潛熱,故在進(jìn)行熱量衡算時可以忽略物流顯熱變化,只考慮物流潛熱對精餾過程的影響[23]。假設(shè)精餾過程中塔板上汽液相焓流保持不變,汽液相流量的改變量除去進(jìn)料與采出的影響外,只取決于組分蒸發(fā)焓變化。據(jù)此,可對如圖1所示的簡單精餾塔的精餾段和提餾段分別進(jìn)行如下分析。

    圖1 簡捷精餾模型

    對于如圖1(a)中的精餾段,從塔頂?shù)降趈塊板做物料衡算,并將塔頂組分的蒸發(fā)焓ΔvapHi,d(J/mol)代入兩側(cè)獲得式(1)。

    式中,Vj和Lj分別代表第j塊塔板的汽相、液相摩爾流量,mol/s;yi,j和xi,j為第j塊塔板組分i的汽相、液相摩爾分率;D與xi,d為塔頂采出摩爾流量以及組分i的塔頂采出摩爾流量,mol/s。

    對所有組分加和式(1),并假設(shè)精餾段各塔板上組分的蒸發(fā)焓ΔvapHi,j均等于塔頂組分的蒸發(fā)焓ΔvapHi,d,且不考慮塔板間顯熱的變化,則有式(2)成立。

    至此,可以得出對于第j塊塔板上的氣液相摩爾流量Vj和Lj,如式(3)和式(4)所示。

    式中,R為回流比;L0為塔頂回流量,mol/s。

    同理,從塔底到第j塊板做物料衡算,將塔頂組分的蒸發(fā)焓ΔvapHi,b代入衡算式并對所有組分加和,并假設(shè)提餾段各塔板上組分的蒸發(fā)焓ΔvapHi,j均等于塔頂組分的蒸發(fā)焓ΔvapHi,b,則對于提餾段塔板氣液相流量,有式(5)和式(6)成立。

    式中,W、xi,b、B和QR分別為塔底采出摩爾流量(mol/s)、組分i的塔底采出組成(mol/s)、再沸比以及再沸器熱負(fù)荷(kW)。

    Rev[23]將該方法應(yīng)用于三元混合物的熱偶精餾計算中,結(jié)果表明CHT 模型可以準(zhǔn)確預(yù)測非關(guān)鍵組分在塔頂塔底中的分布情況,以及熱偶精餾過程中的極限流量。

    1.2 改進(jìn)恒熱傳輸模型

    雖然CHT 模型可以較為準(zhǔn)確地描述塔板流量的變化,但是仍存在如下問題:首先,由于組分的蒸發(fā)焓受溫度的影響并非總是可以忽略,尤其是針對溫度變化較大的精餾,忽略組分蒸發(fā)焓隨溫度的變化將造成塔板間流量計算結(jié)果與實際值偏差過大;其次,上述CHT 方法僅適用于簡單精餾塔,對于多股進(jìn)料或采出的復(fù)雜精餾塔并未考慮。因此,本文考慮了組分蒸發(fā)焓與溫度之間的關(guān)系,并通過物料和能量守恒關(guān)系將CHT 模型適用范圍推廣至圖2(a)所示的復(fù)雜精餾塔。

    圖2 復(fù)雜精餾塔模型

    如圖2(b)所示,從塔頂?shù)降趈塊塔板對該精餾過程做物料以及能量衡算,可得式(7)和式(8)。

    式中,Gj和Uj分別為第j塊塔板汽、液相采出摩爾流量,mol/s;QF,j為第j塊塔板上進(jìn)料所攜帶的熱量,kW;Hj和hj分別為氣相、液相物流的摩爾焓,J/mol。整理式(7)可得第j塊塔板上氣相流量的計算式(9)。

    式中,F(xiàn)k為進(jìn)料摩爾流量,mol/s。

    將式(9)代入式(8)中,可得到第j塊塔板上液相流量的計算式(10)。

    由相同組成的氣相、液相混合物的摩爾焓間關(guān)系,Hj(x)=hj(x)ΔvapHj(x),其中ΔvapHj(x)為各塔板混合物蒸發(fā)焓,又因為塔板間的顯熱變化相對于潛熱可以忽略不記,即在式(10)中,假設(shè)式(11)和式(12)成立。

    式中,下角標(biāo)k為大于0小于j的任意塔板,故對于第j塊塔板的液相流量計算如式(13)所示。

    式中,各塔板蒸發(fā)焓ΔvapHj(x)可根據(jù)塔板溫度和組成進(jìn)行計算。

    最后,將通過式(13)計算得到的液相流量Lj代入式(9)中,即可得到塔板上汽相流量分布。

    由上述可知,在已知塔頂冷凝器熱負(fù)荷的前提下,只需給出塔板溫度以及液相組成,通過式(9)和式(13)即可給出塔板上汽液相流量的分布。

    2 基于改進(jìn)CHT模型的初始化方法

    對于一個塔板數(shù)為N,進(jìn)料組分為M的精餾過程,在已知所有進(jìn)料以及模擬設(shè)置后,通過如下步驟完成初始化賦值,其過程如圖3所示。

    圖3 基于改進(jìn)CHT算法的初始化算法框圖

    2.1 估算各采出塔板的組成(包括側(cè)流)和溫度

    對于一個無共沸物的精餾塔,一般沸點(diǎn)越低的組分其采出位置越靠近塔頂,而隨著沸點(diǎn)增加,其餾出位置也越靠下。對于整個精餾過程,混合所有進(jìn)料,并對所有進(jìn)料組分按照沸點(diǎn)大小進(jìn)行排序,并由小往大依次分配給塔頂?shù)剿椎拿恳粔K采出塔板,即可得到每塊采出塔板的組分流量以及組成的初始估計值。

    如表1所示,對應(yīng)不同的采出情況,將采出板組成賦值各自對應(yīng)塔板上。其中zi,j代表第j塊塔板上有采出時組分i的采出組成。當(dāng)采出相態(tài)為液相時,采出板對應(yīng)塔板的液相組成為采出組成,該組成的泡點(diǎn)溫度為塔板溫度估計值,同時泡點(diǎn)計算求得的汽相組成為塔板汽相組成。同理,當(dāng)采出板相態(tài)為汽相時,則采出板對應(yīng)塔板的汽相組成為采出組成,作露點(diǎn)計算即可得出塔板溫度以及對應(yīng)的液相組成。若塔頂同時采出氣液兩相,采出組成為汽相組成,并在給定的汽相分率或溫度下進(jìn)行閃蒸,即可得出塔頂液相組成。

    表1 不同采出位置所對應(yīng)塔板的組成

    而對于存在共沸物的強(qiáng)非理想體系,當(dāng)壓力和組分一定時,通過共沸物搜尋,可以得出所有進(jìn)料組分中所有可能存在的共沸物組成以及其對應(yīng)的共沸點(diǎn)。Felbab[24]以及Maier等[25]均對混合組分共沸搜尋方法進(jìn)行了研究,結(jié)果表明混合進(jìn)料中組分的所有可能共沸物不難獲得,并且其計算時間基本保持在0.01~0.02s之間。故在平均塔壓下,找出進(jìn)料組分中含有的所有的共沸物以及對應(yīng)組成以及共沸點(diǎn),將每一種共沸物看作是一個虛擬組分。對所有的純組分以及虛擬組分按照沸點(diǎn)(共沸物為其共沸點(diǎn))排序并按溫度大小,依次從塔頂?shù)剿追峙浣o每塊采出塔板。另外,在分配完成后還需重新按每個共沸物的組成計算各自的純組分含量,并去除虛擬組分所代表的共沸物,按純組分重新計算各采出塔板組分的含量及組成。

    2.2 初始組成以及溫度分布

    假設(shè)塔板上氣液相組分組成呈線性變化,并忽略進(jìn)料位置處組分流量的突變。則可由各采出板將精餾塔分為若干段。對每段上分別線性插值,則可給出所有塔板上組分的氣液相組成以及溫度,如式(14)~式(22)所示。

    從塔頂?shù)降谝粔K液相采出塔板位置k,第j塊塔板上的溫度Tj(K)以及對于組分i的汽液相組成yi,j和xi,j,如式(14)~式(16)。

    2.3 基于改進(jìn)的CHT模型計算塔板上流量

    已知塔板上氣液相組成、進(jìn)料條件以及精餾塔模擬設(shè)置,對整個精餾塔作物料衡算,可得到塔頂、塔底采出量。并通過上述條件進(jìn)一步計算出第一塊塔板上的氣液相流量以及塔頂冷凝器熱負(fù)荷QC,如式(23)和式(24)所示。

    全凝器

    部分冷凝器

    然后結(jié)合式(9)和式(13)即可完成對塔板上流量的賦值。

    3 實例驗證

    為評價上述初始化方法,本文分別使用烴類體系以及一個包含共沸物的強(qiáng)非理想體系分離過程作為測試案例,以對比Aspen plus 中RadFrac 模塊自帶初始化算法(即Boston 法)、Venkataraman 法、基于原CHT 模型的初始化方法和本文提出的基于改進(jìn)CHT 模型的初始化方法。分別對照了各方法估算出的溫度和氣液相流量與模擬結(jié)果的相對偏差和各自所耗時間,并將各初始化結(jié)果代入到Aspen plus中進(jìn)行模擬,以對比加入各初值估計后精餾過程的收斂情況。

    3.1 烴類體系分離過程中初始化方法的應(yīng)用與對比

    在本算例中,如圖4所示的為某化工廠氣分裝置中C3、C4 分離塔,C3 和C4 的混合物經(jīng)過精餾塔T1-1 分離后,在塔頂全凝采出大部分C3 組分,而塔底則為C4 組分,其進(jìn)料條件以及精餾塔模擬設(shè)置見表2。選取物性方法PR,分別通過上述3種初始化對該分離過程計算。

    表2 C3、C4分離過程的進(jìn)料情況

    圖4 C3、C4分離塔

    圖5 為3 種初始化結(jié)果與模擬結(jié)果的相對誤差對比。由圖5(a)可知,所提出方法給出的溫度分布較Boston 法和Venkataraman 法更接近最終模擬結(jié)果,而基于原始CHT 模型的初始化方法與基于改進(jìn)CHT方法溫度估計一致。由圖5(b)、(c)可知,基于改進(jìn)CHT 法估算的塔板間汽液相流量更接近模擬結(jié)果,Boston法和Venkataraman法皆采用恒摩爾流假設(shè)進(jìn)行流量估計,所以兩者的流量誤差曲線完全重合;比較原CHT 法和改進(jìn)CHT 法估算的塔板間汽液相流量,原CHT 法中取各塔板上純組分蒸發(fā)焓均等于塔頂組分蒸發(fā)焓,故在越靠近塔底位置其流量估計結(jié)果與實際相差越大,甚至當(dāng)塔板上溫度與塔頂相差過大時,其結(jié)果甚至不如恒摩爾流假設(shè)計算結(jié)果。表3中列出了各自算法的初始化耗時以及代入模擬計算的迭代次數(shù),雖然改進(jìn)CHT 法較前三者耗時稍多,但能夠有效地減少整個模擬計算的迭代次數(shù)。最后,由于該流程非理想性較弱,并不能完全體現(xiàn)新初始化方法的優(yōu)勢。

    表3 C3、C4分離過程中各方法計算對比

    圖5 C3、C4精餾塔初始化結(jié)果分布

    3.2 強(qiáng)非理想體系分離過程中初始化方法的應(yīng)用與對比

    根據(jù)上文3.1節(jié)中的算例,改進(jìn)CHT法計算結(jié)果雖然較其余初始化方法計算結(jié)果有一定改進(jìn),但是在最終模擬計算過程中,新初始化優(yōu)勢并不是很明顯。故在本算例中采用強(qiáng)非理想體系分離過程對三者的結(jié)果進(jìn)行對比。圖6所示為國內(nèi)某化工廠的苯酚丙酮分離裝置中的苯酚、丙酮、水以及部分雜質(zhì)的分離過程,其進(jìn)料條件及精餾模擬設(shè)置見表4,選擇NRTL 活度系數(shù)模型進(jìn)行模擬計算。進(jìn)料為苯酚、丙酮、水以及部分雜質(zhì)的F2-1 經(jīng)過塔T2-1 分離后,從塔頂汽相采出的丙酮和水及其相關(guān)共沸物與進(jìn)料為苯酚、丙酮、水的F2-2 混合形成塔T2-2的進(jìn)料F2-3,在塔T2-2中分離,于T2-2 中側(cè)線采出大部分丙酮,而苯酚和水在塔底采出。由于整個分離過程存在多種共沸物,其模擬計算過程很難收斂,甚至對于塔T2-1 更需要用戶添加合理的估值才能促使流程收斂。通過上述分離過程分別對4種初始化方法進(jìn)行驗證。

    表4 強(qiáng)非理想體系的進(jìn)料情況

    圖6 強(qiáng)非理想體系分離過程

    圖7、圖8 及表5、表6 為分別通過4 種初始化對該流程模擬計算的結(jié)果。同3.1 節(jié)中相同,改進(jìn)CHT 法給出的溫度和流量分布更加貼合于最終的模擬結(jié)果。將各初始化計算結(jié)果分別代入到模擬計算過程中,新初始化方法均能有效地減少計算迭代次數(shù),對于T2-2,使用其他初始化方法不能夠保證流程收斂時原CHT 法和改進(jìn)CHT 法仍然能夠保證計算過程的收斂,并且改進(jìn)CHT 法較原CHT 法能夠有效地促進(jìn)流程收斂。

    表5 T2-1塔中各方法計算對比

    表6 T2-2塔中各方法計算對比

    圖7 T2-1塔初始化結(jié)果分布

    圖8 T2-2塔初始化結(jié)果分布

    4 結(jié)論

    (1)本文通過考慮溫度對組分蒸發(fā)焓的影響,改進(jìn)了精餾模擬的恒熱傳輸模型,將其拓展到復(fù)雜精餾塔,并基此提出了一種新的精餾模擬初始化方法。

    (2)相較于恒摩爾流假設(shè),恒熱傳輸模型在計算過程中考慮了組分間潛熱不同對塔板流量的影響,因此能夠更加準(zhǔn)確地估算塔板上汽液相流量的分布,而改進(jìn)恒熱傳輸模型相較于原模型進(jìn)一步考慮了溫度對于組分蒸發(fā)焓的影響,因此在進(jìn)行流量估值時更加接近于模擬結(jié)果,更加有利于精餾嚴(yán)格計算的迭代收斂。

    (3)兩個精餾過程模擬實例表明,相較于其他初始化方法,所提出的初始化方法能夠更加準(zhǔn)確地估算塔板上溫度、流量以及組成的分布。

    (4)所提出的初始化與其他初始化方法在計算耗時方面相差不大,但能夠有效地減少整個精餾模擬的迭代次數(shù),甚至對于一些難收斂體系,該方法還可以促進(jìn)其計算過程的收斂,增加精餾計算的成功率。

    猜你喜歡
    塔板沸物精餾塔
    國產(chǎn)太陽能光熱發(fā)電有機(jī)熱載體再生工藝研究
    一起精餾塔自燃事故的原因分析及防范措施的改進(jìn)
    氯乙烯精餾高沸物提純系統(tǒng)的改進(jìn)
    聚氯乙烯(2018年2期)2018-07-05 09:49:06
    催化裂解法處理有機(jī)硅高沸物
    固旋閥塔板的流體力學(xué)性能研究及其旋轉(zhuǎn)流場CFD模擬
    石油化工(2016年9期)2017-01-21 02:19:58
    塔器的性能研究及工程應(yīng)用分析
    廣州化工(2016年19期)2016-11-23 05:42:12
    氨水吸收式制冷系統(tǒng)中精餾塔性能模擬與分析
    內(nèi)部熱集成精餾塔分離混合碳五的模擬研究
    聯(lián)合精餾過程分離四氫呋喃-水共沸物的實驗及優(yōu)化模擬
    齒邊導(dǎo)向浮閥塔板流體力學(xué)性能的研究及其工業(yè)應(yīng)用
    石油化工(2015年9期)2015-02-05 03:20:36
    午夜福利视频1000在线观看| 欧美性猛交╳xxx乱大交人| 我要看日韩黄色一级片| 毛片一级片免费看久久久久| 国产精品一区www在线观看| 中文字幕制服av| 欧美日本视频| 永久网站在线| 亚洲综合色惰| 黄色视频,在线免费观看| 能在线免费看毛片的网站| 桃色一区二区三区在线观看| 十八禁国产超污无遮挡网站| 岛国毛片在线播放| 久久99蜜桃精品久久| 3wmmmm亚洲av在线观看| 激情 狠狠 欧美| 国产淫片久久久久久久久| 国产精品久久电影中文字幕| 国产伦理片在线播放av一区 | 最新中文字幕久久久久| 亚洲精品粉嫩美女一区| 一级毛片我不卡| 亚洲久久久久久中文字幕| 在线国产一区二区在线| 国产一级毛片七仙女欲春2| 少妇丰满av| 观看美女的网站| 久久精品夜色国产| 午夜视频国产福利| 男女啪啪激烈高潮av片| 国产一区亚洲一区在线观看| 久久6这里有精品| 成人鲁丝片一二三区免费| 中文字幕久久专区| 1000部很黄的大片| 国产亚洲91精品色在线| 免费观看的影片在线观看| 亚洲欧美精品综合久久99| 亚洲人成网站在线观看播放| 99久久久亚洲精品蜜臀av| 日韩制服骚丝袜av| av卡一久久| 中文字幕免费在线视频6| 日日摸夜夜添夜夜添av毛片| 日韩一本色道免费dvd| 色吧在线观看| 日本一二三区视频观看| 国产乱人偷精品视频| 不卡视频在线观看欧美| 自拍偷自拍亚洲精品老妇| 欧美一区二区国产精品久久精品| 久久午夜亚洲精品久久| 中文在线观看免费www的网站| 国内精品一区二区在线观看| 欧美日韩在线观看h| 亚洲内射少妇av| 欧美色视频一区免费| 精品久久久久久久久久久久久| 亚洲在久久综合| 国产精品久久久久久亚洲av鲁大| 有码 亚洲区| 99在线视频只有这里精品首页| 国产综合懂色| 永久网站在线| 久久精品综合一区二区三区| 给我免费播放毛片高清在线观看| 精品久久久久久久久久久久久| 大香蕉久久网| 久久这里只有精品中国| 欧美zozozo另类| 人妻久久中文字幕网| 中国美白少妇内射xxxbb| 男人舔奶头视频| a级毛片免费高清观看在线播放| 日本一二三区视频观看| 国产在线精品亚洲第一网站| 久久精品91蜜桃| 日韩欧美三级三区| 男人的好看免费观看在线视频| 欧美高清性xxxxhd video| 亚洲在久久综合| 黄片wwwwww| 欧美zozozo另类| 精品99又大又爽又粗少妇毛片| 亚洲成av人片在线播放无| 日韩欧美 国产精品| 黑人高潮一二区| 亚洲四区av| 国产在线精品亚洲第一网站| 色哟哟哟哟哟哟| 长腿黑丝高跟| 在线免费观看不下载黄p国产| 成人综合一区亚洲| 男女啪啪激烈高潮av片| 日韩强制内射视频| 国产单亲对白刺激| 看十八女毛片水多多多| 亚洲欧洲国产日韩| 成人无遮挡网站| 乱码一卡2卡4卡精品| 黄色日韩在线| 国产精品一及| 午夜福利在线观看吧| 中文字幕免费在线视频6| www.av在线官网国产| 欧美日韩在线观看h| 少妇裸体淫交视频免费看高清| 99热6这里只有精品| 少妇的逼水好多| 国产精品福利在线免费观看| а√天堂www在线а√下载| 色视频www国产| 黄色一级大片看看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最近最新中文字幕大全电影3| 在线观看免费视频日本深夜| av免费在线看不卡| 美女cb高潮喷水在线观看| 高清在线视频一区二区三区 | 亚洲一区二区三区色噜噜| 爱豆传媒免费全集在线观看| 在现免费观看毛片| 亚洲人成网站在线观看播放| 久久人人精品亚洲av| 免费观看人在逋| 国产精品久久久久久久久免| 国产毛片a区久久久久| 国产一区二区在线观看日韩| 国产精品久久视频播放| 国产综合懂色| 日产精品乱码卡一卡2卡三| 深爱激情五月婷婷| 九九爱精品视频在线观看| 国产精品.久久久| 亚洲综合色惰| 少妇人妻精品综合一区二区 | 99热只有精品国产| 少妇人妻精品综合一区二区 | 久久久久久伊人网av| 在线观看午夜福利视频| 国产精品麻豆人妻色哟哟久久 | 三级经典国产精品| 久久久久久久久久黄片| 亚洲欧美中文字幕日韩二区| 国产成人午夜福利电影在线观看| 日本与韩国留学比较| 我要搜黄色片| 午夜福利成人在线免费观看| av福利片在线观看| 亚洲无线在线观看| 精品99又大又爽又粗少妇毛片| 日本熟妇午夜| 日韩欧美三级三区| 成年免费大片在线观看| 国产精品一区二区三区四区免费观看| 亚州av有码| 国产美女午夜福利| 黄色一级大片看看| 国产高清视频在线观看网站| 久久精品夜色国产| 黄色一级大片看看| 九色成人免费人妻av| 日韩 亚洲 欧美在线| 麻豆一二三区av精品| 此物有八面人人有两片| 欧美一区二区国产精品久久精品| 国产69精品久久久久777片| 国产午夜精品论理片| www日本黄色视频网| 丝袜美腿在线中文| 国产精品野战在线观看| 村上凉子中文字幕在线| 欧美成人a在线观看| 天堂影院成人在线观看| 亚洲成人久久性| 一级二级三级毛片免费看| 又爽又黄无遮挡网站| 成年版毛片免费区| 精品久久国产蜜桃| 成熟少妇高潮喷水视频| 在线观看免费视频日本深夜| 国产一区二区亚洲精品在线观看| 精品无人区乱码1区二区| 成人永久免费在线观看视频| 国产高清激情床上av| 禁无遮挡网站| 中文资源天堂在线| 国产私拍福利视频在线观看| 非洲黑人性xxxx精品又粗又长| 国产伦在线观看视频一区| 精品久久久久久久久亚洲| 天天躁日日操中文字幕| 久久久精品大字幕| 免费看a级黄色片| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| ponron亚洲| 精品免费久久久久久久清纯| 一本精品99久久精品77| 91午夜精品亚洲一区二区三区| 欧美一区二区亚洲| 国产高清三级在线| eeuss影院久久| 亚洲四区av| 搞女人的毛片| 亚洲精品粉嫩美女一区| 成人亚洲精品av一区二区| 一级黄色大片毛片| 久久综合国产亚洲精品| 国产高清三级在线| 国产不卡一卡二| 国产精品一区www在线观看| 久久久久久久久久成人| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 亚洲国产精品合色在线| av又黄又爽大尺度在线免费看 | 美女高潮的动态| av.在线天堂| 日韩强制内射视频| 亚洲人成网站高清观看| 日韩av不卡免费在线播放| 搡老妇女老女人老熟妇| 国产在线精品亚洲第一网站| 亚洲国产精品sss在线观看| 最后的刺客免费高清国语| a级一级毛片免费在线观看| 国产一区二区三区av在线 | 国产免费一级a男人的天堂| 久久久久久久午夜电影| 国产精品蜜桃在线观看 | 日日摸夜夜添夜夜爱| 99久久人妻综合| 蜜臀久久99精品久久宅男| 女人十人毛片免费观看3o分钟| 久久久久久久久久久丰满| 国产午夜福利久久久久久| a级毛片a级免费在线| 欧美3d第一页| 寂寞人妻少妇视频99o| 我要看日韩黄色一级片| 真实男女啪啪啪动态图| 美女大奶头视频| 午夜精品在线福利| 伦理电影大哥的女人| 国产精品一区二区三区四区久久| 青春草国产在线视频 | av在线播放精品| 联通29元200g的流量卡| 国产一区二区在线观看日韩| 国产精品一区www在线观看| 黄色配什么色好看| 亚洲欧洲日产国产| 亚洲在线自拍视频| 草草在线视频免费看| 亚洲人与动物交配视频| 看十八女毛片水多多多| 久久久精品欧美日韩精品| 亚洲不卡免费看| 亚洲激情五月婷婷啪啪| 在线观看午夜福利视频| 国产伦一二天堂av在线观看| 国产一区亚洲一区在线观看| 亚洲av电影不卡..在线观看| 国产精品一区二区性色av| 成人三级黄色视频| 国产美女午夜福利| 国产精品久久久久久久电影| 少妇裸体淫交视频免费看高清| 国产在线精品亚洲第一网站| 国产不卡一卡二| 麻豆精品久久久久久蜜桃| 国产精品国产高清国产av| 亚洲精品色激情综合| 男人狂女人下面高潮的视频| 国产成人精品久久久久久| 欧美一区二区亚洲| 亚洲欧美成人精品一区二区| 波多野结衣高清作品| 欧美3d第一页| 午夜福利在线在线| 欧美xxxx性猛交bbbb| 欧美激情久久久久久爽电影| 国产精品日韩av在线免费观看| 成熟少妇高潮喷水视频| 99国产极品粉嫩在线观看| 亚洲精品影视一区二区三区av| 美女被艹到高潮喷水动态| 久久亚洲精品不卡| 婷婷色av中文字幕| 日本与韩国留学比较| 一本一本综合久久| 白带黄色成豆腐渣| 日日摸夜夜添夜夜爱| 赤兔流量卡办理| 国产精品免费一区二区三区在线| 一卡2卡三卡四卡精品乱码亚洲| 国内少妇人妻偷人精品xxx网站| 成人午夜高清在线视频| 99热这里只有是精品50| 亚洲真实伦在线观看| 亚洲国产精品成人综合色| 亚洲三级黄色毛片| 97超视频在线观看视频| 亚洲国产精品久久男人天堂| 最好的美女福利视频网| 99在线视频只有这里精品首页| 国产精品.久久久| 久久综合国产亚洲精品| 国产探花在线观看一区二区| 久久精品久久久久久噜噜老黄 | av在线老鸭窝| 男女做爰动态图高潮gif福利片| 久久国内精品自在自线图片| 日韩,欧美,国产一区二区三区 | 国产精品日韩av在线免费观看| 国产精品永久免费网站| av在线观看视频网站免费| 一级毛片我不卡| 伦理电影大哥的女人| 久久久精品欧美日韩精品| 亚洲在线观看片| 亚洲美女视频黄频| 直男gayav资源| 97人妻精品一区二区三区麻豆| 国产老妇伦熟女老妇高清| 国产伦一二天堂av在线观看| 国产高清激情床上av| 国模一区二区三区四区视频| 精品人妻一区二区三区麻豆| 日韩,欧美,国产一区二区三区 | 蜜臀久久99精品久久宅男| a级毛片免费高清观看在线播放| 性色avwww在线观看| 别揉我奶头 嗯啊视频| 一进一出抽搐gif免费好疼| 舔av片在线| 韩国av在线不卡| 久久久久国产网址| 国产一区亚洲一区在线观看| 日本一二三区视频观看| 99热这里只有是精品50| 在线a可以看的网站| 爱豆传媒免费全集在线观看| 亚洲va在线va天堂va国产| 内地一区二区视频在线| 亚洲在久久综合| 丰满人妻一区二区三区视频av| 一个人免费在线观看电影| 丝袜美腿在线中文| 女人被狂操c到高潮| 小蜜桃在线观看免费完整版高清| 男人舔奶头视频| 嫩草影院入口| 中文字幕av在线有码专区| 亚洲精品亚洲一区二区| 色播亚洲综合网| 天堂av国产一区二区熟女人妻| 日产精品乱码卡一卡2卡三| 丝袜美腿在线中文| 亚洲精华国产精华液的使用体验 | 亚洲av熟女| 男女那种视频在线观看| 不卡一级毛片| 乱人视频在线观看| 亚洲不卡免费看| av在线蜜桃| 嫩草影院精品99| 狂野欧美白嫩少妇大欣赏| 亚洲国产日韩欧美精品在线观看| 美女大奶头视频| 久久国内精品自在自线图片| 久久久欧美国产精品| 成人国产麻豆网| 午夜亚洲福利在线播放| 日日撸夜夜添| 久久久久九九精品影院| 国产高清三级在线| 少妇被粗大猛烈的视频| 色播亚洲综合网| 午夜精品一区二区三区免费看| 国产极品天堂在线| 日本欧美国产在线视频| 亚洲人成网站在线观看播放| av免费观看日本| 午夜激情欧美在线| 九九在线视频观看精品| 长腿黑丝高跟| 午夜精品国产一区二区电影 | 精品一区二区三区视频在线| 亚洲国产日韩欧美精品在线观看| 亚洲美女视频黄频| 99久国产av精品国产电影| 国产高清三级在线| 日韩欧美精品v在线| 卡戴珊不雅视频在线播放| 成人二区视频| 插逼视频在线观看| or卡值多少钱| 久久久精品大字幕| 亚洲精品成人久久久久久| 国产亚洲精品久久久com| 精品国内亚洲2022精品成人| 久久精品人妻少妇| av在线播放精品| 免费av观看视频| 色综合站精品国产| 久久精品国产99精品国产亚洲性色| 亚洲在久久综合| 国产成人福利小说| 十八禁国产超污无遮挡网站| 日韩av不卡免费在线播放| 我要搜黄色片| 精品无人区乱码1区二区| 成人毛片60女人毛片免费| 久久久久九九精品影院| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 亚洲国产高清在线一区二区三| 插阴视频在线观看视频| 成年女人看的毛片在线观看| 色播亚洲综合网| 国产精品乱码一区二三区的特点| 97超视频在线观看视频| 偷拍熟女少妇极品色| 日韩中字成人| 日韩一区二区三区影片| 国内精品美女久久久久久| 少妇猛男粗大的猛烈进出视频 | 日日摸夜夜添夜夜添av毛片| 美女黄网站色视频| 中文在线观看免费www的网站| 精华霜和精华液先用哪个| 免费在线观看成人毛片| 全区人妻精品视频| av在线播放精品| av在线亚洲专区| av在线天堂中文字幕| 成人美女网站在线观看视频| 色综合亚洲欧美另类图片| 亚洲欧美精品综合久久99| 91午夜精品亚洲一区二区三区| 午夜老司机福利剧场| 亚洲,欧美,日韩| 又粗又硬又长又爽又黄的视频 | 久久鲁丝午夜福利片| 99在线人妻在线中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 国产爱豆传媒在线观看| 欧美成人一区二区免费高清观看| 赤兔流量卡办理| 中文精品一卡2卡3卡4更新| 国产久久久一区二区三区| av在线观看视频网站免费| 1024手机看黄色片| 国产高潮美女av| 午夜爱爱视频在线播放| 午夜福利在线观看免费完整高清在 | av国产免费在线观看| 亚洲av第一区精品v没综合| 免费看日本二区| 日本五十路高清| 99久国产av精品国产电影| 国产精品野战在线观看| 成人毛片a级毛片在线播放| 校园人妻丝袜中文字幕| 国产成人a区在线观看| 99久久精品热视频| 成人鲁丝片一二三区免费| 女人被狂操c到高潮| 99久久九九国产精品国产免费| 久久久精品94久久精品| 国产精品麻豆人妻色哟哟久久 | 国产精品99久久久久久久久| 亚洲精品色激情综合| 中文在线观看免费www的网站| 久久这里只有精品中国| 国产精品伦人一区二区| 91久久精品电影网| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 日本免费a在线| 91精品一卡2卡3卡4卡| 国产蜜桃级精品一区二区三区| 成年免费大片在线观看| 久久热精品热| 美女 人体艺术 gogo| 国产成人精品久久久久久| 毛片一级片免费看久久久久| 欧美日本视频| 噜噜噜噜噜久久久久久91| 国内久久婷婷六月综合欲色啪| 91av网一区二区| 如何舔出高潮| 色噜噜av男人的天堂激情| 嘟嘟电影网在线观看| 久久热精品热| 午夜精品在线福利| 成人毛片a级毛片在线播放| 毛片女人毛片| 网址你懂的国产日韩在线| 校园人妻丝袜中文字幕| 国产伦精品一区二区三区四那| 亚州av有码| 国产日韩欧美在线精品| 国产午夜精品久久久久久一区二区三区| 能在线免费观看的黄片| 亚洲电影在线观看av| av卡一久久| av天堂中文字幕网| 日本色播在线视频| videossex国产| 国产探花极品一区二区| 国产精品一及| 成年av动漫网址| 亚洲欧美日韩无卡精品| 美女 人体艺术 gogo| 12—13女人毛片做爰片一| 日韩av在线大香蕉| 国语自产精品视频在线第100页| 一夜夜www| 亚洲国产精品合色在线| 美女xxoo啪啪120秒动态图| 波多野结衣巨乳人妻| 日韩三级伦理在线观看| 亚洲久久久久久中文字幕| 亚洲成人中文字幕在线播放| 成年免费大片在线观看| 此物有八面人人有两片| 国产黄色视频一区二区在线观看 | 一个人看的www免费观看视频| 国产黄色小视频在线观看| 麻豆成人av视频| 欧美日韩乱码在线| 亚洲av二区三区四区| 国产精品久久视频播放| 高清午夜精品一区二区三区 | 国产熟女欧美一区二区| 毛片一级片免费看久久久久| 简卡轻食公司| 亚洲,欧美,日韩| 99国产精品一区二区蜜桃av| 国产成人午夜福利电影在线观看| 成年女人永久免费观看视频| 99热这里只有是精品50| 国产精品一及| 国产精品麻豆人妻色哟哟久久 | 国产一级毛片在线| 在线观看一区二区三区| 欧美3d第一页| 国产精品电影一区二区三区| 99久久精品热视频| av在线老鸭窝| 国产黄a三级三级三级人| 精品久久久噜噜| 久久久久久国产a免费观看| 久久这里有精品视频免费| 欧美精品国产亚洲| 成人欧美大片| 综合色av麻豆| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 少妇熟女aⅴ在线视频| av在线播放精品| 欧美最新免费一区二区三区| 成人午夜高清在线视频| 欧美+日韩+精品| 嫩草影院新地址| 国产精品乱码一区二三区的特点| 免费电影在线观看免费观看| 九九热线精品视视频播放| 亚洲第一电影网av| 久久久久久久久久黄片| 老师上课跳d突然被开到最大视频| 国产一区二区在线观看日韩| 女人被狂操c到高潮| 亚洲精品日韩av片在线观看| 国产精品蜜桃在线观看 | a级毛片a级免费在线| 久久精品国产亚洲av涩爱 | 免费观看在线日韩| 久久久久网色| 国产精品福利在线免费观看| 欧美3d第一页| 久久99热这里只有精品18| 内地一区二区视频在线| 国产伦精品一区二区三区视频9| 精品久久久久久成人av| 婷婷六月久久综合丁香| 高清在线视频一区二区三区 | 亚洲欧美日韩无卡精品| 国产精品蜜桃在线观看 | 99在线视频只有这里精品首页| 91久久精品电影网| 桃色一区二区三区在线观看| 亚洲精品日韩av片在线观看| 如何舔出高潮| 国产亚洲91精品色在线| 国内精品一区二区在线观看| 成人av在线播放网站| 国内揄拍国产精品人妻在线| 亚洲在线自拍视频| 少妇猛男粗大的猛烈进出视频 | 亚洲人成网站在线播| 男女边吃奶边做爰视频| 黄色一级大片看看| 麻豆乱淫一区二区| 国产精品一区二区三区四区久久| 久久久久久久久久成人| 国产成人影院久久av| 久久99蜜桃精品久久| 日本一本二区三区精品| 国产黄色小视频在线观看| 两个人的视频大全免费|