
2. 中国科学院 广州能源研究所,广东 广州 510640






2. Guangzhou Institute of Energy Conversion, Chinese Academy of Sciences, Guangzhou 510640, China
如今,地热能的利用已经逐步被世界各国所关注.“世界地热大会WGC2010”统计,截止2009年底,全球已有27个国家实现了地热发电,装机容量和年发电量分别为10 715 MW和672.5亿kWh;地热直接利用装机容量达50 583 MW,年产能达4.38×1014 kJ[1-2].我国地热能分布以中低温为主,地热能较多地被用来采暖、洗浴和热泵等[3].发电方面,国内目前仅存且运行的地热能发电站只有西藏羊八井和广东丰顺两个.丰顺闪蒸地热电站是将91℃地热水抽入闪蒸器闪蒸出蒸汽,从而推动汽轮机发电.调研发现,电站仍有大量未闪蒸的地热水未被进一步利用而直接排弃.而这些直接排放的地热水的温度在68~75℃之间,利用温差较小,利用率较低.
地热发电与地热直接利用相结合的地热梯级综合利用系统能够有效地提高地热资源的利用率,进而提高地热利用系统的热力性能.近年来,国际上有学者在地热的综合利用方面进行研究,如T.A.H. Ratlamwala等[4]将地热水串联通过双级闪蒸发电循环、四效氨水吸收式制冷循环、尾水换热循环和制氢循环,形成地热的梯级综合利用.建模并分析了地热水温度、压力及环境温度对系统性能的影响;Oguz Arslan等[5]对土耳其某地热资源进行了梯级综合利用,地热水依次完成有机朗肯循环发电、居民采暖、温室养殖和温泉洗浴;Duccio Tempesti等[6]研究了一种太阳能地热耦合的有机朗肯循环发电系统,并对不同月份下系统发电效率进行了分析;Can Coskun等[7]对供冷季、供热季不同地热综合利用系统的形式做了热力学的分析,并考察了环境温度对各系统效率的影响;Mehmet Kanoglu等[8]对地热发电、制冷和采暖这3种可能的组合方案分别做了热力学和经济学评价.
因此,根据丰顺夏热冬暖地区的气候特征,利用双级溴化锂吸收式制冷循环(TSARS)在低品位热能回收方面的优势,将地热闪蒸电站未被利用的地热尾水继续级联到TSARS,构成地热梯级综合利用系统,一方面提高地热水的利用率,另一方面则产生冷量以满足制冷需求.但是与传统的电制冷相比,TSARS设备复杂,投资高,且循环COP低.虽然近年来
综上所述,本文提出将TSARS应用到电站地热尾水的回收,建立地热水发电-制冷梯级综合利用系统.以电站运行数据为基准,获得地热水的
广东丰顺地热电站采用的是单级闪蒸发电循环,其流程如图 1所示.从生产井(PW)出来的地热水通过生产井泵(P)输运到电站闪蒸器(FT),进而通过控制闪蒸压力Pf,使闪蒸器内的地热水闪蒸产生蒸汽,蒸汽进入汽轮机(ST)做功而推动发电机(G)发电,经过汽轮机做功后的乏汽进入凝汽器(CON-P)被冷却水(CW)直接冷却.据调研发现,闪蒸器内剩余地热水被直接排弃,流量达60 kg/s,最高温度为75℃.因此,当环境温度为25℃时,该闪蒸电站的地热水利用率仅为24.3%,即热水利用的能量只有其最大做功能力的24.3%,造成了地热能的大量浪费.图 2为2011年1~10月电站的发电功率,平均发电功率约为240 kW,发电效率仅为4.48%.
![]() |
图 1 丰顺地热电站发电流程简图 Figure 1 Flow chart of geothermal power plant in Fengshun |
![]() |
图 2 2011年电站各月发电功率 Figure 2 The generated output of the Plant in 2011 |
基于地热电站较低的地热水利用率,结合地热电站循环排放尾水温位的特点,提出将闪蒸发电后电站排弃的地热尾水级联到TSARS,形成地热梯级综合利用系统,从而提高地热水的利用率.溴化锂吸收式制冷技术在低品位能源利用和工业余热回收方面有较大潜力,而TSARS对热源的温度要求比单效溴化锂吸收式制冷系统更低[12],60℃以上的热水即可驱动,进而可利用其对70℃左右的地热尾水加以利用.TSARS流程如图 3所示,而地热尾水级联到TSARS有串联供入和并联供入两种可能的模式,如图 4所示.
![]() |
图 3 双级溴化锂吸收式制冷循环流程图 Figure 3 Flow chart of the TSARS |
![]() |
图 4 地热尾水供入TSARS的两种模式 Figure 4 The ways of geothermal waste water entering into TSARS |

$ \dot C = c \cdot E = c \cdot \dot m \cdot e, $ | (1) |
$ e = \left( {h - {h_0}} \right) - {T_0} \cdot \left( {s - {s_0}} \right), $ | (2) |
$ {{\dot E}_D} = {{\dot E}_F} - {{\dot E}_P} - {{\dot E}_L}, $ | (3) |
其中,e为单位质量
制冷系统计算模型是建立在如下假设的前提下进行的:
(1) 系统稳态运行且所有的换热器与外界无热交换;
(2) 制冷剂水(水蒸气)在冷凝器进口和蒸发器的出口均为饱和态;
(3) 不考虑系统由于摩擦而带来的管道和换热器的工质压力损失;
(4) 不考虑溴化锂水溶液的动能、势能和化学能;
(5) 溴化锂水溶液在系统各点处于平衡;
(6) 地热水物性采用纯水物性;
(7) 为了保证溴化锂溶液不结晶,进入高低压节流阀的溴化锂浓溶液的温度保证在当前压力下结晶温度8℃以上.
3.1 子系统划分将双级溴化锂制冷系统每一个热力过程划分成一个子系统,即把每个设备当做一个子系统.故每个子系统的“燃料-产品-损失”定义如表 1所示.
![]() |
表 1 各子系统“燃料-产品-损失”的定义 Table 1 The definition of "F-P-L" of each subsystem |

根据
$ \sum {{{\dot C}_{{\rm{in}},k}}} + {{\dot Z}_k} = \sum {{{\dot C}_{{\rm{out}},k}}} , $ | (4) |
其中,
$ {{\dot Z}_k} = {Z_k} \cdot o. $ | (5) |
其中,资金年度化系数[14]
$ o = \frac{{m{{\left( {1 + m} \right)}^n}}}{{{{\left( {1 + m} \right)}^n} - 1}} \cdot \frac{1}{{n \cdot t \cdot 3600}}, $ | (6) |
式(6)中,m为年利率;n为系统运行寿命,年;t为系统年运行小时数,h.
建立各子系统
![]() |
表 2 各子系统![]() |
双级溴化锂吸收式制冷系统的投资包括换热器投资、泵(电机)和节流阀的投资.其中高(低)压发生器、高(低)压吸收器、高(低)压溶液交换器、冷凝器和蒸发器的投资利用式(7)计算,节流阀、泵(电机)的投资利用式(8)计算[15](把节流阀近似为泵的反向过程).
$ {Z_i} = {f_i} \cdot {\left( {\frac{{{A_i}}}{{{A_r}}}} \right)^{0.6}}, $ | (7) |
$ {Z_j} = {f_j} \cdot {\left( {\frac{{{W_i}}}{{{W_r}}}} \right)^{0.5}}, $ | (8) |
其中,换热器的换热面积和泵(电机)功耗分别按照式(9)、(10)确定.
$ {A_i} = \frac{{{Q_i}}}{{{\rm{LMT}}{{\rm{D}}_i} \cdot {K_i}}}, $ | (9) |
$ {W_i} = \frac{{{{\dot m}_i} \cdot g \cdot {H_i}}}{{1000 \cdot {\beta _i}}}, $ | (10) |
其中,Qi为换热器换热量,kW;LMTDi为换热器对数换热温差,℃;Ki为换热器平均换热系数,kW/m2;g为重力加速度,m/s2;Hi为泵的扬程,m;βi为泵的效率,fi、fj、Ar、Wr为修正参数,参考文献[15].
3.3 地热水成本的确定据
$ {{\dot C}_{{\rm{water}}}} = o \cdot {C_{{\rm{water}}}}, $ | (11) |
其中,电站成本Cwater按式(12)计算,各含量取值见表 3.
$ {C_{{\rm{water}}}} = \left( {{C_{{\rm{exp}}}} + {C_{{\rm{per}}}} + {C_d}} \right) \cdot {W_{{\rm{design}}}}, $ | (12) |
![]() |
表 3 地热电站各项目费用[15] ($·kW-1) Table 3 The fees of the geothermal power plant |
所以,得
$ {C_{{\rm{water}}}} = {C_{22}} = \frac{{\dot C}}{{\dot E}} = \frac{{\left( {610 \sim 1265} \right) \cdot o \cdot {W_{{\rm{design}}}}}}{{\dot m \cdot {e_{22}}}}, $ | (13) |
其中,Wdesign为电站设计发电功率,kW; Cexp为勘测费; Cper为许可费; Cd为钻井费.
3.4 系统评价分析指标$ 循环能效比\;{\rm{COP}} = \frac{{{Q_{{\rm{gain}}}}}}{{{Q_{{\rm{pay}}}}}}. $ | (14) |
$ 循环㶲效率\;{\eta _{{\rm{ex}}}} = \frac{{{E_{{\rm{gain}}}}}}{{{E_{{\rm{pay}}}}}}. $ | (15) |
$ 㶲损失率\;{\xi _i} = \frac{{{{\dot C}_{{\rm{L,i}}}} + {{\dot C}_{{\rm{D,i}}}}}}{{{{\dot C}_{{\rm{l}},{\rm{total}}}} + {{\dot C}_{{\rm{D}},{\rm{total}}}}.}} $ | (16) |
$ 㶲经济系数\;f = \frac{{{Z_i}}}{{{Z_i} + {{\dot C}_{L,i}} + {{\dot C}_{D,i}}}}. $ | (17) |
$ 产品单位㶲成本\;{C_p} = \frac{{{{\dot C}_P}}}{{{{\dot E}_P}}} = \frac{{{{\dot C}_{28}} - {{\dot C}_{29}}}}{{{{\dot E}_{28}} - {{\dot E}_{29}}}}. $ | (18) |
表 4给出了TSARS的设计和计算参数,考虑到实际运行条件,设串联、并联供给模式的TSARS利用地热水的温差均为8℃.
![]() |
表 4 双级溴化锂吸收式制冷循环设计及计算参数 Table 4 The designed and calculated parameters of TSARS |
单纯从热力学出发,并联模式TSARS的COP和
![]() |
图 5 地热水成本对TSARS冷冻水![]() |
实地调研发现,由于受到电站运行参数和环境温度变化的影响,闪蒸电站的地热尾水的排放温度T22在68~75℃之间.从图 6得,随着地热水温度的增大,循环的COP不断提高,且COP增大的速率不断减小,在75℃时达到并联TSARS的COP的最大值,约0.417.而从图 7得,随着地热水温度的增大,循环的
![]() |
图 6 地热水进口温度对TSARS的COP的影响 Figure 6 Effect of the waste geothermal water's temperature on the COP of TSARS (T29=10℃,T24=28℃,c22=0.211 $/GJ) |
![]() |
图 7 地热水进口温度对TSARS ![]() |

表 5给出了地热水进口温度为71℃时,并联TSARS各状态点的参数.
![]() |
表 5 并联TSARS各点运行参数和![]() |
在地热水进口温度为71℃,冷却水进口温度为28℃,冷冻水出口温度为10℃,TSARS总输入
![]() |
图 8 并联TSARS各设备![]() ![]() |
![]() |
图 9 并联TSARS各设备![]() |
本文提出将TSARS级联到原闪蒸地热电站,形成地热水的梯级综合利用系统.从电站实际运行出发,确定地热水的成本.采用
(1) 地热水并联级联到TSARS较串联模式具有较好的
(2) 地热水串联TSARS的COP和
(3) 地热水并联TSARS中换热器的总
[1] |
Lund J W, Freeston D H, Boyd T L. Direct utilization of geothermal energy 2010 worldwide review[J].
Geothermics, 2011, 40(03): 159-180.
DOI: 10.1016/j.geothermics.2011.07.004. |
[2] |
Arslan O, Kose R. Exergoeconomic optimization of integrated geothermal system in Simav, Kutahya[J].
Energy Conversion and Management, 2010, 51(4): 663-676.
DOI: 10.1016/j.enconman.2009.11.010. |
[3] |
Luo C, Huang L C, Gong Y L. Thermodynamic comparison of different types of geothermal power plant systems and case studies in China[J].
Renewable Energy, 2012(48): 155-160.
|
[4] |
Ratlamwala T A H, Dincer I, Gadalla M A. Performance analysis of a novel integrated geothermal-based system for multi-generation applications[J].
Applied Thermal Engineering, 2012, 40(04): 71-79.
|
[5] |
Arslan O, Kose R. Exergoeconomic optimization of integrated geothermal system in Simav, Kutahya[J].
Energy Conversion and Management, 2010, 51(04): 663-676.
DOI: 10.1016/j.enconman.2009.11.010. |
[6] |
Tempesti D. Thermodynamic analysis of two micro CHP systems operating with geothermal and solar energy[J].
Applied Energy, 2012, 48(02): 1-9.
|
[7] |
Coskun C, Oktay Z, Dincer I. Thermodynamic analyses and case studies of geothermal based multi-generation systems[J].
Energy, 2012, 32(06): 71-80.
|
[8] |
Kanoglu M, Cengel Y A. Economic evaluation of geothermal power generation, heating, and cooling[J].
Energy, 1999, 24(06): 501-509.
DOI: 10.1016/S0360-5442(99)00016-X. |
[9] |
Sencan A, Yakut K A, Kalogirou S A. Exergy analysis of lithium bromide/water absorption systems[J].
Renewable Energy, 2005, 30(5): 645-657.
DOI: 10.1016/j.renene.2004.07.006. |
[10] |
Misra R D, Sahoo P K, Gupta A. Thermoeconomic evaluation and optimization of a double-effect H2O/LiBr vapour-absorption refrigeration system[J].
International Journal of Refrigeration, 2005, 28(3): 331-343.
DOI: 10.1016/j.ijrefrig.2004.09.006. |
[11] |
Misra R D, Sahoo P K, Gupta A. Thermoeconomic evaluation and optimization of an aqua-ammonia vapour-absorption refrigeration system[J].
International Journal of Refrigeration, 2006, 29(1): 47-59.
DOI: 10.1016/j.ijrefrig.2005.05.015. |
[12] |
Ma W B, Deng S M. Theoretical analysis of low-temperature hot source driven two-stage LiBr/H2O absorption refrigeration system[J].
International Journal of Refrigeration, 1996, 19: 141-146.
DOI: 10.1016/0140-7007(95)00054-2. |
[13] |
Lozano M A, Valero A. Theory of the exergetic cost[J].
Energy, 1993, 18(9): 939-960.
DOI: 10.1016/0360-5442(93)90006-Y. |
[14] |
Tsatsaronis G, Javier P. Exergoeconomics evaluation and optimization of energy systems—application to the CGAM problem[J].
Energy, 1994, 19(3): 287-321.
DOI: 10.1016/0360-5442(94)90113-9. |
[15] |
Bejan A, Tsatsaronis G, Moran M. Thermal design and optimization[M]. New York: John Wiley and Sons Inc, 1996.
|