LNG储存容器非稳态泄漏模型研究.docVIP

LNG储存容器非稳态泄漏模型研究

1

概述

储存LNG的容器由于材料缺陷、操作不当或附属设备仪表失效等原因往往会在使用过程中发生不同程度的泄漏

［1］。LNG泄漏及蒸发产生的可燃气体扩散引发的火灾、爆炸事故时有发生，对周围生命、财产安全和环境都会造成很大程度的威胁

［2-4］。

泄漏率为单位时间内从失效孔口处泄漏的液体或两相流质量。在LNG储存容器定量风险评价中，LNG泄漏率以及某个时段内的泄漏量是预测、分析和评价蒸发扩散后果的重要参数，且容器的实际泄漏过程是一个动态变化的过程，即泄漏率和容器内压力、液面高度等参数是非恒定的，因此，有必要对LNG储存容器非稳态泄漏过程的数学模型进行研究。

由伯努利方程得到的孔泄漏模型仅适用于稳态泄漏的情况，若用于非稳态泄漏计算，得到的是泄漏初始时刻的泄漏率，即最大泄漏率，计算结果会高于实际值。为此，有学者考虑液面下降以及由液面降低引起的气相空间压力降低，建立了新的泄漏模型

［5］，但该模型中没有考虑过热LNG的蒸发现象。

本文将在现有模型基础上，将LNG蒸发过程考虑在内，结合实际泄漏过程中容器内各参数的动态变化规律，建立能够准确表征容器非稳态泄漏过程的模型，结合实例进行分析，并通过分析简化，得到泄漏率的两种简化计算方法。

2

LNG储存容器非稳态泄漏模型

图1

LNG储存容器失效泄漏过程

LNG有较大的比热容，即使在泄漏过程中损失部分热量用于闪蒸，其温度变化也非常小，可忽略不计。蒸发气不断进入气相空间，与原有气体进行混合，也近似认为气相空间的温度不发生变化，均与未发生泄漏时的储存温度一致。因此泄漏率计算式中仅有气相空间压力和液面高度是随时间变化的，其他均为定值。泄漏过程中随着液面的下降，容器内气相空间体积会增加并导致压力降低，但此时液体处于过热状态，而蒸气处于未饱和状态，会有部分液体蒸发扩散至气相空间中，对降低的压力有补充作用。气相空间状态可由理想气体状态方程表示：

图2

泄漏模型计算流程

3

应用实例与分析

某LNG储存容器，其容积为100m3，内部尺寸为14.3m×7m2，充装高度为13m，储存压力（绝对压力）为0.6MPa，容器内处于饱和状态，假定LNG的成分为纯甲烷，可知容器的储存温度为-134

℃，液体的蒸发热为447kJ/kg。因底部连接管道与容器的接口处失效致使大量LNG泄漏，失效裂口近似圆形，当量面积约0.002m2。另有一容器，储存压力为常压101325Pa，储存温度为甲烷常压下的沸点-162

℃，其他参数及失效形式与前述容器均相同。

由于联立的方程数较多且方程形式较为复杂，因此无法得到解析解，但可通过软件计算得到数值解，绘制散点图并连接或拟合成曲线。此处采用MATLAB的ode功能函数对微分方程组进行计算，数值解的计算步长由参数的变化速度决定，最小可至0.1s，可准确地表征实际变化过程。

①容器内气相压力

根据文中提供的非稳态泄漏模型，得到泄漏过程中两个容器内的气相压力随时间变化的散点，其振荡速度非常快，但振荡幅度与储存压力相比却很小，振荡幅度不超过储存压力的0.5%，因此采用线性拟合来表征其变化规律，分别见图3和图4。泄漏开始后极短时间内，压力迅速降低，随后发生快速振荡，而振荡变化总体又呈逐渐上升趋势。

图3

压力容器内气相压力随时间的变化曲线

图4

常压容器内气相压力随时间的变化曲线

②容器内液面高度

泄漏过程中两个容器内液面高度随时间的变化曲线见图5和图6。压力容器泄漏率仅随气相压力变化，由上文可知，气相压力变化幅度很小，因此液面高度随时间可以看作是线性变化的。常压容器中气相压力与外界大气压近似相等，液面产生的位置压头是促使LNG泄漏的动力来源，随着液面的下降，液面高度降低的速度也逐渐减小。当液面高度低于0.5m时，由于气相压力的不断振荡，有时会出现总作用压头为负值的情况，其对应的实际过程为液面短时停留在某一位置，随着蒸发气的产生，气相压力回升，作用压头为正，容器内液体继续向外泄漏，但同时又伴随着液面下降，气相压力降低，总压头转为负值，循环重复上述过程，直至容器内液体泄漏完全。但该过程会持续较长时间并且平均泄漏率与整个泄漏过程的泄漏率相比很小，因此对液面降至0.5m以下的泄漏过程不予考虑，认为此时已泄漏完全。

图5

压力容器内液面高度随时间的变化曲线

图6

常压容器内液面高度随时间的变化曲线

③泄漏率

两个容器任意时刻的泄漏率见图7和图8。压力容器的泄漏率近似为定值，而常压容器的泄漏率随时间逐渐降低，近似呈线性变化。

图7

压力容器内泄漏率随时间的变化曲线

图8

常压容器内泄漏率随时间的变化曲线

4

结论

①LNG储存容器液体泄漏过程是一个非稳态过程，应同时考