238重庆市某酒店水源热泵对长江水环境影响的模拟分析

1重庆市某酒店水源热泵对长江水环境影响的模拟分析重庆大学城市建设与环境工程学院吴浩王勇李文摘要是否具备合适的水源以及排水对水环境的影响有多大，是建造地表水源热泵工程前必须解决的两个问题。本文从水温和水质两方面研究了长江水作为地表水源热泵空调系统冷热源的可行性；对重庆某酒店水源热泵空调系统排水对长江水的影响进行了模拟分析，根据模拟结果得到结论：该工程从长江取、回水是可行的，不但满足工程需要，而且排水只在较小的区域内对江水水体有程度很小的影响，在该区域以外，江水水温都能得到恢复。最后本文展望了地表水源热泵空调系统在长江流域的应用前景。关键词地表水源热泵长江水环境CFD模拟1引言地表水源热泵就是利用地球表面的江、河、湖、海等地表水作为热泵机组的热源和热汇。当建筑物的周围有大量的地表水可以利用时，可通过水泵和输配管路将水体的热量传递给热泵机组或将热泵机组的热量释放到水体中。重庆境内江河纵横，水网密布，水资源十分丰富，全市年平均水资源总量超过5000亿m3，其中地表水资源占绝大部分，由长江、嘉陵江、乌江等流经重庆地区的入境水形成的地表水约4600亿m3，其中长江占80％以上，嘉陵江占9.9％，其它河流约占10％[1]。重庆地区丰富的地表水资源，为水源热泵技术在重庆地区的应用提供了得天独厚的自然资源条件。2工程概况该酒店位于重庆市江津区，建筑总面积为102628m2，地上39层，建筑高度为148.6米，为集商业、办公、酒店客房、休闲会所、会议为一体的超高层五星级酒店综合体，是城市重要的标志。酒店紧靠长江，故拟采用水源热泵空调系统，利用长江水作为低位冷热源通过热泵提升能量为建筑进行供冷和供热。酒店总冷负荷约为12000KW；总热负荷约为6000KW。3长江水利用的自然条件3.1水温参数由文献[1]的统计数据可知：长江水最低温度出现在1月份，为10.30C，最高温度出现在8月份，为24.70C。为进一步摸清长江江津段水温情况，课题组于07年1月份和8月份分别对长江江津段水温进行了测试，测试结果如下图1、图2所示：10.0010.2010.4010.6010.8011.0010:3012:0014:1515:20测试时间水温(℃)23.023.524.024.525.08:309:0010:3511:1014:30测试时间水温(℃)图1冬季水温测试情况(07-1-13)图2夏季水温测试情况(07-8-7)（07-8-7）2测试结果与文献[1]中统计数据非常接近。课题组于06年8月份和08年2月份分别对极端条件下的长江水温进行了测试：06年8月份最高水温为28.50C，08年2月最低水温为7.00C，根据美国制冷学会ARI320标准，即使在百年不遇的连续高温天气和五十年不遇的严寒天气情况下，长江水仍然适合水源热泵，长江水的这种温度特性使其成为水源热泵良好的冷热源。3.2水质参数水质的好坏直接关系到机组的运行效果和使用寿命，地表水的水质指标包括水的浊度、硬度以及藻类和微生物含量等。《采暖通风与空气调节设计规范》（GB50019—2003)中第7.3.3条和《地源热泵系统工程技术规范》（GB50366—2005）中第5.2.8条和第6.2.4条对采用地下水、地表水的水源热泵机组的水质指标进行了规定，表1是长江江津段的各项水质指标测试结果与允许值的对比：序号水质指标测试结果允许值备注1PH值6.42～8.216.5~8.5满足要求2Ca2+(mg/l)38.5～82.23200满足要求3矿化度(g/l)＜1.303满足要求4Cl-(mg/l)8.23～13.16100满足要求5SO42-(mg/l)24.12～32.20200满足要求6Fe2+(mg/l)0.06～0.201满足要求7H2S(mg/l)0.17～0.320.5满足要求8浊度(NTU)620～12706.5不满足要求表1长江江津段水质与水源热泵机组允许的水质要求比较由表1可以看出：除浊度外，其他各项水质指标均能满足机组对水质的要求。由此可见,利用长江水做水源热泵系统冷热源必须解决含沙量大的问题。本工程采用天然河底滤床反向渗滤取水方案，该方案利用长江河床的卵石层作为天然过滤器，通过渗透滤嘴在河底硐室内汇集江水，各个硐室内汇集的江水通过江底集水隧道、导井送至竖井，然后利用水泵直接抽往热泵机组使用。江水通过河床的卵石层过滤、渗透至酮室后，水中含沙量将大幅度降低，经简单处理后就可以直接进入机组，这就解决了江水含沙量大的问题。4水源热泵机组排水对江水影响的模拟分析本文采用CFD通用商业软件PHOENICS来对水源热泵从长江提取冷、热量后，回水对水体的影响范围、影响程度进行模拟分析。4.1江水水体“冷热收支”情况分析水体的“冷热收支”主要取决于以下几个部分：(1)水面热交换量；(2)水体与河底、河床的换热量；(3)水源热泵系统对水体的取热或排热量。4.1.1水面热交换量的计算[2][3][4][5]水面热交换包括辐射热流量、蒸发热流量和对流热流量三个方面。具体地，水表面层的热交换总量n（w/m2）可表示为：nsabec(1)3式中：s为水表面吸收的太阳短波辐射量（w/m2）；a为水表面吸收的大气长波辐射量（w/m2）；b为水面长波辐射（w/m2）；e水表面的净蒸发热通量（w/m2）；c为水表面的热对流通量（w/m2）。(1)水体吸收的太阳短波辐射太阳的日照总辐射I经过水面反射后，除表层吸收一部分外，其余部分透射入水体内部。因此水体吸收的太阳短波辐射量s（w/m2）的计算公式如下，)1(*Is(2)式中：I为平均太阳总辐射强度（w/m2），一般可引用现场或相邻主要气象台所测得的太阳辐射量数据，重庆地区的平均太阳辐射强度按照“中国建筑热环境分析专用气象数据集”里面沙坪坝气象站的数据来取值，夏天取168.2w/m2，冬天取49.9w/m2；为水表面对太阳短波辐射的反射率。太阳辐射在水体中的分布规律是随水深的增加呈指数规律衰减，即：11expzIHz(3)式中：z为高程为z处吸收的太阳辐射热；为水面表层对太阳辐射的吸收率；为太阳辐射在水体中的衰减系数m-1；H为水体的深度（m）。本文在模拟中根据水体的实际情况选定以下数值＝0.4；=0.1；=0.5m-1。(2)水表面吸收的大气长波辐射量a（w/m2）41273aaacaT(4)式中：a为水表面对大气的长波反射率；ac为大气的发射率；是Stefan－Boltzman常数，为-8245.6710/wmk；aT为水面2m处的空气温度(0C)，根据测试数据，aT取值为夏天取32.40C，冬天取7.40C。本模拟中取值a＝0.03；ac＝0.97。(3)水面长波辐射b（w/m2）水体吸收的大气长波辐射能量会向大气进行返回辐射，是水体热损失的很重要的一部分。基于水面为绝对黑体的假设，水面长波辐射量的计算公式如下，4273bwwT(5)式中：wT为水表面的温度(0C)，根据对重庆地区长江水温的测试数据，wT取值为夏天取24.40C，冬天取9.60C；w为水表面的长波发射率，由于水并非绝对黑体，取0.97；是Stefan－Boltzman常4数。(4)水表面的净蒸发热通量e（w/m2）估计蒸发量的方法很多，其中大多数是经验性的，水表面蒸发热通量的计算公式为，ezwafwee(6)式中：zw表示水面上z（m）处的风速，)(zwf是用水面上z（m）处风速表示的风函数（w/m2*mmHg），zfw=29.20.46zw；ae为水面上空气的蒸发压力(mmHg)；we为对应于水面温度wT(0C)的空气饱和蒸发压力,we=exp20.855278/273wT，重庆地区夏季、冬季室外平均风速按以下取值：夏季2.1m/s，冬季0.8m/s。(5)水表面的热对流通量c（w/m2）当气温不等于水面温度时，水汽交接面上会发生对流传热过程，其计算公式为，czwaafwTT(7)式中：a为经验系数，取0.47；)(zwf是用水面上方z（m）风速表示的风函数（w/m2*mmHg）；wT为水表面温度(0C)；aT为水表面上方2m处空气温度(0C)。4.1.2水体与河底、河床的换热量gE(w)根据J.R.Hull[6]得到的半经验公式，水底以及底水层壁面的热损失量gE(w)计算式为gwgEAPTT(8)式中：A为池底面积（m2）；wT为水体的温度；gT为土壤远边界的温度(重庆地区土壤的年平均温度为190C)；p为池底周长（m）；为经验系数（w/m2），与土壤的传热系数K(w/m2.k)成正比，与池底和土壤远边界之间的距离(m)成正比；为经验系数（w/m），包括池底热损失以及底水层壁面热损失的影响，与土壤的传热系数K成正比关系。4.1.3水源热泵系统对水体的取热或排热量制冷工况运行时，水源热泵系统排入到水体中的热量Q为：EEREERQQe1(w)(9)5式中:eQ为夏季空调冷负荷(w)；EER为制冷能效比。制热工况运行时，水源热泵系统从水体中吸取的热量Q为：COPCOPQQc1(w)(10)式中:cQ为冬季空调热负荷(w)；COP为制热性能系数。4.2模型的建立从江水取回水的结构模型如下图3：图3从江水取回水的结构模型模型尺寸：夏季，长800m，宽30m，水深7.5m；冬季，长800m，宽30m，水深4m。取水口位于排水口上游100m，排水管直径1m，排水口位于上游x=30m处，管上边沿与水表面齐。4.3边界条件的确定长江水为流动水体，释放到水体的冷热量会及时被水流带走，不会在当地的水体中聚集，且江水流动速度较大，上下层对流强烈，水体内在不断地进行热量交换，因此沿水深方向不形成温度分层，整个水体在某一横截面处呈等温分布。根据测得的数据，本工程项目长江水体初温设为：夏季23.80C；冬季10.70C。江水平均流速为夏季0.7m/s；冬季0.5m/s。江水表面根据理论计算，简化为常热流边界考虑，通过计算得到热流值：夏季n＝79（w/m2）；冬季n＝-15（w/m2），江水水体与土壤传热相对于水体表面而言较小，设为绝热边界。排水参数设置：夏季排水按极不利情况考虑，取100C的温差，即排水温度为33.80C，排水量为0.8128m3/s，排水管的出水流速1.04m/s。冬季按50C的温差考虑，即排水温度为5.70C，排水量为0.6694m3/s，排水管的出水流速0.853m/s。4.4模拟结果及分析4.4.1夏季模拟结果分析夏季模拟水源热泵机组连续运行1周排水对江水的影响，模拟结果显示：江水表层受影响最大，这是因为江面排水，排水方向与江水流速方向一致，且排水密度比江水密度要小的缘故。江水表层的模拟水温分布图如图4所示：6图4夏季连续运行一周后江水表面的水温分布20.022.024.026.028.030.032.0-2230100250400550700距排水口距离（m）水温（℃）20.022.024.026.028.030.032.00.00.51.01.52.02.5水深(m)水温(℃)上图5是水表面水温分布曲线，从图5可知：表层受影响范围为距排水口上游2m到距排水口下游740m长度区域水体,总共大约742m范围；表层影响区域内平均水温大约为24.10C。根据模拟结果，水温在排水口下游方向1m左右达到最高(为29.920C，温升为6.120C)，该截面处水温随水深变化关系曲线如图6所示，从图6可得出：沿水深方向，受影响的程度是逐层递减的，在2.5m深处，江水水温已不受排水的影响。所以排水在深度方向上影响范围为水面到水面以下2.5m的范围。从模拟计算结果还可得到以下数据：排水在纵向（河宽度方向）上影响范围为从岸边到17m的范围；排水影响区域内平均水温大约为23.970C。可知:即使取、回水温差取100C，