1.本发明属于动力工程技术领域，具体涉及一种实时计算天然气可燃物主成分比例的方法。
背景技术：
2.目前我国燃机机组燃气来源复杂，造成实际运行过程中所使用的天然气主成分比例与设计参数差异较大，天然气主成分比例不同对燃烧稳定性有很大影响，而大多数电厂仅测量有天然气热值数据，没有天然气主成分数据，因此如何利用电厂现有数据实时计算天然气主成分占比有重要的意义。
技术实现要素：
3.为了克服上述现有技术存在的不足，本发明的目的在于提供一种实时计算天然气可燃物主成分比例的方法，该方法能够计算出天然气可燃物主成分比例，为电厂运行人员优化操作、提升燃烧稳定性提供数据支持，对提高机组安全运行水平具有重要的指导意义。
4.为了实现上述目的，本发明采用的技术方案是：
5.一种实时计算天然气可燃物主成分比例的方法，包括以下步骤；
6.步骤一：燃机进出口气体实时参数采集；
7.步骤二：碳元素平衡计算；
8.燃机进口碳元素摩尔量等于空压机进口空气中co2所含碳元素摩尔量与天然气可燃物中所含碳元素摩尔量之和，燃机出口碳元素摩尔量等于燃机出口烟气中含碳元素的气体成分co2摩尔量和co摩尔量之和，分别计算上述参数并根据计算结果进行碳元素平衡计算；
9.步骤三：氢元素平衡计算；
10.燃机进口氢元素摩尔量等于天然气中烷烃部分所含氢元素摩尔量与进口空气中水分子所含氢元素之和，燃机出口氢元素摩尔量等于燃机出口烟气中水分子中氢元素含量，分别计算上述参数并根据计算结果进行氢守恒计算；
11.步骤四：能量守恒计算；
12.天然气热值由其可燃物部分贡献，依据已有天然气高位发热量与低位发热量数据，，对天然气可燃物主成分甲烷、乙烷、丙烷、丁烷进行能量守恒计算；
13.步骤五：天然气可燃物主成分体积分数计算。
14.所述步骤一具体为：
15.依托火电厂厂级监控信息系统平台，实时采集燃机进口空气和出口烟气实时数据，包括空压机进口空气温度t
air
、压力p
air
、体积流量v
air
、co2体积分数相对湿度出口烟气温度t
fg
、压力p
fg
、体积流量v
fg
、co2体积分数co体积分数vol
fg,co
、h2o体积分数以及天然气高位发热量q
hv,fuel
和低位发热量q
lv,fuel
。
16.所述步骤二具体为：
17.1)计算空压机进口空气中co2摩尔量；
18.根据气体状态方程pv＝nrt以及混合气体分体积定律，得出进口空气中co2摩尔量如下式所示：
[0019][0020]
式中，
[0021]
——空压机进口空气中co2摩尔量，kmol/h；
[0022]
p
air
——空压机进口空气压力绝对值，kpa；
[0023]
v
air
——空压机进口空气体积流量，km3/h；
[0024]
——空气中co2体积分数，％；
[0025]
r——摩尔气体常数，kj/(mol
·
k)；
[0026]
t
air
——空压机进口空气温度，k。
[0027]
2)计算烟气中co2摩尔量；
[0028][0029]
式中，
[0030]
——燃气轮机出口烟气中co2摩尔量，kmol/h；
[0031]
p
fg
——燃气轮机出口烟气压力绝对值，kpa；
[0032]
v
fg
——燃气轮机出口烟气体积流量，km3/h；
[0033]
——燃气轮机出口烟气中co2体积分数，％；
[0034]
t
fg
——燃气轮机出口烟气温度，k。
[0035]
3)计算烟气中co摩尔量；
[0036]
n
fg,co
＝f(p
fg
,v
fg
,vol
fg,co
,t
fg
)
ꢀꢀꢀꢀꢀ
(3)
[0037]
式中，
[0038]
n
fg,co
——燃气轮机出口烟气中co摩尔量，kmol/h；
[0039]
vol
fg,co
——燃气轮机出口烟气中co体积分数，％。
[0040]
4)碳元素守恒计算；
[0041]
根据天然气主成分各自燃烧过程的化学反应方程式，得出燃气轮机进出口碳元素守恒公式如下：
[0042][0043]
式中，
[0044]
——天然气中甲烷摩尔量，kmol/h；
[0045]
——天然气中乙烷摩尔量，kmol/h；
[0046]
——天然气中丙烷摩尔量，kmol/h；
[0047]
——天然气中丁烷摩尔量，kmol/h。
[0048]
所述步骤三具体为：
[0049]
1)计算进口空气中h2o摩尔量；
[0050]
根据空气状态参数经验公式可计算得到空气当前状态下的饱和水蒸气压力p
q,b
，kpa；
[0051]
p
q,b
＝f(t
air
)
ꢀꢀꢀꢀꢀ
(5)
[0052]
利用计算出的空气饱和水蒸气压力与相对湿度数据可计算出当前状态下的水蒸气压力p
q
如下：
[0053][0054]
式中：
[0055]
p
q
——进口空气水蒸气压力，kpa；
[0056]
——进口空气相对湿度，％；
[0057]
p
q,b
——进口空气饱和水蒸气压力，kpa；
[0058]
由道尔顿分压定律及气体状态方程可推导出进口空气中水分子摩尔量计算如下：
[0059][0060]
式中：
[0061]
——进口空气中h2o摩尔量，kmol/h。
[0062]
2)计算烟气中h2o摩尔量；
[0063][0064]
式中，
[0065]
——燃气轮机出口烟气中h2o摩尔量，kmol/h；
[0066]
——燃气轮机出口烟气中h2o体积分数，％。
[0067]
3)氢元素守恒计算；
[0068][0069]
所述步骤四具体为：
[0070]
根据气体状态方程pv＝nrt，当压力和温度相同时，体积分数等于摩尔量分数，则天然气总摩尔量可表示为：
[0071][0072]
式中，
[0073]
λ
nc
——不可燃物体积分数，％。
[0074]
天然气热值是由可燃气体热值计算得到，其中丁烷按照正丁烷分析，不考虑其同分异构体，据此进行能量守恒计算；
[0075]
1)低位发热量守恒计算；
[0076][0077]
式中，
[0078]
——甲烷低位发热量，kj/mol；
[0079]
——乙烷低位发热量，kj/mol；
[0080]
——丙烷低位发热量，kj/mol；
[0081]
——丁烷低位发热量，kj/mol；
[0082]
q
lv,fuel
——天然气低位发热量，kj/mol。
[0083]
2)高位发热量守恒计算；
[0084][0085]
式中，
[0086]
——甲烷高位发热量，kj/mol；
[0087]
——乙烷高位发热量，kj/mol；
[0088]
——丙烷高位发热量，kj/mol；
[0089]
——丁烷高位发热量，kj/mol；
[0090]
q
hv,fuel
——天然气高位发热量，kj/mol。
[0091]
所述步骤五具体为：
[0092]
通过前四个步骤得到天然气中可燃物主成分各自的摩尔量通过前四个步骤得到天然气中可燃物主成分各自的摩尔量以及天然气总摩尔量n
fule
；同上，根据气体状态方程pv＝nrt，当压力和温度相同时，体积分数等于摩尔量分数，则可燃物主成分体积分数可表示为：
[0093][0094][0095][0096][0097]
本发明的有益效果：
[0098]
本发明利用电厂厂级监控信息系统(sis)所采集的数据，对燃机进口天然气、空气以及燃机出口烟气中成分分别进行元素守恒计算与能量守恒计算，可实时计算出天然气可燃物主成分中甲烷(ch4)、乙烷(c2h6)、丙烷(c3h8)、丁烷(c4h
10
)的体积分数，方便机组运行人员实时了解天然气主成分参数信息，为提升燃机机组燃烧稳定性提供数据支撑。
具体实施方式
[0099]
下面结合实施例对本发明作进一步详细说明。
[0100]
燃机进出口气体实时参数采集；
[0101]
依托火电厂厂级监控信息系统平台，实时采集燃机进口空气和出口烟气实时数
据，包括空压机进口空气温度t
air
、压力p
air
、体积流量v
air
、co2体积分数相对湿度出口烟气温度t
fg
、压力p
fg
、体积流量v
fg
、co2体积分数co体积分数vol
fg,co
、h2o体积分数以及天然气高位发热量q
hv,fuel
和低位发热量q
lv,fuel
。
[0102]
碳元素平衡计算；
[0103]
燃机进口碳元素来源分别为空压机进口空气中co2所含碳元素及天然气中所含碳元素两部分，燃机出口烟气中含碳元素的气体成分为co2和co，据此进行碳元素守恒计算。
[0104]
1)计算空压机进口空气中co2摩尔量。
[0105]
根据气体状态方程pv＝nrt以及混合气体分体积定律，可以推导出进口空气中co2摩尔量如下式所示：
[0106][0107]
式中，
[0108]
——空压机进口空气中co2摩尔量，kmol/h；
[0109]
p
air
——空压机进口空气压力绝对值，kpa；
[0110]
v
air
——空压机进口空气体积流量，km3/h；
[0111]
——空气中co2体积分数，％；
[0112]
r——摩尔气体常数，kj/(mol
·
k)；
[0113]
t
air
——空压机进口空气温度，k。
[0114]
2)计算烟气中co2摩尔量。
[0115][0116]
式中，
[0117]
——燃气轮机出口烟气中co2摩尔量，kmol/h；
[0118]
p
fg
——燃气轮机出口烟气压力绝对值，kpa；
[0119]
v
fg
——燃气轮机出口烟气体积流量，km3/h；
[0120]
——燃气轮机出口烟气中co2体积分数，％；
[0121]
t
fg
——燃气轮机出口烟气温度，k。
[0122]
3)计算烟气中co摩尔量
[0123]
n
fg,co
＝f(p
fg
,v
fg
,vol
fg,co
,t
fg
)
ꢀꢀꢀꢀꢀ
(3)
[0124]
式中，
[0125]
n
fg,co
——燃气轮机出口烟气中co摩尔量，kmol/h；
[0126]
vol
fg,co
——燃气轮机出口烟气中co体积分数，％。
[0127]
5)碳元素守恒计算
[0128]
根据天然气主成分各自燃烧过程的化学反应方程式，可以得出燃气轮机进出口碳元素守恒公式如下：
[0129][0130]
式中，
[0131]
——天然气中甲烷摩尔量，kmol/h；
[0132]
——天然气中乙烷摩尔量，kmol/h；
[0133]
——天然气中丙烷摩尔量，kmol/h；
[0134]
——天然气中丁烷摩尔量，kmol/h。
[0135]
氢元素摩尔量计算；
[0136]
燃机进口氢元素摩尔量等于天然气中烷烃部分所含氢元素摩尔量与进口空气中水分子所含氢元素之和，燃机出口烟气中含氢元素的气体成分为h2o，据此进行氢元素守恒计算。
[0137]
1)计算进口空气中h2o摩尔量；
[0138]
根据空气状态参数经验公式可计算得到空气当前状态下的饱和水蒸气压力p
q,b
，kpa；
[0139]
p
q,b
＝f(t
air
)
ꢀꢀꢀꢀꢀ
(5)
[0140]
利用计算出的空气饱和水蒸气压力与相对湿度数据可计算出当前状态下的水蒸气压力p
q
如下：
[0141][0142]
式中：
[0143]
p
q
——进口空气水蒸气压力，kpa；
[0144]
——进口空气相对湿度，％；
[0145]
p
q,b
——进口空气饱和水蒸气压力，kpa。
[0146]
由道尔顿分压定律及气体状态方程可推导出进口空气中水分子摩尔量计算如下：
[0147][0148]
式中：
[0149]
——进口空气中h2o摩尔量，kmol/h。
[0150]
2)计算烟气中h2o摩尔量
[0151][0152]
式中，
[0153]
——燃气轮机出口烟气中h2o摩尔量，kmol/h；
[0154]
——燃气轮机出口烟气中h2o体积分数，％。
[0155]
3)氢元素守恒计算
[0156][0157]
能量守恒计算；
[0158]
天然气成分可分为可燃物部分和不可燃物部分，其中不可燃物体积分数λ
nc
可由运行人员根据现场实际进行设置。
[0159]
根据气体状态方程pv＝nrt，当压力和温度相同时，体积分数等于摩尔量分数，则天然气总摩尔量可表示为：
[0160][0161]
天然气热值是由可燃气体热值计算得到，其中丁烷按照正丁烷分析，不考虑其同分异构体，据此进行能量守恒计算。
[0162]
1)低位发热量守恒计算；
[0163][0164]
式中，
[0165]
——甲烷低位发热量，kj/mol；
[0166]
——乙烷低位发热量，kj/mol；
[0167]
——丙烷低位发热量，kj/mol；
[0168]
——丁烷低位发热量，kj/mol；
[0169]
q
lv,fuel
——天然气低位发热量，kj/mol。
[0170]
2)高位发热量守恒计算；
[0171][0172]
式中，
[0173]
——甲烷高位发热量，kj/mol；
[0174]
——乙烷高位发热量，kj/mol；
[0175]
——丙烷高位发热量，kj/mol；
[0176]
——丁烷高位发热量，kj/mol；
[0177]
q
hv,fuel
——天然气高位发热量，kj/mol。
[0178]
天然气可燃物主成分体积分数计算；
[0179]
联立式(4)、(9)、(10)、(11)、(12)可计算得到天然气中可燃物主成分各自的摩尔量以及天然气总摩尔量n
fule
。同上，根据气体状态方程pv＝nrt，当压力和温度相同时，体积分数等于摩尔量分数，则可燃物主成分体积分数可表示为：
[0180][0181][0182][0183]