危化品风险定量分析模型相关整理

重大危险源管理系统危化品风险定量分析模型需求分析2011年6月30日目录1术语 12参考依据 33风险评价程序 34风险评价所需数据 45模型分析 6失效频率分析 6设备设施基础失效（泄漏）频率分析 6泄漏（失效）频率的修正 7火源点火可能性分析 10泄漏后事故场景的频率分析 16失效后果分析 19潜在事故情景的描述（容器破裂，管道破裂，安全阀失灵等） 19危险物质泄漏量的计算（有毒、易燃、爆炸） 22危险物质泄漏后扩散的计算 35火灾、爆炸的计算 39事故后果影响的评估（毒性、热辐射、爆炸冲击波） 49泄漏扩散影响分析 49火灾热辐射影响分析 54冲击波超压影响分析 55个人风险、社会风险分析及等值线 57术语失效：指系统、结构或元件失去其原有包容流体或能量的能力，如泄漏。失效频率：失效事件所发生的频率，单位为/年。单元：具有清晰边界和特定功能的装置、设施或场所，在泄漏时能与其他装置及时切断。事故树：它是通过对可能造成系统失效的各种因素（包括硬件、软件、环境、人为因素）进行分析，画出逻辑框图（失效树），从而确定系统失效原因的各种可能组合方式或其发生概率。事件树分析：是根据一些规则用图形来表示由初因事件可能引起的许多事件链，以追踪事件破坏的过程及各事件链发生的概率。源项：定量风险评价中可能引起急性伤害的触发事件，如危险物质泄漏、火灾、爆炸等。闪火：在不造成超压的情况下物质云团燃烧时所发生的现象。池火灾：所谓池火灾是指贮罐中泄漏后的可燃液体遇火源发生的火灾.发生池火灾时热幅射是主要危害，热辐射效应被限制在液池周围的一个区域内。喷射火：压力易燃气体或液体从小孔释放到大气中被点燃，释放压力使之产生长形火焰而形成的稳定燃烧。火球：大量易燃液体/气体突然释放，并立即被点燃，由于初始动能和热火焰的浮力而上升，形成一个火球迅速燃烧。压力气体和液化烃罐事故会导致火球。闪蒸就是高压的饱和水进入比较低压的容器中后由于压力的突然降低使这些饱和水变成一部分的容器压力下的饱和水蒸气和饱和水两相：指气液混合物。辐射通量：单位时间内以辐射的形式发射、传播或接收的能量，单位为瓦。固体火灾：可燃固体材料的燃烧。蒸气云爆炸（VCE）：当大量可燃液体或蒸气泄漏后与空气混合，在扩散过程中如遇到点火源并发生点火时，由于存在某些特殊原因或条件，火焰加速传播，产生爆炸冲击波，产生蒸气云爆炸。容器瞬间破裂爆炸：如高温高压反应器。沸腾液体膨胀蒸气云爆炸（BLEVE）：压缩液体受火灾影响急剧气化，设备破裂，大量气体轻烃瞬间释放，形成爆炸。锅炉爆炸：压力容器爆炸。低温罐破裂爆炸：低温的液态乙烯丙烯泄漏后蒸发为气体，大量吸热，形成雾团，遇点火源爆炸。粉尘爆炸：细微颗粒积聚，遇点火源爆炸。个人风险：假设个体100%处于某一危险场所且无保护，由于发生事故而导致的死亡频率，单位为次/年。社会风险：能够引起大于等于N人死亡的事故累积频率（F），也即单位时间内（通常每年）的死亡人数。常用社会风险曲线（F-N曲线）表示。F-N曲线：表示累积频率（F）和死亡人数（N）之间关系的曲线图。ERPG：（EmergencyResponsePlanningGuidelines），意思是：紧急响应计划指南。紧急响应计划指南是由美国工业卫生协会（AIHA）所制定的标准。即在紧急情况下，人们持续暴露在有毒环境中1-24h，并完成指定任务所能接受的气体、蒸汽或烟雾的浓度（紧急暴露指导标准）该浓度是ppm浓度。目前美国工业卫生协会已公布了四十多种化学物之ERPG值。例如ERPG125，指限值1，空气中化学物质浓度为25ppm，只有不良反应，不会产生不良伤害。参考依据GB18218-2009危险化学品重大危险源辨识化工企业定量风险评价导则SY/T6714-2008基于风险检验的基础方法风险评价程序准备准备资料数据收集危险辨识失效频率分析失效后果分析风险计算风险标准报告和结论风险评价风险减缓措循环计算图：定量风险评价的基本程序如上图所示：定量风险评价是对危险进行识别、定量评价，并作出全面的、综合的分析。借助于定量风险评价所获得的数据和结论，并综合考虑经济、环境、可靠性和安全性等因素，制定适当的风险管理程序及措施，为设计、运行、安全管理及决策提供技术支持。风险评价所需数据定量风险评价过程涉及大量的数据信息，根据评价的目的和深度确定企业所需收集的资料数据,包括但不限于表中的资料数据:表：定量风险评价收集的资料数据大类小类数据资料企业生产信息危害信息危险物质存量危险物质安全技术说明书（MSDS）现有的工艺危害分析（如危险与可操作性分析（HAZOP））结果点火源信息（点火源的名称、种类、方位、数目以及出现的概率）点源，如加热炉（锅炉）、机车、工厂人员；线源，如公路、铁路、输电线路；面源，如厂区外的化工厂、冶炼厂；设计和运行数据设计说明平面布置图工艺技术规程安全操作规程工艺流程图（PFD）管道和仪表流程图（P&ID）设备相关数据管道数据运行数据减缓控制系统探测和切断系统（气体探测、火焰探测、毒性探测、电视监控、联锁切断等）消防减缓控制系统水幕减缓控制系统管理系统管理制度、操作和维护手册、培训、应急、事故调查、承包商管理、机械完整性管理、变更和作业程序等。这些大指标中有很多小指标是打分指标气候及环境信息气象条件大气参数（气压、温度、湿度、太阳辐射热等）风速大气稳定度周边环境现场周边地形现场建筑物等历史与人口数据历史数据事故案例设备失效统计资料等人口数据企业厂区内的人口分布企业厂区外的人口分布（区别白天与晚上、娱乐场所、体育馆等敏感场所人员的流动性等）模型分析失效频率分析设备设施基础失效（泄漏）频率分析在对风险辨识过程中，需要列出定量风险评价所需的压力容器、常压容器、管线、阀门、泵、压缩机等事故风险点清单。如清单中的设备不包括下表中的典型设备设施中，需要从企业事故统计及工业数据库中获取，并且要获取泄漏场景（4种）。国家安监总局提出的《化工企业定量风险评价导则》中指出SY/T6714-2008典型设备发生完全破裂以及孔泄漏的频率，如下表所示：表：典型设备的泄漏频率设备类型泄漏频率（/年，4种场景）5mm25mm100mm完全破裂单密封离心泵6×10-25×10-41×10-4双密封离心泵6×10-35×10-41×10-4塔器8×10-52×10-42×10-56×10-6离心压缩机1×10-31×10-4往复式压缩机6×10-36×10-4过滤器9×10-41×10-45×10-51×10-5翅片/风扇冷却器2×10-33×10-45×10-82×10-8换热器，壳程4×10-51×10-41×10-56×10-6换热器，管程4×10-51×10-41×10-56×10-619mm直径管道1×10-53×10-725mm直径管道5×10-65×10-751mm直径管道3×10-66×10-2102mm直径管道9×10-76×10-77×10-8152mm直径管道4×10-74×10-77×10-8203mm直径管道3×10-73×10-78×10-82×10-8254mm直径管道2×10-73×10-78×10-82×10-8305mm直径管道1×10-73×10-73×10-82×10-8406mm直径管道1×10-72×10-72×10-82×10-8>406mm直径管道6×10-82×10-72×10-81×10-8压力容器4×10-51×10-41×10-56×10-6反应器1×10-43×10-43×10-52×10-6往复泵7×10-11×10-21×10-31×10-3常压储罐4×10-51×10-41×10-52×10-5泄漏（失效）频率的修正由事故统计及工业数据库得出的泄漏频率称为通用（基础）泄漏频率。针对被评价单元的具体情况，采用设备系数（FE）和管理系统系数（FM）两项对通用（基础）泄漏频率进行修正，修正后的失效频率更能反映被评价单元的实际泄漏频率。如下式：F修正＝F通用（基础）×FE×FM式中：FE是设备系数，FM是管理系数。FE取值：根据设备单元运行环境及运行状况，可为每一设备制定设备修正系数FE。FE包括一般因素和工艺因素，一般因素包括工厂运行周期、气候条件、地质因素；工艺因素包括工艺连续性、稳定性及泄压系统等，参考下图。确定每一模块各个分项的取值，并将所有分项确定值相加，得到该设备项的最终数值。FE的取值范围为，即设备运行环境可造成装置风险水平降低或增加一个数量级，一般企业的分值为1。FE的取值参考SY/T6714-2008详见第部分。图：设备修正因素FMFM取值：FM主要是对企业的管理水平进行评级，一个企业的工艺安全管理系统的有效性会对机械完整性有显著影响。可采用SY/T6714-2008推荐的管理系统评价方法(《化工企业定量风险评价导则》中指出)。此方法采用风险工程学的理论，将管理系统评价分为13大类，总共1000分，见下表，一般美国化工企业管理系统评价分值约为50％。在国内开展评价活动，可参看该评价表的要素，并根据国内情况进行修正，确定每一项指标的分指标评价细则详见SY/T6714-2008附录中描述。表：管理系统评价表项目主题总分值1领导和管理702工艺安全信息803工艺危害性分析1004变更管理805操作规程806安全作业857培训1008机械完整性1209开工前安全审查6010应急措施6511事故调查7512承包商4513安全生产管理系统评价4014总计1000管理系数的评价分值与频率修正系数关系见下图，其中分值为实际分数的10％。将管理系统评价分转换成管理系统修正系数的推荐比例是基于如下的假设：一般的企业装置在管理系统评价上的得分为50%，并且100%分值等于总的装置风险降低一个数量级。图：管理系统分值与修正系数的关系火源点火可能性分析《化工企业定量风险评价导则》中指出应对评价单元的工艺条件、设备设施、平面布局等资料进行分析，结合现场调研，确定最坏事故场景影响范围内的的潜在点火源，并统计点火源的名称、种类、方位、数目以及出现的概率等要素。点火分为立即点火和延迟点火。立即点火的点火概率应考虑泄漏形式（瞬时泄漏或者连续泄漏）、泄漏速率和泄漏化学品活性。立即点火概率下表给出单元内装置的立即点火概率值：表：装置发生泄漏后的立即点火概率泄漏形式及泄漏速率（量）泄漏物质分类连续泄漏瞬时泄漏K1-液体气体-低活性气体-中/高活性<10kg/s<1000kg10kg/s-100kg/s1000kg-10000kg>100kg/s>10000kg根据泄漏物质的活性，可将泄漏物质分成K1-液体、低活性气体和中/高活性气体，一些物质的活性见如各表：表：部分化学品活性分类低中高1-氯2，3-环氧丙烷1，3-二氯丙烷3-氯-1-丙烯氨溴代甲烷一氧化碳氯乙烷氯甲烷氯乙烷甲烷四乙铅1-丁烷1，2-二氨基乙烷乙醛乙腈丁烷氯乙烯二甲胺乙烷乙基乙酰胺甲酸丙烷丙烯丁三醇*乙炔苯*二硫化碳*乙硫醇*环氧乙烷甲酸乙酯*甲醛*甲基丙烯酸酯*甲酸甲酯*甲基环氧乙烷*石脑油，溶剂*四氢噻吩*乙烯基乙酸盐*注1：如果不能获得某物质的详细信息，就认为此物质是高活性物质，这些物质以*符号标出。表：液体化学品分类物质种类限定条件分类极度易燃沸点<35℃K0闪点<0℃易燃非K0且闪点<21℃K1可燃21℃<闪点<K255℃<闪点<K3100℃<K4备注：分类标准参见《TNOPurpleBookGuidelinesforQuantitativeRiskAssessment》（简称《TNO紫皮书》）。对液体化学品立即点火概率的确定方法见下表：表：不同液体点火概率确定物质种类分类点火概率极度易燃K0气体（按气体的泄漏量及活性分类）易燃K1K1-液体可燃K2—K3—K4—延迟点火概率延迟点火的点火概率应考虑点火源特性、泄漏物特性以及泄漏发生时点火源存在的概率，计算公式如下：P(t)=Ppresent(1-e-ωt)式中：P(t)——0~t时间内发生点火的概率；Ppresent——点火源存在的概率；ω——点火效率（s-1），与点火源特性有关；t——时间，单位为s。点火效率可根据点火源在某一段时间内的点火概率计算得出。可参考《TNO紫皮书》中不同点火源在1min内的点火概率，见下表：表：点火源在1min内的点火概率点火源1min内的点火概率点源机动车辆火焰室外燃烧炉室内燃烧炉室外锅炉室内锅炉船危化品船捕鱼船游艇内燃机车电力机车线源输电线路100m公路注1铁路注1面源化工厂/座炼油厂/座重工业区/座轻工业按人口计算人口活动居民/人工人/人注1：发生泄漏事故地点周边的公路或铁路的点火概率与平均交通密度d有关。平均交通密度d的计算公式为：d=N×E/V……（）式中：N——每小时通过的汽车数量，单位为h-1；E——道路或铁路的长度，单位为km；V——汽车平均速度，单位为kmh-1。如果d≤1，则d的数值就是蒸气云通过时点火源存在的概率，此时P(t)=d(1-e-ωt)……（）式中ω为单辆汽车的点火效率，单位为s-1。如果d≥1，则d表示当蒸气云经过时的平均点火源数目；则在0-t时间之间发生点火的概率为：P(t)=1-e-dωt……（）式中ω为单辆汽车的点火效率，单位为s-1。注2：对某个居民区而言，0-t时间内的点火概率可由下式给出：P(t)=1-e-nωt……（）式中ω为每个人的点火效率，单位为s-1，n为居民区中存在的平均人数。注3：如果其他模型中采用不随时间变化的点火概率，则该点火概率等于1min内的点火概率。泄漏后事故场景的频率分析需要根据事件树分析，建立不同事故风险点的事件树，进行量化计算、分析。典型可燃气体泄漏事件树如下图所示：图：典型可燃气体泄漏事件树（WorldBank1985WTP55）（注：无影响是指无火灾、爆炸影响）

图：典型可燃液体泄漏事件树（WorldBank1985WTP55）图：可燃液体或气体瞬时泄漏事件树（SY/T6714-2008）图：可燃液体或气体连续泄漏事件树（SY/T6714-2008）失效后果分析潜在事故情景的描述（容器破裂，管道破裂，安全阀失灵等）在定量风险评价中，泄漏场景应同时满足以下两个条件：发生的概率≥10-8/年；至少导致1％的致死伤害概率泄漏场景可根据泄漏孔径大小分为完全破裂以及孔泄漏两大类，有代表性的泄漏场景如下表所示：表：泄漏场景泄漏场景范围代表值小孔泄漏0mm~5mm5mm中孔泄漏5mm~50mm25mm大孔泄漏50mm~150mm100mm完全破裂>150mm整个设备的直径注：当设备（设施）直径小于150mm时，取小于设备（设施）直径的孔泄漏场景以及完全破裂场景。管道破裂场景泄漏场景见“表：泄漏场景”；对于完全破裂场景，如果泄漏位置严重影响泄漏量或泄漏后果，应至少分别考虑三个位置的完全破裂：管道上游；管道中游；管道下游；对于长管线，应沿管线选择一系列泄漏点，泄漏点的初始间距可取为50m容器破裂场景固定的带压容器和储罐固定的带压容器和储罐分为三类，其泄漏场景见“表：泄漏场景”，分类如下表：表：带压容器和储罐分类类别定义例子带压容器内部绝对压力大于储存容器。分液罐、压力储罐等。工艺容器在容器内，物质发生物理性质变化（如温度或相态）的容器。如果容器仅仅发生液位的变化，则它应作为一个受压容器（不包括具有管程和壳程结构的换热器）。蒸馏塔、过滤器、冷凝器等。反应容器在容器内物质发生了化学变化的容器。如果在一个容器内发生了物质混合放热，则该容器也应作为一个反应容器。通用反应器、釜式反应器、床式反应器等。固定的常压容器和储罐当容器或储罐的内部绝对压力小于或等于MPa时，应考虑为常压容器或储罐，常见的常压容器和储罐的泄漏场景见下表：表：固定的常压容器和储罐的泄漏场景场场景泄漏孔径设备泄漏到大气中泄漏到外罐中5mm孔径泄漏25mm孔径泄漏100mm孔径泄漏完全破裂5mm孔径泄漏25mm孔径泄漏100mm孔径泄漏完全破裂单防罐（单层容器）√√√√双防罐√√√√√全防罐√半地下储罐√√√√注1地下储罐注1注1：对于地下储罐的泄漏场景，周围包围介质（如土壤）应考虑为第二级容器（外罐），如果储罐周围的包围介质允许泄漏物质快速蒸发或包围介质可能失效，则应考虑地下储罐的泄漏场景。对于半地下储罐的地下部分，应遵循地下储罐的原则。注2：低温储罐考虑为常压储罐，低温储罐的泄漏事件和常压储罐一致。安全阀失灵场景当压力释放装置的排放气直接排入大气环境中，应考虑压力释放装置的风险，其场景可取压力释放装置以最大释放速率进行排放。危险物质泄漏量的计算（有毒、易燃、爆炸）当压力释放装置的排放气直接排入大气环境中，应考虑压力释放装置的风险，其场景可取压力释放装置以最大释放速率进行排放。泄漏位置应根据设备（设施）实际情况而确定。在过程或反应容器中，当容器内同时存在气相和液相时，应至少模拟气相泄漏和液相泄漏两种场景。泄漏方向要根据设备安装的实际情况确定。如果没有准确的信息，泄漏方向宜设为水平方向，与风向相同。对于地下管道，泄漏方向宜为垂直向上。液体经管道上的孔流出质量流率为：……（）式中:Qm——质量流率，单位为Kg/s；P——管道内液体压力,单位为Pa；P0——环境压力,单位为Pa；C0——液体泄漏系数；ρ——泄漏液体密度，单位为kg/m3；A——裂口面积，单位为m2；液体泄漏系数C0是雷诺准数和孔洞直径的函数，经验数据如下：对于锋利的孔和雷诺准数大于30000时，液体泄漏系数近似取。对于这种情况，液体的流出速率不依赖于裂口的尺寸。对于圆滑喷嘴，液体泄漏系数可近似取1。对于与容器相连的管嘴（即长度与直径之比不小于3），液体泄漏系数近似取。当液体泄漏系数不知道或不能确定时，取使所计算的流量最大。液体经储罐上的孔流出瞬时质量流率为：……（）式中：Qm——质量流率，单位为Kg/s；P——储罐内液体压力,单位为Pa；P0——环境压力,单位为Pa；C0——液体泄漏系数；g——重力加速度，9.8m/s2；A——裂口面积，单位为m2；ρ——液体密度，单位为Kg/m3；hL——裂口上方液体高度，单位为m；t——泄漏时间，单位为s。液体管道断裂不可压缩液体在管道中流动，能量式如下：……（）式中：ΔP——管道两端压力差，单位为Pa；ρ——液体密度，单位为kg/m3；——液体平均瞬时流速，单位为m/s；α——无量纲速率轮廓修正系数,其取值为:对于层流,α取；对于塞流，α取；对于湍流α→；g——重力加速度，单位为m/s2；ΔZ——终止状态减去初始状态的高度差，单位为m；F——摩擦导致的机械能损失，包括来自流经管道长度的摩擦损失，适用于诸如阀门、弯头、孔、管道的进口和出口，单位kg；Ws——轴功，单位Pa·m；m——质量流速，单位kg/s；Δ函数——终止状态减去初始状态。对于有摩擦阻力的设备，摩擦损失项形式为：……（）式中：Kf——管道或管道配件摩擦导致的压差损失（无量纲）；对于流经管道的液体，压差损失项Kf为：……（）式中：——Fanning（范宁）摩擦系数，无量纲；L——管道长度，单位为m；d——流道直径，单位为m。Fanning（范宁）摩擦系数是雷诺准数Re和管道粗糙度ε的函数。下表给出了各种类型净管道的ε值。表：净管道的粗糙系数ε管道材料ε/mm管道材料ε/mm水泥覆护钢1~10熟铁混凝土~3拉制钢管铸铁玻璃0镀锌铁塑料0型钢对于层流，摩擦系数可按下式计算：……（）对于湍流，可按下式求解：……（）对于粗糙管道中发展完全的湍流，独立于雷诺准数，此时，摩擦系数可按下式计算：……（）对于光滑管道，ε=0，摩擦系数可按下式计算：……（）对于光滑管道，当雷诺准数小于10000时，……（）对于管道附件、阀门及其他流动障碍物，可采用改进的2-K方法来考虑能量损失，2-K方法根据雷诺准数和管道内径定义压差损失。……（）式中：Kf——超压位差损失（无量纲）；K1——常数（无量纲）；K∞——常数（无量纲）；——雷诺准数（无量纲）；IDinches——管道内径，单位为in。2-K常数如下表：表：附件和阀门中损失系数的2-K常数附件附件描述弯头90o标准（=1）,带螺纹800标准（=1）,采用法兰连接/焊接800长半径（=）,所有类型800斜接（=）：1焊缝（90o）10002焊缝（45o）8003焊缝（30o）8004焊缝（）8005焊缝（18o）80045o长半径（=1）,所有类型500长半径（=）500斜接，1焊缝（45o）500斜接，2焊缝（）500180o标准（=1），带螺纹1000标准（=1），采用法兰连接/焊接1000长半径（=），所有类型1000三通管作为弯头使用标准的，带螺纹500长半径，带螺纹800标准的，采用法兰连接/焊接800短分支1000贯通带螺纹200采用法兰连接/焊接150短分支100阀门闸阀、球阀或旋塞阀全尺寸，β=300缩减尺寸，β=500缩减尺寸，β=1000球心阀标准1500斜角或Y型1000隔膜阀Dam（闸坝）类型1000蝶形阀800止回阀提升阀2000回转阀1500倾斜片状阀1000对于管道的入口和出口，可改为：……（）对于管道进口，=160；对于一般的进口，=；对于边界类型的入口，=。对于管道出口，=0；=。对于高雷诺数（Re＞10000），上式中的第一个方程式可忽略，即；对于低雷诺数（Re＜50），小于50，第一项占支配地位，。物质从管道系统中流出,质量流率的求解方程如下：（1）假设：管道长度、直径和类型；沿管道系统的压力和高度变化；来自泵、涡轮等对液体的输入或输出功；管道上附件的数量和类型；液体的特性，包括密度和黏度。（2）指定初始点（点1）和终止点（点2）。（3）确定点1和点2处的压力和高度。确定点1处的初始液体流速。（4）推测点2处的液体流速，如果认为是完全发展的湍流，则不需要这一步。（5）用式（）~式（）确定管道的摩擦系数。（6）确定管道的超压位置损失[式（）]、附件的超压位差损失[式（）]、和进出口效应的超压位差损失[式（）]。将这些压差损失相加，使用式（）计算净摩擦损失项。使用点2处的高度。（7）计算式（）中的所用各项的值，并将其带入到方程中。如果式（）所用项的和等于零，那么计算结束。如果不等于零，返回到第（4）步重新计算。（8）使用方程确定质量流率。如果认为是完全发展的湍流，求解非常简单，将已知项代入到式（）中，将点2处的速度设为变量，直接求解该速速。气体经孔泄漏当式（）成立时，气体流动属音速流动；当式（）成立时，气体流动属亚音速流动。……（）……（）式中：——环境压力，单位为Pa；P——容器内介质压力，单位为Pa；——热容比，=cp/cv。音速流动的气体泄漏质量流率为：……（）亚音速流动的气体泄漏质量流率为：……（）式中：Q——气体泄漏质量流率，单位为kg／s；Cd——气体泄漏系数，与裂口形状有关，裂口形状为圆形时取，三角形时取，长方形时取；A——裂口面积，单位为m2；P——容器内介质压力，单位为Pa；M——泄漏气体或蒸气的分子量，单位为g/mol；Rg——理想气体常数，单位为J/（mol·K）；T——气体温度，单位为K；Y——气体膨胀因子，按式（）计算。……（）气体管道断裂绝热流动对于长管或沿管程有较大压差，气体流速在大部分情况下接近声速。对于涉及塞流绝热流动的情况下，已知管长（L）、内径（d）、上游压力（P1）和温度（T1），计算质量通量G步骤如下：（1）由式（）确定Fanning摩擦系数。假设是高雷诺数的完全发展的湍流。随后将验证这一假设，通常情况下该假设是正确的。（2）确定Ma1：……（）……（）式中：Ma1——马赫数；L——管道长度，单位为m；d——管道内径，单位为m；Y1——气体膨胀系数，无量纲（3）确定质量通量G：……（）……（）……（）式中：Gchoked——质量通量，单位为Kg/m2·s；P1——上游气体压力，单位为Pa；Pchoked——下游气体压力，单位为Pa；T1——上游气体温度，单位为K；Tchoked——下游气体温度，单位为K；（4）由式（）确定Pchoked，以确认处于塞流情况。等温流动对于大多数典型问题，已知管长（L）、内径（d）、上游压力（P1）和温度（T1），质量通量Gchoked计算步骤如下：（1）由式（）确定Fanning摩擦系数。假设是高雷诺数的完全按发展的湍流。随后将验证这一假设，通常情况下该假设是正确的。（2）确定Ma1：……（）（3）确定质量通量G：……（）……（）式中：Gchoked——质量通量，单位为Kg/m2·s；T——上游初始温度，单位为K；Pchoked——下游塞流压力，单位为Pa；P1——上游压力，单位为Pa；Ma1——马赫数；——热容比，=cp/cv；Rg——理想气体常数，单位为J/（mol·K）；M——物质分子量，单位为g/mol。绝热和等温管道方法得到的结果很接近，对于大多数实际情况，并不能很容易地确定热传递特性。因此选择绝热管道方法，通常能够得到较大的计算结果，适合于保守的安全设计。液池蒸发或沸腾过热液体泄漏计算应考虑闪蒸的影响，闪蒸计算标准方程见附录，当闪蒸比例大于时，不考虑形成液池。液池扩展应考虑地面粗糙度、障碍物以及液体收集系统等影响，如果存在围堰、防护堤等拦蓄区，且泄漏的物质不溢出拦蓄区时，液池最大半径为拦蓄区的等效半径。液池内液体蒸发按其机理可分为闪蒸、热量蒸发和质量蒸发三种。（1）闪蒸,见。（2）热量蒸发热量蒸发速度可按下式计算：……（）式中：A1——液池面积，单位为m2；T0——环境温度，单位为K；Tb——液体沸点，单位为K；H——液体蒸发热，单位为J/Kg；L——液池长度，单位为m；α——热扩散系数，单位为m2/s，见表；K——导热系数。单位为J/(m.·K)，见表；t——蒸发时间，单位为s；Nu——努塞尔（Nusselt）数。表：某些地面的热传递（热扩散、导热）系数地面情况导热系数/[J/(m·.K)]热扩散系数/(m2/s)水泥×10-7土地（含水8%）×10-7干涸土地×10-7湿地×10-7沙砾地×10-6（3）质量蒸发蒸发速度可按下式计算：……（）式中：α——分子扩散系数，单位为m2/s；Sh——舍伍德（Sherwood）数；A1——液池面积，单位为m2；L——液池长度，单位为m；ρ1——液体密度，单位为Kg/m3。闪蒸（1）闪蒸比例和闪蒸速度液体中闪蒸部分：……（）闪蒸速度：……（）式中：Fv——泄漏液体的闪蒸比例；TT——储存温度，单位为K；Tb——泄漏液体的沸点，单位为K；Hv——泄漏液体的蒸发热，单位为J/Kg；Cp——泄漏液体的比定压热容，单位为KJ·Kg-1·K-1；Qf——液体闪蒸速度，单位为kg/s；QL——物质泄漏速率，单位为kg/s。（2）闪蒸带走的液体量在液体闪蒸过程中，除了有一部分液体转变成气体外，还有一部分液体以液滴的形式悬浮在气体中，闪蒸带走的液体量的计算如下：a）当Fv≤时：带到空气中的液体量：……（）式中：D——带到空气中的液体量，单位为kg；地面液池内液体量：……（）式中：Ds——地面液池内液体量，单位为kg。b）当Fv＞时，液体全面带走，地面无液池形成。危险物质泄漏后扩散的计算计算扩散时，应至少考虑以下过程：射流扩散；根据扩散初始密度、Richardson数等条件选择重气扩散和非重气扩散。气云释放和扩散过程中的化学过程对气云的危险性影响很大时，应考虑气云中的化学过程。室内的容器、油罐和管道等设备泄漏后，应考虑建筑物对扩散的影响，选择模型时应考虑以下情况：建筑物不能承受物质泄漏带来的压力，可设定物质直接释放到大气中。建筑物可以承受物质泄漏带来的压力，则室外扩散源项应考虑建筑物内的源项以及通风系统的影响，其位置由通风系统的排风口位置来决定。在计算扩散时，应至少选择六种天气类别，代表稳定、中等稳定、不稳定、低风速、中风速和高风速等天气条件。当使用Pasquill大气稳定度时，应至少选择以下六种天气类别，见下表：表:选择的天气条件大气稳定度风速B中风速：3~5m/sD低风速：1~2m/sD中风速：3~5m/sD高风速：8~9m/sE中风速：3~5m/sF低风速：1~2m/s扩散计算时，应至少选择八种风向。气象统计资料宜取自评价单元附近有代表性的气象站，大气稳定度确定大气稳定度确定通常采用Pasquill分类方法,大气稳定度分为A、B、C、D、E和F六类，A类表示气象条件极不稳定，B类表示气象条件中等程度不稳定，C类表示气象条件弱不稳定，D类表示气象条件的稳定性在稳定和不稳定之间，E类表示气象条件弱稳定，F类表示气象条件中等程度稳定。大气稳定度的具体分类见下表：表：Pasquill大气稳定度的确定地面风速（m/s）白天日照夜间条件强中等弱阴天且云层薄，或低空云量为4/8天空云量为3/8<2AA~BB2~3A~BBCEF3~4BB~CCDE4~6CC~DDDD>6CDDDD表：日照强度的确定天空云层情况日照角>60o日照角<60o且>35o日照角>15o且<35o天空云量为4/8,或高空有薄云强中等弱天空云量为5/8~7/8,云层高度为2134m~4877m中等弱弱天空云量为5/8~7/8,云层高度<2134m弱弱弱Pasquill-Gifford模型扩散方程Pasquill-Gifford模型扩散方程位于地面高处的连续稳态源的烟羽……（）式中：——连续排放时，形成稳定的流场后，给定地点的污染物的浓度，单位为kg/m3；——连续排放的物料质量流量，单位为kg/s；u——风速，单位为m/s；，——侧风向和垂直风向的扩散系数，单位为m；x——下风向距离，单位为m；y——侧风向距离，单位为m；z——垂直风向距离，单位为m。位于地面高处的瞬时点源的烟团，坐标系位于地面并随烟团移动烟团中心在x=ut处，平均浓度方程为：……（）式中：——瞬时排放时，给定地点和时间t的污染物的浓度，单位为kg/m3；——瞬时排放的物料质量，单位为kg；u——风速，单位为m/s；，，——下风向，侧风向和垂直风向的扩散系数，单位为m。Pasquill-Gifford模型扩散系数表：推荐的烟羽扩散Pasquill-Gifford模型扩散系数方程（下风向距离x的单位为m）Pasquill-Gifford稳定度等级农村条件ABCDEF城市条件A~BCDE~F表：推荐的烟团扩散Pasquill-Gifford模型扩散系数方程（下风向距离x的单位为m）Pasquill-Gifford稳定度等级Pasquill-Gifford稳定度等级ABCDEF火灾、爆炸的计算火灾、爆炸发生的具体场景与物质特性、储存参数、泄漏类型、点火类型等有关，可采用事件树分析确定各种火灾、爆炸事件发生的类型及概率。池火计算池火灾火焰的几何尺寸及辐射参数按如下步骤计算：=1\*GB3①计算池直径根据泄漏的液体量和地面性质，按下式可计算最大可能的池面积。……（）式中：S——液池面积，单位为m2；W——泄漏液体的质量，单位为kg；——液体的密度，单位为kg/m3；Hmin——最小物料层厚度，单位为m。最小物料层与地面性质对应关系见下表：表：不同性质地面物料层厚度地面性质最小物料层厚度（m）草地粗糙地面平整地面混凝土地面平静的水面=2\*GB3②确定火焰高度计算池火焰高度的经验公式如下：……（）式中：L——火焰高度，单位为m；D——池直径，单位为m；mf——燃烧速率，单位为kg/（m2·s）；ρ0——空气密度，单位为kg/m3；g——重力加速度，单位为m/s2。=3\*GB3③计算火焰表面热通量假定能量由圆柱形火焰侧面和顶部向周围均匀辐射，用下式计算火焰表面的热通量：……（）式中：q0——火焰表面的热通量，单位为kw/m2；ΔHC——燃烧热，单位为kJ/kg；π——圆周率；——热辐射系数，取；——燃烧速率，单位为Kg/（m2·s）；其它符号同前。=4\*GB3④目标接收到的热通量的计算目标接收到的热通量q(r)的计算公式为：……（）式中：q(r)——目标接收到的热通量，单位为kw/m2；q0——按式（）计算的火焰表面的热通量，单位为kw/m2；r——目标到泄漏中心的水平距离，单位为m；V——视角系数。=5\*GB3⑤视角系数的计算视角系数V与目标到火焰垂直轴的距离与火焰半径之比s，火焰高度与直径之比h有关。……（）……（）……（）……（）……（）……（）……（）……（）……（）式中：A、B、J、K、VH、VV——描述方便而引入的中间变量；π——圆周率。火球计算火球几何尺寸及热辐射通量按如下步骤计算：①火球直径的计算火球直径计算公式为：……（）式中：——火球直径，单位为m；——火球中消耗的可燃物质量，单位为kg,对于单罐储存，取罐容量的50%，对于双罐储存，取罐容量的70%；对于多罐储存，取罐容量的90%。②火球持续时间的计算火球持续时间按下式计算：……（）式中：——火球持续时间，单位为s；——火球中消耗的可燃物质量，单位为kg。③目标接收到热辐射通量的计算……（）式中：——火球表面的辐射通量，单位为W/m2；对于柱形罐取270W/m2，对于球形罐取200W/m2；——目标到火球中心的平均距离，单位为m。喷射火计算喷射火热辐射通量计算如下：①火焰高度的计算火焰长度按下式计算：……（）式中：L——火焰长度，单位为m；dj——喷管直径，单位为m；CT——燃料-空气计算化学反应中燃料的摩尔系数；Tf——燃烧火焰的绝热温度，单位为K；Tj——喷射流体的绝热温度，单位为K；αT——燃料-空气计量化学反应中产生每摩尔燃烧产物所需反应物的摩尔数；Ma——空气的摩尔质量，单位为g/mol；Mf——燃料的摩尔质量，单位为g/mol。对于大多数燃料而言，CT远小于1，αT近似等于1，Tf和Tj的比值在7到9之间。②目标接收到热辐射通量的计算……（）式中：q(r)——距离r处目标接收到得热通量，单位为kW/m2；——大气传输率；——热辐射系数；——燃料的质量流速，单位为kJ/s；——燃烧热，单位为kJ/kg；——视角因子。大气传输率可按下式计算：……（）式中:——大气传输率；Pw——大气中水蒸汽的分压，单位为Pa；Xs——目标到火焰表面的距离，单位为m。大气中水蒸汽分压Pw可按下式计算：……（）式中：Pw——大气中水蒸汽的分压，单位为Pa；RH——相对湿度，%；Ta——环境温度，K。视角因子可按下式计算：……（）式中：r——目标到火焰中心的距离，单位为m。蒸气云爆炸（TNO模型）计算使用TNO模型步骤如下：（1）使用扩散模型确定气云的范围。（2）进行区域检查以确定拥挤的空间。（3）在被可燃气云覆盖的范围内，确定引起强烈冲击波的潜在源。强烈冲击波的潜在源包括：拥挤的空间和建筑物；延伸的平行平面之间的距离；管状结构内的空间；由于高压泄放导致的喷射中的燃料—空气混合物的剧烈动荡。（4）通过下列步骤，估算燃料—空气混合物所释放的能量。①认为每一个冲击波源是相互分离的；②假设全部的燃料—空气混合物都存在于部分受限或有障碍物的区域，被确定为气云中的冲击波源有助于冲击波的产生；③估算存在于被确定为冲击波源的单个区域内的燃料—空气混合物体积（估算是基于区域的全部尺寸之上的，注意可燃气云可能没有充满全部的冲击波源体积，以及设备的体积应被认为是它描绘了一个可接受的整个体积的一部分）；④通过将单个混合物的体积同×106J/m3相乘（该值是烃类和空气混合物，在平均化学组成计量下的典型燃料热值），计算每次爆炸的燃料能E(J)；⑤为每个单独冲击波指定一个代表冲击波强度的典型数字。如果假设爆轰的最大强度用数字10来代替，那么对于强烈爆炸的源强的估算是安全和保守的。然而，源强7似乎能更准确地代表真实的爆炸。剩余的未受限制和无障碍的部分气云所产生的爆炸，可通过假设低的初始强度进行模拟。对于延伸的静止的部分，假设最小的强度为1。对于多数不静止的部分，但处于低强度的动荡运动（如由于燃料释放的动量），可假设强度为3。（6）一旦估算出单个的当量燃料—空气混合物所导致的能量E和初始爆炸强度，那么在计算过Sachs比拟距离后，距离爆炸源R处的Sachs比拟爆炸侧向超压和负相持续时间，可从”TNO多能量爆炸模型的Sachs比拟超压与Sachs比拟正相持续时间”所示的爆炸图中查到。……（）式中：——填料的Sachs比拟距离（无量纲）；——距填料的距离，单位为m；E——填料的燃烧能，单位为J；P0——周围环境的大气压，单位为Pa。P0P0—大气压力c0—空气中的声速E—燃烧能量R0—燃烧半径燃烧能量比拟距离R图：TNO多能量爆炸模型的Sachs比拟超压与Sachs比拟正相持续时间爆炸侧向超压峰值和负相持续时间，可根据Sachs比拟超压与Sachs比拟正相持续时间计算。超压由下式计算：……（）负相持续时间，由下式计算：……（）式中：p0——侧向爆炸超压，单位为Pa；——Sachs比拟爆炸超压（无量纲）；——周围环境压力，单位为Pa；——负相持续时间，单位为s；——Sachs比拟负相持续时间，（无量纲）；E——填料的燃烧能，单位为J；——周围环境的声速，单位为m/s。TNT爆炸计算爆炸的发生是瞬间的，所以可以直接计算破坏半径信息。tnt当量法一般是以同心圆表现。计算步骤如下：（1）首先根据容器内所装介质的特性，选用公式计算出其爆破能量。（2）将爆破能量换算成TNT当量。（3）求出爆炸的模拟比。（4）求出在1000kgTNT爆炸试验中的相当距离。（5）.根据相当距离值在表“表：1000kgTNT爆炸时的冲击波超压”中找出距离某处的超压(中间值用插入法)，此即所求距离某处的超压。气体名称高热值气体名称高热值氢气12770乙烯64019氨气17250乙炔58985苯47843丙烷101828一氧化碳17250丙烯94375硫化氨生成SO225708正丁烷134026生成SO330146异丁烷132016甲烷39860丁烯121883乙烷70425表：1000kgTNT爆炸时的冲击波超压距离RO(m)56789101112超压△PO(Mpa)距离△RO(m)1618202530354045超压△PO(Mpa)距离△RO(m)505560657075超压△PO(Mpa)事故后果影响的评估（毒性、热辐射、爆炸冲击波）泄漏扩散影响分析根据章节中的计算结果，得出在某时某点有毒物质浓度值，根据《紧急响应计划指南》中标准ERPG划定致死区、失能区、中毒区，展现方式在结合二、三维地图。《紧急响应计划指南》中如下描述：ERPGs由工业界准备，由美国工业卫生协会出版，给出三个浓度范围：ERPG-1：空气中的最高浓度，低于该值就可以认为，几乎所有人都能暴露于其中1h，除了轻微的短暂的有害于健康的影响，或明显感觉到令人讨厌的气味，而没有其他影响。ERPG-2：空气中的最高浓度，低于该值就可以认为，几乎所有人都能暴露于其中1h，不会逐步显示出不可逆或其他严重的健康影响，或者消弱他们采取防护行动的能力。ERPG-3：空气中的最高浓度，低于该值就可以认为，几乎所有人都能暴露于其中1h，不会逐步显示出危及生命健康的影响。下表给出了应急响应计划指南中各物质：表：应急响应计划指南（ERPGs）（除非注明，所有值的单位均为10-6）化学物质ERPG-1ERPG-2ERPG-3化学物质ERPG-1ERPG-2ERPG-3乙醛102001000氰化氢NA1025丙烯醛3氟化氢542050丙烯酸250750硫化氢30100丙烯腈NA3575异丁腈1050200烯丙基氯3403002-异丙基丙烯酸氰乙酯NA1氨252001000苯501501000氢化锂25μg/m3100μg/m3500μg/m3氯苯11025甲醇20010005000溴15氯甲烷NA40010001,3-丁二烯10505000二氯甲烷2007504000丙烯酸丁酯25250异氰酸甲酯5异氰酸丁酯1甲硫醇25100二硫化碳150500甲基三氯硅烷315四氯化碳20100750一甲胺10100500氯气1320全氟异丁烯NA三氟化氯110苯酚1050200氯乙酰氯110光气NA1三氯硝基甲NA3五氧化二磷5mg/m325mg/m3100mg/m3氯磺酸2mg/m310mg/m330mg/m3环氧丙烷50mg/m3250mg/m3750mg/m3三氟氯乙烯20100300苯乙烯5025010002-丁烯醛21050磺酸21030乙硼烷NA13二氧化硫315双烯酮1550四氟乙烯200100010000二甲胺1100500四氯化钛5mg/m3300mg/m31000mg/m3二甲基氯硅烷525甲苯503001000二甲基二硫醚50250三甲胺100500表氯醇220100六氟化溴5mg/m315mg/m330mg/m3环氧乙烷NA50500乙酸乙烯酯575500甲醛11025乙酸535250六氯丁二烯31030乙酸酐15100六氟丙酮NA1503-氯丙烯340300六氟环丙烷1050500砷化氢NA氯化氢320100苯甲酰氯520铍NA25mg/m3100mg/m3二氯甲醚ID三氟化硼2mg/m330mg/m3100mg/m3无水肼530乙酸丁酯52003000盐酸320150丁基异氰酸酯1氢氰酸NA1525一氧化碳200350500过氧化氢1050100二氧化氯NA3硒化氢NA2一氯二氟乙烷100001500025000氯甲酸异丙酯ID520三氟甲烷NA505000碘5氯甲基甲醚NA10顺丁烯二酸酐220硝基三氯甲烷汞NA氯化氰NA4溴甲烷NA502001,2-二氯乙烷50200300氯甲酸甲酯NA252,4-二氯酚220甲基异氰酸酯二聚环戊二烯575二苯甲撑二异氰酸酯mg/m32mg/m325mg/m31,1-二氟乙烷100001500025000硝酸1678二乙烯酮1520二氧化氮11530N,N-二甲基甲酰胺2100200三氟化氮NA400800二甲硫醚100050001-辛烯40800b20003-氯-1,2-环氧丙烷520100四氯乙烯1002001000丙烯酸乙酯30300磷化氢NA5氯甲酸乙酯ID510三氯化磷315异辛醇100200四氯化硅537氟520氢氧化钠mg/m35mg/m350mg/m3氯磺酸2mg/m310mg/m330mg/m3锑化氢ID呋喃甲醛210100正硅酸乙酯25100300戊二醛15四氢呋喃1005005000六氟-1,3-丁二烯1310正硅酸甲酯NA1020六氟丙烯1050500氯化亚砜2101-己烯NA50050001,1,1-三氯乙烷3507003500三氯乙烯1005005000三氯硅烷1325三甲基氯硅烷320150氯乙烯500500020000乙烯三氯硅烷5501,1-二氯乙烯ID5001000八氧化三铀ID10mg/m350mg/m3二氧化铀ID10mg/m330mg/m3六氟化铀5mg/m315mg/m330mg/m3三氧化铀IDmg/m33mg/m3三乙氧基硅烷410三甲氧基甲硅烷25甲苯-2,4(2,6-)二异氰酸酯甲基丙烯酸异氰基乙酯ID1异戊二烯5100040002,2-二氯-1,1,1-三氟乙烷ID100010000二乙基苯10100500四羟基氢化锢ID1,1,1,2-四氟-2-氯乙烷1000500010000氯乙酰氯510邻氯苄叉缩丙二腈mg/m3mg/m325mg/m3亚乙基降冰片烯100500注1：NA表示尚未分析；ID表示数据不充分。a：表示25%的最低爆炸下限；b：表示10%的最低爆炸下限。上述物质的ERPG数值是2008年1月1日更新的数据。当没有物质的ERPG数据时的确定方法当物质的ERPG值没有时，可利用下面的方法来取得代替值。ERPG-1：(1）使用气味阈值；(2）使用十分之一ERPG-2。ERPG-2：(1）使用STEL（短时间接触容许浓度）或最高浓度限值Cs；(2）使用3倍TWA（时间加权平均容许浓度）值。ERPG-3：(1）LC-50值除以30；(2）使用5倍的ERPG-2。注：ERPG值定期更新，应使用最新的ERPG值。火灾热辐射影响分析火灾一般的影响主要是热辐射引起伤亡。但由于对于人体还是物质的热辐射影响量总是一定的，所以这里使用的标准很规范。根据、、中计算所得值热辐射与标准对照，得出死亡半径、重伤半径、轻伤半径，在二、三维地图上展现。表：稳定火灾下不同临界值热辐射通量的伤害效应临界热通量kW/m2破坏类型临界热通量kW/m2破坏类型加工设备破坏暴露15s的痛阀值木材燃烧（无引火）暴露20s的痛阀值，一度灼伤暴露5s后人严重灼伤PVC绝热电缆破坏木材被引燃暴露1min的痛阀值暴露8s的痛阀值，20s后二度灼伤长时间暴露无不适感冲击波超压影响分析根据从、计算出的超压值，与下表标准比对找出对人员和建筑场的伤害、破坏作用评价。表：基于超压的普通建筑物破坏评价（近似值）压力破坏/psig/kPa令人厌恶的噪声（137dB,或低频10-15Hz）已经处于疲劳状态下的大玻璃突然破碎非常吵的噪声（143dB）、音爆、玻璃破裂处于疲劳状态的小玻璃破裂玻璃破裂的典型压力“安全距离”（低于该值，不造成严重损坏的概率为）；抛射物极限；屋顶出现某些破坏；10%的窗户玻璃被打碎受限较小的建筑物破坏~~大窗户和小窗户通常破碎；窗户框架偶尔遭到破坏房屋建筑物受到较小的破坏房屋部分破坏，不能居住1~2~石棉板粉碎，钢板或铝板起皱，紧固失效，扣件失效，木板固定失效、吹落钢结构的建筑物轻微变形2房屋的墙和屋顶局部坍塌2~3~没有加固的水泥或煤渣石块粉碎低限度的严重结构破坏房屋的砌砖有50%被破坏3工厂建筑物内的重型机械（3000lb）遭到少许破坏；刚结构建筑变形，并离开基础3~4~无框架、自身构架钢面板建筑破坏；原油储罐破裂4轻工业建筑物的覆层破裂5木制的柱折断；建筑物被巨大的水压（40000lb）轻微破坏5~7~房屋几乎完全破坏7装满的火车翻到7~8~未加固的~mm厚的砖板被剪切，或弯曲而失效9装满的火车货车车厢被完全破坏10建筑物可能全部遭到破坏；重型机械工具（7000lb）被移走并遭到破坏，非常重的机械工具（12000lb）幸免3002068有限的爆坑痕迹个人风险、社会风险分析及等值线定量风险评价风险度量分为个体风险和社会风险。个体风险可表现为个人风险等高线，社会风险可表现为F-N曲线和潜在生命损失（PLL）。个体风险和社会风险结果应满足：个体风险应在标准比例尺地理图上以等高线的形式给出，宜表示出频率大于10-8/年的个体风险等高线。社会风险应绘制F-N曲线。在计算个体风险和社会风险时，应对评价区域进行计算网格划分，遵循的原则为：网格单元的尺寸大小取决于当地人口密度和事故影响范围，网格尺寸应尽可能小而不会影响计算结果；确定每个网格单元的人员数量时，假设网格单元内部有相同的人口密度；将点火概率分配到每一个网格单元，如果网格中有多个点火源，则将所有的点火源合并成处于网格单元中心的单个点火源。对于个体风险可只考虑人员处于室外的情况，社会风险应考虑人员处于室外和室内两种情况。在计算个体风险和社会风险时，应按下式对死亡概率进行修正：……（9）……（10）式中：——死亡概率；——个体风险计算时的死亡概率；——社会风险计算时的死亡概率；——个体风险计算时的死亡概率修正因子；——社会风险计算时的死亡概率修正因子。的取值见下表：表：死亡概率修正因子取值场景室外室外室内爆炸爆炸超压≥MPa111爆炸超压~MPa注1爆炸超压≤MPa000闪火范围内111闪火范围外000热辐射通量<KW/h火球10.14a0喷射火10.14a0池火10.14a0热辐射通量≥KW/h火球111喷射火111池火111毒性111b）注1：爆炸超压~MPa半径区域的室外人员的死亡概率为0；在计算社会风险时，室内人员需考虑建筑物破坏的影响，死亡概率为%。a）：当计算社会风险时，通常认为在衣服着火以前，室外人员因受到衣服的保护而减弱了热辐射的影响，与没有衣