近年来,空调水系统变流量控制技术在工程中得到大量应用,其节能效果亦得到了广大专业人士的认同。末端压差控制和冷热源总管温差控制是目前讨论较多的两种变流量控制方式。相比于末端压差控制,温差控制理论上节能效果较好,但是对管网形式的要求较高[1-3],一般建议使用在末端不设调节阀且各用户负荷变化规律相近的系统中,否则容易出现水/热力失调;然而,对于该问题的研究还仅仅是停留在定性分析的层次,温差控制在末端设置通断控制阀的系统中的适用性仍需进一步的讨论。
另一方面,近几年由于变频技术的不断发展,学者们对变流量技术的研究主要还是侧重于提升变流量控制下的机组能效[4-5],或者进行水系统整体控制策略的优化,最大程度降低整个机房的总能耗[6-9];也有少数学者结合风机盘管换热静特性提出了变温差控制法[3, 10],认为该控制方式下末端设备可以获得与实际负荷对应的流量需求,并通过与制冷机房其余设备的联合运行分析,指出变温差控制具有更高的节能潜力[11-12]。然而,上述研究在追求系统低能耗或高效率的同时,却忽略了末端设备在变工况下,尤其是在变温差控制下与房间的动态换热情况,房间最终的舒适性条件能否达到要求,有待进一步的研究。
针对上述问题,首先建立了末端设置通断控制阀的一次泵变流量风机盘管系统的计算模型,分析了其在温差控制下的水力特性,并结合风机盘管与房间的非稳态耦合换热特性,探讨了温差控制在末端设置通断阀的系统中运用的可能性;同时,对现有的变温差控制方法下房间温湿度的稳定性进行了分析,并在此基础上提出了改进建议。
1)研究对象的末端空调器均为采用一般通断控制阀门的风机盘管。
2)调整各支路静态平衡阀,使各末端通断控制阀门全部打开时所有末端流量均等于设计流量。
3)末端用户的两通阀在全开时的阻抗为常数,冷源支路的阻抗为常数,变流量机组蒸发器侧阻抗为常数。
4)为简化分析,不考虑旁通阀动作。
采用一维热流体仿真模拟软件Flowmaster建立模型并进行计算分析。模拟分析的集中空调水系统由4个支路组成,每个支路内有4个用户,为简化分析,假设每个用户的设计流量相同且均为10 m3/h,设计供回水温度为7/12 ℃。为方便水力平衡计算,设末端用户阻力为30 kPa,用户间管段阻力为5 kPa,支路间管段阻力为5 kPa,源侧阻力为100 kPa,PID控制器温差设定值取设计工况下的5 ℃,模拟运行时维持总管温差恒定。支路编号由远端至近端分别为1-1,1-2,…,4-3,4-4。模型如图 1所示。
温差控制通常建议使用在末端不设调节阀且各用户负荷变化规律相近的系统中,这样可以避免由于阀门动作和负荷变化引起的管网水/热力失调[1]。而对于末端设置通断阀的系统,由于部分用户负荷变化将导致末端调节阀的动作,管网负荷分布也具有多样性,势必导致其余用户不同程度的水力失调。因此,有必要对其温差控制下的水力特性进行量化分析,这将有助于直观地评价系统的整体特性,为优化控制提供参考和依据。
为便于分析,假设各用户负荷均为设计负荷,总管温差设定值取设计温差即5 ℃,分别研究负荷集中分布与均匀分布下各用户的流量分配情况。
1)负荷集中分布时各用户流量分配情况。表 1列出了6种不同负荷分布下,维持总管5 ℃温差时各用户的水力失调情况。从中可以看到,当负荷集中分布时,负荷较小的支路水力失调尤为严重,Case2中4-3用户过流度已经达到了34%;对比Case5和Case6以及Case1和Case6,负荷集中在下游或者管网整体负荷率较低时,管网整体的水力失调相对更加严重。
2)负荷均匀分布时各用户流量分配情况。表 2为9种不同负荷分布下各用户的水力失调情况,从中可见,无论管网整体负荷率如何,负荷均匀分布时,各支路水力失调度较小,在±10%以内,且各支路内远端用户的流量普遍大于近端用户的流量,支路内两端的用户水力失调要相对严重,这在设计时应予以注意。
通过建立风机盘管与房间的换热数学模型,可以分析末端设备不同水力失调度和不同负荷率下室内的温湿度变化情况,进而在此基础上探讨温差控制的适用条件以及变温差控制下室内温湿度的控制效果。
以接触系数法[13-14]为基础,建立风机盘管与房间的非稳态换热模型。
分析的假定条件为:
1)空气流量G不变;
2)进水温度tw1不变;
3)风机盘管的进风参数与室内空气的状态参数相同;
4)室内热湿负荷恒定。
根据传热学原理及空气特性可列出方程
式中:t1、ts1分别为盘管进风干湿球温度,℃;t2、ts2分别为盘管出风干湿球温度,℃;tw2为盘管出水温度,℃;cp、c分别为空气和水的比热,J/(kg·℃);E′为通用热交换效率,即接触系数;F为传热面积,m2;W为盘管水流量,kg/s;ξ为盘管瞬时析湿系数;K为盘管传热系数,W/(m2·℃);vy为表冷器迎面风速,m/s;w为盘管内水流速,m/s;p、m、n、r、s分别为实验系数及指数;h1、h2分别为进、出风焓值,kJ/kg。
对于特定的风机盘管,其表冷器结构参数和传热系数公式一般是已知的,在风量G一定的情况下,E’和vy也都是确定的,当水流量W和进水温度tw1已知时,上述5个方程共含有t1、t2、ts1、ts2、tw2、ξ、K等7个未知参数(h值可由干湿球温度确定[13])。一般对于表冷器的校核性计算往往是给定进风参数,以检验出风参数及处理冷量是否满足要求。但是当室内负荷不变,而水流量变化时,风盘与房间的换热过程是一个非稳态的相互耦合的换热过程。初始时刻,风盘进风参数是不变的,若风盘处理冷量由于水量变化而发生改变,室内空气的状态参数也随之改变,即风盘进风参数将逐渐发生变化,因此, 下一时刻风盘的处理冷量也将发生变化,反之, 又将影响风盘的进风参数,如此循环。
因此,要分析整个动态换热过程,还需引入风盘与房间的换热模型,参考文献[15]中空调房间的模型,并在此基础上做了改进,引入湿交换的换热关联式,具体为
式中:ρ为空气密度,取1.2 kg/m3;V为房间容积,m3;qx0为房间显热负荷,W;qxτ为风盘处理的显热冷量,W;qq0为房间潜热负荷,W;qqτ为风盘处理的潜热冷量,W;τ为时间,s;d为含湿量,g/kg;r为水蒸气的汽化潜热,为计算方便,取2 500 kJ/kg。
同时给定初始边界条件t1=t1m、d1=d1m、τ=0(t1m、d1m分别为设计工况下的进风干球温度和含湿量),联立式(1)~(8),进行迭代计算,即可得到各参数随时间变化的函数。
以某品牌FST324风机盘管为例,盘管传热面积F=18.25 m2,盘管通断面积fw=5.06×10-4 m2,额定工况:t1=27 ℃,ts1=19.5 ℃,d1=10.99 g/kg,tw1=7 ℃,tw2=12 ℃,G=1.08 kg/s,vy=2.5 m/s,E’=0.84,W=0.814 kg/s,全热冷量q=17 044 W,显热冷量qx=12 786 W,房间大小为5 m×3 m×3 m。盘管传热系数实验公式为
现假设水流量变为80%即0.642 kg/s,取时间步长Δτ=1 s,边界条件t1=27 ℃,d1=10.99 g/kg,τ=0,风量和进水温度不变,其余参数未知,进行迭代计算。图 2~4为计算后进风干球温度t1、进风湿球温度ts1、进风含湿量d1、风盘处理的全热冷量q以及显热冷量qx随时间变化的曲线。可以看到,由于水流量的突然减小,初始阶段风机盘管的处理冷量急剧下降,但是随着过程的进行,q和qx逐步回升,直至最后与室内热湿负荷相同,此时,室内干湿球温度、含湿量均有上升。从传热学角度不难理解,盘管水流量减少,房间的温湿度是不会无限上升的,理论上房间与风盘的换热最终将达到一个动态平衡的过程。其本质是因为水流量减小后,盘管处理能力下降导致室内干球温度及含湿量增加的同时,室内空气与盘管进行热湿交换的势差也在逐步增大,盘管处理能力也将增加。
由此可见,不管整个换热过程如何进行,在水流量变化后,风机盘管处理冷量的大小是逐步向室内热湿负荷逼近的,只是室内状态参数发生了变化。因此,如果是为了求解室内最终状态参数,可以以室内总负荷和显热负荷作为风机盘管最终所能处理的全热与显热冷量,即q和qx为已知参数,根据传热学原理可列出方程[13]
联立方程(1)~(5)、(10)~(11),所有未知参数的最终收敛值均能求出。
选取某品牌FST324和FST224两种不同型号的风机盘管,计算其不同设计工况在不同水力失调度下室内空气的最终状态参数,计算结果见表 3~4。
由表中可以看到,当水力失调度在±25%范围内,室内的相对湿度变化较小,均在0.5%以内,但是室内干球温度变化较大,最高达2 ℃以上;当水力失调度在±15%内变化时,室内干球温度变化较小,平均在1 ℃以内。而对舒适性空调而言,一定范围内的温湿度波动是可以接受的,对照文献[16]中的相关标准,表 3~4中的几种设计工况下,水力失调度在±15%范围内,室内最终温湿度均仍能满足人体舒适度的要求。
再观察表 2中所列负荷均匀分布时各用户的水力失调情况,均在±10%以内,这种程度的水力失调导致的室内温湿度变化将是很小的。
由此可见,对于末端设置通断阀的风机盘管系统,当管网负荷分布较为均匀时,温差控制法其实是较为适用的,但前提是各开启用户的负荷变化规律较为相近,例如数据机房、对温湿度控制精度要求不高的车间、厂房等。
对于采取末端压差和温差控制相结合的系统,当管网负荷分布集中、变化剧烈时,例如办公楼上下班高峰、午休及加班等时间段内,可以采取以末端压差控制为主的方式;而当管网负荷相对稳定、分布较为均匀时,则可适当以采取温差控制为主的方式,以弥补末端压差控制下管网整体欠流的不足,同时也能取得良好的节能效果和控制效果。另外,对于各支路内两端的用户,由于其在温差控制下的水力失调度相对更加严重,在设计时可以根据房间功能要求调整管道布置,或者在末端加设动态流量调节阀,以满足个别房间较高的舒适度要求。
根据表冷器的换热静特性,当处理冷量发生变化时,其进出口温差值是变化的,因此, 一些学者在此基础上提出了变温差控制的思想,建议以不同的总负荷率分段设置温差值。以FST324风盘为例,根据该思路进行了温差值设定计算,取设计工况为t1=27 ℃,ts1=19.5 ℃(φ1=50%),tw1=7 ℃,tw2=12 ℃,G=1.08 kg/s,E’=0.84,W=0.814 kg/s,全热冷量q=17 044 W,显热冷量qx=12 786 W。计算步骤为:
1)取温度控制精度±1 ℃,计算在额定温差和除湿能力(ξ)下,±1 ℃的温度波动所对应的负荷变化,即已知t1、tw1、tw2、ξ、G,根据式(1)~(5)、式(10)求解ts1、t2、ts2、W、K、q等6个参数。计算结果可见表 5,±1 ℃的温度波动所对应的负荷率为0.93~1.07,即±7%的负荷变化范围,该范围内湿度变化仅为±2.4%,满足文献[16]规定。
2)根据上一步计算结果对负荷进行分段划分,即每个负荷率节点对应于自身±7%的负荷范围为一个负荷区间。接着对各负荷率节点对应的工况进行计算,即已知t1、ts1、tw1、q、G(维持室内空气参数为设计工况),根据式(1)~(5)、式(10)求解t2、ts2、tw2、W、K、ξ等6个参数,即可得出各节点对应的温差设定值(tw2-tw1)。
3)在维持负荷率节点对应温差和除湿能力不变的情况下,计算出对应区间内负荷变化对室内温湿度的影响情况,即已知tw1、tw2、ξ、q、G,根据式(1)~(5)、式(10)求解t1、ts1、t2、ts2、W、K等6个参数。
计算结果见表 5,由表 5可以看到,该负荷分段划分方法是可行的,理论上基本满足室内温湿度的控制要求。但同时也发现,低负荷时表冷器的除湿能力将大幅下降,分析步骤2)计算温差值的过程,ξ值是未知数,意味着表冷器的除湿能力是得不到保证的,最终室内温湿度很可能达不到要求。
因此,为了维持室内温湿度的稳定,仅以总负荷为标准分段设置温差值是不恰当的,因为总负荷并不能反映显热与潜热负荷的权重,尽管风盘处理的全热冷量足够,但如果其中的显热和潜热分量与实际热湿负荷相差较大,室内温湿度仍将出现较大的波动。事实上,部分负荷时,仅仅改变温差是不够的,必须同时改变风量,以增大或减少接触系数,才能在维持室内温湿度不变的同时,保证表冷器的除湿能力。然而,增加风量控制将极大地增加系统的控制成本和难度,工程上也不易实现。
笔者建议,以总负荷为标准分段设置温差值时,可以根据建筑物自身负荷特点,保证表冷器一定的除湿能力,而不必追求室内湿度的稳定,即在步骤2)中,将ξ作为已知参数,而ts1作为未知参数,进行温差设定值的计算。表 6为采用该方法的计算结果,为了分析方便,各负荷率节点对应的ξ均设定为1.333,即保持表冷器除湿能力为设计工况不变。从表中可以看到,维持表冷器除湿能力不变的条件下,各符合率节点对应的室内相对湿度有所上升,但仍符合文献[16]的规定,满足舒适性空调的要求,同时还可保证室内温度的相对稳定。
按照第1种温差值设定方法,若总负荷不变而实际潜热负荷为表 6中对应的潜热负荷,室内最终温湿度变化将如表 7所示。可以看到,当实际潜热负荷所占总负荷比重较大时,室内最终温湿度会有较大程度的变化。对比之下,第2种方法虽然同样导致室内相对湿度发生变化,但优点在于室内干球温度的稳定性可以得到保证。
1)基于温差控制的末端设置通断阀的一次泵变流量系统,当管网负荷集中分布时,将造成低负荷支路较大程度的水利失调;当管网负荷均匀分布时,无论整体负荷率如何,各支路水利失调度均较小,但各支路两端的用户水力失调要相对严重,在设计时应予以注意。
2)在供回水温度不变的情况下,当末端设备水力失调度较小时,室内温湿度的波动也较小。因此,当管网负荷分布均匀时,温差控制法是较为适用的,前提是各开启用户的负荷变化规律比较相近。
3)维持室内温湿度不变的条件下,仅以总负荷为标准分段设置温差值是不恰当的,实际运行时并不能保证室内温湿度的稳定。建议以优先保证风盘除湿能力和室内干球温度为原则设定温差值,而不必追求室内湿度的稳定,相比于现有方法,虽然同样导致室内相对湿度产生较大波动,但室内温度的稳定性可以得到保证。