石油与天然气化工  2018, Vol. 47 Issue (5): 51-57
塔一联含氮天然气对加热炉燃烧的影响
姚丽蓉 1, 赵海洋 1, 赵德银 1, 何金蓬 2     
1. 中国石化西北油田分公司石油工程技术研究院;
2. “油气藏地质及开发工程”国家重点实验室·西南石油大学
摘要:为解决塔河油田注氮开采导致天然气含氮量上升而造成的天然气燃烧炉燃烧过程不稳定的问题,基于数值模拟软件Fluent,对塔一联含氮天然气在加热炉锥形燃烧室中温度场的分布规律及燃烧产物NOx的生成规律进行了研究,并针对塔一联含氮天然气组分的特殊性,分别以燃烧效率、NOx浓度、引燃温度、发热量及沃泊指数为基准,对掺混方案进行计算,并对4种掺混方案进行对比分析。确定选用基于NOx浓度的含氮天然气掺混方案,可保证塔一联含氮天然气燃烧过程燃烧效率、引燃温度、发热量及沃泊指数要求的同时,有效减少氮氧化合物带来的环境污染,具有较强的实际意义。
关键词含氮天然气    加热炉    燃烧    数值模拟    掺混方案    
Effect of natural gas containing nitrogen on combustion of heating furnace in No.1 combined station of Tahe Oilfield
Yao Lirong1 , Zhao Haiyang1 , Zhao Deyin1 , He Jinpeng2     
1. Petroleum Engineering Technology Research Institute of Sinopec Northwest Oilfield Company, Urumqi, Xinjiang, China;
2. State Key Laboratory of Southwest Petroleum University, Chengdu, Sichuan, China
Abstract: In order to solve the unstable combustion process of natural gas combustion furnace due to increase of nitrogen content in natural gas caused by the injection of nitrogen in Tahe Oilfield, based on the numerical simulation software Fluent, which is for the tower-linked natural gas containing nitrogen in the furnace conical combustion chamber, the distribution law of temperature field and the formation law of combustion product NOx are researched in this paper. According to the particularity of the natural gas component of the nitrogen of the tower, the combustion efficiency, NOx content, ignition temperature and heating value and Webbo index used as the benchmark for the blending scheme, four blending schemes are calculated and compared. A nitrogen-containing natural gas blending scheme based on NOx content, which will ensure the combustion efficiency, ignition temperature and calorific value and Webbo index of the tower-linked natural gas containing nitrogen is selected. At the same time, this scheme will effectively reduce the environmental pollution caused by nitrogen oxides, and has strong practical significance.

塔河油田是以缝洞为主要储集体的碳酸盐岩油藏,具有非均质性强、流动区域边界复杂、缝洞尺寸差异大、接触关系复杂、几何形状各异等特征[1-3],注氮提高采收率已成为主要的增产措施,但近年来随着注氮规模的扩大,天然气中氮气含量也不断上升,在实际生产过程中将产出天然气作为塔河油田一号联(以下简称塔一联)加热炉的燃料气,于是随注氮开采的逐渐推进,天然气氮含量波动严重影响了加热炉燃烧过程的稳定性,故需对其加热炉燃烧过程进行分析,研究不同氮含量的天然气燃烧规律,以保证加热炉的稳定运行。

目前,各国学者对加热炉中天然气燃烧过程进行了大量的研究。Jin Kusaka使用多维模型与包含57种化学成分和290个基本反应的化学动力学模型相耦合的方法对双燃料天然气加热炉的燃烧与排放特性进行了研究,在计算中分析了预混合气浓度对燃烧的影响[4]。R.Lauvergne利用三维CFD软件KIVA对火花点火天然气多个方面进行了研究,确定了欧洲3种具有代表性成分的天然气的热力学特性和火焰传播速度。研究表明,三维计算可以指出天然气对各种空气动力学和热力学参数的不同敏感性[5]。窦慧莉采用STAR-CD对电控喷射稀燃天然气加热炉的混合气形成和燃烧过程进行了数值模拟,验证了数值模拟对燃烧过程研究的有效性[6]

虽然以往学者对天然气加热炉有较广泛的研究[7-15],但均未涉及过实际生产过程天然气氮含量变化对燃烧过程的影响,本文针对塔河油田含氮天然气加热炉燃烧过程,基于数值模拟软件Fluent,对塔河油田含氮天然气在加热炉锥形燃烧室中温度场的分布规律及燃烧产物NOx的生成规律进行研究,并获取掺混方案,具有较强的实际意义。

1 塔一联天然气组分波动现状

塔一联2017年1月2日至2017年7月22日天然气全组分分析报告见图 1。分析发现,天然气中C1和N2含量较高,且波动较大,其余各组分含量较少且变化较为平稳。将天然气全组分报告平均化处理得到塔一联天然气平均组分,见表 1。此时,天然气中氮的体积分数为10.17%,发热量为30.6 MJ/m3,不能满足GB 17820-2010《天然气》中的天然气发热量应不小于31.4 MJ/m3的要求。分析塔一联天然气组分波动图发现,随时间的变化,天然气氮含量不断波动,最高可达21.0%,并且随注氮开采的逐渐推进,塔一联氮含量将逐渐提高,而天然气发热量将进一步下降,由此会影响加热炉的热负荷,进而影响塔一联整体工艺运行的稳定性。

图 1     塔一联天然气组分波动图 Figure 1     Natural gas composition of No.1 combined station in Tahe Oilfield

表 1    塔一联2017年1月2日~2017年7月22日天然气平均组分 Table 1    Average composition of natural gas from 2 January 2017 to 22 July 2017 in Tahe Oilfield

2 数学模型
2.1 燃烧模型及网格划分

由于在塔一联的加热炉燃烧过程中,氧化剂和天然气是未经混合直接进入燃烧室进行反应的,属于扩散燃烧,采用Fluent中的PDF燃烧模型(description of the equilibrium mixture fraction/PDF model)进行模拟。该模型通过化学平衡方法或火焰面方法(即混即燃模型)定义反应机理,计算效率很高,且允许预测中间组分和严格的湍流化学反应[18-21],主要是为模拟进行快速化学反应的湍流扩散燃烧而设计,因此,本研究选用此模型来模拟辐射室中的燃烧。

利用Fluent软件研究含氮天然气在加热炉锥形燃烧室中温度场的分布规律及燃烧产物NOx的生成规律。建立Fluent几何模型如图 2所示。塔一联天然气加热炉额定热负荷为1 000 kW。天然气预热至650 K,以5 m/s的速度由入口进入燃烧室,燃烧后的气体经出口流出。选用k-epsilon(2 eqn)湍流模型、有限速率化学反应模型(finite-rate chemistry model)进行计算,并采用自适应网格(self adaptive mesh)提高计算精度。

图 2     加热炉锥形燃烧室几何模型 Figure 2     Geometry model of conical furnace combustion chamber

2.2 天然气燃烧过程评价参数计算模型

为了便于对比分析天然气不同氮含量的燃烧过程,在此引入发热量和沃泊指数。发热量又称为燃料发热量,是指单位质量(指固体或液体)或单位体积的燃料完全燃烧,燃烧产物冷却到燃烧前的温度(一般为环境温度)时所释放出来的热量。根据GB/T 11062-2014《天然气发热量、密度、相对密度和沃泊指数的计算方法》计算出不同氮含量下天然气的发热量,计算方式见式(1):

$ {\bar H^0}\left( {{t_1}} \right) = \sum\limits_{j = 1}^N {{x_j}\cdot} \bar H_j^0({t_1}) $ (1)

式中:H0(t1)为混合物的理想摩尔发热量,kJ/mol;xj为混合物中组分j的摩尔分数;Hj0(t1)为混合物中组分j的理想摩尔发热量,kJ/mol。

沃泊指数是燃气的热负荷指数,等于燃气的高位发热量与相对密度开方的比值,代表燃气性质对热负荷的综合影响。

$ {W_{\rm{S}}} = \frac{{{H_{\rm{S}}}}}{{\sqrt \Delta }} $ (2)

式中:WS为沃泊指数,MJ/m3HS为天然气高位发热量,MJ/m3Δ为天然气相对密度。

沃泊指数是重要的燃气参数之一,其意义在于,具有相同沃泊指数的不同的燃气成分,在相同的燃烧压力下,能释放出相同的热负荷。

3 含氮天然气燃烧规律分析

基于所建立的塔一联加热炉锥形燃烧室数学模型,获得氮含量(摩尔分数, 下同)为0%~10%的天然气完全燃烧时的温度场,发现当天然气氮含量小于10.0%时,天然气燃烧具有相似的温度分布,以不含氮气的天然气为例,完全燃烧时温度场分布如图 3所示,同理获得不同氮含量天然气完全燃烧时的温度场,进而得到天然气的燃烧最高温度随氮含量变化曲线,见图 4。分析发现,当氮含量小于4.4%时,随着天然气氮含量上升,燃烧最高温度基本不变,但当氮含量为4.5%时,燃烧最高温度骤降,当氮含量大于4.5%时,天然气最高燃烧温度下降速度变缓。计算为达到热负荷所需气量、燃烧效率、NOx浓度及过剩空气系数随天然气氮含量变化关系曲线,见图 5~图 8表 2

图 3     不含氮天然气温度场 Figure 3     Temperature field of nitrogen-free natural gas

图 4     不同氮含量天然气燃烧最高温度 Figure 4     Maximum temperature of natural gas combustion with different nitrogen content

图 5     达到热负荷所需天然气气量随天然气氮含量变化 Figure 5     Amount of natural gas required to reach the heat load varies with the natural gas nitrogen content

图 6     天然气燃烧效率随天然气氮含量变化 Figure 6     Natural gas combustion efficiency varies with natural gas nitrogen content

图 7     天然气燃烧后NOx浓度随天然气氮含量变化 Figure 7     NOx content varies with natural gas nitrogen content after combustion

图 8     天然气发热量及沃泊指数随天然气氮含量变化 Figure 8     Natural gas calorific value and Wobbe index varies with natural gas nitrogen content

表 2    天然气燃烧模拟结果 Table 2    Simulation combustion results of natural gas

分析发现,随着天然气氮含量增加,为达到热负荷所需的天然气气量和燃烧后NOx浓度总体上升,燃烧效率及天然气发热量总体下降。当天然气氮含量小于5.0%时,C1含量较多,燃烧后温度较高,N2主要生成热力型NOx,其反应速率较慢,故NOx浓度较小,天然气燃烧效率较高(大于97%),达到热负荷所需天然气气量较少;当含量大于5.0%时,N2反应机理发生改变,主要生成快速型NOx和热力型NOx,反应速率较快,且反应过程吸收大量热量,故天然气燃烧效率迅速下降(小于85%),NOx浓度和为达到热负荷而所需天然气气量急剧上升。

4 塔一联含氮天然气掺混方案

通过研究塔一联含氮天然气燃烧规律发现,天然气氮含量过高将引起天然气发热量下降、天然气无法点燃、天然气燃烧效率下降和燃烧后NOx浓度过高的问题,于是需要对天然气进行掺混处理,掺混用天然气组分如表 3所示。在以往的研究中,所处理的天然气组分较纯净,掺混过程仅需考虑提高加热炉的热负荷及燃烧效率,但对于塔一联含氮天然气的燃烧,还需考虑氮含量对于加热炉引燃温度和NOx浓度的影响,因此,本研究分别以燃烧效率、NOx浓度、引燃温度和发热量为基准对掺混方案进行计算。

表 3    掺混天然气组分 Table 3    Components of blending natural gas

4.1 基于燃烧效率的含氮天然气掺混方案

根据GB 24848-2010《石油工业用加热炉能效限定值及能效等级》,当3级燃气加热炉额定热负荷为1 000 kW,燃烧效率应不小于85%,故考虑用干气进行掺混,得到基于燃烧效率的天然气掺混比例见图 9。分析发现,基于燃烧效率的含氮天然气掺混方案,当燃料气氮含量小于4.8%时,无需掺混天然气。

图 9     基于燃烧效率的含氮天然气掺混比例随其氮含量变化 Figure 9     Blending ratio of nitrogen-containing natural gas based on combustion efficiency varies with its nitrogen content

4.2 基于NOx浓度的含氮天然气掺混方案

根据GB 13271-2014《锅炉大气污染物排放标准》,燃烧后天然气中NOx浓度应小于400 mg/m3,通过计算可知,天然气中氮气摩尔分数应小于2.0%。采用干气(不含氮)进行掺混时,掺混标准为:掺混后氮气摩尔分数不大于2.0%,计算得到基于NOx浓度的含氮天然气掺混比例见图 10。基于NOx浓度的含氮天然气掺混方案,当燃料气中氮气摩尔分数小于2.0%时,无需掺混天然气。

图 10     基于NOx浓度的含氮天然气掺混比例随其氮含量变化 Figure 10     Blending ratio of nitrogen-containing natural gas based on NOx content varies with its nitrogen content

4.3 基于引燃温度的含氮天然气掺混方案

当天然气中氮气含量过高时,天然气发热量将下降,燃烧温度逐渐降低,过低时将无法点燃。基于所建立的塔一联加热炉锥形燃烧室数学模型,获取天然气含氮浓度为30.0%的温度场见图 11。其燃烧最高温度小于C1引燃温度(810 K),故当天然气浓度超过30.0%时,将无法点燃。基于引燃温度的掺混标准为:掺混后氮气摩尔分数小于30.0%,以此计算得到含氮天然气掺混比例见图 12,在此方案中,当燃料气氮含量小于30.0%时,无需掺混天然气。

图 11     天然气氮气摩尔分数为30.0%时温度场 Figure 11     Temperature field of natural gas with 30.0% nitrogen mole fraction

图 12     基于引燃温度的含氮天然气掺混比例随其氮含量变化 Figure 12     Blenging ratio of nitrogen-containing natural gas based on ignition temperature varies with its nitrogen content

4.4 基于发热量及沃泊指数的天然气掺混方案

当天然气中氮气含量过高时,燃烧气发热量及沃泊指数将下降。根据GB 17820-2010《天然气》,天然气发热量应不小于31.4 MJ/m3;同时,根据GB/T 13611—2006《城镇燃气分类和基本特性》,10 T天然气的沃泊指数应不小于41.5 MJ/m3,故考虑用干气进行掺混,掺混标准为:掺混后氮气摩尔分数小于7.0%,计算得到基于发热量的天然气掺混比例见图 13,当燃料气氮含量小于7.5%时,无需掺混天然气。

图 13     基于发热量及沃泊指数的含氮天然气掺混比例随其氮含量变化 Figure 13     Blending ratio of nitrogen-containing natural gas based on calorific value and Wobbe index varies with its nitrogen content

对比塔一联含氮天然气的4种不同基准掺混方案发现,燃烧后NOx浓度对于天然气氮含量变化最为敏感,当塔一联天然气氮气摩尔分数达到2.0%时,已无法满足国家标准GB 13271-2014《锅炉大气污染物排放标准》对燃料产物中NOx浓度应小于400 mg/m3的相关要求,因此,选取塔一联含氮天然气燃烧掺混方案时,不能一味地追求燃烧效率和发热量及沃泊指数的提高,在国家大力推进节能减排的大环境下,须长远考虑氮氧化合物对环境的长远影响;并且,在4种掺混方案中,基于NOx浓度的含氮天然气掺混方案的天然气掺混比例要远高于其他掺混方案,当采用此方案时,可同时满足燃烧效率、引燃温度和发热量及沃泊指数的要求。于是针对塔一联含氮天然气的掺混处理,选用基于NOx浓度的含氮天然气掺混方案。

5 总结

近年来随着塔河油田注氮规模的扩大,天然气中氮气含量也不断上升,为解决塔河油田注氮开采导致的天然气氮含量上升而造成的天然气燃烧炉燃烧过程不稳定的问题,通过建立塔河油田一号联加热炉锥形燃烧室数学模型,进一步研究塔一联天然气氮含量对加热炉燃烧过程关键参数的影响,研究表明,当氮气摩尔分数达到30.0%时,天然气将无法引燃,需对含氮天然气进行掺混处理。

针对塔一联含氮天然气组分的特殊性,掺混方案确定过程增加了氮含量对加热炉引燃温度和NOx浓度影响的考虑,分别以燃烧效率、NOx浓度、引燃温度、发热量及沃泊指数为基准对掺混方案进行计算,对比塔一联含氮天然气的4种不同基准掺混方案,确定选用基于NOx浓度的含氮天然气掺混方案,可保证塔一联含氮天然气燃烧过程的燃烧效率、引燃温度和发热量要求,同时有效地减少了氮氧化合物引发的环境污染。

参考文献
[1]
涂兴万. 碳酸盐岩缝洞型油藏单井注水替油开采的成功实践[J]. 新疆石油地质, 2008, 29(6): 735-736.
[2]
鲁新便. 岩溶缝洞型碳酸盐岩储集层的非均质性[J]. 新疆石油地质, 2003, 24(4): 360-362. DOI:10.3969/j.issn.1001-3873.2003.04.030
[3]
涂兴万, 陈朝晖. 塔河碳酸盐岩缝洞型油藏水动力学模拟新方法[J]. 西南石油学院学报, 2006, 28(5): 53-56. DOI:10.3863/j.issn.1674-5086.2006.05.015
[4]
KUSAKA J, TSUZUKI K I, DAISHO Y, et al. A numerical study on combustion and exhaust gas emissions characteristics of a dual fuel natural gas engine using a multi-dimensional model combined with detailed kinetics[C]//JSAE/SAE International Spring Fuels and Lubricants Meeting. New York: SAE, 2002: 89-93. https://www.sae.org/publications/technical-papers/content/2002-01-1750/
[5]
THOBOIS L, LAUVERGNE R, GIMBRES D, et al. The analysis of natural gas engine combustion specificities in comparison with isooctane through CFD computation[R]. SAE Technical Paper 2003-01-0009, 2003.
[6]
窦慧莉.电控喷射稀燃天然气发动机的关键技术研究[D].长春: 吉林大学, 2006. http://cdmd.cnki.com.cn/Article/CDMD-10183-2006109712.htm
[7]
赵增武, 苍大强, 李保卫, 等.烟气再循环对高温空气平焰燃烧及其NOx排放特性的影响[C]//2007年中国机械工程学会年会论文汇编.长沙: 中国机械工程学会, 2007. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGVE200711001144.htm
[8]
宋少鹏, 卓建坤, 李娜, 等. 燃料分级与烟气再循环对天然气低氮燃烧特性影响机理[J]. 中国电机工程学报, 2016, 36(24): 6849-6858.
[9]
张培理, 杜扬. 油气爆炸的氮气非预混抑制实验[J]. 爆炸与冲击, 2016, 36(3): 347-352.
[10]
SAPKO M J, WEISS E S, CASHDOLLAR K L, et al. Experimental mine and laboratory dust explosion research at NIOSH[J]. Journal of Loss Prevention in the Process Industries, 2000, 13(3/4/5): 229-242.
[11]
胡满银, 乔欢, 杜欣, 等. 烟气再循环对炉内氮氧化物生成影响的数值模拟[J]. 华北电力大学学报, 2007, 34(6): 77-82. DOI:10.3969/j.issn.1007-2691.2007.06.017
[12]
张利琴, 宋蔷, 吴宁, 等. 煤烟气再循环富氧燃烧污染物排放特性研究[J]. 中国电机工程学报, 2009, 29(29): 35-40. DOI:10.3321/j.issn:0258-8013.2009.29.007
[13]
杨伟杰. 烟气再循环燃气燃烧器的试验分析[J]. 工业锅炉, 2008(5): 23-25. DOI:10.3969/j.issn.1004-8774.2008.05.008
[14]
LAM S H, GOUSSIS D A. The CSP method for simplifying kinetics[J]. International Journal of Chemical Kinetics, 1994, 26(4): 461-486. DOI:10.1002/(ISSN)1097-4601
[15]
MASSIAS A, DIAMANTIS D, MASTORAKOS E, et al. An algorithm for the construction of global reduced mechanisms with CSP data[J]. Combustion and Flame, 1999, 117(4): 685-708. DOI:10.1016/S0010-2180(98)00132-1
[16]
SIVATHANU Y R, FAETH G M. Generalized state relationships for scalar properties in nonpremixed hydrocarbon/air flames[J]. Combustion and Flame, 1990, 82(2): 211-230. DOI:10.1016/0010-2180(90)90099-D
[17]
JONES W P, WHITELAW J H. Calculation methods for reacting turbulent flows: a review[J]. Combustion and Flame, 1982, 48: 1-26. DOI:10.1016/0010-2180(82)90112-2
[18]
朱自强. 应用计算流体力学[M]. 北京: 北京航空航天大学出版社, 1998.
[19]
MENON S, YEUNG P K, KIM W W. Effect of subgrid models on the computed interscale energy transfer in isotropic turbulence[J]. Computers and Fluids, 1996, 25(2): 165-180. DOI:10.1016/0045-7930(95)00036-4
[20]
DUKOWICZ J K. A particle-fluid numerical model for liquid sprays[J]. Journal of Computational Physics, 1980, 35(2): 229-253. DOI:10.1016/0021-9991(80)90087-X
[21]
JONES W P, MARQUIS A J, VOGIATZAKI K. Large-eddy simulation of spray combustion in a gas turbine combustor[J]. Combustion and Flame, 2014, 161(1): 222-239. DOI:10.1016/j.combustflame.2013.07.016