天然气放空火炬高度与热辐射设计计算书

天然气放空火炬高度与热辐射设计计算,,,,,,,,,

序号,计算项目,符号表示,参数来源,计算公式,国际单位,计算结果,参数说明,规范依据,其他说明

一、基础设计参数输入,,,,,,,,,

1,火炬气最大质量排放量,W,设计值,工艺专业提供最大泄放量,kg/h,120000.00,紧急放空工况下最大气体泄放量，考虑多套装置安全阀同时起跳泄放的最不利工况,APIRP521-2014第5章,可根据实际工程规模调整该输入值

2,火炬气操作温度,T,设计值,泄放阀后温度+管道温升,K,293.15,考虑节流降温效应及管道环境温升，取20℃即293.15K作为计算温度,石油化工工艺设计手册,夏季高温工况取较高值进行保守计算

3,火炬气平均分子量,M,计算值,Σyi·Mi（组分加权）,kg/kmol,17.50,以甲烷为主的天然气，甲烷分子量16.04，含乙烷、丙烷等重烃时平均分子量增大,天然气工程设计手册,纯甲烷为16.04，典型管输天然气17~19

4,火炬气绝热指数,k,查表值,定压比热/定容比热=Cp/Cv,无量纲,1.32,天然气绝热指数，甲烷约1.32，随重烃组分含量增加略有降低,工程热力学物性手册,一般烃类气体取值范围1.2~1.4

5,气体压缩因子,Z,查表值,实际气体状态方程校正系数,无量纲,0.98,低压高温下实际气体性质接近理想气体，压缩因子接近1，取0.98保守计算,天然气压缩因子计算标准,常压放空条件下近似可取1.0

6,当地大气压力,P,实测值,当地气象站年平均大气压,kPa(A),96.00,项目所在地大气压，海拔较高地区大气压较低，取较低值保守计算,GB50009建筑结构荷载规范,标准大气压为101.325kPa，海拔每升高1000m约降10kPa

7,火炬气低发热值,HL,查表值,按组分燃烧热加权计算,kJ/kg,46500.00,甲烷低发热值约50000kJ/kg，含重烃组分时单位质量热值略低,石油化工设计手册物性数据,典型天然气取值范围42000~50000kJ/kg

8,设计马赫数,Ma,选用值,全装置紧急放空取0.5,无量纲,0.50,APIRP521规定：场/站发生事故造成全装置停工紧急放空时马赫数可取0.5,APIRP521-2014第6.2.1条,马赫数过大易导致火焰吹熄脱火，过小易发生回火

9,火炬筒顶设计风速,v_wind,设计值,当地10m高处50年一遇10min平均最大风速,m/s,25.00,用于计算火焰倾斜角及水平偏移量，大风下火焰倾斜严重，热辐射影响范围增大,GB50009-2012建筑结构荷载规范,沿海及大风地区可取30m/s以上

10,大气相对湿度,RH,实测值,当地年平均相对湿度,%,60.00,影响大气对热辐射的吸收率，湿度越大，大气对热辐射的吸收越强，透过率越低,工业设备及管道绝热设计规范,干燥地区透过率高，热辐射影响距离更远

11,火焰辐射率,F,选用值,气体扩散火焰取值范围,无量纲,0.25,烃类气体扩散火焰辐射系数，火焰越明亮、含碳烟越多取值越大,APIRP521-2014第7.2条,无烟火焰取0.15，含烟雾火焰取0.3

12,受热点距地面高度,h,设计值,按人员操作面或设备顶标高,m,1.50,人员站立操作面按1.5m标高计算，设备顶部按实际标高计算,GB50183-2015石油天然气工程设计防火规范,可根据受热点类型不同分别调整计算

13,火炬中心至受热点水平距离,R,设计值,总平面布置确定的防火间距,m,60.00,根据站场总平面布置，火炬中心与最近需要保护设施的水平距离,GB50183-2015表4.0.4,由总图运输专业根据平面布置确定

14,允许热辐射强度(人员长期暴露),q_allow,规范值,有防护时长时间停留区域,kW/m²,1.58,操作人员着合适防护服可长时间停留的热辐射强度上限,APIRP521-2014表7-1,超过此值需采取隔热措施或限制停留时间

二、火炬气物性参数计算,,,,,,,,,

15,质量流量单位换算,W_s,计算值,W/3600,kg/s,33.33,将小时质量流量换算为秒流量，用于热释放速率等国际单位制公式计算,国际单位制换算标准,工程公式统一采用SI单位制

16,火炬气密度(出口状态),ρ,计算值,(P·M)/(Z·R·T),kg/m³,0.70,采用理想气体状态方程计算火炬出口处气体密度，R为通用气体常数8.314kJ/(kmol·K),工程热力学理想气体状态方程,压力低、温度高时密度小

17,火炬气声速计算,c,计算值,√(k·(R/M·1000)·T),m/s,428.76,理想气体声速公式，声速与气体种类（绝热指数、分子量）和温度有关，与压力无关,流体力学声速计算公式,常温下甲烷声速约430m/s

18,火炬出口允许线速度,v_allow,计算值,Ma·c,m/s,214.38,设计马赫数乘以声速得到允许出口最大流速，是确定火炬直径的核心参数,APIRP521-2014第6.2条,防止回火和脱火的关键控制参数

19,火炬气摩尔流量,n_mol,计算值,W/M,kmol/h,6857.14,质量流量除以平均分子量得到摩尔流量，用于体积流量换算,化工计算基础,

20,标准状态体积流量,Q_std,计算值,n_mol×22.4/1000,Nm³/s,0.04,标准状态（0℃，101.325kPa）下的体积流量，标准摩尔体积22.4Nm³/kmol,天然气流量计量标准,

三、火炬头出口直径计算,,,,,,,,,

21,火炬气工况体积流量,Q,计算值,W_s/ρ,m³/s,47.39,质量流量除以操作条件下密度得到实际工况体积流量,流体力学连续方程,用于计算管道所需流通截面积

22,所需最小流通截面积,A_min,计算值,Q/v_allow,m²,0.22,体积流量除以允许流速得到理论最小流通截面积,APIRP521公式6-1,未考虑安全余量和管径圆整

23,火炬筒理论计算内径,d_calc,计算值,√(4·A_min/π),m,0.53,由圆面积公式反算得到理论内径，需圆整至标准管径,几何公式圆面积计算,连续方程法计算结果

24,API标准公式直径计算,d_api,计算值,3.23×10⁻⁵·√(WTZ/(PMakM)),m,0.01,APIRP521推荐的火炬头直径直接计算公式，两种方法结果应接近,APIRP521-2014公式6.2.1-1,API标准工程简化公式

25,计算结果取大值,d_max,计算值,"MAX(d_calc,d_api)",m,0.53,两种计算方法结果取较大值，确保安全,工程设计保守原则,

26,选用火炬头公称直径,DN,选用值,向上圆整至标准管径,mm,550.00,将计算内径向上圆整至最接近的标准管径（DN50的整数倍）,化工管道设计规范HG/T20592,"常用标准管径：DN200,DN250,DN300,DN350,DN400..."

27,实际火炬筒内径,d_actual,选用值,按钢管规格查取,m,0.55,根据选用的公称直径确定实际内径，按内径等于公称直径保守计算,输送流体用无缝钢管GB/T8163,考虑壁厚后内径略小，流速略高

28,实际出口截面积,A_actual,计算值,π·d_actual²/4,m²,0.24,实际选用管径对应的流通截面积,几何公式,

29,实际出口流速,v_actual,计算值,Q/A_actual,m/s,199.47,选用标准管径后的实际出口气体流速,流体力学连续方程,

30,实际运行马赫数,Ma_actual,计算值,v_actual/c,无量纲,0.47,实际运行工况下的马赫数，应不超过设计马赫数0.5,马赫数定义,

31,马赫数是否满足要求,Ma_check,判断值,"IF(Ma_actual≤0.5,""满足"",""不满足"")",无量纲,满足,校核实际马赫数是否在API标准允许范围内，防止脱火或回火,APIRP521-2014第6.2.1条,马赫数0.2~0.5为火炬稳定燃烧范围

四、火焰长度及倾斜变形计算,,,,,,,,,

32,火炬总热释放速率,Q_heat,计算值,W_s·HL·η,kW,1519000.00,单位时间内火炬燃烧释放的总热量，乘以燃烧效率0.98,APIRP521公式7-1,热辐射计算的核心输入参数，决定火焰大小和辐射强度

33,火焰释放热量(兆瓦),Q_MW,计算值,Q_heat/1000,MW,1519.00,换算为兆瓦单位，便于工程判断火炬规模,单位换算,1MW热释放功率对应火焰长度约3~4m

34,无风条件下火焰长度,L_flame,计算值,0.0148·Q_heat^0.467,m,11.40,APIRP521推荐的烃类气体扩散火焰长度经验公式，适用于Q=1~100000kW范围,APIRP521-2014图7-2及公式,火焰长度随热释放量增大而增长，但不是线性关系

35,火炬出口流速风速比,R_v,计算值,v_wind/v_actual,无量纲,0.13,环境风速与火炬气出口流速的比值，是决定火焰倾斜程度的关键参数,APIRP521火焰变形计算,比值越大，火焰被风吹得越倾斜

36,火焰水平偏移系数,dx_L,查表值,APIRP521图7-3查取,无量纲,0.10,根据流速比查API标准火焰变形图，水平方向火焰投影长度占总火焰长度的比例,APIRP521-2014图7-3火焰变形曲线,风速为0时水平偏移约10%，大风时可达50%以上

37,火焰垂直偏移系数,dy_L,查表值,APIRP521图7-3查取,无量纲,0.90,垂直方向火焰投影长度占总火焰长度的比例，水平与垂直系数平方和约为1,APIRP521-2014图7-3火焰变形曲线,大风时垂直分量显著减小，火焰被吹平

38,火焰中心水平偏移量,ΔX,计算值,dx_L·L_flame,m,1.14,火焰中心相对于火炬头出口中心的水平位移距离，下风方向为正,APIRP521第7.3.2条火焰中心位置,火焰不是从火炬头直接向上，而是向下风方向偏移

39,火焰中心垂直抬升量,ΔY,计算值,dy_L·L_flame,m,10.26,火焰中心相对于火炬头出口的垂直升高距离,APIRP521第7.3.2条,火焰总长度的垂直分量

40,火焰倾斜角,θ,计算值,arctan(ΔX/ΔY),度,6.34,火焰轴线与垂直线的夹角，直观反映火焰倾斜程度,几何三角函数计算,大风时火焰倾角可达45度甚至60度以上

41,火焰总长度校核,L_check,计算值,√(ΔX²+ΔY²),m,10.33,根据水平和垂直偏移量反算火焰总长度，应与L_flame基本一致，验证系数取值正确性,勾股定理,验证计算逻辑自洽性

五、辐射传递系数计算,,,,,,,,,

42,火焰中心X坐标,Xc,计算值,ΔX（下风方向）,m,1.14,火焰中心相对于火炬筒中心轴线的水平距离坐标,坐标系几何计算,下风方向为正方向

43,迭代计算初始假设高度,H_assume,选用值,先假设一高度进行迭代,m,30.00,迭代计算初始假设火炬高度，通过后续热辐射计算验证调整,工程迭代试算法,一般高架火炬高度范围20~60m

44,火焰中心Y坐标(初算),Yc,计算值,H_assume+ΔY,m,40.26,火焰中心距地面的垂直高度，等于火炬头高度加上火焰垂直抬升量,坐标系几何计算,火炬越高，火焰中心越高，地面热辐射越小

45,火焰中心至受热点直线距离,D,计算值,√((R-Xc)²+(Yc-h)²),m,70.48,火焰辐射点源到受热点的空间直线距离，热辐射强度与距离平方成反比,空间两点距离公式勾股定理,距离越远，热辐射衰减越快

46,几何视角因子,τ_view,计算值,点源模型保守取1.0,无量纲,1.0,保守按点源模型取视角因子为1，即受热点面对火焰的最大辐射方向,APIRP521第7.2条,点源模型是偏保守的简化计算方法

47,大气透过率计算,τ_atm,计算值,考虑距离和湿度的衰减,无量纲,0.87,大气中水蒸气和二氧化碳对热辐射的吸收衰减，距离越远、湿度越大，透过率越低,APIRP521公式7-5及图7-5,工程简化计算，精确计算需查API大气透过率图

48,综合辐射传递系数,K,计算值,F·τ_view·τ_atm,无量纲,0.22,综合考虑火焰辐射率、几何视角、大气吸收三项因素的总辐射传递系数,热辐射传输理论,典型取值范围0.12~0.25

六、热辐射强度及火炬高度计算,,,,,,,,,

49,受热点热辐射强度(初算),q_calc,计算值,K·Q_heat/(4·π·D²),kW/m²,5.30,点源模型计算初设高度下受热点的热辐射通量密度,APIRP521公式7-4点源辐射模型,与热释放量成正比，与距离平方成反比

50,所需最小火焰中心高度,Yc_min,计算值,由允许热辐射强度反算,m,116.34,根据规范允许的热辐射强度，反算满足要求所需的火焰中心最低高度,APIRP521第7.4.2条高度计算,由热辐射公式反解得到Y坐标

51,计算所需火炬头最小高度,H_calc,计算值,Yc_min-ΔY,m,106.08,火焰中心最低高度减去火焰垂直抬升量，得到火炬头出口所需最小高度,几何关系推导,尚未考虑安全系数

52,火炬高度安全系数,SF,选用值,工程设计安全余量系数,无量纲,1.20,考虑计算公式误差、气象条件波动、未来扩能等因素预留的安全裕量,石油化工安全设计导则,一般工程设计取值范围1.15~1.3

53,设计选用火炬筒总高度,H_design,选用值,向上圆整至整数,m,128.00,乘以安全系数后向上圆整，得到最终选用的火炬筒实际高度（地面至火炬头出口）,工程设计惯例,此为火炬设计核心结果参数

七、最终工况热辐射强度校核,,,,,,,,,

54,设计工况火焰中心Y坐标,Yc_design,计算值,H_design+ΔY,m,138.26,采用最终设计火炬高度后，实际火焰中心距地面高度,几何关系,

55,设计工况火焰至受热点距离,D_design,计算值,√((R-Xc)²+(Yc_design-h)²),m,148.89,设计高度下，火焰中心至受热点的实际空间直线距离,空间距离公式,

56,设计工况大气透过率,τ_atm_design,计算值,由实际距离重新计算,无量纲,0.87,根据实际设计距离重新计算大气透过率,APIRP521,

57,设计工况辐射传递系数,K_design,计算值,F·τ_view·τ_atm_design,无量纲,0.22,,,

58,设计工况受热点热辐射强度,q_design,计算值,K_design·Q_heat/(4·π·D_design²),kW/m²,1.18,最终设计高度下，受热点实际接收到的热辐射强度,,核心校核指标

59,人员长期暴露区热辐射校核,q_check,判断值,"IF(q_design≤q_allow,""满足"",""不满足"")",无量纲,满足,校核设计热辐射强度是否在人员长期暴露允许范围内,APIRP521-2014第7.4条,核心安全判定指标，必须满足

60,设备区热辐射校核(4.73kW/m²),q_check_eq,判断值,"IF(q_design≤4.73,""满足"",""不满足"")",无量纲,满足,设备区域允许热辐射强度4.73kW/m²，操作人员可穿戴防护服短时间停留（2~3分钟）,APIRP521-2014表7-1,此热辐射下20s可致痛，无灼伤

61,火炬筒底部地面热辐射强度,q_ground,计算值,K·Q_heat/(4π·(Yc_design-h)²),kW/m²,1.40,火炬筒中心点正下方地面处的热辐射强度，用于确定热辐射隔离区范围,GB50183-2015,热辐射大于1.58kW/m²区域应为防火隔离禁入区

62,1.58kW/m²安全热辐射半径,R_safe,计算值,√(KQ/(4π·1.58)-ΔY²)+ΔX,m,1.14,地面热辐射强度1.58kW/m²等值线距火炬中心的水平距离，此范围外人员可长时间安全停留,APIRP521,应以此为依据划定站场防火间距

63,4.73kW/m²热辐射半径,R_eq,计算值,√(KQ/(4π·4.73)-ΔY²)+ΔX,m,1.14,地面热辐射强度4.73kW/m²等值线距火炬中心水平距离,APIRP521,此范围内不宜设置经常有人操作的岗位

64,热辐射安全距离校核,R_check,判断值,"IF(R_safe≤R,""间距满足"",""间距不足"")",无量纲,防火间距满足,将计算得到的安全辐射半径与实际布置距离对比，校核总平面布置安全距离,GB50183-2015第4.0.4条,GB50183规定火炬防火间距应经辐射热计算确定

八、火炬系统其他安全校核,,,,,,,,,

65,火炬头出口当地声速,c_exit,计算值,同前声速公式,m/s,428.76,出口处当地声速，仅与气体组分和温度有关,流体力学,

66,正常检修放空马赫数(30%负荷),Ma_normal,计算值,30%最大流量下马赫数,无量纲,0.14,按30%最大放空量估算正常检修或小流量放空时的实际马赫数,工程经验,低流量工况需防止回火风险

67,低流速回火风险校核,Ma_low_check,判断值,"IF(Ma_normal≥0.05,""基本满足"",""低流量回火风险"")",无量纲,基本满足,低流量放空时马赫数不宜过低，防止火焰传播速度大于气流速度导致火焰回烧火炬头,APIRP521火焰稳定性,低流量工况常设置长明灯伴烧维持火焰稳定

68,火炬头出口气体动压,P_dyn,计算值,0.5·ρ·v_actual²,Pa,13993.09,高速喷射气体的动压，是产生喷注噪声的主要来源,流体力学伯努利方程,

69,火焰吹熄临界风速估算,v_blowoff,计算值,约为出口流速的30%,m/s,59.84,火焰被吹脱熄灭的临界环境风速，大于此值可能发生脱火,燃烧学火焰稳定理论,实际工程配置可靠的自动点火系统保障点火

70,火焰稳定性裕度,stability_margin,计算值,v_blowoff/v_wind,无量纲,2.39,临界吹熄风速与设计风速的比值，大于1表示在设计风速下火焰可自然稳定,APIRP521,

71,火焰稳定性校核结论,flame_stable,判断值,"IF(stability_margin≥1,""满足"",""建议设防风罩"")",无量纲,自然稳定,若自然稳定性不足，需在火炬头设置防风罩或采用强制伴烧措施,高架火炬设计规范,大风地区火炬头常设挡风板/防风罩

72,设计工况燃烧效率,η_comb,选用值,高架扩散火焰燃烧效率,%,98.00,设计良好的高架火炬在正常放空条件下的扩散燃烧效率,APIRP521燃烧性能,高效火炬燃烧效率可达98%以上，少量未燃烃主要是CO和未燃甲烷

九、放空噪声估算及评价,,,,,,,,,

73,喷注噪声声功率级估算,Lw,计算值,经验公式：Lw=L0+10lg(W·P/P0)+...,dB(A),95.56,高速气流喷射产生的空气动力性噪声声功率级，主要与质量流量和压比有关,工业噪声控制手册喷气噪声公式,紧急放空是强噪声源，声功率级可达140~160dB

74,距火炬头30m处声压级,Lp30,计算值,Lw-20lg(r)-8dB,dB(A),58.01,距火炬头水平距离30m处的噪声级，考虑球面波扩散衰减和地面反射效应,声学点声源距离衰减公式,球面波衰减规律：距离加倍声压级降6dB

75,噪声衰减至70dB(A)所需距离,R_noise70,计算值,10^((Lw-8-70)/20),m,7.55,噪声级衰减至工业区厂界标准70dB(A)所需的直线距离,GB12348工业企业厂界环境噪声排放标准,厂界噪声排放标准：昼间70dB(A)，夜间55dB(A)

76,噪声影响达标校核,noise_check,判断值,"IF(R_noise70≤R,""达标"",""噪声影响超标"")",无量纲,厂界噪声达标,紧急放空属于偶发突发噪声，国家标准允许偶发噪声超标不超过15dB(A),环境影响评价技术导则,实际工程中若放空频繁可考虑消声措施

十、未点火工况可燃气体扩散分析,,,,,,,,,

77,气体出口喷射动量,M_gas,计算值,W_s·v_actual,N,6649.04,气体喷射初始动量，决定射流初始段的扩散稀释能力,流体力学射流理论,出口速度越高，初始稀释能力越强

78,可燃下限(LEL)距离估算,R_LEL,计算值,射流扩散经验公式估算,m,55.58,甲烷浓度降至爆炸下限(5%体积比)的下风向水平距离估算,APIRP521气体扩散计算,若未点火，此范围内存在可燃气体云团

79,点火系统可靠性要求,ignition_req,判断值,必须设置可靠点火系统,无量纲,"""必须设置三级自动点火+双长明灯""",为防止大量可燃气体扩散形成可爆炸云团，必须保证放空时及时可靠点燃,GB50183-2015第6.5.7条,高架火炬必须设置高空自动点火、传焰点火、地面内传焰三级点火系统

80,点火安全措施结论,ignition_conclusion,判断值,强制安全措施,无量纲,"""设置双套独立长明灯，24小时常燃备用""",长明灯采用独立燃料气源，设置火焰检测，熄灭自动报警并重新点火,石油化工火炬系统设计规范SH3009,长明灯是防止未点火扩散的关键安全措施

十一、火雨风险及分液校核,,,,,,,,,

81,分液罐出口允许最大液滴直径,d_drop,选用值,APIRP521推荐值,μm,400.00,经分液罐分离后，火炬气中允许携带的最大液滴直径，大于此值易产生火雨,APIRP521第5.4条,大于400μm液滴可能未燃尽落到地面形成火雨

82,液滴重力沉降终端速度,v_settle,计算值,Stokes沉降公式,m/s,3.87,液滴在气体中的重力沉降终端速度，气体流速须小于此值才能有效分离,分离工程Stokes沉降定律,液滴直径越小，沉降速度越低，越难分离

83,火炬筒内气体平均流速,v_pipe,计算值,按400mm直径管道估算,m/s,377.13,火炬竖直筒体内气体平均流速，若流速过高会携带液滴,流体力学,

84,分液效果及火雨风险校核,sep_check,判断值,"IF(v_pipe<v_settle,""液滴分离良好"",""需设高效分离器"")",无量纲,气体流速过高，建议设置高效旋风分离器,校核气体流速是否小于液滴沉降速度，防止可燃液体液滴携带至火炬头燃烧形成火雨,APIRP521第5.4条,天然气放空一般不含重液滴，火雨风险较小

十二、材料选用及结构设计条件,,,,,,,,,

85,火炬筒体下部材质,mat_lower,选用值,压力容器用碳素钢,,"""Q345R正火板""",火炬下部距火焰较远、温度较低段采用优质压力容器碳素钢,GB713锅炉和压力容器用钢板,考虑腐蚀裕量不小于2mm

86,火炬筒体上部材质,mat_upper,选用值,耐热不锈钢复合板,,"""Q345R+304不锈钢复合""",火炬上部靠近火炬头高温段采用不锈钢复合钢板，耐高温腐蚀,HG/T20536化工装置管道材料设计规定,上部约10m范围采用复合板或全不锈钢

87,火炬头材质,mat_tip,选用值,耐高温奥氏体不锈钢,,"""310S(0Cr25Ni20)耐热不锈钢""",火炬头长期处于火焰高温区，采用优质耐热不锈钢，耐温1000℃以上,耐热钢牌号标准,火炬头设计寿命约5年，可更换

88,