一种基于复杂工况的燃气内燃机发电功率及排烟参数的计算
方法
技术领域
[0001] 本
发明涉及一种
热能工程的
旋转机械领域，尤其涉及一种燃气内燃机在复杂工况运行条件下的发电功率及排烟参数的计算方法。
背景技术
[0002] 楼宇式
天然气分布式
能源系统是靠近用户端的能源综合利用系统。它以天然气为
燃料，通过燃气内燃机发电，并将其尾部烟气通过
能量转换设备生产出用户所需要的热能利用形式，如生活热
水、采暖热水或制冷用热能等；按照“
温度对口、
梯级利用”的用能原则，使得整个系统的综合能源利用率大于70％。发展楼宇式天然气分布式能源是降低能源成本、提升能源效率、改善大气环境的一种有效技术途径。
[0003] 燃气内燃机是构建天然气分布式能源系统的主要设备之一，其发电功率、燃机效率、燃气消耗量、排气温度、排气流量、高温
循环水流量、低温循环水流量等参数受到其运行条件下的大气环境的温度、压
力、湿度以及海拔的影响很大。在某些工况条件下，燃气内燃机的发电功率、排气温度及排气流量、等参数很难与系统需求相匹配。确定复杂工况条件下，燃气内燃机的发电功率、排烟参数、高温循环水流量、低温循环水流量对天然气分布式能源系统的方案配置起到了决定性作用。
[0004] 因此，要构建一个适用于复杂工况运行条件下的燃气内燃机发电功率及排烟参数的计算方法，是合理配置楼宇式天然气分布式能源系统的
基础手段，能够为楼宇式天然气分布式能源系统运行方案提供基础数据，具有重要的实用意义。
发明内容
[0005] 针对上述问题，本发明提供一种适用于复杂工况运行条件下的燃气内燃机发电功率及排烟参数的计算方法。
[0006] 为达到上述目的，本发明一种适用于复杂工况运行条件下的燃气内燃机发电功率及排烟参数的计算方法，所述的发电功率和排烟参数的计算方法包括：
[0007] S1，相关参数的采集与测定；
[0008] S2，压缩过程热力计算；
[0010] S4，膨胀过程热力计算；
[0011] S5，燃气内燃机发电功率及排烟参数的计算。
[0012] 进一步地，相关参数的采集与测定具体包括如下步骤：
[0013] S11，根据燃气内燃机项目所在的地点，采集燃气内燃机周围环境的
大气压力、温度、湿度、海拔高度等参数信息；
[0014] S12，根据所
选定的燃气内燃机机型，收集燃气内燃机标准工况下的运行参数，包括：内燃机
气缸数、缸径、行程、标准状态下的发电功率、标准状态下的发电效率、标准状态下燃气内燃机的进燃气量、压缩比、供气压力、确定燃气内燃机的进气方式、油冷器在高低温循环系统的
位置、内燃机排烟温度、排烟流量、高温水供回水温度及流量、中冷水供回水温度及流量；
[0015] S13，根据燃气内燃机项目所在地，收集燃气内燃机所用天然气成分，通过天然气热力计算公式计算得到项目所使用的天然气的低位热值、高位热值、
密度等参数，其具体步骤包括：
[0016] 1)根据第一计算式计算天然气的低位热值，所述第一计算式为：
[0017]
[0018] 式中，QNG，L为天然气低位发热量，MJ/Nm3；mi为天然气各组分的摩尔体积百分数，％；Qi，L为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kJ/mol；PNG为天然气的压力，kPa；TNG为天然气的温度，K；Z为天然气压缩系数；
[0019] 2)根据第二计算式计算天然气的高位热值，所述第二计算式为：
[0020]
[0021] 式中，QNG，H为天然气高位发热量，MJ/Nm3；Qi，H为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kJ/mol；
[0022] 3)根据第三计算式计算天然气的密度，所述第三计算式为：
[0023] ρNG＝(PNG/(R·TNG))×∑mi·Mi
[0024] 式中，ρNG为天然气密度，kg/Nm3；Mi为天然气各组分的摩尔
质量，kg/kmol。
[0025] S14，根据S11和S12所采集的数据，通过湿空气的气体常数计算公式、湿空气的
比热及比热比计算公式等计算大气湿度对空气热物理性质的影响，进而利用相似原理，推导出燃气内燃机性能参数大气湿度修正系数。
[0026] 4)根据第四计算式计算湿空气气体常数，所述第四计算式为：
[0027]
[0028] 式中，Rh为湿空气气体常数；Rk为干空气气体常数；Rs为水
蒸汽气体常数；d为湿空气中相对于每千克质量干空气所含的水蒸汽量；
[0029] 5)根据第五计算式计算湿空气定压比热，所述第五计算式为：
[0030]
[0031] 式中，Cph为湿空气定压比热，J/(kg·K)；Cpk为干空气定压比热，J/(kg·K)；Cps为水蒸汽定压比热，J/(kg·K)；
[0032] 6)根据第六计算式计算湿空气比热比，所述第六计算式为：
[0033]
[0034] 式中，Kh为湿空气比热比；
[0035] 7)根据第七计算式计算干空气比热比，所述第七计算式为：
[0036]
[0037] 式中，Kk为干空气比热比；
[0038] 8)根据第八计算式计算大气湿度变化时，空气流量的修正系数，所述第八计算式为：
[0039]
[0040] 式中，γ为大气湿度变化时，空气流量的修正系数。
[0041] 进一步地，压缩过程热力计算还包括如下步骤：
[0042] 1)根据第九计算式计算实际运行工况下燃气内燃机的进燃气量，所述第九计算式为：
[0043]
[0044] 式中，Mf为实际运行工况下燃气内燃机的进燃气量，kg/s；Mfs为标准状态下燃气内燃机的进燃气量，kg/s；Pa为实际运行工况下的大气压力，MPa；Ta为实际运行工况下的大气温度，℃；Ps为标准状态下的大气压力，MPa；Ts为标准状态下的大气温度，℃；
[0045] 2)根据第十计算式计算天然气在燃气内燃机中燃烧所需要的空气量，所述第十计算式为：
[0046]
[0047] 式中，Vair天然气在燃气内燃机中燃烧所需要的空气量，m3/s；CnHm为天然气成分，％；n为天然气成分的
碳原子个数；m为天然气成分的氢原子个数； 为燃气内燃机运行工况下的环境空气中
氧气的百分含量；αair为空气过量系数；
[0048] 3)根据第十一计算式计算燃气内燃机的
空燃比，所述第十一计算式为：
[0049]
[0050] 式中，f为燃气内燃机的空燃比，m3/m3；
[0051] 4)根据第十二计算式计算燃气内燃机混合气体各成分容积数，所述第十二计算式为：
[0052]
[0053] 式中，Vi为燃气内燃机混合气体各成分容积数，m3/s；ii为燃气内燃机燃烧的天然气的各成分的体积百分含量，％；
[0054] 5)根据第十三计算式计算燃气内燃机混合气体总容积数，所述第十三计算式为：
[0055] V7＝∑Vi
[0056] 式中，V7为燃气内燃机混合气体总容积数，m3/s；
[0057] 6)根据第十四计算式计算燃气内燃机混合气体各成分容积的百分比，所述第十四计算式为：
[0058]
[0059] 式中，γi为燃气内燃机混合气体各成分容积的百分比，％；
[0060] 7)根据第十五计算式计算燃气内燃机混合气体体积热值，所述第十五计算式为：
[0061]
[0062] 式中，Qmix，L为燃气内燃机混合气体体积热值，kJ/Nm3；QNG，L，mol为燃气内燃机混合气体的低位摩尔发热量，kJ/mol；Pmix为燃气内燃机混合气体压力，MPa；tmix为燃气内燃机混合气体温度，℃；Zmix为燃气内燃机混合气体的压缩因子；R为燃气内燃机混合气体的气体常数，8.31451；
[0063] 8)根据第十六计算式计算燃气内燃机混合气体密度，所述第十六计算式为：
[0064]
[0065] 式中，Mi，mol为燃气内燃机混合气体的摩尔质量，kg/kmol；
[0066] 9)根据第十七计算式计算燃气内燃机混合气体完全燃烧后排气中氮气的容积数，所述第十七计算式为：
[0067]
[0068] 式中， 为燃气内燃机混合气体完全燃烧后排气中氮气的容积数，m3/s； 为空气中氮气的百分含量，％； 为天然气中氮气的百分含量，％；
[0069] 10)根据第十八计算式计算燃气内燃机混合气体完全燃烧后排气中二氧化碳的容积数，所述第十八计算式为：
[0070]
[0071] 式中， 为燃气内燃机混合气体完全燃烧后排气中二氧化碳的容积数，m3/s；为空气中氮气的百分含量，％；CnHm为天然气中各碳氢化合物的百分含量，％；n为天然气成分的碳原子个数；
[0072] 11)根据第十九计算式计算燃气内燃机混合气体完全燃烧后排气中水蒸汽的容积数，所述第十九计算式为：
[0073]
[0074] 式中， 为燃气内燃机混合气体完全燃烧后排气中水蒸汽的容积数，m3/s；为空气中氮气的百分含量，％；m为天然气成分的氢原子个数；
[0075] 12)根据第二十计算式计算燃气内燃机混合气体完全燃烧后排气中氧气的容积数，所述第二十计算式为：
[0076]
[0077] 式中， 为燃气内燃机混合气体完全燃烧后排气中氧气的容积数，m3/s； 为空气中氧气的百分含量，％； 为天然气成分中氧气的百分含量，％；
[0078] 13)根据第二十一计算式计算燃气内燃机混合气体完全燃烧后排气中氩气的容积数，所述第二十一计算式为：
[0079]
[0080] 式中，VAr为燃气内燃机混合气体完全燃烧后排气中氩气的容积数，m3/s； 为空气中氩气的百分含量，％；τAr，f为天然气成分中氩气的百分含量，％；
[0081] 14)根据第二十二计算式计算燃气内燃机混合气体完全燃烧后排气中二氧化硫的容积数，所述第二十二计算式为：
[0082]
[0083] 式中， 为燃气内燃机混合气体完全燃烧后排气中二氧化硫的容积数，m3/s；为天然气成分中二氧化硫的百分含量，％； 为天然气成分中
硫化氢的百分含
量，％；
[0084] 15)根据第二十三计算式计算燃气内燃机混合气体完全燃烧后排气的总容积数，所述第二十三计算式为：
[0085]
[0086] 式中，Vpy为燃气内燃机混合气体完全燃烧后排气的总容积数，m3/s；
[0087] 16)根据第二十四计算式计算燃气内燃机混合气体完全燃烧后排气各成分的体积百分数，所述第二十四计算式为：
[0088]
[0089] 式中，γi-py为燃气内燃机混合气体完全燃烧后排气各成分的体积百分数，％；Vi为燃气内燃机混合气体完全燃烧后排气各成分的容积数，m3/s；i代表N2、CO2、H2O、O2、Ar、SO2；
[0090] 17)根据第二十五计算式计算燃气内燃机进气终点压力，所述第二十五计算式为：
[0091]
[0092] 式中，Pin-end为燃气内燃机进气终点压力，MPa；υin为燃气内燃机进气速度，m/s；
[0093] 18)根据第二十六计算式计算燃气内燃机排气终点压力，所述第二十六计算式为：
[0094]
[0095] 式中，Pout-end为燃气内燃机排气终点压力，MPa；υout为燃气内燃机排气速度，m/s；
[0096] 19)根据第二十七计算式计算燃气内燃机剩余废气系数，所述第二十七计算式为：
[0097]
[0098] 式中，γreg为燃气内燃机剩余废气系数，无量纲；γc为燃气内燃机
净化系数，无量纲；ε为内燃机压缩比，无量纲；Pr为燃气内燃机膨胀冲程终点理论压力，MPa；m为气体膨胀指数，无量纲；
[0099] 20)根据第二十八计算式计算燃气内燃机剩余废气量，所述第二十八计算式为：
[0100] Vreg＝γreg×V7
[0101] 式中，Vreg为燃气内燃机剩余废气量，m3/s；
[0102] 21)根据第二十九计算式计算燃气内燃机混合气体量，所述第二十九计算式为：
[0103] Vmix＝γreg+V7
[0104] 式中，Vmix为燃气内燃机混合气体量，m3/s；
[0105] 22)根据第三十计算式计算燃气内燃机混合气体量，所述第三十计算式为：
[0106] Vp＝γreg+Vpy
[0107] 式中，Vp为燃气内燃机排气气体量，m3/s；
[0108] 23)根据第三十一计算式计算燃气内燃机充气终点压力，所述第三十一计算式为：
[0109]
[0110] 式中，Pa为燃气内燃机充气终点压力，MPa：ξ1为补充充气系数；
[0111] 24)根据第三十二计算式计算燃气内燃机充气终点温度，所述第三十二计算式为：
[0112]
[0113] 式中，Ta为燃气内燃机充气终点温度，K；ξr为比热不相等系数；n为多变指数；
[0114] 25)根据第三十三计算式计算燃气内燃机混合压力下的剩余废气温度，所述第三十三计算式为：
[0115]
[0116] 式中，Tr-mix为燃气内燃机混合压力下的剩余废气温度，K；
[0117] 26)根据第三十四计算式计算燃气内燃机排气压力下的废气温度，所述第三十四计算式为：
[0118]
[0119] 式中，Tr-out为燃气内燃机排气压力下的废气温度，K；
[0120] 27)根据第三十五计算式计算燃气内燃机充气系数，所述第三十五计算式为：
[0121]
[0122] 式中，ξc为燃气内燃机充气系数；tair为
环境温度，℃；
[0123] 28)根据第三十六计算式计算燃气内燃机压缩终点温度，所述第三十六计算式为：
[0124]
[0125] 式中，Tc为燃气内燃机压缩终点温度，K；n1为燃气内燃机压缩多变指数；
[0126] 29)根据第三十七计算式计算燃气内燃机压缩终点压力，所述第三十七计算式为：
[0127]
[0128] 式中，Pc为燃气内燃机压缩终点压力，MPa；
[0129] 30)根据第三十八计算式计算燃气内燃机压缩多变指数，所述第三十八计算式为：
[0130] n1＝1.392-0.13×(Tc-aver-273.15)/1000
[0131] 式中，n1为燃气内燃机压缩多变指数；Tc-aver为燃气内燃机压缩终点平均温度，K；
[0132] 进一步地，燃烧热力计算还包括如下步骤：
[0133] 1)根据第三十九计算式计算燃气内燃机最高燃烧温度，所述第三十九计算式为：
[0134]
[0135] 式中，tZ为燃气内燃机最高燃烧温度，℃；hrs为燃烧产物
焓值，kJ/(kg·℃)；
[0136] 2)根据第四十计算式计算燃气内燃机压力升高比，所述第四十计算式为：
[0137]
[0138] 式中，π为燃气内燃机压力升高比；
[0139] 3)根据第四十一计算式计算燃气内燃机最高燃烧压力，所述第四十一计算式为：
[0140] PZ＝Pc×π×κ
[0141] 式中，PZ为燃气内燃机最高燃烧压力，MPa；κ为最高燃烧压力的修正系数；
[0142] 进一步地，膨胀过程热力计算还包括如下步骤：
[0143] 1)根据第四十二计算式计算燃气内燃机膨胀终点温度，所述第四十二计算式为：
[0142]
[0143] 式中，Tb为燃气内燃机膨胀终点温度，K；n2为燃气内燃机膨胀多变指数；
[0144] 2)根据第四十三计算式计算燃气内燃机膨胀终点压力，所述第四十三计算式为：
[0145]
[0146] 式中，Pb为燃气内燃机膨胀终点压力，MPa：
[0147] 3)根据第四十四计算式计算燃气内燃机膨胀多变指数，所述第四十四计算式为：
[0148]
[0149] 式中，n2为燃气内燃机膨胀多变指数；Tb-aver为燃气内燃机膨胀终点平均温度，K；
[0150] 进一步地，燃气内燃机发电功率及排烟参数的计算还包括如下步骤：
[0151] 1)根据第四十五计算式计算燃气内燃机处初膨胀比，所述第四十五计算式为：
[0152]
[0153] 式中，ρ0为燃气内燃机处初膨胀比；
[0154] 2)根据第四十六计算式计算燃气内燃机后膨胀比，所述第四十六计算式为：
[0155]
[0156] 式中，δ为燃气内燃机处后膨胀比；
[0157] 3)根据第四十七计算式计算燃气内燃机平均指示压力，所述第四十七计算式为：
[0158]
[0159] 式中，Pi为燃气内燃机平均指示压力，MPa；
[0160] 4)根据第四十八计算式计算燃气内燃机平均有效压力，所述第四十八计算式为：
[0161] Pmi＝Pi×ηm
[0162] 式中，Pmi为燃气内燃机平均有效压力，MPa；ηm为燃气内燃机机械效率，％；
[0163] 5)根据第四十九计算式计算燃气内燃机有效功率，所述第四十九计算式为：
[0164]
[0165] 式中，Pe为燃气内燃机有效功率，kW；nc为燃气内燃机气缸数；Vw为燃气内燃机的工作容积，m3；ω为燃气内燃机转速，r/min；
[0166] 6)根据第五十计算式计算燃气内燃机发电功率，所述第五十计算式为：
[0167] Pele＝Pe×ηele
[0168] 式中，Pele为燃气内燃机发电功率，kW；ηele为燃气内燃机发
电机效率，％；
[0169] 7)根据第五十一计算式计算燃气内燃机发电效率，所述第五十一计算式为：
[0170]
[0171] 式中，η为燃气内燃机发电效率，％；
[0172] 8)根据第五十二计算式计算燃气内燃机高温循环水在额定工况下带走的热量，所述第五十二计算式为：
[0173] QHTW-in design＝Cpw×MHTW-in design×(tHTW-out-tHTW-in)/3.6
[0174] 式中，QHTW-in design为燃气内燃机高温循环水在额定工况下带走的热量，kJ；Cpw为水的比热，kJ/(kg·℃)；MHTW-in design为燃气内燃机高温循环水在额定工况下的流量，kg/h；tHTW-out为燃气内燃机高温循环水在额定工况下出口处的温度，℃；tHTW-in为燃气内燃机高温循环水在额定工况下进口处的温度，℃；
[0175] 9)根据第五十三计算式计算燃气内燃机高温循环水在额定工况下带走的热量占总热量的百分比，所述第五十三计算式为：
[0176]
[0177] 式中，ηHTM为燃气内燃机高温循环水在额定工况下带走的热量占总热量的百分比，％；Pele-in design为在额定工况下燃气内燃机发电功率，kW；ηele-in design为在额定工况下燃气内燃机发电效率，％；
[0178] 10)根据第五十四计算式计算燃气内燃机高温循环水在运行工况下带走的热量，所述第五十四计算式为：
[0179] QHTW-off design＝Mf×QNG，L×1000×3600×ηHTM
[0180] 式中，QHTW-off design为燃气内燃机高温循环水在运行工况下带走的热量，kJ；
[0181] 11)根据第五十五计算式计算燃气内燃机高温循环水在运行工况下的流量，所述第五十五计算式为：
[0182]
[0183] 式中，MHTW-off design为燃气内燃机高温循环水在运行工况下的流量，kg/s；t′HTW-out为燃气内燃机高温循环水在运行工况下出口处的温度，℃；t′HTW-in为燃气内燃机高温循环水在运行工况下进口处的温度，℃；
[0184] 12)根据第五十六计算式计算燃气内燃机低温循环水在额定工况下带走的热量，所述第五十六计算式为：
[0185] QLTW-in design＝Cpw×MLTW-in design×(tLTW-out-tLTW-in)/3.6
[0186] 式中，QLTW-in design为燃气内燃机低温循环水在额定工况下带走的热量，kJ；MLTW-in design为燃气内燃机低温循环水在额定工况下的流量，kg/h；tLTW-out为燃气内燃机低温循环水在额定工况下出口处的温度，℃；tLTW-in为燃气内燃机低温循环水在额定工况下进口处的温度，℃；
[0187] 13)根据第五十七计算式计算燃气内燃机低温循环水在额定工况下带走的热量占总热量的百分比，所述第五十七计算式为：
[0188]
[0189] 式中，ηLTM为燃气内燃机低温循环水在额定工况下带走的热量占总热量的百分比，％；
[0190] 14)根据第五十八计算式计算燃气内燃机低温循环水在运行工况下带走的热量，所述第五十八计算式为：
[0191] QLTW-off design＝Mf×QNG，L×1000×3600×ηLTM
[0192] 式中，QLTW-off design为燃气内燃机低温循环水在运行工况下带走的热量，kJ；
[0193] 15)根据第五十九计算式计算燃气内燃机高温循环水在运行工况下的流量，所述第五十九计算式为：
[0194]
[0195] 式中，MLTW-off design为燃气内燃机低温循环水在运行工况下的流量，kg/s；t′LTW-out为燃气内燃机低温循环水在运行工况下出口处的温度，℃；t′LTW-in为燃气内燃机低温循环水在运行工况下进口处的温度，℃；
[0196] 16)根据第六十计算式计算燃气内燃机烟气在额定工况下带走的热量，所述第六十计算式为：
[0197] Qg-indesign＝Mg×hg-indesign
[0198] 式中，Qg-indesign为燃气内燃机烟气在额定工况下带走的热量，kJ；Mg为燃气内燃机烟气在额定工况下的流量，m3/s；hg-indesign为燃气内燃机烟气在额定工况下的焓值，kJ/m3；
[0199] 17)根据第六十一计算式计算燃气内燃机烟气在额定工况下带走的热量占总热量的百分比，所述第六十一计算式为：
[0200]
[0201] 式中，ηg为燃气内燃机烟气在额定工况下带走的热量占总热量的百分比，％；
[0202] 18)根据第六十二计算式计算燃气内燃机在额定工况下
散热损失占总热量的百分比，所述第六十二计算式为：
[0203] ηloss＝1-η-ηHTM-ηLTM-ηg
[0204] 式中，ηloss为燃气内燃机在额定工况下散
热损失占总热量的百分比，％；
[0205] 19)根据第六十三计算式计算燃气内燃机在运行工况下散热量，所述第六十三计算式为：
[0206] Qloss-off design＝Mf×QNG，L×1000×3600×ηloss
[0207] 式中，Qloss-off design为燃气内燃机在运行工况下散热量，kJ；
[0208] 20)根据第六十四计算式计算燃气内燃机在运行工况下排烟所携带的热量，所述第六十四计算式为：
[0209] Qg-off design＝Mf×QNG，L×1000×3600-Pele×3600-QHTW-off design-QLTW-off design-Qloss-off design
[0210] 式中，Qg-off design为燃气内燃机在运行工况下排烟所携带的热量，kJ；
[0211] 21)根据第六十五计算式计算燃气内燃机在运行工况下排烟流量，所述第六十五计算式为：
[0212]
[0213] 式中，V7′为核算的燃气内燃机在运行工况下排烟流量，m3/s；hg-offdesign为燃气内燃机烟气在运行工况下的焓值，kJ/m3；
[0214] 本发明一种适用于复杂工况运行条件下的燃气内燃机发电功率及排烟参数的计算方法，通过计算得到复杂工况运行条件下的燃气内燃机的发电功率及排烟参数，具有以下优点：
[0215] 1、克服了
现有技术条件下无法准确计算燃气内燃机所处地理位置的大气环境参数，从而使得变工况运行条件下的燃气内燃机发电功率及排烟参数的计算得以顺利进行；
[0216] 2、本发明的燃气内燃机发电功率及排烟参数的计算方法，能够通过计算结果反映燃气内燃机的运行情况，进而为楼宇式天然气分布式能源系统的配置方案及运行方式、优化运行提供指导，达到节能、降低损耗的目的。
附图简要说明
[0217] 图1是本发明一种适用于复杂工况运行条件下的燃气内燃机发电功率及排烟参数的计算方法示意图。
具体实施方式
[0218] 下面结合
说明书附图对本发明做进一步的描述。
[0219] 如图1所示为一种适用于复杂工况运行条件下的燃气内燃机发电功率及排烟参数的计算方法，其过程包含了S1相关参数的采集与测定，S2压缩过程热力计算，S3燃烧热力计算，S4膨胀过程热力计算，S5燃气内燃机发电功率及排烟参数的计算。所述计算具体如下：
[0220] 首先进行相关参数的采集与测定，通过仪表测量和取样分析获取
燃气轮机的各项输入参数，具体包括：
[0221] 根据燃气内燃机项目所在的地点，采集燃气内燃机周围环境的大气压力、温度、湿度、海拔高度等参数信息；
[0222] 根据所选定的燃气内燃机机型，收集燃气内燃机标准工况下的运行参数，包括：内燃机气缸数、缸径、行程、标准状态下的发电功率、标准状态下的发电效率、标准状态下燃气内燃机的进燃气量、压缩比、供气压力、确定燃气内燃机的进气方式、油冷器在高低温循环系统的位置、内燃机排烟温度、排烟流量、高温水供回水温度及流量、中冷水供回水温度及流量；
[0223] 根据燃气内燃机项目所在地，收集燃气内燃机所用天然气成分，通过天然气热力计算公式计算得到项目所使用的天然气的低位热值、高位热值、密度等参数，其具体步骤包括：
[0224] 根据第一计算式计算天然气的低位热值，所述第一计算式为：
[0225]
[0226] 式中，QNG，L为天然气低位发热量，MJ/Nm3；mi为天然气各组分的摩尔体积百分数，％；Qi，L为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kJ/mol；PNG为天然气的压力，kPa；TNG为天然气的温度，K；Z为天然气压缩系数；
[0227] 根据第二计算式计算天然气的高位热值，所述第二计算式为：
[0228]
[0229] 式中，QNG，H为天然气高位发热量，MJ/Nm3；Qi，H为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kJ/mol；
[0230] 根据第三计算式计算天然气的密度，所述第三计算式为：
[0231] ρNG＝(PNG/(R·TNG))×∑mi·Mi
[0232] 式中，ρNG为天然气密度，kg/Nm3；Mi为天然气各组分的摩尔质量，kg/kmol。
[0233] 根据所采集的数据，通过湿空气的气体常数计算公式、湿空气的比热及比热比计算公式等计算大气湿度对空气热物理性质的影响，进而利用相似原理，推导出燃气内燃机性能参数大气湿度修正系数。
[0234] 根据第四计算式计算湿空气气体常数，所述第四计算式为：
[0235]
[0236] 式中，Rh为湿空气气体常数；Rk为干空气气体常数；Rs为水蒸汽气体常数；d为湿空气中相对于每千克质量干空气所含的水蒸汽量；
[0237] 根据第五计算式计算湿空气定压比热，所述第五计算式为：
[0238]
[0239] 式中，Cph为湿空气定压比热，J/(kg·K)；Cpk为干空气定压比热，J/(kg·K)；Cps为水蒸汽定压比热，J/(kg·K)；
[0240] 根据第六计算式计算湿空气比热比，所述第六计算式为：
[0241]
[0242] 式中，Kh为湿空气比热比；
[0243] 根据第七计算式计算干空气比热比，所述第七计算式为：
[0244]
[0245] 式中，Kk为干空气比热比；
[0246] 根据第八计算式计算大气湿度变化时，空气流量的修正系数，所述第八计算式为：
[0247]
[0248] 式中，γ为大气湿度变化时，空气流量的修正系数。
[0249] 根据第九计算式计算实际运行工况下燃气内燃机的进燃气量，所述第九计算式为：
[0250]
[0251] 式中，Mf为实际运行工况下燃气内燃机的进燃气量，kg/s；Mfs为标准状态下燃气内燃机的进燃气量，kg/s；Pa为实际运行工况下的大气压力，MPa；Ta为实际运行工况下的大气温度，℃；Ps为标准状态下的大气压力，MPa；Ts为标准状态下的大气温度，℃；
[0252] 根据第十计算式计算天然气在燃气内燃机中燃烧所需要的空气量，所述第十计算式为：
[0253]
[0254] 式中，Vair天然气在燃气内燃机中燃烧所需要的空气量，m3/s；CnHm为天然气成分，％；n为天然气成分的碳原子个数；m为天然气成分的氢原子个数； 为燃气内燃机运行工况下的环境空气中氧气的百分含量；αair为空气过量系数；
[0255] 根据第十一计算式计算燃气内燃机的空燃比，所述第十一计算式为：
[0256]
[0257] 式中，f为燃气内燃机的空燃比，m3/m3；
[0258] 根据第十二计算式计算燃气内燃机混合气体各成分容积数，所述第十二计算式为：
[0259]
[0260] 式中，Vi为燃气内燃机混合气体各成分容积数，m3/s；ii为燃气内燃机燃烧的天然气的各成分的体积百分含量，％；
[0261] 根据第十三计算式计算燃气内燃机混合气体总容积数，所述第十三计算式为：
[0262] V7＝∑Vi
[0263] 式中，V7为燃气内燃机混合气体总容积数，m3/s；
[0264] 根据第十四计算式计算燃气内燃机混合气体各成分容积的百分比，所述第十四计算式为：
[0265]
[0266] 式中，γi为燃气内燃机混合气体各成分容积的百分比，％；
[0267] 根据第十五计算式计算燃气内燃机混合气体体积热值，所述第十五计算式为：
[0268]
[0269] 式中，Qmix，L为燃气内燃机混合气体体积热值，kJ/Nm3；QNG，L，mol为燃气内燃机混合气体的低位摩尔发热量，kJ/mol；Pmix为燃气内燃机混合气体压力，MPa；tmix为燃气内燃机混合气体温度，℃；Zmix为燃气内燃机混合气体的压缩因子；R为燃气内燃机混合气体的气体常数，8.31451；
[0270] 根据第十六计算式计算燃气内燃机混合气体密度，所述第十六计算式为：
[0271]
[0272] 式中，Mi，mol为燃气内燃机混合气体的摩尔质量，kg/kmol；
[0273] 根据第十七计算式计算燃气内燃机混合气体完全燃烧后排气中氮气的容积数，所述第十七计算式为：
[0274]
[0275] 式中， 为燃气内燃机混合气体完全燃烧后排气中氮气的容积数，m3/s； 为空气中氮气的百分含量，％； 为天然气中氮气的百分含量，％；
[0276] 根据第十八计算式计算燃气内燃机混合气体完全燃烧后排气中二氧化碳的容积数，所述第十八计算式为：
[0277]
[0278] 式中， 为燃气内燃机混合气体完全燃烧后排气中二氧化碳的容积数，m3/s；为空气中氮气的百分含量，％；CnHm为天然气中各碳氢化合物的百分含量，％；n为天然气成分的碳原子个数；
[0279] 根据第十九计算式计算燃气内燃机混合气体完全燃烧后排气中水蒸汽的容积数，所述第十九计算式为：
[0280]
[0281] 式中， 为燃气内燃机混合气体完全燃烧后排气中水蒸汽的容积数，m3/s；为空气中氮气的百分含量，％；m为天然气成分的氢原子个数；
[0282] 根据第二十计算式计算燃气内燃机混合气体完全燃烧后排气中氧气的容积数，所述第二十计算式为：
[0283]
[0284] 式中， 为燃气内燃机混合气体完全燃烧后排气中氧气的容积数，m3/s； 为空气中氧气的百分含量，％； 为天然气成分中氧气的百分含量，％；
[0285] 根据第二十一计算式计算燃气内燃机混合气体完全燃烧后排气中氩气的容积数，所述第二十一计算式为：
[0286]
[0287] 式中，VAr为燃气内燃机混合气体完全燃烧后排气中氩气的容积数，m3/s；τAr，air为空气中氩气的百分含量，％；τAr，f为天然气成分中氩气的百分含量，％；
[0288] 根据第二十二计算式计算燃气内燃机混合气体完全燃烧后排气中二氧化硫的容积数，所述第二十二计算式为：
[0289]
[0290] 式中， 为燃气内燃机混合气体完全燃烧后排气中二氧化硫的容积数，m3/s；为天然气成分中二氧化硫的百分含量，％； 为天然气成分中硫化氢的百分含
量，％；
[0291] 根据第二十三计算式计算燃气内燃机混合气体完全燃烧后排气的总容积数，所述第二十三计算式为：
[0292]
[0293] 式中，Vpy为燃气内燃机混合气体完全燃烧后排气的总容积数，m3/s；
[0294] 根据第二十四计算式计算燃气内燃机混合气体完全燃烧后排气各成分的体积百分数，所述第二十四计算式为：
[0295]
[0296] 式中，γi-py为燃气内燃机混合气体完全燃烧后排气各成分的体积百分数，％；Vi为3
燃气内燃机混合气体完全燃烧后排气各成分的容积数，m/s；i代表N2、CO2、H2O、O2、Ar、SO2；
[0297] 根据第二十五计算式计算燃气内燃机进气终点压力，所述第二十五计算式为：
[0298]
[0299] 式中，Pin-end为燃气内燃机进气终点压力，MPa；υin为燃气内燃机进气速度，m/s；
[0300] 根据第二十六计算式计算燃气内燃机排气终点压力，所述第二十六计算式为：
[0301]
[0302] 式中，Pout-end为燃气内燃机排气终点压力，MPa；υout为燃气内燃机排气速度，m/s；
[0303] 根据第二十七计算式计算燃气内燃机剩余废气系数，所述第二十七计算式为：
[0304]
[0305] 式中，γreg为燃气内燃机剩余废气系数，无量纲；γc为燃气内燃机净化系数，无量纲；ε为内燃机压缩比，无量纲；Pr为燃气内燃机膨胀冲程终点理论压力，MPa；m为气体膨胀指数，无量纲；
[0306] 根据第二十八计算式计算燃气内燃机剩余废气量，所述第二十八计算式为：
[0307] Vreg＝γreg×V7
[0308] 式中，Vreg为燃气内燃机剩余废气量，m3/s；
[0309] 根据第二十九计算式计算燃气内燃机混合气体量，所述第二十九计算式为：
[0310] Vmix＝γreg+V7
[0311] 式中，Vmix为燃气内燃机混合气体量，m3/s；
[0312] 根据第三十计算式计算燃气内燃机混合气体量，所述第三十计算式为：
[0313] Vp＝γreg+Vpy
[0314] 式中，Vp为燃气内燃机排气气体量，m3/s；
[0315] 根据第三十一计算式计算燃气内燃机充气终点压力，所述第三十一计算式为：
[0316]
[0317] 式中，Pa为燃气内燃机充气终点压力，MPa；ξ1为补充充气系数；
[0318] 根据第三十二计算式计算燃气内燃机充气终点温度，所述第三十二计算式为：
[0319]
[0320] 式中，Ta为燃气内燃机充气终点温度，K；ξr为比热不相等系数；n为多变指数；
[0321] 根据第三十三计算式计算燃气内燃机混合压力下的剩余废气温度，所述第三十三计算式为：
[0322]
[0323] 式中，Tr-mix为燃气内燃机混合压力下的剩余废气温度，K；
[0324] 根据第三十四计算式计算燃气内燃机排气压力下的废气温度，所述第三十四计算式为：
[0325]
[0326] 式中，Tr-out为燃气内燃机排气压力下的废气温度，K；
[0327] 根据第三十五计算式计算燃气内燃机充气系数，所述第三十五计算式为：
[0328]
[0329] 式中，ξc为燃气内燃机充气系数；tair为环境温度，℃；
[0330] 根据
迭代法求解燃气内燃机压缩多变指数n1，具体步骤如下：
[0331] 1)设定一个燃气内燃机压缩多变指数n′1；
[0332] 2)根据第三十六计算式计算燃气内燃机压缩终点温度，所述第三十六计算式为：
[0333]
[0334] 式中，Tc为燃气内燃机压缩终点温度，K；
[0335] 3)根据第三十七计算式计算燃气内燃机压缩终点压力，所述第三十七计算式为：
[0336]
[0337] 式中，Pc为燃气内燃机压缩终点压力，MPa；
[0338] 4)根据第三十八计算式计算燃气内燃机压缩多变指数，所述第三十八计算式为：
[0339] n1＝1.392-0.13×(Tc-aver-273.15)/1000
[0340] 式中，n1为燃气内燃机压缩多变指数；Tc-aver为燃气内燃机压缩终点平均温度，K；
[0341] 5)将n′1与n1做差，得出二者的差值；
[0342] 若差值在预定的误差范围内，则假设的n′1为燃气内燃机压缩多变指数n1；
[0343] 若差值超出预定的误差范围内，则将n′1和n1的平均值作为新的n′1，重新执行上述1)～5)的计算，直到n′1与n1的差值满足设定的误差范围；
[0344] 根据迭代法求解燃气内燃机最高燃烧温度tZ，具体步骤如下：
[0345] 1)设定一个燃气内燃机最高燃烧温度t′Z；
[0346] 2)确定燃气内燃机燃烧产物的焓值hrs；
[0347] 3)根据第三十九计算式计算燃气内燃机最高燃烧温度，所述第三十九计算式为：
[0348]
[0349] 式中，tZ为燃气内燃机最高燃烧温度，℃；hrs为燃烧产物焓值，kJ/(kg·℃)；
[0350] 4)将t′Z与tZ做差，得出二者的差值；
[0351] 若差值在预定的误差范围内，则假设的t′Z为燃气内燃机最高燃烧温度tZ；
[0352] 若差值超出预定的误差范围内，则将t′Z和tZ的平均值作为新的t′Z，重新执行上述1)～4)的计算，直到t′Z与tZ的差值满足设定的误差范围；
[0353] 根据第四十计算式计算燃气内燃机压力升高比，所述第四十计算式为：
[0354]
[0355] 式中，π为燃气内燃机压力升高比；
[0356] 根据第四十一计算式计算燃气内燃机最高燃烧压力，所述第四十一计算式为：
[0357] PZ＝Pc×π×κ
[0358] 式中，PZ为燃气内燃机最高燃烧压力，MPa；κ为最高燃烧压力的修正系数；
[0359] 根据迭代法求解燃气内燃机膨胀多变指数n2，具体步骤如下：
[0360] 1)设定一个燃气内燃机膨胀多变指数n′2；
[0362] 2)根据第四十二计算式计算燃气内燃机膨胀终点温度，所述第四十二计算式为：
[0363]
[0364] 式中，Tb为燃气内燃机膨胀终点温度，K；n2为燃气内燃机膨胀多变指数；
[0365] 3)根据第四十三计算式计算燃气内燃机膨胀终点压力，所述第四十三计算式为：
[0366]
[0367] 式中，Pb为燃气内燃机膨胀终点压力，MPa；
[0368] 4)根据第四十四计算式计算燃气内燃机膨胀多变指数，所述第四十四计算式为：
[0369]
[0370] 式中，n2为燃气内燃机膨胀多变指数；Tb-aver为燃气内燃机膨胀终点平均温度，K；
[0371] 5)将n′2与n2做差，得出二者的差值；
[0372] 若差值在预定的误差范围内，则假设的n′2为燃气内燃机膨胀多变指数n2；
[0373] 若差值超出预定的误差范围内，则将n′2和n2的平均值作为新的n′2，重新执行上述1)～5)的计算，直到n′2与n2的差值满足设定的误差范围；
[0374] 根据第四十五计算式计算燃气内燃机处初膨胀比，所述第四十五计算式为：
[0375]
[0376] 式中，ρ0为燃气内燃机处初膨胀比；
[0377] 根据第四十六计算式计算燃气内燃机后膨胀比，所述第四十六计算式为：
[0378]
[0379] 式中，δ为燃气内燃机处后膨胀比；
[0380] 根据第四十七计算式计算燃气内燃机平均指示压力，所述第四十七计算式为：
[0381]
[0382] 式中，Pi为燃气内燃机平均指示压力，MPa；
[0383] 根据第四十八计算式计算燃气内燃机平均有效压力，所述第四十八计算式为：
[0384] Pmi＝Pi×ηm
[0385] 式中，Pmi为燃气内燃机平均有效压力，MPa；ηm为燃气内燃机机械效率，％；
[0386] 根据第四十九计算式计算燃气内燃机有效功率，所述第四十九计算式为：
[0387]
[0388] 式中，Pe为燃气内燃机有效功率，kW；nc为燃气内燃机气缸数；Vw为燃气内燃机的工作容积，m3；ω为燃气内燃机转速，r/min；
[0389] 根据第五十计算式计算燃气内燃机发电功率，所述第五十计算式为：
[0390] Pele＝Pe×ηele
[0391] 式中，Pele为燃气内燃机发电功率，kW；ηele为燃气内燃机发电机效率，％；
[0392] 根据第五十一计算式计算燃气内燃机发电效率，所述第五十一计算式为：
[0393]
[0394] 式中，η为燃气内燃机发电效率，％；
[0395] 根据第五十二计算式计算燃气内燃机高温循环水在额定工况下带走的热量，所述第五十二计算式为：
[0396] QHTW-in design＝Cpw×MHTW-in design×(tHTW-out-tHTW-in)/3.6
[0397] 式中，QHTW-in design为燃气内燃机高温循环水在额定工况下带走的热量，kJ；Cpw为水的比热，kJ/(kg·℃)；MHTW-in design为燃气内燃机高温循环水在额定工况下的流量，kg/h；tHTW-out为燃气内燃机高温循环水在额定工况下出口处的温度，℃；tHTW-in为燃气内燃机高温循环水在额定工况下进口处的温度，℃；
[0398] 根据第五十三计算式计算燃气内燃机高温循环水在额定工况下带走的热量占总热量的百分比，所述第五十三计算式为：
[0399]
[0400] 式中，ηHTM为燃气内燃机高温循环水在额定工况下带走的热量占总热量的百分比，％；Pele-in design为在额定工况下燃气内燃机发电功率，kW；ηele-in design为在额定工况下燃气内燃机发电效率，％；
[0401] 根据第五十四计算式计算燃气内燃机高温循环水在运行工况下带走的热量，所述第五十四计算式为：
[0402] QHTW-off design＝Mf×QNG，L×1000×3600×ηHTM
[0403] 式中，QHTW-off design为燃气内燃机高温循环水在运行工况下带走的热量，kJ；
[0404] 根据第五十五计算式计算燃气内燃机高温循环水在运行工况下的流量，所述第五十五计算式为：
[0405]
[0406] 式中，MHTW-off design为燃气内燃机高温循环水在运行工况下的流量，kg/s；t′HTW-out为燃气内燃机高温循环水在运行工况下出口处的温度，℃；t′HTW-in为燃气内燃机高温循环水在运行工况下进口处的温度，℃；
[0407] 根据第五十六计算式计算燃气内燃机低温循环水在额定工况下带走的热量，所述第五十六计算式为：
[0408] QLTW-in design＝Cpw×MLTW-in design×(tLTW-out-tLTW-in)/3.6
[0409] 式中，QLTW-in design为燃气内燃机低温循环水在额定工况下带走的热量，kJ；MLTW-in design为燃气内燃机低温循环水在额定工况下的流量，kg/h；tLTW-out为燃气内燃机低温循环水在额定工况下出口处的温度，℃；tLTW-in为燃气内燃机低温循环水在额定工况下进口处的温度，℃；
[0410] 根据第五十七计算式计算燃气内燃机低温循环水在额定工况下带走的热量占总热量的百分比，所述第五十七计算式为：
[0411]
[0412] 式中，ηLTM为燃气内燃机低温循环水在额定工况下带走的热量占总热量的百分比，％；
[0413] 根据第五十八计算式计算燃气内燃机低温循环水在运行工况下带走的热量，所述第五十八计算式为：
[0414] QLTW-off design＝Mf×QNG，L×1000×3600×ηLTM
[0415] 式中，QLTW-off design为燃气内燃机低温循环水在运行工况下带走的热量，kJ；
[0416] 根据第五十九计算式计算燃气内燃机高温循环水在运行工况下的流量，所述第五十九计算式为：
[0417]
[0418] 式中，MLTW-off design为燃气内燃机低温循环水在运行工况下的流量，kg/s；t′LTW-out为燃气内燃机低温循环水在运行工况下出口处的温度，℃；t′LTW-in为燃气内燃机低温循环水在运行工况下进口处的温度，℃；
[0419] 根据第六十计算式计算燃气内燃机烟气在额定工况下带走的热量，所述第六十计算式为：
[0420] Qg-indesign＝Mg×hg-indesign
[0421] 式中，Qg-indesign为燃气内燃机烟气在额定工况下带走的热量，kJ；Mg为燃气内燃机烟气在额定工况下的流量，m3/s；hg-indesign为燃气内燃机烟气在额定工况下的焓值，kJ/m3；
[0422] 根据第六十一计算式计算燃气内燃机烟气在额定工况下带走的热量占总热量的百分比，所述第六十一计算式为：
[0423]
[0424] 式中，ηg为燃气内燃机烟气在额定工况下带走的热量占总热量的百分比，％；
[0425] 根据第六十二计算式计算燃气内燃机在额定工况下散热损失占总热量的百分比，所述第六十二计算式为：
[0426] ηloss＝1-η-ηHTM-ηLTM-ηg
[0427] 式中，ηloss为燃气内燃机在额定工况下散热损失占总热量的百分比，％；
[0428] 根据第六十三计算式计算燃气内燃机在运行工况下散热量，所述第六十三计算式为：
[0429] Qloss-off design＝Mf×QNG，L×1000×3600×ηloss
[0430] 式中，Qloss-off design为燃气内燃机在运行工况下散热量，kJ；
[0431] 根据第六十四计算式计算燃气内燃机在运行工况下排烟所携带的热量，所述第六十四计算式为：
[0432] Qg-off design＝Mf×QNG，L×1000×3600-Pele×3600-QHTW-off design-QLTW-off design-Qloss-off design
[0433] 式中，Qg-off design为燃气内燃机在运行工况下排烟所携带的热量，kJ；
[0434] 根据迭代法求解燃气内燃机排烟温度tpy，具体步骤如下：
[0435] 1)设定一个燃气内燃机排烟温度t′py；
[0436] 2)确定燃气内燃机排烟的焓值hg-offdesign；
[0437] 3)根据第六十五计算式计算燃气内燃机在运行工况下排烟流量，所述第六十五计算式为：
[0438]
[0439] 式中，V7′为核算的燃气内燃机在运行工况下排烟流量，m3/s；hg-offdesign为燃气内燃机烟气在运行工况下的焓值，kJ/m3；
[0440] 4)将V7′与V7做差，得出二者的差值；
[0441] 若差值在预定的误差范围内，则假设的t′py为燃气内燃机排烟温度tpy；
[0442] 若差值超出预定的误差范围内，则将t′py和tpy的平均值作为新的t′py，重新执行上述1)～4)的计算，直到V7′与V7的差值满足设定的误差范围；
[0443] 以上，仅为本发明的较佳实施案例，但本发明的保护范围并不局限于此，任何熟悉
本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替代，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以
权利要求所界定的保护范围为准。
