金春旭，潘振，陳保東，張桂強(qiáng)，黃騰龍，范開峰，李亮

（1. 遼寧石油化工大學(xué) 石油天然氣工程學(xué)院， 遼寧 撫順 113001； 2. 中國石油海洋工程有限公司天津分公司， 天津 300451）

液化天然氣儲(chǔ)罐翻滾事故分析

金春旭1，潘振1，陳保東1，張桂強(qiáng)2，黃騰龍1，范開峰1，李亮1

（1. 遼寧石油化工大學(xué) 石油天然氣工程學(xué)院， 遼寧 撫順 113001； 2. 中國石油海洋工程有限公司天津分公司， 天津 300451）

隨著世界能源緊缺現(xiàn)狀的出現(xiàn)，天然氣工業(yè)在我國的不斷發(fā)展，液化天然氣成為了一種成效卓著的儲(chǔ)運(yùn)方式，那么，隨乊而來的便是液化天然氣的安全問題，分層、翻滾是液化天然氣在儲(chǔ)存過程中有可能出現(xiàn)的一種極其危險(xiǎn)的事故現(xiàn)象。通過FLUENT模擬軟件對(duì)LNG儲(chǔ)罐中不同物性參數(shù)的液化天然氣在罐體漏熱情況下進(jìn)行數(shù)值模擬，在此乊前通過計(jì)算對(duì)大型LNG儲(chǔ)罐漏熱進(jìn)行求解，將計(jì)算得到的熱流密度近似取成10 W/m2，以此作為模擬時(shí)邊界條件的罐壁、罐底漏熱量，模擬得出翻滾發(fā)生的開始和結(jié)束時(shí)間，最后通過翻滾發(fā)生持續(xù)時(shí)間計(jì)算僅翻滾時(shí)所造成的LNG損失，約334 t，翻滾發(fā)生前，儲(chǔ)罐內(nèi)的氣相空間壓力梯度約為32 kPa/h，僅約20 min便會(huì)使儲(chǔ)罐內(nèi)正常操作壓力超出安全閥設(shè)定壓力。由此可見，液化天然氣翻滾事故的發(fā)生不僅會(huì)造成液化天然氣本身的一大浪費(fèi)，同時(shí)也會(huì)對(duì)人身及財(cái)產(chǎn)安全構(gòu)成巨大威脅，因此，深入研究LNG翻滾對(duì)社會(huì)發(fā)展具有深遠(yuǎn)意義。

液化天然氣；儲(chǔ)運(yùn)；事故；模擬；漏熱；分層；翻滾；損失

LNG工業(yè)作為當(dāng)今的新興能源正在蓬勃發(fā)展，被人譽(yù)為清潔、高效而又安全的新型能源。然而在儲(chǔ)存過程中也會(huì)隱藏著某種潛在危險(xiǎn)—翻滾。LNG翻滾是在 LNG儲(chǔ)存時(shí)發(fā)生的一種劇烈蒸發(fā)過程，短時(shí)間內(nèi)儲(chǔ)罐壓力由正常操作壓力上升到壓力安全閥設(shè)定壓力，安全閥跳起向外排放天然氣。如果操作不當(dāng)，翻滾不能及時(shí)得以控制，將引起安全閥頻繁開閉，造成大量低溫天然氣排放到儲(chǔ)罐上空空間，形成爆炸性混合云團(tuán)，或引起沸騰也爭(zhēng)氣膨脹爆炸，對(duì)LNG儲(chǔ)運(yùn)安全構(gòu)成重大威脅[1]。其中最大的兩起翻滾事故：一起發(fā)生在1971年的意大利，另一起發(fā)生在1993年的英國，兩起翻滾事故均在高峰時(shí)不到兩個(gè)小時(shí)的時(shí)間就造成100多噸的LNG損失。

LNG分層是翻滾事故發(fā)生的前提條件，而造成分層現(xiàn)象發(fā)生的原因可以歸納為以下三種：老化；新充注的LNG與儲(chǔ)罐剩余的LNG密度不同，幵超過臨界密度差；受到漏熱影響產(chǎn)生熱分層。LNG分層演化分為四個(gè)階段：固定分界面階段；上層卷攜階段；上層中心射流穿透階段；翻滾階段[2-4]。

1 儲(chǔ)罐罐體漏熱計(jì)算

漏熱是引起儲(chǔ)罐內(nèi)LNG產(chǎn)生分層一種因素，而分層又會(huì)造成翻滾事故發(fā)生，所以儲(chǔ)罐漏熱不可忽視，雖然 LNG儲(chǔ)罐具有良好的保冷性能，但仍會(huì)有漏熱發(fā)生[5,6]。本文以大連LNG接收站為例，大連 LNG接收站站址位于遼寧省大連市大孤山半島的東南面，項(xiàng)目分二期建設(shè)，一期工程設(shè)計(jì)規(guī)模為300×104t/a，計(jì)劃2011年4月1日投產(chǎn)，二期工程的設(shè)計(jì)規(guī)模增加到 600×104t/a，站內(nèi)擁有 3座16×104m3的LNG儲(chǔ)罐，單個(gè)LNG儲(chǔ)罐內(nèi)徑為40 m，罐高36 m，日均蒸發(fā)率為0.05%，充裝率為90%，正常操作壓力18 kPa，安全閥跳起壓力29 kPa。

式中：M—LNG儲(chǔ)罐小時(shí)蒸發(fā)量，kg/h；

r—日均蒸發(fā)率；

ρ—LNG的密度，kg/m3，取445 kg/m3；

V′—充裝率為90%時(shí)的LNG儲(chǔ)罐容積，m3。

式中：hg—單位質(zhì)量氣態(tài)甲烷的焓值，J/kg；

hl—單位質(zhì)量液態(tài)甲烷的焓值，J/kg；

Δ H—?dú)饣療?，J/kg；

Δ H1—甲烷的氣化熱，J/kg；

Tc—環(huán)境溫度，190.55 K；

Tl—標(biāo)準(zhǔn)液體天然氣溫度，111.65 K；

T—液化天然氣溫度，111.15 K；

n—系數(shù)，90.55 K

式中：q—通過LNG儲(chǔ)罐罐體漏熱量，W；

式中： φ—熱流密度，W/m2；

S1─LNG儲(chǔ)罐內(nèi)罐表面積，m2。

根據(jù)以上計(jì)算，選取該儲(chǔ)罐周壁漏熱熱流密度為10 W/m2，幵對(duì)模型罐壁及罐底同時(shí)設(shè)定為漏熱[7]。

2 翻滾模擬

對(duì)大型 LNG儲(chǔ)罐進(jìn)行數(shù)值模擬時(shí)，儲(chǔ)罐內(nèi)分為氣相區(qū)和液相區(qū)兩部分，而在密度差和溫度差的作用下液相區(qū)又分為分層區(qū)和主流區(qū)。密度大的為主流區(qū)在下層，密度小的為分層區(qū)在上層，兩層液體的自然對(duì)流分別被限制在各自的層區(qū)內(nèi)，因此，在層與層乊間會(huì)產(chǎn)生一個(gè)阻礙各自層區(qū)內(nèi) LNG進(jìn)入另一層的液—液界面，它對(duì)兩層液體的傳熱產(chǎn)生很大的影響[8,9]，如圖1所示。

為此，本文在使用FLUENT軟件模擬前對(duì)液化天然氣儲(chǔ)罐液相區(qū)進(jìn)行GAMBIT建模。二維模型中，儲(chǔ)罐壁面及罐底定義為壁面（WALL）邊界，由于該模型定義為軸對(duì)稱圖形，所以，左邊壁面定義為軸（AXIS），上邊由于是氣相區(qū)與液相區(qū)的交界面，故定義為壓力出口（PRESSURE-OUTLET）。

圖1 儲(chǔ)罐分層示意圖Fig.1 Schematic diagram of tank stratification

計(jì)算參數(shù)參考文獻(xiàn)[1]，如表1所示。

表1 主流層和分流層中液化天然氣物性參數(shù)Table 1 Physical parameters of main layer and the shunt layer liquefied natural gas

模擬結(jié)果如下圖所示：

圖 2（a）、（b）為翻滾開始階段的溫度場(chǎng)及密度場(chǎng)示意圖。

圖3（a）、（b）為翻滾結(jié)束階段的溫度場(chǎng)及密度場(chǎng)示意圖。

3 翻滾損失

3.1 翻滾過程的LNG損失

LNG翻滾事故是一種比較復(fù)雜的工況，在沒有具體工程設(shè)計(jì)模型的條件下，根據(jù)經(jīng)驗(yàn)公式（5）計(jì)算：

式中：M翻滾─LNG翻滾事故發(fā)生時(shí)小時(shí)最大蒸發(fā)量，kg/h。

圖2 翻滾開始階段示意圖Fig.2 Schematic diagram of roll beginning

圖3 翻滾結(jié)束階段示意圖Fig.3 Schematic diagram of the end of roll

代入數(shù)據(jù)計(jì)算得翻滾事故發(fā)生時(shí)每小時(shí) LNG最大蒸發(fā)量為133 500 kg。若根據(jù)1971年意大利LNG接收站發(fā)生的LNG翻滾事故高峰時(shí)間計(jì)算，1.25 h便可造成約170 t LNG蒸發(fā)排空[10]。

由模擬結(jié)果可知，翻滾持續(xù)時(shí)間約為2.5 h，因此可以到到翻滾發(fā)生過程所造成的LNG損失約334 t。

3.2 翻滾過程的罐內(nèi)壓力變化

由于儲(chǔ)罐內(nèi)發(fā)生翻滾時(shí)主要表現(xiàn)為物性參數(shù)不同的兩種LNG發(fā)生摻混及罐體漏熱導(dǎo)致的LNG迅速蒸發(fā)，使儲(chǔ)罐上空氣相空間操作壓力上升至安全閥設(shè)定壓力，安全閥跳起向外排放低溫天然氣，翻滾發(fā)生時(shí)，安全閥持續(xù)開啟，大量 LNG蒸發(fā)排空[11,12]。儲(chǔ)罐處于正常情況時(shí)，閃蒸氣量可以根據(jù)HT Hasheim/HR Wesson公式計(jì)算：

式中：M閃蒸—閃蒸氣量，lb/h；

S2—儲(chǔ)罐內(nèi)罐截面積，ft2；

ps—飽和蒸氣壓力，inchH2O。

若大氣壓或罐內(nèi)壓力發(fā)生變化，儲(chǔ)罐內(nèi)壓力為：式中：ps1—50%正常滿罐自然蒸發(fā)氣量時(shí)的飽和蒸氣壓力，inchH2O；

Ps2—LNG儲(chǔ)罐內(nèi)壓力變化后的壓力，inchH2O；—由LNG儲(chǔ)罐內(nèi)壓力變化引起的壓力梯度，inchH2O/h。

將上式中M閃蒸看作翻滾發(fā)生時(shí)的小時(shí)蒸發(fā)量，而由于發(fā)生翻滾時(shí)儲(chǔ)罐內(nèi)壓力急劇上升，此時(shí)儲(chǔ)罐內(nèi)飽和蒸汽壓相對(duì)較小，因此，代入計(jì)算可得 32 kPa/h，以此壓力梯度計(jì)算，僅約20 min便會(huì)使儲(chǔ)罐內(nèi)正常操作壓力超出安全閥設(shè)定壓力。

4 結(jié)束語

由于翻滾發(fā)生時(shí) LNG儲(chǔ)罐及罐內(nèi)液化天然氣所處狀態(tài)各不相同，如：儲(chǔ)罐內(nèi)兩層LNG密度差、溫度差、兩層 LNG厚度及漏熱等條件不同，導(dǎo)致翻滾發(fā)生的開始、結(jié)束時(shí)間、持續(xù)時(shí)間及罐內(nèi)壓力各不相同。通過計(jì)算分析得出處于以上狀態(tài)下的LNG儲(chǔ)罐翻滾結(jié)論如下：

（1）大型LNG儲(chǔ)罐漏熱熱流密度為10 W/m2；

（2）FLUENT模擬軟件模擬得到翻滾持續(xù)時(shí)間約為8 400 s，在此期間所造成的LNG損失約為334 t，幵在翻滾發(fā)生前安全閥的不斷開啟也造成很大一部分損失，因此LNG損失遠(yuǎn)不止334 t；

（3）LNG儲(chǔ)罐在發(fā)生翻滾前，儲(chǔ)罐內(nèi)的氣相空間壓力梯度約為32 kPa/h，僅約20 min便會(huì)使儲(chǔ)罐內(nèi)正常操作壓力超出安全閥設(shè)定壓力。

因此，LNG翻滾事故不容忽視，只有做好預(yù)防幵能夠在事故發(fā)生時(shí)采取適當(dāng)措施才能夠?qū)⑹鹿首钚』?，降低人身及?cái)產(chǎn)損失。

［1］侯金偉．LNG儲(chǔ)罐分層與翻滾的數(shù)值模擬[D]．哈爾濱：哈爾濱工業(yè)大學(xué)，2009．

［2］劉浩，周永春．LNG低溫儲(chǔ)罐壓力安全系統(tǒng)設(shè)計(jì)[J]．化工設(shè)計(jì)，2007，17（1）：7-10．

［3］姚云，張湘鳳，曲順利，等．LNG儲(chǔ)罐內(nèi)翻滾現(xiàn)象臨界判據(jù)的數(shù)值模擬[J]．低溫工程，2012，（6）：37-40．

［4］馬曉茜．LNG儲(chǔ)運(yùn)安全性問題的非線性機(jī)理及其智能評(píng)價(jià)[J]．天然氣工業(yè)，2003，23（6）：123-126．

［5］顧安忠．液化天然氣技術(shù)[M]．北京：機(jī)械工業(yè)出版社，2004：225．

［6］王冰，陳學(xué)東，王國平．大型低溫LNG儲(chǔ)罐設(shè)計(jì)與建造技術(shù)的新進(jìn)展[J]．天然氣工業(yè)，2010，30（5）：108-112．

［7］Q.-S. Chen a,b, J. Wegrzyn c,* , V. Prasad a. Analysis of temperature and pressure changes in liquefied natural gas (LNG) cryogenic tanks[J]．Cryogenics,2004,44:701-709．

［8］李霞，馬貴陽，杜明俊，等．LNG球罐分層翻滾汽化過程的數(shù)值分析[J]．遼寧石油化工大學(xué)學(xué)報(bào)，2010，30（4）：42-44．

［9］程旭東，徐劍，喬國發(fā)．大型液化天然氣儲(chǔ)罐內(nèi)液體內(nèi)分層和翻滾影響因素分析[J]．化工學(xué)報(bào)，2009，60（12）：2990-2996．

［10］康正凌，孫新征．LNG接收站蒸發(fā)氣量計(jì)算方法[J]．油氣儲(chǔ)運(yùn)，2011，30（9）：663-666．

［11］付子航．LNG接收站蒸發(fā)氣處理系統(tǒng)靜態(tài)設(shè)計(jì)計(jì)算模型[J]．天然氣工業(yè)，2011，31（1）：83-85．

［12］付子航．LNG接收站蒸發(fā)氣處理系統(tǒng)動(dòng)態(tài)設(shè)計(jì)計(jì)算模型[J]．天然氣工業(yè)，2011，31（6）：85-88．

Analysis of LNG Storage Tank Rollover Accidents

JIN Chun-xu1，PAN Zhen1，CHEN Bao-dong1，ZHANG Gui-qiang2，HUANG Teng-long1，F(xiàn)ANG Kai-feng1，LI Liang1
（1. College of Petroleum Engineering, Liaoning Shihua University, Liaoning Fushun 113001，China；2. China National Petroleum Offshore Engineering Co., Ltd., Tianjin 300451，China）

With emergence of the world's energy shortage situation, and the natural gas industry obtains continues development in China, liquefying natural gas becomes an effective way of storage and transportation, then safety problems in liquefying natural gas appear, stratification and rollover are extremely dangerous accident phenomena that may occur during LNG storage. In this paper, numerical simulation on storage tanks of liquefied natural gas with different physical parameters under the condition of heat leakage was carried out by FLUENT simulation software. Before the numerical simulation, heat leakage of the large LNG storage tank was calculated; calculated heat flux will be approximately determined 10 W/m2as the heat leakage of the boundary conditions of the tank wall. Then the start time and end time of roll occurring was got through simulation, and finally losses was calculated only under the roll occurrence in the roll time, about 334 t. The pressure gradient of vapor space of the tank is about 32 kPa/h before roll occurs, and about 20 minutes will make the normal operating pressure of the tank beyond the set pressure of safety valve. Thus, LNG rollover not only will result in a great loss of LNG, but also will pose a huge threat to the safety of persons and property, so in-depth study on LNG rollover has a far-reaching significance on the development of society.

Liquefied natural gas; Transportation; Accident; Simulation; Heat leakage; Stratification; Rollover; Losses

TE972

A

1671-0460（2014）10-2065-03

2014-03-05

金春旭（1988-），男，遼寧盤錦人，碩士生，2010年畢業(yè)于遼寧石油化工大學(xué)交通運(yùn)輸專業(yè)，研究方向：LNG衛(wèi)星站工藝設(shè)計(jì)。E-mail：jinchunxu987813@163.com。

指導(dǎo)教師： 潘振（1981-），男，副教授，博士，研究方向：油氣儲(chǔ)運(yùn)工程。E-mail：p6860770@126.com。