基于cfd的井喷失控喷射火温度场和辐射场模拟仿真-火灾科学.pdf

第 l8卷 第 1期 火 灾 科 学 Vo1．18，No．1 2009年 1月 FIRE SAFETY SCm NCE Jan．2009 文章编号：1004．5309(2009)-0034-05 基于 CFD的井喷失控喷射火温度场和辐射场模拟仿真 王其华 ，李 占稳 ，雷文章 ，艾志久 ，陈高松 ，杨 芳 (1中国石油集团安全环保技术研究院，北京 100083；2西南石油大学，四川成都 610500) 摘要：采用CFD方法研究了井喷失控喷射火，阐述了模拟仿真所需的定解条件，论述了模拟仿真采用的数值计算方 法，通过模拟仿真得出了井喷失控喷射火温度场和热辐射场的分布规律。结论表明采用CFD方法研究井喷失控喷 射火是切实可行的。 关键词：井喷失控；喷射火；CFD；模拟仿真 中图分类号：TK121 文献标识码：A 柱形区域；井 口半径 r=0．1m，高度h=4m，几何模 0 引言 型如图 1所示。 在喷射火方面国外研究起步较早，而国内虽然 对喷射火也进行了大量的研究，但仍停留在用经验 公式或评价模型推导流场分布阶段，对井喷失控喷 射火的研究也很匮乏。2006年西南石油大学叶峰、 周开吉等人采用经验公式推导出了井喷失控喷射火 焰温度场、浓度场和速度场分布 ¨J。2007年 中国石 油化工股份有限公司姜巍巍等人引入了对喷射火辐 射通量影响的定量计算模型”THORNTON”，完整地 确立了喷射火评价模型 。 本文把 CFD方法引进井喷失控喷射火研究，采 图 1 几何模型 用CFX软件对其进行模拟仿真，得出了井喷失控喷 Fig．1 Geometricmodel 射火温度场和辐射场的分布规律，并通过实验验证 1．2 定解条件 了结果的正确性 。 边界条件 ： 1 模拟区域和定解条件 (1)进 口边界 (abed圆面和 kli—opn半圆柱 1．1 模拟区域 面) 模拟区域为半径R=lO0m，高度H=100m的圆 abed面是井 口边界，喷射速度 =lOOm／s，方 收稿 日期：2008—07—10；修改 日期 ：2008—10—14 基金项目：中国石油天然气集团公司重点科技项 目”陆上钻井过程重大安全事故井喷失控着火、爆炸和有毒有害气体扩散模 拟仿真研究”，编号：06B7101 作者简介：王其华 (1960一)，男，天津人 ，高级工程师，首批国家注册安全工程师、安全评价师，现在中国石油安全环保技术研 究院从事 HSE体系认证、HSE体系研究工作。 Vo1．18No．1 王其华，等：基于CFD的井喷失控喷射火温度场和辐射场模拟仿真 35 向为 z正方向；压力P=1atm；湍流动能和湍流耗散 (4)在求解联立方程组时，如果不控制变量在 ，】 1．3／2 两次求解过程中的变化幅度，就是可能导致计算过 率分别采用公式．j}=÷(“)和 = 程的发散。为了使变化幅度适当减小，要对因变量 计算。 采用欠松弛法。 kli—opn半圆柱面是风速进 口边界，对于该边 界，T：280K，P=1arm，湍流强度