張永學, 于丹丹, 胡承云, 王元華, 孫 東

(1.中國石油大學機械與儲運工程學院,北京 102249; 2.華東理工大學化學工程聯(lián)合國家重點實驗室,上海 200237;3.中國石化勝利油田技術檢測中心,山東東營 257000)

基于響應面法的油田蓄熱式加熱爐燃燒器優(yōu)化設計

張永學1, 于丹丹1, 胡承云2, 王元華2, 孫 東3

(1.中國石油大學機械與儲運工程學院,北京 102249; 2.華東理工大學化學工程聯(lián)合國家重點實驗室,上海 200237;3.中國石化勝利油田技術檢測中心,山東東營 257000)

為提高水套加熱爐熱效率,降低污染物排放濃度,設計適用于油田水套加熱爐的蓄熱式燃燒器,并確定最優(yōu)結構參數(shù)?；陧憫娣ú捎萌蛩匚逅降闹行膹秃显囼炘O計和嶺嵴分析方法,通過數(shù)值模擬分析加熱爐燃燒器燃料氣噴口與二次空氣噴口之間相對高度、二者之間徑向距離和二次空氣大直徑噴口間夾角三因素對輻射管內NO濃度的影響機制,確定燃燒器最優(yōu)結構參數(shù)。結果表明:二次空氣大直徑噴口間夾角的變化對響應值的影響最為顯著,較小的夾角能夠獲得更均勻的溫度分布和更低的NO排放濃度(體積分數(shù));采用優(yōu)化后燃燒器的加熱爐NO濃度明顯降低,理論熱效率也有一定提高。數(shù)值模擬結果與建立的響應函數(shù)模型預測值基本一致,采用響應面法優(yōu)化水套加熱爐蓄熱式燃燒器結構的方法是可行的。

響應面法; 數(shù)值模擬; 結構優(yōu)化; 蓄熱式燃燒器; NOx排放

在原油和天然氣集輸系統(tǒng)中水套加熱爐廣泛應用于油田聯(lián)合站。長期以來,在役的水套加熱爐熱量利用效率低,NOx等污染物排放量大,不利于節(jié)能減排[1];為提高熱效率,多采用設置空氣預熱器回收余熱的方式,但依靠回收煙氣熱量提高加熱爐效率效果不明顯。新型高效蓄熱式燃燒(high temperature air combustion, HTAC)技術是20世紀90年代初發(fā)展起來的一項新型燃燒技術,在發(fā)達國家得到了普遍的應用和推廣。在提高熱效率的同時大幅節(jié)能,且大量降低CO2、NOx排放,在國際燃燒界被譽為最具發(fā)展?jié)摿Φ男滦腿紵夹g之一[2]。HTAC目前在鋼鐵、玻璃等行業(yè)的應用較為成熟,改造后的工業(yè)爐在達到設計產量的同時節(jié)能效果很顯著[3];但到目前為止該技術在石油領域的應用不多。隨著計算機技術和燃燒學的飛速發(fā)展,數(shù)值模擬方法作為研究高溫復雜燃燒過程的有效手段被廣泛應用到工程中,其為改善燃燒過程,提高加熱爐熱效率,減小煙氣污染物排放等提供理論依據(jù)[4]。Liu 等[5]使用 Fluent 軟件研究了某燃燒器噴嘴出口位置及噴嘴直徑對NOx排放的影響,根據(jù)結果選取了優(yōu)化尺寸。Su 等[6]和陸曉峰等[7]也運用了CFD手段對燃燒器進行了結構優(yōu)化設計。江華等[8]采用數(shù)值模擬方法分析了結構變化對燃燒器性能的影響。響應面法(response surface methodology,RSM)是數(shù)學方法和統(tǒng)計方法結合的產物,用來對受多個變量影響的響應值進行建模和分析,從而優(yōu)化響應值。許多學者對響應面法進行了研究,在生物學、醫(yī)學、工程學領域都得到了廣泛應用[9]。筆者以某油田聯(lián)合站800 kW水套加熱爐為研究對象,對燃燒器進行蓄熱式改造,使用Statistica軟件對新型燃燒器的數(shù)值模擬結果進行響應面分析,研究燃氣噴口與二次空氣噴口的相對高度、徑向距離和二次空氣大直徑噴口間夾角對蓄熱式水套加熱爐NOx排放濃度的影響。

1 數(shù)值模擬模型

1.1 幾何模型及網(wǎng)格劃分

數(shù)值模擬的幾何模型為油田某蓄熱式燃氣水套加熱爐,設計工藝負荷為0.8 MW,爐內輻射管為U型管,長度為6.6 m,管內徑為600 mm(圖1)。改造后的加熱爐在輻射管兩端安裝蓄熱式燃燒器,以一定的頻率切換進氣方向,用以吸收煙氣余熱,預熱助燃空氣。該蓄熱式燃燒器為空氣分級燃燒器,一次助燃空氣為常溫,二次助燃空氣經過燃燒器內部蓄熱體加熱,以高速噴入爐膛;6個一次空氣噴口圍繞燃料噴口均勻布置,二次空氣噴口分為5個大噴口與2個小噴口,在外圍呈偏心分布(圖2)。燃氣噴口與二次空氣噴口之間相對高度h、徑向距離L和二次空氣大直徑噴口之間夾角α為燃燒器待優(yōu)化結構參數(shù)。

圖1 原水套爐爐管結構Fig.1 Fire tube structure of original water jacket furnace

模擬計算區(qū)域為蓄熱式燃燒器(暫不考慮蓄熱體)與U形輻射管。采用混合結構網(wǎng)格處理結構復雜的燃燒器，由于噴嘴和火道等區(qū)域各參數(shù)的梯度較大，故進行局部加密[10]。爐膛內結構簡單，為保證計算精度、節(jié)省計算時間，采用六面體網(wǎng)格。燃燒器模型網(wǎng)格劃分見圖3(a)，整體U形輻射管網(wǎng)格劃分見圖3(b)。最終得到的網(wǎng)格數(shù)為210萬。

圖2 燃燒器結構待優(yōu)化參數(shù)Fig.2 Optimizing parameters of burner structure

圖3 數(shù)值模擬區(qū)域網(wǎng)格劃分Fig.3 Mesh generation of calculation region

1.2 數(shù)學模型及邊界條件

綜合考慮工程數(shù)值模擬精度和計算工作量,輻射管內的湍流計算基于雷諾平均控制方程組,選擇標準k-ε湍流模型和標準壁面函數(shù)。該蓄熱式燃燒器為非預混燃燒器,輻射管內的燃燒為擴散燃燒,因此燃燒模型采用混合分數(shù)PDF模型,爐內輻射選擇離散坐標DO模型。該加熱爐所采用的燃氣不含氮化物,因此NOx的生成選擇熱力型和快速型生成模型,同時由于燃燒過程中NO生成量占NOx生成量的95%以上[11],所以后續(xù)計算只考慮NO的排放。

燃氣和空氣入口選用流量入口邊界條件,出口采用壓力出口邊界條件。燃氣入口流量為0.019 kg/s,溫度為293 K;設置一次、二次空氣入口分級比為1∶9,一次空氣入口流量為0.029 1 kg/s,溫度為293 K,二次空氣入口流量為0.261 9 kg/s,溫度為1 073 K(通過蓄熱體換熱預熱),過?？諝庀禂?shù)為1.05。U型管壁面采用定溫壁面條件。出口壓力為-800 Pa。燃料氣組分主要為CH4,占81.85%;其他烷烴氣體占9.8%,N2和O2體積分數(shù)分別為4.21%和 0.72%。

2 響應面法應用

2.1 響應面理論

在合理的試驗基礎上通過構造一個多元二次回歸方程用于擬合因素與響應值的函數(shù)關系,通過分析該回歸方程,得到最優(yōu)參數(shù)組合,從而解決多因素問題[12-13]。響應值Y與因素X1,X2,…,Xk的函數(shù)關系表示為

Y=f(X1,X2,…,Xk)+ε.

(1)

式中,ε為響應值Y的觀測誤差和隨機誤差;f為響應面,用來表述函數(shù)與因素的關系。

引入一階線性項、二階交互項和平方項,最終得到響應面法的逼近函數(shù)式為

(2)

式中,β0為常數(shù)項;βi為各線性項的系數(shù);βij為交互項的系數(shù);βii為平方項的系數(shù)[14]。

2.2CCD中心復合設計

為了確定式(2)中的所有系數(shù),數(shù)值試驗設計必須保證每個因素至少為三水平。CCD中心復合設計是常用的二階對稱設計中的一種,是響應面法最常用的試驗設計法。以三因素試驗為例,CCD包括:①完全析因設計或者部分析因設計,試驗點位于立方體頂點,稱為立方點;②附加設計,通常是星點設計,試驗點稱為軸向點,位于中心軸線上;③中心點,位于立方體中心[15]。三因素CCD試驗設計點分布見圖4。

圖4 三因素中心復合設計Fig.4 Central composite design for three variables

2.3 試驗設計及響應模擬結果

為了考查燃料噴口與二次空氣噴口相對高度h、二者徑向距離L和二次空氣噴口夾角α的組合對NO排放濃度(體積分數(shù))的影響,對以上三因素進行CCD中心復合設計。對h、L、α三因素進行編碼,編碼水平見表1(h、L的單位為mm,α單位為(°))。

表1 響應面試驗因素和水平Table 1 Experimental factors and levels design

根據(jù)CCD原理,設計了三因素五水平共17個試驗點的響應面分析試驗。其中14個為析因子(立方點與軸向點),3個為零點(中心點)。零點試驗重復3次,用以估計試驗誤差[16]。試驗設計表及模擬結果見表2。

表2 CCD設計表及模擬結果Table 2 CCD design and response results

3 結果及其討論

3.1 回歸模型分析

用Statistica軟件對所有數(shù)據(jù)進行回歸擬合,得到以NO排放濃度為響應值的回歸方程為

Y=24.26+1.97h-2.45h2+2.77L-0.93L2+9.26α+2.16α2-3.34hL-2.87hα+12.84Lα.

(3)

回歸方程的擬合結果見表3。

表3 回歸方程的擬合分析Table 3 Fitting analysis of response surface quadratic model

方差分析結果見表4。如表4所示,模型擬合系數(shù)R2=0.814 5,表明模型中81.45%以上的數(shù)據(jù)變異性均能通過模型進行解釋,說明模擬試驗是可靠的。一般認為,當某項P<0.05時,證明該項指標影響顯著,若某項P<0.01,則表明該項指標的影響極為顯著。模型的P值小于0.001,達到極為顯著的標準,再次證明模擬試驗的可靠性。模型失擬項不顯著(P>0.05),表明回歸模型正確。

表4 回歸方程的方差分析Table 4 Variance analysis for response surface quadratic model

3.2 各因素對NOx生成量的響應面分析

由表3、 4可見,各因素中,L、α乘積項的P值小于0.01,為極顯著的模型項,說明L、α存在明顯的交互作用。α線性項小于0.05,為顯著的模型項。對比P值,各因素對NO排放濃度的影響程度由大到小依次為α>L>h。

為了對結果有更清晰的理解,通過Statistica軟件繪制各因素對NO濃度產生影響的響應面3D圖,如圖5所示。由圖5(a)、(b)可見,當h>0時,響應值變化平緩,梯度值較小,此時主要是h的平方項發(fā)揮作用。當h<0時,響應值梯度增大,說明h的線性項在此區(qū)間占主導地位。同時,當h為定值時,響應值沿L和α方向變化梯度明顯,尤其沿α方向,此時h對響應值影響可忽略不計。由圖5(a)、5(c)可見,當L>0時,響應值變化平緩,梯度值較小,此時主要是L的平方項發(fā)揮作用。當L<0時,響應值梯度增大,說明L的線性項在此區(qū)間占主導地位,同時,在L為定值時,響應值沿h和α方向變化梯度值較大,尤其是沿α方向。由圖5(b)、5(c)可見,響應值沿α方向變化梯度最明顯,說明α對響應值的影響占主導地位,同時由圖5(c)可知,響應值在α方向和L方向上梯度均很大,再次證明α和L存在明顯的交互作用。

圖5 響應面分析圖Fig.5 Analysis of response surface plots

以NO排放濃度為指標,三因素中α取值的變化對指標影響最為明顯。分析認為,較小的α值偏流卷吸的效果更加明顯,降低了燃燒區(qū)域的氧濃度,延緩了預熱后的二次助燃空氣與燃料氣的混合,使燃燒區(qū)位置升高,范圍擴大,局部高溫區(qū)減小,更有助于燃燒火焰的彌散;而較大的α值則使預熱空氣與燃氣的混合區(qū)域降低,燃燒環(huán)境氧含量更高,燃燒速率和強度加大,局部高溫明顯,從而為NOx的生成創(chuàng)造了條件。這一解釋也可以從分布云圖上得到驗證。圖6為α值對溫度分布的影響。

圖6 α值對溫度分布的影響Fig.6 Effect of α on temperature distribution

在其余因素不變的條件下,圖6(a)為CCD第9試驗點的溫度云圖,α編碼尺寸為1;圖6(b)為CCD第4試驗點的溫度云圖,α編碼尺寸為-1。相對于圖6(a),圖6(b)在整個燃燒區(qū)域溫度分布更加彌散,局部高溫區(qū)不明顯。

3.3 嶺嵴分析及驗證

由響應面分析結果可看出,二次響應面是鞍面,沒有唯一最佳值,因此使用SAS軟件嶺嵴分析,以進一步確定燃燒器結構的最優(yōu)值。嶺嵴分析是以原始設計中心點為球心(本試驗取球心小于0,等于0,大于0),在以R為半徑的超球面與響應面的交點(即嵴點)形成的軌跡范圍內找出最佳響應值。嶺嵴分析的結果是對每個坐標從球心開始不斷擴大,半徑R不能超過試驗范圍,本文中選取R在0.1、0.2、…、1.0范圍內計算嵴點。當R>0.8時,標準誤差較大,因此不予考慮。在本試驗水平范圍,最小響應值時R=0.8,此時燃料、二次空氣噴口相對高度h為96 mm,徑向距離L為130 mm,二次空氣噴口夾角α為29°。

按照上述優(yōu)化結構參數(shù)做驗證模擬計算,NO濃度實際模擬結果為7.24×10-6,與預測值7.18×10-6基本一致,驗證了RSM分析的正確性。

圖7為沿中心軸線方向3種結構的NO體積分數(shù)分布對比,三種結構分別為優(yōu)化結構、原始結構(CCD零點)和最差結構(CCD 第 2 試驗點,(h,L,α)編碼尺寸為(-1,1,1))。由圖7可見,各結構燃燒器NO體積分數(shù)在1 ～ 2 m區(qū)域達到峰值。第2組結構輻射管的出口NO體積分數(shù)最高,接近60×10-6,而優(yōu)化結構輻射管出口NO體積分數(shù)僅約為7×10-6。這說明經過 RSM分析和嶺嵴分析得到的優(yōu)化結構有效控制了NOx的生成,實現(xiàn)了 NOx減排效果,達到了預期目的。

圖7 NO濃度云圖對比Fig.7 Comparison of NO concentration distributions

3.4 節(jié)能效果評價

在20 ℃的環(huán)境溫度下,現(xiàn)役水套加熱爐的排煙溫度為167.5 ℃,熱效率為79.70%,其中不完全燃燒損失為3.12%,排煙損失為14.28%,散熱損失為2.9%。同樣條件下優(yōu)化改造后的蓄熱式加熱爐的排煙溫度約為150 ℃,排煙損失為5.43%;采用改造后燃燒器的加熱爐內燃燒過程更加緩慢均勻,燃料能夠完全燃燒,因此不完全燃燒損失為0%;散熱損失為2.9%。計算得到優(yōu)化后的加熱爐理論熱效率為91.67%,比現(xiàn)役加熱爐的熱效率提高了11.97%。

4 結 論

(1)運用響應面分析方法得到的回歸模型中81.45%以上的數(shù)據(jù)變異性均能通過模型進行解釋,模擬試驗有效可靠。

(2)以NO排放濃度為響應值,燃料氣噴口與二次空氣噴口相對高度h、徑向距離L和二次空氣大噴口之間夾角α三個因素對響應值的影響程度依次為α>L>h,其中,α和L存在明顯的交互作用。

(3)基于嶺嵴分析得到燃燒器的優(yōu)化結構:燃料、二次空氣噴口相對高度h為96 mm,徑向距離L為130 mm,二次空氣噴口夾角α為29°。采用最優(yōu)結構燃燒器的輻射管出口NO體積分數(shù)為7.24×10-6,顯著降低了NOx排放;且該模擬結果與預測值吻合良好。

(4)經優(yōu)化改造后的蓄熱式水套加熱爐的理論熱效率為91.67%,比現(xiàn)役加熱爐的熱效率提高了11.97%。

[1] 唐建峰,李玉星,薛吉明,等. 水套爐的參數(shù)優(yōu)化及自耗氣降低技術研究[J]. 中國石油大學學報(自然科學版), 2007,31(5):74-77. TANG Jianfeng, LI Yuxing, XUE Jiming, et al. Research of optimal parameter of fluid-jacket heater and depression of gas consumption[J].Journal of China University of Petroleum(Edition of Natural Science), 2007,31(5):74-77.

[2] 溫治,代朝紅.蓄熱式高溫空氣燃燒技術的研究現(xiàn)狀及應用前景分析[J].河南冶金, 2002(6):3-8,37. WEN Zhi, DAI Zhaohong. Research status and application prospect of high temperature air combustion technology with regenerator [J]. Henan Metallurgy,2002 (6):3-8,37.

[3] 謝國威,于宏,秦文，等. 蓄熱式加熱爐應用與熱態(tài)模擬[J]. 冶金能源,2005(3):30-31,38. XIE Guowei, YU Hong, QIN Wen, et al. Application and thermal simulation of regenerative heating furnace[J]. Energy for Metallurgical Industry, 2005(3):30-31,38.

[4] 歐儉平,吳道洪,肖澤強.蓄熱式加熱爐內流體流動燃燒與傳熱的數(shù)值模擬[J]. 工業(yè)爐,2003,5(1):52-54. OU Jianping, WU Daohong, XIAO Zeqiang. Numerical simulation on the flow, combustion and heat transfer in regenerative reheating furnace [J].Industrial Furnace,2003, 5(1):52-54.

[5] LIU F G, YOU X Y, WANG Q,et al. On optimal design and experimental validation of household appliance burner of low pollutant emission [J]. Energy Conversion and Management, 2013,76:837-845.

[6] SU Y, DENG W, JIANG F, et al. Numerical study of effect of burner jet parameters on high temperature air combustion of coal gas: International Conference on Computer Distributed Control and Intelligent Environmental Monitoring [C].IEEE Computer Society, 2011.

[7] 陸曉峰,李元青. 基于CFD計算的燃燒器結構改進研究[J].石油學報 (石油加工), 2011,27(5):787-791. LU Xiaofeng, LI Yuanqing. Study on structure improvement of the burner based on CFD [J]. Acta Petrolei Sinica (Petroleum Processing Section), 2011,27(5):787-791.

[8] 江華,毛羽,吳德飛,等.結構變化對氣體燃燒器性能影響的數(shù)值模擬[J].中國石油大學學報(自然科學版), 2004,28(6):94-98,102. JIANG Hua, MAO Yu, WU Defei,et al. Numerical simulation of influence of different structureson the combustion of gas burner[J]. Journal of China University of Petroleum(Edition of Natural Science),2004,28(6):94-98,102.

[9] MONTGOMERY D C. Design and analysis of experiments [M]. 8th ed. New York: John Wiley & Sons Inc, 2012.

[11] HILL S C, SMOOT L D. Modeling of nitrogen oxides formation and destruction in combustion systems [J]. Progress in Energy and Combustion Science,2000, 26(4/5/6):417-458.

[12] 隋述會,洪定一,洪東峰.基于響應曲面法的環(huán)管裝置聚丙烯相對分子質量的模擬及優(yōu)化[J]. 石油學報(石油加工), 2013,29(6):1054-1061. SUI Shuhui, HONG Dingyi,HONG Dongfeng. Optimization and simulation of relative molecular mass of polypropylene produced in a loop reactor by response surface methodology[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2013,29(6):1054-1061.

[13] 吳曉磊. 新型低氮燃氣分級燃燒器CFD模擬和熱態(tài)試驗研究[D]. 上海: 華東理工大學, 2015. WU Xiaolei. CFD Research and hot performance test of a new reburning burner with low NOxemission[D]. Shanghai: East China University of Science and Technology, 2015.

[14] 王永菲,王成國.響應面法的理論與應用[J].中央民族大學學報(自然科學版), 2005,14(3):236-240. WANG Yongfei, WANG Chengguo. The application of response surface methodology [J].Journal of the Central University for Nationalities(Natural Science Edition),2005,14(3):236-240.

[15] BEZERRA M A, SANTELLI R E, OLIVEIRA E P, et al. Response surface methodology (RSM) as a tool for optimization in analytical chemistry[J]. Talanta, 2008,76(5):965-977.

[16] OBENG D P, MORREL L S, NAPIER-MUNN T J. Application of central compositer rotatable design to modeling the effect of some operating variables on the performance of the three-product cyclone [J]. International Journal of Mineral Processing, 2005,76(3):181-192.

(編輯 沈玉英)

Optimization design of oilfield regenerative heating furnace burner based on response surface method

ZHANG Yongxue1, YU Dandan1, HU Chengyun2, WANG Yuanhua2, SUN Dong3

(1.CollegeofMachineryandTransportationEngineeringinChinaUniversityofPetroleum,Beijing102249,China;2.State-KeyLaboratoryofChemicalEngineering,EastChinaUniversityofScienceandTechnology,Shanghai200237,China;3.TechnologyTestingCenterinShengliOilfield,SINOPEC,Dongying257000,China)

In order to improve the thermal efficiency of the jacket heater and reduce the pollutant concentration, it is aimed to design a regenerative burner which is suitable for water jacket furnaces, and the purpose is also to determine the optimal structural parameters. Based on the response surface method (RSM) with the central composite design (CCD) of 3 factors and 5 levels and in combination with the ridge analysis method, the numerical simulation was carried out. The relative height, radial distance between the fuel nozzle and secondary air nozzle and the spray angle of secondary air nozzles on the NO emission of the burner were analyzed, then the optimal structural parameters can be obtained. The result shows that the variation of the angle of secondary air nozzles has the most significant effect on the responsive values, and the smaller angle can reach more homogeneous temperature distribution and lower NO emission concentration (volume fraction). It is also found that the NO emission concentration of the burner has been reduced greatly after optimization, and the theoretical efficiency is also enhanced. It indicates that the numerical simulation is consistent with the RMS model result, and the route to optimize the structures of the regenerative burner in the jacket heater is possible to realize.

response surface method; numerical simulation; structural optimization; regenerative burner; NOxemission

2016-04-12

張永學(1977-)，男，教授，博士，博士生導師，研究方向為熱能及流體計算。E-mail:zhyx@cup.edu.cn。

于丹丹(1990-)，女，碩士研究生，研究方向為熱力過程優(yōu)化與節(jié)能。E-mail: yudandan613@163.com。

1673-5005(2017)02-0156-07

10.3969/j.issn.1673-5005.2017.02.019

TK 175

A

張永學,于丹丹,胡承云,等.基于響應面法的油田蓄熱式加熱爐燃燒器優(yōu)化設計[J]. 中國石油大學學報(自然科學版),2017,41(2):156-162.

ZHANG Yongxue, YU Dandan, HU Chengyun, et al. Optimization design of oilfield regenerative heating furnace burner based on response surface method[J]. Journal of China University of Petroleum (Edition of Natural Sciences), 2017,41(2):156-162.