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

    折流式外取熱器內(nèi)傳熱強(qiáng)化特性

    2021-01-27 01:51:22姚秀穎李建濤盧春喜
    關(guān)鍵詞:相區(qū)流板氣速

    劉 璐, 姚秀穎, 李建濤, 盧春喜

    (中國(guó)石油大學(xué)(北京) 重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102249)

    在石油加工、化工等工業(yè)過(guò)程中,為了保證系統(tǒng)溫度的可控性和穩(wěn)定性,通常配置合理的換熱器進(jìn)行調(diào)節(jié)。例如,在催化裂化裝置中,取熱器用于取走再生器內(nèi)多余的熱量,以維持反應(yīng)-再生系統(tǒng)的熱平衡,保證裝置的穩(wěn)定運(yùn)行。早期因取熱負(fù)荷要求不高,催化裂化裝置多采用內(nèi)取熱方式來(lái)取走再生器內(nèi)過(guò)剩的熱量;隨著重油和渣油催化裂化技術(shù)的發(fā)展,內(nèi)取熱方式已難以滿足工業(yè)生產(chǎn)需要,因此外取熱技術(shù)開(kāi)始快速發(fā)展,外取熱器作為獨(dú)立單元安裝在再生器外部,不僅不影響反應(yīng)-再生系統(tǒng)的正常運(yùn)行,而且操作彈性大,取熱負(fù)荷可調(diào)控[1];到20世紀(jì)90年代,已經(jīng)形成多種形式的外取熱器,按照流動(dòng)方式的不同分為4種,即上流式、下流式、氣控式和返混式外取熱器[2-3],其中下流式外取熱器因取熱負(fù)荷大、易于操作等優(yōu)點(diǎn)廣泛地應(yīng)用于工業(yè)生產(chǎn)中。

    但是,隨著中國(guó)催化裂化(FCC)原料不斷重質(zhì)化和劣質(zhì)化,生焦率增加,再生系統(tǒng)的產(chǎn)熱量不斷增加,同時(shí)外取熱器在實(shí)際工業(yè)應(yīng)用中也存在著取熱負(fù)荷不足和換熱管破裂等問(wèn)題,使其難以維持反應(yīng)-再生系統(tǒng)的熱量平衡[4-5]。而外取熱器取熱負(fù)荷不足會(huì)導(dǎo)致裝置劑/油比降低,影響輕質(zhì)油收率和產(chǎn)品分布,降低經(jīng)濟(jì)收益,同時(shí)也影響裝置的長(zhǎng)周期運(yùn)行[6-7]。因此,提高外取熱器取熱負(fù)荷至關(guān)重要。

    在外取熱器內(nèi),從換熱管外高溫催化劑到換熱管內(nèi)低溫流體的傳熱過(guò)程中,顆粒與換熱管表面的對(duì)流傳熱過(guò)程為控制熱阻[8]。外取熱器殼程內(nèi),顆粒與換熱管之間的傳熱過(guò)程不僅受顆粒床層流化性能的影響,而且與換熱管局部的流動(dòng)特性相關(guān)。很多學(xué)者通過(guò)改造換熱管的形狀從而強(qiáng)化傳熱,如:李瑞強(qiáng)等[9]提出在取熱管外側(cè)均布針狀取熱柱,增加換熱管外側(cè)的換熱面積,并通過(guò)釘頭的局部紊流提高傳熱效率;張榮克等[10]提出在換熱管表面焊接釘頭形式的構(gòu)件;唐逢一等[11-12]和顧月章等[13]對(duì)目前工業(yè)上常用的長(zhǎng)翅片管進(jìn)行了優(yōu)化,解決了縱向長(zhǎng)翅片所存在的應(yīng)力集中、易產(chǎn)生裂紋的問(wèn)題。但釘頭管和優(yōu)化后的翅片管對(duì)取熱性能的優(yōu)化效果較小。為了強(qiáng)化殼程內(nèi)傳熱性能,很多學(xué)者通過(guò)改造取熱器殼程從而強(qiáng)化外取熱器的殼程流化性能。王明東等[14-15]針對(duì)催化裂化工業(yè)裝置中出現(xiàn)的外取熱系統(tǒng)流化性能差、取熱器產(chǎn)汽量無(wú)法靈活調(diào)節(jié)等問(wèn)題,對(duì)外取熱器流化風(fēng)進(jìn)行改造,取熱負(fù)荷優(yōu)化明顯。程向峰[16-17]提出將外取熱器內(nèi)取熱管束分為上部和下部?jī)刹糠?,上部管束沿軸向從殼體上端插入,下部管束沿軸向從下部插入,并在下部管束下設(shè)第一流化介質(zhì)分布器,上、下管束間設(shè)第二流化介質(zhì)分布器,在2個(gè)分布器的作用下使顆粒流化。該技術(shù)將取熱管束分為上部和下部?jī)刹糠郑啾扔诠I(yè)上常用的取熱器,其使用的流化介質(zhì)少、傳熱效率高,但對(duì)原設(shè)備改動(dòng)較大。張永民等[18]提出了一種改進(jìn)型的強(qiáng)制內(nèi)混式流化床外取熱器,通過(guò)在外取熱器內(nèi)設(shè)置多個(gè)氣體分布器來(lái)強(qiáng)化傳熱過(guò)程。盧春喜等[19]開(kāi)發(fā)了一種環(huán)流結(jié)構(gòu)外取熱器,該取熱器具有相間接觸好、體積傳熱系數(shù)大等優(yōu)點(diǎn),但該技術(shù)受傳熱強(qiáng)化高度的限制,對(duì)外取熱器上部傳熱強(qiáng)化作用較小。

    為了有效地提高外取熱器的取熱負(fù)荷,筆者設(shè)計(jì)并開(kāi)發(fā)了一種新型折流式外取熱器,即在工業(yè)上常用的下流式外取熱器基礎(chǔ)上安裝若干折流板,通過(guò)強(qiáng)化顆粒徑向混合,提高整個(gè)軸向高度內(nèi)催化劑顆粒與換熱管束間的接觸頻率,進(jìn)而提高傳熱系數(shù)。為了驗(yàn)證其傳熱強(qiáng)化效果,采用一套與工業(yè)裝置具有相似傳熱機(jī)制的大型冷模實(shí)驗(yàn)裝置,分別對(duì)比考察常規(guī)下流式外取熱器(Down-flow catalyst cooler,DCC)和折流式外取熱器(Baffle catalyst cooler,BCC)內(nèi)換熱特性,并根據(jù)顆粒在密相區(qū)和稀相區(qū)的流動(dòng)特點(diǎn),深入分析折流式外取熱器傳熱強(qiáng)化的機(jī)制及流動(dòng)與傳熱特性的耦合關(guān)系,以期為折流式外取熱器的工業(yè)應(yīng)用提供參考。

    1 實(shí)驗(yàn)部分

    1.1 實(shí)驗(yàn)裝置及操作參數(shù)

    采用的外取熱器換熱示意圖如圖1所示。下流式外取熱器和折流式外取熱器具有相同的尺寸,其殼體的內(nèi)徑為0.48 m、高為5.5 m。外取熱器的上部設(shè)有催化劑入口,入口管路與取熱器外壁呈45°,下部設(shè)有垂直向下的催化劑出口。為了保證與工業(yè)外取熱器具有相同的水力學(xué)半徑,8根外徑為 0.076 m 的換熱管沿直徑0.334 m的同心圓周向均勻布置,每根換熱管高為4.4 m,管外均布8個(gè)翅片,每個(gè)翅片分為4段高度為0.45 m的短翅片,以減小長(zhǎng)翅片造成的熱應(yīng)力不均、翅片上易產(chǎn)生小裂紋的問(wèn)題[10]。在換熱管下部,設(shè)有與工業(yè)裝置一致的環(huán)管式氣體分布器。相比于下流式外取熱器,折流式外取熱器在其上部區(qū)域的換熱管內(nèi)部區(qū)域焊接6塊傾斜向下45°的折流板,相鄰折流板間距為0.235 m,共設(shè)計(jì)了2種形式的折流板。其中,由上往下第1塊折流板,其傾斜方向與入口管傾斜方向一致,以增加入口處折流板對(duì)催化劑流股的整流作用,減小外取熱器入口處催化劑對(duì)換熱管的沖刷,同時(shí)折流板上設(shè)8個(gè)直徑為40 mm的開(kāi)孔,以破碎氣泡。由上往下第2~6塊折流板,每個(gè)折流板上設(shè)4個(gè)直徑為40 mm的開(kāi)孔,該形式折流板與外取熱器內(nèi)壁距離較遠(yuǎn),以避免折流板上出現(xiàn)死區(qū),并對(duì)床層氣泡起到一定的破碎作用。

    圖1 外取熱器換熱示意圖Fig.1 Diagram of external catalyst coolers(a) Down-flow catalyst cooler——DCC;(b) Baffle catalyst cooler——BCC

    下流式和折流式外取熱器均采用與工業(yè)裝置中相似的傳熱機(jī)制,實(shí)驗(yàn)過(guò)程中,常溫催化劑顆粒以恒定流量由外取熱器上部的入口管流入,經(jīng)過(guò)氣體流化,與高溫?fù)Q熱管換熱后,從下部出口管流出,入口管和出口管上分別設(shè)有DN100 mm的D71X-10型對(duì)夾式蝶閥,用于控制催化劑循環(huán)量和藏料量。換熱管由高溫水加熱,經(jīng)恒溫水槽加熱的水由換熱管下部流入,上部流出,為換熱管提供熱源。

    實(shí)驗(yàn)流化氣體為常溫常壓空氣,流化顆粒為催化裂化平衡催化劑,其平均粒徑為70.3 μm、顆粒密度為1500 kg/m3、松散堆積密度為856.9 kg/m3。自上部入口流入的催化劑顆粒具有穩(wěn)定的流量,其顆粒質(zhì)量流量(Gs)控制在0.32~44.96 kg/(m2·s)范圍內(nèi),外取熱器藏料量保持為247.80 kg,表觀氣速(Ug)為0.1~0.5 m/s,為滿足實(shí)驗(yàn)測(cè)量需要,換熱管內(nèi)熱水流量恒定。

    1.2 測(cè)量?jī)x器及測(cè)點(diǎn)布置

    為分析外取熱器內(nèi)流動(dòng)和傳熱特性,分別測(cè)量了床層的壓力、床層溫度、換熱管壁溫及其內(nèi)部水溫。其中,壓力采用差壓變送器進(jìn)行測(cè)量,采樣頻率為100 Hz,采樣時(shí)間為60 s;床層溫度、換熱管壁溫及其內(nèi)部水溫均采用鎧裝圓柱形熱敏電阻進(jìn)行測(cè)量,采集時(shí)間1800 s,采樣時(shí)間間隔0.2 s。

    為準(zhǔn)確對(duì)比下流式和折流式外取熱器內(nèi)流動(dòng)和傳熱特性,采用相同的測(cè)點(diǎn)布置。其中,外取熱器內(nèi)氣體分布器出口處所在的軸向高度H=0 m,向上為正,軸向測(cè)點(diǎn)布置如圖2(a)所示,其具體軸向高度如表1所示。外取熱器內(nèi)有壓力p測(cè)點(diǎn)10個(gè),由于靜床高度為2 m,為確定密相床層高度,在H為2~2.5 m區(qū)域內(nèi)壓力測(cè)點(diǎn)分布較密;床層溫度(Tb)測(cè)點(diǎn)4個(gè),如圖2(a)所示,每根換熱管上設(shè)置4段短翅片,4個(gè)Tb測(cè)點(diǎn)沿軸向布置于4段短翅片對(duì)應(yīng)的中間位置;壁面溫度(Tw)測(cè)點(diǎn)14個(gè),其中,12個(gè)位于換熱管外部的翅片間,2個(gè)位于換熱管上部的光管上;水溫度(Twater)測(cè)點(diǎn)8個(gè),其中6個(gè)位于短翅片間的光管處,2個(gè)位于換熱管進(jìn)出口處。p、Tb、Tw和Twater的周向測(cè)點(diǎn)如圖2(b)所示,其中,Tw和Twater測(cè)點(diǎn)均布置在通有熱水的換熱管上。

    圖2 外取熱器內(nèi)實(shí)驗(yàn)測(cè)點(diǎn)布置圖Fig.2 Diagram of experimental measurement pointsin the external catalyst cooler(a) Axial distribution of measurement points;(b) Circumferential distribution of measurement points

    表1 外取熱器內(nèi)實(shí)驗(yàn)測(cè)點(diǎn)軸向高度(H)Table 1 Height (H) of axial measurement points inexternal catalyst cooler

    1.3 數(shù)據(jù)處理方法

    主要考察折流式外取熱器相對(duì)于下流式外取熱器的傳熱強(qiáng)化特性,并結(jié)合床層流動(dòng)特性探討其傳熱強(qiáng)化的原因,分析2種外取熱器內(nèi)傳熱與流動(dòng)特性的耦合關(guān)系。根據(jù)顆粒在不同區(qū)域的流動(dòng)特點(diǎn),將床層分為密相區(qū)和稀相區(qū),床層分界面由床層密度(ρb)確定,ρb由式(1)計(jì)算。床層流動(dòng)特性通過(guò)壓力數(shù)據(jù)進(jìn)行分析,壓力脈動(dòng)的標(biāo)準(zhǔn)偏差(Sd)反映床層脈動(dòng)特性,Sd越大,床層壓力脈動(dòng)越大,氣泡聚并、破碎頻率越高,Sd可由式(2)計(jì)算。

    (1)

    (2)

    在流化床內(nèi),隨著表觀氣速Ug的增加,流化床內(nèi)Sd先增加后減小。當(dāng)Sd達(dá)到最大值時(shí),對(duì)應(yīng)Ug即為起始湍流速度(Uc)[20]。通過(guò)該方法確定Uc,以H=1.330 m,Gs≈18 kg/(m2·s)為例,該位置床層處于穩(wěn)定的密相床層操作。圖3為表觀氣速對(duì)壓力脈動(dòng)標(biāo)準(zhǔn)偏差的影響??梢钥闯觯?種外取熱器內(nèi)Sd均隨著Ug的增加,先增大后減小,下流式和折流式外取熱器內(nèi)的Ug分別為 0.4 m/s 和 0.3 m/s 時(shí)Sd達(dá)到最大值,因此,下流式和折流式外取熱器內(nèi)的Uc分別為0.4 m/s和0.3 m/s。陽(yáng)永榮[21]得到在直徑為800 mm的流化床內(nèi),采用平均粒徑為66.7 μm、顆粒密度為1480 kg/m3的FCC催化劑起始湍流速度為0.35 m/s,與本實(shí)驗(yàn)所用物料參數(shù)相近,起始湍流速度基本一致。

    圖3 表觀氣速(Ug)對(duì)壓力脈動(dòng)標(biāo)準(zhǔn)偏差(Sd)的影響Fig.3 Effects of superficial gas velocity (Ug) on standarddeviation of pressure fluctuation (Sd)

    床層換熱特性通過(guò)傳熱系數(shù)進(jìn)行分析,密相區(qū)和稀相區(qū)的傳熱系數(shù)計(jì)算式分別如式(3)和(4)所示,密相區(qū)和稀相區(qū)的換熱面積計(jì)算式分別如式(5)和(6)所示。

    (3)

    (4)

    Adense=πDHd+2nHdL

    (5)

    Adilute=πD(H-Hd)+2n(H-Hd)L

    (6)

    2 結(jié)果與討論

    2.1 床層分界面的確定

    在實(shí)驗(yàn)操作條件下,外取熱器內(nèi)床層沿軸向高度分為稀相區(qū)和密相區(qū),以Gs≈18 kg/(m2·s)、Ug為0.2和0.4 m/s為例,分析下流式和折流式外取熱器內(nèi)ρb隨H的變化規(guī)律,以確定稀、密相區(qū)的床層分界面。

    圖4為下流式和折流式外取熱器內(nèi)床層密度的軸向分布。如圖4(a)所示,隨H的增加,下流式外取熱器內(nèi)ρb先緩慢增大,后急劇減小。這是因?yàn)樵谙铝魇酵馊崞髦?,?dāng)H<1.1 m時(shí),床層處于氣體分布器影響區(qū),由于分布器射流和射流結(jié)束后產(chǎn)生大量小氣泡的作用,使得ρb隨H的增加逐漸增加;當(dāng)H為1.1~1.8 m時(shí),床層發(fā)展為穩(wěn)定的密相區(qū),ρb波動(dòng)較?。划?dāng)H>1.8 m時(shí),ρb急劇減小,床層由密相區(qū)變?yōu)橄∠鄥^(qū)。如圖4(b)所示,隨H的增加,折流式外取熱器內(nèi)ρb先緩慢增加后降低,再緩慢增加后急劇減小。相比于下流式外取熱器,相同操作條件下,在折流式外取熱器內(nèi),當(dāng)H>1.8 m時(shí),ρb小幅波動(dòng)后急劇減?。磺耶?dāng)H>2.2 m 時(shí),ρb急劇減小。結(jié)合折流板安裝位置分析發(fā)現(xiàn),從上往下第4塊折流板安裝高度為1.90~2.09 m,可能是由于該位置折流板上開(kāi)孔對(duì)氣泡有破碎作用所致,各氣速下折流式外取熱器內(nèi)H=2.2 m附近ρb均較大。

    圖4 外取熱器內(nèi)床層密度的軸向分布Fig.4 Axial profiles of bed density in the external catalyst cooler(a) DCC; (b) BCC

    由上述結(jié)果可知,在外取熱器的不同軸向高度,床層密度分布差異較大,在床層下部為密度較高的密相區(qū),上部則為密度較低的稀相區(qū)。為了獲得不同區(qū)域的傳熱特性,分析其傳熱機(jī)理,首先確定床層稀相區(qū)和密相區(qū)的床層分界面高度Hd。張永民等[22]提出采用相鄰測(cè)點(diǎn)間的壓力脈動(dòng)比Sd,p(i,i+1)來(lái)確定床層分界面,認(rèn)為當(dāng)某兩測(cè)量點(diǎn)之間的Sd,p(i,i+1)最大時(shí),料面便位于第i和第i+1 壓力測(cè)點(diǎn)間,得到Hd的計(jì)算式如式(7)所示。

    (7)

    通過(guò)上述實(shí)驗(yàn)發(fā)現(xiàn),料面所在的第i和第i+1壓力測(cè)點(diǎn)間的壓力差不僅受密相區(qū)密度的影響,當(dāng)氣速較高時(shí),還受到稀相區(qū)密度的影響,因此在張永民等[22]方法的基礎(chǔ)上,進(jìn)一步考慮了稀相區(qū)密度對(duì)第i和第i+1壓力測(cè)點(diǎn)間壓力差的影響,得到改進(jìn)后的Hd計(jì)算式如式(8)所示。

    (8)

    通過(guò)改進(jìn)方法確定的Hd隨著表觀氣速的變化如圖5所示??梢钥闯?,隨著Ug的增加,下流式和折流式外取熱器內(nèi)Hd均先增加后減小。進(jìn)一步分析發(fā)現(xiàn),2種外取熱器內(nèi),當(dāng)表觀氣速達(dá)到Uc時(shí),Hd值最大。

    圖5 表觀氣速(Ug)對(duì)床層分界面高度(Hd)的影響Fig.5 Effects of superficial gas velocity (Ug) on theheight (Hd) of dense phase area

    因此,根據(jù)床層壓力脈動(dòng)特性確定的Hd,將外取熱器內(nèi)床層準(zhǔn)確分為密相區(qū)和稀相區(qū),Hd以下為密相區(qū)換熱,Hd以上為稀相區(qū)換熱。

    2.2 密相區(qū)換熱及流動(dòng)特性

    已有研究[23-24]表明,顆粒與床層之間的換熱主要受顆粒團(tuán)在換熱表面的接觸時(shí)間分率和接觸頻率的影響。筆者通過(guò)換熱管外的傳熱系數(shù)h對(duì)換熱性能進(jìn)行表征。h越大,表明換熱性能越好。密相區(qū)內(nèi)顆粒是連續(xù)相,氣體是分散相,床層換熱主要為顆粒對(duì)流換熱。通過(guò)分析密相區(qū)換熱特性發(fā)現(xiàn),下流式和折流式外取熱器內(nèi)密相區(qū)傳熱系數(shù)hdense均受Ug和Gs的影響。

    圖6為表觀氣速對(duì)hdense的影響。由圖6可見(jiàn):對(duì)于下流式外取熱器,當(dāng)Gs一定時(shí),隨著Ug的增加,hdense逐漸增大;當(dāng)Ug>0.4 m/s時(shí),hdense增幅逐漸減小。這是由于下流式外取熱器內(nèi)Uc=0.4 m/s,當(dāng)Ug

    圖6 表觀氣速(Ug)對(duì)密相區(qū)傳熱系數(shù)(hdense)的影響Fig.6 Effects of superficial gas velocity (Ug) ondense heat transfer coefficient (hdense)

    在H為1.54 m和1.86 m處,床層既不受取熱器下部的分布器影響區(qū)的影響,也不受密相區(qū)和稀相區(qū)交界處的彈濺區(qū)的影響,流動(dòng)狀態(tài)穩(wěn)定,可有效代表密相區(qū)流動(dòng)特性[22]。圖7為表觀氣速對(duì)床層密度的影響。由圖7可知:當(dāng)UgUc時(shí),密相區(qū)內(nèi)床層處于湍動(dòng)流化狀態(tài),床層內(nèi)氣泡尺寸變小但氣泡數(shù)增加,顆粒團(tuán)與換熱管間的接觸頻率較高,但隨Ug的增加,密相區(qū)內(nèi)ρb持續(xù)降低,抑制了接觸時(shí)間分率的進(jìn)一步升高,受氣泡運(yùn)動(dòng)和床層密度的共同影響,hdense的增幅變小。

    對(duì)于折流式外取熱器(如圖6所示),當(dāng)Gs一定時(shí),hdense隨著Ug的增加持續(xù)增大。這是因?yàn)橄啾扔谙铝魇酵馊崞?,在折流式外取熱器?nèi),當(dāng)Ug>Uc時(shí),ρb并未降低,而是幾乎不變(圖7),折流板的存在強(qiáng)化了氣泡的破碎,抑制了氣泡的進(jìn)一步增大,使得折流式外取熱器能在較高氣速下保證較高的床層密度,而顆粒團(tuán)在換熱表面的接觸頻率較高,因此hdense持續(xù)增加。

    圖7 表觀氣速(Ug)對(duì)床層密度(ρb)的影響Fig.7 Effect of superficial gas velocity (Ug)on bed density (ρb)

    相比于表觀氣速,下流式和折流式外取熱器內(nèi)密相區(qū)傳熱系數(shù)受Gs的影響較小,如圖8所示。由圖8可以看出,隨著Gs的增加,2種外取熱器內(nèi)hdense的變化趨勢(shì)基本一致。當(dāng)Ug=0.1 m/s時(shí),hdense均隨著Gs的增大而減??;當(dāng)Ug>0.1 m/s時(shí),hdense均隨著Gs的增大而增大。這是因?yàn)楫?dāng)Ug=0.1 m/s時(shí),ρb較大,床層內(nèi)氣泡數(shù)量較少且運(yùn)動(dòng)不夠劇烈,冷顆粒團(tuán)與熱換熱管之間接觸頻率較低,充分換熱所需時(shí)間較長(zhǎng)。顆粒質(zhì)量流量越小,催化劑在取熱器內(nèi)停留時(shí)間越長(zhǎng),因此換熱效果越好。當(dāng)Ug>0.1 m/s時(shí),床層內(nèi)氣泡數(shù)較多,氣泡運(yùn)動(dòng)較劇烈,冷顆粒團(tuán)與熱換熱管之間接觸頻率較高,充分換熱所需時(shí)間較短。隨著顆粒質(zhì)量流量的增大,催化劑循環(huán)速率加快、顆粒在換熱表面的停留時(shí)間變短,換熱后的熱催化劑很快被流入的新鮮冷催化劑所替代,有利于換熱。

    為確定折流式外取熱器的強(qiáng)化換熱性能,通過(guò)比較2種外取熱器內(nèi)hdense(見(jiàn)圖8)發(fā)現(xiàn):實(shí)驗(yàn)操作條件下,下流式外取熱器內(nèi)hdense為60.2~626.8 W/(m2·℃),折流式外取熱器內(nèi)hdense為155.9~938.2 W/(m2·℃),hdense提高了125.5~345.2 W/(m2·℃)。相比于下流式外取熱器,折流式外取熱器在Ug=0.1 m/s時(shí)的hdense相對(duì)增長(zhǎng)量為159.0%~267.9%;在Ug為 0.2~0.5 m/s時(shí),折流式外取熱器的hdense相對(duì)增長(zhǎng)量為26.9%~79.9%;折流式外取熱器換熱性能優(yōu)化明顯。而且,低氣速下折流式外取熱器內(nèi)密相區(qū)強(qiáng)化換熱性能更強(qiáng)。這是因?yàn)榈蜌馑傧耯dense較小,

    圖8 顆粒質(zhì)量流量(Gs)對(duì)密相區(qū)換熱系數(shù)(hdense)的影響Fig.8 Effect of particle mass flux (Gs) on dense heat transfer coefficient (hdense)(a) DCC; (b) BCC

    使得hdense的相對(duì)增長(zhǎng)量較大。分析發(fā)現(xiàn),折流式外取熱器主要通過(guò)以下幾個(gè)方面強(qiáng)化密相區(qū)換熱:

    (1)強(qiáng)化徑向流動(dòng)。下流式外取熱器中,冷催化劑從上部流入,換熱后的熱催化劑直接從下部流出。折流式外取熱器中,床層內(nèi)催化劑流股在軸向高度2.0 m附近,形成了如圖9所示的2種交替的流動(dòng)狀態(tài),實(shí)驗(yàn)觀察到,短時(shí)間內(nèi)床層呈現(xiàn)向下流動(dòng)(見(jiàn)圖9(a)),大部分時(shí)間內(nèi)則為沖擊流股(見(jiàn)圖9(b)中紅圈附近)。結(jié)合折流板的安裝位置發(fā)現(xiàn),該沖擊區(qū)域位置位于由下往上第4塊折流板的斜下方。即增設(shè)折流板后,增加了床層中催化劑流股的徑向流動(dòng),使催化劑沿折流板形成“折流”式流動(dòng),提高了床層的湍動(dòng)程度,更有利于換熱。

    (2)加劇軸向混合。增設(shè)折流板后,以Gs≈18 kg/(m2·s)為例,各氣速下密相區(qū)床溫的軸向溫差ΔTb均減小,如圖10所示,說(shuō)明增設(shè)折流板后,密相床層的軸向混合程度增加,換熱更充分。這是因?yàn)橄铝魇酵馊崞鲀?nèi)催化劑呈無(wú)序、無(wú)規(guī)律的流動(dòng);增設(shè)折流板后,催化劑流股在折流板的作用下沿折流板呈“折流”式流動(dòng),且折流板上的開(kāi)孔破碎了向上運(yùn)動(dòng)的大氣泡,氣泡與“折流”式流動(dòng)的催化劑多次接觸,使得折流式取熱器內(nèi)催化劑有序流動(dòng),因此強(qiáng)化了密相區(qū)床層的軸向混合。

    (3)加劇氣泡破碎,增加接觸頻率。以H=1.330 m、Gs≈18 kg/(m2·s)為例(圖3),通過(guò)比較下流式和折流式外取熱器床層內(nèi)Sd發(fā)現(xiàn),與下流式外取熱器相比,在折流式外取熱器內(nèi),當(dāng)Ug>0.1 m/s時(shí),密相區(qū)內(nèi)Sd明顯減小。這是因?yàn)檩^高氣速下床層內(nèi)氣泡數(shù)較多,折流式的增設(shè),抑制了氣泡尺寸的增大,同時(shí)折流板上的開(kāi)孔也加劇了氣泡的破碎,使得折流式外取熱器內(nèi)氣泡尺寸減小,氣泡數(shù)量增加,床層壓力脈動(dòng)較小;此時(shí)密相區(qū)床層的氣泡尺寸小且數(shù)量多,不僅增加了顆粒團(tuán)在換熱表面的接觸頻率,而且減小了顆粒的停留時(shí)間,使得換熱性能提高。當(dāng)Ug=0.1 m/s時(shí),密相區(qū)內(nèi)Sd減小的幅度較小,但hdense增幅較大(如圖6所示)。這是因?yàn)楫?dāng)Ug=0.1 m/s時(shí),床層內(nèi)氣泡數(shù)較少,折流板上的開(kāi)孔對(duì)氣泡的破碎作用不明顯,因此造成Sd減小的幅度較小。但是,當(dāng)Ug=0.1 m/s時(shí),增設(shè)折流板后,ΔTb減小(如圖10 所示),密相床層的軸向混合程度增加,利于換熱。此外,當(dāng)Ug=0.1 m/s時(shí),由下往上運(yùn)動(dòng)的氣泡數(shù)較少,催化劑流股主要呈由上往下的“折流”式流動(dòng),顆粒流通當(dāng)量路徑長(zhǎng)度Lp增加,相同換熱時(shí)間內(nèi)催化劑流股流經(jīng)換熱管壁的速率增加,hdense增幅較大。

    圖9 外取熱器內(nèi)床層流動(dòng)直觀圖(H=2.0 m)Fig.9 Pictures of fluidized bed in the external catalyst coolers (H=2.0 m)(a) Condition 1; (b) Condition 2

    圖10 表觀氣速(Ug)對(duì)密相區(qū)床層溫度軸向溫度差(ΔTb)的影響Fig.10 Effect of superficial gas velocity (Ug) onaxial temperature difference (ΔTb) in the dense region

    2.3 稀相區(qū)換熱及流動(dòng)特性

    當(dāng)Ug≥0.4 m/s時(shí),下流式和折流式外取熱器內(nèi)床層均由鼓泡床轉(zhuǎn)變?yōu)橥膭?dòng)床,床層稀相夾帶量增大,稀相區(qū)換熱不可忽略。在稀相區(qū),氣體為連續(xù)相,顆粒為分散相,雖然氣體對(duì)流換熱已不可忽視,但由于顆粒的體積比熱容為氣體的103倍,因此,在較高氣速下,稀相區(qū)傳熱系數(shù)同樣主要受顆粒流動(dòng)特性的影響。分別考察Ug≥0.4 m/s時(shí),下流式和折流式2種外取熱器內(nèi)稀相區(qū)傳熱系數(shù)hdilute隨Gs和Ug的變化規(guī)律。

    圖11為下流式和折流式外取熱器內(nèi)Gs對(duì)稀相區(qū)傳熱系數(shù)hdilute的影響。可以看出,在下流式外取熱器內(nèi),hdilute隨著Gs的增加而逐漸增大。這是因?yàn)殡S著Gs的增加,從外取熱器入口管流入的新鮮催化劑的量增加,單位時(shí)間內(nèi)與換熱管有效接觸的新鮮催化劑的量增加,催化劑在換熱表面的更新頻率加快,因此hdilute增大。在折流式外取熱器內(nèi),hdilute隨著Gs的增加先增大、后基本保持不變。當(dāng)Gs<40 kg/(m2·s)時(shí),hdilute的變化趨勢(shì)與下流式外取熱器基本相同,即隨著Gs的增加,催化劑在換熱表面的更新頻率加快,hdilute增大;當(dāng)Gs>40 kg/(m2·s)時(shí),雖然催化劑在換熱表面的更新頻率進(jìn)一步加快,但是Gs過(guò)高,催化劑流股在床層內(nèi)的循環(huán)速率加快,停留時(shí)間減短,抑制了hdilute的進(jìn)一步增加。

    進(jìn)一步分析發(fā)現(xiàn),顆粒質(zhì)量流量一定時(shí),2種外取熱器內(nèi)hdilute均隨著Ug的增加而增大。這是因?yàn)殡S著Ug的增加,床層稀相夾帶量增加,2種外取熱器內(nèi)稀相區(qū)ρb均增大(如圖4所示),因此顆粒在換熱管外壁的體積分?jǐn)?shù)增大,則傳熱系數(shù)增大。

    圖11 顆粒質(zhì)量流量(Gs)對(duì)稀相區(qū)傳熱系數(shù)(hdilute)的影響Fig.11 Effect of particle mass flux (Gs) on dilute heat transfer coefficient (hdilute)(a) DCC; (b) BCC

    對(duì)比2種外取熱器內(nèi)稀相區(qū)換熱性能(見(jiàn)圖11)發(fā)現(xiàn),在實(shí)驗(yàn)操作條件下,下流式外取熱器內(nèi)hdilute為26.3~152.5 W/(m2·℃),折流式外取熱器內(nèi)hdilute為50.2~214.9 W/(m2·℃),增設(shè)折流板后,hdilute增加了23.9~112.0 W/(m2·℃),hdilute的相對(duì)增長(zhǎng)量為31.6%~152.5%,稀相區(qū)換熱性能得到優(yōu)化。hdilute的相對(duì)增長(zhǎng)量較大是由于下流式外取熱器稀相區(qū)的對(duì)流傳熱系數(shù)較小,相對(duì)增長(zhǎng)量的基數(shù)較小。分析發(fā)現(xiàn):在下流式外取熱器內(nèi),催化劑流股經(jīng)入口管流入后直接以拋物線形式向下流動(dòng),顆粒流僅在單一方向與換熱管進(jìn)行熱量交換;在折流式外取熱器內(nèi),由于多層折流板的存在,催化劑流股經(jīng)入口管流入后沿折流板以“折流”形式向下流動(dòng),即顆粒在向下流動(dòng)的同時(shí),存在徑向流動(dòng),顆粒與換熱管多次進(jìn)行熱量交換,同時(shí)折流板開(kāi)孔使催化劑在取熱器截面上分布更加均勻,增加了由入口管路流入取熱器的新鮮顆粒與換熱管的接觸分率,進(jìn)一步強(qiáng)化了稀相傳熱。此外,稀相區(qū)催化劑的“折流”流動(dòng)增長(zhǎng)了顆粒在稀相區(qū)的過(guò)流路徑,使得稀相換熱更充分,hdilute增大。

    2.4 密相區(qū)與稀相區(qū)強(qiáng)化傳熱特性的比較

    為了詳細(xì)比較下流式和折流式外取熱器在密相區(qū)和稀相區(qū)的傳熱特性,以下流式外取熱器內(nèi)兩區(qū)域的傳熱系數(shù)為相對(duì)強(qiáng)化量,分析折流式外取熱器內(nèi)平均傳熱系數(shù)的絕對(duì)強(qiáng)化度(Δh)和相對(duì)強(qiáng)化比率(δh),其計(jì)算式分別如式(9)和式(10)所示,其結(jié)果如表2所示。

    (9)

    (10)

    由表2可以看出,各操作條件下,Δhdense為186.0~299.9 W/(m2·℃),且Δhdense隨著表觀氣速的增加,先緩慢減小后迅速增大;δh,dense為43.5%~196.0%,其隨著Ug的增加先迅速減小后基本不變。這是因?yàn)樵谙铝魇酵馊崞鲀?nèi),傳熱系數(shù)隨著Ug的增加而增大,使得較高氣速下相對(duì)強(qiáng)化量較大,因此δh,dense有所降低。當(dāng)表觀氣速為 0.4 m/s 和0.5 m/s時(shí),稀相換熱不可忽略,計(jì)算得到Δhdilute分別為47.8和89.3 W/(m2·℃),且Δhdilute隨著Ug的增加而增大。這是因?yàn)殡S著Ug的增加,稀相夾帶量增大,顆粒濃度增大,顆粒團(tuán)與換熱管之間的接觸分率得到提高;δh, dilute分別為140.0%和79.1%,同樣,由于下流式外取熱器傳熱系數(shù)在較高氣速下較大,較高氣速下相對(duì)強(qiáng)化量較大,所以高氣速下δh, dilute較低,低氣速下δh, dilute較高。因此,無(wú)論在密相區(qū)還是在稀相區(qū),折流式外取熱器的傳熱效果均強(qiáng)于下流式外取熱器。

    表2 密相區(qū)和稀相區(qū)內(nèi)傳熱系數(shù)的絕對(duì)強(qiáng)化度(Δh)和相對(duì)強(qiáng)化比率(δh)Table 2 Δh and δh of heat transfer coefficients indense region and dilute region

    比較表觀氣速大于0.4 m/s時(shí)在密相區(qū)和稀相區(qū)的折流式外取熱器的傳熱強(qiáng)化效果發(fā)現(xiàn),雖然δh,dilute為79.1%~140.0%,明顯大于δh,dense(43.5%~51.0%),但是Δhdilute為47.8和89.3 W/(m2·℃),而Δhdense為228.9和299.9 W/(m2·℃),密相區(qū)內(nèi)折流式外取熱器內(nèi)平均傳熱系數(shù)的絕對(duì)強(qiáng)化度更大,因此,折流板的存在,對(duì)密相區(qū)的強(qiáng)化傳熱效果明顯優(yōu)于稀相區(qū)。這主要是因?yàn)槊芟鄥^(qū)床層密度遠(yuǎn)大于稀相區(qū)床層密度,且密相區(qū)顆粒團(tuán)與換熱管的接觸頻率更高;同時(shí),由于受氣體分布器的影響,密相區(qū)下部氣泡的形成、運(yùn)動(dòng)、聚并和破碎也加快了顆粒團(tuán)與換熱管的接觸頻率;而稀相區(qū)床層密度較低,顆粒占據(jù)換熱表面的體積分?jǐn)?shù)較低,因此傳熱系數(shù)較低。

    2.5 密相區(qū)傳熱系數(shù)關(guān)聯(lián)式的建立

    由圖8和圖11可知,當(dāng)表觀氣速大于0.4 m/s時(shí),下流式外取熱器內(nèi)密相區(qū)傳熱系數(shù)在465.6~626.8 W/(m2·℃)范圍內(nèi),稀相區(qū)傳熱系數(shù)在26.3~152.5 W/(m2·℃)范圍內(nèi);而在折流式外取熱器內(nèi)密相區(qū)傳熱系數(shù)在591.0~938.2 W/(m2·℃)范圍內(nèi),稀相區(qū)傳熱系數(shù)在50.2~214.9 W/(m2·℃)范圍內(nèi)。即相同操作條件下,密相區(qū)傳熱系數(shù)是稀相區(qū)的4.4~15.3倍,密相區(qū)為主要換熱區(qū)域。因此,建立了下流式和折流式外取熱器密相區(qū)的無(wú)因次換熱關(guān)聯(lián)式。通過(guò)分析催化劑顆粒的流通路徑和影響流化床內(nèi)床層與壁面之間換熱的因素發(fā)現(xiàn),影響床層與壁面之間換熱的因素主要分為3類,即物性參數(shù)、操作參數(shù)和幾何參數(shù)。物性參數(shù)由顆粒和流化氣體的物性確定,包括氣體密度ρg、氣體導(dǎo)熱系數(shù)kg、氣體比熱容cpg、氣體黏度μg、顆粒直徑dp、顆粒密度ρp和顆粒比熱容cpp;操作參數(shù)包括操作溫度T、操作壓力p、表觀氣速Ug、顆粒質(zhì)量流量Gs、床層密度ρb和稀-密相床層分界面的高度Hd,其中,Hd由筆者所提出的式(8)進(jìn)行計(jì)算;幾何參數(shù)包括換熱管外徑D和床層直徑Db。由于取熱器內(nèi)布置了8根換熱管,因此筆者采用顆粒流通截面當(dāng)量直徑De來(lái)代替Db,De由催化劑流股的流通截面積和流體與壁面的接觸長(zhǎng)度確定。同時(shí)由于折流板的存在,顆粒的流動(dòng)路徑不再能由床層分界面的高度Hd進(jìn)行表述,因此筆者提出顆粒流通當(dāng)量路徑長(zhǎng)度Lp來(lái)代替Hd,在此針對(duì)不同結(jié)構(gòu)的取熱器,Lp采用不同的計(jì)算方法,具體由顆粒流經(jīng)取熱器殼程的路徑確定。

    將各參數(shù)進(jìn)行無(wú)量綱化處理為:

    進(jìn)一步,得到密相區(qū)換熱關(guān)聯(lián)式的形式為:

    (11)

    由于本研究中氣、固兩相物性恒定,因此式(11)可簡(jiǎn)化為式(12)。

    (12)

    通過(guò)擬合,分別得到下流式和折流式外取熱器密相區(qū)傳熱系數(shù)預(yù)測(cè)關(guān)聯(lián)式的擬合系數(shù)如表3所示。圖12為下流式和折流式外取熱器內(nèi)密相區(qū)傳熱系數(shù)的耦合關(guān)聯(lián)式預(yù)測(cè)值與實(shí)驗(yàn)值之間的相對(duì)誤差??梢钥闯觯?jì)算值與實(shí)驗(yàn)值之間的相對(duì)誤差均小于±20%,處于工程允許誤差范圍內(nèi)。

    表3 外取熱器內(nèi)經(jīng)驗(yàn)關(guān)聯(lián)式擬合系數(shù)Table 3 Fitting coefficients of empirical correlationin the external catalyst cooler

    圖12 密相區(qū)傳熱系數(shù)的耦合關(guān)聯(lián)式預(yù)測(cè)值與實(shí)驗(yàn)值之間的相對(duì)誤差Fig.12 Relative error between predicted value from empirical correlation and experimental data in dense phase region(a) DCC; (b) BCC

    3 結(jié) 論

    通過(guò)大型冷模實(shí)驗(yàn)裝置,分別考察了下流式和折流式外取熱器在密相區(qū)和稀相區(qū)的傳熱和流動(dòng)性能。

    (1)在密相區(qū),相比于下流式外取熱器,折流式外取熱器內(nèi)折流板的增設(shè),強(qiáng)化了顆粒的徑向流動(dòng),加劇了密相床層的軸向混合和氣泡破碎,使得折流式外取熱器的密相傳熱系數(shù)明顯提高。

    (2)在稀相區(qū),折流式外取熱器內(nèi)折流板的增設(shè),使催化劑流股經(jīng)入口管流入后,沿折流板以“折流”形式流動(dòng),增長(zhǎng)了顆粒在稀相區(qū)的過(guò)流路徑;折流板開(kāi)孔使稀相區(qū)催化劑分布更加均勻。因此,折流式外取熱器內(nèi)稀相傳熱系數(shù)大于下流式外取熱器。

    (3)在密相區(qū)和稀相區(qū),折流式外取熱器的傳熱效果均明顯優(yōu)于下流式外取熱器,其中密相區(qū)傳熱系數(shù)的相對(duì)增長(zhǎng)量為26.9%~267.9%,稀相區(qū)傳熱系數(shù)的相對(duì)增長(zhǎng)量為31.6%~152.5%。

    (4)為了準(zhǔn)確區(qū)分密相和稀相換熱區(qū),提出了確定床層分界面高度Hd的新方法。在下流式和折流式外取熱器內(nèi)Hd均隨著Ug的增大而先增大后減小,其峰值點(diǎn)出現(xiàn)在Ug=Uc情況下。并引入催化劑流通當(dāng)量直徑和當(dāng)量路徑長(zhǎng)度,建立了下流式和折流式外取熱器的密相區(qū)無(wú)因次換熱關(guān)聯(lián)式,其計(jì)算值與實(shí)驗(yàn)值之間的相對(duì)誤差均小于±20%。

    符號(hào)說(shuō)明:

    Adense——密相區(qū)換熱面積,m2;

    Adilute——稀相區(qū)換熱面積,m2;

    c——顆粒比熱容與氣體比熱容之比;

    cpg——?dú)怏w比熱容,J/(kg·℃);

    cpp——顆粒比熱容,J/(kg·℃);

    cwater——水比熱容,J/(kg·℃);

    dp——顆粒直徑,m;

    D——換熱管外徑,m;

    Db——床層直徑,m;

    De——顆粒流通截面當(dāng)量直徑,m;

    Gs——顆粒質(zhì)量流量,kg/(m2·s);

    g——重力加速度,m/s2;

    H——軸向高度,m;

    Hd——床層分界面高度,m;

    Hi——第i個(gè)壓力測(cè)點(diǎn)對(duì)應(yīng)的軸向高度,m;

    h——傳熱系數(shù),W/(m2·℃);

    hdense——密相區(qū)傳熱系數(shù),W/(m2·℃);

    hdilute——稀相區(qū)傳熱系數(shù),W/(m2·℃);

    K0~K6——擬合系數(shù);

    kg——?dú)怏w導(dǎo)熱系數(shù),W/(m·℃);

    L——翅片高度,m;

    Lp——顆粒流通當(dāng)量路徑長(zhǎng)度,m;

    mwater——水的質(zhì)量流量,kg/s;

    N——數(shù)據(jù)采樣總數(shù);

    Nu——努塞爾準(zhǔn)數(shù);

    n——周向翅片數(shù)目;

    Pr——普朗特準(zhǔn)數(shù);

    p——操作壓力,kPa;

    pi——壓力瞬時(shí)值,kPa;

    Re——雷諾準(zhǔn)數(shù);

    Sd——壓力脈動(dòng)的標(biāo)準(zhǔn)偏差,kPa;

    Sd,p——壓力脈動(dòng)比;

    T——操作溫度,℃;

    Tb——床層溫度,℃;

    Tb,dense——密相區(qū)床層溫度的軸向均值,℃;

    Tb,dilute——稀相區(qū)床層溫度的軸向均值,℃;

    Tw——壁面溫度,℃;

    Tw,dense——密相區(qū)壁面溫度的軸向均值,℃;

    Tw,dilute——稀相區(qū)壁面溫度的軸向均值,℃;

    Twater——水溫度,℃;

    Twater,d——Hd對(duì)應(yīng)的上、下2個(gè)水溫度測(cè)點(diǎn)的均值,℃;

    Twater,in——進(jìn)口處測(cè)點(diǎn)水溫度,℃;

    Twater,out——出口處測(cè)點(diǎn)水溫度,℃;

    Uc——起始湍流速度,m/s;

    Ug——表觀氣速,m/s;

    ΔH——軸向高度差,m;

    Δh——傳熱系數(shù)的絕對(duì)強(qiáng)化度,W/(m2·℃);

    Δp——壓力差,kPa;

    ΔTb——床層溫度的軸向溫度差,℃;

    δh——傳熱系數(shù)的相對(duì)強(qiáng)化比率,%;

    μg——?dú)怏w黏度,Pa·s;

    ρb——床層密度,kg/m3;

    ρg——?dú)怏w密度,kg/m3;

    ρp——顆粒密度,kg/m3。

    猜你喜歡
    相區(qū)流板氣速
    結(jié)構(gòu)參數(shù)對(duì)交錯(cuò)百葉折流板管殼式換熱器性能影響的研究
    新型錐形孔折流板管殼式換熱器數(shù)值模擬研究
    五陽(yáng)熱電廠4×75t/h CFB鍋爐爐型分析及其改造
    同煤科技(2019年1期)2019-05-16 01:46:22
    濕法煙氣脫硫吸收塔阻力特性實(shí)驗(yàn)研究
    二組分溫度組成相圖單雙相區(qū)的經(jīng)驗(yàn)規(guī)則
    化工管理(2017年12期)2017-05-12 08:38:28
    新型折板除霧器的流場(chǎng)和壓降數(shù)值模擬
    螺旋折流板換熱器在高溫高壓環(huán)境中的應(yīng)用研究
    D120 mm流化床冷模實(shí)驗(yàn)研究
    化工科技(2014年5期)2014-06-09 05:17:22
    DSG太陽(yáng)能槽式集熱器的熱性能研究
    三氯氫硅合成爐運(yùn)行總結(jié)
    性欧美人与动物交配| 亚洲精品久久成人aⅴ小说| 九九热线精品视视频播放| 日本 欧美在线| 非洲黑人性xxxx精品又粗又长| 日本一本二区三区精品| 成人特级黄色片久久久久久久| 无遮挡黄片免费观看| 中文字幕人妻丝袜一区二区| 给我免费播放毛片高清在线观看| 精品国产乱子伦一区二区三区| 国产精品国产高清国产av| 国产男靠女视频免费网站| 国产午夜精品久久久久久| 午夜免费激情av| 日韩欧美在线二视频| 亚洲欧美日韩东京热| 亚洲乱码一区二区免费版| videosex国产| x7x7x7水蜜桃| 久久久久免费精品人妻一区二区| 久久国产乱子伦精品免费另类| 视频区欧美日本亚洲| av福利片在线| 91老司机精品| 一区二区三区高清视频在线| 90打野战视频偷拍视频| 女警被强在线播放| 亚洲av片天天在线观看| 国产男靠女视频免费网站| 午夜福利免费观看在线| e午夜精品久久久久久久| 丁香六月欧美| 日本五十路高清| 久久精品综合一区二区三区| x7x7x7水蜜桃| 在线十欧美十亚洲十日本专区| 国产亚洲精品第一综合不卡| 最近视频中文字幕2019在线8| 99久久国产精品久久久| 91九色精品人成在线观看| 嫩草影视91久久| 国产一区二区在线av高清观看| 成年人黄色毛片网站| 91大片在线观看| 黑人巨大精品欧美一区二区mp4| 久久精品夜夜夜夜夜久久蜜豆 | 国产97色在线日韩免费| 欧美绝顶高潮抽搐喷水| 美女大奶头视频| 香蕉丝袜av| 国产日本99.免费观看| 亚洲中文日韩欧美视频| 亚洲成人精品中文字幕电影| 日本a在线网址| 成年女人毛片免费观看观看9| 女人爽到高潮嗷嗷叫在线视频| 婷婷精品国产亚洲av在线| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 久久久久久久久中文| 嫩草影视91久久| 国产精品爽爽va在线观看网站| x7x7x7水蜜桃| 又黄又爽又免费观看的视频| 亚洲美女视频黄频| 久久性视频一级片| 国内精品久久久久精免费| 动漫黄色视频在线观看| 十八禁网站免费在线| 最近在线观看免费完整版| 怎么达到女性高潮| 国产精品久久久久久亚洲av鲁大| 嫩草影院精品99| 亚洲人成伊人成综合网2020| 91九色精品人成在线观看| 亚洲精品一区av在线观看| 天天添夜夜摸| 日韩欧美在线二视频| 在线a可以看的网站| 1024香蕉在线观看| 真人一进一出gif抽搐免费| 午夜精品一区二区三区免费看| 欧美成人午夜精品| 国产亚洲欧美在线一区二区| 99久久综合精品五月天人人| 国产精品免费视频内射| 国产欧美日韩一区二区三| 又大又爽又粗| 看免费av毛片| 少妇人妻一区二区三区视频| 日韩欧美免费精品| 99在线视频只有这里精品首页| 亚洲人成网站高清观看| 激情在线观看视频在线高清| 天天躁狠狠躁夜夜躁狠狠躁| 俺也久久电影网| 亚洲av成人精品一区久久| 色哟哟哟哟哟哟| 日韩三级视频一区二区三区| 国产精品美女特级片免费视频播放器 | 久久久精品大字幕| 久久久久亚洲av毛片大全| 丝袜美腿诱惑在线| 国产av一区在线观看免费| 国产亚洲精品久久久久久毛片| 日韩欧美三级三区| 狂野欧美白嫩少妇大欣赏| 亚洲国产欧美网| 亚洲自拍偷在线| 亚洲七黄色美女视频| 欧美中文综合在线视频| 在线a可以看的网站| 亚洲男人天堂网一区| 久9热在线精品视频| 国产激情偷乱视频一区二区| 中文字幕人成人乱码亚洲影| 757午夜福利合集在线观看| 亚洲精品久久成人aⅴ小说| 亚洲一码二码三码区别大吗| 99热这里只有是精品50| 亚洲欧美日韩无卡精品| 久久性视频一级片| 一级毛片高清免费大全| 一a级毛片在线观看| 一边摸一边抽搐一进一小说| 国产av麻豆久久久久久久| 国产精品久久久久久亚洲av鲁大| 亚洲精品中文字幕一二三四区| 法律面前人人平等表现在哪些方面| 日本成人三级电影网站| 床上黄色一级片| 免费在线观看亚洲国产| 国产精品99久久99久久久不卡| 亚洲精品色激情综合| 成人av在线播放网站| 99久久99久久久精品蜜桃| 无遮挡黄片免费观看| 日本 欧美在线| 精品国产乱子伦一区二区三区| 亚洲av成人不卡在线观看播放网| 老鸭窝网址在线观看| 久久热在线av| 日韩精品青青久久久久久| 欧美日韩乱码在线| 亚洲欧美日韩高清在线视频| 成人三级黄色视频| a在线观看视频网站| 国产探花在线观看一区二区| 国产又色又爽无遮挡免费看| 亚洲成人久久爱视频| 搡老熟女国产l中国老女人| 一区福利在线观看| 1024手机看黄色片| 少妇熟女aⅴ在线视频| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| 色综合婷婷激情| 色噜噜av男人的天堂激情| 十八禁人妻一区二区| 夜夜躁狠狠躁天天躁| 亚洲av成人不卡在线观看播放网| 日韩高清综合在线| 一本大道久久a久久精品| 久久精品国产亚洲av香蕉五月| 午夜精品在线福利| 欧美日韩瑟瑟在线播放| 亚洲熟妇中文字幕五十中出| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区三区四区免费观看 | 精品福利观看| 中文字幕熟女人妻在线| 人妻久久中文字幕网| 欧美成人午夜精品| 亚洲精品中文字幕一二三四区| 一本精品99久久精品77| 午夜影院日韩av| 免费看十八禁软件| 国产亚洲精品一区二区www| 国产主播在线观看一区二区| 三级毛片av免费| 国内精品一区二区在线观看| 女生性感内裤真人,穿戴方法视频| 女人被狂操c到高潮| 日本成人三级电影网站| 男女之事视频高清在线观看| 99精品在免费线老司机午夜| 久久国产乱子伦精品免费另类| 日韩高清综合在线| 欧美乱色亚洲激情| 亚洲精品美女久久av网站| 狂野欧美激情性xxxx| 国产视频一区二区在线看| 久久香蕉国产精品| 伦理电影免费视频| 久久国产乱子伦精品免费另类| 久久久久九九精品影院| 99国产精品一区二区三区| 国产私拍福利视频在线观看| 狂野欧美激情性xxxx| 在线观看www视频免费| 黄色视频,在线免费观看| 成人18禁在线播放| 亚洲国产精品sss在线观看| 嫁个100分男人电影在线观看| 男男h啪啪无遮挡| 亚洲精品久久成人aⅴ小说| 可以在线观看的亚洲视频| 成人高潮视频无遮挡免费网站| 五月玫瑰六月丁香| 国产精品亚洲av一区麻豆| 全区人妻精品视频| 性色av乱码一区二区三区2| 观看免费一级毛片| avwww免费| 别揉我奶头~嗯~啊~动态视频| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 亚洲男人的天堂狠狠| 国产高清视频在线播放一区| cao死你这个sao货| 欧美日本视频| 亚洲精品美女久久久久99蜜臀| 老司机午夜福利在线观看视频| 久久婷婷成人综合色麻豆| 性色av乱码一区二区三区2| 久久久精品大字幕| 香蕉av资源在线| 国产日本99.免费观看| 国产精品久久久久久久电影 | 一二三四社区在线视频社区8| 欧美丝袜亚洲另类 | 午夜免费观看网址| 校园春色视频在线观看| 午夜影院日韩av| 亚洲成av人片免费观看| 变态另类成人亚洲欧美熟女| 欧美日韩乱码在线| 亚洲午夜理论影院| 国产一区二区在线av高清观看| 亚洲国产精品久久男人天堂| 国产精品久久久av美女十八| 亚洲欧美激情综合另类| 夜夜看夜夜爽夜夜摸| 亚洲性夜色夜夜综合| 成人18禁高潮啪啪吃奶动态图| 久久国产乱子伦精品免费另类| 欧美成狂野欧美在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲男人的天堂狠狠| 亚洲av成人不卡在线观看播放网| 欧美人与性动交α欧美精品济南到| 777久久人妻少妇嫩草av网站| 露出奶头的视频| 国产人伦9x9x在线观看| 一级a爱片免费观看的视频| 一级作爱视频免费观看| 亚洲成人免费电影在线观看| 亚洲国产精品成人综合色| 中文字幕av在线有码专区| 国产私拍福利视频在线观看| svipshipincom国产片| av天堂在线播放| 精品欧美国产一区二区三| 国产乱人伦免费视频| 舔av片在线| 两个人看的免费小视频| 亚洲欧美日韩高清专用| 又紧又爽又黄一区二区| 免费人成视频x8x8入口观看| 国产av麻豆久久久久久久| 草草在线视频免费看| 夜夜夜夜夜久久久久| 三级国产精品欧美在线观看 | 精品久久久久久久久久久久久| 无遮挡黄片免费观看| 亚洲人成电影免费在线| 国产一级毛片七仙女欲春2| 啪啪无遮挡十八禁网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产v大片淫在线免费观看| 老熟妇仑乱视频hdxx| 久久精品91无色码中文字幕| 欧美日韩精品网址| 日韩av在线大香蕉| 天天一区二区日本电影三级| 国产黄色小视频在线观看| 看黄色毛片网站| 88av欧美| 99热这里只有是精品50| 在线免费观看的www视频| 校园春色视频在线观看| 美女午夜性视频免费| 亚洲欧美激情综合另类| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 麻豆成人av在线观看| 性色av乱码一区二区三区2| 国产成人精品久久二区二区免费| 久久精品综合一区二区三区| 国产1区2区3区精品| 99精品欧美一区二区三区四区| 熟女少妇亚洲综合色aaa.| 午夜影院日韩av| 一进一出抽搐动态| 一进一出抽搐gif免费好疼| 十八禁人妻一区二区| 老熟妇乱子伦视频在线观看| 国产高清有码在线观看视频 | 19禁男女啪啪无遮挡网站| 久久热在线av| 美女黄网站色视频| 日韩欧美国产在线观看| 国产免费av片在线观看野外av| 特大巨黑吊av在线直播| 精华霜和精华液先用哪个| 久久这里只有精品19| 两个人的视频大全免费| 午夜精品久久久久久毛片777| 一区二区三区激情视频| 青草久久国产| 身体一侧抽搐| 午夜精品久久久久久毛片777| 一区二区三区激情视频| 两性夫妻黄色片| 中文亚洲av片在线观看爽| a级毛片在线看网站| 久久久久国产精品人妻aⅴ院| 中文字幕熟女人妻在线| 中文资源天堂在线| 亚洲性夜色夜夜综合| 亚洲中文日韩欧美视频| or卡值多少钱| 欧美高清成人免费视频www| 亚洲精品美女久久av网站| 久久久久久久久免费视频了| 欧美黄色片欧美黄色片| 精品久久久久久久久久免费视频| 十八禁网站免费在线| 国产又黄又爽又无遮挡在线| 在线免费观看的www视频| 天堂影院成人在线观看| 国产欧美日韩一区二区三| 国产精品 欧美亚洲| av在线播放免费不卡| 可以免费在线观看a视频的电影网站| 午夜免费成人在线视频| 99久久精品国产亚洲精品| 国产三级中文精品| 国产精品美女特级片免费视频播放器 | 99国产极品粉嫩在线观看| 国产成人aa在线观看| 久久久久久人人人人人| 在线观看免费午夜福利视频| 97碰自拍视频| 国内精品久久久久久久电影| 久久久国产欧美日韩av| 国产v大片淫在线免费观看| 国产成人系列免费观看| 欧美日韩亚洲国产一区二区在线观看| 婷婷六月久久综合丁香| 国产片内射在线| 久久99热这里只有精品18| 欧美另类亚洲清纯唯美| 久久中文看片网| 国产午夜精品久久久久久| 精品久久久久久久毛片微露脸| 桃红色精品国产亚洲av| 在线十欧美十亚洲十日本专区| or卡值多少钱| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 日韩三级视频一区二区三区| 久久久精品国产亚洲av高清涩受| 成人国产综合亚洲| 99精品欧美一区二区三区四区| av片东京热男人的天堂| 黄色成人免费大全| 午夜福利免费观看在线| 亚洲天堂国产精品一区在线| 亚洲成人中文字幕在线播放| 亚洲精华国产精华精| 99热这里只有是精品50| 精品国内亚洲2022精品成人| 国产亚洲精品第一综合不卡| 青草久久国产| 亚洲va日本ⅴa欧美va伊人久久| 国产精品免费一区二区三区在线| 国产亚洲欧美98| 午夜两性在线视频| 午夜亚洲福利在线播放| 正在播放国产对白刺激| 无遮挡黄片免费观看| 亚洲成av人片免费观看| 好看av亚洲va欧美ⅴa在| 久久久久久久久久黄片| 99久久综合精品五月天人人| 午夜视频精品福利| 亚洲国产精品成人综合色| 男男h啪啪无遮挡| 黄片小视频在线播放| 在线看三级毛片| 午夜福利在线在线| 亚洲欧美精品综合久久99| 国产一区在线观看成人免费| 久久伊人香网站| 少妇被粗大的猛进出69影院| 亚洲人成网站高清观看| 色综合站精品国产| 搡老熟女国产l中国老女人| 在线观看美女被高潮喷水网站 | 99国产极品粉嫩在线观看| 欧美丝袜亚洲另类 | 午夜亚洲福利在线播放| 777久久人妻少妇嫩草av网站| 久久精品成人免费网站| 免费在线观看影片大全网站| 黄色片一级片一级黄色片| 亚洲真实伦在线观看| 欧美丝袜亚洲另类 | 99国产精品一区二区蜜桃av| av福利片在线观看| 国产亚洲av嫩草精品影院| 99久久久亚洲精品蜜臀av| 国产成人一区二区三区免费视频网站| 欧美av亚洲av综合av国产av| svipshipincom国产片| 亚洲,欧美精品.| 亚洲 欧美一区二区三区| 露出奶头的视频| 国产免费男女视频| 老鸭窝网址在线观看| 日本黄色视频三级网站网址| 制服丝袜大香蕉在线| 好看av亚洲va欧美ⅴa在| 两个人视频免费观看高清| 亚洲专区国产一区二区| 香蕉国产在线看| 国产精品综合久久久久久久免费| 热99re8久久精品国产| 无遮挡黄片免费观看| 午夜成年电影在线免费观看| 亚洲av美国av| 亚洲av中文字字幕乱码综合| 性色av乱码一区二区三区2| 欧美日韩精品网址| 真人一进一出gif抽搐免费| www.自偷自拍.com| 精品久久蜜臀av无| 丰满人妻一区二区三区视频av | 日韩大码丰满熟妇| 999精品在线视频| 很黄的视频免费| 亚洲国产精品sss在线观看| 午夜成年电影在线免费观看| 亚洲熟妇中文字幕五十中出| 亚洲午夜理论影院| 男女视频在线观看网站免费 | 最新在线观看一区二区三区| avwww免费| 国产在线精品亚洲第一网站| 哪里可以看免费的av片| 国产伦一二天堂av在线观看| 久久精品成人免费网站| 亚洲专区中文字幕在线| 色综合欧美亚洲国产小说| 国产私拍福利视频在线观看| 免费一级毛片在线播放高清视频| 90打野战视频偷拍视频| 成人国产综合亚洲| 后天国语完整版免费观看| 97碰自拍视频| www.精华液| 亚洲精品色激情综合| 日韩欧美在线乱码| 老鸭窝网址在线观看| 精品国产超薄肉色丝袜足j| 美女 人体艺术 gogo| av国产免费在线观看| 久久人妻av系列| 九九热线精品视视频播放| 欧美日韩瑟瑟在线播放| 久久午夜亚洲精品久久| av天堂在线播放| www.精华液| 老汉色av国产亚洲站长工具| 在线观看舔阴道视频| 一级片免费观看大全| 色哟哟哟哟哟哟| 亚洲免费av在线视频| 国产主播在线观看一区二区| 人妻丰满熟妇av一区二区三区| 国产精品av久久久久免费| 777久久人妻少妇嫩草av网站| 麻豆国产av国片精品| 免费看a级黄色片| 欧美国产日韩亚洲一区| 国产亚洲欧美在线一区二区| 成人国产综合亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成人精品中文字幕电影| 伦理电影免费视频| 露出奶头的视频| 成人国产综合亚洲| 亚洲色图 男人天堂 中文字幕| 免费看a级黄色片| 欧洲精品卡2卡3卡4卡5卡区| 在线免费观看的www视频| 黑人操中国人逼视频| av视频在线观看入口| 亚洲一区二区三区色噜噜| 国产视频内射| 免费av毛片视频| 国产在线精品亚洲第一网站| 黄色视频,在线免费观看| 国产黄片美女视频| 久99久视频精品免费| av中文乱码字幕在线| 日本三级黄在线观看| 久久久水蜜桃国产精品网| 少妇熟女aⅴ在线视频| 一本一本综合久久| 亚洲国产精品sss在线观看| 两个人看的免费小视频| 亚洲九九香蕉| 国产亚洲av高清不卡| 在线看三级毛片| 亚洲激情在线av| 看黄色毛片网站| 最近最新免费中文字幕在线| 亚洲全国av大片| 欧美极品一区二区三区四区| 国产高清视频在线观看网站| 国产午夜精品论理片| 日本一二三区视频观看| 精品人妻1区二区| 级片在线观看| 国产精品九九99| 老司机在亚洲福利影院| 亚洲成人久久爱视频| 日本 av在线| 99国产精品99久久久久| 老鸭窝网址在线观看| 好男人电影高清在线观看| 熟女少妇亚洲综合色aaa.| 国产三级黄色录像| 欧美精品亚洲一区二区| 国产成人av教育| 久久精品aⅴ一区二区三区四区| 亚洲av中文字字幕乱码综合| 日本 欧美在线| 国产精品影院久久| 级片在线观看| 欧美日韩亚洲综合一区二区三区_| 中文亚洲av片在线观看爽| 欧美黄色淫秽网站| 中文字幕人妻丝袜一区二区| 亚洲va日本ⅴa欧美va伊人久久| 欧美色视频一区免费| 国产不卡一卡二| 麻豆一二三区av精品| 一级片免费观看大全| 久久久久亚洲av毛片大全| 欧美大码av| 我要搜黄色片| 久久婷婷成人综合色麻豆| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区激情短视频| av免费在线观看网站| 中文字幕精品亚洲无线码一区| 亚洲激情在线av| 无限看片的www在线观看| 国产欧美日韩一区二区精品| avwww免费| 色综合婷婷激情| 男插女下体视频免费在线播放| 美女大奶头视频| 亚洲精品av麻豆狂野| 中国美女看黄片| 午夜福利18| 国产99久久九九免费精品| 国产精品久久久久久人妻精品电影| 嫩草影视91久久| 国产又黄又爽又无遮挡在线| 搡老妇女老女人老熟妇| 国内少妇人妻偷人精品xxx网站 | 草草在线视频免费看| 久久精品国产综合久久久| 黄片大片在线免费观看| 精华霜和精华液先用哪个| 麻豆一二三区av精品| 午夜免费激情av| 99re在线观看精品视频| 久久久久久久久中文| 91国产中文字幕| 精华霜和精华液先用哪个| 男男h啪啪无遮挡| 在线观看免费日韩欧美大片| 午夜福利欧美成人| 国产成人精品久久二区二区免费| 在线观看免费日韩欧美大片| av在线播放免费不卡| 黑人操中国人逼视频| 真人一进一出gif抽搐免费| 草草在线视频免费看| 无限看片的www在线观看| 免费在线观看黄色视频的| 亚洲成人国产一区在线观看| 岛国在线观看网站| aaaaa片日本免费| АⅤ资源中文在线天堂| 国语自产精品视频在线第100页| 国产区一区二久久| 国产欧美日韩精品亚洲av| 在线观看66精品国产|