湖水源熱泵溫(冷)排水三維溫度場數(shù)值模擬

?！?韓龍喜

(河海大學(xué)環(huán)境學(xué)院,江蘇 南京　210098)

?

湖水源熱泵溫(冷)排水三維溫度場數(shù)值模擬

祝悅,韓龍喜

(河海大學(xué)環(huán)境學(xué)院,江蘇 南京210098)

摘要：結(jié)合南京某湖水源熱泵工程項目,利用三維數(shù)學(xué)模型對系統(tǒng)運行時溫(冷)排水溫度場進行數(shù)值模擬,預(yù)測溫(冷)排水?dāng)U散范圍及溫升(降)程度。結(jié)果表明:該項目造成的湖泊溫升(降)滿足地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)要求;溫、冷排水受密度影響,從湖泊表層排放后擴散情況不同,溫排水對湖泊表層水溫影響較大,排放口附近區(qū)域水溫層趨水平方向分布,而冷排水對湖泊底層水溫影響較大,排放口附近區(qū)域水溫層趨垂直方向分布;相同排放流量、排水方式情況下,冷排水對湖泊水溫的影響更大。

關(guān)鍵詞：湖水源熱泵;溫(冷)排水;三維溫度場;數(shù)值模擬

地表水源熱泵是一種以江、河、湖、海等地球表面的水體作為熱源進行制冷/制熱循環(huán)的熱泵,在制熱的時候以水作為熱源,在制冷的時候以水作為排熱源[1]。地表水源熱泵技術(shù)是一種可再生能源利用技術(shù),具有運行穩(wěn)定、環(huán)境效益顯著、高效節(jié)能等優(yōu)勢,近年來被大力推廣[2-3]。湖水源熱泵作為地表水源熱泵的一種,其應(yīng)用方案的可行性已經(jīng)在實際應(yīng)用中得到證明[4-6],但其研究主要集中在技術(shù)應(yīng)用上,而針對湖水源熱泵溫(冷)排水在受納水域的溫度擴散規(guī)律的研究很少。

關(guān)于溫排水對自然水體的影響,國內(nèi)外研究多集中于發(fā)電廠冷卻水對海洋熱環(huán)境影響方面[7-9],而有關(guān)冷排水造成受納水域溫降對水環(huán)境的影響的研究較少。對溫排水的數(shù)值模擬,目前多集中在沿垂向積分的平面淺水二維模型上,但二維模型由于深度平均的緣故,無法估計排水口附近水溫層厚度,模擬溫(冷)排水時,可能會導(dǎo)致平面上溫升(降)區(qū)面積偏小[10]。

筆者結(jié)合南京某湖水源熱泵工程項目,利用三維數(shù)學(xué)模型對溫(冷)排水熱擴散情況進行數(shù)值模擬,預(yù)測溫(冷)排水?dāng)U散范圍及溫升(降)程度,分析排水口附近區(qū)域水溫層分布特征,探討溫排水和冷排水對排水口附近區(qū)域水溫的不同影響,旨在為湖水源熱泵系統(tǒng)的設(shè)計和運行提供參考。

1工程概況

1.1　建筑和氣候概況

南京某度假中心項目總建筑面積90 999 m2,是集旅游、度假、休閑、運動、商務(wù)、生活為一體的綜合建筑群。其酒店與部分商務(wù)樓空調(diào)系統(tǒng)設(shè)計采用以地表水為冷(熱)源的水源熱泵系統(tǒng)。該系統(tǒng)就近利用項目附近湖泊水的低溫低位熱能資源,利用熱泵技術(shù),通過少量的高位電能輸入,實現(xiàn)低位熱能向高位熱能的轉(zhuǎn)移,進行夏季制冷,冬季供暖,并兼供熱水。

項目所在區(qū)域常年平均氣溫為15.1℃,較高的1961年平均氣溫為16℃,較低的1969年平均氣溫為14.3℃。2012年最低月平均氣溫為2.2℃(1月),最高月平均氣溫為29.2℃(7月)。全年無霜降期為228 d,平均氣壓為101.55 kPa,略高于海平面標(biāo)準(zhǔn)氣壓。境內(nèi)最多風(fēng)向為東風(fēng),夏季多東風(fēng)和東南風(fēng),多年平均風(fēng)速2.8 m/s;冬季多東北風(fēng)和西北風(fēng),多年平均風(fēng)速為2.3 m/s。多年平均降水量為999.8 mm,多年平均蒸發(fā)量為1 459.6 mm。

1.2　冷熱源方案

對項目取水口附近的3個采樣點進行水質(zhì)監(jiān)測,湖水水質(zhì)為Ⅱ～Ⅲ級,滿足水源熱泵系統(tǒng)對水質(zhì)的要求。湖水由取水口泵送經(jīng)簡單處理后,進入熱泵主機的換熱器,在封閉環(huán)路中循環(huán)后排入湖內(nèi)。取水量和退水量相等,沒有水的損耗,水質(zhì)基本無變化,夏季排水水溫溫升3℃,冬季排水水溫溫降2℃。退水實施分散淺排,可減小對湖泊水體的影響[11]。

湖水源熱泵系統(tǒng)負(fù)擔(dān)夏季冷負(fù)荷4 393.07 kW,冬季熱負(fù)荷2 907.24 kW。夏季空調(diào)制冷期間,設(shè)計取水量為750 m3/h,每天開機12 h,制冷期為120 d;冬季供暖期間,取水量為750 m3/h,每天開機12 h,供暖期為90 d;春秋季節(jié)基本停用,僅啟新風(fēng)和輔助加熱水。

2三維水溫數(shù)值模擬

為實現(xiàn)溫(冷)排水的精細(xì)模擬,筆者建立三維水溫數(shù)學(xué)模型,分別對溫排水和冷排水熱擴散情況進行數(shù)值模擬。

2.1　模型介紹

筆者采用MIKE3建立數(shù)學(xué)模型,對溫(冷)排水熱擴散情況進行數(shù)值模擬。MIKE3軟件由丹麥水資源及水環(huán)境研究所(DHI)研發(fā),可用于模擬港口、河流、湖泊水庫、河口和海洋的水動力、水質(zhì)、泥沙輸移問題,具有先進的前后處理功能和友好的用戶界面。MIKE3模型的數(shù)學(xué)基礎(chǔ)是雷諾平均化的N-S方程,它包括了紊流影響及密度變化,同時包含了鹽度及溫度平衡方程,其水動力學(xué)模塊采用交替方向隱式迭代法(ADI方法)對質(zhì)量及動量守恒方程進行積分,對其產(chǎn)生的數(shù)學(xué)矩陣采用雙精度掃描法(DOUBLE SWEEP)進行求解[12]。

2.2　模擬范圍和網(wǎng)格構(gòu)建

為探討湖水源熱泵溫(冷)排水對受納湖泊水溫結(jié)構(gòu)的影響,以整個湖區(qū)為模擬范圍,但重點關(guān)注排放口附近水溫變化情況。本研究受納湖泊面積約8.8 km2,平均水深約8 m,采用該湖泊2013年實測水下地形作為模型計算條件,離散得到各個節(jié)點的湖底高程,見圖1。為了較精確擬合天然水域的復(fù)雜邊界,水平方向上采用無結(jié)構(gòu)三角網(wǎng)格劃分計算區(qū)域,在排水口附近區(qū)域適當(dāng)加密以提高模型的計算精度,非加密區(qū)域網(wǎng)格邊長約100 m,加密區(qū)域網(wǎng)格邊長約為50m,垂直方向上采用均勻分層模式,即任何兩層網(wǎng)格的間距與總水深的比值相等,垂直方向上分為5層,研究區(qū)域平面網(wǎng)格布置如圖2所示。

圖1　研究區(qū)域地形(單位：m)

圖2　平面網(wǎng)格布置

2.3　定解條件及參數(shù)取值

2.3.1定解條件

a. 初始水位。根據(jù)該項目湖泊水文站1980—2001年水文監(jiān)測資料,該湖泊多年平均水位22.8 m。據(jù)2012年1月至2013年9月逐月水位監(jiān)測數(shù)據(jù),2012年最高水位23.51 m(9月),最低水位22.18 m(7月),平均水位23.03 m,2013年最高水位23.57 m(6月),最低水位22.43 m(9月),平均水位23.17 m。采用多年平均水位作為模擬的初始水位條件。

b. 初始水溫。根據(jù)下關(guān)長江水文站資料,南京地區(qū)夏季7、8月份平均水溫一般在30℃,冬季1、2月份平均水溫一般在7℃,可作為該項目附近湖泊月平均水溫。實測該湖泊水溫,2014年1月4日水溫為6.3℃,2014年1月10日水溫為5.6℃,與南京地區(qū)相同季節(jié)時的地表水溫基本一致。由于湖泊表層水溫和底層水溫相差不大,為著重探討溫(冷)排水對湖泊的水溫結(jié)構(gòu)影響,不考慮湖泊自然水溫分層,采用南京地區(qū)夏季和冬季平均水溫分別作為模擬的初始水溫條件。

c. 流量邊界條件。排水處所在計算單元的流量邊界條件根據(jù)實際退水情況給定;出庫流量邊界采用自由出流邊界條件。

d. 水溫邊界條件。排水處所在計算單元的水溫邊界條件根據(jù)實際退水水溫變化情況給定;出庫水溫邊界條件設(shè)定為湖泊自然水溫。

2.3.2參數(shù)取值

a. 糙率。依據(jù)資料,預(yù)測區(qū)域糙率取值范圍為0.020～0.033,本項目糙率取為0.028。

b. 熱擴散系數(shù)。在MIKE3軟件中,熱擴散系數(shù)是按照比例系數(shù)乘以紊動黏性系數(shù)計算得出,其中,紊動黏性系數(shù)由軟件在進行模擬計算時自動算出,比例系數(shù)由人為給定,取軟件用戶手冊推薦值1.0。

c. 風(fēng)場摩擦力。本次模擬考慮靜風(fēng)條件下溫(冷)排水熱擴散情況。

2.4　模型驗證

水體的水動力特征是分析研究溫(冷)排水熱擴散特征的基礎(chǔ),由于該湖水源熱泵系統(tǒng)還未建成使用,尚無熱泵系統(tǒng)運行時流場實測資料,故對湖泊原流場進行數(shù)值模擬,并利用排放口附近代表性測點實測資料進行水動力模型的率定和驗證。排放口附近,平行于岸邊各測點(離岸80 m)表層流場流速驗證見圖3,垂直岸邊各測點表層流場流速驗證結(jié)果見圖4。從模型流場的驗證結(jié)果來看,代表性測點流速平均誤差小于15%。

圖3　平行于岸邊各測點流速驗證結(jié)果

圖4　垂直岸邊各測點流速驗證結(jié)果

2.5　模擬結(jié)果與分析

2.5.1夏季溫度場

夏季熱泵系統(tǒng)排水量750 m3/h,排水溫升3℃,湖水自然水溫30℃,連續(xù)運行一周后湖體溫升分布情況如圖5所示。模擬結(jié)果是在忽略土壤與湖體換熱基礎(chǔ)上得出的,并且考慮機組連續(xù)運行情況,所以實際湖體熱擴散范圍更小,溫升(降)更小。

由圖5可知,夏季溫排水對湖泊表層水溫影響范圍和影響程度最大,影響隨著水深的增加顯著減小。表層最大溫升區(qū)域溫升約為0.46℃,中層最大溫升區(qū)域溫升約為0.21℃,底層最大溫升區(qū)域溫升約為0.15℃。在排放口附近區(qū)域溫度梯度很大,距離越遠(yuǎn)溫度梯度越小。

退水對湖體造成的溫升均小于1℃,滿足地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)要求:人為造成的環(huán)境水溫變化應(yīng)限制在:周平均最大溫升≤1℃[13]。

2.5.2冬季溫度場

冬季熱泵系統(tǒng)排水量750 m3/h,排水溫降2℃,湖水自然水溫7℃,運行一周后湖體溫降分布情況如圖6所示。

圖5　夏季各層溫升分布

圖6　冬季各層溫降分布

由圖6可知,除極靠近排水口處,冬季冷排水對湖泊底層水溫影響范圍和影響程度最大,影響范圍隨著水深的增加而增加。但由于冷排水從湖泊表層排放,最大溫降區(qū)域在湖泊表層靠近排水口處,最大溫降約為0.63℃,中層最大溫降區(qū)域溫降約為0.33℃,底層最大溫降區(qū)域溫降約為0.30℃,但冷排水對底層的影響范圍最大。與溫排水相似,冷排水在排放口附近區(qū)域溫度梯度很大,距離越遠(yuǎn),溫度梯度越小。

退水對湖體造成的溫降均小于2℃,滿足地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)要求,即人為造成的環(huán)境水溫變化應(yīng)限制在周平均最大溫降≤2℃[13]。

2.5.3夏季溫度場和冬季溫度場比較

圖5和圖6對比可以看出,由于夏季溫排水密度比自然水體密度小,從湖泊表層排放后,溫排水主要分布在湖泊表層,隨水流向四周擴散,熱擴散呈輻射狀,對湖泊表層水溫影響較大,排放口附近區(qū)域水溫層趨水平方向分布;而冬季冷排水密度比自然水體密度大,從湖泊表層排放后,沉降的同時隨水流向壩址處流動,擴散距離較遠(yuǎn),且對湖泊底層水溫影響范圍較大,排放口附近區(qū)域水溫層趨垂直方向分布。

排放流量相同、排水方式相同情況下,溫排水(溫升3℃)造成湖泊最大溫升約為0.46℃,冷排水(溫降2℃)造成湖泊最大溫降約為0.63℃,且冷排水對湖泊底層造成較大影響。可見冷排水對湖泊水溫的影響更為嚴(yán)重，原因在于溫排水分布在湖泊表層,有利于與外界進行熱交換,冷排水的沉降特征導(dǎo)致與外界的熱交換減少。

3結(jié)論

a. 依照該項目水源熱泵空調(diào)系統(tǒng)的規(guī)模,夏季溫排水造成湖泊最大溫升約為0.46℃,冬季冷排水造成湖泊最大溫降約為0.63℃,滿足地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)要求。

b. 溫(冷)排水受密度影響,從湖泊表層排放后擴散情況不同,溫排水對湖泊表層水溫影響較大,排放口附近區(qū)域水溫層趨水平方向分布;冷排水對湖泊底層水溫影響較大,排放口附近區(qū)域水溫層趨垂直方向分布。

c. 排放流量相同、排水方式相同的情況下,冷排水對湖泊水溫的影響更大,在進行湖水源熱泵系統(tǒng)的設(shè)計和運行時,應(yīng)重點考慮冷排水對湖泊的影響。

d. 設(shè)計時,排水口應(yīng)盡量靠近水體表面,分散淺排,有利于溫(冷)排水與外界進行熱交換,減小對湖泊的影響。

參考文獻：

[ 1 ] 張軍.地表水源熱泵的發(fā)展現(xiàn)狀以及面臨問題[J].制冷空調(diào)與電力機械,2007,28(6): 73-77.(ZHANG Jun.Application status and problems of surface water heat pumps for the present[J].Refrigeration Air Conditioning & Electric Power Machinery,2007,28(6): 73-77.(in Chinese))

[ 2 ] 孫強.淺談水源熱泵技術(shù)的國內(nèi)外發(fā)展現(xiàn)狀及趨勢[J].佳木斯大學(xué)學(xué)報: 自然科學(xué)版,2007,25(3): 433-434.(SUN Qiang.On the development of water loop source heat pump technology[J].Journal of Jiamusi University: Natural Science Edition,2007,25(3): 433-434.(in Chinese))

[ 3 ] 狄彥強,王清勤,袁東立,等.水源熱泵的應(yīng)用與發(fā)展[J].制冷與空調(diào),2006,6(5): 1-4.(DI Yanqiang,WANG Qingqin,YUAN Dongli,et al.Discussion about the application and development of water source heat pumps[J].Refrigeration and Air-conditioning,2006,6(5): 1-4.(in Chinese))

[ 4 ] 劉婷婷,彭建國,張國強,等.湖水水源熱泵空調(diào)系統(tǒng)在湖南省的應(yīng)用實例分析[J].建筑熱能通風(fēng)空調(diào),2004,23(6): 40-44.(LIU Tingting,PENG Jianguo,ZHANG Guoqiang,et al.Application of lake water heat pump system in Hunan Province[J].Building Energy & Environment,2004,23(6): 40-44.(in Chinese))

[ 5 ] 陳金華,劉勇,丁勇,等.重慶開縣人民醫(yī)院湖水源熱泵空調(diào)系統(tǒng)實測分析[J].制冷與空調(diào),2008,38(8): 111-114.(CHEN Jinhua,LIU Yong,DING Yong,et al.Test and analysis for lake-water source heat pump air conditioning system in People’s Hospital of Kaixian,Chongqin[J].Refrigeration and Air-conditioning,2008,38(8): 111-114.(in Chinese))

[ 6 ] 劉有勢,馬滿英,李志仁,等.株洲市神農(nóng)城湖水源熱泵空調(diào)系統(tǒng)的應(yīng)用技術(shù)提煉[J].應(yīng)用能源技術(shù),2014(9): 6-10.(LIU Youshi,MA Manying,LI Zhiren,et al.Application technology extraction of lake-water-source heat pump air-conditioning system in Zhuzhou Shennong town[J].Applied Energy Technology,2014(9): 6-10.(in Chinese))

[ 7 ] 石繼香,官云蘭,趙利民,等.國內(nèi)外溫排水熱影響監(jiān)測與評價研究方法探討[J].江西科學(xué),2014,32(1): 73-77.(SHI Jixiang,GUAN Yunlan,ZHAO Limin,et al.The investigation of monitoring and evaluation methods on thermal impact of thermal discharge from NPPs in domestic and abroad[J].Jiangxi Science,2014,32(1): 73-77.(in Chinese))

[ 8 ] ARIELI R N,ALMOGI-LABIN A,ABRAMOVICH S,et al.The effect of thermal pollution on benthic foraminiferal assemblages in the Mediterranean shoreface adjacent to Hadera power plant (Israel)[J].Marine Pollution Bulletin,2011,62: 1002-1012.

[ 9 ] TATIANA P T,LEONARDO M N,FRANCISCO G A.Effects of a nuclear power plant thermal discharge on habitat complexity and fish community structure in Ilha Grande Bay,Brazil[J].Marine Environmental Research,2009,68: 188-195.

[10] 端木琳,李震,蔣爽,等.大連星海灣海水源熱泵空調(diào)系統(tǒng)熱擴散數(shù)值模擬研究 [J].太陽能學(xué)報,2008,29(7): 832-836.(DUANMU Lin,LI Zhen,JIANG Shuang,et al.Numerical simulation of the thermal discharge in Xinghai bay seawater source heat pump project[J].Acta Energiae Solaris Sinica,2008,29(7): 832-836.(in Chinese))

[11] 卿曉霞,崔相楠,周健,等.江水源熱泵尾水排放方式優(yōu)選試驗研究[J].土木建筑與環(huán)境工程,2011,33(7): 105-110.(QING Xiaoxia,CUI Xiangnan,ZHOU Jian,et al.Experimental analysis on optimal tail water discharge method of river water source heat pump[J].Journal of Civil,Architectural & Environmental Engineering,2011,33(7): 105-110.(in Chinese))

[12]馬騰,劉文洪,宋策,等.基于MIKE3的水庫水溫結(jié)構(gòu)模擬研究[J].電網(wǎng)與清潔能源,2009,25(2): 68-71.(MA Teng,LIU Wenhong,SONG Ce,et al.Application of MIKE3 in prediction analysis of the water temperature of reservoir[J].Power System and Clean Energy,2009,25(2): 68-71.(in Chinese))

[13] GB 3838—2002地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)[S].

Numerical simulation of three-dimension thermal environmental impact of lake-water source heat pump warm/cold water discharge

ZHU Yue,HAN Longxi

(CollegeofEnvironment,HohaiUniversity,Nanjing210098,China)

Abstract：In combination with Lake-water source heat pump project in Nanjing, three-dimensional modeling was established to simulate the thermal environmental impact of warm/cold water discharge when the system is operating, and predict the range of warm/cold water discharge diffusion and the extent of temperature rise/drop. The results show that: the lake-water temperature rise/drop caused by the project meets the requirement of environmental quality standards for surface water. Affected by water density, thermal diffusion of discharged warm water is different from cold water and the discharged warm water has a greater impact on the temperature of surface lake-water and water temperature layering appear to be horizontal near the outfall area, while the discharged cold water has a greater impact on bottom water and the water temperature layering appear to be vertical near the outfall area. With the same emission load and the same emissions way, cold water discharge has greater impact on the water temperature of the lake.

Key words：lake-water source heat pump; warm/cold water discharge; three-dimensional temperature field; numerical simulation

(收稿日期：2015-03-26編輯：彭桃英)

中圖分類號：X524

文獻標(biāo)志碼：A

文章編號：1004-6933(2016)01-0130-05

作者簡介:祝悅(1990—),女,碩士研究生,研究方向為水動力及擴散質(zhì)輸運數(shù)值模擬。E-mail:zhuyue5151@sina.com

DOI：10.3880/j.issn.1004-6933.2016.01.023