劉春秀，易麗軍,2，陸小杰

(1 廣西大學(xué)資源與冶金學(xué)院， 廣西　南寧 530004；2 廣西高校礦物加工重點(diǎn)實(shí)驗(yàn)室， 廣西　南寧　530004)

?

某石化企業(yè)LPG球罐區(qū)泄漏擴(kuò)散的數(shù)值模擬

劉春秀1，易麗軍1,2，陸小杰1

(1 廣西大學(xué)資源與冶金學(xué)院， 廣西南寧 530004；2 廣西高校礦物加工重點(diǎn)實(shí)驗(yàn)室， 廣西南寧530004)

以某石油化工企業(yè)的LPG球罐罐區(qū)為研究對象，在對現(xiàn)場實(shí)際情況進(jìn)行簡化和抽象的基礎(chǔ)上，采用ANSYS FLUENT軟件對LPG球罐區(qū)的泄漏擴(kuò)散進(jìn)行三維數(shù)值模擬計算。主要研究了不同風(fēng)速條件下LPG擴(kuò)散的變化規(guī)律，得出了風(fēng)速對LPG球罐泄漏擴(kuò)散的擴(kuò)散方向、距離以及濃度產(chǎn)生的影響。該研究結(jié)果對球罐區(qū)的事故預(yù)防具有指導(dǎo)意義。

LPG；球罐；擴(kuò)散；FLUENT；數(shù)值模擬

隨著社會的發(fā)展，液化石油氣(LPG)作為重要的燃料和工業(yè)原料，已廣泛應(yīng)用于工業(yè)和民用的各個領(lǐng)域[1]。對于石化企業(yè)，儲存LPG的設(shè)備主要為儲罐，而球罐因其儲藏量大、耗材少、占地面積小、基礎(chǔ)工程簡單等特點(diǎn)，被廣泛應(yīng)用于石油、化工、冶金等企業(yè)[2]。LPG具有易燃、易爆等危險特性，在運(yùn)輸、儲存和使用的過程中易引發(fā)重大事故。尤其在LPG儲罐區(qū)，由于儲罐集中且儲存量大，一旦儲罐發(fā)生泄漏，帶有壓力的LPG從儲罐內(nèi)泄漏流出，立即劇烈汽化，發(fā)生閃蒸。若未及時采取有效措施，易引發(fā)火災(zāi)、爆炸事故，造成嚴(yán)重的人員傷亡和巨大的財產(chǎn)損失。由于LPG儲罐泄漏擴(kuò)散具有綜合性和復(fù)雜性，只靠實(shí)驗(yàn)難以全面、細(xì)致的對事故發(fā)生的機(jī)理進(jìn)行研究。因此，在實(shí)驗(yàn)研究的基礎(chǔ)上進(jìn)行數(shù)值模擬研究。

1　計算模型的建立

1.1罐區(qū)基本情況概述

某石化企業(yè)位于廣西欽州港工業(yè)區(qū)內(nèi)。儲罐區(qū)位于廠區(qū)的西北部，包括原料罐區(qū)、產(chǎn)品罐區(qū)和LPG球罐區(qū)。其中LPG罐區(qū)中有8個容積為1000 m3的LPG球罐，球罐主要設(shè)計參數(shù)及結(jié)構(gòu)尺寸見表1。

LPG的泄漏擴(kuò)散受到當(dāng)?shù)貧庀髼l件(風(fēng)速、風(fēng)向、氣溫、濕度以及氣壓等)的影響，該廠區(qū)的氣象條件見表2。

表1　球罐主要設(shè)計參數(shù)及結(jié)構(gòu)尺寸

表2　廠區(qū)的氣象條件

1.2數(shù)學(xué)物理模型的建立

1.2.1幾何模型

為便于模擬與分析，選取X、Y軸方向平行于水平地面，Z軸方向垂直于水平地面，使用Workbench的Geometry模塊建立LPG儲罐泄漏擴(kuò)散的三維幾何模型，如圖1所示。

圖1　LPG球罐區(qū)泄漏擴(kuò)散模擬計算模型

模擬計算區(qū)域?yàn)?25 m×125 m×50 m的長方體區(qū)域，區(qū)域內(nèi)包括8臺容量為1000 m3的LPG球罐(R1～R8)，球罐半徑均為6.15 m。球罐R1的中心位于(0,0,0)， X方向上相鄰兩球罐中心距離均為19.3 m，Y方向上相鄰兩球罐中心距離均為24.3 m，球罐底部距離地面高度為1.85 m。圖1中，面F1為風(fēng)流入口，面F2、F3、F5、F6均為氣體出口，面F4為地面。其中面F1到R1中心的距離為75 m，面F5到R1中心的距離為50 m。泄漏口位于球罐R1的底部，半徑為0.1 m，泄漏速度為250 m/s。氣象參數(shù)分別選?。簹鉁貫?00 K、大氣壓為標(biāo)準(zhǔn)大氣壓(101.325 kPa)、風(fēng)向?yàn)楸憋L(fēng)。設(shè)風(fēng)速方向?yàn)閄軸正方向。

利用ANSYS ICEM CFD軟件將計算域進(jìn)行非均勻的混合網(wǎng)格劃分，球罐周圍采用邊界層網(wǎng)格，并且為了提高模擬計算的精度，對泄漏口附近以及球罐區(qū)域的網(wǎng)格進(jìn)行加密處理。

1.2.2湍流模型

k-ε兩方程模型是用途最廣的湍流數(shù)值模型，包括標(biāo)準(zhǔn)k-ε模型、RNG k-ε模型和Realizable k-ε模型[3]。Realizable k-ε模型不僅比標(biāo)準(zhǔn)k-ε模型在質(zhì)量分?jǐn)?shù)分布上有著更好的精度，而且在預(yù)測由重力和氣體擴(kuò)散引起的質(zhì)量分?jǐn)?shù)波動時更符合試驗(yàn)結(jié)果[4]。因此本文采用Realizable k-ε模型，模型中關(guān)于k和ε的方程如下[5-6]:

(1)

(2)

式中，σk=1.0，σε=1.2，C2=1.9。

2　數(shù)值模擬結(jié)果與分析

圖2　不同風(fēng)速下地面的LPG質(zhì)量分?jǐn)?shù)等值線圖

LPG儲罐的泄漏擴(kuò)散受到當(dāng)?shù)貧庀髼l件的影響，其中風(fēng)速對LPG泄漏后擴(kuò)散具有很大的影響。因此，本文主要研究風(fēng)速對LPG泄漏擴(kuò)散的影響。分別對不同風(fēng)速(0 m/s、5 m/s、10 m/s、20 m/s、30 m/s五種情況)下LPG的泄漏擴(kuò)散情況進(jìn)行三維數(shù)值模擬，分析泄漏擴(kuò)散60 s時LPG氣云的擴(kuò)散規(guī)律，研究風(fēng)速對其擴(kuò)散的影響。圖2為不同風(fēng)速下地面的LPG質(zhì)量分?jǐn)?shù)等值線圖。

由圖2可以看出，無風(fēng)條件下，泄漏出的LPG以泄漏口為中心向四周擴(kuò)散，氣體擴(kuò)散以泄漏口為中心向周圍擴(kuò)散，并且質(zhì)量分?jǐn)?shù)大致呈對稱分布。有風(fēng)條件下，泄漏出的LPG隨著風(fēng)流向下風(fēng)向運(yùn)動，并緩慢的沿著下風(fēng)向擴(kuò)散，并且氣云出現(xiàn)分叉特征。風(fēng)流對擴(kuò)散方向產(chǎn)生了影響，將泄漏氣體向下風(fēng)向輸送。

同時風(fēng)流對擴(kuò)散距離也產(chǎn)生了影響：風(fēng)速越大，氣體擴(kuò)散速度越快，氣云在下風(fēng)方向擴(kuò)散的距離會隨著風(fēng)速的增加而增加。例如，風(fēng)速為10 m/s時，氣云能到達(dá)的最遠(yuǎn)距離明顯比風(fēng)速為5 m/s時的大。隨著風(fēng)速的不斷增大，當(dāng)風(fēng)速達(dá)到某一值時，風(fēng)對LPG氣云的稀釋作用成為影響LPG氣云擴(kuò)散的主要因素，此時LPG氣云在下風(fēng)方向可達(dá)到的距離反而會隨著風(fēng)速的增大而減小，如圖2所示，當(dāng)風(fēng)速為30 m/s時，擴(kuò)散的距離明顯比20 m/s時小。

圖3為不同風(fēng)速下沿X軸上地面LPG質(zhì)量分?jǐn)?shù)的變化曲線。

圖3　不同風(fēng)速下沿X軸上地面LPG質(zhì)量分?jǐn)?shù)的變化曲線

由圖3可以看出，隨著沿X方向距離的增加，LPG的質(zhì)量分?jǐn)?shù)逐漸減小。這是由于隨著風(fēng)速的增大，大氣的湍流程度也在增大，空氣的稀釋作用也增強(qiáng)，泄漏的LPG越容易被稀釋，這樣一來就會降低下風(fēng)向氣體的質(zhì)量分?jǐn)?shù)，且風(fēng)速越大，下風(fēng)向同一位置的質(zhì)量分?jǐn)?shù)越低。例如，下風(fēng)向X=20 m處LPG的質(zhì)量分?jǐn)?shù)，風(fēng)速為5 m/s時約為3.25%， 10 m/s時約為1.25%， 20 m/s時約為0.40%， 30 m/s時約為0.09%。

圖4為不同風(fēng)速下X=8 m處沿Y方向上地面LPG質(zhì)量分?jǐn)?shù)的變化曲線。

圖4　不同風(fēng)速下X=8 m處沿Y方向上地面LPG質(zhì)量分?jǐn)?shù)的變化曲線

由圖4可以看出：在Y方向上，以泄漏孔為中心，隨著兩側(cè)離泄漏口距離的增加，LPG的質(zhì)量分?jǐn)?shù)逐漸減小，并且風(fēng)速越大，質(zhì)量分?jǐn)?shù)減小的速度越快，氣體在Y方向上擴(kuò)散的距離越小。這是由于隨著泄漏的進(jìn)行，氣云變得越來越寬，氣云兩側(cè)受到風(fēng)流的影響越來越大，尤其是氣云最外側(cè)的部分，由于質(zhì)量分?jǐn)?shù)小密度輕，受風(fēng)流的影響最嚴(yán)重，使得氣云兩側(cè)在順風(fēng)向的擴(kuò)散速度明顯高于中間部分，而向橫風(fēng)向的擴(kuò)散速度越小，故擴(kuò)散的距離越小。

圖5為不同風(fēng)速下 X=8 m,Y=0 m處垂直高度上LPG質(zhì)量分?jǐn)?shù)的變化曲線。

圖5　不同風(fēng)速下X=8 m,Y=0 m處垂直高度上LPG質(zhì)量分?jǐn)?shù)的變化曲線

由圖5可知：隨著垂直地面方向上高度的增加，LPG質(zhì)量分?jǐn)?shù)逐漸降低，當(dāng)高度增加到一定值時，LPG質(zhì)量分?jǐn)?shù)變化較小。在靠近地面的位置，LPG的質(zhì)量分?jǐn)?shù)要比遠(yuǎn)離地面位置的質(zhì)量分?jǐn)?shù)高，這說明LPG主要沿地面擴(kuò)散。這是由于LPG比空氣重，泄漏后會向下沉降，擴(kuò)散沿地面移動。因此，在地面附近，人員的活動頻繁，并且焊接作業(yè)、摩擦靜電等點(diǎn)火源出現(xiàn)的可能性高，引起火災(zāi)爆炸的危險性高，因此LPG儲罐泄漏擴(kuò)散的危害范圍主要是在近地面區(qū)域。

3　結(jié)　論

本文運(yùn)用計算流體力學(xué)ANSYS FLUENT軟件對LPG球罐區(qū)泄漏擴(kuò)散情況進(jìn)行了三維數(shù)值模擬，研究了不同風(fēng)速情況下對LPG泄漏擴(kuò)散的影響，結(jié)果表明：

(1)風(fēng)速對LPG泄漏的擴(kuò)散方向、擴(kuò)散距離和濃度分布產(chǎn)生影響。風(fēng)流將泄漏氣體向下風(fēng)向輸送，并緩慢沿下風(fēng)向擴(kuò)散。隨著離泄漏口距離的增加，LPG的質(zhì)量分?jǐn)?shù)逐漸減小。風(fēng)速越大，所形成的氣云向下風(fēng)向擴(kuò)散的距離越大，質(zhì)量分?jǐn)?shù)減小的速度越快。但當(dāng)風(fēng)速達(dá)到一定值時，由于空氣的稀釋，氣云的擴(kuò)散距離反而減小。

(2)LPG屬于重氣，主要沿地面擴(kuò)散。LPG的質(zhì)量分?jǐn)?shù)隨著垂直高度的增加而逐漸降低。

[1]孫文棟,方江敏.LPG儲罐火災(zāi)爆炸事故后果模擬方法研究[J].工業(yè)安全與環(huán)保,2009, 35(12):44-47.

[2]居建龍,張明.LPG球罐的全面檢驗(yàn)和缺陷分析[J].科技傳播,2011(6):105-106.

[3]姜昉.易燃?xì)怏w泄漏擴(kuò)散過程的數(shù)值模擬研究[D].北京:首都經(jīng)濟(jì)貿(mào)易大學(xué),2008.

[4]S.M. Tauseef, D. Rashtchian, S. A. Abbasi. CFD-based simulation of dense gas dispersion in presence of obstacles [J]. Journal of Loss Prevention in the Process Industries, 2011,24: 371-376.

[5]王建.儲罐區(qū)可燃?xì)怏w泄漏擴(kuò)散模擬及爆燃災(zāi)害評估[D].大連:大連理工大學(xué),2013.

[6]寧平,張朝能,沈武燕.危險品泄漏的風(fēng)洞實(shí)驗(yàn)與數(shù)值模擬[M].北京:冶金工業(yè)出版社,2010:80-84.

Numerical Simulation of LPG Spherical Tank’s Leakage Diffusion in A Petrochemical Enterprise

LIUChun-xiu1,YILi-jun1,2,LUXiao-jie1

(1 College of Resources and Metallurgy, Guangxi University, Guangxi Nanning 530004;2 Key Laboratory of Mineral Processing in Guangxi University, Guangxi Nanning 530004, China)

Taking some petrochemical enterprise of LPG spherical tank and tank farm as the research object, based on the actual situation to simplify and abstract, using ANSYS FLUENT software for leakage of LPG spherical tank area spread 3-d numerical simulation, the change rule of LPG diffusion at different wind speeds was mainly studied. It was concluded that the wind speed had an effect on the LPG spherical tank leakage diffusion direction, distance, and concentration. The results of the study have guiding significance for the prevention of tank area accident.

LPG; spherical tank; diffusion; FLUENT; numerical simulation

劉春秀(1989-)，女，廣西大學(xué)安全科學(xué)與工程專業(yè)碩士研究生，主要研究方向?yàn)楣I(yè)安全。

X937;TE88

A

1001-9677(2016)010-0182-03