• <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ī)械控制
    国产精品一区二区在线观看99 | 热99re8久久精品国产| 久久人人爽人人片av| 桃色一区二区三区在线观看| 国产亚洲欧美98| 夜夜夜夜夜久久久久| 精品一区二区三区视频在线| 久久久久网色| 欧美色欧美亚洲另类二区| 成人无遮挡网站| 欧美成人a在线观看| 女人十人毛片免费观看3o分钟| 国产黄色小视频在线观看| 一个人观看的视频www高清免费观看| 性插视频无遮挡在线免费观看| 毛片一级片免费看久久久久| 欧美成人一区二区免费高清观看| 欧美性猛交黑人性爽| 哪里可以看免费的av片| 色综合站精品国产| 午夜免费激情av| 日本在线视频免费播放| 国产老妇女一区| 久久这里只有精品中国| 丝袜美腿在线中文| 亚洲人成网站在线播| 国产免费男女视频| 26uuu在线亚洲综合色| 青青草视频在线视频观看| 国产黄色视频一区二区在线观看 | 国产精品乱码一区二三区的特点| 人人妻人人澡欧美一区二区| 亚洲精品国产成人久久av| 免费观看的影片在线观看| 日本在线视频免费播放| 夜夜看夜夜爽夜夜摸| 久久99热这里只有精品18| а√天堂www在线а√下载| 亚洲一区高清亚洲精品| 91狼人影院| 国产极品天堂在线| 18禁在线播放成人免费| 欧美3d第一页| 免费无遮挡裸体视频| 久久精品国产自在天天线| 一夜夜www| 国产精品伦人一区二区| 亚洲最大成人中文| 国产精品久久久久久精品电影| av卡一久久| 晚上一个人看的免费电影| 不卡视频在线观看欧美| 久久这里只有精品中国| 久久九九热精品免费| 男女那种视频在线观看| 在线观看av片永久免费下载| 欧美一区二区亚洲| 99热这里只有精品一区| 高清午夜精品一区二区三区 | 欧美潮喷喷水| 2021天堂中文幕一二区在线观| 最近视频中文字幕2019在线8| 精品无人区乱码1区二区| 精品久久久久久久人妻蜜臀av| 国产精品综合久久久久久久免费| 听说在线观看完整版免费高清| 最近的中文字幕免费完整| 成人午夜精彩视频在线观看| 国产午夜精品久久久久久一区二区三区| 成人漫画全彩无遮挡| 久久久久性生活片| 老熟妇乱子伦视频在线观看| 色吧在线观看| 男人的好看免费观看在线视频| 丝袜喷水一区| 国内少妇人妻偷人精品xxx网站| 久久精品综合一区二区三区| 一级av片app| 天天一区二区日本电影三级| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲国产色片| 久久久久免费精品人妻一区二区| 看非洲黑人一级黄片| 一区二区三区四区激情视频 | 国产成人精品婷婷| 欧美色欧美亚洲另类二区| 欧美日韩国产亚洲二区| 久久精品夜夜夜夜夜久久蜜豆| av卡一久久| 久久久久久久久久久丰满| 在线观看66精品国产| 免费电影在线观看免费观看| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 色噜噜av男人的天堂激情| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| a级毛片免费高清观看在线播放| 非洲黑人性xxxx精品又粗又长| 日韩强制内射视频| 亚洲av电影不卡..在线观看| 欧美区成人在线视频| 国产精品,欧美在线| 一本一本综合久久| 少妇的逼水好多| 久久人人精品亚洲av| 嫩草影院新地址| 中文字幕免费在线视频6| 黄色日韩在线| 精品久久国产蜜桃| 国产精品久久久久久久久免| 夜夜看夜夜爽夜夜摸| 国产成人福利小说| 欧美3d第一页| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 亚洲久久久久久中文字幕| 成人亚洲欧美一区二区av| 亚洲精品色激情综合| 久久午夜亚洲精品久久| 尾随美女入室| 国产蜜桃级精品一区二区三区| 午夜福利视频1000在线观看| 中国美白少妇内射xxxbb| 哪个播放器可以免费观看大片| 99久久无色码亚洲精品果冻| 亚洲在线自拍视频| av在线观看视频网站免费| 亚洲av第一区精品v没综合| 精品99又大又爽又粗少妇毛片| 老熟妇乱子伦视频在线观看| 中文精品一卡2卡3卡4更新| 狂野欧美白嫩少妇大欣赏| 麻豆成人av视频| 午夜a级毛片| 久久九九热精品免费| 国产在线男女| 日韩制服骚丝袜av| 一级黄片播放器| 久久精品国产亚洲av涩爱 | 日韩一区二区视频免费看| 日产精品乱码卡一卡2卡三| 一区二区三区四区激情视频 | 噜噜噜噜噜久久久久久91| 免费看日本二区| 欧美日韩综合久久久久久| 波多野结衣高清作品| 中文字幕精品亚洲无线码一区| 一个人看视频在线观看www免费| 国产精品av视频在线免费观看| 国产精华一区二区三区| 国产精品99久久久久久久久| 亚洲精品色激情综合| 在线免费观看不下载黄p国产| 伊人久久精品亚洲午夜| 中文精品一卡2卡3卡4更新| 人妻系列 视频| 99久久九九国产精品国产免费| 丰满人妻一区二区三区视频av| 人妻少妇偷人精品九色| 伦精品一区二区三区| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 国产成人午夜福利电影在线观看| 亚洲综合色惰| 又爽又黄无遮挡网站| 免费av不卡在线播放| 久久精品影院6| 99热这里只有是精品在线观看| 久久久久国产网址| 精品久久久久久久久av| 成年av动漫网址| 永久网站在线| 99九九线精品视频在线观看视频| 国产乱人偷精品视频| 中文字幕av成人在线电影| 岛国在线免费视频观看| 成人一区二区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 晚上一个人看的免费电影| 插逼视频在线观看| 欧美一区二区亚洲| 免费观看精品视频网站| 日韩精品有码人妻一区| 久久午夜亚洲精品久久| 网址你懂的国产日韩在线| 亚洲熟妇中文字幕五十中出| 九九爱精品视频在线观看| 边亲边吃奶的免费视频| 亚洲高清免费不卡视频| 伦精品一区二区三区| 免费电影在线观看免费观看| 人妻久久中文字幕网| 我的老师免费观看完整版| 校园人妻丝袜中文字幕| 日韩一区二区视频免费看| 看非洲黑人一级黄片| 欧美色视频一区免费| 伦精品一区二区三区| 在线免费观看的www视频| 少妇熟女欧美另类| 在线天堂最新版资源| 99九九线精品视频在线观看视频| 精品国内亚洲2022精品成人| 美女国产视频在线观看| 国产单亲对白刺激| 在线观看av片永久免费下载| 亚洲欧美清纯卡通| 大香蕉久久网| 99九九线精品视频在线观看视频| 男女那种视频在线观看| 国产午夜精品久久久久久一区二区三区| 一本久久精品| 伦理电影大哥的女人| 有码 亚洲区| 亚洲人成网站在线播| 22中文网久久字幕| 天天躁日日操中文字幕| 久久久久久久久久成人| 久久久久久大精品| 在线播放国产精品三级| 精品久久久久久久久亚洲| 麻豆国产av国片精品| 欧美激情久久久久久爽电影| 变态另类丝袜制服| av在线老鸭窝| 久久精品国产自在天天线| 国产精品美女特级片免费视频播放器| 全区人妻精品视频| 日韩av在线大香蕉| 麻豆成人av视频| 欧美精品国产亚洲| 天堂影院成人在线观看| 婷婷色综合大香蕉| 夜夜夜夜夜久久久久| 免费观看在线日韩| 一区二区三区免费毛片| 亚洲最大成人中文| 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 中文欧美无线码| 99久国产av精品国产电影| 久久久久国产网址| 欧美日韩乱码在线| 国产亚洲精品av在线| 国产一区亚洲一区在线观看| 日本在线视频免费播放| 欧美bdsm另类| 嘟嘟电影网在线观看| 国产高潮美女av| 一区二区三区免费毛片| 欧美+日韩+精品| 国产国拍精品亚洲av在线观看| 欧美高清性xxxxhd video| av视频在线观看入口| 国产av麻豆久久久久久久| 男女做爰动态图高潮gif福利片| 99热这里只有是精品在线观看| 不卡一级毛片| 春色校园在线视频观看| 免费av不卡在线播放| 超碰av人人做人人爽久久| 精品久久久噜噜| 日韩精品有码人妻一区| 特大巨黑吊av在线直播| 国产一级毛片七仙女欲春2| 亚洲七黄色美女视频| 看黄色毛片网站| 免费不卡的大黄色大毛片视频在线观看 | 成人特级黄色片久久久久久久| 嫩草影院新地址| 亚洲人成网站在线观看播放| 97超视频在线观看视频| 国产乱人偷精品视频| 国产 一区精品| 亚洲国产精品成人久久小说 | 午夜激情欧美在线| 一进一出抽搐gif免费好疼| 一区二区三区免费毛片| 亚洲经典国产精华液单| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 狂野欧美激情性xxxx在线观看| 久久久久久久亚洲中文字幕| 久久久久久久久大av| 国产精品蜜桃在线观看 | 日韩亚洲欧美综合| 高清午夜精品一区二区三区 | 99久久久亚洲精品蜜臀av| 好男人视频免费观看在线| 国产精品99久久久久久久久| 一级毛片久久久久久久久女| 亚洲无线观看免费| 一级二级三级毛片免费看| 黄色欧美视频在线观看| 国产成年人精品一区二区| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线 | 搡女人真爽免费视频火全软件| 国产69精品久久久久777片| 天天躁夜夜躁狠狠久久av| 亚洲最大成人手机在线| 日本-黄色视频高清免费观看| 少妇猛男粗大的猛烈进出视频 | 插阴视频在线观看视频| 成人亚洲欧美一区二区av| 日日摸夜夜添夜夜爱| 中国美白少妇内射xxxbb| 亚洲国产欧美人成| 精品久久久久久久久久免费视频| 亚州av有码| 在线播放国产精品三级| 国产精品久久久久久精品电影| 欧美成人精品欧美一级黄| 国产av在哪里看| 乱系列少妇在线播放| 国内精品美女久久久久久| 亚洲精品久久国产高清桃花| 欧美区成人在线视频| 毛片一级片免费看久久久久| 久久鲁丝午夜福利片| 美女xxoo啪啪120秒动态图| 少妇的逼水好多| 精品久久久久久久久av| 国产伦精品一区二区三区四那| 99久久中文字幕三级久久日本| 久久人人爽人人片av| av福利片在线观看| 亚洲av.av天堂| 国产精品久久视频播放| 老女人水多毛片| 日韩人妻高清精品专区| 91久久精品国产一区二区三区| 亚洲一区高清亚洲精品| 国产毛片a区久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久午夜福利片| 精品久久久久久久人妻蜜臀av| 中文欧美无线码| 在线观看午夜福利视频| 亚洲丝袜综合中文字幕| 五月伊人婷婷丁香| 日本一二三区视频观看| 国产伦精品一区二区三区四那| 哪个播放器可以免费观看大片| 国产亚洲精品久久久com| 大香蕉久久网| 国产精品久久视频播放| 在线天堂最新版资源| 最后的刺客免费高清国语| 变态另类成人亚洲欧美熟女| 日韩欧美精品免费久久| 久久精品91蜜桃| 一本久久精品| 国产精品久久久久久久电影| 免费观看的影片在线观看| 夜夜看夜夜爽夜夜摸| 亚洲欧美清纯卡通| 麻豆久久精品国产亚洲av| 男人和女人高潮做爰伦理| 国产成人精品久久久久久| 男人的好看免费观看在线视频| 亚洲人成网站在线播放欧美日韩| 日本欧美国产在线视频| 三级国产精品欧美在线观看| 一进一出抽搐动态| 欧美日本亚洲视频在线播放| 久久99热6这里只有精品| 最好的美女福利视频网| 国产v大片淫在线免费观看| 亚洲精品日韩av片在线观看| 免费黄网站久久成人精品| 午夜a级毛片| 日本三级黄在线观看| 日韩制服骚丝袜av| 又黄又爽又刺激的免费视频.| 丰满人妻一区二区三区视频av| 国产高清三级在线| 长腿黑丝高跟| 精品久久久久久久人妻蜜臀av| 亚洲无线在线观看| 国产免费男女视频| 毛片一级片免费看久久久久| 高清午夜精品一区二区三区 | 欧美日本亚洲视频在线播放| 如何舔出高潮| 国产伦精品一区二区三区视频9| 亚洲精品乱码久久久v下载方式| 精品国产三级普通话版| 最近最新中文字幕大全电影3| 黑人高潮一二区| 联通29元200g的流量卡| 噜噜噜噜噜久久久久久91| 色5月婷婷丁香| 日本av手机在线免费观看| 国产精品乱码一区二三区的特点| 久久6这里有精品| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 国内精品美女久久久久久| 男人和女人高潮做爰伦理| 边亲边吃奶的免费视频| 人妻系列 视频| 少妇人妻一区二区三区视频| 国产精品久久电影中文字幕| 美女国产视频在线观看| 免费观看精品视频网站| or卡值多少钱| 九草在线视频观看| 国产成人影院久久av| 国内精品美女久久久久久| 人妻少妇偷人精品九色| 久久人人精品亚洲av| 精品日产1卡2卡| 男女下面进入的视频免费午夜| 亚洲第一电影网av| 国产一区亚洲一区在线观看| 小蜜桃在线观看免费完整版高清| 亚洲无线在线观看| 亚洲在线观看片| 亚洲欧美成人综合另类久久久 | 国产亚洲精品久久久com| 久久久久免费精品人妻一区二区| 亚洲一级一片aⅴ在线观看| 蜜桃亚洲精品一区二区三区| 91久久精品电影网| 一区二区三区免费毛片| 国产免费男女视频| 久久久久久大精品| 亚洲欧美日韩东京热| 久久久精品94久久精品| 久久人妻av系列| 在线观看一区二区三区| 99热这里只有是精品在线观看| 男的添女的下面高潮视频| 能在线免费看毛片的网站| 晚上一个人看的免费电影| av天堂中文字幕网| 免费人成视频x8x8入口观看| 日本撒尿小便嘘嘘汇集6| 国产成人影院久久av| 在线观看66精品国产| 欧美成人一区二区免费高清观看| 久久精品国产亚洲av涩爱 | 日本黄色片子视频| 免费观看在线日韩| 国产成人影院久久av| 夫妻性生交免费视频一级片| 日本av手机在线免费观看| 中国美白少妇内射xxxbb| 成年女人永久免费观看视频| 国产日韩欧美在线精品| 精品一区二区免费观看| 亚洲久久久久久中文字幕| 久久精品国产自在天天线| 国产精品伦人一区二区| 午夜福利在线在线| 99热全是精品| 国内精品美女久久久久久| 美女xxoo啪啪120秒动态图| 午夜激情欧美在线| 日韩人妻高清精品专区| 最近手机中文字幕大全| 人体艺术视频欧美日本| av卡一久久| 色综合色国产| а√天堂www在线а√下载| 亚洲无线观看免费| 乱人视频在线观看| 欧美人与善性xxx| 午夜亚洲福利在线播放| 久久久久国产网址| 欧美激情在线99| 网址你懂的国产日韩在线| av卡一久久| 人体艺术视频欧美日本| 亚洲精品国产av成人精品| 亚洲欧美日韩东京热| 中国国产av一级| 国产午夜精品一二区理论片| 成人特级av手机在线观看| 午夜精品一区二区三区免费看| 亚洲精品色激情综合| 黄色欧美视频在线观看| 美女国产视频在线观看| 久久精品国产亚洲av天美| 午夜a级毛片| 国产精品久久久久久久久免| 国产成人精品久久久久久| 亚洲精品色激情综合| 日本在线视频免费播放| 久久欧美精品欧美久久欧美| 久久久久网色| 亚洲人成网站在线播放欧美日韩| videossex国产| 少妇熟女aⅴ在线视频| 日韩av在线大香蕉| 亚洲欧美日韩东京热| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 最近视频中文字幕2019在线8| 国产在线男女| 欧美一区二区国产精品久久精品| 亚洲欧美日韩高清在线视频| 日韩欧美一区二区三区在线观看| 午夜精品国产一区二区电影 | 我的女老师完整版在线观看| 一本精品99久久精品77| 久久99热这里只有精品18| 91精品国产九色| 少妇丰满av| 国产91av在线免费观看| 久久这里只有精品中国| 免费看日本二区| av国产免费在线观看| 人妻久久中文字幕网| 国产不卡一卡二| 亚洲精品456在线播放app| 久久国产乱子免费精品| 日本黄色片子视频| 91aial.com中文字幕在线观看| 免费看日本二区| 哪个播放器可以免费观看大片| 成人二区视频| 日日啪夜夜撸| 性欧美人与动物交配| 人人妻人人澡欧美一区二区| 日韩一区二区三区影片| 亚洲精品久久久久久婷婷小说 | 女人十人毛片免费观看3o分钟| 大香蕉久久网| 亚洲高清免费不卡视频| 波野结衣二区三区在线| 亚洲欧美清纯卡通| 欧美三级亚洲精品| 性插视频无遮挡在线免费观看| 女同久久另类99精品国产91| 国产成人影院久久av| 成人毛片60女人毛片免费| 能在线免费观看的黄片| 欧美+亚洲+日韩+国产| 人妻少妇偷人精品九色| 亚洲精品亚洲一区二区| 非洲黑人性xxxx精品又粗又长| 国产三级中文精品| 成人二区视频| 欧美日本亚洲视频在线播放| 最近2019中文字幕mv第一页| 寂寞人妻少妇视频99o| 六月丁香七月| 看免费成人av毛片| 91久久精品国产一区二区三区| 国产色婷婷99| 日韩三级伦理在线观看| 国产精品麻豆人妻色哟哟久久 | 国产激情偷乱视频一区二区| 狂野欧美激情性xxxx在线观看| 国产精品人妻久久久久久| a级毛片a级免费在线| 国产伦一二天堂av在线观看| 中国国产av一级| 国产精品99久久久久久久久| av视频在线观看入口| 亚洲色图av天堂| 亚洲精品久久国产高清桃花| 天堂中文最新版在线下载 | 麻豆成人午夜福利视频| 人妻夜夜爽99麻豆av| 精品无人区乱码1区二区| 亚洲精品自拍成人| 久久人人爽人人片av| 五月伊人婷婷丁香| 亚洲天堂国产精品一区在线| 菩萨蛮人人尽说江南好唐韦庄 | 国产真实乱freesex| 欧美精品一区二区大全| 日韩亚洲欧美综合| 亚洲精品乱码久久久久久按摩| 国产69精品久久久久777片| 18禁黄网站禁片免费观看直播| 99热全是精品| 91aial.com中文字幕在线观看| 亚洲第一电影网av| 国产黄a三级三级三级人| 久久久久久久久久成人| 午夜福利在线观看吧| 成人毛片a级毛片在线播放| 久久精品久久久久久久性| 热99在线观看视频| 99热精品在线国产| 亚洲精品成人久久久久久| 久久久a久久爽久久v久久| 欧美性猛交黑人性爽| 大香蕉久久网| 老司机影院成人| 黄片无遮挡物在线观看| 又爽又黄a免费视频| 久久精品综合一区二区三区| 人妻制服诱惑在线中文字幕| 99久久精品国产国产毛片| 国产精品久久久久久亚洲av鲁大| 亚洲av不卡在线观看| av又黄又爽大尺度在线免费看 | 免费一级毛片在线播放高清视频| 在线观看美女被高潮喷水网站| 亚洲最大成人中文| 成人无遮挡网站| 12—13女人毛片做爰片一| 国产亚洲欧美98| 村上凉子中文字幕在线| 亚洲久久久久久中文字幕| 国产一区亚洲一区在线观看| 日韩在线高清观看一区二区三区| 免费人成视频x8x8入口观看|