包旭,黃宇,黃治軍,段元強,段倫博
(1.能源熱轉(zhuǎn)換及其過程測控教育部重點實驗室,東南大學 能源與環(huán)境學院,江蘇 南京,210096;2.江蘇方天電力技術(shù)有限公司,江蘇 南京,211100)
我國CO2的主要排放來源于能源行業(yè),煤炭資源消耗占比96%,其中約一半煤炭資源用于燃煤發(fā)電。富氧燃燒技術(shù)[1]是最具經(jīng)濟優(yōu)勢的燃煤電站CO2捕集技術(shù)之一,但較高的能耗限制了富氧燃燒技術(shù)進一步發(fā)展。增壓富氧燃燒技術(shù)[2-3]是解決能耗問題的有效手段之一。相比傳統(tǒng)富氧燃燒技術(shù),增壓富氧燃燒系統(tǒng)均在高壓下進行,具有眾多優(yōu)勢[4-6]:
1) 避免漏風,降低純化壓縮系統(tǒng)功耗;
2) 避免壓力變化導致的能量損失;
3) 提高煙氣中水蒸氣露點,汽化潛熱易回收;
4) 降低鍋爐等設(shè)備尺寸。
相關(guān)研究[7-8]表明,流化床鍋爐在高壓下運行,減小爐膛尺寸且變大容積熱負荷,從而導致受熱面布置困難,密相區(qū)布置水平埋管[9-10]是解決受熱面布置面積不足的有效方法之一。
準確測定換熱管傳熱系數(shù)對流化床鍋爐設(shè)計至關(guān)重要。目前,針對流化床密相區(qū)水平埋管傳熱特性的研究較多,但研究大多集中于冷態(tài)工況,而針對增壓富氧條件下開展的熱態(tài)增壓研究報道較少。KIM等[11]在以石英砂為床料的冷態(tài)增壓鼓泡床中研究水平埋管傳熱,發(fā)現(xiàn)壓力由0.10 MPa 提升至0.74 MPa時,最高平均傳熱系數(shù)提升了29.68%,傳熱系數(shù)隨流化風速增加呈現(xiàn)先提升后降低的趨勢,換熱管兩側(cè)傳熱效果優(yōu)于迎風面和背風面;MASOUMIFARD等[12]也得到了相似的結(jié)論。而OLSSON等[13]在結(jié)構(gòu)相似的冷壓增壓流化床中采用電加熱平衡法,利用電加熱埋管研究埋管與床層間的傳熱特性,發(fā)現(xiàn)管束結(jié)構(gòu)對水平管的傳熱特性產(chǎn)生影響,降低水平節(jié)距導致傳熱變差。STENBERG等[14]認為密相區(qū)中床料種類和床料粒徑對于水平管傳熱的影響更為顯著,發(fā)現(xiàn)小粒徑床料對傳熱效果有顯著增強作用,且在所有床料中鈦鐵礦傳熱效果最好,這與MCINTYRE等[15]在冷態(tài)高壓下得到的實驗結(jié)果一致。
目前,人們對流化床密相區(qū)傳熱系數(shù)的研究大多集中于冷態(tài)高壓或熱態(tài)常壓。冷床熱管的熱量傳遞方向與實際工業(yè)鍋爐不同,其所獲得的換熱系數(shù)及換熱關(guān)聯(lián)式不能指導增壓流化床鍋爐埋管設(shè)計。而高壓熱態(tài)實驗由于其實驗臺設(shè)計和操作難度較大,相關(guān)研究報道較少。因此,本文首先研制了熱態(tài)增壓流化床實驗平臺,采用夾套式布置方法克服高溫高壓下材料的選擇困難;其次,選用載氧體輔助燃燒技術(shù)中應(yīng)用廣泛的鈦鐵礦為床料,研究密相區(qū)水平埋管外表面?zhèn)鳠釞C理及不同運行工況對傳熱特性的影響;最后,根據(jù)密相區(qū)埋管實驗結(jié)果,提出適用于不同床料的增壓流化床密相區(qū)水平埋管傳熱關(guān)聯(lián)式。
圖1 所示為換熱平臺結(jié)構(gòu)示意圖。由圖1(a)可見:整個實驗臺由壓力罐、流化床爐體、中頻加熱器、管殼式換熱器、背壓閥和水平換熱管等組成,并配備了氣路系統(tǒng)、電加熱系統(tǒng)、密封系統(tǒng)和數(shù)據(jù)采集系統(tǒng)。實驗臺整體采用夾套式設(shè)計,流化床置于壓力容器內(nèi)部,保證流化床在高溫高壓下穩(wěn)定運行。為增大換熱面積,埋管采用水平U形管結(jié)構(gòu),外徑為8 mm,壁厚為2 mm,在埋管外壁面周向每隔90°開0.7 mm深的凹槽,并放置直徑為0.5 mm 的K 型熱電偶測量管壁溫,熱電偶測點位于爐膛中心處,保證測量的準確性。管內(nèi)工質(zhì)為水,采用浮子流量計控制流量,實驗過程中流量設(shè)為70 L/h。高壓空氣進入壓力罐后,在壓差作用下進入底部中頻加熱器下部封頭,從而進入爐膛作為流化風,冷卻后經(jīng)背壓閥排入外界環(huán)境。背壓閥調(diào)節(jié)范圍為0.1~2.0 MPa,在實驗過程中通過調(diào)節(jié)背壓閥來控制爐膛壓力。有關(guān)實驗裝置的詳細描述見文獻[16]。
圖1 換熱平臺結(jié)構(gòu)示意圖Fig.1 Schematic illustration of heat exchange platform
本工作采用鈦鐵礦作為實驗床料。通過與石英砂進行對比,可以更加形象地展示鈦鐵礦作為床料的埋管傳熱特性。2種床料具有相同的粒徑范圍(450~500 μm),基于氣體置換法測得鈦鐵礦的密度為4 044 kg/m3,石英砂的密度為2 650 kg/m3,2種床料在啟動前的靜床高均為300 mm。臨界流化風速采用李皓宇[17]提出的關(guān)聯(lián)式進行計算。其中壓力為0.1~0.7 MPa、以石英砂為床料的實驗數(shù)據(jù)來自本文作者前期研究[16],壓力大于0.9 MPa的數(shù)據(jù)由本次實驗測得。
在忽略管壁軸向?qū)崆闆r下,水平埋管與床層之間的傳熱分為3 個過程:床層與外管壁間傳熱、埋管壁面徑向一維導熱、水與埋管內(nèi)壁面對流換熱。當傳熱過程達到穩(wěn)定后,水的吸熱量即為床層向埋管外表面?zhèn)鬟f的熱量。首先,根據(jù)實驗測得的水平管進出口水溫Tout和Tin以及流量,可以通過式(1)計算水的吸熱量Q;其次,根據(jù)計算所得熱量Q和實驗測得的水平管進出口管壁周向溫度Ti以及密相區(qū)床層溫度Tb,可以通過式(2)和式(3)分別計算總傳熱系數(shù)ho和周向局部傳熱系數(shù)ho,l。
式中:Q為管內(nèi)工質(zhì)帶走的熱量,W;m為水的質(zhì)量流量,kg/s;cp,w為水的定壓比熱容,J/(kg·℃);Tout和Tin分別為水的出口溫度和進口溫度,℃;Tb為床層溫度,℃;Ti為測點溫度,℃;Ts為埋管外壁面平均溫度,℃;-Ts 為埋管外壁面同周向平均溫度,℃。
在熱穿透模型和接觸熱阻模型的基礎(chǔ)上,結(jié)合增壓流化床內(nèi)乳化相的熱物性參數(shù),BAO等[18]提出了一個適用于熱態(tài)增壓流化床水平埋管的半經(jīng)驗傳熱模型,HUANG等[19-20]通過實驗驗證了該模型的準確性。該傳熱模型如下:
式中:havg為埋管表面平均傳熱系數(shù),W/(m2·K);ξ為管束結(jié)構(gòu)特性系數(shù);δb為鼓泡床中氣相體積分數(shù);hgc為氣體對流傳熱系數(shù),W/(m2·K);hpc為顆粒對流傳熱系數(shù),W/(m2·K);hr為輻射傳熱系數(shù),W/(m2·K);hgc,top和hgc,bottom分別為埋管頂部和底部氣體對流傳熱系數(shù),W/(m2·K);tc為乳化團-管壁接觸均方根時間,s;t0為熱穿透時間,s;σ為Stefan Boltamann常數(shù),取5.67×10-8W/(m2·K4);Teb為有效床溫(床層溫度的0.85 倍),K;eb為床層發(fā)射率,本文取0.9;es為床料發(fā)射率,本文取0.95。
圖2所示為實驗工況下床層與埋管之間傳熱系數(shù)的實驗值和模型預(yù)測值。從圖2可見:床料更換為鈦鐵礦,該模型對床料的水平埋管換熱系數(shù)的預(yù)測相對誤差仍在20%以內(nèi)。本文將結(jié)合該傳熱模型進一步分析爐膛壓力、床層溫度和流化風速對床層與埋管之間傳熱系數(shù)的影響。
圖2 換熱模型對以鈦鐵礦為床料的水平埋管換熱系數(shù)的預(yù)測相對誤差Fig.2 Prediction relative error of heat transfer model with ilmenite as bed material of immersed tube
在床層溫度為800 ℃、流化數(shù)為4 的工況下,不同爐膛壓力下測得以鈦鐵礦為床料的床層與埋管間平均傳熱系數(shù)如圖3(a)所示。從圖3(a)可見:傳熱系數(shù)隨爐膛壓力增加而增加,這與KIM等[11]在高壓冷態(tài)下的實驗結(jié)果保持一致。增加壓力對于傳熱系數(shù)的促進作用可以從2 個方面進行解釋。由于實驗平臺在高溫高壓下運行,對系統(tǒng)的密封性具有嚴格要求,因此,通過實驗分別測量壓力對氣體對流傳熱系數(shù)和顆粒對流傳熱系數(shù)的影響十分困難。然而,結(jié)合式(4)~(7),可以計算和對比實驗中氣體對流傳熱系數(shù)和顆粒對流傳熱系數(shù)。如圖3(b)所示,隨著爐膛壓力增加,氣體對流傳熱系數(shù)和顆粒對流傳熱系數(shù)均單調(diào)遞增,這歸因于氣體密度增加和顆粒流化質(zhì)量改善。
相比于石英砂,以鈦鐵礦為床料的埋管具有更高的傳熱系數(shù)(見圖3(a))。鈦鐵礦更好的傳熱效果仍可以從這2 個方面進行解釋。值得注意的是,為盡可能排除床料種類變化對流化狀態(tài)的影響[21],2 種床料下的換熱實驗均在4 倍流化數(shù)下開展。從圖3(b)可以看出,以鈦鐵礦為床料的氣體對流傳熱系數(shù)比石英砂的略高,這是因為在相同粒徑下,鈦鐵礦具有更大的臨界流化風速。而顆粒對流是2種床料埋管平均傳熱系數(shù)存在差異的主要原因。在相同粒徑下,鈦鐵礦較大的密度導致其具有更大的單位體積熱容(鈦鐵礦和石英砂的單位體積熱容分別為3 640 kJ/m3和2 120 kJ/m3)。因此,在與埋管外表面進行換熱時,相比于石英砂其溫度下降更慢,與埋管之間存在相對更大的溫差,從而增強了顆粒與埋管之間的熱量交換。
圖3 壓力對埋管傳熱系數(shù)的影響Fig.3 Effect of pressure on heat transfer coefficient of immersed tube
而壓力變化也會導致埋管周向熱流密度分布不均,因此,增壓后埋管周向傳熱系數(shù)與常壓下有所區(qū)別。圖4所示為埋管周向局部傳熱系數(shù)在不同壓力下的分布規(guī)律。從圖4可見:以鈦鐵礦為床料的埋管周向傳熱系數(shù)分布特征與以石英砂為原料時大致相同,即水平埋管兩側(cè)的傳熱系數(shù)高于埋管底部和頂部的傳熱系數(shù),驗證了前人基于密相區(qū)埋管表面顆粒運動特性給出預(yù)測結(jié)果[22]的合理性。增加爐膛壓力并沒有明顯改變顆粒在埋管周向的運動特性,埋管周向局部傳熱系數(shù)的分布差異主要取決于顆粒流動特性影響。后續(xù)在床層溫度和流化風速的討論中不再考慮埋管周向傳熱系數(shù)。
圖4 不同壓力下埋管周向傳熱系數(shù)分布規(guī)律Fig.4 Circumferential heat transfer coefficient distribution law under different pressures of immersed tube
圖5 所示為爐膛壓力為0.1 MPa、流化數(shù)為4時,以鈦鐵礦和石英砂為床料的埋管平均傳熱系數(shù)和輻射傳熱系數(shù)隨床層溫度的變化情況。從圖5可見:隨著床層溫度增加,2種床料下的埋管平均傳熱系數(shù)均顯著增加。這一方面得益于氣體熱導率增加,另一方面則是輻射傳熱系數(shù)隨床溫增加而顯著增加。通過對比可以發(fā)現(xiàn),以鈦鐵礦為床料的埋管輻射換熱系數(shù)略大于石英砂的輻射換熱系數(shù)。這主要是由于鈦鐵礦較高的密度導致床層與埋管間的換熱能力更強,埋管表面溫度更高。由輻射換熱模型可知,埋管壁溫增大有助于增加輻射傳熱系數(shù)。盡管輻射傳熱系數(shù)在總傳熱系數(shù)中的占比在高溫下顯著提高,但即使在800 ℃下2種床料的輻射傳熱份額仍未超過13%(鈦鐵礦和石英砂的輻射傳熱份額分別為12.39%和11.81%)。這說明在增壓流化床鍋爐實際運行過程中,床層與埋管之間的傳熱方式依舊以對流傳熱為主。
圖5 床層溫度對埋管平均傳熱系數(shù)和輻射傳熱系數(shù)的影響Fig.5 Effect of bed temperature on average and radiation heat transfer coefficient of immersed tube
在爐膛壓力為0.9 MPa、床層溫度為800 ℃下,流化風速對埋管平均傳熱系數(shù)的影響如圖6所示。從圖6可見:鈦鐵礦和石英砂隨著流化數(shù)增加均呈現(xiàn)了相同的變化趨勢,即低流化數(shù)下,埋管平均傳熱系數(shù)隨著流化數(shù)提高而快速增加;在5倍流化數(shù)時,埋管平均傳熱系數(shù)達到峰值;隨著流化數(shù)進一步增加,埋管平均傳熱系數(shù)卻呈現(xiàn)逐漸下降的趨勢。這是因為隨著流化數(shù)增加,埋管表面顆粒更新頻率逐漸增加,但埋管附近顆粒濃度逐漸減小,因此,存在一個最大值。對比不同流化數(shù)下2種床料的埋管傳熱系數(shù)可以看出,高流化數(shù)下鈦鐵礦的平均傳熱系數(shù)高于石英砂的平均傳熱系數(shù)。而在流化數(shù)等于2時,鈦鐵礦具有與石英砂相似的傳熱系數(shù)。這可能是由于低流化數(shù)下埋管表面的氣泡頻率較低,顆粒與埋管的接觸頻率也較低,導致鈦鐵礦的對流傳熱優(yōu)勢并不顯著。
圖6 流化風速對埋管平均傳熱系數(shù)的影響Fig.6 Effect of fluidization velocity on average heat transfer coefficient of immersed tube
相比于經(jīng)過數(shù)學推導得到的半經(jīng)驗公式,經(jīng)驗公式在工業(yè)流化床密相區(qū)埋管設(shè)計中使用起來更加便捷高效。前人基于各自實驗范圍提出了形式各異的經(jīng)驗關(guān)聯(lián)式。
GENETTI等[23]以玻璃珠為床料,提出了對大粒徑范圍(114~470 μm)適用性較好的埋管傳熱經(jīng)驗關(guān)聯(lián)式:
式中:Nu為努塞爾數(shù);ε為床層孔隙率;ρg為氣體密度,kg/m3;dp為床料顆粒直徑,m;μg為氣體黏度,Pa·s;uo為流化風速,m/s;d0為埋管外徑,m。
ANDEEN[24]在GENETTI等[23]的基礎(chǔ)上,擴大經(jīng)驗關(guān)聯(lián)式粒徑適用范圍(360~710 μm),提出了高流化速下具有較高精度的埋管傳熱經(jīng)驗關(guān)聯(lián)式:
式中:Pr為普朗特數(shù);g為重力加速度,取 9.81 m2/s;ρp為床料顆粒密度,kg/m3。
GREWAL等[25]所提出的經(jīng)驗公式應(yīng)用范圍較廣,在熱態(tài)(400~900 ℃)常壓和冷態(tài)高壓(0.1~ 1.2 MPa)以及鈦鐵礦作為床料下對埋管傳熱系數(shù)的預(yù)測精度進行實驗驗證,經(jīng)驗關(guān)聯(lián)式如下:
式中:cpb為床料比熱容,J/(kg·K);kg為氣體導熱系數(shù),W/(m·K)。
采用上述3種換熱關(guān)聯(lián)式預(yù)測本文以鈦鐵礦為床料的熱態(tài)增壓流化床埋管傳熱系數(shù),實驗結(jié)果與預(yù)測結(jié)果對比見圖7。從圖7 可以看出:上述換熱關(guān)聯(lián)式預(yù)測值與本文實驗值均存在較大差異。這可能是本文實驗工況超出上述關(guān)聯(lián)式的應(yīng)用范圍。對比上述關(guān)聯(lián)式,可以發(fā)現(xiàn)盡管均考慮了埋管直徑對傳熱的影響,但不同關(guān)聯(lián)式認為埋管直徑對傳熱的影響程度不同。從文獻[23-25]可見:管徑變化對于傳熱系數(shù)變化的影響逐漸減小。此外,與GENETTI等[23]和ANDEEN[24]所提出的經(jīng)驗公式相比,GREWAL等[25]所提出的經(jīng)驗公式還額外考慮了顆粒比熱不同帶來的影響。因此GREWAL等[25]所提出的公式有最高的準確性。但對于以鈦鐵礦為床料的埋管表面?zhèn)鳠嵯禂?shù)的預(yù)測,流化介質(zhì)和顆粒導熱系數(shù)對于傳熱系數(shù)的影響便不可忽略。
圖7 前人經(jīng)驗關(guān)聯(lián)式對以鈦鐵礦為床料的埋管傳熱系數(shù)的適用性Fig.7 Applicability of empirical correlation to heat transfer coefficient with ilmenite as bed material of immersed tube
本文基于準確性相對較高的GREWAL等[25]關(guān)聯(lián)式,結(jié)合導熱系數(shù)的影響,擬合出了新的增壓流化床水平埋管換熱關(guān)聯(lián)式:
式中:kp為床料顆粒導熱系數(shù),W/(m·K)。
圖8所示為新經(jīng)驗關(guān)聯(lián)式與本文以及文獻中的實驗結(jié)果的對比。從圖8 可見:94.2%數(shù)據(jù)的預(yù)測結(jié)果相對誤差范圍在25%以內(nèi)。流化床增壓后,流化介質(zhì)物性及眾多流態(tài)化參數(shù)的急劇變化導致大部分傳熱模型適用性較低,相對誤差較大。本文提出的經(jīng)驗?zāi)P拖鄬φ`差范圍在現(xiàn)有模型中較小,具有較高的精度,可用于以鈦鐵礦為床料的流化床密相區(qū)水平埋管設(shè)計。
圖8 以鈦鐵礦為床料的埋管傳熱系數(shù)實驗值與新經(jīng)驗關(guān)聯(lián)式預(yù)測值對比Fig.8 Comparison between experimental value and predicted value of new empirical correlation of immersed tube
1) 隨著爐膛壓力和床層溫度增加,埋管平均傳熱系數(shù)逐漸變大,而隨著流化風速增加,傳熱系數(shù)呈現(xiàn)出先變大后減小的趨勢。
2) 鈦鐵礦整體傳熱效果優(yōu)于石英砂,主要由于鈦鐵礦較高的密度導致其有著更大的體積熱容,與埋管接觸時可以傳遞更多的熱量。
3) 鈦鐵礦埋管周向傳熱系數(shù)的分布規(guī)律與石英砂的周向傳熱系數(shù)分布規(guī)律保持一致,說明鈦鐵礦作為床料并未改變顆粒在埋管周向的分布。
4) 本文在考慮導熱系數(shù)影響下提出了新的經(jīng)驗關(guān)聯(lián)式,其對以鈦鐵礦為床料的埋管傳熱系數(shù)預(yù)測精度相對較高。