泄漏源模型

泄漏源模型

泄漏源模型根据泄漏源位置、形式与特征的不同，可将其分为密封元件的渗漏模型、储罐或管道的泄漏模型和泄压元件的泄放模型3种类型。 1 渗漏模型

化工系统所发生的重大泄漏事故大部分是由于密封失效、密封件设计或安装不合理造成的。流体的密封通常是靠密封面间的相互紧密接触以增加流动阻力来实现的，但由于不可能实现密封面间的完全吻合和密封件毛细孔的完全阻塞，流体就可能通过密封件与被密封件间的间隙或通过密封件本身内部的孔隙渗漏，根据流体的渗漏通道不同，可将渗漏模型分为平行圆板模型、三角沟槽模型和多孔介质模型3种。

1.1 平行圆板模型

平行圆板模型将流体介质通过密封点处的泄漏简化为介质通过间隙高度为h ，由内径r 1处流至外径r 2处的定长、层流流动，其体积泄漏率为： ()3

21216p p h L r In r υπη-=⎛⎫ ⎪⎝⎭ （1）

式中：η为介质粘度，p 2、p 1分别为垫片内、外侧的压力。

1.2 三角沟槽模型

三角沟槽模型认为，在正常的密封情况下，垫片与法兰面的间隙由许多三角沟槽所组成，设H 为三角沟槽的深度，L 为三角沟槽的底宽，b 为流道的长度( 通常为垫片的宽度)，ρ为介质密度，则体积泄漏率为：

对于液体：

3v LH p L C b

η∆= （2） 对于气体：

()2

312v LH p L C p b

η∆= （3） 式中：()2222121,p p p p p p ∆=-∆=-，C 为常数。 1.3 多孔介质模型

多孔介质模型认为非金属垫片可近似看作各向同性的多孔介质，其流道由多个弯弯曲曲、半径大小不等的毛细管组成。气体通过多孔介质可分为层流流动和分子流流动，其气体的总流率为层流流率与分子流流率之和。研究表明毛细管半径r 随垫片残余应力σ的增大而减小，存在

()n r f σ-=的关系。这样就可以得到

气体通过垫片的泄漏率方程： ()()()()()122121nL nM pv L m M L A p p p A b T M p p ησσ--=-+- （4）

式中，L A 、M A 、nL 、nM 为常数，其值可由实验得到，pv L 为PV 泄漏率，()21/2m p p p =+，M 是气体相对分子质量，T 为气体绝对温度。

2 泄漏模型

物质存储形式的多种多样、漏源的大小、形状、位置以及泄漏介质本身物理特性的不同决定了泄漏形式的多样性和复杂性。影响泄漏扩散的因素主要有介质的相态(气态或液态)、储存条件(压力液化储存、冷冻液化储存、常态液体储存和常态气体储存)、弥散限制(泄漏源周围有无防液堤)和泄放形式(连续泄漏、瞬时泄漏和有限时间泄漏)。

2.1 储罐泄漏

罐壁上的腐蚀、疲劳裂纹或孔洞以及碰撞、容器超压都能导致储罐泄漏。

2.1.1 气体泄漏

对于高压(低温)液化气储罐，如果处于满装状态，罐内不存在气相空间，此时即使少许裂缝出现，由于少量液体的泄漏也会引起内压的迅速下降而处于过热状态，液体全部汽化，从而最终导致灾难性破裂(闪蒸)；如果储罐没有满装，当破裂处位于气相空间时，在破裂面积较大的情况下， 高压蒸气通过裂缝或孔洞喷出，储罐内压急剧下降，直到环境压力(常温)。由于内压急剧下降，气液平衡遭到破坏，储罐内流体处于过热状态，过热状态的液体为了再次恢复平衡，内部会均匀地产生沸腾核，同时产生大量气泡，液体体积急剧膨胀，最终也导致蒸气爆炸。对于以上两种情况(闪蒸)，泄漏量可按存储介质瞬间全部泄漏计算。若裂口面积不大，即使有蒸气喷出，但由于储罐内压下降不急剧，液体不会达到过热状态，因此不会发生蒸气爆炸。气体或蒸气不同于液体，它属于可压缩流体。当气体或蒸气以较高速度流动时（>0.3倍音速）或前后压差大于2时，其在流动过程中的动能变化和物理性质的变化（尤其是密度的变化）就必须加以考虑。气体或蒸气的泄放可分为节流泄放（Throttling Release ）和自由泄放（Free Expansion Release ）。对于节流泄放，气体或蒸汽的压缩能绝大部分用来克服摩擦阻力；而对于自由泄放，则气体或蒸汽的压缩能绝大部分转化为动能。节流泄放模型需要裂口的详细物理特征，在这里就不作讨论了，而只考虑较为简单的自由泄放模型。据机械能守恒原理，得到气体或蒸气通过孔洞泄漏的质量流速模型：

{}1/22/(1)/000002/(1)(/)(/)r r r m g Q C AP rM R T r P P P P +⎡⎤⎡⎤=-⨯-⎣⎦⎣⎦ （5）

式中：

m

Q ——质量泄漏率/(kg/s)； 0C ——泄漏系数；

A ——裂口面积/m 2；

0P ——储罐内压/Pa ；

M ——气体或蒸气的摩尔质量（kg/mol ）；

g R ——理想气体常数；

0T ——泄漏源温度/K ；

P ——泄漏处压力/Pa ；

r ——热容比。

泄漏过程在临界状态时，泄漏源流量最大，此时泄漏处于塞压状态。对于理想气体而言，塞压是热容比的函数，见表1。

表1 塞亚和热容比的关系

临界状态下，最大质量流量的计算公式如下： ()()[]{}1/2(1)/(1)000/2/(1)r r m g choked Q C AP rM R T r +-=+ （6）

泄漏系数的确定直接影响气体泄漏速度的计算。一般而言，泄漏系数的取值范围在0.6~1.0之间。按泄漏孔的形状可分为：圆形孔，

0 1.0C =；三角孔，00.95C =；长形孔，00.90C =（棱越多，泄漏系数越小）。孔口为内层腐蚀形

成的渐缩孔，00.9 1.0C <<；孔口为外力机械损伤形成的渐扩孔，00.60.9C <<。

对于气体或蒸气泄漏，必须考虑动力抬升和热力抬升。其中动力抬升是由泄漏方向决定，热力抬升是由介质与环境的温差决定的。目前普遍使用国标GB3840-83推荐的抬升公式，但该公式只适用于泄漏方向竖直向上和泄漏介质的温度大于环境温度。不同的泄漏方向和泄漏介质与环境的温差将产生不同的泄漏扩散效果。

对于常压气体储存，由于交通事故或系统超压导致储罐大面积开裂或超压爆炸，可认为气体瞬间全部泄放到大气中，形成云团。

2.1.2 液体泄漏

对于高压(低温)液化储罐，当裂口处位于液相空间时，尽管液体流出并可能发生闪蒸，但由于液体的流出阻力大，内压下降速度缓慢，储罐内过热液体不会发生蒸气爆炸。闪蒸所需能量来自于过热液体中所储存的能量，即：