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

    共沸精餾隔壁塔與萃取精餾隔壁塔的控制研究

    2017-02-17 12:55:31劉立新陳夢(mèng)琪劉育良王建新孫蘭義
    化工進(jìn)展 2017年2期
    關(guān)鍵詞:體系

    劉立新,陳夢(mèng)琪,劉育良,王建新,孫蘭義

    (1青島理工大學(xué)汽車與交通學(xué)院,山東 青島 266520;2中國(guó)石油大學(xué)化學(xué)(華東)工程學(xué)院,山東 青島 266580;3愛因霍芬理工大學(xué)化學(xué)工程學(xué)院,荷蘭)

    共沸精餾隔壁塔與萃取精餾隔壁塔的控制研究

    劉立新1,陳夢(mèng)琪2,劉育良3,王建新2,孫蘭義2

    (1青島理工大學(xué)汽車與交通學(xué)院,山東 青島 266520;2中國(guó)石油大學(xué)化學(xué)(華東)工程學(xué)院,山東 青島 266580;3愛因霍芬理工大學(xué)化學(xué)工程學(xué)院,荷蘭)

    針對(duì)多元共沸物或近沸點(diǎn)混合物的分離,采用共沸精餾隔壁塔和萃取精餾隔壁塔兩種流程,分別建立穩(wěn)態(tài)模型,并進(jìn)行了溫度靈敏板的選擇。針對(duì)共沸精餾隔壁塔建立若干兩點(diǎn)溫度控制結(jié)構(gòu),針對(duì)萃取精餾隔壁塔建立若干三點(diǎn)(及四點(diǎn))溫度控制結(jié)構(gòu)。通過添加進(jìn)料流量和組成擾動(dòng)進(jìn)行測(cè)試分析,分別為兩種流程挑選了能有效抵抗進(jìn)料擾動(dòng)的溫度控制結(jié)構(gòu)。①共沸精餾隔壁塔最優(yōu)控制結(jié)構(gòu):QMC/F控制TMC,13;QRC/F控制TRC,5。②萃取精餾隔壁塔最優(yōu)控制結(jié)構(gòu):RRM控制TMC1,3;Qr/F控制TMC1,12;RRR控制TRC1,3;αv控制TMC1,9。最后通過分析兩種最優(yōu)控制結(jié)構(gòu)的相似性,總結(jié)得出:帶有再沸器與進(jìn)料量比值(Qr/F)控制的溫度控制結(jié)構(gòu),可有效降低共沸精餾隔壁塔及萃取精餾隔壁塔體系的余差及超調(diào)量。

    共沸物;分離;共沸精餾隔壁塔;萃取精餾隔壁塔;控制

    精餾是化工過程中能耗巨大的單元操作,為了降低精餾過程的能耗可采用兩種技術(shù)手段:一是熱耦合技術(shù)[1],如隔壁塔(dividing wall column,DWC);二是精餾過程與其他過程的集成。DWC技術(shù)[2-3]將用于多組分分離過程的多塔集成在一個(gè)塔中,這種精餾塔集成形式能夠有效降低設(shè)備投資、操作費(fèi)用以及占地面積[4-7]。對(duì)于多元共沸物或近沸點(diǎn)混合物的分離,工業(yè)上常采用共沸精餾和萃取精餾法,通過向混合液中加入第三種組分(夾帶劑),提高待分離組分之間的相對(duì)揮發(fā)度,使其得以分離。

    共沸精餾隔壁塔(azeotropic dividing wall column,ADWC)將共沸精餾過程[8]與DWC的設(shè)計(jì)概念相結(jié)合,如圖1(a)所示。迄今為止,對(duì)ADWC的研究報(bào)道已有不少,WU等[9]利用甲苯分離吡啶和水的體系對(duì)ADWC的完全熱耦合流程進(jìn)行研究,相對(duì)于傳統(tǒng)分離過程,其能耗降低29%,而年總費(fèi)用(TAC)降低31%。BRIONES-RAMíREZ等[10]同樣采用Aspen Plus軟件,選取異丙醇/水/丙酮體系和異丙醇/水/甲醇體系對(duì)ADWC流程進(jìn)行了理論研究,通過引入多目標(biāo)遺傳算法對(duì)體系進(jìn)行優(yōu)化,以獲得最優(yōu)的穩(wěn)態(tài)操作條件,結(jié)果表明,該ADWC流程相對(duì)于傳統(tǒng)的兩塔精餾過程可節(jié)約能耗約50%。而對(duì)于ADWC的工業(yè)應(yīng)用僅有KAIBEL等[11]的文章提及。

    萃取精餾隔壁塔(extractive dividing wall column,EDWC)指將萃取精餾過程的兩個(gè)精餾塔耦合成一個(gè)帶隔板的精餾塔,其集成過程如圖1(b)所示。如今EDWC在工業(yè)應(yīng)用和理論研究方面均取得了較大進(jìn)展:Uhde公司將隔壁塔應(yīng)用于嗎啉法萃取精餾過程,該過程基于不同進(jìn)料純度的苯/甲苯/二甲苯的兩塔分離流程[12-13];MIDORI等[14]利用Aspen Plus軟件從理論上對(duì)EDWC進(jìn)行研究,選取丙酮/甲醇/水體系,結(jié)果顯示與常規(guī)萃取精餾流程相比,其節(jié)省能耗約36%。

    然而,由于DWC將多塔合并為一塔,其內(nèi)部變量增多,且變量間耦合性較強(qiáng)[15],為其控制研究帶來一定困難,其在實(shí)際應(yīng)用中難以操作和控制,成為制約DWC工業(yè)推廣的因素之一。目前,學(xué)者對(duì)ADWC以及EDWC的控制研究還很少,對(duì)于不同的體系得出的結(jié)果也不相同。因此,為了探索不同DWC流程的控制規(guī)律,本文針對(duì)ADWC和EDWC兩種流程,分別建立穩(wěn)態(tài)模型,之后進(jìn)行系統(tǒng)的動(dòng)態(tài)特性分析,為其選擇適宜的控制結(jié)構(gòu),并分析控制結(jié)構(gòu)的相似性。

    圖1 從精餾塔到ADWC及 EDWC的耦合過程

    1 ADWC控制

    1.1 ADWC穩(wěn)態(tài)模型

    異丙醇(IPA)廣泛用作半導(dǎo)體生產(chǎn)過程的清潔介質(zhì),也常用于家庭、醫(yī)院以及工業(yè)中的殺菌消毒。異丙醇廢液通常是IPA和水的混合物,脫水可以實(shí)現(xiàn)IPA的重復(fù)回收利用,但I(xiàn)PA和水易形成共沸物(共沸溫度80℃,87.4%,IPA摩爾分?jǐn)?shù))[16],采用常規(guī)精餾手段分離難度大、能耗高、設(shè)備投資高。目前,分離IPA-水體系的工藝是用苯作為夾帶劑的共沸精餾法,但由于苯對(duì)人體的危害,近期共沸精餾法生產(chǎn)IPA開始采用環(huán)己烷為夾帶劑[17]。

    本文選取環(huán)己烷為共沸劑,將ADWC應(yīng)用于IPA脫水過程,其示意圖見圖2(a)。混合物進(jìn)料流率為100kmol/h;進(jìn)料摩爾組成為50%的IPA和50%的H2O;進(jìn)料為飽和液體,壓力為101.3kPa。此外塔的操作壓力為101.3kPa,每塊理論板壓降為1kPa(文中所述壓降均為理論板壓降)。要求IPA產(chǎn)品中IPA摩爾分率不低于99.99%,廢水中水的摩爾分?jǐn)?shù)不低于99.9%。

    圖2 ADWC示意圖及其在Aspen Plus 中的替代模型

    利用Aspen Plus軟件,采用NRTL物性方法[18],選取一個(gè)完整精餾塔(MC)和一個(gè)提餾塔(RC)組成的熱耦合精餾模型進(jìn)行熱力學(xué)等效模擬,如圖2(b)所示?;旌衔镉蒑C進(jìn)料,IPA產(chǎn)品從MC塔釜采出,IPA、水和環(huán)己烷組成的共沸物由MC塔頂餾出,經(jīng)過冷凝后分為兩相:一相為有機(jī)相,含有較多的環(huán)己烷,作為回流返回MC塔頂;另一相為水相,引入RC進(jìn)行進(jìn)一步分離。廢水從RC塔釜采出。該ADWC包含兩個(gè)再沸器,比常規(guī)精餾過程多一個(gè),因此在控制過程中可以選取的操縱變量比常規(guī)精餾過程也要多一個(gè)。穩(wěn)態(tài)模擬的ADWC操作參數(shù)見表1。

    在建立控制結(jié)構(gòu)之前需要對(duì)共沸精餾隔壁塔的穩(wěn)態(tài)模型添加必要的泵和閥門,在這里閥門壓降設(shè)置為304kPa。由穩(wěn)態(tài)模型可計(jì)算出MC塔徑為2.6m,RC塔徑為1.5m,除此之外,動(dòng)態(tài)模擬時(shí)還需要計(jì)算塔釜和回流罐的尺寸。在計(jì)算過程中我們規(guī)定當(dāng)流體占設(shè)備體積50%時(shí),其停留時(shí)間為5 min。計(jì)算得到的設(shè)備尺寸記錄在表2中。

    表1 穩(wěn)態(tài)模擬的ADWC操作參數(shù)

    表2 塔釜及回流罐參數(shù)

    1.2 ADWC動(dòng)態(tài)控制

    本文選取傳統(tǒng)的PID控制策略,通過建立多回路溫度控制結(jié)構(gòu)對(duì)ADWC進(jìn)行控制研究。從理論上講,ADWC的控制比DWC的控制簡(jiǎn)單,因?yàn)锳DWC不涉及隔壁底部氣相分配比(αv)的穩(wěn)定與調(diào)節(jié)問題,故首先考慮選取兩個(gè)溫度控制板,通過調(diào)節(jié)對(duì)應(yīng)操縱變量來維持被控塔板溫度恒定,進(jìn)而保證產(chǎn)品純度。在建立溫度控制結(jié)構(gòu)之前首先要選取溫度靈敏板。

    采用靈敏板判據(jù)選擇溫度靈敏板,該ADWC中能夠選取的操縱變量有MC再沸器熱負(fù)荷(QMC)、RC再沸器熱負(fù)荷(QRC)以及隔板上方液相分離比[見圖2(a),RLS=LS/L]。然而,在該體系中僅需要控制兩個(gè)溫度靈敏板,因此需要對(duì)不同的操縱變量與溫度靈敏板進(jìn)行匹配,對(duì)比選擇最優(yōu)的匹配策略。首先分別將3個(gè)操縱變量QMC、QRC和RLS在穩(wěn)態(tài)設(shè)計(jì)值基礎(chǔ)上增加、減小0.1%記錄各板溫度開環(huán)增益,如圖3。

    根據(jù)圖3(a)和圖3(b)可以看出,QMC對(duì)應(yīng)的溫度靈敏板為MC的第13塊板(TMC,13),而QRC對(duì)應(yīng)兩塊溫度靈敏板,分別是TMC,13和TRC,5,根據(jù)就近原則,利用QMC來控制TMC,13,利用QRC來控制TRC,5。根據(jù)圖3(c)可以看出,RLS對(duì)應(yīng)兩塊溫度控制板,分別為TMC,11和TRC,5。由于僅需要控制兩個(gè)溫度靈敏板的溫度,且如果兩個(gè)溫度控制板距離較近會(huì)造成兩個(gè)對(duì)應(yīng)的控制回路之間的相互干擾,因此綜合上述條件可以得到3種溫度控制匹配方案,分別為:①RLS控制TMC,11,QRC控制TRC,5;②QMC控制TMC,13,RLS控制TRC,5;③QMC控制TMC,13,QRC控制TRC,5。

    本文根據(jù)3種溫度控制方案,建立了不同的控制結(jié)構(gòu),并分別對(duì)其控制效果進(jìn)行研究,分析對(duì)比不同操縱變量在系統(tǒng)遇到干擾時(shí)的控制效果。

    圖3QMC、QRC、RLS變化±0.1%對(duì)應(yīng)各塔板溫度的靈敏度曲線

    1.2.1 溫度控制結(jié)構(gòu)CS1

    將Aspen Plus中建立的穩(wěn)態(tài)模型導(dǎo)入Aspen Dynamics中,運(yùn)行初始化后利用選定的溫度匹配控制方案一,建立如圖4所示的控制結(jié)構(gòu)CS1?;究刂平Y(jié)構(gòu)如下:①流量控制器控制進(jìn)料流量;②塔頂氣相流率控制MC塔頂壓力;③有機(jī)相采出量控制回流罐中有機(jī)相液位,水相液位由回流罐流出的水相流率控制;④MC、RC塔釜產(chǎn)品采出量分別控制兩塔塔釜液位;⑤補(bǔ)充共沸劑的流量與進(jìn)料量成比例,補(bǔ)充共沸劑的流量控制器設(shè)置為串級(jí)控制,串級(jí)信號(hào)來源于比例控制器的輸出信號(hào)。

    溫度控制結(jié)構(gòu)如下:①Q(mào)RC控制TRC,5;②RLS控制TMC,11。QMC保持恒定,RC上方的液相回流管線設(shè)置流量控制器,將控制器設(shè)置為串級(jí)控制,輸入信號(hào)來自比值控制器的輸出,而比值控制器的輸入信號(hào)來自MC第11塊塔板的溫度控制器輸出。MC塔頂冷凝器溫度會(huì)影響有機(jī)相和水相中的組分分配,從而對(duì)共沸精餾分離效果造成影響,因此對(duì)共沸精餾塔而言,控制塔頂冷凝器出口溫度非常重要,故在冷凝器出口管線設(shè)置溫度控制器,控制器輸出信號(hào)連接冷凝器負(fù)荷。由于溫度控制存在時(shí)間滯后,因此需要在所有溫度控制器的輸入信號(hào)端添加死時(shí)間原件,值為1min,然后運(yùn)行初始化。

    控制結(jié)構(gòu)建立后采用繼電-反饋測(cè)試整定溫度控制器參數(shù),然后通過對(duì)體系施加擾動(dòng)來測(cè)試控制結(jié)構(gòu)的有效性。在閉環(huán)動(dòng)態(tài)平穩(wěn)運(yùn)行0.5h后,對(duì)體系施加±10%進(jìn)料量(F)和±5%水的進(jìn)料摩爾分?jǐn)?shù)(ZFH2O)擾動(dòng),記錄IPA產(chǎn)品純度(ZIPA)和廢水純度(ZH2O)的響應(yīng)曲線,結(jié)果表明,ZH2O在F改變時(shí)均能在短時(shí)間內(nèi)回歸設(shè)定值附近,而在F降低10%時(shí)ZIPA不能得到有效控制,產(chǎn)品純度在20h內(nèi)不能達(dá)到穩(wěn)定更不能回歸設(shè)定值。ZH2O在ZFH2O增加和減少時(shí)下均能在短時(shí)間回歸到設(shè)定值附近,但當(dāng)ZFH2O增加5%時(shí),ZIPA雖然能回歸穩(wěn)定但由0.99995變?yōu)?.992,余差較大。由此可見溫度控制結(jié)構(gòu)CS1控制效果不理想。

    1.2.2 溫度控制結(jié)構(gòu)CS2

    溫度控制結(jié)構(gòu)如下:①Q(mào)MC控制TMC,13;②RLS控制TRC,5。QMC保持恒定,保留RC上方液相回流管線上的流量控制器和比值控制器,將比值控制器的輸入信號(hào)改為TRC,5。其余控制回路與控制結(jié)構(gòu)CS1相同??刂平Y(jié)構(gòu)CS2建立后采用繼電-反饋測(cè)試重新整定溫度控制器參數(shù),然后通過對(duì)體系施加擾動(dòng)來測(cè)試控制結(jié)構(gòu)的有效性。

    結(jié)果顯示,在該控制結(jié)構(gòu)的調(diào)節(jié)作用下,無論F增加或減少,ZH2O和ZIPA均可在較短的時(shí)間內(nèi)達(dá)到穩(wěn)定,但當(dāng)F減少10%時(shí),ZH2O由0.999變?yōu)?.995,余差較大,造成較大的IPA損耗。由此可見根據(jù)溫度控制結(jié)構(gòu)CS2控制效果雖有所改進(jìn),但仍然不能滿足生產(chǎn)過程對(duì)產(chǎn)品及廢水的要求。

    1.2.3 溫度控制結(jié)構(gòu)CS3

    溫度控制如下:①Q(mào)MC控制TMC,13;②QRC控制TRC,5。RLS保持恒定,保留RC上方液相回流管線的流量控制器和比值控制器,但是比例控制器的比例不再是任何控制器的輸出值,而通過手動(dòng)輸入,且其數(shù)值不再隨著系統(tǒng)響應(yīng)發(fā)生改變,其余控制回路與前兩種控制結(jié)構(gòu)相同。平穩(wěn)運(yùn)行0.5h后,對(duì)體系施加相同的擾動(dòng),記錄ZH2O和ZIPA的響應(yīng)曲線。

    結(jié)果顯示,在溫度控制結(jié)構(gòu)CS3調(diào)節(jié)下,無論F擾動(dòng)還是ZFH2O發(fā)生變化,ZH2O和ZIPA均可以在較短的時(shí)間內(nèi)恢復(fù)穩(wěn)定,且余差穩(wěn)定在可接受范圍內(nèi)。由此可見,在以上3種控制方案中,通過調(diào)節(jié)兩塔再沸器負(fù)荷來穩(wěn)定靈敏板溫度的控制方案CS3效果最佳,然而該控制結(jié)構(gòu)也存在一定缺陷,當(dāng)F增加10%時(shí),ZIPA動(dòng)態(tài)響應(yīng)曲線的超調(diào)量較大,高達(dá)1.0%,而當(dāng)F減小10%時(shí),ZH2O動(dòng)態(tài)響應(yīng)曲線的超調(diào)量高達(dá)2.1%。

    圖4 Aspen Dynamics 中建立的控制結(jié)構(gòu)和控制器面板

    1.2.4 改進(jìn)的溫度控制結(jié)構(gòu)CS3-1

    溫度控制結(jié)構(gòu)CS3的響應(yīng)過程中超調(diào)量較大,主要是因?yàn)樵俜衅鳠嶝?fù)荷(Qr)在進(jìn)料量(F)發(fā)生變化時(shí)不能迅速作出響應(yīng)。為克服這一困難,考慮采用兩塔再沸器熱負(fù)荷與進(jìn)料量的比例(Qr/F)控制:將MC第13塊塔板的溫度控制器和RC第5塊塔板溫度控制器的輸出信號(hào)連接到比值控制器,作為比值控制器的比例的設(shè)定值,其余控制回路與溫度控制結(jié)構(gòu)CS3相同,如圖5所示。采用繼電-反饋測(cè)試重新整定溫度控制器參數(shù),然后通過對(duì)體系施加擾動(dòng)來測(cè)試控制結(jié)構(gòu)的有效性。在動(dòng)態(tài)平穩(wěn)運(yùn)行0.5h后,對(duì)體系施加10%的F擾動(dòng)和5%的ZFH2O擾動(dòng),記錄ZIPA、ZH2O、兩個(gè)靈敏板溫度和兩產(chǎn)品流率的響應(yīng)曲線,如圖6和圖7。

    由圖6可以看出,在改進(jìn)控制結(jié)構(gòu)的調(diào)節(jié)作用下,當(dāng)F發(fā)生擾動(dòng)時(shí),兩個(gè)溫度靈敏板的溫度都能在較短的時(shí)間內(nèi)回歸設(shè)定值;IPA流量(FIPA)和水流量(FH2O)能夠隨著F的變化而C呈比例地發(fā)生變化;ZIPA和ZH2O均可穩(wěn)定在較高水平,余差明顯降低,且超調(diào)量相對(duì)于控制結(jié)構(gòu)CS3明顯降低,當(dāng)F增加10%時(shí),ZIPA的超調(diào)量?jī)H為0.03%,而當(dāng)F降低10%時(shí),ZH2O的超調(diào)量?jī)H為0.08%。由圖7可以看出,在改進(jìn)控制結(jié)構(gòu)CS3-1的調(diào)節(jié)作用下,當(dāng)ZFH2O發(fā)生擾動(dòng)時(shí)(進(jìn)料中IPA濃度也對(duì)應(yīng)發(fā)生變化),兩個(gè)溫度靈敏板的溫度都能在較短的時(shí)間內(nèi)回歸設(shè)定值;FIPA、FH2O能夠隨著進(jìn)料中組分含量的的變化而成比例地發(fā)生變化;ZIPA和ZH2O均在短時(shí)間內(nèi)回歸穩(wěn)定,且余差較小。最終結(jié)果顯示:對(duì)于ADWC,帶有Qr/F比值控制的改進(jìn)控制結(jié)構(gòu)CS3-1控制效果最優(yōu)。

    圖5 改進(jìn)的溫度控制結(jié)構(gòu)CS3-1

    圖6 控制結(jié)構(gòu)CS3-1在±10%F擾動(dòng)下的動(dòng)態(tài)響應(yīng)曲線

    圖7 控制結(jié)構(gòu)CS3-1在±5%ZFH2O摩爾分率擾動(dòng)下的動(dòng)態(tài)響應(yīng)曲線

    2 EDWC控制

    2.1 EDWC穩(wěn)態(tài)模型

    萃取精餾分離技術(shù)工業(yè)應(yīng)用廣泛,對(duì)于乙醇-水體系來說一種有效的分離方法,故在萃取精餾的基礎(chǔ)上,本文將EDWC應(yīng)用于乙醇-水的分離過程,乙二醇為萃取劑。EDWC中隔板兩側(cè)共用一個(gè)提餾段和一個(gè)再沸器,隔板的上端與塔壁相連,左右兩側(cè)各有一個(gè)冷凝器,如圖8(a)所示,因此在EDWC中沒有隔板上方的液相分離卻有隔板下方的氣相分離,在精餾塔的設(shè)計(jì)階段可通過調(diào)整隔板下方的氣相分離比[見圖8(a),αv=VS/V]幫助滿足產(chǎn)品純度。其等效的完全熱耦合精餾塔分為一個(gè)主塔(MC1)和一個(gè)只有冷凝器的精餾塔(RC1),利用Aspen Plus軟件對(duì)體系建立穩(wěn)態(tài)模型,如圖8(b)所示。

    在該流程中,混合物由MC1進(jìn)料,進(jìn)料摩爾組成為84%的乙醇和16 %的水,進(jìn)料量45.36 kmol/h,溫度351K,壓力101.3kPa。MC1塔頂采出產(chǎn)品乙醇,要求乙醇產(chǎn)品純度不低于99.5%(摩爾),水從RC1塔頂采出,水的摩爾分?jǐn)?shù)不低于99.5%。冷凝器壓力為101.3kPa,每塊理論板壓降為1kPa。穩(wěn)態(tài)模擬的EDWC操作參數(shù)如表3。

    控制結(jié)構(gòu)建立前的準(zhǔn)備工作與ADWC相似,計(jì)算得到的MC1塔徑為1.2m,RC1塔徑為0.8m,回流罐與塔釜設(shè)備尺寸記錄在表4中。

    圖8 EDWC示意圖及其在Aspen Plus 中的替代模型

    表3 穩(wěn)態(tài)模擬的萃取精餾隔壁塔操作參數(shù)

    表4 塔釜及回流罐參數(shù)

    2.2 EDWC動(dòng)態(tài)控制

    本文選取傳統(tǒng)的PID控制策略,建立溫度控制結(jié)構(gòu)對(duì)EDWC進(jìn)行控制研究。由于塔內(nèi)部結(jié)構(gòu)復(fù)雜,變量間耦合作用嚴(yán)重,采用兩點(diǎn)溫度控制結(jié)構(gòu)來穩(wěn)定兩股產(chǎn)品的純度控制效果往往不佳,因此首先考慮選取3個(gè)溫度控制板,通過調(diào)節(jié)對(duì)應(yīng)操縱變量來維持被控板溫度恒定,最后達(dá)到保證產(chǎn)品純度目的。

    2.2.1 基礎(chǔ)控制結(jié)構(gòu)CS1

    與ADWC類似,建立溫度控制結(jié)構(gòu)之前首先對(duì)溫度靈敏板進(jìn)行選取。該體系中可以選作溫度靈敏板的操縱變量有MC1回流比(RRM)、RC1回流比(RRR)以及再沸器負(fù)荷(Qr)。與ADWC溫度控制板的選取類似,根據(jù)相對(duì)增益矩陣判據(jù)的靈敏板選擇步驟,分別將3個(gè)操縱變量:RRM、RRR和Qr在穩(wěn)態(tài)設(shè)計(jì)值基礎(chǔ)上增加、減小0.1%記錄各板溫度相對(duì)增益,如圖9所示,結(jié)果發(fā)現(xiàn),RRM和Qr所對(duì)應(yīng)的塔板溫度變化峰值位置相同,均為MC1第3、第12塊塔板和RC1第3塊塔板,而RRR對(duì)應(yīng)的塔板溫度變化峰值位置分別為MC1第9、第12塊塔板以及RC1的第3塊塔板。根據(jù)就近控制原則,得到的溫度靈敏板匹配方案為:RRM控制TMC1,3、Qr控制TMC1,12、RRR控制TRC1,3。

    圖9RRM、RRR、Qr變化±0.1%對(duì)應(yīng)各塔板溫度的靈敏度曲線

    將Aspen Plus中建立的穩(wěn)態(tài)模型導(dǎo)入Aspen Dynamics,運(yùn)行初始化后,利用選定的溫度匹配控制方案建立如圖10所示的基礎(chǔ)控制結(jié)構(gòu)CS1?;究刂平Y(jié)構(gòu)如下:乙醇和水混合物進(jìn)料量由流量控制器控制;MC1和RC1塔頂壓力由冷凝器負(fù)荷控制;冷凝器液位由塔頂餾出量控制;RC1塔釜液位由塔釜采出量控制;MC1塔釜液位由補(bǔ)充萃取劑流量進(jìn)行控制。溫度控制結(jié)構(gòu)如下:①RRM控制TMC1,3,回流比比值控制器輸入信號(hào)來自溫度控制器的輸出;②Qr控制TMC1,12;③RRR控制TRC1,3,回流比比值控制器輸入信號(hào)來自溫度控制器的輸出。另外:①由于循環(huán)萃取劑溫度對(duì)EDWC操作影響較大,因此利用萃取劑循環(huán)管線上的冷凝器負(fù)荷穩(wěn)定冷凝器出口溫度;②在循環(huán)萃取劑管線上設(shè)置流量控制器,控制器輸入信號(hào)來自比例控制器的輸出值,以此來保證萃取劑的進(jìn)料量與乙醇/水混合物進(jìn)料量成比例;③保證隔板下方氣相分離比(αv)恒定,在MC1氣相采出管線閥門處設(shè)置流量控制器,控制器輸入信號(hào)來自αv比值控制器的的輸出。所有溫度控制器的輸入信號(hào)端均需添加死時(shí)間原件,輸入死時(shí)間值為1min,然后運(yùn)行初始化。

    圖10 基礎(chǔ)控制結(jié)構(gòu)CS1

    控制結(jié)構(gòu)建立后采用繼電-反饋測(cè)試整定溫度控制器參數(shù),然后通過對(duì)體系施加擾動(dòng)來測(cè)試控制結(jié)構(gòu)的有效性。在閉環(huán)動(dòng)態(tài)平穩(wěn)運(yùn)行0.5h后,對(duì)體系施加±10%F擾動(dòng)和±5%進(jìn)料乙醇(ZFEtOH)擾動(dòng),記錄產(chǎn)品純度、靈敏板溫度以及對(duì)應(yīng)操縱變量的響應(yīng)曲線,如圖11所示,結(jié)果表明,在基礎(chǔ)控制結(jié)構(gòu)CS1的調(diào)節(jié)作用下,F(xiàn)擾動(dòng)時(shí)系統(tǒng)可在短時(shí)間內(nèi)回歸穩(wěn)定狀態(tài),余差不大但水純度(ZH2O)超調(diào)量較大,F(xiàn)增加10%時(shí)為0.1%,降低10%時(shí)為0.08%;而在ZFEtOH發(fā)生±5%擾動(dòng)時(shí),ZH2O的余差稍大。

    圖11 基礎(chǔ)控制結(jié)構(gòu)CS1在±10%F擾動(dòng)下的動(dòng)態(tài)響應(yīng)曲線

    2.2.2 含αv的控制結(jié)構(gòu)CS2

    為了克服基礎(chǔ)控制結(jié)構(gòu)CS1在抵抗±5%的ZFEtOH擾動(dòng)時(shí)ZH2O余差稍大的現(xiàn)象,考慮添加αv作為操縱變量的溫度控制結(jié)構(gòu)。將MC1氣相采出量在穩(wěn)態(tài)設(shè)計(jì)值基礎(chǔ)上增加、減小0.1%記錄各板溫度相對(duì)增益,結(jié)果顯示,MC1氣相采出量變化對(duì)應(yīng)的溫度相對(duì)增益峰值在MC1第10、第12以及RC1第3塊塔板。由于MC1第10塊塔板與再沸器負(fù)荷的溫度控制板位置較近,兩控制回路會(huì)發(fā)生干擾,以此選擇峰值附近的MC1第9塊塔板作為αv的溫度控制板。將αv控制器的輸入設(shè)定為MC1第9塊塔板的溫度控制器的輸出值,其他控制回路保持不變。

    對(duì)體系施加10%F擾動(dòng)和5%ZFEtOH擾動(dòng),記錄產(chǎn)品純度、靈敏板溫度以及對(duì)應(yīng)操縱變量的響應(yīng)曲線,結(jié)果顯示,體系抵抗±10%F擾動(dòng)的能力并未提高,但體系的余差有所降低。但其對(duì)控制效果的改善作用并不顯著。

    2.2.3 帶有再沸器熱負(fù)荷與進(jìn)料量比值(Qr/F)的控制結(jié)構(gòu)CS3

    為了研究Qr/F控制結(jié)構(gòu)在穩(wěn)定EDWC產(chǎn)品純度過程中所發(fā)揮的作用,本節(jié)在含αv的控制結(jié)構(gòu)CS2的基礎(chǔ)上,將Qr控制TMC1,12的控制回路改為Qr/F控制TMC1,12,控制結(jié)構(gòu)如圖12所示。對(duì)體系施加同樣的進(jìn)料擾動(dòng),記錄產(chǎn)品純度的響應(yīng)曲線,由圖13可知,控制結(jié)構(gòu)CS3較基礎(chǔ)控制結(jié)構(gòu)CS1,在抵抗±10%F擾動(dòng)方面有較大優(yōu)勢(shì),而在ZFEtOH發(fā)生±5%擾動(dòng)時(shí),系統(tǒng)可在短時(shí)間內(nèi)回歸穩(wěn)定狀態(tài),余差不大,且ZH2O的超調(diào)量有了較為明顯的降低:當(dāng)F增加10%時(shí),為0.015%,當(dāng)F減少10%時(shí)為0.012%。

    結(jié)果表明,與ADWC的控制規(guī)律相似,利用Qr/F代替Qr直接控制靈敏板溫度,有利降低體系進(jìn)料擾動(dòng)過程中的的超調(diào)量。主要是因?yàn)楫?dāng)F降低時(shí),此類控制結(jié)構(gòu)能迅速降低Qr,防止將大量重組分餾入塔頂;同理當(dāng)F增加時(shí),Qr也迅速做出反應(yīng),防止Qr過小引起輕組分從塔釜流失。

    圖12 帶有Qr/F的控制結(jié)構(gòu)CS3

    圖13Qr/F控制結(jié)構(gòu)CS3在±10%F擾動(dòng)下的動(dòng)態(tài)響應(yīng)曲線

    3 結(jié)論

    針對(duì)多元共沸物或近沸點(diǎn)混合物的分離,本文主要以IPA/水體系的ADWC和乙醇/水體系的EDWC為例進(jìn)行系統(tǒng)的動(dòng)態(tài)特性分析與控制研究,在滿足控制結(jié)構(gòu)簡(jiǎn)單、建設(shè)費(fèi)用低、工業(yè)實(shí)用性強(qiáng)的條件下獲得最優(yōu)的控制效果。本文選取不同的溫度靈敏板選擇方案進(jìn)行溫度靈敏板的確定,并建立了若干種控制結(jié)構(gòu),分析其性能,并通過測(cè)試對(duì)比選擇各過程最優(yōu)的控制結(jié)構(gòu)。

    (1)對(duì)于IPA/水體系的ADWC,調(diào)節(jié)兩塔Qr控制溫度靈敏板的溫度控制效果較好,且利用Qr/F控制靈敏板溫度的控制結(jié)構(gòu)CS3-1,相對(duì)于直接利用Qr作為操縱變量的控制結(jié)構(gòu)CS3可有效減少體系的超調(diào)量:當(dāng)F增加10%時(shí),ZIPA動(dòng)態(tài)響應(yīng)曲線的超調(diào)量由1.0%降為0.03%,當(dāng)F減小10%時(shí),ZH2O動(dòng)態(tài)響應(yīng)曲線的超調(diào)量由2.1%。降為0.08%

    (2)對(duì)于乙醇/水體系的EDWC,增加αv控制回路能夠減低體系調(diào)節(jié)余差,但效果不顯著。與ADWC控制結(jié)構(gòu)相似,利用Qr/F控制靈敏板溫度的控制結(jié)構(gòu)CS3相比于基礎(chǔ)控制結(jié)構(gòu)CS1,在抵抗進(jìn)料擾動(dòng)時(shí)能獲得更小的超調(diào)量:當(dāng)F增加10%時(shí)ZH2O超調(diào)量由0.1%降為0.015%,當(dāng)F減少10%時(shí)ZH2O超調(diào)量由0.08%降為0.012%。

    [1] 孫蘭義,昌興武,譚雅文,等. 熱耦合技術(shù)應(yīng)用于共沸精餾系統(tǒng)的研究[J]. 化工進(jìn)展,2010,29(12):2228-2233. SUN L Y,CHANG X W,TAN Y W,et al. Application of thermally coupled technology to azeotropic distillation system[J]. Chemical Industry and Engineering Progress,2010,29(12):2228-2233.

    [2] 李軍,王純正,馬占華,等. 隔壁塔用于苯、甲苯、二甲苯分離的控制[J]. 化工進(jìn)展,2013,32(4):757-762. LI J,WANG C Z,MA Z H,et al. Control of the DWC for separating the mixture of benzene,toluene,xylene[J]. Chemical Industry and Engineering Progress,2013,32(4):757-762.

    [3] 錢春健,葉青,朱國(guó)彪,等. DWC分離三組分混合物的研究[J]. 化工進(jìn)展,2007,26(8):1174-1177. QIAN C J,YE Q,ZHU G B,et al. Study of using DWC to separate three components mixture[J]. Chemical Industry and Engineering Progress,2007,26(8):1174-1177.

    [4] 馬晨皓,曾愛武. 隔壁塔流程模擬及節(jié)能效益的研究[J]. 化學(xué)工程,2013,41(3):1-5. MA C H,ZENG A W. Simulation of Dividing-Wall-Column process and its energy-saving benefits[J] Chemical Engineering(China),2013,41(3):1-5.

    [5] 李紅海,姜奕. 精餾塔設(shè)備的設(shè)計(jì)與節(jié)能研究進(jìn)展[J]. 化工進(jìn)展,2014,33(s1):14-18. LI H H,JIANG Y. Research progress of design and energy-saving of distillation column[J] Chemical Industry and Engineering Progress,2014,33(s1):14-18.

    [6] SUN L Y,WANG Q Y,LI L M,et al. Design and control of extractive dividing wall column for separating benzene/cyclohexane mixtures[J]. Industrial & Engineering Chemistry Research,2014,53(19):8120-8131.

    [7] 裘兆榮,葉青,李成益. 國(guó)內(nèi)外分隔壁精餾塔現(xiàn)狀與發(fā)展趨勢(shì)[J].江蘇工業(yè)學(xué)院學(xué)報(bào),2005,17(1):58-61. QIU Z R,YE Q,LI C Y. Status and development trends of dividing wall column at home and abroad[J]. Journal of Jiangsu Polytechnic University,2005,17(1):58-61.

    [8] 胡松,張冰劍,佘志鴻,等. 共沸蒸餾在化工生產(chǎn)中的應(yīng)用與研究進(jìn)展[J]. 化工進(jìn)展,2010,29(12):2207-2219. HU S,ZHANG B J,SHE Z H,et al. Advances in application of azeotropic distillation in chemical processes[J]. Chemical Industry and Engineering Progress,2010,29(12):2207-2219.

    [9] WU Y C,LEE H Y,HUANG H P,et al. Energy-saving dividing-wall column design and control for heterogeneous azeotropic distillation systems[J]. Industrial & Engineering Chemistry Research,2014,53(4):1537-1552.

    [10] BRIONES-RAMíREZ A,Gutiérrez-Antonio C. Dividing wall distillation columns for separation of azeotropic mixtures: feasibility procedure and rigorous optimization[J]. Computer Aided Chemical Engineering,2009,26:555-560.

    [11] KAIBEL B,JANSEN H,ZICH E,et al. Unfixed dividing wall technology for packed and tray distillation columns[J]. Distillation Absorption,2006,152:252-266.

    [12] KOLBE B,WENZEL S. Novel distillation concepts using one-shell columns[J]. Chemical Engineering and Processing:Process Intensification,2004,43(3):339-346.

    [13] DIEHL T,KOLBE B,GEHRKE H. Uhde Morphylane? extractive distillation where do we stand?[C]//ERTC Petrochemical Conference,Prague,Chech Republic,2005:1-12.

    [14] MIDORI S,ZHENG S N,YAMADEA K. Analysis of divided-wall column for extractive distillation[J]. Kagaku Kogaku Ronbunshu,2000,26(5):627-632.

    [15] LUYBEN W L. Plantwide control of an isopropyl alcohol dehydration process[J]. AIChE Journal,2006,52(6):2290-2296.

    [16] GMEHLING J. Azeotropic data[M]. Wiley-VCH,1994.

    [17] 王惠媛,許松林. 異丙醇-水分離技術(shù)進(jìn)展[J]. 上?;?,2005,30(6):20-24. WANG H L,XU S L. Separation technologies for isopropyl alcohol-water[J]. Shanghai Chemical Industry,2005,30(6):20-24.

    [18] LIU Y L,ZHAI J,LI L M,et al. Heat pump assisted reactive and azeotropic distillations in dividing wall columns[J]. Chemical Engineering and Processing: Process Intensification,2015,95:289-301.

    Control of azeotropic dividing wall column and extractive dividing wall column

    LIU Lixin1,CHEN Mengqi2,LIU Yuliang3,WANG Jianxin2,SUN Lanyi2
    (1College of Automobile and Transportation,Qingdao University of Technology,Qingdao 266520,Shandong,China;2College of Chemical Engineering,China University of Petroleum,Qingdao 266580,Shandong,China;3Department of Industrial Engineering & Innovation Sciences,Eindhoven University of Technology,Holland)

    Azeotropic and extractive dividing wall column were used for the separation of multiple azeotrope mixture or closing-boiling mixture. Steady-state simulation of the two processes was established by Aspen Plus. The different matching relationships between manipulate variables and control variables were determined through the sensitive analysis of steady state relative gain. Then,some two-point temperature control structures for azeotropic dividing wall column while several three-point (or four-point) temperature control structures for extractive dividing wall column were proposed in Aspen Dynamics. After adding dynamic disturbance,it chose the best one for each process: (1) best control structure for azeotropic dividing wall column:QMC/Fcontrol TMC,13,QRC/Fcontrol TRC,5,and (2) best control structure for extractive dividing wall column: RRM control TMC1,3,Qr/Fcontrol TMC1,12,RRR control TRC1,3,αv control TMC1,9. Finally,through the analysis of the similarity of the best control structures,it showed that the control structure with the ratio of reboiler duty and mixture feed rate (Qr/F) can achieve reasonable control performance,which reduced the overshoot of product purities and residual error of these two systerms.

    azeotrope;separation;azeotropic dividing wall column;extractive dividing wall column;control

    TQ028.1

    :A

    :1000–6613(2017)02–0756–10

    10.16085/j.issn.1000-6613.2017.02.049

    2016-07-11;修改稿日期:2016-09-11。

    山東省自然科學(xué)基金(ZR2013EEM024)及山東省高等學(xué)??萍加?jì)劃(J14LA08)項(xiàng)目。

    及聯(lián)系人:劉立新(1972—),女,碩士,副教授,研究方向?yàn)榈缆沸虏牧?、工藝、設(shè)備及控制。E-mail:llx297403478@163.com。

    猜你喜歡
    體系
    TODGA-TBP-OK體系對(duì)Sr、Ba、Eu的萃取/反萃行為研究
    “三個(gè)體系”助力交通安全百日攻堅(jiān)戰(zhàn)
    杭州(2020年23期)2021-01-11 00:54:42
    構(gòu)建體系,舉一反三
    探索自由貿(mào)易賬戶體系創(chuàng)新應(yīng)用
    常熟:構(gòu)建新型分級(jí)診療體系
    如何建立長(zhǎng)期有效的培訓(xùn)體系
    E-MA-GMA改善PC/PBT共混體系相容性的研究
    汽車零部件(2014年5期)2014-11-11 12:24:28
    “曲線運(yùn)動(dòng)”知識(shí)體系和方法指導(dǎo)
    加強(qiáng)立法工作 完善治理體系
    浙江人大(2014年1期)2014-03-20 16:19:53
    日本終身學(xué)習(xí)體系構(gòu)建的保障及其啟示
    五月伊人婷婷丁香| 男女边吃奶边做爰视频| 久久久色成人| 亚洲欧美精品自产自拍| 人人妻人人看人人澡| 久久久久久久久久久丰满| 麻豆精品久久久久久蜜桃| 成人无遮挡网站| 一级片'在线观看视频| 久久精品人妻少妇| 日日摸夜夜添夜夜添av毛片| 国产欧美日韩精品一区二区| 国产久久久一区二区三区| 欧美成人精品欧美一级黄| av在线播放精品| 又粗又硬又长又爽又黄的视频| 亚洲一区二区三区欧美精品 | 亚洲图色成人| 精品一区二区三区视频在线| 亚洲精品色激情综合| 午夜福利网站1000一区二区三区| 国产爱豆传媒在线观看| 午夜福利在线在线| 国产伦理片在线播放av一区| 不卡视频在线观看欧美| 人妻少妇偷人精品九色| xxx大片免费视频| 黄片wwwwww| 久久99精品国语久久久| 久久久久久久久久人人人人人人| 欧美97在线视频| 日本三级黄在线观看| 赤兔流量卡办理| 一级毛片黄色毛片免费观看视频| 最近最新中文字幕免费大全7| 免费观看无遮挡的男女| 又爽又黄无遮挡网站| 少妇的逼水好多| 人妻夜夜爽99麻豆av| 七月丁香在线播放| 91精品伊人久久大香线蕉| 超碰97精品在线观看| 97超碰精品成人国产| 中文天堂在线官网| 亚州av有码| 97热精品久久久久久| 蜜桃久久精品国产亚洲av| 日韩,欧美,国产一区二区三区| 欧美亚洲 丝袜 人妻 在线| 少妇的逼好多水| 午夜福利视频1000在线观看| 国产成人a区在线观看| 又黄又爽又刺激的免费视频.| 五月开心婷婷网| 天天躁日日操中文字幕| 午夜日本视频在线| 五月伊人婷婷丁香| 国产国拍精品亚洲av在线观看| 亚洲av电影在线观看一区二区三区 | 久久99蜜桃精品久久| 日韩精品有码人妻一区| 国产成年人精品一区二区| 狂野欧美激情性xxxx在线观看| 亚洲国产精品成人综合色| 国产老妇伦熟女老妇高清| 亚洲欧美清纯卡通| 午夜福利视频1000在线观看| 免费黄网站久久成人精品| 亚洲精华国产精华液的使用体验| 国产极品天堂在线| 欧美国产精品一级二级三级 | 欧美一区二区亚洲| 免费黄色在线免费观看| 可以在线观看毛片的网站| 亚洲av日韩在线播放| 国产精品国产三级国产av玫瑰| 亚洲va在线va天堂va国产| 国产91av在线免费观看| 内射极品少妇av片p| 天天一区二区日本电影三级| 国产色爽女视频免费观看| 亚洲aⅴ乱码一区二区在线播放| 夜夜看夜夜爽夜夜摸| 国产片特级美女逼逼视频| 亚洲怡红院男人天堂| 国产爽快片一区二区三区| 日韩大片免费观看网站| av专区在线播放| 18禁在线无遮挡免费观看视频| 夫妻午夜视频| 青青草视频在线视频观看| 亚洲国产精品成人综合色| 国产免费一区二区三区四区乱码| 大片免费播放器 马上看| 欧美国产精品一级二级三级 | 日韩视频在线欧美| 亚洲成人一二三区av| 一本一本综合久久| 亚洲精品国产色婷婷电影| 十八禁网站网址无遮挡 | 别揉我奶头 嗯啊视频| 免费电影在线观看免费观看| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美一区二区三区黑人 | 一区二区三区四区激情视频| 99精国产麻豆久久婷婷| 成人欧美大片| 人人妻人人爽人人添夜夜欢视频 | 日本色播在线视频| 国产精品秋霞免费鲁丝片| 亚洲综合色惰| 亚洲国产欧美人成| 色婷婷久久久亚洲欧美| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 久久人人爽av亚洲精品天堂 | 97超碰精品成人国产| 久久99精品国语久久久| 黄色怎么调成土黄色| 久久久久久久国产电影| 五月开心婷婷网| 精品一区二区免费观看| 精品久久久久久久末码| av线在线观看网站| 大陆偷拍与自拍| 免费电影在线观看免费观看| 久久精品夜色国产| 亚洲人成网站在线观看播放| 少妇人妻一区二区三区视频| 成人无遮挡网站| 国产 一区 欧美 日韩| 午夜福利在线观看免费完整高清在| xxx大片免费视频| 在线免费十八禁| 国产在线男女| 深夜a级毛片| 午夜亚洲福利在线播放| 亚洲成人一二三区av| 少妇高潮的动态图| 日日撸夜夜添| 亚洲成人久久爱视频| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区欧美精品 | 九九在线视频观看精品| 爱豆传媒免费全集在线观看| 秋霞在线观看毛片| 人妻一区二区av| 少妇丰满av| 日韩免费高清中文字幕av| 国产成人精品一,二区| 大话2 男鬼变身卡| 欧美国产精品一级二级三级 | 国产精品久久久久久精品电影小说 | 久久久欧美国产精品| av在线观看视频网站免费| 国产精品女同一区二区软件| 欧美日韩国产mv在线观看视频 | 免费观看性生交大片5| 精品国产三级普通话版| 网址你懂的国产日韩在线| 亚洲婷婷狠狠爱综合网| 亚洲av免费高清在线观看| 一级av片app| av在线播放精品| 寂寞人妻少妇视频99o| 国产 精品1| 日韩不卡一区二区三区视频在线| 成年版毛片免费区| 亚洲第一区二区三区不卡| 亚洲天堂av无毛| 国产一区二区亚洲精品在线观看| 在线免费十八禁| 久久久久久久国产电影| 日本三级黄在线观看| 久久久a久久爽久久v久久| 亚洲在久久综合| 尾随美女入室| 小蜜桃在线观看免费完整版高清| 国产大屁股一区二区在线视频| 啦啦啦啦在线视频资源| 免费播放大片免费观看视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 久久热精品热| 国产精品久久久久久精品古装| 久久久久久久久久久丰满| 97在线人人人人妻| 亚洲av电影在线观看一区二区三区 | 亚洲熟女精品中文字幕| 国产永久视频网站| 成人国产麻豆网| 激情 狠狠 欧美| 一本久久精品| 免费黄网站久久成人精品| 人体艺术视频欧美日本| 国产男人的电影天堂91| 麻豆乱淫一区二区| 黄色一级大片看看| 亚洲精品乱久久久久久| 少妇熟女欧美另类| 少妇高潮的动态图| 一边亲一边摸免费视频| 亚洲精品中文字幕在线视频 | 免费黄网站久久成人精品| 国内精品美女久久久久久| 偷拍熟女少妇极品色| 高清在线视频一区二区三区| 国产极品天堂在线| 久久午夜福利片| 免费人成在线观看视频色| 日本与韩国留学比较| 嫩草影院入口| 国产又色又爽无遮挡免| 国产精品国产三级专区第一集| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美日韩卡通动漫| 国产欧美亚洲国产| 九九在线视频观看精品| 日本wwww免费看| 欧美精品一区二区大全| 边亲边吃奶的免费视频| 婷婷色综合大香蕉| 国产精品秋霞免费鲁丝片| 免费观看无遮挡的男女| 亚洲成人一二三区av| 欧美日本视频| 亚洲精品日韩在线中文字幕| 色吧在线观看| 日韩 亚洲 欧美在线| 高清视频免费观看一区二区| 亚洲内射少妇av| 亚洲精品456在线播放app| 插逼视频在线观看| 美女xxoo啪啪120秒动态图| av线在线观看网站| 国产黄色免费在线视频| 欧美日韩国产mv在线观看视频 | 五月开心婷婷网| 精品久久久久久久久av| 免费黄网站久久成人精品| 久久久久久九九精品二区国产| 一二三四中文在线观看免费高清| 亚洲精华国产精华液的使用体验| 自拍欧美九色日韩亚洲蝌蚪91 | 乱系列少妇在线播放| 成人高潮视频无遮挡免费网站| 日韩三级伦理在线观看| 干丝袜人妻中文字幕| 全区人妻精品视频| av在线天堂中文字幕| 亚洲av在线观看美女高潮| 精品人妻偷拍中文字幕| 国产成人免费无遮挡视频| 在线观看一区二区三区激情| 日本午夜av视频| 毛片女人毛片| 免费av不卡在线播放| 97在线视频观看| 精品亚洲乱码少妇综合久久| 欧美日韩在线观看h| 亚洲精品国产av成人精品| 亚洲,欧美,日韩| 又粗又硬又长又爽又黄的视频| 欧美 日韩 精品 国产| 在线看a的网站| 91在线精品国自产拍蜜月| 欧美潮喷喷水| 亚洲一级一片aⅴ在线观看| 欧美丝袜亚洲另类| 99视频精品全部免费 在线| 国产高清国产精品国产三级 | 欧美日韩综合久久久久久| 欧美性感艳星| 国产欧美日韩精品一区二区| 亚洲精品色激情综合| 亚洲国产成人一精品久久久| 中文欧美无线码| 在线 av 中文字幕| www.av在线官网国产| 成人特级av手机在线观看| 又爽又黄a免费视频| 只有这里有精品99| 精品一区二区三区视频在线| av.在线天堂| 日本av手机在线免费观看| 可以在线观看毛片的网站| 亚洲国产av新网站| 青春草国产在线视频| 三级国产精品欧美在线观看| 久久久久九九精品影院| 欧美激情久久久久久爽电影| 国产精品精品国产色婷婷| 国产精品麻豆人妻色哟哟久久| 美女主播在线视频| 99热全是精品| 男女啪啪激烈高潮av片| 国产亚洲av嫩草精品影院| 国产精品一区二区三区四区免费观看| 成人漫画全彩无遮挡| 日韩一区二区视频免费看| 精品久久久久久电影网| 欧美丝袜亚洲另类| 交换朋友夫妻互换小说| 日韩欧美精品免费久久| 看免费成人av毛片| 精品一区在线观看国产| tube8黄色片| 美女主播在线视频| 久久久久性生活片| 国产高潮美女av| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久成人| 一级a做视频免费观看| 乱码一卡2卡4卡精品| 深夜a级毛片| 亚洲欧美日韩另类电影网站 | 亚洲成人精品中文字幕电影| 在线免费观看不下载黄p国产| 精品国产一区二区三区久久久樱花 | 观看免费一级毛片| 日本色播在线视频| 国产精品国产三级国产专区5o| 成年女人看的毛片在线观看| 日韩免费高清中文字幕av| 99久久九九国产精品国产免费| 欧美性感艳星| 日日摸夜夜添夜夜添av毛片| 久久久精品欧美日韩精品| 国产探花在线观看一区二区| 各种免费的搞黄视频| 狠狠精品人妻久久久久久综合| av在线观看视频网站免费| 下体分泌物呈黄色| 夜夜看夜夜爽夜夜摸| 大片免费播放器 马上看| 国产一区二区三区综合在线观看 | 亚洲欧美日韩另类电影网站 | 国产精品成人在线| 最近的中文字幕免费完整| 精品熟女少妇av免费看| 日韩欧美精品免费久久| 九草在线视频观看| 观看美女的网站| 国产在线男女| 成人特级av手机在线观看| 联通29元200g的流量卡| 亚洲精品色激情综合| 极品少妇高潮喷水抽搐| 亚洲怡红院男人天堂| 在线亚洲精品国产二区图片欧美 | 久久人人爽av亚洲精品天堂 | 免费黄频网站在线观看国产| 水蜜桃什么品种好| 国产伦精品一区二区三区四那| 国产免费福利视频在线观看| 成人特级av手机在线观看| 在线免费观看不下载黄p国产| 观看免费一级毛片| 日本-黄色视频高清免费观看| 国产伦精品一区二区三区视频9| 欧美少妇被猛烈插入视频| 最近的中文字幕免费完整| 国产在线一区二区三区精| 在线观看av片永久免费下载| 午夜亚洲福利在线播放| 国产综合懂色| 亚洲精品视频女| 久久97久久精品| 成人美女网站在线观看视频| av网站免费在线观看视频| 丝袜美腿在线中文| 我要看日韩黄色一级片| 日韩欧美精品v在线| 国产免费一区二区三区四区乱码| 18+在线观看网站| www.色视频.com| 激情 狠狠 欧美| 中文字幕亚洲精品专区| 一级毛片黄色毛片免费观看视频| 美女xxoo啪啪120秒动态图| 青春草视频在线免费观看| 免费看a级黄色片| 一区二区三区四区激情视频| 亚洲人成网站在线观看播放| 亚洲av福利一区| 嘟嘟电影网在线观看| 亚洲欧美一区二区三区黑人 | 欧美成人午夜免费资源| 精品人妻偷拍中文字幕| 久久精品综合一区二区三区| 男人狂女人下面高潮的视频| 亚洲av欧美aⅴ国产| 国产精品久久久久久av不卡| 国产乱人偷精品视频| 大话2 男鬼变身卡| 免费大片18禁| 成人美女网站在线观看视频| 一区二区三区精品91| 国产在视频线精品| 日韩成人伦理影院| 国产一区亚洲一区在线观看| 欧美 日韩 精品 国产| 中文字幕av成人在线电影| 最近手机中文字幕大全| 国产毛片a区久久久久| 男女边吃奶边做爰视频| av线在线观看网站| 高清午夜精品一区二区三区| 黄色怎么调成土黄色| 小蜜桃在线观看免费完整版高清| 大码成人一级视频| 欧美 日韩 精品 国产| freevideosex欧美| 国产精品熟女久久久久浪| 人妻一区二区av| 国产国拍精品亚洲av在线观看| 蜜桃久久精品国产亚洲av| 日日啪夜夜撸| 日本欧美国产在线视频| 少妇高潮的动态图| 中国国产av一级| 国产日韩欧美在线精品| 亚洲av成人精品一二三区| 色视频www国产| 91久久精品电影网| 国语对白做爰xxxⅹ性视频网站| 少妇高潮的动态图| 99九九线精品视频在线观看视频| 精品久久久噜噜| 免费播放大片免费观看视频在线观看| 国产视频内射| 男女啪啪激烈高潮av片| 91精品伊人久久大香线蕉| 人妻一区二区av| 高清av免费在线| 亚洲不卡免费看| 国产亚洲5aaaaa淫片| 男人和女人高潮做爰伦理| 搞女人的毛片| 久久久久久国产a免费观看| 蜜桃亚洲精品一区二区三区| 天堂中文最新版在线下载 | 国产 一区 欧美 日韩| av播播在线观看一区| 在线 av 中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 最近最新中文字幕免费大全7| 国产视频内射| 婷婷色麻豆天堂久久| 2021少妇久久久久久久久久久| 夫妻性生交免费视频一级片| 日韩强制内射视频| 国产大屁股一区二区在线视频| 日韩欧美 国产精品| 国精品久久久久久国模美| 色5月婷婷丁香| 欧美性猛交╳xxx乱大交人| 在线观看一区二区三区激情| 亚洲伊人久久精品综合| 国产美女午夜福利| 九九在线视频观看精品| 国产在线男女| 国产精品国产三级国产av玫瑰| 久久热精品热| 99久久精品热视频| 亚洲精品乱码久久久v下载方式| 一区二区av电影网| 亚洲成人一二三区av| 在线看a的网站| 禁无遮挡网站| 精品国产露脸久久av麻豆| 男男h啪啪无遮挡| 伦精品一区二区三区| 亚洲精品国产色婷婷电影| 秋霞伦理黄片| 高清日韩中文字幕在线| 国产v大片淫在线免费观看| 日产精品乱码卡一卡2卡三| 亚洲精品第二区| 91午夜精品亚洲一区二区三区| 国内精品宾馆在线| 精品国产一区二区三区久久久樱花 | 激情五月婷婷亚洲| av在线播放精品| 一区二区三区乱码不卡18| 蜜桃久久精品国产亚洲av| 国产一区二区三区av在线| 亚洲精品色激情综合| 在线观看人妻少妇| 日韩精品有码人妻一区| 建设人人有责人人尽责人人享有的 | 国产午夜福利久久久久久| 久久ye,这里只有精品| 国产精品久久久久久久久免| 大香蕉久久网| 日本av手机在线免费观看| 中文在线观看免费www的网站| 亚洲成人精品中文字幕电影| 久久精品综合一区二区三区| 别揉我奶头 嗯啊视频| 人体艺术视频欧美日本| 丰满人妻一区二区三区视频av| 色视频www国产| 人人妻人人澡人人爽人人夜夜| 嫩草影院新地址| 99热这里只有是精品50| 欧美高清成人免费视频www| 七月丁香在线播放| 麻豆精品久久久久久蜜桃| 国产成人精品久久久久久| 免费电影在线观看免费观看| av女优亚洲男人天堂| 国产精品久久久久久久电影| 午夜免费鲁丝| 亚洲在久久综合| 少妇的逼好多水| 嘟嘟电影网在线观看| 99视频精品全部免费 在线| 大片免费播放器 马上看| 午夜精品一区二区三区免费看| 22中文网久久字幕| 成年人午夜在线观看视频| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 极品教师在线视频| 纵有疾风起免费观看全集完整版| 午夜亚洲福利在线播放| 免费av观看视频| 亚洲国产av新网站| 极品少妇高潮喷水抽搐| 亚洲av中文av极速乱| 卡戴珊不雅视频在线播放| 欧美激情久久久久久爽电影| 少妇人妻久久综合中文| 国产精品一区www在线观看| 成人一区二区视频在线观看| 亚洲国产精品专区欧美| 国产 精品1| 日韩国内少妇激情av| 免费看a级黄色片| 最近手机中文字幕大全| 美女脱内裤让男人舔精品视频| 18禁在线播放成人免费| 看黄色毛片网站| 18禁裸乳无遮挡动漫免费视频 | 免费观看在线日韩| 男插女下体视频免费在线播放| 久久久欧美国产精品| 观看美女的网站| 国产伦精品一区二区三区四那| 国产精品人妻久久久影院| 麻豆精品久久久久久蜜桃| 国产极品天堂在线| 久久国内精品自在自线图片| 国产免费福利视频在线观看| 久久综合国产亚洲精品| 亚洲不卡免费看| 插逼视频在线观看| 久久99热这里只频精品6学生| 久久久久久久久久久丰满| 又大又黄又爽视频免费| 色婷婷久久久亚洲欧美| 日本黄大片高清| 91久久精品国产一区二区三区| 国产伦精品一区二区三区四那| 精品国产三级普通话版| 蜜桃久久精品国产亚洲av| 国内精品美女久久久久久| 高清欧美精品videossex| 国产精品一区二区三区四区免费观看| 免费大片18禁| 精品人妻一区二区三区麻豆| 久久久久久伊人网av| 性色av一级| 国产一区二区三区av在线| 午夜爱爱视频在线播放| 亚洲欧美一区二区三区黑人 | 欧美高清性xxxxhd video| 亚洲欧美成人综合另类久久久| 亚洲av中文字字幕乱码综合| 天天躁日日操中文字幕| 国产av码专区亚洲av| 国内精品宾馆在线| 少妇的逼水好多| 婷婷色综合www| 久久久亚洲精品成人影院| 最近手机中文字幕大全| av国产免费在线观看| 久久久精品免费免费高清| 国产男女内射视频| 久久久精品免费免费高清| 最后的刺客免费高清国语| 伊人久久精品亚洲午夜| 日韩 亚洲 欧美在线| 国产精品一及| 欧美一级a爱片免费观看看| 卡戴珊不雅视频在线播放| 亚洲真实伦在线观看| 国产精品蜜桃在线观看| 成人特级av手机在线观看| 国产成人精品一,二区| 国产精品女同一区二区软件| 纵有疾风起免费观看全集完整版| 男人和女人高潮做爰伦理| 亚洲av成人精品一二三区| 中国美白少妇内射xxxbb| av在线天堂中文字幕| 午夜福利网站1000一区二区三区| 亚洲图色成人| 97超视频在线观看视频| 日韩,欧美,国产一区二区三区| 最近的中文字幕免费完整| 亚洲熟女精品中文字幕| 午夜视频国产福利| 毛片女人毛片| 亚洲精品乱码久久久久久按摩|