基于FLACS的水上LNG泄漏事故數值模擬分析

程 康，金全洲，范洪軍

(中國船級社武漢規范研究所，武漢430022)

基于FLACS的水上LNG泄漏事故數值模擬分析

程 康，金全洲，范洪軍

(中國船級社武漢規范研究所，武漢430022)

為準確評估水上LNG泄漏事故的危害以降低事故造成的損失,采用挪威三維計算流體力學(CFD)軟件FLACS進行水上LNG泄漏事故數值模擬,分析數值模擬中相關場景及參數設定,為水上LNG泄漏事故的仿真計算提供參考。

數值模擬;液化天然氣;計算流體力學;泄漏擴散;FLACS

液化天然氣(LNG)替代傳統燃油作為船舶動力燃料已引起各界關注,并成為業內公認的未來船舶綠色能源的首選。鑒于LNG的深冷和可燃特性,水上LNG安全應用與發展的首要問題就是LNG泄漏事故的風險控制。早期的研究基于唯象模型(經驗模型)進行無因次處理擬合出解析公式,該方法簡單易行但精度較差［1-2］。箱及相似模型［3］具有概念清晰、計算量較小等優點,適合于數值模擬分析,但模型需假定速度和濃度的自相似分布,并且通常涉及不連續界面,具有很大的不確定性,不能模擬障礙物周圍或者復雜地形上的流動擴散。近期比較成熟的方法有美國海岸警衛隊開發的DEGADIS大氣擴散模型［4］和計算流體力學(CFD)模型,三維CFD模型［5］是基于納維-斯托克斯方程組(N-S方程)的流體守恒定律模型,具有較強的科學依據,可以提供最詳細的關于LNG蒸氣云在大氣中流動和擴散的描述,具有模擬復雜地形情形氣云擴散的能力。也正因為如此,基于N-S方程的CFD模型被視為研究LNG泄漏蒸氣云擴散的最好工具,代表了該領域研究的方向。三維CFD模型計算通常采用Fluent軟件完成,LNG為低溫液體,一旦發生泄漏,液態LNG急劇蒸發形成蒸氣,同時吸收大量的熱,使周圍空氣降溫,空氣和蒸氣混合形成重氣云團,這時蒸氣云擴散形式為重氣擴散行為,因此必須要考慮重力對擴散行為的影響。當LNG泄漏時間較長或泄漏量較大時,泄漏出的LNG射流來不及全部蒸發,將會在地上或水面上形成液池,液池的擴展、蒸發、消失等行為將會使LNG泄漏擴散過程非常復雜［6］,采用通用軟件在參數設定上極為不便。挪威三維CFD軟件FLACS在LNG泄漏模擬方面具有權威性,并且是美國聯邦法規CFR惟一認可用于模擬LNG泄漏擴散的CFD軟件。該軟件將大量參數內置,在進行水上LNG泄漏事故仿真計算時具有良好的操作性和準確性,但目前國內相關的研究較少。為此,介紹采用FLACS軟件進行水上LNG泄漏擴散模擬的過程,以期推進國內對該方法的參考和應用。

1 事故分析

L.Vandebroek和J.Berghmans對水上LNG泄漏事故進行了研究,給出泄漏事故樹［7］,見圖1。

圖1 LNG泄漏事故樹

水上LNG泄漏后,會對人體造成低溫灼傷、對船體結構造成冷脆,擴散LNG蒸氣云的爆炸極限為5%～15%［8］,遇到點火源后,會形成閃火、池火,在圍壁處所內點燃則會出現蒸氣云爆炸。因此,水上LNG泄漏事故數值模擬的重點在于確定泄漏后形成的液池直徑和濃度處于爆炸極限內的蒸氣云范圍。

2 仿真計算

2.1 網格選取

圖2 網格大小與計算精度、時間的關系

當網格大于1 m時,計算精度下降很快,同時選取1 m和1 m以下網格時精度基本沒有差異,計算時間卻相差很大,綜合可知進行LNG泄漏擴散分析時選擇核心區域網格大小為1 m為宜。

2.2 泄漏文件設置

通過FLACS外置液池泄漏文件*.POOL可以方便地在其液池擴展模塊中設定泄漏速率和泄漏時間,進而給出泄漏量。該泄漏速率可以是定制速率,也可以是時間的變化量。例如,對于一個前10 s為定值泄漏速率(300 kg/s),后10 s線性減小到零的泄漏場景,可編輯液池泄漏文件見圖3。

另外還可通過FLACS內置的leak wizard來設定泄漏量,見圖4。

在指定泄漏形式、初始泄漏體積、初始壓力、噴嘴直徑、液體溫度等參數后,leak wizard會自動生成相關的泄漏文件。使用中發現leak wizard切斷泄漏控制的3個條件(迭代次數,切斷壓力和泄漏總量)并列,即滿足其中任意一項則泄漏被終止。

圖3 利用POOL文件生成液池效果

圖4 FLACS中的leak w izard

用leak wizard模擬不同方向噴射泄漏的效果,見圖5。

利用POOL文件設置方便,計算效率高,但前提需已知泄漏量,且噴射方向無法設置;利用leak wizard設置,需要準確把握初始參數和切斷參數,計算時間較長,但計算精度和模擬效果較POOL文件好。

3 實例分析

以內河21 000 m3LNG運輸船為例,其單罐容積為7 500 m3,初始壓力為0.4 MPa,液體溫度為-138℃［10］,泄漏時間1 800 s,風速3 m/s［11］,穩態D,核心區域網格尺寸1 m,利用POOL文件設置泄漏量。

圖5 不同方向的泄漏噴射模擬

計算結果與比利時大型LNG基礎設施公司Fluxys2012年的研究報告對比見表1。

表1 兩種結果對比

文中計算結果與Fluxys研究報告基本吻合。

4 結論

1)水上LNG泄漏事故主要造成的危害有低溫損傷、池火、閃火、噴射火和蒸氣云爆炸,為將事故后續傷亡降至最低,需對其進行數值模擬,確定泄漏形成的液池直徑和可燃蒸氣云范圍,以制定相應保護措施;

2)采用FLACS進行水上LNG泄漏事故分析時,為同時保證計算的精度和效率,建議核心區域網格大小取為1 m;

3)FLACS對于泄漏量的設定有兩種方式, POOL文件簡捷、計算速度快,但需提供泄漏速率,且不能模擬復雜泄漏情況;leak wizard能依據輸入參數自動計算泄漏速率,且可模擬多種形式的泄漏,計算準確度高,但設置較復雜,計算速度慢,建議結合實際需要選擇;

4)FLACS軟件能夠很好地對水上LNG泄漏事故進行模擬和分析,并能將評價結果以圖表形式直觀地表達出來。相較于傳統方法,具有參數設置便捷,計算速度快,精度高等優點。但是要合理正確地使用該軟件須結合工程實際,對模擬過程中每一參數都要加以推敲,以確保每一輸入項的合理性,從而使模型能更準確地反映工程事故實際。

［1］KOOPMAN,RONALD PE,DONALD L.Lessons learned from LNG safety research［J］.Journal of Hazardous Materials,2007,140(3):412-428.

［2］VDI.Dispersion of heavy gas emissions by accidental release-safety study［M］.Germany:Verein Dutscher Ingeniere,1990.

［3］KUNSCH J P,WEBBER D M.Simple box model for dense-gas dispersion in a straight sloping channel［J］. Journal of Hazardous Materials,2000,75(1):29-46.

［4］潘旭海,華 敏,蔣成軍.環境條件對LNG泄漏擴散影響的模擬研究［J］.天然氣工業,2009,29(1):117 -119.

［5］RICCIARDI L P R,VOST C,BOUILLOUX L.Experimental and numerical study of heavy gas dispersion in a ventilated room［J］.Journal of Hazardous Materials, 2008,152(2):493-505.

［6］莊學強.大型液化天然氣儲罐泄漏擴散數值模擬［D］.武漢:武漢理工大學,2012.

［7］VANDEBROEK L,BERGHMANS J.Safety aspects of the use of LNG formarine propulsion［J］.Procedia Engineering,2012,45:21-26.

［8］JOHN L,Woodward,Robin M,Pitblado.LNG risk based safety modeling and consequence analysis［M］.Hoboken:John Wiley&Sons,2010.

［9］GexCon.FLACS manual［M］.Norway:GexCon AS, 2014.

［10］Rue Guimard.Safety study for LNG supply in Flemish ports［R］.Brussels:Fluxys LNG sa,2012.

［11］段春鋒,繆啟龍.中國風速變化特征及其影響因素［C］∥中國科協常委會青年工作專門委員會、國務院學位委員會辦公室、中國科協組織人事部.第八屆博士生學術年會論文摘要集.中國科協常委會青年工作專門委員會、國務院學位委員會辦公室、中國科協組織人事部,2010.

Numerical Simulation for LNG Leakage Accident of Ship Based on FLACS

CHENG Kang,JIN Quan-zhou,FAN Hong-jun
(Wuhan Rules and Research Institute of China Classification Society,Wuhan 430022,China)

To assess accurately the consequences thatmay be led by liquefied naturalgas(LNG)leakage accidentof ships, numerically simulation for LNG leakage accidentofa ship is implemented based on the computational fluid dynamics(CFD)software FLACS.It is intended as a reference for simulation research relating to LNG leakage accident overwater by analyzing in detail the scenario and parameter setting in process of simulation.

numerical simulation;liquefied natural gas;computational fluid dynamics;leakage;FLACS

10.3963/j.issn.1671-7953.2015.03.040

U698

A

1671-7953(2015)03-0169-04

2015-01-30

修回日期:2015-03-05

程 康(1986-),男,碩士,工程師

研究方向:船舶與海洋結構物設計制造

E-m ail:kcheng@ccs.org.cn