柳樹河油頁巖異步旋轉干燥技術

楊 寧， 周云龍， 齊天宇， 朱孝宇

(1.東北電力大學 能源與動力工程學院， 吉林 吉林 132012； 2.清華大學 核能與新能源技術研究院， 北京 100084)

柳樹河油頁巖異步旋轉干燥技術

楊 寧1， 周云龍1， 齊天宇2， 朱孝宇1

(1.東北電力大學 能源與動力工程學院， 吉林 吉林 132012； 2.清華大學 核能與新能源技術研究院， 北京 100084)

為充分挖掘流化床的最大干燥效力進而對油頁巖流化干燥過程進行強化，將常規布風板改造為異步旋轉布風裝置，調節內嵌式中心轉盤布風裝置與外套式行星轉環布風裝置的轉速，實現不同的布風工況。實驗研究了不同干燥工況對柳樹河油頁巖干燥特性的影響，并選取9種常用干燥模型對實驗數據進行干燥動力學分析。結果表明：適當地降低中心轉盤布風裝置的轉速，增加外套式行星轉環布風裝置的轉速，可以改善物料在流化床內的分布狀態，降低干燥所需的時間；物料的流化干燥過程主要分為升速干燥、恒速干燥和降速干燥3個階段，合理的調整布風板轉速會促使物料由升速干燥階段直接進入降速干燥階段；Two-term干燥模型對異步旋轉工況下物料的干燥實驗數據進行擬合的效果最好，但由于物料空隙內部吸附的氣體會隨水分一同擴散，且每次實驗結束后物料仍殘留少量結合水，導致模擬結果大于實驗結果。

流化床； 異步旋轉； 油頁巖； 干燥技術

油頁巖作為一種能源資源，既可以低溫干餾煉油，又可以燃燒發電，由于其特殊的組成和結構決定了它在能源、礦產、化工、環保方面具有廣闊的潛在用途[1-5]。油頁巖是一種多孔物質，在開采及儲運過程中，由于其內在或外在因素使得油頁巖含有一定的水分，若油頁巖含水率過高，直接干餾，爐中很大一部分熱量將用于油頁巖水分的蒸發，造成輸送氣體熱介質的管道發生形變，設備投資加大，同時大量的水蒸氣進入干餾爐中，加大了回收和油水分離系統的負荷，使油頁巖的回收率降低，因此，油頁巖在熱解干餾前必先經歷一個干燥脫水過程。

流化床干燥器作為一類常見的工業操作設備，在生產中的應用已十分廣泛[6-10]。然而傳統氣-固流化床存在床面波動大、顆粒橫向混合差等明顯缺點，造成了干燥效率降低等缺陷。近年來，大量學者致力于提升流化床干燥效率的相關研究，李凡等[11]通過在流化床內部加裝攪拌器，有效地防止了顆粒之間的黏附、結塊，以及顆粒在床層壁面的附著，明顯地改善了流化床內濕物料及易凝聚物料顆粒間的傳熱效果。宮國清等[12]改變傳統流化床恒定送風為周期性送風，通過調節氣流的脈沖頻率或脈沖氣流導通率，使通過孔板的氣體流量或流化區發生周期性變化，實現了物料的流化干燥。Kuipers等[13]將特定要求的振動源施加于常規流化床上的新型干燥裝置，裝置克服了溝流、騰涌的缺點，在很大程度上強化了傳熱效率和傳質速率。2006年，Madhiyanon等[14]設計研發旋轉流化布風裝置，將布風板中心位置與電機轉軸連接，通過調節電機的轉速間接帶動布風板隨之發生旋轉，在一定程度上極大地改善了流化床內物料的分布形態，至此強化流化床干燥能力實現階躍性突破。Sobrino等[15]進一步研究了流化床旋轉布風板轉速在0～100 r/min的范圍內變化時對物料干燥效果造成的影響，研究發現隨著旋轉轉速的增加，顆粒的橫向分布狀態趨于均勻，流化床的干燥能力也隨之增加。然而，布風板的旋轉速率控制在100 r/min是否是最佳的干燥流化工況，是否可以通過改變旋轉布風板局部風帽的旋轉速率進一步提升流化床的干燥能力，這一問題仍值得研究。筆者創造性地研發出一套異步轉速布風裝置，通過調整內嵌式中心轉盤布風板與外套式行星轉環布風板的轉速，對比異步旋轉布風裝置在不同布風工況下的干燥效率，尋求最佳的旋轉工況，采用9種常用干燥數學模型對油頁巖干燥實驗結果進行數學模擬，建立適用于異步旋轉流化狀態下的干燥方程，為油頁巖有效干燥提供重要參考依據。

1 實驗部分

1.1 實驗樣品

實驗所用樣品為吉林省柳樹河油頁巖，收到基含水率為31.99%，將油頁巖樣品破碎篩分至1～3 mm，置于恒溫箱內備用。樣品的工業分析及元素分析結果見表1。

表1 柳樹河油液巖工業分析及元素分析結果Table 1 Proximate and elemental analysis results of Liushuhe oil shale

M—Moisture; V—Volatile matter; A—Ash; FC—Fixed carbon; C—Carbon; H—Hydrogen; N—Nitrogen; S—Sulphur;Qnet—Net calorific value; ad—Air-dried basic

1.2 實驗設備及方法

1.2.1 實驗設備

空氣壓縮機(W-0.9/7C)，沈陽市東陵空壓機廠產品；空氣加熱器(WT-JQR)，鹽城市沃特機械設備有限公司產品；溫度傳感器(3144P)，美國艾默生有限公司產品；渦街流量計(LUGB-231106FA)，上海帆揚機電有限公司產品；變頻電機A(112M1-4)，上海通太電機有限公司產品；電容層析成像裝置(ERT)，英國工業層析成像公司產品；電子天平(AUW) ，日本島津公司產品。

1.2.2 實驗方法

油頁巖異步旋轉干燥實驗裝置如圖1所示。實驗在由樹脂玻璃加工而成的圓柱形截面流化床中進行，流化床高為1.5 m，直徑為0.21 m，由空氣壓縮機(0.9 m3/min)供給的空氣經過空氣加熱器(50～250℃)升溫后再經過閥門進入流化床風室，異步旋轉布風裝置被安裝在流化床的底部位置。流化床入口位置安裝有溫度傳感器(0～300℃)、渦街流量計(0～5 m3/min)，可以實現氣流溫度及速率的實時監測，本實驗中選取干燥氣流溫度為(200±0.5)℃，流化氣流速率為1.6 m/s，所采用的工況參見文獻[16]。

異步旋轉布風裝置如圖2所示，裝置主要由內嵌式中心轉盤布風板與外套式行星轉環布風板組成，其中，中心轉盤布風板表面均勻布置有13個半球形風帽，行星轉環布風板均勻布置有20個半球形風帽，直徑為1 cm的風帽上均勻地開有8個直徑為 1 mm 的小孔。異步旋轉布風板的旋轉原理構造圖如圖3所示，變頻電機A(0～1400 r/min)通過固定在布風板軸心位置的傳動轉軸帶動中心轉盤布風板轉動，變頻電機B通過連接皮帶間接帶動固定在布風板外環的傳動支架，實現行星轉環布風板的轉動。實驗選用變頻電機額定轉速0～1400 r/min，轉軸通過聯軸器連接速速比為1/7的減速器，實現布風裝置在0～200 r/min的范圍內轉速可調。分別調節變頻電機A、B轉速，即可實現旋轉布風裝置的異步轉動，異步旋轉布風裝置的布風工況如表2所示，其中N1為中心轉盤布風板的轉速，N2為行星轉環布風板的轉速。當N1=N2時，中心轉盤布風板與行星轉環布風板同步轉動。

圖1 油頁巖異步旋轉干燥實驗裝置示意圖Fig.1 Flow diagram of the experimental setup of the oil shale asynchronous rotation drying system1—Computer; 2—Data acquisition; 3—Fluidized bed; 4—Sampling aperture; 5—ERT; 6—Temperature sensor; 7—Vortex flow meter; 8—Roots blower; 9—A synchronous rotating air distributor; 10—Valve; 11—Air heater

圖2 異步旋轉布風板示意圖Fig.2 Asynchronous rotating air distributor

圖3 異步旋轉布風板旋轉原理構造圖 Fig.3 Schematic diagram of the asynchronous rotating air distributor

表2 異步旋轉布風裝置布風工況Table 2 Different air supply conditions of the asynchronous rotating air distributor

Cond—Condition

顆粒在流化內的分布狀態通過安裝在高度為0.5 m的電容層析成像裝置(ERT)進行測量，ERT由數據采集頻率為2幀/s、激勵頻率為9.6 Hz的電容陣列傳感器組成，多個傳感器的組合可以反映顆粒在橫截面的分布狀態。

1.3 水分測定實驗

油頁巖干燥過程的含水率通過電子天平(320 g/0.1 mg)進行讀取。每次取制備好的油頁巖樣品 5 kg 置于流化床內，在設定的布風工況下干燥，干燥時間為90 min，每隔2 min從取料口提取料樣，每次取樣的質量為5 g，將樣品放入恒溫干燥箱內進行干燥，直至前后兩次稱量的質量之差達到0.001 g，即干燥完畢。在整個干燥過程中取樣的總質量不超過干燥樣品質量的5%，因此對油頁巖的干燥過程不會產生影響。在油頁巖干燥過程中的質量損失，即為不同時刻油頁巖含水率(wR)，其計算式見式(1)。

(1)

式(1)中，wt為t時刻物料干基含水量，g/g；w0為物料原始含水量，g/g；we為物料平衡含水量，g/g；t為物料干燥時間，min。

2 結果與討論

油頁巖在布風板同步旋轉布風工況(N1=N2)下的干燥情況見圖4。由圖4可以看到，隨著轉速的逐漸增加，所需干燥時間不斷縮短，干燥效率有所提高，這是由物料在流化床內的分布狀態所決定的。經過研究[17-18]表明，當轉速N1=N2=0時，物料在流化床內主要呈現環-核態流動結構。圖5為同步旋轉布風工況(N1=N2)下油頁巖顆粒的電容層析成像結果，其中電容層析成像裝置測量的顆粒流動圖像顏色越深，代表顆粒的濃度越高。由圖5(a)可以看到，流化床中心區域物料濃度的分布較均勻；而靠近壁面的位置，由于物料之間的相互碰撞以及壁面的阻礙效應，局部氣速明顯低于中心區域，造成壁面位置發生物料的堆積。如圖5(b)所示，適當地增加布風板轉速，可以有效地混合中心區域與壁面區域的氣流，改善氣流的不均勻分布狀態，減輕壁面區域物料的堆積效應，增加物料與干燥氣流之間的接觸面積，因此干燥效率有所增加，這一點與Sobrino等[15]的研究結果相一致。然而轉速增加到N1=N2=100 r/min時，干燥效率并沒有達到極大值，繼續增加旋轉速率，干燥所需時間繼續縮短。當轉速增加到N1=N2=150 r/min時，干燥效率達到極大值；隨著轉速的繼續增加，干燥效率開始下降，這主要是由于布風板轉速過高，旋轉離心力對物料的作用明顯超過布風板旋轉均流的作用，此時大量物料重新聚集在流化床的壁面位置，干燥氣流無法實現對物料的有效干燥。

圖4 同步旋轉布風工況(N1=N2)下油頁巖顆粒含水率隨干燥時間的變化Fig.4 Moisture content of oil shale particles changing with time at N1=N2

圖5 同步旋轉布風工況(N1=N2)下油頁巖顆粒的電容層析成像結果Fig.5 The image of ERT measurement result of oil shale particles at N1=N2N1,N2/(r·min-1): (a) N1=N2=0; (b) N1=N2=100; (c) N1=N2=150; (d) N1=N2=175R—Radius of fluidized bed’s cross section; r—Current position of fluidized bed’s cross section

分別調整內嵌式中心轉盤與外套式轉盤布風板的轉速，進行布風板異步旋轉布風工況下的干燥特性研究，由于布風工況數量較大，對實驗結果進行對比篩分選取典型工況，在布風板異步旋轉布風工況下干燥的油頁巖含水率隨著時間的變化如圖6所示，圖7為異步旋轉布風工況下油頁巖顆粒的電容層析成像結果。由圖6可以看出，與同步轉動布風的最優工況N1=N2=150 r/min相比，轉速N1=150 r/min、N2=125 r/min及轉速N1=175 r/min、N2=150 r/min的干燥時間均有所延長，即減小外套式轉環布風裝置的轉速，提高內嵌式轉盤布風裝置的轉速都不能起到提高流化床干燥效率的作用，這一現象可以通過圖7(a)和(b)加以解釋。可以發現，上述2種布風方式會促使中心區域的旋流強度增加，導致物料向壁面位置遷移，破壞物料的干燥效果。相比之下，N1=125 r/min、N2=150 r/min對應的干燥時間開始縮短，這說明相比于最優的同步旋轉干燥工況，干燥效率在異步旋轉干燥狀態下可以被進一步提高，如圖7(c)所示，降低內嵌式中心轉盤布風裝置的轉速可以有效地降低流化床中心區域的旋轉強度，從而促使邊壁區域的物料向中心區域遷移，物料的分布狀態逐漸趨于均勻，隨著內嵌式中心轉盤布風裝置轉速的降低，N1=75 r/min、N2=150 r/min對應的干燥時間開始延長，此時中心轉盤布風裝置的轉速明顯低于行星轉環布風裝置的轉速，中心區域的旋轉強度較低，如圖7(d)所示，邊壁區域的物料在流化床中心區域發生了堆積，阻礙了干燥氣流與物料的有效接觸。增加外套式行星轉環布風裝置的轉速，N1=150 r/min，N2=175 r/min對應的干燥時間亦有所縮短，但干燥時間明顯大于降低內嵌式轉盤布風裝置轉速所需的干燥時間，而且繼續增加外套式行星轉環布風裝置的轉速至N1=150 r/min、N2=200 r/min，干燥效率急劇降低，說明此時物料向中心區域遷移狀態提前發生。

圖6 異步旋轉布風工況下油頁巖顆粒含水率隨干燥時間的變化Fig.6 Moisture content of oil shale particles changing with time at the asynchronous rotation condition

圖7 異步旋轉布風工況下油頁巖顆粒的電容層析成像結果Fig.7 The image of ERT measurement result of oil shale particles at the asynchronous rotation conditionN1,N2/(r·min-1): (a) N1=150, N2=125; (b) N1=175, N2=150; (c) N1=125, N2=150; (d) N1=150, N2=175

油頁巖樣品的干燥速率(RD)計算公式如式(2)所示。

(2)

式(2)中，wt+dt為時間間隔為dt的物料干基含水率，g/g。

實驗計算得到的布風板同步旋轉布風工況下油頁巖干燥速率(RD)隨含水率(wR)的變化如圖8所示。從圖8可以看出，油頁巖在流化床內的干燥過程主要分為3個干燥階段，即升速干燥階段、恒速干燥階段和降速干燥階段。在升速干燥階段，物料表面的水分首先受熱蒸發，導致表面溫度低于內部溫度，推動水分向外表面遷移。隨著轉速的增加，干燥曲線的斜率逐漸增加，這是由于在旋轉過程中物料的分布逐漸趨于均勻，物料的表面與干燥氣流的接觸面積隨之增加，較大的溫度梯度產生較大的水分推動力，干燥速率隨之增加。當干燥曲線達到最大值后，曲線的斜率在一定時間內保持不變，這一階段稱為恒速干燥階段，這時物料表面水分散失速率等于內部水分的擴散速率。曲線斜率減小的階段為降速干燥階段，該階段以脫去結合水為主，由于此時物料內部水分向表面移動的速率低于表面水分的蒸發速率，因此干燥速率逐漸減小，且油頁巖在降速干燥階段干燥速率呈現由慢到快再到慢的降低趨勢。原因可能是剛進入降速階段時，油頁巖中仍存在部分非結合水分，到了中期，由于油頁巖粒度較小，水分的內部遷移阻力較大，水分來不及向表面擴散，致使干燥速率下降較快，到了后期，由于油頁巖的干燥速率已經很低，所以呈現下降幅度變緩的趨勢，如果停留時間足夠長，物料達到平衡含水量。當轉速增加到N1=N2=150 r/min時，干燥速率達到最大值0.0615 gH2O/(g·min)，此時為同步轉動布風的最優干燥工況；繼續增加轉速，物料聚集在流化床的壁面區域，阻礙了干燥氣流對物料的有效干燥，干燥速率隨之降低。

圖8 同步旋轉布風工況(N1=N2)下油頁巖顆粒干燥速率(RD)隨含水率(wR)的變化Fig.8 Drying rate (RD) changing with moisture content (wR) of oil shale particles at N1=N2

相比之下，布風板異步轉動布風工況的干燥速率明顯大于同步轉動布風工況的干燥速率，如圖9所示。布風板異步轉動布風最優工況N1=100 r/min、N2=150 r/min的最大干燥速率達到0.1145 gH2O/(g·min)，為最佳同步轉動布風干燥速率的1.86倍。結合圖8、圖9可以看到，布風板異步轉動布風干燥與同步轉動布風干燥的曲線形狀明顯不同，異步轉動布風的干燥速率經歷升速干燥階段后幾乎不具備恒速干燥階段直接進入降速干燥階段，這是因為異步轉動布風干燥不同于一般的干燥方式，通過調節內嵌式中心轉盤與外套式轉環的轉速會使堆積在壁面區域的物料引向流化床的中心區域，明顯提高升速干燥階段物料的干燥強度，在升速干燥階段物料內部的非結合水分已經全部散失殆盡，從而直接進入以結合水在顆粒內部擴散為主要特征的降速干燥階段。進一步降低內嵌式中心轉盤布風板的轉速，由于流化床的干燥能力迅速下降，恒速干燥階段隨之開始產生。

圖9 異步旋轉布風工況下油頁巖顆粒干燥速率(RD)隨含水率(wR)的變化Fig.9 Drying rate (RD) changing with moisture content (wR) of oil shale particles at the asynchronous rotation conditionN1, N2/(r·min-1): N1=150, N2=150; N1=75， N2=150； N1=125, N2=150; N1=50， N2=150； N1=100, N2=150; N1=25， N2=150

3 實驗數據干燥動力學分析

為尋求適于描述異步轉動布風工況下柳樹河油頁巖的理想干燥模型，通過對現有干燥模型進行總結，得到如表3所示的9種常用干燥模型，對物料在典型異步轉動布風配比工況的干燥實驗數據進行動力學模擬。以相關系數R2、卡方誤差χ2作為指標，評價不同模型與實驗數據的擬合效果，R2越接近于1、χ2越小，數學模型的匹配程度越好，其中χ2和R2可按式(3)、式(4)計算。

(3)

(4)

表3 干燥模型方程Table 3 Drying mathematical models

a,b,k,n—Coefficients to be determined in the models equation

干燥模型模擬結果如表4所示。從表4可知，隨著內嵌式中心轉盤布風裝置轉速的降低，即外套式轉環與中心轉盤布風裝置轉速差的持續增加，其相關系數先增大后減小，其卡方誤差先減小后增加，盡管多個模型模擬的相關系數都達到0.9以上，卡方誤差都接近于0。但比較可見，Two-term模型在對所有異步轉動布風工況下的實驗數據進行模擬時均取得了較好的擬合效果，其相關系數在0.997以上，大于其他模型擬合值；卡方誤差均在0.0003以下，小于其他模型的擬合值。

表4 干燥模型擬合結果Table 4 Fitting results of drying mathematical models

Con—Condition

圖10為異步轉動布風工況下物料含水率實驗值與Two-term模擬值的對比。其數據基本在斜率為45°的直線周圍，從另一個方面說明Two-term模型較適用于異步轉動布風工況下油頁巖干燥過程的模擬。模擬值與實驗值的平均相對誤差為6.6%，隨著內嵌式中心轉盤布風裝置轉速的降低，模擬效果逐漸降低。這一現象歸咎于兩方面原因：一方面由于相應工況下的干燥強度有所提高，物料內部空隙中吸附的氣體隨水分一起擴散出來；另一方面是由于每次實驗結束后物料仍存在少量結合水。物料質量減小部分包括吸附氣體質量、水分蒸發及殘留水分的質量之和，而模型中所計算的物料減小質量僅為水分蒸發質量，導致模擬值小于實驗值。當中心轉盤布風裝置的轉速工況降低為N1=75 r/min、N2=150 r/min時，干燥強度開始降低，模擬效果有所改善。

圖10 異步旋轉布風工況下油頁巖顆粒含水率(wR)的Two-term模型驗證結果Fig.10 Moisture content (wR) of oil shale particles verification results by the two-term drying model at asynchronous rotation condition

4 結 論

(1)隨著布風板同步旋轉布風轉速適當的增加，物料與干燥氣流之間的接觸面積得到改善，干燥時間不斷縮短。當轉速增加到N1=N2=150 r/min時，物料達到最佳干燥狀態；繼續增加轉速，會導致大量物料聚集在流化床的壁面位置，破壞流化床的有效干燥能力。

(2)降低內嵌式中心轉盤布風裝置的轉速可以降低流化床中心區域的氣流強度，進一步改善物料的分布狀態。N1=75 r/min、N2=150 r/min為異步旋轉布風工況的最佳干燥狀態，相比于布風板同步旋轉，干燥效率被進一步提高。繼續降低中心轉盤轉速，會導致物料在流化床中心區域發生堆積，阻礙干燥氣流與物料的有效接觸。

(3)物料的流化干燥過程主要分為升速干燥、恒速干燥和降速干燥3個階段，布風板異步旋轉布風最佳工況的干燥速率為同步旋轉布風最佳工況的1.7倍。合理的調節內嵌式中心轉盤與外套式轉環的轉速會強化升速干燥階段物料內部非結合水的散失殆盡，幾乎不存在恒速干燥階段，直接進入降速干燥階段。

(4)Two-term干燥模型對典型異步旋轉工況下物料的干燥實驗數據進行擬合的效果最好。隨著流化床干燥能力的逐漸強化，導致物料空隙內部吸附的氣體隨水分一同擴散，顆粒質量損失包括吸收氣體及水分蒸發的質量兩部分，模擬結果小于實驗結果。

[1] 王擎， 張宏喜， 隋義， 等. 蒙脫石熱反應性及對干酪根熱解的影響[J].石油學報(石油加工)， 2016， 32(1)： 50-55. (WANG Qing, ZHANG Hongxi, SUI Yi, et al. Thermal properties of smectite and its effect on kerogen pyrolysis[J].Acta Petrolei Sinica (Petroleum and Processing Section), 2016, 32(1): 50-55.)

[2] 柏靜儒， 白章， 王擎， 等. 油頁巖固體熱載體綜合利用系統工藝模擬[J].石油學報(石油加工)， 2014, 30(5): 902-908. (BAI Jingru, BAI Zhang, WANG Qing, et al. Process simulation of oil shale comprehensive utilization system on solid heat carrier technology[J].Acta Petrolei Sinica (Petroleum and Processing Section), 2014, 30(5): 902-908.)

[3] 錢家麟, 尹亮. 油頁巖： 石油的補充能源[M].北京： 中國石油出版社， 2008: 56-70.

[4] LI W H, LU S F, XUE H T, et al. Oil content in argillaceous dolomite from the Jianghan Basin, China: Application of new grading evaluation criteria to study shale oil potential[J].Fuel, 2015, 143(3): 424-429.

[5] 周新航， 王擎， 柏靜儒， 等. 樺甸油頁巖微波干餾及其升溫特性的研究[J].東北電力大學學報， 2014， 34(5)： 43-47. (ZHOU Xinhang, WANG Qing, BAI Jingru, et al. Study of Huadian oil shale’s microwave dry distillation and temperature rising characteristic[J].Journal of Northeast Dianli University, 2014，34(5): 43-47.)

[6] LAW C L, SITI M T, WAN R W D. The effect of vertical internal baffles on fluidization hydrodynamics and grain drying characteristics[J].Chinese Journal of Chemical Engineering, 2004, 12(6): 801-808.

[7] 孫佰仲， 宮殿臣. 流化床油頁巖流化特性的數值模擬研究[J].東北電力大學學報， 2013, 33(5): 10-14. (SUN Baizhong, GONG Dianchen. Numerical simulation on the flow characteristic of oil shale in fluidized bed[J].Journal of Northeast Dianli University, 2013, 33(5): 10-14.)

[8] 于德平， 張玉明，楊運泉. 流化床中造紙黑液催化石油焦氣化特性[J].石油學報(石油加工)， 2013， 29(3)： 438-446. (YU Deping, ZHANG Yuming, YANG Yunquan, et al. Gasification characteristic of petroleum coke catalyzed by biack liquor in a fluidized bed[J].Acta Petrolei Sinica (Petroleum and Processing Section), 2013, 29(3): 438-446.)

[9] 金涌， 祝京旭， 汪展文， 等. 流態化工程原理[M].北京： 清華大學出版社， 2001: 36-40.

[10] 周云龍， 楊寧. 流化床顆粒團聚物水分測量新技術[J].化工學報， 2014, 65(10): 3878-3883. (ZHOU Yunlong, YANG Ning. A new method for measuring moisture distribution in fluidized bed[J].CIESC Journal, 2014, 65(10): 3878-3883.)

[11] 李凡, 馮連芳, 顧雪萍, 等. 氣固攪拌流化床的床層壓降[J].高校化學工程學報, 2002, 16(4): 384-388.(LI Fan, FENG Lianfang, GU Xueping, et al. Press drop in the gas-solid stirred fluidized bed[J].Journal of Chemical Engineering of Chinese Universities, 2002, 16(4): 384-388.)

[12] 宮國清， 呂岳琴， 余華瑞， 等. 脈沖式流化床干燥器操作性能的研究[J].南京林業大學學報， 1997, 21(S1): 43-45. (GONG Guoqing, Lü Yueqin, YU Huarui, et al. A study operation characteristics of pulse fluidized bed dryer[J].Journal of Nanjing Forestry University, 1997, 21(S1): 43-45.)

[13] KUIPERS N J M, STAMHUIS E J, BEENACKERS A A C M. Fluidization of potato starch in a stirred vibrating fluidized bed[J].Chemical Engineering Science, 1996, 51(11): 2727-2732.

[14] MADHIYANON T, LAPIRATTANAKUN A, SATHITRUANGSAK P, et al. A novel cyclonic fluidized bed combustor: Combustion and thermal efficiency, temperature distributions, combustion intensity, and emission of pollutants[J].Combust & Flame, 2006, 146(1-2): 232-245.

[15] SOBRINO C, ALMENDROS-IBANEZ J A, SANTANA D, et al. Fluidization of group B particles with rotating distributor[J].Powder Technology, 2008, 181(3): 273-280.

[16] 范新欣， 王寶和， 于才淵， 等. 基于氣固流態化原理的油頁巖干燥動力學[J].化學工程， 2010, 38(10): 210-214. (FAN Xinxin, WANG Baohe, YU Caiyuan, et al. Drying kinetics of oil shale based on gas-solid fluidization principles[J].Chemical Engineering, 2010, 38(10): 210-214.)

[17] SAEIDEH H, FARAMARZ H. Drying kinetics of coated sodium percarbonate particles in a conical fluidized bed dryer[J].Powder Technology, 2015, 269(8): 30-37.

[18] SRINIVAS G, SETTY Y P. Drying behavior of uniform and binary mixture of solids in a batch fluidized bed dryer[J].Powder Technology, 2013, 241(3): 181-197.

[19] HENDERSON S M, PABIS S. Grain drying theory Ⅱ Temperature effects on drying coefficients[J].Journal of Agricultural Engineering Research, 1961, 6(3): 169-174.

[20] OVERHULES D G, WHITE G M, HAMILTION M E, et al. Drying soybeans with heated air[J].Amer Soc Agr Eng Trans Asae, 1973, 16(1): 112-113.

[21] VERMALR, BCKLIN R A, EDAN J B, et al. Effects of drying air temperature on rice drying models[J].Transactions of the ASAE, 1985, 28: 296-301.

[22] AYENSU A. Dehydration of food crops using a solar dryer with convective heat flow[J].Solar Energy, 1997, 59(4-6): 121-126.

[23] AKGUN A, DOYMAZ I. Modeling of olive cake thin-layer drying process[J].Journal of Food Engineering, 2005, 68(4): 455-461.

[24] TOGRUL I T, PEHLIVAN D. Mathematical modeling of solar drying of apricots in thin layers[J].Journal of Food Engineering, 2002, 55(3): 209-216.

[25] MIDILLI A, KUCUK H, YAPAR Z. A new model for single layer drying[J].Drying Technology, 1997, 20(7): 1503-1513.

[26] BABALIS S J, PAPANICOLAOU E, KYRIAKIS N, et al. Evaluation of thin-layer drying models for describing drying kinetics of figs (Ficus carica)[J].Journal of Food Engineering, 2006,75(2): 205-214.

Asynchronous Rotating Drying Technology of Liushuhe Oil Shale

YANG Ning1, ZHOU Yunlong1, QI Tianyu2, ZHU Xiaoyu1

(1.CollegeofEnergyResourceandPowerEngineering,NortheastDianliUniversity,Jilin132012,China;2.InstituteofUnclearandNewEnergyTechnology,TsinghuaUniversity,Beijing100084,China)

A novel asynchronous rotating air distributor was designed to replace the conventional distributor in this study to optimize the drying ability of fluidized bed and strengthen the drying performance of oil shale particles. The rotating speed of the asynchronous rotating air distributor with embedded center disk and encircled disk was regulated to achieve the various air supply conditions. The impacts of different drying conditions on the drying characteristics of Liushuhe oil shale particles were studied with the help of Electrical Capacitance Tomography (ERT). The dynamics of experimental data was analyzed with 9 common drying models. The results indicate that decrease of embedded center disk and increase of the encircled disk air distributor rotating speed can improve particles distribution in the fluidized bed, and reduce the drying time. The drying process of oil shale particles consists of three periods, i.e. raising drying rate period, constant drying rate period and falling drying rate period. The method that regulating air distributor rotating speed reasonably leads to accelerate the particles from raising drying rate period to the falling drying rate period directly. The two-term model fits the oil shale particles drying simulation within 9 drying models at different air supply conditions. Yet the air absorbed in the particles pores is diffused together with moisture evaporation, and a small amount of bound moisture remains after each experiment; consequently, the results of drying simulation are less than the experimental values.

fluidized bed; asynchronous rotation; oil shale; drying technology

2016-03-15

國家自然科學基金項目(51276033，51541608)資助

楊寧，男，講師，從事油頁巖干燥技術及氣固參數檢測研究；E-mail:m13620752851@163.com

周云龍，男，教授，從事多相流動參數檢測研究；Tel:0432-64807495;E-mail:1020219438@qq.com

1001-8719(2017)03-0497-10

TK121

A

10.3969/j.issn.1001-8719.2017.03.014