周軍 李傳錢 蒲鶴 梁光川

1.西南石油大學(xué)石油與天然氣工程學(xué)院 2.中國(guó)石油化工股份有限公司青島安全工程研究院

LNG主要成分為甲烷，儲(chǔ)存運(yùn)輸?shù)臏囟葹?162 ℃，空氣中的爆炸濃度極限約為5%～15%，具有低溫、易燃易爆等特性。槽車運(yùn)輸為L(zhǎng)NG主要的陸地運(yùn)輸方式，在槽車LNG裝卸作業(yè)過程中，通過裝車棧臺(tái)的鶴管與槽車通過快速接頭進(jìn)行鏈接，快速接頭的安全和穩(wěn)定性極大程度地影響了整個(gè)LNG裝卸作業(yè)的安全和效率。由于長(zhǎng)時(shí)間低溫高壓下工作，接頭的密封圈極易損壞，導(dǎo)致LNG泄漏，LNG裝卸作業(yè)時(shí)，相關(guān)人員距離危險(xiǎn)源近，一旦發(fā)生LNG泄漏，后果將不堪設(shè)想。在LNG儲(chǔ)存運(yùn)輸環(huán)節(jié)中，王志寰等[1]對(duì)LNG接收站大型儲(chǔ)罐低溫出液管(150 mm)泄漏爆炸進(jìn)行模擬，采用相對(duì)偏差率對(duì)不同環(huán)境下LNG泄漏擴(kuò)散進(jìn)行分析和研究，預(yù)測(cè)了災(zāi)害效應(yīng)。而凌曉東[2]對(duì)接收站中槽車裝卸作業(yè)時(shí)管線(80 mm)斷裂造成的LNG泄漏的事故場(chǎng)景進(jìn)行模擬，分析了不同環(huán)境條件下可燃?xì)庠频挠绊懛秶?jì)算爆炸產(chǎn)生的超壓最大值。國(guó)外學(xué)者研究LNG船LNG泄漏事故較多，F(xiàn)ilippo Gavelli等[3]對(duì)裝卸碼頭LNG船在LNG裝卸作業(yè)過程中發(fā)生LNG泄漏爆炸事故進(jìn)行研究，分析部分受限和部分擁擠的障礙物中蒸氣云爆炸情形，對(duì)事故潛在后果進(jìn)行了評(píng)價(jià)。Sydney D Ryan等[4]模擬了LNG船的LNG泄漏(孔徑12 mm)擴(kuò)散情形，對(duì)比分析LNG蒸氣云在有無實(shí)體障礙物的環(huán)境下的危險(xiǎn)區(qū)域范圍，并指出障礙物的存在，降低了可燃范圍。

國(guó)內(nèi)外關(guān)于LNG接收站LNG泄漏的研究，多數(shù)都是LNG快速接頭脫落或者管道斷裂等大型LNG泄漏事故的研究和風(fēng)險(xiǎn)評(píng)價(jià)[5-10]，少有對(duì)小孔徑或少量LNG泄漏事故進(jìn)行研究和評(píng)價(jià)。為此，研究采用CFD技術(shù)，以LNG裝卸作業(yè)快速接頭密封圈損壞或失效為泄漏工況，通過FLACS軟件模擬該泄漏工況下，LNG迅速泄漏和汽化向周圍擴(kuò)散的過程，深入分析LNG泄漏擴(kuò)散規(guī)律以及燃燒爆炸風(fēng)險(xiǎn)，為現(xiàn)場(chǎng)LNG裝卸作業(yè)LNG泄漏事故應(yīng)急預(yù)案編制提供理論依據(jù)，對(duì)裝車區(qū)域LNG泄漏風(fēng)險(xiǎn)安全監(jiān)控和應(yīng)急響應(yīng)管理具有實(shí)際意義。

1 氣體擴(kuò)散數(shù)學(xué)模型

LNG泄漏到大氣時(shí)，會(huì)迅速汽化形成氣云，向四周擴(kuò)散。對(duì)于氣體擴(kuò)散，F(xiàn)LACS數(shù)值模擬采用有限體積法在三維笛卡爾坐標(biāo)系下，求解可壓縮流體流動(dòng)的質(zhì)量、動(dòng)量、能量和輸運(yùn)方程組成的守恒方程，見式(1)[11]。

(1)

式中：t為泄漏時(shí)間，s；ρ為氣體密度，kg/m3；xi為在i方向的坐標(biāo)，m；ui為i方向上的速度矢量，m/s；φ為通用求解變量，單位為對(duì)應(yīng)求解變量的國(guó)際制單位；Гφ為擴(kuò)散系數(shù)，m2/s；Sφ為定義的相應(yīng)求解變量φ源項(xiàng)。

相對(duì)于穩(wěn)定的層流，LNG泄漏形成的氣云擴(kuò)散流動(dòng)是一種復(fù)雜的不定常隨機(jī)流動(dòng)[12]，湍流模型更加適用于該情形下的流動(dòng)。湍流模型使用目前廣泛應(yīng)用的標(biāo)準(zhǔn)k-ε模型，使得守恒方程組封閉可解，方程如式(2)和式(3)所示。

Gk+Gb-ρε-YM+Sk

(2)

式中：k為湍動(dòng)能，m2/s2；ε為湍動(dòng)能耗散率，m2/s3；μ為動(dòng)力黏度，Pa·s；μt為動(dòng)力湍流黏度，Pa·s；σk為湍動(dòng)能k對(duì)應(yīng)的普朗特?cái)?shù)，取1.0；σε為湍動(dòng)能耗散率ε對(duì)應(yīng)的普朗特?cái)?shù)，取1.3；Gk為平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng)，kg/(m·s3)；Gb為浮力引起的湍動(dòng)能k的產(chǎn)生項(xiàng)，kg/(m·s3)；C1ε、C2ε和C3ε為經(jīng)驗(yàn)常數(shù)，分別取1.44、1.92和0.80；Sk為定義的湍動(dòng)能源項(xiàng)，kg/(m·s3)；Sε為定義的湍動(dòng)能耗散率源項(xiàng)，kg/(m·s4)；YM為可壓湍流中脈動(dòng)擴(kuò)張，kg/(m·s3)。

2 LNG泄漏擴(kuò)散CFD模型建立與驗(yàn)證

2.1 Falcon系列實(shí)驗(yàn)簡(jiǎn)述

Falcon系列實(shí)驗(yàn)是美國(guó)國(guó)家重點(diǎn)實(shí)驗(yàn)室-勞倫斯利物莫實(shí)驗(yàn)室(LLNL)，以研究障礙物對(duì)LNG泄漏擴(kuò)散的阻擋及LNG泄漏氣云在復(fù)雜障礙物下的擴(kuò)散行為為目的所進(jìn)行的系列現(xiàn)場(chǎng)大型實(shí)驗(yàn)。該系列實(shí)驗(yàn)總共進(jìn)行5次LNG泄漏試驗(yàn)，LNG通過1個(gè)多出口溢出的“蜘蛛”式溢流管線，溢出到矩形水池(60 m×40 m×0.76 m)中?！爸┲搿笔揭缌鞴芫€由4只臂組成，每只臂與相鄰的臂相距90°，單臂長(zhǎng)11.6 m，出口直徑為0.114 m和0.038 m(僅4號(hào)實(shí)驗(yàn))，LNG在水面上方直接向下釋放(見圖1)。水池中的水循環(huán)系統(tǒng)保持水溫穩(wěn)定，使LNG蒸發(fā)率近似等于溢出率。該試驗(yàn)區(qū)域?yàn)橐粋€(gè)圍欄封閉區(qū)域，圍欄長(zhǎng)88 m、寬44 m、高8.7 m，圍欄中高13.3 m、寬17.1 m的屏障位于水池的上風(fēng)。該現(xiàn)場(chǎng)實(shí)驗(yàn)中部署多個(gè)氣體含量傳感器，傳感器分別位于圍欄下風(fēng)50 m、150 m和250 m處，高度分別為1 m、5 m、11 m和17 m。Falcon系列實(shí)驗(yàn)的LNG溢流速率為8.7～30.3 m3/min，溢出量為20.6～66.4 m3，溢流時(shí)間為78～301 s?，F(xiàn)場(chǎng)在2 m高度測(cè)得的風(fēng)速范圍為1.7～5.2 m/s，大氣條件為中性至穩(wěn)定。

2.2 LNG泄漏擴(kuò)散模型

通過FLACS-CASD，以Falcon系列實(shí)驗(yàn)現(xiàn)場(chǎng)布置為參考，1∶1建立如圖2所示障礙物下的LNG泄漏擴(kuò)散模型。阻礙LNG擴(kuò)散的障礙物為圍欄和屏障，4個(gè)LNG泄漏孔分別位于水池中央的4段溢流管線的末端，視LNG泄漏速率與汽化速率相同，6個(gè)含量監(jiān)測(cè)點(diǎn)布置如圖3所示。

因Falcon系列實(shí)驗(yàn)中2號(hào)和5號(hào)實(shí)驗(yàn)存在不同程度的數(shù)據(jù)缺失，所以本次驗(yàn)證選取實(shí)驗(yàn)報(bào)告中Falcon系列1號(hào)、3號(hào)和4號(hào)實(shí)驗(yàn)相關(guān)數(shù)據(jù)[14]，對(duì)有障礙物的情況下LNG泄漏擴(kuò)散模型進(jìn)行定性和定量評(píng)價(jià)驗(yàn)證，相關(guān)現(xiàn)場(chǎng)實(shí)驗(yàn)參數(shù)見表1。本次建立的LNG泄漏擴(kuò)散模型，所有泄漏及環(huán)境參數(shù)均參照表1設(shè)置。

表1 Falcon系列1、3、4號(hào)實(shí)驗(yàn)參數(shù)表實(shí)驗(yàn)編號(hào)單孔泄漏直徑/mm管線壓力/MPa單孔泄漏體積流量/(m3·min-1)總泄漏體積流量/(m3·min-1)總泄漏質(zhì)量流量/(kg·s-1)泄漏時(shí)間/s平均風(fēng)速/(m·s-1)環(huán)境溫度/℃大氣穩(wěn)定度風(fēng)向泄漏孔坐標(biāo)1號(hào)114.30.4487.17528.72101311.732.8G(6.9, -2.4)3號(hào)114.30.2764.72518.91381544.134.9D(9.0, 0.7)4號(hào)38.10.8622.1758.7643015.231.1D/E(9.2, -1.3)(-37.5, -7.5, 0.5)(-22.5, -7.5, 0.5)(-22.5, 7.5, 0.5)(-37.5, 7.5, 0.5)

2.3 擴(kuò)散模型有效性定性評(píng)價(jià)

經(jīng)過模擬，得到3組數(shù)值模擬結(jié)果，因篇幅原因，僅取3組實(shí)驗(yàn)的不同點(diǎn)位的模擬預(yù)測(cè)值和Falcon實(shí)驗(yàn)值進(jìn)行有效性定性評(píng)價(jià)，分別取Falcon 1號(hào)實(shí)驗(yàn)的M3和M4點(diǎn)、Falcon 3號(hào)實(shí)驗(yàn)的M5和M6點(diǎn)，以及Falcon 4號(hào)實(shí)驗(yàn)的M2和M3點(diǎn)的時(shí)間-含量數(shù)據(jù)[13]，繪制數(shù)值模擬預(yù)測(cè)值與Falcon實(shí)驗(yàn)值曲線進(jìn)行對(duì)比，結(jié)果分別如圖4～圖6所示。分析時(shí)間-含量曲線可以看出，在有障礙物的情況下，1號(hào)和3號(hào)實(shí)驗(yàn)LNG擴(kuò)散時(shí)間-含量曲線的模擬值與現(xiàn)場(chǎng)實(shí)驗(yàn)吻合，4號(hào)實(shí)驗(yàn)泄漏量較小，泄漏時(shí)間最長(zhǎng)，LNG擴(kuò)散趨勢(shì)和含量也完全和現(xiàn)場(chǎng)實(shí)驗(yàn)結(jié)果相吻合，定性說明了建立的模型的有效性。

2.4 擴(kuò)散模型有效性定量評(píng)價(jià)

一般定量評(píng)價(jià)需對(duì)各監(jiān)測(cè)點(diǎn)的每個(gè)時(shí)間點(diǎn)的含量數(shù)據(jù)進(jìn)行統(tǒng)計(jì)偏差分析，但由于現(xiàn)場(chǎng)實(shí)驗(yàn)與數(shù)值模擬無法精確統(tǒng)計(jì)兩者間在各個(gè)相同時(shí)間點(diǎn)的含量值，此方法可行性低。為了定量評(píng)價(jià)該氣體擴(kuò)散模型的有效性，本研究采用LNG泄漏擴(kuò)散模型偏差評(píng)價(jià)標(biāo)準(zhǔn)進(jìn)行評(píng)價(jià)[14]。該標(biāo)準(zhǔn)是通過對(duì)相對(duì)平均偏差(MRB)、幾何平均偏差(MG)、幾何方差(VG)、相對(duì)均方誤差(MRSE)、歸一化均方誤差(NMSE)和FAC2參數(shù)進(jìn)行計(jì)算和判斷，達(dá)到評(píng)價(jià)目的。評(píng)價(jià)中分別統(tǒng)計(jì)各監(jiān)測(cè)點(diǎn)的最大濃度值，然后計(jì)算3組實(shí)驗(yàn)的數(shù)值模擬值和Falcon實(shí)驗(yàn)值的各個(gè)偏差[13]，得到偏差統(tǒng)計(jì)評(píng)價(jià)結(jié)果，如表2所列。結(jié)果表明，F(xiàn)alcon系列實(shí)驗(yàn)與模擬預(yù)測(cè)值的各項(xiàng)評(píng)價(jià)偏差值均在合理有效的范圍之內(nèi)，因此可認(rèn)為建立的模型是有效的。經(jīng)過定性和定量評(píng)價(jià)，證明建立的氣體擴(kuò)散模型的有效性。

表2 模擬預(yù)測(cè)值與Falcon系列實(shí)驗(yàn)實(shí)測(cè)偏差值統(tǒng)計(jì)結(jié)果評(píng)價(jià)參數(shù)NMSEMRBMGVGMRSEFAC2評(píng)價(jià)標(biāo)準(zhǔn)<4.00-0.40～0.400.67～1.50<3.30<2.300.50～2.00理想值001101Falcon 10.080.271.321.080.080.76Falcon 30.040.171.181.030.030.85Falcon 40.01-0.110.901.010.011.11有效性合理合理合理合理合理合理

3 裝卸作業(yè)LNG泄漏擴(kuò)散數(shù)值模擬與分析

3.1 裝車區(qū)三維模型

以現(xiàn)場(chǎng)實(shí)際情況為參考，使用FLACS-CASD建立裝車區(qū)三維模型，如圖7所示。裝車區(qū)由5輛槽車、裝車罩棚等組成，罩棚長(zhǎng)40 m、寬22 m、高8 m。裝卸作業(yè)時(shí)，鶴管與槽車連接的部位位于車后，當(dāng)LNG快速接頭發(fā)生密封圈失效時(shí)，泄漏點(diǎn)位于槽車后方，泄漏點(diǎn)坐標(biāo)為(0, 0, 0.8)，如圖7所示。為監(jiān)測(cè)LNG泄漏后不同位置的可燃?xì)怏w的含量變化規(guī)律，在離地面高度0.65 m，沿槽車左右和車后依次共設(shè)置15個(gè)氣體含量監(jiān)測(cè)點(diǎn)，M1-M5位于泄漏點(diǎn)上風(fēng)向，M6-M10位于泄漏點(diǎn)下風(fēng)向，M11-M15位于槽車后方(見圖7)，坐標(biāo)具體參數(shù)如表3所列。

表3 監(jiān)測(cè)點(diǎn)坐標(biāo)參數(shù)m監(jiān)測(cè)點(diǎn)坐標(biāo)監(jiān)測(cè)點(diǎn)坐標(biāo)監(jiān)測(cè)點(diǎn)坐標(biāo)M1(0, -1.5, 0.65)M6(0, 1.5, 0.65)M11(1.5, 0, 0.65)M2(0, -3, 0.65)M7(0, 3, 0.65)M12(3, 0, 0.65)M3(0, -4.5, 0.65)M8(0, 4.5, 0.65)M13(4.5, 0, 0.65)M4(0, -6, 0.65)M9(0, 6, 0.65)M14(6, 0, 0.65)M5(0, -7.5, 0.65)M10(0, 7.5, 0.65)M15(7.5, 0, 0.65)

3.2 泄漏參數(shù)

LNG裝卸作業(yè)時(shí)，LNG接頭密封圈失效形式有密封圈缺失、密封圈斷裂、密封圈磨損等[15]。在本次LNG泄漏數(shù)值模擬中，考慮最嚴(yán)重的泄漏情況為密封圈完全缺失，其次為密封圈部分缺失。

LNG裝卸作業(yè)過程中，介質(zhì)溫度為-162 ℃，鶴管接頭處管內(nèi)流壓為0.7 MPa，管徑為DN50，泄漏時(shí)間為120 s，對(duì)于不同密封圈失效工況，分析計(jì)算得出相應(yīng)的泄漏速度，取不同的環(huán)境參數(shù)，得到不同工況下泄漏基礎(chǔ)參數(shù)(見表4)。

表4 裝卸作業(yè)LNG泄漏參數(shù)項(xiàng)目泄漏工況等效泄漏孔徑/mm風(fēng)速/(m·s-1)泄漏質(zhì)量流量/(kg·s-1)離地高度/m大氣穩(wěn)定度TEST 1密封圈部分缺失4.73無風(fēng)0.200.8FTEST 2密封圈完全缺失6.69無風(fēng)0.40 0.8FTEST 3密封圈部分缺失4.732 0.20 0.8FTEST 4密封圈完全缺失6.692 0.400.8F

3.3 數(shù)值模擬結(jié)果及分析

3.3.1監(jiān)測(cè)點(diǎn)含量

通過模擬，得到4組15個(gè)點(diǎn)位不同泄漏工況下可燃?xì)怏w含量隨時(shí)間變化曲線，分別如圖8～圖11所示。在無風(fēng)環(huán)境下，LNG泄漏擴(kuò)散最初20 s，因空氣阻礙，擴(kuò)散前緣形成LNG氣云堆積，含量高。隨著LNG持續(xù)泄漏，氣云持續(xù)向周圍擴(kuò)散，泄漏后30 s，監(jiān)測(cè)點(diǎn)含量達(dá)到穩(wěn)定。120 s泄漏停止，在大氣作用下，氣體繼續(xù)擴(kuò)散，含量迅速下降。在有風(fēng)環(huán)境下，氣體擴(kuò)散的不定常性和隨機(jī)性明顯，含量波動(dòng)性大。受到風(fēng)場(chǎng)影響，上風(fēng)向可燃?xì)怏w含量極低，下風(fēng)向可燃?xì)怏w含量很高。

相同位置處可燃?xì)怏w含量與密封圈缺失程度和所處風(fēng)場(chǎng)位置有關(guān)，密封圈缺失程度越嚴(yán)重，位于下風(fēng)向，距離泄漏點(diǎn)越近，可燃?xì)怏w含量越高。

在LNG裝卸作業(yè)時(shí)，發(fā)生快速接頭密封圈缺失造成泄漏，以5%(體積分?jǐn)?shù)，下同)為爆炸下限(LEL)，在0.65 m的高度上，各個(gè)不同泄漏工況下造成的危險(xiǎn)氣體LEL擴(kuò)散最遠(yuǎn)距離也不同。分析時(shí)間-含量曲線，得到如圖12所示的0.65 m高度不同泄漏工況下氣體擴(kuò)散最遠(yuǎn)距離。泄漏量越大，環(huán)境風(fēng)速越大，氣體LEL擴(kuò)散距離越遠(yuǎn)。

3.3.2擴(kuò)散含量分布

分別取4組數(shù)值模擬泄漏擴(kuò)散過程中的40 s、80 s、120 s和140 s時(shí)刻，地面高度為0.5 m時(shí)可燃?xì)怏w體積分?jǐn)?shù)為5%～15%的分布圖，進(jìn)行密封圈失效工況下LNG泄漏擴(kuò)散含量分布規(guī)律分析。

無風(fēng)環(huán)境下，密封圈部分缺失和完全缺失LNG泄漏擴(kuò)散的規(guī)律相同，所以僅對(duì)最危險(xiǎn)的泄漏工況，即密封圈完全缺失情形下LNG泄漏模擬結(jié)果進(jìn)行分析，可燃?xì)怏w含量隨時(shí)間變化結(jié)果如圖13所示。由圖13可以看出，在無風(fēng)條件下，泄漏過程前40 s，LNG泄漏后，氣云向四周以環(huán)狀均勻擴(kuò)散。80 s后，由于槽車和罩棚區(qū)等障礙物的存在，LNG泄漏后形成的重氣云沿著槽車底部擴(kuò)散時(shí)，氣云運(yùn)移方向受到限制，導(dǎo)致槽車底部、車頭和車尾可燃?xì)怏w集中，含量高，危險(xiǎn)程度高。

有風(fēng)環(huán)境下，密封圈部分缺失和完全缺失LNG泄漏擴(kuò)散的規(guī)律相同，所以僅對(duì)密封圈完全缺失數(shù)值模擬結(jié)果進(jìn)行分析，數(shù)值模擬結(jié)果如圖14所示。因?yàn)榇嬖陲L(fēng)場(chǎng)，40 s時(shí)，泄漏點(diǎn)左邊裝卸車區(qū)域，可燃?xì)怏w擴(kuò)散受到阻礙，擴(kuò)散速度慢，右邊無障礙，擴(kuò)散速度快，危險(xiǎn)區(qū)域范圍右邊大于左邊。80 s后，在泄漏點(diǎn)左邊槽車和裝卸車區(qū)域，由于障礙物的存在，使得風(fēng)速減緩，大氣湍流復(fù)雜，裝車區(qū)域可燃?xì)怏w被稀釋速度小于右邊無障礙區(qū)，使得裝車區(qū)的危險(xiǎn)區(qū)域比右邊范圍更廣。

對(duì)比分析不同環(huán)境因素下的模擬結(jié)果可以發(fā)現(xiàn)，在有風(fēng)條件下，大氣對(duì)氣云的稀釋非常明顯，有效降低了整體危險(xiǎn)區(qū)域范圍。風(fēng)在通過障礙物時(shí)，其稀釋和運(yùn)移氣云的能力降低，使得裝車區(qū)局部可燃?xì)怏w含量較高，危險(xiǎn)程度高。

通過模擬結(jié)果可以看出，在0.5 m的高度上，不同泄漏工況下的危險(xiǎn)含量氣體擴(kuò)散范圍不同?？偨Y(jié)擴(kuò)散規(guī)律分布，得到不同泄漏工況下氣體LEL擴(kuò)散最大范圍，如圖15所示。在密封圈完全缺失，無風(fēng)條件下，可燃?xì)怏w擴(kuò)散的范圍最廣，LNG泄漏后擴(kuò)散形成一個(gè)直徑約為58.0 m可燃?xì)怏w危險(xiǎn)燃燒爆炸區(qū)域。

3.3.3可燃?xì)怏w云團(tuán)及TNT當(dāng)量

通過數(shù)值模擬，得到4種不同工況下的等量化學(xué)式可燃?xì)怏w云團(tuán)體積變化規(guī)律和最大值。從圖16可以看出，在無風(fēng)環(huán)境下，從泄漏開始可燃?xì)怏w云團(tuán)量直線增長(zhǎng)，約127 s后停止增長(zhǎng)，迅速下降。在有風(fēng)環(huán)境下，氣云的增加速度較慢，約125 s后開始下降。

在4種不同泄漏工況下，于120 s泄漏停止后，可燃?xì)怏w云團(tuán)量均在繼續(xù)增長(zhǎng)，形成最大可燃?xì)怏w云團(tuán)量的時(shí)間均在120 s以后的一段時(shí)間內(nèi)。原因在于，120 s后擴(kuò)散仍在持續(xù)，大氣稀釋的速度小于云團(tuán)自身擴(kuò)散的增加速度，氣云量繼續(xù)增加，一段時(shí)間后，可燃?xì)怏w云團(tuán)量增加速度開始小于大氣稀釋速度，氣云量開始下降。

LNG裝卸作業(yè)時(shí)，因密封圈損壞導(dǎo)致LNG泄漏，形成的可燃?xì)怏w云團(tuán)量會(huì)受到環(huán)境風(fēng)速和密封圈損壞狀態(tài)的影響。環(huán)境越穩(wěn)定，密封圈損失越嚴(yán)重，形成的可燃?xì)怏w云團(tuán)量越大。其中，在無風(fēng)環(huán)境下，密封圈缺失工況形成的可燃?xì)怏w云團(tuán)量最大，為278.7 m3。按照常用的TNT當(dāng)量法計(jì)算[16]，計(jì)算得到不同泄漏工況下的可燃?xì)怏w云團(tuán)TNT當(dāng)量(見表5)，最大可燃?xì)怏w云團(tuán)(甲烷體積分?jǐn)?shù)為5%～15%)TNT當(dāng)量為7.99～23.98 kg。

表5 可燃?xì)怏w云團(tuán)體積和TNT爆炸當(dāng)量項(xiàng)目LNG泄漏總量/kg可燃?xì)怏w云團(tuán)總量/m3TNT爆炸當(dāng)量/kgTEST 124 93.52.68～8.05TEST 248278.77.99～23.98TEST 324 21.60.62～1.86TEST 448 64.6 1.85～5.56

4 結(jié)論

(1) 通過Falcon現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)與數(shù)值模擬結(jié)果的定性及定量評(píng)價(jià)，驗(yàn)證了建立的LNG泄漏擴(kuò)散模型的有效性。

(2) LNG裝卸作業(yè)時(shí)，接頭處密封圈缺失導(dǎo)致LNG發(fā)生泄漏后，相同位置處可燃?xì)怏w含量與密封圈缺失程度和環(huán)境因素有關(guān)，密封圈缺失越嚴(yán)重，泄漏點(diǎn)同側(cè)距離泄漏點(diǎn)越近，可燃?xì)怏w含量越高。在0.65 m高度上，有風(fēng)環(huán)境下，密封圈完全缺失，氣體LEL擴(kuò)散距離最遠(yuǎn)，其距離超過8 m。

(3) 受到泄漏量和環(huán)境因素影響，密封圈缺失越嚴(yán)重，大氣條件越穩(wěn)定，氣體擴(kuò)散范圍越廣。0.5 m高度上，密封圈完全缺失，無風(fēng)環(huán)境下，氣體LEL擴(kuò)散范圍最大，區(qū)域直徑約為58 m。

(4) 形成的可燃?xì)怏w云團(tuán)體積受到密封圈缺失程度和環(huán)境因素影響，無風(fēng)環(huán)境下，密封圈完全缺失，形成的可燃?xì)怏w云團(tuán)體積最大，體積約為178.7 m3，TNT爆炸當(dāng)量最大約為24 kg。

(5) LNG泄漏擴(kuò)散，具有重氣擴(kuò)散特征。高度越低，擴(kuò)散范圍越廣；在有風(fēng)條件下，氣體擴(kuò)散濃度波動(dòng)加強(qiáng)，擴(kuò)散方向不定常性和隨機(jī)性越明顯，有效降低了危險(xiǎn)區(qū)域范圍。裝車區(qū)障礙物較多，空間狹窄流動(dòng)性差，氣體聚集，燃燒爆炸的危險(xiǎn)程度最高。