障碍物投影长度对海上可燃气云爆燃特性影响*
高占胜1，刘文鹏1，2*，金良安1，蒋晓刚1
（1.海军大连舰艇学院航海系，辽宁 大连 116018；2.解放军91404 部队，河北 秦皇岛 066001）
摘 要： 针对海上运输和军事行动中可燃气泄漏形成气云等问题，利用涡耗散模型，对6 种障碍物投影长度（0 m～180 m）下海上可燃气云的形成和爆燃过程进行了数值模拟，重点研究障碍物表面及周围气云爆燃压力和温度特性的具体变化规律。结果表明，存在障碍物时气云湍流性增强，可燃气浓度分布更加不均匀；随着投影长度的增加，障碍物前方最大爆燃压力略增，但前表面无法形成持续高温；上方最大压力略降，投影长度小于100 m 时高温持续时间显著增加；后方最大压力骤降，投影长度大于100 m 时高温持续时间大幅缩短。结果可用于预测和评估不同投影长度下气云爆燃对船舶等障碍物的危害，为交通运输和军事行动中的安全防护等实际应用提供参考。
关键词： 海上可燃气云，爆燃特性，涡耗散模型，数值模拟，障碍物投影长度
0 引言
随着现代石油化工技术的快速发展和规模的不断扩大，海洋天然气开采、LNG 海上运输等也在不断增加。一旦发生意外导致可燃气泄漏，将在海上形成极易发生爆燃的可燃气云，可能导致海上平台和运输工具等遭到破坏［1-3］。同时，海上可燃气云爆燃也被广泛应用于军事攻击或防御等相关领域［4-6］，气云爆燃产生的高压和高温将对身处其中的舰船和人员造成严重伤害，因此，对海上可燃气云爆燃特性的研究具有重要的意义。
创新意味着突破传统，走不同寻常的道路，与创新相伴左右的便是风险，风险管理是考量一个新创企业面对市场环境的动荡所体现出来的能力，风险管理即营销职能寻求降低企业对外部环境的脆弱性和依赖性，通过协同营销方案提高企业的灵活性。“三只松鼠”已经走过了6个年头，取得了巨大的成就，成为中国互联网休闲食品的第一品牌，但在其成长过程中面临着众多风险与挑战，如食品质量问题、来自于同行的模仿、资金链运转等，其中章燎原认为“三只松鼠”最大的风险是食品安全。为了将风险控制到最低程度，“三只松鼠”利用互联网信息化、数据综合处理等手段，制定未来发展的四大战略：大健康、大娱乐、大品类、大消费，这构成了松鼠未来新消费主义。
目前大多数对可燃气云爆燃的研究文献中［7-11］，其实验或模拟的环境主要为陆地，在爆燃特性研究方面多以压力为主，气云形状也都假设为理想的对称形状，并认为气云在爆燃时已经混合均匀。而海上属于开敞空间，相对于产生的高压，可燃气云爆燃造成的高温伤害更加不容忽视，且海上可燃气云的形状往往并不规则［12］，可燃气的不断泄漏也会导致气云内部存在浓度差异［13-14］。重要的是，可燃气扩散主要受风速和风向影响［15］，因此，气云中的障碍物在下风向的投影长度会对爆燃过程产生很大影响，而文献中对此却鲜有研究。
鉴于以上现状，为突出湍流混合对燃烧速率的控制作用，本文基于k-ε 模型并结合涡耗散模型，以船舶作为障碍物原型，将其简化为长方体，针对障碍物在气云下风向不同投影长度下的爆燃过程进行数值模拟，重点研究投影长度对障碍物附近气云爆燃压力和温度特性的具体影响，以探讨气云爆燃可能给船舶等障碍物带来的危害。
1 数值模拟
1.1 物理模型
物理模型如图1 所示，x 轴正向为风向，y 轴代表高度，计算区域长400 m，高80 m；可燃气泄漏源位于海面，入口宽10 m，中心坐标（40，0）；点火区域位于泄漏源处，点火源中心坐标（40，0.5）；根据一般船舶尺寸，假设障碍物长度为180 m，宽度为20 m，高度为16 m，根据船舶首尾连线与下风向所成夹角，依次设定其在下风向的投影长度为20（首尾连线垂直风向，投影为船宽）、60、100、140、180 m（首尾连线平行风向，投影为船长）5 个等差值，并添加无障碍物组（0 m）作为对比。为了全面研究障碍物前方（靠近泄漏源，1、2 号监测点）、上方（3、4 号监测点）、后方（远离泄漏源，5、6 号监测点）等位置的压力和温度变化，分别设置6 个监测点，其中1、3、5号监测点位于障碍物各表面中间，同时在水平（垂直）距离各表面点20 m 处设置2、4、6 号监测点作为对比，无障碍物组的监测点位置与投影长度20 m 组相同。
整个流场采用四边形网格进行划分，由于可燃气入口位置初始浓度较大且为中心点火区域，障碍物附近流场参数变化较为剧烈，因此，对这两个区域的网格进行了局部加密，加密网格大小为0.4 m。
根据人体临界伤害温度391 K、人体严重伤害温度453 K、钢结构部分失效温度673 K 和钢结构全部失效温度873 K 这4 个温度界限［18］，表2 分别对6 个监测点超出最高临界温度873 K 的时间进行了统计，同时建立图6 所示6 种投影长度下气云爆燃300 s 时的临界温度等值面。
1.2 数值计算方法
可燃气扩散流动的连续方程、动量方程、能量方程和组分方程可以统一表示如下［1］：
“李老师，不好了，小刘用厕所门把小邹的手指夹破了。”班长小冯急匆匆地来向我报告。我急忙跟着小冯来到了教室，只见小邹捂着手指痛哭不止。我一面安抚小邹，带她去医务室简单处理，一面打电话给她妈妈告知情况。而一旁的小刘似乎也很委屈，一直强调自己是无意的，同行的几名学生也连连点头。赶来的小邹妈妈看到小邹手指并无大碍，悬着的心也放下了。这件事就过去了。
式（1）中，ρ 为密度，φ 为通用求解变量，t 为时间，ui为i 方向上的速度分量，Γφ 为通量φ 的扩散系数，Sφ为源项。
本文基于标准k-ε 模型，并利用涡耗散模型模拟可燃气云的爆燃过程：
由图6，从气云爆燃的空间分布来看，障碍物前方始终处于453 K～673 K 的人体严重伤害区域；上方则随着投影长度的增加，爆燃温度逐渐升高，进入873 K～2 000 K 的钢结构全部失效区域，甚至部分区域的温度已经超过2 000 K；后方则随着投影长度的增加，爆燃温度逐渐降低，由钢结构全部失效区域进入低于391 K 的安全区域。
图1 物理模型
乙炔-空气混合气的爆炸上、下限分别为80%和2.5%（体积分数），因此，以2.5%等值面（转换为质量分数为2.27%）作为扩散危险区域。图2 为扩散300 s 时6 种投影长度下可燃气在海面上方的危险区域分布，相对于无障碍物组，障碍物组前方的可燃气浓度较高，各组障碍物均处于气云包围中。随着障碍物投影长度的增加，气云波动性增加，湍流性增强，这是因为障碍物阻碍了可燃气向外扩散，当可燃气翻越障碍物与空气形成对流时，会在障碍物上方形成小涡团，引起气云的不稳定；同时随着投影长度的增加，障碍物后方的可燃气浓度逐渐降低，当投影长度为140 m～180 m 时，后方可燃气的质量分数已经低于0.225。
图5 为不同测点处最大爆燃压力随投影长度的变化关系，随着投影长度的增加，障碍物前方的最大爆燃压力略有增加，上方则略有减小，而后方的最大爆燃压力显著减小。因为障碍物前方的压力冲击波遭到阻挡，导致爆燃产生的瞬间压力更大；而随着投影长度的增加，可燃气绕过障碍物扩散的难度加大，不易扩散至后方海面，导致该处可燃气较为稀薄，无法产生爆燃高压。
1.3 初始条件及边界条件
计算环境初始压力为1.01×105 Pa，初始温度为300 K，氧气质量分数为0.23，二氧化碳及水蒸气的质量分数均为0。为了尽可能突出气云爆燃的危险性，可燃气选择爆炸极限范围大、爆炸能较高的乙炔气体，泄漏源喷出质量流40 kg/s，垂直从海面进入大气，并考虑到可燃气入口处也存在部分氧气，设入口处乙炔和空气的质量比为9：1。气云形成后，采用瞬间高温进行点火，点火温度2 000 K。
计算域顶部流体变量梯度为零，来风口采用速度入口边界条件，出风口选择自由出流边界条件，风向平行于海面，不随时间、位置变化，y 方向的风速则采用幂次律计算［16］：
表1 为不同投影长度下6 个监测点的最高爆燃温度，相对于无障碍物组，障碍物组前表面测点1处温度普遍较低，前方测点2 处最高温度相差不大；随着投影长度的增加，上方测点3、4 处最高温度略有增加，但当投影长度为100 m～180 m 时温差并不大，后方测点5、6 处的趋势与上方类似。
2 结果分析
2.1 投影长度对可燃气云形成的影响
式（2）～式（4）中，μeff 为有效粘性系数，G 为阻力影响因素，δk和σε为输运特性的湍流普朗特数，C1和C2为常数；Ri，r为燃烧反应速率，v'i，r和v''i，r 分别为反应物和生成物的化学计量数，M 为分子量，m 为质量分数，A 和B 为常数。
图2 可燃气扩散300 s 时形成的危险区域
2.2 投影长度对爆燃压力特性的影响
图3 投影长度为100 m 的爆燃模拟压力
图4 文献［11］的爆燃实验压力
图5 不同测点处的最大爆燃压力
对扩散300 s 时的气云点火开始爆燃数值模拟，爆燃持续时间300 s。图3 是投影长度为100 m时不同测点处气云爆燃压力随时间的变化关系，在靠近泄漏源的测点2 处有最大正超压和最小负压，而在远离泄漏源的测点6 处正负超压则较小，这是因为爆燃中心位于泄漏源处，爆燃产生的压力冲击波从中心逐渐向周围传播，随着传播距离的增加，波阵面面积不断增大，加之过程中有能量损耗，单位体积的冲击波能量越来越小。图4 为气云爆燃实验压力随时间的变化关系［11］，爆燃过程中压力变化的总趋势为：在气云爆燃初期，超压迅速增加，达到最大正超压后又急剧减小到负压，并且随着时间的推移，达到二次正超压。可见本文的数值模拟结果与文献［11］的实验压力变化过程基本吻合，说明该数值模型能够较好地预测内置障碍物的气云爆燃过程。
利用FLUENT 进行基本控制方程的离散和求解。其中，方程的离散采用有限体积法，压力-速度耦合方程求解采用SIMPLE 算法，梯度差分采用基于单元体的Least-Squares 方法，对源项进行局部线性化处理以使方程封闭，时间采用二阶隐式积分，压力采用二阶格式，其余皆采用二阶迎风格式。
根据压力伤害准则［17］，人体轻微伤害临界值为20 kPa、人体致命伤害临界值为100 kPa、钢结构轻微毁坏临界值为3 kPa、钢结构严重毁坏临界值为30 kPa，可见海上可燃气云爆燃产生的压力较小，高压作用时间不长，至多使障碍物产生轻微损坏。这是因为海面为完全开敞空间，可燃气云爆燃后产生的压力很容易泄放，且海面可燃气扩散区域面积较大，障碍物的存在又使得可燃气浓度分布十分不均匀，很难恰好处于燃烧反应当量浓度，自然无法产生很高的超压。
我国对问题银行进行早期干预、接管、重组、破产、撤销等相关法律规定散见于《存款保险条例》《银行业监督管理法》《企业破产法》等法律法规。相关具体条款主要是明确各监管机构相应的风险处置职责，且非常原则性，没有清晰地界定出问题银行是否需要救助、是否需要关闭，怎样重组、怎样清算退出等一整套处置规则，且有的法规之间缺少有效衔接。
2.3 投影长度对爆燃温度特性的影响
式中，uz 为距海面高度z 处的风速，m/s；u10 为海面上方10 m 处的风速，m/s；海上可取幂指数n =0.107。式（5）由UDF 文件建立。
表1 不同投影长度下各测点的最高温度
在他的带领下，全局工作热情高涨，纪律作风严明，涌现出了一批先进典型，3名同志分别被评为市“十大杰出青年”“十大道德模范”和省局优秀党员，1名同志被评为济南市基层监管新闻人物，成为干部群众学习的楷模。
秀容月明对皇上说，梨友共带了八万人马，你只让我带两千秀容兵，再把驻守桂州的军士算上，也不过一万七千人，太少了。
表2 各测点超出临界温度873 K 的累计时间
根据表2 的统计数据，相对于无障碍物组，障碍物组测点1 处的爆燃温度低于873 K，测点2 处的爆燃高温持续时间更长；随着投影长度的增加，当投影长度为0 m～60 m 时测点3、4 处的高温持续时间迅速增加，100 m～180 m 时高温持续时间的变化不大，表面测点3 比上方测点4 的高温持续时间更长；而在测点5、6 处，0 m～60 m 时高温持续时间的变化不大，100 m～180 m 时随投影长度的增加高温持续时间迅速减少。
图6 可燃气云爆燃300 s 时临界损伤温度区域分布
此前毫无战争经验的卫青被刘彻任命为车骑将军，率一万骑兵从上谷出塞。卫青率部突袭匈奴的祭天圣地龙城，斩杀了七百余匈奴人，战胜而归。卫青第一仗就取得胜利，被封为关内侯。
综合表1、表2 和图5、图6 可见，障碍物的存在提高了爆燃火焰传播的湍流程度，湍流又加速了火焰的传播，从而进一步增强了燃烧的速率和强度，产生了更大的破坏力。由于障碍物前表面处可燃气扩散受到阻碍，导致该处可燃气浓度过高，可供燃烧的氧气较少，因而燃烧并不剧烈，但该处仍会受到前方压力冲击波的伤害；同时，随着投影长度的增加，障碍物上方依附的可燃气体量逐渐增多，导致爆燃温度升高。而在障碍物后方，由于障碍物对扩散和火焰的阻碍作用，随着投影长度的增加，爆燃火焰基本无法蔓延至该区域。
由于船舶等障碍物的两侧都是较为坚固的钢结构，能够耐受一定程度的冲击和高温伤害，但其甲板以上部分的防护往往较为薄弱，因此，当投影长度较大时，障碍物上方将持续处于高温中，成为遭到爆燃火焰波及的重灾区，可能导致处于该区域的装备设施完全失效，造成作业人员严重伤亡。
最后，还需要引导学生将全域旅游理念运用到导游实践工作中。在当今时代，随着全球可持续发展的发展和旅游业的迅速发展，全域旅游的概念引起了人们的关注。全域旅游理念的本质是协调和平衡旅游与文化、自然、人的生活环境之间的关系，并将其作为一个整体，在全球范围内实现社会发展与经济发展目标的统一。因此，在全域旅游发展过程中，需要引导学生将全域旅游理念运用到导游实践工作中，开发城市绿色旅游，促进旅游绿色发展。
3 结论
利用k-ε 湍流模型和涡耗散模型，对海上可燃气云的形成和爆燃过程进行了数值模拟，通过对比数值模拟压力与文献实验压力的变化过程，验证了模拟的可靠性，由数值模拟结果揭示了投影长度对气云爆燃压力和温度特性的具体影响，得到以下结论：
1）存在障碍物时气云的湍流性增强，可燃气的浓度分布更加不均匀，在部分区域障碍物加强了爆燃造成的伤害。
2）随着投影长度的增加，障碍物前方爆燃最大压力略增，但在前表面无法形成持续高温；障碍物上方爆燃最大压力略降，投影长度小于100 m 时爆燃高温的持续时间显著增加；障碍物后方爆燃最大压力骤降，投影长度大于100 m 时爆燃高温的持续时间大幅缩短。
3）较大的投影长度将使船舶防护较弱的甲板以上区域处于持续高温中，从而加重了气云爆燃对该区域内的装备设施和作业人员的伤害；当投影长度大于180 m 时，障碍物后方将成为安全区。
参考文献：
［1］杨冬冬，陈国明，师吉浩.海洋平台井喷含硫天然气扩散危险区域研究［J］. 中国安全生产科学技术，2017，13（8）：114-120.
［2］聂细亮，刘强，李国帅，等.LNG 船舶海上运输安全研究现状及进展［J］.中国安全科学学报，2017，27（5）：87-92.
［3］LI X H，CHEN G M，ZHU H W. Quantitative risk analysis on leakage failure of submarine oil and gas pipelines using Bayesian network ［J］. Process Safety and Environmental Protection，2016（103）：163-173.
［4］周楠，詹秀秀，朱燕燕，等.舰艇编队防空能力仿真和可疑目标识别模型［J］. 火力与指挥控制，2017，42（7）：39-43.
［5］刘文鹏，金良安，迟卫，等.海上水基燃烧技术的高可靠性自致燃方法［J］.科学技术与工程，2017，17（24）：132-137.
［6］刘文鹏，金良安，高占胜，等.海上水基燃烧技术多燃烧模式的实现研究［J］.火工品，2017，39（2）：19-23.
［7］毕明树，王淑兰，丁信伟.无约束气云爆燃压力场的计算［J］.爆炸与冲击，2002，22（2）：148-151.
［8］毕明树，丁信伟，王淑兰，等.半球形气云爆炸的数值模拟［J］.天然气工业，2003，23（4）：86-89.
［9］杨国刚，岳丹婷，毕明树.圆柱形可燃气云爆炸实验研究与数值模拟［J］.化工学报，2008，59（11）：2954-2959.
［10］陈文瑛，柴建设.有障碍物开敞空间可燃气云爆炸超压场的数值模拟［J］.中国安全科学学报，2009，19（6）：35-40.
［11］党福辉.开敞空间可燃气云爆炸数值模拟研究［D］.天津：天津理工大学，2017.
［12］艾志久，吴光武，王其华，等.井喷失控气云爆炸的3维数值模拟［J］.爆炸与冲击，2010，30（1）：96-100.
［13］刘文鹏，金良安，蒋晓刚，等.海面稀氧区的人工形成及其战术效能研究［J］.安防技术，2017，5（4）：51-58..
［14］金良安，刘文鹏，高占胜，等.大气稳定度对海面人工稀氧区形成的影响［J］.交通运输工程学报，2016，16（6）：99-106.
［15］刘文鹏，金良安，高占胜.风速对海面稀氧区人工形成的影响［J］.海军大连舰艇学院学报，2017，40（3）：42-46.
［16］HOLTSLAG M C，BIERBOOMS W A A M，BUSSEL G J W V.Extending the diabatic surface layer wind shear profile for offshore wind energy［J］.Renewable Energy，2017（101）：96-110.
［17］王建.储罐区可燃气体泄漏扩散模拟及爆燃灾害评估［D］.大连：大连理工大学，2013.
［18］魏超南，陈国明，刘康.浮式生产系统泄漏天然气爆燃特性与安全区域［J］.石油学报，2014，35（4）：786-794.
Influence of Obstacle Projection Length on Offshore Flammable Gas Cloud Deflagration Characteristics
GAO Zhan-sheng1，LIU Wen-peng1，2*，JIN Liang-an1，JIANG Xiao-gang1
（1.Department of Navigation，Dalian Naval Academy，Dalian 116018，China；2.Unit 91404 of PLA，Qinhuangdao 066001，China）
Abstract： For the problem of flammable gas leakage in the process of marine transportation and military action，such as the formation of flammable gas cloud，by eddy-dissipation model，the formation and deflagration process of offshore flammable gas cloud are numerically simulated with six kinds of obstacle projection length （0～180 m），the variation of gas cloud deflagration pressure and temperature characteristic in obstacle’s surface and surrounding is studied specifically. The results show that when the obstacle exists，the gas cloud is more turbulent and the flammable gas concentration distribution is more uneven. With the increase of projection length，the maximum deflagration pressure in front of the obstacle is slightly increased，but the front surface cannot form.The maximum pressure above is slightly lower，and the continuous time of high temperature significantly increases when projection length is less than 100 m. The maximum pressure rear is sudden drawdown，and the continuous time of high temperature significantly shortens when projection length is greater than 100 m. The results can be used to predict and evaluate the hazards of gas cloud deflagration on the ships or other obstacles under different projection lengths，and provide reference for practical applications such as safety protection in transportation and military operations.
Key words： offshore flammable gas cloud，deflagration characteristics，eddy-dissipation model，numerical simulation，obstacle projection length
中图分类号： TJ99；O383
文献标识码： A
DOI： 10.3969/j.issn.1002-0640.2019.11.018
引用格式： 高占胜，刘文鹏，金良安，等.障碍物投影长度对海上可燃气云爆燃特性影响［J］.火力与指挥控制，2019，44（11）：88-92.
文章编号： 1002-0640（2019）11-0088-05
收稿日期： 2018-06-08
修回日期： 2018-09-12
*基金项目： “ 十三五”国防预研基金资助项目（30203010303）
作者简介： 高占胜（1965- ），男，河北昌黎人，副教授。研究方向：船舶操纵的理论和方法、军事海洋环境。
*通信作者： 刘文鹏（1993- ），男，辽宁本溪人，博士研究生。研究方向：军事航海安全保障与防护技术。
Citation format： GAO Z S，LIU W P，JIN L A，et al.Influence of obstacle projection length on offshore flammable gas cloud deflagration characteristics［J］.Fire Control & Command Control，2019，44（11）：88-92.
标签：海上可燃气云论文; 爆燃特性论文; 涡耗散模型论文; 数值模拟论文; 障碍物投影长度论文; 海军大连舰艇学院航海系论文; 解放军91404部队论文;