泄漏工况下LNG燃爆事故时空演化特征

张洋杰　周瑶　何华刚　丁金金　张方舟　佟国强

摘要：為研究泄漏工况下液化天然气的时空演化规律，基于CFD和FVM方法，通过FLACS软件，对不同风速条件下LNG接收站泄漏的扩散及蒸气云爆炸发展规律进行了研究。另外，在距火焰中心不同距离下，采用FDS火灾仿真软件对LNG接收站持续泄漏时引发的池火灾进行了模拟。结果表明，无风条件下初期LNG扩散区域大部分为低于爆炸浓度下限（5%）、但高于0.94%的泄漏浓度。因此在距离泄漏点较近的敏感区内必须考虑无风情况下泄漏危险浓度对于敏感区的影响;无风工况下，LNG泄漏危险浓度区域和爆炸蒸气云浓度横向扩散范围更大;有风工况下，泄漏的LNG形成的危险浓度区域向敏感区扩散速度更快，垂直天空扩散的高度更高，爆炸蒸气云浓度沿风速方向蔓延距离更远;另外无风条件下，175 s左右时距火焰中心30 m的地方辐射热流密度达到最大90 kW/m²，火焰中心距离与辐射热流密度成反比，揭示了泄漏工况下LNG燃爆事故时空演化特征，为LNG接收站的选址、站内布局和储罐泄漏应急处置提供理论依据和支撑。

关键词：泄漏;液化天然气;时空演化;蒸气云爆炸;池火灾

中图分类号：TE 972;TE 88文献标志码：A

文章编号：1672-9315（2022）04-0637-10

DOI：10.13800/j.cnki.xakjdxxb.2022.0403

Spatial and temporal evolution characteristics of LNG

explosion accident under leakage conditionsZHANG Yangjie ZHOU Yao HE Huagang DING Jinjin ZHANG Fangzhou TONG Guoqiang

（1.Research Center of Emergency Management Department，Beijing 100013，China;

2.Faculty of Engineering，China University of Geosciences，Wuhan 430074 China;

3.Shenzhen GDT Security Technology Co.，Shenzhen 518000，China;

4.Saifeite Engineering Technology Group Co.，Ltd.，Qingdao 266061，China）Abstract：The spatial and temporal evolution law of LNG under leakage conditions was studied by CFD and FVM methods.And the diffusion and vapor cloud explosion development law of LNG receiving station leaks under different wind speed conditions were examined by FLACS software.At different distances from the center of the flame，pool fires caused by the continuous leakage of the LNG receving stations were simulated by FDS fire simulation software.The results indicated that most of the initial LNG diffusion area under windless conditions was below the lower explosive concentration limit（5%）but above the leakage concentration of 0.94%.Therefore，in the sensitive areas close to the leakage point，the impact of the hazardous concentration of the leak on the sensitive area under windless conditions must be taken into consideration.Under windless conditions，the dangerous concentration area of LNG leakage and the lateral diffusion range of explosion vapor cloud concentration are larger;under windy conditions，the area of dangerous concentration formed by the leaking LNG spreads faster to the sensitive areas，the height of vertical sky diffusion is higher，and the concentration of explosive vapor clouds spreads farther along the wind speed direction.In addition，under windless conditions，the radiant heat flux density at a distance of 30 m from the center of the flame reaches a maximum of 90 kW/m2 at about 175 s，and the distance from the center of the flame is inversely proportional to the radiant heat flow.The results reveal the spatial and temporal evolution characteristics of LNG flaring and explosion accidents under leakage conditions，which provides a theoretical basis and support for the site selection of LNG receiving stations，station layout and emergency disposal of storage tank leakage.

Key words：leakage;Liquefied Natural Gas;spatial and temporal evolution;vapor cloud explosions;pool fires

0引言

随着中国“碳达峰”和“碳中和”目标的实施，天然气将成为其他化石能源的最佳替代品之一。而目前，中国天然气需求主要来源于进口，进口天然气的途径共有2种：一种是通过跨国天然气管道输送，以中俄东线天然气管道为代表，另一种则为进口液化天然气[1-2]（liquefied natural gas，LNG）。而液化天然气储气站作为管道沿线的重要组成，在储运调峰、燃气发电等方面发挥着重要的作用。其中LNG储罐是储气站的主要储存容器，LNG的储存方式为低压、低温，但由于LNG液化生产的过程中不可避免地会混入H₂O，H₂S，HCL等气体，而这些气体在露点温度下会变成腐蚀性液体，从而形成酸性腐蚀环境，对罐体产生腐蚀作用[3-4]。一旦发生泄漏，根据泄漏情形的不同，会引发蒸气云爆炸、池火灾、喷射火等燃爆事故，易形成多米诺效应，造成大面积火灾，导致严重的经济损失和人员伤亡[5-6]。因此，研究LNG泄漏扩散的时空演化规律，准确预测LNG泄漏对人员财产的影响范围，对应急处置和救援具有十分重要的意义。

国内外学者针对LNG泄漏的研究，目前主要集中在理论、相似模拟实验、数值模拟等方面。其中理论方面主要以模型为主，包含有高斯、唯象（即经验）、BM（bench mark）、箱型相似、浅层理论模型等[7-9]。相似模拟实验主要有Maplin Sands系列实验、Burro系列实验、Coyote系列实验、Falcon系列实验、Meroney实验以及阿肯色大学危险化学品研究中心进行的一系列风洞试验等[10-12]。数值模拟方面主要利用Fluent，FLACS（flame acceleratioin simulation）、FDS（fire dynamic simulation）等软件进行模拟研究。如：ZHANG等通过数值模拟不同风向LNG泄漏情况，发现LNG在顺风水平方向扩散最远[13]。MARSEGAN等基于数值模拟发现主动设置屏障可有效有效减小LNG泄漏后的扩散范围[14]。NGUYEN等在2017年首次通过试验手段发现LNG的泄漏率与液池扩展速度和蒸发速率呈正相关，同时根据一维热传导模型建立了蒸发速率与泄漏率和时间的函数关系[15]。

综上所述，目前对LNG泄漏扩散的研究主要侧重于不同大气环境对LNG扩散距离的影响[16-20]，但针对因LNG泄漏导致引发连锁反应，发生蒸气云爆炸、池火灾等燃爆事故的演化特征研究较少[21-24]。因此，文中结合相关风险报告和标准，以中国广州深圳大鹏接收站为研究对象，综合考虑其泄漏点位置、泄漏速率、风速、风向以及大气环境温度等多种因素，对LNG接收站泄漏的扩散过程、蒸气云爆炸和油池火灾进行了模拟，深入分析了因LNG接收站泄漏引发连锁反应导致蒸气云爆炸与油池火灾的时空演化特征，为LNG接收站的选址、站内布局和储罐泄漏应急处置提供理论依据和支撑。

1模型构建

1.1FLACS数值模型

LNG接收站泄漏与蒸气云爆炸的数值模拟均基于CFD（computational fluid dynamics）和FVM（finite volume method）方法，采用FLACS软件进行模拟。FLACS对气相模型的计算主要基于以下3個流体力学基本公式[25]。

流体力学连续性方程

1.2FDS火灾数值模型

1.3工况参数

LNG接收站泄漏扩散模型环境参数参考当地的气象资料，泄漏管道、LNG的理化性质、场区防火堤的高度等参数设置参考《深圳市大鹏新区油气库企业和高压管道安全风险咨询报告》《深圳市空港油料有限公司成品油经营单位安全评价报告》《空港油料公司重大危险源评估报告》《空港油库改造说明》等。参数设置见表1。

点火源位置如图1白色圆点所示，具体坐标为（-47.5，153.5，41.1），分别对应XYZ轴。由图2可知，广东大鹏LNG接收站的敏感区包括迭福路、东部电厂和码头卸油区，因此针对该区域，对该场区内1号罐顶部阀门处，分别进行了无风条件与风速为5 m/s的模拟研究，总网格数为789 276，最小网格尺寸为1 m，场区模型如图3所示;考虑因泄漏引发连锁反应导致池油火灾、蒸气云爆炸事故的发生，为人员疏散以及应急预案的编写提供依据，因此将风速设置1.3 m/s，泄漏速率为145.92 m/s，风向为西北风，大气稳定度F级。

2计算实例

2.1LNG接收站泄漏模拟分析

在不同风速条件下通过对LNG接收站泄漏模拟，发现LNG持续泄漏扩散过程具有明显的分阶段特征，同时甲烷扩散范围在不同风速条件下随时间的变化也有明显的不同。

从图4可以看出，泄漏点在无风条件下扩散初期呈椭圆形扩散，图中蓝色区域大部分为低于爆炸浓度下限（5%），但高于0.94%的泄漏浓度。此浓度范围以上，人员穿戴和佩戴防护用具时可不受伤害，反之亦然，因此该区域范围为危险区域范围。随时间的增长椭圆形扩散的径向尺寸逐渐增大，当泄漏时间到达200 s左右时椭圆形危险区域范围开始发生破裂，并向四周蔓延扩散，此时甲烷燃烧爆炸上限（UFL）距离、爆炸下限（LFL）距离和容许接触浓度（0.5LFL）安全距离增加。

在无风条件下，泄漏口设置在1号罐顶部，泄漏口方向正对东部电厂敏感区位置，且储罐内外存在较大压力差，LNG泄漏浓度较大，在气压的作用下泄漏出来的LNG首先会向东部电厂敏感区方向径向扩散，释放气压，在气压得到释放后，再继续横向扩散。泄漏发生第50 s时，泄漏的LNG已经蔓延至场区边缘;在泄漏发生第100 s时，泄漏的LNG造成的危险区域范围继续沿东部电厂方向扩大;泄漏发生第150 s时，危险区域范围区域已蔓延至东部电厂区域;泄漏发生第200 s时，危险区域范围沿着东部电厂方向蔓延，并且有将整个电厂全部覆盖的趋势。

从图5可以看出，泄漏点在风速5 m/s条件下扩散初期呈现径向性扩散。当设置风向为东部电厂方向时，泄漏发生第50 s时，LNG泄漏产生的危险区域范围已经扩散到东部电厂敏感区边缘;在泄漏发生第100 s时，LNG泄漏扩散范围明显扩大，泄漏的LNG逐渐向空中扩散;泄漏发生第150 s时，LNG扩散至空中高度达到将近100 m，并且整个东部电厂敏感区已经开始被危险区域范围所覆盖;泄漏发生第200 s时，危险区域范围沿着东部电厂方向进一步扩大，并且有将整个场区完全覆盖的趋势。

通过将无风条件与有风条件进行对比，发现在无风条件下，LNG泄漏危险浓度区域横向扩散范围更大;而在有风情况下，泄漏的LNG形成的危险浓度区域向敏感区扩散速度更快，垂直向天空扩散的高度更高，但扩散区域内LNG浓度范围更小。因此，在距离泄漏点较近的敏感区内必须设置无风情况来考虑泄漏危险浓度对于敏感区的影响。造成这样不同的原因为，风速加快了大气的湍流，使得LNG泄漏产生的有害浓度被稀释消散。

2.2蒸气云爆炸模拟分析

为了研究LNG接收站持续泄漏时引发爆炸，不同风速对爆炸蒸气云浓度扩散的影响，对编号3的工况进行了模拟。

通过对比图6（a）和图6（b），发现无风条件下LNG接收站持续泄漏时引发爆炸蒸气云浓度形成规律为：以泄漏口为起点，向东部电厂方向扩散至罐区边缘处（浓度在5%～15%的甲烷浓度遇火极容易发生爆炸）;风速5 m/s条件下爆炸蒸气云浓度形成规律为：以泄漏口为起点，向东部电厂方向扩散至罐区边缘，并蔓延至东部电厂与储罐区间的小山坡地区。在此区域范围内遇明火则会引发爆炸事故，其爆炸范围最远波及至电厂与储罐区之间的小山坡地区。同时其爆炸蒸气云浓度分布相比无风条件，纵向（沿风速方向）蔓延的距离更远，而横向范围，无风条件下的爆炸蒸气云浓度的范围更广。

由图7可知，爆炸的燃烧区域并未覆盖整个爆炸极限浓度蒸汽云分布范围，这是因为爆炸火焰产生的高温热辐射使周围环境温度升高，间接加热了火焰周围的可燃气云，促使部分可燃气云的浓度迅速降至可燃下限以下，导致无法被燃烧。而由于爆炸高温气流所产生的压力差，促使大量LNG气云迅速爬升。此时近地面的不在可燃烧范围内的LNG气云也同时向燃烧爆炸的中心涌动，随着爆炸高温气流一起上升。在爆炸瞬间产生的超压冲击波的配合作用下，导致高温燃烧形成抽吸作用，造成附近的LNG气云大量减少，因此云爆最终发生在东部电厂附近，而并未发生在小山坡地区。

爆炸测点由图8绿色圆点位置M1处所示，M1的坐标为（0.125，-400.5，10.5）。由图9可知，爆炸初期，测点处超压冲击波的压强变化剧烈，随时间推移压强变化频率与冲击波压强的峰值逐渐减小。造成压强值反复变化的原因可归咎于，超压冲击波在传播过程中受到了周围障碍物的影响出现反弹现象，在反弹的过程中，受到空气阻力等影响，导致冲击波的能量逐渐衰减，因此反弹回来的爆炸冲击波压力峰值在逐渐减小。爆炸发生360 s后，超压冲击波压强值逐渐恢复到爆炸发生前状态。

由图10可知，爆炸开始时，超压冲击波的值短暂上升，此后迅速下降至约2.5×10-4 bar（1 bar=1标准大气压），之后在0.2 s内超压冲击波又迅速到达峰值，在0.3～0.8 s时超压冲击波压强值又开始迅速下降。造成超压冲击波压强值波动的原因可归咎于，前期爆炸蒸气云受到高温冲击后累计释放以及火焰运动趋势共同作用导致曲线上升;随后由于受到泄放慣性的影响，引发负压膨胀波的形成，从而导致超压冲击波压强值迅速下降;随着LNG气云与氧气混合经火焰射流引燃后会进一步燃烧，促使压力快速回升形成外部爆燃峰值，之后外场爆燃逐渐减缓压力值也随之下降。

2.3油池火灾模拟分析

为研究LNG接收站持续泄漏时引发油池火灾，其发生发展过程对厂区的影响，对编号4工况进行模拟，油池尺寸为289 m×50 m×0.005 m。

该场景油池火灾模拟地点为LNG接收站1号罐发生泄漏后在地面聚集形成液池，遇到火源后形成油池火灾，如图11、图12所示。

由图13、图14可知，火焰温度从液面到火焰底部存在一段迅速增加的过程，越过这一过渡区到达火焰底部后火焰的温度逐渐进入稳定阶段。火焰底部部分温度相对较低，是因为火焰底部到油面存在一个蒸气带，会吸收部分热量，且该处存在空气的吸卷作用，对流传热加强，使得温度略有降低。在整个模拟的时间里热释放速率最大值约为7.35×107 kW，空间内最大温度约为1 200 ℃。

由图15可知，距离火焰中心30 m的地方既1号储罐附近，在t=0到t=175 s左右辐射热流密度逐渐增大到45 kW/m²，在t=175 s左右时辐射热流密度达到最大90 kW/m²，这主要是由于油池燃烧过程中，在与燃烧壁面接触位置会产生气泡，随着持续燃烧，气泡产生量不断增加，当气泡量足够多时，就会跃离燃料与油池接触壁面，转移到油池液面，气泡扰动作用使得油池燃料流动性增强，形成了核态沸腾形象，因此辐射热流密度会迅速增加。最后在t=175到t=200 s的时间段内辐射热流密度维持在30 kW/m²上下波动。

由图16可知，距离火焰中心110 m的地方既靠近电厂附近，在t=0到t=150 s左右辐射热流密度呈上升趋势，在t=150到t=175 s左右时，辐射热流密度有下降的趋势，最后在t=175到t=200 s的时间段内辐射热流密度突然快速增大。通过对比图15与图16可得，火焰中心距离与辐射热流密度成反比，相邻油池表面各点的热辐射通量值与其距挡火墙距离成反比，距挡火墙越远热辐射通量值越小。

3对策措施

1）风速对LNG泄漏速度及浓度有较大影响，因此针对LNG接收站选址问题，在考虑经济效益的同时，应尽量选择通风良好的开阔地带。另外，可在泄漏点上风侧采取强制加风的措施，如防爆风机等，可加快大气湍流，使得有害浓度更快被稀释消散。

2）定期检查并维修罐体、阀门或焊缝等位置的功能完整性，加强日常巡检，避免因裂缝的出现导致LNG泄漏;同时降低地面粗糙度，避免地表液池的形成，可在罐区附近人为设置坡度或封闭的阻墙，在阻隔的同时引导气流方向。

3）重点关注泄漏点下风侧的危险区域，抢险人员进行救援时，尽量站于上风侧，并合理配戴防护用具。

4结论

1）基于CFD方法，采用FLACS软件对不同风速条件下LNG接收站泄漏进行了模拟，发现LNG持续泄漏扩散过程具有明显的分阶段特征。在无风条件下，LNG泄漏危险浓度区域横向扩散范围更大;当风速为5 m/s时，泄漏的LNG向敏感区扩散速度更快，垂直于天空扩散的高度更高，而扩散区域内LNG浓度较小。因此，在距离泄漏点较近的敏感区内必须考虑无风情况下泄漏危险浓度对于敏感区的影响。

2）通过研究不同风速对爆炸蒸气云浓度扩散的影响，发现5 m/s风速条件下，爆炸蒸气云浓度分布相比无风条件，纵向（沿风速方向）蔓延的距离更远，而横向范围，无风条件下的爆炸蒸气云浓度的范围更广;爆炸初期超压冲击波短暂上升后迅速下降，但在0.2 s内再次达到峰值。

3）采用FDS软件，对LNG接收站持续泄漏时引发油池火灾进行了模拟，掌握了其發生发展过程，发现火焰中心距离与辐射热流密度成反比;揭示了因LNG接收站泄漏引发连锁反应导致蒸气云爆炸与油池火灾的时空演化特征，为LNG接收站的选址、站内布局和储罐泄漏应急处置提供理论依据和支撑。

参考文献（References）：

[1]贾承造，张永峰，赵霞.中国天然气工业发展前景与挑战[J].天然气工业，2014，34（2）：8-18.JIA Chengzao，ZHANG Yongfeng，ZHAO Xia.Prospects of and challenges to natural gas industry development in China[J].Natural Gas Industry，2014，34（2）：8-18.

[2]邓军，梁天天，程方明.液化天然气厂区风险量化计算与分级管控[J].西安科技大学学报，2019，39（4）：571-580.DENG Jun，LIANG Tiantian，CHENG Fangming.Quantitative risk calculation and grading control of LNG plant area[J].Journal of Xian University of Science and Technology，2019，39（4）：571-580.

[3]王志赢，李成兵，周宁.大型LNG接收站泄漏事故灾害效应分析与预测[J].天然气工业，2019，39（5）：145-153.WANG Zhiying，LI Chengbin，ZHOU Ning.Analysis and prediction on the disaster effect of leakage accidents at large LNG receiving stations[J].Natural Gas Industry，2019，39（5）：145-153.

[4]凌晓东.槽车装卸区泄漏爆炸CFD模拟研究[J].消防科学与技术，2018，37（9）：1282-1286.LING Xiaodong.CFD simulation on leakage and explosion accident of truck loading station[J].Fire Science and Technology，2018，37（9）：1282-1286.

[5]曾岳梅，凌晓东.LNG接收站蒸气云爆炸数值模拟分析[J].消防科学与技术，2013，32（8）：834-837.ZENG Yuemei，LING Xiaodong.Numerical simulation of vapor cloud explosion in LNG receiving terminal[J].Fire Science and Technology，2013，32（8）：834-837.

[6]徐大用，蒋会春，姜威，等.基于FLACS的汽油槽车运输泄漏爆炸事故数值模拟研究[J].常州大学学报（自然科学版），2019，31（4）：16-25.XU Dayong，JIANG Huichun，JIANG Wei，et al.CFD simulation of gasoline tanker transportation explosion using FLACS[J].Journal of Changzhou University（Natural Science Edition），2019，31（4）：16-25.

[7]周宁，张倩，李雪，等.风速对LNG泄漏扩散过程影响的数值模拟[J].安全与环境学报，2021，21（1）：285-294.ZHOU Ning，ZHANG Qian，LI Xue，et al.Numerical simulation of the effect of wind speed on LNG leakage and the diffusion process[J].Journal of Safety and Environmen，2021，21（1）：285-294.

[8]杨磊，邹冀翼，寇丽颖，等.液化天然气储罐泄漏扩散模拟研究[J].石油化工自动化，2021，57（1）：48-50，62.YANG Lei，ZOU Jiyi，KOU Liying，et al.Simulation study of leakage diffusion of LNG tank[J].Automation in Petro-Chemical Industry，2021，57（1）：48-50，62.

[9]康泉胜，陆守香，陈兵.小尺度庚烷池火燃烧速率实验研究[J].科学通报，2010，55（1）：87-93.KANG Quansheng，LU Shouxiang，CHEN Bin.Experimental study on fire combustion rate of small-scale heptane ponds[J].Chinese Science Bulletin，2010，55（1）：87-93.

[10]黄坤，李芹，徐洁，等.风速对成品油库储罐池火灾热辐射分布的影响研究[J].热科学与技术，2019，18（5）：386-391.HUANG Kun，LI Qin，XU Jie，et al.Influence of wind velocity on heat radiation distribution in tank pool fire of refined oil depot[J].Journal of Thermal Science and Technology，2019，18（5）：386-391.

[11]杨帆，李康.两相流传热模型在LNG加热气化装置中的研究[J].农业装备与车辆工程，2021，59（9）：86-89.YANG Fan，LI Kang.Research on Two-phase circulating heat model in LNG heating gasification device[J].Agricultural Equipment & Vehicle Engineering，2021，59（9）：86-89.

[12]劉珊珊，焦文玲.LNG空温式气化器管内气化传热与流型数值模拟[J].油气储运，2017，36（8）：918-923.LIU Shanshan，JIAO Wenling.Numerical simulation on the vaporization，heat transfer and flow pattern inside the tube of LNG ambient air vaporizer[J].Oil & Gas Storage and Transportation，2017，36（8）：918-923.

[13]ZHANG Q X，LIANG D.Numerical simulations of LNG vapor dispersion from LNG jetting in different directions[J].Procedia Engineering，2016，13（5）：316-321.

[14]MARSEGAN C，BUSINI V，ROTA R.Influence of active mitigation barriers on LNG dispersion[J].Journal of Loss Prevention in the Process Industries，2016，4（4）：380-389.

[15]NGUYEN LD，KIM M，CHOI B.An experimental investigation of the evaporation of cryogenic liquid pool spreading on concrete ground[J].Applied Thermal Engineering，2017，12（3）：196-204.

[16]樊玉光，豆少刚，高琳.基于遗传算法的LNG气化混合工质配比研究[J].低温与超导，2022，50（1）：95-102.FAN Yuguang，DOU Shaogang，GAO Lin.Research on LNG gasification mixed working medium ratio based on genetic algorithsm[J].Cryogenics & Superconductivity，2022，50（1）：95-102.

[17]姜波，苗东升，刘岗.一起LNG储罐泄露原因分析[J].内蒙古石油化工，2018，44（4）：53-54.JIANG Bo，MIAO Dongsheng，LIU Gang.Analysis of leakage reasons from a LNG storage tank[J].Inner Mongolia Petrochemical Industry，2018，44（4）：53-54.

[18]李锐锋.国内LNG接收站设备设施完整性管理研究[J].化工管理，2021（34）：146-148.Ll Ruifeng.Research on management of the domestic LNG terminal equipment and facilities integrity[J].Chemical Management，2021（34）：146-148.

[19]张乙.基于FDS的LNG储备站火灾安全分析[D].北京：北方工业大学，2016.ZHANG Yi.Fire safety analysis of LNG reserve station based on FDS[D].Beijing：North China University of Technology，2016.

[20]高晋.大型LNG储罐大跨度穹顶结构稳定性分析[J].石油和化工设备，2021，24（11）：147-151.GAO Jin.Stability analysis of large span dome structure of large LNG storage tank[J].Petro & Chemical Equipment，2021，24（11）：147-151.

[21]鲁盈利，吕植勇，庄学强.LNG加注趸船蒸气云爆炸危害后果分析[J].武汉理工大学学报（交通科学与工程版），2015，39（4）：871-875.LU Yingli，LYU Zhiyong，ZHUANG Xueqiang.Harmful consequences analysis of vapor cloud explosion in LNG filling barges[J].Journal of Wuhan University of Technology（Transportation Science & Engineering），2015，39（4）：871-875.

[22]唐振东.LNG罐箱安全监督技术研究[D].大连：大连海事大学，2015.TANG Zhendong.Study on safe practice of portable LNG container operations[D].Dalian：Dalian Maritime University，2015.

[23]姚緒飞.海运LNG泄漏事故危害性研究[D].大连：大连海事大学，2015.YAO Xufei.Study on hazards caused by leakage from LNG tankers[D].Dalian：Dalian Maritime University，2015.

[24]杜彦.LNG加气站几个关键安全技术的研究[D].济南：山东建筑大学，2013.DU Yan.Research on several security technologies in LNG filling station[D].Jinan：Shandong Jianzhu University，2013.

[25]王志寰.大型LNG接收站场泄漏扩散与爆炸的灾害效应数值模拟研究[D].成都：西南石油大学，2019.WANG Zhihuan.Numerical simulation of the disaster effects of leakage dispersion and explosion at a large LNG receiving site[D].Chengdu：Southwest Petroleum University，2019.