劉偉，李越

中國(guó)礦業(yè)大學(xué)(北京)應(yīng)急管理與安全工程學(xué)院，北京 100083

可燃性油品在裝卸、儲(chǔ)運(yùn)過(guò)程中可能泄漏，當(dāng)油品泄漏到地面后受地形限制，將在某一限定區(qū)域內(nèi)積聚，若遇到火源將引發(fā)油池火災(zāi)，造成人員傷亡和財(cái)產(chǎn)損失[1]。因此，有必要對(duì)油池火開(kāi)展深入全面的研究，以準(zhǔn)確認(rèn)識(shí)其致災(zāi)機(jī)理。

近幾十年來(lái)，油池火災(zāi)被廣泛地研究。在研究?jī)?nèi)容上，主要包括基本燃燒特性、溫度場(chǎng)和熱輻射場(chǎng)?？椎门蟮萚2]為探討不同直徑及初始油層厚度對(duì)沸溢火災(zāi)燃燒速率的影響，進(jìn)行了小尺度沸溢油池火災(zāi)燃燒速率特性實(shí)驗(yàn)研究，并在此基礎(chǔ)上建立了沸溢強(qiáng)度模型。吳向君等[3]在封閉空間開(kāi)展了油池火實(shí)驗(yàn)，研究不同羽流模型的適用性。陳陣等[4]對(duì)油池火經(jīng)典熱輻射模型的適用性進(jìn)行了總結(jié)分析。在研究思路上，主要從實(shí)驗(yàn)環(huán)境和燃油兩個(gè)方面設(shè)置實(shí)驗(yàn)工況。張培紅等[5]進(jìn)行了不同通風(fēng)條件下柴油池火實(shí)驗(yàn)，對(duì)比分析不同通風(fēng)條件下燃燒及引燃特性。張?jiān)平艿萚6]通過(guò)不同尺寸油池火災(zāi)實(shí)驗(yàn)，主要研究了臨界參數(shù)。在研究方式上，主要有全尺寸實(shí)驗(yàn)研究、小尺寸實(shí)驗(yàn)研究和數(shù)值模擬三種。康泉?jiǎng)賉7]開(kāi)展了小尺度油池火實(shí)驗(yàn)，對(duì)其非穩(wěn)態(tài)燃燒特性和熱力學(xué)過(guò)程進(jìn)行了細(xì)致的研究。Sikanen等[8]使用FDS研究了火焰對(duì)液體燃料表面輻射和對(duì)流熱反饋問(wèn)題。周健楠等[9]用FDS對(duì)丙烷儲(chǔ)罐罐內(nèi)池火災(zāi)進(jìn)行數(shù)值模擬，確定其熱輻射分布規(guī)律及煙氣影響范圍，進(jìn)而分析了事故后果。然而，油池火研究還不夠多樣、大尺寸實(shí)驗(yàn)研究不夠充分[10]、火羽流外側(cè)溫度研究還不夠系統(tǒng)、數(shù)值模擬不夠準(zhǔn)確。

針對(duì)以上研究現(xiàn)狀和問(wèn)題，本文選擇低硫原油、高硫原油和重石腦油在室外開(kāi)放空間進(jìn)行不同燃燒面積的自由燃燒實(shí)驗(yàn)，并將實(shí)驗(yàn)得到的燃燒速率自定義輸入FDS中，對(duì)每組實(shí)驗(yàn)進(jìn)行模擬。首先，通過(guò)比較各實(shí)驗(yàn)燃燒速率變化曲線，探究不同種類油料在不同燃燒面積時(shí)各燃燒階段的變化特征與聯(lián)系；其次，依據(jù)油盤正上方溫度測(cè)量數(shù)據(jù)，分析在不同燃燒階段油盤中心線和外側(cè)溫度分布及其變化機(jī)理，并與FDS的模擬結(jié)果進(jìn)行比較；再次，選擇典型的輻射模型對(duì)測(cè)點(diǎn)接收熱輻射通量進(jìn)行計(jì)算，將其與模擬值和實(shí)驗(yàn)值進(jìn)行對(duì)比，分析適用性；最后，對(duì)于大功率油池火災(zāi)，確定其“安全距離-時(shí)間”臨界曲線，得到距火源中心不同位置處人員可以暴露的時(shí)間。通過(guò)實(shí)驗(yàn)和數(shù)值模擬研究，從燃燒速率、火羽流溫度和測(cè)點(diǎn)接收熱輻射通量三個(gè)方面認(rèn)識(shí)油池火災(zāi)的燃燒特性，為科學(xué)防治提供理論基礎(chǔ)。

1 實(shí)驗(yàn)設(shè)置與模擬

1.1 實(shí)驗(yàn)設(shè)置

選擇低硫原油、高硫原油和重石腦油在室外開(kāi)放空間進(jìn)行自由燃燒實(shí)驗(yàn)，實(shí)驗(yàn)編號(hào)和燃料的基本性質(zhì)見(jiàn)表1。

為降低自然風(fēng)對(duì)池火影響，在油盤來(lái)風(fēng)向和垂直于風(fēng)向搭建擋風(fēng)板，燃燒過(guò)程中燃燒區(qū)未出現(xiàn)明顯的風(fēng)流。實(shí)驗(yàn)選用的單個(gè)油盤內(nèi)部長(zhǎng)、寬、高分別為0.841 m、0.595 m、0.13 m，油盤數(shù)量按照2、4、6遞增，油盤面積依次為1 m2、2 m2、3 m2；在其中一個(gè)油盤底部布置天平，記錄質(zhì)量隨時(shí)間的變化；在與油盤上邊沿相同高度水平距火源中心2.5 m、4.5 m、6.5 m處布置輻射熱流計(jì)，記錄熱輻射通量隨時(shí)間的變化；在油盤正上方距地面0.5 m處布置一束熱電偶，間距為0.3 m，對(duì)于Test1—Test5和Test8，熱電偶樹(shù)布置在油盤中心位置，記錄火羽流內(nèi)部溫度隨時(shí)間的變化；對(duì)于Test6和Test7熱電偶樹(shù)布置在油盤外側(cè)，記錄火羽流外部溫度隨時(shí)間的變化。實(shí)驗(yàn)裝置布置如圖1所示。

表1 實(shí)驗(yàn)編號(hào)和油料的基本性質(zhì)

圖1 實(shí)驗(yàn)布置示意圖Fig.1 Schematic diagram of experimental arrangement

1.2 模擬簡(jiǎn)介

本文依據(jù)實(shí)驗(yàn)環(huán)境，確定模擬區(qū)域的長(zhǎng)、寬、高分別為9 m、5 m、5 m，經(jīng)計(jì)算和網(wǎng)格獨(dú)立性驗(yàn)證之后，選用0.15 m×0.15 m×0.15 m的方形網(wǎng)格。網(wǎng)格四周和頂部的邊界條件設(shè)置為開(kāi)放，網(wǎng)格底部設(shè)置為0.5 m厚水泥，環(huán)境溫度為295 K，立體角數(shù)設(shè)定為800，煙灰產(chǎn)率0.1。油盤擺放、熱電偶和熱輻射計(jì)布置與實(shí)驗(yàn)情況保持一致。

在FDS中，定義燃燒反應(yīng)的方式有兩種：一種是指定單位面積熱釋放速率，另一種是指定燃燒熱與其他熱參數(shù)。這兩種情況下，燃燒模型都可被使用，本文選擇前者，將實(shí)驗(yàn)得到的燃燒速率曲線自定義輸入到FDS中，選用LES湍流模型進(jìn)行數(shù)值模擬。FDS技術(shù)參考指南[11]中詳細(xì)地介紹了FDS的求解原理。其中，能量守恒方程如下：

2 結(jié)果分析與討論

2.1 燃燒速率

燃燒速率作為油池火的基本參數(shù)，是研究燃燒特性的基礎(chǔ)。低硫原油、高硫原油和重石腦油在油盤面積分別為1 m2、2 m2、3 m2時(shí)，燃燒速率隨時(shí)間的變化曲線如圖2所示。文獻(xiàn)[12]已經(jīng)簡(jiǎn)要提及本文各實(shí)驗(yàn)燃燒速率隨時(shí)間的變化規(guī)律，本文將對(duì)其變化機(jī)理進(jìn)行詳細(xì)討論。

由圖2(a)(b)可知，低硫原油和高硫原油的燃燒過(guò)程均可以劃分為三個(gè)燃燒階段，即初期發(fā)展階段、準(zhǔn)穩(wěn)態(tài)燃燒階段和衰減熄滅階段。在燃燒的初期發(fā)展階段，燃燒速率近似線性趨勢(shì)增大，迅速達(dá)到峰值，且油盤面積越大，這種趨勢(shì)越明顯，燃燒速率峰值越早到達(dá)。由此可知，在防治原油池火災(zāi)時(shí)，應(yīng)著重控制燃燒面積，避免形成過(guò)大的原油燃燒面積，以最大程度減小池火災(zāi)帶來(lái)的危害。燃燒速率達(dá)到峰值之后，在一定范圍內(nèi)保持較小波動(dòng)趨勢(shì)，稱之為準(zhǔn)穩(wěn)態(tài)燃燒階段。由低硫原油和高硫原油燃燒速率變化曲線可以看出其穩(wěn)定燃

圖2 燃燒速率隨時(shí)間的變化Fig.2 Change of combustion rate with time

燒階段持續(xù)時(shí)間很短，這主要是因?yàn)樵褪嵌嘟M分混合物，沒(méi)有固定的沸點(diǎn)，燃料不能以恒定的速率吸熱蒸發(fā)。在準(zhǔn)穩(wěn)態(tài)燃燒階段之后燃燒速率呈快速下降趨勢(shì)，直至燃燒結(jié)束，此過(guò)程為衰減熄滅階段。在整個(gè)燃燒過(guò)程中并未出現(xiàn)明顯的多個(gè)穩(wěn)定燃燒階段，此現(xiàn)象與隧道內(nèi)小尺寸原油池火的燃燒速率出現(xiàn)多個(gè)穩(wěn)定燃燒階段顯著不同[12]。一方面是因?yàn)楸菊撐乃M(jìn)行的實(shí)驗(yàn)油盤面積較大，屬于大尺寸油池火燃燒，原油升溫迅速，多組分的沸騰蒸發(fā)過(guò)程發(fā)生在短時(shí)間內(nèi)，彼此分界不明顯，基本呈現(xiàn)出一個(gè)準(zhǔn)穩(wěn)態(tài)燃燒階段，并且隨著油盤面積的增大這種現(xiàn)象更為明顯；另一方面是由于燃燒環(huán)境的影響，隧道屬于狹長(zhǎng)的受限空間，燃燒過(guò)程中會(huì)受到氧氣濃度和空氣流動(dòng)等多重因素的影響。

由圖2(c)可以看出，在重石腦油燃燒的準(zhǔn)穩(wěn)態(tài)階段燃燒速率迅速增至最大，考慮本文所進(jìn)行的實(shí)驗(yàn)油料厚度均為0.01 m，綜合實(shí)驗(yàn)現(xiàn)場(chǎng)觀察到油池壁出現(xiàn)的氣泡和火焰高度的變化特征，可以斷定此階段為沸騰燃燒階段。在沸騰燃燒階段油池壁產(chǎn)生了沸騰現(xiàn)象，液面和池壁的燃料大量氣化，致使燃料蒸氣在短時(shí)間內(nèi)迅速增加，燃燒速率快速增大。該結(jié)果與康泉?jiǎng)俚萚13]對(duì)油池火沸騰燃燒階段研究時(shí)出現(xiàn)的現(xiàn)象類似，這說(shuō)明油池火的沸騰燃燒與油品性質(zhì)和厚度關(guān)系緊密。

在整個(gè)燃燒過(guò)程中每個(gè)實(shí)驗(yàn)燃燒速率均出現(xiàn)小幅度的波動(dòng)。這一方面是由于環(huán)境對(duì)稱重系統(tǒng)的影響；另一方面是因?yàn)樵蛯儆诙嘟M分油品，燃燒時(shí)產(chǎn)生的燃料蒸氣不是單一組分，不利于燃燒速率的穩(wěn)定。在整個(gè)燃燒過(guò)程中出現(xiàn)的燃燒速率波動(dòng)，一定程度上影響著火焰脈動(dòng)和跳躍現(xiàn)象的產(chǎn)生，這種影響在火羽流溫度和測(cè)點(diǎn)接收熱輻射通量隨時(shí)間變化的趨勢(shì)上都有體現(xiàn)。

比較圖2可知，低硫原油和高硫原油燃燒速率的變化趨勢(shì)相似，相同油盤面積下燃燒速率相近，這說(shuō)明硫元素的存在對(duì)燃燒速率幾乎沒(méi)有影響。而同等條件下重石腦油燃燒速率遠(yuǎn)遠(yuǎn)大于低硫原油和高硫原油，其燃燒更為劇烈。

2.2 火羽流溫度

火羽流是火災(zāi)的初始基本形態(tài)，可分為連續(xù)火焰區(qū)、間歇火焰區(qū)和浮力羽流區(qū)，溫度是其重要參數(shù)[14]。本文對(duì)各測(cè)點(diǎn)穩(wěn)定燃燒階段的溫度數(shù)據(jù)求平均，得到火羽流溫度垂直分布，如圖3所示。

圖3 火羽流溫度隨高度的變化Fig.3 Variation of plume temperature with height

由圖3(a)可知，油盤中心溫度隨高度的增加呈先增大后減小的變化趨勢(shì)。這主要是在靠近液面區(qū)域燃料蒸氣濃度比較高，而氧氣濃度不足以使可燃蒸氣完全燃燒，溫度相對(duì)較低。隨著高度的增加，由于體積膨脹和浮升力的作用燃料蒸氣向上流動(dòng)，與此同時(shí)羽流不斷地卷吸新鮮空氣，在某一高度燃料蒸氣與氧氣的比值幾乎達(dá)到當(dāng)量比，在此區(qū)域燃燒最為劇烈，溫度達(dá)到最高；當(dāng)高度進(jìn)一步增加時(shí)，燃料蒸氣濃度相對(duì)較小，不能進(jìn)行劇烈燃燒，溫度趨于變低。對(duì)于同一種燃油，油盤面積越大，火羽流溫度越高。對(duì)比三種實(shí)驗(yàn)油料可以看出，重石腦油火羽流溫度高于低硫原油和高硫原油，這主要是因?yàn)橹厥X油屬于原油蒸餾或石油二次加工切取相應(yīng)餾分而得，以生產(chǎn)高辛烷值汽油為目的，熱值更高，燃燒更為劇烈，一旦發(fā)生重石腦油池火災(zāi)，會(huì)造成更大的危害。

由圖3(b)可以看出，火羽流外部溫度隨高度的增加而降低，且這種下降呈現(xiàn)出先急后緩的趨勢(shì)。這主要是在靠近液面區(qū)域燃料蒸氣充足，并且火羽流外側(cè)有充足的氧氣供給，為充分燃燒提供了條件，這與火羽流中心溫度的變化機(jī)理相印證。通過(guò)對(duì)比低硫原油和高硫原油火羽流外部溫度變化曲線，也能說(shuō)明兩種原油中硫元素的差異對(duì)溫度影響較小。

由以上分析可知，低硫原油、高硫原油和重石腦油的火羽流在油盤中心垂直方向溫度的變化趨勢(shì)相似，而對(duì)于相同面積的池火，重石腦油燃燒時(shí)溫度更高，產(chǎn)生的威脅更大。以重石腦油油盤面積1 m2(Test3)為例，通過(guò)FDS對(duì)其模擬得到不同燃燒時(shí)間段溫度切片，如圖4(a)所示。對(duì)于低硫原油和高硫原油的大尺寸油池火，其火羽流外側(cè)垂直方向溫度的變化情況相似，與含硫量相關(guān)性不大。本節(jié)以高硫原油油盤面積3 m2(Test7)為例，通過(guò)FDS對(duì)其模擬得到不行燃燒時(shí)間段溫度切片，如圖4(b)所示。

由圖4可以直觀地看出，火羽流溫度在垂直方向上呈不連續(xù)分布，這說(shuō)明其不是連續(xù)的火焰柱。從溫度的空間分布特征來(lái)看，油盤中心線區(qū)域的溫度是整個(gè)火羽流溫度最高的部分，基本呈現(xiàn)出自下而上先增加后減小的變化，而火羽流外部垂直方向的溫度自下而上逐漸減小。從溫度的大小來(lái)看，重石腦油油盤面積為1 m2時(shí)，油盤的中心線溫度達(dá)400 ℃～900 ℃，而單個(gè)油盤中心線溫度也達(dá) 300 ℃～600 ℃；高硫原油油盤面積為3 m2時(shí)，火羽流外側(cè)溫度距地面0.5 m處大約750 ℃，隨高度的上升逐漸下降到200 ℃左右。從火羽流區(qū)域來(lái)看，在連續(xù)火焰區(qū)高溫趨于集中，在間歇火焰區(qū)和浮力羽流區(qū)徑向溫度的變化較大。由此可知，溫度的模擬結(jié)果基本與實(shí)驗(yàn)測(cè)量結(jié)果一致。

由圖4(a)可知，重石腦油燃燒在沸騰燃燒階段火羽流的高度和寬度達(dá)到最大，這說(shuō)明在此階段燃燒最為劇烈，產(chǎn)生的危害最大，應(yīng)采取措施避免形成薄層重石腦油池火。由圖4(b)可知，對(duì)于大尺寸油池火，在其燃燒的三個(gè)階段中，火羽流外側(cè)溫度始終很高，影響范圍很大，有直接灼傷附近人員、引燃周圍物體的危險(xiǎn)，應(yīng)該采取措施避免形成大面積油池火。

圖4 不同燃燒階段溫度切片F(xiàn)ig.4 Temperature slices at different combustion stages

2.3 測(cè)點(diǎn)接收熱輻射通量

對(duì)于大尺寸油池火，無(wú)風(fēng)或者較小風(fēng)速時(shí)只有少量熱氣往下風(fēng)向傳遞，即熱對(duì)流相對(duì)較弱，熱輻射是主要的傳熱方式。本節(jié)根據(jù)大型開(kāi)放空間碳?xì)浠衔锘馂?zāi)危險(xiǎn)性計(jì)算公式[15]得到輻射模型計(jì)算值，并與FDS模擬值和實(shí)驗(yàn)測(cè)量值相比較，如圖5所示。

圖5 測(cè)點(diǎn)接收熱輻射通量比較Fig.5 Comparison of heat radiation flux received at measuring points

由圖5可知，在靠近火源位置時(shí)點(diǎn)源模型和圓柱模型的計(jì)算值偏大，模擬值與測(cè)量值幾乎相等；隨著離火源距離的增大，點(diǎn)源模型和圓柱模型計(jì)算值與測(cè)量值逐漸接近，在距火源較遠(yuǎn)的位置計(jì)算值、模擬值和測(cè)量值基本保持一致。這主要是因?yàn)辄c(diǎn)源模型將火焰看作一個(gè)對(duì)外輻射為各向同性的點(diǎn)，距離火源越近，這種假設(shè)的局限性越大；圓柱模型認(rèn)為火焰表面均勻向外產(chǎn)生熱輻射，雖然可以將火焰輻射通量分為羽流底部發(fā)光區(qū)域和上方煙羽流兩部分進(jìn)行計(jì)算，但距離火源很近時(shí)，計(jì)算結(jié)果仍較測(cè)量值偏高。

通過(guò)對(duì)圖5橫向和縱向比較可知，油盤等效直徑越大，測(cè)點(diǎn)接受的熱輻射通量越大；燃燒面積相同時(shí)，低硫原油和高硫原油相同位置接收的熱輻射通量相差不大，而重石腦油遠(yuǎn)遠(yuǎn)高于兩者，這種變化趨勢(shì)與前文燃燒速率和火羽流溫度的變化特征在本質(zhì)上相一致。從火源功率來(lái)看，當(dāng)火源功率在0.747 0～1.910 8 MW時(shí)，模擬值與實(shí)驗(yàn)值幾乎相等，而點(diǎn)源模型和圓柱模型遠(yuǎn)遠(yuǎn)高于實(shí)驗(yàn)值；火源功率在2.270 4～6.918 8 MW時(shí)，F(xiàn)DS的模擬結(jié)果相對(duì)于測(cè)量值偏小，點(diǎn)源模型計(jì)算值高出測(cè)量值很多，而圓柱模型計(jì)算結(jié)果與實(shí)驗(yàn)值的逐漸靠近，這主要是因?yàn)镕DS模擬時(shí)，室外燃燒環(huán)境存在瞬時(shí)變化，而邊界條件設(shè)置不能與其完全保持一致。

總體來(lái)看，通過(guò)將實(shí)驗(yàn)燃燒速率曲線輸入到FDS中，對(duì)測(cè)點(diǎn)接收熱輻射通量的模擬是比較準(zhǔn)確的，這為不同場(chǎng)景大尺寸油池火災(zāi)的研究提供了可靠的選擇。同時(shí)，也驗(yàn)證了圓柱模型對(duì)于大功率油池火災(zāi)不同位置接收熱輻射通量計(jì)算的科學(xué)性。本文在此研究基礎(chǔ)上，以Test8為例，基于熱通量-時(shí)間準(zhǔn)則、聯(lián)立圓柱模型方程和比率方程[16-17]，確定其“安全距離-時(shí)間”臨界曲線，如圖6所示。

通過(guò)“安全距離-時(shí)間”臨界曲線可以直觀地看出在距火源中心不同位置處人員可以暴露的時(shí)間，這對(duì)于科學(xué)地預(yù)防次生災(zāi)害發(fā)生和消防救援意義重大。

3 結(jié) 論

本文通過(guò)對(duì)室外開(kāi)放空間油池火的實(shí)驗(yàn)研究與數(shù)值模擬，得到以下結(jié)論：

(1) 低硫原油和高硫原油的燃燒過(guò)程均出現(xiàn)三個(gè)燃燒階段。在燃燒初期發(fā)展階段，燃燒速率近似線性趨勢(shì)增大，穩(wěn)定燃燒階段不明顯，且未出現(xiàn)明顯的多個(gè)穩(wěn)定燃燒階段，重石腦油出現(xiàn)沸騰燃燒現(xiàn)象。

(2) 油盤中心溫度隨高度的增加先增大后減小，火羽流外部溫度隨高度的增加而減小，且下降趨勢(shì)先急后緩。在連續(xù)火焰區(qū)高溫集中，在間歇火焰區(qū)和浮力羽流區(qū)溫度變化較大。

(3) 將實(shí)驗(yàn)得到的燃燒速率曲線自定義到FDS中，可以實(shí)現(xiàn)對(duì)測(cè)點(diǎn)接收熱輻射通量的準(zhǔn)確模擬。火源功率越大，圓柱模型的計(jì)算值與實(shí)驗(yàn)值越接近。

(4) 同等條件下，重石腦油產(chǎn)生的危害大于低硫原油和高硫原油，大尺寸油池火產(chǎn)生的危險(xiǎn)遠(yuǎn)遠(yuǎn)大于小尺寸油池火。

(5) 對(duì)于大功率原油池火災(zāi)，可以確定“安全距離-時(shí)間”臨界曲線，以直觀地掌握距火源中心不同位置處人員可以暴露的時(shí)間。