楊兆晶 侯 磊 朱 淼

中國(guó)石油大學(xué)(北京)機(jī)械與儲(chǔ)運(yùn)工程學(xué)院, 北京 102249

0 前言

隨著國(guó)家能源結(jié)構(gòu)的調(diào)整及LNG接收站的逐步建成投產(chǎn),LNG進(jìn)口量迅速攀升,目前我國(guó)已成為世界第二大LNG進(jìn)口國(guó)[1]。大型LNG儲(chǔ)罐作為接收站的核心設(shè)備,其安全運(yùn)行是整個(gè)站場(chǎng)安全生產(chǎn)的重中之重。一旦LNG發(fā)生泄漏,低溫氣云在近地面聚積,可能對(duì)周圍環(huán)境、設(shè)備和人員造成嚴(yán)重危害[2-4]。國(guó)外從20世紀(jì)60年代末開始對(duì)LNG泄漏擴(kuò)散過程進(jìn)行分析研究,主要研究方法包括大型現(xiàn)場(chǎng)試驗(yàn)[5-6]、風(fēng)洞實(shí)驗(yàn)研究[7-8]和數(shù)值模擬,目前應(yīng)用最為廣泛的是數(shù)值模擬。常用數(shù)學(xué)模型主要包括高斯修正模型、箱及相似模型、經(jīng)驗(yàn)唯象模型、淺層模型和計(jì)算流體力學(xué)(CFD)模型,其中,CFD模型最為真實(shí)可靠[9]。Gavelli F等[10-13]將現(xiàn)場(chǎng)試驗(yàn)結(jié)果與模擬結(jié)果進(jìn)行對(duì)比,驗(yàn)證了CFD軟件模擬LNG泄漏擴(kuò)散的有效性。莊學(xué)強(qiáng)、李清等[14-15]選用CFD軟件中的離散相模型對(duì)LNG泄漏擴(kuò)散進(jìn)行模擬,秦雅琦等[16]將泄漏源從罐壁小孔轉(zhuǎn)化為液池再進(jìn)行氣體泄漏擴(kuò)散模擬,均未對(duì)LNG連續(xù)泄漏后的相變過程進(jìn)行充分研究。針對(duì)此問題,本文在利用CFD軟件模擬LNG泄漏擴(kuò)散時(shí)未對(duì)泄漏源進(jìn)行轉(zhuǎn)化,而是加入相變模型的編寫設(shè)置,并在此基礎(chǔ)上分析環(huán)境條件對(duì)大型LNG儲(chǔ)罐泄漏擴(kuò)散的影響規(guī)律,對(duì)比分析單罐及罐區(qū)泄漏的氣云分布,為L(zhǎng)NG儲(chǔ)罐布局及泄漏防控提供參考。

1 研究場(chǎng)景和理論模型

1.1 場(chǎng)景建立

以我國(guó)南方某LNG接收站內(nèi)16×104m3大型儲(chǔ)罐為研究對(duì)象,假設(shè)地面平坦且周圍無(wú)其他建筑,利用ICEM軟件建立簡(jiǎn)化后的單罐及罐區(qū)模型,在外流場(chǎng)區(qū)域進(jìn)行網(wǎng)格劃分。單個(gè)LNG儲(chǔ)罐模型見圖1,罐體高38.45 m,直徑82 m,穹頂高11 m,LNG罐區(qū)模型見圖2,由4座容量相同的儲(chǔ)罐組成。

圖1 單罐模型圖Fig.1 Single tank model

圖2 罐區(qū)模型圖Fig.2 Tank farm model

LNG儲(chǔ)罐泄漏源分為氣相和液相兩種,氣相泄漏主要發(fā)生在儲(chǔ)罐頂部,離地面較遠(yuǎn)且泄漏后擴(kuò)散較快,與液相泄漏相比發(fā)生危險(xiǎn)的可能性較低。因此,主要研究LNG儲(chǔ)罐液相連續(xù)泄漏的擴(kuò)散過程以及環(huán)境因素可能對(duì)其造成的影響。根據(jù)《化工裝備事故分析與預(yù)防》[17]中1949～1988年全國(guó)化工行業(yè)事故發(fā)生情況相關(guān)資料,持續(xù)泄漏小孔徑為100 mm,大孔徑為200 mm,為了模擬LNG泄漏的不利情況,將泄漏定為大孔徑泄漏,為顯著觀察LNG重氣云擴(kuò)散情況,將泄漏位置定在儲(chǔ)罐背風(fēng)中心高處。研究風(fēng)速影響時(shí),將環(huán)境風(fēng)速設(shè)為輕風(fēng)3 m/s、微風(fēng)5 m/s、和風(fēng)7 m/s、勁風(fēng)10 m/s四種情況。

1.2 理論模型

1.2.1 泄漏源強(qiáng)計(jì)算模型

儲(chǔ)罐泄漏口面積遠(yuǎn)小于儲(chǔ)罐內(nèi)LNG液面積,泄漏不會(huì)引起液位明顯下降,可將泄漏源看成連續(xù)泄漏源。根據(jù)伯努利方程計(jì)算液相泄漏速率,其表達(dá)式為:

(1)

式中:QL為液相泄漏速率,kg/s;Chole為流量系數(shù),無(wú)量綱常數(shù);Ahole為泄漏孔的面積,m2;ρL為儲(chǔ)罐中液體的密度,kg/m3;pT為儲(chǔ)罐內(nèi)液體上方壓力,Pa;patm為大氣壓力,Pa;g為重力加速度,9.8 m2/s;H為泄漏源與儲(chǔ)罐內(nèi)液面的高度差,m。

模擬計(jì)算采用的LNG密度為430 kg/m3,比熱為2 055 J/(kg·K),導(dǎo)熱率為0.21 W/(m·K)，黏度為0.000 118 3 kg/(m·s),相對(duì)分子質(zhì)量為16.41 kg/kmol,蒸發(fā)潛熱為509 332 J/kg,沸騰溫度為111.66 K。

1.2.2 流動(dòng)模型

LNG儲(chǔ)罐液相泄漏擴(kuò)散過程涉及相變、質(zhì)量傳遞、熱量傳遞、湍流流動(dòng)等復(fù)雜物理問題,需要根據(jù)連續(xù)性方程、動(dòng)量守恒方程、能量守恒方程對(duì)質(zhì)量傳遞、動(dòng)量傳遞和能量傳遞等過程進(jìn)行細(xì)致描述。LNG儲(chǔ)罐液相泄漏擴(kuò)散過程涉及LNG、天然氣和空氣三相,可采用多相流模型來描述流動(dòng)過程,其中,混合模型適用于流動(dòng)中存在相混合或分離的情況[18]。

目前標(biāo)準(zhǔn)k-ε模型已經(jīng)得到廣泛應(yīng)用,但是標(biāo)準(zhǔn)k-ε模型假設(shè)湍動(dòng)黏度系數(shù)各向同性,對(duì)于儲(chǔ)罐這種彎曲壁面流動(dòng)的模擬誤差較大,而Realizable k-ε模型在湍動(dòng)黏度系數(shù)計(jì)算公式中引入了與旋轉(zhuǎn)和曲率有關(guān)的變量,考慮了湍流流動(dòng)各向異性的情況,能夠有效用于彎曲壁面流動(dòng)、彎曲流線流動(dòng)、邊界層流動(dòng)。Tauseef S M等[19-20]利用CFD方法將不同湍流模型模擬結(jié)果與實(shí)驗(yàn)值對(duì)比,結(jié)果表明,Realizable k-ε模型在時(shí)間及空間濃度分布方面與實(shí)際結(jié)果最為吻合,因此,選擇Realizable k-ε模型作為湍流模型。

1.2.3 蒸發(fā)相變模型

泄漏出的LNG為液態(tài),經(jīng)過與空氣的質(zhì)量熱量傳遞氣化為天然氣,利用FLUENT自帶的UDF功能編寫LNG蒸發(fā)相變模型。在已有的相變模型中,Lee W H[21]提出的模型應(yīng)用最為廣泛,兩相間的質(zhì)量傳遞公式為:

(2)

(3)

式中:T為混合區(qū)單元溫度,K;Tsat為飽和溫度,K;ml→v為液相轉(zhuǎn)化為氣相的相變率,kg/(m3·s);mv→l為氣相轉(zhuǎn)化為液相的相變率,kg/(m3·s);r為控制相變強(qiáng)度的因子,根據(jù)Schepper S C K D等[22]對(duì)烴原料流動(dòng)蒸發(fā)過程的模擬將r取為0.1 s-1;αl、αv分別為液相、氣相的表面張力,N/m;ρl、ρv分別為液相、氣相的密度,kg/m3。

2 結(jié)果與分析

2.1 各因素對(duì)LNG儲(chǔ)罐泄漏擴(kuò)散的影響

2.1.1 風(fēng)速

風(fēng)對(duì)LNG泄漏擴(kuò)散會(huì)產(chǎn)生兩方面的影響:一方面,風(fēng)的平流輸送作用會(huì)使LNG氣云隨風(fēng)向遠(yuǎn)處運(yùn)動(dòng);另一方面,風(fēng)會(huì)對(duì)LNG氣云起到混合稀釋作用,有助于高濃度氣云的擴(kuò)散。通常天然氣的爆炸上限(UFL)為15 %,爆炸下限(LFL)為5 %,安全區(qū)域分界線為1/2 LFL。對(duì)比四種風(fēng)速下氣云隨時(shí)間到達(dá)下風(fēng)向的最遠(yuǎn)距離見圖3～4,大風(fēng)速可加快下風(fēng)向氣云的擴(kuò)散速度,使得1/2 LFL氣云到達(dá)更遠(yuǎn)距離,而UFL氣云則會(huì)因稀釋速度的加快呈現(xiàn)與之相反的趨勢(shì)。由此可見,風(fēng)對(duì)高濃度氣云產(chǎn)生的平流輸送作用較弱,對(duì)其產(chǎn)生的混合稀釋作用較強(qiáng)。

圖3 不同風(fēng)速下UFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離曲線圖Fig.3 The farthest downwind distance UFL gas cloud can reach under different wind speeds

圖4 不同風(fēng)速下1/2 LFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離曲線圖Fig.4 The farthest downwind distance 1/2 LFL gas cloud can reach under different wind speeds

在連續(xù)泄漏80 s內(nèi),UFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離總體呈現(xiàn)出隨時(shí)間小波動(dòng)上升的趨勢(shì),但其隨風(fēng)速的變化趨勢(shì)最初卻并不穩(wěn)定。這是由于此時(shí)的氣云仍處于儲(chǔ)罐背風(fēng)接近罐壁處,儲(chǔ)罐的阻擋使得此處風(fēng)速整體較低,差異性不明顯,同時(shí)此處的紊流作用也使得氣云擴(kuò)散變得更加復(fù)雜,擾亂固有的流動(dòng)趨勢(shì)。但隨著連續(xù)泄漏的進(jìn)行,風(fēng)速的影響逐漸凸顯,直到泄漏80 s時(shí),3 m/s風(fēng)速下的UFL最遠(yuǎn)距離達(dá)到82 m,而10 m/s風(fēng)速下則剛至53.5 m。高風(fēng)速使得UFL氣云在未達(dá)到較遠(yuǎn)距離時(shí)即稀釋為低濃度氣云,縮短了高濃度氣云的存在時(shí)間。就1/2 LFL氣云而言,風(fēng)對(duì)氣云的平流輸送作用非常顯著,80 s時(shí)的最遠(yuǎn)距離已經(jīng)有很大差距,3 m/s風(fēng)速下為156 m,10 m/s風(fēng)速下已達(dá)到367 m。因此,在氣云分布未達(dá)到穩(wěn)定狀態(tài)時(shí),高風(fēng)速促進(jìn)了氣云的順風(fēng)擴(kuò)散,同一時(shí)刻的1/2 LFL氣云隨風(fēng)速增大能夠到達(dá)更遠(yuǎn)的距離,變化趨勢(shì)顯著且變化幅度均勻。

LNG的泄漏與擴(kuò)散在達(dá)到平衡時(shí),氣云分布不再出現(xiàn)大幅度變化,基本趨于穩(wěn)定。表1顯示,3 m/s風(fēng)速下的1/2 LFL氣云在泄漏580 s時(shí)能夠到達(dá)下風(fēng)向最遠(yuǎn)距離554 m;風(fēng)速大于3 m/s時(shí),1/2 LFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離均低于554 m,且隨風(fēng)速增大而減小,同時(shí),1/2 LFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離時(shí)間縮短,表明大風(fēng)速不僅能夠促進(jìn)氣云擴(kuò)散,減小其分布范圍,還能夠使氣云更快達(dá)到穩(wěn)定狀態(tài)。

表1 不同風(fēng)速下1/2 LFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離及時(shí)間表

Tab.1 The farthest downwind distance and time 1/2 LFL gas cloud can reach under different wind speeds

風(fēng)速/(m·s-1)1/2LFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離/m1/2 LFL氣云到達(dá)下風(fēng)向最遠(yuǎn)距離時(shí)間/s35545805492500744947010405450

2.1.2 風(fēng)向

上述研究場(chǎng)景均是基于泄漏孔在背風(fēng)側(cè)的假設(shè),為詳細(xì)分析風(fēng)向?qū)Υ笮蚅NG儲(chǔ)罐泄漏擴(kuò)散的影響,選取順風(fēng)、側(cè)風(fēng)和逆風(fēng)三個(gè)典型風(fēng)向進(jìn)行模擬對(duì)比。

在LNG連續(xù)泄漏至80 s時(shí),順風(fēng)泄漏產(chǎn)生的1/2 LFL氣云體積為24 229.9 m3,側(cè)風(fēng)泄漏產(chǎn)生的1/2 LFL氣云體積為22 680.9 m3,逆風(fēng)泄漏產(chǎn)生的1/2 LFL氣云體積為21 905.4 m3。根據(jù)圖5順風(fēng)泄漏到逆風(fēng)泄漏的結(jié)果,同一時(shí)刻的氣云體積逐漸減小,這是由于儲(chǔ)罐本身對(duì)氣云擴(kuò)散的阻礙作用逐漸加強(qiáng),風(fēng)對(duì)氣云的輸送作用逐漸減弱,導(dǎo)致氣云擴(kuò)散速度減慢。對(duì)長(zhǎng)時(shí)間泄漏而言,順風(fēng)泄漏時(shí)泄漏速率與蒸發(fā)速率能最快達(dá)到穩(wěn)定,而逆風(fēng)泄漏時(shí)則最遲達(dá)到穩(wěn)定,同時(shí)緩慢的擴(kuò)散速度會(huì)使得高濃度氣云的存在時(shí)間延長(zhǎng)進(jìn)而導(dǎo)致危險(xiǎn)氣云體積增加,危險(xiǎn)區(qū)域增大。

a)順風(fēng)泄漏a)Downwind leakage

b)側(cè)風(fēng)泄漏b)Crosswind leakage

c)逆風(fēng)泄漏c)Upwind leakage

2.1.3 障礙物

儲(chǔ)罐周圍還有一些如辦公樓、控制室等較大的建筑物,為了說明儲(chǔ)罐周邊建筑物對(duì)儲(chǔ)罐泄漏擴(kuò)散的影響,在儲(chǔ)罐下風(fēng)向50 m處設(shè)置障礙物(擋板),寬100 m,高25 m。

在LNG重氣云翻越障礙物之前,受擴(kuò)散及復(fù)雜風(fēng)場(chǎng)影響,氣云內(nèi)部濃度分布非常不均,會(huì)在儲(chǔ)罐與障礙物之間短期聚積,形成高濃度氣云,這是由于儲(chǔ)罐與障礙物的雙向阻擋限制了氣云的流動(dòng),導(dǎo)致氣云稀釋作用減弱。在LNG氣云翻越障礙物之后,大部分重氣云仍會(huì)受重力影響向下沉降,但其擴(kuò)散速度大大加快,云團(tuán)迅速向下風(fēng)向擴(kuò)展。根據(jù)圖6可知,障礙物除了能夠使中間部分重氣云沿障礙物向上攀爬,直至翻越障礙物,還能夠使得重氣云在順風(fēng)擴(kuò)散的過程中分流,呈現(xiàn)出“鉗”狀,“鉗”代表的兩個(gè)分支云團(tuán)會(huì)隨著LNG的泄漏擴(kuò)散不斷增大,最終融合為一體。

因此,重氣云繞過障礙物存在攀爬和分流兩種形式,為了防止泄漏形成的蒸汽云集聚在辦公行政樓、員工宿舍、控制室等引起危險(xiǎn),接收站建筑物建設(shè)布局必須經(jīng)過嚴(yán)格計(jì)算得出站內(nèi)各設(shè)施的最小安全距離,再進(jìn)行總平面布置設(shè)計(jì)。

a)50 s

b)70 s

2.2 單罐及罐區(qū)泄漏擴(kuò)散分析及對(duì)比

LNG泄漏擴(kuò)散在大氣環(huán)境中進(jìn)行,受到大氣流動(dòng)的影響,在單個(gè)儲(chǔ)罐周圍沒有其他障礙物的情況下,空氣從X負(fù)方向朝正方向流動(dòng),圖7顯示風(fēng)場(chǎng)分布僅受儲(chǔ)罐自身影響。儲(chǔ)罐迎風(fēng)面氣流部分沿儲(chǔ)罐輪廓向上抬升,其余向儲(chǔ)罐兩側(cè)分流,流線相對(duì)規(guī)則;儲(chǔ)罐背風(fēng)面氣流受到儲(chǔ)罐自身阻擋出現(xiàn)明顯的漩渦及回流,流線紊亂。

圖7 單罐風(fēng)場(chǎng)流線圖Fig.7 Wind field streamline figure of single tank

常溫常壓下空氣密度約為1.225 kg/m3,LNG從罐壁泄漏后與空氣進(jìn)行強(qiáng)烈快速的熱交換,部分氣化形成高濃度低溫蒸汽云,其初始密度約為空氣的1.5倍,從圖8顯示的氣云分布可見,在擴(kuò)散的同時(shí)受重力影響向地面沉降。由于儲(chǔ)罐背風(fēng)側(cè)氣流的漩渦及回流作用,部分氣云逆風(fēng)向罐壁處擴(kuò)散,導(dǎo)致罐壁壁面形成大量高濃度氣云,罐壁溫度降低。而大部分氣云受順流風(fēng)場(chǎng)的影響向遠(yuǎn)離儲(chǔ)罐的方向擴(kuò)散,因氣云對(duì)空氣有卷吸作用,高濃度氣云逐步擴(kuò)散為低濃度氣云,因而隨著時(shí)間的進(jìn)行,低濃度氣云的面積逐步增大,爆炸濃度區(qū)域增大。氣云濃度降低會(huì)導(dǎo)致其密度降低,當(dāng)與空氣密度相差較小甚至低于空氣密度時(shí),氣云受大氣浮力影響,在順風(fēng)擴(kuò)散的同時(shí)向上部空間流動(dòng)。

a)10 s

b)20 s

c)30 s

d)40 s

目前國(guó)內(nèi)LNG接收站罐區(qū)大多具備2座及2座以上的大型LNG儲(chǔ)罐,LNG泄漏擴(kuò)散過程中遇到的罐區(qū)內(nèi)部最大障礙物即是其他未發(fā)生泄漏的大型儲(chǔ)罐。以南方沿海某罐區(qū)分布為例,除發(fā)生泄漏的儲(chǔ)罐外還有3座相同儲(chǔ)罐。根據(jù)圖9的風(fēng)場(chǎng)分布情況可見,每座儲(chǔ)罐背風(fēng)側(cè)均會(huì)產(chǎn)生漩渦及回流,兩罐之間的風(fēng)場(chǎng)更因雙重阻擋導(dǎo)致紊流作用加強(qiáng),回流增多,流動(dòng)更加復(fù)雜。

圖9 罐區(qū)風(fēng)場(chǎng)流線圖Fig.9 Wind field streamline figure of tank farm

罐間空氣回流導(dǎo)致罐區(qū)1/2 LFL氣云在未到達(dá)下風(fēng)向儲(chǔ)罐時(shí)超前抬升,此后大量氣云沿罐壁被動(dòng)向上擴(kuò)散并受風(fēng)場(chǎng)影響繼續(xù)向下風(fēng)向運(yùn)動(dòng)。從圖10中兩個(gè)時(shí)刻的氣云分布來看,順風(fēng)向氣云由于儲(chǔ)罐的阻擋未到達(dá)罐區(qū)外部,但豎直向氣云卻因此高達(dá)一百多米,使得罐區(qū)危險(xiǎn)氣云體積增多,增加了發(fā)生火災(zāi)爆炸的危險(xiǎn),但從另一方面看,儲(chǔ)罐的阻擋延長(zhǎng)了危險(xiǎn)氣云到達(dá)罐區(qū)外部的時(shí)間。

泄漏80 s時(shí)1/2 LFL氣云分布見圖11,對(duì)比可見兩種幾何模型的模擬結(jié)果差異巨大,開敞空間單罐泄漏后的氣云體積遠(yuǎn)小于罐區(qū)氣云。一方面,罐區(qū)風(fēng)場(chǎng)比單罐風(fēng)場(chǎng)復(fù)雜得多,存在更多的回流,高濃度氣云停滯時(shí)間增加,氣云向下風(fēng)向擴(kuò)散的速度減慢;另一方面,下風(fēng)向儲(chǔ)罐對(duì)氣云的阻擋導(dǎo)致氣云沿儲(chǔ)罐向兩側(cè)擴(kuò)散,風(fēng)對(duì)橫向氣云擴(kuò)散的抑制作用減弱,氣云寬度不斷增加。以此可見,儲(chǔ)罐對(duì)氣云擴(kuò)散的阻擋作用非常大,罐區(qū)應(yīng)盡量遠(yuǎn)離周圍建筑,內(nèi)部安裝高敏感性報(bào)警器,以便能夠快速采取安全防護(hù)措施,在堵漏的同時(shí)利用水幕稀釋高濃度氣云。

a)40 s

b)60 s

a)單罐a)Single tank

b)罐區(qū)b)Tank farm

3 結(jié)論

1)風(fēng)能有效促進(jìn)氣云的擴(kuò)散,且對(duì)高濃度氣云產(chǎn)生的混合稀釋作用較強(qiáng),而對(duì)低濃度氣云產(chǎn)生的平流輸送作用較強(qiáng),在LNG泄漏及擴(kuò)散達(dá)到穩(wěn)定之前,風(fēng)速越大,UFL氣云到達(dá)下風(fēng)向距離越短,1/2 LFL氣云到達(dá)下風(fēng)向距離越遠(yuǎn)。在3 m/s低風(fēng)速下,1/2 LFL氣云在泄漏580 s時(shí)能夠到達(dá)下風(fēng)向最遠(yuǎn)距離554 m。

2)儲(chǔ)罐本身對(duì)氣云擴(kuò)散起阻擋作用,泄漏80 s時(shí)的1/2 LFL氣云體積從順風(fēng)24 229.9 m3,側(cè)風(fēng)22 680.9 m3,直至逆風(fēng)21 905.4 m3,阻礙作用逐漸加強(qiáng),氣云擴(kuò)散速度逐步減慢。因此,逆風(fēng)泄漏最不利于氣云擴(kuò)散,穩(wěn)定后形成的危險(xiǎn)區(qū)域最大。

3)障礙物阻擋重氣云擴(kuò)散,使得高濃度氣云短期聚積,增大火災(zāi)爆炸危險(xiǎn)性;LNG重氣云繞過障礙物存在攀爬和分流兩種形式,氣云短時(shí)間內(nèi)會(huì)呈現(xiàn)“鉗”狀。

4)LNG泄漏后初始?xì)庠泼芏却笥诳諝?呈重氣擴(kuò)散特征,同時(shí)風(fēng)場(chǎng)分布影響氣云擴(kuò)散,建筑分布越復(fù)雜,出現(xiàn)漩渦和回流的情況越多,高濃度氣云停滯時(shí)間越長(zhǎng),氣云擴(kuò)散減慢,因此,同一時(shí)刻單罐泄漏形成的危險(xiǎn)區(qū)域遠(yuǎn)小于罐區(qū)儲(chǔ)罐泄漏形成的危險(xiǎn)區(qū)域。