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

    不同溫度反應(yīng)與精餾集成生產(chǎn)醋酸叔丁酯的過程模擬

    2015-02-14 09:34:34黃玉鑫湯吉海陳獻(xiàn)宋建國崔咪芬費(fèi)兆陽喬旭
    化工學(xué)報(bào) 2015年10期
    關(guān)鍵詞:異丁烯板數(shù)丁酯

    黃玉鑫,湯吉海,陳獻(xiàn),宋建國,崔咪芬,費(fèi)兆陽,喬旭

    (1南京工業(yè)大學(xué)化學(xué)化工學(xué)院,江蘇省工業(yè)節(jié)水減排重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210009;2南京工業(yè)大學(xué)材料化學(xué)工程國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210009)

    不同溫度反應(yīng)與精餾集成生產(chǎn)醋酸叔丁酯的過程模擬

    黃玉鑫1,湯吉海1,陳獻(xiàn)1,宋建國1,崔咪芬1,費(fèi)兆陽2,喬旭2

    (1南京工業(yè)大學(xué)化學(xué)化工學(xué)院,江蘇省工業(yè)節(jié)水減排重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210009;2南京工業(yè)大學(xué)材料化學(xué)工程國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210009)

    針對醋酸與異丁烯加成酯化可逆反應(yīng)溫度低、精餾分離溫度高的特點(diǎn),采用帶側(cè)反應(yīng)器的反應(yīng)精餾集成過程(SRC)建立了低溫反應(yīng)與高溫精餾集成的醋酸叔丁酯生產(chǎn)新工藝。固定塔釜上升汽化量100 kmol·h?1,規(guī)定新鮮醋酸進(jìn)料的轉(zhuǎn)化率達(dá)到99.9%、醋酸叔丁酯選擇性達(dá)到97.0%,采用過程模擬考察了進(jìn)入側(cè)反應(yīng)器的精餾塔采出量、精餾段塔板數(shù)、側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)和側(cè)反應(yīng)器臺數(shù)等參數(shù)對合成醋酸叔丁酯的SRC過程的影響。模擬結(jié)果表明,醋酸與異丁烯加成酯化生產(chǎn)醋酸叔丁酯的SRC過程中只有反應(yīng)能力與分離能力達(dá)到最佳匹配才能使單位產(chǎn)品的生產(chǎn)成本最小。研究結(jié)果為醋酸叔丁酯生產(chǎn)新工藝的放大設(shè)計(jì)與優(yōu)化奠定了基礎(chǔ)。

    反應(yīng)精餾;反應(yīng)器;酯化;醋酸叔丁酯;集成

    Key words:reactive distillation;reactors;esterification;tert-butyl acetate;integration

    引 言

    醋酸叔丁酯是一種重要的含氧有機(jī)溶劑,被美國國家環(huán)保局列為非揮發(fā)性有機(jī)溶劑和非空氣有害污染物,是一種可用于涂料、油墨、黏合劑和工業(yè)清洗的綠色環(huán)保溶劑[1-3]。由于叔丁醇的空間位阻大,傳統(tǒng)醋酸叔丁酯生產(chǎn)采用醋酐和叔丁醇為原料進(jìn)行乙?;磻?yīng)生產(chǎn)[4],該方法副產(chǎn)醋酸導(dǎo)致后續(xù)分離能耗較大。而醋酸與異丁烯直接加成酯化合成醋酸叔丁酯的反應(yīng)屬于100%的原子經(jīng)濟(jì)反應(yīng)[5],是一種綠色原子經(jīng)濟(jì)反應(yīng)工藝[6]。

    然而,醋酸與異丁烯的加成酯化反應(yīng)是一個可逆反應(yīng),反應(yīng)不僅受熱力學(xué)平衡限制,而且在酯化反應(yīng)的同時伴隨異丁烯二聚副反應(yīng)。為了抑制異丁烯二聚副反應(yīng),不僅需要加入叔丁醇等極性溶劑,而且要求反應(yīng)在較低溫度下進(jìn)行(約30℃)[7]。雖然反應(yīng)精餾集成技術(shù)[8]在許多羧酸與醇酯化合成有機(jī)酯的反應(yīng)中得到應(yīng)用[9-11],可提高轉(zhuǎn)化率和選擇性,降低分離能耗,但是,醋酸叔丁酯和醋酸的正常沸點(diǎn)分別為96℃和118℃,反應(yīng)物和產(chǎn)物可以采用常壓精餾分離,精餾塔內(nèi)的溫度遠(yuǎn)遠(yuǎn)高于醋酸與異丁烯加成酯化的反應(yīng)溫度,在一個精餾塔內(nèi)耦合的傳統(tǒng)反應(yīng)精餾過程不適用于羧酸與異丁烯的酯化反應(yīng)[12]。帶側(cè)反應(yīng)器的反應(yīng)精餾集成過程(SRC)中,側(cè)反應(yīng)器與精餾塔相互有質(zhì)量和能量交換,可以達(dá)到與傳統(tǒng)反應(yīng)精餾集成過程相同的效果[13-19]。而且,SRC過程的側(cè)反應(yīng)器和精餾塔又處于相對獨(dú)立的空間,不僅具有放大容易的優(yōu)勢,更可以使反應(yīng)器和精餾塔分別在不同的溫度或壓力工況下操作,大大拓展了反應(yīng)精餾集成過程的應(yīng)用范圍[20]。

    已有的SRC過程研究主要針對連串反應(yīng)或者可逆反應(yīng),而醋酸與異丁烯的反應(yīng)過程既存在加成酯化可逆反應(yīng),也存在異丁烯二聚的副反應(yīng),屬于可逆-平行復(fù)雜反應(yīng)體系,而且酯化反應(yīng)和精餾溫度差異大,對該過程的反應(yīng)精餾強(qiáng)化工藝研究還處于空白。為此,本文提出了較低溫度反應(yīng)與較高溫度精餾集成生產(chǎn)醋酸叔丁酯的新型反應(yīng)精餾集成工藝,采用Aspen Plus軟件對該過程進(jìn)行模擬研究,為不同溫度反應(yīng)與精餾集成生產(chǎn)醋酸叔丁酯的過程設(shè)計(jì)與優(yōu)化提供參考。

    1 工藝流程

    較低溫度反應(yīng)與較高溫度精餾集成生產(chǎn)醋酸叔丁酯的工藝流程如圖1所示。精餾塔內(nèi)進(jìn)行組分的分離過程,醋酸叔丁酯沸點(diǎn)低于醋酸,在塔頂冷凝器冷凝后與少量副產(chǎn)品二異丁烯作為粗產(chǎn)品采出,醋酸在精餾塔底部富集,進(jìn)入各臺側(cè)反應(yīng)器反應(yīng)。側(cè)反應(yīng)器和精餾塔均為常壓操作,側(cè)反應(yīng)器的操作溫度維持在30℃,精餾塔的操作溫度范圍為91~117℃。新鮮醋酸原料從側(cè)反應(yīng)器1中加入,異丁烯按一定比例加入各臺側(cè)反應(yīng)器中,側(cè)反應(yīng)器出口物料返回上方塔板。因?yàn)閭?cè)反應(yīng)器溫度遠(yuǎn)低于精餾塔,為了提高熱量利用率,每臺側(cè)反應(yīng)器配置一臺換熱器,對進(jìn)反應(yīng)器熱物料與出反應(yīng)器冷物料進(jìn)行換熱。其中側(cè)反應(yīng)器、換熱器與精餾塔塔板的編號均從下到上計(jì)數(shù)。

    圖1 不同溫度反應(yīng)與精餾集成裝置Fig. 1 Structure of distillation column integrated with side reactors under different temperatures

    2 模擬方法

    本文采用化工流程模擬軟件Aspen Plus進(jìn)行過程模擬,因?yàn)镾RC過程的側(cè)反應(yīng)器和精餾塔獨(dú)立,所以在Aspen Plus流程模擬中,精餾塔模擬采用嚴(yán)格精餾塔模型RadFrac,側(cè)反應(yīng)器采用化學(xué)計(jì)量反應(yīng)器模型RStoic,側(cè)反應(yīng)器和精餾塔之間的流股換熱模型采用HeatX。當(dāng)采用多臺側(cè)反應(yīng)器與精餾塔集成時,本文均對異丁烯在側(cè)反應(yīng)器中的分配進(jìn)行優(yōu)化。為了提高優(yōu)化計(jì)算效率,采用了基于獨(dú)立反應(yīng)量[21]的優(yōu)化設(shè)計(jì)方法。這種方法在給定了每臺側(cè)反應(yīng)器的異丁烯反應(yīng)量的情況下,需要根據(jù)側(cè)反應(yīng)器的進(jìn)料組成、反應(yīng)器大小和反應(yīng)動力學(xué)(表1)計(jì)算側(cè)反應(yīng)器中醋酸叔丁酯的選擇性。

    醋酸與異丁烯的反應(yīng)既存在加成酯化可逆反應(yīng),也存在異丁烯二聚副反應(yīng),其反應(yīng)方程式和反應(yīng)動力學(xué)模型見表1。在穩(wěn)定狀態(tài)下,醋酸叔丁酯和二異丁烯的選擇性關(guān)系式可用式(1)和式(2)表示[21]。

    將表1的動力學(xué)參數(shù)代入式(1),得到醋酸叔丁酯選擇性與側(cè)反應(yīng)器出口各物質(zhì)液相濃度的關(guān)系式(3)

    式中,cAc、cIB、cTBAc分別代表側(cè)反應(yīng)器出口物料中醋酸、異丁烯、醋酸叔丁酯的液相濃度,kmol·m?3。根據(jù)式(2)和式(3)編寫Fortran程序,與Aspen Plus中的RStoic反應(yīng)器模塊相連接,實(shí)現(xiàn)在給定的側(cè)反應(yīng)器進(jìn)口組成條件下,直接迭代求出側(cè)反應(yīng)器出口組成。

    由于醋酸和醋酸叔丁酯屬于強(qiáng)液相非理想體系,且存在氣相中醋酸的締合現(xiàn)象[22],故在Aspen Plus中選擇UNIFAC-HOC模型,液相活度采用UNIFAC預(yù)測,氣相逸度采用狀態(tài)方程Hayden-O?Connell預(yù)測。

    本文選用的優(yōu)化目標(biāo)為單位產(chǎn)品成本最低,用年度總費(fèi)用與年產(chǎn)量的比值表示,年度總費(fèi)用為設(shè)備固定投資費(fèi)用與能源費(fèi)用之和。其中,設(shè)備固定投資費(fèi)用包括精餾塔、側(cè)反應(yīng)器和換熱器的費(fèi)用,采用Douglas[23]對化工過程的估算方法;能源費(fèi)用包括再沸器的蒸汽費(fèi)用和冷凝器與維持側(cè)反應(yīng)器恒溫的循環(huán)冷卻水費(fèi)用,蒸汽和循環(huán)水單價分別按250 元/噸與0.5 元/噸計(jì)算。

    表1 醋酸與異丁烯加成酯化合成醋酸叔丁酯動力學(xué)特性Table 1 Kinetic characteristics oftert-butyl acetate synthesis with acetic acid and isobutylene

    3 結(jié)果與討論

    規(guī)定新鮮醋酸進(jìn)料的轉(zhuǎn)化率達(dá)到99.9%、醋酸叔丁酯選擇性達(dá)到97.0%,固定塔釜上升汽化量(V)為100 kmol·h?1,優(yōu)化異丁烯進(jìn)料量(F)和異丁烯分配系數(shù)(fk)使單位產(chǎn)品成本最小。在此最優(yōu)條件下,考察進(jìn)入側(cè)反應(yīng)器的精餾塔采出量(B)、精餾段塔板數(shù)(NS)、側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)(NSR)和反應(yīng)器臺數(shù)(NR)對SRC過程的影響。

    3.1 進(jìn)入側(cè)反應(yīng)器的精餾塔采出量的影響

    因?yàn)轷セ磻?yīng)產(chǎn)物醋酸叔丁酯從反應(yīng)精餾塔頂采出,原料醋酸從塔釜采出進(jìn)入第1臺側(cè)反應(yīng)器,反應(yīng)后返回至第8塊塔板,從第9塊塔板側(cè)線采出進(jìn)入第2臺側(cè)反應(yīng)器,反應(yīng)后返回至第16塊塔板。全文均假設(shè)進(jìn)入兩臺側(cè)反應(yīng)器的精餾塔采出量相同。在V=100 kmol·h?1,NS=22,NSR=7,NR=2的條件下,考察進(jìn)入側(cè)反應(yīng)器的精餾塔采出量對SRC過程的影響,結(jié)果如表2所示。

    從表2可以看出,在精餾塔上升汽化量固定(即精餾塔塔徑固定)的情況下,隨著進(jìn)入側(cè)反應(yīng)器的精餾塔采出量從10 kmol·h?1增加到40 kmol·h?1,醋酸叔丁酯年產(chǎn)量從5552 t·a?1增加到13106 t·a?1,增加了1.4倍,并且單位產(chǎn)品成本從626元/噸顯著減小到294 元/噸,下降了53.0%。當(dāng)進(jìn)入側(cè)反應(yīng)器的精餾塔采出量超過40 kmol·h?1時,年產(chǎn)量增加趨勢變緩,單位產(chǎn)品成本由下降趨勢轉(zhuǎn)為上升趨勢。這是因?yàn)樵黾舆M(jìn)入側(cè)反應(yīng)器的精餾塔采出量實(shí)際是增加了進(jìn)入側(cè)反應(yīng)器的進(jìn)料流率,使得側(cè)反應(yīng)器的反應(yīng)能力增加。而在精餾塔的直徑和上升汽化量不變的情況下,精餾塔的設(shè)備投資費(fèi)用和能源費(fèi)用不變,側(cè)反應(yīng)器的進(jìn)料流率增加導(dǎo)致反應(yīng)器投資費(fèi)用和換熱費(fèi)用增加。當(dāng)反應(yīng)能力、投資費(fèi)用和能源費(fèi)用同時增加時,單位產(chǎn)品的成本下降、增加或不變的可能性都存在。當(dāng)進(jìn)入側(cè)反應(yīng)器的精餾塔采出量在40 kmol·h?1以下時,增加進(jìn)入側(cè)反應(yīng)器的精餾塔采出量,醋酸叔丁酯產(chǎn)量的增加幅度大于生產(chǎn)成本增加幅度,所以單位產(chǎn)品的總成本下降。但是當(dāng)進(jìn)入側(cè)反應(yīng)器的精餾塔采出量大于40 kmol·h?1時,增加進(jìn)入側(cè)反應(yīng)器的精餾塔采出量后,醋酸叔丁酯產(chǎn)量基本不增加,而反應(yīng)器投資費(fèi)用和反應(yīng)物料流股換熱能源費(fèi)用增加,所以單位產(chǎn)品的成本開始增加。以上結(jié)果說明在SRC集成過程中,反應(yīng)能力和分離能力存在最佳匹配,只有當(dāng)反應(yīng)能力和分離能力達(dá)到最佳匹配,才能使過程達(dá)到最優(yōu),反應(yīng)能力或者分離能力某一方面過剩均使過程的成本增加。

    3.2 精餾段塔板數(shù)的影響

    在V=100 kmol·h?1,B=40 kmol·h?1,NSR=7,NR=2的條件下,考察精餾段塔板數(shù)對SRC過程的影響,結(jié)果如表3和圖2所示。

    從表3可以看出,精餾段塔板數(shù)從14增加到22,醋酸叔丁酯年產(chǎn)量從11107 t·a?1增加到13106 t·a?1,增加了17.2%,當(dāng)精餾段塔板數(shù)超過22塊之后,繼續(xù)增加精餾段塔板數(shù),醋酸叔丁酯年產(chǎn)量增加很小。精餾段塔板數(shù)從14增加到22,單位產(chǎn)品能源費(fèi)用從256 元/噸下降至215 元/噸,減小16.0%,當(dāng)精餾段塔板數(shù)超過22塊之后,繼續(xù)增加精餾段塔板數(shù),能源費(fèi)用基本不變。精餾段塔板數(shù)從14增加到30,雖然單位產(chǎn)品反應(yīng)器投資費(fèi)用從35元/噸下降至32元/噸,換熱器投資費(fèi)用從25元/噸下降至21 元/噸,兩者分別減小8.6%和16.0%,但是隨著精餾段塔板數(shù)增加到30,單位產(chǎn)品精餾塔投資費(fèi)用逐漸增加,最終導(dǎo)致單位產(chǎn)品設(shè)備費(fèi)用先減小后增加。

    圖2 精餾段塔板數(shù)對塔釜醋酸液相摩爾組成和再沸器能耗的影響Fig. 2 Effects of NS on liquid mole fraction of acetic acid in bottom and reboiler duty (V=100 kmol·h?1,B=40 kmol·h?1, NSR=7, NR=2)

    表2 進(jìn)入側(cè)反應(yīng)器的精餾塔采出量對SRC過程的影響Table 2 Effect ofBon SRC process

    表3 精餾段塔板數(shù)對SRC過程的影響Table 3 Effects of NS on SRC process

    精餾段塔板數(shù)的增加提高了系統(tǒng)的分離能力,使醋酸叔丁酯從精餾段移出能力增強(qiáng),從精餾塔采出液相物料流股中醋酸的摩爾分?jǐn)?shù)增大(如圖2所示),從而促進(jìn)反應(yīng)平衡向有利于生成醋酸叔丁酯的方向移動,提高了SRC過程的反應(yīng)能力和醋酸叔丁酯產(chǎn)量。另外,因?yàn)榇姿岬哪柶剩?3.93 MJ·kmol?1)比醋酸叔丁酯(31.61 MJ·kmol?1)小24.3%,所以在塔釜上升汽化量一定的情況下,塔釜中醋酸的摩爾分?jǐn)?shù)增加還可以有效降低塔釜加熱負(fù)荷,減少能源費(fèi)用。從圖2還可以看出在塔釜上升汽化量固定的情況下,當(dāng)精餾段塔板數(shù)小于26時,隨著精餾段塔板數(shù)增加,塔釜醋酸的液相組成增大,且再沸器熱負(fù)荷下降。當(dāng)精餾段塔板數(shù)大于26時,增加精餾段塔板數(shù),塔釜醋酸的液相組成變化已不明顯,說明此時精餾段塔板數(shù)的增加對分離效率的增加不明顯,故反應(yīng)能力的增加很小,即醋酸叔丁酯年產(chǎn)量增加趨勢變緩。同時,當(dāng)精餾段塔板數(shù)大于26時,塔釜組成隨精餾段塔板數(shù)變化不明顯,再沸器蒸發(fā)量基本不變,也就是能源費(fèi)用基本不變。同時,單位產(chǎn)品的反應(yīng)器費(fèi)用和換熱器費(fèi)用基本不變。在此情況下,繼續(xù)增加精餾段塔板數(shù),單位產(chǎn)品精餾塔費(fèi)用仍然增加,因此,單位產(chǎn)品成本呈現(xiàn)增加趨勢。

    表4 NSR對SRC過程的影響Table 4 Effect of NSR on SRC process

    3.3 側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)的影響

    在V=100 kmol·h?1,B=40 kmol·h?1,NS=22,NR=2的條件下,考察側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)對SRC過程的影響,模擬結(jié)果如表4和圖3所示。

    從表4可以看出,當(dāng)NSR從1增加到7時,醋酸叔丁酯的年產(chǎn)量從10961 t·a?1增加到13106 t·a?1,增幅為19.6%,單位產(chǎn)品能源費(fèi)用顯著降低,從356 元/噸下降至294 元/噸,減小了21.5%,而單位產(chǎn)品設(shè)備費(fèi)用變化并不明顯,單位產(chǎn)品成本逐漸降低。當(dāng)NSR超過7之后,醋酸叔丁酯年產(chǎn)量增加與單位產(chǎn)品能源費(fèi)用下降緩慢,而單位產(chǎn)品設(shè)備投資中的精餾塔費(fèi)用仍然增加,導(dǎo)致單位產(chǎn)品成本降低趨勢放緩。

    圖3顯示了不同NSR下,反應(yīng)精餾段醋酸叔丁酯的液相摩爾分?jǐn)?shù)的分布,圖中的向右箭頭表示進(jìn)入側(cè)反應(yīng)器的反應(yīng)精餾塔采出流股塔板位置,向左箭頭表示側(cè)反應(yīng)器進(jìn)入反應(yīng)精餾塔的塔板位置。

    圖3 側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)對連接側(cè)反應(yīng)器精餾段內(nèi)醋酸叔丁酯液相摩爾分?jǐn)?shù)的影響Fig. 3 Effects of NSR on liquid mole fraction oftert-butyl acetate in distillation zone connected with side reactors(V=100 kmol·h?1,B=40 kmol·h?1, NS=22, NR=2)

    從圖3可以看出,隨著NSR從1增加到9時,進(jìn)入側(cè)反應(yīng)器1的反應(yīng)精餾塔采出流股中醋酸叔丁酯含量顯著減小,而側(cè)反應(yīng)器2進(jìn)入反應(yīng)精餾塔的塔板上醋酸叔丁酯含量顯著增加;側(cè)反應(yīng)器1進(jìn)入反應(yīng)精餾塔的塔板上和進(jìn)入側(cè)反應(yīng)器2的反應(yīng)精餾塔采出流股中的醋酸叔丁酯含量增加很小。這是由于增加NSR,增加了反應(yīng)精餾段的分離能力,使得反應(yīng)精餾段上部塔板的醋酸叔丁酯含量增加。增加側(cè)反應(yīng)器2的NSR,不僅增加了側(cè)反應(yīng)器2生成醋酸叔丁酯產(chǎn)物的分離能力,而且增加了從側(cè)反應(yīng)器1進(jìn)入精餾塔的反應(yīng)產(chǎn)物分離的精餾段塔板數(shù),進(jìn)一步提升了側(cè)反應(yīng)器1生成產(chǎn)物的分離能力。酸烯酯化反應(yīng)是受反應(yīng)熱力學(xué)平衡限制的可逆反應(yīng),從精餾塔采出進(jìn)入側(cè)反應(yīng)器的反應(yīng)物濃度增加和產(chǎn)物濃度減小有利于酯化反應(yīng)。因此,從精餾塔釜采出進(jìn)入側(cè)反應(yīng)器1的進(jìn)料流股中醋酸叔丁酯摩爾分?jǐn)?shù)隨著NSR增加而降低,也就是進(jìn)入側(cè)反應(yīng)器1的產(chǎn)物濃度低,反應(yīng)物濃度高,從而增加了側(cè)反應(yīng)器1的反應(yīng)推動力,提高側(cè)反應(yīng)器1的反應(yīng)能力。

    3.4 側(cè)反應(yīng)器臺數(shù)的影響

    增加側(cè)反應(yīng)器臺數(shù),不僅增加了反應(yīng)能力,同時也增加了反應(yīng)精餾段塔板數(shù),提升分離能力,有利于反應(yīng)能力的提升。在V=100 kmol·h?1,B=40 kmol·h?1, NS=22, NSR=7的條件下,考察了側(cè)反應(yīng)器臺數(shù)的影響,模擬結(jié)果如表5所示。

    表5 側(cè)反應(yīng)器臺數(shù)對SRC過程的影響Table 5 Effect of NR on SRC process

    從表5可以看出,2臺側(cè)反應(yīng)器與1臺相比,醋酸叔丁酯年產(chǎn)量從11043 t·a?1增加到13106 t·a?1,增加了18.7%,單位產(chǎn)品能源費(fèi)用從253元/噸下降至215 元/噸,降低15.0%,單位產(chǎn)品設(shè)備費(fèi)用雖然略有增加,但是單位產(chǎn)品成本還是從325元/噸下降至294 元/噸,降低9.5%。當(dāng)側(cè)反應(yīng)器臺數(shù)增加到3臺時,醋酸叔丁酯年產(chǎn)量與單位產(chǎn)品能源費(fèi)用的變化并不明顯,但是單位產(chǎn)品設(shè)備費(fèi)用增加。與2臺側(cè)反應(yīng)器相比,3臺側(cè)反應(yīng)器的單位產(chǎn)品成本增加9.2%。從表5還可以看出,當(dāng)側(cè)反應(yīng)器臺數(shù)為2臺或3臺時,異丁烯對各臺側(cè)反應(yīng)器的分配系數(shù)從精餾塔釜向上依次逐漸減小,第1臺側(cè)反應(yīng)器的異丁烯分配系數(shù)達(dá)到0.651以上,表明大部分的反應(yīng)能力集中在下部第1臺側(cè)反應(yīng)器中進(jìn)行。3臺側(cè)反應(yīng)器與2臺側(cè)反應(yīng)器相比,f1的變化并不顯著,這表明增加的側(cè)反應(yīng)器主要分擔(dān)了上部第2臺側(cè)反應(yīng)器的反應(yīng)能力。根據(jù)圖3的分析可以知道,第2臺側(cè)反應(yīng)器的進(jìn)料流股中醋酸叔丁酯的摩爾分?jǐn)?shù)達(dá)到0.2,產(chǎn)物濃度較高,對反應(yīng)的推動力較小,因而反應(yīng)能力的增加有限。而f3小于0.1也說明了第3臺側(cè)反應(yīng)器提升產(chǎn)能較小,而設(shè)備費(fèi)用與能源費(fèi)用的增加成為主要影響因素。

    3.5 最優(yōu)結(jié)構(gòu)下SRC過程的精餾塔內(nèi)參數(shù)分布

    綜合3.1~3.4節(jié)的結(jié)果與討論,固定精餾塔釜汽化量100 kmol·h?1時,得到醋酸與異丁烯加成酯化SRC過程的最優(yōu)結(jié)構(gòu)參數(shù)為:V=100 kmol·h?1,B=40 kmol·h?1,NS=22,NSR=7,NR=2,f1=0.675,f2=0.325。在此最優(yōu)參數(shù)下SRC過程的精餾塔塔板溫度分布與塔內(nèi)液相摩爾組成分布如圖4所示。

    圖4 最優(yōu)結(jié)構(gòu)下的塔板溫度分布與塔內(nèi)液相摩爾組成分布Fig. 4 Temperature profiles and liquid composition profiles with optimum design parameters(V=100 kmol·h?1,B=40 kmol·h?1, NS=22, NSR=7, NR=2,f1=0.675,f2=0.325)

    從圖4可以看出,在精餾的分離作用下,產(chǎn)品醋酸叔丁酯不斷從反應(yīng)精餾段移出,在精餾塔頂部富集。反應(yīng)精餾段以及側(cè)反應(yīng)器的進(jìn)料流股保持醋酸的濃度高、醋酸叔丁酯的濃度低,有利于異丁烯與醋酸發(fā)生酯化反應(yīng),降低異丁烯二聚副反應(yīng)的選擇性。反應(yīng)精餾段的醋酸和醋酸叔丁酯的組成分布、溫度分布呈現(xiàn)單調(diào)變化規(guī)律,說明側(cè)反應(yīng)器沒有帶來逆向精餾副作用[24]。說明本文固定對應(yīng)反應(yīng)器2的精餾塔采出與反應(yīng)器1的返塔位置相鄰、每臺側(cè)反應(yīng)器的進(jìn)入側(cè)反應(yīng)器的精餾塔采出量(B)以及側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)(NSR)假設(shè)為相等是合理的。

    4 結(jié) 論

    (1)針對醋酸與異丁烯加成酯化反應(yīng)可逆平行的反應(yīng)特點(diǎn),建立了較低溫度反應(yīng)與較高溫度精餾的帶側(cè)反應(yīng)器反應(yīng)精餾集成過程,模擬結(jié)果表明該過程是可行的。

    (2)通過模擬考察了進(jìn)入側(cè)反應(yīng)器的精餾塔采出量、精餾段塔板數(shù)、側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)和側(cè)反應(yīng)器臺數(shù)等參數(shù)對合成醋酸叔丁酯SRC過程的影響。模擬結(jié)果表明SRC過程的反應(yīng)能力與分離能力存在最佳匹配,否則將導(dǎo)致過程的經(jīng)濟(jì)性下降。

    (3)固定精餾塔釜汽化量100 kmol·h?1,規(guī)定新鮮醋酸進(jìn)料的轉(zhuǎn)化率達(dá)到99.9%、醋酸叔丁酯選擇性達(dá)到97.0%,得到了最佳的參數(shù):塔釜采出量為40 kmol·h?1、精餾段塔板數(shù)為22、側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)為7、側(cè)反應(yīng)器臺數(shù)為2臺、異丁烯在兩臺側(cè)反應(yīng)器中的分配系數(shù)分別為0.675和0.325時,SRC過程生產(chǎn)能力達(dá)到13106 t·a?1,此時單位產(chǎn)品成本達(dá)到最小為294 元/噸。

    符 號 說 明

    B——塔底采出量,kmol·h?1

    F——異丁烯進(jìn)料流率,kmol·h?1

    fk——異丁烯在各臺反應(yīng)器中的分配系數(shù)

    Ke——平衡常數(shù),L·mol?1

    k1——主反應(yīng)正反應(yīng)速率常數(shù),L·mol?1·h?1

    k2——副反應(yīng)速率常數(shù),L·mol?1·h?1

    NR ——側(cè)反應(yīng)器臺數(shù)

    NS ——精餾段塔板數(shù)

    NSR ——側(cè)反應(yīng)器進(jìn)出口間隔塔板數(shù)

    r1——主反應(yīng)醋酸叔丁酯生成速率,mol·L?1·h?1

    r2——副反應(yīng)二異丁烯生成速率,mol·L?1·h?1

    S1——主反應(yīng)醋酸叔丁酯的選擇性

    S2——副反應(yīng)二異丁烯的選擇性

    V——塔釜上升汽化量,kmol·h?1

    下角標(biāo)

    k——反應(yīng)器編號

    [1] Faber W, Kirkpatrick D, Coder P, Li A, Borghoff S, Banton M. Subchronic, reproductive, and maternal toxicity studies with tertiary butyl acetate (TBAC) [J].Regulatory Toxicology and Pharmacology, 2014, 68(3): 332-342.

    [2] Luo D, Corey R, Propper R, Collins J, Komorniczak A, Davis M, Berger N, Lum S. Comprehensive environmental impact assessment of exempt volatile organic compounds in California [J].Environmental Science & Policy, 2011, 14(6): 585-593.

    [3] Montón J B, Mu?oz R, Burguet M C, Torre J. Isobaric vapor-liquid equilibria for the binary systems isobutyl alcohol+isobutyl acetate andtert-butyl alcohol+tert -butyl acetate at 20 and 101.3kPa [J].Fluid Phase Equilibria, 2005, 227(1): 19-25.

    [4] Shyamprasad K, Shamshuddin S Z M, Vasantha V T. Effective acetylation of alcohols, phenols and amines over mesoporous aluminophosphate solid acids under solvent free conditions [J].Journal of Porous Materials, 2014, 21(6): 1079-1090.

    [5] Gu Y L, Shi F, Deng Y Q. Esterifcation of aliphatic acids with olefn promoted by Br?nsted acidic ionic liquids [J].Journal of Molecular Catalysis A:Chemical, 2004, 212(1/2): 71-75.

    [6] Trost B M. The atom economy - a search for synthetic efficiency [J].Science, 1991, 254(5037): 1471-1477.

    [7] Karas L J. Ester preparation[P]: US, 5994578. 1999-11-30.

    [8] Yu Guocong(余國琮), Yuan Xigang(袁希鋼), Li Genhao(李根浩). Review of papers on distillation research published in “CIESC Journal” in past 60 years [J].CIESC Journal(化工學(xué)報(bào)), 2013, 64(1): 11-27.

    [9] Malone M F, Doherty M F. Reactive distillation [J].Industrial & Engineering Chemistry Research, 2000, 39(11): 3953-3957.

    [10] Niesbach A, Fuhrmeister R, Keller T, Lutze P, Gorak A. Esterification of acrylic acid andn-butanol in a pilot-scale reactive distillation column-experimental investigation, model validation, and process analysis [J].Industrial & Engineering Chemistry Research, 2012, 51(50): 16444-16456.

    [11] Xiao Wu(肖武), Li Mingyue(李明月), Ruan Xuehua(阮雪華), He Gaohong(賀高紅), Du Jian(都健). Optimization of ethyl acetate process conditions for sodium bisulfate fluidized catalytic distillation using response surface methodology [J].CIESC Journal(化工學(xué)報(bào)), 2014, 65(11): 4465-4471.

    [12] Degarmo J L, Parulekar V N, Pinjala V. Consider reactive distillation [J].Chemical Engineering Progress, 1992, 88(3): 43-50.

    [13] Ding Lianghui(丁良輝), Tang Jihai(湯吉海), Cui Mifen(崔咪芬), Chen Xian(陳獻(xiàn)), Bo Cuimei(薄翠梅), Qiao Xu(喬旭). Analysis and comparison of RD and SRC involving the consecutive reaction of chlorination of toluene [J].CIESC Journal(化工學(xué)報(bào)), 2013, 64(9): 3277-3284.

    [14] Baur R, Krishna R. Distillation column with reactive pump arounds: an alternative to reactive distillation [J].Chemical Engineeringand Processing, 2004, 43: 435-445.

    [15] Ye J C, Li J, Sha Y, Lin H D, Zhou D W. Evaluation of reactive distillation and side reactor configuration for direct hydration of cyclohexene to cyclohexanol [J].Industrial & Engineering Chemistry Research, 2014, 53: 1461-1469.

    [16] Bisowarno B H, Tian Y C, Tade M O. Application of side reactors on ETBE reactive distillation [J].Chemical Engineering Journal, 2004, 99(1): 35-43.

    [17] Kaymak D B, Luyben W L. Design of distillation columns with external side reactors [J].Industrial & Engineering Chemistry Research, 2004, 43: 8049-8056.

    [18] Kaymak D B, Luyben W L. Optimum design of a column/side reactor process [J].Industrial & Engineering Chemistry Research, 2007, 46(15): 5175-5185.

    [19] Gangadwala J, Kienle A. MINLP optimization of butyl acetate synthetsis [J].Chemical Engineeringand Processing, 2007, 46: 107-118.

    [20] Ding Lianghui(丁良輝), Tang Jihai(湯吉海), Cui Mifen(崔咪芬), Chen Xian(陳獻(xiàn)), Bo Cuimei(薄翠梅), Qiao Xu(喬旭). Optimum design of coupled atmospheric reaction—vacuum distillationtechnology for benzyl chloride production [J].CIESC Journal(化工學(xué)報(bào)), 2011, 62(8): 2323-2327.

    [21] Ding L H, Tang J H, Cui M F, Bo C M, Chen Xian, Qiao X. Optimum design and analysis based on independent reaction amount for distillation column with side reactors: production of benzyl chloride [J].Industrial & Engineering Chemistry Research, 2011, 50(19): 11143-11152.

    [22] Kim Y, Pokki J P, Rehak K, Matous J, Aittamaa J. Excess enthalpy and vapor-liquid equilibrium at atmospheric pressure for the binary systems 2-methyl-2-propanol+2,4,4-trimethyl-1-pentene and 2-butanol+ 2,4,4-trimethyl-1-pentene [J].Journal of Chemical and Engineering Data, 2004, 49(2): 246-250.

    [23] Douglas J M. Conceptual Design of Chemical Processes [M]. New York: McGraw-Hill, 1988: 345-350.

    [24] Lu Enxi(陸恩錫), Li Xiaoling(李小玲). Retrograde distillation [J].Chemical Engineering(China) (化學(xué)工程), 2010, 38(4): 97-102.

    Process simulation oftert-butyl acetate production by distillation column integrated with side reactors under different temperatures

    HUANG Yuxin1, TANG Jihai1, CHEN Xian1, SONG Jianguo1,CUI Mifen1, FEI Zhaoyang2, QIAO Xu2
    (1Jiangsu Key Laboratory of Industrial Water-conservation & Emission Reduction,College of Chemistry and Chemical Engineering,Nanjing Tech University,Nanjing210009,Jiangsu,China;2State Key Laboratory of Materials-Oriented Chemical Engineering,Nanjing Tech University,Nanjing210009,Jiangsu,China)

    An advanced technology integrating a high-temperature distillation column with low-temperature side reactors(SRC) was proposed fortert-butyl acetate productionviathe additive esterification of acetic acid and isobutylene, since the reaction and separation required completely different operations. With a vapor boilup rate of 100 kmol·h?1, acetic acid conversion of 99.9% andtert-butyl acetate selectivity of 97.0%, the effects of flow rate from distillation column to reactors, stage number of rectification zone, side-reactor number and interval stage number between inlet and outlet of the reactor were systemically investigated. The results demonstrated that a higher economy process could be achieved with a more optimized match between the reaction capability and the separation capability. The obtained results laid a foundation for design and optimization of the industrial SRC process.

    Prof. TANG Jihai, jhtang@njtech.edu.cn; Prof. QIAO Xu, qct@njtech.edu.cn

    10.11949/j.issn.0438-1157.20150207

    TQ 021.8

    :A

    :0438—1157(2015)10—4039—08

    2015-02-09收到初稿,2015-06-07收到修改稿。

    聯(lián)系人:湯吉海,喬旭。

    :黃玉鑫(1990—),男,碩士研究生。

    國家自然科學(xué)基金項(xiàng)目(21276126,61203020);江蘇省普通高校研究生科研創(chuàng)新計(jì)劃項(xiàng)目(SJZZ_0098);江蘇高校優(yōu)勢學(xué)科建設(shè)工程項(xiàng)目。

    Received date: 2015-02-09.

    Foundation item: supported by the National Natural Science Foundation of China (21276126, 61203020), the Graduate Research and Innovation Projects in Jiangsu Province (SJZZ_0098) and the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

    猜你喜歡
    異丁烯板數(shù)丁酯
    二硫化碳-四氯化碳精餾工藝改進(jìn)與模擬研究
    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
    Global Bioenergies公司的生物異丁烯純度已達(dá)99.77%
    3,5-二氨基對氯苯甲酸異丁酯的合成研究
    乙酸仲丁酯的催化合成及分析
    山東魯清石化17萬t UOP異丁烯生產(chǎn)裝置順利投產(chǎn)
    上海化工(2015年12期)2015-04-07 17:07:14
    鄰苯二甲酸二丁酯的收縮血管作用及其機(jī)制
    2,4-D丁酯對棉花造成的藥害及補(bǔ)救措施
    土茯苓含量測定分離度和理論板數(shù)的調(diào)整
    激情视频va一区二区三区| 亚洲,欧美,日韩| 18禁国产床啪视频网站| 一本一本久久a久久精品综合妖精| 一二三四在线观看免费中文在| 国产成人免费无遮挡视频| 日韩av不卡免费在线播放| www.999成人在线观看| 精品少妇黑人巨大在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 晚上一个人看的免费电影| 好男人电影高清在线观看| 亚洲少妇的诱惑av| 午夜福利,免费看| kizo精华| 久久女婷五月综合色啪小说| 国产精品国产三级专区第一集| 操出白浆在线播放| 一本—道久久a久久精品蜜桃钙片| 可以免费在线观看a视频的电影网站| 精品少妇内射三级| 日韩中文字幕欧美一区二区 | 久久久久精品人妻al黑| 久久中文字幕一级| 日韩制服骚丝袜av| 天天躁日日躁夜夜躁夜夜| 亚洲av成人不卡在线观看播放网 | 99国产精品一区二区蜜桃av | 欧美97在线视频| 90打野战视频偷拍视频| 性少妇av在线| 一本综合久久免费| 搡老岳熟女国产| 成人亚洲欧美一区二区av| 国产精品久久久久久人妻精品电影 | 亚洲欧美一区二区三区国产| 天天影视国产精品| 欧美精品人与动牲交sv欧美| av片东京热男人的天堂| 我的亚洲天堂| 在线 av 中文字幕| 午夜福利一区二区在线看| 大香蕉久久网| 色精品久久人妻99蜜桃| 最新在线观看一区二区三区 | 亚洲成人免费av在线播放| h视频一区二区三区| 欧美日本中文国产一区发布| www.av在线官网国产| 午夜av观看不卡| 18禁裸乳无遮挡动漫免费视频| 成人影院久久| 在线 av 中文字幕| 国产在线一区二区三区精| 五月开心婷婷网| 欧美在线一区亚洲| 一本色道久久久久久精品综合| 中文欧美无线码| 大片电影免费在线观看免费| 欧美国产精品va在线观看不卡| 亚洲精品国产区一区二| 巨乳人妻的诱惑在线观看| 国产伦人伦偷精品视频| 在线观看免费午夜福利视频| 久久热在线av| 亚洲欧美清纯卡通| 一个人免费看片子| 涩涩av久久男人的天堂| 在线 av 中文字幕| 欧美性长视频在线观看| 丁香六月天网| 久久久久精品人妻al黑| 日本av免费视频播放| 黄色片一级片一级黄色片| 精品欧美一区二区三区在线| 国语对白做爰xxxⅹ性视频网站| 久久久久久免费高清国产稀缺| 无遮挡黄片免费观看| 日韩欧美一区视频在线观看| 大香蕉久久成人网| 国产亚洲午夜精品一区二区久久| 婷婷色综合www| 亚洲三区欧美一区| 亚洲男人天堂网一区| 满18在线观看网站| 性少妇av在线| 亚洲色图综合在线观看| 亚洲午夜精品一区,二区,三区| 国产精品香港三级国产av潘金莲 | 亚洲av电影在线进入| 精品久久久久久久毛片微露脸 | 亚洲伊人色综图| 国产成人精品在线电影| 国产国语露脸激情在线看| 搡老乐熟女国产| 一本综合久久免费| 国产无遮挡羞羞视频在线观看| 亚洲成av片中文字幕在线观看| 国产成人精品无人区| 日本wwww免费看| 青春草亚洲视频在线观看| 啦啦啦啦在线视频资源| 1024香蕉在线观看| 久久av网站| 欧美xxⅹ黑人| 操美女的视频在线观看| 国产1区2区3区精品| 每晚都被弄得嗷嗷叫到高潮| 91精品三级在线观看| 电影成人av| 日本vs欧美在线观看视频| 欧美国产精品一级二级三级| 丝袜在线中文字幕| 欧美变态另类bdsm刘玥| 男女下面插进去视频免费观看| 久9热在线精品视频| 夫妻午夜视频| 亚洲欧美清纯卡通| 日本91视频免费播放| 免费看不卡的av| 久久精品亚洲av国产电影网| 中文字幕最新亚洲高清| 亚洲精品日本国产第一区| 婷婷成人精品国产| av有码第一页| av网站免费在线观看视频| 国产一区二区三区av在线| 中文字幕人妻丝袜制服| 欧美人与性动交α欧美精品济南到| 午夜激情久久久久久久| av电影中文网址| av天堂在线播放| 亚洲av日韩精品久久久久久密 | 各种免费的搞黄视频| 国产男女内射视频| 久热爱精品视频在线9| 久久久国产欧美日韩av| 丝袜美腿诱惑在线| 美女主播在线视频| 丝袜美足系列| 黑人巨大精品欧美一区二区蜜桃| 成人免费观看视频高清| 嫁个100分男人电影在线观看 | 国产精品久久久人人做人人爽| 国产av精品麻豆| 这个男人来自地球电影免费观看| 国产高清国产精品国产三级| 建设人人有责人人尽责人人享有的| 精品一区二区三卡| 肉色欧美久久久久久久蜜桃| 99热网站在线观看| 丰满少妇做爰视频| h视频一区二区三区| 大型av网站在线播放| 99热网站在线观看| 精品久久久精品久久久| 女人高潮潮喷娇喘18禁视频| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 我要看黄色一级片免费的| 999精品在线视频| 叶爱在线成人免费视频播放| 一级毛片 在线播放| 丁香六月欧美| 国产在线一区二区三区精| 永久免费av网站大全| 最新在线观看一区二区三区 | 性色av一级| 成人免费观看视频高清| av片东京热男人的天堂| 一本久久精品| 在线精品无人区一区二区三| 在线观看免费午夜福利视频| 色播在线永久视频| 国产又爽黄色视频| 免费在线观看完整版高清| 永久免费av网站大全| 久久精品国产综合久久久| 日韩视频在线欧美| 免费观看av网站的网址| 女人精品久久久久毛片| 中文字幕av电影在线播放| 精品少妇一区二区三区视频日本电影| 男女床上黄色一级片免费看| 操出白浆在线播放| 久久久久网色| 波多野结衣av一区二区av| 久久亚洲国产成人精品v| h视频一区二区三区| 日本欧美视频一区| 国产爽快片一区二区三区| 久久精品国产a三级三级三级| 在线观看人妻少妇| 久久狼人影院| 亚洲人成网站在线观看播放| 啦啦啦 在线观看视频| 国产有黄有色有爽视频| 久久精品国产亚洲av高清一级| 中国美女看黄片| 性色av乱码一区二区三区2| 欧美精品一区二区大全| 成人影院久久| 晚上一个人看的免费电影| 99国产综合亚洲精品| 日韩制服骚丝袜av| 亚洲精品日本国产第一区| 建设人人有责人人尽责人人享有的| 成人亚洲欧美一区二区av| 日韩一卡2卡3卡4卡2021年| 亚洲国产最新在线播放| 国产成人系列免费观看| 中文字幕人妻熟女乱码| 真人做人爱边吃奶动态| 亚洲图色成人| 亚洲av国产av综合av卡| 久久人人爽av亚洲精品天堂| 久久久久久久国产电影| 亚洲国产欧美一区二区综合| 伦理电影免费视频| 亚洲精品成人av观看孕妇| 另类亚洲欧美激情| 久久久亚洲精品成人影院| 中文字幕av电影在线播放| 亚洲专区中文字幕在线| 欧美激情高清一区二区三区| 一级毛片女人18水好多 | 性色av乱码一区二区三区2| 国产麻豆69| 亚洲精品一卡2卡三卡4卡5卡 | 又黄又粗又硬又大视频| 少妇精品久久久久久久| 国产男女超爽视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲中文日韩欧美视频| 国产精品麻豆人妻色哟哟久久| 嫁个100分男人电影在线观看 | av线在线观看网站| 天天操日日干夜夜撸| 欧美国产精品一级二级三级| 免费在线观看影片大全网站 | 99热国产这里只有精品6| 日本av手机在线免费观看| 老汉色∧v一级毛片| 久久久久精品人妻al黑| 免费看十八禁软件| 9191精品国产免费久久| 日韩电影二区| 日韩大码丰满熟妇| 2021少妇久久久久久久久久久| av在线老鸭窝| 久久精品成人免费网站| 我要看黄色一级片免费的| avwww免费| 1024视频免费在线观看| 777米奇影视久久| 老汉色∧v一级毛片| 视频区欧美日本亚洲| 亚洲伊人久久精品综合| 美女国产高潮福利片在线看| 高清视频免费观看一区二区| 视频区图区小说| 在线 av 中文字幕| 高清视频免费观看一区二区| 极品人妻少妇av视频| 免费一级毛片在线播放高清视频 | 久久免费观看电影| 成人影院久久| 国产又爽黄色视频| 亚洲精品久久久久久婷婷小说| 亚洲av成人精品一二三区| 亚洲五月色婷婷综合| 欧美黄色片欧美黄色片| 七月丁香在线播放| 我要看黄色一级片免费的| 国产亚洲欧美在线一区二区| 国产成人免费观看mmmm| 欧美日韩av久久| 一本一本久久a久久精品综合妖精| 在线观看免费午夜福利视频| 亚洲av片天天在线观看| 国产又爽黄色视频| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品成人av观看孕妇| 人人妻人人爽人人添夜夜欢视频| 日本wwww免费看| 婷婷成人精品国产| 欧美xxⅹ黑人| 人人妻人人添人人爽欧美一区卜| 午夜免费观看性视频| 国产一区二区三区av在线| 黑人欧美特级aaaaaa片| 午夜福利一区二区在线看| 欧美日韩视频精品一区| 黄色毛片三级朝国网站| 1024香蕉在线观看| 亚洲精品乱久久久久久| 久久人人爽人人片av| 天堂中文最新版在线下载| 男女高潮啪啪啪动态图| 三上悠亚av全集在线观看| 最新的欧美精品一区二区| 日日摸夜夜添夜夜爱| 好男人电影高清在线观看| 香蕉国产在线看| av电影中文网址| 80岁老熟妇乱子伦牲交| 久久精品亚洲av国产电影网| a级毛片黄视频| 色94色欧美一区二区| 久久综合国产亚洲精品| 欧美+亚洲+日韩+国产| 国产色视频综合| 久久精品国产综合久久久| 精品少妇黑人巨大在线播放| 国产精品免费大片| 亚洲专区中文字幕在线| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产最新在线播放| 另类亚洲欧美激情| 七月丁香在线播放| 亚洲成人手机| 9热在线视频观看99| 王馨瑶露胸无遮挡在线观看| 国产成人a∨麻豆精品| 超碰成人久久| 男女午夜视频在线观看| 国产免费现黄频在线看| 性色av乱码一区二区三区2| 国产国语露脸激情在线看| 一区二区三区四区激情视频| 九色亚洲精品在线播放| 另类精品久久| 天天添夜夜摸| 黄色片一级片一级黄色片| 好男人视频免费观看在线| 男人添女人高潮全过程视频| 男人操女人黄网站| 婷婷色av中文字幕| 欧美在线一区亚洲| 免费在线观看视频国产中文字幕亚洲 | www.熟女人妻精品国产| 中文字幕亚洲精品专区| 国产伦人伦偷精品视频| 久久影院123| 久久天堂一区二区三区四区| 在线av久久热| 国产视频一区二区在线看| 国产亚洲欧美在线一区二区| 欧美精品av麻豆av| 精品一区二区三卡| 国产成人91sexporn| 777久久人妻少妇嫩草av网站| 国产成人精品久久二区二区91| 午夜久久久在线观看| 亚洲午夜精品一区,二区,三区| 男女国产视频网站| 欧美性长视频在线观看| 精品一区二区三卡| 欧美精品一区二区大全| 国产高清videossex| 黄色视频在线播放观看不卡| 亚洲国产成人一精品久久久| 欧美激情极品国产一区二区三区| svipshipincom国产片| 丝瓜视频免费看黄片| 一级毛片黄色毛片免费观看视频| 丰满少妇做爰视频| 十八禁网站网址无遮挡| 丰满人妻熟妇乱又伦精品不卡| av视频免费观看在线观看| bbb黄色大片| 又紧又爽又黄一区二区| 夫妻性生交免费视频一级片| 中文字幕高清在线视频| 免费在线观看完整版高清| 日韩一卡2卡3卡4卡2021年| 日本五十路高清| 青草久久国产| 在线 av 中文字幕| 亚洲av欧美aⅴ国产| 亚洲成人手机| 欧美激情 高清一区二区三区| 国产片内射在线| 亚洲av在线观看美女高潮| 亚洲一区二区三区欧美精品| 亚洲专区国产一区二区| 满18在线观看网站| 咕卡用的链子| 免费观看av网站的网址| 国精品久久久久久国模美| 青春草视频在线免费观看| 亚洲国产中文字幕在线视频| 亚洲精品一区蜜桃| e午夜精品久久久久久久| 免费看av在线观看网站| 国产免费福利视频在线观看| 国产精品99久久99久久久不卡| 看十八女毛片水多多多| 各种免费的搞黄视频| 国产免费又黄又爽又色| 国产极品粉嫩免费观看在线| 国产黄色视频一区二区在线观看| 麻豆av在线久日| 婷婷色综合www| 亚洲人成网站在线观看播放| 欧美久久黑人一区二区| 久热这里只有精品99| 后天国语完整版免费观看| 丝袜在线中文字幕| 十八禁人妻一区二区| 天天躁日日躁夜夜躁夜夜| 精品免费久久久久久久清纯 | 又黄又粗又硬又大视频| 成人免费观看视频高清| 欧美大码av| 国产伦理片在线播放av一区| 国产亚洲午夜精品一区二区久久| 女人被躁到高潮嗷嗷叫费观| 只有这里有精品99| 亚洲精品久久午夜乱码| 国产成人精品在线电影| xxxhd国产人妻xxx| 人妻 亚洲 视频| 亚洲国产精品一区三区| 精品视频人人做人人爽| 日本猛色少妇xxxxx猛交久久| 精品人妻1区二区| 亚洲熟女精品中文字幕| 天堂俺去俺来也www色官网| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 免费观看av网站的网址| 午夜福利影视在线免费观看| 90打野战视频偷拍视频| 国产精品久久久久成人av| 国产亚洲欧美在线一区二区| 精品一区在线观看国产| 国产一区二区三区综合在线观看| 久久久精品94久久精品| 国产精品久久久av美女十八| 亚洲精品国产色婷婷电影| 少妇精品久久久久久久| 黄色片一级片一级黄色片| 丝袜喷水一区| 久久这里只有精品19| 久久久国产欧美日韩av| 国产97色在线日韩免费| 9热在线视频观看99| 91国产中文字幕| 亚洲午夜精品一区,二区,三区| 午夜福利免费观看在线| 亚洲中文日韩欧美视频| 亚洲av综合色区一区| 亚洲成av片中文字幕在线观看| 国产视频首页在线观看| 日韩欧美一区视频在线观看| 一级毛片黄色毛片免费观看视频| 国产激情久久老熟女| 国产不卡av网站在线观看| 亚洲伊人久久精品综合| 成年动漫av网址| 最新在线观看一区二区三区 | 性色av乱码一区二区三区2| 国产成人免费无遮挡视频| 国产一卡二卡三卡精品| 91精品国产国语对白视频| 尾随美女入室| 日韩一卡2卡3卡4卡2021年| 18禁观看日本| 亚洲国产精品一区二区三区在线| 国产成人精品在线电影| 超色免费av| 久久影院123| 中文字幕制服av| 免费女性裸体啪啪无遮挡网站| 丰满迷人的少妇在线观看| 欧美+亚洲+日韩+国产| 日本a在线网址| 国产成人啪精品午夜网站| 精品亚洲成国产av| av网站免费在线观看视频| 人人妻人人添人人爽欧美一区卜| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 一级毛片女人18水好多 | 国产欧美日韩一区二区三 | 国产熟女欧美一区二区| 国产视频首页在线观看| 国产成人精品在线电影| 欧美精品亚洲一区二区| 久9热在线精品视频| 一个人免费看片子| av天堂久久9| 久久久久久久大尺度免费视频| 女人精品久久久久毛片| 国产亚洲精品第一综合不卡| 免费观看a级毛片全部| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 欧美变态另类bdsm刘玥| 欧美另类一区| 久久久精品免费免费高清| 欧美日韩亚洲国产一区二区在线观看 | 大香蕉久久网| 一级毛片电影观看| 精品少妇黑人巨大在线播放| 两个人看的免费小视频| 午夜福利在线免费观看网站| 国产1区2区3区精品| 男女床上黄色一级片免费看| 亚洲自偷自拍图片 自拍| 一级毛片黄色毛片免费观看视频| 啦啦啦视频在线资源免费观看| 精品福利观看| 男女之事视频高清在线观看 | 亚洲精品久久午夜乱码| 免费高清在线观看日韩| 丝袜美足系列| 国产熟女欧美一区二区| 久久国产亚洲av麻豆专区| 亚洲,欧美,日韩| 久久青草综合色| 中文字幕人妻丝袜一区二区| 无限看片的www在线观看| 欧美日韩综合久久久久久| av国产精品久久久久影院| 最近最新中文字幕大全免费视频 | 亚洲av成人不卡在线观看播放网 | 老司机影院毛片| 中文字幕av电影在线播放| 欧美在线一区亚洲| 亚洲欧美一区二区三区国产| 伊人亚洲综合成人网| 每晚都被弄得嗷嗷叫到高潮| 日韩伦理黄色片| 秋霞在线观看毛片| 久久精品成人免费网站| 亚洲中文日韩欧美视频| 亚洲激情五月婷婷啪啪| 精品少妇一区二区三区视频日本电影| videosex国产| 91精品伊人久久大香线蕉| 男女边摸边吃奶| 一级a爱视频在线免费观看| 男的添女的下面高潮视频| 丝袜脚勾引网站| 国产精品 国内视频| 夜夜骑夜夜射夜夜干| 大型av网站在线播放| 又黄又粗又硬又大视频| 麻豆国产av国片精品| 中文字幕制服av| 久久久国产一区二区| 欧美黑人精品巨大| 午夜福利视频精品| 这个男人来自地球电影免费观看| 精品亚洲成a人片在线观看| 国产高清国产精品国产三级| 国产熟女欧美一区二区| 蜜桃国产av成人99| 高清av免费在线| 视频区欧美日本亚洲| 亚洲精品成人av观看孕妇| 女人精品久久久久毛片| 中文字幕人妻丝袜一区二区| 亚洲一区中文字幕在线| 少妇人妻 视频| 18禁国产床啪视频网站| 欧美人与性动交α欧美软件| 别揉我奶头~嗯~啊~动态视频 | 国产伦人伦偷精品视频| 国产高清视频在线播放一区 | 亚洲欧美激情在线| 中文字幕色久视频| 国产三级黄色录像| 日韩 亚洲 欧美在线| 国产精品秋霞免费鲁丝片| 久久亚洲国产成人精品v| 超碰成人久久| 日韩伦理黄色片| 久久天堂一区二区三区四区| 我要看黄色一级片免费的| av福利片在线| 蜜桃在线观看..| 欧美性长视频在线观看| 女警被强在线播放| a级片在线免费高清观看视频| 最近手机中文字幕大全| 99精品久久久久人妻精品| 国产午夜精品一二区理论片| 中文字幕精品免费在线观看视频| 香蕉国产在线看| 午夜福利乱码中文字幕| 天天躁日日躁夜夜躁夜夜| 成人黄色视频免费在线看| 国产又爽黄色视频| 制服人妻中文乱码| 满18在线观看网站| 黄色视频在线播放观看不卡| 两人在一起打扑克的视频| 晚上一个人看的免费电影| 午夜激情久久久久久久| 免费在线观看完整版高清| 青青草视频在线视频观看| 欧美精品人与动牲交sv欧美| 亚洲国产最新在线播放| 免费人妻精品一区二区三区视频| 国产成人系列免费观看| 日本一区二区免费在线视频| 尾随美女入室| 亚洲国产精品成人久久小说| 亚洲av欧美aⅴ国产| 91成人精品电影| 免费一级毛片在线播放高清视频 | 亚洲欧美一区二区三区久久|