中国石油大学(北京)《计算机辅助安全工程》第4章_安全模拟与仿真

安全工程信息化技术基础第四章安全模拟与仿真叶迎春机械与储运工程学院2(1)模型参数任意调整(2)系统模型快速求解(3)运算结果准确可靠(4)实景仿真形象直观数字仿真的特点数字仿真的作用：(1)优化系统设计(2)降低实验成本(3)减少失败风险(4)提高预测能力3I.安全评价II.安全应急III.安全训练（虚拟现实）IV.寿命预测(ANSYS，ABAQUS提供结构疲劳分析功能)V.故障诊断数字仿真在安全工程中的应用本章目录4.5基于CFD的事故灾害模拟4.6虚拟现实技术4CFD计算流体动力学：基本思想：将时间和空间域上的连续物理场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的方程组，然后求解代数方程组，获得场变量的近似值。流动基本方程（质量守恒方程、动量守恒方程、能量守恒方程、本构方程及状态方程）控制下对流动的数值模拟。54.5基于CFD的事故灾害模拟事故灾害模拟软件•Fluent•PHAST和SAFETI（DNV开发）•FLACS（基于Unix或Linux平台）•AutoReaGas（美国CenturyDynamics开发）•AUTOHYN（美国CenturyDynamics开发）•CHARM（连续地或瞬时地释放气体或液体的模拟）•WHAZAN（定量评价石油、化工等工业危险性）64.5基于CFD的事故灾害模拟SAFETI软件2020/2/1174.5基于CFD的事故灾害模拟PHAST软件功能•泄漏模块：计算物料泄漏到大气环境中的流速和状态•扩散模块：对泄漏模块得到的结果以及天气情况进行计算来得到云团的传播扩散情况。•燃烧、爆炸模块•毒性模块2020/2/1184.5基于CFD的事故灾害模拟PHAST软件实例2020/2/119某进站管线破裂形成喷火最大影响范围图4.5基于CFD的事故灾害模拟某化肥厂发生氯气泄漏的实时扩散图2020/2/11104.5基于CFD的事故灾害模拟•Fluent软件系统目前处于世界领先地位的CFD软件之一，可用于模拟和分析具有复杂几何结构的流体流动与热传导问题。软件主要由三部分组成：•Fluent解算器•prePDF（用于计算燃烧过程的预处理软件）•Gambit(前处理程序，用于几何图形模拟和网格生成）应用：航空航天，汽车，空调供热，流体机械，化学反应，电子元件冷却，石油开采，换热器等众多领域114.5.2Fluent数值模拟Fluent的强大功能•典型应用:–汽车的内流和外流及气缸内的流动–高速航天空气动力学–火箭的内部流动–涡轮机械–反应容器的设计–台风或飓风，建筑物风环境–气泡动力学及两相流–混合容器–流化床–流体诱导噪声的预测–煤粉或者油，燃气的燃烧4.5.2Fluent数值模拟12•流动问题几乎都能解决•流体和固体的共轭传热问题•流体和固体的耦合问题（水流冲击一板条或者鸟类的翅膀扇动的过程，利用动网格）发动机的冷却水套的温度分布催化剂固体颗粒的浓度分布2020/2/11134.5.2Fluent数值模拟Velocitycontourplotu=u(t),sinusoidalinletvelocityrigidwallrigidwallcompliantwallinitialmesh动脉血管脉动流动的模拟2020/2/11144.5.2Fluent数值模拟4.5.2Fluent数值模拟电子计算机机箱的热分析与流体的粒子轨迹图搅拌机的速度流线图2020/2/1115Fluent软件模拟过程164.5.2Fluent数值模拟用fluent软件如何作论文和应掌握那些方面的知识•熟练的建模能力和划分网格的技巧•扎实的流体力学和传热学知识，包括计算传热学和计算流体力学知识•数据处理和分析的能力•编程能力•研究的方法和思路要清晰，尽可能快的吸收需要的知识为我所用2020/2/11174.5基于CFD的事故灾害模拟研究目的和意义含硫气田的钻探开发是石油工业的一个世界性难题。高含硫气井井喷失控是石油及天然气勘探开发中灾难性事故。严重危及人民的生命财产，污染良田，破坏生态环境，造成地下资源损失。开展高含硫气井井喷扩散规律研究具有重大社会意义和经济意义。2003年12月23日，重庆市开县高桥镇罗家寨16H井井喷243人死亡，数千人不同程度中毒，近十万名群众紧急疏散。直接经济损失近亿元；土壤硫化严重。•2000年10月，•天津大港油田板深七0三井井喷在气井上空形成雾状对附近大气环境造成了严重污染。1993年9月28日,华北油田一口预探井，在试油射孔作业中发生井喷村民6人死亡，24人中等中毒，440余人轻度中毒。污染庄稼近700亩；22.6万人紧急疏散。2020/2/11184.5基于CFD的事故灾害模拟研究目的和意义目前在川渝、新疆等地区发现多个大型天然气藏。该地区天然气田地质构造复杂，人口稠密，气藏埋藏深度大，具有高含H2S、高压力和高产能的“三高”特点，天然气勘探开发生产存在重大安全风险。0102030405060708090100蓬莱镇组遂宁组上沙溪庙下沙溪庙千佛崖组自流井组须家河组雷口坡组嘉陵江组飞仙关组长兴组龙潭组茅口组栖霞组梁山组黄龙组韩家店组发生次数川东地区钻完井事故及事故地层分布统计分析事故统计结果表明：对有毒有害气体扩散规律的认识欠缺是导致人员疏散延误和盲目疏散的直接原因。研究泄漏气体的扩散规律，准确预测有毒有害气体的浓度和蔓延方向，是制定应急预案和逃生路线的重要依据。2020/2/11194.5基于CFD的事故灾害模拟高含硫井喷气体扩散动态机理研究射流极值点D自由射流的流动结构与速度分布YysyTR外边界面射流核心区内边界面）射流边界层Vy﹐mVy0Vy0喷管出口R0θ初始段基本段转折截面S0井喷气体喷射速度云图射流核心区三个阶段：初始射流阶段；重力下沉阶段；自由扩散阶段。2020/2/11204.5基于CFD的事故灾害模拟数值计算软件实体模型生成2D/3D网格组分运输方程质量动量能量状态方程所支持的物理模型网格的调用和修改材料属性边界条件初始条件计算模型：紊流燃烧辐射多相流相转换动区域动网格前处理Gambit后处理TecplotFluent求解Fluent模拟计算流程2020/2/11214.5基于CFD的事故灾害模拟研究方案高含硫气体扩散模型模型控制参数整定地形模型建立组分输运方程边界条件的设置简化及算法研究网格划分技术研究复杂山区气体扩散数值模拟计算多种工况不同风速不同硫化氢含量不同大气湿度不同障碍物形式数值计算结果分析2020/2/11224.5基于CFD的事故灾害模拟边界条件计算参数及边界条件的确定网格设计解算器选择：耦合隐式求解器混合物模型：组分输运模型初始条件：计算域内各点的泄放物质浓度为零壁面条件：标准壁面函数湍流模型：标准k-ε模型计算类型：非定常计算（unsteady）（a）空间结构斜视图（b）空间结构正视图风速入口GXGXGyGZGZ井喷气体出口2020/2/11234.5基于CFD的事故灾害模拟喷射速度等值线图（硫化氢含量50%）喷射速度等值线图（硫化氢含量90%）喷射速度等值线图（硫化氢含量10%）不同硫化氢含量对高含硫气体扩散的影响（1）不同硫化氢含量（无风）2020/2/11244.5基于CFD的事故灾害模拟井喷发生稳定后模拟范围内H2S浓度小于100ppm的区域（a）硫化氢含量10%；（b）硫化氢含量50%；（c）硫化氢含量90%硫化氢含量越大，模拟区域内的安全空间越小硫化氢含量越高，含硫天然气的下部硫化氢的积聚越明显。（1）不同硫化氢含量（无风）2020/2/11254.5基于CFD的事故灾害模拟VelocityVectorsColoredByVelocityMagnitude(m/s)(Time=7.5000e+01)FLUENT6.3(2d,pbns,spe,ske,unsteady)Mar12,20102.80e+012.66e+012.52e+012.38e+012.24e+012.10e+011.96e+011.82e+011.68e+011.54e+011.40e+011.26e+011.12e+019.80e+008.40e+007.00e+005.60e+004.20e+002.80e+001.40e+000.00e+00硫化氢初始质量分数10%扩散75s后气体分布图（a）速度云图；（b）硫化氢浓度云图（c）速度矢量图（2）不同硫化氢含量（风速2m/s）2020/2/11264.5基于CFD的事故灾害模拟扩散75s后气体分布图（a）速度云图；（b）硫化氢浓度云图（黑色区域为超过100ppm的危险区域）硫化氢初始质量分数90%（2）不同硫化氢含量（风速2m/s）2020/2/11274.5基于CFD的事故灾害模拟风速5m/s风速1m/s3.2风速对高含硫气体扩散的影响2020/2/11284.5基于CFD的事故灾害模拟同一湿度，不同温度下气体扩散规律湿度对高含硫气体扩散的影响253K273K293K2020/2/11294.5基于CFD的事故灾害模拟湿度90%湿度10%湿度对高含硫气体扩散的影响同一温度，不同湿度下气体扩散规律2020/2/11304.5基于CFD的事故灾害模拟泄漏源与障碍物距离170m，障碍物尺寸为40×50m，风速2m/s硫化氢90%硫化氢10%有毒气体在障碍物背风面沉积。2020/2/11314.5基于CFD的事故灾害模拟典型高含硫气井井喷事故数值模拟井喷气体成分组成CH4H2SCO2体积分数%84.19%9.02%6.79%质量分数%68%16%16%2003年12月23日发生在重庆市开县高桥镇罗家寨16H井。国内乃至世界气井井喷史上罕见的特大井喷事故建国以来重庆历史上死亡人数最多、损失最重的一次特大事故。243人死亡，数千人不同程度中毒，数万群众紧急撤离。对周围土壤、大气等环境因素造成不同程度的污染。4.5典型高含硫气井井喷事故数值模拟2020/2/1132井喷气体成分组成CH4H2SCO2体积分数%84.19%9.02%6.79%质量分数%68%16%16%气候温暖湿润，雨量充足，日照偏少气温4.0~8.00C，相对湿度为94%~99%。无阻流量为400×104~1000×104m3/d平均风速为1m/s大气稳定度为D级或E级稳定度外界气压取101325pa。事故现场环境条件4.5典型高含硫气井井喷事故数值模拟2020/2/1133名称高度h（m）R1（m）R3（m）底面中心坐标地形形状计算区域50025002500（0，0）六棱柱A400900300（-1100，440）六棱台B250800600（-800，-1100）四棱台C100300150（1050，-1250）四棱台D200600400（1370，-230）六棱台E400550250（1650，900）四棱台F4501900100（0，-3400）六棱台重庆开县罗家16H井周边三维地形模型4.5典型高含硫气井井喷事故数值模拟2020/2/11344.5典型高含硫气井井喷事故数值模拟2020/2/11364.5典型高含硫气井井喷事故数值模拟2020/2/11374.5典型高含硫气井井喷事故数值模拟2020/2/11384.5典型高含硫气井井喷事故数值模拟2020/2/1139距离井口附近500m处井喷初期，硫化氢浓度50ppm；30min后,硫化氢浓度达到200ppm.死亡人员主要分布在距离井口500m以内的范围内;死亡的最远距离出现在距离井口1200m处。不同扩散时间的硫化氢浓度横切面图4.5典型高含硫气井井喷事故数值模拟2020/2/1140井口上方100m和200m，气体的扩散范围比地面更大。由于16H井的喷射流量非常大，导致高含硫气体的初始动量非常大，带动硫化氢向高处运动。整个毒气气云将周围的几个村庄覆盖，即使住在山顶也不能幸免。扩散30分钟不同高度硫化氢浓度横切面图4.5典型高含硫气井井喷事故数值模拟2020/2/1141本章目录4.5基于CFD的事故灾