dnv-phast在危险化学品泄漏后果分析与定量风险评估中的应用

简介这是一份DNV PHAST软件讲解完整版PPT资料面向石油天然气、化工与能源领域的工艺安全工程师、风险分析及安全管理人员系统梳理PHAST工艺危险源分析工具的核心功能、界面结构和操作流程。资料共1个PPT文件压缩包大小约10MB。内容围绕闪火、喷火、池火、火球、爆炸及有毒气体扩散等事故后果模拟展开详细介绍了窗体结构、主操作菜单、分析树窗口、工作窗口和日志运行记录窗口并结合建模过程、参数输入、混合物料配置、后果情景选择、爆炸超压影响等关键环节进行说明。现有561人学习/下载。借助这份完整讲解读者可以快速掌握PHAST软件的建模思路与分析方法理解泄漏量、系统隔离时间、热辐射通量、爆炸下限LEL等关键参数的实际意义同时还能了解软件内含61种物质及理化特性数据库的使用方式减少自行摸索界面与流程的时间。无论是开展HAZOP后续定量评估、编写安全分析报告还是用于教学培训和项目汇报都具有较高的参考价值。1. dnv-phast 是什么从安全评价到定量后果分析Phast是DNV旗下用于危险化学品泄漏后果分析和定量风险评估QRA的工程软件中文项目里常把“dnv-phast”当成液烃罐区、长输管道和化工装置安全专篇的标准计算工具。它解决的是一个非常具体的工程问题当容器或管道发生泄漏后可燃或有毒物质如何在大气中扩散遇点火源会形成池火、喷射火、蒸气云爆炸中的哪一种在多大距离内会产生致死、受伤或设备损坏。很多工程师做HAZOP和重大危险源评估时最缺的不是风险意识而是能拿得出手的泄漏后果定量数据Phast恰好是用来补这一环的。它不追求CFD级别的流场精细度而是用成熟工程模型把几十个工况批量计算输出热辐射、超压和毒性等值线产出安全距离划定和风险决策所用的一系列后果曲线。EHS工程师、石化设计人员、安评机构审查人员和事故调查工程师都会在两条线上用到它一是具体泄漏事故的后果快速估算二是QRA中把事件频率和后果概率组合成风险曲线。2. dnv-phast 泄漏与扩散模型从基础参数到UDM模拟2.1 泄漏场景建模孔口、容器与管道在dnv-phast 里建后果分析第一步不是选池火还是爆炸而是把泄漏源定义准确。物料、储存温度、储存压力和泄漏孔直径这四项输入会先经过内部的泄漏模型算出泄漏速率和质量流量再把泄漏物交给下游的扩散模型。孔口直径、泄漏相态和容器压力三个参数中任何一个填错后面的扩散和火灾结果都会系统性偏移而且这种偏移在后期很难排查。泄漏速率的理论基础是孔口流动方程纯气相按等熵膨胀和临界流处理液体按Bernoulli方程带压液化气体则要算闪蒸比例与两相临界流。Phast对这三种相态分别建模界面上由“Release Phase”选项控制。项目中最常见的做法是按孔径分级建多个泄漏场景例如1mm、3mm、5mm、10mm、25mm、50mm、100mm、200mm直至完全断裂。孔径等级的选择一般参考项目安全专篇和业主风险标准而不是随意取一个管径值。下表是工程上常用的孔径分级和场景用途对照孔径等级 / mm典型对应部件后果用途1阀门盘根微漏长期低浓度泄漏、毒性边界5法兰垫片失效可燃气体泄漏、LFL边界25仪表管断裂中等泄漏QRA主要贡献场景100管线破裂喷射火、VCE主要场景200以上/全断容器整体破裂最坏事故工况用于应急预案孔径选得越大后果半径越大但对应频率越低。若做QRA每个孔径必须搭配独立的年泄漏频率不能把100mm的后果和1mm的频率放在一起算。2.1.1 孔口泄漏计算的温度与压力条件温度与压力直接决定相态选择。对气相泄漏容器压力决定出口马赫数压力接近临界压力比时泄漏速率不再随下游压力降低而继续增大。对带压液化气储存温度下的饱和蒸气压是判断相态的关键如果容器压力高于饱和蒸气压按气相泄漏如果等于饱和蒸气压按两相泄漏建模更符合实际。我一般会先在Phast的物料库里查看该物料在储存温度下的饱和蒸气压和沸点再对比设计压力。一个常见的坑是LPG储罐按“纯气体泄漏”建模忽略了闪蒸产生的液滴夹带。液滴在空气中继续蒸发会形成比纯气云持续时间更长的重气云团下风向可燃气浓度反而比纯气体模型高。2.2 扩散计算与UDM模型的关键气象参数泄漏出的物料进入大气后扩散计算在Phast中默认走UDMUnified Dispersion Model统一弥散模型。UDM把高斯型连续烟羽、重气沉降、浮力抬升和边界层湍流放在同一个框架内求解避免了过去重气到轻气切换点浓度突变的问题。计算输出包括沿下风向的中心线浓度、云团宽度和高度以及用于后果评估的LFL等值线、1/2 LFL等值线或毒性浓度等值线。气象条件对扩散结果的影响在很多案例里超过泄漏量本身。Phast的Weather界面需要设定风速、大气稳定度、环境温度、相对湿度和地表粗糙度。最常用的是把风速按1.5、3、5、10 m/s四档稳定度按B、D、F三类组合成12种气象工况逐一计算。2.2.1 大气稳定度与风速的组合Pasquill-Gifford稳定度分类中B类代表白天强日照下的强湍流D类为中性条件F类代表夜间稳定层结。同样泄漏量下F类1.5 m/s时烟羽横向扩散最弱下风向浓度下降最慢可燃或毒性边界最远B类5 m/s以上时湍流强浓度快速稀释影响距离明显缩短。所以单个“最大泄漏后果半径”通常一定来自稳定F类或静风场景而不是大风天。地面粗糙度也要按场地实际情况改。开阔农田可取0.03 m厂区内取0.5 m到1 m。Phast默认值通常偏开阔地形直接用于化工园区会高估影响距离。改完粗糙度后要重跑一遍扩散不要只改气象组合。2.3 一个可复现的泄漏扩散输入的JSON配置示例从Phast 8.x开始工况可以用JSON文件保存和交换这对批量建模很有用。下面是一份25mm孔径丙烷气相泄漏的配置骨架{ chemical: PROPANE, scenario: { type: LEAK, source: { phase: VAPOUR, pressure: 1000000, temperature: 298.0, hole: { diameter_mm: 25.0, orientation: HORIZONTAL } }, atmospheric: { wind_speed_ms: 3.0, stability_class: D, ambient_temperature_k: 298.0, relative_humidity_percent: 70.0 } } }pressure字段用Pa作单位1000000代表1MPa表压diameter_mm为泄漏孔直径stability_class只接受B/D/F等Pasquill字符。保存成case.json后在Phast Manager中导入即可复现同一场景。批量做气象敏感性时只需循环改wind_speed_ms和stability_class两个字段软件会自动重算全部下游后果。要注意的是贴这份JSON骨架前先检查物料库里的压缩因子Z是否随压力和温度自动更新。气体泄漏速率对Z很敏感Z偏差0.05就会让质量流量偏差接近5%高压场景下更明显。3. dnv-phast 火灾与爆炸后果计算模型选择与参数整定泄漏云团遇火源后具体走哪条后果路径取决于点火时刻、泄漏相态和空间约束度。很多计算偏差并不是软件算错而是在池火、喷射火、闪火和蒸气云爆炸四种模型之间选错了对象。这里把四类模型的关键参数和适用边界分开讲。3.1 池火计算燃料类型与地面条件池火针对液体泄漏到地面形成液池后持续燃烧的场景。Phast中需要输入的三个核心量是液池面积、燃料燃烧热和火焰热辐射分数。液池面积可以由围堰尺寸直接给定也可以由泄漏速率和蒸发模型自动推算。实际项目中围堰内的池火应按围堰等效直径建模如果按无边界自由扩散的液池算火灾半径会明显偏大。热辐射分数代表火焰总热量中转化为辐射的比例是整条计算链里不确定性最大的参数。常见烃类取0.25至0.35醇类取0.15左右烟尘大的重质燃料辐射分数反而偏低。Phast物料库内置了常用值但自定义混合物料时务必手动核对否则软件会沿用默认值导致热辐射通量偏离实际。下表是四类火焰/爆炸模型的适用边界速查模型适用场景关键输入主要输出池火液池稳定燃烧液池直径、燃烧热、辐射分数热辐射通量、火焰高度喷射火带压孔口喷出立即点燃孔径、压力、出口速度火焰长度、热辐射包络闪火可燃气云遇火无约束爆燃LFL等值线云团火焰扫过范围无超压VCE爆炸受限或部分受限空间爆燃云团体积、爆炸强度超压等值线选模型的判断点很简单液体在地面燃烧用池火气体从破口喷出被立即点燃用喷射火云团延迟点燃且无明显约束用闪火云团进入管廊或设备密集区后点燃才考虑VCE。3.2 喷射火与闪火的参数边界喷射火是带压流体从破口喷出后被点火源立即引燃的扩散火焰Phast按API 521和Chamberlain方法计算火焰长度和表面辐射功率。泄漏口直径和上游压力决定出口马赫数出口马赫数越高火焰被拉得越长。对应地地面目标点接收到的热辐射也显著增大。实际工况里仪表管断裂和阀门填料失效这类小孔径泄漏如果物料压力高喷射火长度也能达到几十米不能因为孔径小就忽略风险。闪火则完全不同。它不计算超压只输出可燃气云在LFL范围内的火灾边界。Phast在闪火模块中直接读取扩散阶段生成的LFL等值线并按点火模式给出火焰到达范围。3.2.1 喷射火热辐射的视角因子模拟在完整版Phast的喷射火结果中可以指定目标点的空间坐标软件会结合火焰形状、视角因子和大气透射率计算该点的热辐射通量。大气透射率默认取1.0即忽略大气吸收。实际工程中水蒸气和二氧化碳会吸收一部分辐射能湿热地区建议把透射率调到0.85至0.95再用于安全距离判定否则偏保守。拿到软件输出的辐射距离后可以用点源辐射模型快速复核数量级。把火焰总辐射功率按点源向四周均匀辐射反推某一热辐射阈值对应的距离fire_power_kw 50000.0 # 池火/喷射火总辐射功率来自泄漏速率与燃烧热 threshold_kw_m2 12.5 # 12.5 kW/m²人员可短暂停留的参考阈值 distance_m (fire_power_kw / (4 * 3.14159 * threshold_kw_m2)) ** 0.5 print(f点源估算距离: {distance_m:.1f} m)fire_power_kw是火焰总辐射功率的估算值threshold_kw_m2取12.5 kW/m²作为热辐射参考阈值。点源模型忽略火焰体积和视角因子只用于复核数量级Phast详细模型会给出更真实的方向性辐射场二者相差3倍以上时优先检查泄漏速率单位。3.3 蒸气云爆炸VCE与多能法参数蒸气云爆炸是所有泄漏后果中破坏力最大的路径也是评审最关注的内容。Phast的VCE模块提供TNT当量法和多能法两种算法。TNT当量法把可燃气云的爆炸能量折算成TNT质量再用TNT爆炸超压曲线查距离简单但保守。多能法把云团按空间约束度划分成不同爆炸强度的子区域更接近工程实际我一般优先用多能法。多能法需要云团体积和爆炸强度两个核心参数。Phast会自动从扩散结果中提取LFL范围内的云团体积但工程师要判断这个云团是瞬时泄漏形成的瞬时云还是连续泄漏形成的稳态云。瞬时云体积大但持续时间短连续云体积按释放时间累积泄漏时间越长参与爆炸的质量越大。3.3.1 爆炸强度与置信水平爆炸强度等级在Phast里通常取1到10。开阔农田取4到6管廊下方和设备密集的化工装置区取7到10。等级越高同一超压等值线覆盖的距离越远。0.3 bar超压是评价建筑物结构性破坏的常用阈值这条等值线会随爆炸强度从6级升到10级而显著外扩。做安全专篇时我会同时输出多能法6级和10级两条0.3 bar等值线让评审看到约束度假设带来的不确定性范围。只报一个强度等级等于把不确定性藏起来了审查时反而更容易被挑战。4. dnv-phast 毒性分析与风险计算从概率值到可控QRA可燃物质泄漏之后要考虑点火后果有毒物质泄漏则要单独处理毒性暴露。Phast在毒性模块中用概率值模型把浓度场随时间变化折算成人员伤害概率再在风险模块中叠加事件频率得到个体风险和社会风险。4.1 毒性负荷计算与Probit模型Probit模型的通用形式是Yab·ln(C^n·t)。其中C为暴露浓度t为暴露时间a、b、n为物质专属参数。Phast对持续泄漏按固定浓度和固定暴露时间计算对瞬时云团按浓度随时间积分计算完成后输出1%致死率、50%致死率等浓度-距离关系。暴露时间默认取600秒相当于工业事故中人员从发现到撤离的常用假设。如果场景设定为现场有自动喷淋和报警人员可以在3分钟内撤离可以改为180秒。改短暴露时间会显著降低致死概率但需要拿出疏散能力证据否则评审不会接受。4.1.1 常用毒性概率值参数表下面是Phast物料库中几种常见有毒物的常用Probit参数与参考暴露时间供建模时核对物质Probit n值暴露时间/s常见参考浓度/ppm氨2.060025氯2.756001硫化氢4.060010二氧化硫2.06002表格里的数值适合用来核对软件默认参数是否被意外修改不能直接当作最终安全标准引用。毒性参数的权威性直接决定安全距离项目报告中要注明Probit参数来源例如参考荷兰TNO或美国环保署的公开数据集Phast各版本的物料库也随版本更新。4.2 个体风险与社会风险的输出参数当泄漏场景配置好后果模型后QRA还要挂上事件频率。每个泄漏孔径对应一个年频率频率乘以点火概率、暴露概率和致死概率才能得到某一坐标点的个体风险值。Phast输出常见为10^-4、10^-5、10^-6三条个体风险等值线以及社会风险的FN曲线。如果某条10^-6等值线画到了厂界外先检查该场景的年泄漏频率和点火概率是否合理而不是急着改泄漏孔径。频率差一个数量级风险等值线位置变化比后果模型参数改动更明显。用一段简单脚本可以快速验证频率和风险的乘积关系leak_freq_per_meter_year 1.2e-6 # 每米每年50mm孔径泄漏频率 pipe_length_m 100 # 管线长度 ignition_prob 0.09 # 点火概率 exposure_prob 1.0 # 人员在场概率 annual_risk leak_freq_per_meter_year * pipe_length_m * ignition_prob * exposure_prob print(f该场景年风险贡献: {annual_risk:.2e})这段代码的逻辑是把泄漏频率与管线长度相乘得到该管段年泄漏次数再乘以点火概率和人员暴露概率。输出结果如果比项目风险标准高说明这条管线在QRA中是主导场景需要进一步细化孔径分布和点火源位置而不是盲目改后果模型。4.2.1 风险管理矩阵与风险等值线导入Phast风险模块支持把个体风险等值线导出为DXF或shapefile格式。把10^-5和10^-6两条等值线叠加到厂区总图上能直观看出风险是否越界。导出前先确认坐标原点设置Phast把泄漏点作为相对原点要换算到厂区大地坐标后再叠加地形图否则等值线会整体偏移甚至错到厂外道路上去。5. dnv-phast 结果工程化灵敏度分析、校准与批量导出的关键操作5.1 用灵敏度分析找出主导场景Phast的灵敏度分析功能允许同时变化孔径、风速、稳定度等多个变量自动批量重算几十个组合。不要一开始就全组合铺开先固定最极端气象工况只变化孔径的2到3个等级观察影响距离是否出现量级跳变。如果25mm孔径的热辐射等值线已经越界50mm只是把边界再推远十几米那主导场景就锁定在25mm孔径与会影响频率的那些小孔径上把精力放在频率数据核查而不是继续放大孔径。5.2 用公开实验数据校准扩散结果扩散模型的可靠性可以用大规模泄漏实验数据验证。国际上常被引用的LNG泄漏试验Burro、Coyote系列和重气扩散试验Thorney Island都有公开的浓度-距离数据适合用来校验UDM在重气阶段的浓度衰减趋势。Phast安装目录下的验证案例库通常自带实验对比文件把软件输出和试验峰值浓度放在同一张图上比较。偏差超过30%时先检查气象条件是否和试验记录一致再看泄漏源高度和地面粗糙度。源高度设置错误是校准失败最常见的原因。5.3 批量导出结果到Excel和GIS批量计算完成后建议按场景名、气象条件和距离阈值三列组织导出数据。Excel导出时把最大影响距离和对应面积做成数据透视表快速识别各场景的包络。给GIS出图时用shapefile导出并确认坐标已换算为大地坐标再拖入QGIS叠加卫星底图比CAD总图更接近真实周边环境。最后留一个排错技巧拿到一组后果数据后用点源热辐射经验公式或近似的蒸气云爆炸超压公式复核数量级。如果软件给出的辐射半径和经验公式相差3倍以上先回泄漏模型中看质量流率的单位是kg/s还是kg/minPhast在单位制转换上不会自动换算多数异常结果都出在这一步。本文还有配套的精品资源点击获取