湖水源热泵系统排热工况下散点排水对湖体温度场的影响.docVIP

湖水源热泵系统排热工况下散点排水对湖体温度场的影响 摘要：采用流体计算软件建立三维湖体模型，模拟 设计日下湖体自然水温，用实测水温进行修正作为初始条 件；根据别墅群全年逐时负荷的模拟结果，采用动态取排水 温和水量作为边界条件；以典型气象年数据为基础，通过用 户自定函数输入水面换热边界；分别对某别墅群水库水源热 泵系统的排放水管网采用与雨水排水管网相结合的散点排 水以及采用传统的集中排水方式进行了动态数值模拟，得到 供冷季逐时排热工况下湖体温度场的分布情况，分散排水优 势明显。分析模拟结果表明：集中排水方案水体周温升 0. 33 °C,分散排水方案水体周温升0. 16 °C；水体水道狭长 体积有限的区域在连续排水的情况下，有明显热累积现象。 连续运行条件下，散点排水取水口水温的模拟结果表明：7? 10月，位于常年水位4. 6 m水深下的取水口取水温度在27? 23.5 ￡之间，可实现热泵系统高效运行。 关键词：水源热泵；散点排水；数值模拟；动态负荷； 温度分布 中图分类号：TU831文献标志码：A 文章编号(2013)湖水水体作为一种天然冷热源，实现了热泵系统高能效 比运行，产生良好节能效益［1］。以重庆为例，地表水水源 热泵系统得到了大力推广和广泛应用，依托辖区内的水系及 湖库塘资源，已建和在建地表水水源热泵项目有25个［2］。 目前，水源水质［34］、系统性能研究［58］和系统排水对水 体的影响［911］是本领域研究重点。温排水对于自然水体的 影响研究，国外文献多集中于发电厂类研究［1213］,国内学 者多针对二维水温计算方法［9］、特定排水方式下的温度场 特性研究［10］及不同排水方式的评价研究［11］。针对滞留水 体（湖库塘水）水源热泵系统在动态负荷边界条件下的研究 较少。 陈金华，等：湖水源热泵系统排热工况下散点排水对湖 体温度场的影响 地表水系可分河流和湖泊（水库、池塘），前者为流动 水体，而后者为滞留水体。两者的热承载能力大不相同［14］。 前者热承载能力大，流动性好，可以依靠流动的水较快地带 走热量，因此排放水对于水体的局部和总体温升影响一般都 不大；而水库水流动性差，热容量有限，对于空调系统持续 不断的排放水，长期运行后，湖体的承载能力是否会超出其 自恢复能力的极限就要另当别论了，这直接关系到系统的运 行效率以及对周围环境的影响。随着系统的运行，温排水排 入源水水体后，水体温度会发生改变。建筑负荷是全年动态 变化的，导致取水温度、热泵机组能效比、系统实际排热量 相互影响也是动态变化的。而目前的文献研究中采用的边界 条件多是在设计负荷下的，较少考虑多因素的影响动态过 程。 笔者采用以实测水温修正后的设计日下湖体自然水温 为初始条件，根据某别墅群全年逐时负荷的模拟结果，采用 动态取排水作为数值模拟的边界条件；以典型气象年数据为 基础，通过用户自定函数输入水面换热边界；对别墅群水库 水源热泵系统的排放水管网采用与雨水排水管网相结合的 散点排水方式进行了动态数值模拟，得到了供冷季排热工况 下湖体温度场的分布情况。 1控制方程及数值计算模拟流程图 非稳态、常物性、不可压缩流体湍流的湍流基本控制方 程组包括： 上述方程组中：u、v、w分别为x、y、z方向的速度瞬 时值，是时均值和脉动值之和；P为水体的密度；cp为水 体的定压比热容；P为压力；T为热力学温度；动量方程中 忽略体积力，能量方程忽略耗散项；ST为源项；qxyzt为排 水口的质量源项；en为水体表面的热交换量；ATh为热泵 系统的取排水温差。 △ Th 由式 ATh=QlhmxyztcpEERh+lEERh 确定，式中：Qlh 为建筑的逐时冷负荷；mxyzt为各排水口的设计质量流量； EERh为计算时刻下机组的能效比时均值。当△Th253C时， mxyzt为各排水口的设计流量；当ATh 考虑到水体模 型尺寸较大，而管径较小，实际划分网格时分块划分，并进 行网格过渡，最远处网格尺寸参照文献［16］的计算值，以保 证计算质量。 4定解条件 如图1所示，排水工况下的温度场模拟以自然水温为初 始条件，将建筑的逐时负荷作为能量源项加入数学模型，由 于模型假设水面高度不变，取排水量保持平衡，因此无需添 加质量源项。 水面边界条件：以典型气象年数据为准，按照文献［14］ 的方法动态输入水面边界。 水底土壤水体边界条件：以该地区的空调季节的土壤平 均温度为准，设为定温边界，参照文献［17］中的3.2 m以下 土壤层的温度取值23.25 °Co 水体水体交界面处的边界条件：即利用湖体与未利用湖 体的交界面。在此界面上，存在着复杂的热量和水体的交换, 结合水体模型建立过程中是以整个水体的中间界线作为 水体水体”的界面。再考虑到后期发展，假定左