儲罐池火特性研究進(jìn)展*

張日鵬，趙祥迪，王 正，袁紀(jì)武，馬浩然

(中國石化青島安全工程研究院化學(xué)品安全控制國家重點實驗室，山東青島 266071)

隨著石油化工行業(yè)的蓬勃發(fā)展,儲罐的數(shù)量和體積急劇增加，儲罐內(nèi)的燃料都具有易燃易爆的特性，而明火、雷擊、靜電和硫化亞鐵自燃都有可能成為火源[1]，因此，儲罐火災(zāi)事故頻頻發(fā)生；而在眾多的儲罐火災(zāi)事故當(dāng)中，儲罐池火災(zāi)事故又占了很大比重。為了能更好地對儲罐區(qū)進(jìn)行安全評估、制定罐區(qū)池火災(zāi)事故應(yīng)急預(yù)案、確定儲罐之間的安全距離以及確保滅火作業(yè)時人員安全和設(shè)備的可操作性，需對儲罐池火的特性進(jìn)行研究。

國內(nèi)外專家對儲罐池火特性進(jìn)行了大量研究，主要集中在燃燒速率、火焰形態(tài)(包括火焰高度、火焰傾斜和拖曳)、火焰脈動頻率、熱輻射通量[2]等參數(shù)。本文對儲罐池火特性參數(shù)的半經(jīng)驗公式進(jìn)行總結(jié)和對比研究。

1 燃燒速率

儲罐油品的燃燒由油品蒸發(fā)和蒸氣燃燒兩部分組成[3]。油品燃燒，釋放出熱量，未燃油品通過輻射或?qū)α鲝幕鹧嫖諢崃浚ㄟ^傳導(dǎo)與儲罐進(jìn)行熱量交換；開始階段未燃油品吸收的熱量大于損失的熱量，未燃油品溫度升高并蒸發(fā)加劇，燃燒速率逐漸增加；一段時間后，未燃油品吸收的熱量等于損失的熱量和蒸發(fā)需要的熱量，蒸發(fā)速率趨于穩(wěn)定，此時為穩(wěn)定燃燒階段，燃燒速率基本保持不變；當(dāng)燃料產(chǎn)生的熱量不足以維持該燃燒狀態(tài)進(jìn)行下去時，燃燒進(jìn)入衰減階段，燃燒速率和溫度逐漸下降。

Blinov和Khudiakov通過實驗得出：在D<0.03 m時，燃料為層流燃燒，隨著儲罐直徑的增加，燃燒速率下降；在D>1 m時,燃料為湍流燃燒，燃燒速率是常數(shù)，與儲罐直徑和燃料類型無關(guān)；在0.03 m

直徑>1 m罐，Burgess等認(rèn)為線性燃燒速率為：

(1)

質(zhì)量燃燒速率m等于線性燃燒速率與密度乘積,如密度不可得，則可通過式(2)估算:

(2)

式中：y——線性燃燒速率，m/s；

ΔH——燃料燃燒熱，kJ/kg；

ΔH*——改進(jìn)的燃料蒸發(fā)吸收熱，kJ/kg。

公式(1)和(2)都是建立在“池火穩(wěn)定燃燒時燃燒速率等于蒸發(fā)速率”的基礎(chǔ)之上，然后通過實驗數(shù)據(jù)確定前面系數(shù)得出；其中，公式(1)與實驗結(jié)果吻合得更好[5,6]。

如果燃料特性參數(shù)未知并且要考慮直徑對燃燒速率的影響，可以使用Burgess和Zabetakis提出的經(jīng)驗公式來計算燃燒速率：

m=m∞(1-e-kD)

(3)

式中：m——質(zhì)量燃燒速率，kg/(m2·s)；

m∞——無限大池的質(zhì)量燃燒速率，kg/(m2·s)；

D——池火直徑，m；

k——常數(shù)，m-1。

許多學(xué)者通過實驗確定了不同燃料的m∞和k值[5,6]。對于汽油，使用Babrauskas推薦參數(shù)，做出質(zhì)量燃燒速率隨儲罐直徑變化的關(guān)系圖，見圖1。

圖1 汽油燃燒速率隨直徑的變化規(guī)律

從圖1可以看出，隨著儲罐直徑的增加，質(zhì)量燃燒速率先逐漸增加，然后維持不變(其最大值與公式(2)計算結(jié)果一致[5])，該變化規(guī)律與池火湍流燃燒時燃燒速率隨直徑變化規(guī)律相一致。因此，公式(3)適用于直徑大于1 m池火湍流穩(wěn)定燃燒時，質(zhì)量燃燒速率的計算；公式(2)適用于儲罐直徑大于2 m時穩(wěn)定燃燒的質(zhì)量燃燒速率計算。

關(guān)于風(fēng)速對燃燒速率的影響，由于風(fēng)對燃燒過程的影響極為復(fù)雜，因此，大多數(shù)熱輻射模型在估算燃燒速率時都沒有考慮風(fēng)的影響，而用風(fēng)對火焰高度的影響來彌補(bǔ)[5,9]。

2 火焰形態(tài)

目前，一般認(rèn)為儲罐池火火焰為圓柱體?；鹧嫘螒B(tài)由火焰直徑、火焰高度、火焰傾斜和拖曳構(gòu)成。

2.1 火焰高度

燃燒過程中火焰高度變化很大，一般提到的火焰高度均指平均火焰高度[11]。關(guān)于火焰高度的定義有很多，本文的火焰高度采用Zukoski利用間歇率提出的平均火焰高度[4]。

2.1.1不考慮風(fēng)速

Heskestad[12,13]、Brotz[14]在池火災(zāi)實驗基礎(chǔ)上以及Thomas[15]在木垛實驗和量綱分析的基礎(chǔ)上提出不考慮風(fēng)速影響的火焰高度計算公式(4)～(6)：

(4)

(5)

(6)

式中：H——平均火焰高度，m；

Q——總放熱量，kW；

ρa(bǔ)——空氣密度，kg/m3。

利用上述公式計算輕柴油池火火焰高度，做出火焰高度與直徑比(H/D)隨直徑變化的曲線，見圖2。

圖2 H/D隨儲罐直徑變化規(guī)律

從圖2可看出，隨著直徑的增加，Brotz公式計算的H/D值先減小然后維持不變；Heskestad和Thomas公式計算的H/D逐漸減小。

筆者曾對11.5 m和20.3 m油池進(jìn)行柴油池火實驗，測得火焰高度與直徑的比值分別為1.5和1.4，結(jié)合文獻(xiàn)[3]中池火高度數(shù)據(jù)，可以得出：Thomas公式與實驗結(jié)果更符合；同時,在點源輻射模型和Shokri-Beyler輻射模型中多使用Heskestad公式，而在表面輻射模型中多使用Thomas公式；表面輻射模型更廣泛地應(yīng)用于熱輻射計算，這也說明了Thomas公式應(yīng)用更為普遍。

2.1.2考慮風(fēng)速

Thomas提出有風(fēng)條件下火焰長度(L)修正公式：

(7)

Moorhouse[17]基于大尺度LNG池火實驗提出了考慮風(fēng)速的火焰長度公式：

(8)

考慮風(fēng)作用時，輕柴油池火火焰長度隨儲罐直徑的變化規(guī)律如圖3。

圖3 考慮風(fēng)作用時L/D值隨儲罐直徑變化規(guī)律

從圖3可以看出，L/D值隨著直徑的增加先稍有增加然后逐漸減小，最后趨于定值；針對Thomas公式，考慮風(fēng)作用時計算得出的火焰長度低于不考慮風(fēng)作用時，并且風(fēng)速越大，火焰長度越??；不同風(fēng)速時，L/D值隨直徑增加而減小的速率基本相等;針對Moorhouse公式，風(fēng)速越大，火焰長度越??；不同風(fēng)速時，L/D值隨直徑增加而減小的速率基本相等，但減小速率小于Thomas公式的減小速率。

考慮風(fēng)作用時，Thomas公式和Moorhouse公式應(yīng)用都較廣泛，但Thomas公式應(yīng)用相對更多。

2.2 火焰傾斜

在風(fēng)的作用下，池火火焰將會傾斜。1966年，Sliepcevich[18]基于小尺度池火實驗提出了火焰傾角與弗勞德數(shù)Fr以及雷諾數(shù)Re的關(guān)系式：

(9)

同年，Sliepcevich基于實驗數(shù)據(jù)分析得出：a′=3.3,b′=0.8,c′=0.07；1992年，Blinding基于大尺寸LNG池火實驗數(shù)據(jù)分析得出：a′=0.666,b′=0.333,c′=0.117。

Thomas也提出了計算火焰傾角的表達(dá)式：

cosθ=0.7×(u*)-0.49

(10)

AGA(American Gas Association)[19]針對火焰傾斜提出了計算公式：

cosθ=1u*≤1

(11)

對于儲罐直徑20 m的汽油池火，計算風(fēng)速從1到10 m/s時，火焰傾角變化規(guī)律，如圖4所示。

圖4 火焰傾角隨風(fēng)速的變化規(guī)律

由圖4可見，隨著風(fēng)速增加，各公式計算得出的火焰傾角都先增加后趨于穩(wěn)定；其中Thomas公式計算得出的火焰傾角最大；其余3個公式在風(fēng)速大于3 m/s時，計算結(jié)果相近；在風(fēng)速小于3 m/s時，Sliepcevich和Blinding公式認(rèn)為，只要有風(fēng)，火焰就會傾斜，Blinding公式計算得出的火焰傾斜角較大；而AGA公式認(rèn)為，存在一個臨界風(fēng)速，小于該風(fēng)速時，火焰不發(fā)生傾斜。文獻(xiàn)[20]指出，AGA公式計算結(jié)果與實際更符合，應(yīng)用更廣[20]。

2.3 火焰拖曳

由于風(fēng)的作用,火焰基底向下風(fēng)向移動，而上風(fēng)向的火焰邊緣以及火焰寬度均保持不變,這種現(xiàn)象稱為火焰拖曳?，F(xiàn)有的熱輻射計算模型更多的只考慮了火焰傾斜而沒有考慮火焰拖曳。

Moorhouse認(rèn)為火焰發(fā)生拖曳之后，風(fēng)的方向上直徑計算公式如式(12)所示：

(12)

式中：D′——火焰在風(fēng)向上的直徑，m；

θ——火焰傾角，°。

國外專家通過實驗確定了公式(12)中的參數(shù)值，見表1。

表1 公式(12)參數(shù)取值

計算當(dāng)風(fēng)速為5 m/s時，汽油池火的D′/D值隨D的變化規(guī)律，如圖5所示：

圖5 D′/D隨直徑的變化規(guī)律

從圖5可以看出，隨著儲罐直徑的增加，D′/D值逐漸減小，最后維持為定值；宋雪飛等認(rèn)為，Moorhouse計算結(jié)果與實驗更接近，使用廣泛[21,22]。

3 火焰脈動

湍流火焰都有脈動的特性，脈動的主要原因是周圍空氣對火焰的卷吸效應(yīng)[23]，目前常用的脈動頻率計算公式主要有以下幾個。

Cetegen[24]基于無量綱特征數(shù)分析得出脈動頻率計算式(13),Pagni[25]通過分析實驗數(shù)據(jù)得出脈動頻率的經(jīng)驗計算式(14),Byram[9]通過對函數(shù)進(jìn)行解析求解得出火焰脈動周期與直徑的關(guān)系(15):

(13)

(14)

π1=t(gD)0.5

(15)

式中：f——火焰脈動頻率，s-1；

Ta——環(huán)境溫度，K；

π1——火焰脈動周期，s；

t——時間常數(shù)。

可見，火焰脈動頻率與燃料無關(guān)，與D-0.5成正比。

4 熱輻射模型

4.1 點源模型

假定池火火焰可以用一點源表示，該點源位于火焰的幾何中心，并朝所有方向輻射熱量。點源模型熱輻射通量計算公式如(16)所示：

(16)

式中：q——目標(biāo)點接收到的熱輻射通量，kW/m2；

l——點源到目標(biāo)點距離，m；

η——輻射分?jǐn)?shù)；

θ——目標(biāo)物體法線方向與點源和目標(biāo)物體連線的夾角,°。

4.2 Shokri-Beyler模型

Shokri-Beyler模型主要基于實驗數(shù)據(jù)得出：

q=EF

(17)

E=58(10-0.00823D)

(18)

式中：E——火焰表面發(fā)射功率，kW/m2；

F——視角系數(shù)。

4.3 表面輻射模型

表面輻射模型為了便于建模和計算，假設(shè)熱輻射是由火焰表面釋放，事實上，火焰釋放出的熱輻射是由整個火焰(包括燃料氣體、燃燒產(chǎn)物、煙氣顆粒等)產(chǎn)生的，而不僅僅由火焰表面產(chǎn)生，因此表面輻射模型是對復(fù)雜三維熱輻射問題的二維簡化。該模型認(rèn)為火焰是靜止的灰體，火焰之上的煙羽(不可見火焰)輻射也部分考慮。

Mudan首先提出表面輻射模型的計算公式(19)：

q=τFE

(19)

式中：E——火焰表面單位面積的熱輻射釋放速率，kW/m2。

發(fā)射功率計算主要有以下2種方法。

4.3.1定義求解法

發(fā)射功率是單位面積熱輻射釋放速率，采用總熱輻射釋放速率除以總面積，如公式(20)表示：

(20)

4.3.2分層法

火焰分為上下兩部分，下部火焰不被煙氣遮擋，上部火焰存在被煙氣遮擋的情況。Mudan和Croce認(rèn)為火焰發(fā)射功率Em=140 kW/m2，煙氣發(fā)射功率Es=20 kW/m2，通過對汽油、煤油和JP-5進(jìn)行池火實驗，建議采用公式(21)來計算高分子量烴類有煙氣池火的發(fā)射功率(其中S=0.12 m-1)：

Eav=Eme-SD+Es(1-e-SD)

(21)

Binding和Pritchard提出如下公式：

Eav=xlumElum+(1-xlum)Esoot

(22)

式中：xlum——火焰表面中發(fā)光火焰占的比例；

Elum——發(fā)光火焰區(qū)域最大發(fā)射功率，kW/m2；

Esoot——火焰非發(fā)光區(qū)域的發(fā)射功率，kW/m2。

文獻(xiàn)[4]認(rèn)為：對于柴/汽油，Esoot=40 kW/m2；在直徑小于5 m時，Elum隨直徑的增加而增加，關(guān)系式為：Elumgasline=53.64D0.474，Elumdiesel=28.03D0.877，xlum gasoline=0.45,xlum diesel oil=0.30；當(dāng)直徑大于5 m時，Elum=115 kW/m2，xlum隨直徑增加而線性減小，直到池直徑大于20 m時，xlum=0。

利用上述方法，計算不同直徑柴油池火的發(fā)射功率，見圖6。

圖6 發(fā)射功率隨直徑的變化規(guī)律

從圖6可以看出，隨儲罐直徑的增加，定義法和分層法(Binding&Pritchard)計算得出的發(fā)射功率先增加然后減小，分層法(Mudan&Croce)計算得出的發(fā)射功率一直減小。當(dāng)儲罐直徑大于7 m時，三者相差不大。因此，在計算大型儲罐池火災(zāi)的發(fā)射功率時，三者均廣泛應(yīng)用。

不考慮風(fēng)的影響，利用上述三種輻射模型，計算直徑為1 m和20 m的柴油儲罐池火燃燒時，儲罐周圍的熱輻射通量，并作出熱輻射通量與距儲罐中心距離的關(guān)系，見圖7和圖8。

圖7 直徑1 m儲罐周圍的熱輻射通量

從圖7和圖8可以看出，在儲罐直徑為20 m時，對于與儲罐底面處于同一平面的目標(biāo)物體，點源模型、Shokri-Beyler模型和表面輻射模型的計算值非常相近，三種模型可以通用。

圖8 直徑20 m儲罐周圍的熱輻射通量

在儲罐直徑為1 m且目標(biāo)物體距儲罐中心較近時，Shokri-Beyler計算值最大，表面輻射模型計算值居中，點源模型計算值最??；在儲罐直徑為1 m且目標(biāo)物體距儲罐中心較遠(yuǎn)時，Shokri-Beyler計算值依然最大，表面輻射模型和點源模型計算值相近，這說明在儲罐直徑較小時，Shokri-Beyler模型高估了熱輻射通量，儲罐直徑較小且靠近火源時，點源模型由于沒有考慮火焰形狀而低估了熱輻射通量。

5 結(jié)論

a)基于“燃料蒸發(fā)速率等于燃燒速率”的燃燒速率計算公式適合于計算直徑大于2 m，池火穩(wěn)定燃燒時的燃燒速率；考慮直徑影響的燃燒速率計算公式適合于計算池火湍流燃燒的燃燒速率；風(fēng)對燃燒速率的影響還沒有完全研究清楚，在計算熱輻射時一般也沒有考慮，還有待于進(jìn)一步研究。

b)無風(fēng)作用時，Thomas公式更適合用于火焰高度計算；有風(fēng)作用時，Thomas公式和Moorhouse公式都可用來計算火焰高度;AGA公式更適合火焰傾斜角的計算；現(xiàn)有的熱輻射計算中一般不考慮火焰拖曳，火焰拖曳直徑的計算常使用Moorhouse公式。

c)火焰脈動頻率與燃料無關(guān)，與D-0.5成正比。

d)在儲罐直徑為20 m時，對于與儲罐底面處于同一平面的目標(biāo)物體，點源模型、Shokri-Beyler模型和表面輻射模型的計算值非常相近，三種模型可以通用。在儲罐直徑為1 m且目標(biāo)物體距儲罐中心較近時，Shokri-Beyler計算值最大，表面輻射模型計算值居中，點源模型計算值最小；在儲罐直徑為1 m且目標(biāo)物體距儲罐中心較遠(yuǎn)時，Shokri-Beyler計算值依然最大，表面輻射模型和點源模型計算值相近，這說明在儲罐直徑較小時，Shokri-Beyler模型高估了熱輻射通量，儲罐直徑較小且靠近火源時，點源模型由于沒有考慮火焰形狀而低估了熱輻射通量。

e)現(xiàn)有的池火熱輻射模型沒有考慮風(fēng)對燃燒速率的影響、火焰脈動、火焰拖曳以及火焰的三維形態(tài)，因此，建議通過實驗或模擬，進(jìn)一步探究池火在這些方面的規(guī)律，建立更精確的熱輻射模型。