事故泄漏模型分析及其在环境风险评价中的应用

事故泄漏模型分析及其在环境风险评价中的应用

——化工石化医药类环境影响评价登记培训论文

摘要：化学危险品事故泄漏风险评价是整个环评过程中的一个有机组成部分，是进行风险评价预测和模拟的前提。本文以《建设项目环境风险评价技术导则》(HJ/T169－2004)为基础，分析了可能出现事故泄漏的设备和事故泄漏的后果。根据不同事故情况，对导则中事故泄漏模型及预测模型的应用进行了分析。

关键词：风险评价；事故泄漏；源模型；临界流；两相泄漏；气体泄漏；液体泄漏；多烟团模型

1 前言

在化工、石油化工及相关行业中，易燃、易爆及有毒有害物质在生产、储存和运输过程中经常发生泄漏事故。事故的发生不仅会导致巨大的经济损失，而且还会造成严重的人员及环境生态的毒性伤害和污染。更为严重的是可能会继而发生火灾或爆炸等灾害，使得灾害损失与破坏进一步加剧。《建设项目环境风险评价技术导则》(HJ/T169-2004)要求对建设项目建设和运行期间发生的可预测突发性事件或事故引起有毒有害、易燃易爆等物质泄漏，或突发事故产生的新的有毒有害物质，所造成的对人身安全与环境的影响和损害，进行评估，提出防范、应急与减缓措施［1］。对事故泄漏源进行分析，主要是根据项目所涉及的危险物品的化学性质、事故下设备情况，采取相应的数学模型来估算泄漏物的排放量、排放时间等。在计算得到事故泄漏源强参数后即可采用扩散模型进一步对事故泄漏对环境的影响进行预测分析。

2 泄漏情况

根据各种设备泄漏情况分析，可将工厂(特别是化工厂)中易发生泄漏的设备分类，通常归纳为：管道、挠性连接器、过滤器、阀门、压力容器或反应器、泵、压缩机、储罐、加压或冷冻气体容器及火炬燃烧装置或放散管等10类［2］。每一种设备的典型损坏类型及其典型的损坏尺寸不同，一般可按设备大小的20%～100%计算。泄漏一旦出现，其后果不单与物质的数量、易燃性、毒性有关，而且与泄漏物的相态、压力、温度等状态有关。［3］这些状态可有多种不同的结合，在后果分析中，常见的可能结合有常压液体、加压液化气体、低温液化气体、加压气体4种。无论是气体泄漏还是液体泄漏，泄漏量的多少都是决定泄漏后果严重程度的主要因素，而泄漏量又与泄漏时间长短有

关。

3 泄漏模型

物质存储形式的多种多样、漏源的大小、形状、位置以及泄漏介质本身物理特性的不同决定了泄漏形式的多样性和复杂性[4]

。影响泄漏扩散的因素主要有介质的相态(气态或液态)、储存条件(压力液化储存、冷冻液化储存、常态液体储存和常态气体储存)、弥散限制(泄漏源周围有无防液堤)和泄放形式(连续泄漏、瞬时泄漏和有限时间泄漏)。

当发生泄漏的设备的裂口是规则的，而且裂口尺寸及泄漏物质的有关热力学、物理化学性质及参数已知时，可根据流体力学中的有关方程式计算泄漏量。当裂口不规则时，可采取等效尺寸代替；当遇到泄漏过程中压力变化等情况时，往往采用经验公式计算。

3.1 气体泄漏

以储罐为例，对于高压(低温)液化气储罐，如果处于满装状态，罐内不存在气相空间，此时即使少许裂缝出现，由于少量液体的泄漏也会引起内压的迅速下降而处于过热状态，液体全部汽化，从而最终导致灾难性破裂；如果储罐没有满装，当破裂处位于气相空间时，在破裂面积较大的情况下，高压蒸气通过裂缝或孔洞喷出，储罐内压急剧下降，直到环境压力(常温)。由于内压急剧下降，气液平衡遭到破坏，储罐内流体处于过热状态，过热状态的液体为了再次恢复平衡，内部会均匀地产生沸腾核，同时产生大量气泡，液体体积急剧膨胀，最终也导致蒸气爆炸。对于以上两种情况，泄漏量可按存储介质瞬间全部泄漏计算。若裂口面积不大，即使有蒸气喷出，但由于储罐内压下降不急剧，液体不会达到过热状态，因此不会发生蒸气爆炸。气体或蒸气的泄放可分为节流泄放(ThrottlingRelease)和自由泄放(FreeExpansionRelease)。对于节流泄放，气体或蒸汽的压缩能绝大部分用来克服摩擦阻力;而对于自由泄放，则绝大部分转化为动能。节流泄放模型需要裂口的详细物理特征，因而只考虑较为简单的

自由泄放模型。［5］ 气体从裂口泄漏的速率与其流动状态有关。因此，计算泄漏量时首先要判断泄漏时气体流动属于音速还是亚音速流动，前者称为临界流，后者称为次临界流。

当下式成立时，气体流动属于音速流动：

102-⎪⎭⎫ ⎝⎛1+≤k k P P κ

当下式成立时，气体流动属于亚音速流动：

102-⎪⎭⎫ ⎝⎛1+>k k

P P κ

式中：P ——容器内介质压力，Pa ；P 0——环境压力，Pa ；κ——气体的绝热指数（热容比），即定压热容C p 与定容热容C V 之比。

气体泄漏量计算公式为：

11

12-+⎪⎭⎫ ⎝⎛+=κκκκG d G RT M AP YC Q

式中：Q G ——气体泄漏速度，kg/s ；P ——容器压力，Pa ；C d ——气体泄漏系数；A ——裂口面积，m 2；M ——分子量；R ——气体常数，J/(mol ·k)；T G ——气体温度，K ；Y ——流出系数，对于临界流Y=1.0对于次临界流按下式计算：

()()211121*********⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡+⨯⎥⎦⎤⎢⎣⎡-⨯⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡-⨯⎥⎦⎤⎢⎣⎡=-+-κκκκκκκp p P P Y

气体或蒸气泄漏系数的确定直接影响气体泄漏速度的计算。一般而言，泄漏系数的取值范围在0.6～1.0之间。按泄漏孔的形状可分[6]:圆形孔，Ｃd =1.0;三角孔，Ｃd =0.95；长形孔，Ｃd =0.90。孔口为内层腐蚀形成的渐缩孔，0.9<Ｃd <1.0；孔口为外力机械损伤形成的渐扩孔，0.6<Ｃd <0.9。环境风险评价中常用到的是前三种情况。

容器压力P 的大小同样直接影响气体的泄漏速速。对于压力液化气储罐气相小孔泄漏而言，储罐的压力可采用该物质在该温度下的蒸气压作为储罐压力。由于蒸气压是温度的函数，只要知道储罐的温度就可通过查阅物性数据手册得到该物质蒸气压。目前，最常用的蒸气压计算公式是安托因公式。对于压缩气体储罐(未液化)则不能采用上述方法进行计算，必须以储罐内的实际压力作计算的依据。

气体的绝热指数κ也是温度的函数，是定压热容C p 与定容热容C V 之比。一般情况下，物性手册中不直接提供绝热指数与温度的关系，只提供C p 与温度的关系值。为了求得绝热指数κ还需进一步知道C V 的值。对于理想气体而言，C P 与C V 之间有如下关系

[7]：

C P -C V =nR 或C P ，m -C V,m =R

C P ，m , C V,m 分别表示1mol 物质的等压热容和等容热容。

对于实际气体，在缺少等容热容C V 的实测资料时也可采用上式进行估算。

当容器内物质随泄漏而减少或压力降低而影响泄漏速度时，泄漏速率的计算比较