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

    循環(huán)流化床顆粒輸送斜管的壓力脈動(dòng)特性

    2016-10-20 03:39:54曹曉陽周發(fā)戚張慧敏魏志剛魏耀東
    關(guān)鍵詞:斜管流率蝶閥

    曹曉陽,周發(fā)戚,陳 勇,張慧敏,魏志剛,3,魏耀東

    (1.中國石油大學(xué) 重質(zhì)油國家重點(diǎn)實(shí)驗(yàn)室,北京 102249; 2.浙江海洋學(xué)院 石化與能源工程學(xué)院,浙江 舟山 316000; 3.中國石油 撫順石化公司,遼寧 撫順 113008)

    ?

    循環(huán)流化床顆粒輸送斜管的壓力脈動(dòng)特性

    曹曉陽1,周發(fā)戚1,陳勇2,張慧敏1,魏志剛1,3,魏耀東1

    (1.中國石油大學(xué) 重質(zhì)油國家重點(diǎn)實(shí)驗(yàn)室,北京 102249; 2.浙江海洋學(xué)院 石化與能源工程學(xué)院,浙江 舟山 316000; 3.中國石油 撫順石化公司,遼寧 撫順 113008)

    以循環(huán)流化床輸送斜管為研究對象,通過改變顆粒質(zhì)量流率,測量斜管內(nèi)蝶閥上、下兩處的動(dòng)態(tài)壓力。結(jié)果表明,隨著顆粒質(zhì)量流率增加,蝶閥上、下方顆粒的流動(dòng)形式均發(fā)生變化,斜管內(nèi)的顆粒流態(tài)發(fā)生變化,動(dòng)態(tài)壓力的波動(dòng)幅度逐漸增加。蝶閥下方的動(dòng)態(tài)壓力標(biāo)準(zhǔn)偏差與顆粒質(zhì)量流率呈線性關(guān)系,能夠用于表征循環(huán)流化床的顆粒質(zhì)量流率。壓力脈動(dòng)的小波分析表明,斜管內(nèi)蝶閥上、下方均存在一個(gè)主頻,是由于斜管自身的流動(dòng)結(jié)構(gòu)所致,但蝶閥下方還存在一個(gè)次頻,是由于氣體脈動(dòng)以及氣、固兩相的相互作用所致,且次頻所占能量隨顆粒質(zhì)量流率增加而降低。斜管下料產(chǎn)生的脈動(dòng)壓力是一種低頻壓力,是斜管振動(dòng)的激振源。

    循環(huán)流化床;斜管;兩相流;動(dòng)態(tài)壓力;小波分析

    斜管是循環(huán)流化床顆粒循環(huán)回路中的一個(gè)關(guān)鍵部件,主要用于輸送顆粒返回流化床,同時(shí)維持系統(tǒng)的壓力平衡。例如,催化裂化裝置的再生斜管和待生斜管用于再生器和沉降器之間的待生催化劑和再生催化劑的輸送,同時(shí)維持再生器的燒焦和提升管反應(yīng)器的裂化反應(yīng)過程。到目前為止,針對循環(huán)流化床的研究主要集中在流化床和提升管上,考察其氣、固兩相流的流態(tài)及轉(zhuǎn)變條件,建立氣、固流動(dòng)參數(shù)的數(shù)學(xué)模型等,而針對斜管的研究報(bào)道非常有限。斜管內(nèi)的顆粒流動(dòng)是一種順重力的下行流動(dòng),顆粒受到垂直重力的作用和斜管器壁的約束,會(huì)沉積在器壁上,導(dǎo)致斜管內(nèi)的氣、固兩相流動(dòng)的顆粒濃度和顆粒速度在橫截面上分布的不均勻性和對稱性,氣體流動(dòng)可能是下行也可能是上行,氣、固流動(dòng)比垂直立管的更復(fù)雜。盧春喜等[1]將斜管內(nèi)氣、固兩相流態(tài)劃分為黏附滑移流動(dòng)、過渡流和充氣流動(dòng)。斜管通常連接提升管反應(yīng)器,結(jié)合部位的流動(dòng)是一個(gè)受到關(guān)注的區(qū)域。Bai等[2]對連接流化床至提升管底部的斜管不同入口結(jié)構(gòu)進(jìn)行了實(shí)驗(yàn)研究,認(rèn)為斜管上蝶閥的開度變化對提升管施加了不同強(qiáng)弱程度的入口約束;蝶閥開度越小,對入口顆粒流動(dòng)的阻力越大。Engelandt等[3]實(shí)驗(yàn)測量并分析了松動(dòng)風(fēng)對提升管Y型35°斜管入口區(qū)域的流動(dòng)的影響。Arastoopour等[4]對水平0°和傾斜45°提升管入口結(jié)構(gòu)的數(shù)值模擬結(jié)果表明,水平管底部易于產(chǎn)生顆粒的堆積,而45°斜管內(nèi)顆粒堆積減少,同時(shí)顆粒分布更加均勻。斜管的傾角對斜管內(nèi)的氣、固兩相流態(tài)有很大的影響。O’dea等[5]、Sarkar等[6-7]、Levy等[8]分別對不同傾角斜管內(nèi)的氣、固流動(dòng)進(jìn)行了實(shí)驗(yàn)研究,考察了斜管內(nèi)的流態(tài)變化,建立了斜管壓降的計(jì)算模型。斜管內(nèi)的氣、固兩相流具有很強(qiáng)的動(dòng)態(tài)特性,表現(xiàn)出壓力脈動(dòng)的變化和流動(dòng)的不穩(wěn)定性。這種不穩(wěn)定性表現(xiàn)為下行顆粒流動(dòng)形成陣發(fā)式股流,使通過斜管截面的顆粒濃度發(fā)生周期性的變化,并表現(xiàn)為壓力的低頻脈動(dòng)[9-10]。Martin等[11]測量了45°斜管的壓力脈動(dòng),并建立了顆粒循環(huán)流率與壓力脈動(dòng)的標(biāo)準(zhǔn)偏差的關(guān)系式,可實(shí)現(xiàn)顆粒循環(huán)流率的在線測量。

    循環(huán)流化床的斜管上通常安裝控制閥,用以調(diào)節(jié)顆粒循環(huán)量。閥門的開度對斜管內(nèi)的流態(tài)有重要影響,使閥前和閥后的流態(tài)存在很大的變化,導(dǎo)致了斜管內(nèi)流態(tài)的多樣性和多變性。為此,筆者采用大型循環(huán)流化床裝置,以45°斜管為研究對象,用流化催化裂化(FCC)平衡催化劑顆粒作為實(shí)驗(yàn)物料,通過改變顆粒循環(huán)流率,測量不同操作條件下斜管內(nèi)的動(dòng)態(tài)壓力,并觀察其內(nèi)部流態(tài),運(yùn)用小波分析方法分析斜管內(nèi)壓力脈動(dòng)特性,探討壓力脈動(dòng)的特點(diǎn)和組成,提高對斜管內(nèi)復(fù)雜氣、固兩相流的認(rèn)識。

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

    1.1實(shí)驗(yàn)物料

    FCC平衡催化劑為GeldartA類顆粒,顆粒物性及粒徑分布見表1。

    表1 實(shí)驗(yàn)用FCC平衡催化劑的物性及粒徑分布

    1.2實(shí)驗(yàn)裝置與分析儀器

    圖1為循環(huán)流化床實(shí)驗(yàn)裝置示意圖。流化床尺寸為φ800 mm×12000 mm,提升管尺寸為φ200 mm×12500 mm,預(yù)提升段尺寸為φ600 mm×1500 mm,斜管尺寸為φ250 mm×1300 mm,斜管傾角45°。為便于觀察斜管內(nèi)的流動(dòng)狀態(tài),采用有機(jī)玻璃制造。一級旋風(fēng)分離器筒體直徑為φ400 mm,料腿為φ150 mm×9000 mm,料腿出口為直口結(jié)構(gòu)插入流化床的密相床層內(nèi)。

    圖1 循環(huán)流化床實(shí)驗(yàn)裝置

    采用英國Gems Sensors Ltd.公司動(dòng)態(tài)壓力傳感器測量斜管內(nèi)動(dòng)態(tài)壓力,量程0~0.03 MPa,靈敏度20 Pa/mV。測得的壓力由壓力變送器轉(zhuǎn)換為1~5 V的標(biāo)準(zhǔn)電壓信號,用多功能數(shù)據(jù)采集板采樣,采樣頻率500 Hz,采樣時(shí)間120 s。斜管上壓力測量點(diǎn)示于圖2。圖2中的C、D兩測量點(diǎn)分別距離蝶閥400 mm。

    圖2 斜管測壓點(diǎn)位置

    1.3實(shí)驗(yàn)方法

    如圖1所示,催化劑顆粒循環(huán)從流化床6經(jīng)45°顆粒輸送斜管16進(jìn)入提升管下部的預(yù)提升器17,預(yù)提升器17內(nèi)提升風(fēng)攜帶催化劑顆粒垂直向上通過提升管15,在旋風(fēng)分離器10內(nèi)進(jìn)行氣、固分離后,催化劑顆粒通過其下方的料腿7下行返回流化床,從而完成一個(gè)循環(huán)。通過旋風(fēng)分離器料腿上的測量床9測量顆粒質(zhì)量流率,即在穩(wěn)定操作條件下關(guān)閉插板閥,測量一定時(shí)間內(nèi)顆粒堆積量,計(jì)算顆粒質(zhì)量流率。由斜管16上的蝶閥調(diào)節(jié)顆粒質(zhì)量流率。顆粒循環(huán)流率(Gs)范圍為0~140 kg/(m2·s),提升管入口表觀氣速(ug)為8.2 m/s。

    2 結(jié)果與討論

    2.1循環(huán)流化床顆粒輸送斜管實(shí)驗(yàn)現(xiàn)象

    實(shí)驗(yàn)中隨蝶閥開度變化,可觀察到斜管內(nèi)顆粒的流態(tài)發(fā)生不同的變化,結(jié)果示于圖3。當(dāng)?shù)y開度非常小(即催化劑顆粒質(zhì)量流率較小)時(shí),蝶閥前部的催化劑顆粒會(huì)堆滿整個(gè)斜管,形成移動(dòng)床流態(tài),顆粒密度接近堆積密度,在斜管的頂部存在向上運(yùn)動(dòng)的小氣泡,其直徑較小但個(gè)數(shù)較多,蝶閥后部顆粒因重力作用沿斜管底部形成顆粒股流動(dòng),如圖3(a)所示。隨蝶閥開度增大,顆粒質(zhì)量流率增加,蝶閥前部的顆粒流態(tài)與圖3(a)相同,但氣泡直徑變大,氣泡頻率增加,蝶閥后部沉積在斜管底部的顆粒流所占體積增加,如圖3(b)所示。當(dāng)?shù)y開度繼續(xù)增大,整個(gè)斜管內(nèi)出現(xiàn)顆粒的分層流動(dòng),顆粒沿斜管底部下行流動(dòng),斜管頂部是氣體通道,氣體沿斜管向上運(yùn)動(dòng),如圖3(c)所示。當(dāng)?shù)y開度達(dá)到某個(gè)值時(shí),顆粒質(zhì)量流率增到一定值,斜管內(nèi)顆粒流動(dòng)表現(xiàn)為滿管流化流動(dòng),斜管內(nèi)顆粒濃度均勻,顆粒和氣體整體沿斜管向下流動(dòng),如圖3(d)所示。

    2.2循環(huán)流化床顆粒輸送斜管動(dòng)態(tài)壓力

    圖4為ug=8.2 m/s時(shí),在斜管內(nèi)蝶閥下方D點(diǎn)處測量的不同顆粒質(zhì)量流率下的動(dòng)態(tài)壓力。從圖4 可見,當(dāng)Gs=0時(shí),斜管內(nèi)為純氣相,其壓力波動(dòng)的幅值很小,頻率較高,此壓力脈動(dòng)屬于提升風(fēng)的脈動(dòng),約為0.5 kPa,平均值是提升風(fēng)的壓力7.1 kPa;當(dāng)加入顆粒后,隨Gs增大,動(dòng)態(tài)壓力的平均值逐漸增加,壓力脈動(dòng)的波動(dòng)幅度也逐漸增大。加入顆粒使壓力脈動(dòng)由兩種不同成分的脈動(dòng)疊加而成,一種為脈動(dòng)頻率較高(周期可能為毫秒級)、幅值相對較小的脈動(dòng),即高頻低幅脈動(dòng);另一種是脈動(dòng)頻率較低(周期可達(dá)十幾秒)、幅值相對較大的脈動(dòng),即低頻高幅脈動(dòng)。高頻低幅脈動(dòng)主要是由于氣體本身的脈動(dòng)和顆粒簇的時(shí)聚時(shí)散產(chǎn)生,而低頻高幅脈動(dòng)由顆粒流動(dòng)過程不穩(wěn)定引起,來自蝶閥閥口的排料不穩(wěn)定。

    圖3 顆粒和氣體在斜管中的流動(dòng)狀態(tài)

    圖4 蝶閥下方D點(diǎn)處測量的不同顆粒質(zhì)量流率(Gs)下的動(dòng)態(tài)壓力

    圖5為ug=8.2 m/s時(shí),在斜管內(nèi)蝶閥上方C點(diǎn)處測量的不同顆粒質(zhì)量流率下的動(dòng)態(tài)壓力。從圖5 可見,當(dāng)Gs=0時(shí),由于蝶閥上方堆滿催化劑顆粒,壓力平均值是9.0 kPa,大于蝶閥下方的壓力平均值;隨蝶閥開度不斷增大,即顆粒質(zhì)量流率變大,有一定量的氣體竄入上行,從而使堆積的顆粒充氣形成流化流動(dòng),此處的氣體主要是從預(yù)提升器中竄入斜管的流化氣體。壓力脈動(dòng)也是由高頻低幅部分和低頻高幅部分構(gòu)成,壓力脈動(dòng)的幅度隨著顆粒質(zhì)量流率增加而增大。

    圖6為ug=8.2 m/s時(shí),斜管內(nèi)蝶閥上、下兩個(gè)測點(diǎn)在不同顆粒質(zhì)量流率下的壓力平均值。從圖6 能夠看出,隨著Gs的增加,閥下壓力從小于閥上壓力到逐漸轉(zhuǎn)變?yōu)榇笥陂y上壓力;隨著開度的增大,壓力逐漸由正壓差演變?yōu)樨?fù)壓差;轉(zhuǎn)變?yōu)樨?fù)壓差后,即Gs超過70 kg/(m2·s)時(shí),動(dòng)態(tài)壓力的波動(dòng)幅度也顯著增大,閥下的壓力波動(dòng)相比于閥上增加更為明顯(見圖4和圖5)。

    圖5 蝶閥上方C點(diǎn)處測量的不同顆粒質(zhì)量流率(Gs)下的動(dòng)態(tài)壓力

    圖6 斜管內(nèi)蝶閥上、下兩個(gè)測點(diǎn)在不同顆粒質(zhì)量流率下的壓力平均值

    2.3循環(huán)流化床壓力脈動(dòng)的標(biāo)準(zhǔn)偏差分析

    (1)

    Sd可用于表征壓力脈動(dòng)強(qiáng)度。圖7為圖4、5所示數(shù)據(jù)的Sd隨Gs的變化。從圖7可見,在純氣流條件下,碟閥上、下兩個(gè)測點(diǎn)的壓力脈動(dòng)的幅值非常??;隨Gs增大,兩個(gè)測點(diǎn)的標(biāo)準(zhǔn)偏差值增加,但增加的幅度不同。當(dāng)Gs增大到一定程度后,蝶閥下面的壓力脈動(dòng)強(qiáng)度大于上面的壓力脈動(dòng)強(qiáng)度。在實(shí)驗(yàn)范圍內(nèi),蝶閥下方的動(dòng)態(tài)壓力標(biāo)準(zhǔn)偏差與顆粒質(zhì)量流率近似呈線性關(guān)系,可以用于表征循環(huán)流化床的顆粒質(zhì)量流率并進(jìn)行在線監(jiān)測,與Martin等[11]的實(shí)驗(yàn)結(jié)果一致。但Martin實(shí)驗(yàn)測量的是閥前脈動(dòng)壓力,動(dòng)態(tài)壓力標(biāo)準(zhǔn)偏差與顆粒質(zhì)量流率呈非性線關(guān)系,并不利于建立合理的公式來預(yù)測顆粒質(zhì)量流率。由于受到蝶閥的約束作用,閥上顆粒流動(dòng)的壓力脈動(dòng)強(qiáng)度不能很好地反映管內(nèi)顆粒質(zhì)量流率,因此顆粒質(zhì)量流率與標(biāo)準(zhǔn)偏差并不具有線性關(guān)系,而閥下的顆粒由于擺脫了蝶閥的約束,其壓力脈動(dòng)強(qiáng)度能很好地反映斜管內(nèi)顆粒質(zhì)量流率,顆粒質(zhì)量流率與動(dòng)態(tài)壓力標(biāo)準(zhǔn)偏差呈線性關(guān)系說明,隨著Gs的增大,壓力脈動(dòng)強(qiáng)度隨著增強(qiáng)。此外,斜管上蝶閥前后的壓力隨著開度的增大逐漸由正壓差演變?yōu)樨?fù)壓差。正是由于蝶閥孔口變截面的約束作用,造成了斜管上負(fù)壓差操作的條件,從而使斜管內(nèi)排料具有一定的不穩(wěn)定特性[13],表現(xiàn)為下料的顆粒流量的陣發(fā)性變化。因此,蝶閥下面的壓力脈動(dòng)增加主要是由于蝶閥閥口約束形成負(fù)壓差的操作條件所致。

    圖7 斜管內(nèi)蝶閥上、下兩個(gè)測點(diǎn)壓力的標(biāo)準(zhǔn)偏差(Sd)隨顆粒循環(huán)流率(Gs)的變化

    2.4循環(huán)流化床壓力脈動(dòng)的小波分析

    小波分析是信號處理的一種方法。小波分析方法中的Daubechies小波由于具有支集緊、正交性的優(yōu)點(diǎn)[14],已被應(yīng)用于流化床壓力脈動(dòng)信號的分析[15-16]。通過對壓力信號的Daubechies小波分解誤差比較,選用誤差相對較小的db5小波基對圖4和圖5中的動(dòng)態(tài)壓力信號進(jìn)行14層多尺度分解。原始信號被分解成細(xì)節(jié)信號D1和近似信號A1。近似信號A1再被分解成A2和D2,如此往復(fù)進(jìn)行分解。重構(gòu)小波的傳遞是從細(xì)節(jié)尺度D1到近似尺度D14和A14。細(xì)節(jié)尺度對應(yīng)的是原始信號的高頻部分,而近似尺度對應(yīng)的是原始信號的低頻部分[17]。本研究中,筆者采用細(xì)節(jié)能量分率[18-19]表征催化劑顆粒在斜管內(nèi)的壓力脈動(dòng)特性。

    圖8為ug=8.2 m/s時(shí),不同顆粒循環(huán)流率下蝶閥下方D點(diǎn)處動(dòng)態(tài)壓力信號的小波能量百分比。從圖8可看出,當(dāng)Gs=0時(shí),即在純氣體工況下,動(dòng)態(tài)壓力信號的小波能量只有一個(gè)峰值,即只有一個(gè)主頻,來自于氣體的脈動(dòng);當(dāng)有催化劑顆粒加入時(shí),動(dòng)態(tài)壓力信號的小波能量存在兩個(gè)峰值。較低峰值能量對應(yīng)的頻率為次頻,較高峰值能量對應(yīng)的頻率為主頻,主、次兩個(gè)頻段均屬于低頻脈動(dòng)。次頻是D8~D10(0.4883~1.9531 Hz),由斜管內(nèi)氣、固兩相間的相互作用產(chǎn)生,基本不隨顆粒質(zhì)量流率的變化而變化,其小波能量峰值隨著顆粒質(zhì)量流率增加而降低;主頻是D12~D13(0.0610~0.2441 Hz),基本不隨顆粒質(zhì)量流率的變化而變化,其小波能量峰值隨顆粒質(zhì)量流率增加而增大。主頻由斜管閥控機(jī)制引起的顆粒不穩(wěn)定下料所產(chǎn)生,隨蝶閥開度增大,顆粒質(zhì)量流率增加,壓力脈動(dòng)強(qiáng)度會(huì)逐漸增強(qiáng),因此其小波能量分率有所加強(qiáng)。

    圖8 不同顆粒質(zhì)量流率下蝶閥下方動(dòng)態(tài)壓力信號的小波能量百分比

    圖9為ug=8.2 m/s時(shí),不同顆粒循環(huán)質(zhì)量流率下蝶閥上方C點(diǎn)處動(dòng)態(tài)壓力信號的小波能量百分比。圖9表明,動(dòng)態(tài)壓力信號的小波能量只有一個(gè)峰值,即只有一個(gè)主頻,沒有次頻,而主頻段是D11~D14(0.0305~0.4883 Hz),隨著顆粒質(zhì)量流率的增大,主頻從0.0305 Hz逐漸增大到0.4883 Hz,但是小波能量分率降低。主頻主要是斜管閥控機(jī)制引起的顆粒不穩(wěn)定下料所致,隨蝶閥開度增大,孔口約束作用減小,因此其小波能量分率有所降低,說明蝶閥形成的孔口排料作用在降低,此主頻脈動(dòng)屬低頻脈動(dòng)。

    圖9 不同顆粒質(zhì)量流率下蝶閥上方動(dòng)態(tài)壓力信號的小波能量百分比

    2.5循環(huán)流化床顆粒輸送斜管壓力脈動(dòng)

    輸送斜管的動(dòng)態(tài)壓力測量和分析結(jié)果表明,隨著Gs的增加,斜管上兩個(gè)測點(diǎn)的動(dòng)態(tài)壓力波動(dòng)幅度增大,尤其是閥下的壓力脈動(dòng)幅度增長較顯著,最大波動(dòng)幅值可達(dá)1.5 kPa,會(huì)造成斜管較大程度的振動(dòng),且其振動(dòng)強(qiáng)度會(huì)隨著壓力的波動(dòng)幅度的增大而加強(qiáng)。

    斜管內(nèi)顆粒下行流動(dòng)的不穩(wěn)定導(dǎo)致了壓力脈動(dòng),這種壓力脈動(dòng)對工業(yè)裝置有兩方面的影響。一方面是,壓力脈動(dòng)的幅度與顆粒循環(huán)流率密切相關(guān),脈動(dòng)壓力的標(biāo)準(zhǔn)偏差Sd與Gs基本呈線性關(guān)系,因此,在結(jié)構(gòu)與入口速度一定的條件下,標(biāo)準(zhǔn)偏差Sd是顆粒質(zhì)量流率Gs的單值函數(shù)。通過對脈動(dòng)壓力信號的標(biāo)準(zhǔn)偏差分析可以獲得循環(huán)流化床內(nèi)顆粒循環(huán)流率數(shù)據(jù),用于工業(yè)裝置運(yùn)行參數(shù)的測量和監(jiān)視。

    另一方面,脈動(dòng)壓力的小波分析表明,閥前、后的主、次兩個(gè)頻段均屬于低頻脈動(dòng)。在工業(yè)裝置中,顆粒輸送斜管一般采用柔性設(shè)計(jì)結(jié)構(gòu),斜管比較細(xì)長,目的是吸收溫差產(chǎn)生的應(yīng)力,結(jié)果斜管本身的固有頻率比較低。由于斜管下料產(chǎn)生的脈動(dòng)壓力是一種低頻壓力變化,形成了斜管振動(dòng)的激振源,易于造成斜管的振動(dòng)斷裂[20-24]。因此,為了防止斜管發(fā)生共振響應(yīng),應(yīng)將斜管的固有頻率盡量遠(yuǎn)離斜管下料時(shí)脈動(dòng)壓力的主頻和次頻。可以通過在斜管上增加支撐結(jié)構(gòu),改變其固有頻率或改變操作參數(shù)的范圍,規(guī)避斜管發(fā)生共振。

    3 結(jié) 論

    (1)隨斜管內(nèi)顆粒質(zhì)量流率Gs的增大,蝶閥上方催化劑顆粒流態(tài)依次呈現(xiàn)移動(dòng)床流動(dòng),氣、固分層流動(dòng),滿管流化流動(dòng);而蝶閥下方的催化劑顆粒態(tài)依次為顆粒股流,氣、固分層流動(dòng),滿管流化流動(dòng)。

    (2)隨顆粒質(zhì)量流率Gs的增加,斜管內(nèi)壓力脈動(dòng)的波動(dòng)幅度增大,蝶閥下方的動(dòng)態(tài)壓力的標(biāo)準(zhǔn)偏差與Gs呈較好的線性關(guān)系,可以用壓力信號的標(biāo)準(zhǔn)偏差來在線檢測循環(huán)流化床內(nèi)的顆粒質(zhì)量流率。

    (3)小波分析表明,斜管內(nèi)蝶閥上方和下方均存在一個(gè)主頻,由蝶閥及斜管自身結(jié)構(gòu)引起的壓力脈動(dòng)所致;蝶閥下方還有一個(gè)次頻,由斜管內(nèi)氣體脈動(dòng)以及氣、固兩相的相互作用所致,其小波能量分率隨顆粒質(zhì)量流率的增加而降低。

    (4)斜管下料時(shí)產(chǎn)生的脈動(dòng)壓力均屬于低頻壓力,是造成斜管振動(dòng)的激振源。工業(yè)中斜管的固有頻率比較低,應(yīng)將斜管的固有頻率盡量遠(yuǎn)離斜管下料時(shí)脈動(dòng)壓力的主頻和次頻。

    [1]盧春喜,王祝安.催化裂化流態(tài)化技術(shù)[M].北京:中國石化出版社,2002:226-235.

    [2]BAI D R,JIN Y,YU Z Q,et al.The axial distribution of the cross-sectionally averaged voidage in fast fluidized beds[J].Powder Technology,1992,71(1):51-58.

    [3]ENGELANDT G V,WILDE J D,HEYNDERICKX G J,et al.Experimental study of inlet phenomena of inclined non-aerated and aerated Y-inlets in a dilute cold-flow riser[J].Chemical Engineering Science,2007,62(1/2):339-355.

    [4]ARASTOOPOUR H.Numerical simulation and experimental analysis of gas/solid flow system:1999 Fluor-daniel plenary lecture[J].Powder Technology,2001,119(2/3):59-67.

    [5]O’DEA D P,RUDOLPH V,CHONG Y O.Gas-solids flow through the bottom restriction of an inclined standpipe[J].Powder Technology,1990,62(3):291-297.

    [6]SARKAR M,GUPTA S K,SARKAR M K.Experimental investigation on gravity flow of solids through inclined pipes[J].Chemical Engineering Science,1991,46(4):1137-1144.

    [7]SARKAT M,GUPTA S K,SARKAR M K.An experimental investigation of the flow of solids from a fluidized bed through an inclined pipe[J].Powder Technology,1991,64(3):221-231.

    [8]LEVY A,MOONEY T,MARJANOVIC P,et al.A comparison of analytical and numerical models with experimental data for gas-solid flow through a straight pipe at different inclinations[J].Powder Technology,1997,93(3):253-260.

    [9]張毅,彭園園,魏耀東,等.循環(huán)流化床下料立管內(nèi)氣固兩相流動(dòng)狀態(tài)與壓力脈動(dòng)的關(guān)系[J].過程工程學(xué)報(bào),2008,8(1):23-27.(ZHANG Yi,PENG Yuanyuan,WEI Yaodong,et al.Relationships between pressure fluctuation and flow patterns in standpipes of circulating fluidized bed[J].The Chinese Journal of Process Engineering,2008,8(1):23-27.)

    [10]張毅,魏耀東,時(shí)銘顯.氣固循環(huán)流化床負(fù)壓差下料立管的壓力脈動(dòng)特性[J].化工學(xué)報(bào),2007,58(6):1417-1420.(ZHANG Yi,WEI Yaodong,SHI Mingxian.Characteristics of pressure fluctuation in standpipe at negative pressure gradient in circulating fluidized bed[J].CIESC Journal,2007,58(6):1417-1420.)

    [11]MARTIN L D,OMMEN J R V.Estimation of the overall mass flux in inclined standpipes by means of pressure fluctuation measurements[J].Chemical Engineering Journal,2012,204-206:125-130.

    [13]NAGASHIMA H,ISHIKURA T,IDE M.Flow characteristics of a small moving bed downcomer with an orifice under negative pressure gradient[J].Powder Technology,2009,192(1):110-115.

    [14]趙貴兵,陳紀(jì)忠,陽永榮.Daubechies小波對流化床壓力波動(dòng)的分解研究[J].高?;瘜W(xué)工程學(xué)報(bào),2001,17(3):272-278.(ZHAO Guibing,CHEN Jizhong,YANG Yongrong.Study on decomposition of pressure fluctuations using a series of daubechies wavelets in a fluidized bed[J].Journal of Chemical Engineering of Chinese Universities,2001,17(3):272-278.)

    [15]冀海峰,黃志堯,吳賢國.基于小波變換的氣固流化床壓力波動(dòng)信號的分析[J].高校化學(xué)工程學(xué)報(bào),2000,14(6):553-557.(JI Haifeng,HUANG Zhiyao,WU Xianguo.Analysis of pressure fluctuation signal of gas-solidfluidized bed based on wavelets transform[J].Journal of Chemical Engineering of Chinese Universities,2000,14(6):553-557.)

    [16]黃海,黃軼倫,張衛(wèi)東.氣固流化床壓力脈動(dòng)信號的相關(guān)結(jié)構(gòu)模型與分析[J].化工學(xué)報(bào),1999,50(6):812-817.(HUANG Hai,HUANG Yilun,ZHANG Weidong.Modeling and analysis of pressure fluctuations in a gas-solid fluidized bed[J].CIESC Journal,1999,50(6):812-817.)

    [17]SHOU M C,LEU L P.Energy of power spectral density function and wavelet analysis of absolute pressure fluctuation measurements in fluidized beds[J].Chemical Engineering Research and Design,2005,83(5):478-491.

    [18]陳恒志,陳小翠,李洪鐘.提升管內(nèi)氣固流動(dòng)的小波分析[J].化學(xué)工程,2009,37(8):20-23.(CHEN Hengzhi,CHEN Xiaocui,LI Hongzhong.Study on gas-solids flow in riser based on wavelets analysis[J].Chemical Engineering (China),2009,37 (8):20-23.)

    [19]周發(fā)戚,陳勇,魏志剛,等.循環(huán)流化床提升管T形彎頭動(dòng)態(tài)壓力的小波分析[J].化工學(xué)報(bào),2015,66(5):1697-1703.(ZHOU Faqi,CHEN Yong,WEI Zhigang,et al.Wavelet analysis of dynamic pressure in T-abrupt of CFB riser[J].CIESC Journal,2015,66(5):1697-1703.)

    [20]許緘濤,王軍峰,羅萬明.FCC裝置半再生斜管流化不穩(wěn)定的原因分析[J].齊魯石油化工,2008,36(4):297-300.(XU Jiantao,WANG Junfeng,LUO Wanming.Cause analysis of unstable fluidization in semi-regeneration chute of FCC unit[J].Qilu Petrochemical Technology,2008,36(4):297-300.)

    [21]李健.催化裂化反應(yīng)再生系統(tǒng)斜管上松動(dòng)點(diǎn)的合理設(shè)置[J].煉油技術(shù)與工程,2003,33(9):16-18.(LI Jian.Reasonable aeration tap setting in the slope tube of FCC reaction-regeneration system[J].Petroleum Refinery Engineering,2003,33(9):16-18.)

    [22]劉瑞豐.催化裂化再生斜管安定性分析[J].煉油設(shè)計(jì),1996,26(1):57-60.(LIU Ruifeng.Stability analysis of regeneration sloped pipe in FCCU[J].Petroleum Refinery Engineering,1996,26(1):57-60.)

    [23]王恒,耿興東.催化裂化裝置待生斜管流化異常原因分析[J].煉油與化工,2013,24(2):29-30.(WANG Heng,GENG Xingdong.Cause for fluidizing abnormity of to-be-generated sloping tube in FCCU[J].Refining and Chemical Industry,2013,24(2):29-30.)

    [24]劉春貴,馬俊,李慶文,等.重油催化裂化裝置再生斜管流化效果不好的技術(shù)改造[J].石化技術(shù)與應(yīng)用,2013,31(1):44-46.(LIU Chungui,MA Jun,LI Qingwen,et al.Revamping of regenerator sloped tube in a residue fluidized-bed catalytic cracking unit for improvement of fluidizing effectiveness[J].Petrochemical Technology & Application,2013,31(1):44-46.)

    Characteristics of Pressure Fluctuations in the Particle-TransportInclined Standpipe of a Circulating Fluidized Bed

    CAO Xiaoyang1,ZHOU Faqi1,CHEN Yong2,ZHANG Huimin1,WEI Zhigang1,3,WEI Yaodong1

    (1.State Key Laboratory of Heavy Oil Processing,China University of Petroleum,Beijing 102249,China; 2.School of Petrochemical and Energy Engineering,Zhejiang Ocean University,Zhoushan 316000,China; 3.Fushun Petrochemical Company,CNPC,Fushun 113008,China)

    In large circulating fluidized bed device,the characteristics of pressure fluctuation in a particle-transport inclined standpipe under different solid circulating flux were investigated based on measuring the dynamic pressure of a gas-solid two-phase flow by using FCC equilibrium catalyst as the particles.The experimental results showed that with increase of particle mass flux,the flow patterns of FCC particles above and below the butterfly valve in the inclined standpipe both successively changed,and the fluctuation amplitude of dynamic pressures increased.The standard deviation of the dynamic pressures below the butterfly valve possessed a linear relation with the particle mass flux,which could be used to characterize the particle mass flux in the inclined standpipe.By wavelet analysis,it is indicated that there existed a main frequency caused by inclined standpipe flow structures in the pressure fluctuations above and below the butterfly valve,besides,there was a secondary frequency in the pressure fluctuation below the butterfly valve caused by gas turbulence and gas-solid interactions.With the increase of particle mass flux,the wavelet energy ratio of the secondary frequency gradually decreased.The fluctuation pressures caused by discharge in the inclined standpipe is a low frequency pressure and the excitation source of inclined standpipe vibration.

    circulating fluidized bed;inclined standpipe;two-phase flow; pressure fluctuations;wavelet analysis

    2015-09-16

    國家自然科學(xué)基金項(xiàng)目(21176250)資助

    曹曉陽,男,碩士研究生,從事石油化工過程裝備方面的研究

    魏耀東,男,教授,博士,從事氣固分離、流態(tài)化及化工過程裝備優(yōu)化等方面的研究;Tel:010-89733939;E-mail:weiyd@cup.edu.cn

    1001-8719(2016)05-0913-08

    TQ022.3

    Adoi:10.3969/j.issn.1001-8719.2016.05.007

    猜你喜歡
    斜管流率蝶閥
    CDN 緩存空間調(diào)整優(yōu)化方法
    恒壓排氣制動(dòng)蝶閥的結(jié)構(gòu)及應(yīng)用簡析
    帶交換和跳躍的一維雙向自驅(qū)動(dòng)系統(tǒng)的仿真研究
    催化裂化裝置輸送斜管內(nèi)催化劑流化狀態(tài)分析
    Predictive factors of histological response of colorectal liver metastases after neoadjuvant chemotherapy
    多束光在皮膚組織中傳輸?shù)拿商乜_模擬
    聚合物氣輔共擠中熔體體積流率對擠出脹大的影響
    基于蝶閥的優(yōu)化設(shè)計(jì)研究
    斜管沉淀池處理轉(zhuǎn)爐除塵水斜管內(nèi)積泥原因分析及對策
    河南冶金(2015年2期)2015-12-22 06:21:08
    風(fēng)機(jī)進(jìn)口蝶閥改用電器、機(jī)械控制
    亚洲精品中文字幕在线视频| 国产黄色免费在线视频| 免费高清在线观看日韩| 久久人妻福利社区极品人妻图片| 中国美女看黄片| 在线观看免费午夜福利视频| 国产区一区二久久| 久久影院123| 亚洲国产毛片av蜜桃av| 每晚都被弄得嗷嗷叫到高潮| 99国产极品粉嫩在线观看| 欧美午夜高清在线| 老汉色av国产亚洲站长工具| 欧美激情高清一区二区三区| 亚洲欧美一区二区三区黑人| 欧美日韩福利视频一区二区| 久久久精品94久久精品| 王馨瑶露胸无遮挡在线观看| 国产精品久久久久成人av| 九色亚洲精品在线播放| 欧美 日韩 精品 国产| 国产成人欧美| 国产精品久久久av美女十八| 久久中文看片网| 香蕉丝袜av| 两性夫妻黄色片| 精品久久蜜臀av无| 黑人欧美特级aaaaaa片| 热99国产精品久久久久久7| 国产无遮挡羞羞视频在线观看| 国产成人系列免费观看| 国产一区二区三区综合在线观看| 最近最新中文字幕大全免费视频| 欧美日韩视频精品一区| 欧美黄色淫秽网站| 久久ye,这里只有精品| 久久午夜综合久久蜜桃| 啦啦啦中文免费视频观看日本| 黑人操中国人逼视频| 岛国在线观看网站| 又紧又爽又黄一区二区| 黄网站色视频无遮挡免费观看| 久久久久久久国产电影| 97在线人人人人妻| 亚洲专区中文字幕在线| 日韩制服丝袜自拍偷拍| 在线观看一区二区三区激情| 男女无遮挡免费网站观看| 国产精品免费一区二区三区在线 | 国产成+人综合+亚洲专区| 亚洲中文字幕日韩| 人人妻人人爽人人添夜夜欢视频| 国产黄频视频在线观看| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久 | 中文字幕色久视频| 国产精品99久久99久久久不卡| 一区福利在线观看| 99riav亚洲国产免费| 成人影院久久| 日本黄色视频三级网站网址 | 欧美乱妇无乱码| 最新美女视频免费是黄的| 麻豆成人av在线观看| 国产av一区二区精品久久| 精品福利永久在线观看| 国产日韩一区二区三区精品不卡| 成人av一区二区三区在线看| 免费高清在线观看日韩| 亚洲精品国产一区二区精华液| 激情在线观看视频在线高清 | 日韩人妻精品一区2区三区| 国产精品久久电影中文字幕 | 午夜视频精品福利| 久久中文字幕人妻熟女| 免费av中文字幕在线| 久久国产精品影院| 99re6热这里在线精品视频| 成年人黄色毛片网站| 人人澡人人妻人| 久久婷婷成人综合色麻豆| 99国产精品免费福利视频| 波多野结衣一区麻豆| 日日夜夜操网爽| 999久久久国产精品视频| 国产精品.久久久| 无人区码免费观看不卡 | 大型av网站在线播放| 国产亚洲欧美精品永久| 一级a爱视频在线免费观看| 少妇被粗大的猛进出69影院| 欧美人与性动交α欧美精品济南到| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 一个人免费在线观看的高清视频| 色综合欧美亚洲国产小说| 国产精品国产高清国产av | 一边摸一边抽搐一进一小说 | 又紧又爽又黄一区二区| 亚洲中文av在线| 视频区图区小说| 99久久99久久久精品蜜桃| 丁香六月天网| 十八禁人妻一区二区| 50天的宝宝边吃奶边哭怎么回事| 国产男靠女视频免费网站| 51午夜福利影视在线观看| 免费久久久久久久精品成人欧美视频| 婷婷成人精品国产| 精品国产乱码久久久久久男人| 可以免费在线观看a视频的电影网站| 一边摸一边抽搐一进一出视频| 国产不卡一卡二| 亚洲一区二区三区欧美精品| 啦啦啦在线免费观看视频4| 久久人人97超碰香蕉20202| 国产精品av久久久久免费| 中文字幕色久视频| 在线观看www视频免费| 日本精品一区二区三区蜜桃| 99热网站在线观看| 欧美成人午夜精品| 欧美日韩中文字幕国产精品一区二区三区 | 午夜视频精品福利| 免费在线观看完整版高清| 久久性视频一级片| 久久精品亚洲精品国产色婷小说| 国产淫语在线视频| 69精品国产乱码久久久| 亚洲国产欧美网| 国产日韩欧美在线精品| 久久午夜亚洲精品久久| 又黄又粗又硬又大视频| 亚洲精品一二三| 国产av又大| 美女福利国产在线| 国产成人系列免费观看| tocl精华| 女性被躁到高潮视频| 91麻豆精品激情在线观看国产 | 高清毛片免费观看视频网站 | 免费在线观看黄色视频的| 真人做人爱边吃奶动态| 久久狼人影院| 久久国产精品大桥未久av| 久久香蕉激情| 夜夜骑夜夜射夜夜干| 每晚都被弄得嗷嗷叫到高潮| 免费观看人在逋| 午夜激情久久久久久久| 99久久99久久久精品蜜桃| 一二三四在线观看免费中文在| 我的亚洲天堂| av不卡在线播放| 老司机靠b影院| 国精品久久久久久国模美| 国产免费视频播放在线视频| 色老头精品视频在线观看| 热re99久久精品国产66热6| 乱人伦中国视频| 午夜福利欧美成人| 成人免费观看视频高清| a级毛片在线看网站| 午夜精品久久久久久毛片777| 亚洲成国产人片在线观看| 一级毛片电影观看| 亚洲av片天天在线观看| 精品国产超薄肉色丝袜足j| 国产精品久久久av美女十八| 国产精品av久久久久免费| 亚洲一卡2卡3卡4卡5卡精品中文| 香蕉久久夜色| 精品少妇黑人巨大在线播放| 一本色道久久久久久精品综合| 久久国产精品影院| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品久久午夜乱码| 亚洲精品久久成人aⅴ小说| 成年人黄色毛片网站| 咕卡用的链子| 捣出白浆h1v1| 夫妻午夜视频| 99国产极品粉嫩在线观看| 日韩人妻精品一区2区三区| 日日爽夜夜爽网站| 在线永久观看黄色视频| 亚洲成a人片在线一区二区| 黄色视频在线播放观看不卡| 亚洲九九香蕉| 亚洲色图综合在线观看| h视频一区二区三区| 亚洲熟女精品中文字幕| 国产精品1区2区在线观看. | 亚洲精品国产精品久久久不卡| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三区在线| 又大又爽又粗| av不卡在线播放| 国产精品二区激情视频| 夫妻午夜视频| 日韩制服丝袜自拍偷拍| 香蕉国产在线看| 国产亚洲精品一区二区www | 热99re8久久精品国产| a级毛片黄视频| 国产在线免费精品| 亚洲avbb在线观看| 午夜免费成人在线视频| 国产av又大| 免费日韩欧美在线观看| 精品少妇黑人巨大在线播放| 黄色丝袜av网址大全| av天堂在线播放| 久久精品成人免费网站| 午夜日韩欧美国产| 天堂动漫精品| 亚洲欧洲精品一区二区精品久久久| 国产在线一区二区三区精| 丝袜人妻中文字幕| a级毛片黄视频| 国产三级黄色录像| 国产精品亚洲av一区麻豆| 午夜免费鲁丝| 精品午夜福利视频在线观看一区 | 欧美激情久久久久久爽电影 | 操美女的视频在线观看| 老熟妇乱子伦视频在线观看| 国产一卡二卡三卡精品| 制服诱惑二区| 热re99久久国产66热| 国产精品久久久久久精品电影小说| 巨乳人妻的诱惑在线观看| 夜夜爽天天搞| 91大片在线观看| 激情在线观看视频在线高清 | 午夜福利一区二区在线看| 岛国毛片在线播放| 91国产中文字幕| 欧美中文综合在线视频| 亚洲va日本ⅴa欧美va伊人久久| 欧美+亚洲+日韩+国产| 欧美激情极品国产一区二区三区| 久久午夜亚洲精品久久| 正在播放国产对白刺激| 曰老女人黄片| 午夜免费成人在线视频| 亚洲人成77777在线视频| 极品人妻少妇av视频| 黑人操中国人逼视频| 午夜福利欧美成人| 免费观看人在逋| 2018国产大陆天天弄谢| 欧美日韩精品网址| 好男人电影高清在线观看| 19禁男女啪啪无遮挡网站| 99精品久久久久人妻精品| 国产野战对白在线观看| 十八禁网站网址无遮挡| 亚洲熟女精品中文字幕| 国产精品一区二区免费欧美| 欧美成人午夜精品| 热99久久久久精品小说推荐| 久久久久久久久久久久大奶| 国产精品久久电影中文字幕 | 午夜免费鲁丝| 久久久久久久国产电影| 亚洲九九香蕉| 高清av免费在线| 丰满饥渴人妻一区二区三| 啦啦啦 在线观看视频| 久久精品亚洲av国产电影网| 法律面前人人平等表现在哪些方面| 免费在线观看黄色视频的| 欧美 亚洲 国产 日韩一| 久久av网站| 两人在一起打扑克的视频| 一级黄色大片毛片| 午夜福利免费观看在线| 国产亚洲午夜精品一区二区久久| 亚洲欧美日韩另类电影网站| 露出奶头的视频| 亚洲欧洲日产国产| 不卡一级毛片| 狠狠狠狠99中文字幕| 啦啦啦视频在线资源免费观看| 无遮挡黄片免费观看| 日韩大码丰满熟妇| 天天躁日日躁夜夜躁夜夜| 黄片小视频在线播放| 国产成人一区二区三区免费视频网站| 中文字幕精品免费在线观看视频| 少妇被粗大的猛进出69影院| 精品久久久久久电影网| 中文亚洲av片在线观看爽 | 亚洲精品美女久久久久99蜜臀| 在线av久久热| 久久久精品区二区三区| 日韩欧美免费精品| 淫妇啪啪啪对白视频| 无遮挡黄片免费观看| 91国产中文字幕| 中文字幕高清在线视频| 日韩欧美一区视频在线观看| 日本av免费视频播放| 久久久久网色| 十八禁网站免费在线| 国产成人免费观看mmmm| av免费在线观看网站| 欧美老熟妇乱子伦牲交| 后天国语完整版免费观看| 国产一区二区 视频在线| 久久婷婷成人综合色麻豆| 在线天堂中文资源库| 日韩一卡2卡3卡4卡2021年| 高清欧美精品videossex| av在线播放免费不卡| 色精品久久人妻99蜜桃| 人人妻人人添人人爽欧美一区卜| 国产高清视频在线播放一区| 成人免费观看视频高清| 国产成人欧美在线观看 | 建设人人有责人人尽责人人享有的| 99在线人妻在线中文字幕 | 搡老乐熟女国产| 美女午夜性视频免费| 首页视频小说图片口味搜索| 丁香六月天网| 两人在一起打扑克的视频| www.自偷自拍.com| 老司机午夜十八禁免费视频| 免费看a级黄色片| 18禁黄网站禁片午夜丰满| 国产一区二区三区视频了| 嫁个100分男人电影在线观看| 黄色视频在线播放观看不卡| 一级,二级,三级黄色视频| 成年版毛片免费区| 亚洲av欧美aⅴ国产| 五月天丁香电影| 黄频高清免费视频| 亚洲自偷自拍图片 自拍| 极品教师在线免费播放| 9191精品国产免费久久| 国产av精品麻豆| 王馨瑶露胸无遮挡在线观看| 99国产综合亚洲精品| 精品高清国产在线一区| 亚洲性夜色夜夜综合| 色播在线永久视频| 国内毛片毛片毛片毛片毛片| av不卡在线播放| 性色av乱码一区二区三区2| netflix在线观看网站| av有码第一页| netflix在线观看网站| 国产一区二区三区在线臀色熟女 | 欧美久久黑人一区二区| www.999成人在线观看| 黄片播放在线免费| 黄色视频不卡| 亚洲成人手机| 免费观看人在逋| 又黄又粗又硬又大视频| 大型av网站在线播放| 成年人午夜在线观看视频| 高清黄色对白视频在线免费看| 又黄又粗又硬又大视频| 久久久久久久精品吃奶| 水蜜桃什么品种好| 18在线观看网站| 亚洲精品成人av观看孕妇| 麻豆av在线久日| 国产精品亚洲一级av第二区| 国产黄频视频在线观看| 午夜视频精品福利| 日韩一卡2卡3卡4卡2021年| av天堂在线播放| 夜夜夜夜夜久久久久| 午夜福利视频在线观看免费| 精品欧美一区二区三区在线| 久久久精品免费免费高清| 成人18禁在线播放| 黄片小视频在线播放| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区| 亚洲午夜理论影院| 日本a在线网址| 五月天丁香电影| 国产精品久久久久久人妻精品电影 | 日本五十路高清| 性少妇av在线| 欧美精品高潮呻吟av久久| 自拍欧美九色日韩亚洲蝌蚪91| 精品亚洲成国产av| 久久久国产一区二区| 久热这里只有精品99| 国产精品 欧美亚洲| 一区二区日韩欧美中文字幕| 欧美久久黑人一区二区| 美女高潮到喷水免费观看| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 亚洲色图av天堂| 久久精品成人免费网站| 丝袜在线中文字幕| aaaaa片日本免费| 男女午夜视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品一卡2卡三卡4卡5卡| a级毛片在线看网站| 欧美成人午夜精品| 99精品久久久久人妻精品| 日韩欧美国产一区二区入口| 亚洲第一欧美日韩一区二区三区 | 性少妇av在线| 国产不卡一卡二| 国产片内射在线| 欧美日韩一级在线毛片| 免费在线观看影片大全网站| 日韩一卡2卡3卡4卡2021年| 免费黄频网站在线观看国产| 国产成人啪精品午夜网站| 国产一卡二卡三卡精品| 亚洲精品av麻豆狂野| 成人亚洲精品一区在线观看| 久久精品亚洲熟妇少妇任你| 亚洲精华国产精华精| 亚洲成人免费电影在线观看| 国产精品美女特级片免费视频播放器 | 久久亚洲真实| 精品国产亚洲在线| 国产麻豆69| 色在线成人网| 午夜成年电影在线免费观看| 91精品三级在线观看| 中文字幕高清在线视频| 男女免费视频国产| 水蜜桃什么品种好| 亚洲色图 男人天堂 中文字幕| 久久人人97超碰香蕉20202| 久久久久久亚洲精品国产蜜桃av| 精品人妻在线不人妻| 多毛熟女@视频| 1024香蕉在线观看| 三上悠亚av全集在线观看| 黄色 视频免费看| 国产免费福利视频在线观看| 99精品在免费线老司机午夜| 日韩视频在线欧美| 亚洲第一av免费看| 青草久久国产| 国产精品99久久99久久久不卡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图综合在线观看| 黄片小视频在线播放| 国产高清videossex| 欧美在线黄色| 国产成人免费观看mmmm| 久久人人97超碰香蕉20202| 日本av免费视频播放| 亚洲精品成人av观看孕妇| 夜夜爽天天搞| 久久久久网色| 免费观看av网站的网址| 9191精品国产免费久久| 国产在线免费精品| 精品少妇久久久久久888优播| 国产精品免费一区二区三区在线 | 一级片免费观看大全| 757午夜福利合集在线观看| 色精品久久人妻99蜜桃| 亚洲av日韩精品久久久久久密| 日日摸夜夜添夜夜添小说| 亚洲精品久久成人aⅴ小说| 大香蕉久久成人网| 中文亚洲av片在线观看爽 | 久久久久国内视频| 一二三四在线观看免费中文在| 久久天堂一区二区三区四区| 国产人伦9x9x在线观看| a级片在线免费高清观看视频| 国产一区二区 视频在线| 久久久久精品人妻al黑| 啦啦啦视频在线资源免费观看| 午夜福利在线免费观看网站| 看免费av毛片| 亚洲熟妇熟女久久| xxxhd国产人妻xxx| 国产午夜精品久久久久久| 王馨瑶露胸无遮挡在线观看| 一区二区日韩欧美中文字幕| 国产高清国产精品国产三级| 在线观看人妻少妇| 欧美黄色淫秽网站| 91成人精品电影| 精品欧美一区二区三区在线| 少妇被粗大的猛进出69影院| 欧美久久黑人一区二区| 19禁男女啪啪无遮挡网站| 亚洲av成人一区二区三| 精品久久久久久久毛片微露脸| 午夜福利影视在线免费观看| 国产亚洲精品久久久久5区| 一本综合久久免费| 人成视频在线观看免费观看| 欧美激情久久久久久爽电影 | 日日夜夜操网爽| 一个人免费在线观看的高清视频| 一区二区三区国产精品乱码| 国产精品.久久久| 香蕉久久夜色| 亚洲精品中文字幕一二三四区 | 高清黄色对白视频在线免费看| 在线十欧美十亚洲十日本专区| 男女下面插进去视频免费观看| 亚洲精品一卡2卡三卡4卡5卡| 黄片大片在线免费观看| 日本黄色视频三级网站网址 | 亚洲综合色网址| 久久国产精品男人的天堂亚洲| 欧美成人免费av一区二区三区 | 韩国精品一区二区三区| 国产一区二区在线观看av| 日本五十路高清| 久久久久国产一级毛片高清牌| 人妻 亚洲 视频| 欧美大码av| 欧美激情极品国产一区二区三区| 大片免费播放器 马上看| 久久中文字幕人妻熟女| 无限看片的www在线观看| 欧美精品一区二区大全| 国产亚洲av高清不卡| 欧美成狂野欧美在线观看| 中文字幕制服av| av福利片在线| 亚洲少妇的诱惑av| 一级片'在线观看视频| 免费少妇av软件| 一级毛片电影观看| 国产伦理片在线播放av一区| 青青草视频在线视频观看| 日韩熟女老妇一区二区性免费视频| 欧美乱码精品一区二区三区| 18禁美女被吸乳视频| 老熟妇乱子伦视频在线观看| 十八禁人妻一区二区| 飞空精品影院首页| 色婷婷av一区二区三区视频| 午夜福利欧美成人| 久久ye,这里只有精品| 国产欧美日韩一区二区三| 国产99久久九九免费精品| 天天添夜夜摸| 国产老妇伦熟女老妇高清| 日韩人妻精品一区2区三区| 黑人巨大精品欧美一区二区mp4| 亚洲九九香蕉| 国产精品98久久久久久宅男小说| 国产伦人伦偷精品视频| 亚洲成人免费av在线播放| 国产国语露脸激情在线看| 亚洲精品国产区一区二| 国产精品秋霞免费鲁丝片| 亚洲全国av大片| 三级毛片av免费| 久久久久久人人人人人| 桃花免费在线播放| 午夜福利一区二区在线看| 黄片播放在线免费| 女性生殖器流出的白浆| 亚洲一区二区三区欧美精品| 久久久精品区二区三区| 国产成人一区二区三区免费视频网站| 国产精品麻豆人妻色哟哟久久| www.精华液| 日本黄色日本黄色录像| 一区二区三区乱码不卡18| 亚洲av第一区精品v没综合| 在线 av 中文字幕| cao死你这个sao货| 一边摸一边做爽爽视频免费| 午夜精品久久久久久毛片777| 99国产极品粉嫩在线观看| 一级片'在线观看视频| 80岁老熟妇乱子伦牲交| 国产精品一区二区免费欧美| 岛国毛片在线播放| 久久精品亚洲精品国产色婷小说| 电影成人av| 蜜桃国产av成人99| 一区二区av电影网| 日韩视频在线欧美| 丰满迷人的少妇在线观看| 女人久久www免费人成看片| 高清欧美精品videossex| 男女边摸边吃奶| 中文字幕制服av| 亚洲第一欧美日韩一区二区三区 | 亚洲国产精品一区二区三区在线| 精品少妇一区二区三区视频日本电影| 欧美人与性动交α欧美精品济南到| 免费看a级黄色片| 国产精品 欧美亚洲| 99在线人妻在线中文字幕 | 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 国产野战对白在线观看| 丰满人妻熟妇乱又伦精品不卡| av有码第一页| 水蜜桃什么品种好| 精品国产乱码久久久久久小说| 天天躁狠狠躁夜夜躁狠狠躁| 高清在线国产一区| 亚洲中文av在线|