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

    碘硫循環(huán)制氫中HI 濃縮精餾的模擬

    2020-03-03 11:11:58黃秀輝崔國(guó)民
    關(guān)鍵詞:板數(shù)精餾塔模擬計(jì)算

    王 俊, 黃秀輝, 崔國(guó)民

    (上海理工大學(xué)能源與動(dòng)力工程學(xué)院,上海 200093)

    氫能具有資源豐富、無污染、高揮發(fā)性、高能量等特點(diǎn),在工業(yè)生產(chǎn)中有廣泛應(yīng)用。傳統(tǒng)的制氫方法主要是用化學(xué)方法從天然氣、石油和煤等化石燃料中制氫,但從水中制氫的方法能解決氫的節(jié)能、高效、可持續(xù)大規(guī)模生產(chǎn)的問題[1]。

    在標(biāo)準(zhǔn)壓力下,要利用高于2 000 ℃的熱源才能直接熱解水制氫。而采用熱化學(xué)循環(huán)的方法能夠較大程度地降低操作條件[2],熱化學(xué)循環(huán)分解水制氫是在不同階段、不同溫度下把熱量供給含有中間介質(zhì)的水的分解系統(tǒng)最終分解為氫氣和氧氣,而其余產(chǎn)物作為中間產(chǎn)物或中間反應(yīng)物循環(huán)利用[3]。當(dāng)今碘硫(I-S)循環(huán)被認(rèn)為是眾多熱化學(xué)循環(huán)分解水制氫路線中的優(yōu)選方案之一,它可以在溫度低于1 000 ℃時(shí)進(jìn)行,制氫率達(dá)到52%左右[4]。但由于碘化氫(HI)與水的共沸現(xiàn)象,導(dǎo)致了過程的復(fù)雜性增加,因此HI 的濃縮分離成為了碘硫循環(huán)能否高效運(yùn)行的瓶頸步驟。美國(guó)GA 公司于20 世紀(jì)70 年代提出磷酸萃取的技術(shù)路線[5],向HI-H20-I2中加入H3PO4作為萃取劑,難溶于水的I2會(huì)首先從料液中析出,H3PO4與H2O 的強(qiáng)締合作用減少水與HI 的共沸現(xiàn)象,對(duì)該體系進(jìn)行精餾即可在塔頂?shù)玫礁邼舛鹊腍I,而精餾塔釜底的稀磷酸經(jīng)濃縮除水后可循環(huán)利用。德國(guó)亞琛工業(yè)大學(xué)(RWTH)于1987 年提出了反應(yīng)精餾方法[6],在高溫高壓下的反應(yīng)精餾塔內(nèi)進(jìn)行精餾和HI 分解的耦合操作,但對(duì)設(shè)備要求過高,大量生產(chǎn)難度較大。日本原子力機(jī)構(gòu)(JAEA)在1997 年提出了一個(gè)更為新穎的研究方法,利用電解電滲析預(yù)濃縮-精餾,也稱為膜電解工藝(Electro-Electrodialysis,EED)對(duì)HI-H20-I2體系進(jìn)行預(yù)濃縮處理[7],再進(jìn)行精餾,相對(duì)而言,該方法效率高,條件更為合適。

    目前,HI 濃縮精餾環(huán)節(jié)的設(shè)計(jì)模擬都停留在簡(jiǎn)化模型基礎(chǔ)上,對(duì)于該體系內(nèi)三元物系的物性及二元交互作用參數(shù)都還沒有系統(tǒng)完整的研究。本文重點(diǎn)從最基礎(chǔ)的物性研究出發(fā),針對(duì)HI-I2-H2O 的三元共沸體系,選用合適的熱力學(xué)模型來修正氣液平衡的非理想性,同時(shí)結(jié)合相平衡的實(shí)驗(yàn)數(shù)據(jù),采用極大似然法關(guān)聯(lián)得到二元交互作用參數(shù)。本文通過Aspen Plus 平臺(tái)和相圖分析所得體系的熱力學(xué)精餾特性對(duì)該精餾塔進(jìn)行模擬設(shè)計(jì),對(duì)關(guān)鍵操作參數(shù)對(duì)于過程分離效果和能耗的影響進(jìn)行了分析,并得到了相關(guān)參數(shù)的優(yōu)化設(shè)計(jì)結(jié)果。

    1 碘硫循環(huán)及HI 精餾的熱力學(xué)模型

    碘硫循環(huán)包括3 個(gè)部分:本生(Bunsen)反應(yīng)[8]、硫酸分解(SA 部分)、氫碘酸分解(HI 部分):

    (1)本生反應(yīng)(放熱,反應(yīng)溫度20~120 ℃)

    (2)硫酸分解反應(yīng)(吸熱,反應(yīng)溫度800~900 ℃)

    (3)氫碘酸分解反應(yīng)(吸熱,反應(yīng)溫度300~500 ℃)

    圖1 所示為碘硫循環(huán)流程圖[9],由圖可見,硫酸和氫碘酸是由水、碘和二氧化硫反應(yīng)而得。在反應(yīng)中當(dāng)存在過量的碘時(shí)這兩種酸將自發(fā)分離成兩層不混合的溶液相,其上層(H2SO4相)含有大量H2SO4和少量的HI 和I2,而下層(HI-I2-H2O)中含有大量的HI 和I2,但是也有少量硫酸的存在[10]。用本生反應(yīng)的逆反應(yīng)分別對(duì)兩層液體分離和提純以消除雜質(zhì)(即在H2SO4中的HI、I2,HI-I2-H2O 相中的H2SO4)。進(jìn)行分離和純化時(shí)兩種酸離開本生反應(yīng)部分,硫酸支流進(jìn)入到SA 部分被濃縮并分解;純化的HI-I2-H2O 進(jìn)入到HI 部分,被分離,整流,并分解產(chǎn)生氫。在SA 部分產(chǎn)生SO2和H2以及在HI 部分產(chǎn)生的I2循環(huán)回到本生反應(yīng)。碘硫循環(huán)過程的凈投入為水和用來分解H2SO4和HI 的高溫?zé)幔瑑糨敵鰹闅?、氧和本生反?yīng)產(chǎn)生的低溫?zé)?,整個(gè)過程的效率主要取決于H2SO4和HI 的分解和濃縮需求的熱量。HI 分解部分被認(rèn)為是整個(gè)過程優(yōu)化的關(guān)鍵部分,因?yàn)镠II2-H2O 中來自本生反應(yīng)的混合物很難分離且HI 分解反應(yīng)的轉(zhuǎn)化率低(約20%),是整個(gè)過程中高耗能的部分[11]。所以對(duì)于整個(gè)提純體系,本文采用精餾塔分離得到高濃度HI 這一環(huán)節(jié)有至關(guān)重要的影響。

    精餾過程的實(shí)質(zhì)是精餾塔內(nèi)所有物質(zhì)的氣液相之間通過氣液相平衡(VLE)實(shí)現(xiàn)輕重組分間的傳質(zhì)過程,因此通過混合物相平衡的研究獲得整個(gè)精餾過程所需的熱力學(xué)數(shù)據(jù)及模型是所有精餾過程操作分析和模擬計(jì)算的重要基礎(chǔ)。針對(duì)HI-I2-H2O 三元體系的物系性質(zhì)及特點(diǎn),文獻(xiàn)[12]表明NRTL(Non-Random Two-Liquids)模型能較好地修正氣液相平衡的非理想性,因此本研究選用NRTL 模型,該模型中未知的二元交互作用參數(shù)表達(dá)示如式(4)所示:

    式中:ia,ib 分別代表兩種組分;二元交互作用參數(shù)Aia,ib和 Bia,ib可以通過已有相平衡實(shí)驗(yàn)數(shù)據(jù)進(jìn)行回歸擬合,得到活度系數(shù)方程中的能量參數(shù) τia,ib。本研究基于Larousse 等[13]給出的氣液平衡實(shí)驗(yàn)數(shù)據(jù),采用極大似然法關(guān)聯(lián)參數(shù):

    圖 1 碘硫循環(huán)流程圖Fig. 1 Iodine sulfur cycle flow chart

    式中: Q 為參數(shù)關(guān)聯(lián)所要求取最小值的目標(biāo)函數(shù);NP為數(shù)據(jù)點(diǎn)數(shù)目; Nc為組分個(gè)數(shù);T 為溫度;y 為氣相摩爾分?jǐn)?shù); σ 為相應(yīng)關(guān)聯(lián)數(shù)據(jù)的標(biāo)準(zhǔn)偏差( σT=0.05 ℃, σy=0.003),上標(biāo) e 和m 分別表示預(yù)計(jì)數(shù)據(jù)(Estimated data)和實(shí)驗(yàn)測(cè)量數(shù)據(jù)(Measured data)。

    基于Aspen Plus 軟件平臺(tái)對(duì)NRTL 模型中的二元交互作用參數(shù) Aia,ib和 Bia,ib進(jìn)行回歸擬合所得的二元交互作用參數(shù)如表1 所示, Cia,ib為混合有序特性常數(shù),根據(jù)文獻(xiàn)資料[14]取值為0.300。

    表 1 NRTL 方程的二元交互作用參數(shù)Table 1 Binary interaction parameters of NRTL equation

    圖2 示出了采用二元交互作用參數(shù)預(yù)測(cè)的H2OHI 氣液平衡相圖與實(shí)驗(yàn)數(shù)據(jù)[15]的對(duì)比結(jié)果。從圖中可以看出,預(yù)測(cè)結(jié)果與實(shí)驗(yàn)數(shù)據(jù)無明顯差異,趨勢(shì)一致,吻合良好,驗(yàn)證了所得二元交互作用參數(shù)的準(zhǔn)確性及可靠性,可以進(jìn)一步用于預(yù)測(cè)該三元體系中其他的二元相平衡及三元相平衡。

    圖 2 H2O-HI 氣液平衡相圖與實(shí)驗(yàn)點(diǎn)的對(duì)比Fig. 2 Comparison of H2O-HI vapor-liquid equilibrium phase diagram and experimental point

    圖3 示出了HI-I2-H2O 體系基于所得二元交互作用參數(shù)的相平衡預(yù)測(cè)圖,其中圖3(a)、圖3(b)和圖3(c)分別為H2O-HI、H2O-I2和HI-I2三組二元體系的氣液平衡相圖。由圖3(a)可以看出,H2O-HI 二元組分在126.7 ℃附近兩曲線相交于一點(diǎn),此點(diǎn)處水的氣相摩爾分?jǐn)?shù)和液相摩爾分?jǐn)?shù)相等,均為0.81,表明HI和H2O 發(fā)生了共沸,共沸物中含HI 的物質(zhì)的量分?jǐn)?shù)為0.19,這對(duì)于HI 的濃縮精餾有著關(guān)鍵的影響。圖3(d)示出了為基于擬合參數(shù)所預(yù)測(cè)的HI-H2O-I2剩余曲線圖,即精餾曲線圖。該三元體系中,HI 沸點(diǎn)為-35.38 ℃,H2O 沸點(diǎn)為100.02 ℃,I2沸點(diǎn)為184.39 ℃,從圖中可以看出,HI 和H2O 的共沸點(diǎn)位于D 點(diǎn)(126.74 ℃),AD 為精餾邊界線,即當(dāng)精餾塔進(jìn)料組分處于ABD 區(qū)域內(nèi)時(shí),即進(jìn)料混合物中HI 濃度小于0.19 時(shí),精餾塔內(nèi)塔頂最多得到該區(qū)域內(nèi)沸點(diǎn)最低的純組分H2O,塔底最多得到沸點(diǎn)最高的純組分I2,無論如何也不可能得到高純度HI。而大多數(shù)碘硫循環(huán)工藝中本生反應(yīng)階段出來的HI-I2-H2O 混合物中HI 的濃度低于0.19,因此大多數(shù)工藝在本生反應(yīng)和精餾濃縮之間增加了電解電滲析預(yù)濃縮-精餾(EED)環(huán)節(jié),就是為了讓HI 濃縮精餾塔的進(jìn)料組分處于圖3(d)中的ADC 區(qū)域內(nèi),以便在精餾塔頂部獲得高濃度低沸點(diǎn)物質(zhì)HI。圖3(d)中E 點(diǎn)為文獻(xiàn)[16]中經(jīng)EED 預(yù)濃縮后的HI 濃縮精餾塔進(jìn)料組成,經(jīng)過該點(diǎn)的曲線CEA 即為該進(jìn)料條件下精餾塔內(nèi)塔頂?shù)剿慕M分分布預(yù)測(cè)曲線,可以根據(jù)該曲線來進(jìn)行后期精餾塔模擬計(jì)算的初值預(yù)設(shè)。

    2 HI 濃縮精餾的模擬與設(shè)計(jì)

    某個(gè)EED 工藝的HI 濃縮精餾的流程如圖4 所示。從本生反應(yīng)得到的原料液F1,被分成兩個(gè)流股(F2、F3)分別進(jìn)入EED 的陽極和陰極,而通過EED 工藝處理后的流股,因?yàn)镠I 富集在陰極側(cè)和出口流體的陰極側(cè),流股F9 被返回到本生反應(yīng)中作為原料液,而流股F8 進(jìn)入到換熱器中進(jìn)行加熱,達(dá)到泡點(diǎn)溫度,流股F10 進(jìn)入精餾塔,塔頂餾出高純度的HI,高純度HI 的F11 流股給分解反應(yīng)提供原料,產(chǎn)出大量氫氣,而塔底產(chǎn)品F12通過換熱器和F13 流股,分別與F3 和F2 流股混合,返回到EED 工藝中補(bǔ)充流失的HI。

    對(duì)于該HI-H20-I2體系的精餾塔的模擬計(jì)算,其數(shù)學(xué)穩(wěn)態(tài)平衡級(jí)模型主要是由物料平衡方程(M)、相平衡方程(E)、濃度加和方程(S)以及熱量平衡方程(H)四大基本方程式構(gòu)成的方程組,通常也被稱為MESH 方程[17],如表2 所示。

    圖 3 HI-I2-H2O 體系相平衡預(yù)測(cè)圖(p=101.325 kPa)Fig. 3 Phase equilibrium prediction chart of HI-I2-H2O system (p=101.325 kPa)

    圖 4 EED 工藝的HI 精餾流程圖Fig. 4 HI distillation flow chart of EED process

    表 2 氣液兩相MESH 方程Table 2 MESH equations of vapor-liquid(VL) phase

    將各塔板的MESH 方程與冷凝器、再沸器方程以及相應(yīng)的物性方程等聯(lián)合求解即可得到每塊塔板上的氣液相摩爾分?jǐn)?shù)、塔板溫度以及氣液相流率等變量。本文擬設(shè)計(jì)負(fù)荷為1 000 mol/s 的HI 濃縮精餾塔進(jìn)行模擬優(yōu)化,將精餾塔看作理想塔,塔內(nèi)壓降為10 kPa,塔板效率為100%。對(duì)于精餾塔的進(jìn)料采用文獻(xiàn)[16]的進(jìn)料組分n(HI)∶n (I2)∶n (H2O)=0.129∶0.342∶0.529,進(jìn)料設(shè)為泡點(diǎn)進(jìn)料,在Aspen Plus 溫度模塊設(shè)Vapor fraction 為0;壓力設(shè)為1.17 MPa;流量設(shè)為1 000 mol/s。根據(jù)精餾塔快捷模擬計(jì)算及物料衡算,理論塔板數(shù)暫時(shí)設(shè)為12、冷凝器為部分冷凝、塔頂產(chǎn)品流量為380 mol/s、回流比為1.5,其他項(xiàng)目采用缺省值。物流的進(jìn)料設(shè)為第8 塊塔板,塔頂壓力設(shè)定為1.17 MPa。參數(shù)設(shè)置完畢,運(yùn)行模擬計(jì)算,結(jié)果如表3 所示,從表可以看出,塔頂出料F11 中,HI 的質(zhì)量流量為1.748×105kg/h,從冷凝器精餾出的HI 的質(zhì)量分?jǐn)?shù)為0.997,由輸出結(jié)果可以看到,塔底出料F12 的HI 的質(zhì)量分?jǐn)?shù)雖然只有0.330,但由于流量較大,所以精餾不是很徹底,說明模擬相關(guān)的參數(shù)設(shè)定輸入不是最優(yōu)值,還可以進(jìn)一步優(yōu)化。

    如圖5 所示為精餾塔內(nèi)的模擬計(jì)算結(jié)果分布,模擬計(jì)算結(jié)果中溫度在塔內(nèi)的分布如圖5(a)所示,氣相和液相在塔內(nèi)的組分分布分別如圖5(b)和5(c)所示。從圖5(a)可以看出,第1 塊板和第8 塊塔板溫差較大,可以考慮設(shè)置為靈敏板。從圖5(b)和5(c)可以看出,3 種組分在氣相和液相在塔內(nèi)的分布趨勢(shì)相似,在塔頂?shù)臍庀嗪鸵合喽嫉玫搅烁邼舛鹊腍I,并且在第7、第8 塊塔板數(shù)組分分布迅速變化,其中很大一部分影響應(yīng)該來自于進(jìn)料口塔板位置的設(shè)置。

    再沸器熱負(fù)荷是影響精餾過程費(fèi)用的主要因素,在其他參數(shù)保持不變的情況下,本文比較了不同塔板數(shù)下精餾出HI 的質(zhì)量分?jǐn)?shù)(w)及其對(duì)應(yīng)的再沸器的熱負(fù)荷(Q),整理數(shù)據(jù)繪制成如圖6 所示的塔板數(shù)對(duì)精餾效果的影響曲線圖。圖6(c)示出了單位質(zhì)量分?jǐn)?shù)對(duì)應(yīng)的再沸器熱負(fù)荷(Qm, 下同)與塔板數(shù)之間的變化關(guān)系。由曲線的變換趨勢(shì)可以看出,當(dāng)塔板數(shù)增加到13 時(shí),對(duì)應(yīng)的Qm值驟降,之后隨著塔板數(shù)的增加,Qm的降低趨于平緩??紤]到每增加一個(gè)塔板,相應(yīng)的設(shè)備費(fèi)用增加,所以本文推薦塔板數(shù)13 為精餾塔最佳值。

    表 3 精餾塔進(jìn)出料模擬計(jì)算結(jié)果Table 3 Simulated calculation results of feed and discharge in distillation column

    圖 5 精餾塔內(nèi)的模擬計(jì)算結(jié)果分布圖Fig. 5 Distillation column simulation calculation results distribution map

    圖 6 塔板數(shù)對(duì)精餾效果的影響曲線圖Fig. 6 Effect of the number of plates on the rectification performance

    回流比也是影響精餾塔分離效果及能耗的主要因素,不同回流比是通過改變回流量來決定的,本文通過設(shè)置不同回流量來比較精餾結(jié)果,整理數(shù)據(jù)繪制成如圖7 所示的回流比對(duì)精餾效果的影響曲線圖。由圖7(a)可以發(fā)現(xiàn),隨著回流比的增加,餾出HI 質(zhì)量分?jǐn)?shù)也在增加,但由圖7(b)可以發(fā)現(xiàn),回流比增加的同時(shí)也帶來了熱負(fù)荷的增加。單位質(zhì)量分?jǐn)?shù)對(duì)應(yīng)的再沸器熱負(fù)荷Qm與回流比之間的關(guān)系如圖7(c)所示,由圖可以看出Qm隨著回流比的增大而均勻增大,所以想得到更高純度的HI,再沸器的熱負(fù)荷需要增加。

    圖 7 回流比對(duì)精餾效果的影響曲線圖Fig. 7 Effect of reflux ratio on rectification performance

    不同進(jìn)料組分比例也會(huì)對(duì)精餾塔的分離效果和能耗產(chǎn)生影響[18],本文分別模擬計(jì)算了4 組不同EED 方法中預(yù)濃縮得到的HI-I2-H2O 進(jìn)料組分,結(jié)果見表4。從表4 中可以發(fā)現(xiàn)當(dāng)進(jìn)料采用第2 種組分占比時(shí),HI 濃度高,且耗能少。所以并不是HI 組分物質(zhì)的量之比越大,精餾出HI 的質(zhì)量分?jǐn)?shù)就越大,塔釜加熱量也不是隨著HI 組分物質(zhì)的量之比進(jìn)行單調(diào)變化,所以并不能盲目提高預(yù)濃縮的濃度,需要通過相關(guān)的模擬分析來得出最適合的濃度,以達(dá)到高濃度的分離指標(biāo)的同時(shí)實(shí)現(xiàn)低能耗。

    表 4 不同進(jìn)料組分下HI 精餾塔模擬結(jié)果Table 4 Simulation results of HI distillation column of different feeds

    3 結(jié) 論

    本文針對(duì)HI-I2-H2O 三元混合物系選用NRTL方程來修正氣液平衡的非理想性,結(jié)合已有的實(shí)驗(yàn)數(shù)據(jù),回歸擬合出NRTL 方程中的二元交互作用參數(shù),并基于精餾原理,進(jìn)料采用實(shí)際進(jìn)料組分,借助化工模擬軟件Aspen Plus 對(duì)HI 濃縮精餾過程進(jìn)行了模擬、分析及優(yōu)化,得到以下結(jié)論:

    (1)隨著塔板數(shù)的增加,單位質(zhì)量對(duì)應(yīng)的再沸器熱負(fù)荷Qm整體是降低的,但是當(dāng)塔板數(shù)增加超過13 時(shí),Qm的降低很緩慢,所以選擇13 為最佳塔板數(shù)。

    (2)隨著回流比的增加,單位質(zhì)量對(duì)應(yīng)的再沸器熱負(fù)荷Qm整體呈勻速增加。當(dāng)回流比為1.3 時(shí),餾出HI 的質(zhì)量分?jǐn)?shù)已經(jīng)達(dá)到0.96,之后隨著回流比的增加,餾出HI 的質(zhì)量分?jǐn)?shù)增加不明顯,再考慮到回流比增大會(huì)導(dǎo)致操作費(fèi)用和設(shè)備費(fèi)用也會(huì)相應(yīng)增加,所以本文回流比以1.3 為最佳。

    (3)進(jìn)料中的HI 質(zhì)量分?jǐn)?shù)需要通過具體模擬分析才能得出最優(yōu)值,與出料的HI 質(zhì)量分?jǐn)?shù)和塔內(nèi)能耗并不是單調(diào)變化的關(guān)系。

    符號(hào)說明

    i——關(guān)鍵組分,i = 1, 2, 3

    j——塔板數(shù),0 ≤j ≤N+1

    Vj, Lj, Fj——第j 塊塔板氣相、液相、進(jìn)料流率

    xj,i, yj,i, zj,i——組分i 在第j 塊塔板上的液相、氣相、進(jìn)料物質(zhì)的量分?jǐn)?shù)

    猜你喜歡
    板數(shù)精餾塔模擬計(jì)算
    二硫化碳-四氯化碳精餾工藝改進(jìn)與模擬研究
    一起精餾塔自燃事故的原因分析及防范措施的改進(jìn)
    R1234ze PVTx熱物性模擬計(jì)算
    能源工程(2022年1期)2022-03-29 01:06:26
    Aspen Plus 軟件在低溫甲醇洗工藝中的應(yīng)用
    化工管理(2020年19期)2020-07-28 02:59:02
    用Excel進(jìn)行精餾塔理論板數(shù)的計(jì)算
    山東化工(2018年19期)2018-10-29 08:12:20
    氨水吸收式制冷系統(tǒng)中精餾塔性能模擬與分析
    擠出發(fā)泡片材褶皺分析及模擬計(jì)算
    內(nèi)部熱集成精餾塔分離混合碳五的模擬研究
    實(shí)際發(fā)射工況下底排藥柱結(jié)構(gòu)完整性的模擬計(jì)算
    一種用于處理冷氫化料的低能耗精餾工藝
    化工管理(2014年12期)2014-03-16 02:03:12
    内射极品少妇av片p| 亚洲精品美女久久久久99蜜臀| 一区二区三区国产精品乱码| 国产一区二区三区在线臀色熟女| av天堂在线播放| 国产亚洲精品综合一区在线观看| 老司机午夜十八禁免费视频| 岛国在线免费视频观看| eeuss影院久久| 黄色成人免费大全| 久久伊人香网站| av视频在线观看入口| 美女黄网站色视频| 在线观看av片永久免费下载| 亚洲av五月六月丁香网| 怎么达到女性高潮| 久久草成人影院| 欧美性感艳星| 九九久久精品国产亚洲av麻豆| 亚洲一区二区三区色噜噜| 欧美乱码精品一区二区三区| 成人鲁丝片一二三区免费| 桃红色精品国产亚洲av| 欧美一区二区国产精品久久精品| www日本黄色视频网| 白带黄色成豆腐渣| 久久国产精品人妻蜜桃| 欧美乱妇无乱码| 亚洲精品色激情综合| 日韩成人在线观看一区二区三区| 在线a可以看的网站| av专区在线播放| 成熟少妇高潮喷水视频| xxxwww97欧美| a级毛片a级免费在线| 成人国产一区最新在线观看| 少妇高潮的动态图| 可以在线观看毛片的网站| 可以在线观看的亚洲视频| 久久久久久久精品吃奶| 国产成年人精品一区二区| 五月伊人婷婷丁香| 欧美+亚洲+日韩+国产| 国产亚洲精品一区二区www| 国产精品久久电影中文字幕| 18禁裸乳无遮挡免费网站照片| 久久久久久大精品| 亚洲自拍偷在线| 女同久久另类99精品国产91| 国产成+人综合+亚洲专区| 国产 一区 欧美 日韩| 青草久久国产| 欧美性猛交黑人性爽| 久久精品国产自在天天线| 伊人久久大香线蕉亚洲五| 熟妇人妻久久中文字幕3abv| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 麻豆成人午夜福利视频| 超碰av人人做人人爽久久 | 日本 av在线| 日韩人妻高清精品专区| 少妇高潮的动态图| 亚洲不卡免费看| 嫁个100分男人电影在线观看| 宅男免费午夜| 一个人看的www免费观看视频| 日本在线视频免费播放| 女人被狂操c到高潮| 免费观看的影片在线观看| 搡老妇女老女人老熟妇| 床上黄色一级片| 九九久久精品国产亚洲av麻豆| 亚洲精品日韩av片在线观看 | www.熟女人妻精品国产| 国产欧美日韩一区二区精品| 最后的刺客免费高清国语| 日韩人妻高清精品专区| 久久久久精品国产欧美久久久| 欧美zozozo另类| 成人亚洲精品av一区二区| 亚洲精品日韩av片在线观看 | www日本黄色视频网| 午夜福利在线在线| 亚洲电影在线观看av| 岛国在线观看网站| 欧美高清成人免费视频www| 国产日本99.免费观看| 夜夜爽天天搞| 国产一区二区亚洲精品在线观看| 国产野战对白在线观看| 国产精品99久久久久久久久| 尤物成人国产欧美一区二区三区| 国产97色在线日韩免费| 真人一进一出gif抽搐免费| 日日干狠狠操夜夜爽| 日本精品一区二区三区蜜桃| 男女下面进入的视频免费午夜| 免费在线观看影片大全网站| 日韩精品中文字幕看吧| 久久国产精品人妻蜜桃| 亚洲av电影不卡..在线观看| 99久久精品国产亚洲精品| 99久久精品国产亚洲精品| 久久天躁狠狠躁夜夜2o2o| 天天添夜夜摸| 欧美一区二区国产精品久久精品| 热99在线观看视频| 欧美在线黄色| 亚洲电影在线观看av| 国产精品99久久99久久久不卡| 一级毛片高清免费大全| 亚洲av不卡在线观看| 亚洲人成电影免费在线| 母亲3免费完整高清在线观看| 成人特级av手机在线观看| 五月玫瑰六月丁香| 国产欧美日韩一区二区三| 色噜噜av男人的天堂激情| 国产精品久久久久久久电影 | 3wmmmm亚洲av在线观看| 久久久久亚洲av毛片大全| 久久国产精品影院| 亚洲成av人片在线播放无| 黄色丝袜av网址大全| 欧美丝袜亚洲另类 | 欧美绝顶高潮抽搐喷水| 99精品在免费线老司机午夜| 两个人看的免费小视频| 国产精品永久免费网站| 国产精品嫩草影院av在线观看 | 久久久久久久久大av| 国产三级在线视频| 日韩欧美 国产精品| 1000部很黄的大片| 国产一区二区三区视频了| 伊人久久精品亚洲午夜| 久久久久免费精品人妻一区二区| 真实男女啪啪啪动态图| 国内精品久久久久精免费| 制服人妻中文乱码| 真实男女啪啪啪动态图| 淫妇啪啪啪对白视频| 亚洲人成网站高清观看| 看免费av毛片| 国产激情欧美一区二区| a级一级毛片免费在线观看| 一级毛片高清免费大全| 欧美日本视频| 一个人看视频在线观看www免费 | 男人舔奶头视频| 国内少妇人妻偷人精品xxx网站| 免费看a级黄色片| 欧美又色又爽又黄视频| 日韩欧美精品v在线| 亚洲精品日韩av片在线观看 | 香蕉丝袜av| 免费看美女性在线毛片视频| 精品免费久久久久久久清纯| 搡老妇女老女人老熟妇| 亚洲黑人精品在线| 午夜老司机福利剧场| 欧美成人a在线观看| 69人妻影院| 最近最新中文字幕大全免费视频| 91在线精品国自产拍蜜月 | 18禁裸乳无遮挡免费网站照片| 亚洲美女视频黄频| 久久久久亚洲av毛片大全| 日韩欧美 国产精品| 国产99白浆流出| 我的老师免费观看完整版| 999久久久精品免费观看国产| 欧美日韩综合久久久久久 | 97碰自拍视频| 亚洲精品亚洲一区二区| 人人妻人人澡欧美一区二区| 最新中文字幕久久久久| 久久精品影院6| 一二三四社区在线视频社区8| 国产黄a三级三级三级人| a级毛片a级免费在线| 黄色视频,在线免费观看| 亚洲av一区综合| 91麻豆精品激情在线观看国产| 黄色女人牲交| 久久久久久久久大av| 中文字幕人妻丝袜一区二区| 亚洲精品美女久久久久99蜜臀| 天天添夜夜摸| 丝袜美腿在线中文| 高清毛片免费观看视频网站| 久久亚洲真实| 日本a在线网址| 欧美xxxx黑人xx丫x性爽| 日韩高清综合在线| 精品一区二区三区视频在线 | 手机成人av网站| 手机成人av网站| 在线观看舔阴道视频| 久久久久国产精品人妻aⅴ院| h日本视频在线播放| 欧美三级亚洲精品| 九九在线视频观看精品| 国产老妇女一区| 国产精品久久久久久久电影 | 欧美成人a在线观看| 亚洲最大成人手机在线| 中文字幕精品亚洲无线码一区| 午夜激情欧美在线| 欧美3d第一页| 国产久久久一区二区三区| 制服人妻中文乱码| 波多野结衣巨乳人妻| 网址你懂的国产日韩在线| 人妻夜夜爽99麻豆av| 五月伊人婷婷丁香| 国产av不卡久久| 偷拍熟女少妇极品色| 亚洲精品美女久久久久99蜜臀| 淫秽高清视频在线观看| 观看美女的网站| 51午夜福利影视在线观看| 日韩欧美精品免费久久 | 夜夜爽天天搞| 亚洲欧美日韩高清专用| 又爽又黄无遮挡网站| 级片在线观看| 国产亚洲精品综合一区在线观看| 欧美激情久久久久久爽电影| 全区人妻精品视频| 午夜福利在线在线| 精品国产美女av久久久久小说| 男女床上黄色一级片免费看| 国产精品女同一区二区软件 | 又黄又爽又免费观看的视频| 黄色日韩在线| 欧美又色又爽又黄视频| 99久久精品热视频| 色综合站精品国产| 99久久无色码亚洲精品果冻| 1024手机看黄色片| 亚洲欧美激情综合另类| 中出人妻视频一区二区| 亚洲av成人精品一区久久| 女同久久另类99精品国产91| 亚洲专区中文字幕在线| 五月伊人婷婷丁香| 床上黄色一级片| 18美女黄网站色大片免费观看| 精品99又大又爽又粗少妇毛片 | 欧美成人性av电影在线观看| 亚洲欧美激情综合另类| 国产精品亚洲美女久久久| 欧美日韩黄片免| 欧美午夜高清在线| 在线观看午夜福利视频| 亚洲国产精品久久男人天堂| 搡老熟女国产l中国老女人| 国产一区二区在线观看日韩 | 国产成年人精品一区二区| 男人的好看免费观看在线视频| 国产69精品久久久久777片| 丰满人妻熟妇乱又伦精品不卡| 婷婷精品国产亚洲av| 一级黄片播放器| 久99久视频精品免费| 国产激情偷乱视频一区二区| 国产成人系列免费观看| 天天添夜夜摸| 成人av在线播放网站| 精品久久久久久成人av| 国产熟女xx| 亚洲一区二区三区色噜噜| 99久久精品热视频| 日本免费a在线| 99精品久久久久人妻精品| 亚洲天堂国产精品一区在线| 国产精品1区2区在线观看.| 免费看a级黄色片| 老汉色av国产亚洲站长工具| 黄片大片在线免费观看| 狂野欧美激情性xxxx| 国产精品电影一区二区三区| 熟女电影av网| 美女免费视频网站| 久久精品综合一区二区三区| 毛片女人毛片| 老司机福利观看| 久久精品国产综合久久久| 一个人看的www免费观看视频| 亚洲 国产 在线| 欧美黄色片欧美黄色片| 99久久无色码亚洲精品果冻| 成人鲁丝片一二三区免费| 欧美最黄视频在线播放免费| 精品不卡国产一区二区三区| 97人妻精品一区二区三区麻豆| 久久精品影院6| 日韩高清综合在线| 中出人妻视频一区二区| 亚洲中文日韩欧美视频| 男女午夜视频在线观看| 国产精品女同一区二区软件 | 国内少妇人妻偷人精品xxx网站| 听说在线观看完整版免费高清| 精品99又大又爽又粗少妇毛片 | bbb黄色大片| 国产亚洲精品综合一区在线观看| av中文乱码字幕在线| 久久久精品欧美日韩精品| 国产中年淑女户外野战色| 久久久久久久久中文| 一进一出抽搐gif免费好疼| 嫩草影院精品99| 国产美女午夜福利| 小蜜桃在线观看免费完整版高清| 老鸭窝网址在线观看| 国产精品亚洲一级av第二区| tocl精华| 日韩欧美在线二视频| 99久久99久久久精品蜜桃| 99久久久亚洲精品蜜臀av| 一进一出好大好爽视频| 久久6这里有精品| 黄色视频,在线免费观看| 最近在线观看免费完整版| 亚洲在线观看片| 中文字幕久久专区| 国产成人a区在线观看| 一级毛片女人18水好多| 亚洲专区国产一区二区| 男女床上黄色一级片免费看| 日本一二三区视频观看| 亚洲人与动物交配视频| 午夜福利成人在线免费观看| 脱女人内裤的视频| 亚洲一区二区三区不卡视频| 欧美日韩瑟瑟在线播放| 亚洲一区二区三区色噜噜| 神马国产精品三级电影在线观看| 90打野战视频偷拍视频| 国产麻豆成人av免费视频| 国产精品久久久久久久久免 | 国产精品99久久99久久久不卡| 午夜福利高清视频| 男女下面进入的视频免费午夜| 最近视频中文字幕2019在线8| 国产麻豆成人av免费视频| 一个人看视频在线观看www免费 | 又爽又黄无遮挡网站| 国产毛片a区久久久久| 国产欧美日韩一区二区三| 日本一二三区视频观看| 国产熟女xx| 国产欧美日韩一区二区三| 免费观看的影片在线观看| 成熟少妇高潮喷水视频| 国产伦在线观看视频一区| 在线观看av片永久免费下载| 男女床上黄色一级片免费看| 精品久久久久久成人av| 国产伦精品一区二区三区视频9 | 国产精品久久久久久精品电影| 一个人免费在线观看电影| 亚洲最大成人手机在线| 午夜福利成人在线免费观看| 一区二区三区国产精品乱码| 久久国产精品影院| 国产一区在线观看成人免费| 极品教师在线免费播放| 男人和女人高潮做爰伦理| 最近最新中文字幕大全电影3| 少妇的丰满在线观看| 日本三级黄在线观看| 色av中文字幕| 天堂av国产一区二区熟女人妻| 白带黄色成豆腐渣| 国产精品99久久99久久久不卡| 国语自产精品视频在线第100页| 一个人看视频在线观看www免费 | 草草在线视频免费看| 九九热线精品视视频播放| 中国美女看黄片| 99久久综合精品五月天人人| 中文在线观看免费www的网站| 又黄又粗又硬又大视频| 日韩免费av在线播放| 久久久久久久久久黄片| 亚洲精品影视一区二区三区av| 国产97色在线日韩免费| 精品国内亚洲2022精品成人| 一个人看的www免费观看视频| 97碰自拍视频| 亚洲18禁久久av| 国模一区二区三区四区视频| 久久草成人影院| 人人妻人人看人人澡| 国产极品精品免费视频能看的| 婷婷六月久久综合丁香| 成人午夜高清在线视频| 日本黄色片子视频| 亚洲国产日韩欧美精品在线观看 | 亚洲狠狠婷婷综合久久图片| 99精品久久久久人妻精品| 国产 一区 欧美 日韩| 国产黄片美女视频| 一级a爱片免费观看的视频| 国产精品三级大全| av女优亚洲男人天堂| 99久久99久久久精品蜜桃| 欧美日韩黄片免| 久久香蕉国产精品| 法律面前人人平等表现在哪些方面| 欧美最新免费一区二区三区 | 丰满人妻熟妇乱又伦精品不卡| 亚洲在线观看片| 欧美乱色亚洲激情| 老汉色∧v一级毛片| 丰满的人妻完整版| 久久亚洲真实| 精品午夜福利视频在线观看一区| 真人做人爱边吃奶动态| 精品人妻偷拍中文字幕| 午夜免费观看网址| 中国美女看黄片| 色尼玛亚洲综合影院| 国产精品98久久久久久宅男小说| av中文乱码字幕在线| 免费高清视频大片| 亚洲欧美一区二区三区黑人| 午夜激情福利司机影院| 国产成人av激情在线播放| h日本视频在线播放| 18禁黄网站禁片免费观看直播| 91久久精品国产一区二区成人 | 一本综合久久免费| 午夜久久久久精精品| 欧美乱码精品一区二区三区| 在线天堂最新版资源| 精品久久久久久,| 麻豆国产97在线/欧美| e午夜精品久久久久久久| 亚洲成人免费电影在线观看| 中文字幕久久专区| 亚洲av电影在线进入| 99热这里只有精品一区| www日本在线高清视频| 亚洲精品在线观看二区| 欧美性猛交黑人性爽| 国产成人aa在线观看| 小蜜桃在线观看免费完整版高清| av在线蜜桃| 俄罗斯特黄特色一大片| 国产在线精品亚洲第一网站| 一区二区三区激情视频| 久久久久久久久大av| 欧美激情在线99| 中文字幕人妻熟人妻熟丝袜美 | 国内少妇人妻偷人精品xxx网站| 亚洲av成人av| 国产中年淑女户外野战色| 搡老熟女国产l中国老女人| 国产野战对白在线观看| 天美传媒精品一区二区| 亚洲av不卡在线观看| 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 久久久久九九精品影院| 99久久精品一区二区三区| 欧美又色又爽又黄视频| 国产高清有码在线观看视频| 国产精品永久免费网站| 亚洲成人久久性| 久久久精品欧美日韩精品| 欧美中文综合在线视频| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 麻豆一二三区av精品| 亚洲成人中文字幕在线播放| 国产亚洲精品久久久com| 日日摸夜夜添夜夜添小说| 97人妻精品一区二区三区麻豆| 很黄的视频免费| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 国产在线精品亚洲第一网站| 亚洲人成网站在线播| 中文资源天堂在线| 日本 av在线| 欧美丝袜亚洲另类 | 欧美xxxx黑人xx丫x性爽| 在线观看免费午夜福利视频| 麻豆成人午夜福利视频| 国产 一区 欧美 日韩| 亚洲自拍偷在线| www.www免费av| 亚洲av中文字字幕乱码综合| 国产一区在线观看成人免费| 亚洲精品影视一区二区三区av| 99久久99久久久精品蜜桃| 婷婷六月久久综合丁香| 操出白浆在线播放| 99久久久亚洲精品蜜臀av| 亚洲成av人片免费观看| 亚洲av中文字字幕乱码综合| av国产免费在线观看| 99精品久久久久人妻精品| 少妇熟女aⅴ在线视频| 色吧在线观看| xxxwww97欧美| 国产99白浆流出| 久久亚洲真实| 午夜日韩欧美国产| 国产免费av片在线观看野外av| 国产精品嫩草影院av在线观看 | 搞女人的毛片| 他把我摸到了高潮在线观看| 别揉我奶头~嗯~啊~动态视频| 韩国av一区二区三区四区| 99国产极品粉嫩在线观看| 香蕉丝袜av| 国产成人av激情在线播放| 国产成人系列免费观看| 亚洲精华国产精华精| 国产午夜精品论理片| 国产伦精品一区二区三区视频9 | 亚洲精品影视一区二区三区av| 国产高清视频在线播放一区| 全区人妻精品视频| 动漫黄色视频在线观看| 国产高清视频在线观看网站| 99在线视频只有这里精品首页| 蜜桃亚洲精品一区二区三区| 动漫黄色视频在线观看| av天堂在线播放| 我要搜黄色片| 国产精品永久免费网站| 性欧美人与动物交配| 亚洲中文字幕一区二区三区有码在线看| 男女视频在线观看网站免费| 精品久久久久久成人av| www日本黄色视频网| 中亚洲国语对白在线视频| 国产毛片a区久久久久| 中文字幕熟女人妻在线| 亚洲第一电影网av| 国内少妇人妻偷人精品xxx网站| a在线观看视频网站| 精品熟女少妇八av免费久了| 欧美绝顶高潮抽搐喷水| 欧美黄色片欧美黄色片| 久9热在线精品视频| 亚洲人成网站在线播| www日本在线高清视频| 他把我摸到了高潮在线观看| 亚洲精华国产精华精| 国产亚洲精品av在线| 亚洲av免费在线观看| 亚洲五月天丁香| 窝窝影院91人妻| 亚洲专区中文字幕在线| 欧美日韩乱码在线| svipshipincom国产片| 国产伦一二天堂av在线观看| 丁香欧美五月| 高清毛片免费观看视频网站| 精品人妻一区二区三区麻豆 | 色综合欧美亚洲国产小说| 日本免费一区二区三区高清不卡| 亚洲国产欧美人成| 国产精品久久电影中文字幕| 午夜福利在线观看免费完整高清在 | 天堂网av新在线| 久久国产乱子伦精品免费另类| 国产欧美日韩精品亚洲av| 国产在视频线在精品| 精品人妻1区二区| 在线观看午夜福利视频| 99久久99久久久精品蜜桃| 国产av在哪里看| 免费看十八禁软件| 一级毛片高清免费大全| 一级a爱片免费观看的视频| 久久九九热精品免费| 久久久精品大字幕| 国产三级中文精品| 久久久久久久午夜电影| 午夜福利视频1000在线观看| 亚洲av成人精品一区久久| 91字幕亚洲| 99国产精品一区二区蜜桃av| 18禁美女被吸乳视频| 免费在线观看日本一区| 国语自产精品视频在线第100页| 乱人视频在线观看| 精品国内亚洲2022精品成人| 一进一出抽搐gif免费好疼| 国产精品一及| 亚洲成av人片免费观看| 欧美一级a爱片免费观看看| 久久欧美精品欧美久久欧美| 观看美女的网站| 小蜜桃在线观看免费完整版高清| 长腿黑丝高跟| 美女cb高潮喷水在线观看| 熟妇人妻久久中文字幕3abv| 麻豆久久精品国产亚洲av| 男女视频在线观看网站免费| 国产综合懂色| 久久久久久人人人人人| 欧美乱妇无乱码| 一个人看视频在线观看www免费 | 激情在线观看视频在线高清| 最近视频中文字幕2019在线8| 精品久久久久久久久久免费视频|