一种再生水源热泵供暖系统协同优化运行方法

摘要

本发明公开了一种再生水源热泵供暖系统协同优化运行方法，以不同负荷率下再生水源热泵供暖系统能耗最低为目标函数，建立再生水源热泵供暖系统协同优化模型，并确定协同因子；带入已知参数，根据协同优化模型进行计算，得出不同负荷率下的再生水源热泵供暖系统优化解及相应运行方法。该协同优化运行方法能很好满足实际工程的需要，有效减少热泵系统能耗，节省运行费用，便于管理人员操作，具有很好的应用前景。

权利要求

1.一种再生水源热泵供暖系统协同优化运行方法，其特征在于，包括下述步骤：

1)根据水源侧供水流量mG，以再生水提升泵能耗Ppump,re最低为目标，确定非设计工况下水源侧输配系统运行方法，得到再生水提升泵能耗：Ppump,re＝f1(mG)  (1)；

2)根据地辐管散热量Q和用户末端负荷需求Qdemand，得出不同用户末端负荷需求Qdemand下相应的用户侧供水温度tL,in；

3)根据热泵机组制热量Qh、热泵机组输入功率Powerh、修正系数ηQ-t、ηP-t及蒸发器换热量Qe，得出不同负荷率x下热泵机组的蒸发器换热量Qe与蒸发器进口温度tE,in的关系式：Qe＝f2(tE,in)  (2)；

4)根据水源侧换热量Qre、中介水换热量Qmc、蒸发器换热量Qe及蒸发器进口温度tE,in，得到水源侧供水流量：mG＝f3(tE,in)  (3)；

5)根据得到的公式(1)、(2)和(3)，以不同负荷率x下再生水源热泵供暖系统能耗P(x)最低为目标函数，分别以蒸发器进口温度tE,in、热泵机组制热量Qh、水源侧供水流量mG范围为约束条件，建立再生水源热泵供暖系统协同优化模型；

6)根据步骤5)的再生水源热泵供暖系统协同优化模型，确定协同因子为蒸发器进口温度tE,in；

7)根据再生水提升泵样本参数，拟合出再生水提升泵相应关系式，然后根据再生水提升泵相应关系式结合热泵机组参数，代入再生水源热泵供暖系统协同优化模型，确定协同因子数值，由协同因子数值得出水源侧输配系统再生水提升泵的开启台数和变频大小，及热泵机组的开启台数和运行负荷大小。

2.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤1)中，再生水提升泵能耗Ppump,re具体计算式如下：其中：HP为再生水提升泵扬程，m；γ为流体的容重，kN/m3；ηp为再生水提升泵效率，％；ηm为电动机效率，％；ηvfd为变频器效率，％。

3.根据权利要求2所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，电动机效率ηm通过下式计算得到：ηm＝94.187(1-e-0.0904ω)

变频器效率ηvfd通过下式计算得到：

ηvfd＝50.87+1.263ω-0.0142ω2+5.834×10-5ω3其中：ω为再生水提升泵转速比，％；

所述非设计工况下水源侧输配系统运行方法包括：阀门节流调节和水泵变速调节。

4.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤2)中，地辐管散热量Q具体表达式如下：Q＝Qu+Qd＝Qc+Qr＝Qdemand

其中：Qu为地辐管向上的传热量，W；Qd为地辐管向下的热损失，W；Qc为地面的对流散热量，W；Qr为地面的辐射散热量，W；Qdemand为用户末端负荷需求，W；

用户末端负荷需求Qdemand计算式具体为：

Qdemand＝cpmu(tL,in-tL,out)其中：cp为水的定压比热容；mu为用户侧循环水流量，m3/s；tL,in、tL,out分别为用户侧供水、回水温度，℃。

5.根据权利要求4所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，地辐管向上的传热量Qu和地辐管向下的热损失Qd分别通过下式计算得到：其中：F为地板面积，m2；tr为室内空气温度，℃；Ku为地面向上传热系数，W/(m2·K)；Kd为2

地面向下传热系数，W/(m·K)；

地面的对流散热量Qc和辐射散热量Qr分别通过下式计算得到：Qc＝2.17F(tf-tr)1.31

Qr＝5×10-8F[(tf+273)4-(UMRT+273)4]其中：tf为地板表面温度，℃；UMRT为非加热面的平均辐射温度，℃。

6.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤3)中，热泵机组制热量Qh、热泵机组输入功率Powerh及修正系数ηQ-t、ηP-t分别为：Qh＝Qh,ref·ηQ-t·PLR

Powerh＝Powerh,ref·ηP-t·ηPLR其中：Qh,ref为热泵机组额定制热量，W；Powerh,ref为热泵机组额定输入功率，W；ηQ-t、ηP-t分别为热泵机组制热量、热泵机组输入功率关于机组进出口温度变化的修正系数；PLR为热泵部分负荷容量；ηPLR为不同PLR下的热泵机组输入功率修正系数；C1—C6、D1—D6均为拟合系数；

蒸发器换热量Qe关系式为：

Qe＝Qh-Powerh。

7.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤4)中，水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe之间关系为：Qre＝Qmc＝Qe

水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe计算公式分别为：Qre＝mGcp(t1,in-t1,out)

Qmc＝kexFexΔt

Qe＝mZcp(tE,in-tE,out)

其中：Fex、kex分别为宽流道换热器换热面积和传热系数，m2、kW/m2·℃；t1,in、t1,out分别为水源侧供水、回水温度，℃；Δt为对数平均温差，Δt＝(Δt'-Δt”)/(lnΔt'/Δt”)，Δt'＝t1,in-tE,in，Δt”＝t1,out-tE,out，℃；tE,out为蒸发器出口温度，℃；mG、mZ分别为水源侧供水、中介水流量，m3/s。

8.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤5)中，再生水源热泵供暖系统协同优化模型具体为：minP(x)＝Pheatpump(x,tE,in)+Ppump,re(x,mG)+P′pump式中：Pheatpump(x,tE,in)、Ppump,re(x,mG)分别为负荷率x时对应的热泵机组能耗和再生水提升泵能耗，P′pump为定频水泵能耗，包含中介水循环泵Ppump,mc、用户侧循环泵Ppump,user，两者通过台数调节来满足流量需求，某一负荷率x下的P′pump可视为定值，k为热泵机组额定制热量修正系数。

9.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤7)中，再生水提升泵样本参数包括再生水提升泵在某一流量下对应的扬程及效率；热泵机组参数包括某一水源侧进出口温度与对应的制热量及能耗。

10.根据权利要求1所述的再生水源热泵供暖系统协同优化运行方法，其特征在于，所述步骤7)中，拟合出相应关系式包括：再生水提升泵扬程与水源侧供水流量关系式：

再生水提升泵效率与水源侧供水流量关系式：

其中，a1、a2、a3、b1、b2、b3均为拟合系数。

说明书

一种再生水源热泵供暖系统协同优化运行方法

技术领域

[0001] 本发明属于水源热泵系统工程技术领域，特别涉及一种再生水源热泵供暖系统协同优化运行方法。

背景技术

[0002] 再生水源热泵是一种高效节能系统。再生水中蕴含了大量的低品位热能，环境友好，水量稳定，在冬季提取再生水中的热量用于建筑物供暖，能够有效减少环境污染并适当缓解能源缺乏的现状。

[0003] 随着再生水源热泵系统应用数量的增多，由于缺乏对热泵系统协同优化的认识，导致运行管理不当和能源浪费的现象也逐渐引起重视。目前，在应用再生水源热泵系统时，运行人员往往只考虑优化输配系统能耗或热泵主机能耗，但再生水源热泵供暖系统是一个由多个子系统组成的复杂系统；同时，大多数运行时间均在非设计工况下工作，单纯考虑某一设备或子系统对其进行能耗优化并不能表示热泵系统能耗达到最低，因此，为了降低热泵系统能耗，必须对各子系统进行协同优化。

[0004] 在现有应用中，为了在满足不同用户需求的同时降低系统能耗，需要对热泵系统的多个参数进行控制调节，但操作难度较大且节能效果不明显。热泵系统中的水源侧输配系统和热泵机组耗能较大，且两者能耗存在此消彼长的关系，故水源侧输配系统和热泵机组之间存在一个协同因子，再生水源热泵供暖系统协同优化运行的目标就是找到协同因子，并以降低热泵系统能耗为目标，确定不同负荷率下协同因子的数值。

发明内容

[0005] 针对目前再生水源热泵供暖系统难以实现协同优化、运行能耗较高的问题，本发明的目的是提供一种通用的再生水源热泵供暖系统协同优化运行方法，在满足不同用户末端负荷需求的同时，使整个再生水源热泵供暖系统能耗达到最低，实现系统的整体最优化运行。

[0006] 为实现上述目的，本发明采用如下技术方案：

[0007] 一种再生水源热泵供暖系统协同优化运行方法，包括下述步骤：

[0008] 1)根据水源侧供水流量mG，确定非设计工况下水源侧输配系统运行方法，得到再生水提升泵能耗：

[0009] Ppump,re＝f1(mG)  (1)；

[0010] 2)根据地辐管散热量Q表达式和用户末端负荷需求Qdemand，得出不同用户末端负荷需求Qdemand下的用户侧供水温度tL,in；

[0011] 3)根据热泵机组制热量Qh、热泵机组输入功率Powerh、修正系数ηQ-t、ηP-t及蒸发器换热量Qe，得出不同负荷率x下热泵机组的蒸发器换热量Qe与蒸发器进口温度tE,in的关系式：

[0012] Qe＝f2(tE,in)  (2)；

[0013] 4)根据水源侧换热量Qre、中介水换热量Qmc、蒸发器换热量Qe及蒸发器进口温度，得到水源侧供水流量：

[0014] mG＝f3(tE,in)  (3)；

[0015] 5)根据得到的公式(1)、(2)和(3)，以不同负荷率x下再生水源热泵供暖系统能耗P最低为目标函数，建立再生水源热泵供暖系统协同优化模型；

[0016] 6)根据步骤5)的再生水源热泵供暖系统协同优化模型，确定协同因子为蒸发器进口温度tE,in；

[0017] 7)根据再生水提升泵样本参数，拟合出再生水提升泵相应关系式，然后根据再生水提升泵相应关系式结合热泵机组参数，代入再生水源热泵供暖系统协同优化模型，确定协同因子数值，由协同因子数值得出水源侧输配系统再生水提升泵的开启台数和变频大小，及热泵机组的开启台数和运行负荷大小。

[0018] 进一步，所述步骤1)中，再生水提升泵能耗Ppump,re具体计算式如下：

[0019]

[0020] 其中：HP为再生水提升泵扬程，m；γ为流体的容重，kN/m3；ηp为再生水提升泵效率，％；ηm为电动机效率，％；ηvfd为变频器效率，％。

[0021] 进一步，所述的电动机效率ηm通过下式计算得到：

[0022] ηm＝94.187(1-e-0.0904ω)

[0023] 变频器效率ηvfd通过下式计算得到：

[0024] ηvfd＝50.87+1.263ω-0.0142ω2+5.834×10-5ω3

[0025] 其中：ω为再生水提升泵转速比，％。

[0026] 进一步，所述步骤1)中，非设计工况下水源侧输配系统运行方法包括：阀门节流调节和水泵变速调节。

[0027] 进一步，所述步骤2)中，地辐管散热量Q具体表达式如下：

[0028] Q＝Qu+Qd＝Qc+Qr＝Qdemand

[0029] 其中：Qu为地辐管向上的传热量，W；Qd为地辐管向下的热损失，W；Qc为地面的对流散热量，W；Qr为地面的辐射散热量，W；Qdemand为用户末端负荷需求，W。

[0030] 进一步，所述步骤2)中，用户末端负荷需求Qdemand计算式具体为：

[0031] Qdemand＝cpmu(tL,in-tL,out)

[0032] 其中：cp为水的定压比热容；mu为用户侧循环水流量，m3/s；tL,in、tL,out分别为用户侧供水、回水温度，℃。

[0033] 进一步，所述的地辐管向上的传热量Qu和地辐管向下的热损失Qd分别通过下式计算得到：

[0034]

[0035]

[0036] 其中：F为地板面积，m2；tr为室内空气温度，℃；Ku为地面向上传热系数，W/(m2·K)；Kd为地面向下传热系数，W/(m2·K)；

[0037] 地面的对流散热量Qc和辐射散热量Qr分别通过下式计算得到：

[0038] Qc＝2.17F(tf-tr)1.31

[0039] Qr＝5×10-8F[(tf+273)4-(UMRT+273)4]

[0040] 其中：tf为地板表面温度，℃；UMRT为非加热面的平均辐射温度，℃。

[0041] 进一步，所述步骤3)中，热泵机组制热量Qh、热泵机组输入功率Powerh及修正系数ηQ-t、ηP-t分别为：

[0042] Qh＝Qh,ref·ηQ-t·PLR

[0043] Powerh＝Powerh,ref·ηP-t·ηPLR

[0044]

[0045]

[0046] 其中：Qh,ref为热泵机组额定制热量，W；Powerh,ref为热泵机组额定输入功率，W；ηQ-t、ηP-t分别为热泵机组制热量、热泵机组输入功率关于机组进出口温度变化的修正系数；
PLR为热泵部分负荷容量；ηPLR为不同PLR下的功率修正系数；C1—C6、D1—D6均为拟合系数；

[0047] 蒸发器换热量Qe关系式为：

[0048] Qe＝Qh-Powerh。

[0049] 进一步，所述步骤4)中，水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe之间关系为：

[0050] Qre＝Qmc＝Qe

[0051] 水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe计算公式分别为：

[0052] Qre＝mGcp(t1,in-t1,out)

[0053] Qmc＝kexFexΔt

[0054] Qe＝mZcp(tE,in-tE,out)

[0055] 其中：Fex、kex分别为宽流道换热器换热面积和传热系数，m2、kW/m2·℃；t1,in、t1,out分别为水源侧供水、回水温度，℃；Δt为对数平均温差，Δt＝(Δt'-Δt”)/(lnΔt'/Δt”)，Δt'＝t1,in-tE,in，Δt”＝t1,out-tE,out，℃；tE,out为蒸发器出口温度，℃；mG、mZ分别为水源侧供水、中介水流量，m3/s。

[0056] 进一步，所述步骤5)中，再生水源热泵供暖系统协同优化模型具体为：

[0057] minP(x)＝Pheatpump(x,tE,in)+Ppump,re(x,mG)+P′pump

[0058]

[0059] 式中：Pheatpump(x,tE,in)、Ppump,re(x,mG)分别为负荷率x时对应的热泵机组能耗和再生水提升泵能耗，P′pump为定频水泵能耗，包含中介水循环泵Ppump,mc、用户侧循环泵Ppump,user，两者通过台数调节来满足流量需求，某一负荷率x下的P′pump可视为定值，k为热泵机组额定制热量修正系数。

[0060] 进一步，所述步骤7)中，再生水提升泵样本参数包括再生水提升泵在某一流量下对应的扬程及效率；热泵机组参数包括某一水源侧进出口温度与对应的制热量及能耗。

[0061] 所述步骤7)中，拟合出相应关系式包括：

[0062] 再生水提升泵扬程与水源侧供水流量关系式

[0063] 再生水提升泵效率与水源侧供水流量关系式

[0064] 其中，a1、a2、a3、b1、b2、b3均为拟合系数。

[0065] 本发明的再生水源热泵供暖系统协同优化运行方法的优点在于：该方法对热泵系统整体能耗做出协同优化，而非现有技术中仅针对某一子系统进行能耗降低；同时，发现了影响热泵系统能耗的协同因子——蒸发器进口温度，基于此建立了再生水源热泵供暖系统协同优化模型，通过改变蒸发器进口温度值，就可实现不同负荷率下热泵系统的协同优化，大大降低了现有技术中对热泵系统协同优化控制的难度。因此，该协同优化运行方法能很好满足实际工程的需要，有效减少热泵系统能耗，节省运行费用，便于管理人员操作，具有很好的应用前景。

附图说明

[0066] 图1为再生水源热泵供暖系统协同优化运行方法流程图；

[0067] 图2为再生水源热泵供暖系统示意图；

[0068] 图3不同负荷率下确定协同因子流程图。

[0069] 图中的标记分别表示：1.再生水管道；2.再生水提升泵；3.管道调节阀；4.再生水换热器；5.中介水循环泵；6.用户侧循环泵；7.用户末端。

具体实施方式

[0070] 下面结合附图及具体实施例对本发明做进一步的说明。

[0071] 如图1所示，本发明再生水源热泵供暖系统协同优化运行方法，包括下述步骤：

[0072] 1)根据水源侧供水流量mG，以再生水提升泵能耗Ppump,re最低为目标，确定非设计工况下水源侧输配系统运行方法，得到水源侧供水流量mG与再生水提升泵能耗Ppump,re的关系式Ppump,re＝f1(mG)。

[0073] 再生水提升泵能耗Ppump,re具体计算式如下：

[0074]

[0075] 其中：HP为再生水提升泵扬程，m；γ为流体的容重，kN/m3；ηp为再生水提升泵效率，％；ηm为电动机效率，％；ηvfd为变频器效率，％。

[0076] 电动机效率ηm通过下式计算得到：

[0077] ηm＝94.187(1-e-0.0904ω)

[0078] 变频器效率ηvfd通过下式计算得到：

[0079] ηvfd＝50.87+1.263ω-0.0142ω2+5.834×10-5ω3

[0080] 其中：ω为再生水提升泵转速比，％；

[0081] 非设计工况下水源侧输配系统运行方法包括：阀门节流调节、水泵变速调节。

[0082] 假设水源侧输配系统采用阀门节流调节，再生水提升泵工作点为A，A点再生水提升泵扬程Hp,A通过下式计算得到：

[0083]

[0084] A点再生水提升泵效率ηp,A通过下式计算得到：

[0085]

[0086] 其中：Gpr,A为A点单台水泵额定转速下的流量，m3/h， nA为A点再生水提升泵运行台数，台；a1～b3为单台水泵的设备性能常数，a1＜0，b1＜0，a3＞0；

[0087] 阀门节流调节的再生水提升泵转速比ω＝100％；

[0088] 假设水源侧输配系统采用水泵变频调节，再生水提升泵工作点为B，B点再生水提升泵扬程Hp,B通过输配系统管网阻力表达式计算得到：

[0089]

[0090] 其中：H0为输配系统静扬程；S为输配系统管道阻抗，s2/m5；

[0091] B点再生水提升泵效率ηp,B等于其相似工况点B'效率ηp,B'，ηp,B'通过下式计算得到：

[0092]

[0093] 其中：Gpr,B'为B点单台水泵额定转速下的流量，m3/h， nB为B点再生水提升泵运行台数，台；mG,B'为相似工况点B'的流量，m3/h；

[0094] 水泵变频调节的再生水提升泵转速比；

[0095] 2)根据地辐管散热量Q表达式和用户末端负荷需求Qdemand计算式，得出不同用户末端负荷需求Qdemand下的用户侧供水温度tL,in。

[0096] 地辐管散热量Q具体表达式如下：

[0097] Q＝Qu+Qd＝Qc+Qr＝Qdemand

[0098] 其中：Qu为地辐管向上的传热量，W；Qd为地辐管向下的热损失，W；Qc为地面的对流散热量，W；Qr为地面的辐射散热量，W；Qdemand为用户末端负荷需求，W。

[0099] 地辐管向上的传热量Qu和地辐管向下的热损失Qd分别通过下式计算得到：

[0100]

[0101]

[0102] 其中：F为地板面积，m2；tr为室内空气温度，℃；Ku为地面向上传热系数，W/(m2·K)；Kd为地面向下传热系数，W/(m2·K)；tL,out为用户侧回水温度，℃；

[0103] 地面的对流散热量Qc和辐射散热量Qr分别通过下式计算得到：

[0104] Qc＝2.17F(tf-tr)1.31

[0105] Qr＝5×10-8F[(tf+273)4-(UMRT+273)4]

[0106] 其中：tf为地板表面温度，℃；UMRT为非加热面的平均辐射温度，℃；

[0107] 用户末端负荷需求Qdemand计算式具体为：

[0108] Qdemand＝cpmu(tL,in-tL,out)

[0109] 其中：cp为水的定压比热容，cp＝4.18kJ/kg℃；mu为用户侧循环水流量，m3/s；tL,in、tL,out分别为用户侧供水、回水温度，℃。

[0110] 3)根据热泵机组制热量Qh表达式、热泵机组输入功率Powerh表达式、修正系数ηQ-t、ηP-t表达式及蒸发器换热量Qe与热泵机组制热量Qh关系式，得出不同负荷率x下热泵机组的蒸发器换热量Qe与蒸发器进口温度tE,in的关系式Qe＝f2(tE,in)。

[0111] 热泵机组制热量Qh表达式、热泵机组输入功率Powerh表达式及修正系数ηQ-t、ηP-t表达式分别为：

[0112] Qh＝Qh,ref·ηQ-t·PLR

[0113] Powerh＝Powerh,ref·ηP-t·ηPLR

[0114]

[0115]

[0116] 其中：Qh,ref为热泵机组额定制热量，W；Powerh,ref为热泵机组额定输入功率，W；ηQ-t、ηP-t分别为热泵机组制热量、热泵机组输入功率关于机组进出口温度变化的修正系数；
PLR为热泵部分负荷容量；ηPLR为不同PLR下的功率修正系数。C1—C6、D1—D6均为拟合系数；

[0117] 蒸发器换热量Qe关系式为：

[0118] Qe＝Qh-Powerh。

[0119] 4)根据水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe之间关系及计算公式，得出水源侧供水流量mG与蒸发器进口温度tE,in的关系式mG＝3(tE,in)。

[0120] 水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe之间关系为：

[0121] Qre＝Qmc＝Qe

[0122] 水源侧换热量Qre、中介水换热量Qmc及蒸发器换热量Qe计算公式分别为：

[0123] Qre＝mGcp(t1,in-t1,out)

[0124] Qmc＝kexFexΔt

[0125] Qe＝mZcp(tE,in-tE,out)

[0126] 其中：Fex、kex分别为宽流道换热器换热面积和传热系数，m2、kW/m2·℃；t1,in、t1,out分别为水源侧供水、回水温度，℃；Δt为对数平均温差，Δt＝(Δt'-Δt”)/(lnΔt'/Δt”)，Δt'＝t1,in-tE,in，Δt”＝t1,out-tE,out，℃；tE,out为蒸发器出口温度，℃；mG、mZ分别为水源侧供水、中介水流量，m3/s。

[0127] 5)根据步骤1)到步骤4)得出的关系式Ppump,re＝f1(mG)、Qe＝f2(tE,in)、mG＝f3(tE,in)，以不同负荷率x下再生水源热泵供暖系统能耗P最低为目标函数，建立再生水源热泵供暖系统协同优化模型。

[0128] 再生水源热泵供暖系统协同优化模型具体为：

[0129] minP(x)＝Pheatpump(x,tE,in)+Ppump,re(x,mG)+P′pump

[0130]

[0131] 式中：Pheatpump(x,tE,in)、Ppump,re(x,mG)分别为负荷率x时对应的热泵机组能耗和再生水提升泵能耗，P′pump为定频水泵能耗，包含中介水循环泵Ppump,mc、用户侧循环泵Ppump,user，两者通过台数调节来满足流量需求，某一负荷率x下的P′pump可视为定值，k为热泵机组额定制热量修正系数。

[0132] 6)根据再生水源热泵供暖系统协同优化模型，确定协同因子。

[0133] 确定的为协同因子为蒸发器进口温度tE,in，原因如下：在确定的用户末端负荷需求Qdemand下，热泵机组制热量Qh一定，随着蒸发器进口温度tE,in的改变，蒸发器换热量Qe随之变化，进而影响热泵机组输入功率Powerh；同时，在水源侧供水温度t1,in一定的条件下，不同的蒸发器换热量Qe导致水源侧供水流量mG不同，而水源侧供水流量mG又直接影响再生水提升泵能耗Ppump,re，因此，蒸发器进口温度tE,in是再生水源热泵供暖系统的协同因子。

[0134] 7)根据再生水提升泵样本参数，拟合出再生水提升泵相应关系式，然后根据再生水提升泵相应关系式结合热泵机组参数，代入再生水源热泵供暖系统协同优化模型，确定协同因子数值，由协同因子数值得出水源侧输配系统再生水提升泵的开启台数和变频大小，及热泵机组的开启台数和运行负荷大小。

[0135] 其中，再生水提升泵样本参数包括：再生水提升泵在某一流量下对应的扬程及效率，热泵机组参数包括某一水源侧进出口温度与对应的制热量及能耗。

[0136] 拟合出相应关系式包括：

[0137] 再生水提升泵扬程与水源侧供水流量关系式

[0138] 再生水提升泵效率与水源侧供水流量关系式

[0139] 热泵机组制热量关于机组进出口温度变化的修正系数ηQ-t、热泵机组输入功率关于机组进出口温度变化的修正系数ηP-t，其中a1-a3、b1-b3均为拟合系数。

[0140] 热泵机组参数还包括：热泵机组额定制热量Qh,ref、热泵机组额定输入功率Powerh,ref、中介水流量mZ、宽流道换热器面积Fex、宽流道换热器传热系数kex、水源侧供水温度t1,in、地板面积F、室内空气温度tr、非加热面的平均辐射温度UMRT、地面向上传热系数Ku、地面向下传热系数Kd。

[0141] 水源侧输配系统的运行参数包括：水源侧供水流量mG、水源侧回水温度t1,out。

[0142] 热泵机组的运行参数为：蒸发器出口温度tE,out。

[0143] 下面通过具体实施例对本发明方法做进一步详细说明。

[0144] 实施例：某再生水源热泵供暖系统协同优化运行方法。

[0145] 某再生水源热泵供暖系统以污水处理厂排出的二级出水为热源，由取水点经再生水管道1流入蓄水池内，再经再生水提升泵2加压后，经管道调节阀3二级出水流入各个热泵机房的宽流道再生水换热器4内，并与中介水进行换热，换热后的二级出水再经管道流至退水点；中介水通过中介水循环泵5，将热量传递给热泵机组，热泵机组通过用户侧循环泵6，将热量传递给用户末端7，如此循环，具体见图2。用户末端采用低温地板辐射采暖。

[0146] 项目冬季供暖负荷需求为70MW，建筑面积为200万m2，建筑热负荷指标为35W/m2，统一设置提升泵房，分设多个热泵机房。

[0147] 水源侧输配系统管网阻力表达式为：H＝10+5.536×10-7Q2。

[0148] 由于再生水需求量大，单台水泵无法满足其流量需求；同时，考虑到实际用户需求基本在部分流量下，为了节约能耗，便于输配系统在不同工况下运行，采用不同型号的水泵。因此，水源侧的输配系统的水泵配置由五台同型号大流量水泵并联及三台同型号小流量水泵并联构成，大流量再生水提升泵型号为WQ1700-50-355，小流量再生水提升泵型号为WQ800-38-132，根据再生水提升泵的详细样本参数，拟合出再生水提升泵的设备性能常数，得到单台大流量再生水提升泵扬程及水泵效率关系式分别为：

[0149]

[0150]

[0151] 单台小流量再生水提升泵扬程及水泵效率关系式分别为：

[0152]

[0153]

[0154] 结合热泵机组参数，通过分析计算，得到水源侧供水流量mG与再生水提升泵能耗Ppump,re的关系式为：

[0155] Ppump,re＝56.18333-0.02616mG+2.48034×10-5mG2

[0156] 根据地辐管散热量Q表达式和用户末端负荷需求Qdemand计算式，得出不同用户末端负荷需求Qdemand下的用户侧供水温度tL,in。考虑到在地面构造时会铺设一层绝热层，因此，忽略地辐管向下的热损失Qd，假设地辐管的散热量Q等于向上的传热量Qu；由于能量守恒定律，地辐管传给室内空气的热量等于地辐管传给地板表面的热量，因此，Q又可表示为：

[0157]

[0158] 其中：K'u为填充层以上传热系数，W/(m2·K)，K'u＝1.40W/(m2·K)。

[0159] 以该再生水源热泵供暖系统的某一片区的热泵机房为例进行计算，已知在满负荷情况下用户末端的负荷需求Qdemand＝2820kW，由此得到不同负荷率x下的用户末端负荷需求；同时，在不同负荷率x下用户侧供水流量mu也有所不同，具体见表1。

[0160] 表1不同负荷率下的热负荷及用户侧供水流量

[0161]

[0162] 已知参数为：F＝8×105m2、tr＝16℃、UMRT＝15.6℃、cp＝4.18kJ/kg℃；同时，满足规范要求：用户侧供水温度tL,in≥35℃。结合已知参数计算得出结果见表2。

[0163] 表2不同负荷率下的用户侧供水温度

[0164]

[0165] 该机房采用三台热泵机组运行，热泵机组型号为KLSH-270S-F。根据热泵机组样本参数，拟合出ηQ-t、ηP-t分别为：

[0166]

[0167]

[0168] 同时，根据样本数据拟合出ηPLR为：ηPLR＝0.13321+0.64618PLR+0.19528PLR2。

[0169] 根据Qh、Powerh、ηQ-t、ηP-t表达式，Qe与Qh关系式，得出不同负荷率x下蒸发器换热量Qe与蒸发器进口温度tE,in的关系式Qe＝f2(tE,in)。

[0170] 已知参数为：Qh,ref＝1142kW、Powerh,ref＝235kW。

[0171] 假设在负荷率x∈[100％,60％)时，选择三台热泵机组运行，中介水流量mZ＝426m3/h；在负荷率x∈[60％,30％)时，选择两台热泵机组运行，mZ＝284m3/h；在负荷率x∈[30％,0％)时，选择一台热泵机组运行，mZ＝142m3/h。

[0172] 结合已知参数，得出不同负荷率x下蒸发器换热量Qe与蒸发器进口温度tE,in的关系式。具体见表3。

[0173] 表3不同负荷率下的单台热泵机组蒸发器换热量

[0174]

[0175] 根据Qre、Qmc及Qe之间关系及计算公式，得出水源侧供水流量mG与蒸发器进口温度tE,in的关系式。

[0176] 已知参数包括：Fex＝200m2、kex＝1.2kW/m2·℃、t1,in＝16℃。

[0177] 结合已知参数，得出不同负荷率x下，水源侧供水流量mG与蒸发器进口温度tE,in的关系mG＝f3(tE,in)。具体见表4。

[0178] 表4不同负荷率下的水源侧供水流量

[0179]

[0180]

[0181] 结合实际工程，得出再生水源热泵供暖系统协同优化模型具体为：

[0182] minP(x)＝Pheatpump(x,tE,in)+Ppump,re(x,mG)+P′pump

[0183]

[0184] 其中：min mG,L＝mG,Lref×50％、min mG,S＝mG,Sref×50％。mG,L为采用大流量水泵时的水源侧供水流量，mG,S为采用小流量水泵时的水源侧供水流量；mG,Lref为大流量水泵的额定流量值，m3/h；mG,Sref为小流量水泵的额定流量值，m3/h。

[0185] 根据再生水源热泵供暖系统协同优化模型进行协同优化计算，首先计算出不同负荷率x下的协同因子值，具体步骤见图3，在根据协同因子值得出不同负荷率x下的温度最优解，具体数值见表5。

[0186] 表5不同负荷率下热泵系统的温度最优解

[0187]

[0188]

[0189] 水源侧输配系统仅设置一个提升泵房，而该片区热负荷占总项目热负荷的4％，因此，以该片区水源侧需求量为基准，计算出水源侧输配系统的供水流量。不同负荷率x下的水源侧供水流量优化值见表6。

[0190] 表6不同负荷率下输配系统的最优流量值

[0191]

[0192] 这些优化解在实际应用中对应着设备的运行方法，因此，给出不同负荷率x下的再生水提升泵及热泵机组运行方法，如表7所示。

[0193] 表7不同负荷率下的设备运行方法

[0194]

[0195] 为体现协同优化方法的节能效果，分别采用两种现有技术作为协同优化方法的比较对象，在相同的约束条件下，进行不同负荷率x下热泵系统能耗的计算。

[0196] 现有技术一：在满足用户负荷需求的条件下，对输配系统进行运行优化，即以输配系统能耗最低为目标，对水源侧输配系统进行非设计工况下运行调节方法的优化；而热泵机组按设计工况运行，即蒸发器进口温度按照设计值运行。

[0197] 现有技术二：在满足用户负荷需求的条件下，对热泵机组进行运行优化，即以热泵机组能耗最低为目标，优化蒸发器进口温度值；而输配系统的再生水提升泵采用定频方式运行，即仅通过改变水泵台数来满足流量需求。

[0198] 经计算后，得出本方法与现有技术的对比，具体见表8。

[0199] 表8本方法与现有技术对比

[0200]

[0201] 由表8可以看出，在负荷率较低的条件下，采用本方法有更好的节能效果(最高达到32.68％)；在设计负荷需求附近即负荷率在90％～100％范围内，采用本方法也有一定的节能空间(2％左右)。

[0202] 上述实例只为说明本发明的技术构思和特点，并不能以此限制本发明的保护范围。对于本领域的技术人员来说，凡是根据本发明精神实质所做的等效变换或修饰改进，都应涵盖在本发明的保护范围之内。