重大危险源分级标准.docx
- 文档编号:10426745
- 上传时间:2023-02-11
- 格式:DOCX
- 页数:41
- 大小:397.74KB
重大危险源分级标准.docx
《重大危险源分级标准.docx》由会员分享,可在线阅读,更多相关《重大危险源分级标准.docx(41页珍藏版)》请在冰豆网上搜索。
重大危险源分级标准
重大危险源分级标准
(征求意见稿)
1适用范围
本规范规定了重大危险源评估分级的方法和程序。
本规范为重大危险源评估分级技术规范,适用于包括储罐区、库区、生产场所等重大危险源。
2规范性引用文件
下列文件中的条款,通过本规范的引用而成为本标准的条款。
凡是标注日期的引用文件,其随后所有的修改单(不包括勘误的内容)或修订版均不适用于本规范,然而,鼓励根据本规范达成协议的各方研究是否可使用这些文件的最新版本。
凡是不注日期的引用文件,其最新版本适用于本规范。
《中华人民共和国安全生产法》
《危险化学品安全管理条例》
《安全生产许可证条例》
《重大危险源辨识》(GB18218)
《安全评价通则》
《关于规范重大危险源监督与管理工作的通知》(安监总协调字[2005]125号)
3术语和定义
下列术语和定义适用于本规范。
3.1重大危险源majorhazardinstallations
重大危险源是指长期地或者临时地生产、搬运、使用或者储存危险物品,且危险物品的数量等于或超过临界量的单元(包括场所和设施)。
4重大危险源分级判据
重大危险源分级判据如表1所示。
表1重大危险源分级判据
危险源等级
分级判据
死亡人数
一级重大危险源
可能造成30人(含30人)以上
二级重大危险源
可能造成10一29人
三级重大危险源
可能造成3—9人
四级重大危险源
可能造成1-2人
具体判别的依据如下:
①一级重大危险源:
可能造成死亡30人(含30人)以上的重大危险源;
②二级重大危险源:
可能造成死亡10-29人的重大危险源;
③三级重大危险源:
可能造成死亡3-9人的重大危险源;
④四级重大危险源:
可能造成死亡1-2人的重大危险源。
5重大危险源死亡人数及财产损失计算方法
可能造成的死亡人数评价程序为:
1将重大危险源的周边区域划分成等间隔的网格区,用一笛卡尔坐标体系的网格覆盖城市的区域地图(如图1所示),网格间距大小取决于当地人口密度,以不影响计算结果为准。
2确定每一网格内的人员数量,通过火灾(室内火灾除外)、爆炸、毒物泄漏扩散事故后果模型计算重大危险源事故在每一网格中心处产生的热辐射、超压或
毒物浓度的数值,然后通过热辐射、冲击波超压、中毒概率函数将其其转化为造成死亡的概率。
3将每一网格中心的死亡率与人口数量相乘,即得到死亡的人数。
4将所有网格的死亡人数求和,即得到总的死亡人数。
具体用下式表示:
(1)
NDiSvi
i1
式中,N为总的死亡人数,D为第i个网格的人口密度,S为网格面积,Vi为第i个网格的个人死亡率,n为网格的数目。
图1死亡人数计算原理示意图
采用财产损失半径的方法评估事故后果造成的损失,并假定此半径内没有损失的财产与此半径外损失的财产相互抵消,或者说此半径内的财产完全损失,半径外的财产完全无损失。
财产损失半径通过火灾、爆炸事故后果模型确定。
财产损失半径按下式计算:
1/3
Ri
KiWTNT
21/6
WTNT
3175
(2)
热辐射对建筑物的影响直接取决于热辐射强度的大小及作用时间的长短,
引燃木材的热通量作为对建筑物破坏财产损失半径,计算公式如下:
tW/M
(4)
式中,q为引燃木材的热通量(W/m2),t为热辐射作用时间,即火灾持续时间(s)。
6重大危险源评价分级程序
源具体事故情景选择、后果计算及死亡概率计算过程参见附录Ao
图2重大危险源评价分级程序
重大危险源的评价分级程序如下图所示。
如果一种危险物质具有多种事故形
态,按照后果最严重的事故形态考虑,即遵循“最大危险原则”。
各类重大危险
附录A:
重大危险源事故后果模型
A.1储罐区重大事故后果分析
A.1.1储罐区的主要事故后果类型
A.1.1.1池火灾
易燃液体如汽油、苯、甲醇、乙酸乙酯等,一旦从储罐及管路中泄漏到地面后,将向四周流淌、扩展,形成一定厚度的液池,若受到防火堤、隔堤的阻挡,液体将在限定区域(相当于围堰)内得以积聚,形成一定范围的液池。
这时,若遇到火源,液池可能被点燃,发生地面池火灾。
A.1.1.2蒸气云爆炸
易燃易爆气体如H2、天然气等,泄漏后随着风向扩散,与周围空气混合成易燃易爆混合物,在扩散扩过程中如遇到点火源,延迟点火,由于存在某些特殊原因和条件,火焰加速传播,产生爆炸冲击波超压,发生蒸气云爆炸。
易燃易爆的液化气体如液化石油气、液化丙烷、液化丁烷等,其沸点远小于环境温度,泄漏后将会由于自身的热量、地面传热、太阳辐射、气流运动等迅速蒸发,在液池上面形成蒸气云,与周围空气混合成易燃易爆混合物,并且随着风向扩散,扩散扩过程中如遇到点火源,也会发生蒸气云爆炸。
A.1.1.3喷射火
对于易燃易爆气体如H2、天然气,以及易燃易爆的液化气体来说,泄漏后
可能因摩擦产生的静电立即点火,产生喷射火。
A.1.1.4沸腾液体扩展蒸气云爆炸
易燃易爆的液化气体容器在外部火焰的烘烤下可能发生突然破裂,压力平衡被破坏,液体急剧气化,并随即被火焰点燃而发生爆炸,产生巨大的火球。
这种事故被称为沸腾液体扩展为蒸气云爆炸。
A.1.1.5中毒事故
毒性的液化气体如液氯、液氨等,由于沸点小于环境温度,泄漏后会因自身
热量、地面传热、太阳辐射、气流运动等迅速蒸发,生成有毒蒸气云,密集在泄漏源周围,随后由于环境温度、地形、风力和湍流等因素影响产生漂移、扩散,
范围变大,浓度减小。
A.1.2储罐区主要事故后果模型
A.1.2.1池火灾事故后果模型
池火灾火焰的几何尺寸及辐射参数按如下步骤计算。
1计算池直径
根据泄漏的液体量和地面性质,按下式可计算最大可能的池面积。
SW/Hmin
(1)
式中,S为液池面积(m2),W为泄漏液体的质量(kg),为液体的密度(kg/m3)
Hmin为最小油层厚度(m)。
最小物料层厚度与地面性质对应关系见表1。
表1不同性质地面物料层厚度表
地面性质
最小物料层厚度(m)
草地
0.020
粗糙地面
0.025
平整地面
0.010
混凝土地面
0.005
平静的水面
0.0018
2确定火焰高度
计算池火焰高度的经验公式如下:
hLD42[m」(o.gD)]0.61
(2)
式中:
L为火焰高度(m),D为池直径(m),mf为燃烧速率(kg/m2s),po为空气密度(kg/m3),g为引力常数。
3计算火焰表面热通量
假定能量由圆柱形火焰侧面和顶部向周围均匀辐射,用下式计算火焰表面的
热通量:
周率,f为热辐射系数(可取为0.15),mf为燃烧速率(kg/m2s),其它符号同前
4目标接收到的热通量的计算
目标接收到的热通量q(r)的计算公式为:
q(r)qo(10.0581nr)V(4)
式中,q(r)为目标接收到的热通量(kw/m2),qo为由式(3)计算的火焰表面的热通量(kw/m2),r为目标到油区中心的水平距离(m),V为视角系数。
5视角系数的计算
角系数V与目标到火焰垂直轴的距离与火焰半径之比s,火焰高度与直径之
比h有关。
V,.2厂VHb
1n05
Vvtan(h/(s1)-)/sh(JK)/s
0.5
2
b(1s)/(2s)
其中A、B、J、K、Vh、Vv是为了描述方便而引入的中间变量,n为圆周率。
A.1.2.2蒸气云爆炸事故后果模型
蒸气云爆炸产生的冲击波超压是其主要危害。
冲击波超压可通过传统的TNT
当量系数法进行计算,将事故爆炸产生的爆炸能量等同于一定当量的TNT也可
根据爆炸能量直接计算。
(1)TNT当量法
1确定闪蒸系数
在热力学数据资料的基础上,用下式估算燃料的闪蒸部分
CpT
F1exp——
L
L(14)
式中,F为蒸发系数,Cp为燃料的平均比热(J/kgK),△T为环境压力下容器内温度与沸点的温差(K),L为汽化热(kJ/kg)。
2计算云团中燃料的质量:
(15)
Wf2FW
式中,W为云团中燃料的质量(kg),W为泄漏的燃料的质量(kg),F为闪蒸系数。
3计算TNT当量:
(16)
WTNTeWfHf/HTNT
式中,Wtnt为燃料的TNT当量(kg),Wf为云团中燃料的质量(kg),Hf为燃料的燃烧热(MJ/kg),Htnt为TNT的爆热(MJ/kg),a为TNT当量系数,推荐a=0.03。
4
将实际距离转化为无因次距
离:
Rr/WtNt(17)
式中,R为离爆炸点的实际距离(m),R为无因次距离(m)。
在离爆炸点距离为R处,根据相应的R值,查图1得到超压,进而预测人员受伤害和建筑受破坏的情况。
图1△PsR曲线
(2)直接计算法
在得到云团中燃料的质量的情况下,可按下式直接计算爆炸冲击波超压△p。
23
ln(Ps/Pa)0.91261.5058(1nZ)0.1675(1nZ)20.0320(1nZ)3(18)
(0.3WZ<12)
zR/^E)'3(19)
pa
E1.8WQC(20)
式中,△Ps为冲击波正相最大超压(Pa),Z为无量纲距离,Pa为环境压力,R为目标到爆源的水平距离(m),E为爆源总能量(J),a为蒸气云当量系数,一般取0.04,W为蒸气云中对爆炸冲击波有实际贡献的燃料质量(Kg),Q为燃料的燃烧热(J/Kg)。
A.1.2.3喷射火事故后果模型
加压的可燃物泄漏时形成射流,如果在泄漏裂口处被点燃,则形成喷射火。
假定火焰为圆锥形,并用从泄漏处到火焰长度4/5处的点源模型来表示。
1火焰长度计算
喷射火的火焰长度可用如下方程得到:
(21)
0444(Hem).
161.66
式中,L为火焰长度(m),He为燃烧热(J/kg),m为质量流速(kg/s)
2热辐射的通量计算
距离火焰点源为X(m)处接收到的热辐射通量可用下式表示:
(22)
fHCm
4X21000
式中,q为距离X处接收的热辐射的通量(KW/m),f为热辐射率,t为大气传输率。
大气传输率T按下式计算:
A.1.2.4沸腾液体扩展为蒸气云爆炸事故后果模型
计算主要包括如下步骤。
1火球直径
D2.665W0.327(24)
式中,D为火球直径(m),W为火球中消耗的可燃物质量(Kg)。
对单罐储存,W取罐容量的50%;对双罐储存,W取罐容量的70%;对多罐储存,W取罐容量的90%。
2火球持续时间
t1.089W0.327(25)
式中,t为火球持续时间(s),W同式(24)。
3火球抬升高度火球在燃烧时,将抬升到一定高度。
火球中心距离地面的高度H由下式估
计:
HD(26)
4火球表面热辐射能量假设火球表面热辐射能量是均匀扩散的。
火球表面热辐射能量SEP由下式
计算:
2
SEPFsmHa/(D2t)(27)
式中,Fs为火球表面辐射的能量比,Ha为火球的有效燃烧热(J/Kg)。
Fs与储罐破裂瞬间储存物料的饱和蒸气压力P(MPa)有关:
Fs0.27P0.32(28)
对于因外部火灾引起的BLEVE事故,上式中的P值可取储罐安全阀启动压力Pv(MPa)的1.21倍,即:
P1.21Pv(29)
30)
Ha由下式求得:
HaHcHvCpT
式中,He为燃烧热(J/kg),Hv为常沸点下的蒸发热(J/kg),Cp为恒压比热(J/(kg.K)),T为火球表面火焰温度与环境温度之差(K),一般来说T=1700K。
5视角系数
视角系数F的计算公式如下:
2
F((D/2)/r)2(31)
式中,r为目标到火球中心的距离(m)。
令目标与储罐的水平距离为X(m),贝U:
r(X2H2)0.5(32)
6大气热传递系数
火球表面辐射的热能在大气中传输时,由于空气的吸收及散射作用,一部分能量损失掉了。
假定能量损失比为a,则大气热传递系数a1。
a和大气中的CO2和H20的含量、热传输距离及辐射光谱的特性等因素有关。
a可由以下的经验公式来求取:
a2.02(pwr')0.09(33)
式中,pw为环境温度下空气中的水蒸气压(N/m2),r为目标到火球表面的距离(m)。
pwpW0RH(34)
0
式中,Pw为环境温度下的饱和水蒸气压(N/m2),RH为相对湿度。
r'rD/2(35)
7火球热辐射强度分布函数在不考虑障碍物对火球热辐射产生阻挡作用的条件下,距离储罐X处的热
辐射强度q(W/m2)可由下式计算:
qSEPFa(36)
A.1.2.5中毒事故后果模型
(1)泄漏模型
①液体泄漏速率模型
液体泄漏可根据流体力学中的柏努力方程计算泄漏量。
当裂口不规则时,可采取等效尺寸代替;当泄漏过程中压力变化时,则往往采用经验公式。
柏努力方程如下:
QCdA.2(PPo)2gh丫(37)
式中,Q为液体泄漏速率(kg/s),Cd为无量纲泄漏系数,是液体密度(kg/m3),A是泄漏孔面积(m2),P为罐压(Pa),Po为大气压力(Pa),g为引力常数(9.8m/s2),h为液压高度(m)。
液体出口速度可按下式计算:
Q
u
(38)
CdA
式中,u为液体出口速度(m/s),其他符号如前。
持续时间按下式计算:
ts[u°/(Cdg)](AJA)(39)
式中,U0为初始流速(m/s),At为罐内液面积(m2)。
泄漏系数Cd的取值通常可从标准化学工程手册中查到。
对于管道破裂,Cd
的典型取值为0.8。
表2为常用的液体泄漏系数数据。
表2液体泄漏系数Cd
雷诺数Re
裂口形状
圆形(多边形)
三角形
长方形
>100
0.65
0.60
0.55
<100
0.50
0.45
0.40
这个方法没有考虑泄漏速率对时间的依赖关系(压力随时间而降低以及液压
高度下降)。
因此,计算出的泄漏速率是保守的最大可能泄漏速率。
②气体泄漏模型
压力气体泄漏通常以射流的方式发生,泄漏的速度与其流动的状态有关,其
特征可用临界流(最大出口速度等于声速)或亚临界流来描述。
Perry等人用如下的关系式作为临界流的判断准则:
当式(40)成立时,气
体流动属音速流动;当式(41)成立时,气体流动属亚音速流动
k
2Flk1
k
2n
式中,Po为环境大气压力(Pa),P为容器压力(Pa),k为气体的绝热指数,
即定压比热CP和定容比热G之比。
对于很多气体,临界比值(P/Po)cr近似为2,也就是说储压近似等于大气
压力的两倍,此时流体泄漏的出口速度近似等于声速。
临界流的质量泄漏速率可
按下式计算:
k1
Mk2k1
气体呈亚音速流动时,其泄漏量为:
(43)
QYCdAP.
VRTk1
式中,Q是气体泄漏速率(kg/s),Cd为气体泄漏稀疏,A为裂口面积(ml),
M是气体相对分子质量,R是普适气体常数(8.31436Jmol-1K、,T是气体的储
存温度(K,丫为气体膨胀因子,按式(44)计算。
上述考虑的为理想气体的不可逆绝热扩散过程。
此外,没有考虑气体泄漏速率随时间的变化,因此使用初始储存条件必然导致保守的结果。
3两相流泄漏模型
Cude在1975年建议了两相流泄漏关系式。
假设源容器和泄漏点之间的管道
长度和管道直径之比L/D>12泄漏点压力与泄漏点上流压力之比Pc/P=0.55
具体计算方法如下:
第一步,按下式计算两相流的质量分数:
(45)
(TTJCp
Hv
式中,M为蒸发的液体占液体总量的比例,Tc是对应于泄漏点压力R的平衡
温度(K),T是对应于泄漏点上流压力P的平衡温度(K),CP是液体的定压比热
[J/(kg•K)〕,Hz是液体的蒸发热(J/kg)。
第二步,按下式计算两相流的平均密度:
(46)
1
MV1MV
V1
式中,、v和i分别是两相流、蒸气和液体的密度(kg/m3)。
第三步,按下式计算两相流的质量泄漏速率Q(kg/S):
QACd(p—巳)(47)
式中,G为泄漏系数,多数情况下,取G二0.8,A为裂口面积(吊),P为两相混合物的压力(Pa),Pc为临界压力(Pa)。
如果L/D<12先按前面介绍的方法计算纯液体泄漏速率和两相流泄漏速
率,再用内插法加以修正。
两相流实际泄漏速率的计算公式为:
流泄漏速率和纯液体泄漏速率(kg/S)。
(3)非重气云扩散模型
①瞬间泄漏扩散模型
②连续泄漏扩散模型
源瞬间泄漏量(kg),Q'为源连续泄漏速率(kg/s),V为风速(my7s),t为泄漏后的时间(s),x、y和z分别为x、y和z方向的扩散系数(m)。
对于连续泄漏,平均时间取10min。
其中exyz与地面的有效粗糙度
有关。
地面有效粗糙度长度如下表所示。
表3地面有效粗糙度长度表
有效粗糙度Z)<0.1m地区的扩散参数按下表选取。
表4Zo<0.1m地区的扩散参数
大气稳定度
ey/m
ez/m
A
-1/2
0.22x(1+0.0001x)
0.20x
B
0.16x(1+0.0001x)-1/2
0.12x
C
-1/2
0.11x(1+0.0001x)
-1/2
0.08x(1+0.0002x)
D
-1/2
0.08x(1+0.0001x)
-1/2
0.06x(1+0.0015x)
E
-1/2
0.06x(1+0.0001x)
-1/2
0.03x(1+0.0003x)
F
-1/2
0.04x(1+0.0001x)
-1/2
0.016x(1+0.0003x)
有效粗糙度Z)>0.1m的粗糙地形扩散系数为:
y0y
zz0z
fy(Z0)1a0Z0
式中,CyO、CzO按表4中的数值取值。
其他系数按表5取值
表5不同大气稳定度下的系数值
稳定度
A
B
C
D
E
F
a0
0.042
0.115—
0.15
0.38
0.3
0.57
b0
1.10
1.5
1.49
2.53
2.4
2.913
C0
0.0364
0.045
0.0182
0.13
0.11
0.0944
d0
0.4364
0.853
0.87
0.55
0.86
0.753
e°
0.05
0.0128
0.01046
0.042
0.01682
0.0228
f。
0.273
0.156
0.089
0.35
0.27
0.29
g0
0.024
0.0136J
0.0071
0.03
0.022
0.023
式(49)和式(50)中泄漏源有效高度是指泄漏气体形成的气云基本上变成水平状时气云中心的离地高度。
在大多数问题中,泄漏源有效高度难以与泄漏源实际高度相一致。
事实上,它等于泄漏源实际高度加泄漏源抬升高度。
泄漏源抬升高度可以用下面的公式近似计算:
1
HVsd[1.50.268Fa(TsTa)Tsd]/V(51)
或
H2.4Vsd/V(52)
式中:
△H是泄漏源抬升高度(m).Vs是气云出口速度(m/s),d是出口直径(m),V是环境风速(m/s),pa是环境大气压力(Pa),Ts是气云出口温度(K),Ta是环境大气温度(K)。
计算出泄漏源抬升高度以后,将泄漏源抬升高度与泄漏源实际高度相加就得到了泄漏源有效高度。
(3)重气云扩散模型
常用模型有盒子模型和平板模型两类。
盒子模型用来描述瞬间泄漏形成的重气云团的运动,平板模型用来描述连续泄漏形成的重气云羽的运动。
这两类模型
的核心是因空气进入而引起的气云质量增加速率方程。
1盒子模型
盒子模型使用如下假设:
I、重气云团为正立的坍塌圆柱体,圆柱体初始高度等于初始半径的一半。
II、在重气云团内部,温度、密度和危险气体浓度等参数均匀分布。
III、重气云团中心的移动速度等于风速。
重气扩散的盒子模型示意图如下图所示。
777^777777777777^77777
二云团向外的魏向运动
图2重气云团盒子模型
坍塌圆柱体的径向蔓延速度由下式确定:
1/2
Vfdr/dtg(a)/ah(53)
式中,Vf为圆柱体的径向蔓延速度(m/s),r为圆柱体半径(m),h为圆柱体
高度(m),t为泄漏后时间(s)。
等式两边同时乘以2r,上式变成下面的形式:
由于假设重气云团和环境之间没有热量交换,重气云团的浮力将守恒,即:
将上式代入式(53),积分后得到:
r2r。
22g(oa)/aVo/1/2t(56)
式中,ro为重气云团的初始半径(m),Vo为重气云团的初始体积(m3),。
为重气云团的初始密度(kg/m3)。
由于假设重气云团是圆柱体,初始高度等于初始半径的一半,因此重气云团初始半径的计算公式为:
ro(2Vo/)1/3(57)
随着空气的不断进入,云团的高度和体积也将不断变化。
云团体积随时间的变化速率由下式确定:
2
dV/dt(R)VT(2Rh)VE(58)
式中,重气云团体积VR2h,Vt和VP分别为空气从顶部和边缘进入重气云
团的速率(m/s)o
由于重气云团内部危险气体质量守恒,因此,在重气云团扩散过程中,下式存立:
22
C/COVO/V(horo2)/(hr2)(59)
式中,CO和C分别为初始时刻和t时刻重气云团内部危险物质浓度(kg/m3)任意时刻重气云团的半径按式(56)计算。
如果知道任意时刻重气云团高度的计算公式,利用上式就可计算任意时刻重气云团内部危险物质浓度。
但这里不准备采用先推导重气云团高度的计算公式,然后计算重气云团体积和危险物质浓度的方法。
而是先采用量纲分析法求重气云团的体积和浓度,然后利用上式反推重气云团的高度。
无量纲量V/Vo与x/vO/3之间存在如下函数关系:
VVo(x/Vo1/3)1.5,xVo1/3(60)
式中,x为下风向距离(m)。
它与时间、风速之间的关系为:
xVt(61)
将上式代入式(52),得到:
CCo(x/Vo1/3)1.5,xVo1/3(62)
将圆柱形重气云团的体积Vr2h代入式(60),可推导出
hVo(x/Vo1/3)1.5/(r2),xVo1/3(63)
随着空气的不断进入,重气云团的密度将不断减小,重气坍塌引起的扩散将逐步让位于环境湍流引起的扩散。
目前,判断重气坍塌过程终止的准则主要有:
I、&准则
定义(pa)/a。
&准则
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 重大 危险源 分级 标准