趙裕 呂植勇 劉俊雯 昌雪玲 沈繼超 陶瑋凡

摘要：為進一步提高液化天然氣（liquefied natural gas，LNG）罐柜內河運輸?shù)陌踩?，研究LNG罐柜在船上裝載位置的合理性。根據爆炸產生的火球模型和熱輻射模型，考慮人員受到不同程度傷害的概率，建立LNG罐柜爆炸危害的系統(tǒng)動力學模型。利用系統(tǒng)動力學仿真軟件Vensim_PLE模擬2個罐柜爆炸，得出結論：爆炸中心的熱輻射強度最高，但在距爆炸中心一定距離處人員受到傷害的概率最大。為保證人員安全，船舶罐柜裝載區(qū)與人員活動區(qū)應保持足夠的安全距離。

關鍵詞：? 安全管理工程; 安全距離; 系統(tǒng)動力學; LNG罐柜

中圖分類號：? U698

文獻標志碼：? A

Meeting of the Waterborne Transport Division， World Transport Convention 2021 （WTC 2021）

Safety study on loading position of LNG tanks on ships

ZHAO Yu1， LYU Zhiyong1， LIU Junwen1， CHANG Xueling2， SHEN Jichao1， TAO Weifan1

（1. National Engineering Research Center for Water Transportation Safety， Wuhan University of Technology， Wuhan 430000， China;

2. Zhejiang Scientific Research Institute of Transport， Hangzhou 310000， China）

Abstract： In order to further improve the safety of LNG tank transportation in inland rivers， the rationality of loading position of LNG tanks on ships is studied. According to the fireball model and the thermal radiation model， where the fireball and thermal radiation are produced by the explosion， the system dynamics model of the explosion hazard of LNG tanks is established considering the probability of different degrees of injury to the crew. Using the system dynamics simulation software Vensim_PLE to simulate the explosion of two tanks， it is concluded that the thermal radiation intensity of the explosion center is the highest， but the probability of crew injury is maximum at a certain distance from the explosion center. In order to ensure crew safety， a sufficient safe distance should be maintained between the tank loading area and the crew activity area on a ship.

Key words： safety management engineering; safety distance; system dynamics; LNG tank

收稿日期： 2021-03-29

修回日期： 2021-06-23

基金項目： 浙江省交通運輸廳科技計劃（2020043）

作者簡介：

趙裕（1995—），男，山東濟寧人，碩士研究生，研究方向為危化品運輸，（E-mail）895110093@qq.com;

呂植勇（1964—），男，湖北武漢人，教授，博士，研究方向為載運工具及運用，（E-mail）banlzy@whut.edu.cn

0 引 言

目前，我國液化天然氣（liquefied natural gas，LNG）罐柜水路運輸主要是海運，內河運輸尚未實施。2020年7月6日交通運輸部正式印發(fā)《整船載運液化天然氣可移動罐柜安全運輸要求（試行）》，使得推廣LNG罐柜內河運輸成為可能。在我國LNG罐柜內河運輸營運之前，首先需要解決LNG罐柜運輸?shù)陌踩詥栴}。LNG罐柜的水路運輸條件與陸路運輸條件存在較大差別，LNG罐柜水路運輸會更大程度地對載運船舶及其航行水域內的基礎設施、通航建筑的安全造成影響：一方面，載運船舶本身在航行、停泊、作業(yè)等營運過程中存在一定的風險;另一方面，若罐柜質量差則極易造成LNG泄漏，從而導致火災、爆炸、中毒、人身傷害等事故的發(fā)生，造成嚴重的生命財產損失。

我國LNG罐柜運輸?shù)难芯慷嗉性贚NG罐柜鐵路和公路運輸?shù)陌踩芯?。鄭偉[1]從運輸通道的比較入手，分析LNG罐柜江海聯(lián)運的必要性與可行性，認為相對于傳統(tǒng)運輸方式，罐柜將成為新形勢下的發(fā)展趨勢之一。船運LNG罐柜的運動是一種三維運動，LNG在罐柜內翻滾、晃動會使罐內壓力處在迅速上升與下降的變化過程中，存在一定的安全隱患[2-3]。2007年，彭建華[4]進行了我國第一次LNG罐柜水路運輸試驗，試驗結果為：在水路運輸過程中，罐內壓力緩慢增長，處于安全范圍內，罐體未出現(xiàn)損壞、腐蝕、變形等現(xiàn)象，然而該試驗時間僅為4 h，與實際船舶運輸時間存在較大差距。LNG罐柜水路運輸?shù)木嚯x遠、時間長，經長時間使用后其強度會逐漸降低，當降低至某一極限時，罐柜會破裂;由于罐內儲存有高壓氣體，罐內外的壓力差會使LNG迅速汽化;當存在點火源時，有可能燃燒，從而引起爆炸。LUKETA-HANLIN[5]在試驗中研究了儲存以及水路運輸LNG的過程中LNG泄漏引發(fā)火災的危險性，總結了LNG的特點及危害。對于LNG罐柜爆炸的研究，付珊珊[6]通過數(shù)值模擬分析了風速、大氣穩(wěn)定度、溫度等不同環(huán)境條件對LNG罐柜泄漏擴散的影響，得到：隨著風速的增大，LNG泄漏擴散的橫向影響范圍呈下降趨勢;隨著大氣穩(wěn)定度的上升，LNG泄漏擴散的橫向影響范圍呈上升趨勢;溫度變化對LNG泄漏擴散范圍的影響不是特別明顯。姜永勝等[7]和單彤文等[8]研究發(fā)現(xiàn)爆炸荷載能讓LNG罐體和穹頂產生位移變化。紀純明[9]使用TNT（trinitrotoluene）當量法，估算出LNG低溫罐柜全部泄漏引發(fā)爆炸所產生的沖擊波危害后果，認為LNG低溫儲罐的爆炸危害范圍較廣，且化學爆炸沖擊波的危害遠大于物理爆炸沖擊波的危害。管義鋒等[10]采用TNT當量法，對在LNG船加注過程中被加注氣罐泄漏進行了研究，預測了發(fā)生蒸氣云爆炸的危害范圍。伍蒙等[11]將液池蒸發(fā)過程考慮到蒸氣云計算模型中，使計算結果更加準確。以往的研究表明，LNG罐柜運輸安全性較高，但一旦發(fā)生LNG罐柜泄漏，引發(fā)爆炸，其危害巨大。前人對LNG泄漏后可能發(fā)生的爆炸危害研究大部分基于爆炸荷載以及TNT當量法計算爆炸危害，但是LNG蒸氣燃燒產生的火球和熱輻射帶來的破壞和危害更大。本文在前人研究的基礎上，改進火球和熱輻射模型，并從 LNG蒸氣燃燒產生的熱輻射強度和人員受到傷害的概率兩個角度出發(fā)，以目標（人員）與罐柜之間的距離為方程變量，采用系統(tǒng)動力學方法研究LNG罐柜爆炸的危害，對確定LNG罐柜在船上的裝載位置，提高運輸安全有重要的作用。

1 系統(tǒng)動力學評價模型

LNG罐柜使用周期長，但在船舶運輸中其腐蝕率和銹蝕率高，且罐柜之間發(fā)生碰撞的可能性較高，從而造成LNG罐柜的破損和泄漏。LNG爆炸會產生很強的熱輻射能量，而且一旦發(fā)生爆炸，船員無法迅速逃離船舶，會處于燃燒環(huán)境中，在離爆炸中心的不同距離處會受到不同程度的傷害。本文建立LNG罐柜爆炸危害的系統(tǒng)動力學模型（見圖1），從整體上判定安全距離，減少人員傷害。該模型可為合理制定LNG罐柜的裝載位置和設定隔離區(qū)域提供參考。

為便于理解各要素間的關系，使用FreeFta軟件，以死亡傷害概率為例，對各要素之間的關系進行分析。引發(fā)死亡傷害概率的事故樹見圖2。

2 系統(tǒng)動力學方程的建立及參數(shù)值的確定

LNG是深冷液體，泄漏后很容易受熱升溫成天然氣。LNG從周圍環(huán)境中吸收熱量，當溫度高于-107 ℃時，會擴散在空氣中。當蒸氣云的體積分數(shù)在燃燒范圍內時，如果遇到火源，哪怕是靜電，都會發(fā)生燃燒，產生巨大的危害。本文假設一艘運輸LNG罐柜的船在航行過程中由于碰撞等因素其2個LNG罐柜發(fā)生泄漏，并引發(fā)LNG燃燒。假設罐內LNG的密度ρ=426 kg/m3，罐柜體積V=250 m3，罐柜充裝系數(shù)k=95%。參考文獻[12-15]，將相關參數(shù)以系統(tǒng)的角度聯(lián)系起來，建立相關的系統(tǒng)動力學方程如下：

D=2.665M0.327

M=kρV

t=1.089M0.327

式中：D為LNG蒸氣燃燒產生的火球半徑，m;M為LNG蒸氣燃燒產生的火球中天然氣的質量，kg;t為LNG蒸氣燃燒產生的火球能夠持續(xù)燃燒的時間，s。

對于只有1個罐柜泄漏并引發(fā)LNG燃燒的情況，M可以取罐柜容量的50%;當有2個罐柜泄漏并引發(fā)LNG燃燒時，M取2個罐柜容量之和的70%;當存在多個罐柜泄漏并引發(fā)LNG燃燒時，M取所有罐柜總容量的90%。

火球在燃燒過程中，產生的熱量會向上升起一定的高度，本文假定火球中心距地面的高度（即火球抬升高度）h與LNG蒸氣燃燒產生的火球半徑D近似相等。

ESR=FSMHaπD2t

FS=0.27p0.32， p=1.21pV

Ha=HC-HV-CpT

式中：ESR為火球表面的熱輻射能量，J/kg;FS為火球表面的輻射能量比;Ha為火球的有效燃燒熱，J/kg;p為LNG的飽和蒸氣壓，MPa;pV為LNG罐柜安全閥的啟動壓力，本文取pV=1.25 MPa;HC為LNG的燃燒熱，J/kg;HV為沸點下LNG的蒸發(fā)熱，J/kg;Cp為恒壓狀態(tài)下LNG的比熱容，J/（kg·K）;T為火球表面與周圍環(huán)境間的溫差，一般取T=1 700 K。

kT=2.02（pWr1）-0.99

pW=pW0Hr

r1=r-D/2， r=（X2+h2）0.5

式中：kT為大氣熱傳遞系數(shù)，Pa·m;pW為環(huán)境溫度下空氣中的飽和水蒸氣壓，Pa;pW0為環(huán)境溫度下的飽和水蒸氣壓，Pa;Hr為相對濕度;r1為目標到火球表面的距離，m;r為目標到火球中心的距離，m;X為目標到罐柜的水平距離。設定初始值0為爆炸中心。

q=ESR4 153.8X2+16 615.2kT

式中：q為距離罐柜X處的熱輻射強度，J/（m2·s）。

qc=μq

式中：qc為有服裝防護下的熱輻射強度，J/（m2·s）;μ為有服裝防護下的熱接受率，一般取μ=0.4。

PDr=-43.14+3.018 8ln（L/W）

P（2）r=-36.38+2.56ln（L/W）

P（1）r=-39.83+3.018 6ln（L/W）

L=t*q4/3c

式中：PDr、P（1）r和P（2）r分別為人員受到死亡傷害、一度傷害、二度傷害的概率;L為熱負荷，W;t*為人員在熱輻射條件下暴露的時間，s。

3 數(shù)值模擬及數(shù)據分析

將各參數(shù)值輸入Vensim_PLE軟件進行模擬，取時間步長為1 s，時間持續(xù)100 s，研究各參數(shù)隨時間的變化趨勢，結果見圖3。

由圖3可以看出，隨著時間的增加，大氣熱傳遞系數(shù)和熱輻射強度會逐漸降低，人員受到傷害的概率和熱負荷先增加后降低。以死亡傷害概率為例，由圖2可以看出，人員受到死亡傷害的概率與熱負荷有直接關系。根據熱負荷公式L=t*q4/3c可得出，熱負荷變化與人員在熱輻射條件下暴露的時間和有服裝防護下的熱輻射強度成正比，熱負荷與有服裝防護下的熱輻射強度變化趨勢應相同。模擬中采用的是步進時間，故可認為t*是單位時間1 s，然而，圖3顯示，有服裝防護下的熱輻射強度與熱負荷的變化趨勢不同。根據圖2分析產生這種不同的原因，由事故樹可以看出，影響熱負荷的直接原因有時間和有服裝防護下的熱輻射強度，這與公式得到的結果相同，無法解釋產生這種不同的原因。利用Vensim_PLE中自帶的原因樹功能分析熱輻射強度變化的原因，見圖4。

由圖4可以看出，影響熱輻射強度變化的原因有大氣熱傳遞系數(shù)、目標到罐柜的水平距離（以下簡稱“距離”）以及火球表面的熱輻射能量。3種影響因素中，距離是存在的唯一變量。因此，距離通過影響熱輻射強度進而在一定程度上影響熱負荷大小，最終影響人員受到傷害的概率。距離與時間存在對應關系，人員受到傷害的概率也與時間存在對應關系。為便于研究人員受到傷害的概率與距離之間的關系，建立關于時間、距離和人員受到傷害的概率的三維坐標系，見圖5。

在這個三維坐標系中，x軸為時間，y軸為距離，z軸為人員受到傷害的概率。將三維坐標系轉動到y(tǒng)-z視圖下得到人員受到傷害的概率與距離的關系，見圖6。采用相同的分析方法，可以得到熱輻射強度與距離的關系，見圖7。

由圖6可知，人員受到傷害的概率隨距離的增加先增加后減小。對于雙罐爆炸模型來說，在距罐柜100 m左右處人員受到傷害的概率最大。在區(qū)間（0，100 m）與區(qū)間[100 m，400 m）內，在不同位置處人員受到低傷害的概率相同。然而，由圖7可見，隨著距離的增加，熱輻射強度降低，因此為確保安全，不同位置處危險概率相同時取較遠位置作為LNG罐柜裝載在船上的安全位置。

4 結束語

本文結合爆炸引起的火球模型、熱輻射模型和系統(tǒng)動力學方法，建立LNG罐柜爆炸危害系統(tǒng)動力學模型。利用MATLAB分析了時間、目標（人員）到罐柜的水平距離、人員受到傷害的概率、熱負荷和熱輻射強度之間的關系，為LNG罐柜在船上的裝載位置提供了新的理論和計算方法，結論如下：

（1）系統(tǒng)動力學建模方法將系統(tǒng)內不同參數(shù)之間的關系和相互作用通過相關反饋回路聯(lián)系起來，全面考慮各影響因素，分析過程簡單直接，提高了分析效率。

（2）距離爆炸中心越遠，LNG罐柜爆炸產生的熱輻射強度越低，穿防護服對有效降低傷害起重要作用。

（3）距爆炸中心一定距離處，人員受到傷害的概率最大。在不同位置處人員受到低傷害的概率相同，但距離較近的位置，熱輻射強度更高，因此應選擇距離較遠的地方進行人員疏散或設置船員居住區(qū)。

本文提供一種研究船運LNG罐柜裝載位置的方法，沒有根據船舶長度給出罐柜存儲位置建議。在以后的研究中將結合船舶實際長度對模型進行修改。

參考文獻：

[1]鄭偉. 液化天然氣罐式集裝箱江海聯(lián)運初探[J].? 四川水泥， 2017（10）： 280.

[2]FERRIN J L， PREZ-PREZ L J. Numerical simulation of natural convection and boil-off in a small size pressurized LNG storage tank[J]. Computers and Chemical Engineering， 2020， 138： 106840. DOI： 10.1016/J.COMPCHEMENG.2020.106840.

[3]LUO Min， WANG Xin， JIN Xin， et al. Three-dimensional sloshing in a scaled membrane LNG tank under combined roll and pitch excitations[J]. Ocean Engineering， 2020， 211： 107578. DOI： 10.1016/j.oceaneng.2020.107578.

[4]彭建華. LNG罐式集裝箱水運安全性探討[J]. 世界海運， 2007， 30（3）： 27-29. DOI： 10.16176/j.cnki.21-1284.2007.03.011.

[5]LUKETA-HANLIN A. A review of large-scale LNG spills： experiments and modeling[J]. Journal of Hazardous Materials， 2006， 132（2/3）： 119-140. DOI： 10.1016/j.jhazmat.2005.10.008.

[6]付姍姍. LNG罐式集裝箱內河運輸泄漏擴散模擬研究[D]. 武漢： 武漢理工大學， 2013.

[7]姜永勝， 蘇娟， 蘇龍龍， 等. 大型LNG儲罐爆炸載荷效應數(shù)值分析[J]. 石油和化工設備， 2020， 23（1）： 23-26.

[8]單彤文， 陳團海， 張超， 等. 爆炸荷載作用下LNG全容罐安全性優(yōu)化設計[J]. 油氣儲運， 2020， 39（3）： 334-341. DOI： 10.6047/j.issn.1000-8241.2020.03.013.

[9]紀純明. LNG低溫儲罐爆炸沖擊波危害后果估算[J]. 中國設備工程， 2020（9）： 84-85.

[10]管義鋒， 史騰飛， 朱培培， 等. LNG雙燃料船舶加注過程氣罐泄漏爆炸后果分析[J]. 船舶工程， 2016， 38（1）： 31-34. DOI： 10.13788/j.cnki.cbgc.2016.01.031.

[11]伍蒙， 許淵， 李左， 等. 改進的蒸氣云爆炸模型在LNG儲罐爆炸模擬中的應用[J]. 工業(yè)安全與環(huán)保， 2019， 45（10）： 44-47.

[12]范洪軍， 張暉， 徐建勇. LNG加注躉船的池火危險距離分析[J]. 中國造船， 2013， 54（4）： 186-195.

[13]LEE C， KIM B G， KO S Y. A new Lagrangian stochastic model for the gravity-slumping/spreading motion of a dense gas[J]. Atmospheric Environment， 2007， 41： 7874-7886. DOI： 10.1016/j.atmosenv.2007.06.007.

[14]別麗麗. LNG儲罐泄漏蒸氣擴散及爆炸危害研究[D]. 濟南： 山東建筑大學， 2019.

[15]GONG Maoqiong， MA Jia， WU Jianfeng， et al. Nucleate pool boiling of liquid methane and its natural gas mixtures[J]. International Journal of Heat and Mass Transfer， 2009， 52（11/12）： 2733-2739. DOI： 10.1016/j.ijheatmasstransfer.2008.12.011.

（編輯 趙勉）