朱麗云，趙文斌，仲理科，張玉春，王振波，金有海

（1中國石油大學(xué)（華東）重質(zhì)油國家重點實驗室，山東 青島 266580；2山東理工大學(xué)農(nóng)業(yè)工程與食品科學(xué)學(xué)院，山東 淄博 255049）

催化裂化用短接觸旋流反應(yīng)器內(nèi)氣固滑移特性

朱麗云1，趙文斌1，仲理科1，張玉春2，王振波1，金有海1

（1中國石油大學(xué)（華東）重質(zhì)油國家重點實驗室，山東 青島 266580；2山東理工大學(xué)農(nóng)業(yè)工程與食品科學(xué)學(xué)院，山東 淄博 255049）

利用歐拉-歐拉雙流體模型對短接觸旋流反應(yīng)器分離腔內(nèi)氣固滑移特性進行了數(shù)值模擬，主要研究了切向氣固滑移速度的分布規(guī)律，并考察了操作參數(shù)和物性參數(shù)對分離腔內(nèi)切向滑移速度的影響。計算結(jié)果表明，分離腔內(nèi)切向氣固滑移速度沿徑向呈“駝峰”分布；當(dāng)氣相入口速度增大或者劑氣比減小時，切向氣固滑移速度變小，顆粒切向速度增大，離心力增大，有利于提高氣固分離效率；顆粒密度對切向滑移速度分布影響不大；顆粒粒徑較大時，在排塵口易出現(xiàn)堵塞，不利于長周期運行；建立了截面平均切向氣固滑移速度計算模型，計算值與模擬結(jié)果誤差在±7.0%以內(nèi)。

短接觸旋流反應(yīng)器；分離腔；滑移速度；模擬；分離；模型

引　言

以重質(zhì)油為原料生產(chǎn)汽、柴油兼產(chǎn)丙烯等化工原料的催化裂化工藝（FCC）是世界上最重要的重油輕質(zhì)化的工藝之一［1-2］。在原料劣質(zhì)化和環(huán)保要求嚴(yán)格化的推動下，中國石油大學(xué)（華東）提出了短接觸旋流反應(yīng)器［3］，由混合腔和分離腔組成。混合腔內(nèi)切向入口進氣與軸向入口進劑接觸混合反應(yīng)并進行預(yù)旋，經(jīng)導(dǎo)向葉片旋轉(zhuǎn)加速后進入分離腔形成三維剪切流場，催化劑在離心力和重力作用下沿壁面旋轉(zhuǎn)下行，產(chǎn)物則向上流出，完成快速分離，從而實現(xiàn)反應(yīng)與分離一體化過程?？梢姡眯鞣磻?yīng)器內(nèi)各相混合強度高、停留時間短、油氣與催化劑易快速分離等特點，可以有效改善系統(tǒng)設(shè)備結(jié)焦、過裂化等問題，同時可以提高反應(yīng)選擇性、改善產(chǎn)品分布和提高目標(biāo)產(chǎn)品收率。目前，該旋流反應(yīng)器處于基礎(chǔ)研究階段［4-12］，僅對反應(yīng)器內(nèi)氣相流場及氣固兩相流動開展了初步研究。

在氣-固非均勻體系中，氣固之間的速度存在一定的差異［13］，即滑移速度，它反映了氣固間相互作用，是影響氣固兩相的質(zhì)量傳遞、動量傳遞和化學(xué)反應(yīng)的重要因素。對催化裂化提升管反應(yīng)器［14-17］來說，器內(nèi)氣固滑移速度越大，顆粒返混越嚴(yán)重，氣固接觸效率越低，原油轉(zhuǎn)化深度和產(chǎn)品收率降低。下行床反應(yīng)器［18-19］內(nèi)滑移速度越大，顆粒聚集程度越大，同樣不利于反應(yīng)的進行。提升管和下行床內(nèi)均只在軸向上存在氣固滑移現(xiàn)象。

而短接觸旋流反應(yīng)器，由于器內(nèi)存在反應(yīng)與分離的耦合過程，滑移速度對旋流反應(yīng)器內(nèi)氣固反應(yīng)過程的影響較復(fù)雜。在混合腔內(nèi)，主要為軸向滑移速度，滑移速度越大，氣固接觸效率越低，越不利于油劑接觸反應(yīng)；在分離腔內(nèi)，氣固兩相流為三維剪切流場，徑向遷移和螺旋下行同時存在，分離效率不僅取決于軸向和徑向上氣固滑移速度的大小，更與氣固切向滑移速度密切相關(guān)，切向滑移速度越大，顆粒跟隨性越差，顆粒切向速度減小，所受離心力也減小，不利于顆粒沿徑向向邊壁移動，則徑向滑移速度變小，顆粒易被氣體攜帶向中心區(qū)流動，并向上流出分離腔，從而導(dǎo)致氣固分離效率降低。因此，對旋流反應(yīng)器分離腔內(nèi)切向滑移速度開展研究，有助于理解旋流反應(yīng)器的分離機制。

本文利用計算流體力學(xué)軟件（CFD）對旋流反應(yīng)器分離腔內(nèi)切向氣固滑移特性進行數(shù)值模擬，并考察操作參數(shù)和物性參數(shù)對滑移速度的影響，為短接觸旋流反應(yīng)器的優(yōu)化設(shè)計提供理論指導(dǎo)。

1　數(shù)值模擬方法

1.1模型的選擇

結(jié)合前期結(jié)果，多相流模型選用歐拉-歐拉雙流體模型［7］，對催化劑顆粒采用擬流體的處理方法，將其視為連續(xù)相介質(zhì)。湍流模型選用RNG k-ε湍流模型［4-6，8］。兩個模型涉及的連續(xù)方程、動量方程、能量方程、κ-ε運輸方程等見表1。方程中的各系數(shù)參見文獻［7］。

1.2幾何建模與網(wǎng)格劃分

模擬過程中所用的短接觸旋流反應(yīng)器如圖1所示，分離腔直徑為100 mm，混合腔直徑為250 mm，計算坐標(biāo)采用直角坐標(biāo)系，原點設(shè)在混合腔頂部中心，Z軸沿軸向向下，X軸與切向入口方向平行?？紤]到計算精度，網(wǎng)格劃分本著結(jié)構(gòu)化網(wǎng)格優(yōu)于非結(jié)構(gòu)化的原則［20］，采用分區(qū)劃分網(wǎng)格的方法，將短接觸旋流反應(yīng)器劃分為26萬個網(wǎng)格。

表1　模型中涉及的方程Table 1　Equations in models

圖1　短接觸旋流反應(yīng)器幾何模型Fig.1　Geometric model of short-contact cyclone reactor

1.3數(shù)值解法及邊界條件

求解方法采用非耦合隱式，時間依賴采用非穩(wěn)態(tài)模型［21］。在離散化格式中，壓力插值格式采用PRESTO算法，壓力速度耦合插值采用SIMPLE算法［22］。為了計算結(jié)果的真實性，采用非穩(wěn)態(tài)計3～6 s的數(shù)據(jù)平均值作比較。

1.3.1進口邊界采用速度入口，根據(jù)已知流量以及入口直徑，直接得到氣相入口速度并計算出其他相應(yīng)的湍流參數(shù)，包括湍流強度和水力直徑。固相濃度和入口速度根據(jù)顆粒循環(huán)量計算獲得。具體操作參數(shù)：切向入口氣速ug1=2.5 m·s-1，軸向入口氣速ug2=5 m·s-1，催化劑入口速度us=1.2 m·s-1（對應(yīng)劑氣比為16），催化劑中位粒徑70 μm，催化劑顆粒密度為2000 kg·m-3。

1.3.2出口邊界出口按照湍流流動充分發(fā)展處理，采用自由出流邊界。

1.3.3固壁邊界壁面為無滑移邊界條件，壁面粗糙度為0.5，采用標(biāo)準(zhǔn)壁面函數(shù)法處理邊界湍流。其他參數(shù)設(shè)置見文獻［7］。

1.4切向氣固滑移速度

本文定義切向氣固滑移速度為氣相切向速度與固相切向速度之差，即

2　模擬結(jié)果與討論

2.1模擬結(jié)果可靠性驗證

圖2　模擬值與實驗值對比Fig.2　Comparison between experimental data and simulation results

氣固在旋流反應(yīng)器分離腔內(nèi)形成三維剪切流場，尚未有精確的測量儀器可以準(zhǔn)確測量出分離腔內(nèi)固含率和顆粒速度的分布情況，因此，通過比較混合腔內(nèi)z=200 mm截面上實驗測量的固含率分布與模擬結(jié)果進行對比，以此驗證模擬結(jié)果的可靠性。實驗裝置、測量方法等見文獻［7］。圖2為在 ug1= 2.5 m·s-1，ug2=5 m·s-1的條件下，實驗測量結(jié)果與模擬結(jié)果的對比。在z=200 mm截面上，固含率沿徑向呈現(xiàn)中心小、邊壁大的分布，且實驗值與模擬值吻合較好。因此，選用歐拉-歐拉雙流體模型來模擬短接觸旋流反應(yīng)器內(nèi)氣固兩相流動分布，并研究氣固兩相滑移速度分布情況是可行的。

2.2操作參數(shù)對切向氣固滑移速度的影響

2.2.1氣相入口速度（處理量）對切向氣固滑移速度的影響為考察氣相入口速度（處理量）對反應(yīng)器分離腔內(nèi)滑移速度分布的影響，分別設(shè)置：①ug1=1.5 m·s-1，ug2=3.0 m·s-1（Qg=53 m3·h-1）；②ug1=2.0 m·s-1，ug2=4.0 m·s-1（Qg=70 m3·h-1）；③ug1=2.5 m·s-1，ug2=5.0 m·s-1（Qg=88 m3·h-1）；④ug1=3.0 m·s-1，ug2=6.0 m·s-1（Qg=106 m3·h-1）。

圖3為不同氣相入口速度下旋流反應(yīng)器分離腔內(nèi)不同截面上切向滑移速度的分布曲線。由圖3（a）可知，在z=600 mm截面上，不同氣相入口速度下切向滑移速度沿徑向的分布與切向速度曲線變化規(guī)律基本一致，近似為“駝峰”形分布［8-9，23］，在中心區(qū)切向滑移速度出現(xiàn)最小值，說明顆粒在中心區(qū)跟隨性較好，有利于顆粒在氣體攜帶下獲得切向速度，并在離心力作用下向邊壁區(qū)聚集；切向滑移速度較大值出現(xiàn)在近邊壁區(qū)（r/R=±0.8），說明顆粒在該區(qū)域的跟隨性變差，分析其原因可能是在邊壁區(qū)顆粒由于邊壁效應(yīng)以及顆粒之間相互作用使得顆粒切向速度減小，增大了切向氣固滑移速度。隨著氣相入口速度增大，由ug1=2.5 m·s-1，ug2=5.0 m·s-1和ug1=3.0 m·s-1，ug2=6.0 m·s-1曲線可知，滑移速度在逐漸變小，數(shù)值在0～1.5 m·s-1之間，且大多數(shù)區(qū)域滑移速度小于0.5 m·s-1，說明隨著氣相入口速度的不斷增大，氣固之間速度越接近，顆粒的跟隨性越好，顆粒所受的離心力越大，越利于顆粒向邊壁區(qū)移動，從而更有利于氣固的分離，這是因為氣相入口速度增大時，氣體進入反應(yīng)器的能量增大，即加速和輸送顆粒的能量增大，顆粒跟隨氣體做切向運動。

圖3　氣相入口速度對切向氣固滑移速度分布的影響Fig.3　Effect of gas velocity at inlets on tangential slip velocity

如圖3（b）所示，在z=1000 mm截面，隨著氣相入口速度的增大，滑移速度不斷減小，當(dāng)氣相入口速度達到ug1=3.0 m·s-1，ug2=6.0 m·s-1時，滑移速度為零，說明氣相入口速度小時顆粒跟隨性差，而氣相入口速度大時顆粒獲得的輸送能量最大，顆粒跟隨性也最好，最有利于顆粒在氣體攜帶下做離心運動。當(dāng)氣相入口速度較小時，滑移速度方向不同，大小各異，曲線波動較大，甚至出現(xiàn)負(fù)值，分析其原因可能是z=1000 mm截面位于排塵段區(qū)域，由于結(jié)構(gòu)的變化使得上、下行氣流在此處速度較大，造成強烈的動量交換，使得該區(qū)域內(nèi)存在二次渦流，當(dāng)氣體氣相入口速度較小時，顆粒獲得的輸送能力較小，顆粒跟隨性變差，則顆粒容易受局部渦流區(qū)的影響，使其運動軌跡發(fā)生偏轉(zhuǎn)，即出現(xiàn)所謂的“龍擺尾”［24］，部分顆粒被卷入渦流區(qū)，并在中心上行流的夾帶下被走，從而不利于氣固的分離。

2.2.2劑氣比對切向氣固滑移速度的影響本文模擬是在冷態(tài)下進行，定義催化劑與空氣質(zhì)量比為劑氣比（catalyst to air，CTA）來對應(yīng)反應(yīng)條件下的劑油比。在保持氣相入口速度不變（88 m3·h-1）時，考察不同劑氣比（分別為16、22、26、30）對分離腔內(nèi)切向氣固滑移速度的影響，如圖4所示。

由圖4（a）可看出，在z=600 mm截面上，不同劑氣比下切向氣固滑移速度沿徑向分布趨勢變化不大，仍近似為“駝峰”形分布，但是隨著劑氣比的增大，切向滑移速度逐漸變大，這是因為劑氣比的增大導(dǎo)致旋流反應(yīng)器內(nèi)顆粒體積分?jǐn)?shù)增大，其所需的輸送能量也隨之增大，在氣相入口速度（ug1=2.5 m·s-1，ug2=5.0 m·s-1）不變的情況下，顆粒跟隨性變差。

由圖4（b）可知，在z=1000 mm截面滑移速度則沒有明顯變化規(guī)律，劑氣比最大時曲線波動最大，滑移速度最大值可達3.0 m·s-1，氣固相速度差增大，顆粒跟隨性變差，則顆粒所受的離心力會變小，不利于顆粒做離心運動，從而會影響到氣固的分離效率。然而，對于重油催化裂化來說，增大劑油比可以提高氣固接觸效率，從而有效提高產(chǎn)品的收率。另外，短接觸旋流反應(yīng)器內(nèi)催化裂化反應(yīng)是一個短時接觸的過程，需要在較大的劑油比下才能達到較高的產(chǎn)品收率。雖然增大劑氣比理論上不利于提高氣固分離效率，但是考慮到旋流反應(yīng)器是反應(yīng)分離一體化設(shè)備，較大的劑氣比會促進反應(yīng)器內(nèi)油劑接觸混合反應(yīng)，因此，需要綜合考慮各個因素選擇合適的劑氣比。

圖4　劑氣比對切向氣固滑移速度分布的影響Fig.4　Effect of CTA on tangential slip velocity

2.3物性參數(shù)對切向氣固滑移速度的影響

2.3.1顆粒粒徑對切向氣固滑移速度的影響一般地，相對小的顆粒通過曳力的作用會降低氣相湍流強度，而由于顆粒尾渦的作用，大顆粒的加入則有可能增強氣相的湍流強度。因此，分別采用平均粒徑為50、60和70 μm的催化劑顆粒進行數(shù)值計算，考察催化劑顆粒粒徑對旋流反應(yīng)器分離腔內(nèi)切向氣固滑移速度的影響規(guī)律，計算結(jié)果如圖5所示。

由圖5可知，在z=600 mm截面上，粒徑為50 和 60 μm顆粒切向氣固滑移速度（最大值為 4.0 m·s-1）大于70 μm顆粒（2.0 m·s-1左右），這是因為在離開導(dǎo)向葉片時，顆粒與氣相切向速度較接近，在旋轉(zhuǎn)向下運動時，相較于粒徑小的顆粒，粒徑大的顆粒因自身慣性力較大，受短路流、徑向氣流等的影響較小，顆粒切向速度變化小，切向氣固滑移速度要小，跟隨性好。由圖5（b）可知，與z=600 mm截面相比，在z=1000 mm截面處，粒徑為50 和 60 μm 的顆粒最大切向氣固滑移速度由 4.0 m·s-1左右降低到1.0 m·s-1左右，而70 μm顆粒切向滑移速度變化不大，且沿徑向呈現(xiàn)不規(guī)則分布，最大值出現(xiàn)在r/R≤±0.3附近，而在邊壁區(qū)切向滑移速度較小，且出現(xiàn)負(fù)值，說明此處產(chǎn)生了次級流動［25］，其擾流作用不容忽視。z=1000 mm截面處于分離腔與灰斗連接的重要部位，次級流動的產(chǎn)生會導(dǎo)致顆粒返混夾帶的發(fā)生，對于油氣與催化劑的分離效果以及催化劑顆粒的順利排出都會產(chǎn)生較大影響，進而影響整個循環(huán)系統(tǒng)的長周期運行。

圖5　顆粒粒徑對切向氣固滑移速度分布的影響Fig.5　Effect of dpon tangential slip velocity

2.3.2顆粒密度對切向氣固滑移速度的影響圖6為不同催化劑顆粒密度（1800、2000和2200 kg·m-3）對旋流反應(yīng)器內(nèi)切向滑移速度的影響規(guī)律。由圖 6可看出，在z=600 mm截面切向滑移速度近似呈對稱分布，r/R=-1～0三條曲線重合度很高，在r/R=0～1范圍內(nèi)則存在差別，ρp=1800， 2000kg·m-3沒有出現(xiàn)較為明顯的切向滑移速度最大值。在z=1000 mm截面，ρp=1800 kg·m-3時切向滑移速度在0.1 m·s-1以下，氣固兩相速度較接近，顆粒在氣體攜帶下做螺旋向下運動，利于顆粒從排塵口排出而實現(xiàn)氣固分離。圖中在 0.4＜r/R＜1范圍內(nèi)出現(xiàn)不對稱現(xiàn)象，分析原因是氣固兩相在分離腔內(nèi)主要以螺旋的流動形態(tài)向下流動，則在同一軸向截面處沿徑向方向兩側(cè)處在不同的螺旋帶上，并不是完全對稱的，因此出現(xiàn)了如圖所示的不對稱現(xiàn)象。

圖6　顆粒密度對切向氣固滑移速度分布的影響Fig.6　Effect of ρpon tangential slip velocity

2.4切向氣固滑移速度模型的建立

計算結(jié)果表明z=1000 mm位置處于旋流反應(yīng)器的排塵錐段，旋流強度已明顯減弱且由于存在各種次級流動，使其流動狀況變得更為隨機和復(fù)雜，規(guī)律性差，因此選取z=600 mm截面為特征截面，研究截面平均切向氣固滑移速度（vmean，slip）與各操作參數(shù)（氣相處理量 Qg，劑氣比 CTA，催化劑粒徑dp，顆粒密度ρp）的變化關(guān)系。根據(jù)局部切向氣固滑移速度計算得到截面平均切向氣固滑移速度，假設(shè)其與主要影響因素的關(guān)聯(lián)式為

通過將模擬數(shù)據(jù)進行回歸可得到

對比模擬值和理論模型計算值，如圖7所示。兩者間相對誤差為±7.0%。由關(guān)聯(lián)模型可知，對氣固滑移速度影響最大的是催化劑顆粒粒徑，后續(xù)將系統(tǒng)研究顆粒粒徑對氣固滑移速度的影響。

圖7　旋流反應(yīng)器vmean，slip模擬值與計算值的比較Fig 7　Comparison of simulated and calculated values of vmean，slip

3　結(jié)　論

（1）切向氣固滑移速度沿徑向呈“駝峰”形分布，在中心區(qū)顆粒跟隨性較好，而在近邊壁區(qū)由于顆粒團聚的影響氣固滑移速度增大。

（2）隨著氣相入口速度增大，顆粒加速和輸送顆粒的能量增大，切向氣固滑移速度減小，顆粒切向速度增大，離心力增大，有利于顆粒做離心運動，提高氣固分離效率。

（3）當(dāng)氣相入口速度不變時，降低劑氣比，切向氣固滑移速度變小，顆粒跟隨性變好，有利于顆粒在氣體攜帶下做離心運動而提高分離效率。但過低的劑氣比將引起產(chǎn)品收率的降低，因此要選擇合適的劑氣比，匹配催化裂化反應(yīng)過程和分離過程。

（4）根據(jù)不同操作參數(shù)和物性參數(shù)下切向氣固滑移速度分布規(guī)律，建立了截面平均切向滑移速度關(guān)聯(lián)模型，計算值與模擬值誤差為±7.0%。從模擬結(jié)果和關(guān)聯(lián)模型中可知，顆粒粒徑對切向滑移速度影響最大。

符號說明

CTA——劑氣比

c1，c2——根據(jù)經(jīng)驗確定的常數(shù)

dp——顆粒直徑，μm

g——重力加速度，m·s-2

pk——壓力，Pa

Qg——處理量，m3·h-1

r/R——量綱1半徑

u slip——滑移速度，m·s-1

ut——切向速度，m·s-1

vmean，slip——截面平均切向滑移速度，m·s-1

β ——相間動量轉(zhuǎn)換系數(shù)，kg·m-3·s-1

γs——能量碰撞耗散，kg·m-3·s-1

εg——氣含率

εp——固含率

θ ——顆粒溫度，m2·s-2

κs——脈動能傳導(dǎo)率，kg·m-3·s-1

ρ ——密度，kg·m-3

φ ——湍能耗散率，m2·s-3

References

［1］王剛， 高金森， 徐春明. 高溫短接觸催化裂化反應(yīng)技術(shù)的研究與進展 ［J］. 石油與天然氣化工， 2004， 33 （1）: 37-39.

WANG G， GAO J S， XU C M. Advances in the high-temperature and short-contact time technology of FCC ［J］. Chem ical Engineering of Oil & Gas， 2004， 33 （1）: 37-39.

［2］ZHU C， JUN Y， PATEL R， et al. Interactions of flow and reaction in fluid catalytic cracking risers ［J］. AIChE Journal， 2011， 57 （11）: 3122-3131.

［3］王振波， 王建， 任相軍， 等. 一種超短接觸旋流反應(yīng)-再生系統(tǒng): 201010615284.7 ［P］. 2015-01-07.

WANG Z B， WANG J， REN X J， et al. A kind of short contact cyclone reaction/regeneration system: 201010615284.7 ［P］. 2015-01-07.

［4］張玉春， 王振波， 金有海. 入口結(jié)構(gòu)對超短接觸旋流反應(yīng)器流場影響的數(shù)值模擬 ［J］. 中國石油大學(xué)學(xué)報（自然科學(xué)版）， 2012， 36 （5）: 169-174.

ZHANG Y C， WANG Z B， JIN Y H. Numerical simulation of effects of inlet structure on flow field in quick-contact cyclone reactor ［J］. Journal of China University of Petroleum （Edition of Natural Science），2012， 36 （5）: 169-174.

［5］張玉春， 王振波， 金有海， 等. 超短接觸旋流反應(yīng)器混合腔內(nèi)氣固混合特性的數(shù)值模擬 ［J］. 化工學(xué)報， 2012， 63 （2）: 424-432.

ZHANG Y C， WANG Z B， JIN Y H， et al. Numerical simulation on gas-solid m ixing characteristics in m ixing cavity of quick-contact cyclone reactor ［J］. CIESC Journal， 2012， 63 （2）: 424-432.

［6］張玉春， 王振波， 金有海， 等. 超短接觸旋流反應(yīng)器分離腔氣固分離特性的數(shù)值模擬 ［J］. 流體機械， 2012， 40 （2）: 24-29， 45.

ZHANG Y C， WANG Z B， JIN Y H， et al. Numerical simulation of gas-solid separation characteristics in separation space of quick-contact cyclone reactor ［J］. Fluid Machinery， 2012， 40 （2）: 24-29， 45.

［7］ZHANG Y C， WANG Z B， JIN Y H. Simulation and experiment of gas-solid flow field in short-contact cyclone reactors ［J］. Chem ical Engineering Research & Design， 2013， 91 （9）: 1768-1776.

［8］張玉春， 王振波， 金有海. 工藝參數(shù)對短接觸旋流反應(yīng)器內(nèi)顆粒流動特性影響的數(shù)值研究 ［J］. 石油學(xué)報（石油加工）， 2013， 29 （1）: 102-109.

ZHANG Y C， WANG Z B， JIN Y H. Numerical simulation of effect of operation parameters on particle flow characteristics in the quick-contact cyclone reactor ［J］. Acta Petrolei Sinica （Petroleum Processing Section）， 2013， 29 （1）: 102-109.

［9］黃峰， 張玉春， 王振波， 等. 短接觸旋流反應(yīng)器導(dǎo)葉位置對氣相流動的影響 ［J］. 化工學(xué)報， 2013， 64 （9）: 3153-3160.

HUANG F， ZHANG Y C， WANG Z B， et al. Effect of vane-guide location on gas flow behavior in short-contact cyclone reactor ［J］. CIESC Journal， 2013， 64 （9）: 3153-3160.

［10］ZHANG Y， WANG Z， JIN Y， et al. CFD simulation and experiment of residence time distribution in short-contact cyclone reactors ［J］. Advanced Powder Technology， 2015， 26 （4）: 1134-1142.

［11］王振波， 張玉春， 徐春明. 不同曳力模型及顆粒碰撞恢復(fù)系數(shù)對短接觸旋流反應(yīng)器內(nèi)氣固流場的影響 ［J］. 化工學(xué)報， 2014， 65 （6）: 2034-2041.

WANG Z B， ZHANG Y C， XU C M. Effect of drag model and restitution coefficient on gas-solids flow field in quick-contact cyclone reactors ［J］. CIESC Journal， 2014， 65 （6）: 2034-2041.

［12］黃峰， 孫治謙， 王振波， 等. 旋流反應(yīng)器內(nèi)顆粒停留時間分布的CFD模擬 ［J］. 石油化工設(shè)備， 2014， 43 （1）: 9-13.

HUANG F， SUN Z Q， WANG Z B， et al. CFD simulation of particle residence time distribution in cyclone reactor ［J］. Petro-Chem ical Equipment， 2014， 43 （1）: 9-13.

［13］白丁榮， 金涌， 俞芷青， 等. 快速流化床中平均滑落速度及絮狀物的特性 ［J］. 化學(xué)工程， 1989， 17 （6）: 44-49.

BAI D R， JIN Y， YU Z Q， et al. The average slip velocity and characteristics of clusters in fast fluidized bed ［J］. Chemical Engineering （China）， 1989， 17 （6）: 44-49.

［14］YANG Y L， JIN Y， YU Z Q， et al. Investigation on slip velocity distributions in the riser of dilute circulating fluidized bed ［J］. Powder Technology， 1992， 73 （1）: 67-73.

［15］楊勇林， 蔣大洲， 金涌， 等. 提升管循環(huán)流化床氣、固局部滑落速度的分布 ［J］. 高校化學(xué)工程學(xué)報， 1994， 8 （2）: 157-164.

YANG Y L， JIANG D Z， JIN Y， et al. Local gas-solid slip velocity in the riser of the circulating fluidized bed ［J］. Journal of Chem ical Engineering of Chinese Universities， 1994， 8 （2）: 157-164.

［16］李松庚， 林偉剛， 姚建中. 提升管中氣固相的局部滑移速度 ［J］.化工學(xué)報， 2003， 54 （1）: 119-123.

LI S G， LIN W G， YAO J Z. Local gas-solid slip velocity in riser ［J］. Journal of Chem ical Industry and Engineering （China）， 2003， 54 （1）: 119-123.

［17］CHIAN W C， SEVILLE J P K， PARKER D J， et al. Particle velocities and their residence time distribution in the riser of a CFB ［J］. Powder Technology， 2010， 203 （2）: 187-197.

［18］LI S， LIN W， YAO J. Local slip velocity in a downer ［J］. China Particuology， 2003， 1 （6）: 258-261.

［19］曹春社， 金涌， 俞芷青， 等. 并流下行循環(huán)流化床中氣固速度剖面以及滑落現(xiàn)象的研究 ［J］. 化工學(xué)報， 1995， 46 （2）: 180-187.

CAO C S， JIN Y， YU Z Q， et al. Gas-solids velocity profiles and slip phenomenon in a cocurrent downflow circulating fluidized bed ［J］. Journal of Chemical Industry and Engineering （China）， 1995， 46 （2）: 180-187.

［20］王振波， 馬藝， 金有海. 切流式旋流器內(nèi)兩相流場的模擬 ［J］. 中國石油大學(xué)學(xué)報（自然科學(xué)版）， 2010， 34 （4）: 136-140.

WANG Z B， MA Y， JIN Y H. Simulation of two-phase flow field in tangential hydrocyclone ［J］. Journal of China University of Petroleum （Edition of Natural Science）， 2010， 34 （4）: 136-140.

［21］倪建軍， 梁欽鋒， 代正華， 等. 撞擊流氣化爐內(nèi)氣粒兩相流動的數(shù)值模擬 ［J］. 化工學(xué)報， 2009， 60 （4）: 864-871.

NI J J， LIANG Q F， DAI Z H， et al. Numerical simulation of gas-particle two-phase flow in impinging streams gasifier ［J］. CIESCJournal， 2009， 60 （4）: 864-871.

［22］金有海， 姬廣勤， 曹晴云， 等. 旋風(fēng)分離器排氣管內(nèi)氣相流場的數(shù)值模擬 ［J］. 中國石油大學(xué)學(xué)報（自然科學(xué)版）， 2008， 32 （6）: 109-112.

JIN Y H， JI G Q， CAO Q Y， et al. Numerical simulation of gas-phase flow field in vortex finder of cyclone separators ［J］. Journal of China University of Petroleum （Edition of Natural Science）， 2008， 32 （6）: 109-112.

［23］CORTES C， GIL A. Modeling the gas and particle flow inside cyclone separators ［J］. Progress in Energy and Combustion Science， 2007， 33 （5）: 409-452.

［24］吳小林， 熊至宜， 姬忠禮， 等. 旋風(fēng)分離器旋進渦核的數(shù)值模擬［J］. 化工學(xué)報， 2007， 58 （2）: 383-390.

WU X L， XIONG Z Y， JI Z L， et al. Numerical simulation of precessing vortex core in cyclone separator ［J］. Journal of Chemical Industry and Engineering （China）， 2007， 58 （2）: 383-390.

［25］WANG B， XU D L， CHU K W， et al. Numerical study of gas-solid flow in a cyclone separator ［J］. Applied Mathematical Modelling，2006， 30 （11）: 1326-1342.

Characteristics of slip velocity between gas and solid in short-contact cyclone FCC reactor

ZHU Liyun1， ZHAO Wenbin1， ZHONG Like1， ZHANG Yuchun2， WANG Zhenbo1， JIN Youhai1
（1State Key Laboratory of Heavy Oil Processing， China University of Petroleum， Qingdao 266580， Shandong， China；2School of Agricultural Engineering and Food Science， Shandong University of Technology， Zibo 255049， Shandong， China）

The characteristics of slip velocity between gas and solid in a short-contact cyclone FCC reactor was simulated using the Eulerian-Eulerian two fluid model， especially about radial distribution of the tangential slip velocity between gas and solid as well as its dependence on operational and physical parameters in the separation chamber. Simulation results showed that the tangential slip velocity between gas and solid had a two-humped distribution in the radial direction. With increase of the inlet gas velocity or decrease of the ratio of catalyst to gas，the tangential slip velocity between gas and solid weakened but the tangential velocity of solid particles enhanced，which would strengthen the centrifugal force on particles and improve the gas-solid separation efficiency. The particle density exhibited little influence on radial distribution of the tangential slip velocity between gas and solid. However， when the particle diameters were increased， the dust outlet would easily be clogged which was not conducive to running the reactor for a long cycle of time. Finally， a model of the average cross-section tangential slip velocity between gas and solid was established w ith no more than ±7.0% difference among model simulations and calculated results.

short-contact cyclone reactor； separation chamber； slip velocity； simulation； separation； model

date: 2016-03-29.

Prof. WANG Zhenbo， wangzhb@upc.edu.cn

supported by the National Basic Research Program of China （2012CB215000） and the National Natural Science Foundation of China （21276281）.

TQ 051.8

A

0438—1157（2016）08—3468—08

10.11949/j.issn.0438-1157.20160369

2016-03-29收到初稿，2016-05-13收到修改稿。

聯(lián)系人：王振波。第一作者：朱麗云（1987—），女，博士后。

國家重點基礎(chǔ)研究發(fā)展計劃項目（2012CB215000）；國家自然科學(xué)基金項目（21276281）；中國博士后基金項目（2016M 592266）。