• <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
    无限看片的www在线观看| 九色成人免费人妻av| 日日干狠狠操夜夜爽| 在线观看日韩欧美| 国产在线精品亚洲第一网站| 久久精品91蜜桃| 国产真人三级小视频在线观看| 亚洲性夜色夜夜综合| 亚洲欧美日韩高清专用| 免费人成在线观看视频色| 午夜两性在线视频| 国产精品一区二区三区四区久久| 一本久久中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国产伦在线观看视频一区| 精品人妻偷拍中文字幕| 久久久久久大精品| xxxwww97欧美| 欧美乱色亚洲激情| 国产真人三级小视频在线观看| 少妇裸体淫交视频免费看高清| 欧美日韩一级在线毛片| 一级毛片高清免费大全| 国产高清视频在线播放一区| 久久精品影院6| 亚洲五月婷婷丁香| 啦啦啦观看免费观看视频高清| 国产精品亚洲一级av第二区| 国产伦精品一区二区三区视频9 | 在线播放无遮挡| 久久久久久人人人人人| 九九热线精品视视频播放| 99在线视频只有这里精品首页| 亚洲人成网站高清观看| 一区二区三区激情视频| 少妇的逼好多水| 国产精品嫩草影院av在线观看 | 给我免费播放毛片高清在线观看| 国产色婷婷99| 两个人的视频大全免费| 欧美在线一区亚洲| 亚洲一区高清亚洲精品| 欧美日韩瑟瑟在线播放| 久久99热这里只有精品18| 免费人成视频x8x8入口观看| 亚洲第一电影网av| 免费看光身美女| 色在线成人网| 叶爱在线成人免费视频播放| 亚洲av电影不卡..在线观看| 亚洲色图av天堂| 少妇丰满av| 狠狠狠狠99中文字幕| www国产在线视频色| 中文资源天堂在线| 国产一区在线观看成人免费| 一区二区三区免费毛片| 中文字幕久久专区| 国产精品影院久久| а√天堂www在线а√下载| 国产单亲对白刺激| 在线十欧美十亚洲十日本专区| 色在线成人网| 亚洲精华国产精华精| 欧洲精品卡2卡3卡4卡5卡区| 国产一级毛片七仙女欲春2| 亚洲av美国av| 国产麻豆成人av免费视频| 真人做人爱边吃奶动态| 日韩成人在线观看一区二区三区| 天美传媒精品一区二区| 国内久久婷婷六月综合欲色啪| 黄色日韩在线| 国产精品,欧美在线| 久久久久九九精品影院| 99热这里只有是精品50| 97碰自拍视频| 岛国在线免费视频观看| 国产成人啪精品午夜网站| 波野结衣二区三区在线 | 国产精品亚洲av一区麻豆| 欧美绝顶高潮抽搐喷水| 99久久精品热视频| 桃红色精品国产亚洲av| 国产探花在线观看一区二区| 日本五十路高清| av在线天堂中文字幕| 亚洲成a人片在线一区二区| 久久久久免费精品人妻一区二区| 亚洲美女黄片视频| 色精品久久人妻99蜜桃| 日日夜夜操网爽| 日本免费一区二区三区高清不卡| 久久久国产精品麻豆| 免费在线观看成人毛片| 国产毛片a区久久久久| 久久这里只有精品中国| 十八禁人妻一区二区| 久久国产精品人妻蜜桃| 91字幕亚洲| 舔av片在线| 又爽又黄无遮挡网站| 99久久九九国产精品国产免费| 国产精品香港三级国产av潘金莲| 欧美乱妇无乱码| 99热精品在线国产| 欧美在线黄色| 亚洲精华国产精华精| 男女那种视频在线观看| 日韩人妻高清精品专区| av视频在线观看入口| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人免费电影在线观看| 午夜福利免费观看在线| 国产精品久久电影中文字幕| 午夜两性在线视频| 亚洲国产欧美人成| 欧美另类亚洲清纯唯美| 免费在线观看影片大全网站| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av电影不卡..在线观看| 一进一出抽搐gif免费好疼| 成人特级av手机在线观看| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品一区二区www| 亚洲精品成人久久久久久| 夜夜爽天天搞| av专区在线播放| 神马国产精品三级电影在线观看| 欧美黑人巨大hd| 久久久久久久亚洲中文字幕 | 亚洲av日韩精品久久久久久密| 欧美黄色淫秽网站| 在线天堂最新版资源| 亚洲国产欧洲综合997久久,| 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av涩爱 | 欧美成人一区二区免费高清观看| 看片在线看免费视频| 最好的美女福利视频网| 18禁裸乳无遮挡免费网站照片| 免费观看精品视频网站| 一区福利在线观看| 青草久久国产| 少妇的逼水好多| 国产真实伦视频高清在线观看 | 午夜亚洲福利在线播放| 啦啦啦观看免费观看视频高清| av中文乱码字幕在线| 欧美成人性av电影在线观看| 中文字幕久久专区| avwww免费| 禁无遮挡网站| 18+在线观看网站| 18禁黄网站禁片免费观看直播| 中国美女看黄片| 成人欧美大片| 国产精品av视频在线免费观看| 免费观看的影片在线观看| 麻豆一二三区av精品| 亚洲乱码一区二区免费版| 欧美在线一区亚洲| 免费av观看视频| 一区福利在线观看| 老司机午夜十八禁免费视频| av福利片在线观看| 最近视频中文字幕2019在线8| 精品午夜福利视频在线观看一区| 身体一侧抽搐| 亚洲性夜色夜夜综合| 日韩欧美在线二视频| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久午夜电影| 欧美国产日韩亚洲一区| 1000部很黄的大片| 最近最新免费中文字幕在线| 真人做人爱边吃奶动态| 色噜噜av男人的天堂激情| 欧美成狂野欧美在线观看| 18+在线观看网站| 国产主播在线观看一区二区| 中文字幕人成人乱码亚洲影| 成人av在线播放网站| 国产成+人综合+亚洲专区| www.色视频.com| 国产淫片久久久久久久久 | 久久精品91蜜桃| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 悠悠久久av| 男女之事视频高清在线观看| 成年女人毛片免费观看观看9| 黄色日韩在线| 免费大片18禁| 99热精品在线国产| 男女做爰动态图高潮gif福利片| 99精品久久久久人妻精品| 内射极品少妇av片p| 久久久久久久久久黄片| 色综合欧美亚洲国产小说| 午夜两性在线视频| 最近最新中文字幕大全免费视频| 免费人成在线观看视频色| 亚洲欧美日韩无卡精品| 免费观看的影片在线观看| 欧美日韩国产亚洲二区| 草草在线视频免费看| 小蜜桃在线观看免费完整版高清| 国产成人aa在线观看| 亚洲av电影不卡..在线观看| 一级毛片高清免费大全| 欧美一区二区精品小视频在线| 欧美av亚洲av综合av国产av| 乱人视频在线观看| 18禁黄网站禁片午夜丰满| 天堂影院成人在线观看| 动漫黄色视频在线观看| 91麻豆精品激情在线观看国产| 午夜福利免费观看在线| 欧美中文日本在线观看视频| 欧美激情在线99| 欧美三级亚洲精品| 欧美午夜高清在线| 国产精华一区二区三区| 精品一区二区三区av网在线观看| 欧美最新免费一区二区三区 | 欧美绝顶高潮抽搐喷水| 男女做爰动态图高潮gif福利片| 一本久久中文字幕| 欧美大码av| 久久精品人妻少妇| 日本免费a在线| 日本三级黄在线观看| 又黄又爽又免费观看的视频| 亚洲美女黄片视频| 欧美成人a在线观看| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩东京热| 国产主播在线观看一区二区| 青草久久国产| 91九色精品人成在线观看| 国产一区二区在线观看日韩 | 最近最新中文字幕大全电影3| 成人av一区二区三区在线看| 免费大片18禁| 国产精品亚洲av一区麻豆| 亚洲中文日韩欧美视频| 国产乱人视频| 国产一区二区在线观看日韩 | 国产亚洲精品久久久com| 国产中年淑女户外野战色| 99热这里只有是精品50| 亚洲国产日韩欧美精品在线观看 | 色播亚洲综合网| 色综合欧美亚洲国产小说| 一级毛片高清免费大全| 国产老妇女一区| 波多野结衣高清无吗| 欧美黑人欧美精品刺激| 听说在线观看完整版免费高清| 午夜福利在线观看吧| 变态另类成人亚洲欧美熟女| 亚洲专区国产一区二区| 神马国产精品三级电影在线观看| a级毛片a级免费在线| 男人和女人高潮做爰伦理| 久久国产精品人妻蜜桃| 中文字幕高清在线视频| 欧美+日韩+精品| 黄色成人免费大全| 三级国产精品欧美在线观看| 国产男靠女视频免费网站| 97人妻精品一区二区三区麻豆| 五月伊人婷婷丁香| 高清日韩中文字幕在线| or卡值多少钱| av视频在线观看入口| 最后的刺客免费高清国语| 精品熟女少妇八av免费久了| 色视频www国产| 午夜精品一区二区三区免费看| 久久久久久九九精品二区国产| 超碰av人人做人人爽久久 | 日韩国内少妇激情av| 国产免费av片在线观看野外av| 九九热线精品视视频播放| 黄片大片在线免费观看| 女人被狂操c到高潮| 可以在线观看毛片的网站| 久久精品亚洲精品国产色婷小说| 日韩国内少妇激情av| 国内精品美女久久久久久| 精华霜和精华液先用哪个| 免费人成视频x8x8入口观看| 99久久综合精品五月天人人| 亚洲一区二区三区不卡视频| 变态另类成人亚洲欧美熟女| 国产乱人视频| 久久亚洲精品不卡| 国产三级黄色录像| 欧美zozozo另类| 禁无遮挡网站| 国内精品久久久久久久电影| 精品日产1卡2卡| 淫妇啪啪啪对白视频| 搡老熟女国产l中国老女人| 在线观看66精品国产| 伊人久久大香线蕉亚洲五| www.色视频.com| 偷拍熟女少妇极品色| 在线观看午夜福利视频| xxxwww97欧美| 淫妇啪啪啪对白视频| 岛国在线观看网站| 老汉色∧v一级毛片| 18禁黄网站禁片免费观看直播| 国产高清视频在线观看网站| 欧美性感艳星| 久久久久九九精品影院| tocl精华| 国产高清三级在线| 欧美最新免费一区二区三区 | 免费av观看视频| 91av网一区二区| av在线天堂中文字幕| 国内精品久久久久久久电影| 真人做人爱边吃奶动态| 老熟妇仑乱视频hdxx| 亚洲熟妇熟女久久| 国产亚洲精品久久久com| 久久国产精品人妻蜜桃| 欧美区成人在线视频| 久久精品国产亚洲av香蕉五月| 啪啪无遮挡十八禁网站| 精品久久久久久成人av| 激情在线观看视频在线高清| 国产真人三级小视频在线观看| 精品免费久久久久久久清纯| 久久草成人影院| 99国产精品一区二区三区| 日韩欧美精品v在线| 一个人免费在线观看电影| 狂野欧美激情性xxxx| 99国产精品一区二区三区| av国产免费在线观看| 日韩av在线大香蕉| 久久伊人香网站| 久久香蕉精品热| 99热只有精品国产| 免费观看的影片在线观看| 国产成人a区在线观看| 国产午夜福利久久久久久| 亚洲精品美女久久久久99蜜臀| 国产淫片久久久久久久久 | x7x7x7水蜜桃| 香蕉久久夜色| 精品久久久久久,| 两个人视频免费观看高清| 蜜桃久久精品国产亚洲av| 麻豆成人av在线观看| 国产蜜桃级精品一区二区三区| 亚洲精品久久国产高清桃花| 国产又黄又爽又无遮挡在线| 国产一级毛片七仙女欲春2| 九色成人免费人妻av| 少妇人妻精品综合一区二区 | 天堂√8在线中文| bbb黄色大片| 国产亚洲精品一区二区www| 91av网一区二区| 欧美黑人巨大hd| xxx96com| 在线观看av片永久免费下载| 亚洲无线在线观看| 久久久国产成人精品二区| 欧美精品啪啪一区二区三区| 草草在线视频免费看| 黄片大片在线免费观看| 免费人成在线观看视频色| 国产视频一区二区在线看| 波多野结衣高清作品| 亚洲人与动物交配视频| 欧美日韩国产亚洲二区| www日本黄色视频网| 欧美成人a在线观看| 国产在线精品亚洲第一网站| 日本免费一区二区三区高清不卡| 精品久久久久久久毛片微露脸| 亚洲第一欧美日韩一区二区三区| 日韩精品中文字幕看吧| 久久精品夜夜夜夜夜久久蜜豆| 久久久久九九精品影院| 美女被艹到高潮喷水动态| 神马国产精品三级电影在线观看| 一个人免费在线观看电影| 色尼玛亚洲综合影院| 亚洲一区二区三区色噜噜| 亚洲av电影在线进入| 欧美成人免费av一区二区三区| 女同久久另类99精品国产91| 亚洲熟妇熟女久久| 国产69精品久久久久777片| 国产乱人伦免费视频| 日本精品一区二区三区蜜桃| 99热6这里只有精品| 国产毛片a区久久久久| 岛国视频午夜一区免费看| 午夜激情欧美在线| 岛国在线观看网站| 美女高潮喷水抽搐中文字幕| 三级毛片av免费| 亚洲最大成人中文| 日韩欧美三级三区| 免费在线观看成人毛片| 老熟妇乱子伦视频在线观看| 天堂动漫精品| 国产视频内射| 在线a可以看的网站| 九色成人免费人妻av| 少妇的丰满在线观看| 国产亚洲精品久久久久久毛片| 久久精品影院6| 老司机午夜福利在线观看视频| 国产视频一区二区在线看| 一个人免费在线观看电影| 久久国产精品人妻蜜桃| av天堂在线播放| 成人无遮挡网站| a级一级毛片免费在线观看| 久久精品国产自在天天线| 国产精华一区二区三区| 国产熟女xx| 91麻豆精品激情在线观看国产| 午夜激情福利司机影院| h日本视频在线播放| 国产精品三级大全| 欧美极品一区二区三区四区| 日本黄色视频三级网站网址| 黄色视频,在线免费观看| 99久久成人亚洲精品观看| 午夜福利高清视频| 99久久无色码亚洲精品果冻| 可以在线观看毛片的网站| 少妇丰满av| 操出白浆在线播放| 成人永久免费在线观看视频| 欧美午夜高清在线| 日本撒尿小便嘘嘘汇集6| 中文亚洲av片在线观看爽| or卡值多少钱| 午夜久久久久精精品| 亚洲av成人不卡在线观看播放网| 亚洲aⅴ乱码一区二区在线播放| 欧美最黄视频在线播放免费| 国产一区在线观看成人免费| 久久久久久久久久黄片| 精品无人区乱码1区二区| h日本视频在线播放| 成人特级av手机在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人不卡在线观看播放网| 欧美xxxx黑人xx丫x性爽| 9191精品国产免费久久| 国产私拍福利视频在线观看| 夜夜夜夜夜久久久久| 午夜福利在线观看免费完整高清在 | 国产欧美日韩一区二区精品| 欧美日韩国产亚洲二区| 亚洲18禁久久av| 欧美日韩亚洲国产一区二区在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美一区二区三区黑人| 淫妇啪啪啪对白视频| 欧美bdsm另类| 一进一出抽搐动态| 女同久久另类99精品国产91| av在线天堂中文字幕| 夜夜躁狠狠躁天天躁| 国内揄拍国产精品人妻在线| x7x7x7水蜜桃| 国产99白浆流出| 欧美成人性av电影在线观看| a级一级毛片免费在线观看| 精品国产美女av久久久久小说| aaaaa片日本免费| 国产亚洲精品久久久久久毛片| 一个人免费在线观看电影| 可以在线观看毛片的网站| 精品国产美女av久久久久小说| 欧美3d第一页| 亚洲专区中文字幕在线| 99riav亚洲国产免费| 一个人看视频在线观看www免费 | 夜夜躁狠狠躁天天躁| 一区福利在线观看| e午夜精品久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 欧美成人一区二区免费高清观看| 国产精品亚洲美女久久久| 长腿黑丝高跟| 日韩欧美国产在线观看| 久久精品国产自在天天线| aaaaa片日本免费| 日韩欧美 国产精品| 一区二区三区免费毛片| 日本熟妇午夜| 嫩草影院精品99| 成人欧美大片| 深夜精品福利| 在线播放国产精品三级| 啦啦啦免费观看视频1| 少妇人妻精品综合一区二区 | 久久天躁狠狠躁夜夜2o2o| 天天一区二区日本电影三级| 岛国在线免费视频观看| 一级毛片高清免费大全| 欧美黑人欧美精品刺激| 亚洲第一电影网av| 日韩 欧美 亚洲 中文字幕| 一级黄片播放器| 我的老师免费观看完整版| 亚洲精品影视一区二区三区av| 日韩国内少妇激情av| 国产精品自产拍在线观看55亚洲| 美女被艹到高潮喷水动态| 亚洲av免费在线观看| 99久久成人亚洲精品观看| 少妇丰满av| 精品午夜福利视频在线观看一区| 亚洲真实伦在线观看| 国产成人av教育| 欧美一区二区精品小视频在线| 午夜精品在线福利| 午夜两性在线视频| 中国美女看黄片| 国产精品1区2区在线观看.| 岛国在线免费视频观看| 无限看片的www在线观看| 极品教师在线免费播放| 午夜日韩欧美国产| 麻豆国产97在线/欧美| 夜夜躁狠狠躁天天躁| 黄片小视频在线播放| svipshipincom国产片| 一区福利在线观看| 色av中文字幕| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 亚洲欧美日韩无卡精品| 成人鲁丝片一二三区免费| 国产三级在线视频| 夜夜夜夜夜久久久久| 午夜免费成人在线视频| 免费在线观看亚洲国产| 男人舔奶头视频| 在线观看免费午夜福利视频| 国产高潮美女av| 久99久视频精品免费| 国产爱豆传媒在线观看| 欧美日韩福利视频一区二区| 老汉色∧v一级毛片| 最近在线观看免费完整版| 悠悠久久av| 欧美日韩一级在线毛片| 老熟妇仑乱视频hdxx| 女警被强在线播放| 亚洲av成人不卡在线观看播放网| 亚洲人与动物交配视频| 又粗又爽又猛毛片免费看| 中文字幕精品亚洲无线码一区| 97碰自拍视频| aaaaa片日本免费| 久久人人精品亚洲av| 国产单亲对白刺激| www日本黄色视频网| 熟女少妇亚洲综合色aaa.| 国产精品av视频在线免费观看| 欧美国产日韩亚洲一区| www.色视频.com| 观看免费一级毛片| 国产精品一区二区三区四区久久| 免费在线观看日本一区| 久久久久久国产a免费观看| 波多野结衣巨乳人妻| a级毛片a级免费在线| 搡老熟女国产l中国老女人| 制服丝袜大香蕉在线| 国内精品一区二区在线观看| 99久久无色码亚洲精品果冻| 变态另类成人亚洲欧美熟女| 男插女下体视频免费在线播放| 亚洲中文字幕一区二区三区有码在线看| 亚洲av免费高清在线观看| 欧美一区二区精品小视频在线| 中文字幕精品亚洲无线码一区| 亚洲成人精品中文字幕电影| 精品午夜福利视频在线观看一区| 亚洲片人在线观看| 在线播放国产精品三级| 88av欧美| 波多野结衣巨乳人妻| 在线观看美女被高潮喷水网站 | 亚洲avbb在线观看| 欧美午夜高清在线| 日韩大尺度精品在线看网址| 女生性感内裤真人,穿戴方法视频| 国产亚洲精品综合一区在线观看| 国产伦人伦偷精品视频| 99久久九九国产精品国产免费| 最后的刺客免费高清国语| 麻豆久久精品国产亚洲av| 无人区码免费观看不卡| 欧美av亚洲av综合av国产av| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,|