本发明涉及一种热能工程领域的制冷设备，尤其涉及一种直燃型溴化锂制冷机热力特性的计算方法。
背景技术：
天然气分布式能源是分布在用户端的能源综合利用系统。它是以气体燃料为主、可再生能源为辅，通过燃气轮机或者内燃机发电，将其尾部烟气通过能量转换设备生产出用户所需要的热能利用形式，如蒸汽、生活热水、采暖热水或制冷用热能等；按照“温度对口、梯级利用”的用能原则，使得天然气分布式能源系统的综合能源利用率大于70％。发展天然气分布式能源是降低能源成本、提升能源效率、改善大气环境的一种有效技术途径。
直燃型溴化锂制冷机是天然气分布式能源系统的主要调峰设备之一，其制冷量极易受到燃气成分、环境参数等因素的影响。在某些工况条件下，直燃型溴化锂制冷机的制冷量很难与项目负荷需求相匹配。确定变工况条件下，直燃型溴化锂制冷机的制冷量对天然气分布式能源系统的方案配置起到了决定性作用。
因此，要构建一个适用于直燃型溴化锂制冷机热力特性的计算方法，是合理配置天然气分布式能源系统的基础手段，能够为天然气分布式能源系统运行方案提供基础数据，具有重要的实用意义。
技术实现要素：
针对上述问题，本发明提供一种适用于直燃型溴化锂制冷机热力特性的计算方法。
为达到上述目的，本发明一种确定直燃型溴化锂制冷机热力特性的计算方法，所述的热力特性的计算方法包括：
s1，燃气燃烧计算及冷却介质参数的确定；
s2，制冷机设计参数的选定；
s3，热力循环点参数的计算；
s4，设备负荷、传热面积计算及制冷量的确定。
进一步地，燃气燃烧计算及冷却介质参数的确定还包括如下步骤：
s11，根据直燃型溴化锂制冷机项目所在的地点，采集直燃型溴化锂制冷机周围环境的大气压力、温度、湿度或者海拔高度参数信息；
s12，根据s11的参数信息，计算空气中各成分的体积百分数；根据其温度及成分，确定冷空气的比热；
s13，根据直燃型溴化锂制冷机项目所在地，收集直燃型溴化锂制冷机所用天然气成分，通过天然气热力计算公式计算得到项目所使用的天然气的低位热值、高位热值、密度等参数，其具体步骤包括：
(1)根据第一计算式计算天然气的低位热值，所述第一计算式为：
式中，qng，l为天然气低位发热量，mj/nm3；mi为天然气各组分的摩尔体积百分数，％；qi，l为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kj/mol；png为天然气的压力，kpa；tng为天然气的温度，k；z为天然气压缩系数；
(2)根据第二计算式计算天然气的高位热值，所述第二计算式为：
式中，qng，h为天然气高位发热量，mj/nm3；qi，h为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kj/mol；
(3)根据第三计算式计算天然气的密度，所述第三计算式为：
ρng＝(png/(r·tng))×∑mi·mi
式中，ρng为天然气密度，kg/nm3；mi为天然气各组分的摩尔质量，kg/kmol。
s14，根据第四计算式计算单位体积天然气在直燃型溴化锂制冷机中燃烧所需要的空气量，所述第四计算式为：
式中，vair为单位体积天然气在直燃型溴化锂制冷机中燃烧所需要的空气量，m3/s；cnhm为天然气成分，％；n为天然气成分的碳原子个数；m为天然气成分的氢原子个数；
s15，确定天然气燃烧时的空气过量系数α；按以下方法计算烟气各成分体积百分数：
(1)通过第五计算式计算单位体积天然气燃烧产生的烟气中氮气的容积数，所述第五计算式为：
式中，
(2)通过第六计算式计算单位体积天然气燃烧产生的烟气中二氧化碳的容积数，所述第六计算式为：
式中，
(3)通过第七计算式计算单位体积天然气燃烧产生的烟气中水蒸汽的容积数，所述第七计算式为：
式中，
(4)通过第八计算式计算单位体积天然气燃烧产生的烟气中二氧化硫的容积数，所述第八计算式为：
式中，
(5)通过第九计算式计算单位体积天然气燃烧产生的烟气中剩余氧气的容积数，所述第九计算式为：
式中，
(6)通过第十计算式计算单位体积天然气燃烧产生的烟气中氩气的容积数，所述第十计算式为：
mar＝α×vair×τar/100+τar-gas/100
式中，mar为单位体积天然气燃烧产生的烟气中氩气的容积数，m3/s；τar为环境空气中氩气的容积百分数，％；τar-gas为天然气中氩气的容积百分数，％；
(7)通过第十一计算式计算单位体积天然气燃烧生成的烟气总容积数，所述第十一计算式为：
式中，mfg为单位体积天然气燃烧生成的烟气总容积数，m3/s；
(8)通过第十二计算式计算单位体积天然气燃烧生成的烟气各成分体积百分数，所述第十二计算式为：
式中，γi为单位体积天然气燃烧生成的烟气各成分体积百分数，％；其中，i为n2、co2、h2o、so2、o2、ar；
s16，根据直燃型溴化锂制冷机进口及出口的烟气温度及成分，确定其进口及出口的烟气的比热；
s17，判断直燃型溴化锂制冷机是否有空气预热装置，确定空气进入燃烧器的温度，进而确定进入直燃型溴化锂制冷机的燃烧器的空气比热；
s18，根据第十三计算式计算空气预热器的平均温差，所述第十三计算式为：
δtm1＝(tpy-tair-in)-0.35×min{(tpy-tpy-preh)，(tair-preh-tair-in)}-0.65×max{(tpy-tpy-preh)，(tair-preh-tair-in)}
式中，δtm1为空气预热器的平均温差，℃；tpy为直燃型溴化锂制冷机的排烟温度，℃；tpy-preh为直燃型溴化锂制冷机空气预热器出口处的排烟温度，℃；tair-in为直燃型溴化锂制冷机燃烧器进口空气温度，℃；tair-preh为直燃型溴化锂制冷机空气有预热器时燃烧器进口空气温度，℃；
s19，确定加热燃气损失的附加系数，确定冷媒水的进、出口温度及冷却水进、出口温度；
进一步地，制冷机设计参数的选定还包括如下步骤：
1)设定直燃型溴化锂制冷机流程结构，是串联还是并联；
2)设定直燃型溴化锂制冷机吸收器与冷凝器的连接方式，是串联还是并联；
3)设定冷却水总温升，并对吸收器出口和冷凝器出口的冷却水温升进行分配；
4)设定吸收器出口处稀溶液与吸收器冷却水温差；
5)设定溴化锂溶液冷凝温度与冷却水出口处温度的差值；
6)设定溴化锂溶液蒸发温度与冷媒水出口处温度的差值；
7)设定吸收器压损及再循环倍数；
8)设定溴化锂溶液从吸收器中流到低压发生器后浓度变化量；
9)设定溴化锂溶液在高压发生器中浓度变化量；
10)设定凝结水离开凝结水换热器的温度；
11)设定蒸发器的再循环倍数。
进一步地，热力循环点参数的计算还包括如下步骤：
1)根据第十四计算式计算蒸发器出口处冷剂蒸气的温度，所述第十四计算式为：
t0＝tx′-δt
式中，t0为蒸发器出口处冷剂蒸气的温度，℃；tx′为冷媒水出口温度，℃；δt为蒸发温度与冷媒水出口温度的差值，℃；
2)根据蒸发器出口处冷剂蒸气的温度确定其压力p0及焓值h0；
3)根据第十五计算式计算从蒸发器流入到吸收器内的饱和水的压力，所述第十五计算式为：
p1＝p0-δp0
式中，p1为从蒸发器流入到吸收器内的饱和水的压力，mpa；δp0为吸收器进口处与出口处稀溶液的压力差，mpa；
4)根据从蒸发器流入到吸收器内的饱和水的压力确定其温度t1及焓值h1；
5)根据第十六计算式计算从蒸发器流入到吸收器内的蒸汽压力，所述第十六计算式为：
p1a＝p1
式中，p1a为从蒸发器流入到吸收器内的蒸汽压力，mpa；
6)根据从蒸发器流入到吸收器内的蒸汽压力确定其温度t1a及焓值h1a；
7)根据第十七计算式计算吸收器出口处稀溶液的温度，所述第十七计算式为：
t2＝tw+δtw
式中，t2为吸收器出口处稀溶液的温度，℃；tw为吸收器冷却水进口温度，℃；δtw为吸收器出口处稀溶液与吸收器冷却水温差，℃；
8)根据第十八计算式计算吸收器出口处稀溶液的压力，所述第十八计算式为：
p2＝p1
式中，p2为吸收器出口处稀溶液的压力，mpa；
9)根据第十九计算式计算吸收器出口处稀溶液的浓度，所述第十九计算式为：
式中，ξ2为吸收器出口处稀溶液的浓度，％；a、b、c、d为吸收器出口处稀溶液浓度的计算系数，其取值根据吸收器出口处稀溶液压力下的饱和温度计算得到；
10)根据第二十计算式计算吸收器出口处稀溶液的焓值，所述第二十计算式为
式中，h2为吸收器出口处稀溶液的焓值，kj/kg；e、f、g为吸收器出口处稀溶液焓值的计算系数，其取值根据吸收器出口处稀溶液的温度和浓度计算得到；
11)根据第二十一计算式计算冷凝器出口处冷剂水的温度，所述第二十一计算式为：
t3＝tcw+δtk
式中，t3为冷凝器出口处冷剂水的温度，℃；tcw为冷却水出口温度，℃；δtk为冷凝温度与冷却水出口温度的差值，℃；
12)根据冷凝器出口处冷剂水的温度确定其压力p3及焓值h3；
13)根据第二十二计算式计算低压发生器产生的水蒸汽的压力，所述第二十二计算式为：
p3a＝p3
式中，p3a为低压发生器产生的水蒸汽的压力，mpa；
14)根据低压发生器产生的水蒸汽的压力，确定其焓值h3a；
15)根据第二十三计算式计算高压发生器内的饱和水的压力，所述第二十三计算式为：
p3b＝pr
式中，p3b为高压发生器内的饱和水的压力，mpa；pr为为高压发生器的压力，mpa；
16)根据高压发生器内的饱和水的压力确定其温度t3b及焓值h3b；
17)根据第二十四计算式计算高压发生器内水蒸汽的压力，所述第二十四计算式为：
p3c＝pr
式中，p3c为高压发生器内水蒸汽的压力，mpa；
18)根据高压发生器内水蒸汽的压力确定其焓值h3c；
19)根据第二十五计算式计算低压发生器出口处浓溶液的压力，所述第二十五计算式为：
p4＝p3
式中，p4为低压发生器出口处浓溶液的压力，mpa；
20)根据第二十六计算式计算低压发生器出口处浓溶液的浓度，所述第二十六计算式为：
ξ4＝ξ2+δξ1(直燃吸收式溴化锂制冷机流程结构为并联流程时)
ξ4＝ξ5+δξ1(直燃吸收式溴化锂制冷机流程结构为串联流程时)
式中，ξ4为低压发生器出口处浓溶液的浓度，％；δξ1为溴化锂溶液从吸收器中流到低压发生器后浓度变化量，％；ξ5为高温热交换器出口处的浓溶液(饱和溶液)的浓度，％；
21)根据第二十七计算式和第二十计算式分别计算低压发生器出口处浓溶液的温度t4和焓值h4；
22)根据第二十八计算式计算高温热交换器出口处的浓溶液(饱和溶液)的压力，所述第二十八计算式为：
p5＝p3(直燃吸收式溴化锂制冷机流程结构为串联流程时)
式中，p5为高温热交换器出口处的浓溶液(饱和溶液)的压力，mpa；
23)根据第二十九计算式计算高温热交换器出口处的浓溶液(饱和溶液)的浓度，所述第二十九算式为：
ξ5＝ξ12(直燃吸收式溴化锂制冷机流程结构为串联流程时)
式中，ξ12为高压发生器出口处浓溶液的浓度，％；
24)根据第十九计算式和第二十计算式分别计算高温热交换器出口处的浓溶液(饱和溶液)的温度t5和焓值h5(直燃吸收式溴化锂制冷机流程结构为串联流程时)；
25)根据第三十计算式计算低温热交换器的稀溶液的压力，所述第三十计算式为：
p7＝p3
式中，p7为低温热交换器的稀溶液的压力，mpa；
26)根据第三十一计算式计算低温热交换器的稀溶液的浓度，所述第三十一计算式为：
ξ7＝ξ2
式中，ξ7为低温热交换器的稀溶液的浓度，％；
27)根据第十九计算式计算低温热交换器的稀溶液的温度t7；
28)根据第三十二计算式计算低温热交换器的稀溶液的焓值，所述第三十二计算式为：
式中，h7为低温热交换器的稀溶液的焓值，kj/kg；a1为低位发生器循环倍率；h8为低温热交换器出口处的浓溶液(过冷溶液)的焓值，kj/kg；
29)根据第三十三计算式计算低温热交换器出口处的浓溶液(过冷溶液)的压力，所述第三十三计算式为：
p8＝p3
式中，p8为低温热交换器出口处的浓溶液(过冷溶液)的压力，mpa；
30)根据第三十四计算式计算低温热交换器出口处的浓溶液(过冷溶液)的浓度，所述第三十四计算式为：
ξ8＝ξ4
式中，ξ8为低温热交换器出口处的浓溶液(过冷溶液)的浓度，％；
31)根据第三十五计算式计算低温热交换器出口处的浓溶液(过冷溶液)的温度，所述第三十五计算式为：
t8＝t2+δth
式中，t8为低温热交换器出口处的浓溶液(过冷溶液)的温度，℃；δth为低温热交换器出口处的浓溶液与入口处的温差，℃；
32)根据第二十计算式计算低温热交换器出口处的浓溶液(过冷溶液)的焓值h8；
33)根据第三十六计算式计算吸收器的喷淋溶液的压力，所述第三十六计算式为：
p9＝p2
式中，p9为吸收器的喷淋溶液的压力，mpa；
34)根据第三十七计算式计算吸收器的喷淋溶液的浓度，所述第三十七计算式为：
式中，ξ9为吸收器的喷淋溶液的浓度，％；ξ13为高温热交换器出口处的浓溶液(过冷溶液)的浓度，％；qmh为单位时间进入高位发生器的稀溶液质量流量，kg/s；qmdh为高压发生器中产生的水蒸汽量，kg/s；qml为单位时间进入低位发生器的稀溶液质量流量，kg/s；qmdl为低压发生器中产生的水蒸汽量，kg/s；qmd为冷剂水质量流量，kg/s；f为吸收器的再循环倍数；
35)根据第十九计算式计算吸收器的喷淋溶液的温度t9；
36)根据第三十八计算式计算吸收器的喷淋溶液的焓值，所述第三十八计算式为：
式中，h9为吸收器的喷淋溶液的焓值，kj/kg；h13为高温热交换器出口处的浓溶液(过冷溶液)的焓值，kj/kg；
37)根据第三十九计算式计算高温热交换器的稀溶液的压力，所述第三十九计算式为：
p10＝pr
式中，p10为高温热交换器的稀溶液的压力，mpa；
38)根据第四十计算式计算高温热交换器的稀溶液的浓度，所述第四十计算式为：
ξ10＝ξ2
式中，ξ10为高温热交换器的稀溶液的浓度，％；
39)根据第十九计算式计算高温热交换器的稀溶液的温度t10；
40)根据第四十一计算式计算高温热交换器的稀溶液的焓值，所述第四十一计算式为：
式中，h10为高温热交换器的稀溶液的焓值，kj/kg；ah为高压发生器循环倍率；h12为高压发生器出口处浓溶液的焓值，kj/kg；
41)根据第四十二计算式计算高压发生器出口处浓溶液的压力，所述第四十二计算式为：
p12＝pr
式中，p12为高压发生器出口处浓溶液的压力，mpa；
42)根据第四十三计算式计算高压发生器出口处浓溶液的浓度，所述第四十三计算式为：
ξ12＝ξ2+δξ2
式中，ξ12为高压发生器出口处浓溶液的浓度，％；δξ2为溴化锂溶液从吸收器中流到高压发生器后浓度变化量，％；
43)根据第十九计算式和第二十计算式分别计算高压发生器出口处浓溶液的温度t12和焓值h12；
44)根据第四十四计算式计算高温热交换器出口处的浓溶液(过冷溶液)的压力，所述第四十四计算式为：
p13＝pr(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，p13为高温热交换器出口处的浓溶液(过冷溶液)的压力，mpa；
由于高温热交换器入口处浓溶液的温度较高，所以在热交换器中温降大于低温热交换器中浓溶液的温升，通常浓溶液的出口温度在60～70℃范围内(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
45)根据第四十五计算式计算高温热交换器出口处的浓溶液(过冷溶液)的浓度，所述第四十五计算式为：
ξ13＝ξ12(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，ξ13为高温热交换器出口处的浓溶液(过冷溶液)的浓度，％：
46)根据第二十计算式计算高温热交换器出口处的浓溶液(过冷溶液)的焓值h10(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
47)根据第四十六计算式计算低压发生器循环倍率，所述第四十六计算式为：
式中，al为低压发生器循环倍率；
48)根据第四十七计算式计算高压发生器循环倍率，所述第四十七计算式为：
式中，ah为高压发生器循环倍率；
进一步地，设备负荷、传热面积计算及制冷量的确定还包括如下步骤：
1)根据第四十八计算式计算制冷机中冷剂水的流量，所述第四十八计算式为：
式中，qmd为制冷机中冷剂水的流量，kg/s；q0为直燃型溴化锂制冷机的制冷量，kw；
2)根据第四十九计算式计算低压发生器中产生的水蒸汽量，所述第四十九计算式为：
式中，qmdl为低压发生器中产生的水蒸汽量，kg/s；
3)根据第五十计算式计算高压发生器中产生的水蒸汽量，所述第五十计算式为：
qmdh＝qmd-qmdl
式中，qmdh为高压发生器中产生的水蒸汽量，kg/s；
4)根据第五十一计算式计算单位时间进入高压发生器的稀溶液质量，所述第五十一计算式为：
qmh＝ah×qmdh
式中，qmh为单位时间进入高压发生器的稀溶液质量，kg/s；
5)根据第五十二计算式计算单位时间进入低压发生器的稀溶液质量，所述第五十二计算式为：
qml＝al×qmdl
式中，qml为单位时间进入高低压发生器的稀溶液质量，kg/s；
6)根据第五十三计算式计算高压发生器中产生1kg水蒸气时，在吸收器内产生的热负荷，所述第五十三计算式为：
qah＝(ah-1)×h13-ah×h2+h1a(仅直燃吸收式溴化锂制冷机流程结构为并联流程时)；
式中，qah为高压发生器中产生1kg水蒸气时在吸收器内产生的热负荷，kj/kg；
7)根据第五十四计算式计算低压发生器中产生1kg水蒸气时在吸收器内产生的热负荷，所述第五十四计算式为：
qal＝(al-1)×h8-al×h2+h1a
式中，qal为低压发生器中产生1kg水蒸气时在吸收器内产生的热负荷，kj/kg；
8)根据第五十五计算式计算高压发生器中产生1kg水蒸气时的热负荷，所述第五十五计算式为：
qgh＝(ah-1)×h12-ah×h10+h3c
式中，qgh为高压发生器中产生1kg水蒸气时的热负荷，kj/kg；
9)根据第五十六计算式计算低压发生器中产生1kg水蒸气时的热负荷，所述第五十六计算式为：
式中，qlh为低压发生器中产生1kg水蒸气时的热负荷，kj/kg；
10)根据第五十七计算式计算吸收器热负荷，所述第五十七计算式为：
qa＝qmdh×qh+qmdl×ql(直燃吸收式溴化锂制冷机流程结构为并联流程时)；
qa＝qmdl×ql(直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，qa为吸收器热负荷，kw；
11)根据第五十八计算式计算高压发生器热负荷，所述第五十八计算式为：
qgh＝qmdh×qgh
式中，qgh为高压发生器热负荷，kw；
12)根据第五十九计算式计算低压发生器热负荷，所述第五十九计算式为：
qgl＝qmdl×qgl
式中，qgl为低压发生器热负荷，kw；
13)根据第六十计算式计算冷凝器热负荷，所述第六十计算式为：
qk＝qmdh×(h3c-h3)+qmdl×(h3a-h3)(直燃吸收式溴化锂制冷机流程结构为并联流程时)；
qk＝qmdh×(h3c-h3)-qgh+qmdl×(h3a-h3)(直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，qk为冷凝器热负荷，kw；
14)蒸发器热负荷等于溴化锂制冷机的制冷量q0；
15)根据第六十一计算式计算高压溶液热交换器的总热负荷，所述第六十一计算式为：
qexh＝qmh×(h10-h2)
式中，qexh为高压溶液热交换器的总热负荷，kw；
16)根据第六十二计算式计算低压溶液热交换器的总热负荷，所述第六十二计算式为：
qexl＝qml×(h7-h2)
式中，qexl为低压溶液热交换器的总热负荷，kw；
17)凝结水热交换器的总热负荷qexv等于低温水总放热量荷qv；
18)根据第六十三计算式计算直燃型溴化锂制冷机的热力系数，所述第六十三计算式为：
式中，ζ为直燃型溴化锂制冷机的热力系数；
19)根据第六十四计算式计算直燃型溴化锂制冷机的加热工作热水流量，所述第六十四计算式为：
式中，qmv2为直燃型溴化锂制冷机的加热工作热水流量，kg/s；
20)根据第六十五计算式计算直燃型溴化锂制冷机的冷媒水流量，所述第六十五计算式为：
式中，qv0为直燃型溴化锂制冷机的冷媒水流量，kg/s；tx″为冷媒水进口温度，℃；tx′为冷媒水出口温度，℃；
21)根据第六十六计算式计算直燃型溴化锂制冷机的冷却水泵的流量，所述第六十六计算式为：
式中，
22)根据第六十七计算式计算直燃型溴化锂制冷机的蒸发器泵的流量，所述第六十七计算式为：
式中，qvd为直燃型溴化锂制冷机的蒸发器泵的流量，m3/h；f0为蒸发器的再循环倍数；
23)根据第六十八计算式计算直燃型溴化锂制冷机进入发生器的溶液密度，所述第六十八计算式为：
式中，ρa为直燃型溴化锂制冷机的吸收器溶液密度，kg/l；a0、a1、a2、a3、a4、a5、a6为直燃型溴化锂制冷机的吸收器溶液密度的计算系数；
24)根据第六十九计算式计算直燃型溴化锂制冷机吸收器的溶液密度，所述第六十九计算式为：
式中，ρ9为直燃型溴化锂制冷机的吸收器溶液密度，kg/l；a0、a1、a2、a3、a4、a5、a6为直燃型溴化锂制冷机的吸收器溶液密度的计算系数；
25)根据第七十计算式计算直燃型溴化锂制冷机发生器泵的流量，所述第七十计算式为：
式中，qvg为直燃型溴化锂制冷机发生器泵的流量，m3/h；
26)根据第七十一计算式计算直燃型溴化锂制冷机吸收器泵的流量，所述第七十一计算式为：
式中，qva为直燃型溴化锂制冷机吸收器泵的流量，m3/h；f为吸收器的再循环倍数；
27)根据第七十二计算式计算直燃型溴化锂制冷机的高压发生器的传热面积，所述第七十二计算式为：
式中，fgh为直燃型溴化锂制冷机的高压发生器的传热面积，m2；kgh为直燃型溴化锂制冷机高压发生器的传热系数，w/(m2·℃)；tf-in为加热工作热水进口温度，℃；bg1为发生器换热过程中，冷流体换热系数，无量纲；
28)根据第七十三计算式计算直燃型溴化锂制冷机的低压发生器的传热面积，所述第七十三计算式为：
式中，fgl为直燃型溴化锂制冷机的低压发生器的传热面积，m2；kgl为直燃型溴化锂制冷机低压发生器的传热系数，w/(m2·℃)；
29)根据第七十四计算式计算直燃型溴化锂制冷机的冷凝器的传热面积，所述第七十四计算式为：
式中，fk为直燃型溴化锂制冷机的冷凝器的传热面积，m2；kk为直燃型溴化锂制冷机冷凝器的传热系数，w/(m2·℃)；bk1为冷凝器换热过程中，冷流体换热系数，无量纲；
30)根据第七十五计算式计算直燃型溴化锂制冷机的吸收器的传热面积，所述第七十五计算式为：
式中，fa为直燃型溴化锂制冷机的吸收器的传热面积，m2；ka为直燃型溴化锂制冷机吸收器的传热系数，w/(m2·℃)；aa1为吸收器换热过程中，热流体换热系数，无量纲；ba1为吸收器换热过程中，冷流体换热系数，无量纲；
31)根据第七十六计算式计算直燃型溴化锂制冷机的蒸发器的传热面积，所述第七十六计算式为：
式中，f0为直燃型溴化锂制冷机的蒸发器的传热面积，m2；k0为直燃型溴化锂制冷机蒸发器的传热系数，w/(m2·℃)；b01为蒸发器换热过程中，冷流体换热系数，无量纲；
32)根据第七十七计算式计算直燃型溴化锂制冷机的高温溶液热交换器的传热面积，所述第七十七计算式为：
式中，fexh为直燃型溴化锂制冷机的高温溶液热交换器的传热面积，m2；kexh为直燃型溴化锂制冷机高温溶液热交换器的传热系数，w/(m2·℃)；aex1为高温溶液热交换器换热过程中，热流体换热系数，无量纲；bex1为高温溶液热交换器换热过程中，冷流体换热系数，无量纲；
33)根据第七十八计算式计算直燃型溴化锂制冷机低温溶液热交换器的传热面积，所述第七十八计算式为：
式中，fexl为直燃型溴化锂制冷机的低温溶液热交换器的传热面积，m2；kexl为直燃型溴化锂制冷机低温溶液热交换器的传热系数，w/(m2·℃)；aex2为低温溶液热交换器换热过程中，热流体换热系数，无量纲；bex2为低温溶液热交换器换热过程中，冷流体换热系数，无量纲；
34)根据第七十九计算式计算直燃型溴化锂制冷机的高压发生器的对数平均温差，所述第七十九计算式为：
式中，δtm为直燃型溴化锂制冷机的高压发生器的对数平均温差，℃；tf-out为加热工作热水出口温度，℃；
35)根据第八十计算式计算直燃型溴化锂制冷机的高压发生器的核算对数平均温差，所述第八十计算式为：
式中，δt′m为直燃型溴化锂制冷机的高压发生器的核算对数平均温差，℃；
本发明一种确定直燃型溴化锂制冷机热力特性的计算方法，通过计算得到直燃型溴化锂制冷机的制冷量，具有以下优点：
1、克服了现有技术条件下因加热工作热水温度、压力及流量的变化，而无法准确计算直燃型溴化锂制冷机制冷量的计算得以顺利进行；
2、本发明的直燃型溴化锂制冷机热力特性的计算方法，能够通过计算结果反映直燃型溴化锂制冷机的运行情况，进而为天然气分布式能源系统的配置方案及运行方式、优化运行提供指导，达到节能、降低损耗的目的。
附图简要说明
图1是本发明一种确定直燃型溴化锂制冷机热力特性的计算方法示意图。
具体实施方式
下面结合说明书附图对本发明做进一步的描述。
如图1所示为一种确定直燃型溴化锂制冷机热力特性的计算方法，其过程包含了s1燃气燃烧计算及冷却介质参数的确定，s2制冷机设计参数的选定，s3热力循环点参数的计算，s4设备负荷、传热面积计算及制冷量的确定。所述计算具体如下：
首先根据直燃型溴化锂制冷机项目所在的地点，采集直燃型溴化锂制冷机周围环境的大气压力、温度、湿度或者海拔高度参数信息；
计算空气中各成分的体积百分数；根据其温度及成分，确定冷空气的比热；
根据直燃型溴化锂制冷机项目所在地，收集直燃型溴化锂制冷机所用天然气成分，通过天然气热力计算公式计算得到项目所使用的天然气的低位热值、高位热值、密度等参数，其具体步骤包括：
根据第一计算式计算天然气的低位热值，所述第一计算式为：
式中，qng，l为天然气低位发热量，mj/nm3；mi为天然气各组分的摩尔体积百分数，％；qi，l为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kj/mol；png为天然气的压力，kpa；tng为天然气的温度，k；z为天然气压缩系数；
根据第二计算式计算天然气的高位热值，所述第二计算式为：
式中，qng，h为天然气高位发热量，mj/nm3；qi，h为天然气各组分在不同的燃烧和计量参比条件下的理想气体摩尔低位发热量，kj/mol；
根据第三计算式计算天然气的密度，所述第三计算式为：
ρng＝(png/(r·tng))×∑mi·mi
式中，ρng为天然气密度，kg/nm3；mi为天然气各组分的摩尔质量，kg/kmol。
根据第四计算式计算单位体积天然气在直燃型溴化锂制冷机中燃烧所需要的空气量，所述第四计算式为：
式中，vair为单位体积天然气在直燃型溴化锂制冷机中燃烧所需要的空气量，m3/s；cnhm为天然气成分，％；n为天然气成分的碳原子个数；m为天然气成分的氢原子个数；
确定天然气燃烧时的空气过量系数α；按以下方法计算烟气各成分体积百分数：
通过第五计算式计算单位体积天然气燃烧产生的烟气中氮气的容积数，所述第五计算式为：
式中，
通过第六计算式计算单位体积天然气燃烧产生的烟气中二氧化碳的容积数，所述第六计算式为：
式中，
通过第七计算式计算单位体积天然气燃烧产生的烟气中水蒸汽的容积数，所述第七计算式为：
式中，
通过第八计算式计算单位体积天然气燃烧产生的烟气中二氧化硫的容积数，所述第八计算式为：
式中，
通过第九计算式计算单位体积天然气燃烧产生的烟气中剩余氧气的容积数，所述第九计算式为：
式中，
通过第十计算式计算单位体积天然气燃烧产生的烟气中氩气的容积数，所述第十计算式为：
mar＝α×vair×τar/100+τar-gas/100
式中，mar为单位体积天然气燃烧产生的烟气中氩气的容积数，m3/s；τar为环境空气中氩气的容积百分数，％；τar-gas为天然气中氩气的容积百分数，％；
通过第十一计算式计算单位体积天然气燃烧生成的烟气总容积数，所述第十一计算式为：
式中，mfg为单位体积天然气燃烧生成的烟气总容积数，m3/s；
通过第十二计算式计算单位体积天然气燃烧生成的烟气各成分体积百分数，所述第十二计算式为：
式中，γi为单位体积天然气燃烧生成的烟气各成分体积百分数，％；其中，i为n2、co2、h2o、so2、02、ar；
根据直燃型溴化锂制冷机进口及出口的烟气温度及成分，确定其进口及出口的烟气的比热；
输入直燃型溴化锂制冷机的燃料流量、冷空气温度，计算冷空气焓值；
根据迭代法求解直燃型溴化锂制冷机高压发生器对数平均温差δtm，具体步骤如下：
1)设定一个直燃型溴化锂制冷机的排烟出口温度tf-out，计算排烟比热；
判断直燃型溴化锂制冷机是否有空气预热装置，确定空气进入燃烧器的温度，进而确定进入直燃型溴化锂制冷机的燃烧器的空气比热；
根据第十三计算式计算空气预热器的平均温差，所述第十三计算式为：
δtm1＝(tpy-tair-in)-0.35×min{(tpy-tpy-preh)，(tair-preh-tair-in)}-0.65×max{(tpy-tpy-preh)，(tair-preh-tair-in)}
式中，δtm1为空气预热器的平均温差，℃；tpy为直燃型溴化锂制冷机的排烟温度，℃；tpy-preh为直燃型溴化锂制冷机空气预热器出口处的排烟温度，℃；tair-in为直燃型溴化锂制冷机燃烧器进口空气温度，℃；tair-preh为直燃型溴化锂制冷机空气有预热器时燃烧器进口空气温度，℃；
确定加热燃气损失的附加系数，确定冷媒水的进、出口温度及冷却水进、出口温度；
其次，制冷机设计参数的选定还包括如下步骤：
设定直燃型溴化锂制冷机流程结构，是串联还是并联；
设定直燃型溴化锂制冷机吸收器与冷凝器的连接方式，是串联还是并联；
设定冷却水总温升，并对吸收器出口和冷凝器出口的冷却水温升进行分配；
设定吸收器出口处稀溶液与吸收器冷却水温差；
设定溴化锂溶液冷凝温度与冷却水出口处温度的差值；
设定溴化锂溶液蒸发温度与冷媒水出口处温度的差值；
设定吸收器压损及再循环倍数；
设定溴化锂溶液从吸收器中流到低压发生器后浓度变化量；
设定溴化锂溶液在高压发生器中浓度变化量；
设定凝结水离开凝结水换热器的温度；
设定蒸发器的再循环倍数。
根据第十四计算式计算蒸发器出口处冷剂蒸气的温度，所述第十四计算式为：
t0＝tx′-δt
式中，t0为蒸发器出口处冷剂蒸气的温度，℃；tx′为冷媒水出口温度，℃；δt为蒸发温度与冷媒水出口温度的差值，℃；
2)根据蒸发器出口处冷剂蒸气的温度确定其压力p0及焓值h0；
3)根据第十五计算式计算从蒸发器流入到吸收器内的饱和水的压力，所述第十五计算式为：
p1＝p0-δp0
式中，p1为从蒸发器流入到吸收器内的饱和水的压力，mpa；δp0为吸收器进口处与出口处稀溶液的压力差，mpa；
4)根据从蒸发器流入到吸收器内的饱和水的压力确定其温度t1及焓值h1；
5)根据第十六计算式计算从蒸发器流入到吸收器内的蒸汽压力，所述第十六计算式为：
p1a＝p1
式中，p1a为从蒸发器流入到吸收器内的蒸汽压力，mpa；
6)根据从蒸发器流入到吸收器内的蒸汽压力确定其温度t1a及焓值h1a；
7)根据第十七计算式计算吸收器出口处稀溶液的温度，所述第十七计算式为：
t2＝tw+δtw
式中，t2为吸收器出口处稀溶液的温度，℃；tw为吸收器冷却水进口温度，℃；δtw为吸收器出口处稀溶液与吸收器冷却水温差，℃；
8)根据第十八计算式计算吸收器出口处稀溶液的压力，所述第十八计算式为：
p2＝p1
式中，p2为吸收器出口处稀溶液的压力，mpa；
9)根据第十九计算式计算吸收器出口处稀溶液的浓度，所述第十九计算式为：
式中，ξ2为吸收器出口处稀溶液的浓度，％；a、b、c、d为吸收器出口处稀溶液浓度的计算系数，其取值根据吸收器出口处稀溶液压力下的饱和温度计算得到；
10)根据第二十计算式计算吸收器出口处稀溶液的焓值，所述第二十计算式为
式中，h2为吸收器出口处稀溶液的焓值，kj/kg；e、f、g为吸收器出口处稀溶液焓值的计算系数，其取值根据吸收器出口处稀溶液的温度和浓度计算得到；
11)根据第二十一计算式计算冷凝器出口处冷剂水的温度，所述第二十一计算式为：
t3＝tcw+δtk
式中，t3为冷凝器出口处冷剂水的温度，℃；tcw为冷却水出口温度，℃；δtk为冷凝温度与冷却水出口温度的差值，℃；
12)根据冷凝器出口处冷剂水的温度确定其压力p3及焓值h3；
13)根据第二十二计算式计算低压发生器产生的水蒸汽的压力，所述第二十二计算式为：
p3a＝p3
式中，p3a为低压发生器产生的水蒸汽的压力，mpa；
14)根据低压发生器产生的水蒸汽的压力，确定其焓值h3a；
15)根据第二十三计算式计算高压发生器内的饱和水的压力，所述第二十三计算式为：
p3b＝pr
式中，p3b为高压发生器内的饱和水的压力，mpa；pr为为高压发生器的压力，mpa；
16)根据高压发生器内的饱和水的压力确定其温度t3b及焓值h3b；
17)根据第二十四计算式计算高压发生器内水蒸汽的压力，所述第二十四计算式为：
p3c＝pr
式中，p3c为高压发生器内水蒸汽的压力，mpa；
18)根据高压发生器内水蒸汽的压力确定其焓值h3c；
19)根据第二十五计算式计算低压发生器出口处浓溶液的压力，所述第二十五计算式为：
p4＝p3
式中，p4为低压发生器出口处浓溶液的压力，mpa；
20)根据第二十六计算式计算低压发生器出口处浓溶液的浓度，所述第二十六计算式为：
ξ4＝ξ2+δξ1(直燃吸收式溴化锂制冷机流程结构为并联流程时)
ξ4＝ξ5+δξ1(直燃吸收式溴化锂制冷机流程结构为串联流程时)
式中，ξ4为低压发生器出口处浓溶液的浓度，％；δξ1为溴化锂溶液从吸收器中流到低压发生器后浓度变化量，％；ξ5为高温热交换器出口处的浓溶液(饱和溶液)的浓度，％；
21)根据第二十七计算式和第二十计算式分别计算低压发生器出口处浓溶液的温度t4和焓值h4；
22)根据第二十八计算式计算高温热交换器出口处的浓溶液(饱和溶液)的压力，所述第二十八计算式为：
p5＝p3(直燃吸收式溴化锂制冷机流程结构为串联流程时)
式中，p5为高温热交换器出口处的浓溶液(饱和溶液)的压力，mpa；
23)根据第二十九计算式计算高温热交换器出口处的浓溶液(饱和溶液)的浓度，所述第二十九算式为：
ξ5＝ξ12(直燃吸收式溴化锂制冷机流程结构为串联流程时)
式中，ξ12为高压发生器出口处浓溶液的浓度，％；
24)根据第十九计算式和第二十计算式分别计算高温热交换器出口处的浓溶液(饱和溶液)的温度t5和焓值h5(直燃吸收式溴化锂制冷机流程结构为串联流程时)；
25)根据第三十计算式计算低温热交换器的稀溶液的压力，所述第三十计算式为：
p7＝p3
式中，p7为低温热交换器的稀溶液的压力，mpa；
26)根据第三十一计算式计算低温热交换器的稀溶液的浓度，所述第三十一计算式为：
ξ7＝ξ2
式中，ξ7为低温热交换器的稀溶液的浓度，％；
27)根据第十九计算式计算低温热交换器的稀溶液的温度t7；
28)根据第三十二计算式计算低温热交换器的稀溶液的焓值，所述第三十二计算式为：
式中，h7为低温热交换器的稀溶液的焓值，kj/kg；a1为低位发生器循环倍率；h8为低温热交换器出口处的浓溶液(过冷溶液)的焓值，kj/kg；
29)根据第三十三计算式计算低温热交换器出口处的浓溶液(过冷溶液)的压力，所述第三十三计算式为：
p8＝p3
式中，p8为低温热交换器出口处的浓溶液(过冷溶液)的压力，mpa；
30)根据第三十四计算式计算低温热交换器出口处的浓溶液(过冷溶液)的浓度，所述第三十四计算式为：
ξ8＝ξ4
式中，ξ8为低温热交换器出口处的浓溶液(过冷溶液)的浓度，％；
31)根据第三十五计算式计算低温热交换器出口处的浓溶液(过冷溶液)的温度，所述第三十五计算式为：
t8＝t2+δth
式中，t8为低温热交换器出口处的浓溶液(过冷溶液)的温度，℃；δth为低温热交换器出口处的浓溶液与入口处的温差，℃；
32)根据第二十计算式计算低温热交换器出口处的浓溶液(过冷溶液)的焓值h8；
33)根据第三十六计算式计算吸收器的喷淋溶液的压力，所述第三十六计算式为：
p9＝p2
式中，p9为吸收器的喷淋溶液的压力，mpa；
34)根据第三十七计算式计算吸收器的喷淋溶液的浓度，所述第三十七计算式为：
式中，ξ9为吸收器的喷淋溶液的浓度，％；ξ13为高温热交换器出口处的浓溶液(过冷溶液)的浓度，％；qmh为单位时间进入高位发生器的稀溶液质量流量，kg/s；qmdh为高压发生器中产生的水蒸汽量，kg/s；qml为单位时间进入低位发生器的稀溶液质量流量，kg/s；qmdl为低压发生器中产生的水蒸汽量，kg/s；qmd为冷剂水质量流量，kg/s；f为吸收器的再循环倍数；
35)根据第十九计算式计算吸收器的喷淋溶液的温度t9；
36)根据第三十八计算式计算吸收器的喷淋溶液的焓值，所述第三十八计算式为：
式中，h9为吸收器的喷淋溶液的焓值，kj/kg；h13为高温热交换器出口处的浓溶液(过冷溶液)的焓值，kj/kg；
37)根据第三十九计算式计算高温热交换器的稀溶液的压力，所述第三十九计算式为：
p10＝pr
式中，p10为高温热交换器的稀溶液的压力，mpa；
38)根据第四十计算式计算高温热交换器的稀溶液的浓度，所述第四十计算式为：
ξ10＝ξ2
式中，ξ10为高温热交换器的稀溶液的浓度，％；
39)根据第十九计算式计算高温热交换器的稀溶液的温度t10；
40)根据第四十一计算式计算高温热交换器的稀溶液的焓值，所述第四十一计算式为：
式中，h10为高温热交换器的稀溶液的焓值，kj/kg；ah为高压发生器循环倍率；h12为高压发生器出口处浓溶液的焓值，kj/kg；
41)根据第四十二计算式计算高压发生器出口处浓溶液的压力，所述第四十二计算式为：
p12＝pr
式中，p12为高压发生器出口处浓溶液的压力，mpa；
42)根据第四十三计算式计算高压发生器出口处浓溶液的浓度，所述第四十三计算式为：
ξ12＝ξ2+δξ2
式中，ξ12为高压发生器出口处浓溶液的浓度，％；δξ2为溴化锂溶液从吸收器中流到高压发生器后浓度变化量，％；
43)根据第十九计算式和第二十计算式分别计算高压发生器出口处浓溶液的温度t12和焓值h12；
44)根据第四十四计算式计算高温热交换器出口处的浓溶液(过冷溶液)的压力，所述第四十四计算式为：
p13＝pr(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，p13为高温热交换器出口处的浓溶液(过冷溶液)的压力，mpa；
由于高温热交换器入口处浓溶液的温度较高，所以在热交换器中温降大于低温热交换器中浓溶液的温升，通常浓溶液的出口温度在60～70℃范围内(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
45)根据第四十五计算式计算高温热交换器出口处的浓溶液(过冷溶液)的浓度，所述第四十五计算式为：
ξ13＝ξ12(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，ξ13为高温热交换器出口处的浓溶液(过冷溶液)的浓度，％；
46)根据第二十计算式计算高温热交换器出口处的浓溶液(过冷溶液)的焓值h10(仅直燃吸收式溴化锂制冷机流程结构为串联流程时)；
47)根据第四十六计算式计算低压发生器循环倍率，所述第四十六计算式为：
式中，al为低压发生器循环倍率；
48)根据第四十七计算式计算高压发生器循环倍率，所述第四十七计算式为：
式中，ah为高压发生器循环倍率；
根据迭代法求解直燃型溴化锂制冷机的加热工作热水流量qmv2″，具体步骤如下：
(1)设定一个直燃型溴化锂制冷机的制冷量q0；
(2)根据第四十八计算式计算制冷机中冷剂水的流量，所述第四十八计算式为：
式中，qmd为制冷机中冷剂水的流量，kg/s；q0为直燃型溴化锂制冷机的制冷量，kw；
(3)根据第四十九计算式计算低压发生器中产生的水蒸汽量，所述第四十九计算式为：
式中，qmdl为低压发生器中产生的水蒸汽量，kg/s；
(4)根据第五十计算式计算高压发生器中产生的水蒸汽量，所述第五十计算式为：
qmdh＝qmd-qmdl
式中，qmdh为高压发生器中产生的水蒸汽量，kg/s；
(5)根据第五十一计算式计算单位时间进入高压发生器的稀溶液质量，所述第五十一计算式为：
qmh＝ah×qmdh
式中，qmh为单位时间进入高压发生器的稀溶液质量，kg/s；
(6)根据第五十二计算式计算单位时间进入低压发生器的稀溶液质量，所述第五十二计算式为：
qml＝al×qmdl
式中，qml为单位时间进入高低压发生器的稀溶液质量，kg/s；
(7)根据第五十三计算式计算高压发生器中产生1kg水蒸气时，在吸收器内产生的热负荷，所述第五十三计算式为：
qah＝(ah-1)×h13-ah×h2+h1a(仅直燃吸收式溴化锂制冷机流程结构为并联流程时)；
(8)根据第五十四计算式计算低压发生器中产生1kg水蒸气时在吸收器内产生的热负荷，所述第五十四计算式为：
qal＝(al-1)×h8-al×h2+h1a
式中，qal为低压发生器中产生1kg水蒸气时在吸收器内产生的热负荷，kj/kg；
(9)根据第五十五计算式计算高压发生器中产生1kg水蒸气时的热负荷，所述第五十五计算式为：
qgh＝(ah-1)×h12-ah×h10+h3c
式中，qgh为高压发生器中产生1kg水蒸气时的热负荷，kj/kg；
(10)根据第五十六计算式计算低压发生器中产生1kg水蒸气时的热负荷，所述第五十六计算式为：
式中，qlh为低压发生器中产生1kg水蒸气时的热负荷，kj/kg；
(11)根据第五十七计算式计算吸收器热负荷，所述第五十七计算式为：
qa＝qmdh×qh+qmdl×ql(直燃吸收式溴化锂制冷机流程结构为并联流程时)；
qa＝qmdl×ql(直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，qa为吸收器热负荷，kw；
(12)根据第五十八计算式计算高压发生器热负荷，所述第五十八计算式为：
qgh＝qmdh×qgh
式中，qgh为高压发生器热负荷，kw；
(13)根据第五十九计算式计算低压发生器热负荷，所述第五十九计算式为：
qgl＝qmdl×qgl
式中，qgl为低压发生器热负荷，kw；
(14)根据第六十计算式计算冷凝器热负荷，所述第六十计算式为：
qk＝qmdh×(h3c-h3)+qmdl×(h3a-h3)(直燃吸收式溴化锂制冷机流程结构为并联流程时)；
qk＝qmdh×(h3c-h3)-qgh+qmdl×(h3a-h3)(直燃吸收式溴化锂制冷机流程结构为串联流程时)；
式中，qk为冷凝器热负荷，kw；
(15)蒸发器热负荷等于溴化锂制冷机的制冷量q0；
(16)根据第六十一计算式计算高压溶液热交换器的总热负荷，所述第六十一计算式为：
qexh＝qmh×(h10-h2)
式中，qexh为高压溶液热交换器的总热负荷，kw；
(17)根据第六十二计算式计算低压溶液热交换器的总热负荷，所述第六十二计算式为：
qexl＝qml×(h7-h2)
式中，qexl为低压溶液热交换器的总热负荷，kw；
(18)凝结水热交换器的总热负荷qexv等于低温水总放热量荷qv；
(19)根据第六十三计算式计算直燃型溴化锂制冷机的热力系数，所述第六十三计算式为：
式中，ζ为直燃型溴化锂制冷机的热力系数；
(20)根据第六十四计算式计算直燃型溴化锂制冷机的加热工作热水流量，所述第六十四计算式为：
式中，qmv2为直燃型溴化锂制冷机的加热工作热水流量，kg/s；
(21)将qmv2″与qmv2做差，得出二者的差值；
若差值在预定的误差范围内，则假设的q0为直燃型溴化锂制冷机的制冷量；
若差值超出预定的误差范围内，则将qmv2″和qmv2的平均值作为新的qmv2″，重新执行上述(1)～(21)的计算，直到qmv2″与qmv2的差值满足设定的误差范围；
根据第六十五计算式计算直燃型溴化锂制冷机的冷媒水流量，所述第六十五计算式为：
式中，qv0为直燃型溴化锂制冷机的冷媒水流量，kg/s；tx″为冷媒水进口温度，℃；tx′为冷媒水出口温度，℃；
根据第六十六计算式计算直燃型溴化锂制冷机的冷却水泵的流量，所述第六十六计算式为：
式中，
根据第六十七计算式计算直燃型溴化锂制冷机的蒸发器泵的流量，所述第六十七计算式为：
式中，qvd为直燃型溴化锂制冷机的蒸发器泵的流量，m3/h；f0为蒸发器的再循环倍数；
根据第六十八计算式计算直燃型溴化锂制冷机进入发生器的溶液密度，所述第六十八计算式为：
式中，ρa为直燃型溴化锂制冷机的吸收器溶液密度，kg/l；a0、a1、a2、a3、a4、a5、a6为直燃型溴化锂制冷机的吸收器溶液密度的计算系数；
根据第六十九计算式计算直燃型溴化锂制冷机吸收器的溶液密度，所述第六十九计算式为：
式中，ρ9为直燃型溴化锂制冷机的吸收器溶液密度，kg/l；a0、a1、a2、a3、a4、a5、a6为直燃型溴化锂制冷机的吸收器溶液密度的计算系数；
根据第七十计算式计算直燃型溴化锂制冷机发生器泵的流量，所述第七十计算式为：
式中，qvg为直燃型溴化锂制冷机发生器泵的流量，m3/h；
根据第七十一计算式计算直燃型溴化锂制冷机吸收器泵的流量，所述第七十一计算式为：
式中，qva为直燃型溴化锂制冷机吸收器泵的流量，m3/h；f为吸收器的再循环倍数；
根据第七十二计算式计算直燃型溴化锂制冷机的高压发生器的传热面积，所述第七十二计算式为：
式中，fgh为直燃型溴化锂制冷机的高压发生器的传热面积，m2；kgh为直燃型溴化锂制冷机高压发生器的传热系数，w/(m2·℃)；tf-in为加热工作热水进口温度，℃；bg1为发生器换热过程中，冷流体换热系数，无量纲；
根据第七十三计算式计算直燃型溴化锂制冷机的低压发生器的传热面积，所述第七十三计算式为：
式中，fgl为直燃型溴化锂制冷机的低压发生器的传热面积，m2；kgl为直燃型溴化锂制冷机低压发生器的传热系数，w/(m2·℃)；
根据第七十四计算式计算直燃型溴化锂制冷机的冷凝器的传热面积，所述第七十四计算式为：
式中，fk为直燃型溴化锂制冷机的冷凝器的传热面积，m2；kk为直燃型溴化锂制冷机冷凝器的传热系数，w/(m2·℃)；bk1为冷凝器换热过程中，冷流体换热系数，无量纲；
根据第七十五计算式计算直燃型溴化锂制冷机的吸收器的传热面积，所述第七十五计算式为：
式中，fa为直燃型溴化锂制冷机的吸收器的传热面积，m2；ka为直燃型溴化锂制冷机吸收器的传热系数，w/(m2·℃)；aa1为吸收器换热过程中，热流体换热系数，无量纲；ba1为吸收器换热过程中，冷流体换热系数，无量纲；
根据第七十六计算式计算直燃型溴化锂制冷机的蒸发器的传热面积，所述第七十六计算式为：
式中，f0为直燃型溴化锂制冷机的蒸发器的传热面积，m2；k0为直燃型溴化锂制冷机蒸发器的传热系数，w/(m2·℃)；b01为蒸发器换热过程中，冷流体换热系数，无量纲；
根据第七十七计算式计算直燃型溴化锂制冷机的高温溶液热交换器的传热面积，所述第七十七计算式为：
式中，fexh为直燃型溴化锂制冷机的高温溶液热交换器的传热面积，m2；kexh为直燃型溴化锂制冷机高温溶液热交换器的传热系数，w/(m2·℃)；aex1为高温溶液热交换器换热过程中，热流体换热系数，无量纲；bex1为高温溶液热交换器换热过程中，冷流体换热系数，无量纲；
根据第七十八计算式计算直燃型溴化锂制冷机低温溶液热交换器的传热面积，所述第七十八计算式为：
式中，fexl为直燃型溴化锂制冷机的低温溶液热交换器的传热面积，m2；kexl为直燃型溴化锂制冷机低温溶液热交换器的传热系数，w/(m2·℃)；aex2为低温溶液热交换器换热过程中，热流体换热系数，无量纲；bex2为低温溶液热交换器换热过程中，冷流体换热系数，无量纲；
2)根据第七十九计算式计算直燃型溴化锂制冷机的高压发生器的对数平均温差，所述第七十九计算式为：
式中，δtm为直燃型溴化锂制冷机的高压发生器的对数平均温差，℃；tf-out为加热工作热水出口温度，℃；
3)根据第八十计算式计算直燃型溴化锂制冷机的高压发生器的核算对数平均温差，所述第八十计算式为：
式中，δt′m为直燃型溴化锂制冷机的高压发生器的核算对数平均温差，℃；
4)将δtm与δt′m做差，得出二者的差值；
若差值在预定的误差范围内，则假设的tf-out为直燃型溴化锂制冷机工作热水的出口温度；
若差值超出预定的误差范围内，则重新假设tf-out，重新执行上述1)～4)的计算，直到δtm与δt′m的差值满足设定的误差范围；
以上，仅为本发明的较佳实施案例，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替代，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求所界定的保护范围为准。