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

    FCC提升管反應(yīng)器中終止劑注入對裂化反應(yīng)的影響

    2016-09-18 09:56:08王敏吳迎亞藍(lán)興英高金森
    化工學(xué)報 2016年8期
    關(guān)鍵詞:催化裂化管內(nèi)油氣

    王敏,吳迎亞,藍(lán)興英,高金森

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

    FCC提升管反應(yīng)器中終止劑注入對裂化反應(yīng)的影響

    王敏,吳迎亞,藍(lán)興英,高金森

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

    通過對催化裂化提升管注入終止劑前后的工況進(jìn)行數(shù)值模擬,研究了終止劑注入對提升管內(nèi)速度分布、催化劑顆粒濃度分布、溫度分布以及組分濃度分布的影響,考察了不同注入量以及注入高度的終止劑在提升管內(nèi)的作用區(qū)域及其對裂化反應(yīng)的影響。研究表明,終止劑的注入大幅提升了提升管內(nèi)的油氣速度,降低了催化劑濃度、油氣和催化劑的溫度,使得提升管內(nèi)原料的裂化程度降低,二次反應(yīng)減少。且不同注入量和注入高度的作用區(qū)域不同,對裂化反應(yīng)的影響不同,應(yīng)根據(jù)實(shí)際工況進(jìn)行分析。

    提升管;催化裂化;終止劑;計算流體力學(xué)

    引 言

    催化裂化在煉油工業(yè)中占據(jù)著至關(guān)重要的地位,其全球日產(chǎn)量可達(dá)1500萬桶[1],約占全球成品油總產(chǎn)量的18%。在我國,催化裂化作為重質(zhì)油輕質(zhì)化的主要手段,其2014年催化裂化原油加工能力占總原油加工能力的28.21%[2]。作為催化裂化裝置的核心部分,提升管內(nèi)發(fā)生著復(fù)雜的平行-順序反應(yīng),其目標(biāo)產(chǎn)物汽油和柴油為中間產(chǎn)物,因此控制反應(yīng)深度尤為重要。目前隨著原料重質(zhì)化、劣質(zhì)化,為了保證提升管進(jìn)料段原料油能夠在短時間內(nèi)充分汽化,大多數(shù)煉廠采用提高油氣和催化劑的混合溫度的方法。然而,進(jìn)料段溫度的提高會導(dǎo)致整個提升管內(nèi)部溫度的上升,使得油氣過裂化,進(jìn)而導(dǎo)致目標(biāo)產(chǎn)物收率的降低。為獨(dú)立控制提升管進(jìn)料段油劑混合溫度以及中上部油氣反應(yīng)溫度,國外采用的技術(shù)有 IFP公司提出的 MTC混合控溫技術(shù)[3-4]、Kellogg公司設(shè)計的提升管急冷技術(shù)[5]等。目前我國采用的一種靈活、簡單、高效的方法為終止劑技術(shù),通過在提升管的某一部位注入冷卻介質(zhì),優(yōu)化提升管內(nèi)部的溫度分布,從而優(yōu)化產(chǎn)物分布。

    在我國80%以上的FCC裝置具有使用終止劑的能力或正在使用終止劑技術(shù)。影響終止劑注入效果的工藝參數(shù)主要有終止劑的注入位置、注入量以及終止劑的類型等[6]。對于這些參數(shù)的確定,高金森等[7-9]通過對采用了反應(yīng)終止劑技術(shù)的提升管反應(yīng)器進(jìn)行數(shù)值模擬研究,確定了相應(yīng)的終止劑注入位置,選取了終止劑類型和注入量。在理論分析及模擬計算的基礎(chǔ)上,董偉等[10]根據(jù)勝利石油化工總廠重油催化裂化裝置提升管反應(yīng)器的實(shí)際情況,選定輕污油及直餾汽油作為終止劑進(jìn)行了工業(yè)試驗(yàn),發(fā)現(xiàn)終止劑的使用提高了輕質(zhì)油的收率,降低了干氣、焦炭等二次產(chǎn)物產(chǎn)率。此外多家煉廠[11-20]也通過在裝置上探索,對比了終止劑的注入位置、注入量以及終止劑的類型等對提升管內(nèi)溫度分布、各產(chǎn)物產(chǎn)率等方面的影響。然而,大部分對于終止劑的探索的關(guān)注點(diǎn)都在出口處的溫度以及產(chǎn)物濃度,實(shí)際上,不同注入條件下終止劑的流動區(qū)域不同,終止劑與油氣和催化劑進(jìn)行質(zhì)量傳遞、動量傳遞、熱量傳遞的作用范圍也不盡相同,進(jìn)而對反應(yīng)過程產(chǎn)生了不同影響,最終導(dǎo)致了出口處溫度以及產(chǎn)物濃度的變化。因此本文通過數(shù)值模擬,研究提升管終止劑的不同注入條件對系統(tǒng)內(nèi)流動、傳熱以及反應(yīng)行為的影響,揭示終止劑的作用機(jī)理,為其實(shí)際工業(yè)應(yīng)用提供理論依據(jù)。

    1 模擬方法

    1.1模擬對象

    以某煉廠催化裂化裝置提升管反應(yīng)器為模擬對象,取預(yù)提升段5 m,提升管直管段33.12 m,出口橫管段7.176 m。具體的操作條件見表1。終止劑選用120℃的水蒸氣,注入量分別為原料的0%、5%、10%,終止劑入口高度分別選取距原料油噴嘴15、24 m。

    表1 模擬條件Table 1 Simulation conditions

    1.2提升管綜合數(shù)學(xué)模型

    提升管內(nèi)低溫油氣和高溫催化劑接觸瞬間汽化,氣液固三相湍流流動,流動、傳熱、反應(yīng)高度耦合[21]。目前針對提升管的模型主要為耦合 CFD流動模型與集總動力學(xué)模型的綜合數(shù)學(xué)模型。Theologos等[22-25]最早將三維CFD流動模型與集總動力學(xué)模型結(jié)合來模擬提升管內(nèi)流動與反應(yīng)過程,但未考慮各相湍流流動特性。Gao等引入 k-ε-kp湍流模型并結(jié)合描述我國重油催化裂化的十三集總動力學(xué)模型對工業(yè)提升管進(jìn)行模擬[26],之后考慮原料油霧滴汽化,建立了提升管三維三相流動反應(yīng)耦合模型[27]。對于提升管內(nèi)的流動,目前多采用歐拉-歐拉雙流體模型,氣相與固相具有各自的流動特性且相間發(fā)生著動量、質(zhì)量、熱量等傳遞過程。對于反應(yīng),主要采用集總動力學(xué)模型。Weekman等首先提出集總理論,開發(fā)了三集總動力學(xué)模型[28],之后對原料油集總進(jìn)行進(jìn)一步劃分,提出催化裂化十集總反應(yīng)動力學(xué)模型[29]。洛陽石化工程公司和華東理工大學(xué)[30]結(jié)合我國催化裂化特點(diǎn),建立了催化裂化十一集總動力學(xué)模型。本文流動模型采用歐拉-歐拉雙流體模型,裂化反應(yīng)動力學(xué)模型為渣油十四集總動力學(xué)模型[31]。由于提升管反應(yīng)器體系新增終止劑組分,在原有提升管模型[31]的基礎(chǔ)上增加終止劑組分方程,見表2,其中i為終止劑組分。

    表2 終止劑組分方程Table 2 Component equation of term inator

    1.3模擬方法

    根據(jù)工業(yè)裝置的實(shí)際情況,構(gòu)建該提升管反應(yīng)器的三維幾何模型,采用ICEM軟件對該幾何模型進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分,并對噴嘴入口和壁面附近的網(wǎng)格進(jìn)行加密處理,網(wǎng)格數(shù)量總共約為50萬個。之后以FLUENT軟件為平臺,結(jié)合某煉廠裝置所測數(shù)據(jù),對模型進(jìn)行驗(yàn)證,表3為工業(yè)裝置所測數(shù)據(jù)與模擬值對比情況,由表3可見模型誤差較小。之后采用驗(yàn)證后的模型對提升管內(nèi)注入終止劑前后的流動、傳熱以及反應(yīng)情況進(jìn)行考察。

    表3 工業(yè)裝置數(shù)據(jù)與模擬數(shù)據(jù)對比情況Table 3 Com parisons between industrial data and simulation data

    2 結(jié)果與討論

    2.1終止劑注入對裂化反應(yīng)的影響

    終止劑注入后,提升管反應(yīng)器內(nèi)流動狀態(tài)發(fā)生明顯變化。圖1為終止劑注入前后提升管內(nèi)氣相速度矢量的變化情況,由圖可見,終止劑注入后,其入口附近油氣速度大幅提升,部分油氣向管中心區(qū)域流動。該油氣速度的變化一方面受注入終止劑氣流的影響,另一方面由于低溫終止劑與高溫油氣接觸換熱后升溫汽化而引起。終止劑的注入對油氣速度矢量的影響程度與終止劑注入量、入射角度、終止劑類型等有關(guān)。

    提升管內(nèi)顆粒體積分?jǐn)?shù)分布特點(diǎn)為邊壁大中心小,中心區(qū)域與邊壁區(qū)域的顆粒體積分?jǐn)?shù)相差較大。提升管反應(yīng)器從其邊壁處注入終止劑后,必然會對其邊壁處的顆粒體積分?jǐn)?shù)分布產(chǎn)生影響。注入終止劑后提升管內(nèi)不同截面上催化劑濃度分布如圖2所示,在終止劑流過的區(qū)域,催化劑濃度明顯降低,而終止劑未流過的區(qū)域,催化劑濃度變化不大。這是由于終止劑的注入,使得該區(qū)域氣相體積增加,相對而言催化劑體積分?jǐn)?shù)就降低。沿提升管高度方向,隨著終止劑的徑向擴(kuò)散作用,催化濃度差異減小。

    由于終止劑的汽化吸熱和升溫作用,終止劑的注入會引起提升管反應(yīng)器內(nèi)溫度的變化。圖 3為注入終止劑后提升管各截面油氣溫度的變化情況。由圖可看出,在終止劑流過的區(qū)域,溫度明顯降低,而終止劑未流過的區(qū)域,溫度變化不大,沿著提升管高度方向,油氣溫度逐漸分布均勻。這是由于與油氣相比終止劑的溫度較低,終止劑流經(jīng)的區(qū)域,終止劑所占比例大,而且還來不及與油氣進(jìn)行熱量交換,使得該區(qū)域溫度較低;沿提升管高度方向,終止劑不斷與油氣進(jìn)行熱量交換,溫度逐漸升高。

    圖1 終止劑注入前后氣相速度矢量圖Fig.1 Gas velocity vector w ith injection of term inator

    圖2 終止劑注入后提升管內(nèi)催化劑顆粒體積分?jǐn)?shù)分布Fig.2 Contour of solid volume fraction w ith injection of terminator

    圖3 終止劑注入后提升管內(nèi)油氣溫度分布Fig.3 Contour of temperature in riser with injection of terminator

    注入終止劑后,提升管內(nèi)各組分的濃度受到的影響有:終止劑的注入造成流動不均勻性,使得各組分濃度在不同區(qū)域產(chǎn)生不同變化;終止劑的注入影響了停留時間、反應(yīng)組分濃度、催化劑濃度、反應(yīng)溫度,進(jìn)而對裂化反應(yīng)造成影響,使得反應(yīng)生成的各產(chǎn)物的量發(fā)生改變。因此終止劑的注入對組分濃度的影響為多因素共同作用的結(jié)果。

    表4 終止劑注入后提升管內(nèi)產(chǎn)物收率分布Table 4 Change in yield of products w ith injection of term inator

    表4為終止劑注入前后提升管反應(yīng)器出口處產(chǎn)物的收率。與不注入終止劑相比,終止劑注入后,提升管內(nèi)原料的裂化程度降低,二次反應(yīng)也相應(yīng)減少,汽油收率提高,干氣、液化氣和焦炭的收率都有所降低??梢?,所注入的終止劑的確起到了終止部分二次裂化反應(yīng)的作用,且終止劑對二次反應(yīng)的終止作用隨終止劑的注入量的增加而增大,隨注入高度的升高而減小。

    圖4 終止劑注入后汽油濃度分布Fig.4 Contour of mass fraction of gasoline w ith injection of term inator

    為進(jìn)一步考察終止劑注入對組分濃度的影響,圖4給出了終止劑注入后提升管內(nèi)汽油濃度分布。終止劑注入后,在終止劑流過的區(qū)域,汽油濃度明顯降低;而終止劑未流過的區(qū)域,汽油濃度變化不大。沿提升管高度方向,由于汽油組分的擴(kuò)散作用以及終止劑對二次反應(yīng)的抑制,汽油濃度逐漸增加。

    2.2終止劑注入量對裂化反應(yīng)的影響

    終止劑的注入使得一定范圍內(nèi)的分子數(shù)增加,流速明顯提高,這將縮短油氣在提升管的反應(yīng)時間,自然也減少了二次裂化反應(yīng)的發(fā)生,有利于產(chǎn)物分布。然而不同注入量下終止劑的主要作用區(qū)域不同,對提升管內(nèi)不同位置的流速造成的影響不同,進(jìn)而對各處反應(yīng)的影響也不相同。

    圖5為不同終止劑注入量下的終止劑濃度分布,由圖可以看出,當(dāng)終止劑注入量為 5%時,由于終止劑注入量較小,射流沖擊作用小,終止劑主要貼著提升管管壁流動,因此邊壁處的油氣速度大幅增加,停留時間大幅縮短。當(dāng)終止劑注入量為10%時,終止劑注入量較大,射流沖擊作用強(qiáng),大部分終止劑向管中心區(qū)域流動,因此中心處的油氣速度大幅增加,停留時間大幅縮短。

    終止劑在沿高度方向流動的同時,在徑向逐漸向中心區(qū)域擴(kuò)散,不同注入量的終止劑擴(kuò)散情況也不盡相同。表5統(tǒng)計了注入終止劑后提升管不同位置處終止劑的質(zhì)量分?jǐn)?shù),選取位置如圖6所示。當(dāng)終止劑注入量為5%時,注入10 m后,終止劑質(zhì)量分?jǐn)?shù)在噴嘴上方邊壁附近B處為15.7%,而在中心區(qū)域A處為2.5%,邊壁C處僅0.5%,說明注入10 m后,終止劑逐漸由邊壁B處擴(kuò)散到了中心區(qū)域,但還未擴(kuò)散至邊壁C處。終止劑注入15 m后,終止劑逐漸擴(kuò)散至邊壁C處,但質(zhì)量分?jǐn)?shù)僅為0.7%,遠(yuǎn)低于預(yù)期。當(dāng)終止劑注入量為10%時,注入5 m后,終止劑質(zhì)量分?jǐn)?shù)在中心區(qū)域A處為11.5%,在噴嘴上方邊壁附近B處為3.5%,邊壁C處僅0.5%,說明注入5 m后,終止劑逐漸由中心區(qū)域擴(kuò)散到了邊壁B處,但還未擴(kuò)散至邊壁C處。終止劑注入10 m后,終止劑逐漸擴(kuò)散至邊壁C處,質(zhì)量分?jǐn)?shù)為2.4%,低于邊壁B處。終止劑注入15 m后,邊壁B處和C處的終止劑濃度基本相等。

    表5 不同終止劑注入量下提升管內(nèi)不同位置處終止劑濃度分布Table 5 M ass fraction of term inator in different position w ith different mass flux of term inator

    綜合以上分析,終止劑注入量較少時,終止劑的主要作用區(qū)域?yàn)閲娮焐戏竭叡谔帲⑷胍欢尉嚯x后終止劑才擴(kuò)散到中心區(qū)域且濃度始終低于噴嘴上方邊壁處,除了噴嘴上方附近區(qū)域的邊壁處外,其他邊壁處終止劑濃度始終偏低,終止劑起不到應(yīng)有的作用。終止劑注入量較多時,終止劑的主要作用區(qū)域?yàn)橹行膮^(qū)域,注入后逐漸擴(kuò)散至其他區(qū)域,各邊壁處的終止劑濃度基本相等,但仍略小于中心區(qū)域。

    圖5 不同終止劑注入量下的終止劑濃度分布Fig.5 Mass fraction of term inator w ith different mass flux of term inator

    圖6 選取位置分布Fig. 6 Distribution of different position in cross section

    為考察不同終止劑注入量對催化劑濃度分布的影響,對終止劑噴嘴上方的區(qū)域的催化劑濃度分布進(jìn)行重點(diǎn)分析,如圖7所示。不同終止劑注入量下催化劑濃度降低的區(qū)域不同,終止劑注入量較大時,催化劑濃度降低的范圍主要在提升管中心區(qū)域,而終止劑注入量較小時,催化劑濃度降低的范圍主要為提升管管壁區(qū)域。當(dāng)終止劑注入量為10%時,由于射流作用較大,噴嘴附近出現(xiàn)負(fù)壓區(qū),使得附近的催化劑往負(fù)壓區(qū)流動,出現(xiàn)了圖中所示的在噴嘴正上方的區(qū)域催化劑濃度反而高的情況。

    圖7 不同終止劑注入量下催化劑濃度分布Fig.7 Contour of solid volume fraction in riser w ith different mass flux of term inator

    對比不同終止劑注入量下的提升管內(nèi)溫度分布情況,從圖8可以看出,隨著終止劑注入量的增加,終止劑流過的區(qū)域面積增加,使得提升管內(nèi)溫度低的區(qū)域增加。沿提升管高度方向,隨著終止劑不斷與油氣發(fā)生熱量交換,溫度降低的幅度逐漸減弱。當(dāng)終止劑注入量較少時,終止劑主要是與噴嘴上方的管壁附近的油氣發(fā)生熱量交換,該區(qū)域附近溫度大幅降低;當(dāng)終止劑注入量較大時,終止劑進(jìn)入了提升管中心,終止劑主要是與管中心的油氣發(fā)生熱量交換,中心區(qū)域溫度大幅降低。注入終止劑的主要目的就是降低油氣溫度,減少二次裂化反應(yīng)發(fā)生。終止劑注入量影響管內(nèi)油氣溫度分布,從而也將影響二次反應(yīng)的發(fā)生情況,最終影響產(chǎn)物組成情況。

    圖8 不同終止劑注入量下提升管不同高度截面溫度分布Fig.8 Contour of temperature in riser w ith different mass flux of terminator

    為了便于對比分析終止劑注入量對組分濃度的影響,對終止劑噴嘴上方區(qū)域的汽油濃度分布情況進(jìn)行重點(diǎn)分析,如圖9所示。對比不同終止劑注入量下提升管內(nèi)組分濃度分布情況,當(dāng)終止劑注入量較少時,噴嘴上方管壁附近組分濃度明顯下降,而管中心濃度變化不大。當(dāng)終止劑注入量較大時,終止劑進(jìn)入了提升管中心,管中心的組分濃度明顯下降。隨著終止劑濃度的增加,該區(qū)域內(nèi)一次反應(yīng)產(chǎn)物汽油、柴油收率提高,二次反應(yīng)產(chǎn)物干氣、液化氣、焦炭收率降低。

    在提升管中,邊壁處催化劑濃度比中心處高,二次反應(yīng)速率也比中心處高,注入終止劑時,為抑制二次反應(yīng),理想的終止劑分布應(yīng)為邊壁處高于中心處。然而,當(dāng)終止劑注入量較高時,雖能有效擴(kuò)散至各區(qū)域,但中心區(qū)域濃度大于邊壁處,造成了終止劑的浪費(fèi);當(dāng)終止劑注入量較低時,雖然邊壁處終止劑濃度大于中心處,但由于目前大多數(shù)提升管裝置的終止劑入口只有兩個,終止劑無法擴(kuò)散至全區(qū)域,遠(yuǎn)離噴嘴上方的邊壁區(qū)域終止劑濃度過低,無法形成理想的終止劑濃度分布。因此建議增加終止劑噴嘴,即可在終止劑注入量較小的條件下增加終止劑在邊壁處的作用區(qū)域,從而更為有效地減少二次反應(yīng)的發(fā)生。

    圖9 不同終止劑注入量下汽油濃度分布Fig. 9 Contour of mass fraction of gasoline w ith different mass flux of term inator

    2.3終止劑注入高度對裂化反應(yīng)的影響

    對于終止劑注入高度對終止劑作用區(qū)域的影響,從圖10不同終止劑注入高度下的終止劑濃度分布可看出,其主要影響的是終止劑在軸向上的作用范圍。由圖 10(a)可以看出,當(dāng)終止劑注入高度為15 m時,終止劑注入約17 m進(jìn)入出口橫管段。而當(dāng)終止劑注入高度為24 m時,終止劑注入約8 m就進(jìn)入出口橫管段。

    圖10 不同終止劑注入高度下的終止劑濃度分布Fig.10 Mass fraction of term inator w ith different teem ing height

    表6統(tǒng)計了終止劑注入量為10%時,不同注入高度下提升管不同位置處終止劑的質(zhì)量分?jǐn)?shù)。當(dāng)終止劑注入高度為24 m時,出口橫管段前的截面處,邊壁B處終止劑濃度為5.9%,低于注入高度15 m的提升管的該處濃度,而邊壁C處終止劑濃度只有0.8%,說明終止劑在非噴嘴上方的邊壁處始終未達(dá)到足夠濃度。終止劑注入后,需要空間進(jìn)行充分?jǐn)U散,因此注入位置的升高會導(dǎo)致在高度方向沒有足夠的區(qū)域使得終止劑充分?jǐn)U散,進(jìn)而增加終止劑對速度影響的不均勻性。

    表6 不同終止劑注入高度下提升管內(nèi)不同位置處終止劑濃度分布Table 6 Mass fraction of term inator in different position w ith different teem ing height

    注入高度對終止劑徑向作用范圍的影響較弱,由圖10(b)可看出,對于不同注入高度的提升管,當(dāng)終止劑注入8 m后,橫截面上的終止劑濃度分布區(qū)別不大。注入高度為24 m的提升管終止劑濃度高的區(qū)域分布略為分散,這是由于在提升管中隨著高度的增加,裂化生成的氣體體積增加,油氣氣速增大,終止劑注入時射流沖擊作用變?nèi)酢?/p>

    圖11為終止劑注入量10%時不同注入高度下提升管內(nèi)溫度的分布情況。從圖中可看出終止劑注入后,提升管內(nèi)溫度急劇降低,且不同終止劑注入高度使得提升管內(nèi)溫度降低的范圍不同,終止劑注入高度越高,提升管內(nèi)溫度降低的范圍越小,終止劑未能與油氣和催化劑進(jìn)行充分的熱量交換,使得提升管內(nèi)終止劑入口以上區(qū)域溫度分布不均勻,進(jìn)而影響到反應(yīng)過程。

    圖11 不同注入高度下提升管溫度變化Fig.11 Temperature in riser w ith different teem ing height

    終止劑的注入高度對組分濃度的影響見表4,隨著注入位置的升高,汽油收率不斷降低,而干氣、液化氣和焦炭收率不斷增加。提高終止劑的注入高度,高溫條件下油氣與催化劑的接觸時間增加,因而終止裂化反應(yīng)的效果逐漸減弱,此時原料油的轉(zhuǎn)化率逐漸增加,汽油收率降低,干氣、液化氣及焦炭的收率則逐漸增加。可見,終止劑注入位置的選擇十分重要,應(yīng)在保證提升管反應(yīng)深度的同時盡可能減少發(fā)生二次裂化反應(yīng)。隨著終止劑注入量的增加,對裂化反應(yīng)的終止效果越明顯,因此汽油收率不斷增加,而干氣、液化氣和焦炭收率降低。

    3 結(jié) 論

    通過對某廠催化裂化裝置的提升管反應(yīng)器終止劑技術(shù)實(shí)施方案進(jìn)行模擬研究,分析了終止劑的注入對提升管內(nèi)速度分布、催化劑顆粒濃度分布、溫度分布以及組分濃度的影響,分析了不同終止劑注入量和注入高度的不同作用區(qū)域及其影響,得到如下結(jié)論。

    (1)終止劑的注入使得作用區(qū)域內(nèi)入口附近油氣速度大幅提升,油氣停留時間縮短,催化劑濃度有所降低。低溫終止劑與油氣和催化劑接觸換熱使得作用區(qū)域內(nèi)溫度降低。受以上變化的綜合影響,終止劑的注入減少了二次反應(yīng),提高了中間產(chǎn)物汽油、柴油的產(chǎn)率。

    (2)不同終止劑的注入量對于終止劑在徑向上的分布情況影響較大。當(dāng)終止劑注入量較少時,終止劑濃度在徑向上的分布為噴嘴上方的邊壁區(qū)域>中心區(qū)域>遠(yuǎn)離噴嘴上方的邊壁區(qū)域;注入量較多時,終止劑濃度在徑向上的分布為中心區(qū)域>噴嘴上方的邊壁區(qū)域≥遠(yuǎn)離噴嘴上方的邊壁區(qū)域。

    (3)終止劑注入高度對于終止劑在軸向上的分布情況影響較大,對終止劑徑向分布情況影響較小。隨著注入高度的增加,終止劑的作用范圍減小,對二次反應(yīng)的抑制作用減弱。

    (4)在提升管中注入終止劑時,為抑制二次反應(yīng),理想的終止劑分布應(yīng)為邊壁處高于中心處。然而目前大多數(shù)提升管裝置的終止劑入口只有兩個,無法形成理想的終止劑濃度分布,因此建議增加終止劑噴嘴,增加終止劑在邊壁處的作用區(qū)域,從而更為有效地減少二次反應(yīng)的發(fā)生。

    符號說明

    Pr ——湍流擴(kuò)散Prandtl數(shù)

    S ——源項

    Y ——質(zhì)量分?jǐn)?shù)

    Γ ——擴(kuò)散系數(shù),kg·m-2·s

    μ ——黏度,N·m-2·s

    Φ ——一般變量

    下角標(biāo)

    i ——組分

    g ——?dú)庀?/p>

    p ——顆粒相

    t ——時間

    References

    [1] VASILEIOS K, LYNNE X T, MELISSA C, et al. Zeolites in sustainable chemistry (FCC) [C]//XIAO F S, MENG X J. Zeolites in Sustainable Chemistry. Berlin: Springer Berlin Heidelberg, 2016: 271-297.

    [2] 金云, 朱和. 調(diào)整轉(zhuǎn)型升級 適應(yīng)新常態(tài) 引領(lǐng)新常態(tài)——中國煉油工業(yè)發(fā)展現(xiàn)狀及“十三五”發(fā)展趨勢[J]. 國際石油經(jīng)濟(jì), 2015, (5): 14-21.

    JIN Y, ZHU H. Development strategy and innovation in China’s petroleum and chemical industry during the 13th5-Year Plan [J]. International Petroleum Economics, 2015, (5): 14-21.

    [3] DEAN R R, MAULEON J L, LETZSCH W S. New resid cracker (part 1) [J]. Oil Gas Journal, 1982, 80(40): 75-80.

    [4] DEAN R R, MAULEON J L, LETZSCH W S. New resid cracker (part 2) [J]. Oil Gas Journal, 1982, 80(41): 168-176.

    [5] 鈕根林, 王新元. 渣油催化裂化工藝反應(yīng)技術(shù)新進(jìn)展(Ⅰ): 圍繞提升管反應(yīng)器的技術(shù)開發(fā)[J]. 石化技術(shù)與應(yīng)用, 2001, (3): 169-173.

    NIU G L, WANG X Y. Progress of residuum catalytic cracking reaction technology (Ⅰ ): Technology development around commercial riser reactor [J]. Petrochemical Technology & Application, 2001, (3): 169-173.

    [6] 毛安國, 常學(xué)良, 顧潔. 對催化裂化裝置使用終止劑技術(shù)的認(rèn)識[J]. 石化技術(shù)與應(yīng)用, 2003, (4): 276-278.

    MAO A G, CHANG X L, GU J. Analysis of applying quench technology in FCCU [J]. Petrochemical Technology & Application,2003, (4): 276-278.

    [7] 高金森, 徐春明, 林世雄, 等. 提升管反應(yīng)器氣固兩相流動反應(yīng)模型及數(shù)值模擬(Ⅲ): 注終止劑技術(shù)的數(shù)值模擬 [J]. 石油學(xué)報(石油加工), 1998, (3): 41-49.

    GAO J S, XU C M, LIN S X, et al. A gas-solid two phase flow-reaction model of FCC riser reactors and its numerical simulation (Ⅲ): Numerical modeling on the reaction-term inating technique [J]. Acta Petrolei Sinica (Petroleum Processing Section),1998, (3): 41-49.

    [8] 王艷花, 張紅梅, 高金森, 等. RFCC提升管注反應(yīng)終止劑方案的數(shù)值模擬[J]. 大慶石油學(xué)院學(xué)報, 1999, (1): 37-41.

    WANG Y H, ZHANG H M, GAO J S, et al. Numerical modeling on the reaction-terminating technique in commercial RFCC riser reactor [J]. Journal of Daqing Petroleum Institute, 1999, (1): 37-41.

    [9] 高金森, 鄭曉軍, 董偉, 等. 催化裂化提升管反應(yīng)器終止劑注入位置的確定[J]. 石油大學(xué)學(xué)報(自然科學(xué)版), 2000, (6): 1-3.

    GAO J S, ZHENG X J, DONG W, et al. Determination of location of terminator in FCC riser reactor [J]. Journal of the University of Petroleum, China (Edition of Natural Science), 2000, (6): 1-3.

    [10] 董偉, 徐春明, 高金森, 等. 催化裂化提升管反應(yīng)器終止劑技術(shù)的工業(yè)化試驗(yàn)[J]. 石油大學(xué)學(xué)報(自然科學(xué)版), 2000, (6): 4-6.

    DONG W, XU C M, GAO J S, et al. Industry experiment of reaction-terminating technique in FCC reactor [J]. Journal of the University of Petroleum, China (Edition of Natural Science), 2000,(6): 4-6.

    [11] 解亞平, 萬勝林. 提升管注終止劑技術(shù)總結(jié)[J]. 催化裂化, 1996,15(4): 9-11.

    XIE Y P, WAN S L. Summary of terminating technique in riser [J]. Catalytic Cracking, 1996, 15(4): 9-11.

    [12] 李洪, 馮羅明, 仲偉萍, 等. 終止劑技術(shù)在大慶重油催化裂化裝置上的應(yīng)用[J]. 石油煉制與化工, 2000, (11): 55-56.

    LI H, FENG L M, ZHONG W P, et al. Application of terminating technique in Daqing RFCC riser [J]. Petroleum Processing and Petrochemicals, 2000, (11): 55-56.

    [13] 張培平, 景春明, 魏俊文. 淺議反應(yīng)終止劑技術(shù)的應(yīng)用[J]. 浙江化工, 2004, (2): 20-21.

    ZHANG P P, JING C M, WEI J W. Application of reaction terminating technique [J]. Zhejiang Chemical Industry, 2004, (2): 20-21.

    [14] 曾毅. 重油催化裂化裝置反應(yīng)終止劑技術(shù)的應(yīng)用[J]. 河南石油,2004, (2): 66-67.

    ZENG Y. Application of reaction terminating technique in RFCC unit [J]. Henan Petroleum, 2004, (2): 66-67.

    [15] 張志芳, 韓劍敏. 提升管注終止劑技術(shù)的應(yīng)用與探討[J]. 煉油設(shè)計, 1994, (1): 6-9.

    ZHANG Z F, HAN J M. Application and discussion of terminating technique in riser [J]. Petroleum Refinery Engineering, 1994, (1): 6-9.

    [16] 李躍民, 王建軍, 云宏俊. FCC原料重質(zhì)化與終止劑技術(shù)的應(yīng)用[J].內(nèi)蒙古石油化工, 1997, (3): 98-101.

    LI Y M, WANG J J, YUN H J. FCC heavy feedstock and application of terminating technique [J]. Inner Mongolia Petrochemical Industry,1997, (3): 98-101.

    [17] 李躍民, 王建軍. 提升管注終止劑技術(shù)的應(yīng)用[J]. 催化裂化, 1998,17(5): 8-10.

    LI Y M, WANG J J. Application of terminating technique in riser [J]. Catalytic Cracking, 1998, 17(5): 8-10.

    [18] 邢穎春, 徐世泰. 提升管注反應(yīng)終止劑技術(shù)的應(yīng)用[J]. 催化裂化,1996, 15(2): 34-37.

    XING Y C, XU S T. Application of reaction terminating technique in riser [J]. Catalytic Cracking, 1996, 15(2): 34-37.

    [19] 莫國紅. 重油催化提高產(chǎn)品收率的工藝技術(shù)探索[J]. 杭州化工,2000, 30(3): 21-23.

    MO G H. Process to increase yield of RFCC product [J]. Hangzhou Chemical Industry, 2000, 30(3): 21-23.

    [20] 仲偉萍, 李洪. 終止劑技術(shù)在RFCCU上的應(yīng)用[J]. 黑龍江石油化工, 2001, 12(1): 8-11.

    ZHONG W P, LI H. Application of terminating technique in RFCCU [J]. Heilongjiang Petrochemical Technology, 2001, 12(1): 8-11.

    [21] 高金森, 徐春明, 楊光華, 等. 提升管反應(yīng)器氣固兩相流動反應(yīng)模型及數(shù)值模擬(Ⅰ): 氣固兩相流動反應(yīng)模型的建立[J]. 石油學(xué)報(石油加工), 1998, (1): 29-35.

    GAO J S, XU C M, YANG G H, et al. A gas-solid two-phase flow-reaction model of FCC riser reactors and numerical simulation for them (Ⅰ): Development of the gas-solid two-phase flow-reaction model [J]. Acta Petrolei Sinica (Petroleum Processing Section), 1998,(1): 29-35.

    [22] THEOLOGOS K N, MARKATOS N C. Advanced modeling of fluid catalytic cracking riser-type reactors [J]. AIChE Journal, 1993, 39: 1007-1017.

    [23] THEOLOGOS K N, NIKOU I D, LYHEROS A I, et al. Simulation and design of fluid catalytic-cracking riser-type reactors [J]. AIChE Journal, 1997, 43: 486-494.

    [24] THEOLOGOS K N, LYGEROS A I, MARKATOS N C. Feedstock atomization effects on FCC riser reactors selectivity [J]. Chemical Engineering Science, 1999, 54: 5617-5625.

    [25] THEOLOGOS K N, NIKOU D, LYGEROS A I. Simulation and design of fluid catalytic cracking riser type reactor [J]. Computers & Chemical Engineering, 1996, 43: 486-494.

    [26] GAO J, XU C, LIN S, et al. Advanced model for turbulent gas-solid flow and reaction in FCC riser reactors [J]. AIChE Journal, 1999, 45: 1095-113.

    [27] GAO J, LIN S, XU C, et al. Simulations of gas-liquid-solid 3-phase flow and reaction in FCC riser reactors [J]. AIChE Journal, 2001, 47: 677-692.

    [28] WEEKMAN JR V W. Model of catalytic cracking conversion in fixed, moving, and fluid bed reactors [J]. Industrial & Engineering Chemistry Process Design and Development, 1968, 7(1): 90-95.

    [29] GROSS B, JACOB S M, NACE D M, et al. Simulation of catalytic cracking process: US3960707 [P]. 1976-6-1.

    [30] 翁惠新, 歐陽福生, 馬軍. 重油催化裂化反應(yīng)集總動力學(xué)模型(Ⅰ):模型的建立[J]. 化工學(xué)報, 1995, 46(6): 662-668.

    WENG H X, OUYANG F S, MA J. Lumped model for heavy oil catalytic cracking reaction (Ⅰ): Establishment of the model [J]. Journal of Chemical Industry and Engineering (China), 1995, 46(6): 662-668.

    [31] LAN X Y, XU C M, WANG G, et al. CFD modeling of gas-solid flow and cracking reaction in two-stage riser FCC reactors [J]. Chemical Engineering Science, 2009, 64: 3847-3858.

    Effect of term inator injection on catalytic cracking reactions in FCC riser

    WANG M in, WU Yingya, LAN Xingying, GAO Jinsen
    (State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249, China)

    The flow, heat transfer and reaction behaviors in the fluid catalytic cracking (FCC) riser before and after the injection of terminator were simulated by the computational fluid dynam ics (CFD) method. Distributions of velocity, solid catalyst volume fraction, temperature, and species concentration in FCC riser were obtained after the injection of term inator. Besides, the action zone of the terminator and the effect on cracking reactions were analyzed w ith different mass fraction and teem ing height of the term inator. The simulation results showed that w ith injection of the terminator, the velocity of oil gas was increased but volume fraction and temperature of the catalyst as well as temperature of the oil gas were decreased. The change in flow and the heat transfer caused the reduction in degree of the cracking reactions and the secondary reaction of oil gas in the riser, which resulted in a higher yield of gasoline and lower yields of dry gas, LPG and coke. The action zone and the effect on cracking reaction varied w ith different mass fraction and teem ing height of the terminator, which should be evaluated by actual process conditions.

    riser; fluid catalytic cracking; term inator; computational fluid dynam ics

    date: 2016-04-07.

    LAN Xingying, lanxy@cup.edu.cn

    supported by the National Basic Research Program of China(2012CB215003).

    TQ 021.1

    A

    0438—1157(2016)08—3191—11

    10.11949/j.issn.0438-1157.20160438

    2016-04-07收到初稿,2016-06-02收到修改稿。

    聯(lián)系人:藍(lán)興英。第一作者:王敏(1992—),女,碩士研究生。

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃項目(2012CB215003)。

    猜你喜歡
    催化裂化管內(nèi)油氣
    平?jīng)?,油氣雙破2萬噸
    “峰中”提前 油氣轉(zhuǎn)舵
    《非常規(guī)油氣》第二屆青年編委征集通知
    上傾管內(nèi)油水兩相流流型實(shí)驗(yàn)研究
    提高催化裂化C4和C5/C6餾分價值的新工藝
    催化裂化裝置摻渣比改造后的運(yùn)行優(yōu)化
    油氣體制改革迷局
    能源(2017年5期)2017-07-06 09:25:55
    催化裂化汽油脫硫工藝淺析
    揚(yáng)礦管內(nèi)高速螺旋流的數(shù)值模擬與仿真
    催化裂化多產(chǎn)丙烯
    国产精品久久久久久精品电影| 亚洲一区二区三区不卡视频| 亚洲人与动物交配视频| 1024手机看黄色片| 国内精品久久久久久久电影| av片东京热男人的天堂| 91麻豆av在线| 亚洲第一电影网av| av黄色大香蕉| 久久精品国产综合久久久| 一级毛片精品| 老熟妇仑乱视频hdxx| 国产成年人精品一区二区| 制服丝袜大香蕉在线| 免费观看人在逋| 一本久久中文字幕| svipshipincom国产片| 亚洲成人久久爱视频| 色av中文字幕| 亚洲无线观看免费| 成人av一区二区三区在线看| 97超视频在线观看视频| 夜夜爽天天搞| 熟女少妇亚洲综合色aaa.| 99热这里只有是精品50| 一个人看的www免费观看视频| 熟女电影av网| 草草在线视频免费看| 女生性感内裤真人,穿戴方法视频| 最新中文字幕久久久久 | 色综合站精品国产| 中文字幕精品亚洲无线码一区| 久久久久久久久免费视频了| 国产99白浆流出| 免费电影在线观看免费观看| 成人特级av手机在线观看| 色综合亚洲欧美另类图片| 久久午夜综合久久蜜桃| 亚洲人成网站高清观看| 观看免费一级毛片| 国产在线精品亚洲第一网站| 少妇的逼水好多| 757午夜福利合集在线观看| 亚洲国产高清在线一区二区三| 国产精品久久久久久亚洲av鲁大| 色噜噜av男人的天堂激情| АⅤ资源中文在线天堂| 一个人免费在线观看电影 | 丁香六月欧美| www日本在线高清视频| 亚洲真实伦在线观看| 国产精品爽爽va在线观看网站| 国产成人aa在线观看| 脱女人内裤的视频| av天堂在线播放| 久久久国产精品麻豆| 国产成人一区二区三区免费视频网站| 在线观看舔阴道视频| 99国产极品粉嫩在线观看| 久久久久久久午夜电影| 少妇的丰满在线观看| 搡老熟女国产l中国老女人| 午夜福利在线在线| 露出奶头的视频| 99热6这里只有精品| 一区二区三区国产精品乱码| 久久人人精品亚洲av| 天天添夜夜摸| 成年女人毛片免费观看观看9| 亚洲自偷自拍图片 自拍| 少妇人妻一区二区三区视频| 欧美性猛交╳xxx乱大交人| 亚洲人成伊人成综合网2020| 国产亚洲精品一区二区www| 午夜免费成人在线视频| 久久中文字幕人妻熟女| 国产激情久久老熟女| 国产乱人伦免费视频| 国产蜜桃级精品一区二区三区| www.www免费av| 日韩欧美国产在线观看| 欧美一区二区精品小视频在线| 久久精品综合一区二区三区| 国产精品亚洲av一区麻豆| 亚洲专区字幕在线| 一区福利在线观看| 性色avwww在线观看| 久久人妻av系列| 国产午夜福利久久久久久| 天堂动漫精品| 国产日本99.免费观看| 久久精品亚洲精品国产色婷小说| 在线观看一区二区三区| 高清毛片免费观看视频网站| 国产精品电影一区二区三区| 午夜影院日韩av| 亚洲精品中文字幕一二三四区| 欧美成人免费av一区二区三区| www日本在线高清视频| 久久这里只有精品中国| 久久久色成人| 在线观看舔阴道视频| 国产精品综合久久久久久久免费| 又紧又爽又黄一区二区| 老熟妇仑乱视频hdxx| 国产真人三级小视频在线观看| 亚洲七黄色美女视频| 美女 人体艺术 gogo| 免费在线观看成人毛片| 99久久精品热视频| 日韩欧美三级三区| 久久久国产精品麻豆| 18美女黄网站色大片免费观看| 国产97色在线日韩免费| 国内精品久久久久精免费| 欧美日韩瑟瑟在线播放| 欧美+亚洲+日韩+国产| 国产伦精品一区二区三区视频9 | 久久久久久人人人人人| 久久天堂一区二区三区四区| av女优亚洲男人天堂 | 日本免费a在线| 欧美成人免费av一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 美女黄网站色视频| 99久久精品国产亚洲精品| 少妇熟女aⅴ在线视频| 色在线成人网| 男人的好看免费观看在线视频| 高潮久久久久久久久久久不卡| 亚洲av成人不卡在线观看播放网| 免费看美女性在线毛片视频| 日本在线视频免费播放| 大型黄色视频在线免费观看| 免费观看人在逋| 法律面前人人平等表现在哪些方面| 亚洲成人免费电影在线观看| a级毛片a级免费在线| 成人鲁丝片一二三区免费| 老熟妇乱子伦视频在线观看| 国产午夜精品久久久久久| 99久久精品国产亚洲精品| 亚洲国产精品合色在线| 天天躁狠狠躁夜夜躁狠狠躁| 又大又爽又粗| 国产1区2区3区精品| 国产三级在线视频| 嫁个100分男人电影在线观看| 成人欧美大片| avwww免费| 午夜免费激情av| 真人做人爱边吃奶动态| 精品国产三级普通话版| 午夜精品在线福利| 免费在线观看视频国产中文字幕亚洲| 久久久久久久久久黄片| 久久久国产成人精品二区| 手机成人av网站| 国产高清视频在线观看网站| 99久久综合精品五月天人人| 久久精品91蜜桃| 久久久国产精品麻豆| 色av中文字幕| 国内精品美女久久久久久| 国产精品一及| 久久久久久大精品| 91av网一区二区| 久久亚洲精品不卡| 亚洲精品国产精品久久久不卡| 精品久久久久久久末码| 成人性生交大片免费视频hd| www.熟女人妻精品国产| 免费在线观看视频国产中文字幕亚洲| 伊人久久大香线蕉亚洲五| 国产精品av视频在线免费观看| 日日干狠狠操夜夜爽| 老司机午夜十八禁免费视频| 窝窝影院91人妻| 精品国产亚洲在线| 国产精品av视频在线免费观看| 亚洲性夜色夜夜综合| 一级作爱视频免费观看| 中文字幕最新亚洲高清| 亚洲欧美精品综合一区二区三区| 91麻豆精品激情在线观看国产| 可以在线观看的亚洲视频| 国产午夜精品论理片| 久久婷婷人人爽人人干人人爱| 麻豆成人午夜福利视频| 国产精品99久久久久久久久| 哪里可以看免费的av片| 欧洲精品卡2卡3卡4卡5卡区| 国产av在哪里看| 日本 欧美在线| 成人欧美大片| 国产欧美日韩一区二区三| 成人性生交大片免费视频hd| 99视频精品全部免费 在线 | 搞女人的毛片| 国产欧美日韩精品亚洲av| 最新在线观看一区二区三区| 免费在线观看影片大全网站| 最近最新免费中文字幕在线| 国产成+人综合+亚洲专区| 色吧在线观看| 特大巨黑吊av在线直播| 亚洲 国产 在线| 免费人成视频x8x8入口观看| 小蜜桃在线观看免费完整版高清| 国产精品久久久人人做人人爽| 99久久国产精品久久久| 夜夜看夜夜爽夜夜摸| 极品教师在线免费播放| 啦啦啦韩国在线观看视频| 巨乳人妻的诱惑在线观看| 青草久久国产| 搡老熟女国产l中国老女人| 美女黄网站色视频| 国产成人啪精品午夜网站| 精品电影一区二区在线| 午夜福利在线观看吧| 看黄色毛片网站| 热99在线观看视频| 精品久久久久久久末码| 69av精品久久久久久| 国产午夜精品久久久久久| 真人一进一出gif抽搐免费| 国产精品亚洲av一区麻豆| 可以在线观看毛片的网站| 久久精品91蜜桃| 日韩欧美三级三区| 午夜福利视频1000在线观看| 男女那种视频在线观看| 国产精品久久久久久久电影 | 国产精品香港三级国产av潘金莲| 国产高清有码在线观看视频| 99热这里只有精品一区 | 国产高清视频在线观看网站| 最近在线观看免费完整版| 啦啦啦观看免费观看视频高清| 每晚都被弄得嗷嗷叫到高潮| 嫩草影院精品99| 亚洲avbb在线观看| 国产精品国产高清国产av| 亚洲中文字幕一区二区三区有码在线看 | 国产成人影院久久av| 国产成人欧美在线观看| 国产精品99久久久久久久久| 丁香欧美五月| 国产成人精品久久二区二区91| 欧美黄色淫秽网站| 天天一区二区日本电影三级| 男女下面进入的视频免费午夜| 中国美女看黄片| 日本一本二区三区精品| 色哟哟哟哟哟哟| 国产一区二区在线av高清观看| 午夜激情福利司机影院| 国产免费男女视频| 国产免费av片在线观看野外av| 久久九九热精品免费| 一级毛片高清免费大全| 2021天堂中文幕一二区在线观| 午夜福利在线观看吧| 久久精品亚洲精品国产色婷小说| 精品久久久久久久毛片微露脸| 久久人妻av系列| 国产不卡一卡二| 99久久国产精品久久久| 久久久久久九九精品二区国产| 嫩草影院入口| 国产毛片a区久久久久| 九九热线精品视视频播放| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 婷婷丁香在线五月| 少妇人妻一区二区三区视频| 日本黄色片子视频| 女人高潮潮喷娇喘18禁视频| 91麻豆精品激情在线观看国产| 成年人黄色毛片网站| 欧美绝顶高潮抽搐喷水| 亚洲狠狠婷婷综合久久图片| 国产1区2区3区精品| 黄片小视频在线播放| 欧美黄色淫秽网站| 国产野战对白在线观看| 久久久久性生活片| 精品国产亚洲在线| 男人的好看免费观看在线视频| 国产毛片a区久久久久| 俄罗斯特黄特色一大片| 韩国av一区二区三区四区| 国产欧美日韩精品亚洲av| 伦理电影免费视频| 热99在线观看视频| www.精华液| 丁香欧美五月| 亚洲第一欧美日韩一区二区三区| 精品久久久久久,| 久久性视频一级片| 欧美高清成人免费视频www| 亚洲熟女毛片儿| 91在线观看av| 丁香欧美五月| 亚洲七黄色美女视频| 女人被狂操c到高潮| 看片在线看免费视频| 夜夜看夜夜爽夜夜摸| 欧美黄色片欧美黄色片| 国产成人av教育| 成人高潮视频无遮挡免费网站| 久久婷婷人人爽人人干人人爱| 亚洲成av人片免费观看| 中文字幕高清在线视频| www日本黄色视频网| 午夜福利视频1000在线观看| 在线国产一区二区在线| 99re在线观看精品视频| 日韩高清综合在线| 国产精品电影一区二区三区| 少妇人妻一区二区三区视频| www.www免费av| 色综合婷婷激情| 久久婷婷人人爽人人干人人爱| av女优亚洲男人天堂 | 99在线人妻在线中文字幕| 国产人伦9x9x在线观看| 久久精品国产清高在天天线| 亚洲精品乱码久久久v下载方式 | 色精品久久人妻99蜜桃| 又紧又爽又黄一区二区| 十八禁人妻一区二区| av福利片在线观看| www.www免费av| 一二三四在线观看免费中文在| 狠狠狠狠99中文字幕| 国产熟女xx| 久久精品综合一区二区三区| 国产激情欧美一区二区| 麻豆国产97在线/欧美| 欧美一级a爱片免费观看看| 午夜a级毛片| 观看美女的网站| 日韩三级视频一区二区三区| 亚洲人与动物交配视频| 国产伦在线观看视频一区| 又紧又爽又黄一区二区| 国产精品一区二区三区四区久久| 欧美中文日本在线观看视频| 国产亚洲精品久久久久久毛片| 国产欧美日韩一区二区三| 国产97色在线日韩免费| 国产精品久久久久久久电影 | 一级a爱片免费观看的视频| 国产精品 国内视频| 国产精品野战在线观看| 日韩 欧美 亚洲 中文字幕| 嫁个100分男人电影在线观看| 亚洲午夜理论影院| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 老司机深夜福利视频在线观看| 宅男免费午夜| 老司机午夜十八禁免费视频| 黄色成人免费大全| 变态另类丝袜制服| 热99re8久久精品国产| av天堂中文字幕网| 偷拍熟女少妇极品色| 夜夜躁狠狠躁天天躁| 99久久精品热视频| e午夜精品久久久久久久| 白带黄色成豆腐渣| 色综合站精品国产| 午夜激情欧美在线| 69av精品久久久久久| 亚洲av电影不卡..在线观看| 精品国产亚洲在线| 久久中文字幕一级| 日韩欧美国产在线观看| 国产精品一区二区三区四区免费观看 | 国产伦人伦偷精品视频| 欧美大码av| 亚洲 国产 在线| 欧美日韩乱码在线| 欧美成人免费av一区二区三区| 欧美三级亚洲精品| 午夜亚洲福利在线播放| 亚洲av成人精品一区久久| 美女大奶头视频| 精品一区二区三区视频在线 | 欧美三级亚洲精品| 午夜亚洲福利在线播放| 日韩大尺度精品在线看网址| 男人舔女人下体高潮全视频| 午夜激情欧美在线| 成人特级黄色片久久久久久久| 我要搜黄色片| 窝窝影院91人妻| 亚洲激情在线av| 在线免费观看不下载黄p国产 | 亚洲精品色激情综合| 欧美中文日本在线观看视频| 成人三级黄色视频| 色尼玛亚洲综合影院| 国产私拍福利视频在线观看| 免费在线观看日本一区| 亚洲国产欧美人成| 欧美乱色亚洲激情| 色播亚洲综合网| 婷婷精品国产亚洲av| 久久精品91无色码中文字幕| 国产视频一区二区在线看| 青草久久国产| 91字幕亚洲| 成年免费大片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 伦理电影免费视频| 亚洲精品美女久久av网站| 窝窝影院91人妻| 国产一区二区三区在线臀色熟女| www日本在线高清视频| 嫩草影院入口| 999精品在线视频| 一级作爱视频免费观看| 久久人人精品亚洲av| 人妻久久中文字幕网| 国产一区二区在线观看日韩 | 人妻丰满熟妇av一区二区三区| 久久久水蜜桃国产精品网| 国产1区2区3区精品| 国产精品av久久久久免费| 一个人看视频在线观看www免费 | 伦理电影免费视频| 99在线人妻在线中文字幕| 亚洲av片天天在线观看| 亚洲七黄色美女视频| 一个人看的www免费观看视频| 俺也久久电影网| 国产高潮美女av| 欧美色欧美亚洲另类二区| 啦啦啦免费观看视频1| 欧美中文综合在线视频| 久久久精品大字幕| 亚洲精品色激情综合| 国产欧美日韩精品一区二区| 免费看美女性在线毛片视频| 亚洲欧美精品综合一区二区三区| 麻豆国产av国片精品| 婷婷精品国产亚洲av在线| 欧美中文日本在线观看视频| 丰满的人妻完整版| 美女被艹到高潮喷水动态| www.999成人在线观看| 最新美女视频免费是黄的| 网址你懂的国产日韩在线| 久久人人精品亚洲av| 成人av在线播放网站| 亚洲色图av天堂| 成在线人永久免费视频| 窝窝影院91人妻| 久久久久久久久久黄片| 久久伊人香网站| 丁香欧美五月| 久久久久性生活片| 身体一侧抽搐| 这个男人来自地球电影免费观看| 欧美日韩瑟瑟在线播放| 俄罗斯特黄特色一大片| 亚洲av中文字字幕乱码综合| 1024香蕉在线观看| 亚洲国产精品久久男人天堂| 欧美av亚洲av综合av国产av| 国产日本99.免费观看| 国产乱人视频| 亚洲欧美日韩高清专用| 国产视频一区二区在线看| 99国产精品一区二区蜜桃av| 不卡av一区二区三区| 国产一区二区三区在线臀色熟女| 小蜜桃在线观看免费完整版高清| 亚洲av电影在线进入| 国产熟女xx| 国内毛片毛片毛片毛片毛片| 国产熟女xx| 俄罗斯特黄特色一大片| 熟妇人妻久久中文字幕3abv| 亚洲av第一区精品v没综合| 国产亚洲精品av在线| 久久香蕉国产精品| 日日夜夜操网爽| 国产高清视频在线观看网站| 可以在线观看毛片的网站| 女生性感内裤真人,穿戴方法视频| 啦啦啦观看免费观看视频高清| 免费高清视频大片| 国产欧美日韩精品一区二区| 噜噜噜噜噜久久久久久91| 国产亚洲精品一区二区www| 亚洲av免费在线观看| 欧美成人性av电影在线观看| 特级一级黄色大片| 亚洲性夜色夜夜综合| 中文字幕av在线有码专区| 女生性感内裤真人,穿戴方法视频| 国产精品电影一区二区三区| 欧美性猛交黑人性爽| 国产精品野战在线观看| 老汉色av国产亚洲站长工具| 精品久久久久久久久久免费视频| www.999成人在线观看| 国产探花在线观看一区二区| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 久久精品综合一区二区三区| 亚洲无线观看免费| 最近最新中文字幕大全免费视频| 岛国在线观看网站| 国产伦人伦偷精品视频| 757午夜福利合集在线观看| 亚洲七黄色美女视频| 国产精品1区2区在线观看.| 日本一二三区视频观看| 亚洲成人久久性| 一进一出抽搐动态| 成人av一区二区三区在线看| 别揉我奶头~嗯~啊~动态视频| 两个人的视频大全免费| 精品无人区乱码1区二区| 人人妻人人澡欧美一区二区| 两人在一起打扑克的视频| 久久人人精品亚洲av| 国产v大片淫在线免费观看| 久久国产乱子伦精品免费另类| 成人鲁丝片一二三区免费| 日韩欧美精品v在线| 日韩欧美在线二视频| 国产私拍福利视频在线观看| 久久午夜综合久久蜜桃| 精品国产美女av久久久久小说| 国产在线精品亚洲第一网站| 亚洲精品粉嫩美女一区| 久久精品国产综合久久久| 一个人免费在线观看电影 | 亚洲熟妇中文字幕五十中出| 怎么达到女性高潮| 国产成年人精品一区二区| 免费观看精品视频网站| 亚洲av成人av| 丁香六月欧美| 中文字幕熟女人妻在线| 亚洲成av人片在线播放无| 亚洲狠狠婷婷综合久久图片| 日韩精品中文字幕看吧| 国产精品 国内视频| 亚洲欧美精品综合一区二区三区| 亚洲精品久久国产高清桃花| 一区二区三区激情视频| tocl精华| 欧美三级亚洲精品| 99精品久久久久人妻精品| 国产精品综合久久久久久久免费| 一进一出抽搐动态| 2021天堂中文幕一二区在线观| 特大巨黑吊av在线直播| 国产伦精品一区二区三区视频9 | 很黄的视频免费| 亚洲成人中文字幕在线播放| 熟女少妇亚洲综合色aaa.| 99精品欧美一区二区三区四区| 人人妻人人看人人澡| 日本免费a在线| 国产亚洲精品久久久久久毛片| 给我免费播放毛片高清在线观看| 老熟妇仑乱视频hdxx| 国产69精品久久久久777片 | 色综合欧美亚洲国产小说| avwww免费| 丰满人妻熟妇乱又伦精品不卡| 精品乱码久久久久久99久播| 国产激情欧美一区二区| 91av网站免费观看| 亚洲成人免费电影在线观看| 美女 人体艺术 gogo| 人人妻人人澡欧美一区二区| 午夜福利成人在线免费观看| 1024香蕉在线观看| 脱女人内裤的视频| 亚洲av免费在线观看| 天堂av国产一区二区熟女人妻| 色综合站精品国产| 国产精品乱码一区二三区的特点| 色播亚洲综合网| 欧美三级亚洲精品| 一本一本综合久久| 亚洲av成人不卡在线观看播放网| 欧美极品一区二区三区四区| 久久精品国产清高在天天线| 淫妇啪啪啪对白视频| 欧美一区二区国产精品久久精品| 黑人操中国人逼视频| 色精品久久人妻99蜜桃| 老汉色∧v一级毛片| 久久精品国产99精品国产亚洲性色| 久久婷婷人人爽人人干人人爱| 精品一区二区三区av网在线观看| 欧美性猛交╳xxx乱大交人| 在线十欧美十亚洲十日本专区| 99国产精品一区二区蜜桃av| 变态另类丝袜制服| 一级毛片精品| 男女之事视频高清在线观看| 999久久久精品免费观看国产| e午夜精品久久久久久久|