段慶昊， 余建星， 梁 峰， 趙羿羽， 樊志遠(yuǎn)

(1. 天津大學(xué) 建筑工程學(xué)院 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072；2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240；3.海洋石油工程股份有限公司設(shè)計(jì)公司，天津 300451)

?

海上溢油池火災(zāi)災(zāi)害的計(jì)算方法研究

段慶昊1,2， 余建星1,2， 梁 峰3， 趙羿羽1,2， 樊志遠(yuǎn)1,2

(1. 天津大學(xué) 建筑工程學(xué)院 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072；2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240；3.海洋石油工程股份有限公司設(shè)計(jì)公司，天津 300451)

基于國(guó)內(nèi)外池火災(zāi)研究的理論,建立海面上有風(fēng)情況下溢油火災(zāi)模型,并提出對(duì)海上作業(yè)溢油事故引發(fā)的池火災(zāi)新計(jì)算方法。其驗(yàn)證表明,風(fēng)對(duì)池火災(zāi)的影響主要體現(xiàn)在火焰傾角上,風(fēng)速越大,火焰傾角也越大,隨著池火災(zāi)直徑的增加,風(fēng)的影響作用逐漸減?。伙L(fēng)的影響還體現(xiàn)在火災(zāi)發(fā)生區(qū)域上風(fēng)向及下風(fēng)向的熱輻射分布上,導(dǎo)致同一熱輻射強(qiáng)度下下風(fēng)向的距離要大于上風(fēng)向的距離。研究了給定風(fēng)速下不同直徑池火災(zāi)熱輻射分布,確定不同傷害準(zhǔn)則下的上風(fēng)向和下風(fēng)向的安全距離公式,可以為溢油火災(zāi)現(xiàn)場(chǎng)人員疏散和安全技術(shù)管理提供有效依據(jù)。

海上溢油；池火災(zāi)；安全距離

0 引言

在海洋油氣開發(fā)、生產(chǎn)以及運(yùn)輸作業(yè)中，易燃、易爆、有毒的危險(xiǎn)化學(xué)品應(yīng)用十分廣泛，使得火災(zāi)、爆炸等惡性事件屢屢發(fā)生，這類事件的破壞作用一般并不僅僅作用于事故發(fā)生區(qū)域的人員和設(shè)施，其周圍一定范圍內(nèi)的人員和設(shè)施都會(huì)受到不同程度的傷害：一方面會(huì)導(dǎo)致更大范圍內(nèi)人員死亡、受傷概率增加；另一方面，周圍的一些特殊設(shè)施在熱輻射、爆炸波等作用下可能會(huì)引發(fā)更為嚴(yán)重的事故，其造成的損失是難以估計(jì)的。

國(guó)外對(duì)石油化工行業(yè)火災(zāi)事故的研究起步較早，主要有實(shí)驗(yàn)研究和數(shù)值模擬研究?jī)蓚€(gè)方面。在實(shí)驗(yàn)研究方面，Thomas[1]、Heskestad[2]等人根據(jù)實(shí)驗(yàn)數(shù)據(jù)，提出了池火災(zāi)相關(guān)的數(shù)學(xué)模型。數(shù)值模擬方面，Sinai從大量的數(shù)值模擬結(jié)果同實(shí)驗(yàn)結(jié)果相吻合得出結(jié)論，利用CFD 技術(shù)對(duì)火災(zāi)進(jìn)行數(shù)值模擬是可行的[3]。國(guó)內(nèi)主要有大連海事大學(xué)、上海海事大學(xué)、天津大學(xué)等幾所院校[4]和科研機(jī)構(gòu)[5]在積極開展海上安全風(fēng)險(xiǎn)評(píng)估技術(shù)的研究，目前我國(guó)在溢油事故引起的危害及后果方面還缺乏系統(tǒng)全面的研究工作。現(xiàn)階段海洋油氣開發(fā)事業(yè)高速發(fā)展，對(duì)海洋油氣泄漏事故提出有效的預(yù)防和救助措施顯的尤其重要。

該文在國(guó)內(nèi)外學(xué)者研究的基礎(chǔ)上，利用Fluent軟件進(jìn)行海上溢油池火災(zāi)事故模擬，確定不同直徑池火災(zāi)對(duì)應(yīng)的安全距離，為海上溢油火災(zāi)事故的處理工作提供安全保障。

1 池火災(zāi)理論基礎(chǔ)

1.1 池火災(zāi)概念

池火災(zāi)[6]一般是指可燃液體泄漏后遇到點(diǎn)火源成為固定形狀或者不定形狀的液池火災(zāi)，其對(duì)周圍環(huán)境的主要傷害來(lái)源為熱輻射。

1.2 熱通量準(zhǔn)則

熱通量準(zhǔn)則[7]主要依據(jù)目標(biāo)接受的熱通量來(lái)衡量目標(biāo)是否被傷害破壞，該準(zhǔn)則主要適用于熱通量的作用時(shí)間比目標(biāo)達(dá)到熱平衡所需要的時(shí)間長(zhǎng)的情況。其中的熱通量臨界值需要通過對(duì)不同物質(zhì)的大量實(shí)驗(yàn)來(lái)確定。如果目標(biāo)接受的熱通量處于被傷害破壞的臨界值之上，就會(huì)受到傷害破壞作用，否則，目標(biāo)就不會(huì)受到傷害破壞。在較長(zhǎng)時(shí)間內(nèi)能夠比較穩(wěn)定地燃燒的穩(wěn)態(tài)火災(zāi)熱輻射作用下，人員傷害和設(shè)備破壞的臨界值見表1。

表1 熱輻射傷害準(zhǔn)則

1.3 物理幾何模型

海上溢油引發(fā)的池火災(zāi)，屬于非預(yù)混燃燒，燃燒產(chǎn)物中主要的輻射物質(zhì)為CO2、H2O 等氣體，以及炭黑微粒、焦炭微粒和灰粒?；鹧嬷袣怏w成分的輻射能力比固體微粒的輻射能力低，所以，可以將氣體成分的輻射疊加到固體微粒的輻射上。

(1) 物理模型[8]

為了方便計(jì)算，對(duì)池火災(zāi)燃燒模型做出如下簡(jiǎn)化：

(a) 池火災(zāi)火焰簡(jiǎn)化為圓柱形，有風(fēng)情況下為傾斜圓柱。

(b) 將火焰對(duì)周圍的熱輻射歸結(jié)為火焰外表面對(duì)周圍的輻射，并假定為灰體輻射，黑度取0.96炭黑的黑度。

(c) 假定火災(zāi)為穩(wěn)態(tài)，研究區(qū)域外圍為黑體。

(2) 幾何模型

研究區(qū)域?yàn)槔锩鎴A柱到外面圓柱的中間部分，里面圓柱的尺寸為假定火焰的尺寸，外面圓柱的尺寸根據(jù)經(jīng)驗(yàn)取得，幾何模型示意圖如圖1所示。

圖1 幾何模型示意圖

1.4 火焰模型的尺寸計(jì)算

1.4.1 燃燒速率v

液池中液體沸點(diǎn)高于環(huán)境溫度時(shí)其燃燒速率為：

(1)

式中：v為單位表面積燃燒速度；m為常，取0.001；Hc為液體燃燒熱；Cp為液體的定壓比熱；Tb為液體的沸點(diǎn)；Ta為環(huán)境溫度；H為液體的氣化熱。

液池中液體沸點(diǎn)低于環(huán)境溫度時(shí)，其燃燒速率為：

(2)

1.4.2 有風(fēng)情況下火焰高度h

池火災(zāi)火焰高度的關(guān)系模型是基于實(shí)驗(yàn)數(shù)據(jù)及量綱分析等所得到的經(jīng)驗(yàn)公式，應(yīng)用較多的是Thomas關(guān)系模型。

(3)

式中：h為火焰高度；D為液池直徑；g為重力加速度；ρ0為空氣密度；uc為特征風(fēng)速，由下式計(jì)算：

(4)

1.4.3 火焰傾角θ

火焰受到風(fēng)的作用就會(huì)產(chǎn)生一定的傾角，該角度大小與風(fēng)速有直接關(guān)系，通常定義火焰傾角為火焰流在垂直方向上的傾斜角度。

關(guān)于火焰傾角θ的預(yù)測(cè)模型，一般采用SFPE(2002)的規(guī)定，即：

當(dāng)計(jì)算得uc≥1時(shí)，

(5)

當(dāng)計(jì)算得uc≤1時(shí)，

(6)

總結(jié)大量實(shí)驗(yàn)，得出不同的預(yù)測(cè)模型下B、C的取值見表2。

表2 火焰傾角預(yù)測(cè)模型B與C取值

1.4.4 研究區(qū)域的確定

模型研究區(qū)域是火焰與外圓柱中間部分，經(jīng)過多次試算，研究區(qū)域高度取2.5H，半徑取液池半徑的20倍，既保證了研究區(qū)域足夠大，也考慮到了模型的計(jì)算速度。

2 新計(jì)算方法

為了達(dá)到研究的目的，選取丙烷(C3H8)作為研究對(duì)象，已知燃燒熱為49 036 kJ/kg，氣化熱為410.5 kJ/kg、沸點(diǎn)為-42.1℃。初步擬定取10組不同的池火災(zāi)直徑進(jìn)行計(jì)算分析，分別為D=5 m，10 m，15 m，20 m，25 m，30 m，35 m，40 m，45 m，50 m。

以D=30 m直徑池火災(zāi)為例介紹計(jì)算過程，由式(1)～式(6)計(jì)算得出h=45 m，θ=40.8℃。

2.1 網(wǎng)格劃分

網(wǎng)格劃分選用ICEM軟件，最低網(wǎng)格質(zhì)量0.45，網(wǎng)格數(shù)量146 851。模型底面和火焰表面網(wǎng)格劃分情況如圖2所示。

圖2 模型底面和火焰網(wǎng)格

2.2 邊界條件

四周以及頂部壁面黑度為1，溫度為0。地面黑度為1，溫度為0。輻射熱源黑度取0.96，溫度取1 300 K，大氣介質(zhì)吸收系數(shù)為0.261/m。

2.3 結(jié)果分析

將模型運(yùn)用Fluent求解，熱輻射強(qiáng)度在垂直面的分布如圖3所示?？紤]到安全的因素，取海面上5 m處平面為參考平面，熱輻射強(qiáng)度δ的分布如圖4所示。

圖3 垂直面熱輻射強(qiáng)度分布

圖4 h=5 m水平面熱輻射強(qiáng)度分布

分別計(jì)算D=5 m，10 m，15 m，20 m，25 m，30 m，35 m，40 m，45 m，50 m的池火災(zāi)模型，將每個(gè)模型上風(fēng)向和下風(fēng)向?qū)?yīng)各個(gè)傷害準(zhǔn)則的安全距離整理見表3。

表3 結(jié)果統(tǒng)計(jì)各傷害準(zhǔn)則對(duì)應(yīng)的安全距離

將計(jì)算結(jié)果做統(tǒng)計(jì)學(xué)分析，可以看出，池直徑(m)與安全距離(m)的線性關(guān)系良好，以池直徑為自變量，安全距離為因變量做線性擬合。

圖5 q=4 000計(jì)算結(jié)果線性擬合

圖6 q=12 500計(jì)算結(jié)果線性擬合

上風(fēng)向： y=12.27+3.9·x

下風(fēng)向： y=23.47+4.0·x

上風(fēng)向： y=5.8+2.6·x(9)

下風(fēng)向： y=18.07+2.6·x

圖7 q=25 000計(jì)算結(jié)果線性擬合

圖8 q=375 000計(jì)算結(jié)果線性擬合

上風(fēng)向： y=3.8+1.94·x

下風(fēng)向： y=14.53+1.93·x

上風(fēng)向： y=-0.73+1.75·x

下風(fēng)向： y=11.27+1.73·x

3 結(jié)論

該文研究了有風(fēng)情況下池火災(zāi)模型，模型相關(guān)參數(shù)選取主要依據(jù)國(guó)內(nèi)外普遍使用的經(jīng)驗(yàn)公式，通過實(shí)例驗(yàn)證，提出了上風(fēng)向以及下風(fēng)向安全距離的數(shù)學(xué)經(jīng)驗(yàn)公式，具有可靠性和適用性。

由計(jì)算分析，可以得出下述規(guī)律：

(1) 相同傷害準(zhǔn)則下，隨著池直徑的增加，對(duì)應(yīng)的安全距離增大。

(2) 上風(fēng)向安全距離直線與下風(fēng)向安全距離擬合直線近似平行。

(3) 不同的傷害準(zhǔn)則下，擬合直線的斜率不同，隨著熱輻射強(qiáng)度臨界值的增加，直線的斜率減小。

[1] Thomas P H.The Size of Flames from Natural Fires[J]. Symposium(International) on Combustion,1963,9(1):844-859.

[2] Heskestad G. Engineering Relations for Fire Plumes[J].Fire Safety,1987(7):25-32.

[3] Sinai Y L, Owens M P. Validation of CFD Modeling of Unconfined Pool Fires with Cross-Wind:Flame Geometry[J].Fire Safety,1995(24):1-34.

[4] Fan P, Liu Q Q, Li J C, et al.Numerical simulations of the infiltration of rain on the crack slope[J]. Sciences in China, 2005：107-120.

[5] Li X S, Yang B, Li G, et al. Numerical Simulation of Gas Production from Natural Gas Hydrate Using a Single Horizontal Well by Depressurization in Qilian Mountain Permafrost[J]. Ind Eng Chem Res， 2012，51(11):24-32.

[6] 王曉寧，李麗霞，王明賢.有風(fēng)情況下池火災(zāi)的數(shù)值模擬[J].江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,27(4):328-331.

[7] 宇德明,馮長(zhǎng)根,曾慶軒,等.熱輻射的破壞準(zhǔn)則和池火災(zāi)的破壞半徑[J].中國(guó)安全科學(xué)報(bào),1996,6(2):5-10.

[8] 李慧,蔣軍成,王若菌.池火災(zāi)熱輻射的數(shù)值研究[J].中國(guó)安全科學(xué)報(bào),2005,15(10):7-10.

The New Calculation Method on the Fire Hazards by Oil Spill Pool Fire at Sea

DUAN Qing-hao1,2, YU Jian-xing1,2, LIANG Feng3，ZHAO Yi-yu1,2, Fan Zhi-yuan1,2

(1.State Key Laboratory of Hydraulic Engineering Simulation and Safety ,Tianjin University, Tianjin 300072,China; 2. Collaborative Innovation Center for Advanced Ship and Deep-sea Exploration, Shanghai 200240,China; 3. Offshore Oil Engineering Co., Ltd, Designing Company, Tianjin 300451,China)

Based on the theory of pool fire research at home and abroad, the model of oil spill was set up, and the method of CFD simulation was adopted to research pool fire caused by offshore oil spill accidents.The results show that the effect of wind is mainly reflected in the flame angle. The higher the wind speeds, the greater the flame angle was. With the increasing of pool fire diameter, the influence of wind decreased;also the influence of wind is the distribution of thermal radiation in the fire area, at the same thermal radiation intensity ,the distance of down wind direction is greater than the up winddirection.Studying the different diameter of pool fires for a given wind speed, and determining the safety distance formula, can provide effective basis for the fire evacuation and the technology management of oil spill.

oil spill at sea; pool fire; safety distance

2015-07-23

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃資助項(xiàng)目(2014CB046803)，國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào)：51239008)資助項(xiàng)目，國(guó)家自然科學(xué)基金(批準(zhǔn)號(hào)：51379145)資助項(xiàng)目，“海洋工程作業(yè)安全模擬系統(tǒng)及工程應(yīng)用研究(首期)”項(xiàng)目資助。

段慶昊(1989-)，男，碩士研究生。

1001-4500(2016)05-0034-06

TE56

A