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

    分離混合二甲苯的對二甲苯-間二甲苯混合生產(chǎn)工藝的模擬

    2015-06-06 11:54:03劉肖肖許新望王麗軍成有為
    石油化工 2015年2期

    劉肖肖,許新望,王麗軍,成有為,李 希

    (浙江大學(xué) 化學(xué)工程與生物工程學(xué)系,浙江 杭州 310027)

    分離混合二甲苯的對二甲苯-間二甲苯混合生產(chǎn)工藝的模擬

    劉肖肖,許新望,王麗軍,成有為,李 希

    (浙江大學(xué) 化學(xué)工程與生物工程學(xué)系,浙江 杭州 310027)

    在結(jié)晶熱力學(xué)分析的基礎(chǔ)上,提出用于分離混合二甲苯的對二甲苯-間二甲苯(PX-MX)混合生產(chǎn)工藝,并建立數(shù)學(xué)模型進(jìn)行模擬??疾炝司s塔理論塔板數(shù)、精餾塔頂采出量與進(jìn)料量的摩爾比(精餾塔采出比)和結(jié)晶溫度對PX-MX混合產(chǎn)品回收率(YMX/PX)及產(chǎn)品中MX與PX摩爾比的影響,并對4種混合C8芳烴原料進(jìn)行全流程模擬。模擬結(jié)果表明,在精餾塔理論塔板數(shù)為50、精餾塔采出比為0.8、結(jié)晶溫度為高于各進(jìn)料條件下PX-MX-OX三組分共晶溫度5 K時,4種混合C8芳烴為原料通過PX-MX混合生產(chǎn)工藝分離,均可獲得含量大于99.9%(x)的PX-MX混合產(chǎn)品,YMX/PX最高可達(dá)61.55%,PX回收率最高可達(dá)84.77%;同時省去MX單組分提純過程,可降低生產(chǎn)成本。

    過程模擬;精餾-結(jié)晶;對二甲苯;間二甲苯;混合二甲苯

    對苯二甲酸(PTA)和間苯二甲酸(PIA)為重要的化工基礎(chǔ)原料,主要用于生產(chǎn)聚酯纖維和瓶片[1-2]。目前,PTA和PIA均采用高純度的單體原料通過分別氧化獲得。國內(nèi)外廣泛采用吸附法制備高純度的對二甲苯(PX)和間二甲苯(MX),但因?yàn)槲絼┑膿p耗及更換使得生產(chǎn)成本較高,對于生產(chǎn)規(guī)模較小的MX裝置顯得尤為不經(jīng)濟(jì)[3-4]。聚酯生產(chǎn)過程中PIA的主要用途是與PTA進(jìn)行共聚,使用時要求二者充分混合[5-6]。因此,如果能直接利用混合二甲苯原料而不是其高純度單組分物質(zhì)進(jìn)行氧化,得到混合苯二羧酸,則生產(chǎn)工藝可大幅簡化,生產(chǎn)成本也會顯著降低。張永昭[7]提出了一種PX-MX的共氧化技術(shù),研究結(jié)果表明,使用PX-MX混合原料生產(chǎn)PTA-PIA混合苯二羧酸是完全可行的。

    基于PX-MX共氧化考慮,提出一種新的混合二甲苯純化思路,即不以單組分產(chǎn)品為目標(biāo),而是以PX-MX混合物為最終產(chǎn)品,省去PX與MX的分離過程。由于結(jié)晶分離法操作彈性更大,通過調(diào)整結(jié)晶溫度即可獲得單一組分產(chǎn)品,又可得到固定組成的混合物。同時精餾分離可使結(jié)晶器進(jìn)料遠(yuǎn)離最低共熔點(diǎn)組成,提高PX-MX混合產(chǎn)品收率。因此將精餾分離與結(jié)晶分離相結(jié)合,可從混合二甲苯中經(jīng)濟(jì)地分離出PX-MX混合物,從而解決苯二羧酸原料難以廉價獲得的難題。

    本工作在結(jié)晶熱力學(xué)分析的基礎(chǔ)上,提出用于分離混合二甲苯的PX-MX混合生產(chǎn)工藝,并建立數(shù)學(xué)模型進(jìn)行模擬,考察精餾塔理論塔板數(shù)、精餾塔頂采出量與進(jìn)料量的摩爾比(精餾塔采出比)和結(jié)晶溫度對分離效果的影響。

    1 混合二甲苯物系的固液相平衡

    1.1 結(jié)晶熱力學(xué)計(jì)算

    混合二甲苯主要包括MX、PX、鄰二甲苯(OX)、乙苯(EB),各物質(zhì)的物性參數(shù)見表1[8-9]。

    由表1可知,混合二甲苯不同組分間的最大沸點(diǎn)差僅為5.25 K,采用單一的精餾方法分離混合二甲苯很困難。而不同二甲苯之間的熔點(diǎn)差比較大,因而可采用結(jié)晶法分離混合二甲苯?;旌隙妆街械?個二組分物系均為簡單低共熔物系,其理想溶液模型固液相平衡基本關(guān)系式可表示為van’t Hoff方程[8,10-13],見式(1)。

    表1 混合二甲苯中各組分的物性參數(shù)Table 1 Properties of the components in mixed xylenes

    混合二甲苯結(jié)晶過程中的最大回收率受到低共熔點(diǎn)組成的限制,以四組分混合二甲苯平衡組成為基礎(chǔ),分離其中的單組分PX最大回收率為65%[14]。采用式(1)計(jì)算得到的混合二甲苯中二組分、三組分、四組分共熔點(diǎn)溫度及組成及文獻(xiàn)值[10,14]見表2。

    表2 混合二甲苯物系的低共熔點(diǎn)溫度及組成Table 2 Eutectic point temperature and composition of mixed xylenes systems

    從表2可看出,式(1)的計(jì)算結(jié)果與文獻(xiàn)值符合良好。因此,四組分混合二甲苯物系結(jié)晶過程的固液相平衡數(shù)據(jù)采用式(1)計(jì)算,并采用該模型預(yù)測四組分混合二甲苯物系的析晶過程。

    1.2 析晶過程分析

    混合二甲苯物系的析晶過程與溫度和原料組成有關(guān)。隨溫度的降低,混合二甲苯物系結(jié)晶析出的先后順序分別涉及初始結(jié)晶、二組分固熔面結(jié)晶、三組分固熔線結(jié)晶以及四組分最低共熔點(diǎn)結(jié)晶過程。

    混合二甲苯四組分固液相圖見圖1。該相圖為正四面體,4個頂點(diǎn)為4個純組分,6條邊為6個二組分物系,4個正三角形為4個三組分物系,正四面體內(nèi)部為四組分物系。以典型催化重整獲得的混合二甲苯為原料(摩爾組成為:MX 45%,PX 20%,OX 20%,EB 15%)對析晶過程進(jìn)行分析,進(jìn)料組成M落在組分PX初始結(jié)晶區(qū)域內(nèi),溫度降低時PX單晶體首先析出,隨PX析出其濃度逐漸降低,而MX,OX,EB濃度相應(yīng)增加,液相組成沿PX-M延長線變化,至L點(diǎn)。L點(diǎn)為PX-MX固熔面e1e2e3E上的點(diǎn),溫度進(jìn)一步降低時PX-MX同時析出,各組分濃度變化在固熔面e1e2e3E內(nèi)。當(dāng)達(dá)到三組分固熔線PE上的P點(diǎn)時,OX隨MX和PX析出。溫度進(jìn)一步降低,液相濃度沿PE變化,至四組分最低共熔E點(diǎn),達(dá)五相固液平衡狀態(tài),液相組成不再變化,直至析晶結(jié)束[15]。

    由圖1可知,只要合理控制結(jié)晶原料的組成和溫度,就可獲得PX-MX混合結(jié)晶產(chǎn)物,且保證OX 和EB不析出。

    2 PX-MX混合生產(chǎn)工藝及模擬

    2.1 流程描述

    PX-MX混合生產(chǎn)工藝的流程見圖2,主要由精餾、結(jié)晶、晶體提純3個單元組成?;旌螩8芳烴原料首先進(jìn)入精餾塔進(jìn)行精餾分離,在塔底得到富含OX的物流,在塔頂?shù)玫礁患疨X和MX的物流。精餾塔頂物流進(jìn)入結(jié)晶器進(jìn)行結(jié)晶,PX和MX混合晶體從液相中析出。漿料在固液分離器中通過離心分離獲得PX和MX混合晶體,混合晶體進(jìn)入逆流晶體洗滌塔,通過熔融液與晶體的逆流接觸去除雜質(zhì)[16-17]。

    混合C8芳烴的組成主要為二甲苯,其來源包括催化重整油、裂解汽油、甲苯歧化和煤焦油。不同來源的混合C8芳烴的組成見表3[18]。

    圖1 混合二甲苯四組分固液相圖Fig.1 Phase diagram for the crystallization of mixed xylenes.e1:eutectic point of PX-MX;e2:eutectic point of PX-MX-OX;e3:eutectic point of PX-MX-EB;E:eutectic point of PX-MX-OX-EB;M:feed composition;L:point on the solid fusing surface of PX-MX;P:point on the solid fusing line of PX-MX-OX.

    圖2 PX-MX混合生產(chǎn)工藝流程Fig.2 Process for the production of mixed PX and MX.

    表3 不同來源的混合C8芳烴的組成Table 3 Compositions of mixed C8aromatics from different sources

    2.2 數(shù)學(xué)模型

    利用數(shù)學(xué)工具M(jìn)atlab對PX-MX混合生產(chǎn)工藝流程建立數(shù)學(xué)模型。

    2.2.1 精餾塔模型

    精餾塔模型采用MESH方程組[19]:

    物料衡算方程

    相平衡方程

    歸一化方程

    能量衡算方程

    利用Wang-Henke法求解精餾塔模型,Newton-Raphson迭代過程采用高斯消元法求解。

    2.2.2 結(jié)晶器模型

    結(jié)晶器是PX-MX混合生產(chǎn)工藝流程中的核心設(shè)備。圖3為單級結(jié)晶理論板的示意圖。結(jié)晶器模型主要包括兩部分:首先按照理想溶液模型計(jì)算給定進(jìn)料組成下,i組分的結(jié)晶溫度(i=1,2,…,N);然后計(jì)算在給定結(jié)晶溫度和進(jìn)料條件下的固液平衡,得出固體量、母液量及相應(yīng)組成。

    圖3 單級結(jié)晶理論板的示意圖Fig.3 Schematic diagram of crystallizing stage.Fc:feed;S:solid at solid-liquid equilibrium;Lc:liquid attached to solid;W:mother liquid.

    結(jié)晶器模型包括物料平衡方程、固液相平衡約束式和歸一化方程等:

    物料平衡方程

    結(jié)晶的組分固液相平衡約束式

    非結(jié)晶的組分固液相平衡約束式

    歸一化方程

    其他約束方程

    結(jié)晶器模型的求解步驟見圖4。

    圖4 結(jié)晶器模型的求解步驟Fig.4 Flow diagram of calculation for crystallizer.Tm:melting point;Δ hm:enthalpy of fusion;x0:feed composition;Tc:crystallization temperature.

    2.2.3 逆流晶體洗滌塔

    賴家鳳等[20]認(rèn)為,通過降溫結(jié)晶形成的固液非均相體系的分離純化過程可類似于精餾塔的精餾過程,以熔融結(jié)晶洗滌塔內(nèi)的液固相平衡類比精餾塔內(nèi)的汽液相平衡。因此,以連續(xù)多級逆流結(jié)晶器模型來模擬多級固液相平衡,每一級均調(diào)用結(jié)晶器模型進(jìn)行固液相平衡的計(jì)算,模型求解過程與精餾塔類似。

    3 結(jié)果與討論

    以催化重整油的混合C8芳烴原料為例,考察精餾塔理論塔板數(shù)、精餾塔采出比、結(jié)晶器溫度對PX-MX混合產(chǎn)品回收率(YMX/PX)及產(chǎn)品中MX與PX摩爾比(rMX/PX)的影響。YMX/PX定義為經(jīng)提純后的產(chǎn)品中總的PX和MX與進(jìn)料中總的PX和MX的摩爾比。PX-MX混合生產(chǎn)工藝相關(guān)參數(shù)的基準(zhǔn)值見表4。

    表4 PX-MX混合生產(chǎn)工藝相關(guān)參數(shù)的基準(zhǔn)值Table 4 Reference values of parameters in the process for the production of mixed PX and MX

    3.1 精餾塔理論塔板數(shù)的影響

    精餾塔理論塔板數(shù)對YMX/PX和rMX/PX的影響見圖5。由圖5可見,隨精餾塔理論塔板數(shù)從10塊增至50塊,YMX/PX和rMX/PX均增大;當(dāng)精餾塔理論塔板數(shù)達(dá)50塊以上時,YMX/PX和rMX/PX趨于穩(wěn)定。這是因?yàn)殡S精餾塔理論塔板數(shù)的增大,對OX切割的效果變好,在采出量不變的情況下塔頂物流中PX-MX含量增加,使得進(jìn)入結(jié)晶器的物料組成逐漸遠(yuǎn)離PXMX-OX三組分共晶點(diǎn)的組成,當(dāng)結(jié)晶溫度不變時YMX/PX和rMX/PX均增大;當(dāng)精餾塔理論塔板數(shù)超過50塊,在其他條件不變時繼續(xù)增大理論塔板數(shù),逐漸接近分離的極限,對精餾分離效果影響不大。

    圖5 精餾塔理論塔板數(shù)對YMX/PX和rMX/PX的影響Fig.5 Effects of number of theoretical plates of distillation column on YMX/PXand rMX/PX.YMX/PX:mole ratio of MX and PX in products to MX and PX in feed;rMX/PX:mole ratio of MX to PX in products.

    3.2 精餾塔采出比的影響

    精餾塔采出比對YMX/PX和rMX/PX的影響見圖6。由圖6可以看出,隨精餾塔采出比的增大,YMX/PX和rMX/PX均先增大后減小,精餾塔采出比為0.85時達(dá)到最大。這是因?yàn)?,在其他條件不變時,當(dāng)精餾塔采出比較小時,隨精餾塔采出比的增大,一方面精餾塔對OX的切割可達(dá)到要求,另一方面塔頂采出量增大,即結(jié)晶器的進(jìn)料量增大,二者共同作用使得YMX/PX和rMX/PX增大。當(dāng)精餾塔采出比增大至0.85時,繼續(xù)增大精餾塔采出比,由于精餾塔頂過度采出,使得對OX的切割達(dá)不到要求,進(jìn)而使YMX/PX和rMX/PX降低。

    圖6 精餾塔采出比對YMX/PX和rMX/PX的影響Fig.6 Effects of mole ratio of distillate to feed in distillation column on YMX/PXand rMX/PX.

    3.3 結(jié)晶溫度對析晶過程的影響

    結(jié)晶溫度對析晶過程的影響見圖7。由圖7可見,當(dāng)結(jié)晶溫度高于208 K時,結(jié)晶產(chǎn)物為PX單組分;結(jié)晶溫度處于208~191.6 K時,結(jié)晶產(chǎn)物為PX-MX二組分;結(jié)晶溫度低于191.6 K時,得到含有OX和EB的三組分或四組分結(jié)晶物。為獲得較多的PX和MX物質(zhì),適宜的結(jié)晶溫度應(yīng)略高于191.6 K。

    圖7 結(jié)晶溫度對析晶過程的影響Fig.7 Effect of crystallizing temperature on the number of crystals.

    結(jié)晶溫度對YMX/PX和rMX/PX的影響見圖8。由圖8可見,隨結(jié)晶溫度的降低,YMX/PX和rMX/PX均呈現(xiàn)逐漸增大的趨勢;當(dāng)結(jié)晶溫度為208 K時,由于第二結(jié)晶物質(zhì)MX的析出出現(xiàn)拐點(diǎn)。這可由圖1固液相圖進(jìn)行解釋,隨溫度的降低,更多的PX-MX結(jié)晶析出,且MX析晶量大于PX,故rMX/PX逐漸增大。

    圖8 結(jié)晶溫度對YMX/PX和rMX/PX的影響Fig.8 Effects of crystallizing temperature on YMX/PXand rMX/PX.

    3.4 全流程的模擬結(jié)果

    對4種混合C8芳烴原料,在精餾塔理論塔板數(shù)為50、采出比為0.8、結(jié)晶溫度高于相應(yīng)進(jìn)料條件下PX-MX-OX三元共晶溫度5 K且不低于194.3 K[21]時,進(jìn)行PX-MX混合生產(chǎn)工藝的全流程模擬,模擬結(jié)果見表5。由表5可看出,PX的回收率最低可達(dá)到77.45%,突破了單一組分結(jié)晶制備PX時由于共熔點(diǎn)造成的PX產(chǎn)品回收率為65%的限制。同時結(jié)晶制備PX-MX混合產(chǎn)品時,省去MX的單組分提純過程,可將原料混合二甲苯中的MX與PX一起進(jìn)行回收,可降低生產(chǎn)成本。

    表5 PX-MX混合生產(chǎn)全流程計(jì)算結(jié)果Table 5 Calculation for the production process of mixed MX-PX

    4 結(jié)論

    1)四組分混合二甲苯物系結(jié)晶過程的固液相平衡數(shù)據(jù)可采用van’t Hoff 方程計(jì)算,通過控制結(jié)晶溫度和原料組成可得到PX-MX二組分混合物。

    2)針對四組分混合二甲苯物系的分離,提出包括精餾單元、結(jié)晶單元和晶體提純單元的PX-MX混合生產(chǎn)工藝流程。建立數(shù)學(xué)模型對該工藝進(jìn)行模擬計(jì)算。

    3)分別對4種混合C8芳烴原料進(jìn)行流程模擬,在精餾塔理論塔板數(shù)為50、精餾塔采出比為0.8、結(jié)晶溫度不低于194.3 K時,均可得到純度大于99.9%(x)的混合PX-MX產(chǎn)品,且PX的回收率突破了單一組分結(jié)晶制備PX時回收率為65%的限制,最低可達(dá)到77.45%。

    符 號 說 明

    F 進(jìn)料摩爾流量,kmol/h

    H 焓值,kJ/kmol

    Δ hm標(biāo)準(zhǔn)摩爾生成焓,kJ/kmol

    L 液相流量,kmol/h

    Lc晶體附著液量,kmol/h

    L0結(jié)晶器固液平衡液相量,kmol/h

    p 壓力,MPa

    R 氣體常數(shù),J/(mol·K)

    S 結(jié)晶器固液平衡固相量,kmol/h

    T 溫度,K

    Tm標(biāo)準(zhǔn)熔點(diǎn),K

    V 汽相流量,kmol/h

    W 結(jié)晶器母液量,kmol/h

    x 液相摩爾分?jǐn)?shù)

    y 汽相摩爾分?jǐn)?shù)

    z 精餾塔進(jìn)料摩爾分?jǐn)?shù)

    φ 逸度系數(shù)

    γ 活度系數(shù)

    上角標(biāo)

    L 液相

    V 汽相

    下角標(biāo)

    c 結(jié)晶器

    d 精餾塔

    F 進(jìn)料

    i 組分,i=1,2,…,N

    l 精餾塔塔板位置

    [1] 陳亮,肖劍,謝在庫. 對二甲苯結(jié)晶分離技術(shù)進(jìn)展[J]. 現(xiàn)代化工,2009,29(2):10 - 14.

    [2] 嚴(yán)明,張曉東. 影響間二甲苯工藝的因素分析及調(diào)整措施[J]. 化學(xué)工業(yè),2009,27(6):43 - 45.

    [3] Bárcia P S,Nicolau M P M,Gallegos J M,et al. Modeling Adsorption Equilibria of Xylene Isomers in a Microporous Metal-Organic Framework[J]. Microporous Mesoporous Mater,2012,155:220-226.

    [4] Marta S P S,José P B M,Alírio E R. Fixed-Bed Adsorption of Aromatic C8Isomers:Breakthrough Experiments,Modeling and Simulation[J]. Sep Purif Technol,2012,90:246 - 256.

    [5] 林生兵,姚峰,瞿中凱. 低熔點(diǎn)聚酯的合成與性能研究[J].合成纖維工業(yè),2005,28(2):13 - 16.

    [6] Rasouli M,Yaghobi N,Chitsazan S,et al. Adsorptive Separation of Meta-Xylene from C8Aromatics[J]. Chem Eng Res Des,2012,90(9):1407 - 1415.

    [7] 張永昭. 間二甲苯氧化動力學(xué)與混合氧化新工藝研究[D]. 杭州:浙江大學(xué),2006.

    [8] Jakob A,Joh R,Rose C,et al. Solid-Liquid Equilibria in Binary Mixtures of Organic Compounds[J]. Fluid Phase Equilib,1995,113(1/2):117 - 126.

    [9] 趙紅玲,王鳳坤,陳圣坤,等. 氣液物性估算手冊[M]. 5版.北京:化學(xué)工業(yè)出版社,2004:515 - 516.

    [10] 陳亮, 肖劍, 謝在庫. 對二甲苯結(jié)晶過程的固液相平衡研究[J]. 聚酯工業(yè),2009,22(1):7 - 11.

    [11] Roger S,Julian F. Extended Xylene Solubility Studies[J]. J Chem Eng Data,1967,12(3):392 - 394.

    [12] William F,Julian F. Solubility Data for p-Xylene[J]. J Chem Eng Data,1964,9(1):158 - 159.

    [13] Daniel B P,James B R,Hazim A M. Crystallization of para-Xylene in Scraped-Surface Crystallizers[J]. AIChE J,2001,47(11):2441 - 2451.

    [14] 趙仁殿,金彰禮,陶志華,等. 芳烴工學(xué)[M]. 北京:化學(xué)工業(yè)出版社,2001:239 - 244.

    [15] 陳樹江,田鳳仁,李國華. 相圖分析及應(yīng)用[M]. 冶金工業(yè)出版社,2007:91 - 96.

    [16] 李群生,張遠(yuǎn)謀. 低共熔型物系連續(xù)多級逆流分步結(jié)晶提純塔的研究[J]. 化工學(xué)報(bào),1993,44(4):472 - 479.

    [17] 李群生,張遠(yuǎn)謀. 連續(xù)多級逆流分步結(jié)晶的研究[J]. 北京化工大學(xué)學(xué)報(bào),1997,24(1):1 - 8.

    [18] 成有為. MC工藝烴類液相催化氧化過程的研究[D]. 杭州:浙江大學(xué),2004.

    [19] 于志家,趙宗昌,王寶和,等. 化工過程物料平衡與能量平衡[M]. 大連:大連理工大學(xué)出版社,2008:167 - 175.

    [20] 賴家鳳,宗弘元,徐瑤,等. 熔融結(jié)晶洗滌分離過程的傳遞行為[J]. 過濾與分離,2007,17(4):15 - 18.

    [21] ExxonMobil Chemical Company. Process and Apparatus for para-Xylene Production:US,20080249345 A1[P]. 2008-10-09.

    (編輯 李治泉)

    Simulation of the Production Process of Mixed p-Xylene and m-Xylene from Mixed Xylenes

    Liu Xiaoxiao,Xu Xinwang,Wang Lijun,Cheng Youwei,Li Xi
    (Department of Chemical and Biochemcal Engineering,Zhejiang University,Hangzhou Zhejiang 310027,China)

    A production process for mixed p-xylene(PX) and m-xylene(MX) was proposed and simulated based on crystallization thermodynamics. The effects of number of theoretical plates,mole ratio of distillate to feed in distillation column and crystallization temperature on the yield of the mixed PX-MX(YMX/PX) and ratio of MX/PX were investigated. The process simulation was carried out for 4 kinds of C8raw materials. The results showed that,under the conditions of number of theoretical plates of distillation column 50,mole ratio of distillate to feed 0.8 and crystallization temperature 5 K higher than the eutectic temperature of PX-MX-OX,the purity of the mixed PX and MX products could reach more than 99.9%(x) for the C8raw materials,and both YMX/PXand the yield of PX in the products could reach the maximum,61.55% and 84.77%,respectively. Furthermore,the recovery of MX from mixed xylenes can be realized,which omits the MX purification process to reduce the production cost.

    process simulation;distillation-melt crystallization;p-xylene;m-xylene;mixed xylenes

    1000 - 8144(2015)02 - 0168 - 07

    TQ 021.8

    A

    2014 - 08 - 19;[修改稿日期] 2014 - 11 - 03。

    劉肖肖(1990—),女,河南省駐馬店市人,碩士生,電話 0571 - 87952210, 電郵 21228060@zju.edu.cn。聯(lián)系人:王麗軍,電話 0571 - 87952210,電郵 wang_lijun@zju.edu.cn。

    日本免费a在线| 久久精品亚洲精品国产色婷小说| 国产精品 欧美亚洲| 国产高清videossex| 欧美中文综合在线视频| 久久久久精品国产欧美久久久| 中文在线观看免费www的网站| 特大巨黑吊av在线直播| 国产成人一区二区三区免费视频网站| 久久午夜亚洲精品久久| 欧美乱妇无乱码| 一本久久中文字幕| 禁无遮挡网站| 国产亚洲av嫩草精品影院| 午夜福利在线观看免费完整高清在 | 琪琪午夜伦伦电影理论片6080| 精品国产美女av久久久久小说| 草草在线视频免费看| 最新中文字幕久久久久 | 长腿黑丝高跟| www.999成人在线观看| 亚洲国产欧美人成| 噜噜噜噜噜久久久久久91| 亚洲九九香蕉| 亚洲狠狠婷婷综合久久图片| 日韩欧美免费精品| 久久久国产精品麻豆| 午夜亚洲福利在线播放| 国产精品久久久久久人妻精品电影| 性色avwww在线观看| 99在线视频只有这里精品首页| 国产精品美女特级片免费视频播放器 | 亚洲avbb在线观看| 特大巨黑吊av在线直播| 中文字幕最新亚洲高清| 精品国产乱子伦一区二区三区| 久久久久久久久免费视频了| 国产探花在线观看一区二区| 国产伦精品一区二区三区四那| 欧美高清成人免费视频www| 国产1区2区3区精品| 好看av亚洲va欧美ⅴa在| 色综合婷婷激情| 成人高潮视频无遮挡免费网站| 免费无遮挡裸体视频| 亚洲国产高清在线一区二区三| 亚洲 欧美 日韩 在线 免费| 久久久久久人人人人人| 日本五十路高清| 美女免费视频网站| 在线观看一区二区三区| 五月玫瑰六月丁香| 国产高清视频在线观看网站| a级毛片a级免费在线| 黄色视频,在线免费观看| 99久久无色码亚洲精品果冻| 欧美成人一区二区免费高清观看 | 91九色精品人成在线观看| 免费看光身美女| 成年版毛片免费区| 亚洲人成网站在线播放欧美日韩| 免费在线观看影片大全网站| 18禁观看日本| 国产三级在线视频| 一a级毛片在线观看| 亚洲,欧美精品.| 90打野战视频偷拍视频| 一个人免费在线观看电影 | 成人欧美大片| 国产激情偷乱视频一区二区| 亚洲自偷自拍图片 自拍| 国产毛片a区久久久久| 成人特级黄色片久久久久久久| 精品国产亚洲在线| 日韩免费av在线播放| 一级a爱片免费观看的视频| 他把我摸到了高潮在线观看| 久久久久久久久久黄片| 久久久国产精品麻豆| 亚洲七黄色美女视频| 禁无遮挡网站| 亚洲人成伊人成综合网2020| 99久国产av精品| 天堂影院成人在线观看| 国产欧美日韩精品一区二区| h日本视频在线播放| ponron亚洲| 在线视频色国产色| 亚洲国产欧洲综合997久久,| 我要搜黄色片| 久久天躁狠狠躁夜夜2o2o| 给我免费播放毛片高清在线观看| 精品欧美国产一区二区三| 国产日本99.免费观看| 性欧美人与动物交配| 蜜桃久久精品国产亚洲av| 欧美黄色淫秽网站| 亚洲天堂国产精品一区在线| 国产熟女xx| 国产成+人综合+亚洲专区| 麻豆成人av在线观看| 久久精品91蜜桃| 精品不卡国产一区二区三区| 老熟妇仑乱视频hdxx| 亚洲国产看品久久| 午夜a级毛片| 在线免费观看的www视频| 久久久国产精品麻豆| 一个人观看的视频www高清免费观看 | 黄色视频,在线免费观看| 欧美激情在线99| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品成人综合色| 久久性视频一级片| 午夜激情欧美在线| 国产成人啪精品午夜网站| 可以在线观看毛片的网站| 在线十欧美十亚洲十日本专区| 国产99白浆流出| 日本与韩国留学比较| 亚洲熟妇中文字幕五十中出| 日韩欧美一区二区三区在线观看| 午夜免费观看网址| 久久久久久久午夜电影| 亚洲欧美日韩卡通动漫| 精品一区二区三区视频在线 | 亚洲av电影不卡..在线观看| 亚洲av成人精品一区久久| 欧美激情久久久久久爽电影| 身体一侧抽搐| 女人高潮潮喷娇喘18禁视频| 日本精品一区二区三区蜜桃| 男女之事视频高清在线观看| 欧美成人一区二区免费高清观看 | 国产亚洲av高清不卡| 国产99白浆流出| 日韩欧美精品v在线| 变态另类丝袜制服| 午夜福利在线观看吧| 亚洲成人久久爱视频| 真实男女啪啪啪动态图| 两性午夜刺激爽爽歪歪视频在线观看| 婷婷亚洲欧美| 久久久久久九九精品二区国产| 国产一区在线观看成人免费| 亚洲熟妇熟女久久| 俄罗斯特黄特色一大片| 久久久久国内视频| 欧美中文日本在线观看视频| 天堂网av新在线| 日韩精品中文字幕看吧| av中文乱码字幕在线| 久久久久久久久久黄片| 国产成人精品久久二区二区91| 久久久色成人| 欧美日韩亚洲国产一区二区在线观看| 国产激情偷乱视频一区二区| 国产精品1区2区在线观看.| 男女午夜视频在线观看| 日本精品一区二区三区蜜桃| 亚洲aⅴ乱码一区二区在线播放| 日韩有码中文字幕| av黄色大香蕉| 亚洲av成人一区二区三| 不卡一级毛片| 18禁黄网站禁片免费观看直播| 岛国视频午夜一区免费看| 亚洲av熟女| 午夜激情欧美在线| 天堂影院成人在线观看| 特级一级黄色大片| 成人av一区二区三区在线看| 精品国产三级普通话版| 亚洲av电影在线进入| 母亲3免费完整高清在线观看| 久久久色成人| 精品久久久久久成人av| 亚洲精品美女久久av网站| 午夜影院日韩av| 欧美黑人欧美精品刺激| 51午夜福利影视在线观看| 热99re8久久精品国产| 亚洲天堂国产精品一区在线| 国产伦一二天堂av在线观看| 老汉色∧v一级毛片| 美女免费视频网站| 日韩av在线大香蕉| 最新在线观看一区二区三区| 国产欧美日韩精品一区二区| 黑人巨大精品欧美一区二区mp4| 午夜影院日韩av| 中文字幕熟女人妻在线| 午夜精品久久久久久毛片777| www.熟女人妻精品国产| 国内精品一区二区在线观看| 日本免费一区二区三区高清不卡| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕熟女人妻在线| 最近在线观看免费完整版| 蜜桃久久精品国产亚洲av| 国产亚洲精品一区二区www| 999久久久国产精品视频| 国产精品免费一区二区三区在线| 免费大片18禁| 在线看三级毛片| 亚洲av片天天在线观看| 中文亚洲av片在线观看爽| 免费av毛片视频| 观看美女的网站| www国产在线视频色| 一级毛片高清免费大全| 亚洲中文字幕日韩| 一a级毛片在线观看| 国内毛片毛片毛片毛片毛片| 很黄的视频免费| 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影 | 精华霜和精华液先用哪个| 极品教师在线免费播放| 精品无人区乱码1区二区| 国产激情欧美一区二区| 亚洲中文字幕日韩| 香蕉av资源在线| 日日干狠狠操夜夜爽| 婷婷丁香在线五月| 久久久久免费精品人妻一区二区| 欧美一区二区精品小视频在线| 十八禁网站免费在线| 香蕉国产在线看| 亚洲人成网站高清观看| 亚洲国产欧美一区二区综合| 欧美丝袜亚洲另类 | 国产免费男女视频| 免费看a级黄色片| 琪琪午夜伦伦电影理论片6080| 麻豆成人午夜福利视频| 国产aⅴ精品一区二区三区波| 精品福利观看| 国内揄拍国产精品人妻在线| 国产精品久久久久久人妻精品电影| 免费看日本二区| 亚洲va日本ⅴa欧美va伊人久久| 日本与韩国留学比较| 在线观看美女被高潮喷水网站 | 国产精品一区二区三区四区免费观看 | 亚洲国产精品成人综合色| 99热6这里只有精品| 一区二区三区国产精品乱码| 亚洲成人久久性| 一二三四社区在线视频社区8| 99riav亚洲国产免费| 黄色丝袜av网址大全| 又紧又爽又黄一区二区| 欧美在线黄色| 欧美成人性av电影在线观看| 成熟少妇高潮喷水视频| 欧美乱色亚洲激情| 人人妻,人人澡人人爽秒播| 亚洲精品国产精品久久久不卡| av女优亚洲男人天堂 | 国产不卡一卡二| 夜夜爽天天搞| 午夜久久久久精精品| 麻豆成人av在线观看| 五月玫瑰六月丁香| 日韩有码中文字幕| 国产亚洲精品综合一区在线观看| 国产精品影院久久| 精品国产乱码久久久久久男人| 免费看美女性在线毛片视频| 精品国产乱码久久久久久男人| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 精品久久久久久久毛片微露脸| 中文字幕高清在线视频| 免费观看精品视频网站| 在线观看舔阴道视频| 久久欧美精品欧美久久欧美| 国产精品日韩av在线免费观看| 午夜激情福利司机影院| 国产av麻豆久久久久久久| 99久久久亚洲精品蜜臀av| 午夜福利免费观看在线| 99精品久久久久人妻精品| 可以在线观看毛片的网站| www.精华液| 精品福利观看| 久久99热这里只有精品18| 最近最新免费中文字幕在线| 黄色丝袜av网址大全| 1000部很黄的大片| 中文字幕熟女人妻在线| 精品久久久久久成人av| 亚洲精品在线美女| 视频区欧美日本亚洲| 精品日产1卡2卡| 色老头精品视频在线观看| 国产精品久久久人人做人人爽| 午夜精品在线福利| 亚洲av成人一区二区三| 久久久久国内视频| 日韩有码中文字幕| 亚洲国产日韩欧美精品在线观看 | 精品电影一区二区在线| 欧美日韩乱码在线| 国产成人欧美在线观看| h日本视频在线播放| 丰满的人妻完整版| 成年人黄色毛片网站| 成人一区二区视频在线观看| 国产高清视频在线播放一区| 亚洲中文av在线| 夜夜夜夜夜久久久久| 国产精品乱码一区二三区的特点| 五月伊人婷婷丁香| 在线观看午夜福利视频| 一级毛片女人18水好多| 国产精品久久久久久亚洲av鲁大| 久久香蕉国产精品| 97超级碰碰碰精品色视频在线观看| 在线观看舔阴道视频| avwww免费| 午夜精品在线福利| 看黄色毛片网站| 亚洲自拍偷在线| 久久久久久久精品吃奶| 精品久久久久久成人av| 亚洲中文字幕日韩| 久久久精品欧美日韩精品| 熟女少妇亚洲综合色aaa.| av天堂中文字幕网| 美女高潮喷水抽搐中文字幕| tocl精华| 亚洲一区二区三区不卡视频| 婷婷六月久久综合丁香| 给我免费播放毛片高清在线观看| 亚洲欧美日韩高清在线视频| 国产三级在线视频| 国产精品久久视频播放| 国产 一区 欧美 日韩| 变态另类成人亚洲欧美熟女| a级毛片a级免费在线| tocl精华| 一区二区三区国产精品乱码| 午夜a级毛片| 麻豆国产av国片精品| 9191精品国产免费久久| 久久久久亚洲av毛片大全| 岛国在线免费视频观看| 久久久精品大字幕| av天堂在线播放| 免费人成视频x8x8入口观看| 一级毛片女人18水好多| 亚洲乱码一区二区免费版| 亚洲精华国产精华精| 中文字幕久久专区| 欧美日本亚洲视频在线播放| 美女cb高潮喷水在线观看 | 变态另类丝袜制服| 欧美一级毛片孕妇| 免费在线观看影片大全网站| 黄频高清免费视频| 欧美大码av| 精品人妻1区二区| 免费搜索国产男女视频| 可以在线观看毛片的网站| 99re在线观看精品视频| 一本久久中文字幕| 亚洲国产日韩欧美精品在线观看 | 一进一出好大好爽视频| 久9热在线精品视频| 无限看片的www在线观看| 黄片小视频在线播放| 国内精品久久久久久久电影| 国产精品一区二区精品视频观看| 哪里可以看免费的av片| 99久国产av精品| av在线蜜桃| 成人午夜高清在线视频| 精品一区二区三区av网在线观看| 最好的美女福利视频网| 国模一区二区三区四区视频 | 国产欧美日韩精品亚洲av| 欧美丝袜亚洲另类 | 2021天堂中文幕一二区在线观| 欧美成人一区二区免费高清观看 | 国产黄a三级三级三级人| 色av中文字幕| 一进一出抽搐gif免费好疼| 欧美激情在线99| 99国产精品99久久久久| 男女做爰动态图高潮gif福利片| 狂野欧美白嫩少妇大欣赏| 国产精品99久久99久久久不卡| 噜噜噜噜噜久久久久久91| 伦理电影免费视频| 国产精品久久久久久亚洲av鲁大| 麻豆av在线久日| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 99热这里只有精品一区 | 最近视频中文字幕2019在线8| 国产久久久一区二区三区| 婷婷精品国产亚洲av在线| 欧美xxxx黑人xx丫x性爽| 一区二区三区国产精品乱码| 18禁国产床啪视频网站| 色综合欧美亚洲国产小说| 欧美又色又爽又黄视频| 中国美女看黄片| 一个人看视频在线观看www免费 | 中文字幕熟女人妻在线| 午夜精品一区二区三区免费看| 久久精品夜夜夜夜夜久久蜜豆| 精品无人区乱码1区二区| 悠悠久久av| 又大又爽又粗| 一本综合久久免费| 十八禁网站免费在线| 亚洲男人的天堂狠狠| 国产av在哪里看| 亚洲欧美日韩高清专用| 手机成人av网站| 久久久精品欧美日韩精品| 国产 一区 欧美 日韩| 国产三级在线视频| 最近最新免费中文字幕在线| 日韩欧美在线乱码| 久久午夜综合久久蜜桃| 国产野战对白在线观看| 男人的好看免费观看在线视频| 欧美乱码精品一区二区三区| 久久久久久久久久黄片| 99久国产av精品| 亚洲美女黄片视频| 国产真实乱freesex| 一级黄色大片毛片| 一个人观看的视频www高清免费观看 | 在线观看午夜福利视频| 国产日本99.免费观看| 一级黄色大片毛片| 九九热线精品视视频播放| 日本三级黄在线观看| 亚洲 国产 在线| 这个男人来自地球电影免费观看| 亚洲国产日韩欧美精品在线观看 | 国产精品久久久久久精品电影| 欧美3d第一页| 国产亚洲精品久久久久久毛片| 观看美女的网站| 99久久99久久久精品蜜桃| 十八禁网站免费在线| 国产高清videossex| 18禁黄网站禁片免费观看直播| 国产成人av教育| 日韩中文字幕欧美一区二区| 国产黄色小视频在线观看| 高潮久久久久久久久久久不卡| 美女免费视频网站| 很黄的视频免费| 999精品在线视频| 天天躁日日操中文字幕| 99久久99久久久精品蜜桃| 久久精品影院6| 一本精品99久久精品77| 人人妻人人澡欧美一区二区| 亚洲精品中文字幕一二三四区| 国产1区2区3区精品| 欧美黑人欧美精品刺激| 日韩欧美在线二视频| 日韩欧美精品v在线| 一进一出抽搐动态| 中文亚洲av片在线观看爽| 欧美zozozo另类| 亚洲午夜理论影院| 19禁男女啪啪无遮挡网站| 亚洲人成网站高清观看| 啦啦啦免费观看视频1| 午夜激情福利司机影院| 午夜福利在线观看免费完整高清在 | 极品教师在线免费播放| 人人妻人人看人人澡| 亚洲精品粉嫩美女一区| 国产又黄又爽又无遮挡在线| 国产一区在线观看成人免费| 色视频www国产| 国内精品美女久久久久久| 欧美黑人巨大hd| 香蕉久久夜色| a级毛片a级免费在线| 免费一级毛片在线播放高清视频| 在线免费观看的www视频| 丰满人妻熟妇乱又伦精品不卡| 99re在线观看精品视频| e午夜精品久久久久久久| 好男人电影高清在线观看| 欧美色视频一区免费| 亚洲 欧美 日韩 在线 免费| 最近最新免费中文字幕在线| 欧美不卡视频在线免费观看| 国产精品女同一区二区软件 | 国产精品一区二区三区四区免费观看 | 亚洲欧洲精品一区二区精品久久久| 97人妻精品一区二区三区麻豆| 又爽又黄无遮挡网站| 国产精品影院久久| 日日干狠狠操夜夜爽| 国产伦精品一区二区三区四那| 91麻豆精品激情在线观看国产| 精品日产1卡2卡| 搡老岳熟女国产| 欧美成人免费av一区二区三区| 男女午夜视频在线观看| xxxwww97欧美| 男人舔奶头视频| 精品无人区乱码1区二区| 大型黄色视频在线免费观看| 18禁黄网站禁片午夜丰满| 久久久久久人人人人人| 99在线视频只有这里精品首页| 精品国产亚洲在线| 国产精品久久视频播放| 国产免费男女视频| 久久中文字幕人妻熟女| 中文字幕高清在线视频| 淫秽高清视频在线观看| 精品久久久久久久人妻蜜臀av| 一二三四社区在线视频社区8| 99精品在免费线老司机午夜| 床上黄色一级片| 免费一级毛片在线播放高清视频| 亚洲精品粉嫩美女一区| 小蜜桃在线观看免费完整版高清| 久久久久久久久久黄片| 午夜成年电影在线免费观看| 亚洲成人中文字幕在线播放| 日本黄色片子视频| 亚洲成人精品中文字幕电影| 可以在线观看的亚洲视频| 91字幕亚洲| 国产视频内射| 亚洲国产精品999在线| 夜夜夜夜夜久久久久| 国产三级在线视频| 欧美色欧美亚洲另类二区| 伊人久久大香线蕉亚洲五| 熟女人妻精品中文字幕| 免费在线观看影片大全网站| 好男人在线观看高清免费视频| 12—13女人毛片做爰片一| 午夜精品在线福利| 我的老师免费观看完整版| 香蕉国产在线看| 身体一侧抽搐| 精品熟女少妇八av免费久了| 国产精品影院久久| 亚洲精品456在线播放app | 日本黄大片高清| 麻豆成人午夜福利视频| 伊人久久大香线蕉亚洲五| 在线免费观看不下载黄p国产 | 亚洲,欧美精品.| 国内精品久久久久精免费| 久久精品影院6| 国产1区2区3区精品| 国产精品日韩av在线免费观看| 成人精品一区二区免费| 欧美性猛交╳xxx乱大交人| 亚洲国产日韩欧美精品在线观看 | 成在线人永久免费视频| 香蕉久久夜色| 精品国产乱码久久久久久男人| 熟妇人妻久久中文字幕3abv| 亚洲精品中文字幕一二三四区| 白带黄色成豆腐渣| 啪啪无遮挡十八禁网站| 又紧又爽又黄一区二区| 美女cb高潮喷水在线观看 | 1024手机看黄色片| 男人舔奶头视频| 国产97色在线日韩免费| 日本五十路高清| 精品久久蜜臀av无| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 亚洲av美国av| 99国产精品99久久久久| 国产精品av久久久久免费| 99国产精品一区二区三区| 色播亚洲综合网| 国产视频一区二区在线看| 真人一进一出gif抽搐免费| 999久久久精品免费观看国产| 麻豆成人午夜福利视频| 亚洲av成人不卡在线观看播放网| 午夜激情欧美在线| 国产av在哪里看| 蜜桃久久精品国产亚洲av| 欧美黄色片欧美黄色片| 中文字幕av在线有码专区| 美女黄网站色视频| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久久黄片| 国内精品久久久久精免费| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 99国产综合亚洲精品| www日本在线高清视频| 国模一区二区三区四区视频 | 丝袜人妻中文字幕| 色哟哟哟哟哟哟| 在线免费观看不下载黄p国产 | 国产一区二区激情短视频| 熟妇人妻久久中文字幕3abv| 国产精品 国内视频|